7

Flexión de placas planas

7.1. INTRODUCCIÓN
Las placas son estructuras continuas, formadas por un dominio material plano, de
espesor constante o variable h mucho menor que las dimensiones transversales del
dominio, y sobre el que actúan unas cargas exteriores que pueden ser fuerzas
perpendiculares al plano de la placa, o bien momentos contenidos en dicho plano
(figura 7.1). Estas cargas son absorbidas por el sólido por medio de un efecto de
flexión y deformación transversal. En sentido estricto, no se consideran en la flexión
de placas las fuerzas contenidas en el plano de la placa, las cuales harían trabajar al
material en su plano, y por lo tanto su tratamiento corresponde a un problema de
elasticidad en dos dimensiones, ya estudiado.

Figura 7.1 Placa plana a flexión

En los problemas de elasticidad estudiados hasta ahora no se ha introducido ninguna

simplificación relativa al estado de deformaciones unitarias ε del sólido: éste es
 Flexión de placas planas 121

introducido de forma exacta en la formulación, y únicamente se introduce la

hipótesis de discretización del campo de desplazamientos, inherente al MEF. Sin
embargo en la teoría de flexión de placas se introducen ciertas hipótesis para
simplificar el problema, y reducirlo a dos dimensiones. Estas simplificaciones son:
• Las tensiones en la dirección perpendicular a la placa son nulas, es decir σz=0.
Esto sólo ocurre si el espesor es muy pequeño frente a las dimensiones
transversales de la placa.
• Las secciones rectas perpendiculares al plano medio de la placa permanecen
rectas durante la deformación, aunque no siempre se supone que siguen siendo
perpendiculares al plano medio. Más adelante se estudiarán ambas situaciones.
En el planteamiento por el método de los elementos finitos, estas dos hipótesis se
siguen manteniendo, por lo que las soluciones que se obtengan serán exactas en la
medida que se cumplan dichas hipótesis.

7.2. TEORÍA CLÁSICA DE FLEXIÓN DE PLACAS

La teoría clásica de flexión de placas se basa en la suposición de que las secciones
rectas perpendiculares al plano medio de la placa permanecen rectas y
perpendiculares a dicho plano medio en el estado deformado de la placa. La teoría
completa de la flexión de placas usando esta hipótesis fue enunciada por G. R.
Kirchhoff en 1850, aunque deben mencionarse asimismo otros trabajos previos como
los de Sophie Germain (1816), Navier (1823), Lagrange (1828) o Poisson (1829).
Como consecuencia de la hipótesis de deformación efectuada, el giro que sufre una
sección recta perpendicular al plano medio de la placa es igual a la pendiente de
dicho plano medio. Con el sistema de ejes adoptado, esto implica las siguientes
relaciones entre los dos giros de la sección recta y las derivadas parciales de la flecha
w (figura 7.2):
∂w ∂w
θx = θy = − (7.1)
∂y dx
En consecuencia, la deformación de la placa queda perfectamente definida
conociendo únicamente el valor de la deformación w en dirección z transversal a ella,
que es una función w(x,y).
Como se verá más adelante, esta suposición respecto a la deformación equivale a
despreciar la deformación de cortadura que existe transversalmente en la placa y por
lo tanto a despreciar la energía debida al esfuerzo cortante. Este planteamiento es
únicamente válido cuando el espesor de la placa es muy pequeño, de tal manera que
realmente pueda despreciarse la energía debida a las deformaciones de cortadura.
122 Flexión de placas planas

Este planteamiento se corresponde con la teoría clásica de flexión de vigas, pues de

hecho la hipótesis de deformación efectuada se corresponde con la hipótesis de
Navier.

w
Y

Figura 7.2 Deformaciones en el plano medio de una placa

7.2.1 Estado de deformación

Un punto situado en el plano medio de la placa (z=0), sólo tiene una deformación
transversal w en dirección z. Sin embargo, los puntos que están fuera del plano
medio sufren además deformaciones laterales u, v en las direcciones x, y, que son
debidas al giro de la normal a la placa (figura 7.3).

Figura 7.3 Deformaciones interiores en una placa

 Flexión de placas planas 123

Para un punto cualquiera P, situado a una distancia z del plano medio las
deformaciones son:
u = z θy v = −z θx (7.2)

En base a la hipótesis de Kirchhoff los giros de la sección recta son iguales a las
derivadas de la flecha w, respecto a x e y. Por lo tanto las deformaciones u, v del
punto P son:
∂w ∂w
u = −z v = −z (7.3)
dx ∂y

7.2.2 Deformaciones unitarias

Las deformaciones unitarias en el punto P son:
∂u ∂ 2w
εx = = −z
∂x ∂x 2
∂v ∂ 2w
εy = = −z
∂y ∂y 2
∂w
εz = =0 (7.4)
∂z
∂u ∂v ∂ 2w
γxy = + = −2z
∂y ∂ x ∂ x ∂y
∂v ∂w ∂w ∂w
γyz = + =− + =0
∂z ∂y ∂y ∂y
∂w ∂u ∂w ∂w
γzx = + =− + =0
∂x ∂z ∂x ∂x
Se observa que no aparecen distorsiones de cortante en dirección vertical, como
consecuencia de la hipótesis de deformación que se ha efectuado de que los giros
coinciden con las derivadas de la flecha. Las tres únicas deformaciones unitarias no
nulas tienen una distribución lineal en el espesor de la placa, y se pueden poner
agrupadas en forma de vector:
⎧
⎪ ∂ 2w ⎫ ⎪
⎪
⎪ − ⎪
⎪
⎧ ⎫ ⎪
⎪ ∂ x 2 ⎪
⎪
⎪
⎪ εx ⎪
⎪ ⎪ ⎪
⎪ ⎪ ⎪
⎪ ∂w ⎪
2
⎪
⎪ ⎪
ε = ⎨ εy ⎬ = z ⎨ − 2 ⎬ = z b (7.5)
⎪
⎪ ⎪ ⎪ ∂y ⎪
⎪γxy ⎪
⎪ ⎪
⎪
⎪
⎪
⎪
⎪ ⎭
⎩ ⎪ ⎪ ∂w ⎪
2
⎪
⎪
⎪−2 ⎪
⎪
⎪
⎪
⎩ ∂ x ∂ y ⎪
⎪
⎭
124 Flexión de placas planas

En esta expresión se ha definido el vector de curvaturas b, habitual en la teoría de

placas, que contiene las tres derivadas parciales segundas de la superficie deformada
de la placa. Este vector de curvaturas no depende de z, sino sólo de la posición en la
superficie media de la placa. Para los desarrollos posteriores conviene poner la
expresión anterior en la forma siguiente:
⎡ 2 ⎤
⎢ − ∂ ⎥
⎧ε ⎪ ⎫ ⎢ ∂ x 2 ⎥⎥
⎪
⎪ ⎪ ⎢
⎪ x
⎪ ⎢ ∂ 2 ⎥⎥
⎪ ⎪ ⎢
⎨ εy ⎬ = z ⎢ − 2 ⎥ w = z ∂ w (7.6)
⎪
⎪ ⎪ ∂y ⎥
⎪γxy ⎪
⎪ ⎢
⎢ 2 ⎥
⎪ ⎭
⎩ ⎪
⎢−2 ∂ ⎥
⎢⎢ ∂ x ∂y ⎥⎥
⎣ ⎦
donde se identifica al operador ∂ que en este caso es de orden 2.

Y Y

XY
Z
X
Y

Figura 7.4 Deformaciones unitarias en una placa. Teoría clásica

7.2.3 Estado de tensiones

El estado de tensiones en la placa tiene cinco componentes: σx σy τxy τyz τzx . Para
facilitar el estudio posterior, estas tensiones se agrupan en dos vectores diferentes
(figura 7.5): el vector σ contiene las tres tensiones contenidas en el plano de la placa,
y el vector τ contiene las dos tensiones cortantes perpendiculares a ella.
⎧
⎪ σx ⎫⎪
⎪ ⎪
⎪ ⎪ ⎪τyz ⎫
⎧ ⎪
σ = ⎨ σy ⎪⎬
⎪ τ = ⎪⎨τ ⎪⎬ (7.7)
⎪
⎪ ⎪
⎪ ⎪
⎪ zx ⎪
⎩ ⎪
⎭
⎪
⎪τ ⎪
⎩ xy ⎪⎭
 Flexión de placas planas 125

Y YZ

XY XZ
Z X

Figura 7.5 Estado de tensiones en una placa

Dado que σz=0 y que sólo existen las tres deformaciones unitarias situadas en el
plano XY, la relación entre las tensiones y dichas deformaciones unitarias es:
⎡ ⎤
⎢ ⎥
⎪⎧⎪ σx ⎪⎫⎪ ⎢1 ν 0 ⎥ ⎛⎜⎪⎧⎪ εx ⎪⎫⎪ ⎪⎧⎪ ε ⎪⎫⎪⎞⎟
⎪⎪ ⎪⎪ E ⎢ ⎥ ⎜⎜⎪⎪ ⎪⎪ ⎪⎪ 0x ⎪⎪⎟ ⎟⎟
σ = ⎨ σy ⎬ = ⎢ ν 1 ⎥ ⎜
0 ⎜⎨ εy ⎬ − ⎨ ε0y ⎬⎟⎟ ⎪ (7.8)
⎪⎪ ⎪⎪ 1 − ν 2 ⎢ ⎥ ⎜⎜⎪⎪ ⎪⎪ ⎪⎪ ⎪⎟
⎪τxy ⎪ ⎢ 1 − ν ⎥ ⎜⎪γxy ⎪ ⎪γ 0xy ⎪⎪⎟⎟
⎩⎪ ⎭⎪ ⎢0 0 ⎥ ⎝⎩⎪ ⎭⎪ ⎩⎪ ⎭⎪⎠⎟
⎢⎣ 2 ⎦ ⎥

