Universidad de Concepción

Departamento de Ingeniería Mecánica

Estructuras Aeroespaciales
UNIDAD 3: Elasticidad 2D

Luis Quiroz Larrea

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 1

 Elasticidad 2D
Temario

• Esfuerzo plano
• Deformación plana
• Axisimetría
• Soluciones analíticas
• Función esfuerzo de Airy
– Métodos inverso y semi-inverso
– Principio de Saint Venant
• Efectos Ambientales

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 2

 Elasticidad 2D
Esfuerzo plano

Físicamente trata de placas con componentes de esfuerzos

normales y de corte sólo en el plano de la placa

El espesor h es mucho más pequeño que las otras dimensiones

Por lo tanto el esfuerzo en dicha dirección es pequeño o nulo.

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 3

 Elasticidad 2D
Esfuerzo plano

Así: σ 33 = σ 13 = σ 23 = 0

Con lo que las ecuaciones de elasticidad se reducen:

Equilibrio
∂ σ 11 ∂ σ 12
+ + f1 = 0
∂ x1 ∂ x2 ∂σ ij
⇒ + fi = 0 i , j = 1,2
∂ σ 21 ∂ σ 22 ∂ xj
+ + f2 = 0
∂ x1 ∂ x2

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 4

 Elasticidad 2D
Esfuerzo plano

Deformaciones – Esfuerzo (Material Ortotrópico)

1
ε 11 = (σ 11 − ν 12σ 22 )
E11
1
ε 22 = (σ 22 − ν 21σ 11 )
E 22
1
ε 12 = σ 12
2G12
1
ε 33 = − (ν 31σ 11 + ν 32σ 22 )≠ 0
E33
El estado de esfuerzo plano es una aproximación y desprecia el
esfuerzo en la dirección 3:
∂ u3
ε 33 ≠ 0 ⇒ ≠ 0 ¡inconsistencia!
∂ x3

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 5

 Elasticidad 2D
Esfuerzo plano

Deformaciones -Desplazamientos

∂ u1 
ε 11 = 
∂ x1

∂ u2 
ε 22 =
∂ x2 
∂ u3 
ε 33 = 
∂ x3 
1  ∂ u1 ∂ u2  
ε 12 = ε 21 =  +  
2  ∂ x 2 ∂ x1  

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 6

 Elasticidad 2D
Esfuerzo plano

Incógnitas - Ecuaciones

Incógnitas:
•Esfuerzos (3) σ x ,σ y ,τ xy
•Deformaciones (4) ε x ,ε y ,ε z ,γ xy
•Desplazamientos (3) u x ,uy ,uz

Ecuaciones:
•Equilibrio (2)
•Deformación – desplazamientos (4)
•Esfuerzo – deformación (4)

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 7

 Elasticidad 2D
Deformación plana

Físicamente trata de cuerpos prismáticos con una

dimensión (x3) mucho mayor que las otras dos (x1
y x2).

Suponemos que el cuerpo es “infinito” en

dirección x3, y la cargas importantes están en el
plano x1 x2 y no cambian con x3:
∂ ∂
= = 0
∂ x3 ∂ z
Esto implica que no hay gradiente de
desplazamiento a lo largo de z

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 8

 Elasticidad 2D
Deformación plana

Así las ecuaciones de elasticidad se reducen:

Equilibrio
∂ σ 11 ∂ σ 12
+ + f1 = 0 primaria
∂ x1 ∂ x2
∂ σ 21 ∂ σ 22
+ + f2 = 0 primaria
∂ x1 ∂ x2
∂ σ 31 ∂ σ 32
+ + f3 = 0 secundaria
∂ x1 ∂ x2

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 9

 Elasticidad 2D
Deformación plana

Deformaciones -Desplazamientos

