lunes, 9 de septiembre de 2024

Problemas seleccionados de cálculo Integral infinita y repaso de la función W de Lambert

 

Medellín, octubre 2024

 

Problemas seleccionados de cálculo

Integral infinita y repaso de la función W de Lambert

 

1.    1. Integral de un integrando infinito


2.    Algunas funciones discretas de x

Función y = x= piso(x); y = x = cielo(x), función parte fraccionaria o decimal de x; y = parte decimal (x) = {x}

Función y = piso(x) = x            y es igual al menor entero mayor que x

Ejemplo            y = piso (2,8) = 2,         y = Piso(-101,3) = -102    y = piso(-0,7) = -1


Figura 1    Función Piso de x

Función y = techo(x) = cielo(x) = ceiling(x) =x         y es igual al mayor entero mayor que x


 Figura 2, función techo de x (ceiling(x)

 

3,4= 4

 -2,4= -2

 -0,7= 0

 -100,6= -100

 Función parte decimal     y = {x} = x - x

Ej. {3,5} = 3,5 – piso(3,5) = 3,5 – 3 = 0,5

{-3,6} = -3,6 - piso(-3,6) = -3,6 –(- 4) = 0,4

{-1,3} = -1,3- piso(-1,3) = -1,3 –(- 2) = 0,7

Resolvamos una ecuación que involucre una de estas 3 funciones

x/x + 2x = 8                                          (1)   Alto grado de dificultad!

 Lo primero que observamos es que x⏌ debe ser  ≠0       y   además    x (0, 1)

Veamos qué pasa cuando x<0

Si x es negativo, k<= x <= k+1   y la ecuación (1) se convertirá en (k+1)/k + 2x = 8

(k +1/k) > 0 y tiende a 1, cuando x tiende a -∞  (Es positivo, porque ambos k +1 y k son negativos.

Veamos qué pasa con la ecuación (1) cuando x = -4,3      x/x = (k+1)/k =

( -4+1)/(-5) = (-3)/(-5) = 3/5

La ecuación queda 2x = 8 – 3/5 = 37/5         x = 3,7; lo que es absurdo, porque habíamos partido de x = -4,3

Por lo anterior excluimos todos los reales negativos, de la solución de x.

Para encontrar una solución, suponemos que x es un real positivo >1, tal que

k <= x <= k + 1 y la ecuación original         (2)

(k+1)/k + 2x = 8     de la cual deducimos que x = (7k -1) / (2k)         (3)

(3) en (2)

k (7k -1) / (2k) k +1                (4

La ecuación (3) nos provee dos desigualdades, que deben ser satisfechas por los valores de k que sean solución

2k2 – 7k +1 0                      (5)        y            7k – 1 - 2k2 -2k 0                              (6)

K2 – (7/2) k + ½ 0            (5a)       y                k2 –(5/2)k + ½ 0                              (6a)

Resolviendo y factorizando las ecuaciones (5a) y (6a) con la ayuda de Derive, obtenemos:

La (5a)     (k – 0,142)(k – 3,3508)       que es negativa para k [0,142, 3,3508]

La (6a)     (k – 0,2192)(k – 2,2808)     que es positiva para k[-∞, 0,2192]  cuyo intervalo no cumple las condiciones identificadas anteriormente.

Igualmente (6a) es positiva en el intervalo k (2,2808, ∞)

Los enteros que son intersección de soluciones a (5a) y (6a) son la solución para k

Sólo k = 3 sirve como solución. Ahora encontremos x en la ecuación original

4/3 + 2x = 8                      x = (8 – 1,33333…)/2        x = 3,3333333.. = 10/3

Remplacemos en la ecuación original:

4/3 + 2(10/3) = 24/3 = 8

 

3.    Problema de solución de una ecuación, aplicando la función de Lambert.

 xx =12

Refresquemos el concepto de la función de W de Lambert. (Hay varios blogs del autor Maths, dedicados a esta función)

La función xex tiene un dominio (-∞, ∞) = Reales     rango = [-1/e, ∞)

No es biyectiva, en el sentido que para dos valores diferentes de x, podemos obtener un mismo valor de y


Figura 3           gráfico de f(x) = xex

Estrictamente f(x) no tiene una inversa única, pero dada la importancia de la función se ha obtenido una inversa para la parte decreciente, que hemos llamado W-1y otra para la parte creciente, que hemos llamado Wo

Es obvio que no es posible encontrar fórmulas algebraicas para las funciones Wo y W-1, ya que no es posible despejar x en la ecuación y = xex

 


Figura 4, funciones Wo y W-1

 

Limitemonos a la rama principal de Wo

 

. Cuyo dominio va des [-1/e, ) y cuyo rango va de [-1, )

 

Aunque la inversa de xex no tiene fórmula, pero hay algoritmos que nos permiten encontrar valores de específicos de W.

 

Si tenemos un valor específico de xex, digamos 12,5 podemos encontrar el valor de x que produjo ese valor, simplemente hallando W(12,5) = 1.889445

 

Lo obtuve de una calculadora libre en la red, que tiene la función W y se encuentra en el siguiente link:

 

https://www.had2know.org/academics/lambert-w-function-calculator.html

 

Problema:

 

Resolver la ecuación xx = 12 , utilizando el método de Lambert

 

xx = 12      x lnx = ln12     todavía no se parece a Lambert

 

Cambiando x por e lnx

 

elnx lnx = lnx e lnx = ln12;  esto si se parece a la aplicación de Lambert

 

Sea u = lnx                        ueu = ln12          y

 

W(ueu) = u = W(ln12)

 

Ln12 = 2,484991

 

W(2,484991) = u     Cómo hallamos el valor de W(2,484991)? R_ Con el link que indiqué anteriormente:

 

u = 0,9556415 = lnx

 

e 0,9556415 = x        x = 2,600338

 

Verifiquemos la respuesta:

 

2,600338 2,600338 = 12

 

Otro ejemplo

 

Resolver por medio de la inversa W de Lambert, la ecuación x + ex = 5

 

Para que se parezca a una ecuación que se pueda resolver por Lambert hacemos los siguientes cambios.

ex = 5 - x

1= (5 .-  x) e -x

 e5 = (5 – x) e -x e5

(5 – x) e5 - x = e5

u = 5 - x

ueu = e5

W(ueu)  = u = W(e5)

W(e5) = 3,693441      Con el link sugerido.

u = 5 – x = 3,693441                x = 5 - 3,693441 = 1,306559

Reemplacemos en la ecuación:

1,306559 + e 1,306559 = 5,000001      la millonésima se debe a los decimales, pero en realidad no existe,

 

Juan Fernando Sanin E

juanfernando.sanin@gmail.com

 





 

viernes, 12 de julio de 2024

Curvatura

Jardín, octubre de 2024

 

Curvatura en 2D y 3D

 

Definición

Si en el plano 2D, tenemos una curva C, tal como se indica en la figura 1; sea P un punto de esa curva, cuya tangente hace un ángulo θ (radianes) con el eje x y et el vector tangente unitario a C, en P, se define la curvatura de C en P así:



Figura 1

k(x) = det/ds                 curvatura de C en el punto P       (1)

La curvatura escalar será:                          k = │k(x)│= │det/ds│

Si R(t) es el vector posición de un punto P, que recorre una curva C (2D o 3D), en el caso de que la variable t sea el tiempo, R’(t) es el vector velocidad del punto P y este vector tiene la dirección de et (vector tangente unitario), que ahora estará en un plano, que no tiene que ser plano xy. Ese plano se llama plano osculador. Más adelante describiremos como se determina ese plano.

La definición se repite para 3D:         k(t) = det/ds     et y ds en el plano osculador.





Figura 2

k = det/ds       Curvatura de una curva 2D o 3D, en un punto P es la razón de cambio del vector tangente unitario et, respecto del cambio en el arco s

k = det/ds = (det/dt) (dt/ds)                                                                  (2)

Si t fuese la variable tiempo, ds/dt sería la rapidez v de la partícula, entonces, ds/dt = v, y es igual a R’(t).       Aunque la variable t no sea el tiempo, esta relación se cumple, siendo t cualquier variable real.

ds/dt = v = │R’(t)│                                                                              (3)

La (2) se convierte en         k = det/ds = (det/dt) (dt/ds) = (det/dt) / │R’(t)│     (4)

que ya es una definición práctica, para curvatura vectorial o escalar, según se aplique.

La fórmula (4) se puede convertir en la (5) con un poquito de algebra t combinando las todas las relaciones anteriores. El libro de Cálculo de Leithold, 7ª edición, coloca esta demostración como un ejercicio para resolver.

k=│[dR(t)/dt) X d2R(t)/dt2]│ / │dR(t)/dt│3                                                           (5)

Vector normal

Observemos la figura 3-2 y 3-1 en este orden:



Figura 3

Si el vector et hace un ángulo θ con el eje x, los vectores et y e’t hacen un ángulo Δθ. Los vectores et y e’t son los vectores tangentes unitario en los puntos P y P’.

Si llevamos estos vectores a un sistema cartesiano con origen O’, ambos tienen radio 1 y hacen un ángulo Δθ. Figura 3-1

El vector Δet es tal que                                                     et + Δet = e’t

Y │Δet │= 2 (1) sen (Δθ/2) =

Y │Δet │/ Δθ = sen (Δθ/2) / (Δθ/2)

Cuando Δθ→0, entonces [│Δet │/ Δθ = sen (Δθ/2) / (Δθ/2)] →1

Δet es un vector unitario y perpendicular a et (Sólo en el plano x, y)

Por lo tanto, det │/ dθ; es un vector unitario normal, perpendicular a et    (6       en = det/dθ)

Veamos que es la curvatura, en forma práctica, en el plano x y

La parametrizamos así:

y= f(x)                                               (7)

x = t                                                  (7a)

f(x) = f(t)                                           (7b)  

R(t) = (t, f(t), 0)                                (8)

R’(t) = (1, f’(t), 0)                              (9)

R”(t)= (0, f”(t), 0)                               (10)

et = [(1, f’(t), 0)]/√(1 + (f(t))2)              (11)

Continuamos con la curvatura en y = f(x)



Recordemos que x = t

Toda curva en el plano o en el espacio, en un punto dado, tiene curvatura, siempre y cuando, en ese punto P(x, y, z), definido por el vector posición R(t), existan  R’(t) y R”(t).

En los puntos de inflexión de esas curvas, la curvatura es igual a 0

En un punto donde hay curvatura, existen infinitas circunferencias tangentes. La pregunta es: ¿Cuál de esas circunferencias es la que mejor se adapta a la curva y a la curvatura en ese punto?

Para responder esto, busquemos la curvatura en una circunferencia

x2 + y2 = r2

Utilicemos sólo la parte superior           y = √ (r2 -x2)

y’ = (1/2) (-2x) / √ (r2 -x2) = -x/√ (r2 -x2)

y” = - [ √ (r2 -x2) (1) – x ((1/2) (-2x) / √ (r2 -x2)]/ (r2 -x2)

   

 = r2 / (r2 – x2) 3/2

k(x) = [ r2 / (r2 - x2) 3/2] / [1 + x2/ (r2 – x2)] 3/2

k(x) = r2 / (r2 )3/2= r2 / r = 1/r

En un punto dado de cualquier circunferencia, la curvatura es igual a 1/Radio. Todos los puntos de la misma circunferencia tienen la misma curvatura y en el caso de varias circunferencias a mayor curvatura menor radio y viceversa.

Volviendo al caso de una curva 2D o 3D, en un punto P, donde exista curvatura, se ha definido que la circunferencia tangente a la curva en P, que mejor se acomoda a la curva es aquella cuya radio sea 1/k y su centro esté dentro de la curva y tenga la dirección del vector normal.

Por tanto, no es gratuito afirmar que la circunferencia de curvatura, de una curva, en un punto dado P, es aquella, que contiene ese punto, que el radio sigue la dirección interior del vector normal y cuyo radio es

ρ= 1/k

Si la curva es en el plano x, y, el plano de la circunferencia es el x, y.

Si la curva es en 3D, la circunferencia está en el plano osculador, definido así:

1.Contiene el punto P de la curva,

2.Su normal es el vector binormal eb = et x en

A veces es difícil y laborioso hallar et y en

Por tanto, el plano osculador, queda definido por el punto P(xo, yo, zo) y el vector:

et x en


Ejercicio 1

R(t) = (t-1, 2ln(t), 2t)                               (1)

El punto P queda en t=2

Hallar el plano osculador en ese punto, la curvatura escalar y el vector curvatura en ese punto

Para t = 2, el punto P es P(1/2, 2ln2, 4)

R’ = (-t-2, 2t-1, 2)                │R’│ = (2t2 + 1)     (2)     Se deja como ejercicio para el lector

et = (-1/(1 + 2t2), 2t/(1 +2 t2), 2t2/(1 +2 t2))

El punto P para t = 2        et = (-1/9, 4/9, 8/9)     (3)               Verificar que es unitario

R” = (2t-3, -2t-2, 0)           

N(t) = det/dt = [d(-1/(1 + t2), 2t/(1 + t2), 2t2/(1 + t2)) / dt]/        Que es muy laborioso. Lo hice por

derive.                             en = N(t) /│N(t)│

N(t) = (1/(2t2 + 1)2)(4t, 2 – 4t2, 4t)

en= (1/2)(1/(2t2 + 1)2)(4t, 2 – 4t2, 4t)

Para comprobar que et y en son perpendiculares hacemos el producto escalar y se ve que nos da igual a 0

En t = 2             en =(1/18) (8, -14, 8)

El vector normal, que define el plano osculador será:

(-1, 4, 8) X (8, -14, 8)        Podemos ignorar el denominador 9 en et  y coeficiente que precede a en

 (80, -72, -18)                           el vector eb =  (80, -72, -18)/ √(802 +722+182) = (80, -72, -18/2√2977

La ecuación del plano osculador es 80x -72y -18z + C = 0

Aplicado a P(1/2, 2ln2, 4)

40 -144ln2 -112 + C = 0

C = 72 + 144ln2

80x – 72y -18z +144ln2 = 0                   (4)           Ecuación del plano osculador.

Para la curvatura escalar k utilizaremos la fórmula (15) para t = 2  formulas (2) y (3) del ejercicio

En t = 2 R’ y R” son, respectivamente, los vectores:

R’ = (-1/4, 1, 2)

R” = (1/4, -1/2,0)

El producto cruz o vectorial de R’ X R” es el vector          (1, ½, -1/8)

Cuya magnitud es 9/8

La magnitud del vector R’ es 9/4

k= (9/8)/(9/4)3 = 8/81

Y el radio del círculo osculador será ρ = 81/8

El vector curvatura k = (8/1458) (8, -14, 8)

 

Juan Fernando Sanin E

juanfernando.sanin@gmail.com


lunes, 8 de julio de 2024

Movimiento curvilíneo de partículas

 

Jardín, septiembre de 2024

 

Movimiento curvilíneo de partículas

 

1.   1.Vector de posición, velocidad y aceleración

Si una partícula se mueve sobre una curva, diferente a una recta, entonces describe un movimiento curvilíneo. Para definir la posición P, ocupada por la partícula en un tiempo determinado t, en un sistema x, y, z que se muestran en la figura a), y se dibuja el vector r que une al origen en O y el punto P. Como el vector r define la magnitud y su dirección con respecto a los ejes referidos, éste vector r define la posición de la partícula con respecto a esos ejes; el vector r se conoce como el vector de posición de la partícula en el tiempo t.

Veamos el vector r’ que define la posición P’ ocupada por la partícula en un tiempo posterior t + Δt. El vector Δr que une a P y a P’ representa el cambio en el vector de posición durante el intervalo del tiempo Δt, pues, como se puede verificar en la figura a) el vector r’ se obtiene al sumar los vectores r y Δr. Δr representa un cambio de dirección y magnitud del vector r. La velocidad promedio de la partícula sobre este intervalo, se define como el cociente de Δr y Δt. Puesto que Δr es un vector y Δt es un escalar, el cociente de Δr/Δt es un vector unido a P, de la misma dirección que Δr y de magnitud igual a la magnitud de Δr dividida entre Δt. figura b)

La velocidad instantánea de la partícula en el tiempo t se obtiene al elegir intervalos de tiempo Δt cada vez más cortos y, de manera correspondiente, incrementos vectoriales Δr cada vez menores. La velocidad instantánea se representa en consecuencia mediante el vector:

v = dr/dt                                                                                                      (1)



          Figura 1

El ABS v del vector v se llama la rapidez de la partícula y se consigue al sustituir, en vez del vector r en la fórmula (1), la magnitud de este vector, representado por el segmento de línea recta PP’. La longitud del segmento PP’ es igual a Δs del arco PP’ en el límite, cuando Δt tienda a 0.

ABS(v)= ds/dt              rapidez de la partícula = ds/dt                                  (2)

Ahora, si el móvil se mueve a lo largo de una curva, la velocidad cambiara´ también, no sólo de rapidez, sino de dirección. La aceleración de la partícula no necesariamente, coincide con la dirección de la velocidad y por tanto no tiene porque tener la dirección de dr/dt = v

                                                                                     



Figura 2. El vector aceleración, no necesariamente, tiene la dirección de dr/dt

Con un razonamiento similar al que utilizamos con la velocidad v

a = dv/dt                                                                                            (3)

Y en muchos casos es importante trabajar con ABS(a), sin que este valor tenga un nombre especial.

2.    Derivadas de funciones vectoriales

Si P(u) es el vector de posición de una partícula en función de una variable u, que podría ser el tiempo t u otra.

Sabemos que ΔP = P (u + Δu) – P(u)

Y

dP/dt = Lim cuando Δu →0 de ΔP/Δu = Lim cuando Δu→0 de (P (u + Δu) – P(u)) /Δu    (4)

Si P y Q son funciones de u

d (P + Q) /du = Lim cuando Δu→0 de Δ (P + Q) / Δu = Lim Δu→0 de (ΔP/Δu + ΔQ/Δu)

y por tanto:

d (P + Q)/du = dP/du + dQ/du)                                                                                           (5)

