02.06.2015 Views

geometría de masas - MecFunNet

geometría de masas - MecFunNet

geometría de masas - MecFunNet

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!