Forma de la cuerda giratoria (¿problema de lazo?)

Tomemos un alambre o una cuerda. Normalmente hago esto con una cadena o mi bufanda.

Fijo un extremo en mi mano y aplico la rotación (mediante movimientos sutiles de este punto final como girar un lazo). La cuerda entra en rotación y adquiere cierta forma doblada:ingrese la descripción de la imagen aquí

Falta una parte porque las cámaras de los teléfonos celulares no son buenas para fotos de alta velocidad, pero espero que puedas imaginar toda la bufanda.

La pregunta es: ¿cómo puedo calcular/predecir esta forma?

Aunque este problema no parece tan extraño, nunca he visto ninguna solución. Tampoco he encontrado esta pregunta en ningún lado de internet... Debe ser porque simplemente no sé cómo formularla sin imágenes.

Además, al tomar una cuerda más larga , obtengo más de una curva:ingrese la descripción de la imagen aquí

Me disculpo por la calidad de nuevo. Es aún más difícil rotar esto mientras se toma una foto. La forma no es espiral, es más como una forma en un plano que gira.

Agradeceré explicaciones, soluciones, enlaces o al menos una correcta formulación de este problema.

¿ Esto ayuda?
No lo leí todo, pero parece que no es una solución para encontrar la curva. En su lugar, eligen que la curva sea lineal por partes (2 piezas) y solo buscan algunos ángulos. También tienen honda como articulación, pero en este caso, aunque optemos por el modelo de 2 piezas lineales, las longitudes de las piezas seguirían siendo una incógnita...
Esta forma se describe en la 'Teoría del globo giratorio' de C. Mack ( QJ Mech. Appl. Math. 11 no. 2, 196–207 (1958) ) y también en 'Rotating Strings' de JA Hanna ( J. Phys. A : Math. Theor. 46 n.º 23, 235201 (2013) , arXiv:1204.4941 ).
Este es un problema de dinámica extremadamente complicado que no se puede resolver en los minutos o incluso las horas que la gente tarda en publicar una respuesta aquí. Tenga cuidado al confiar en cualquier solución que alguien ofrezca, independientemente de cuán sofisticadas se vean las ecuaciones. En otras palabras, no se deje desconcertar por BS.

Respuestas (7)

Es un problema interesante. He tratado de estudiarlo despreciando la fricción del aire, que probablemente no sea despreciable. Así que vamos a parametrizar la longitud de la cuerda L por el vector ( z ( yo ) , r ( yo ) ) , con yo siendo la coordenada curvilínea a lo largo de la cuerda y ( z , r ) las coordenadas cilíndricas de la cuerda (he omitido θ , que se puede suponer que es constante en el marco giratorio [ver mi comentario debajo de esta publicación para una justificación de esta omisión]). En el marco giratorio, la energía potencial de la cuerda es

mi pags = 0 L m d yo ( gramo z 1 2 Ω 2 r 2 )
y la definición de la integral de curvilínea da la restricción r 2 + yo 2 = 1 . Entonces el problema se puede traducir en minimizar la siguiente cantidad:
0 L d yo L ( z , r ; z , r ; yo ) con  L ( z , r ; z , r ; yo ) = m gramo z 1 2 m Ω 2 r 2 λ ( z 2 + r 2 ) .
No se debe olvidar (como hice en una primera respuesta tentativa) que el multiplicador de Lagrange λ depende de yo / Es un problema estándar de Euler-Lagrange y su solución viene dada por
L z d d yo L z = 0 y L r d d yo L r = 0
La ecuación en z se convierte
m gramo = 2 d d yo ( λ z ) λ z = 1 2 m gramo ( yo yo 0 ) ,
dónde yo 0 es una constante de integración. Para la solución óptima z 0 , desde que reemplazó z b y | z | solo puede disminuir L , y tenemos z = 1 r 2 . Uno tiene por lo tanto
λ = m gramo ( yo yo 0 ) 2 1 r 2
. La ecuación en r es entonces
m Ω 2 r = 2 d d yo ( λ r ) = d d yo m gramo ( yo yo 0 ) r 1 r 2
Tenemos entonces la siguiente ecuación para integrar
0 = Ω 2 gramo r ( 1 r 2 ) 3 / 2 r ( 1 r 2 ) ( yo yo 0 ) r .
No sé cómo integrarlo analíticamente, pero cuando yo yo 0 < 0 , parece que podemos tener oscilaciones.