σ = D(ε − ε0 ) (7.9)
Siendo D la matriz elástica y ε0 el vector de deformaciones unitarias iniciales.
Cuando está originado por una variación de temperatura T este vector tiene el valor
siguiente:
⎧
⎪αT ⎫
⎪
⎪
⎪ ⎪ ⎪
ε0T = ⎨αT ⎪
⎪ ⎬ (7.10)
⎪
⎪ ⎪
⎪
⎪
⎪ 0 ⎪
⎪ ⎪
⎩ ⎪
⎭
Al ser la variación de las deformaciones unitarias lineal en el espesor, también las
tensiones varían de forma lineal en él. Cada punto del material está sometido a un
estado plano de tensiones σx σy τxy situado en el plano XY, que varía de forma
lineal a lo largo del espesor de la placa.
Las tensiones cortantes verticales τ son proporcionales a las deformaciones unitarias
correspondientes γ, que según la hipótesis de deformación efectuada son ambas
nulas:
⎧
⎪τyz ⎪⎫ ⎪⎧⎪γyz ⎪⎫⎪ ⎧⎪0⎫⎪
⎪ ⎪ = G = G ⎪ ⎪
⎨τ ⎬ ⎨γ ⎬ ⎨ ⎬ (7.11)
⎪ ⎪
⎪ zx ⎭⎪
⎩ ⎪ ⎪
⎩⎪ zx ⎭⎪ ⎪⎪0⎪⎪
⎩ ⎭
126 Flexión de placas planas

Pero en realidad deben existir tensiones cortantes verticales τ, para que pueda existir
esfuerzo cortante vertical. Esta incoherencia se resuelve considerando que el módulo
de elasticidad de cortadura G es infinito, es decir que la placa no se deforma ante la
cortadura vertical, aunque existan tensiones en dicha dirección. Más adelante se
obtendrá una expresión de las tensiones cortantes verticales, basándose en criterios
de equilibrio en lugar de en la ecuación constitutiva.

7.2.4 Esfuerzos internos

En lugar de las tensiones, se suelen utilizar en el estudio de las placas los esfuerzos
internos (fuerzas y momentos) correspondientes a dichas tensiones, que son sus
integrales a lo largo del espesor h de la placa.
Momentos flectores. Son el momento estático de las tensiones respecto al plano
medio de la placa, por unidad de anchura de la misma (figura 7.6).
+h / 2 +h / 2

Mx = ∫ σx z dz My = ∫ σy z dz (7.12)
−h / 2 −h / 2

Momento torsor. Es el momento estático de las tensiones cortantes situadas en el

plano de la placa, respecto al plano medio de la misma, y por unidad de anchura
(figura 7.7).
+h / 2

M xy = ∫ τxy z dz (7.13)
−h / 2

Estos tres momentos se suelen agrupar en un vector de momentos M:

⎧
⎪ Mx ⎫⎪ ⎧ ⎪
⎪ ⎫
⎪
⎪ ⎪
⎪ +h / 2 ⎪ σx ⎪ +h / 2
⎪
⎪ ⎪ ⎪ ⎪
⎨ My ⎪⎬= ∫ ⎪ ⎨ σy ⎪
⎬ z dz M= ∫ σ z dz (7.14)
⎪
⎪ ⎪
⎪ −h/2 ⎪
⎪ ⎪
⎪
⎪M xy ⎪
⎪ ⎪τxy ⎪
⎪ −h/2
⎪
⎩ ⎪
⎭ ⎪ ⎭
⎩ ⎪
Esfuerzos cortantes. Son la resultante de las tensiones cortantes transversales a la
placa (figura 7.7).
⎧Qzx ⎪
⎪ ⎫ +h / 2 ⎧
⎪τzx ⎫
⎪
+h / 2
⎪
⎨ ⎪ ⎬= ∫ ⎪ ⎨τ ⎪
⎪Qzy ⎪ −h/2 ⎩ ⎪
⎬ dz
zy ⎪
⎪ ⎭ ⎪
Q= ∫ τ dz (7.15)
⎪ ⎭
⎩ ⎪ −h/2
 Flexión de placas planas 127

Figura 7.6 Esfuerzos de flexión en una placa

YZ
XY XZ

QYZ
QXZ

MXY
MXY

Figura 7.7 Esfuerzos cortantes y de torsión en una placa

7.2.5 Relación fuerza deformación

Sustituyendo la expresión de σ en la definición de los esfuerzos internos M, y a
continuación la expresión de las deformaciones unitarias en función de las curvaturas
se obtiene:

∫ σzdz = ∫ D(ε − ε )zdz =Db∫ z dz − D∫ ε zdz

2
M= 0 0
(7.16)
128 Flexión de placas planas

La primera integral proporciona el momento de inercia por unidad de anchura de la

placa, y la segunda corresponde a los momentos producidos por las deformaciones
iniciales, que se definen como:

M0 =− D ∫ ε0zdz (7.17)

Por lo tanto la expresión de los momentos puede ponerse de forma más compacta
como:
h3
M= Db + M0 (7.18)
12
Desarrollando esta ecuación se obtiene:
⎧
⎪ ⎛ ∂ 2w ∂ 2w ⎞⎟⎪⎫⎪
⎪ ⎜
⎪−D ⎜⎜ 2 + ν 2 ⎟⎟⎪⎪
⎧⎪ M ⎫⎪ ⎪⎪⎪ ⎝ ∂x ∂y ⎠⎪ ⎧⎪
⎪⎪ ⎪ M 0x ⎫⎪⎪
⎪⎪ x ⎪⎪ ⎪
⎪ ⎪ ⎪ ⎛∂ w2
∂ w ⎞⎟⎪⎪ ⎪⎪⎪
2 ⎪⎪
M = ⎪⎨ M y ⎪⎬ = ⎪⎨−D ⎜⎜ 2 + ν ⎟ ⎬ + ⎨ M ⎪⎬ (7.19)
⎪⎪ ⎪⎪ ⎪⎪ ⎝⎜ ∂y ∂ x 2 ⎠⎟⎪⎪ ⎪⎪ 0y ⎪⎪
⎪⎪M xy ⎪⎪ ⎪⎪ ⎪⎪ ⎪⎪M ⎪⎪
⎩⎪ ⎭⎪ ⎪ ∂ 2
w ⎪ ⎩⎪ 0xy ⎭⎪
⎪⎪ −D(1 − ν ) ⎪⎪
⎪ ∂ x ∂y ⎪⎪⎭
⎪⎩
siendo D la rigidez equivalente de la placa:
Eh 3
D= (7.20)
12(1 − ν 2 )
La ecuación (7.19) relaciona los momentos interiores de flexión y torsión M en la
placa, con sus deformaciones, a través de las curvaturas b: es el equivalente para
placas de la ecuación de la elástica de las vigas con cargas térmicas. El segundo
término de dicha ecuación corresponde a los momentos que se originan en la placa
debido a las curvaturas iniciales producidas por la variación de temperatura.

7.2.6 Deformaciones unitarias de origen térmico

Se supone un campo de temperatura lineal en el espesor de la placa, variando entre
un valor +TS en la cara superior y uno ‐TS en la cara inferior. En un punto cualquiera
situado a una distancia z del plano medio, la temperatura vale:
T = zTg (7.21)

siendo Tg = 2T / h el gradiente de temperatura a lo largo del espesor h. Las

deformaciones unitarias iniciales producidas por este campo de temperaturas son:
 Flexión de placas planas 129

⎧αT z ⎪
⎪ ⎫ ⎧1⎪
⎪ ⎫
⎪
⎪ g ⎪ ⎪ ⎪
⎪ ⎪⎪
ε0 = ⎪⎨αTg z ⎪⎬ = αTg z ⎪⎨1⎪⎬ (7.22)
⎪
⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ 0 ⎪ ⎪ ⎪
⎪0⎪
⎪
⎩ ⎪
⎭ ⎪ ⎪
⎩ ⎪
⎭
Los momentos iniciales debidos a esta temperatura son:
⎧
⎪1⎫
⎪
⎪
⎪ ⎪
⎪ ⎪
3
h ⎪
M0 = −D∫ ε0zdz = −α Tg D ⎨1⎬ (7.23)
12 ⎪⎪ ⎪
⎪
⎪
⎪0⎪
⎪ ⎪
⎩ ⎪
⎭
Cuyas componentes son:
⎪⎧⎪ M ⎪ ⎫ ⎧
⎪1⎫⎪
⎪⎪ 0x ⎪ ⎪
⎪ E h 3
α T
⎪
⎪ ⎪
⎪
⎪
M 0 = ⎨ M 0y ⎬ = − ⎪ g ⎪ ⎪
⎨1⎬ (7.24)
⎪⎪ ⎪
⎪ 12(1 − ν ) ⎪
⎪ ⎪
⎪
⎪⎪M ⎪ ⎪ ⎪
⎪0⎪
⎩⎪ 0 xy ⎪
⎭ ⎪ ⎪⎭
⎩ ⎪

7.2.7 Ecuaciones de equilibrio

Para obtener las ecuaciones de equilibrio de una placa se considera el equilibrio
estático de un trozo diferencial de placa sobre el que actúan los esfuerzos internos
de flexión y cortadura y la carga exterior en dirección z (figura 7.8).

M xy Qxz
Qyz
M xy Mx ∂Qxz Qyz +
∂Qyz
dy
My My +
∂M y
dy Qxz + dx ∂y
∂y ∂x

