On 13 Nov., 10:30, Antonio González <
gonfe...@gmail.com> wrote:
> El 12/11/2012 20:03, Dr. Wolfgang Hintze escribió:
>
>
>
>
>
> > On 12 Nov., 18:59, Antonio González <
gonfe...@gmail.com> wrote:
> >> El 12/11/2012 17:21, Dr. Wolfgang Hintze escribi :
>
> >>> (no s porque mi mensaje aparec as , otra vez!)
> >>> El resultado aproximado tiene una interpretaci n muy simple: puesto
> >>> que la frecuencia de las deviaciones es la misma que la b sica, la
> >>> curva es muy bien aproximada (incluso hasta eps = 1) por un c nico,
> >>> es
> >>> decir una elipse.
>
> >> Necesariamente doblada, ya que una elipse no tiene curvatura constante.
>
> >> --
>
> >> Antonio
>
> > Por supuesto una elipse doblada. La curva en el espacio ya he dado:
>
> > vr[phi] = {Cos[phi]*(R + eps*Cos(phi)), Sin(phi)*(R + eps*Cos(phi)),
> > a*(R + eps*Cos(phi)}
>
> > Pero lo interesante es que la curva aproximada es cerrada.
> > ¿Se puede probarlo en general?, es decir ¿es la curva de curvatura
> > constante en un cono recto una curva cerrada?
>
> No tiene por qué. Ya te digo yo que me imagino una solución de la
> siguiente forma: comienza cerca del vértice siendo casi paralela a una
> generatriz (para compensar la gran curvatura de las secciones con un
> "enderezamiento" a lo largo). A medida que va pasando por secciones de
> radio creciente, se va haciendo cada vez más horizontal, tendiendo a una
> circunferencia transversal de curvatura la dada, más allá de la cual no
> puede pasar.
>
> --
>
> Antonio
>
>
Sí, lo que dices ya he dicho el 3 de noviembre "Por lo tanto la curva
buscada es una espiral que aproxima gradualmente el círculo r = rho0."
pero, sorprendentemente, parece que no sea así.
Además de la solución lineal (en la vecinidad de r=R) también he
integrado la ecuación exacta numéricamente con el resultado de un
comportamiento oscilatorico sin atenuación.
También, para ver si un comportamineto de atenuación sea válido, he
puesto r(phi) = R(1-Exp(-A phi)). Pero la curvatura no sale muy
constante.
Tal vez sea útil presentar las ecuaciones para las cuantidades
relevantes de la curva en el cono recto:
Tomamos el ángulo phi como parametro u de la curva:
In[512]:= \[Phi][u_] := u;
El vector de posición vr
In[513]:= vr[u]
Out[513]= {Cos[u]*r[u], r[u]*Sin[u], a*r[u]}
El vector tangencial vt = d/ds vr
In[514]:= Simplify[vt[u]]
Out[514]= {((-r[u])*Sin[u] + Cos[u]*Derivative[1][r][u])/
Sqrt[r[u]^2 + (1 + a^2)*Derivative[1][r][u]^2], (Cos[u]*r[u] +
Sin[u]*Derivative[1][r][u])/
Sqrt[r[u]^2 + (1 + a^2)*Derivative[1][r][u]^2],
(a*Derivative[1][r][u])/
Sqrt[r[u]^2 + (1 + a^2)*Derivative[1][r][u]^2]}
El vector normal (d/ds)^2 vr
In[516]:= FullSimplify[vn[u]]
Out[516]= {((-Cos[u])*
r[u]*(r[u]^2 + (2 + a^2)*Derivative[1][r][u]^2 -
r[u]*Derivative[2][r][u]) +
Sin[u]*Derivative[1][r][
u]*(-r[u]^2 -
2*(1 + a^2)*Derivative[1][r][u]^2 + (1 + a^2)*r[u]*
Derivative[2][r][u]))/
(r[u]^2 + (1 + a^2)*Derivative[1][r][u]^2)^2, ((-r[u])*
Sin[u]*(r[u]^2 + (2 + a^2)*Derivative[1][r][u]^2 -
r[u]*Derivative[2][r][u]) +
Cos[u]*
Derivative[1][r][
u]*(r[u]^2 +
2*(1 + a^2)*Derivative[1][r][u]^2 - (1 + a^2)*r[u]*
Derivative[2][r][u]))/(r[u]^2 + (1 + a^2)*
Derivative[1][r][u]^2)^2,
(a*r[u]*(-Derivative[1][r][u]^2 +
r[u]*Derivative[2][r][u]))/(r[u]^2 + (1 + a^2)*
Derivative[1][r][u]^2)^2}
La curvatura k[u] (= el inverso del radio de curvatura)
k[u_] = Simplify[ Sqrt[vn[u].vn[u]] ]
Sqrt[
(r[u]^4
+ 4*(1 + a^2)*r'[u]^4
- 2*r[u]^3*r''[u]
- 4*(1 + a^2)*r[u]*r'[u]^2*r''[u]
+ r[u]^2*((4 + a^2)*r'[u]^2
+ (1 + a^2)*r''[u]^2))
/(r[u]^2 + (1 + a^2)*r'[u]^2)^3
]
La ecucación diferencial para r[u] es k[u] = 1/R = cte.
Saludos,
Wolfgang