ELASTICIDAD BIDIMENSIONAL

EN COORDENADAS POLARES
1. ESTADO PLANO DE DEFORMACIÓN
En un estado plano de deformación las componentes del vector corrimiento deben ser:

Las deformaciones se obtienen de forma inmediata particularizando las expresiones de coordenadas cilíndricas:

Quedando:

Teniendo en cuenta la Ley Hooke generalizada y que

Las tensiones queda:
1. ESTADO PLANO DE DEFORMACIÓN
Las ecuaciones de equilibrio se obtienen particularizando las ecuaciones vistas en coordenadas cilíndricas:

Si f v F , F ,0  son las fuerzas de masa, estas ecuaciones se reducen a:



Expresando estas relaciones en términos de desplazamientos dan lugar a las ecuaciones de Navier:

  u 1 v u 
F    G         G   2u  0
       u 1 v u
ya que: div     
   
1   u 1 v u 
F    G         G   2v  0
       
1. ESTADO PLANO DE DEFORMACIÓN

La ecuación del estado plano  2  x   y   

1
 div f v se convierte en:
1 


2        1
1 
 div f v

2 1  1 2
2    2 2
 2    
F 1 F F
div f v    
   


2      0 
2. ESTADO PLANO DE TENSIÓN
2. ESTADO PLANO DE TENSIÓN

Las ecuaciones de equilibrio interno son las

mismas que en el caso de deformación plana. O sea:

Expresando estas relaciones en términos de

desplazamientos dan lugar a las ecuaciones de Navier:

E   u 1 v u 
F         G   2u  0
2  1           u 1 v u
ya que: div     
E 1   u 1 v u     
F          G   2v  0
2  1           
2. ESTADO PLANO DE TENSIÓN
La ecuación del estado plano 2  x   y   1     div fv se convierte en:

2      1     div fv

2 1 2
1
  2 
2
 
    2  2
F 1 F F
div f v    
   

 
2      0
3. FUNCIÓN DE TENSIONES EN COORDENADAS
POLARES
Las relaciones entre la función de Airy y las componentes de las tensiones en coordenadas cartesianas, se
pueden transformar fácilmente en coordenadas polares:
1  1  2  2   1  
             
   2  2  2     

Puede también verificarse que estas relaciones satisfacen trivialmente a las ecuaciones de equilibrio, y en
ausencia de fuerzas de volumen, las ecuaciones de compatibilidad se traducen a:
 2 1  1  2   2 1  1 2 
 4   2    2  2  2    2  2   0
            

La forma final de la solución, conocida como la solución de Michell, se escribe como:
 
  ,   a0  a1  log   a2   2  a3   2  log   a4  a5  log   a6   2  a7   2  log   
 a 
  a11    a12    log   13  a14   3  a15      a16      log    cos
  
 b 
  b11    b12    log   13  b14   3  b15      b16      log    sin 
  

 an1  r n  an 2  r 2 n  an3  r  n  an 4  r 2 n  cosn  



n2

 bn1  r n  bn2  r 2 n  bn3  r  n  bn4  r 2 n  sinn  



n2

Donde an, anm y bnm son constantes a determinar.

Notar que esta solución es estrictamente periódica, lo que la hace muy práctica para aplicar el método de
Fourier para tratar con las condiciones de borde.
4. DISTRIBUCIÓN SIMÉTRICA DE TENSIONES CON
RESPECTO A UN EJE EN CASOS DE DEFORMACIÓN O
TESIÓN PLANAS, SIN FUERZA DE MASA

4 4

Su solución general es:

  A  log    B   2  log    C   2  D

donde A, B, C y D son constantes de integración. Luego, las tensiones están dadas por:
1 
 B  1  2  log    2  C
A
   
   2

 2
 B  3  2  log    2  C
A
  
 2 2
   0
4. DISTRIBUCIÓN SIMÉTRICA DE TENSIONES CON
RESPECTO A UN EJE EN CASOS DE DEFORMACIÓN O
TESIÓN PLANAS, SIN FUERZA DE MASA
Los corrimientos correspondientes a estas tensiones se pueden determinar por integración de las
relaciones cinemáticas. Para el caso de un estado plano de tensión:
du
d
1

  r       
E
 dv
d
1

     u            u
E