∂M xy
∂M x M xy + dy
Mx + dx ∂y
∂x

∂M xy
M xy + dx
∂x

Figura 7.8. Equilibrio de un elemento diferencial de placa

130 Flexión de placas planas

• El equilibrio de momentos respecto del eje X implica que:

∂M y ∂ M xy
Qyzdxdy + M ydx − M ydx − dydx + M xydy − M xydy − dxdy = 0 (7.25)
∂y ∂x
Simplificando se obtiene:
∂M y ∂ M xy
+ = Qyz (7.26)
∂y ∂x
• El equilibrio de momentos respecto del eje Y implica que :
∂M x ∂M xy
−Qxzdydx + dxdy + dydx = 0 (7.27)
∂x ∂y

Simplificando se obtiene:
∂M x ∂M xy
+ = Qxz (7.28)
∂x ∂y

• El equilibrio en la dirección z es:

∂Qzx ∂Q
dxdy + zy dydx + qzdxdy = 0 (7.29)
∂x ∂y
Simplificando se obtiene:
∂Qzx ∂Qzy
+ + qz = 0 (7.30)
∂x ∂y
Donde qz es la fuerza distribuida por unidad de superficie de la placa.
Sustituyendo los valores de los esfuerzos cortantes dados por las ecuaciones (7.28) y
(7.26) en la (7.29) se obtiene una ecuación que relaciona sólo a los tres momentos:
∂ 2M x ∂ 2M xy ∂ 2M y
+2 + + qz = 0 (7.31)
∂x 2 ∂ x ∂y ∂y 2
Sustituyendo los valores de los momentos en función de la deformación w dados por
(7.19) (no incluyendo el término debido a la temperatura), se obtiene una única
ecuación:
∂ 4w ∂ 4w ∂ 4w qz
+2 2 2 + 4 = (7.32)
∂x 4 ∂ x ∂y ∂y D
que es la ecuación diferencial que controla el problema de la flexión de placas en la
teoría clásica, en ausencia de temperaturas. Obsérvese la similitud con la ecuación
diferencial de la flexión de una viga, pero extendida a un problema de dos
dimensiones. Se trata de la ecuación bi‐armónica, que puede ponerse en la forma
compacta siguiente empleando el operador laplaciano ∇2 :
 Flexión de placas planas 131

qz
∇2 (∇2w ) = (7.33)
D
La relación entre los esfuerzos cortantes y la deformación puede obtenerse a base de
sustituir la expresión de los momentos (7.19) en las ecuaciones (7.26) y (7.28):
∂ ⎛⎜ ∂ 2w ∂ 2w ⎞⎟ ∂
Qyz = −D ⎜⎜ 2 + 2 ⎟⎟ = −D (∇2w ) (7.34)
∂y ⎝ ∂ x ∂y ⎠ ∂y

∂ ⎛⎜ ∂ 2w ∂ 2w ⎞⎟ ∂
Qxz = −D ⎜⎜ 2 + 2 ⎟⎟ = −D (∇2w )
∂x ⎝ ∂x ∂y ⎠ ∂x

7.2.8 Relación entre tensiones y esfuerzos

Para obtener una relación entre las tensiones y los momentos en la placa, resulta
conveniente despejar de la expresión (7.18) el valor de las curvaturas b en función de
los momentos:
12 −1
b= D (M − M 0 ) (7.35)
h3
Por otra parte la expresión de las tensiones en función de las curvaturas es:
σ = D(ε − ε0 ) = z Db − Dε0

Sustituyendo el valor de las curvaturas se obtiene:

12z 12z
σ= 3
M − 3 M0 − Dε0
h h
Sustituyendo las expresiones detalladas de M0 y de ε0 , se anulan los dos últimos
sumandos debidos a las temperaturas, y se obtiene la expresión:
12z
σ= M (7.36)
h3
z Mx z My z M xy
σx = σy = σxy = (7.37)
(h 3 /12) (h 3 /12) (h 3 /12)
En la que se identifica el momento de inercia de la sección recta de la placa, por
unidad de anchura.
• Para obtener la expresión de las tensiones cortantes resulta ventajoso emplear las
ecuaciones de equilibrio del elemento diferencial expresadas en función de las
tensiones. Así la ecuación de equilibrio en la dirección x es:
∂σx ∂τxy ∂τxz
+ + =0 (7.38)
∂x ∂y ∂z
132 Flexión de placas planas

Empleando las expresiones de las tensiones dadas por (7.37) y la expresión del
esfuerzo cortante dada por (7.28) se obtiene:
∂τxz 12z ⎛ ∂ M x ∂ M xy ⎞⎟ 12z
= − 3 ⎜⎜ + ⎟⎟ = − 3 Qxz
∂z h ⎜⎝ ∂ x ∂y ⎠ h
Integrando esta ecuación en z e imponiendo la condición de que la tensión cortante
en la cara superior (z=h/2) es nula se obtiene:
3h 2 − 12z 2
τxz = Qxz (7.39)
2h 3
Que corresponde a una distribución parabólica, igual que la existente en una sección
rectangular de canto h y anchura unidad.
De la misma forma, la otra tensión cortante es:
3h 2 − 12z 2
τyz = Qyz (7.40)
2h 3

7.2.9 Energía elástica de deformación

La energía elástica acumulada en toda la placa es:

∫ε D ε dv − ∫ εT D ε0dv
T
U = 1
2 (7.41)

• Energía elástica en función de las deformaciones

Sustituyendo las deformaciones unitarias en función de las curvaturas se obtiene:

∫b D b z 2 dz dA − ∫ z bT D ε0 dz dA
T
U = 1
2 (7.42)

En el primer término se puede integrar en la coordenada z, y en el segundo se

identifica el valor de los momentos debidos a las deformaciones unitarias iniciales:
3
1 h
∫ 12 b D b dA + ∫ bT M0 dA
T
U = (7.43)
2
Desarrollando el valor de las curvaturas en el primer sumando se obtiene:
⎛ 2 2 2 2⎞
⎜⎜⎜⎛ ∂ w ⎞⎟ + ⎜⎛ ∂ w ⎞⎟ + 2ν ∂ w ∂ w + 2 (1 − ν ) ⎛⎜ ∂ w ⎞⎟ ⎟⎟ dA + U
2 2 2 2
D
U =
2 ∫ ⎜ ⎟ ⎜ ⎟
⎜⎜⎝⎜⎝ ∂ x 2 ⎠⎟ ⎜⎝ ∂y 2 ⎠⎟ ∂ x 2 ∂y 2
⎟
⎜⎝⎜ ∂ x ∂y ⎠⎟ ⎟⎟⎠⎟ T (7.44)

donde se ha denominado UT al segundo sumando, que corresponde a la energía

debida a las deformaciones unitarias iniciales. Si éstas se deben a una distribución
lineal de temperatura, su valor es:
 Flexión de placas planas 133

⎧1⎪
⎪ ⎫
⎪
⎪ ⎪
⎪ ⎛ E h 3 αTg ⎞⎟
T ⎪ ⎪⎜ ⎟⎟ dA
UT = ∫ b ⎨1⎬ ⎜⎜− (7.45)
⎪
⎪ ⎪⎜
⎪⎝⎜ 12(1 − ν )⎠⎟⎟
⎪
⎪0⎪
⎪ ⎪
⎩ ⎪
⎭
E h 3 α Tg ⎛ ∂ 2w ∂ 2w ⎞
UT = −∫ ⎜⎜ + 2 ⎟⎟⎟dA
⎜
12(1 − ν ) ⎝ ∂ x 2
∂y ⎠

• Energía elástica en función de los esfuerzos internos.

El valor de las curvaturas en función de los momentos es:
12 −1
3 (
b= D M − D−1M0 ) (7.46)
h
Sustituyendo en (7.43) y operando se obtiene:
12
∫ h (M )
1
U = 3
T
D−1M − MT0 D−1M0 dA (7.47)
2

7.2.10 Condiciones de contorno

Para profundizar en el estudio de las condiciones de contorno sin complicar
excesivamente el desarrollo algebraico, se considera una placa rectangular de
dimensiones (a × b), sometida a una carga uniforme. Aplicando el principio del
trabajo virtual se obtienen las ecuaciones de equilibrio y las condiciones de contorno
asociadas a ellas:

∫ q δw dA = δW ∫σ
T
δWEX = IN = δε dv (7.48)
A

Figura 7.9 Placa rectangular

Sustituyendo el valor de las deformaciones unitarias en función de las curvaturas:

∫ q δw dA = ∫ σ
T
δ b z dzdA
A V
134 Flexión de placas planas

Integrado en la variable z el término de la derecha, se identifican los momentos en la

placa:

∫ q δw dA = ∫ M
T
δ b dA
A A

Sustituyendo los valores y agrupando las dos integrales se obtiene:

⎛ 2 2 2 ⎞
⎜⎜M x ∂ δw + M y ∂ δw + 2M xy ∂ δw ⎟⎟ dA + q δw dA = 0
∫ ⎝⎜ ∂x 2 ∂y 2 ∂x ∂y ⎠⎟ ∫
A A

El tercer sumando del integrando de la izquierda se puede descomponer en dos

sumandos:
⎛ 2 2 2 2 ⎞
⎜⎜M x ∂ δw + M y ∂ δw + M xy ∂ δw + M xy ∂ δw ⎟⎟ dA + q δw dA = 0
∫ ⎜⎝ ∂x 2 ∂y 2 ∂ x ∂y ∂y ∂x ⎠⎟ ∫
A A

Integrando por partes la integral de la izquierda se obtiene:

⎛ ∂2M ∂2M y ∂2M xy ⎞⎟ ⎡ ∂δw ⎤
a
⎡ ∂δw ⎤
b
⎜⎜
∫ ⎜⎝ ∂x 2
x
+
∂y 2
+ 2
∂ x ∂y
+ q ⎟⎟ δw dA + ∫
⎠⎟
⎢M x
⎢⎣
⎥ dy + ∫
∂x ⎥⎦ 0
⎢M y
⎢⎣
⎥ dx
∂y ⎥⎦ 0
A A A
a b
⎡⎛ ∂M x ∂M xy ⎞ ⎤ ⎡ ⎤
−∫ ⎢⎜⎜ ⎟⎟ δw ⎥ dy − ⎢⎛⎜⎜ ∂M y + ∂M xy ⎞⎟⎟ δw ⎥ dx
⎢⎝⎜ ∂x
+
∂y ⎠⎟ ⎥⎦ 0 ∫ ⎢⎝⎜ ∂y ∂x ⎠⎟ ⎥⎦ 0
A ⎣ A ⎣
a
⎡ ∂δw ⎤
b
⎡ ∂δw ⎤
+∫ ⎢M xy ⎥ dx + ∫ ⎢M xy ⎥ dy = 0
⎣⎢ ∂x ⎦⎥ 0 ⎢⎣ ∂ y ⎥⎦ 0
A A

En las integrales 4º y 5º se identifican los esfuerzos cortantes:

⎛ ∂2M ∂2M y ∂2M xy ⎞⎟ ⎡ ∂δw ⎤
a
⎡
b
∂δw ⎤
⎜⎜
∫ ⎜⎝ ∂x 2
x
+
∂y 2
+ 2
∂ x ∂y
+ q ⎟⎟ δw dA + ∫
⎠⎟
⎢M x
⎢⎣
⎥ dy + ∫
∂x ⎥⎦ 0
⎢M y
⎢⎣
⎥ dx
∂y ⎥⎦ 0
A A A
a
b ⎡ ∂δw ⎤
b
⎡ ∂δw ⎤
−∫ [Qxz δw ]0 dy − ∫ ⎣⎡Qyz δw ⎦⎤ 0 dx + ∫ ⎢M xy
a
⎥ dx + ∫ ⎢M ⎥ dy = 0
⎢
⎣ ∂ x ⎥
⎦ 0
⎢⎣ xy ∂y ⎦⎥ 0
A A A A

Las dos últimas integrales se pueden volver a integrar por partes:

⎛ ∂2M ∂2M y ∂2M xy ⎞⎟ ⎡ ∂δw ⎤
a
⎡ ∂δw ⎤
b
⎜⎜
∫ ⎜⎝ ∂x 2
x
+
∂y 2
+ 2
∂ x ∂y
+ q ⎟⎟ δw dA + ∫
⎠⎟
⎢M x
⎢⎣
⎥ dy + ∫
∂x ⎥⎦ 0
⎢M y
⎢⎣
⎥ dx
∂y ⎥⎦ 0
A A A
b
b ⎡ ∂M xy ⎤ ⎡ b ⎤a
a
−∫ [Qxz δw ]0 dy − ∫ ⎡Qyz δw ⎤ dx − ⎢
⎣ ⎦0 ∫ ⎢⎣ ∂x ⎥⎦ 0
δw ⎥ dx + ⎢
⎣
(M xy δw ) ⎥
0 ⎦0
A A A
a
⎡ ∂M xy ⎤ a b
−∫ ⎢ δw ⎥ dy + ⎡⎢(M xy δw ) ⎤⎥ = 0
⎢ ∂y ⎦⎥ 0 ⎣ 0 ⎦0
A ⎣

Agrupando:
 Flexión de placas planas 135

⎛ 2 2 2 ⎞ a
⎡
b
∂δw ⎤
⎜⎜ ∂ M x + 2 ∂ M xy + ∂ M y + q ⎟⎟ δw dA + ⎡ ∂δw ⎤
∫ ⎜⎝ ∂x 2 ∂x ∂y ∂y 2 ⎟
⎠⎟
∫ ⎢M x
⎢⎣
⎥ dy + ∫
∂x ⎥⎦ 0
⎢M y
⎢⎣
⎥ dx
∂y ⎥⎦ 0
A A A
a b
⎡⎛ ∂M xy ⎞⎟ ⎤ ⎡⎛ ∂M xy ⎞⎟ ⎤
⎢⎜
−∫ ⎜Qxz + ⎥ ⎜
⎟⎟ δw dy − ∫ ⎢⎜Qyz + ⎟⎟ δw ⎥ dx
⎢⎜⎝ ∂ y ⎠ ⎥ ⎢⎜⎝ ∂x ⎠ ⎥
A ⎣ ⎦0 A ⎣ ⎦0
a
+ ⎡⎢(M xy δw ) ⎤⎥ + ⎡⎢(M xy δw ) ⎤⎥ = 0
b a b

⎣ 0 ⎦0 ⎣ 0 ⎦0

Intercambiando la derivada y la variación en las integrales 2ª y 3ª se obtiene:

⎛ ∂2M ∂ 2M xy ∂2M y ⎞⎟
⎜⎜
∫ ⎜⎜⎜ ∂x 2 + 2 ∂x ∂y + ∂y 2 + q ⎟⎟⎟⎟ δw dA
x

A ⎝ ⎠
a b
⎡ ⎛ ∂w ⎟⎞⎤ ⎡ ⎛ ∂w ⎞⎟⎤
+∫ ⎢⎢M x δ ⎜⎜ ⎟⎟⎥⎥ dy + ∫ ⎢⎢M y δ ⎜⎜ ⎟⎟⎥⎥ dx
⎜⎝ ∂x ⎟⎠ ⎜⎝ ∂y ⎠⎟
A ⎢⎣ ⎥⎦ 0 A ⎢⎣ ⎥⎦ 0 (7.49)
a b
⎡⎛ ∂ M ⎞ ⎤ ⎡⎛ ∂M ⎞⎟ ⎤
⎜ ⎟ ⎜
−∫ ⎢⎢⎜⎜Qxz + xy ⎟ ⎥
⎟⎟ δw ⎥ dy − ∫ ⎢⎜⎜Qyz +
⎢ xy ⎟ ⎥
⎟⎟ δw ⎥ dx
A ⎢⎣⎝
⎜ ∂y ⎟⎠ ⎥ A ⎢⎣⎝
⎜ ∂x ⎠⎟ ⎥
⎦0 ⎦0
a b
⎡⎡ b⎤ ⎡ a⎤
+ ⎢ ⎢M xy δw ⎤⎥ ⎥ + ⎢ ⎡⎢M xy δw ⎤⎥ ⎥ = 0
⎣⎢ ⎣ ⎦ 0 ⎦⎥ 0 ⎣⎢ ⎣ ⎦ 0 ⎦⎥ 0

Esta ecuación debe satisfacerse para cualquier variación arbitraria de la flecha δw , y

por lo tanto deben ser nulos todos los factores que afectan a la δw o a su derivada.
La anulación del integrando de la primera integral proporciona la ecuación de
equilibrio, ya obtenida anteriormente por otro método:
∂2M x ∂2M xy ∂2M y
+2 + +q = 0
∂x 2 ∂x ∂y ∂y 2
Además, la anulación de los integrandos de las restantes integrales proporcionan las
condiciones de contorno que deben cumplirse en los 4 lados. Éstas pueden ser de
distintos tipos:
• Lados paralelos al eje Y (x=0 y x=a):
a. Mx = 0 o bien
⎛ ∂w ⎞ ∂w
b. δ ⎜⎜ ⎟⎟⎟ = 0 → = Cte
⎝ ∂x ⎠ ∂x
∂M xy
c. Qx + =0 o bien
∂y
d. δw = 0 → w = Cte

• Lados paralelos al eje X (y=0 y y=b):

136 Flexión de placas planas

a. My = 0 o bien
⎛ ∂w ⎞ ∂w
b. δ ⎜⎜ ⎟⎟⎟ = 0 → = Cte
⎝⎜ ∂y ⎠ ∂y
∂M xy
c. Qy + =0 o bien
∂x
d. δw = 0 → w = Cte

Las condiciones de apoyo encontradas habitualmente en la práctica se consiguen

mediante combinación de dos de estas condiciones de contorno, una del tipo (a) o (b)
y la otra del tipo (c) o (d). Así un empotramiento implica las condiciones (b) y (d); un
apoyo simple consiste en las condiciones (a) y (d) y un lado libre corresponde a las
condiciones (a) y (c).
Resulta a veces sorprendente la condición de tipo (c) que corresponde en realidad a
un lado libre de tensiones cortantes en él, como es el caso de un lado libre
(añadiendo por otra parte la condición (a) de momento flector nulo).
Una explicación más intuitiva a esta condición fue presentada por Thomson y Tait,
analizando para ello un lado de la placa paralelo al eje X, sobre el que actúa el
momento torsor Mxy. Considerando un elemento diferencial dx, el momento total
actuante sobre él vale Mxy dx. Este momento se puede sustituir de forma
estáticamente equivalente por un par de fuerzas iguales y de sentido contrario de
módulo Mxy, separadas una distancia dx.
z z
dx dx

dx dx
x

M xy

⎛ ⎞
M xydx ⎜⎜M + ∂M xy ⎟⎟dx ∂M xy
M xy +
⎝⎜ ∂x ⎠⎟⎟
xy
∂x

z
dx dx

∂M xy x
qyz* =
∂x

Figura 7.10 Esfuerzos sobre el lado de una placa.

Considerando ahora otro elemento diferencial contiguo al anterior, el valor del

