Placas Mindlin PDF

Mtodo de los Elementos Finitos
para Anlisis Estructural
Flexin de placas planas

Teora de Mindlin-Reissner
Estado de deformacin
Q
Hiptesis:
Secciones rectas perpendiculares al plano medio se
mantienen rectas pero no perpendiculares a dicho plano
medio.
w
x
Plano medio:
Punto P:
z y
u=
v = z x
w w
v
u
y
x
X
Flexin de placas - Teora de Mindlin Reissner
Deformaciones unitarias (I)

derivada de los giros = curvatura
y
u
x =
=z
x
x
x
y
x
v
y =
= z
y
y
xy
xy
u v
=
+
=z
z x
y x
y
x
Deformaciones unitarias de flexin

lineales en z
2
Deformaciones unitarias (II)

v w
w
yz =
+
= x +
z
y
y
zx =
w u
w
+
=
+ y
x
z
x
zx
x
yz
Existe distorsin de cortante:

diferencia entre giro y derivada
xz
yz
Y
X
Deformaciones unitarias de cortadura

constantes en el espesor
3
Agrupando las deformaciones unitarias

b curvatura = derivada de los giros
Flexin
= y = z
=
b
=
0
z
z
y
y
xy
x
y
x
x

x w

0 x = z 1 w

y
y
x
Cortadura
y
x
yz
y
= =
=
w

xz
+ y
1 0
x = 2 w
0 1
xy
x
xz
yz
Tensiones
Q
Flexin: tensin plana, variable linealmente con z
1
0
x 0x
= y =
1
0
y
2
0y
xy 0xy
xy
0 0

Cortadura vertical constante
G zx
zx
G
zy
zy
= D( 0 )
xy
x
=G
5
xz
yz
Relacin fuerza-deformacin
M=
zdz =
D( 0 )zdz =Db z 2dz D 0zdz
h3
M=
Db + M0
12
=z b
y
x
M
M 0x
y x
M = M y = D
+ M 0y
x
M
M xy
D(1 ) y x 0xy
x
2
+h / 2
Cortantes
Q=
dz = h G
h /2
6
Rigidez de placa
Eh 3
D=
12(1 2 )
Ecuaciones de equilibrio
Idem. a la teora clsica
M x
x
M xy
y
= Qxz
M xy
x
My
y
= Qyz
Qzx Qzy
+
+ qz = 0
x
y
Sustituyendo M y Q en funcin de las 3 deformaciones
x 1 y x G t w
y
+ y
+
=
x x
y
x
D x
2 y y
x
y 1 y x G t w
+
x
=
+
y y
x
x
D y
2 x y
qz
w
w
+ y +
x +
=0
y y
x x
Gt
7
3 ecuaciones
Orden m = 2
Energa elstica
Densidad
U0 =
T d + T d =
0
U0 =
Energa
U =
U =
U =
1
2
1
2
1
2
1
2
T D d T0 D d + G T d
0
T D T0 D + G T
2
T D dv T D 0dv + 21 G T dv
bT D b z 2dz dA z bT D 0dz dA + 21 G T dz dA
h3
12
b D b dA + b M0 dA +
T
1
2
T
G
k
t
dA
Orden n = 2 funciones C0
8
Elementos finitos para placas Mindlin

Q
Interpolacin C0 para las tres incgnitas

)
w = N i Wi
w

x =

y

x = N i xi
y = N i yi
w3
x3
w4
N
0 ... ... N n
1 0
0 N 1 0 ... ... 0
0
0 N 1 ... ... 0
x2
Y w1
x
Y4
Nn
Y2
3
w = N e
w2
Y3
x4
W1
x1
0 y 1
...
Wn
N n
xn
yn
Elemento normalizado de lado 2 en XY.
Y1
x1
Formulacin isoparamtrica
Q
Interpolacin de coordenadas:
x
y =

Q
Q
N 1 0 N 2
0 N
0
1
0
N2
x1
1
...
x 2
...
y2