Y siendo:
 B  3  2  log   2  C
A
 
A
 B  1  2  log   2  C   
 2 2
queda:
1  1    
u    A  2  1     B    log   1     B    2  C  1        A  sin   B  cos
E  
4   
v B  A  sin   B  cos  C  
E
Las constantes A’, B’ y C’ se asocian al movimiento como cuerpo rígido.
5. APLICACIONES
i- TUBERÍA DE PARED GRUESA SOMETIDA A PRESIÓN

p2
Se considera una tubería de pared gruesa, cuya sección recta es una
corona circular de radios R1 y R2, sometida a presión interior
uniforme p1 y presión exterior, también uniforme, p2
Es evidente que por la simetría de la forma y solicitación externa
que presenta la tubería en estas condiciones, el estado elástico de la
misma será de deformación plana en el caso que los extremos
estuvieran fijos, y de tensión plana, en el caso que estuvieran libres.
En ambos casos, existe un plano director perpendicular a la tubería.
Las tensiones tienen la forma:

 B  1  2  log    2  C
A
 
2
 B  3  2  log    2  C
A
  
 2
   0

Como los corrimientos son independientes de , se demuestra que la constante B debe ser nula. O sea:
A
   2 C
 2
A
    2 C
2
5. APLICACIONES
i- TUBERÍA DE PARED GRUESA SOMETIDA A PRESIÓN
A A
Deben determinarse ahora las constantes de integración A y C de    2C     2 C
 2  2

Las condiciones de contorno son:

Que se traducen en:

Cuya solución es:

Quedando las tensiones:

Los valores máximos de ambas componentes ocurren en la pared interior de la tubería, siendo:
5. APLICACIONES
i- TUBERÍA DE PARED GRUESA SOMETIDA A PRESIÓN

A A
Con y el par de ecuaciones     2  C     2  C se llega:
2  2


Análogamente, pero utilizando ahora s se llega a que:

La variación del espesor es:

con
5. APLICACIONES
ii- HUECO PRESURIZADO EN UN MEDIO INFINITO
Partiendo de la solución para tuberías sometidas a presión interna:

du 1     p2  p1   R12  R2 2 1 p1  R12  p2  R2 2 
       1  2     
dr E  R2 2  R12 2 R2 2  R12 
Se plantea que la condición p2  0 cuando R2   . Así:

   lim 
 p1  p2  R12  R2 2  1

p1  R12  p2  R2 2
    p1 
R12
R2   R2 2  R12 2 R2 2  R12 2
p 2 0

   lim
 p1  p2  R12  R2 2  1

p1  R12  p2  R2 2    p1 
R12
R2   R2  R1
2 2
 2
R2  R1
2 2 2
p 2 0

du 1    p1  R12  1    p1  R12 
     u  
dr E   2  E   

Tanto las tensiones como las deformaciones se reducen a medida que la distancia al orificio es mayor
tendiendo a cero cuando r tiende a infinito, pero las tensiones lo hacen en una forma más rápida.
 5. APLICACIONES
iii- DISCO SOMETIDO A PRESIÓN UNIFORME p
Partiendo de la solución:
1 
 2  B  1  2  log   2  C
A
  
   2

   0   2   2  B  3  2  log   2  C
A
 
Dado que en el centro, donde  = 0, las tensiones no pueden infinitas. Esto conduce a: A  B  0
Sobre la periferia, donde  = R, se tiene que     p , lo que conduce a 2  C   p
Sustituyendo los valore obtenidos:
      p    0

Tratándose de un estado plano de tensión, se tiene que:

1
E
   E

          1    
p du
d
El corrimiento radial en cualquier punto, teniendo en cuenta que es nulo en el centro, está dado por:

p
 1       1     
p
u  
0
E E
p
 1     R
Luego la variación de radio del disco es: R  u  R   
E
 2 p   dw
Para un estado plano de tensión, se tiene también que: z  
E

      
E
dz
2p 2p
z
El corrimiento en la dirección z de cualquier punto del disco está dado por: w  dz  z
0
E E
2 p  e 2 p 
Luego la variación del espesor del disco es: e  2  wz  e   2   e
2 E 2 E
5. APLICACIONES
iv- DISCO MACIZO GIRATORIO

Para calcular el estado tensional que se induce en un