∂ u1 
ε 11 = 
∂ x1 
∂ u2   ∂ 
ε 22 =  ε 33 = ε 13 = ε 23 = 0  = 0 y u3 = 0 
∂ x2   ∂ x3 
1  ∂ u1 ∂ u2  
ε 12 = ε 21 =  +  
2  ∂ x 2 ∂ x1  

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 10

 Elasticidad 2D
Deformación plana

Deformaciones – Esfuerzo (Material Ortotrópico)

1
ε 33 = (σ 33 − ν 31σ 11 − ν 32σ 22 )= 0⇒ σ 33 = ν 31σ 11 + ν 32σ 22 ≠ 0
E33
1
ε 11 = ( (1 − ν 13ν 31 )σ 11 − (ν 12 + ν 13ν 32 )σ 22 )
E11
1
ε 22 = ( (1 − ν 23ν 32 )σ 22 − (ν 21 + ν 23ν 31 )σ 11 )
E 22
1
ε 12 = σ 12
2G12
ε 13 = ε 23 = 0
El estado de deformación plana es una aproximación y hay
esfuerzos en la dirección 3 a pesar de no existir carga en dicha
dirección.

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 11

 Elasticidad 2D
Deformación plana

Incógnitas - Ecuaciones

Incógnitas:
•Esfuerzos (4) σ x ,σ y ,σ z ,τ xy

•Deformaciones (3) ε x ,ε y ,γ xy
•Desplazamientos (3) u x ,uy ,uz

Ecuaciones:
•Equilibrio (3)
•Deformación – desplazamientos (3)
•Esfuerzo – deformación (4)

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 12

 Elasticidad 2D
Axisimetría

Físicamente trata de cuerpos axisimétricos cuyo eje de simetría es

x3 ó z y su sección transversal está definida por r(z) y θ .

Suponemos que las cargas importantes están en dirección radial y

z y no cambian con θ:
∂
= 0 y uθ = 0
∂θ
Esto implica que no hay gradiente de desplazamientos,
deformaciones o esfuerzos a lo largo de θ

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 13

 Elasticidad 2D
Axisimetría

Deformaciones -Desplazamientos
∂ ur
ε rr =
∂r
1 ∂ uθ ur ur
ε θθ = + =
r ∂θ r r
∂ uz
ε zz =
∂z
∂ uθ 1 ∂ ur uθ
ε rθ = ε θ r = + − = 0
∂r r ∂θ r
1 ∂ u z ∂ uθ
ε θ z = ε zθ = + = 0
r ∂θ ∂z
∂ ur ∂ u z
ε rz = ε zr = +
∂z ∂r

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 14

 Elasticidad 2D
Axisimetría

Deformaciones – Esfuerzo (Material Isotrópico)

1
ε rr = [σ rr − ν (σ θ θ + σ zz ) ]
E
1
ε θ θ = [σ θ θ − ν (σ rr + σ zz ) ]
E
1
ε zz = [σ zz − ν ( σ rr + σ θ θ ) ]
E
2(1 + ν )
ε rθ = σ rθ = 0 ⇒ σ rθ = 0
E
2(1 + ν )
ε θz = σ θz = 0 ⇒ σ θz = 0
E
2(1 + ν )
ε rz = σ rz
E
.

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 15

 Elasticidad 2D
Axisimetría

Así las ecuaciones de elasticidad se reducen:

Equilibrio

∂ σ rr ∂ σ rz σ rr − σ θ θ
+ + + fr = 0
∂r ∂z r
∂ σ zr ∂ σ zz σ zr
+ + + fz = 0
∂r ∂z r

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 16

 Elasticidad 2D
Axisimetría

Incógnitas - Ecuaciones

Incógnitas:
•Esfuerzos (4) σ r ,σ θ ,σ z ,τ rz
•Deformaciones (4) ε r , ε θ , ε z , γ rz
•Desplazamientos (2) ur , u z

Ecuaciones:
•Equilibrio (2)
•Deformación – desplazamientos (4)
•Esfuerzo – deformación (4)

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 17

 Elasticidad 2D