momento en él es M xy + (∂M xy / ∂x )dx , y su resultante sobre todo el elemento vale:
⎡M xy + (∂M xy / ∂x )dx ⎤ dx . Este momento se puede a su vez considerar como formado
⎣ ⎦
por un par de fuerzas iguales de valor M xy + (∂M xy / ∂x )dx separadas la distancia dx.
 Flexión de placas planas 137

En la frontera m‐n entre ambos elementos diferenciales, la fuerza neta resultante es

vertical en la dirección z y de valor M xy + (∂M xy / ∂x )dx − M xy .
Así pues, con este artificio se ha sustituido el momento torsor en el lado Mxy por una
fuerza transversal distribuida sobre dicho lado de valor ∂M xy / ∂x . La fuerza lateral
total se denomina fuerza cortante efectiva Vy, y es la suma de la fuerza equivalente
al torsor, más el esfuerzo cortante en el lado Qy. Esta fuerza efectiva total es la que
debe ser nula si el lado está libre, es decir:
∂M xy
Vy ≡ Qy + =0
∂x
Todo este desarrollo se ha efectuado para una placa rectangular, a fin de simplificar
la formulación, pero de la misma manera puede hacerse para un lado curvo,
obteniéndose expresiones similares. Si se denominan n y s a la normal y tangente al
lado respectivamente, las condiciones resultan ser:
a. Mn = 0 o bien
⎛ ∂w ⎞ ∂w
b. δ ⎜⎜ ⎟⎟⎟ = 0 → = Cte
⎝ ∂n ⎠ ∂n
∂M ns
c. Vn ≡ Qn + =0 o bien
∂s
d. δw = 0 → w = Cte

Figura 7.11 Normal y tangente a un lado curvo

Fuerzas en las esquinas. Los dos últimos términos de la ecuación de los trabajos
virtuales introducen las condiciones en las cuatro esquinas de la placa. Todas estas
condiciones son de la forma: 2M xy δw Para que se satisfagan ante cualquier variación
de w, debe cumplirse bien que w=cte es decir que la esquina esté apoyada o bien
que el momento torsor sea nulo.
En cualquier caso estos términos 2M xy δw representan el trabajo virtual efectuado
por una fuerza de valor 2 Mxy cuando la esquina se mueve un valor δw, con lo que
puede concluirse que realmente existe una fuerza en la esquina de valor Mxy. De
138 Flexión de placas planas

hecho estas fuerzas en las esquinas son la consecuencia de la transformación del

momento torsor del lado en la fuerza efectiva cortante equivalente: en los extremos
del lado las fuerzas debidas a los distintos pares de fuerzas no se cancelan entre sí y
producen la fuerza en la esquina.

Mxy Mxy

Figura 7.12 Fuerzas en las esquinas debidas al momento torsor.

7.3. ELEMENTO PLACA RECTANGULAR DE CUATRO NUDOS

Como ejemplo de un elemento finito basado en la teoría clásica de placas se estudia
a continuación un elemento rectangular de cuatro nudos, que fue uno de los
primeros elementos placa desarrollados (Melosh, 1963, Zienkiewicz, Cheung, 1964).
Se supone que el elemento es rectangular, con lados 2a y 2b (figura 7.13).
2a

Figura 7.13 Elemento placa rectangular

7.3.1 Campo de desplazamientos

Se utilizan como parámetros nodales la flecha y los dos giros en cada uno de sus
cuatro nudos, lo cual da un total de doce grados de libertad para el elemento.
T
δe = ⎡⎢W1 θx 1 θy 1 W2 θx 2 θy 2 W3 θx 3 θy 3 W4 θx 4 θy 4 ⎤⎥ (7.50)
⎣ ⎦
 Flexión de placas planas 139

Figura 7.14 Grados de libertad del elemento rectangular

Esto permite emplear como función de interpolación de w un polinomio de 12

términos. De todos los posibles se utiliza uno incompleto de grado cuatro, en el que
faltan los términos x 4 , y 4 , x 2y 2 . Por lo tanto, el polinomio interpolante es:

w = α1 + α2x + α3y + α4x 2 + α5xy + α6y 2

(7.51)
+α7x 3 + α8x 2y + α9xy 2 + α10y 3 + α11x 3y + α12xy 3

que se puede poner en la forma:

w = Rα (7.52)

siendo:
• α un vector con todos los coeficientes del polinomio αi
• R un vector fila que contiene todas las potencias de x e y del polinomio.
Particularizando la expresión de w, y sus dos derivadas parciales respecto a x,y para
los cuatro nudos se obtienen doce expresiones del tipo:
Wi = α1 + α2x i + α3yi + α4x i2 + α5x iyi + α6yi2 +

+α7x i3 + α8x i2yi + α9x iyi2 + α10yi3 + α11x i3yi + α12x iyi3 (7.53)
θxi = α3 + α5x i + 2α6yi + α8x i2 + 2α9x iyi + 3α10yi2 + α11x i3 + 3α12x iyi2

θyi = −α2 − 2α4x i − α5yi − 3α7x i2 − 2α8x i yi − α9yi2 − 3α11x i2yi − α12yi3
donde i=1,4. Estas doce ecuaciones se pueden agrupar en la forma:
δe = C α (7.54)

siendo C una matriz constante de dimensión 12x12, cuyos términos son función
únicamente de las coordenadas de los nudos. Se puede demostrar que la matriz C
tiene inversa siempre que las coordenadas de los nudos definan un rectángulo
correctamente.
Despejando α y sustituyendo en la expresión de w se obtiene:
w = R C−1 δe = N δe (7.55)
140 Flexión de placas planas

con lo que se define la matriz de funciones de interpolación N para la flecha w. La

forma explícita de N se puede calcular analíticamente, resultando un conjunto de 12
funciones, que pueden agruparse para los cuatro nudos en la forma:

N = ⎡⎢N1 N2 N3 N4 ⎤⎥ (7.56)
⎣ ⎦
donde Ni son matrices 1x3 que contienen las tres funciones de interpolación del
nudo i. Para simplificar el cálculo de las funciones de interpolación se emplea un
sistema de ejes ξ,η, local al elemento y situado en su centro (cuyas coordenadas son
xc, yc) de tal manera que las coordenadas locales varían entre –1 y +1. Por lo tanto:
x − xc y − yc
ξ= η= (7.57)
a b
El valor de las funciones de interpolación del nudo i es:
1⎡
Ni = ⎢ (ξ ξi + 1)(η ηi + 1)(2 + ξ ξi + η ηi − ξ 2 − η 2 ) ...
2 ⎣

... aξi (ξ ξi + 1)2 (ξ ξi − 1)(η ηi + 1) aηi (ξ ξi + 1)(η ηi + 1)2 (η ηi − 1) ⎤⎥ (7.58)

⎦

7.3.2 Deformaciones unitarias

Sustituyendo el valor de w en la expresión de las deformaciones unitarias se obtiene
la matriz B:
ε = z ∂ w = z ∂ N δe = z B δe = B δe (7.59)

que está compuesta por 4 submatrices y cuyo valor es:

⎡ ∂ 2N ∂ 2 N2 ∂ 2N3 ∂ 2N4 ⎤⎥
⎢ 1
⎢ ∂x 2 ∂x 2 ∂x 2 ∂ x 2 ⎥⎥
⎢
⎢ ∂ 2N ∂ 2 N2 ∂ 2N3 ∂ 2N4 ⎥⎥
B = −z ⎢⎢ 1
(7.60)
⎢ ∂y
2
∂y 2 ∂y 2 ∂y 2 ⎥⎥
⎢ ∂ 2N ∂ 2 N2 ∂ 2N3 ∂ 2N4 ⎥⎥
⎢2 1
2 2 2
⎢ ⎥
⎢⎣ ∂ x ∂y ∂ x ∂y ∂ x ∂y ∂ x ∂y ⎥⎦
También se puede poner como:
B = z ∂ N = z ∂ R C−1 = z Q C−1 (7.61)

donde la matriz Q puede calcularse explícitamente, y vale

⎡ ⎤
⎢ 0 0 0 −2 0 0 −6x −2y 0 0 −6xy 0 ⎥
⎢ ⎥
Q = ∂ R = ⎢0 0 0 0 0 −2 0 0 −2x −6y 0 −6xy ⎥ (7.62)
⎢ ⎥
⎢ 0 0 0 0 −2 0 0 −4x −4y 0 −6x 2 2⎥
−6y ⎥
⎢⎣ ⎦
 Flexión de placas planas 141

En ella se observa que se puede alcanzar un estado de ε constante si α7 a α12 son

cero con lo que desaparecen todos los términos en x,y.

7.3.3 Matriz de rigidez

La expresión general para la matriz de rigidez de un elemento finito es:

K = ∫ BT D B dx dy dz (7.63)
v

Sustituyendo B mediante (7.61) y sacando fuera de la integral los términos que no

dependen de x,y se obtiene:
h3
K = ∫ z 2 BT D B dx dy dz = ∫ BT D B dx dy
12
(7.64)
h 3 −T ⎡
K= C ⎢ ∫ QT D Q dx dy ⎤⎥ C−1
12 ⎣ ⎦
Tras laboriosos cálculos, es posible obtener expresiones explícitas de esta matriz de
rigidez, que están publicadas en la literatura especializada.

7.3.4 Vector de fuerzas nodales equivalentes

El vector de fuerzas nodales equivalentes para una carga distribuida por unidad de
superficie q normal a la placa es:

∫ N q dxdy = C ∫ R
T −T T
Ps = q dxdy (7.65)
s s

cuyo valor puede ser asimismo evaluado analíticamente y en el que se observa la

