Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

LM(ndofe) arreglo que relaciona cada grado de libertad local con las ecuación global correspondiente,
si es un valor positivo corresponde a un grado de libertad efectivo y si es un valor
negativo corresponde a un grado de libertad restringido asociado con el correspondiente valor
en KNOWN.
Rutina STIFFE
Esta rutina calcula la matriz de rigidez de un elemento por integración numérica usando la
expresión
K e =
NGaus
∑
G=1
Las variables definidas localmente son:
B T (ξ G , η G ) D B (ξ G , η G ) w G |J G |
JAC determinante jacobiano de la transformación, luego multiplicado por el peso del punto de
integración.
ELCOD(dimen,nnode) coordenadas nodales del elemento considerado
D(nstre,nstre) matriz de elasticidad (o la correspondiente para problemas de campo) evaluada
en la rutina DMATRX en función de las características del material. Constante en cada elemento.
B(ndofe,nstre) traspuesta de la matriz B en el punto de integración considerado, evaluada en
la rutina BMATRX en función del gradiente cartesiano de las funciones de forma (y del tipo de
problema)
Rutina BMATRX
Esta rutina evalúa en cada punto de integración el operador B T que relaciona las variables
nodales con las variables derivadas asociadas al flujo mediante las ecuaciones constitutivas.
Ecuación de Laplace
[ N
B I I
′ 1
=
N I
′ 2
]
Estado plano de tensiones y deformaciones
⎡
B I = ⎣ N ⎤
I
′ 1
0
0 N I ⎦
′ 2
N I
′ 2 N I
′ 1
Estado de deformación axilsimétrica
⎡
B I =
⎢
⎣
N I
′ 1 0
0 N I
′ 2
N I
′ 2 N I
′ 1
N 1
x 1
0
⎤
⎥
⎦
Flexión de placas incluyendo el corte
170
⎡
B I =
⎢
⎣
⎤
0 N I
′ 1
0
0 0 N I
′ 2
0 N I
′ 2 N I
′ 1 ⎥
′ 1
N I 0 ⎦
′ 2 0 N I
N I
N I