1

20 – Dinámica + elementos finitos(caso lineal)

Diego Andrés Alvarez MarínProfesor Asociado

Universidad Nacional de ColombiaSede Manizales

2

Ecuaciones de la elástodinámicaLas ecuaciones diferenciales de equilibrio son:

Las ecuaciones dinámicas de Cauchy-Navier son:

3

Interpolación de los desplazamientos, velocidades y

aceleraciones en el caso dinámico

4

Energía cinética

Matriz de masaconsistente

5

Energía potencial elástica

6

Trabajo realizado por las fuerzas externas Energía externa ejercida

respectivamente por las fuerzas másicas, superficiales y puntuales actuando sobre el sólido

Vector de fuerzas nodales de equilibrio

Vector de fuerzas nodales equivalentes

7

Funcional lagrangiano

8

Fuerzas conservativas

9

Propiedades del variacional primero

10

Principio de Hamilton

11

Las condiciones iniciales y finales son conocidas (condición 3 del principio de Hamilton)

12

13

Principio de D'Alembert

Este principio permite extender la segunda ley de Newton sobre equilibrio estático al caso dinámico por medio de la simple consideración que la fuerza de inercia (que para un cuerpo de masa constante es igual al producto de la masa por la aceleración) es una fuerza de sentido opuesto al sentido positivo de la aceleración.

https://en.wikipedia.org/wiki/D%27Alembert%27s_principle

15

Deducción alternativa (e incorrecta, pero popular):modificación del principio de los trabajos virtuales

Utilizando el principio de D'Alembert se puede incluir las fuerzas de inercia en el principio del trabajo virtual como términos adicionales a las fuerzas másicas:

donde:

y ρ es la masa por unidad de volumen (densidad).

DEDUCCION INCORRECTA!Recuerde que la deducción correcta se hace usando el principio de Hamilton

16

Por lo tanto, haciendo:

y reemplazando en el PTV obtenemos (aquí por cuestión de espacio solo consideramos las fuerzas de inercia y las otras fuerzas másicas b):

17

Ecuación de movimiento de un sistema no amortiguado

18

Matriz de masa consistente (consistent mass matrix)

La matriz M que hemos encontrado es la llamada matriz de masa consistente:

donde ρ es la masa por unidad de volumen (densidad).

19

Matriz de masa concentrada (lumped mass matrix)

En los algoritmos de dinámica se requiere frecuentemente invertir la matriz M; por dicha razón, se emplea frecuentemente la matriz de masa concentrada.

Este esquema crea una matriz diagonal empleando así menos recursos computacionales. Existen varias formas de generar esta matriz. Un esquema correcto debe preservar la masa total correcta para el elemento y preservar correctamente el centro de gravedad del mismo. Sería deseable también preservar el primer y segundo momento de inercia, pero esto raramente se logra.

20

Matriz de masa concentrada (lumped mass matrix)

No existe un procedimiento único para calcular la matriz de masa concentrada. Existen en la literatura varios métodos para calcular dicha matriz, por ejemplo:

● Método de Archer (Direct Mass Lumping)● Método HRZ● Otros métodos

21

Método de Archer (1963, 1965)

Matriz de masa concentrada (lumped) de una viga. Observe que la masa soloestá dispuesta en el grado de libertad traslacional

Fue el primer método para crear la matriz de masa concentrada.

Este método solo asigna masas a los grados e libertad traslacionales

Ejemplo:

22

DLMM=DirectLumpedMassMatrix

23

Método HRZ(Hinton-Rock-Zienkiewicks, 1976)

24

25

Condensación de HRZ para un elemento de viga

26

Otros métodos

● Isaac Fried, David S. Malkus (1975). Finite element mass matrix lumping by numerical integration with no convergence rate loss, International Journal of Solids and Structures, Volume 11, Issue 4, Pages 461-466, http://dx.doi.org/10.1016/0020-7683(75)90081-5

● Shen R. Wu (2006). Lumped mass matrix in explicit finite element method for transient dynamics of elasticity, Computer Methods in Applied Mechanics and Engineering, Volume 195, Issues 44–47, , Pages 5983-5994. http://dx.doi.org/10.1016/j.cma.2005.10.008

27

¿Cuál tipo de formulación utilizar?La matriz de masa consistente en general dará resultados más precisos (pero no mucho más). Algunos otros autores no afirman esto.

El cálculo con la matriz de masa concentrada es más rápido ya que la inversa de una matriz diagonal también es diagonal.

Si una malla es lo suficientemente refinada, los resultados serán aproximadamente los mismos independientemente si se usa matriz de masa consistente o concentrada.

28

La matriz de masa consistente hace una representación más precisa de las propiedades inerciales de la estructura y además produce frecuencias naturales mayores o iguales que las exactas; además utilizando un número suficiente de elementos se pueden obtener resultados muy precisos.

La representación de masa concentrada no es tan precisa (algunos autores no afirman esto) y puede producir frecuencias naturales más altas o más bajas que las exactas. Sin embargo, para problemas grandes la representación de masa concentrada supone un ahorro considerable en el cálculo ya que da lugar a una matriz de masa diagonal.

29

En conclusión:

No hay una respuesta definitiva a que tipo de matriz utilizar. La matriz de

masa consistente no es necesariamente la que brinda mejores

resultados.

30

Particularizaciones de la matriz de masa consistente

Caso tridimensional:

Caso bidimensional (tensión/deformación plana):

Caso axisimétrico:

Caso barra unidimensional:

31

Formulación para barras sometidas a fuerzas axiales

32

b

33

34

Matriz de rigidez del elemento cercha

35

Matriz de masa del elemento cercha

36

37

Formulación para vigas

38

Formulación para el elemento viga de dos nodos

Matriz de masa consistente

Matriz de masa concentrada (lumped)Observe que la masa solo está dispuestaen el grado de libertad translacional

39

Formulación para elementos de pórtico

40

Matriz de rigidez de un elemento prismático sometido en sus extremos a carga axial,

flexión y cortante

41

Matriz de masa para un elemento bajo movimiento axial

Matriz de masa para el elemento viga

Y superponiendo las contribuciones:

Matriz masa consistente para elemento pórtico 2D

42

Matriz de masa condensada para elemento pórtico 2D

43

Matriz de transformación

45

46

Elemento triangular de tres nodos

t

t

47

Elemento rectangular de cuatro nodos

t

t

48

49

Matriz de masa consistente para un elemento hexahédrico de 8 nodos

50

Análisis modal de sistemas no amortiguados

Dicha ecuación se puede resolver ya sea con:● Análisis modal● Integración directa

– Método de Newmark

– Método de Wilson

– Método de las diferencias centrales

– etc.

51

Análisis modal

Métodos numéricos para calcular valores y vectores propios:

● Método de Jacobi

● Método de Given

● Método de Housholder

● Método de la bisección (usando secuencias de Sturm)

● Iteración inversa

● Método QR

● Método de Lanczos

52

Análisis modal con MATLAB

53

Análisis modal de sistemas amortiguados

54

Ecuación de movimiento para vibración forzada para sistemas de

múltiples grados de libertad