existencia de componentes tanto según w como según los giros.

7.4. REQUERIMIENTOS DE CONVERGENCIA

El elemento rectangular de placa delgada desarrollado satisface los dos primeros
criterios de convergencia, pero no cumple el tercero, es decir que es incompatible.
Sin embargo este elemento se comporta bastante bien y fue muy utilizado en la
práctica.
Para poner de manifiesto la incompatibilidad de deformaciones, considérese la
situación mostrada en la figura 7.15, en la que dos elementos A y B comparten un
lado 2‐3. Para mayor sencillez, se supone que este lado es la recta x=0. La flecha y los
giros en este lado son:
142 Flexión de placas planas

wx =0 = α1 + α3y + α6y 2 + α10y 3

⎛ ∂w ⎞
θx (x =0) = ⎜⎜ ⎟⎟⎟ = α3 + 2α6y + 3α10y 2 (7.66)
⎜⎝ ∂y ⎠
x =0

⎛ ∂w ⎞
θy (x =0) = − ⎜⎜ ⎟⎟⎟ = −α2 − α5y − α9y 2 − α12y 3
⎝ ∂ x ⎠x =0
Particularizando en el nudo 2 (y=0):
w2 = α1 (7.67)
θx 2 = α3 (7.68)
θy 2 = −α2 (7.69)

Particularizando en el nudo 3 (y=b):

w 3 = α1 + α3b + α6b 2 + α10b 3 (7.70)
θx 3 = α3 + 2α6b + 3α10b 2 (7.71)

θy 3 = −α2 − α5b − α9b 2 − α12b 3 (7.72)

Las 4 ecuaciones (7.67), (7.68), (7.70) y (7.71) permiten hallar los 4 parámetros
α1, α3, α6, α10 en función de las 4 deformaciones nodales w2, w3, θx2, θx3. Por lo tanto
la deformación w y el giro θx quedan unívocamente definidos en el lado x=0
mediante las deformaciones nodales existentes sólo en dicho lado (figura 7.15).

Figura 7.15 Deformación de la frontera entre dos elementos placa

Para definir el giro θy en el lado 2‐3 se necesitan los cuatro parámetros α2, α5, α9, α12
(ver la tercera ecuación de (7.66)). Sin embargo para ello sólo se dispone de las dos
 Flexión de placas planas 143

ecuaciones (7.69) y (7.72), por lo que dicho giro no queda unívocamente definido
empleando únicamente los grados de libertad del lado 2‐3. En consecuencia, cada
uno de los dos elementos A, o B definirá un valor diferente del giro θy, con lo que no
queda garantizada la continuidad del giro θy en el lado. De esta manera se incumple
el tercer criterio de convergencia que exige la continuidad hasta la derivada de orden
1 en los lados entre los elementos.
La figura 7.16 ilustra la situación de incompatibilidad en el giro θy.

4
Y

Y3
A
Y
3
A
B
1 Y

Y
Y2

B
2
X

Figura 7.16 Incompatibilidad del giro normal a la frontera entre 2 elementos placa

Esta conclusión ha sido deducida para una situación particular, pero el fenómeno
puesto de manifiesto es general en elementos placa que empleen la formulación aquí
descrita: existen discontinuidades entre elementos en el giro normal al lado, sea cual
sea su orientación (figura 7.17). A los elementos que muestran este comportamiento
se les denomina no conformes.

tB

tA

Figura 7.17 Continuidad de giros en la frontera entre dos elementos placa

144 Flexión de placas planas

Una primera solución a este problema es adoptar otro polinomio interpolante

diferente, que garantice la compatibilidad del giro normal. En esta línea, Melosh
presentó otro elemento en el que, utilizando como grados de libertad la flecha y los
dos giros, el polinomio interpolante garantiza la compatibilidad del giro normal al
lado. Sin embargo, dicho polinomio no tiene el término (xy ), por lo que no es capaz
de representar estados de curvatura de torsión constante (que es proporcional a la
derivada segunda cruzada), y por lo tanto no satisface los criterios 1 y 2 de
convergencia.
En realidad se puede demostrar que es imposible definir, mediante polinomios
sencillos, unas funciones de interpolación que aseguren la compatibilidad completa,
cuando sólo se usan como parámetros nodales, la flecha w y sus derivadas primeras.
O bien se incumple la compatibilidad del giro normal, o bien el elemento es incapaz
de representar algún estado de tensión constante.
Una solución más adecuada al problema consiste en adoptar otros parámetros,
además de la flecha y el giro, como grados de libertad en los nudos. De esta manera
se pueden generar elementos compatibles, como se describe a continuación.

7.5. ELEMENTOS TRIANGULARES INCOMPATIBLES

De la misma forma que el elemento rectangular ya explicado, se han desarrollado
muchos elementos placa de tres lados. El principal problema para su desarrollo está
en el hecho de que si se emplean como grados de libertad la flecha y los giros en los
tres nudos, se dispone de 9 parámetros para definir la función de interpolación de la
flecha, lo cual no permite emplear un polinomio completo de grado 3, que requiere
un total de 10 parámetros. En función de cuál sea el término del polinomio no
incluido en la función de interpolación, pueden formularse diversos elementos,
aunque muchos de ellos no tienen la calidad suficiente para su empleo práctico.
Un elemento triangular que sí muestra un comportamiento aceptable es el
desarrollado por Cheung, King y Zienkiewicz en 1968. Emplea un polinomio cúbico
incompleto, que le permite representar campos de curvatura constante (es decir
momentos flectores y torsor) de valor arbitrario. Es incompatible, pues no garantiza
la continuidad del giro según la normal a las fronteras, pero sin embargo garantiza la
convergencia. Además se han efectuado sobre su formulación inicial varias
modificaciones, añadiendo términos de grados más elevados, lo cual ha permitido
mejorar la calidad de los resultados obtenidos con él.
Como ejemplo de un elemento incompatible curioso debe citarse el propuesto por
Morley en 1971. Este elemento emplea como grados de libertad la flecha lateral en
cada uno de los tres nudos (sin los giros), y el giro normal en el centro de cada uno de
sus lados. Estos 6 grados de libertad permiten emplear un polinomio completo de
grado 2, el cual permite representar cualquier estado de curvatura constante dentro
 Flexión de placas planas 145

del elemento. Sin embargo el elemento no permite una variación cúbica de la flecha
en sus lados, con lo que no se garantiza la continuidad de los giros en los mismos,
incumpliendo el tercer criterio. A pesar de ello, el elemento converge.

7.6. ELEMENTOS CONFORMES

Para tratar de resolver el problema de los elementos incompatibles se han
desarrollado otros, que se llaman conformes, y que básicamente lo que hacen es
añadir nuevos términos a las funciones de interpolación, de tal manera que se
satisfagan a la vez la compatibilidad de giros en los bordes del elemento, y se
representen estados de deformación unitaria constante.
Para poder añadir estos nuevos términos a las funciones de interpolación, se requiere
introducir grados de libertad adicionales en el elemento. Entre ellos se suelen incluir
los giros según la tangente al punto medio de los lados, o incluso derivadas segundas
de la flecha. Estos elementos pueden tener distinto número de grados de libertad en
cada nudo, y se han desarrollado gran cantidad de ellos. Los resultados obtenidos en
casos comparativos entre distintos tipos de ellos son bastante dispares y no puede
decirse cuál es el mejor de todos, incluyendo en la comparación elementos
conformes y no conformes.
Entre los elementos conformes de cuatro lados, debe mencionarse el elemento
hermítico de 16 grados de libertad desarrollado por Bogner, Fox y Schmidt en 1967.
Tiene cuatro variables en cada nudo: la flecha, los dos giros x e y y la derivada
segunda cruzada de la flecha (figura 7.18). El elemento utiliza como función
interpolante de la flecha w un producto de dos polinomios de Hermite de grado 3,
uno en cada dirección. Este elemento cumple los tres criterios de convergencia, y
resulta ser más preciso que el de Melosh, debido a que utiliza un polinomio
interpolante de mayor grado. En todo caso, no es posible generalizarlo a formas no
rectangulares, lo cual limita su aplicación.
⎧
⎪ wi ⎫ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪ θ xi ⎪
⎪
⎪ ⎪
δi = ⎪
⎨ θyi ⎪ ⎬
⎪
⎪ ⎪
⎪
⎪
⎪ ∂ wi ⎪
2
⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎩ ∂ x ∂ y ⎪
⎪
⎭

Figura 7.18 Elemento rectangular hermítico

En lo que respecta a los elementos triangulares compatibles, se han desarrollado

varios de ellos, con distintos tipos de grados de libertad. La figura 7.19 muestra
algunos de ellos.
146 Flexión de placas planas

⎪⎧⎪ ∂2w ⎪⎫⎪

{w θx θy } {∂∂wn } ⎨w
⎪⎪
⎩⎪
∂w
∂n
⎬
∂s ∂n ⎪⎪
⎭⎪
⎧
⎪
⎪ ∂w ∂w ∂ 2w ∂ 2w ∂ 2w ⎫⎪
⎪
⎨w ⎬
⎪
⎪
⎩ ∂x ∂y ∂x 2 ∂y 2 ∂x ∂ y ⎪
⎪
⎭

Figura 7.19 Elementos triangulares compatibles

 8
Flexión de placas con energía de
esfuerzo cortante

