13.05.2013 Views

Apuntes de Mecánica - Grupo de Mecánica Computacional

Apuntes de Mecánica - Grupo de Mecánica Computacional

Apuntes de Mecánica - Grupo de Mecánica Computacional

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.

Aptdo. 11.1. Ecuaciones <strong>de</strong>l Movimiento 11.3<br />

0). Esto equivale a consi<strong>de</strong>rar un sistema autónomo, en el que las únicas<br />

fuerzas actuantes provienen <strong>de</strong> cambios en la propia configuración <strong>de</strong>l sistema<br />

(fuerzas interiores), o <strong>de</strong> la velocidad (resistencias viscosas), no existiendo<br />

fuerzas <strong>de</strong>bidas a agentes exteriores variables con el tiempo. Con estas premisas<br />

realizamos un <strong>de</strong>sarrollo en serie <strong>de</strong> primer or<strong>de</strong>n alre<strong>de</strong>dor <strong>de</strong> la posición<br />

<strong>de</strong> equilibrio, que nos permitrá linealizar la expresión <strong>de</strong> las fuerzas generalizadas<br />

en función <strong>de</strong> coor<strong>de</strong>nadas y velocida<strong>de</strong>s:<br />

Qj(q, ˙q) = Qi| 0<br />

<br />

=0<br />

+ ∂Qi<br />

∂qj<br />

<br />

<br />

0<br />

qj + ∂Qi<br />

∂ ˙qj<br />

<br />

<br />

<br />

0<br />

˙qj + O(q 2 ) + O( ˙q 2 )<br />

Emplearemos la terminología siguiente para los coeficientes que aparecen:<br />

kij = − ∂Qi<br />

∂qj<br />

cij = − ∂Qi<br />

∂ ˙qj<br />

Coeficientes <strong>de</strong> rigi<strong>de</strong>z<br />

Coeficientes <strong>de</strong> amortiguamiento viscoso<br />

Supondremos aquí que las fuerzas provienen <strong>de</strong> un potencial (Qi =<br />

−∂V /∂qi), por lo que los coeficientes <strong>de</strong> rigi<strong>de</strong>z serán:<br />

kij = ∂2V . (11.4)<br />

∂qi∂qj<br />

Es fácil comprobar la simetría <strong>de</strong> estos coeficientes, heredada <strong>de</strong> las propieda<strong>de</strong>s<br />

<strong>de</strong> las <strong>de</strong>rivadas: kij = kji.<br />

Análogamente, para las fuerzas <strong>de</strong>pendientes <strong>de</strong> la velocidad, se pue<strong>de</strong><br />

adoptar la hipótesis <strong>de</strong> que provienen <strong>de</strong> una función R, <strong>de</strong>nominada función<br />

<strong>de</strong> disipación <strong>de</strong> Rayleigh, <strong>de</strong>finida a partir <strong>de</strong> los coeficientes simétricos cij<br />

como:<br />

R = 1<br />

2 cij ˙qi ˙qj, (11.5)<br />

<strong>de</strong> forma que<br />

cij = ∂2R . (11.6)<br />

∂ ˙qi∂ ˙qj<br />

El significado <strong>de</strong> R pue<strong>de</strong> establecerse calculando la tasa <strong>de</strong> energía disipada<br />

por unidad <strong>de</strong> tiempo por las fuerzas viscosas no conservativas:<br />

D = −Q N i ˙qi = (cij ˙qj) ˙qi = 2R<br />

Por el segundo principio <strong>de</strong> la termodinámica esta disipación <strong>de</strong>be ser positiva<br />

o nula, lo que conduce a la condición R ≥ 0, o <strong>de</strong> forma equivalente, a que<br />

los coeficientes cij <strong>de</strong>finen una forma cuadrática semi<strong>de</strong>finida positiva.<br />

Las ecuaciones <strong>de</strong>l movimiento (11.2) quedan pues expresadas como:<br />

ajk ¨qk + [kl, j] ˙qk ˙ql = −kijqj − cij ˙qj + O(q 2 ) + O( ˙q 2 )

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

Saved successfully!

Ooh no, something went wrong!