...
x = N xe
Elementos con lados curvos en el plano XY

Cualquier nmero de nudos por lado
Z
Y
3
10
Deformaciones unitarias de flexin

Q
Matriz BF (Deformaciones de flexin)
=
=
z
y
=zb
xy
y
BF = BF 1
11
BF 2 ... BFn
= z 1 w = z 1 N e = z BF e
0
0
N i
BFi = 0
N i
0
N i
x
N i
Deformaciones unitarias de cortadura

Q
Matriz BC (Deformaciones de cortadura)

w
yz y
= =
=
xz w +
y
x

1 0 w

x = 2 w

0 1
y
= 2 w = 2 N e = BC e
BC = BC 1
12
BC 2 ... BCn
N i
y
BCi =
N i
N i
0
Ni
Ecuacin de equilibrio
Potencial total
= U +V =
1
2
T D dv T D 0dv + 21 G T dv wT qs ds
Flexin
Cortante
Hiptesis MEF
= 21 eT
z 2 BTF D BF dv e eT
z BTF D 0 dv +
+ 21 eT G BTC BC dv e eT NT qs ds
Estacionario respecto de las deformaciones nodales e
z 2 BTF D BF dv e + G BTC BC dv e = z BTF D 0 dv + NT qs ds

K flexin
K cortante
13
KF + KC e = PT + Ps
Matriz de rigidez (I)

Q
Flexin
KF =
z 2 BTF D BF dv
Integrando en z
KF =
Q
Cortante
h3 T
BF D BF dx dy =
12
KC =
+1
h3 T
12 BF D BF J d d
1
T
G
B
C BC dv
Teor.
+1
KC =
h k G BTC BC J d d
Real
Factor de correccin
de cortante (0.8)
14
Matriz de rigidez (II)
+1
+1
h T
KF =
BF D BF J d d
12
1
KC =
h k G BTC BC J d d
(0,0)
(-1,-1)
BF BC contienen las derivadas de N respecto de x,y

Relacin entre las derivadas:
N i

=
N i
N i
y N i

x
x
= J
N i
y N i
Necesario para BF
N i
N i

x
1
= J
N i
N i
Jacobiana de la transformacin x,y / ,

15
(+1,+1)
Matriz de rigidez (III)

Matriz Jacobiana J
x

J=
x
Conocidos de N(,)
Ni
y
xi

=
y
Ni
xi

Ni
yi
Ni
yi
Coordenadas de los nudos
Espesor variable interpolado
16
h = N i hi
Fuerzas de superficie
q qz
z
qs = q x = 0

q y 0

Slo fuerza vertical

en cada nudo
qs = N qes
Variacin permitida
Ps =
Ms =
17
qz
NT qs ds =
T
N1

...

NT
n
N
1
NT N ds qes = Ms qes
11
1n
Ms ... Ms
... Nn dx dy = ... ... ...
Mn 1 ... Mnn
s
s
Fuerzas de superficie (cont.)
s1

... =

Psn
11
1n
M
...
M
q
s

... ... ... ...

Mn 1 ... Mnn q n
s
s

N N
j
+1 i
Msij = 0
1
0
18
0
Ni N j
q q
iz iz
q i = qix = 0

qiy 0

0 J d d
N i N j
qz
Ps
Fuerzas de lnea sobre los lados
q Lz
qL =
q
L x
qLy
qLz
qLY
q L = N qeL
Variacin permitida
PL =
ML =
19
qLx
NT N dl =
NT N dl qeL = ML qeL
T
N1

...

NT
n
N
1
11
ML ... M1Ln
... Nn dl = ... ... ...
Mn 1 ... Mnn
L
L
Fuerzas trmicas
PT = z BTF D 0 dv
Temperatura lineal en el canto T=Tm+zTg

PT =
Tm + zTg
z BTF D
Tm + zTg dv =
1 +h / 2
BTF D 1 (z Tm + z 2Tg )dz
h / 2
0