Si tenemos el vector fP, donde f es una función escalar de u y P una función vectorial de u

Sea el vector de posición P(t) = f(t) i + g(t) j + h(t) k = (f(t), g(t), h(t))

Derivemos el vector f(f) i

d f(t) i/ d(t) = f(t) d i/dt + if’(t) = f’(t) i           ya que di/dt = 0 por ser i un vector constante, por lo que:

P’(t) = (f’(t), g’(t), h’(t))                                                          y

P”(t) = (f”(t), g”(t), h”(t))                                                         

Si R(t) es el vector de posición de una partícula, R(t) = (f(t), g(t), h(t)), entonces:

R’(t) = (f’(t), g’(t), h’(t))         Es la velocidad de la partícula y tiene una dirección tangencial

R”(t) = (f”(t), g”(t), h”(t))       Es la aceleración de la partícula y no tiene la dirección de la velocidad o la tangente, tal como se indica en la figura 2.

Inclusive, es posible y eso lo veremos más adelante, descomponerla en una componente tangencial y otra radial; radial en la dirección del radio de curvatura, cuyo concepto también revisaremos.

 

Movimiento relativo a un sistema de referencia en translación.


 Figura 3

Si el punto A, centro del sistema que se traslada (En forma paralela a x, y, z) tiene vector de posición RA, el cual tiene velocidad y aceleración vA y aA. A su vez el punto B tiene dos vectores de posición, RB, respecto al eje fijo y RB/A, respecto al sistema que se está trasladando, se dan las siguientes ecuaciones:

RB = RA + RB/A

vB = vA + vB/A

aB = aA + aB/A

3.    UNIDADES

Sistema Internacional de Unidades (unidades del SI).

En SI, las unidades básicas son: longitud (m) metro, masa (kg) kilogramo y tiempo (t) segundos. La unidad de fuerza es una unidad derivada. Se denomina newton (N) y se define como la fuerza que produce una aceleración de 1 m/s2 a una masa de 1 kg

1 N = (1 kg) (1 m/s2) = 1 kg -m/s2

1kgf (fuerza) produce a un kg (masa), una aceleración de -9,8 m/s2 = g m/s2    g = -9,8 m/s2

En el sistema inglés:

Los matemáticos estadounidenses siguen utilizando de forma común longitud, fuerza y tiempo como lo medían los ingleses; estas unidades corresponden, respectivamente, al pie (ft), la libra (lb) y el segundo (s). El segundo es el mismo que la unidad correspondiente del SI. El pie se define como 0.3048 m. La libra se define como 0.453kg. (kg fuerza)

En este sistema, el valor normal de la gravedad es, g = 32,2 ft/s2. La unidad de masa consistente con el pie, la libra fuerza y el segundo es la masa, que recibe una aceleración de 1 ft/s2 , cuando se le aplica una fuerza de 1 lb