disco de radio R cuando gira con una velocidad
angular w, se considera un elemento diferencial de
espesor unitario en coordenadas polares.
Sean r y  las tensiones radial y circunferencial,
respectivamente, así como dfc la fuerza centrífuga que
actúa sobre el elemento, como consecuencia de la
rotación.
Por razones de simetría  = 0
Siendo  la densidad de material, la expresión de la
fuerza centrífuga es:

El equilibrio de fuerzas en la dirección radial resulta:

5. APLICACIONES
iv- DISCO MACIZO GIRATORIO

El equilibrio de fuerzas en la dirección lugar resulta:

Simplificando y operando queda:

Si u y v son las componentes del corrimiento, la deformaciones

unitarias radial y circunferencial están dadas respectivamente por:
5. APLICACIONES
iv- DISCO MACIZO GIRATORIO

Igualando y

queda:

Sustituyendo :

Integrando:

Donde A es una constate de integración

5. APLICACIONES
iv- DISCO MACIZO GIRATORIO
Restándole la ecuación queda:

Luego
 5. APLICACIONES
iv- DISCO MACIZO GIRATORIO
Siendo:
5. APLICACIONES
v- DISCO GIRATORIO CON ORIFICIO CENTRAL
Se parte de:

De donde:
5. APLICACIONES
v- DISCO GIRATORIO CON ORIFICIO CENTRAL
Siendo:
5. APLICACIONES
vi- CHAPA PERFORADA TRACCIONADA

Sometida a tracción uniforme t en la dirección del eje x, en virtud del principio de Saint Venant, las
tensiones en los puntos de una circunferencia concéntrica con taladro de radio r2, grande en
comparación con r1, son las mismas que existirían en la placa sin taladrar. Las expresiones de estas
tensiones en coordenadas polares tomando el centro de la placa como polo y eje x como eje polar, se
obtienen de forma inmediata utilizando la circunferencia de Mohr
5. APLICACIONES
vi- CHAPA PERFORADA TRACCIONADA
5. APLICACIONES
vi- CHAPA PERFORADA TRACCIONADA

Al estado I podemos aplicarle las solución

obtenida en el estudio elástico de la tubería de
pared gruesa expuesto anteriormente, que
verifica las ecuaciones de equilibrio y
compatibilidad

Para el cálculo de las constantes A y C haremos p1=0, p2=-t/2,

despreciando el cociente r1/r2. Así:

A t  r12 1 t t  r12 
   2  2C  
I
    1
 2  2 2 2   2 

Sustituyendo, que para el estado 1: A t  r12 1 t t  r12 

   2  2  C 
I
    1
 2  2 2 2   2 
 
I
0
 5. APLICACIONES
vi- CHAPA PERFORADA TRACCIONADA

Para obtener la solución de tensiones del estado II, se determina la función f()
 
imponiendo la condición de compatibilidad  2  II   II  0 o, lo que es lo mismo,
la condición de biarmonicidad de la función  de Airy. El laplaciando de  es:

 2 1  1  2  d 2 f 1 df 4 f 
 2           cos 2
 2    2  2  d 2  d  2 
La ecuación de compatibilidad será:
 d2 1 d 4   d 2 f 1 df 4 f 
        0
 d 2  d  2   d 2  d  2 
   
Cuyo desarrollo da:
d4 f 1 d3 f 9 d2 f 9 df
      0
d 4  d 3
 d
2 2
 3 d
Su solución es
 C 
   A  2  B   4   D   cos 2 Donde A, B, C y D son constantes a determinar
 2 
 
 5. APLICACIONES
vi- CHAPA PERFORADA TRACCIONADA

Luego, las tensiones son: 

 1  2 6 C 4 D 
 cos2 
1
   
II
    2 A  
   2  2   4  2 

 2  6 C 
 II     A  12B 2   cos2 
 2   4 
     6 C 2 D 
 
II
     2 A  12B 2    sin 2 
      4  2 
Sumando las soluciones de los estados I y II queda:
5. APLICACIONES
vi- CHAPA PERFORADA TRACCIONADA
5. APLICACIONES
vi- CHAPA PERFORADA TRACCIONADA
Reemplazando el valor de las constantes:

Los valores extremos son:

    t  2  t  3 t