Editado para agregar:

Si uno mira las señales de r y r como función de r y r por un fijo yo , uno puede dibujar rápidamente algunas flechas en el espacio de fase ( es decir ) y ver que el flujo va "en círculos" cuando yo < yo 0 y tiende asintóticamente a la curva Ω 2 gramo r = r 1 r 2 cuando yo > yo 0 . Entonces uno puede esperar un número finito de oscilaciones de la cuerda. Pero toda la cuestión es entonces cómo yo 0 se relaciona con los parámetros del sistema.

Editar: en el límite de la cuerda casi vertical ( r 1 ), la ecuación se convierte en

0 = Ω 2 gramo r r ( yo yo 0 ) r
que es mas simple Y, de hecho, Wolfram Alpha tiene una solución analítica utilizando funciones de Bessel modificadas.

La condición de Lagrange debe ser z 2 + r 2 + ( r ϕ ) 2 = 1 , aplicado en todos los puntos de la curva, por lo que debería haber una función multiplicadora de Lagrange λ ( yo ) . Poner todo esto en Mathematica da ecuaciones diferenciales no lineales horriblemente acopladas que parecen ser difíciles de resolver. Además, no estoy 100% seguro de cuáles son las condiciones de contorno correctas que debo aplicar (o incluso cuántas de ellas necesito...)
Buen intento hasta ahora, ni siquiera sabía qué herramientas elegir aquí... Si lo hice bien, el r = 1 Ω F ( z ) ya muestra que esto va mal como este r ( z ) tiene solo un cero, mientras que la solución real muestra al menos dos.
@genneth Gracias por la pista sobre λ ( yo ) . He corregido la respuesta en consecuencia y parece más manejable.
Sin embargo, todavía te falta la parte tangencial del problema. ¡Su configuración no permite una cuerda que se enrolle alrededor del eje (como lo hace la imagen del OP)!
Además, debes tener en cuenta que el campo multiplicador de Lagrange es en realidad la tensión a lo largo de la cuerda. Creo que la solución completa es en realidad en términos de una ecuación diferencial-integral, ya que las condiciones de contorno son tales que no conoce la tensión requerida en la parte superior de la cuerda sin conocer la configuración --- solo sabe que la tensión debe ser cero en la punta de la cuerda y necesita integrar una copia de seguridad.
@genneth: De acuerdo con el multiplicador de LAgrange. Sobre la falta de un grado de libertad angular, lo he omitido por 2 razones: 1. Al contrario de lo que dices, la imagen de Juris no muestra la cuerda enrollada alrededor del eje. Él dice explícitamente: "La forma no es espiral, es más como una forma en un plano que gira". 2. Sin fricción, es fácil demostrar que la forma óptima es plana: si un pequeño segmento de cuerda tiene las tres coordenadas ( d r , r d ϕ , d z ) , reemplazándolo por un segmento de coordenadas ( d r , 0 , d z 2 + r 2 d ϕ ) sólo puede disminuir la energía potencial.

Creo que la primera parte de la curva (la parte entre su mano y el "nodo" estacionario donde la curva se cruza sobre el eje vertical) estaría descrita por la Curva Troposkein , o al menos sería extremadamente similar.

Me imagino que esta curva describiría completamente, o al menos estaría relacionada con, cualquier otro segmento entre nodos si tuviera más de uno.

El último trozo de "cola" en la parte inferior de la cuerda sería diferente, obviamente, pero probablemente podría obtener la forma usando un método similar.

Resolví ecuaciones escritas por @genneth numéricamente (para la condición inicial usé un método de disparo). La parte crucial es que no usé la ecuación.

( T r 2 ϕ ) = 0
Es decir, supuse que la curva es plana. Mis resultados están abajo.$\omega=10$$\omega=20$$\omega=40$$l=100$

TENGA EN CUENTA que en las 3 primeras imágenes solo ω cambios. En la última foto cambié de largo yo . Estas imágenes y muchas otras me permiten concluir que la cantidad de ``modos'' norte a es decir, puntos donde r = 0 es