Esta unidad, llamada en ocasiones un slug, puede deducirse de la ecuación

F = ma, después de sustituir

1lb = 1 Slug x 1pie/s2

1 Slug = 1 lb /(1 pie/s2)     El slug es la unidad de masa en el sistema inglés. Cada que se escriba lb, se refiere a una libra fuerza.

Ejemplo 1

El automóvil A viaja hacia el oeste, con una rapidez constante de 36 km/h- Cuando el automóvil A cruza por la intersección, el automóvil B parte del reposo desde una distancia 35 m al norte y se mueve hacia el sur con una aceleración constante de 1,2 m/s2.

Determinar la posición, velocidad y aceleración relativa de B respecto de A, 5 s después que A pasa por la intersección.


Figura ejercicio 1

Se eligen los ejes x e y con el origen en la intersección de las dos calles. El sistema x’ y’ no dibujado va con el carro A.

vA = 36 i km/h     Constante                                 vBt=0 = 0        aB = -1,2 j m/s2

vA = 36 i km/h = 10 m/s       para tener unidades consistentes

aA            vA          xA       en  t = 0  ;     en  t = 5      aB       vB          yBt=0       yBt = 5

Movimiento del automóvil A.               Trabajemos con el concepto escalar de rapidez

aA = 0                                           vA  = 10 i m/s           xA = (xA)o + vA t = 0 + 10t (m)

