Ver el código de «Dinámica impulsiva analítica (CMR)»

==Introducción==
Al estudiar la mecánica en su formulación vectorial introducimos las [[Din%C3%A1mica_impulsiva_vectorial_(CMR)|fuerzas impulsivas]] como aquellas que era de una gran intensidad pero actuaban durante un breve periodo de tiempo. Idealmente, una fuerza impulsiva sería una de intensidad infinita que actúa durante un intervalo infinitesimal. Cuando actúa una fuerza de este tipo, se dice que el sistema ha experimentado una ''percusión''.

Para caracterizar una fuerza impulsiva que actúa entre <math>t_1</math> y <math>t_1+\Delta t</math>, definimos el ''impulso''

<center><math>\vec{P}=\int_{t_1}^{t_1+\Delta t}\vec{F}\,\mathrm{d}t</math></center>

Esta cantidad mide la magnitud de la fuerza impulsiva, pero no es una fuerza, ya que se mide en N·s.

El efecto de una percusión sobre una partícula es modificar su velocidad, pero no su posición

<center><math>m\,\Delta \vec{v}=m\left(\vec{v}_+-\vec{v}_-\right)=\vec{P}\qquad\qquad \Delta\vec{r}=\vec{r}_+-\vec{r}_-=\vec{0}</math></center>

Aquí el incremento se calcula como la diferencia entre el estado justo tras la percusión menos el estado inmediatamente previo a ella. 

La posición no cambia porque se trata de la integral de la velocidad en el intervalo que dura la percusión. Dado que la velocidad no se hace infinita (aunque sea casi discontínua) y el intervalo tiende a 0, resulta que esta integral se anula.

Para el caso de un sistema de partículas con percusiones en los puntos <math>A_i</math>, tenemos los análogos al teorema de la cantidad de movimiento (TCM)

<center><math>\Delta \vec{p}=m\,\Delta \vec{v}_G=\sum_i \vec{P}_i^\mathrm{ext}</math></center>

y del teorema del momento cinético

<center><math>\Delta \vec{L}_O=\sum_i \overrightarrow{OA}_i\times\vec{P}^{ext}_i</math></center>

Aquí O es cualquier punto, no necesariamente inmóvil o el CM, aunque suele aplicarse a estos dos últimos casos.

Las percusiones internas se anulan como consecuencia de la tercera ley de Newton.

Las percusiones, tanto internas como externas, pueden ser aplicadas o pueden ser de reacción vincular.
==Impulsos generalizados==
Al escribir las ecuaciones de la mecánica en términos de coordenadas generalizadas, aparecen las fuerzas generalizadas

<center><math>Q_k = \sum_i F_i \frac{\partial x_i}{\partial q_k}</math></center>

La generalización a fuerzas impulsivas es inmediata. Si consideramos que <math>F_i</math> es una fuerza impulsiva e integramos sobre la duración de la percusión resulta

<center><math>P_k = \int_{t_1}^{t_1+\Delta t}Q_k\,\mathrm{d}t=\sum_i \left(\int_{t_1}^{t_1+\Delta t}\vec{F}\,\mathrm{d}t\right) \frac{\partial x_i}{\partial q_k}</math></center>

A esta cantidad se la conoce como ''impulso generalizado''. Puede corresponder a un impulso en una dirección o ser un par impulsivo que produce un giro.

==Ecuaciones de Lagrange impulsivas==
Una vez definidos los impulsos generalizados, las ecuaciones de la dinámica analítica en presencia de fuerzas impulsivas pueden escribirse de diferentes formas. Aquí solo consideraremos el caso más simple.

Las ecuacionesde Lagrange para un sistema pueden escribirse en la forma

<center><math>\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial K}{\partial\dot{q}_k}\right)-\frac{\partial K}{\partial q_k}=Q_k</math></center>

Aquí <math>Q_k</math> incluye, en principio, todas las fuerzas: conservativas, no conservativas y de reacción vincular. En cada categoría, algunas pueden ser no impulsivas (actúan de forma distribuida en el tiempo, como el peso) o ser impulsiva (por ejemplo, colisiones). Separando en estas dos categorías

<center><math>\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial K}{\partial\dot{q}_k}\right)-\frac{\partial K}{\partial q_k}=Q^\mathrm{imp}_k+Q^\mathrm{no\,imp}_k</math></center>

Si ahora integramos esta ecuación respecto al tiempo sobre el breve  intervalo que dura la percusión, obtenemos los siguientes cuatro términos:

Para las fuerzas impulsivas

<center><math>\int_{t_1}^{t_1+\Delta t}Q^\mathrm{imp}_k\,\mathrm{d}t=P_k</math></center>

Para las no impulsivas