Soluciones Analíticas

El problema general de elasticidad posee 15 ecuaciones y 15

incógnitas (3D). Este debe ser resuelto para un cuerpo genérico
bajo ciertas cargas genéricas sujeto a condiciones de borde (CB)
prescritas.
Se encuentran dos tipos de condiciones de borde:
3.Normal (esfuerzos prescritos)
4.Geométrica (desplazamientos prescritos)
Sólo una de ellas debe tenerse a la vez en el contorno.
Para resolver el sistema de ecuaciones existen básicamente dos
procedimientos de solución:
•Exacta (satisface todas las ecuaciones y CB)
•Numérica (satisface tanto como sea posible ecuaciones y CB)
A continuación se verán procedimientos exactos para problemas
2D. Numéricos serán vistos en el curso Elementos Finitos

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 18

 Elasticidad 2D
Funciones esfuerzo

Relaciona 6 funciones esfuerzo (3 en caso 2D) definidas de forma

que satisfagan las ecuaciones de equilibrio.

Consideremos:
•Esfuerzo plano
•Material isotrópico
8 ecuaciones primarias 8 incógnitas primarias
- 2 equilibrio - 3 deformaciones
- 3 deformación – desplazamiento - 2 desplazamientos
-3 esfuerzo – deformación - 3 esfuerzos

Definamos la función esfuerzo de Airy: φ ( x , y )

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 19

 Elasticidad 2D
Funciones esfuerzo

Tal que:
∂ 2φ ∂ 2φ ∂ 2φ
σ xx = +V σ yy = +V σ xy = −
∂y 2
∂x 2
∂ x∂ y

donde: V es una función potencial para las fuerzas de cuerpo

∂V ∂V
fx = − fy = −
∂x ∂y
∂ f x ∂ fy
V existe si ∇ × f = 0 ó = (fuerzas conservativas)
∂y ∂x

Reemplazando en las ecuaciones de equilibrio

∂ σ 11 ∂ σ 12
+ + f1 = 0
∂ x1 ∂ x2
∂ σ 21 ∂ σ 22
+ + f2 = 0
∂ x1 ∂ x2
Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 20
 Elasticidad 2D
Funciones esfuerzo

Se obtiene para la primera ecuación:

∂  ∂ 2φ  ∂  ∂ 2φ  ∂ V
 2 + V  +  −  − = 0
∂x ∂y  ∂ y  ∂ x∂ y  ∂ x
∂ 3φ ∂V ∂ 3φ ∂V
⇒ + − − = 0
∂ x∂ y 2
∂ x ∂ x∂ y 2
∂x

y la segunda ecuación:
∂  ∂ 2φ  ∂  ∂ 2φ  ∂V
 −  +  2 + V  − = 0
∂ x  ∂ x∂ y  ∂ y  ∂ x  ∂y
∂ 3φ ∂ 3φ ∂V ∂V
⇒ − 2 + + − = 0
∂x ∂y ∂x ∂y ∂y ∂y
2

Se satisfacen automáticamente usando la función esfuerzo de Airy

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 21

 Elasticidad 2D
Funciones esfuerzo
¿Significa que cualquier función φ ( x , y ) es válida?
No. Debe satisfacer todavía las relaciones deformación – desplazamiento y
esfuerzo – deformación.

 xx E (σ xx − ν σ yy )
1
ε =

Esfuerzo - Deformación ε
 yy =
1
(σ yy − ν σ xx )
 E
 ε xy = 1 σ xy = (1 + ν ) σ
 2G E
xy

tal que:
1  ∂ 2φ ∂ 2φ  (1 − ν )
ε xx =  2 − ν +
2 
V
E  ∂y ∂x  E
1  ∂ 2φ ∂ 2φ  (1 − ν )
ε yy =  2 − ν + V
E ∂x ∂ y 2  E