Para t = 5

aA = 0                                           vA  = 10 i m/s           xA = vA x5 = 50 i m

Movimiento del automóvil B

aB = -1,2 j m/s2        vB = (vB)o + at = 0 – 1,2t        yB = (yB)o + (vB)ot + aBt2/2 = 35 – 0,6t2

Para t = 5

aB = -1,2 j m/s2    vB = 5(-1,2) = -6 j m/s      yB = 35 – 0,6(25) = 20 j m ;  todavía no llega al cruce.

Movimiento de B relativo a A

RB = RA + RB/A          20 j = 50 i + RB/A        RB/A = 50 i – 20 j    ABS RB/A = √(400 + 2500) = 

53,9

α= tan-1 (-50/20) = -21,8014 deg

aB = aA + aB/A          aB/A = aB + 0 = -1,2 j m/s2        

Ejercicio 2

Un proyectil se lanza con una velocidad inicial de 800 m/s, a un blanco ubicado a 2000 m de altura, por arriba del cañón y a una distancia horizontal de 12000 m. La resistencia del aire es despreciable. Calcular el valor del ángulo de disparo α.

 Figura ejercicio 2

El movimiento horizontal es uniforme ax = 0         vx= 800cosα        Uniforme

x= vx t = (800 cosα)t

El tiempo que se requiere para que el proyectil recorre x = 12000 m      t =12000/ (800 cosα) =15/cosα

