10 – Análisis de láminas con elementos planos

10 – Análisis de láminas con
elementos planos
Diego Andrés Alvarez Marín
Profesor Asociado
Universidad Nacional de Colombia
Sede Manizales
1
Clasificación de los cascarones
●
●
De acuerdo a la forma de su superficie media:
–
Cascarones planos
–
Cascarones curvos
–
Cascarones axisimétricos
–
Cascarones prismáticos
De acuerdo con la teoría usada:
–
Reissner-Mindlin (cascarones gruesos y delgados)
–
Kirchhoff (cascarones delgados)
2
Resistencia de los cascarones
●
La resistencia de los cascarones está dado por
la mezcla de
–
El estado resistente típico de flexión (soporta
fuerzas de flexión y fuerzas cortantes)
–
El estado resistente de membrana (soporta las
fuerzas axiales)
3
4
5
Elemento finito rectangular de
lámina plana
Las direcciones de
los ejes x' y y'
coinciden
usualmente con las
de los lados del
elemento finito, sin
embargo, pueden
ser arbitrarias
6
Desplazamientos (en coordenadas locales)
7
Desplazamientos
Desplazamientos de un punto de un elemento de cascarón plano en los planos
locales x'z' y y'z', de acuerdo con la teoría de RM.
8
Campo vectorial de movimientos en
coordenadas locales
9
Campo vectorial de deformaciones
(en coordenadas locales)
10
(x',y',z')
Vector de deformaciones
generalizadas de membrana (alargamiento)
Vector de deformaciones
generalizadas de flexión (curvaturas)
Vector de deformaciones
generalizadas de cortante
11
Vector de deformaciones generalizadas
(en coordenadas locales)
Tenga en cuenta que:
donde
12
Nota: esta no es la representación
verdadera de los esfuerzos (ver más adelante)
Campo de
esfuerzos
(en coord.
locales)
Observe que
aquí no se está
teniendo en
cuenta σ'z ya que
según la tercera
hipótesis de RM
su valor es
despreciable.
vector de esfuerzos debidos
a efectos de flexión
vector de esfuerzos debidos
13
a efectos de cortante
transversal
Ley de Hooke
(relación entre esfuerzos y deformaciones)
14
Matrices constitutivas
Esta es la misma matriz
constitutiva utilizada en
tensión plana y en la teoría de
Kirchhoff
Esta es la misma matriz
constitutiva la teoría de R-M, la
cual incluye el coeficiente de
distorsión transversal α=5/6
15
Ley de Hooke
(relación entre esfuerzos y deformaciones)
16
Distribución
de
esfuerzos
en el
espesor del
cascarón
17
Vector de esfuerzos generalizados locales
fuerza
Los momentos y fuerzas que aquí
mostradas son por unidad de longitud
momento
cortante
y
NOTA: El libro de Oñate en
inglés usa una convención de
momentos flectores y
torsores contraria a la aquí
mostrada
x
18
Los momentos y fuerzas que aquí
mostradas son por unidad de longitud
En el caso que se tengan esfuerzos iniciales
debidos a temperatura, el incremento de
temperatura se asume que varía linealmente a
lo largo del espesor de la cáscara a partir de su
definición en ambas caras de la cáscara. Error
de Oñate? Serían deformaciones iniciales.
19
Relaciones entre esfuerzos y
deformaciones generalizadas locales
MCG de acoplamiento membrana-flexión
Todas las
matrices D
son
simétricas
20
MCG de acoplamiento membrana-flexión
Si las propiedades del material no varían con el espesor o si existe simetría de
dichas propiedades con respecto a z'=0, tenemos que:
Ya que:
Se anula para un
material homogéneo
21
Son las mismas matrices que en la teoría de RM
PTV para láminas planas
Aquí las cargas están expresadas en coordenadas locales:
22
PTV para láminas planas
Aquí las cargas están expresadas en coordenadas locales.
De donde se deduce que el trabajo de deformación virtual puede obtenerse
como suma directa de las contribuciones de membrana, flexión y cortante.
23
PTV para láminas planas
Observe que todas las derivadas que aparecen en los integrandos son
de primer grado, lo que permite utilizar elementos finitos de clase C0
24
Formulación de EFs de cascarón
planos de Reissner-Mindlin
Considere un elemento finito isoparamétrico de clase C0 de n nodos.
El campo de movimientos en coordenadas locales se discretiza así:
Vector de movimientos locales del nodo i
25
Convenio
de signos
Este es el mismo
convenio de signos
para los giros que
se empleó en el
capítulo de placas.
26
Discretización del campo de
deformaciones generalizadas
Recuerde que:
(x',y',z')
Vector de deformaciones
generalizadas de membrana (alargamiento)
Vector de deformaciones
generalizadas de flexión (curvaturas)
Vector de deformaciones
generalizadas de cortante
27
Vector de
deformaciones
generalizadas
locales del
elemento
Matriz de deformaciones
generalizadas locales del
nodo i
Matriz de deformaciones
generalizadas locales 28
del elemento
Matriz de deformaciones
generalizadas locales del
nodo i
Matriz de deformaciones
generalizadas locales de
membrana del nodo i
Matriz de deformaciones
generalizadas locales de
flexión del nodo i
Matriz de deformaciones
generalizadas locales de
cortante del nodo i
29
Obtención de la matriz de rigidez y del
vector de fuerzas nodales equivalentes local
Operando de la forma usual se obtiene que:
donde la matriz de rigidez y el vector de fuerzas nodales
equivalentes expresados en coordenadas locales está dado por:
30
Matriz de rigidez local
Matriz de rigidez
de membrana
Matriz de rigidez
de flexión
Matriz de rigidez
de cortante
MdR de acoplamiento flexiónmembrana
31
32
Es posible de mostrar que la matriz de rigidez local del elemento
puede escribirse como:
donde K'PS y K'PB son las matrices de rigidez correspondientes al
problema de "tensión plana" y de “análisis de losas por ReissnerMindlin” respectivamente; se asume que ambas matrices K'PS y K'PB se
calcularon para el mismo número de nodos y la misma tipología que el
elemento de lámina plana utilizado.
Por consiguiente, si no existe acoplamiento membrana-flexión, la matriz
de rigidez local de un elemento de lámina plana puede obtenerse
directamente ampliando la matriz de rigidez para el caso de flexión de
33
placas con la del elemento de tensión plana correspondiente.
Se podría decir que a nivel local los esfuerzos de
membrana equilibran las acciones contenidas en
el plano de elemento, mientras que las acciones
normales a dicho plano provocan un estado de
flexión independiente, pudiendo obtenerse,
siempre a nivel local, los movimientos,
deformaciones y esfuerzos de ambos estados de
manera totalmente desacoplada. El acoplamiento
entre los estados de membrana y flexión se
produce al ensamblar en ejes globales la matriz
(e)
(e)
K y el vector f .
34
Cambio de ejes de coordenadas
5 variables
6 variables: incluye θz
Convención de signos para las
rotaciones locales y globales
35
Recuerdo de cambio de base
36
Se incluyen θz y Mz
Matriz de transformación
de desplazamientos nodales
Matriz de transformación
de giros
Y como el elemento es plano:
37
Cálculo de los
cosenos directores
locales
Producto cruz
38
39
Cálculo de los cosenos directores
locales
40
Matriz de rigidez y vector de fuerzas
nodales equivalentes en coord. globales
donde,
El uso de la matriz B reduce
significativamente el número de cálculos
41
(e)
Cálculo de K y f
(e)
Primero que todo se debe definir la coordenadas
de los nodos del elemento con respecto a los ejes
locales x' y' z'. Suponiendo que los orígenes de
los sistemas local y global coinciden:
Basta esta transformación por que K(e) es
independiente del origen de coordenadas del
elemento.
A partir de aquí se procede como un elemento
isoparamétrico 2D.
42
(e)
Cálculo de K y f
(e)
43
Cálculo de fuerzas de membrana,
momentos flectores y fuerzas cortantes
en el elemento
Estas se realizan en coordenadas locales al
elemento, por lo que:
Tenga en cuenta que se pueden utilizar las
matrices:
44
Elementos de cascarón plano de
RM más usuales
Un elemento de cascarón puede considerarse como una
simple superposición de un elemento de tensión plana y otro
de flexión de placas, que en el caso más usual (D'mf = 0)
contribuyen de forma totalmente desacoplada a la matriz de
rigidez local del elemento. Desde este punto de vista, se
puede afirmar que cualquiera de los elementos de tensión
plana y de placa de RM podrían combinarse para formar un
elemento de cascarón plano.
Como regla general, es conveniente seleccionar elementos
de la misma familia y con el mismo número de nodos. Así
mismo el elemento de placa no puede tener bloqueo de la
solución por efecto del cortante y ni puede tener mecanismos
propagables.
45
Bloqueo de la solución
46
47
48
Métodos para evitar el bloqueo de la
solución
1.Métodos de integración reducida y selectiva:
son métodos que subintegran la matriz K'c +
K'm.
2.Métodos que utilizan campos de deformación
por cortante impuestos.
49
Integración con cuadraturas de GaussLegendre y singularidad de la matriz K
50
Singularidad de la matriz de rigidez
Cuando K es singular se tiene que j-kp>0. Esta
es una condición necesaria pero no suficiente.
Si j-kp>0, muy probablemente K es singular
Si j-kp≤0, K es invertible
El criterio j-kp>0 es aplicable a cualquier tipo de
elemento finito y también es aplicable a la
estructura en su totalidad. Es aplicable
individualmente a la matriz K, a la matriz Kr o
donde r = m, f, c.
51
Puntos de integración de Gauss-Legendre
Ejemplo:
subintegrando Kf
Num #gld/nodo
nodos
En este caso en
particular se debe
usar la estrategia de
integración c, para
subintegrar la matriz
Kf
#gdl
restringidos
El criterio j-pk>0
es aplicable a
cualquier tipo de
elemento finito y
también
es
aplicable a la
estructura en su
totalidad.
probablemente es
52
29 nodos
j = 29x2 – 3 = 55 gdl libres
k = 3 componentes deformación (ex, ey, gxy)
p = 6 (puntos de integración)
j – kp = 55 – 3x6 = 27 > 0 (Kdd probablemente
es singular)
29 nodos
j = 29x2 – 3 = 55 gdl libres
k = 3 componentes deformación (ex, ey, gxy)
p = 24 (puntos de integración)
53
j – kp = 55 – 3x24 = -17 < 0 (Kdd es invertible)
Métodos de integración
reducida/selectiva
54
Prueba de
singularidad
Borde
apoyado
FALTA ENTENDER DE DONDE SALEN ESTOS NUMEROS
Borde
empotrado
Se deduce que el EF de 8 nodos no satisface la prueba de singularidad
para ciertas configuraciones de malla, lo que pone en entredicho su
buen comportamiento como elemento de lámina delgada. De otro lado
55
los EFs de 4 y 9 nodos siempre satisfacen la prueba de singularidad.
Teorías de láminas planas de
Kirchhoff
●
●
●
Requiere elementos de continuidad C0 para discretizar los
movimientos de membrana (u', v') y elementos de continuidad
C1 para discretizar los movimientos de flexión. Esto obliga a
formulaciones más complicadas.
Presentan incompatibilidades de desplazamientos a lo largo
de lados comunes de elementos no coplanares. Dicha
incompatiblidad se traduce en una rigidización de la estructura
y sólo puede evitarse utilizando mallas más tupidas o
utilizando elementos de tensión plana especiales en lo que el
campo de desplazamientos tenga una variación polinómica
del mismo grado que las flechas en el EF de placa.
Las dificultades anteriores no se presentan en elementos de
láminas planas de Reissner-Mindlin, por lo que es preferible la
segunda teoría.
56
Tratamiento de nodos coplanares:
evitando la singularidad de K
●
Ver Sección 10.8
57
Problemas de cuasi-coplanaridad
●
Ver Sección 10.8.5
58
Elementos de lámina plana
rebajada
59