ε xy = −
1 ∂ 2φ
= −
( 1 + ν ) ∂ 2φ
2G ∂ x∂ y E ∂ x∂ y
Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 22
 Elasticidad 2D
Funciones esfuerzo
Introduciendo éstas en la ecuación de compatibilidad
∂ 2ε xx ∂ ε yy ∂ 2ε xy
2
+ = 2
∂y 2
∂x 2
∂ x∂ y

y arreglando términos se obtiene

∂ 4φ ∂ 4φ ∂ 4φ  ∂ 2(∆ T ) ∂ 2(∆ T )   ∂ 2V ∂ 2V 
+ 2 2 2+ = − Eα  +  − (1 − ν )  2 + 2
∂x 4
∂x ∂y ∂y 4
∂ x 2
  ∂ y 2
 ∂ x ∂ y 

Término térmico aún no considerado
α = coeficiente de expansión térmica
∆ T = diferencial de temperatura

Ecuación básica para esfuerzo plano de una placa isotrópica en forma de función
esfuerzo ( φ ( x , y ) es un escalar)

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 23

 Elasticidad 2D
Funciones esfuerzo
El operador de Laplace en 2D es
∂2 ∂2 ∂4 ∂4 ∂4
∇ =2
+ ⇒ ∇ ∇ = ∇ =
2 2 4
+ 2 2 2+
∂x 2
∂y 2
∂x 4
∂x ∂y ∂ y4

que es el operador biarmónico. Por lo tanto la ecuación anterior se puede escribir:

∇ 4φ = − Eα ∇ 2
( ∆ T ) − (1 − ν ) ∇ 2V
En ausencia de efectos térmicos y de fuerzas de cuerpo se obtiene la forma
homogénea: 4
∇ φ = 0

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 24

 Elasticidad 2D
Funciones esfuerzo
METODO INVERSO
Para casos de esfuerzo plano sin efectos térmicos ni fuerzas de cuerpo.
3. Se asume una función esfuerzo φ ( x , y ) que satisface la ecuación biarmónica
4. Se determinan los esfuerzos a partir de la función esfuerzo
5. Se imponen las condiciones de borde de esfuerzos
6. Se encuentra el problema estructural que se satisface.

• Los matemáticos ya crearon muchas soluciones por lo que hay funciones

esfuerzo disponibles que resuelven problemas estructurales específicos
• Como son soluciones lineales se puede aplicar el principio de superposición
para combinar soluciones que resuelven otros problemas estructurales
• El método inverso lleva a una solución exacta. En la realidad difícilmente se
puede obtener una solución exacta.

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 25

 Elasticidad 2D
Funciones esfuerzo
METODO SEMI-INVERSO
Básicamente es el mismo método inverso excepto que:
• Se hacen simplificaciones para obtener ecuaciones más simples de resolver.
Estas simplifcaciones pueden ser en relación a:
– Componentes de esfuerzos
– Componentes de desplazamientos
• Se hacen suposiciones basadas en la intuición física, evidencia experimental,
experiencia previa, etc.
• Suposiciones pueden deberse a:
– Condiciones de borde no completamente cumplidas
– Otras ... (Como asumir estado de esfuerzos plano)

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 26

 Elasticidad 2D
Funciones esfuerzo
Ejemplo 1: Sin efectos térmicos ni fuerzas de cuerpo.

Sea φ = C1y que satisface ∇ φ = 0

2 4

∂ 2φ ∂ 2φ ∂ 2φ
σ xx = = 2C1 σ yy = = 0 σ xy = − = 0
∂y 2
∂x 2
∂ x∂ y
¿Que problema resuelve?
¡Evidentemente un estado de esfuerzos uniaxial!.

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 27

 Elasticidad 2D
Funciones esfuerzo
Condiciones de borde:
En x =  ,0
σ0
σ xx = σ 0 = 2C1 ⇒ C1 =
2
σ xy = 0

w
En y = ±
2
σ yy = 0
σ xy = 0
σ0 2
Así: φ = y
2