El movimiento vertical tiene aceleración g.           g = -9,8 m/s2

(vy)o = 800 senα

y = (vy)ot + gt2/2 = (800 senα)t -9,8t2/2                2000 = (800 senα)t – 4,9t2

2000 = (800 senα) (15/cosα) – 4,9 (15/cosα)2 = 12000 tanα - 1102,5 sec2α

Sec2α = 1 + tan2α

2000 = 12000 tanα - 1102,5 (1 + tan2α)

3102,5 = 12000 tanα - 1102,5 tan2α    una ecuación de segundo grado para tanα

tanα1 = 0,264993         α1 = 14,84 deg

tan α2 = 10,619360       α2 = 84,62 deg

Con cualquiera de los dos ángulos, el proyectil da en el blanco.

4.    Componentes tangencial y normal

Como se vio en la figura 2 del título 1, si en un punto dado y un tiempo dado t, la partícula cambia de rapidez y de dirección, la aceleración no es tangente a la curva descrita por R(t) y se puede descomponer en componentes x, y, z o simplemente x, y si se trata de una curva plana. Pero igualmente, en 2D o 3D, a(t) se puede descomponer en una componente tangencial y otra normal a la tangente.

Si una partícula se mueve sobre una curva (2D o 3D) con rapidez constante, la aceleración sólo tiene componente radial o normal.

En un punto dado de la curva P (x, y, z) o simplemente P (x, y) existen los vectores unitarios et tangencial y en, normal y un plano, llamado por los matemáticos: plano osculador, cuya normal es el producto vectorial

et vector unitario tangencial

en vector unitario normal

eb = et X en         el vector unitario que define el plano osculador

Es demostrable y se encuentra en casi todos los libros de cálculo que:

en = det/dθ                                                                                    (1)

Las figuras, en la figura (4), ayudan a hacer las demostraciones correspondientes a et, en, v, a

En la primera hemos colocado


Figura 4

 

En la figura 1 de (4), colocamos los dos vectores unitarios et y et que forman el ángulo Δθ, en la (4) 2.  

Siendo unitarios, la curva que subtiende al vector Δet será una circunferencia de radio =1.          ABS Δet =2sen(Δθ/2)         y Lim cuando Δθ→0 de Δet/Δθ será un vector perpendicular a et

y, teniendo en cuenta que:

Lim cuando Δθ→0 de sen(Δθ)/2) /(Δθ/2) = 1

Entonces Lim cuando Δθ→0, Lim de Δet/Δθ será un vector perpendicular a et = en   y unitario

Y por tanto                                                    det/dθ = en               que es la ecuación (1)

El vector velocidad v tiene la dirección de et         y

en=det/dθ                              vector unitario normal

v = v et                                                                         (2)  v = ABS v = rapidez

Sabemos que:

a =dv/dt ;     Recordemos que    v = vet       

a = dvet/dt = vdet /dt + et dv/dt                                                       (3)

det/dt = (det/dθ) (dθ/ds)(ds/dt)

