Ver el código de «Partícula en barra horizontal con rozamiento (MRGIC)»

= Enunciado =

==[[Partícula en barra horizontal con rozamiento (MRGIC)| Partícula en barra horizontal con rozamiento]]==
[[Archivo:Rozamiento.png|sinmarco|derecha]]
Una partícula <math>P</math> de masa <math>m</math> se ecuentra en un plano vertical y desliza libremente sobre una barra horizontal fija y rugosa (coeficiente dinámico <math>\mu</math>) situada en la posición <math>y=h(cte)</math>. Un resorte ideal de constante recuperadora <math>k</math> y longitud natural nula tiene uno de sus extremos en <math>P</math> y el otro en el origen de coordenadas <math>O</math>. Considera que la partícula se está moviendo hacia la derecha <math>(\dot{x}>0)</math>. Aplica el principio de liberación y, usando los métodos de la Mecánica Analítica, obtén la ecuación diferencial de movimiento de <math>P</math> y el valor de la fuerza normal.

= Solución =
[[Archivo:ParticulaRozamientoMRFuerzas.png|sinmarco|derecha]]
Una partícula libre tiene tres grados de libertad. En este problema, la partícula debe cumplir dos vínculos
<center>
<math>
y = h, \qquad z = 0.
</math>
</center>
Por tanto, el sistema tiene un grado de libertad. Podemos escoger la coordenada <math> \{ x \} </math>. 

Sin embargo, la presencia del rozamiento hace que no podamos tratar el problema con un sólo grado de libertad. Para poder incluir la fuerza de rozamiento, hay que usar el Principio de liberación, y tratar las fuerzas vinculares (incluida la de rozamiento), como fuerzas activas. Por cada vínculo que rompamos hay que añadir una componente de reacción vincular en la dirección del movimiento afectado por el vínculo. En este caso
<center>
<math>
\begin{array}{lcl}
y = h & \to & N_y\,\vec{\jmath}, \\
z = 0 & \to & N_z\,\vec{k}.
\end{array}
</math>
</center>
Agrupamos estas dos fuerzas en una sóla fuerza normal
<center>
<math>
\vec{N} = N_y\,\vec{\jmath} + N_z\,\vec{k}.
</math>
</center>
Como la partícula desliza sobre la barra, estamos en régimen de rozamiento dinámico. En este régimen el módulo de la fuerza de rozamiento es conocido. Su dirección y sentido es opuesto a la velocidad de deslizamiento sobre la barra. Como el problema específica que estudiemos el caso <math>\dot{x}>0 </math> (movimento hacia la derecha), la fuerza de rozamiento es
<center>
<math>
\vec{F}_R = -\mu|\vec{N}|\,\vec{\imath},
</math>
</center>
donde <math>\mu </math> es el coeficiente de rozamiento dinámico entre la partícula y la barra.

El efecto del peso y de la fuerza del muelle se incluyen en la función de Lagrange a través de la energía potencial.

Al aplicar el Principio de liberación, consideramos que el sistema tiene tres grados de libertad <math>\{ x, y, z \} </math>. Los vectores de posición y velocidad son
<center>
<math>
\begin{array}{l}
\vec{r} = x\,\vec{\imath} + y\,\vec{\jmath} + z\,\vec{k},
\\
\vec{v} = \dot{x}\,\vec{\imath} + \dot{y}\,\vec{\jmath} + \dot{z}\,\vec{k},
\end{array}
</math>
</center>
La energía cinética es
<center>
<math>
T = \dfrac{1}{2}m|\vec{v}|^2 = \dfrac{1}{2}m\,(\dot{x}^2 + \dot{y}^2 + \dot{z}^2).
</math>
</center>
Tomamos como referencia la altura del eje <math>X </math> para evaluar la energía pontencial gravitatoria
<center>
<math>
U_g = mgy
</math>
</center>
El muelle tiene longitud natural nula. En la situación liberada, su energía potencial es
<center>
<math>
U_m = \dfrac{1}{2}k\,|\overrightarrow{OP}|^2 = \dfrac{1}{2}\,k\,({x}^2 + {y}^2 + {z}^2).
</math>
</center>
La función de Lagrange es
<center>
<math>
L = T - (U_g + U_m) = \dfrac{1}{2}m\,(\dot{x}^2 + \dot{y}^2 + \dot{z}^2) - mgy -  \dfrac{1}{2}\,k\,({x}^2 + {y}^2 + {z}^2).
</math>
</center>
Las ecuaciones de Lagrange son
<center>
<math>
\begin{array}{lcl}
\dfrac{\mathrm{d}}{\mathrm{d}t} \left(\dfrac{\partial L}{\partial\dot{x}}\right) - \dfrac{\partial L}{\partial x} = Q_x^{NC},\\
\dfrac{\mathrm{d}}{\mathrm{d}t} \left(\dfrac{\partial L}{\partial\dot{y}}\right) - \dfrac{\partial L}{\partial y} = Q_y^{NC},\\
\dfrac{\mathrm{d}}{\mathrm{d}t} \left(\dfrac{\partial L}{\partial\dot{z}}\right) - \dfrac{\partial L}{\partial z} = Q_z^{NC}.
\end{array}
</math>
</center>
Las fuerzas generalizadas tiene contibuciones de las fuerzas vincular y de rozamiento. 
<center>
<math>
\begin{array}{l}
Q^{NC}_x = (\vec{N} + \vec{F}_R)\cdot\dfrac{\partial \vec{v}}{\partial \dot{x}} = -\mu\,\sqrt{N_y^2 + N_z^2},\\
Q^{NC}_y = (\vec{N} + \vec{F}_R)\cdot\dfrac{\partial \vec{v}}{\partial \dot{y}} = N_y,\\
Q^{NC}_z = (\vec{N} + \vec{F}_R)\cdot\dfrac{\partial \vec{v}}{\partial \dot{z}} = N_z.\\
\end{array}
</math>
</center>
Después de hacer las derivadas de la función de Lagrange, las ecuaciones de Lagrange quedan
<center>
<math>
\begin{array}{lr}
m\ddot{x} + kx = -\mu\sqrt{N_y^2 + N_z^2}, & (1)\\
m\ddot{y} + ky + mg = N_y, & (2)\\
m\ddot{z} + kz = N_z. & (3)
\end{array}
</math>
</center>
Tenemos cinco incógnitas <math>\{ x, y, z, N_y, N_z \} </math>. Las dos ecuaciones que faltan son las relaciones vinclares, así como sus derivadas y/o integrales.
<center>
<math>
\begin{array}{lllr}
y = h, & \dot{y} = 0, & \ddot{y}=0, & (4)\\
z = 0, & \dot{z} = 0, & \ddot{z}=0. & (5)
\end{array}
</math>
</center>
Aplicando las ecuaciones (4)-(5) en (1)-(3) obtenemos 
<center>
<math>
\begin{array}{l}
m\ddot{x} + kx = -\mu\,(mg + kh),\\
N_y = mg + kh,\\
N_z = 0.
\end{array}
</math>
</center>