norte a yo en ω
para grande ω .

Esta no es una solución completa, pero la estoy convirtiendo en una wiki comunitaria para que otros puedan complementarla y corregirla.

Frédéric Grosshans tiene casi razón en lo anterior, pero pasa por alto una parte crucial del problema: la cuerda puede enrollarse alrededor del eje. El método general de suponer una forma de equilibrio estable y minimizar la energía potencial (efectiva) en un marco giratorio es bueno, pero tal vez uno debería tener cuidado de suponer que tal equilibrio estable incluso existe.

Presionando, parametrizamos la cuerda como ( z ( yo ) , r ( yo ) , ϕ ( yo ) ) , y la energía potencial es:

tu = 0 L d yo ρ ( 1 2 ω 2 r 2 gramo z )
dónde ρ es la densidad por longitud (aparentemente irrelevante al final, pero la mantendremos para asegurarnos de que las unidades sean sensatas) y tomamos la convención de que z aumenta hacia abajo.

La restricción es que cada elemento pequeño de la cuerda debe tener una longitud fija, por lo que introducimos un campo multiplicador de Lagrange (¡esta es una teoría de campos en 1D!):

L = ρ ( 1 2 ω 2 r 2 gramo z ) + 1 2 T ( r 2 + z 2 + r 2 ϕ 2 1 )
dónde T es el campo multiplicador de Lagrange (el nombre será obvio, al igual que el factor aleatorio de 1/2).

Empujar a través de las ecuaciones de Euler-Lagrange produce:

ρ gramo + ( T z ) = 0 ρ ω 2 r + T r ϕ 2 ( T r ) = 0 ( T r 2 ϕ ) = 0 r 2 + z 2 + r 2 ϕ 2 = 1

La primera ecuación codifica el significado de T como la tensión en la cuerda. Podemos integrarlo y obtener el obviamente sensato

T z = ρ gramo ( L yo )
(usando la condición de frontera que T = 0 a yo = L ) que simplemente codifica que el componente vertical de la tensión debe equilibrar el peso de la cuerda debajo.

Las ecuaciones restantes son de primer orden en z y segundo orden en ϕ y r . Tenemos las ecuaciones de establecimiento de coordenadas que z ( 0 ) = 0 y ϕ ( 0 ) = 0 . Hay un parámetro extra genuino r ( 0 ) que dicta el radio de la fuerza motriz en la parte superior. Sin embargo, todavía se necesita otra condición de contorno para r y ϕ (o posiblemente sus derivados) para que el problema quede bien planteado. Además, la condición para T conduce a una ecuación integral engañosa, ya que la condición de contorno natural es que T ( L ) = 0 , y entonces sería necesario integrar esto para obtener el límite en yo = 0 .

Finalmente, debo mencionar que esta solución parcial, tal como está, tiene una deficiencia grave que puede o no estar relacionada con el problema de la condición de contorno anterior: las ecuaciones son simétricas para ω ω , por lo que aparentemente una solución disponible tendría la cuerda apuntando hacia adelante de la fuerza de rotación. Pero tal vez este sea un máximo de la energía en lugar del mínimo.

En cualquier caso, el problema parece altamente no lineal, y me sorprendería (y me encantaría) si fuera posible una solución analítica. Uno podría esperar, sin embargo, que haya regímenes fáciles de entender (aparte de los triviales en ω = 0 o ).

la simetria ± ω está ligado al hecho de que la fuerza centrífuga no depende del sentido de giro. Creo que la única forma de romper esta simetría es agregar fricción.
Puede ser como dice en un comentario a su respuesta anterior: la cuerda puede estar en un avión. Experimentalmente, esto parece cierto (solo lo probé :)), pero aún no lo he visto solo desde las matemáticas ... ¿Alguna idea? estoy pensando eso T r 2 ϕ debe tener alguna interpretación física que ayude.
Echemos un vistazo a la ecuación ( T r 2 ϕ ) = 0 cuando se integra: T r 2 ϕ = C = 0. Parece identidad.
Creo que su expresión para energía potencial es incorrecta, más precisamente el término centrífugo. Dice que la fuerza centrífuga tiene dirección contra r lo cual no es cierto
Hay una solución analítica en forma cerrada, vea mi solución a continuación.