2
     0t  2  t  0
5. APLICACIONES
vi- CHAPA PERFORADA TRACCIONADA
Reemplazando el valor de las constantes:
5. APLICACIONES
vii- CHAPA PERFORADA COMPRIMIDA
Teniendo en cuenta lo anterior y que:
5. APLICACIONES
viii- CHAPA PERFORADA TRACCIONADA Y COMPRIMIDA
Superponiendo ambas soluciones:

+
5. APLICACIONES
viii- CHAPA PERFORADA TRACCIONADA Y COMPRIMIDA
5. APLICACIONES
ix- CHAPA SOMETIDA A TENSIÓN TANGENCIAL

Es decir, la solución de tensiones de una placa cuadrada de grandes dimensiones en cuyo centro se
realiza un taladro y se la somete a una tensión tangencial uniforme t en sus lados, estará dada por las
ecuaciones vistas justo anteriormente, tomando como eje polar una de las diagonales. En el bore del
agujero de esta placa el valor máximo de  se presenta en los puntos de intersección del mismo
con las diagonales y su valor es cuatro veces mayor que la tensión tangencial aplicada
5. APLICACIONES
A h
x- BARRAS DE EJE CURVO
Re  1
Se considera el caso de una pieza curva de sección Ri
rectangular estrecha, de directriz circular, sometida a M R M
flexión en el plano de curvatura, mediante un par M
constante a lo largo de la pieza
En cualquier sección radial, el momento flexor es el
mismo por lo que la distribución de las tensiones O
será la misma. Por lo tanto, se utilizarán coordenadas
polares con centro en O.

Como el problema es independiente de la coordenada radial, se utiliza la solución general para

problemas axilmente simétricos. O sea:

 B  1  2  log   2  C
A
 
2
 B  3  2  log   2  C
A
  
2
   0

Aplicando la condiciones de contorno:

  Re   Re
    R   0     R   0   1 d  0     1 d  M
i e
  Ri   Ri
5. APLICACIONES A h
x- BARRAS DE EJE CURVO
Re  1
Se llega a que: Ri
M R M

 B  1  2  logRi   2  C
A
0
Ri 2

 B  1  2  logRe   2  C
A
0
Re 2 O
   
0  Re  2  B  logRe   2  B  2  C   Ri  2  B  logRi   2  B  2  C 
A A
 Re   Ri 
R 
  
 M   A  log e   B  Re2  logRe   Ri2  logRi   C  Re2  Ri2 
 Ri 
La tercera ecuación es una combinación lineal de las dos primeras. Resolviendo, queda que:

4 M 2 2 R 
A  Re  Ri  log e 
N  Ri 
B
2 M
N

 Re2  Ri2 
C
M
N
 
 Re2  Ri2  2  Re2  logRe   Ri2  logRi  

Siendo  2
 R 
N  Re 2  Ri 2  4  Re 2  Ri 2  log 2  e 
 Ri 
5. APLICACIONES A h
x- BARRAS DE EJE CURVO
Re  1
Reemplazando las constantes de integración en la Ri
M R M
solución general, queda:

4 M  2  Ri     Re 2  Ri 2  Re 
     Ri  log   Re  log  
2
 log  
N    Re  2  Ri 
  O
4 M R     R 2 R 2 R 
     Ri 2  log i   Re 2  log   e 2 i  log e   Re 2  Ri 2 
N    Re    Ri  
   0

En la figura siguiente se muestra un diagrama de  en una sección para un relación Re/Ri=4

Re/Ri=4
 Ri / M

Eje neutro

/Ri
 R   
M
  
Si la altura de la pieza es pequeña frente al radio de curvatura R distribución de J
tensiones queda: r  0
5. APLICACIONES
xi- SOLIDO SEMI-INDEFINIDO SOMETIDO A CARGA UNIFORMEMENTE
DISTRBUIDA NORMAL AL PLANO. Problema de Flamant-Boussinesq

Se plantea encontrar el tensor de tensiones en cada punto de una placa, elástica y semi-infinita,
sometida a una carga uniformemente distribuida normal al plano que lo limita.
5. APLICACIONES
xi- SOLIDO SEMI-INDEFINIDO SOMETIDO A CARGA UNIFORMEMENTE
DISTRBUIDA NORMAL AL PLANO. Problema de Flamant-Boussinesq

