Ver el código de «Estudio de la velocidad de tres puntos»

==Enunciado==
En un hipotético sólido rígido, las posiciones y velocidades de tres puntos son respectivamente:
<center><math>\begin{array}{rclcrcl}
\overrightarrow{OA}&=&\vec{\imath}&\qquad &
\vec{v}^A & = & 4\vec{\imath}+2\vec{k}\\
\overrightarrow{OB}&=&\vec{\jmath}&\qquad &
\vec{v}^B& = &2\vec{\imath}-2\vec{\jmath}-\vec{k}\\
\overrightarrow{OC}&=&\vec{k}&\qquad &
\vec{v}^C&=&2\vec{\imath}-\vec{\jmath}
\end{array}
</math></center>
# Demuestre que estas velocidades son compatibles con la condición de rigidez.
# Halle la velocidad del punto <math>O(0,0,0)</math>.
# Calcule la velocidad del punto <math>\overrightarrow{OP}=2\vec{\imath}-2\vec{\jmath}-\vec{k}</math>.
# ¿Existe algún punto que tenga velocidad nula? ¿Dónde estaría situado?
Todas las cantidades están expresadas en las unidades del SI.

==Condición de rigidez==
Para que un movimiento sea posible en un sólido rígido debe cumplirse, para todos y cada uno de sus pares de puntos la condición cinemática de rigidez, consistente en la equiproyectividad del campo de velocidades:

<center><math>\vec{v}^P\cdot\overrightarrow{PQ}=\vec{v}^Q\cdot\overrightarrow{PQ}</math></center>

Aplicando esta condición para cada par de partículas resulta:

;Partículas A y B:

<center><math>\overrightarrow{AB}=-\vec{\imath}+\vec{\jmath}\qquad\left\{\begin{matrix}\vec{v}^A\cdot\overrightarrow{AB} & = & (4\vec{\imath}+2\vec{k})\cdot(-\vec{\imath}+\vec{\jmath}) & = & -4 \\ \vec{v}^B\cdot\overrightarrow{AB} & = & (2\vec{\imath}-2\vec{\jmath}-\vec{k})\cdot(-\vec{\imath}+\vec{\jmath}) & = & -4 \end{matrix}\right.</math></center>

;Partículas A y C:

<center><math>\overrightarrow{AC}=-\vec{\imath}+\vec{k}\qquad\left\{\begin{matrix}\vec{v}^A\cdot\overrightarrow{AC} & = & (4\vec{\imath}+2\vec{k})\cdot(-\vec{\imath}+\vec{k}) & = & -2 \\ \vec{v}^C\cdot\overrightarrow{AC} & = & (2\vec{\imath}-\vec{\jmath})\cdot(-\vec{\imath}+\vec{k}) & = & -2 \end{matrix}\right.</math></center>

;Partículas B y C:

<center><math>\overrightarrow{BC}=-\vec{\jmath}+\vec{k}\qquad\left\{\begin{matrix}\vec{v}^B\cdot\overrightarrow{BC} & = & (2\vec{\imath}-2\vec{\jmath}-\vec{k})\cdot(-\vec{\jmath}+\vec{k}) & = & 1 \\ \vec{v}^C\cdot\overrightarrow{BC} & = & (2\vec{\imath}-\vec{\jmath})\cdot(-\vec{\jmath}+\vec{k}) & = & 1 \end{matrix}\right.</math></center>

Puesto que los tres pares verifican la condición, estas velocidades son compatibles con un movimiento de un sólido rígido.

==Velocidad del origen==
La velocidad de cualquier otro punto se puede obtener aplicando la equiproyectividad con los tres puntos que ya tenemos. Si la velocidad de O es

<center><math>\vec{v}^O = a\vec{\imath}+b\vec{\jmath}+c\vec{k}</math></center>

la condición cinemática de rigidez nos da, para el par formado por O y A

<center><math>\vec{v}^O\cdot\overrightarrow{OA}= \vec{v}^A\cdot\overrightarrow{OA}</math>{{tose}}<math>(a\vec{\imath}+b\vec{\jmath}+c\vec{k})\cdot\vec{\imath}=(4\vec{\imath}+2\vec{k})\cdot\vec{\imath}</math>{{tose}}<math>a=4\,</math></center>

Para O y B

<center><math>\vec{v}^O\cdot\overrightarrow{OB}=\vec{v}^B\cdot\overrightarrow{OB}</math>{{tose}}<math>(a\vec{\imath}+b\vec{\jmath}+c\vec{k})\cdot\vec{\jmath}= (2\vec{\imath}-2\vec{\jmath}-\vec{k})\cdot\vec{\jmath}</math>{{tose}}<math>b=-2\,</math></center>

Y para O y C

<center><math>\vec{v}^O\cdot\overrightarrow{OC}=\vec{v}^C\cdot\overrightarrow{OC}</math>{{tose}}<math>(a\vec{\imath}+b\vec{\jmath}+c\vec{k})\cdot\vec{k}= (2\vec{\imath}-\vec{\jmath})\cdot\vec{k}</math>{{tose}}<math>c=0\,</math></center>

Reuniendo los tres resultados obtenemos la velocidad

<center><math>\vec{v}^O = (4\vec{\imath}-2\vec{\jmath})\,\mathrm{m/s}</math></center>

==Velocidad del punto P==
El mismo procedimiento puede emplearse para cualquier otro punto, si bien en general resultará un sistema de ecuaciones lineales. Si la velocidad del punto P es de la forma

<center><math>\vec{v}^P = d\vec{\imath}+e\vec{\jmath}+f\vec{k}</math></center>

resultan las ecuaciones

Para P y A

<center><math>\overrightarrow{AP}=\overrightarrow{OP}-\overrightarrow{OA}=\vec{\imath}-2\vec{\jmath}-\vec{k}</math>{{qquad}}{{qquad}}<math>\vec{v}^P\cdot\overrightarrow{AP}= \vec{v}^A\cdot\overrightarrow{AP}</math>{{tose}}<math>d-2e-f=2\,</math></center>

Para P y B

<center><math>\overrightarrow{BP}=\overrightarrow{OP}-\overrightarrow{OB}=2\vec{\imath}-3\vec{\jmath}-\vec{k}</math>{{qquad}}{{qquad}}<math>\vec{v}^P\cdot\overrightarrow{BP}= \vec{v}^B\cdot\overrightarrow{BP}</math>{{tose}}<math>2d-3e-f=11\,</math></center>

Para P y C

<center><math>\overrightarrow{CP}=\overrightarrow{OP}-\overrightarrow{OC}=2\vec{\imath}-2\vec{\jmath}-2\vec{k}</math>{{qquad}}{{qquad}}<math>\vec{v}^P\cdot\overrightarrow{CP}= \vec{v}^C\cdot\overrightarrow{CP}</math>{{tose}}<math>2d-2e-2f=6\,</math></center>

Tenemos entonces el sistema de ecuaciones

<center><math>\left\{\begin{matrix}d & - & 2e & - & f & = & 2 \\ 2d& - & 3e & - & f & = & 11 \\ 2d & - &2e & -& 2f &= & 6\end{matrix}\right.</math></center>

con solución

<center><math>d = 10\qquad e = 1\qquad f = 6\qquad\vec{v}^P = (10\,\vec{\imath}+\vec{\jmath}+6\,\vec{k})\,\mathrm{m/s}</math></center>

Un procedimiento alternativo para calcular la velocidad de un cuarto punto del sólido rígido, conocidas las velocidades de tres puntos no alineados, consiste en determinar primero el vector velocidad angular <math>\vec{\omega}\,</math>. Las velocidades conocidas están obligadas a satisfacer la ecuación del campo de velocidades. Así, por ejemplo:
<center><math>
\vec{v}^{B}=\vec{v}^{A}+\,\vec{\omega}\times\overrightarrow{AB}\,\,\,\Longrightarrow\,\,\,2\vec{\imath}\, -\, 2\,\vec{\jmath}\, -\,\vec{k}=4\,\vec{\imath}\, +\,2\,\vec{k}\,+\left|\begin{array}{ccc} \vec{\imath} & \vec{\jmath} & \vec{k} \\ \omega_x & \omega_y & \omega_z \\ -1 & 1 & 0 \end{array}\right|\,\,\,\Longrightarrow\,\,\,\left\{\begin{array}{lllll} 2 & = & 4-\omega_z & & \\ -2 & = & -\omega_z &\,\,\,\longrightarrow\,\,\,& \omega_z=2 \\ -1 & = & 2+\omega_x+\omega_y & \,\,\,\longrightarrow\,\,\, & \omega_y=-3-\omega_x \end{array} \right.
</math></center>

Y por otra parte:
<center><math>
\vec{v}^{C}=\vec{v}^{A}+\,\vec{\omega}\times\overrightarrow{AC}\,\,\,\Longrightarrow\,\,\,2\vec{\imath}\, -\, \vec{\jmath}=4\,\vec{\imath}\, +\,2\,\vec{k}\,+\left|\begin{array}{ccc} \vec{\imath} & \vec{\jmath} & \vec{k} \\ \omega_x & -3-\omega_x & 2 \\ -1 & 0 & 1 \end{array}\right|\,\,\,\Longrightarrow\,\,\,\left\{\begin{array}{lllll} 2 & = & 4-3-\omega_x & \,\,\,\longrightarrow\,\,\, & \omega_x=-1 \\ ... & = & ... & &  \\ ... & = & ... &  &  \end{array} \right.
</math></center>
Ya tenemos información suficiente para escribir el vector velocidad angular:
<center><math>
\vec{\omega}=(-\vec{\imath}-2\,\vec{\jmath}+\,2\vec{k})\,\mathrm{rad/s}
</math></center>
Y ahora, una vez conocido <math>\vec{\omega}\,</math>, se puede calcular la velocidad de cualquier punto utilizando de nuevo la ecuación del campo de velocidades:
<center><math>
\vec{v}^{P}=\vec{v}^{O}+\,\vec{\omega}\times\overrightarrow{OP}\,\,\,\Longrightarrow\,\,\,\vec{v}^{P}=4\,\vec{\imath}\, -\,2\,\vec{\jmath}\,+\left|\begin{array}{ccc} \vec{\imath} & \vec{\jmath} & \vec{k} \\ -1 & -2 & 2 \\ 2 & -2 & -1 \end{array}\right|\,\,\,\Longrightarrow\,\,\,\vec{v}^{P}=(10\,\vec{\imath}+\vec{\jmath}+6\,\vec{k})\,\mathrm{m/s}
</math></center>

==Puntos con velocidad instantánea nula==
===A partir de la ecuación del campo de velocidades===
Se trata de buscar un punto (o un conjunto de ellos), I, tales que:
<center><math>\overrightarrow{OI}=x\vec{\imath}+y\vec{\jmath}+z\vec{k}\,\,;</math>{{qquad}}{{qquad}}<math>\vec{v}^I = \vec{0}</math></center>

Apoyándonos en la velocidad ya conocida del punto <math>O\,</math>, podemos escribir:
<center><math>
\vec{v}^I=\vec{v}^O+\,\vec{\omega}\times\overrightarrow{OI}\,\,\,\Rightarrow\,\,\,\vec{0}=4\,\vec{\imath}\, -\,2\,\vec{\jmath}\,+\left|\begin{array}{ccc} \vec{\imath} & \vec{\jmath} & \vec{k} \\ -1 & -2 & 2 \\ x & y & z \end{array}\right|\,\,\,\Rightarrow\,\,\,\left\{\begin{array}{lllll} 0 & = & 4-2z-2y & \,\rightarrow\, & 2y+2z=4 \\ 0 & = & -2+2x+z &\,\rightarrow\,& 2x+z=2 \\ 0 & = & -y+2x & \,\rightarrow\, & 2x-y=0 \end{array} \right.
</math></center>
Por tanto, el sistema de ecuaciones cuya solución nos da la posición de los puntos de velocidad instantánea nula es:

<center><math>\left\{\begin{matrix} &  & 2y & + & 2z & = & 4 \\ 2x&  &  & + & z & = & 2 \\ 2x & - &y & &  &= & 0\end{matrix}\right.</math></center>

Este sistema es compatible indeterminado. Podemos prescindir de la primera ecuación porque es linealmente dependiente de las otras dos (se obtiene restándole la tercera ecuación multiplicada por 2 a la segunda ecuación multiplicada por 2). Podemos quedarnos, pues, con el sistema de dos ecuaciones:

<center><math>\left\{\begin{matrix}  2x+z & = & 2 & \,\,\,\longrightarrow\,\,\, & z & = & 2-2x\\ 2x-y & = & 0 & \,\,\,\longrightarrow\,\,\, & y & = & 2x \end{matrix}\right.</math></center>

Por tanto, obtenemos un conjunto infinito de soluciones, situadas sobre la recta

<center><math>\left\{\begin{matrix}x & = & \mu \\ y & = & 2\mu \\ z & = & 2-2\mu\end{matrix}\right.</math>{{tose}}<math>\overrightarrow{OI}=[2\vec{k}+\mu(\vec{\imath}+2\vec{\jmath}-2\vec{k})]\,\mathrm{m}</math></center>

===Cálculo del eje central===
Otra forma más directa de llegar a la solución consiste en caer en la cuenta de que, si existen puntos de velocidad instantánea nula, seguro que esos puntos pertenecen al eje central del campo de velocidades (lugar geométrico de los puntos cuya velocidad tiene módulo mínimo). Por tanto, lo primero que hacemos es confirmar que efectivamente, en el ejercicio que nos ocupa, la velocidad en los puntos del eje central es nula:
<center><math>\vec{v}^{\,\mathrm{min}}=\vec{v}_d=v_d\,\vec{u}_{\omega}=v_d\,\frac{\vec{\omega}}{|\vec{\omega}|}=\frac{\vec{v}_A\cdot\vec{\omega}}{|\vec{\omega}|}\,\frac{\vec{\omega}}{|\vec{\omega}|}=\vec{0}
</math></center>
Y localizamos dichos puntos mediante la ecuación vectorial del eje central:
<center><math>
\overrightarrow{OI}=\frac{\vec{\omega}\times\vec{v}^{\, O}}{|\vec{\omega}|^2}+\lambda\,\vec{\omega}=\left[\frac{4}{9}\,\vec{\imath}+\frac{8}{9}\,\vec{\jmath}+\frac{10}{9}\,\vec{k}+\lambda\,(-\vec{\imath}-2\,\vec{\jmath}+2\,\vec{k})\right]\,\mathrm{m}
</math></center>
Esta ecuación parece diferente de la obtenida por el procedimiento anterior, pero en realidad ambas corresponden a la misma recta. En efecto, si hacemos el cambio de parámetro:

<center><math>\mu = \frac{4}{9}-\lambda</math></center>

podemos comprobar cómo la anterior ecuación se convierte en ésta de ahora.

<!--
===A partir del campo de velocidades===
La condición cinemática de rigidez también puede emplearse para hallar posiciones, en lugar de velocidades. En este caso buscamos un punto (o un conjunto de ellos), I, tales que

<center><math>\overrightarrow{OI}=x\vec{\imath}+y\vec{\jmath}+z\vec{k}</math>{{qquad}}{{qquad}}<math>\vec{v}^I = \vec{0}</math></center>

Imponiendo la equiproyectividad con las velocidades de A, B y C tenemos

Para A e I

<center><math>\overrightarrow{AI}=\overrightarrow{OI}-\overrightarrow{OA}=(x-1)\vec{\imath}+y\vec{\jmath}+z\vec{k}</math>{{qquad}}{{qquad}}<math>\vec{v}^I\cdot\overrightarrow{AI}= \vec{v}^A\cdot\overrightarrow{AI}</math>{{tose}}<math>0 =4(x-1)+2z\,</math></center>

Para B e I

<center><math>\overrightarrow{BI}=\overrightarrow{OI}-\overrightarrow{OB}=x\vec{\imath}+(y-1)\vec{\jmath}+z\vec{k}</math>{{qquad}}{{qquad}}<math>\vec{v}^I\cdot\overrightarrow{BI}= \vec{v}^B\cdot\overrightarrow{BI}</math>{{tose}}<math>0 =2x-2(y-1)-z\,</math></center>

Para C e I

<center><math>\overrightarrow{CI}=\overrightarrow{OI}-\overrightarrow{OC}=x\vec{\imath}+y\vec{\jmath}+(z-1)\vec{k}</math>{{qquad}}{{qquad}}<math>\vec{v}^I\cdot\overrightarrow{CI}= \vec{v}^C\cdot\overrightarrow{CI}</math>{{tose}}<math>0 =2x-y\,</math></center>

El sistema de ecuaciones para la posición es

<center><math>\left\{\begin{matrix}4x & &  & + & 2z & = & 4 \\ 2x& - & 2y & - & z & = & -2 \\ 2x & - &y & &  &= & 0\end{matrix}\right.</math></center>

Este sistema no siempre tiene solución, ya que si hallamos el rango de los coeficientes, resulta menor que el número de ecuaciones, ya que

<center><math>\left|\begin{matrix} 4 & 0 & 2 \\ 2 & -2 & -1 \\ 2 & -1 & 0 \end{matrix}\right| = \vec{v}^A\cdot(\vec{v}^B\times\vec{v}^C) = 0</math></center>

por lo que, o no obtenemos solución alguna o resultan infinitas soluciones. 

Despejando en la última ecuación

<center><math>
y = 2x\,</math>{{tose}}<math>\left\{\begin{matrix} 4x & + & 2z & = & 4 \\ -2x & - & z & = & -2\end{matrix}\right.</math>{{tose}}<math>z = 2-2x\,</math></center>

Por tanto obtenemos un conjunto infinito de soluciones, situadas sobre la recta

<center><math>\left\{\begin{matrix}x & = & \lambda \\ y & = & 2\lambda \\ z & = & 2-2\lambda\end{matrix}\right.</math>{{tose}}<math>\overrightarrow{OI}=2\vec{k}+\lambda(\vec{\imath}+2\vec{\jmath}-2\vec{k})</math></center>

===Cálculo del EIR===

Un procedimiento sistemático para hallar los puntos de velocidad nula es determinar la posición del eje instantáneo de rotación y mínimo deslizamiento (EIRMD) empleando las fórmulas generales.

Para ello precisamos calcular en primer lugar la reducción <math>\{\vec{v}^O,\vec{\omega}\}</math> del campo de velocidades, que nos permite hallar la velocidad de cualquier otro punto como

<center><math>\vec{v}^P = \vec{v}^O + \vec{\omega}\times\overrightarrow{OP}</math></center>

La velocidad del punto O (que es un punto arbitrario, pero para el que por comodidad tomamos el origen de coordenadas) ya la conocemos. Debemos determinar la velocidad angular del movimiento. Hacemos esto, observando que las velocidades de O, A y B se relacionan como

<center><math>\vec{v}^A - \vec{v}^O = \vec{\omega}\times\overrightarrow{OA}</math>{{qquad}}{{qquad}}<math>\vec{v}^B - \vec{v}^O = \vec{\omega}\times\overrightarrow{OB}</math></center>

y por tanto la velocidad angular es ortogonal a las dos diferencias de velocidades, por lo que puede escribirse como

<center><math>\vec{\omega}=\lambda(\vec{v}^A - \vec{v}^O)\times(\vec{v}^B - \vec{v}^O)=\lambda\left|\begin{matrix}\vec{\imath} & \vec{\jmath} & \vec{k} \\ 0 & 2 & 2 \\ -2 & 0 & -1\end{matrix}\right|=-2\lambda(\vec{\imath}+2\vec{\jmath}-2\vec{k})</math></center>

El valor del parámetro <math>\lambda</math> lo obtenemos a partir de alguno de los pares concretos

<center><math>\vec{v}^A -\vec{v}^O = 2\vec{\jmath}+2\vec{k}=\vec{\omega}\times\overrightarrow{OA}=-2\lambda(\vec{\imath}+2\vec{\jmath}-2\vec{k})\times\vec{\imath}= 4(\vec{\jmath}+\vec{k})</math></center>

de donde

<center><math>\lambda = \frac{1}{2}\qquad\Rightarrow\qquad \vec{\omega}=-\vec{\imath}-2\vec{\jmath}+2\vec{k}</math></center>

Una vez que conocemos la reducción en el punto O, podemos identificar el movimiento hallando la velocidad de deslizamiento

<center><math>v_d = \frac{\vec{v}^O\cdot\vec{\omega}}{\omega}=\frac{(4\vec{\imath}-2\vec{\jmath})\cdot(-\vec{\imath}-2\vec{\jmath}+2\vec{k})}{3}=0</math></center>

Por ser nula, se trata de un movimiento de rotación pura. La posición del eje instantáneo de rotación, EIR, se halla mediante la [[Cinemática_del_sólido_rígido_(G.I.T.I.)#An.C3.A1lisis_del_caso_general|fórmula]]

<center><math>\overrightarrow{OI}=\frac{\vec{\omega}\times\vec{v}^O}{\omega^2} + \mu\vec{\omega} = \frac{1}{9}\left|\begin{matrix}\vec{\imath} & \vec{\jmath} & \vec{k} \\ -1 & -2 & 2 \\ 4 & -2 & 0\end{matrix}\right|+\mu(-\vec{\imath}-2\vec{\jmath}+2\vec{k})=\left(\frac{4}{9}-\mu\right)\vec{\imath}+\left(\frac{8}{9}-2\mu\right)\vec{\jmath}+\left(\frac{10}{9}+2\mu\right)\vec{k}</math></center>

Esta ecuación parece diferente de la obtenida anteriormente para los puntos de velocidad nula, pero representa a la misma recta. Si hacemos

<center><math>\mu = \frac{4}{9}-\lambda</math></center>

se convierte en

<center><math>\overrightarrow{OI}=\lambda\vec{\imath}+2\lambda\vec{\jmath}+(2-2\lambda)\vec{k}=2\vec{k}+\lambda\left(\vec{\imath}+2\vec{\jmath}-2\vec{k}\right)</math></center>

que es la misma a la que llegamos antes.
-->

[[Categoría:Problemas de Cinemática del Sólido Rígido (GITI)]]