Ver código fuente

== Enunciado ==
[[Imagen:F1A_GIA_muelle_esfera_enunciado.png|right]]
Un punto material <math>M</math> de peso <math>P</math> está obligado a permanecer en la superficie de una esfera de radio <math>R</math> y centro <math>O</math>. Además, <math>M</math> es atraído por un punto fijo <math>A</math> del ecuador de la superficie esférica, debido a la existencia de un resorte elástico ideal, de longitud natural nula y de constante recuperadora <math>k=P/\sqrt{3}R</math>, que conecta ambos puntos. Determina las posiciones de equilibrio del punto material <math>M</math>, y la fuerza de reacción vincular en ellas.

==Solución==

=== Ecuaciones de equilibrio ===
Se trata de determinar la posición de equilibrio de un punto material
que está obligado a moverse en una superficie dada. Para una partícula
libre, el vector que indica su posición en el espacio se escribe, en
general,
<center><math>
  \vec{r}=x\,\vec{\imath}+y\,\vec{\jmath}+z\,\vec{k}
</math></center>
En este caso, la superficie impone una ligadura bilateral: la
ecuación que cumplen los puntos situados sobre una superficie esférica
de radio <math>R</math>
<center><math>
  x^2+y^2+z^2=R^2\,\Longrightarrow\,z=\pm\sqrt{R^2-x^2-y^2}
</math></center>

La condición de equilibrio mecánico para un punto material es que la
suma de fuerzas sobre él sea nula. Aquí tenemos tres fuerzas,
el peso <math>\vec{P}</math>, la fuerza ejercida por el muelle, <math>\vec{F}_k</math> y la
fuerza de reacción vincular que impide que el punto penetre en la
esfera, <math>\vec{\Phi}</math>.  Veamos la expresión de cada una de ellas.

El peso apunta en la dirección negativa del eje <math>OZ</math> que hemos
elegido, es decir
<center><math>
  \vec{P} = -P\,\vec{k}
</math></center>


La fuerza <math>\vec{F}_k</math> es proporcional a la distancia entre el punto y el
punto de sujeción del muelle, pues el enunciado dice que su longitud
natural es cero. La posición de este punto es <math>\vec{r}_A=R\,\vec{\imath}</math>.  Por
tanto
<center><math>
  \vec{F}_k=-k\,|\vec{r}-\vec{r}_A| = -k\,(x-R)\,\vec{\imath}-k\,y\,\vec{\jmath}-k\,z(x,y)\,\vec{k}
</math></center>
Recordemos que <math>z</math> viene dada en función de <math>x</math> e <math>y</math> a través de la
ecuación de la esfera.

Por ultimo la fuerza de reacción vincular <math>\vec{\Phi}</math> es normal a
la superficie, pues suponemos que el contacto es liso. Tenemos
<center><math>
  \vec{\Phi} = \Phi\dfrac{\vec{r}}{|\vec{r}|}=
  \dfrac{\Phi}{R}\,\vec{r}=\dfrac{\Phi}{R}\,(x\,\vec{\imath}+y\,\vec{\jmath}+z(x,y)\,\vec{k})
</math></center>
Hemos usado que en cada punto de la esfera podemos construir un vector
unitario perpendicular a ella dividiendo el vector de posición de ese
punto por su módulo, <math>\vec{r}/|\vec{r}|</math>.  Y para todos los puntos en la
superficie <math>|\vec{r}|=R</math>.  Para resolver el problema vamos a suponer que el vínculo es no rugoso y que la
partícula se mueve '''sobre''' la superficie de la esfera. Por tanto, la f.r.v. sólo puede estar
dirigida hacia el exterior de la superficie. Si estuviera dirigida hacia dentro querría decir que habría
una fuerza adhesiva entre la superficie y la partícula, con lo que el vinculo no sería liso.


Ahora podemos plantear la condición de equilibrio para el punto
<center><math>
  \vec{P}+\vec{F}_K+\vec{\Phi}=\vec{0}