1
h Tg T
BF D 1 dA =

12
0

3
PT =
1
E h Tg T
BF 1 dA
12(1 )

0
3
Slo el gradiente produce

momentos nodales
20
dA
Fuerzas trmicas (cont.)
PT =
21
N
1
y
N
1
E h Tg
... dA
12(1 )
0
N n
y
N n
M0x
Z
M0y
Y
M0y
M0x
Slo se producen momentos nodales,

pero no fuerza z
Bloqueo por cortante

Q
Q
Q
Rigidez total: K = KF + KC
Elementos de pequeo espesor: rigidez Kc sobre estimada.
Relacin entre KF y KC
=
Q
Q
Q
Rig. cortante
Rig. flexin
Gk
E h2
Si h 0, . Rigidez a cortante sobrestimada

Elementos con pocos nudos (4, 6).
Soluciones:
X Teora de Kirchhoff discreta (DKT). Imponer Ucort = 0 en una
serie de puntos. Otros elementos.
X Integracin reducida de KC Fcil de implementar.
22

Placas Mindlin PDF

Cargado por

Copyright:

Formatos disponibles

Placas Mindlin PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Placas Mindlin PDF

Cargado por

Copyright:

Formatos disponibles

Mtodo de los Elementos Finitos

para Anlisis Estructural

Flexin de placas planas

Flexin de placas - Teora de Mindlin Reissner

Deformaciones unitarias (I)

Deformaciones unitarias de flexin

Flexin de placas - Teora de Mindlin Reissner

Deformaciones unitarias (II)

Existe distorsin de cortante:

Deformaciones unitarias de cortadura

Flexin de placas - Teora de Mindlin Reissner

Agrupando las deformaciones unitarias

Flexin de placas - Teora de Mindlin Reissner

Flexin: tensin plana, variable linealmente con z

Cortadura vertical constante

Flexin de placas - Teora de Mindlin Reissner

D( 0 )zdz =Db z 2dz D 0zdz

Flexin de placas - Teora de Mindlin Reissner

Sustituyendo M y Q en funcin de las 3 deformaciones

Flexin de placas - Teora de Mindlin Reissner

Flexin de placas - Teora de Mindlin Reissner

Elementos finitos para placas Mindlin

Interpolacin C0 para las tres incgnitas

Elemento normalizado de lado 2 en XY.

Elementos con lados curvos en el plano XY

Flexin de placas - Teora de Mindlin Reissner

Deformaciones unitarias de flexin

Matriz BF (Deformaciones de flexin)

Flexin de placas - Teora de Mindlin Reissner

Deformaciones unitarias de cortadura

Matriz BC (Deformaciones de cortadura)

Flexin de placas - Teora de Mindlin Reissner

Estacionario respecto de las deformaciones nodales e

z 2 BTF D BF dv e + G BTC BC dv e = z BTF D 0 dv + NT qs ds

Matriz de rigidez (I)

Flexin de placas - Teora de Mindlin Reissner

Matriz de rigidez (II)

BF BC contienen las derivadas de N respecto de x,y

Jacobiana de la transformacin x,y / ,

Flexin de placas - Teora de Mindlin Reissner

Matriz de rigidez (III)

Coordenadas de los nudos

Espesor variable interpolado

Flexin de placas - Teora de Mindlin Reissner

Slo fuerza vertical

... Nn dx dy = ... ... ...

Flexin de placas - Teora de Mindlin Reissner

Fuerzas de superficie (cont.)

Flexin de placas - Teora de Mindlin Reissner

Fuerzas de lnea sobre los lados

... Nn dl = ... ... ...

Flexin de placas - Teora de Mindlin Reissner

Temperatura lineal en el canto T=Tm+zTg

Slo el gradiente produce

Flexin de placas - Teora de Mindlin Reissner

Fuerzas trmicas (cont.)

Slo se producen momentos nodales,

Bloqueo por cortante

Si h 0, . Rigidez a cortante sobrestimada