Las deformaciones pueden ser evaluadas usando las relaciones esfuerzo –

deformación.
Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 28
 Elasticidad 2D
Funciones esfuerzo
Los desplazamientos pueden ser determinados utilizando las relaciones
deformaciones – desplazamientos e integrando y aplicando las CB de
desplazamiento. σ
u= 0x
En este caso se obtiene: E
ν
v= − σ 0y
E
En un ensayo de tracción:
x

Las mordazas no permiten desplazamiento en la dirección y v = 0 en x =  ,0

ν
Sin embargo la solución anterior da v = − σ 0 y en x =  ,0
E
Lejos de los efectos de introducción de la carga, la solución anterior es válida.
=>Principio de Saint-Venant

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 29

 Elasticidad 2D
Principio de Saint Venant
“Si la distribución de carga en una pequeña sección de la superficie de un
cuerpo elástico es reemplazada por otra carga con la misma fuerza y
momento resultante que la carga original, no hay entonces cambios
apreciables en los esfuerzos del cuerpo excepto en la región cerca de la
superficie donde la carga ha sido alterada”

Discusión:
•¿Cuán lejos de la zona de aplicación de carga este principio es válido?
•¿Qué importancia tiene la zona de aplicación de carga en la resistencia de una
estructura?

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 30

 Elasticidad 2D
Funciones esfuerzo
Ejemplo 2: Distribución de esfuerzos alrededor de un agujero.
Placa infinita sujeta a un esfuerzo de tracción uniforme lejos del agujero.

Como la zona de interés es un círculo parece conveniente usar coordenadas

polares: ∂ 2φ 
σ = + V 1∂φ2
1 ∂ φ
xx σ = + 2 +V
∂y 2 rr
r ∂r r ∂θ 2

∂ φ
2
 ∂ 2φ
σ yy = + V  ⇒ σ θθ = +V
∂x 2
 ∂r 2

∂ φ 
2
1 ∂ 2φ 1 ∂φ
σ xy = − σ = − + 2

∂ x∂ y 
rθ
r ∂ r∂ θ r ∂ θ
Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 31
 Elasticidad 2D
Funciones esfuerzo
Tal como en el caso cartesiano, estos esfuerzos están definidos para satisfacer
automáticamente las ecuaciones de equilibrio, donde
∂V 1∂V
fr = − fθ = −
∂r r ∂θ

∂ fθ fθ 1 ∂ fr
y V existe si: ∇ × f = 0 ó + =
∂r r r ∂θ
La ecuación resultante es nuevamente:

∇ 4φ = − Eα ∇ 2
( ∆ T ) − (1 − ν ) ∇ 2V
pero el operador de Laplace en coordenadas polares es:

∂2 1 ∂ 1 ∂2
∇ = 2
+ + 2
∂r 2
r ∂r r ∂θ 2

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 32

 Elasticidad 2D
Funciones esfuerzo
Tal como en el ejemplo anterior suponemos primero una función esfuerzo
φ ( r ,θ ) que satisface ∇ 4φ = 0

[ 
] B 
φ ( r ,θ ) = A0 + B0 ln r + C0r 2 + D0r 2 ln r +  A2r 2 + 22 + C2r 4 + D2  cos 2θ
r
 
satisface las ecuaciones de equilibrio. También se demuestra (otros ya lo hicieron)
que valores no nulos de D0 resultan en desplazamientos v no únicos por lo que
dicho coeficiente debe ser anulado.
Determinamos a continuación los esfuerzos:
B0  6B2 4D2 
σ rr = 2
+ 2C0 +  − 2 A2 − 4
− 2  cos 2θ
r  r r 
B  6B 
σ θθ = − 20 + 2C0 +  2A2 + 42 + 12C2r 2  cos 2θ
r  r 
 6B 2D 
σ rθ =  2 A2 − 42 + 6C2r 2 − 2 2  sin 2θ
 r r 

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 33

 Elasticidad 2D
