geometrÃa de masas - MecFunNet
geometrÃa de masas - MecFunNet
geometrÃa de masas - MecFunNet
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
UPM<br />
DFAII<br />
ETSII, UPM<br />
Geometría <strong>de</strong> <strong>masas</strong><br />
M ɛ cFunN ɛ t<br />
Resumen<br />
Contenido<br />
Página Inicial<br />
◭◭<br />
◭<br />
◮◮<br />
◮<br />
Volver<br />
Cerrar<br />
Salir
Contenido<br />
1 Características <strong>de</strong> inercia <strong>de</strong> un sólido 3<br />
2 Distribuciones <strong>de</strong> <strong>masas</strong> 4<br />
3 Momentos estáticos 6<br />
4 Momentos <strong>de</strong> inercia 11<br />
5 Tensor y elipsoi<strong>de</strong> <strong>de</strong> inercia 17<br />
6 Simetrías en la distribución 22<br />
7 Ejes <strong>de</strong> los elipsoi<strong>de</strong>s <strong>de</strong> inercia 24<br />
8 Clasificación <strong>de</strong> las rectas 26<br />
9 Clasificación <strong>de</strong> los planos 32<br />
10 Clasificación <strong>de</strong> las direcciones 36<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
◮◮<br />
◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 2 <strong>de</strong> 41
1. Características <strong>de</strong> inercia <strong>de</strong> un sólido<br />
Consi<strong>de</strong>rando la mecánica como la ciencia que se ocupa <strong>de</strong>l estudio <strong>de</strong> la<br />
evolución <strong>de</strong> los sistemas materiales y las causas que la producen, po<strong>de</strong>mos<br />
preguntarnos sobre los aspectos o parámetros <strong>de</strong>l sólido que tienen transcen<strong>de</strong>ncia<br />
en el ámbito <strong>de</strong> la mecánica. Si centramos el objeto <strong>de</strong> nuestro<br />
estudio en el sólido rígido, entonces su evolución viene <strong>de</strong>terminada por la<br />
cinemática <strong>de</strong> los sistemas in<strong>de</strong>formables. A<strong>de</strong>más, a<strong>de</strong>lantando un resultado<br />
<strong>de</strong> Dinámica, las causas que producen la evolución <strong>de</strong>l sólido rígido<br />
pue<strong>de</strong>n ser encajadas en los sistemas <strong>de</strong> vectores <strong>de</strong>slizantes. Queda, por<br />
tanto, el estudio <strong>de</strong> la relación entre el movimiento <strong>de</strong>l sólido y sus causas.<br />
En esta relación intervienen las propias características <strong>de</strong> cada sólido y, <strong>de</strong><br />
éstas, sólo los momentos <strong>de</strong> ór<strong>de</strong>nes cero, uno y dos. Dicho <strong>de</strong> otra forma,<br />
la única información necesaria para relacionar los sistemas <strong>de</strong> fuerzas que<br />
actúan sobre un sólido con la evolución posterior <strong>de</strong> éste son los momentos<br />
<strong>de</strong> or<strong>de</strong>n hasta el segundo <strong>de</strong> dicho sólido. En [1] se analiza ampliamente<br />
esta relación y en estas notas se presenta y <strong>de</strong>sarrolla el tratamiento <strong>de</strong><br />
estos momentos, que es lo que se conoce tradicionalmente como Geometría<br />
<strong>de</strong> <strong>masas</strong>.<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
◮◮<br />
◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 3 <strong>de</strong> 41
2. Distribuciones <strong>de</strong> <strong>masas</strong><br />
Des<strong>de</strong> el punto <strong>de</strong> vista <strong>de</strong> la mecánica racional, un sólido rígido es un ente<br />
i<strong>de</strong>al compuesto por un conjunto <strong>de</strong> puntos materiales que mantienen sus<br />
distancias mutuas durante su evolución. Dicho conjunto <strong>de</strong> puntos materiales<br />
pue<strong>de</strong> estar compuesto por <strong>masas</strong> finitas, constituyendo un conjunto<br />
discreto, por distribuciones volúmicas, superficiales o lineales, en cuyo caso<br />
recibe el nombre <strong>de</strong> continuo, o bien, pue<strong>de</strong>n existir puntos <strong>de</strong> masa finita<br />
y distribuciones continuas en un mismo sólido.<br />
En sistemas discretos, pue<strong>de</strong> suponerse la existencia <strong>de</strong> un conjunto<br />
finito <strong>de</strong> n puntos P i <strong>de</strong> masa m i posicionados respecto a un origen O por<br />
su vector <strong>de</strong> posición OP i = r i i ∈ {1 . . . n}. La masa total, M <strong>de</strong>l sólido<br />
es<br />
n∑<br />
M =<br />
i=1<br />
En sistemas continuos no pue<strong>de</strong> hablarse <strong>de</strong> la existencia <strong>de</strong> una masa<br />
finita en cada punto, sino <strong>de</strong> una <strong>de</strong>nsidad <strong>de</strong> masa, ya sea volúmica (ρ),<br />
superficial (σ) o lineal (λ). Según sea el caso, la masa total vendrá dada<br />
por un integral volúmica, <strong>de</strong> superficie o <strong>de</strong> línea:<br />
m i<br />
M = ∫ ρdV<br />
M = ∫ σdS<br />
M = ∫ λdl<br />
La <strong>de</strong>finición <strong>de</strong> momentos <strong>de</strong> distribuciones <strong>de</strong> <strong>masas</strong> utiliza sumatorios<br />
o integrales <strong>de</strong> los tipos anteriores para exten<strong>de</strong>r una suma a un conjunto<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 4 <strong>de</strong> 41
discreto o continuo <strong>de</strong> puntos. Con objeto <strong>de</strong> no repetir <strong>de</strong>finiciones en cada<br />
uno <strong>de</strong> los cuatro casos anteriores, se asumirá una distribución volúmica.<br />
Los sumatorios y las integrales <strong>de</strong> superficie y <strong>de</strong> línea pue<strong>de</strong>n reducirse a<br />
integrales volúmicas con la introducción <strong>de</strong> funciones basadas en la δ <strong>de</strong><br />
Dirac.<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
Título<br />
Contenido<br />
◭◭<br />
◭<br />
◮◮<br />
◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 5 <strong>de</strong> 41
3. Momentos estáticos<br />
Dado un punto O, que se toma como origen <strong>de</strong> coor<strong>de</strong>nadas, se <strong>de</strong>fine como<br />
momento estático M O respecto a O al vector:<br />
∫<br />
∫<br />
M O = OP ρdV = rρdV<br />
Una distribución cualquiera <strong>de</strong> <strong>masas</strong>, por tanto, <strong>de</strong>fine un campo <strong>de</strong><br />
vectores M O en cada punto O <strong>de</strong>l espacio. Veamos a continuación la<br />
relación entre los momentos estáticos <strong>de</strong> una distribución <strong>de</strong> <strong>masas</strong> respecto<br />
a dos puntos O y O ′ cualesquiera.<br />
∫<br />
∫<br />
M O ′ = O ′ P ρdV = (O ′ O + OP )ρdV =<br />
∫<br />
−<br />
∫<br />
OO ′ ρdV +<br />
OP ρdV = M O − MOO ′ (1)<br />
es <strong>de</strong>cir, la diferencia entre los momentos estáticos <strong>de</strong> una distribución<br />
respecto a dos puntos O y O ′ cualesquiera es el producto <strong>de</strong> la masa por<br />
el vector OO ′ .<br />
Se <strong>de</strong>fine centro <strong>de</strong> <strong>masas</strong> <strong>de</strong> una distribución al punto respecto al cual<br />
el momento estático es nulo. Llamaremos C a este punto, cuya existencia<br />
y unicidad está garantizada , pues, por la ecuación (1), se obtiene:<br />
M C = 0 = M O − MOC ⇐⇒ OC = M O<br />
M<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 6 <strong>de</strong> 41
lo que permite posicionar el centro <strong>de</strong> <strong>masas</strong> C <strong>de</strong> una distribución conocido<br />
el momento estático respecto a un punto cualquiera y la masa.<br />
Ya se ha citado que una distribución <strong>de</strong> <strong>masas</strong> <strong>de</strong>fine un campo <strong>de</strong><br />
momentos estáticos. Pue<strong>de</strong> utilizarse la ecuación (1) para visualizar la<br />
forma <strong>de</strong> este campo. En efecto, M O = MOC por lo que se trata <strong>de</strong><br />
un campo central emergente con simetría radial, lineal con la distancia al<br />
centro <strong>de</strong> <strong>masas</strong>.<br />
Al igual que se ha <strong>de</strong>finido el momento estático respecto a un punto,<br />
se pue<strong>de</strong> <strong>de</strong>finir el momento estático respecto a una recta. La posición<br />
relativa entre dos puntos viene dada por el vector que los une, sin embargo<br />
la posición relativa <strong>de</strong> un punto P respecto a una recta viene <strong>de</strong>terminada<br />
por el vector perpendicular a esta última que tiene como origen uno <strong>de</strong> sus<br />
puntos y como extremo el punto P . La expresión analítica <strong>de</strong> este vector<br />
es, según se explica en [2], u × (OP × u), siendo O un punto <strong>de</strong> la recta y<br />
u un vector unitario <strong>de</strong> su misma dirección. El momento estático respecto<br />
a una recta δ se <strong>de</strong>fine,<br />
∫<br />
M δ = u × (OP × u)ρdV<br />
Para que esta <strong>de</strong>finición sea correcta, el momento estático no <strong>de</strong>be <strong>de</strong>pen<strong>de</strong>r<br />
ni <strong>de</strong>l sentido elegido para el vector u ni <strong>de</strong> la elección <strong>de</strong>l punto<br />
O sobre la recta. Evi<strong>de</strong>ntemente, la elección <strong>de</strong>l sentido <strong>de</strong> u, al aparecer<br />
éste multiplicando dos veces el integrando, no influye en el resultado <strong>de</strong><br />
la integral. En cuanto a la elección <strong>de</strong>l punto O, supondremos que se ha<br />
elegido un punto O ′ y veremos que el resultado obtenido es el mismo. Dado<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 7 <strong>de</strong> 41
que O, O ′ ∈ δ ⇒ O ′ O = λu ⇒ O ′ P = O ′ O + OP Por lo que<br />
∫<br />
M ′ δ = u × ((O ′ O + OP ) × u)dV =<br />
∫<br />
∫<br />
u × ((λu + OP ) × u)dV = u × (OP × u)dV = M δ<br />
Pue<strong>de</strong> relacionarse el momento estático respecto a una recta con el momento<br />
estático respecto a uno cualquiera <strong>de</strong> sus puntos.<br />
∫<br />
∫<br />
M δ = u × (OP × u)dV = u × [( OP dV ) × u] = u × [M O × u]<br />
Es <strong>de</strong>cir, el momento estático respecto a una recta es la proyección sobre<br />
el plano perpendicular a dicha recta <strong>de</strong>l momento estático respecto<br />
a cualquiera <strong>de</strong> sus puntos. A<strong>de</strong>más, ya que OC = M O<br />
, se concluye<br />
M<br />
que el centro <strong>de</strong> <strong>masas</strong> se encuentra sobre una recta paralela a δ <strong>de</strong>splazada<br />
respecto a ésta el vector M δ<br />
. Si la recta δ contiene a C, entonces,<br />
M<br />
obviamente, M δ = 0.<br />
De la misma forma que se ha <strong>de</strong>finido el momento estático para rectas,<br />
se pue<strong>de</strong> <strong>de</strong>finir para planos. Sea π un plano, u un vector unitario perpendicular<br />
y O un punto <strong>de</strong> π. Se <strong>de</strong>fine el momento estático respecto al<br />
plano π como:<br />
∫<br />
M π = (u · OP )uρdV<br />
Esta <strong>de</strong>finición no <strong>de</strong>be <strong>de</strong>pen<strong>de</strong>r <strong>de</strong> la elección <strong>de</strong>l sentido <strong>de</strong> u ni <strong>de</strong> la<br />
<strong>de</strong>l punto O. En efecto, la aparición <strong>de</strong>l vector multiplicando dos veces en<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 8 <strong>de</strong> 41
el integrando hace que éste no <strong>de</strong>penda <strong>de</strong>l sentido <strong>de</strong> u. En cuanto a la<br />
elección <strong>de</strong> O, sea O ′ ∈ π:<br />
∫<br />
∫<br />
M ′ π = (u · O ′ P )uρdV = (u · [O ′ O + OP ])uρdV<br />
∫<br />
M ′ π = (u · OP )uρdV = M π<br />
ya que u · O ′ O = 0 pues O, O ′ pertenecen a π.<br />
Se comprueba que<br />
∫<br />
∫<br />
M π = (u · OP )uρdV = (u · [<br />
OP ρdV ])u<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
es <strong>de</strong>cir, que el momento estático respecto a un plano es la proyección<br />
perpendicular al mismo <strong>de</strong>l momento estático respecto a cualquiera <strong>de</strong> sus<br />
puntos. Igualmente, el momento estático respecto a un plano proporciona<br />
información sobre la posición <strong>de</strong>l centro <strong>de</strong> <strong>masas</strong> respecto a dicho plano.<br />
En efecto, si el momento estático respecto a π <strong>de</strong> una distribución <strong>de</strong><br />
<strong>masas</strong> es M π , entonces C se encuentra en un plano paralelo al anterior,<br />
<strong>de</strong>splazado M π<br />
M y, si C ∈ π ⇒ M π = 0.<br />
Dado que si i, j, k son tres vectores unitarios perpendiculares, se tiene:<br />
se tiene que<br />
OP = (OP · i)i + (OP · j)j + (OP · k)k<br />
Título<br />
Contenido<br />
◭◭ ◮◮<br />
◭ ◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 9 <strong>de</strong> 41<br />
M O = (M yz ·i)i+(M zx·j)j +(M xy ·k)k ⇒ M O = M yz +M zx +M xy
es <strong>de</strong>cir, el momento estático respecto a un punto es la suma <strong>de</strong> los momentos<br />
estáticos respecto a tres planos mutuamente perpendiculares que<br />
pasen por dicho punto.<br />
Por otra parte, dado que<br />
OP = 1 [i × (OP × i) + j × (OP × j) + k × (OP × k)]<br />
2<br />
se tiene que<br />
M O = 1 2 [M x + M y + M z ]<br />
es <strong>de</strong>cir, el momento estático respecto a un punto es la semisuma <strong>de</strong> los<br />
momentos estáticos respecto a tres rectas mutuamente perpendiculares que<br />
pasen por dicho punto.<br />
Se pue<strong>de</strong> completar la relación anterior aña diendo la relación entre los<br />
momentos estáticos respecto a rectas y planos. En efecto, dado que<br />
se tiene que<br />
i × (OP × i) = (j · OP )j + (k · OP )k<br />
M x = M zx + M xy<br />
es <strong>de</strong>cir, el momento estático respecto a una recta es la suma <strong>de</strong> los momentos<br />
estáticos respecto a dos planos mutuamente perpendiculares que la<br />
contengan.<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 10 <strong>de</strong> 41
4. Momentos <strong>de</strong> inercia<br />
Los momentos estáticos tratados en la sección anterior se conocen también<br />
como momentos <strong>de</strong> primer or<strong>de</strong>n, ya que las distancias <strong>de</strong> los puntos<br />
másicos a los elementos respecto a los que están <strong>de</strong>finidos intervienen con<br />
exponente uno (el momento <strong>de</strong> or<strong>de</strong>n cero sería la masa total <strong>de</strong> la distribución).<br />
Como hemos dicho en la introducción, nos interesan los momentos<br />
<strong>de</strong> hasta el segundo or<strong>de</strong>n. Estos se <strong>de</strong>nominan momentos <strong>de</strong> inercia e<br />
intervienen en casi todas las ecuaciones <strong>de</strong> la dinámica <strong>de</strong>l sólido rígido.<br />
Al igual que los momentos estáticos, se <strong>de</strong>finen respecto a puntos, rectas<br />
y planos y, a diferencia <strong>de</strong> aquéllos, también se <strong>de</strong>finen respecto a pares<br />
<strong>de</strong> planos. Comencemos por los momentos <strong>de</strong> inercia respecto a puntos o<br />
momentos <strong>de</strong> inercia centrales.<br />
Dado un punto O, se <strong>de</strong>fine el momento <strong>de</strong> inercia respecto a O, y se<br />
<strong>de</strong>notará I O , a:<br />
∫<br />
I O = nor(OP )ρdV<br />
El momento <strong>de</strong> inercia será siempre, por su <strong>de</strong>finición, mayor que cero<br />
(excepto en el caso <strong>de</strong> una masa puntual respecto al punto en que esté<br />
concentrada, en cuyo caso sería nulo. Pue<strong>de</strong> relacionarse el momento <strong>de</strong><br />
inercia respecto a un punto O cualquiera con el momento <strong>de</strong> inercia respecto<br />
al centro <strong>de</strong> <strong>masas</strong> C.<br />
∫<br />
∫<br />
I O = nor(OC+CP )ρdV = [nor(OC)+2OC·CP +nor(CP )]ρdV =<br />
∫<br />
= nor(OC) ρdV + I C + 2OC · M C<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 11 <strong>de</strong> 41
Pero como M C = 0, se tiene<br />
I O = M nor(OC) + I C<br />
lo que indica que el momento <strong>de</strong> inercia central es mínimo en el centro <strong>de</strong><br />
<strong>masas</strong> y se incrementa respecto a éste en el producto <strong>de</strong> la masa por el<br />
cuadrado <strong>de</strong> la distancia entre C y O.<br />
Los momentos <strong>de</strong> inercia más útiles son los que se <strong>de</strong>finen respecto a<br />
rectas. Sean una recta δ, un vector unitario paralelo a ella u y uno <strong>de</strong> sus<br />
puntos O. Se <strong>de</strong>fine el momento <strong>de</strong> inercia I δ respecto a δ como<br />
∫<br />
I δ = nor(OP × u)ρdV<br />
La <strong>de</strong>finición anterior es correcta, pues no <strong>de</strong>pen<strong>de</strong> ni <strong>de</strong>l sentido asignado a<br />
u (neutralizado por la norma), ni <strong>de</strong> la elección <strong>de</strong> O, ya que cualquier otro<br />
punto <strong>de</strong>finiría un vector O ′ P = λu+OP que multiplicado vectorialmente<br />
por u daría el mismo resultado.<br />
El momento <strong>de</strong> inercia será siempre positivo, excepto si la distribución<br />
<strong>de</strong> <strong>masas</strong> está contenida en la recta δ, en cuyo caso es nulo.<br />
Se pue<strong>de</strong> relacionar el momento <strong>de</strong> inercia <strong>de</strong> una distribución cualquiera<br />
<strong>de</strong> <strong>masas</strong> respecto a una recta δ con el correspondiente respecto a una recta<br />
δ C paralela a δ que contenga al centro <strong>de</strong> <strong>masas</strong>.<br />
∫<br />
I δ = nor[(OC + CP ) × u]ρdV =<br />
∫<br />
= [nor(OC × u) + 2(OC × u) · (CP × u) + nor(CP × u)]ρdV =<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 12 <strong>de</strong> 41
∫<br />
= nor(OC × u)<br />
ρdV + I δC + 2(OC × u) · (M C × u)<br />
Pero como M C = 0, se tiene<br />
I δ = M nor(OC × u) + I δC<br />
El primer sumando <strong>de</strong>l segundo miembro <strong>de</strong> la expresión anterior es el<br />
producto <strong>de</strong> la masa por el cuadrado <strong>de</strong> la distancia <strong>de</strong>l centro <strong>de</strong> <strong>masas</strong> a<br />
δ (o, si se prefiere, la distancia entre δ y δ C ). Esta ecuación se conoce como<br />
fórmula <strong>de</strong> Steiner. De ella se <strong>de</strong>duce que <strong>de</strong> todas las rectas paralelas a<br />
una dada, la que menor momento <strong>de</strong> inercia presenta es la que pasa por el<br />
centro <strong>de</strong> <strong>masas</strong>.<br />
Se suele <strong>de</strong>finir el radio <strong>de</strong> giro k δ <strong>de</strong> una distribución <strong>de</strong> <strong>masas</strong> respecto<br />
a la recta δ como el radio <strong>de</strong> un cilindro <strong>de</strong> revolución <strong>de</strong> eje δ, homogéneo y<br />
<strong>de</strong> masa igual a la <strong>de</strong> la distribución original que <strong>de</strong>fina el mismo momento<br />
<strong>de</strong> inercia respecto a δ. Por lo tanto,<br />
√<br />
Iδ<br />
k δ =<br />
M<br />
Sean un plano π , un vector unitario perpendicular u y un punto O<br />
contenido en π. Se <strong>de</strong>fine el momento <strong>de</strong> inercia respecto a π como:<br />
∫<br />
I π = (OP · u) 2 ρdV<br />
que es, obviamente, in<strong>de</strong>pendiente <strong>de</strong>l sentido <strong>de</strong> u y <strong>de</strong> la elección <strong>de</strong> O, ya<br />
que otro punto O ′ haría O ′ P = O ′ O + OP y O ′ O · OP = 0. El momento<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 13 <strong>de</strong> 41
I π será positivo excepto si la distribución <strong>de</strong> <strong>masas</strong> está contenida en el<br />
plano π en cuyo caso sería nulo.<br />
Se pue<strong>de</strong> relacionar el momento <strong>de</strong> inercia <strong>de</strong> una distribución cualquiera<br />
<strong>de</strong> <strong>masas</strong> respecto a un plano π con el correspondiente respecto a un plano<br />
π C paralelo a π que contenga al centro <strong>de</strong> <strong>masas</strong>.<br />
∫<br />
∫<br />
I π = [(OC+CP )·u] 2 ρdV = [(OC·u) 2 +2(OC·u)(CP ·u)+(CP ·u) 2 ]ρdV =<br />
= (OC · u) 2 ∫<br />
Pero como M C = 0, se tiene<br />
ρdV + I πC + 2(OC · u)(M C · u)<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
I π = M(OC · u) 2 + I πC<br />
El primer sumando <strong>de</strong>l segundo miembro <strong>de</strong> la expresión anterior es el<br />
producto <strong>de</strong> la masa por el cuadrado <strong>de</strong> la distancia <strong>de</strong>l centro <strong>de</strong> <strong>masas</strong> a<br />
π (o, si se prefiere, la distancia entre π y π C ). De esta ecuación se <strong>de</strong>duce<br />
que <strong>de</strong> todos los planos paralelos a uno dado, el que menor momento <strong>de</strong><br />
inercia presenta es el que pasa por el centro <strong>de</strong> <strong>masas</strong>.<br />
Finalmente, se <strong>de</strong>fine el producto <strong>de</strong> inercia respecto al par <strong>de</strong> planos<br />
orientados π, π ′ . Sean u, u ′ dos vectores unitarios normales a π y π ′ congruentes<br />
con sus respectivas orientaciones. Sean asímismo O, O ′ dos puntos<br />
<strong>de</strong> π y π ′ respectivamente. Se <strong>de</strong>fine el producto <strong>de</strong> inercia<br />
∫<br />
P ππ ′ = (OP · u)(O ′ P · u ′ )ρdV<br />
Título<br />
Contenido<br />
◭◭ ◮◮<br />
◭ ◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 14 <strong>de</strong> 41
que no <strong>de</strong>pen<strong>de</strong> <strong>de</strong> la elección <strong>de</strong> O, O ′ ya que OP , O ′ P están multiplicados<br />
escalarmente por U, u ′ respectivamente. Evi<strong>de</strong>ntemente, el producto <strong>de</strong><br />
inercia es simétrico respecto a los planos que lo <strong>de</strong>finen. A diferencia <strong>de</strong><br />
los momentos <strong>de</strong> inercia respecto a puntos, rectas o planos, el producto <strong>de</strong><br />
inercia pue<strong>de</strong> ser positivo, negativo o nulo.<br />
Se pue<strong>de</strong> relacionar el producto <strong>de</strong> inercia <strong>de</strong> una distribución cualquiera<br />
<strong>de</strong> <strong>masas</strong> respecto a un par <strong>de</strong> planos π, π ′ con el correspondiente a un par<br />
<strong>de</strong> planoa paralelos π C , π C ′ que pasen por el centro <strong>de</strong> <strong>masas</strong>.<br />
∫<br />
P π,π ′ = (OP · u)(O ′ P · u ′ )ρdV =<br />
∫<br />
= (OC + CP · u)(O ′ C + CP · u ′ )ρdV =<br />
∫<br />
[(OC · u)(O ′ C · u ′ ) + (OC · u)(CP · u ′ )+<br />
+(CP · u)(O ′ P · u ′ ) + (CP · u)(CP · u ′ )]ρdV<br />
Pero como M C = 0, se tiene<br />
P ππ ′ = M(OC · u)(O ′ C · u ′ ) + P πC π ′ C<br />
El primer sumando <strong>de</strong>l segundo miembro <strong>de</strong> la expresión anterior es el<br />
producto <strong>de</strong> la masa por el producto <strong>de</strong> las distancias con signo <strong>de</strong>s<strong>de</strong> los<br />
planos orientados al centro <strong>de</strong> <strong>masas</strong> (o, si se prefiere, el producto <strong>de</strong> las<br />
distancias <strong>de</strong> π a π C y <strong>de</strong> π ′ a π ′ C ).<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 15 <strong>de</strong> 41
Si i, j, k son tres vectores unitarios perpendiculares, se tiene:<br />
nor OP = (OP · i) 2 + (OP · j) 2 + (OP · k) 2<br />
nor OP = 1 [nor(OP × i) + nor(OP × j) + nor(OP × k)]<br />
2<br />
se tiene que<br />
y que<br />
I O = I yz + I zx + I xy = I x + I y + I z<br />
2<br />
I x = I zx + I xy<br />
2<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
Título<br />
Contenido<br />
◭◭<br />
◭<br />
◮◮<br />
◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 16 <strong>de</strong> 41
5. Tensor y elipsoi<strong>de</strong> <strong>de</strong> inercia<br />
Si el vector unitario u y el vector OP se expresan por sus componentes en<br />
una base ortonormal u = cos αi + cos βj + cos γk, OP = xi + yj + zk, se<br />
tiene:<br />
nor(OP × u) = nor OP nor u − (OP · u) 2 =<br />
= (x 2 + y 2 + z 2 ) − (x cos α + y cos β + z cos γ) 2 =<br />
x 2 (1 − cos 2 α) + y 2 (1 − cos 2 β) + z 2 (1 − cos 2 γ)−<br />
−2yz cos β cos γ − 2zx cos γ cos α − 2xy cos α cos β =<br />
= x 2 (cos 2 β + cos 2 γ) + y 2 (cos 2 γ + cos 2 α) + z 2 (cos 2 α + cos 2 β)−<br />
−2yz cos β cos γ − 2zx cos γ cos α − 2xy cos α cos β =<br />
= cos 2 α(y 2 + z 2 ) + cos 2 β(z 2 + x 2 ) + cos 2 γ(x 2 + y 2 )−<br />
−2yz cos β cos γ − 2zx cos γ cos α − 2xy cos α cos β<br />
Por lo que, si δ es una recta paralela a u que pasa por O, se tiene<br />
∫<br />
I δ = [cos 2 α(y 2 + z 2 ) + cos 2 β(z 2 + x 2 ) + cos 2 γ(x 2 + y 2 )−<br />
−2yz cos β cos γ − 2zx cos γ cos α − 2xy cos α cos β]ρdV =<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
◮◮<br />
◮<br />
Salir<br />
I δ = cos 2 αI x + cos 2 βI y + cos 2 γI z − 2 cos β cos γP yz − 2 cos γ cos αP zx − 2 cos α cos βP xy<br />
Página 17 <strong>de</strong> 41<br />
(2)
La expresión (2) representa una función cuadrática en cos α, cos β, cos γ<br />
que se pue<strong>de</strong> representar <strong>de</strong> la forma:<br />
⎛<br />
⎞ ⎛ ⎞<br />
I δ = ( cos α cos β cos γ ) I x −P xy −P zx cos α<br />
⎝ −P xy I y −P yz<br />
⎠ ⎝ cos β ⎠ (3)<br />
−P zx −P yz I z cos γ<br />
⎞<br />
I δ es un escalar, ( cos α cos β cos γ ) son las componentes <strong>de</strong>l vector u<br />
y la matriz<br />
⎛<br />
I x −P xy −P zx<br />
−P zx −P yz I z<br />
⎝ −P xy I y −P yz<br />
⎠<br />
representa las componentes en forma matricial simétrica <strong>de</strong> la función<br />
cuadrática I δ (u). Por lo tanto 1 , la expresión<br />
⎛<br />
⎞ ⎛ ⎞<br />
I x −P xy −P zx cos α<br />
l(u) = ⎝ −P xy I y −P yz<br />
⎠ ⎝ cos β ⎠<br />
−P zx −P yz I z cos γ<br />
representa una función vectorial <strong>de</strong> variable vectorial. Por consiguiente, las<br />
matriz M que representa esta función cambia sus componentes al cambiar<br />
<strong>de</strong> base <strong>de</strong> referencia según las fórmulas <strong>de</strong> cambio <strong>de</strong> base ortonormal:<br />
M ′ = T MT ′<br />
1 según los resultados <strong>de</strong>l Algebra Lineal, una forma cuadrática <strong>de</strong>termina un tensor<br />
o diádica simétrica, la cual, a su vez, <strong>de</strong>fine una función vectorial lineal<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 18 <strong>de</strong> 41
y representa las componentes <strong>de</strong>l tensor o diádica <strong>de</strong> inercia.<br />
El concepto <strong>de</strong> tensor traduce la i<strong>de</strong>a <strong>de</strong> una función vectorial <strong>de</strong> variable<br />
vectorial con existencia intrínseca, es <strong>de</strong>cir, in<strong>de</strong>pendiente <strong>de</strong> la base<br />
o forma elegida para su representación. La función l(u) representa una<br />
auténtica función vectorial, que no pue<strong>de</strong> representarse meramente por una<br />
matriz, pues frente a un cambio <strong>de</strong> base, <strong>de</strong>be cambiar <strong>de</strong> componentes.<br />
En nuestro <strong>de</strong>sarrollo hemos partido <strong>de</strong> la radiación <strong>de</strong> rectas que pasaba<br />
por un punto O. Queda, pues, <strong>de</strong>finido un tensor <strong>de</strong> inercia en cada punto<br />
<strong>de</strong>l espacio. Utilizaremos una doble barra sobre una letra mayúscula para<br />
<strong>de</strong>notar un tensor. El tensor <strong>de</strong> inercia <strong>de</strong>l punto O se escribirá I O y en la<br />
base i, j, k queda <strong>de</strong>finido por la matriz<br />
⎛<br />
⎞<br />
I x −P xy −P zx<br />
(I 0 ) x,y,z = ⎝ −P xy I y −P yz<br />
⎠<br />
−P zx −P yz I z<br />
Frente a un cambio <strong>de</strong> base a la i ′ , j ′ , k ′ se tiene:<br />
⎛<br />
I x ′ −P x ′ y ′ −P ⎞ ⎛<br />
⎞<br />
z ′ x ′<br />
I x −P xy −P zx<br />
⎝ −P x ′ y ′ I y ′ −P ⎠<br />
y ′ z ′ = T ⎝ −P xy I y −P yz<br />
⎠ T ′<br />
−P z ′ x ′ −P y ′ z ′ I z ′ −P zx −P yz I z<br />
y el momento <strong>de</strong> inercia respecto a δ se pue<strong>de</strong> escribir I δ = u · I O u = u · l<br />
En algunas ocasiones es útil disponer <strong>de</strong> una representación geométrica<br />
<strong>de</strong> las características <strong>de</strong> inercia <strong>de</strong> un sólido en un punto O. Para ello se<br />
<strong>de</strong>fine el elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> O <strong>de</strong> la siguiente forma:<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 19 <strong>de</strong> 41
Sobre cada una <strong>de</strong> las rectas que pasan por O, se lleva a ambos lados la<br />
distancia √ k<br />
Iδ<br />
, don<strong>de</strong> k es una constante arbitraria positiva con dimensiones<br />
<strong>de</strong> [m] 1 2 [L] 2 . Se tiene entonces la superficie:<br />
r = ± k √<br />
Iδ<br />
u ⇒ I δ nor r = k 2 =<br />
= (|r|u) · (I) O (|r|u) = k 2<br />
= r · (I) O r = k 2<br />
I x x 2 + I y y 2 + I z z 2 − 2P yz yz − 2P zx zx − 2P xy xy = k 2<br />
ecuación <strong>de</strong> un elipsoi<strong>de</strong> referido a su centro. El tensor <strong>de</strong> inercia viene<br />
representado por una matriz simétrica que, dado que I δ ≥ 0 <strong>de</strong>be tener tres<br />
valores propios reales no negativos y , al menos, tres direcciones invariantes<br />
asociadas mutuamente perpendiculares. Si se eligen estas direcciones como<br />
base para representar el tensor, la matriz será diagonal con los momentos<br />
<strong>de</strong> inercia (que serán los valores propios) respectivos en dicha diagonal. La<br />
representación <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia tensor en esta referencia sería<br />
I x x 2 + I y y 2 + I z z 2 = k 2<br />
en la que se aprecia que está referido a sus ejes.<br />
Cuando el elipsoi<strong>de</strong> <strong>de</strong> inercia es <strong>de</strong> revolución, entonces existe un eje<br />
<strong>de</strong> revolución y todos los perpendiculares a este eje que pasan por el centro<br />
son asímismo ejes <strong>de</strong>l elipsoi<strong>de</strong>, todos ellos tienen el mismo momento <strong>de</strong><br />
inercia y <strong>de</strong>terminan lo que se conoce como plano ecuatorial. En este caso el<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 20 <strong>de</strong> 41
tensor <strong>de</strong> inercia tiene un valor propio doble y otro simple, correspondiendo<br />
el autovector <strong>de</strong> este último a la dirección <strong>de</strong>l eje <strong>de</strong> revolución, mientras<br />
que el subespacio propio asociado al autovalor doble es el plano ecuatorial.<br />
Si el elipsoi<strong>de</strong> es esférico, entonces todas las rectas que pasan por el centro<br />
son ejes <strong>de</strong>l mismo, comparten el mismo momento <strong>de</strong> inercia y el tensor <strong>de</strong><br />
inercia se pue<strong>de</strong> escribir, en cualquier base, como el producto <strong>de</strong>l momento<br />
<strong>de</strong> inercia común por el tensor i<strong>de</strong>ntidad U.<br />
La ecuación recién obtenida evi<strong>de</strong>ncia que el mayor y el menor momento<br />
<strong>de</strong> inercia <strong>de</strong>l sólido se da, respectivamente sobre el menor y mayor eje <strong>de</strong>l<br />
elipsoi<strong>de</strong>.<br />
Las fórmulas <strong>de</strong> Steiner afectan a las componentes <strong>de</strong>l tensor <strong>de</strong> inercia.Escribimos<br />
su expresión matricial :<br />
(I) = (I ∗ ) + M nor(r c )(U) − M(r c ) T (r c ) (4)<br />
don<strong>de</strong> U es la matriz i<strong>de</strong>ntidad 3 × 3 y I ∗ es el tensor <strong>de</strong> inercia <strong>de</strong>l centro<br />
<strong>de</strong> <strong>masas</strong>. Esta fórmula resume todas las <strong>de</strong> Steiner.<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
Título<br />
Contenido<br />
◮◮<br />
◭<br />
◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 21 <strong>de</strong> 41
6. Simetrías en la distribución<br />
Cuando la distribución <strong>de</strong> <strong>masas</strong> <strong>de</strong> partida presenta algunas características<br />
especiales <strong>de</strong> simetría, entonces los momentos estáticos y los <strong>de</strong> inercia<br />
presentan algunas relaciones pre<strong>de</strong>cibles, que se explicitan a continuación:<br />
• si π es un plano <strong>de</strong> simetría <strong>de</strong>l sólido, entonces:<br />
– C ∈ π es <strong>de</strong>cir, el momentyo estático ⃗ M π = ⃗0<br />
– el producto <strong>de</strong> inercia <strong>de</strong> π con cualquier otro plano π ′ perpendicular<br />
a π es nulo P ππ ′ = 0.<br />
– todas las rectas perpendiculares a π son ejes <strong>de</strong> los elipsoi<strong>de</strong>s <strong>de</strong><br />
inercia <strong>de</strong> los puntos en que cortan a π.<br />
• si δ es un eje <strong>de</strong> simetría <strong>de</strong>l sólido, entonces:<br />
– C ∈ δ es <strong>de</strong>cir, el momentyo estático ⃗ M δ = ⃗0<br />
– el producto <strong>de</strong> inercia <strong>de</strong> cualquier plano perpendicular a δ π con<br />
cualquier otro plano π ′ que contenga a δ es nulo P ππ ′ = 0.<br />
– δ es eje <strong>de</strong> los elipsoi<strong>de</strong>s <strong>de</strong> inercia <strong>de</strong> todos sus puntos.<br />
• si δ es un eje <strong>de</strong> revolución <strong>de</strong>l sólido, entonces:<br />
– C ∈ δ es <strong>de</strong>cir, el momentyo estático ⃗ M δ = ⃗0<br />
– el producto <strong>de</strong> inercia <strong>de</strong> cualquier plano perpendicular a δ π con<br />
cualquier otro plano π ′ que contenga a δ es nulo P ππ ′ = 0.<br />
– δ es eje <strong>de</strong> los elipsoi<strong>de</strong>s <strong>de</strong> inercia <strong>de</strong> todos sus puntos.<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 22 <strong>de</strong> 41
– todos los elipsoi<strong>de</strong>s <strong>de</strong> inercia <strong>de</strong> los puntos <strong>de</strong> δ son <strong>de</strong> revolución<br />
en torno a δ.<br />
– las rectas <strong>de</strong>l plano ecuatorial <strong>de</strong> dichos elipsoi<strong>de</strong>s son ejes <strong>de</strong> los<br />
mismos y tienen todas el mismo momento <strong>de</strong> inercia.<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
Título<br />
Contenido<br />
◭◭<br />
◭<br />
◮◮<br />
◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 23 <strong>de</strong> 41
7. Ejes <strong>de</strong> los elipsoi<strong>de</strong>s <strong>de</strong> inercia<br />
Las características <strong>de</strong> inercia <strong>de</strong> un sólido rígido respecto a un punto <strong>de</strong>l<br />
espacio están con<strong>de</strong>nsadas en su elipsoi<strong>de</strong> <strong>de</strong> inercia.<br />
La ecuación <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> un punto O, si I es su tensor <strong>de</strong><br />
inercia, es<br />
r · I · r = k 2<br />
don<strong>de</strong> r es un vector <strong>de</strong> posición con origen en O y k es una constante<br />
arbitraria con dimensiones.<br />
Dado un punto O <strong>de</strong>l espacio, po<strong>de</strong>mos clasificar el conjunto <strong>de</strong> rectas<br />
que pasan por O en dos grupos: las que son ejes <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia<br />
<strong>de</strong> O, y las que no lo son. La condición para que una recta δ, paralela<br />
al vector unitario u, sea eje <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia anterior, es que corte<br />
perpendicularmente a dicho elipsoi<strong>de</strong>.<br />
La normal al elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> O en un punto P posicionado por<br />
el vector r es paralela al vector<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
Título<br />
Contenido<br />
◮◮<br />
grad(r · I · r) = 2I · r<br />
◭<br />
◮<br />
pero, dado que r es paralelo a u, un vector normal al elipsoi<strong>de</strong> <strong>de</strong> inercia<br />
es:<br />
I · u<br />
y la condición para que δ sea eje <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> O es que<br />
I o · u = λu<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 24 <strong>de</strong> 41
para algún λ real, lo cual expresa la ortogonalidad entre recta y elipsoi<strong>de</strong>.<br />
El primer miembro <strong>de</strong> la ecuación anterior representa una función vectorial<br />
<strong>de</strong>finida para cada punto O <strong>de</strong>l espacio y para cada recta que pase<br />
por ese punto. Llamemos a dicha función vector <strong>de</strong> inercia l(O, u), ya<br />
documentada en [3]:<br />
l(O, u) = I o · u<br />
y la condición anterior para que δ sea eje <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> O es:<br />
l(O, u) × u = 0<br />
Escribimos la expresión matricial <strong>de</strong>l teorema <strong>de</strong> Steiner:<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
(I) = (I ∗ ) + M nor(r c )(U) − M(r c ) T (r c )<br />
don<strong>de</strong> U es la matriz i<strong>de</strong>ntidad 3 × 3 y I ∗ es el tensor <strong>de</strong> inercia <strong>de</strong>l centro<br />
<strong>de</strong> <strong>masas</strong>. La función l(O, u pue<strong>de</strong> escribirse, a la vista <strong>de</strong> la ecuación<br />
anterior:<br />
l(O, u) = l(C, u) + M(nor r c )u − Mr c (r c · u)<br />
lo que constituye la fórmula <strong>de</strong> Steiner para el vector <strong>de</strong> inercia.<br />
La condición para que δ sea eje <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> O queda<br />
l(C, u) × u = M(r c · u)r c × u (5)<br />
Título<br />
Contenido<br />
◭◭ ◮◮<br />
◭ ◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 25 <strong>de</strong> 41
8. Clasificación <strong>de</strong> las rectas<br />
En la sección anterior se ha estudiado una propiedad <strong>de</strong>finida entre una recta<br />
y el elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> un punto. Proce<strong>de</strong>remos a estudiar el número<br />
<strong>de</strong> puntos <strong>de</strong> cuyos elipsoi<strong>de</strong>s pue<strong>de</strong> ser eje una recta dada. El resultado<br />
<strong>de</strong> este análisis permitirá realizar una clasificación <strong>de</strong>l conjunto <strong>de</strong> rectas<br />
<strong>de</strong>l espacio en función <strong>de</strong> unas características <strong>de</strong> inercia cualitativamente<br />
distintas.<br />
Realizaremos un análisis por casos, según el centro <strong>de</strong> <strong>masas</strong> pertenezca<br />
a la recta δ o sea exterior a la misma.<br />
8.1. La recta pasa por el centro <strong>de</strong> <strong>masas</strong><br />
En este caso, r c × u = 0. El vector <strong>de</strong> inercia es:<br />
y la condición (5) queda:<br />
l(O, u) = l(C, u) + M nor r c u − M(r c · u)r } {{ } c<br />
=0<br />
l(O, u) × u = l(C, u) × u + M nor r c u<br />
} {{<br />
× u<br />
}<br />
⇒<br />
=0<br />
l(O, u) × u = l(C, u) × u<br />
es <strong>de</strong>cir, la recta δ será eje <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> O si y sólo si la recta<br />
δ es eje <strong>de</strong>l elipsoi<strong>de</strong> central <strong>de</strong> inercia. En este caso, dado que O es un<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 26 <strong>de</strong> 41
punto cualquiera <strong>de</strong> la recta δ, ésta será eje <strong>de</strong> los elipsoi<strong>de</strong>s <strong>de</strong> inercia <strong>de</strong><br />
todos sus puntos. Este resultado se pue<strong>de</strong> enunciar solemnemente en el<br />
siguiente:<br />
Teorema 1 Entre todas las rectas que pasan por el centro <strong>de</strong> <strong>masas</strong> se<br />
distinguen dos únicas clases: las que son ejes <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong><br />
todos sus puntos y las que no lo son <strong>de</strong> ninguno <strong>de</strong> ellos.<br />
8.2. La recta no pasa por el centro <strong>de</strong> <strong>masas</strong><br />
En este caso, r c × u ≠ 0, el vector <strong>de</strong> inercia en un punto O <strong>de</strong> la recta δ<br />
es<br />
l(O, u) = l(C, u) + M nor r c u − M(r c · u)r } {{ }<br />
′ c (6)<br />
=0<br />
buscamos un punto O ′ <strong>de</strong> dicha recta cuyo elipsoi<strong>de</strong> <strong>de</strong> inercia la tenga por<br />
eje. Es <strong>de</strong>cir, un punto O ′ tal que:<br />
Su vector <strong>de</strong> inercia es:<br />
l(O ′ , u)u = 0<br />
l(O ′ , u) = l(C, u) + M nor r ′ cu − M(r ′ c · u)r ′ c (7)<br />
como O, P ∈ δ entonces el vector O ′ O es paralelo a u y se pue<strong>de</strong> escribir<br />
O ′ O = λu don<strong>de</strong> λ es un parámetro real que posiciona P sobre la recta,<br />
una vez elegido el punto O, que actúa como origen <strong>de</strong> δ. De esta forma:<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 27 <strong>de</strong> 41<br />
r ′ c = O ′ C = O ′ O + OC = O ′ O + r c = λu + r c
La ecuación (7)se escribe<br />
l(O ′ , u) = l(C, u) + M nor(λu + r c )u − M((λu + r c ) · u)(λu + r c )<br />
multiplicando vectorialmente por u<br />
l(O ′ , u)×u = l(C, u)×u+M nor(λu+r c ) u<br />
} {{<br />
× u<br />
}<br />
−M((λu+r c )·u)(λu+r c )×u<br />
=0<br />
l(O ′ , u)×u = l(C, u)×u−M((λu+r c )·u)λ u<br />
} {{<br />
× u<br />
}<br />
−M((λu+r c )·u)r c ×u<br />
=0<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
l(O ′ , u) × u = l(C, u) × u − λMr c × u − M(r c · u)r c × u (8)<br />
multiplicando la ecuación (6) vectorialmente por u<br />
l(O, u) × u = l(C, u) × u − M(r c · u)r c × u (9)<br />
Haciendo l(O ′ , u) × u = 0 y restando (8) <strong>de</strong> (9) queda:<br />
Título<br />
Contenido<br />
◭◭ ◮◮<br />
l(O, u) × u = λMr c × u (10)<br />
◭<br />
◮<br />
Si existe algún valor <strong>de</strong> λ que satisface la ecuación (10) entonces existe<br />
un punto O ′ posicionado por OO ′ = −λu cuyo elipsoi<strong>de</strong> <strong>de</strong> inercia tiene a<br />
la recta δ como eje. Pero la ecuación anterior sólo pue<strong>de</strong> satisfacerse si los<br />
vectores l(O, u)×u y Mr c ×u son paralelos pues, en caso contrario, ningún<br />
escalar λ, multiplicado por el segundo, pue<strong>de</strong> dar el primero. La colinealidad<br />
<strong>de</strong> estos dos vectores exige el paralelismo <strong>de</strong> los planos (l(O, u), u)<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 28 <strong>de</strong> 41
y (r c , u) lo que matemáticamente se expresa por la nulidad <strong>de</strong>l producto<br />
mixto:<br />
(l(O, u), r c , u) = 0 (11)<br />
que pue<strong>de</strong> escribirse, aplicando la fórmula <strong>de</strong> Steiner para el vector <strong>de</strong><br />
inercia:<br />
(l(C, u) + M nor r c u − M(u · r c )r c , r c , u) = 0<br />
don<strong>de</strong> los dos últimos sumandos <strong>de</strong>l primer factor son paralelos al tercer<br />
y segundo factores respectivamente, por lo que pue<strong>de</strong>n simplificarse, obteniendo<br />
la condición<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
(l(C, u), r c , u) = 0 (12)<br />
que representa la consición analítica necesaria y suficiente para que una<br />
recta que no pase por C sea eje <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> alguno <strong>de</strong> sus<br />
puntos. Llamaremos a la ecuación anterior test <strong>de</strong> inercia o condición <strong>de</strong><br />
inercialidad.<br />
Supuesto que (12) es satisfecha, entonces sólo hay un valor <strong>de</strong> λ que<br />
satisfaga la ecuación (10), por lo que sólo existe un punto O ′ <strong>de</strong> cuyo<br />
elipsoi<strong>de</strong> <strong>de</strong> inercia sea eje la recta δ. Dicho punto está posicionado por:<br />
OO ′ = − (l(C, u), u, u π)<br />
M(r c , u, u π ) u<br />
don<strong>de</strong> u π es un vector perpendicular al plano r c , r c y el <strong>de</strong>nominador no<br />
pue<strong>de</strong> anularse, pues si lo hiciese, δ contendría C, contra lo supuesto. Las<br />
conclusiones obtenidas pue<strong>de</strong>n con<strong>de</strong>nsarse en el siguiente:<br />
Título<br />
Contenido<br />
◭◭ ◮◮<br />
◭ ◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 29 <strong>de</strong> 41
Teorema 2 Entre todas las rectas que no pasan por el centro <strong>de</strong> <strong>masas</strong> se<br />
distinguen dos únicas clases: las que son ejes <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong><br />
uno y sólo uno <strong>de</strong> sus puntos y las que no lo son <strong>de</strong> ninguno <strong>de</strong> ellos.<br />
Las consi<strong>de</strong>raciones anteriores permiten efectuar una partición <strong>de</strong>l conjunto<br />
<strong>de</strong> rectas <strong>de</strong>l espacio en las siguientes clases :<br />
Rectas principales <strong>de</strong> inercia las rectas que pasan por el centro <strong>de</strong><br />
<strong>masas</strong> C, y son ejes <strong>de</strong>l elipsoi<strong>de</strong> central <strong>de</strong> inercia.<br />
Rectas centrales <strong>de</strong> inercia las rectas que pasan por el centro <strong>de</strong> <strong>masas</strong><br />
y no son ejes <strong>de</strong>l elipsoi<strong>de</strong> central <strong>de</strong> inercia.<br />
Rectas secundarias <strong>de</strong> inercia las rectas que no pasan por el centro <strong>de</strong><br />
<strong>masas</strong> y son ejes <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> alguno <strong>de</strong> sus puntos.<br />
Rectas terciarias <strong>de</strong> inercia las rectas no incluidas en las clases anteriores,<br />
es <strong>de</strong>cir las que no pasan por el centro <strong>de</strong> <strong>masas</strong> y no son eje<br />
<strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> ninguno <strong>de</strong> sus puntos.<br />
De las rectas <strong>de</strong> cada una <strong>de</strong> estas clases, po<strong>de</strong>mos enunciar las siguientes<br />
propieda<strong>de</strong>s:<br />
• Las rectas principales <strong>de</strong> inercia son ejes <strong>de</strong> los elipsoi<strong>de</strong>s <strong>de</strong> inercia<br />
<strong>de</strong> todos sus puntos.<br />
• Las rectas centrales <strong>de</strong> inercia no son ejes <strong>de</strong> los elipsoi<strong>de</strong>s <strong>de</strong> inercia<br />
<strong>de</strong> ninguno <strong>de</strong> sus puntos.<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 30 <strong>de</strong> 41
• Las rectas secundarias <strong>de</strong> inercia son ejes <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong><br />
uno y sólo uno <strong>de</strong> sus puntos. Este punto se llamará punto característico<br />
<strong>de</strong> la recta. A<strong>de</strong>más, todas cumplen la condición expresada<br />
en (12).<br />
• Las rectas terciarias <strong>de</strong> inercia no son ejes <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong><br />
ninguno <strong>de</strong> sus puntos, y no cumplen (12). Algunos <strong>de</strong> estos resultados<br />
aunque más reducidos, son mencionados, aunque <strong>de</strong>ducidos en<br />
un contexto dinámico en [2],[3], quedando <strong>de</strong>mostrado en las líneas<br />
anteriores, que no es necesario partir <strong>de</strong> éste.<br />
Po<strong>de</strong>mos sintetizar las conclusiones anteriores en el siguiente cuadro que<br />
proporciona el procedimiento <strong>de</strong> clasificación <strong>de</strong> una recta cualquiera:<br />
⎧ { l(C, u) × u = 0 ⇐⇒ Rectas principales<br />
⎪⎨ C ∈ δ<br />
δ<br />
{<br />
l(C, u) × u ≠ 0 ⇐⇒ Rectas centrales<br />
(l(C, u), u, rc ) = 0 ⇐⇒ Rectas secundarias<br />
⎪⎩ C ∉ δ<br />
(l(C, u), u, r c ) = 0 ⇐⇒ Rectas terciarias<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
Título<br />
Contenido<br />
◮◮<br />
◭<br />
◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 31 <strong>de</strong> 41
9. Clasificación <strong>de</strong> los planos<br />
Del mismo modo que hemos clasificado las rectas <strong>de</strong>l espacio, po<strong>de</strong>mos<br />
utilizar el campo <strong>de</strong> elipsoi<strong>de</strong>s <strong>de</strong> inercia <strong>de</strong> una distribución arbitraria<br />
<strong>de</strong> <strong>masas</strong> para clasificar el conjunto <strong>de</strong> planos <strong>de</strong>l espacio. Definamos<br />
una propiedad <strong>de</strong> un plano π en un punto O perteneciente a π. Esta<br />
propiedad es la <strong>de</strong> ser plano <strong>de</strong> simetría <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> O.<br />
Geométricamente esta propiedad equivale a la siguiente: la normal a π por<br />
O es eje <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> O. Para esta segunda condición ya<br />
hemos <strong>de</strong>sarrollado una expresión analítica, que po<strong>de</strong>mos volver a utilizar.<br />
Si u es un vector unitario perpendicular a π, y r es el vector <strong>de</strong> posición<br />
<strong>de</strong> O <strong>de</strong>s<strong>de</strong> C, entonces π es un plano <strong>de</strong> simetría <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia<br />
<strong>de</strong> O si y sólo si :<br />
l(O, u) × u = 0<br />
es <strong>de</strong>cir,<br />
l(C, u) × u = m(r · u)(r × u)<br />
Procedamos al estudio separadamente <strong>de</strong> planos que pasen por C, y planos<br />
que no pasen por C.<br />
9.1. π contiene a C<br />
en este caso r c · u π = 0. Aplicando la fórmula <strong>de</strong> Steiner para el vector <strong>de</strong><br />
inercia, se tiene:<br />
l(O, u) = l(C, u) + M nor r c u − M(r c · u)r } {{ } c<br />
=0<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 32 <strong>de</strong> 41
y, por lo tanto,<br />
l(O, u) × u = l(C, u) × u + M nor r c u × u<br />
} {{ }<br />
=0<br />
l(O, u) × u = l(C, u) × u<br />
y π será plano <strong>de</strong> simetría <strong>de</strong>l elipsoi<strong>de</strong> <strong>de</strong> inercia <strong>de</strong> O si y sólo si lo es <strong>de</strong>l<br />
elipsoi<strong>de</strong> central. En este caso, a<strong>de</strong>más, lo sería <strong>de</strong> los <strong>de</strong> todos sus puntos,<br />
pues O es un punto genérico. Esto conduce al<br />
Teorema 3 Entre todos los planos que pasan por el centro <strong>de</strong> <strong>masas</strong> se<br />
distinguen dos únicas clases: los que son planos <strong>de</strong> simetría <strong>de</strong>l elipsoi<strong>de</strong><br />
<strong>de</strong> inercia <strong>de</strong> todos sus puntos y los que no lo son <strong>de</strong> ninguno <strong>de</strong> ellos.<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
9.2. π no contiene a C<br />
En este caso, r c · u π ≠ 0 .Tomaremos un origen O ∈ π e investigaremos<br />
la existencia <strong>de</strong> algún o ,algunos otros puntos O ′ ∈ π cuyos elipsoi<strong>de</strong>s <strong>de</strong><br />
inercia tengan a π como plano <strong>de</strong> simetría.<br />
El vector <strong>de</strong> inercia en O es, por la fórmula <strong>de</strong> Steiner,<br />
l(O, u) = l(C, u) + M nor r c u − M(r c · u)r c<br />
El vector <strong>de</strong> inercia en O ′ será<br />
l(O ′ , u) = l(C, u) + M nor r ′ cu − M(r ′ c · u)r ′ c<br />
Título<br />
Contenido<br />
◭◭ ◮◮<br />
◭ ◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 33 <strong>de</strong> 41
don<strong>de</strong> r ′ c = O ′ O + OC = O ′ O + r c<br />
Multiplicando las dos expresiones anteriores vectorialmente por u queda<br />
{ l(O, u) × u = l(C, u) × u − M(rc · u)r c × u<br />
l(O ′ , u) × u = l(C, u) − M(r ′ c · u)r ′ c × u<br />
igualando l(O ′ , u) × u = 0 y restando, queda<br />
l(O, u) × u = M[(r ′ c · u)r ′ c × u − (r c · u)r c × u]<br />
Sustituyendo<br />
l(O, u) × u = M[((O ′ O + r c ) · u)(O ′ O + r c ) × u − (r c · u)r c × u]<br />
l(O, u) × u = M[(O } ′ O {{· u}<br />
) + (r c · u)O ′ O × u + r c × u − (r c · u)r c × u]<br />
=0<br />
l(O, u) × u = M[(r c · u)O ′ O × u + (r c · u)r c × u − (r c · u)r c × u]<br />
} {{ }<br />
=0<br />
l(O, u) × u = M(r c · u)O ′ O × u<br />
El primer miembro es un vector paralelo al plano π, por lo que, dado que<br />
M(r c · u) ≠ 0, siempre existirá un vector O ′ O <strong>de</strong>l plano π que satisfaga<br />
la ecuación anterior. Este vector pue<strong>de</strong> <strong>de</strong>spejarse premultiplicando ambas<br />
expresiones por u<br />
u × (l(O, u) × u) = M(r c · u)u × (O ′ O × u)<br />
u × (l(O, u) × u) = M(r c · u)O ′ O<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 34 <strong>de</strong> 41
y, por tanto, siempre existirá una y sólo una solución para O ′ ∈ π, dada<br />
por la expresión:<br />
O ′ u × (l(O, u) × u)<br />
O =<br />
M(r c · u)<br />
De acuerdo con estas consi<strong>de</strong>raciones, po<strong>de</strong>mos efectuar una partición<br />
<strong>de</strong>l conjunto <strong>de</strong> los planos según sus características <strong>de</strong> inercia en las siguientes<br />
clases:<br />
planos principales <strong>de</strong> inercia los planos <strong>de</strong> simetría <strong>de</strong>l elipsoi<strong>de</strong> central<br />
<strong>de</strong> inercia.<br />
planos centrales <strong>de</strong> inercia los planos que pasan por el ceentro <strong>de</strong> <strong>masas</strong><br />
y no son planos <strong>de</strong> simetría <strong>de</strong>l elipsoi<strong>de</strong> central <strong>de</strong> inercia.<br />
planos secundarios <strong>de</strong> inercia los planos que no pasan por el centro<br />
<strong>de</strong> <strong>masas</strong>.<br />
Po<strong>de</strong>mos enunciar las propieda<strong>de</strong>s <strong>de</strong> inercia <strong>de</strong> los planos <strong>de</strong> cada una<br />
<strong>de</strong> estas clases a la luz <strong>de</strong> los resultados anteriores.<br />
• los planos principales <strong>de</strong> inercia son planos <strong>de</strong> simetriaa <strong>de</strong> los elipsoi<strong>de</strong>s<br />
<strong>de</strong> inercia <strong>de</strong> todos sus puntos.<br />
• los planos centrales <strong>de</strong> inercia no son planos <strong>de</strong> simetría <strong>de</strong>l elipsoi<strong>de</strong><br />
<strong>de</strong> inercia <strong>de</strong> ninguno <strong>de</strong> sus puntos.<br />
• los planos secundarios <strong>de</strong> inercia son planos <strong>de</strong> simetría <strong>de</strong>l elipsoi<strong>de</strong><br />
<strong>de</strong> inercia <strong>de</strong> uno y sólo uno <strong>de</strong> sus puntos.<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 35 <strong>de</strong> 41
10. Clasificación <strong>de</strong> las direcciones<br />
Una vez clasificados el conjunto <strong>de</strong> rectas y el <strong>de</strong> planos <strong>de</strong>l espacio, hacemos<br />
una abstracción y queremos analizar qué tipos <strong>de</strong> direcciones <strong>de</strong>l espacio<br />
nos <strong>de</strong>fine las características <strong>de</strong> inercia <strong>de</strong> una distribución dada <strong>de</strong><br />
<strong>masas</strong>.<br />
El punto <strong>de</strong> partida lo tenemos resuelto, pues la clasificación <strong>de</strong> rectas<br />
y planos nos proporciona una clasificación natural <strong>de</strong> las direcciones <strong>de</strong>l<br />
espacio. Según esto, <strong>de</strong>finamos la partición <strong>de</strong>l conjunto <strong>de</strong> direcciones <strong>de</strong>l<br />
espacio en las siguientes clases:<br />
ortodirecciones planarias <strong>de</strong> inercia son las direcciones perpendiculares<br />
a los planos principales <strong>de</strong> inercia, es <strong>de</strong>cir, las paralelas a alguna<br />
recta principal <strong>de</strong> inercia.<br />
ortodirecciones axiales <strong>de</strong> inercia son las direcciones perpendiculares<br />
a las rectas principales <strong>de</strong> inercia, o lo que es lo mismo, paralelas a<br />
algún plano principal <strong>de</strong> inercia, no incluidas en la clase anterior.<br />
clinodirecciones <strong>de</strong> inercia las direcciones no incluidas en las clases anteriores.<br />
Una vez clasificadas las direcciones, nos hacemos las siguientes preguntas:<br />
dada un dirección <strong>de</strong>l espacio . . .<br />
• ¿ en qué puntos <strong>de</strong>l espacio existirá un elipsoi<strong>de</strong> <strong>de</strong> inercia con un eje<br />
paralelo a esa dirección ?<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 36 <strong>de</strong> 41
• ¿ dón<strong>de</strong> estarán, si las hay, las rectas principales, secundarias, centrales,<br />
terciarias <strong>de</strong> inercia, paralelas a la misma ?<br />
Contestaremos separadamente a estas preguntas, según el tipo <strong>de</strong> dirección<br />
<strong>de</strong> que se trate.<br />
En los tres casos (ortodirecciones planarias, ortodirecciones axiales, clinodirecciones<br />
<strong>de</strong> inercia), tomaremos un sistema <strong>de</strong> coor<strong>de</strong>nadas con origen<br />
en el centro <strong>de</strong> <strong>masas</strong> C, y con eje x, paralelo a la dirección en estudio.<br />
Los ejes y, z se <strong>de</strong>finirán <strong>de</strong> forma distinta en cada situación.<br />
10.1. Ortodirecciones planarias <strong>de</strong> inercia<br />
En este caso po<strong>de</strong>mos tomar los ejes y, z según dos rectas principales <strong>de</strong><br />
inercia. La expresión matricial <strong>de</strong>l tensor central <strong>de</strong> inercia en esta base<br />
tiene forma diagonal:<br />
⎛<br />
(I) = ⎝<br />
I ∗ x 0 0<br />
0 I ∗ y 0<br />
0 0 I ∗ z<br />
El tensor <strong>de</strong> inercia <strong>de</strong>l punto O, posicionado por sus coor<strong>de</strong>nadas (ξ, η, ζ),<br />
aplicando la expresión matricial <strong>de</strong> las fórmulas <strong>de</strong> Steiner queda:<br />
⎛<br />
⎞<br />
⎛ ⎞ ⎛<br />
⎞<br />
Ix ∗ 0 0<br />
1 0 0<br />
ξ 2 ξη ζξ<br />
(I) = ⎝ 0 Iy ∗ 0 ⎠+M(ξ 2 +η 2 +ζ 2 ) ⎝ 0 1 0 ⎠−M ⎝ ξη η 2 ηζ ⎠<br />
0 0 Iz<br />
∗ 0 0 1 −ζξ ηζ ζ 2<br />
⎞<br />
⎠<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 37 <strong>de</strong> 41
⎛<br />
⎞<br />
Ix ∗ 0 0<br />
(I) = ⎝ 0 Iy ∗ 0 ⎠ + M<br />
0 0 Iz<br />
∗<br />
buscamos los puntos O en los que<br />
es <strong>de</strong>cir<br />
I ·u ‖ u ⇒ ⎝<br />
Por lo tanto<br />
⎛<br />
y es evi<strong>de</strong>nte que:<br />
I ∗ x 0 0<br />
0 I ∗ y 0<br />
0 0 I ∗ z<br />
⎞<br />
⎛<br />
⎝<br />
l(O, i) × i = 0<br />
⎠+M ⎝<br />
⎛<br />
= ⎝<br />
⎛<br />
⎞<br />
η 2 + ζ 2 −ξη −ζξ<br />
−ξη ζ 2 + ξ 2 −ηζ ⎠<br />
−ζξ −ηζ ξ 2 + η 2<br />
⎞ ⎛<br />
η 2 + ζ 2 −ξη −ζξ<br />
−ξη ζ 2 + ξ 2 −ηζ ⎠ ⎝<br />
−ζξ −ηζ ξ 2 + η 2<br />
I ∗ x + M(η 2 + ζ 2 )<br />
−Mξη<br />
−Mζξ<br />
{ { ξη = 0 ξ = 0 plano ⊥ x<br />
ξζ = 0 ⇒ y, z = 0 eje x<br />
1<br />
0<br />
0<br />
⎞<br />
⎠ =<br />
Teorema 4 el lugar geométrico <strong>de</strong> puntos cuyo elipsoi<strong>de</strong> <strong>de</strong> inercia tiene<br />
un eje con ortodirección planaria está compuesto por los <strong>de</strong> la recta principal<br />
paralela y los <strong>de</strong>l plano principal perpendicular.<br />
lo que conduce al<br />
Corolario 1 el conjunto <strong>de</strong> rectas secundarias <strong>de</strong> inercia paralelas a esa<br />
dirección es tá compuesto por todas ellas excepto la principal.<br />
⎞<br />
⎠<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 38 <strong>de</strong> 41
10.2. Ortodirecciones axiales <strong>de</strong> inercia<br />
En este caso, suponemos que el eje x es perpendicular a una recta principal<br />
<strong>de</strong> inercia y siempre se podrá tomar el eje z según una recta principal,<br />
aunque nunca podremos tomar el eje y según otra. La expresión matricial<br />
<strong>de</strong>l tensor <strong>de</strong> inercia <strong>de</strong>l punto O(ξ, η, ζ) queda:<br />
⎛<br />
(I) = ⎝<br />
⎛<br />
(I)(i) = [ ⎝<br />
Ix ∗ −Pxy ∗ 0<br />
−Pxy ∗ Iy ∗ 0<br />
0 0 Iz<br />
∗<br />
I ∗ x −P ∗ xy 0<br />
−P ∗ xy I ∗ y 0<br />
0 0 I ∗ z<br />
⎞<br />
⎠ + M ⎝<br />
⎞<br />
(I)(i) = ⎝<br />
Por lo que será necesario que<br />
{<br />
⎠+M ⎝<br />
⎛<br />
⎛<br />
⎛<br />
⎞<br />
η 2 + ζ 2 −ξη −ζξ<br />
−ξη ζ 2 + ξ 2 −ηζ ⎠<br />
−ζξ −ηζ ξ 2 + η 2<br />
⎞<br />
η 2 + ζ 2 −ξη −ζξ<br />
−ξη ζ 2 + ξ 2 −ηζ ⎠]<br />
−ζξ −ηζ ξ 2 + η 2<br />
I x + M(η 2 + ζ 2 )<br />
−P xy − Mξη<br />
−Mζξ<br />
ξη = − P xy<br />
M<br />
ζ = 0<br />
Por lo tanto, se pue<strong>de</strong> enunciar el siguiente:<br />
Teorema 5 el conjunto <strong>de</strong> puntos cuyo elipsoi<strong>de</strong> <strong>de</strong> inercia tiene un eje<br />
según una ortodirección axial dada es una hipérbola contenida en el plano<br />
principal paralelo a la ortodirección estudiada. Esta hipérbola se conoce<br />
como hipérbola <strong>de</strong> inercia <strong>de</strong> la dirección.<br />
⎞<br />
⎠<br />
⎛<br />
⎝<br />
1<br />
0<br />
0<br />
⎞<br />
⎠<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 39 <strong>de</strong> 41
<strong>de</strong>l que se <strong>de</strong>duce el<br />
Corolario 2 El conjunto <strong>de</strong> rectas secundarias paralelas a una ortodirección<br />
axial dada es el conjunto <strong>de</strong> rectas contenidas en el plano <strong>de</strong> la<br />
hipérbola, excepto la central.<br />
10.3. Clinodirecciones <strong>de</strong> inercia<br />
Para éstas, nunca podremos elegir ejes y, z según rectas principales, y los<br />
tomaremos según dos ejes cualesquiera perpendiculares entre sí y al eje x.<br />
La expresión matricial <strong>de</strong>l tensor <strong>de</strong> inercia <strong>de</strong> O(ξ, η, ζ) queda:<br />
⎛<br />
(I) = ⎝<br />
Ix ∗ −Pxy ∗ −Pzx<br />
∗<br />
−Pxy ∗ Iyy ∗ −Pyz<br />
∗<br />
−Pzx ∗ −Pyz ∗ Iz<br />
∗<br />
⎞<br />
⎠ + M ⎝<br />
con lo que <strong>de</strong>be cumplirse que:<br />
{<br />
ξη = −<br />
P ∗ xy<br />
M<br />
ζξ = − Pzx<br />
M<br />
⎛<br />
⎞<br />
η 2 + ζ 2 −ξη −ζξ<br />
−ξη ζ 2 + ξ 2 −ηζ ⎠<br />
−ζξ −ηζ ξ 2 + η 2<br />
ecuación <strong>de</strong> una hipérbola equilátera una <strong>de</strong> cuyas asíntotas es la recta central<br />
<strong>de</strong> inercia paralela a la clinodirección <strong>de</strong> partida. Pue<strong>de</strong> <strong>de</strong>mostrarse<br />
que el plano <strong>de</strong> la hipérbola es el <strong>de</strong>terminado por la dirección y la normal<br />
al elipsoi<strong>de</strong> central <strong>de</strong> inercia en el punto en que éste es cortado por la<br />
recta central paralela a la dirección original.<br />
Por lo tanto, se enuncia un teorema similar al <strong>de</strong> los casos anteriores:<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
◭◭<br />
◭<br />
Título<br />
Contenido<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
◮◮<br />
◮<br />
Página 40 <strong>de</strong> 41
Teorema 6 el conjunto <strong>de</strong> puntos cuyo elipsoi<strong>de</strong> <strong>de</strong> inercia tiene un eje<br />
según una clinodirección dada es una hipérbola contenida en el plano central<br />
paralelo a la dirección y a la normal al elipsoi<strong>de</strong> central. Esta hipérbola se<br />
conoce como hipérbola <strong>de</strong> inercia <strong>de</strong> la dirección.<br />
y, como consecuencia, se llega al:<br />
Corolario 3 El conjunto <strong>de</strong> rectas secundarias paralelas a una clinodirección<br />
dada es el conjunto <strong>de</strong> rectas contenidas en el plano <strong>de</strong> la hipérbola,<br />
excepto la central.<br />
En la Tabla I se recogen los resultados anteriores.<br />
Clase <strong>de</strong> dirección Puntos característicos<br />
Ortoplanaria recta ppal ‖ y plano ppal ⊥<br />
Ortoaxial hipérbola equilátera en plano ppal ‖<br />
Clinodirección<br />
hipérbola equilátera<br />
DFAII<br />
M ɛ cFunN ɛ t<br />
Título<br />
Contenido<br />
◭◭<br />
◭<br />
◮◮<br />
◮<br />
Volver Atrás<br />
Cerrar<br />
Salir<br />
Página 41 <strong>de</strong> 41