Scribd
Cargar
69 vistas2 páginas

Matrices Esfuerzo-Deformación en Elementos Finitos

El documento presenta el cálculo de las matrices de esfuerzo-deformación para dos elementos triangulares. Se dan las coordenadas de los nodos de cada elemento y se calculan el jacobiano y su determinante. Luego, se establece la forma general de la matriz B y se resuelve para cada elemento triangular.

Cargado por

Angel Ferreyra
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
69 vistas2 páginas

Matrices Esfuerzo-Deformación en Elementos Finitos

El documento presenta el cálculo de las matrices de esfuerzo-deformación para dos elementos triangulares. Se dan las coordenadas de los nodos de cada elemento y se calculan el jacobiano y su determinante. Luego, se establece la forma general de la matriz B y se resuelve para cada elemento triangular.

Cargado por

Angel Ferreyra
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

José Ángel Ferreyra Ríos

Actividad 7
IAE
Análisis de elemento finito
M1
26/05/2021

Ejemplo 6.3
Encuentra las matrices esfuerzo-deformación nodal para cada elemento mostrado
en la Fig.E6.3. Utiliza los números locales dados en las esquinas.

Podemos saber las coordenadas en los ejes x y y a partir de las medidas que
se muestran en la figura, tomando como referencia de origen al nodo 1.
Se tomará como caso de estudio el elemento e=1.

x1 ≔ 0 x2 ≔ 0 x3 ≔ -3
y1 ≔ 0 y2 ≔ 2 y3 ≔ 0

El Jacobiano y su determinante:
⎡ (x1 - x3)) ((y1 - y3)) ⎤ ⎡ 3 0 ⎤
J1 ≔ ⎢ ( =
⎣ ((x2 - x3)) ((y2 - y3)) ⎥⎦ ⎢⎣ 3 2 ⎥⎦

detJ1 ≔ ⎛J1 ⋅ J1 ⎞ - ⎛J1 ⋅ J1 ⎞ = 6


⎝ 0,0 1 , 1⎠ ⎝ 1,0 0 , 1⎠

Si sabemos que la matriz B es de la forma

⎡ y2 - y3 0 y3 - y1 0 y1 - y2 0 ⎤
1 ⎢
B1 ≔ ―― 0 x3 - x2 0 x1 - x3 0 x2 - x1 ⎥
detJ1 ⎢ x3 - x2 y2 - y3 x1 - x3 y1 - y3 x2 - x1 y1 - y2 ⎥
⎣ ⎦

Entonces:

⎡ 0.333 0 0 0 -0.333 0 ⎤
⎢
B1 = 0 -0.5 0 0.5 0 0 ⎥
⎢ ⎥
⎣ -0.5 0.333 0.5 0 0 -0.333 ⎦

Prof. Dr. Francisco Ramírez Cruz


Facultad de Ingeniería Mecánica y Eléctrica
Universidad Autónoma de Nuevo León
José Ángel Ferreyra Ríos
Actividad 7
IAE
Análisis de elemento finito
M1
26/05/2021

Coordenadas para el elemento triangular número 2, tomando como referencia


el nodo 1 para el origen:

x1 ≔ 0 x2 ≔ 0 x3 ≔ 3
y1 ≔ 0 y2 ≔ -2 y3 ≔ 0

Obtenemos el Jacobiano y su determinante.

⎡ (x1 - x3)) ((y1 - y3)) ⎤ ⎡ -3 0 ⎤


J2 ≔ ⎢ ( =
⎣ ((x2 - x3)) ((y2 - y3)) ⎥⎦ ⎢⎣ -3 -2 ⎥⎦

detJ2 ≔ ⎛J1 ⋅ J1 ⎞ - ⎛J1 ⋅ J1 ⎞ = 6


⎝ 0,0 1 , 1⎠ ⎝ 1,0 0 , 1⎠

De la misma manera se establece la forma y solución de la matriz de


esfuerzo-deformación nodal pero ahora para el segundo elemento:
⎡ y2 - y3 0 y3 - y1 0 y1 - y2 0 ⎤
1 ⎢
B2 ≔ ―― 0 x3 - x2 0 x1 - x3 0 x2 - x1 ⎥
detJ2 ⎢ x3 - x2 y2 - y3 x1 - x3 y1 - y3 x2 - x1 y1 - y2 ⎥
⎣ ⎦

⎡ -0.333 0 0 0 0.333 0 ⎤
⎢
B2 = 0 0.5 0 -0.5 0 0 ⎥
⎢ ⎥
⎣ 0.5 -0.333 -0.5 0 0 0.333 ⎦

Prof. Dr. Francisco Ramírez Cruz


Facultad de Ingeniería Mecánica y Eléctrica
Universidad Autónoma de Nuevo León

También podría gustarte