Funciones esfuerzo
El término que contiene r2 se hace cada vez mayor a medida que r se hace mayor.
El esfuerzo podría llegar a ser infinito. Esto no es posible por lo que el coeficiente
C2 debe ser nulo.
Con esto quedan 5 constantes por determinar: A2 ,B0 ,B2 ,C0 ,D2 mediante las
condiciones de borde.
•En el borde del agujero no hay esfuerzos:
⇒ r = a : σ rr = 0, σ θ r = 0

•En el borde y no hay esfuerzos:

⇒ y = ±∞ : σ yy = 0, σ xy = 0

•En los bordes x el esfuerzo es igual al esfuerzo aplicado

⇒ x = ±∞ : σ xx = σ 0, σ xy = 0

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 34

 Elasticidad 2D
Funciones esfuerzo
Como estamos tratando con coordenadas polares necesitamos cambiar las dos
últimas CB.

~x   r    ~xx  ~xy   x   cos θ sin θ   x 

 ~ =   =    =  
 y   θ    y~x  y~y   y   − sin θ cos θ   y 

[σ~] = [ ][σ ][ ]T ⇒ σ~ij =  ~i k  ~j lσ kl

 σ~xx σ~xy   σ σ rθ    ~xx  ~xy   σ xx σ xy    ~xx  y~x 

⇒  = rr
 = 
σ~xy σ~yy   σ σ θθ   y~x  y~y   σ σ 
yy    ~  y~y 
  rθ xy xy

σ σ rθ   cos θ sin θ   σ xx σ xy   cos θ − sin θ 

 =  − sin θ
rr
σ  
 rθ σ θθ   cos θ   σ xy σ yy   sin θ cos θ 
Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 35
 Elasticidad 2D
Funciones esfuerzo
En r = ∞ ⇒ x = ∞ : y = ∞
σ rr σ rθ   cos θ sin θ   σ 0 0  cos θ − sin θ 
σ σ θθ  =  − sin θ cos θ   0 0  sin θ cos θ 
 rθ  
σ rr σ rθ   cos θ sin θ   σ 0 cos θ − σ 0 sin θ 
σ σ θθ  =  − sin θ cos θ   0 0 
 rθ   
σ rr = σ 0 cos 2 θ , σ θθ = σ 0 sin 2 θ , σ rθ = − σ 0 sin θ cos θ
σ0 σ0 σ0
σ rr = + cos 2θ , σ rθ =− sin 2θ
2 2 2

Resumiendo, hay 4 CB:

σ0 σ0
σ rr = + cos 2θ en r = ∞
2 2
σ
σ rθ = − 0 sin 2θ en r = ∞
2
σ rr = 0 en r = a
σ rθ = 0 en r = a
Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 36
 Elasticidad 2D
Funciones esfuerzo
Hay 4 CB y 5 constantes. En realidad la primera CB posee dos partes por lo que
son dos CB en una. Obviando el trabajo matemático, se obtiene:

σ0 a2  σ 0  a2 a4 
σ =  1 − 2  +  1 − 4 2 + 3 4  cos 2θ
2  r  2 
rr
r r 
σ0 a2  σ 0  a4 
σ θθ =  1 + 2  −  1 + 3 4  cos 2θ
2  r  2  r 
σ0 a2 a4 
σ =−  1 + 2 2 − 3 4  sin 2θ
2 
rθ
r r 

El valor máximo se obtiene en θ = 90º y r = a

σ0 a2  σ 0  a4 
σ = σ =  1 + 2  −  1 + 3 4  ( − 1) = 3σ
2  a  2 
θθ xx 0
a 

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 37

 Elasticidad 2D
Funciones esfuerzo
Definiendo:
Esfuerzo local
Factor de concentración de esfuerzos (FCE) =
Esfuerzo nominal
FCE = 3 en agujero de una placa infinita isotrópica

