1
01 – Repaso de análisis estructural matricial
Diego Andrés Alvarez MarínProfesor Asociado
Universidad Nacional de ColombiaSede Manizales
2
Contenido
● Solución al sistema de ecuaciones q=Ka-f● Conceptos básicos del análisis matricial de
estructuras de barras● Cerchas● Marcos● Etapas básicas del análisis de un sistema de
barras
3
El método de los elementos finitos
Es una técnica numérica para la solución aproximada de ecuaciones diferenciales e integrales.
En el ámbito del análisis estructural (estructuras, geotecnia, pavimentos) es un poderoso método para la estimación de los esfuerzos, deformaciones y desplazamientos de una estructura bajo la acción de un conjunto de cargas.
4
Historia del método de los EFs● 1934-1938: Collar, Duncan: primer artículo sobre análisis matricial de
estructuras
● 1943: Richard Courant (New York University): primer artículo matemático
● 1954: Argyris: desarrolla el concepto de ensamblaje matricial
● 1956: Turner, Clough, Martin, Topp (Boeing): primer artículo en ingeniería
● 1959: Turner: presenta un artículo sobre el método matricial de la rigidez, tal y como lo conocemos hoy
● 1960: Clough (Boeing, Berkeley) define el término “elemento finito”.
● 1965: La NASA invierte 4E6 USD (30E6 USD de hoy) para el desarrollo del software NASTRAN
● 1965: Se empieza a aplicar a otros campos de la ingeniería
● 1967: Zienkiewicz (Swansea): primer libro sobre EFs
● 1970: Se empieza a desarrollar el software ANSYS
● 1973: Strang, Fix (M.I.T.): primer libro sobre EFs desde el punto de vista matemático
5
Solución al sistema de ecuaciones
6
Solución al sistema de ecuaciones
7
Resolviendo
● Cholesky● Métodos que tienen en cuenta las matrices
ralas (matrices sparse). Muchos de estos métodos son iterativos:– PARDISO http://www.pardiso-project.org/
– Método del gradiente conjugado
– Comando linsolve de MATLAB
– Comando \ de MATLAB– http://www.mathcom.com/corpdir/techinfo.mdir/q207.html
8
Ancho de
banda
La numeración de los nodos debe hacerse de modo tal que el ancho de banda sea tan pequeño como sea posible. Existen algoritmos especializados que hacen esta labor (como el algoritmo invertido de Cuthill-McKee). En algunos programas de elementos finitos se puede hacer click en un comando que renumera los nodos de modo que el ancho de banda se reduzca.
9
Ensamblaje matricial
10
11
12
13
Ejemplo 1.1 Oñate
14
Coordenadas locales y globalesMatrices LaG y GaL
11 12 3 1 2
21 2 1 1 32 3 2 2 3
31 3 3 3 42 4
LaG LaG(e, nudo local)barra nudo local nudo global
nudo local
barra (e)
1 2 3 41 1 x 2 x2 x 1 2 x3 x x 1 2
GaL(e, nudo global)
nudo global
barra (e)
15
Solución en MATLAB 1/2
16
Solución en MATLAB 2/2
17
18
19
Cerchas
20
Ejemplo 11.3 Uribe Escamilla
21
22
23
24
Resultados
25
Matriz de rigidez de un elemento prismático sometido en sus extremos a carga axial,
flexión y cortante
26
Matriz de transformación
27
Algunas fuerzas nodales equivalentes
28
Fuerzas nodales equivalentes
29
Cálculo de fuerzas nodales equivalentes
Obteniendo:
30
Ejemplo 11.23 Uribe Escamilla
Mirar solución en la sección código de la WIKI
31
32
idx
Ejemplo: idx para la barra 2:
33
Rótulas interiores a un
pórticoRótula a la derecha del elemento:
Rótula a la izquierda del elemento:
Ver programa:
c1_K_elemento_empotrado_rodillo.m
34
Ver: c1_ejemplo_rotula1.m
35
36
37
38
39
40
41
EJEMPLO
42
a(e)
q
q
Elemento de barra 3D
43
Elemento de barra 3D
44
Matriz de transformación de coordenadas T(e)
45
46
x,y,z = ejes locales X,Y,Z = ejes globales
Tomado de Liu&Quek, p118
47
48
49
50
51
Fórmula de rotación de Rodrigues
52
Etapas básicas del análisis de un sistema de barras
● Etapa de preproceso: – definición de geometría del problema
– definición de las propiedades del material
– definición de las condiciones de carga
– definición de la malla de elementos finitos (numeración local y global de los grados de libertad, de los nodos y de los elementos)
53
Etapas básicas del análisis de un sistema de barras
● Etapa de cálculo
– Cálculo de las matrices de rigidez K(e) y los vectores de fuerza nodales f(e) de cada elemento del sistema
– Ensamblaje de la matriz de rigidez global K y vector de fuerzas global f
54
Etapas básicas del análisis de un sistema de barras
● Etapa de cálculo– Solución del sistema de ecuaciones
– para calcular los desplazamientos en los nudo a (Cholesky, métodos especiales que manejan matrices ralas)
55
Etapas básicas del análisis de un sistema de barras
● Etapa de cálculo– A partir de los valores de las incógnitas en los
nudos obtener información sobre los esfuerzos y deformaciones en las barras
● Etapa de postproceso– Presentación gráfica de la información
Top Related