1

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

Diego Andrés Alvarez MarínProfesor Asociado

Universidad Nacional de ColombiaSede Manizales

2

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)

3

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)

4

5

6

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

7

Desplazamientos

8Desplazamientos 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.

Desplazamientos

9

Campo vectorial de movimientos en coordenadas locales

10

Campo vectorial de deformaciones (en coordenadas locales)

11

Vector de deformacionesgeneralizadas de membrana (alargamiento)

Vector de deformacionesgeneralizadas de flexión (curvaturas)

Vector de deformacionesgeneralizadas de cortante

(x',y',z')

12

Vector de deformaciones generalizadas (en coordenadas locales)

Tenga en cuenta que:

donde

13

Campo de esfuerzos

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 a efectos de cortantetransversal

Nota: esta no es la representación verdadera de los esfuerzos (ver más adelante)

14

Ley de Hooke (relación entre esfuerzos y deformaciones)

15

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

16

Ley de Hooke (relación entre esfuerzos y deformaciones)

17

Distribución de

esfuerzos en el

espesor del cascarón

18

Vector de esfuerzos generalizados locales

x

y

19

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.

20

Relaciones entre esfuerzos y deformaciones generalizadas locales

MCG de acoplamiento membrana-flexión

21

Si las propiedades del material no varían con el espesor, tenemos que:

Ya que:

Se anula para un material homogéneo

MCG de acoplamiento membrana-flexión

22

PTV para láminas planas

Aquí las cargas están expresadas en coordenadas locales:

23

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.

24

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

25

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

26

Conveniode signosEste es el mismo convenio de signos para los giros que se empleó en el capítulo de placas.

27

Discretización del campo de deformaciones generalizadas

Vector de deformacionesgeneralizadas de membrana (alargamiento)

Vector de deformacionesgeneralizadas de flexión (curvaturas)

Vector de deformacionesgeneralizadas de cortante

(x',y',z')

Recuerde que:

28Matriz de deformaciones generalizadas locales delnodo i

Matriz de deformaciones generalizadas localesdel elemento

Vector de deformaciones generalizadas locales del elemento

29

Matriz de deformaciones generalizadas locales delnodo i

Matriz de deformaciones generalizadas locales demembrana del nodo i

Matriz de deformaciones generalizadas locales deflexión del nodo i

Matriz de deformaciones generalizadas locales decortante del nodo i

30

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:

31

Matriz de rigidezde membrana

Matriz de rigidezde flexión

Matriz de rigidezde cortante

MdR de acopla-miento flexión-membrana

Matriz de rigidez local

32

33

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 Reissner-Mindlin” 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 placas con la del elemento de tensión plana correspondiente.

34

Cambio de ejes de coordenadas

Convención de signos para las rotaciones locales y globales

35

Y como el elemento es plano:

36

Cálculo de los cosenos directores locales

Producto cruz

37

Cálculo de los cosenos directores locales

38

39donde,

Ensamblaje de las ecuaciones de rigidez

El uso de la matriz B reduce significativamente el número de cálculos

40

Cálculo de K(e) 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.

41

Cálculo de la matriz K y el vector f

42

Cálculo de fuerzas de membrana, momentos flectores y fuerzas cortantes

en el elementoEstas se realizan en coordenadas locales al elemento, por lo que:

donde,

43

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.