Este es un concepto general. Mientras más fuerte la discontinuidad mayor es el

FCE.
El FCE depende del material. Mientras mayor es la relación Ex/Ey mayor es el FCE.
En compuestos unidireccionales puede existir FCE= 7.

Las funciones esfuerzo pueden aplicarse a materiales ortotrópicos pero es

necesario utilizar cálculo en variable compleja.

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 38

 Elasticidad 2D
Efectos ambientales (temperatura)
Los materiales se contraen y expanden cuando la temperatura cambia. Así:
ε T = α (∆ T )
ε T : deformación témica
α : coeficiente de expansión térmica (CET) unidades (1/grados)
∆ T : variación de temperatura
Si la expansión/contracción térmica es resistida por algún medio, se producen
“esfuerzos térmicos”. Este término está mal utilizado ya que lo que se produce son
“esfuerzos debido a efectos térmicos”. Los esfuerzos son siempre “mecánicos”
Un material 3-D genérico:
ε ij = α ij ( ∆ T ) i, j = 1,2,3
T

La deformación total en una estructura es la suma de la deformación mecánica y la

deformación térmica:
M T
ε ij = ε ij + ε ij
  
total mecánica térmica

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 39

 Elasticidad 2D
Efectos ambientales (temperatura)
La relación entre la deformación mecánica y los esfuerzos es:
M
ε ij = Sijkl σ kl

matriz de flexibilidad

ε ij = Sijklσ kl + α ij ∆ T ⇒ Sijklσ kl = ε ij − α ij ∆ T

Multiplicando a ambos lados por la inversa de la matriz de flexibilidad (matriz de

rigidez o tensor elasticidad):
−1
Sijkl = Eijkl ⇒ σ kl = Eijklε ij − Eijklα ij ∆ T

Que es la misma ecuación que se tenía antes más los esfuerzos debido a efectos
térmicos.

ε M = −ε T ⇒ ε total
= 0

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 40

 Elasticidad 2D
Efectos ambientales (temperatura)
Valores de CET:
Materiales anisotrópicos 6 posibilidades (como esfuerzos y deformaciones es un
tensor): α ,α ,α ,α ,α ,α α = α
11 22 33 12 13 23 ij ji

⇒ ∆ T puede causar deformaciones angulares

Esto último no es posible en materiales reales.

Materiales ortotrópicos 3 posibilidades:
α 11 ,α 22 ,α 33
⇒ ∆ T sólo causa deformaciones normales

Si el sistema de referencia es diferente de los ejes principales del material, puede

tenerse un α12. Las transformaciones operan igual que para las deformaciones (es
un tensor).
α~ij =  ~i k  ~j lα kl

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 41

 Elasticidad 2D
Efectos ambientales (temperatura)
En 2-D:
α~11 = α 11 cos 2 θ + α 22 sin 2 θ
α~22 = α 11 sin 2 θ + α 22 cos 2 θ
α~12 = (α 22 − α 11 ) cos θ sin θ

Por lo tanto en un material isotrópico no se producen componentes de cizalle (α11= α22.)

Valores Típicos de CET:

Material CET(x10-6/ºF)
Acero 6
Aluminio 12.5
Titanio 5
Uni Gr/Ep (dirección fibras) -0.2
Uni Gr/Ep (perpendicular fibras) 16

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 42

 Elasticidad 2D
Efectos ambientales (humedad)
Algunos materiales absorben humedad. Esto es caracterizado por un factor de “hinchado” o
de expansión:
S
ε ij = β ij c i, j = 1,2,3

M T S
Así: ε ij = ε ij + ε ij + ε ij
   
total mecánica térmica humedad

Esto puede ser generalizado para cualquier efecto ambiental sobre la deformación como:
E
ε ij = χ Χ
 ij escalar ambiental
deformación ambiental operador ambiental

Estructuras Aeroespaciales Unidad 3: Elasticidad 2D 43