A la solución del problema planteado se llega teniendo en cuenta las siguientes consideraciones sobre las
tensiones que actúan sobre las caras del elemento indicado en la figura.

2  P  cos 
  
 
      0
5. APLICACIONES
xi- SOLIDO SEMI-INDEFINIDO SOMETIDO A CARGA UNIFORMEMENTE
DISTRBUIDA NORMAL AL PLANO. Problema de Flamant-Boussinesq

P     sin
Puede concebirse que esta solución deriva de la función de Airy:  

2  P  cos 
  
 

En efecto:   0

   0

Las ecuaciones de equilibrio, reducidas a , se satisfacen trivialmente

La ecuación de compatibilidad, también se verifica

También se verifica que  es biarmónica:

 2

 4

Luego, esta es la solución del problema, y es única

5. APLICACIONES
xi- SOLIDO SEMI-INDEFINIDO SOMETIDO A CARGA UNIFORMEMENTE
DISTRBUIDA NORMAL AL PLANO. Problema de Flamant-Boussinesq

Por los puntos de una circunferencia de diámetro d que pasa por el origen y tiene su centro en el eje polar,
la tensión radial es:

Utilizando la circunferencia de Mohr, pueden obtenerse las componentes del tensor de tensiones referidas
a las coordenadas cartesianas:
5. APLICACIONES
xi- SOLIDO SEMI-INDEFINIDO SOMETIDO A CARGA UNIFORMEMENTE
DISTRBUIDA NORMAL AL PLANO. Problema de Flamant-Boussinesq

En los puntos distantes a del borde

superior de la placa, haciendo
a

cos 
queda:

2 P
x    cos 4 
 a

2 P
y   cos 2   sin 2 
 a

2  P sin  cos 3 
 xy   
 a 
5. APLICACIONES
xii- SOLIDO SEMI-INDEFINIDO SOMETIDO A CARGA UNIFORMEMENTE
DISTRBUIDA NORMAL AL PLANO. Problema de Flamant-Boussinesq

Lo analizado anteriormente, puede ser

aplicado al estudio del estado de tensión
creado cuando la carga es continua,
sustituyendo en las expresiones anteriores:
  d
P p  dy  p 
cos 
En efecto:

Para la carga total: Si la carga p es constante:

5. APLICACIONES
xii - SOLIDO SEMI-INDEFINIDO SOMETIDO A CARGA UNIFORMEMENTE
DISTRBUIDA NORMAL AL PLANO. Problema de Flamant-Boussinesq

verifica las ecuaciones de equilibrio y compatibilidad

5. APLICACIONES
xiii. Placa semi-infinita sometida a una cupla en un punto de su borde

Sea (x,y) la función de Airy correspondiente a la fuerza P,

aplicada en O
Sea (x,y+a) la función de Airy correspondiente a la fuerza
P’, aplicada en O’
La función de Airy correspondiente al par de fuerzas es:

Siendo:
5. APLICACIONES
xiii- Placa semi-infinita sometida a una cupla en un punto de su borde
a

Siendo:
5. APLICACIONES
xiv- Disco comprimido diametralmente
Se investiga la solución del problema de un disco cilíndrico comprimido diametralmente. Este problema
es de especial interés dado que esta geometría se utiliza para ensayar (ASTM D4123) materiales
bituminosos y frágiles, como el hormigón, asfalto, rocas y cerámicos. Normalmente se conoce como el
ensayo de tracción indirecta o brasileño (ya que fue desarrollado allí por Lobo Carneiro). La compresión
genera un campo de tensiones de tracción sobre la mayor parte del diámetro cargado, lo que permite
determinar muy simplemente la resistencia a tracción del material.
El problema se resuelve por superposición:

La dos soluciones de Flamant proporcionan el comportamiento singular en los extremos cargados,

mientras que la carga radial quita las tensiones de tracción remanentes en los bordes del disco.
5. APLICACIONES
xiv- Disco comprimido diametralmente
Para combinar las dos soluciones de Flamant, conviene re-definir D

la angular, así:
2 P 2 P
 1    cos3 1  2     cos3  2
x   r1 x   r2
2 P 2 P
 1    cos1  sin 2 1  2     cos 2  sin 2  2