</math></center>
Igualando a cero las componentes de esta ecuación vectorial obtenemos
tres ecuaciones escalares, a las que añadimos la ecuación que deben cumplir las coordenadas de los puntos sobre una esfera
de radio <math>R </math>
<center><math>
  \begin{array}{rcl}
    -k(x-R)+\dfrac{\Phi}{R}x &=&0\\ &&\\
    -k\,y+\dfrac{\Phi}{R}y &=&0\\ &&\\
    -k\,z+\dfrac{\Phi}{R}z-P &=&0 \\ &&\\
    z^2 &=& R^2 - x^2 - y^2 
  \end{array}
  \quad\Longrightarrow\quad
  \begin{array}{rcl}
    \left(k-\dfrac{\Phi}{R}\right)\,x&=&k\,R\\ &&\\
    \left(k-\dfrac{\Phi}{R}\right)\,y&=&0\\ &&\\
    \left(k-\dfrac{\Phi}{R}\right)\,z&=&-P\\ &&\\
    z &=&\pm\sqrt{ R^2 - x^2 - y^2 }
  \end{array}
</math></center>

=== Análisis de las soluciones ===
Las incógnitas del problema son <math>\Phi</math>, <math>x</math>, <math>y</math> y <math>z </math>.  En la segunda ecuación
<math>(k-\Phi/R)</math> no puede ser cero, pues entonces, de la primera ecuación
tendríamos <math>x\to\infty</math>.  Así pues 
<center><math>
  y=0
</math></center>
Dividiendo la primera ecuación por la tercera y teniendo en cuenta
la ecuación de la esfera obtenemos
<center><math>
  x=\pm \dfrac{R}{\sqrt{1+(P/k\,R)^2}} = \pm\dfrac{R}{\sqrt{1+C^2}}
</math></center>
siendo el parámetro <math>C </math>
<center>
<math>
C = \dfrac{P}{KR}
</math>
</center>
Este parámetro mide la importancia relativa del peso y la fuerza ejercida por el muelle. Si <math>C\gg 1 </math> la fuerza del 
muelle es despreciable frente a la del peso. Si <math>C\ll1 </math> es el peso el que es despreciable frente a la fuerza del 
muelle.

Para cada valor de <math>x </math> tenemos dos posibles valores de <math>z </math> (recordemos que <math>y=0 </math>) 
<center>
<math>
z = \pm\sqrt{R^2-x^2} = \pm\dfrac{C\,R}{\sqrt{1+C^2}}
</math>
</center>
Es decir, hay cuatro posibles puntos de 
equilibrio. 

A su vez, la componente de la fuerza de reacción vincular puede despejarse de la primera y de la tercera ecuación de equilibrio.
Obtenemos entonces dos posibles valores de <math>\Phi </math> para cada par de valores de <math>x </math> y <math>z </math>.
Para que el punto sea solución de equilibrio estos dos valores deben coincidir. Resumimos las posibilidades en la siguiente
tabla
<center>
{| border="1" style="text-align:center;background:#ffffff" align="center" 
! width="60 px" style="background:Lavender; color:Black"|Punto 
! width="80 px" style="background:Lavender; color:Black"|<math>x </math>
! width="100 px" style="background:Lavender; color:Black"|<math>z </math>
! width="100 px" style="background:Lavender; color:Black"|<math>\Phi </math> de <math>x </math>
! width="100 px" style="background:Lavender; color:Black"|<math>\Phi </math> de <math>z </math>
! width="100 px" style="background:Lavender; color:Black"|Posible
|-
|<math>P_1 </math>||<math>+\dfrac{R}{\sqrt{1+C^2}} </math>||<math>+\dfrac{C\,R}{\sqrt{1+C^2}} </math>
||<math>KR(1-\sqrt{1+C^2})</math>||<math>KR(1+\sqrt{1+C^2})</math>
||NO
|-
|<math>P_2 </math>||<math>+\dfrac{R}{\sqrt{1+C^2}} </math>||<math>-\dfrac{C\,R}{\sqrt{1+C^2}} </math>
||<math>KR(1-\sqrt{1+C^2})</math>||<math>KR(1-\sqrt{1+C^2})</math>
||SÍ
|-
|<math>P_3 </math>||<math>-\dfrac{R}{\sqrt{1+C^2}} </math>||<math>+\dfrac{C\,R}{\sqrt{1+C^2}} </math>
||<math>KR(1+\sqrt{1+C^2})</math>||<math>KR(1+\sqrt{1+C^2})</math>
||SÍ
|-
|<math>P_4 </math>||<math>-\dfrac{R}{\sqrt{1+C^2}} </math>||<math>-\dfrac{C\,R}{\sqrt{1+C^2}} </math>
||<math>KR(1+\sqrt{1+C^2})</math>||<math>KR(1-\sqrt{1+C^2})</math>
||NO
|}
</center>
Sólo los puntos <math>P_2 </math> y <math>P_3 </math> dan el mismo valor de <math>\Phi </math> al despejarla
de la primera y tercera ecuaciones de equilibrio. 