<center><math>\int_{t_1}^{t_1+\Delta t}Q^\mathrm{no\,imp}_k\,\mathrm{d}t=0</math></center>

La razón de que se anule esta integral es que estamos considerando un integrando finito (ya que las fuerzas no impulsivas, como el peso, no son infinitamente intensas, sino mucho más débiles que las impulsivas) sobre un intervalo infinitesimal. 

De los términos del primer miembro, tenemos la para la derivada temporal de los momentos generalizados

<center><math>\int_{t_1}^{t_1+\Delta t}\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial K}{\partial\dot{q}_k}\right)\,\mathrm{d}t=\Delta\left(\frac{\partial K}{\partial\dot{q}_k}\right)</math></center>

ya que la integral de una derivada da el incremento de la función original.

Por último tenemos que

<center><math>\int_{t_1}^{t_1+\Delta t}\frac{\partial K}{\partial q_k}\,\mathrm{d}t=0</math></center>

Este término se anula por la misma razón que el de las fuerzas no impulsivas. La energía cinética es una función de las coordenadas <math>q_k</math> y velocidades <math>\dot{q}_k</math> (no de las aceleraciones). Su derivada respecto a una coordenada es otra función de las <math>q_k</math> y <math>\dot{q}_k</math>. Ninguna de estas cantidades es infinita, aunque pueda ser discontinua. Por tanto su integral sobre un intervalo infinitesimal se anula.

Reuniendo todo esto queda finalmente

<center><math>\Delta\left(\frac{\partial K}{\partial\dot{q}_k}\right)=P_k</math></center>

es decir, lo que varía el momento generalizado es igual al impulso generalizado.

==Fuerzas de reacción impulsivas==
La presencia de vínculos provoca la aparición de fuerzas impulsivas de reacción. Éstas pueden deberse a varias causas, según el tipo de vínculo. Pueden clasificarse cuatro casos según los vínculos existan antes, durante y después de la percusión:

;Vínculos preexistentes persistentes: son aquellos que existen antes, durante y después de la percusión. 
:Llevan aparejados fuerzas impulsivas en respuesta a otras percusiones. Así, por ejemplo, una viga empotrada que es golpeada en algún punto provoca la aparición d euna fuerza impulsiva y de un par impulsivo en el mpunto de empotramiento.
;Vínculos impuestos persistentes: son los que existen durante y después de la percusión, pero no antes de ella. Es lo que ocurre, por ejemplo al tensarse el hilo de un péndulo flexible.
:Esta situación implica fuerzas impulsivas porque las velocidades previas a la imposición pueden no cumplir el vínculo y por tanto deben cambiar bruscamente. Esto se manifiesta en que el hilo pega un tirón al tensarse.
;Vínculos anulados: son los que existen antes y durante la percusión, pero no después de ella. Siguiendo con el ejemplo del péndulo, sería el caso de que se rompa la cuerda.
:A diferencia de la imposición, la anulación no tiene por qué provocar fuerzas impulsivas. La razón es que la eliminación del vinculo significa una mayor generalidad en las velocidades posibles y por tanto las velocidades previas ya cumplen las restricciones posteriores, que son menores.
;Vinculos impulsivos: son los que solo existen durante el tiempo que dura la percusión. Así, por ejemplo, en una colisión entre dos sólidos, durante el tiempo que dura ésta, se cumple la condición de impenetrabilidad y por tanto la velocidad normal relativa debe ser nula.
:Dado que en los análisis de percusiones lo que se estudia es la relación entre el estado previo y el posterior, y no lo que ocurre durante, los vínculos impulsivos requieren algún tipo de modelo de cómo es esa colisión, ya que no es lo mismo que los sólidos se comporten como cuerpos rígidos a que sean deformables. Es por ello, que se habla de colisiones elásticas o inelásticas, como forma empírica de caracterizar la colisión sin entrar en el detalle de cómo se produce.

<center>[[Archivo:clasificacion-percusiones.png|776px]]</center>
Cuando tenemos un vínculo escrito en la forma

<center><math>\sum_k B_k \dot{q}_k + B_0=0\,</math></center>

admitimos que los coeficientes <math>B_k</math> no dependen de las velocidades, sino solo de las coordenadas. Sabemos que en este caso, las componentes de las fuerzas de reacción son de la forma

<center><math>Q_k = \lambda B_k\,</math></center>

siendo <math>\lambda</math> el multiplicador de Lagrange correspondiente.

Si integramos esta ecuación respecto al tiempo, obtenemos los impulsos generalizados

<center><math>P_k = \Lambda B_k\,</math></center>