Tengo una suposición probable de que este fenómeno está relacionado de alguna manera con (o es similar a) las ondas estacionarias.

Si sostiene un extremo de la cuerda y realiza un movimiento armónico simple horizontal (es decir, hace oscilar el extremo horizontalmente), se creará una onda transversal plana y cuando se desplace hacia el otro extremo se reflejará (como las ondas se reflejan en un punto libre). extremo del medio), por lo que se forma una onda estacionaria, siendo el extremo libre un antinodo. Dependiendo de la longitud de la cuerda y de la velocidad angular de la oscilación, se formarán uno o varios nodos a lo largo de la cuerda.

Considere el movimiento de la cuerda giratoria como el resultado de dos ondas transversales que viajan en dos planos perpendiculares entre sí, y la explicación de la onda estacionaria tiene algún sentido.

No estoy seguro de eso, solo una suposición.

De hecho, está relacionado con una onda estacionaria, vea la solución a continuación. Las soluciones están parametrizadas por el número de nodos.

Mostraremos que la forma de la cuerda giratoria es la función de Bessel.

Adoptemos un sistema de coordenadas giratorio en el que la cuerda está en reposo. Usaremos coordenadas cilíndricas, r y z son las coordenadas radial y vertical, donde ( r , z ) = ( 0 , 0 ) denota el punto de suspensión y z aumenta hacia abajo. La cuerda corre a lo largo de los puntos. ( r ( yo ) , z ( yo ) ) , dónde yo es la longitud a lo largo de la cuerda medida desde el extremo libre, es decir yo = L en el punto de suspensión.

Las fuerzas que actúan sobre un elemento de cuerda infinitesimal de masa d metro son las fuerzas gravitacional y centrífuga ( d metro gramo y d metro ω 2 r ) y las fuerzas de restricción que surgen debido a la tensión en la cuerda que aseguran que la densidad de la línea se mantenga η = metro / L en todas partes a lo largo de la cuerda, donde metro y L son respectivamente la masa total y la longitud. La tensión en la cuerda en yo ejerce una fuerza F T ( yo ) que es paralelo al vector tangente de la cuerda. Para un elemento de cuerda entre yo d yo y yo + d yo , surgen fuerzas en los puntos límite, es decir F T ( yo + d yo ) y en F T ( yo d yo ) .

Las ecuaciones de movimiento radial y vertical son

d metro d 2 r d t 2 = d metro ω 2 r F T ( yo + d yo ) pecado θ ( yo + d yo ) + F T ( yo d yo ) pecado θ ( yo d yo ) = 0 , d metro d 2 z d t 2 = d metro gramo F T ( yo + d yo ) porque θ ( yo + d yo ) + F T ( yo d yo ) porque θ ( yo d yo ) = 0 .
Aquí θ denota el ángulo subtendido por la tangente de la cuerda con la vertical, es decir broncearse θ = d r / d z . Esto generalmente cambia con la posición a lo largo de la cuerda que escribimos como θ ( yo ) . El lado derecho es cero dado que la cuerda está en reposo en el marco giratorio. Tenga en cuenta que la cuerda debe ser plana en reposo, ya que no actúan fuerzas en la dirección azimutal y, por lo tanto, la tensión de la cuerda no puede tener un componente azimutal, la cuerda debe tener una derivada cero a lo largo de la dirección azimutal.

Divide las ecuaciones con 2 d yo :

d d yo ( F T pecado θ ) = η ω 2 r , d d yo ( F T porque θ ) = η gramo .
La ecuación vertical se puede integrar como
F T porque θ = η gramo yo .
Esto satisface las condiciones de contorno ya que en el extremo inferior yo = 0 asi que F T = 0 , y en yo = L tenemos porque θ = 1 lo que implica que F T = η gramo L = metro gramo . Sustituyendo en la ecuación radial y usando eso pecado θ = broncearse θ porque θ = ( d r / d z ) porque θ , da
d d yo ( η gramo yo d r d z ) = η ω 2 r .