En el caso del punto <math>P_2 </math>, como el parámetro <math>C </math> es siempre positivo tenemos
<center>
<math>
\Phi^{P_2} = KR(1-\sqrt{1+C^2}) < 0
</math>
</center>
Este signo significa que la fuerza de reacción vincular apuntaría hacia el centro de la esfera. Pero si el vínculo
es liso esto no es posible. La partícula esta '''sobre''' la superficie de la esfera, no pegada a ella. Una fuerza
hacia dentro supondría que hay algún tipo de pegamento entre la partícula y la esfera. Por tanto, esta solución no es posible.

[[Imagen:F1A_GIA_muelle_esfera_fuerzas.png|right]]
Así pues, el único punto de equilibrio es el <math>P_3 </math>. Las coordenadas y la componente de la fuerza de reacción
vincular son
<center>
{| border="1" style="text-align:center;background:#ffffff" align="center" 
! width="80 px" style="background:Lavender; color:Black"|<math>x </math>
! width="80 px" style="background:Lavender; color:Black"|<math>y </math>
! width="100 px" style="background:Lavender; color:Black"|<math>z </math>
! width="100 px" style="background:Lavender; color:Black"|<math>\Phi </math>
|-
|<math>-\dfrac{R}{\sqrt{1+C^2}} </math>||0 ||<math>+\dfrac{C\,R}{\sqrt{1+C^2}} </math>
||<math>KR(1+\sqrt{1+C^2})</math>
|}
</center>
con <math>C = P/KR </math>.


Puede verse gráficamente que el único punto posible de equilibrio es el <math>P_3 </math>. La figura muestra
la dirección de las fuerzas en cada uno de los puntos de equilibrio, teniendo en cuenta que la fuerza de reacción
vincular sólo puede apuntar hacia el exterior de las esferas. El único punto en el que la suma vectorial de 
las tres fuerzas puede dar cero es en el <math>P_3</math>. En el <math>P_1 </math> y el <math>P_4 </math>, si proyectamos
las fuerzas en las direcciones paralela y perpendicular a <math>\vec{\Phi} </math>, vemos que siempre queda una 
componente no nula perpendicular a <math>\vec{\Phi} </math>. Lo mismo ocurre en <math>P_2 </math> si proyectamos en las direcciones 
paralela y perpendicular a <math>\vec{F}_k </math>.

===Caso límite cuando el muelle es muy fuerte===
Vamos a considerar la situación en que se cumple
<center><math>
  k\,R\gg P \Longrightarrow \dfrac{P}{k\,R}\ll 1
</math></center>
Esto ocurre si la constante del muelle es muy grande, es decir, cuando
tenemos un muelle muy fuerte. Para obtener el valor aproximado
usaremos el desarrollo de Taylor
<center><math>
  (1+\delta^n)\simeq 1 + n\,\delta\qquad\mathrm{(cuando\quad\delta\to0)}
</math></center>
En las soluciones para <math>x </math> y <math>\Phi </math>  consideramos <math>\delta=(P/k\,R)^2</math>,
con lo cual  los valores aproximados de <math>x</math> y <math>\Phi</math> en este caso límite
<center><math>
  x\simeq
  -R\left(1-\dfrac{1}{2}\left(\dfrac{P}{k\,R}\right)^2\right)\to -R
  \qquad\qquad
  \Phi\simeq
  k\,R\,\left(2+\dfrac{1}{2}\left(\dfrac{P}{k\,R}\right)^2\right)\to 2\,k\,R
</math></center>
Es decir, si el muelle tira muy fuerte la posición de equilibrio de la
masa es el punto opuesto de la esfera respecto al punto de sujeción
del muelle, y la fuerza de reacción vincular depende responde sólo a
la fuerza del muelle, pues entonces la contribución del peso es despreciable.


[[Categoría:Dinámica del punto material|1]]
[[Categoría:Física I (G.I.A.)]]
[[Categoría:Física I (G.I.T.I.)]]
[[Categoría:Física I (G.I.C.)]]