Si s es el arco, entonces         v = ds/dt                                  y              det/dθ = en

det/dt = ven /ρ        (xx)                                              

1/ρ(dθ/ds)           (yy)

  

det/dθ = en       y        1/ρ(dθ/ds)       

Nota.    Estas ecuaciones, (xx) y (yy), que vienen del cálculo diferencial e integral, ameritan una explicación adicional, que se hará en un próximo blog.  Por hora, las utilizaremos, tanto en los conceptos faltantes, como en los ejercicios.                                                                      

a = (dv/dt) et + (v2/ ρ) en                                                                    (4)

at = dv/dt                                                                                             (5)

an = v2/ρ                                                                                              (6)

La fórmula (6) es una fórmula de la física elemental y describe la aceleración centrípeta de una partícula que se mueve sobre una circunferencia de radio ρ, con rapidez uniforme.

¿Qué es el radio de curvatura?  En cada punto de una curva P(x, y, z) hay una curvatura k(t) y un radio de curvatura ρ (todavía no definidos)

1.El radio de curvatura ρ, en un punto en la curva es, a grandes rasgos, el radio del círculo que mejor se ajusta la curva es ese punto.

2.La curvatura, denotada como k, es 1 dividido entre el radio de la curvatura ρ

Respecto de las dos proposiciones anteriores se discutirán brevemente más adelante, aunque la posible respuesta, no es del todo satisfactoria.

5.Componente tangencial y normal es una curva 3D


Figura 5, tomada del libro Mecánica de Beer and Johnston, Dinámica.

 

Movimiento de una partícula en el espacio.

Las fórmulas (1) a (6) se cumplen también para una partícula que se mueve a lo largo de una curva en el espacio. Sin embargo, puesto que hay un número infinito de líneas rectas que son perpendiculares a la tangente en un punto dado P de una curva en el espacio, es necesario definir con más precisión la dirección del vector unitario en.

Se considerarán de nuevo los vectores unitarios et y e’t, tangentes a la trayectoria de la partícula en dos puntos vecinos P y P’ (figura 11.24a) y el vector Δet que representa la diferencia entre et y e’t (figura 5-a). Imagine ahora un plano que pasa por P (figura 11.24a) paralelo al plano definido por los vectores et, e’t y Δet (figura 5-b). Este plano contiene la tangente a la curva en P y es paralelo a la tangente en P’.  Si se deja que P’ se acerque a P, se obtiene en el límite el plano que mejor se ajuste a la curva en la vecindad de P.

Este plano recibe el nombre de plano osculador en P. De esta definición se deduce que el plano osculador contiene al vector unitario en, puesto que este vector representa el límite del vector Δet/Δθ.

La normal definida por en está contenida en el plano osculador; ésta recibe el nombre de normal principal en P. El vector unitario eb = et X en, que completa la tríada derecha et, en, eb (figura 5-c) define la binormal en P. En consecuencia, la binormal es perpendicular al plano osculador.

Por tanto, la aceleración de la partícula en P puede descomponerse en dos componentes, una a lo largo de la tangente, y la otra a lo largo de la normal principal en P. Hay que observar que la aceleración no tiene componente a lo largo de la binormal.

Cuando se trabaja en el plano x, y el vector binormal es el eje z (vector k) y el plano osculador es el propio plano x, y.

Nota. Este Párrafo está tomado casi literalmente del libro de Mecánica de Beer, Johnston and  Cornwell, 9ª edición,pag 667

6.Curvatura y radio de curvatura.

Sea una gráfica en 2D o en 3D.

Si R(t) =(f(t), g(t), h(t))              el vector posición de una partícula P, en términos del parámetro t, que podría ser el  tiempo.

v(t) = R’(t) = (f’(t), g’(t), h’(t))                           Velocidad tangencial de la partícula en P 