y   r1 y   r2
2 P 2 P
 1    cos12  sin 1  2     cos 2 2  sin  2
xy   r1 xy   r2
La solución general de problema de Flamant produce sólo tensiones
radiales constantes en un contorno circular de diámetro d:
2 P      0
  
 D
Para puntos del borde se cumple que: D2  r12  r22   1   2
2 P

Dado que se trata de un estado de tensión hidrostático, la

tensión resultante debe ser la misma, o sea: r1
1
2 P
 x  y    0
 D xy D

Luego la tensión radial debe compensar esos valores. O sea: r2

2

2 P
 3   3   3    3  0
 x y  D xy

P
5. APLICACIONES
xiv- Disco comprimido diametralmente
Por ello queda:

2 P 2 P 2 P
 1    cos3 1  2     cos3  2  3  
x   r1 x   r2 x  D
2 P
 1  
2 P
 cos1  sin 2 1  2   
2 P
 cos 2  sin 2  2  3  
y   r1 y   r2
y  D
2 P 2 P  3  0
 1    cos12  sin 1  2     cos 2 2  sin  2 xy
xy   r1 xy   r2

x  

2  P cos 1 cos  2 1

3

3 
 
2 P  2


 d x3 d x3  
 24  
1
 
  r1 r2 d    r1 4
r2 D
 

y 
2  P  cos1  sin 2 1 cos 2  sin 2  2 1 
     

2 P  2  x  y
d 2
 2

d  x  y2
1
 
 
  r1 r2 d    r14 4
r2 D
 

 xy 

2  P cos1  sin 1 cos 2  sin  2

2

2  2 P  2  2 4

 d  x 2  y d  x 2  y
  
 
 
   r1 r2     r14
r2
 
5. APLICACIONES
xiv- Disco comprimido diametralmente
D 
2 P  2  x
3
D x3  
1   xy  
D 2

2 P  2  x  y
 2 4

D  x 2  y


x     2   
  r14 r24 D  r14
r2
   

   
D
2 P  2  x  y 1
D 2 D  x  y2
y  
 

r14
 2
r24
 
D

r1,2  y 2  D  x
2
2

 

Sobre el eje y, donde x=0, se tiene que

r1,2  y 2  d  2 2
 

 x 0, y   
2 P 

 D 3
2
  
 
D 3
2 
1


2  P  4  D4
 

 1
y  2   
   2 d 2 2  2 d 2 2 D 
 y       
  D  D2  4 y 2 2  
    2  
 

 y 0, y   
2 P 

 
 D y
2
2

 
D y
2
2
1
 

2  P 
 2
D 2
 4 y 4 2
2
  2 D 2
        D  4 y 
2 2
 2  D D
y  2   y  2 
2 D 
  

 xy 0, y  
 d 2
2 P 
 2
 
y

d 2 2  y   0
2
  2 d 2

y 
2
  y 2  d  
2 
5. APLICACIONES
xiv- Disco comprimido diametralmente
D
2 P  2  x 
3
 
D x3 
1   xy  
D 2
2 P  2  x  y 
 2 4

D  x 2  y

x     2   
  r14 r24 D  r14
r2
   

   
D
2 P  2  x  y 1
D 2 D  x  y2
y  
 

r14
 2
r24
 
D

r1,2  y 2  D  x
2
2

 

Sobre el eje x, donde y=0, se tiene que

r1,2  D  x
2

 x  x ,0    

2 P  2  
 D x3 D x3
 2 
1


2 P  2
 
2 1
 
 
 D x
 2
4
 D x
2
4 D 


 D  2 x D  2 x D

 

 y  x ,0   
2 P 

  
D  y2
2  2
 
D  y2


1  2 P



  
 
   2 D 2 2  2 D 2 2 D    D
y  2   y  2  


 xy  x,0   0
5. APLICACIONES
xiv- Disco comprimido diametralmente

x  

D
2 P  2  x

 
3


D x3
2 
 
1
 xy  
D 2
2 P  2  x  y 
 2 4

D  x 2  y

  r14 r24 D   r14
r2 
   

   
D
2 P  2  x  y 1
D 2 D  x  y2
y  
 

r14
 2
r24
 
D

r1,2  y 2  D  x
2
2

 