donde los coeficientes son los mismos ya que si solo dependen de las coordenadas y estas no cambian durante la percusión, los coeficientes son constantes para esta integral. Por tanto solo habría que integrar el multiplicador de Lagrange lo que nos da una cierta cantidad, también desconocida, que habrá que hallar en la solución del problema.

Esto es una aproximación, ya que supone que la ecuación del vínculo no se ve afectada por la percusión, lo cual es discutible. En la colisión de dos esferas de acero, por ejemplo, las fuerzas pueden ser tan importantes como para producir una deformación apreciable, incumpliéndose la condición de rigidez.

Admitiendo que las fuerzas impulsivas de reacción cumplen esta condición, la ecuación de Lagrange para la dinámica impulsiva con diversos vínculos queda en la forma

<center><math>\Delta\left(\frac{\partial K}{\partial \dot{q}_k}\right)=P_k+\sum_j \Lambda_j B_{jk}</math></center>

donde aquí las <math>P_k</math> ya solo incluyen las fuerzas impulsivas aplicadas.

La estrategia más efectiva para abordar la imposición de vínculos suele ser el elegir las coordenadas de manera que el nuevo vínculo se convierta en trivial de manera que el impulso generalizado aparezca en una sola de las ecuaciones de Lagrange.

{{ejemplo|'''El péndulo que se tensa''' 

[[Archivo:Pendulo-que-tensa.png|300px|right]] 
Supongamos una masa que cuelga de un hilo flexible de longitud <math>\ell</math>. Inicialmente la masa se halla a la misma altura que el punto al que se encuentra atado el hilo, y a una distancia d de este punto. Se deja caer la masa. ¿Cuánto vale la tensión impulsiva en el hilo?

Si empleamos coordenadas polares, la ecuación del vínculo es 

<center><math>\rho\leq \ell</math></center>

Antes de que se tense se cumple la desigualdad y es como si este vínculo no existiera. Una vez tensado se cumple la igualdad.

La energía cinética de la partícula es

<center><math>K=\frac{1}{2}m(\dot{\rho}^2+\rho^2\dot{\theta}^2)</math></center>

Cuando el hilo se tensa se impone el vínculo

<center><math>\rho=\ell\qquad\Rightarrow\qquad \dot{\rho}=0\qquad\Rightarrow\qquad B_\rho=1,\ B_\theta=0</math></center>

lo que nos da las ecuaciones para el impulso

<center><math>\Delta\left(\frac{\partial K}{\partial \dot{\rho}}\right)=\Lambda\qquad\qquad \Delta\left(\frac{\partial K}{\partial \dot{\theta}}\right)=0</math></center>

Obsérvese que el peso no aparece en estas ecuaciones, por ser una fuerza no impulsiva, su efecto es despreciable en este breve intervalo de tiempo.

Calculando los momentos conjugados nos queda

<center><math>m\left(\dot{\rho}_+-\dot{\rho}_-\right)=\Lambda\qquad\qquad m\ell^2\left(\dot{\theta}_+-\dot{\theta}_-\right)=0</math></center>

La segunda de las ecuaciones nos permite calcular la velocidad justo después del tirón, que es la misma que antes de él. La velocidad anterior se halla resolviendo las ecuaciones de movimiento previas. En este caso es simple por conservación de la energía y proyectando posteriormente en las direcciones radial y acimutal

<center><math>\dot{\rho}_-=\sqrt{2gh}\cos(\theta_0)\qquad\qquad \dot{\theta}_-=-\frac{\sqrt{2gh}}{\ell}\mathrm{sen}(\theta_0)</math></center>

Aquí h es lo que desciende y <math>\theta_0</math> el ángulo que forma con la vertical en el momento en que se tensa

<center><math>h=\sqrt{\ell^2-d^2}\qquad\qquad \mathrm{sen}(\theta_0)=\frac{d}{\ell}</math></center>

Por tanto, tras el tirón

<center><math>\dot{\theta}_+=\dot{\theta}_-=-\frac{d\sqrt{2gh}}{\ell^2}</math></center>

La magnitud de la percusión la obtenemos del salto en la velocidad radial

<center><math>\Lambda = m(\overbrace{\dot{\rho}_+}^{=0}-\dot{\rho}_-)=-\frac{mh\sqrt{2gh}}{\ell}</math></center>

Es negativo porque va en la dirección radial hacia adentro. Si el tirón se produce en un tiempo <math>\Delta t</math> la tensión impulsiva tendrá el valor

<center><math>F_t\simeq \frac{\Lambda}{\Delta t}= -\frac{mh\sqrt{2gh}}{\ell\,\Delta t}</math></center>

Esta cantidad puede ser muy elevada, por lo que podría llegar a romperse el hilo (o escaparse de las manos). 


}}

[[Categoría:Mecánica analítica (CMR)]]