a(t) = R”(t) = ((f”(t), g”(t), h”(t)                          Aceleració tangencial de la partícula en P 

Históricamente, se ha definido la curvatura para curvas 2D y 3D así:

k(t) = d (et)/ds                                                                                             (7)

Es la rata de variación del vector unitario et (también nominado T(t)) respecto de la variación del arco. La curvatura escalar será:

ABS(kt) = k(t) > 0

Se define radio de curvatura como el inverso de k(t)        ρ = 1/k(t)          >0       (8)

En otro blog se ampliarán ambos conceptos, curvatura y radio de curvatura.

Las fórmulas que salen de estas definiciones son:

k(t) = ABS[R’(t)X R”(t)]/[ABS(R’(t)]3 = │R’(t)XR”(t)│/│R’(t)│3                                 (9)

y el radio de curvatura ρ = 1/k(t)                                                                                 (10)

Para una curva en x, y tenemos y = f(x)

La parametrisamos así: x = t              y = f(t)           R(t) = (t, f(t),0)

R’(t) = (1, f’(t), 0)          Velocidad

R”(t) = (0, f”(t), 0)          Aceleración

Encontramos el producto vectorial R’(t) X R”(t) = f”(t) k   y     R’(t) X R”(t) │= f”(t)

Encontramos además, │R’(t)│ = √(1 + [f”(t)]2)

k = f”(t) / [√(1 + (f”(t))2]3 = f”(t) / [(1 + (f”(t))2]3/2                                                     (11)

y recordando que x = t

k = f”(x) / [(1 + (f”(x))2]3/2             ρ = [(1 + (f”(x))2]3/2/ f”(x)                                    (12)

 

Ejercicio 3

Un avión va en trayectoria horizontal a velocidad constante. En cierto punto comienza a descender en una trayectoria circular y conserva la velocidad tangencial. Los pilotos saben que, si llegan a obtener una aceleración de 5,5g (g = 9,8 m/s2) se desmayan. El radio de la circunferencia es 325 m (Como lo indica la figura).

a.A qué velocidad se desmaya el piloto?

b. En otro caso similar, si el piloto va con una velocidad de 150 m/s y desciende por una trayectoria circular, manteniendo la velocidad de 150 m/s, cual será el radio, para que el piloto se desmaye.

Figura ejercicio 3

a. Sabemos que, aunque el avión va con una velocidad constante, existe una aceleración normal, que le produce el cambio de dirección al avión.  At = 0 y an = v2/r = 5,5g

v = √(5,5x9,8x325) = 132,42 m/s

b. Sabemos que, el avión lleva una velocidad tangencial constante igual a 150 m/s.

at = 0 y an = v2/r                                                 r = v2/ an = (159)2/(5,5x9,8) = 417,01 m

Ejercicio 4

El automóvil mostrado en la figura, entra a la curva, en el punto O, a una velocidad tangencial desconocida y el conductor aplica allí, una aceleración tangencial de at =.-0,25m/s2, .Cuál es magnitud de la aceleración, en el punto A, de x = 200 m, si en ese punto lleva una velocidad tangencial de 10 m/s

La gráfica del ejercicio es:



Figura ejercicio 4

El cambio de dirección en la velocidad, implica la existencia de una aceleración normal, que será igual a                        an = v2/r =                               r = ρ: radio de curvatura en el punto indicado.

at = -0,25 m/s2    aceleración tangencial

El valor del radio de curvatura ρ = r = ((1+ (y’(x))2) 3/2) /y” (x)             (1)

Para lo cual necesitamos y’ y y”           de la función y = 100e x/500     

y’ = (1/5) ex/500   

y en x = 200    y’ = (1/5) e 200/500

y” = (1/2500) e x/500                        y en x = 200           y" = (1/2500) e 200/500       

Llevando estos valores a la fórmula (1) (de este ejercicio) r (200) = 1904,48 m

Ahora, la componente tangencial de la aceleración es igual at = -0,25 m/s2

La componente normal será:     an = v2/r= (100m2/s2)/1904,48 m = 0,0525 m/s2

│a│=√(-0,25)2 + (0,0525)2) = 0,255 m/s2

Ejercicio 5

Tal como lo indica la figura del ejercicio, si el pasador sube, también sube la bolita que está en el canal. El pasador está subiendo desde el origen, a una velocidad constante de 3 m/s. Calcular la magnitud de la aceleración de la bolita, cuando la bolita está en x= 6

 


Figura problema 5  

vy = 3 m/s       Uniforme              x e y son funciones del tiempo

x = y2/6                    vx = dx/dt = (2y dy/dt) /6 = (1/3) ydy/dt = (1/3) y vx = y        ya que vy = 3 

Esto indica que las magnitudes de vy de y son iguales, cada una en las unidades consecuentes.

vx = y                                           y= √(6x)    y cuando x = 6 y = 6

La velocidad vx = 6 m/s

La magnitud de la velocidad total será v = √ (32 + 62) = √45 = 3√5 m/s

La aceleración según x será:    ax = dvx/dt = dy/dt = vy = 3 m/s2     la aceleración ay = 0

│a│ = √ (02 + 32) = 3 m/s2

 

 

 

Juan Fernando Sanín E

juanfernando.sanin@gmail.com