8.1. INTRODUCCIÓN
Cuando la energía debida al esfuerzo cortante no es despreciable, debe emplearse
una teoría adecuada que la tenga en cuenta. Esto ocurre en placas cuyo espesor no
es totalmente despreciable frente a sus dimensiones transversales, aunque si debe
seguir siendo lo suficientemente pequeño para que la tensión en la dirección z sea
despreciable.
Los primeros desarrollos de la teoría de flexión de placas incluyendo la energía de
esfuerzo cortante son debidos a Reissner (1945) y Mindlin (1950). En esta teoría no
se emplea la hipótesis de Kirchhoff, por lo que las secciones rectas normales al plano
medio de la placa no se mantienen normales a dicho plano medio en el estado
deformado. Es decir que el ángulo que gira una sección plana cualquiera no tiene por
qué coincidir con la tangente a la deformada.
Como consecuencia de la diferencia entre giro y tangente a la deformada aparece en
el material una deformación de cortadura que, asociada a las tensiones cortantes
verticales existentes, hace que se acumule energía de cortadura.

8.2. ESTADO DE DEFORMACIÓN

En un punto del plano medio de la placa, las deformaciones son la flecha w y los dos
giros θx y θy. Estas deformaciones se agrupan en el vector de deformaciones w, que
es función de la situación del punto (x,y) dentro de la placa:
148 Flexión de placas con energía de esfuerzo cortante

⎧w ⎪
⎪ ⎫
⎪
⎪ ⎪ ⎪
w = ⎨θx ⎪
⎪ ⎬ (8.1)
⎪
⎪ ⎪
⎪
⎪
⎪θy⎪
⎪
⎪ ⎭
⎩ ⎪
Las deformaciones de un punto cualquiera de la placa, situado a una altura z respecto
al plano medio, son (figura 8.1):
⎧u ⎪
⎪ ⎫ ⎪ ⎧ z θy ⎪
⎫
⎪
⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪ ⎪
u = ⎨ v ⎬ = ⎨−z θx ⎪
⎪ ⎬ (8.2)
⎪
⎪ ⎪ ⎪ ⎪
⎪w ⎪ ⎪ w ⎪
⎪ ⎪ ⎪
⎪ ⎭
⎩ ⎪ ⎩ ⎪ ⎪
⎭
Con este planteamiento el estado de deformación de la placa ya no queda definido
únicamente por la flecha w, sino que intervienen también los giros de cada sección
recta, que ya no coinciden con las derivadas de w.

v
u
w Z

Y Y

Figura 8.1 Estado de deformación interior en placas Mindlin

Las deformaciones unitarias producidas por estas deformaciones son:

∂u ∂θ
εx = =z y
∂x ∂x
∂v ∂θ
εy = = −z x
∂y ∂y
∂u ∂v ∂θ ∂θ
γxy = + = z y −z x (8.3)
∂y ∂ x ∂y ∂x
∂w ∂u ∂w
γzx = + = + θy
∂x ∂z ∂x
 Flexión de placas con energía de esfuerzo cortante 149

∂v ∂w ∂w
γyz = + = −θx +
∂z ∂y ∂y
En esta expresión se observa que:
• Las deformaciones unitarias varían linealmente en el espesor de la placa, como
en la teoría clásica.
• No aparecen las derivadas segundas de la flecha w como en las placas delgadas,
sino las derivadas primeras de los giros de la sección.
• Aparecen dos deformaciones de cortadura verticales que se deben a la diferencia
entre el giro de la sección y la pendiente a la deformada (derivada), y que no
existían en la teoría clásica.

Figura 8.2 Rotaciones en placas Mindlin

• Las tres primeras deformaciones unitarias se pueden poner en la forma:

⎧
⎪ ∂θy ⎫ ⎪
⎪
⎪ ⎪
⎪
⎧ ⎫ ⎪
⎪ ∂ x ⎪
⎪
⎪
⎪ εx ⎪
⎪ ⎪ ⎪
⎪ ⎪ ⎪
⎪ ∂θx ⎪ ⎪
⎪ ⎪
εF = ⎨ εy ⎬ = z ⎨ − ⎬=zκ (8.4)
⎪
⎪ ⎪
⎪ ⎪
⎪ ∂ y ⎪
⎪
⎪
⎪ γ ⎪ ⎪ ⎪
⎩ xy ⎪
⎭ ⎪
⎪ ∂θ ∂θ ⎪
x⎪
⎪
⎪
y
− ⎪
⎪
⎪ ∂y
⎩ ∂x ⎪
⎪
⎪
⎭
Donde se ha introducido el vector de curvaturas κ propio de esta teoría, que es
diferente del de la teoría clásica pues ahora contiene las derivadas primeras de los
giros. Introduciendo el operador de derivación ∂1 de orden 1, se obtiene:
150 Flexión de placas con energía de esfuerzo cortante

⎡ ∂ ⎤⎥
⎢0 0
⎢ ∂ x ⎥⎥ ⎪
⎧w ⎪
⎫
⎢ ⎪ ⎪
⎢ ∂ ⎥⎪ ⎪
⎪
εF = z κ = ⎢ 0 − 0 ⎥ ⎨θx ⎪⎬ = z ∂1 w (8.5)
⎢ ∂y ⎥⎪ ⎪
⎢ ⎥ ⎪θ ⎪
⎪ ⎪
⎢ ∂ ∂ ⎥⎪ ⎪ y⎪
⎪
⎢0 − ⎥ ⎭
⎩
⎢⎣ ∂x ∂y ⎥⎦
• Las deformaciones de cortadura valen:
⎧ ∂w
⎪ ⎫ ⎡∂ ⎤ ⎧w ⎪ ⎫
γ ⎪
⎪ + θy ⎪ ⎪
⎪ ⎢ 0 1⎥ ⎪⎪ ⎪ ⎪⎪
⎧
⎪
⎪ ⎫
xz ⎪
⎪ ⎪ ∂x
⎪ ⎪
⎪ ⎢ ∂x ⎥⎪ ⎪
γ = ⎨γ ⎬ = ⎨ ⎬= ⎢ ⎥ ⎨⎪θx ⎬⎪ = ∂ 2 w (8.6)
⎪
⎪ yz ⎪
⎩ ⎪ ⎪⎪⎪ ∂w − θx ⎪⎪⎪
⎭ ⎢∂
⎢
⎪ ⎪
−1 o ⎥⎥ ⎪⎪θ ⎪ ⎪
⎪⎪ ∂y
⎩ ⎭⎪
⎪ ⎣⎢ ∂y ⎥⎦ ⎩⎪⎪ y ⎭
⎪⎪
Se observa que estas deformaciones son constantes en todo el espesor de la placa,
mientras que en la solución real varían de forma parabólica. Esto es debido a que la
hipótesis de deformación efectuada no es perfecta, ya que en la realidad las
secciones rectas perpendiculares al plano medio no se mantienen rectas en el estado
deformado, sino que se distorsionan debido a la deformación de cortadura.

8.3. RELACIÓN TENSIÓN DEFORMACIÓN

La relación entre las tensiones y las deformaciones unitarias es la correspondiente al
estado plano de tensiones, ya que la tensión en dirección z es nula. Suponiendo un
material isótropo dicha relación es:
⎡ ⎤
⎢ ⎥
⎧
⎪⎪ σx ⎫
⎪ ⎢1 ν 0 ⎥ ⎛⎜⎧ ⎪ εx ⎫ ⎪ ⎧
⎪ ε0x ⎫ ⎪⎞⎟
⎪⎪ ⎪⎪⎪ E ⎢ ⎜⎪
⎥ ⎜⎪ ⎪ ⎪ ⎪ ⎪ ⎪
⎪⎟⎟
σ F = ⎨ σy ⎬ = ⎢ν 1 0 ⎥ ⎜⎜⎨⎪ εy ⎬⎪ − ⎨⎪ ε0y ⎬⎪⎟⎟⎟ (8.7)
⎪ ⎪⎪ 1 − ν 2
⎪ ⎢ ⎥ ⎜⎜⎪ ⎪ ⎪ ⎪⎟⎟⎟
⎪τxy ⎪
⎪ ⎪
⎢ 1 − ν ⎥ ⎜⎝⎪
⎪ γ
⎪ ⎭
⎪
⎪ ⎪
⎪
⎪ ⎩
⎪ γ ⎪
⎪
⎪⎠
⎩ ⎭ ⎢0 0 ⎥ ⎩ xy 0 xy
⎭
⎢⎣ 2 ⎦⎥
σ F = D(εF − ε0 ) (8.8)
que es la misma que se emplea en la teoría clásica. Además, las dos tensiones
cortantes verticales se relacionan con las correspondientes deformaciones de
cortadura, a través de la parte correspondiente de la ecuación constitutiva:

⎪τxz ⎫
⎧ ⎧G γxz ⎪
⎪ ⎪ ⎫
τ=⎪ ⎪ ⎪ =G γ
⎨τ ⎪ ⎬=⎨ ⎬ (8.9)
⎪ yx ⎪
⎪ ⎭
⎩ ⎪
⎪ ⎩
⎪G γ ⎪⎪
yx
⎭
donde G es el módulo de elasticidad en cortadura.
 Flexión de placas con energía de esfuerzo cortante 151

8.4. ESFUERZOS INTERNOS

