05 Problemas de elasticidad bidimensional
Diego Andrs Alvarez Marn Profesor Asistente Universidad Nacional de Colombia Sede Manizales
Convencin para los esfuerzos positivos
Deformaciones
Ley de Hooke (relacin esfuerzos deformaciones)
Ley de Hooke para materiales anisotrpicos (relacin esfuerzos-deformaciones)
D
1=x, 2=y, 3=z
5
Tensin plana
Deformacin plana
Ley de Hooke para tensin plana
Ley de Hooke para deformacin plana
Deformaciones iniciales
10
Deformaciones iniciales
11
Esfuerzos iniciales
12
Esfuerzos iniciales
13
Esfuerzos iniciales
14
Malla de elementos finitos
15
Numeracin local vs numeracin global de los nodos de la malla
16
17
18
Reglas para la creacin de la malla de elementos finitos
Es importante reconocer que la malla de elementos finitos representa una idealizacin de la geometra real. Por consiguiente, el anlisis por elementos finitos reproduce el comportamiento de la malla escogida, y no el de la estructura real. Solamente comprobando la convergencia de la solucin podemos estimar el grado de aproximacin de la solucin de elementos finitos a la exacta.
19
20
21
Seleccin del tipo de elemento
En caso que se tenga una cierta idea de la forma polinmica de la solucin, conviene utilizar elementos con funciones de forma del mismo grado que la solucin conocida (rara vez ocurre en la prctica) En zonas donde se intuya que pueden existir gradientes de esfuerzos elevados es ms adecuado utilizar elementos de mayor orden (mtodo p) o mallas ms tupidas (mtodo h).
22
Seleccin del tipo de elemento
Debe evitarse colocar un elemento pequeo contiguo a uno grande. La transicin en tamao debe ser gradual Se recomienda utilizar elementos finitos de pocos nodos (pero no tan pocos!) En el caso de elementos Lagrangianos, tener cuidado con el problema de Runge. Por lo tanto no es bueno escojer tantos nodos.
23
24
Convergencia de la solucin
En lo posible, se deben hacer anlisis con mallas cada vez ms tupidas, de modo que podamos observar si la solucin ha convergido.
25
26
Funciones de forma globales
27
Funciones de forma locales
28
Elemento triangular de tres nodos
29
30
31
Discretizacin del campo de deformaciones
32
33
Discretizacin del campo de deformaciones
34
Discretizacin del campo de tensiones
35
Fuerzas sobre un elemento triangular de tres nodos
36
Las fuerzas de superficie pueden ser de dos tipos: a) Debidas a fuerzas exteriores que actuan sobre los lados del elemento que forman parte del contorno exterior de la estructura b) Debidas a las fuerzas de interaccin entre elementos que se transmiten a travs de lados comunes. Estas ltimas se ignoran desde un principio pues se anulan en el ensamblaje (ya que tienen igual magnitud y direccin, pero sentidos 37 opuestos).
PTV aplicado a un elemento
38
39
40
Hay que destacar que estas expresiones son totalmente generales y, por consiguiente, aplicables a cualquier elemento bidimensional
41
Matriz de rigidez para un elemento triangular de tres nodos
42
43
Vectores de fuerzas nodales equivalentes para un elemento triangular de tres nodos
44
45
46
47
48
Ejercicio de programacin
Considere la viga mostrada, suponiendo que el peso del material es 7.8 kg/m3, E = 200GPa, el coeficiente de Poisson es 0.30 y el espesor de la viga es 10 cm. Calcule los campos de esfuerzos, desplazamientos y deformaciones de la viga
49
Elemento rectangular de 4 nodos
50
Elemento rectangular de 4 nodos
51
Elemento rectangular de 4 nodos
52
Elemento rectangular de 4 nodos
Nota: la matriz de rigidez que aparece en el libro de Oate est mala. Esta es la correcta:
53
54
Elemento rectangular de 4 nodos
Este elemento es muy bueno para problemas de traccin/compresin pura, pero es malo para problemas de flexin debido a su incapacidad natural de adoptar formas curvas. Por esta razn se necesitan mallas muy tupidas 55 para obtener resultados mnimamente aceptables.
Ejercicio de programacin
56
El tringulo de Pascal
57
Tringulo de Pascal
58
Funciones de forma de un elemento rectangular de clase C0 y lados rectos
Estos elementos estn expresados en las llamadas coordenadas naturales o intrnsecas
59
Elemento rectangular lagrangiano
vs
Elemento rectangular serendpito
60
Polinomios de Lagrange
61
Funciones de forma 1D (2 nodos)
62
Funciones de forma 1D (3 nodos)
63
Funciones de forma 1D (4 nodos)
64
Elemento rectangular lagrangiano de 4 nodos
65
66
Elemento rectangular lagrangiano de 9 nodos
67
68
Elemento rectangular lagrangiano de 16 nodos
69
Mostrar programa de MATLAB
70
Elemento rectangular curtico lagrangiano
71
Otros elementos rectangulares de la familia de Lagrange
72
Intercontinuidad elemental
Despus de la deformacin:
Esto implica que si se hace una transicin en el orden de los elementos finitos, se deben utilizar elementos finitos Lagrangianos con diferente nmero de nodos en cada 73 lado para hacer la transicin.
Serendipia (chiripa)
Una serendipia es un descubrimiento o un hallazgo afortunado e inesperado. Se puede denominar as tambin a la casualidad, coincidencia o accidente. El trmino serendipia deriva del ingls serendipity, neologismo acuado por Horace Walpole en 1754 a partir de un cuento persa del siglo XVIII llamado Los tres prncipes de Serendip, en el que los protagonistas, unos prncipes de la isla Serendip (que era el nombre rabe de la isla de Ceiln, la actual Sri Lanka), solucionaban sus problemas a travs de increbles casualidades. NOTA: chiripa si est en el diccionario, serendipia no lo est. Serendipity si existe en el diccionario ingls. 74
Elementos serendpitos rectangulares
Se obtienen de la siguiente manera:
Se selecciona el nmero de nodos de cada lado para definir una variacin lineal, cuadrtica, cbica, etc., sobre dichos lados que garantice la continuidad interelemental. Se escoge el mnimo nmero de nodos en su interior de modo que se obtenga una variacin polinmica de xi y eta completa y simtrica, del mismo grado que la variacin sobre los lados.
75
Elemento rectangular serendpito de 4 nodos
Este elemento pertenece a ambas familias: Lagrangiana y Serendpita
76
Elemento rectangular serendpito de 8 nodos
77
78
Elemento rectangular serendpito de 12 nodos
79
Elemento rectangular serendpito de 17 nodos
80
Elemento rectangular lagrangiano
vs
Elemento rectangular serendpito (GANADORES!)
81
Funciones de forma de elementos triangulares de lados rectos
Estas funciones de forma se caracterizan porque sus funciones de forma contienen exactamente todos los trminos de un polinomio completo de un determinado grado.
1 trmino 3 trminos (lineal) 6 trminos (cuadrtico) 10 trminos (cbico)
82
Coordenadas de rea
83
Coordenadas de rea
Interpolacin paramtrica de la geometra
84
Elemento triangular de 3 nodos
85
Elemento triangular de 6 nodos
86
Elemento triangular de 10 nodos
Mostrar programa de MATLAB
87
Coordenadas naturales del tringulo
88
Cul elemento finito tiene ms precisin?
Los elementos rectangulares son ms precisos que los triangulares para el mismo nmero de grados de libertad. No obstante, los elementos triangulares son mucho ms verstiles que los rectangulares en la discretizacin de geometras complejas. Los elementos de bajo orden son ms sencillos de utilizar aunque en problemas con altos gradientes de esfuerzos la precisin slo se alcanza a cambio de introducir un gran nmero de elementos sencillos, lo que puede hacer obligatorio, e incluso ms rentable en ocasiones, el utilizar elementos de orden ms elevado. 89
La matriz Jacobiana
90
91
El teorema de la funcin inversa
92
El Jacobiano (determinante de la matriz Jacobiana)
El Jacobiano se puede entender como la candidad de estiramiento que una impone una transformacin de variables.
93
94
Cambios de variable en integrales mltiples
95
Ver: http://www.wikimatematica.org/index.php?title=Cambio_de_variables_en_integrales_mltiples
96
La transformacin de coordenadas rectangulares a polares. Se puede notar que el rea de la regin polar es distinta que la de la regin rectangular, lo que justifica la necesidad del jacobiano. 97
Elementos isoparamtricos bidimensionales
98
Elementos cuadrilteros isoparamtricos bidimensionales
99
100
x y
Si se utilizan funciones de forma lineales ningn ngulo interior entre dos lados del elemento sea mayor de 180o. Si las funciones de forma son cuadrticas es necesario adems que los nodos sobre los lados se encuentre en el tercio central de la distancia entre los nodos esquina adyacentes. Para funciones de forma de rdenes superiores es necesario comprobar 101 el signo del Jacobiano.
102
103
El integrando es una funcin racional por lo que debe hacerse uso de la integracin numrica
104
Elementos triangulares isoparamtricos bidimensionales
105
106
Integracin numrica utilizando las cuadraturas de Gauss Legendre sobre dominios cuadrilteros
107
Cuadraturas de Gauss Legendre
108
Cuadraturas de Gauss Legendre
Recuerde que una cuadratura de orden n en cada direccin natural integra exactamente un polinomio de grado 2n-1 o menor en la correspondiente coordenada natural
109
Integracin numrica utilizando las cuadraturas de Gauss Legendre sobre dominios triangulares
David Dunavant, High Degree Efficient Symmetrical Gaussian Quadrature Rules for the Triangle, International Journal for Numerical Methods in Engineering,Volume 21, 1985, pages 1129-1148.
http://people.sc.fsu.edu/~jburkardt/m_src/dunavant/dunavant.html
110
En la tabla la precisin indica el grado del polinomio que se integra exactamente. En los artculos cientficos usualmente se tabulan los Wi de modo que sumen 1. Sin embargo en la frmula se requiere dividir por 1/2. Aqu los pesos ya se han 111 dividido por 1/2.
Seleccin del orden de integracin
En nuestro caso las integrales son funciones racionales y la integracin exacta no es posible. Escoja una nmero de puntos de integracin que integre exactamente los trminos de K correspondientes al polinomio completo contenido en las funciones de forma esta estrategia se llama la cuadratura mnima para obtener la convergencia.
112
113
Si se escojen menos puntos de integracin podran aparecer mecanismos internos.
114
EJEMPLO MATLAB
115
116
117
Integracin numrica de la matriz de rigidez del elemento
Elemento rectangular:
Elemento triangular:
118
Integracin numrica del vector de fuerzas msicas
119
120
Clculo de los esfuerzos principales
121
Teoras de falla
En el mbito de la teora de la elasticidad la falla se produce cuando se produce fluencia en el material. Para calcular el esfuerzo de fluencia las dos teoras de falla ms populares son:
Criterio de falla de Tresca (teora del mximo esfuerzo cortante) Criterio de falla de Von Mises (teora de la mxima energa de deformacin)
122
Tresca
123
Von Mises
124
Comparacin de las superficies de fluencia para los criterios de Von Mises y Tresca en usando las tensiones principales como coordenadas. Observe 125 que el criterio de Tresca es ms conservador
126
127