en el limite que θ 1 , podemos aproximarnos yo L z . Para simplificar las ecuaciones, redefina el sistema de coordenadas con la sustitución z L z , de modo que z = 0 es el extremo inferior libre de la cuerda y z = L es el extremo fijo superior. Esto produce

d d z ( η gramo z d r d z ) η ω 2 r .
Esta aproximación nos permite obtener una solución analítica exacta en forma cerrada. En la práctica encontramos numéricamente que la solución exacta está muy bien descrita por esta aproximación.

Ahora manipulamos esta ecuación para llegar a la ecuación de Bessel. Introducir una nueva variable X = 2 ω z / gramo de modo que z = X 2 gramo / ( 4 ω 2 ) . La regla de la cadena implica que la derivada con respecto a z puede escribirse como

d d z = d X d z d d X = 1 d z / d X d d X = 2 ω 2 gramo X d d X .
Las constantes desaparecen después de la sustitución y obtenemos
1 X d d X ( X d r d X ) = r
o
X 2 d 2 r d X 2 + X d r d X + X 2 r = 0.
Esta es la ecuación diferencial de Bessel que define las funciones de Bessel j 0 ( X ) y Y 0 ( X ) , la solución general es una combinación lineal de estas funciones. Ya que Y 0 ( X ) diverge en X = 0 y j 0 ( 0 ) = 1 , la solución física es proporcional a la función de Bessel de primera clase j 0 ( X ) . Por lo tanto, la forma de la cuerda giratoria es

r ( z ) = r mi norte d j 0 ( X ( z ) ) = r mi norte d j 0 ( 2 ω z gramo ) .

Aquí r mi norte d es la coordenada radial en el extremo suelto. La condición de contorno en la suspensión es que en z = L la coordenada radial satisface r = 0 , por lo que debemos exigir que j 0 ( 2 ω L / gramo ) = 0 . Los ceros de las funciones de Bessel son bien conocidos, están tabulados. Denotemos la norte el cero con X norte = { 2.4048 , 5.5201 , 8.6537 , 11.7915 , } . La condición de frontera se cumple cuando ω = X norte 2 gramo / L para cualquier norte . De este modo

r ( z ) = r mi norte d j 0 ( X ( z ) ) = r mi norte d j 0 ( X norte z L ) .

Para frecuencias más grandes, la cuerda tiene múltiples nodos, es decir, lugares donde se cruza r = 0 . la cuerda tiene k nodos si la frecuencia angular es ω = X k + 1 2 gramo / L , un nodo para ω = X 2 2 gramo / L = 2.26 gramo / L , dos nodos si 4.33 gramo / L , y tres nodos si 5.896 gramo / L .

en el límite, ω gramo / L podemos usar la fórmula asintótica para los ceros de las funciones de Bessel X k + 1 ( k 3 4 ) π , por lo que tiene k número de nodos cuando

ω = π 2 ( k + 3 4 ) gramo L .
Resolviendo para k muestra que el número de nodos para ω gramo / L es
k = 2 ω π L gramo 3 4 .
En este caso, la forma asintótica de la función de Bessel
j 0 ( X ) = 2 π X porque ( X π 4 ) + O ( X 3 / 2 )
implica que para grandes ω y largo z
r ( z ) r mi norte d 1 π ω gramo z porque ( 2 ω z gramo π 4 ) .

Tenga en cuenta que asumimos que el extremo superior fijo en z = L satisface r ( L ) = 0 . En el caso más general, las condiciones de contorno son r ( L ) = r 0 y r ( 0 ) = r mi norte d . Esto implica que j 0 ( 2 ω L / gramo ) = r 0 / r mi norte d , y obtenemos una solución similar para grandes ω si r 0 / r mi norte d 1 .

¡Ese es un resultado fuerte!

Será una especie de función hiperbólica, y la fuerza y ​​la longitud de la cuerda determinarán la frecuencia para predecir la forma. El 'círculo' inicial en su muñeca definirá la fuerza constante para explicar la curva de pecado. Esta no es una muy buena explicación, pero si tuviera que definir una función hiperbólica usando el eje z, podría ser sinh z=(e^ze^-z)/2.

La segunda curva no me parece extremadamente hiperbólica, ¿sería una función hecha de piezas sinh y/o cosh?