La definición de los esfuerzos internos en la placa en función de las tensiones es la
misma que en la teoría clásica. Sustituyendo las tensiones en función de las
deformaciones unitarias y éstas en función de las curvaturas e integrando en z, se
obtiene, para los momentos de flexión y torsión:
h3
M= D κ + M0 (8.10)
12
Sustituyendo las curvaturas κ se obtiene:
⎧⎪ ⎛ ⎞ ⎫
⎪⎪ D ⎜⎜ ∂θy − ν ∂θx ⎟⎟ ⎪⎪⎪
⎪⎧⎪ M x ⎪⎫⎪ ⎪⎪
⎪ ⎜⎝ ∂ x ∂y ⎠⎟ ⎪⎪ ⎪⎧
⎪⎪ ⎪ M 0x ⎪⎫⎪
⎪⎪ ⎪⎪ ⎪⎪ ⎛ ∂θ ∂θ ⎞ ⎪ ⎪⎪ ⎪⎪
M = ⎨⎪ M y ⎬⎪ = ⎪⎨ D ⎜⎜ν y − x ⎟⎟⎟ ⎪⎬ + ⎨⎪ M 0y ⎬⎪ (8.11)
⎪⎪ ⎪ ⎪ ⎜⎝ ∂ x ∂y ⎠ ⎪⎪ ⎪⎪ ⎪
⎪⎪M xy ⎪⎪⎪ ⎪⎪⎪ ⎪⎪ ⎪⎪M ⎪⎪⎪
⎩⎪ ⎭⎪ ⎪ D(1 − ν ) ⎛ ∂θy ∂θx ⎞⎟⎪ ⎩⎪ 0xy ⎭⎪
⎪ ⎪ ⎜⎜⎜ − ⎟⎪⎪
⎪ 2 ⎝ ∂y ⎟
∂ x ⎠⎪⎭⎪
⎪
⎩
Los esfuerzos cortantes son la resultante de las tensiones cortantes verticales:
+h / 2

Q= ∫ τ dz = h G γ (8.12)
−h/2

8.5. ECUACIONES DE EQUILIBRIO

La obtención de las ecuaciones de equilibrio en esta teoría sigue los mismos pasos
que en la teoría clásica, con lo que obtienen las mismas ecuaciones de equilibrio
expresadas en función de los esfuerzos internos:
∂M x ∂M xy
+ = Qxz
∂x ∂y

∂ M xy ∂M y
+ = Qyz (8.13)
∂x ∂y
∂Qzx ∂Qzy
+ + qz = 0
∂x ∂y
Sustituyendo en ellas los valores de los momentos y los cortantes en función de las
tres deformaciones, mediante (8.11) y (8.12), se obtiene:
∂ ⎜⎛ ∂θy ∂θ ⎞ 1 − ν ∂ ⎛⎜ ∂θy ∂θx ⎞⎟ G t ⎛ ∂w ⎞⎟
⎜⎜ + ν x ⎟⎟ + ⎜ − ⎟= ⎜ + θ ⎟
∂x ⎝ ∂x ∂y ⎠⎟ 2 ∂y ⎝⎜ ∂y ∂ x ⎠⎟ D ⎜⎝ ∂ x
y
⎠⎟
152 Flexión de placas con energía de esfuerzo cortante

∂ ⎛⎜ ∂θx ∂θy ⎞⎟ 1 − ν ∂ ⎛⎜ ∂θy ∂θx ⎞⎟ G t ⎛⎜ ∂w ⎞

⎜⎜ +ν ⎟+ ⎜ − ⎟= ⎜ − θx ⎟⎟ (8.14)
∂y ⎝ ∂y ∂ x ⎠⎟ 2 ∂ x ⎝⎜ ∂y ∂ x ⎠⎟ D ⎝⎜ ∂y ⎠⎟

∂ ⎛ ∂w ⎞ ∂ ⎛ ∂w ⎞ q
⎜⎜ + θy ⎟⎟⎟ + ⎜⎜ − θx ⎟⎟⎟ + z = 0
∂x ⎝ ∂x ⎠ ∂y ⎝⎜ ∂y ⎠ Gt

Se trata de un sistema de tres ecuaciones con tres incógnitas, que son las
deformaciones del plano medio de la placa. Este sistema de ecuaciones es de orden
2, en comparación con el orden 4 obtenido en la teoría clásica. Para resolverlo por el
MEF será necesario emplear funciones de interpolación con continuidad C0, en lugar
de las funciones de tipo C1 necesarias en la teoría clásica.

8.6. EXPRESIÓN DE LA ENERGÍA ELÁSTICA

La densidad de energía elástica acumulada en un punto de la placa incluye dos
sumandos: uno debido a las tres tensiones incluidas en el vector σ y otro debido a las
dos tensiones cortantes τ.
ε γ ε ε γ

∫σ d εF + ∫ τ d γ = ∫ε D d εF − ∫ ε D d εF + ∫ G γT d γ (8.15)
T T T T
U0 = F F 0
0 0 0 0 0

1 1
U0 = εTF D εF − εT0 D εF + G γT γ (8.16)
2 2
La energía elástica acumulada en toda la placa es:

∫ε D εF dv − ∫ εTF D ε0 dv + 21 ∫ G γT γ dv
T
U = 1
2 F (8.17)

Sustituyendo las deformaciones unitarias εF en función de las curvaturas se obtiene:

∫b D b z 2dz dA − ∫ z bT D ε0dz dA + 21 ∫ G γT γ dz dA
T
U = 1
2

En el primer término se puede integrar en la coordenada z, y en el segundo se

identifican los momentos debidos a las deformaciones unitarias iniciales:
3
1 h 1
∫ 12 b D b dA + ∫ bT M0 dA + ∫ Gkt
T
U = γT γ dA (8.18)
2 2
En esta ecuación se ha incluido, en el término del cortante, un factor k denominado
factor de corrección de la energía de cortante. Su misión es corregir el valor de la
energía calculada con la distribución uniforme de tensiones cortantes e igualarla al
valor real de la energía, que está producida por una distribución parabólica de las
tensiones cortantes.
 Flexión de placas con energía de esfuerzo cortante 153

8.7. ELEMENTOS FINITOS PARA PLACAS CON ENERGÍA DE CORTANTE

8.7.1 Funciones de interpolación

Para resolver por el MEF un problema de flexión de placas empleando la teoría de
Mindlin‐Reissner se puede usar cualquier tipo de función de interpolación que tenga
continuidad C0. En particular todas las utilizadas para los problemas de elasticidad
plana son aplicables, tanto en formulación de Lagrange como Serendipity.
Estas funciones de interpolación se usan para interpolar las tres incógnitas del
problema, que son la flecha w y los dos giros:
w= ∑N W i i
θx = ∑ N i θxi θy = ∑ N i θyi (8.19)

Las tres leyes de interpolación se pueden agrupar como:

⎧
⎪ W1 ⎫ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪ θx 1 ⎪⎪
⎪
⎪ ⎪
⎪
⎪⎧w ⎪⎫ ⎡ ⎤
0 ⎥ ⎪ θy 1 ⎪⎪
⎪
⎪⎪⎪ ⎪⎪⎪ ⎢N 1 0 0 ... ... N n 0
⎪ ⎪
⎢ ⎥⎪ ⎪
⎨θx ⎬ = ⎢ 0 N 1 0 ... ... 0 Nn 0 ⎥ ⎪⎨ ... ⎪⎬ (8.20)
⎪⎪ ⎪⎪ ⎢ ⎥⎪ ⎪ ⎪
⎪
⎪⎪θy ⎪⎪ ⎢0 0 N 1 ... ... 0 0 ⎥ W
Nn ⎥ ⎪ n ⎪
⎪
⎩⎪ ⎭⎪ ⎢⎣ ⎦⎪ ⎪ ⎪
⎪
⎪ θxn ⎪⎪
⎪
⎪ ⎪
⎪
⎪
⎪ θyn ⎪⎪
⎪ ⎭
⎩ ⎪
Que puede ponerse en la forma matricial habitual:
w = N δe (8.21)

donde δe es el vector de todas las deformaciones nodales del elemento (figura 8.3):
T
δe = ⎡⎢W1 θx 1 θy 1 W2 θx 2 θy 2 ... Wn θxn θyn ⎤⎥ (8.22)
⎣ ⎦

w3 w2
Z
Y2
Y
Y3
X X3 X2
w4 w1

Y1

Y4
X1
X4

Figura 8.3 Grados de libertad del elemento placa rectangular

154 Flexión de placas con energía de esfuerzo cortante

Las funciones de interpolación se definen habitualmente en un sistema de

coordenadas local al elemento. A este efecto se emplean los mismos sistemas usados
para los elementos de elasticidad plana: para elementos de cuatro lados se emplean
las coordenadas ξ,η normalizadas tal que el elemento es un cuadrado de lado 2
(figura 8.4) y para triángulos de emplean las tres coordenadas de área Li .

Figura 8.4 Coordenadas locales para elementos placa

8.7.2 Interpolación de coordenadas.

Para definir la forma del elemento en el plano XY se emplea la interpolación de
coordenadas habitual en el MEF. Pueden emplearse elementos iso, sub o super
paramétricos, tal y como se ha explicado para el problema plano. De hecho la forma
de un elemento placa y de un elemento para elasticidad plana es la misma. Por lo
tanto la ley de interpolación es del tipo:
⎧⎪x 1 ⎫⎪
⎪⎪ ⎪⎪
⎪y1 ⎪
⎧⎪⎪x ⎫⎪⎪ ⎡N 1 0 N 2 0 ...⎤ ⎪⎪ ⎪⎪
⎨y ⎬ = ⎢⎢ ⎥ ⎪⎨x 2 ⎪⎬ (8.23)
⎪⎪⎩ ⎪⎪⎭ ⎢ 0 N 1 0 N 2 ...⎥⎥ ⎪⎪ ⎪⎪
⎣ ⎦ ⎪y2 ⎪
⎪⎪ ⎪⎪
⎪⎪...⎪⎪
⎩ ⎭
x = N xe (8.24)

El vector xe agrupa a las coordenadas (x,y) de todos los nudos del elemento. Esta ley
permite obtener elementos con lados curvos en el plano de la placa.

Figura 8.5 Elementos placa con lados curvos