Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Apuntes Mecanica Medio Continuo Tesis PDF

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 166

UNIVERSIDAD NACIONAL AUTÓNOMA

DE MÉXICO

FACULTAD DE INGENIERÍA

"APUNTES DE LA ASIGNATURA ELEMENTOS


DE MECÁNICA DEL MEDIO CONTINUO"

TESIS
PRESENTADA A LA DMSIÓN DE INGENIERÍA
CIVIL TOPOGRÁFICA Y GEODÉSICA
DE LA FACULTAD DE INGENIERÍA

POR
ESTÉVEZ CA YETANO JOSÉ RAÚL
PARA OBTENER EL TÍTULO
DE LICENCIADO EN
INGENIERÍA CIVIL

DIRIGIDA POR: DR. RIGOBERTO RIVERA CONSTANTINO

MÉXICO, D.F., JUNIO DEL 2005


~~;-:--~F;~~~~~::.- __'::l FACULTAD DE INGENIERíA
--;:::';"'~~~!...~~!-._-""'1tJ G/DCTG/SEAC/UTIT/049/05
DIRECCiÓN

-----
Señor
JOSÉ RAÚL ESTÉVEZ CAYETANO
Presente

En atención a su solicitud me es grato hacer de su conocimiento el tema que propuso el


profesor DR. RIGOBERTO RIVERA CONSTANTINO, que aprobó esta Dirección, para que lo
desarrolle usted como tesis de su examen profesional de INGENIERO CIVIL.

"APUNTES DE LA ASIGNATURA ELEMENTOS DE MECÁNICA DEL MEDIO CONTINUO"


PRÓLOGO
INTRODUCCiÓN
1. ESTADO DE ESFUERZO
11. ESTADO DE DEFORMACiÓN
111. PRINCIPIOS GENERALES DE LA MÉCANICA
IV. ELASTICIDAD LINEAL
V. TEORIAS DE FALLA Y RUPTURA
APÉNDICES
BIBLlOGRAFIA

Ruego a usted cumplir con la disposición de la Dirección General de la Administración Escolar


en el sentido de que se imprima en lugar visible de cada ejemplar de la tesis el Titulo de ésta.

Asimismo le recuerdo que la Ley de Profesiones estipula que deberá prestar servicio social
durante un tiempo minimo de seis meses como requisito para sustentar Examen Profesional.

Atentamen
"POR MI HABLARÁ EL EspIRITU"

'i., '¡L .
15 In lo f:>

~.
, / .• r1~ G.
AGRADECIMIENTOS

A la Universidad Nacional Autónoma de México por haberme


permitido mi formación profesional.

A la Facultad de Ingeniería: En especial a todos mis Profesores


con respeto y dedicatoria.

Al Dr. Rigoberto Rivera; mi Asesor: Por el apoyo, esfuerzo


y sinceridad que me brindó para lograr los objetivos que me
propongo.

A Juana; mi Madre: Con profundo cariño y como tributo


a sus incontables sacrificios para lograr uno de mis grandes
anhelos.

A Irene, Dolores, Arturo y José; mis Hermanos: Que son la


motivación de mi vida.

A mis Amigos, Compañeros, por hacer mis días mas agradables


ya todas aquellas personas que me brindaron su cariño, amistad,
comprensión, apoyo, tolerancia y conocimientos necesarios para
llevar a cabo este trabajo.
11
,
Indice General

Prólogo Xl

Introducción Xln

1 Estado de Esfuerzo 1
1.1 Teoría del estado de esfuerzo . . . . . . . . . . . . . . . . . 2
1.2 Componentes normal y tangencial del vector esfuerzo. 3
1.3 Determinación de las ecuaciones de Cauchy 5
1.4 Esfuerzos principales . . . . . . 10
1.5 Elipsoide de Lamé . . . . . . . . . 17
1.6 Solución gráfica de Mohr . .. .. . 18
1. 7 Procedimiento gráfico de Mohr. . . 24
1.8 Casos particulares de estados de esfuerzo 25
1.9 Descomposición de un estado general de esfuerzos en
sus componentes volumétrica y desviadora 28
1.10 Estado de esfuerzo plano . . . . . . . . . . . . . . . 30
1.11 Planos y ejes principales de esfuerzo. . . . . . . . . 33
1.12 Esfuerzos octaédricos . . . 37
1.13 Equilibrio de partículas en ·u~ . nied.i~ . c"o~tin~¿ 39
1.14 Ecuaciones de equilibrio de momentos en partículas de
un continuo . . . . . . . . 45
2 Estado de Deformación 49
2.1 Características de rotación y cambio de forma de la partícula. 52
2.2 Interpretación física del tensor rotacional . . . . . . . . . . .. 54
2.3 Interpretación física del tensor deformación . . . . . . 57
2.4 Componentes normal y cortante del vector deformación. . 63
2.5 Deformaciones principales y direcciones principales
de deformación . . . . . . . . . . . . . . . . . . . . . . . . . . 64
2.6 Significado físico del invariante lineal JI . . . . . . . . . . . . . 65
2.7 Componentes volumétrica y desviadora del tensor deformación 66
2.8 Ecuaciones de compatibilidad de deformaciones. 67

3 Principios Generales de la Mecánica 73


3.1 Princi pio de la conservación de la masa o ecuación de continuidad 74
3.2 Principio de conservación de la cantidad de movimiento . . .. 76

III
IV ÍNDICE GENERAL

3.3 Principio de conservaclOn de energía; pnmera ley de


la termodinámica . . . . . . . 78
3.4 Segunda ley de la termodinámica, de¿igu81~f¡id d~
Clausius-Duhem . . . . . . . . 81
4 Elasticidad Lineal 83
4.1 Planteamiento matemático para definir las relaciones
constitutivas en un continuo cualquiera . . . . . . . . . .. 83
4.2 Ecuaciones constitutivas de los materiales elásticos lineales
en un marco de referencia principal . . . . . . . . . . . . .. 88
4.3 Ecuaciones constitutivas de los materiales elásticos lineales
homogéneos e isótropos en un marco de referencia cartesiano. 89
4.4 Solución de problemas elásticos haciendo uso de funciones de
esfuerzo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
4.5 Función de esfuerzos de Airy . . . . . . . . . . . . . . . . . . . 97
4.6 Funciones de esfuerzos en coordenadas cilíndricas . . . . . 105
4.7 Energía de deformación elástica para un estado uniaxial de
esfuerzos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ll3
4.8 Energía de deformación elástica para un estado triaxial de
esfuerzos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ll5
5 Teorías de Falla y Ruptura 119
5.1 Teoría de Rankine . . . . . ll9
5.2 Teoría de Coulomb-Tresca 121
5.3 Teoría de Saint Venant . . 124
5.4 Teoría de Nadai . . . . . . . . . . . 125
5.5 Teoría de Von Mises Hencky(VMH) 127

A Apéndice A: Análisis Tensorial 139

B Apéndice B: Rotación de Ejes Coordenados 149

Bibliografía 153
,
Indice de Figuras

1.1: Cuerpo deformable sometido a un sistema de cargas cualesquiera. 2


1.2: Definición de vector esfuerw en el entorno a un punto de un
cuerpo deformable. . . . . . . . . . . . . . . . . . . . . . . .. 3
1.3: Componentes normal y tangencial del vector esfuerw In. . . .. 4
1.4: Componentes del vector esfuerzo cortante ;:¡:n en un sistema de
referencia cartesiano. . . . . . . . . . . . . 5
1.5: Determinación de las ecuaciones de Cauchy. 6
1.6: Equilibrio estático del tetraedro de Cauchy. . . . . . 7
1. 7: Vector esfuerzo asociado a un plano. . . . . . . . . . 8
1.8: Esfuerzos principales. . . . . . . . . . . . . . . . . . 11
1.9: Vector esfuerzo en un marco de referencia principal. 17
1.10: Elipsoide de Lamé. . . . . 18
1.11: Plano de Mohr. . . . . . . . . . . . . . . . . . . . . 19
1.12: Región de Mohr. . . . . . . . . . . . . . . . . . . . 23
1.13: Solución gráfica de Mohr. . . . . . . . . . . . . . . 25
1.14: Representación gráfica de un estado de esfuerzo nulo. 26
1.15: Estado de esfuerw uniaxial, espacio de Lamé. . . . 26
1.16: Estado de esfuerw uniaxial, plano de Mohr. . . . . 26
1.17: Estado de esfuerw plano, espacio de Lamé. 27
1.18: Estado de esfuerw plano, región de Mohr, ()3 = o. . 27
1.19: Estado de esfuerw hidrostático. . . . . . . . . . . . . . . . .. 28
1.20: Estado de esfuerzo plano, en un sistema de referencia cartesiano. 30
1.21: Derivación de los esfuerzos ()n YT n para un estado de esfuerzos
plano. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
1.22: Relación geométrica entre ()xx, ()yy, T xy Y (j . . . . . . . . . 34
1.23: Representación del estado de esfuerzo para el ejemplo 2.5. 35
1.24: Solución gráfica de Mohr, ejemplo 2.5. . . . . . . . . . . . 37
1.25: Definición de esfuerws octaédricos. . . . . . . . . . . . . . 38
1.26: Equilibrio estático de una partícula interior del continuo. . 40
1.27: Fuerzas de cuerpo y de superficie en una partícula elemental
de un medio continuo. .. . . . . . .. . . . . . . . . . . .. . 40
1.28: Viga de sección prismática sometida a flexión simple. . . . .. 42
1.29: Distribución de esfuerws cortantes en el plano x = o. . . . . . 43
1.30: Distribución de esfuerws normales en el plano x = 1. . . . .. 44
1.31: Ecuaciones de equilibrio de una partícula de un medio continuo. 45
2.1: Definición de vector desplazamiento. . . . . . . . . . . . . . .. 49

v
VI ÍNDICE DE FIGURAS

2.2: Definición de vector desplazamiento relativo. 51


2.3: Interpretación física del término TotS x e. . 56
2.4: Rotación pura de una barra rígida. . . . . . 56
" f "ISlca dI
2.5: I nterpretaclOn e térmmo 8v
. -8u + -. 58
8y 8x
2.6: Barra cilíndrica sometida a un campo de desplazamientos. 60
2.7: Componentes del vector desplazamiento en la sección z = l. 62
2.8: Componentes normal y cortante del vector desplazamiento
relativo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
2.9: Cambio volumétrico de una partícula de un medio continuo. 66
3.1: Volumen de referencia. . . . . . . . . . . . . . . . . . . . . 74
3.2: Conservación de la cantidad de movimiento. . . . . . . . . . 77
4.1: Relaciones entre esfuerzos y deformaciones en un medio continuo
elástico lineal. . . . . . . . . . . . . . . . . . . . . . . . . . 83
4.2: Interpretación física de las constantes E y /l. • • • . • . . . 95
4.3: Barra prismática sometida a un estado de esfuerzo uniaxial. 98
4.4: Barra prismática sometida a un estado de esfuerzos biaxial. 99
4.5: Estado de cortante puro. .. . . . . . . . . . . . . . . . . . 99
4.6: Barra prismática sometida a flexión pura. . . . . . . . . .. 102
4.7: Volumen de esfuerzos normales en x = O. . . . . . . . . .. 103
4.8: Ménsula sometida a una carga puntual en su extremo libre. 105
4.9: Estado radial simple. . . . .. . . . . . . . . . . . . . . . . 108
4.10: Variación de O"rr para () = cte y r variable. . . . . . . .. 109
4.11: Lugar geométrico del semiespacio con el mismo valor de O"rr. 110
4.12: Isobaras para el estado radial simple. . . . . . . . . . . . . . . 111
4.13: Roseta de deformación en un punto de un cuerpo deformable. 111
4.14: Definición de energía de deformación elástica. . . . . . . . . . 113
4.15: Definición de módulo de Resiliencia. . . . . . . . . . . . . . . 115
4.16: Energía de deformación elástica para un estado general de
esfuerzos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
5.1: Volumen de fluencia. Teoría de Rankine. . . . . . . . . . . . . 120
5.2: Superficie de fluencia para el estado de esfuerzo plano. Teoría
de Rankine. . . . . . . . . . . . . . . . . . . . . . . . 121
5.3: Superficie de fluencia para el estado de esfuerzos plano. T~o~í~ .
de Coulomb-Tresca. . . . . . . . . . . . . . . . 122
5.4: Estado de esfuerzo plano; tensión. .. . . . . . 123
5.5: Estado de esfuerzo plano; compresión. . . . . . . . . . 123
5.6: Superficie de fluencia, criterio de Saint Venant. . . . . 125
5.7: Superficie de fluencia; teoría de Nadai. . . . . . . . . . 127
5.8: Barra prismática sometida a una carga puntual en su extremo
libre. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 130
5.9: Viga sometida a flexión pura. . . . . . . . . . . . . . . . . .. 132
5.10: Relación Mp/ MI para diferentes secciones. . . . . . . . . .. 134
5.11: Muestra cilíndrica en un molde rígido, a) alzado, b) planta. 136
ÍNDICE DE FIGURAS VII

A.1: Coordenadas de un punto "P" en un sistema de referencia


cartesiano. . . . . . . . . . . . . 140
A.2: Sistema de referencia cartesiano. . . . . . . . . . . . . . . . . . 141
B.1: Rotación de ejes coordenados. . . . . . . . . . . . . . . . . .. 149
8.2: a) Transformación de un sistema de coordenadas x], X2, X3 en
x;, x; y x~ . .... ....... .. . ... .. ......... 151
'"
Indice de Tablas
A.l : Transformación de coordenadas . . . . . 139
8.1: Tabla de deformación de coordenadas .. 151

lX
Prólogo

La elaboración de la Tesis "Apuntes de la Asignatura Elementos de Mecánica


del Medio . Continuo", se preparó arduamente al poner en forma el material
preparado. .

Por otra parte la presentación del material se hace en notación tensorial,


como es común en muchas universidades del país, sin embargo, esto dificulta
el aprendizaje de dicha asignatura, sobre todo cuando se trata de un primer
curso a nivel licenciatura.

El objetivo de esta Tesis es de contar con un material que cubra el temario


actual de la asignatura Elementos de Mecánica del Medio Continuo que forma
parte del programa de estudios de la carrera de Ingeniero Civil que se imparte
en la Facultad de Ingeniería de la UNAM.

Los Elementos de Mecánica del Medio Continuo es una asignatura básica


en la carrera de Ingeniero Civil, pues la hipótesis de medio continuo que se
adopta en esta mecánica permite abordar, sin ninguna distinción, problemas
de sólidos, líquidos y gases.

Adoptando la hipótesis de medio continuo se establecen los conceptos de


estado de esfuerzo y estado de deformación y mediante las ecuaciones
constitutivas correspondientes se hace la liga entre esfuerzos y deformaciones,
dando lugar a las diferentes teorías que utiliza el Ingeniero Civil en el diseño
de sus obras, tal es el caso de las teorías de la elasticidad, plasticidad,
viscosidad y combinaciones de éstas.

Esta Tesis, tiene como finalidad poner al alcance de los alumnos de la carrera
de Ingeniería Civil, de la Facultad de Ingeniería de la UNAM, los primeros
conceptos de esta herramienta que les servirán como base para cursos
posteriores que tienen que ver con la ciencia de la ingeniería y la ingeniería
aplicada.

Xl
Introducción

La mecánica del medio continuo es una rama de la física que estudia el estado
de esfuerzo y deformación (o flujo) de sólidos, líquidos y gases.

Esta mecánica no toma en cuenta la estructura molecular de la materia,


sólo la considera como carente de vacíos, huecos o discontinuidades. Se da
por hecho que la materia se halla distribuida en forma continua en todo su
volumen, llenando por completo el espacio que ocupa.

Desde un punto de vista físico un medio continuo es cualquier cuerpo


deformable que se estudia macroscópicamente, en el cual cada partícula del
cuerpo siempre permanece en contacto con sus vecinas.

Este concepto de medio continuo constituye el postulado fundamental de la


mecánica del medio continuo, lo que permite estudiar de manera semejante
el comportamiento de sólidos, líquidos y gases.

La presentación del material se hace sobre la base del análisis vectorial, sin
embargo; muchas de las ecuaciones que surgen en esta mecánica se expresan
tambien en notación indicial, con el objeto de que el alumno se familiarice
con este enfoque, que es el que se emplea en la literatura técnica relacionada
con esta rama de la física.

Xlll
xiv o. INTRODUCCiÓN

El concepto de estado de esfuerzo (estát.ica del medio cont.inuo) se establece en


el capít.ulo 1, mientras que el concepto de estado de deformación (cinemática
del medio continuo) se trat.a en el capítulo 2, para lo cual sólo fue necesario
adopt.ar la hipótesis de medio continuo. En estos dos capítulos se demuestra
que la estructura matemática para estudiar el estado de esfuerzo y el estado
de deformación es enterament.e similar. En el capítulo 3 se estudian las
leyes generales de la mecánica, como son los principios de conservación de
la masa, de la rapidez de la cantidad de movimiento, de la conservación de
energía o primera ley de la termodinámica y el de aumento de entropía o
desigualdad de Clausius-Duhem, los cuales son aplicables a cualquier medio
cont.inuo. Para ligar los esfuerzos con l~ deformaciones es necesario hablar
de las propiedades del material, lo cual se hace a través de lo que se conoce
como ecuaciones constitutivas de los materiales. En este primer curso sólo
se aborda, en el capítulo 4, la teoría de los materiales elásticos lineales
homogéneos e isótropos. La solución de ciertos problemas elásticos se hace
mediante funciones de esfuerzo, por considerar que estos métodos de solución
siguen siendo muy didácticos, sin embargo, no debernos olvidar que en la
actualidad la solución de problemas complejos se lleva a cabo mediante
herramientas numéricas muy poderosas, como es el caso de los métodos
de diferencias finitas, elemento finito y otros. Finalmente en el capítulo
5, se tratan las teorías de falla y ruptura de los materiales, las cuales están
Íntimamente ligadas con el diseño de las obras civiles que proyecta el ingeniero.
Capítulo 1
Estado de Esfuerzo

Los tipos de fuerzas que se aceptan en el estudio de los continuos son:

a) Fuerzas de superficie, aplicadas en las fronteras del continuo por la acción


de otros cuerpos que se encuentran en contacto con el medio.

La fuerza resultante de todas las fuerzas de superficie que actuan sobre un


área A, de un medio continuo está dada por:

1 In dA = ik 1
tk dA (1.1)

El subíndice k que aparece en la ecuación 1.1 se usa para expresar el mismo


concepto en notación indicial.
donde:
In : vector fuerza de superficie
dA: elemento diferencial de área
Adoptando un sistema de referencia cartesiano la ecuación 1.1 queda corno:

(1.2)

Las unidades de t x, t y Y tz son [F , L -2]

b) Fuerzas de cuerpo, provocadas por la acción de cuerpos distantes que


generan campos gravitacionales, de temperatura y electromagnéticos.

La fuerza resultante de todas las fuerzas de cuerpo actuantes en un volumen


finito dV, está dada por:

(1.3)

1
2 1. ESTADO DE ESFUERZO

doude:
7 : vector fuerza de cuerpo
dV: elemento diferencial de volumen
p : densidad del medio

En notación cart.esiana la ecuación 1.3 se puede escribir como;

Las unidades de fx, fy y fz son [F, M- 1 ] .

1.1 Teoría del estado de esfuerzo

Cuando un cuerpo deformable se somete a solicitaciones de cualquier tipo,


éste se deforma hasta cierto límite. Esto se debe a que las fuerzas cohesivas .
han entrado en juego, tomando un valor tal que permiten equilibr¡u- a las
fuerzas externas aplicadas.

Para describir las acciones entre todas las partículas de un continuo,


imaginemos un sistema de fuerzas aplicadas a un medio continuo, tal como
se ilustra en la Fig. 1.1.

\ Plano de corte i,!!aginario


M1v-./cuya normal es n.
Figura 1.1: Cuerpo deformable sometido a un sistema de cargas
cualesquiera.

Al hacer un corte imaginario a través de un plano cualesquiera, cuya normal


está definida por el vector unitario Ti, se obtienen los cuerpos 1 y II, mostrados
en la Fig. 1.2.
1.2. COMPONENTES NORMAL Y TANGENCIAL DEL VECTOR ESFUERZO 3

M1
Figura 1.2: Definición de vector esfuerzo en el entorno a un punto de un
cuerpo defonnable.

Consideremos un punto "P" de dicha sección de corte y tomemos un entorno


de área del mismo, óA, y sea óF la parte de la acción que el cuerpo II
comunica al cuerpo 1, únicamente a través de óA.

DEFINICIÓN: El vector esfuerzo, In, en un punto "P", asociado a un plano


de corte cualesquiera, está dado por:

- lím ÓF
(1.5)
tn = óA -t üóA

1.2 Componentes normal ' y tangencial del


vector esfuerzo

El vector esfuerzo In puede descomponerse en un vector esfuerzo normal Un


y en un vector esfuerzo cortante T" tal como se muestra en la Fig. 1.3. De
esta manera se obtiene:

In = Un + Tn
4 1. ESTADO DE ESFUERZO

Figura 1.3: Componentes normal y tangencial del vector esfuerzo I,..

El módulo del vector esfuerzo normal (fn, se puede calcular como la proyección
del vector esfuerzo In sobre la dirección n.

(1.6)
El vector esfuerzo cortante T n se puede calcular como la diferencia vectorial:

(1. 7)
La normal n que define el plano donde actua tn, se puede expresar en función
de sus cosenos directores como:

n = cos a + cos (3 + cos,


Siendo a, (3 y ,los ángulos que forman los ejes de referencia cartesianos con
la dirección de dicha normal.
Llamando a:

cosa = 1Lx cos(3 = 1Ly y cas, = 1Lz


La normal n se puede expresar como:

El vector esfuerzo (fn es positivo si tiene la misma dirección que n; siempre


será normal al plano de corte.

El vector esfuerzo cortante T n siempre estará alojado en el plano de corte,


adoptando cualquier dirección.

En un marco de referencia cartesiano, es usual descomponer a Tn en dos


componentes paralelas a los ejes del marco de referencia.
1.3. DETERMINACiÓN DE LAS ECUACIONES DE CAUCHY 5

Figura 1.4: Componentes del vector esfuerzo cortante T n en un sistema de


referencia cartesiano.

Por lo tanto el vector esfuerzo In se puede expresar como:

(1.8)

En esta ecuación, el primer subíndice indica la dirección de la normal al plano


donde actua dicho esfuerzo, mientras que el segundo indica la dirección del
eje cartesiano al cual es paralelo.

En caras positivas las componentes de esfuerzos serán positivas, cuando


tengan la dirección de los ejes del marco de referencia.

En cara..c; negativas los componentes de esfuerzo serán positivas, cuando sean


de dirección contraria a la de los ejes del marco de referencia. En lo que
respecta (f n, éste será de tensión si tiene signo positivo y de compresión si
tiene signo negativo.

1.3 Determinación de las ecuaCIones de


.
Cauchy

Por lo expuesto anteriormente podemos decir que para describir las acciones
que se generan en un plano de corte de un continuo, será necesario describir
al conjunto de vectores In, asociados a todos los puntos del plano de corte.
6 1. ESTADO DE ESFUERZO

Ya que en un punt.o "P" del continuo ~e pueden establecer una infinidad de


planos de corte, se puede afirmar que en cada punto del continuo se debe
est.ablecer a la infinidad de vectores In, asociados a todos los planos de corte
que pasan por "P". La envolvente de todos los vectores In debe ser una
función continua.

Para describir a esa función eauchy propuso el siguiente procedimiento,


basado en consideraciones de equilibrio.

Consideremos un punto "P" de un medio continuo y supongamos conocido


el esfuerzo en dicho punto, asociado a tres planos mutuamente ortogonales,
mismos que tomaremos como coordenados. Supongamos ahora que se desea
calcular el esfuerzo en "P", asociado al plano J que pasa por él y que está
definido mediante la normal n.

Para encontrar el esfuerzo citado se procede de la siguiente manera (Fig. 1.5).

z
z
8'
n

y ---~-~y

a) b)
Figura 1.5: Determinación de las ecuaciones de eauchy.

Se elige un plano J' paralelo a J y separado una distancia "h" del origen. En
estas condiciones analicemos el equilibrio estático del tetraedro así obtenido,
utilizando como esfuerzo actuante en la cara e, B del tetraedro, al medio de
los realmente actuantes en dicha cara. Al tomar límite cuando h --. 0, se
obtiene:

a) El esfuerzo medio In, se reduce al realmente actuante en "P".

b) El plano J' coincidirá con 8.

Analicemos el equilibrio del tetraedro que se muestra en la Fig. 1.6. En este


análisis no se est.án incluyendo fuerzas de cuerpo, causadas por efectos de
campos inerciales, gravitacionales, magnéticos o eléctricos.
1.3. DETERMINACION DE LAS ECUACIONES ... 7

t nz

au
I
y..
Figura 1.6: Equilibrio estático del tetraedro de Cauchy.

Definiendo:
AABC = A = área comprendida entre ABe.
p= masa específica del medio; [M, L -3].

El área ABe, al proyectarse en cada uno de los ejes coordenados queda:

Acosa = Anx
A cos¡3 = Any
Acos, = Anz

Estableciendo el equilibrio del tetraedro, se tiene.


Por L Fx = O se obtiene:

Por lo tanto,

Tomando límite, cuando h ---t O, se obtiene:


8 1. ESTADO DE ESFUERZO

-O"xxAnx - 'TyxAny - 'TzxAn z + txA = O


t x = O"xxnx + 'Tyxny + 'Tzxnz (1.10)

Análogamente se puede obt.ener por ¿ Fy = O Y ¿ Fz = O:

'Txynx + O"yyny + 'Tzynz (1.11)


'Txznx + 'TyzTLy + O"zzn z ( 1.12)

Estas últ.imas tres ecuaciones pueden ser escritas en forma matricial corno
sigue:

(1.13)

donde:
Tij = Tensor esfuerzo de orden 2 (Apéndice A).

Físicamente los elementos del tensor esfuerzo representan los esfuerzos


actuantes en tres planos mutuamente ortogonales. Matemáticamente el
tensor esfuerzo es un operador que al actuar sobre n, produce el vector In.

Ejemplo 1.1. En un punto P de un medio continuo el tensor esfuerzo es:

14 7 -7]
Tij = 7 21 O
[ -7 en MPa
O 35

Determine el vector esfuerzo en un plano que contiene al punto P y es paralelo


a los planos a) BGE, b) BGFC, del paralelepípedo de la Fig. 1.7.

6cm

4cm I

F
#'
X
Figura 1. 7: Vector esfuerzo asociado a un plano.
1.3. DETERMINACION DE LAS ECUACIONES ... 9

Dado que

a) Es necesario calcular la normal al plano BG E, para lo cual procedemos


como sigue:

Coordenadas del punto B (0,0,4) BG =0,= 2I -4k


Coordenadas del punto G (2, 0, O) BE = b= 6} -4k
E (0,6, O)

El producto cruz a x {j resulta:

ax b [~ 6 !4] = I (24) - 2 (-4) - 6k)


=
°6 -4 .
ax b 24I + 8} + 12k

la x {jI J24
2 2
+ 8 + 12
2
= 28

Por lo tanto, la normal al plano BG E es:

rt
a x b = -1
= --_- (- -
24i + 8j
6- 2- 3-
+ 12k-) = -i + -j + -k
la x bl 28 7 7 7

Operando para calcular el vector esfuerzo normal, tenemos:

ll11
[ -7 (}][~j+]
° 35 3/7
nI + 12) + 9k M Pa
b) Para el plano BGFC, procedemos de manera similar.

Punto B (0,0,4) a = 2I -
BG = 4k
Punto G (2, 0, O) Be = e = 6}
Punto e (O, 6, 4)
10 1. ESTADO DE ESFUERZO

[o~ ~ !4] = I
6 O
(+24) - 2(-6k)
axe 24I + 12k
la x el J720 = 12vÍS

La normal al plano BGFC resulta:

ax e 1 (- -)
Tí la x el = vÍS 2i +k

In [1.(
-7
dIo el]
35
[ 2~1/5v's]
tn = ~ (21I + 14] + 21k) MPa

Ejemplo 1.2. Determinar las componentes normal y cortante del vector


esfuerzo tn en el plano BGFC del ejercicio 1.1.

In ·Tí
1 1
v's (21, 14, 21) . (2 , 0,1) v's
1 63
(J" n = 5 (42 + 21) = 5
63 MPa
5

Tn [rq2 - (J"!] ~
ITnl 2 ~ (21 2 + 142 + 212) = 215.60
1
Tn [215.60 - 158.76]2 = 7.54
T" 7.54 MPa

1.4 Esfuerzos principales


Para determinar los valores máximo y mínimo del conjunto de vectores In
que definen el estado de esfuerzos en el punto, es necesario encontrar el plano
que pasa por el punto cuya normal coincida con la dirección del vector In Fig.
1.8.
1.4. ESFUERZOS PRINCIPALES 11

x
Figura 1.8: Esfuerzos principales.

Si la normal al plano coincide con la dirección del vector esfuerzo In, entonces
T n =0, por lo que:

Por otra parte, el vector esfuerzo normal se puede expresar como:


Un = Iunl n = un
Por lo tanto

desarrollando:

uxxnx + 7 yxny + 7 zxn z (1.14)


7 xynx + Uyyny + 7 zy n z (1.15)
7 xz n x + 7 yzny + uzzn z (1.16)

En este sistema de ecuaciones las incógnitas son: u, n x , ny, n z •


Dado que se tienen tres ecuaciones y cuatro incógnitas, es necesario introducir
una ecuación adicional, ella es:
n x2 +n2y + n z2 = 1 (1.17)
Las ecuaciones 1.14, 1.15 Y 1.16 se pueden representar como:

(u xx - u)nx + 7 yx ny + 7 zx n z o (1.18)
7 xynx + (uyy - u)ny + 7 zy n z O (1.19)
7 xz n x + 7 y;ny + (u zz - u)nz O (1.20)
12 1. ESTADO DE ESFUERZO

La..') ecuaciones 1.18, 1.19 Y 1. 20, t.ienen corno modelo matemático:


(1.21)

Siendo:

<T valor característico


n vector característico

Es decir, la solución de lru:; ecuaciones 1.18, 1.19 Y 1.20 es un problema de


valores y vectores característicos, por lo tanto, la expresión 1.21 puede ser
escrita como: .
(1.22)

siendo 1 = matriz identidad


Para que exista solución diferente a la trivial para n x , ny, nz(n =1= O), es
necesario que el determinante de la matriz [7ij - <T 1] sea igual a cero. Esto
es:

<Txx - <T
T xy
T yx
<T yy - <T
T zx
T zy
I= O
I T xz T yz <Tzz - <T

Desarrollando el determinante se obtiene la siguiente ecuación característica:

1I <T 2 + 12 <T - 13 = O
<T
3
- (1.23)
Donde los coeficientes de la eco 1.23 están dados por:

+ 1 T=y ~x 1;
1
<T xx
T yx
T xy 1
<Ty y
<Tyy T y=
<T zz 1+1 xx
invariante cuadrático
I[Tij ] I; invariante cúbico

Resolviendo la ecuaClOn característica 1.23 se obtienen los esfuerzos


principales <TI, <T2 Y <T3 ·

Sustituyendo los valores de <TI, <T2 Y <T3 en las ecuaciones 1.18, 1.19 y 1.20, en
forma sucesiva, y haciendo uso de la ident.idad fundamental n; + n~ + n~ =
1, más dos de estas ecuaciones así obtenidas, se determinan los vectores
característicos respectivos, nI, n2, y n3'

Se puede demostrar que las raíces de la ecuación cúbica son siempre reales,
si el tensor T ij es simétrico.
1.4. ESFUERZOS PRINCIPALES 13

Las raíces de la ecuación cúbica reciben el nombre de esfuerzos principales,


siendo (T1, el mayor, (T2, el intermedio y (T3 el menor, por lo que:

Donde los subíndices 1, 2,3 nos indican tres ejes ortogonales que pasan por
el punto, a los cuales se les pueden asociar tres planos ortogonales en los que
únicamente existe esfuerzo normal.
Cuando el tensor [Tij ] es un tensor simétrico es posible demostrar que los
productos escalares n1 . n2 = O; n2 . n3 = O; n3 . n1 = 0, lo cual indica la
ortogonalidad de las tres direcciones principales de esfuerzo.
Ya que los planos principales son ortogonales entre sí, sus normales pueden
ser utilizadas como ejes principales de referencia, el tensor referido a ese
marco, queda representado por:

Ejemplo 1.3. En un punto de un medio continuo se establece el siguiente


tensor de esfuerzo.

Obtener:
a) El vector esfuerzo en el plano que pasa por el punto cuya normal queda
definida por:

b) Los esfuerzos principales (T}, (T2, (T3

c) Las direcciones de los ejes principales definidas por nI, n2, n3


d) Los productos escalares nI . n2, n2 . n3 y n3 . n]
a) El vector esfuerzo se calcula como:
14 1. ESTADO DE ESFUERZO

['Ji;] Inl
[(4) (~) + (3) (-~) + (O) (~)] ~
+ [(3) (~) + (-4) (-~) + (O) (~)] J
+ [(O) (~) + (O) (-~) + (2) (~)] k
2- 14- 2-
-i+-j+-k MPa
3 3 3

El esfuerzo normal quedará definido por:

El esfuerzo cortante se definirá como:

b) Los e>fuerzos principale> (JI , (J2, (J3 , se calculan resolviendo la ecuación


característica.
De>arrollando:
1.4. ESFUERZOS PRINCIPALES 15

[
4 "3 u
3
_4_ u g ]= O
O 2-uO
4 - u [( - 4 - u) (2 - u)]- 3 [(3) (2 - u)] = O
4 - U (u 2 + 2u - 8) - 3 (6 - 3u) = O
-32 + 16u + 2u 2 - u 3 - 18 + 9u = O
-50 + 25u + 2u 2 - u 3 = O
u 3 - 2u 2 - 25u + 50 = O; Ecuación característica
Resolviendo la ecuación característica, se obtiene:
u} = 5 MPa
U2 = 2 MPa
U3 = -5 MPa
Por lo tanto, la matriz de esfuerzos principales queda:

5 O O]
T... =
IJ [ OO O2 -5O [MPa]

Dado que los invariantes del tensor esfuerzo en los sistemas de referencia
cartesiano y principal son iguales, entonces dichos tensores representan el
mismo estado de esfuerzos en el punto "P" del continuo.
c) Las direcciones de los ejes principales definidas por nI> n2, n3, se calculan
como:
Para u} = 5 [M Pa], tenemos:
nn + 3n21 + On31
3n n 9n 21 + On31
Onn + On21 + (-3) n31 O
Pero, sabemos que

(nn)2 + (n21)2 + (n3d 2 = 1


Por lo tanto, sustituyendo

(3n21)2 + (n2d 2 + (0)2 = 1


9n~} + n~l = 1
10n~1 = 1
n21 = 1jJ1O , nll = 3jJ1O , n31 = O

Para U2 = 2, [M Pa]
2n}2 + 3n22 + O O
3n12 6n22 + O O
O + O + O O
16 1. ESTADO DE ESFUERZO

. 2 2 2
HacH~nd() (nJ2) + (n22) + ('/1,32) =1
sllstituyemlo, nos queda:

(n32)2 = 1
:. 1L32 = 1
1LJ2 = O , 1L22 = O , n32 =1

Para 0"3 = -5 [M Pa], tenemos:

9nJ3 + 3n23 + O = O
3nOI3 ++ n,23 ++ O O
U 7n33 O

Haciendo (n13)2 + (n23)2 + (n33)2 = 1


sustituyendo:

(n13)2 + (-3n13)2 = 1
n~3 + 9n~3 = 1
10n~3 = 1
n13 = l/JiO , n23 = -3/JiO , n33 = O
Por lo tanto, la matriz de cosenos directores de los ejes principales, queda:

d) Productos escalares

Con lo cual se concluye que las direcciones nI , n2 y n3 son ortogonales.


1.5. ELIPSOIDE DE LAMÉ 17

1.5 Elipsoide de Lamé


Se dice que dos tensores representan el mismo estado de esfuerzos en un
punto cuando sus invariantes son ig;uales.
Cuando se seleccionan a los ejes principales como marco de referencia, en un
plano cualesquiera que pase por el punto, definido por 11 = nlt + n2J + n3k;
las expresiones que definen a las componentes de In en ese plano, paralelas a
los ejes 1,2 y 3, están dados por (Fig. 1.9).

1
Figura 1.9: Vector esfuerzo en un marco de referencia principal.

In' 111 = (TI nI (1.24)


In' 112 (T2 n 2 (1.25)
In' 113 (T3 n 3 (1.26)

Elevando al cuadrado las ecs. 1.24, 1.25 Y 1.26, despejando ni, n~ y n~, y
sumando miembro a miembro se obtiene:

etn '11d 2 (Tl n l)2


- 2
(t n . 112) ( (T2 n 2)2

(In ' 113)2 ( (T3 n 3)2

(In' 111)2 2 t;,¡


(T2
nI =-2 (1.27)
1 (TI

(In' 112)2 2 t;'2


n
2- -
- (T2 (1.28)
(T2
2 2
(In' 113)2 2 t;3
n3 =2 (1.29)
(T2 (T3
3
18 1. ESTADO DE ESFUERZO

Sumando 1.27, 1.28 Y 1.29 se llega a:

( 1.30)

El conjunto de vectores In que describe el estado de esfuerzo en un punto,


queda envuelto por un elipsoide (Lamé), tal corno se muestra en la Fig. 1.10.
Por lo tanto, el máximo valor de In debe ser UI Y el mínimo U3'

Figura 1.10: Elipsoide de Lamé.

1.6 Solución gráfica de Mohr


Puesto que conocido el tensor esfuerzo y elegido un vector n,
podernos determinar a la pareja de valores (un, 7 n), mediante:

(1.31 )

(1.32)

Al investigador Otto Mohr se le ocurrió inventar el plano coordenado (un, 7 n),


con lo cual queda establecida la correspondencia entre las normales a los
planos que pasan por el punto "P" del material en estudio y el punto (un, 7,,)
de dicho plano (Fig. 1.11) .
1.6. SOLUCiÓN GRÁFICA DE MOHR 19

,~ ____.Q (crn;tn)
'tn r-

Figura 1.11: Plano de Mohr.

Sin embargo, la situación inversa no es verdadera, es decir, existen puntos


del plano de Mohr, cuyas coordenadas no representan a esfuerzos actuantes
en ninguno de los planos que pasan por el punto "P" . Por lo tanto, se plantea
la siguiente pregunta.
~ ¿Cuál es la región del plano de Mohr, cuyos puntos representan a esfuerzos
reabnente actuantes en el punto "P"?
Para responder a la pregunta planteada se procede como sigue:
Sea Ti la normal a un plano que pasa por un punto "P" de un continuo, cuyos
cosenos directores sean n x , n y , n z .
Asociado a dicho plano existirá un vector esfuerzo In, tal que:
(1.33)
Aceptemos además que en el punto existan las direcciones principales 1, 2 Y
3, en las cuales se definan a los planos principales, con esfuerzos (TI, (T2 Y (T3.
Tomando como base el marco de referencia anterior, la ecuación 1.33 puede
ser escrita como:

In = [81i 2
~3] I ~~ I (1.34)

El vector esfuerzo resulta:

(1.35)

El esfuerzo normal es igual a:

(Tn = -t n en
-
= (Tlnx2 + (T2ny2 + (T3nz2 (1.36)
Haciendo uso de estas dos últimas ecuaciones, conjuntamente con:
2
nx + n 2y + n z2 = 1
20 1. ESTADO DE ESFUERZO

resolvemos el problema plant.eado, es decir, determinamos al Vé'Ctor Ti que


define al plano en el cual actúan a n y Tu.
El sistema. de ecuaciones planteado es:

IInl 2 a 2lx 2n 2 + a2n~


n 2 + a 2y l~
(1.37)
222
an alnx + a2ny + aln;; (1.38)
1 n x2 +n2y +n ;:;2 ( 1.39)

Las ecuaciones 1.36, 1.37 Y 1.38 pueden representarse de manera. matricial


como:

(1.40)

Resolviendo el sistema de ecuaciones planteado para n;, n~ y n~,mediante el


método de Kramer, se tiene:
Cálculo del determinante de la matriz de coeficientes.

Para n!, se tiene:

(a! + r!) (a2 - (3) - a~ (an - (3)

+a~ (a n - (2)
(a! + r!) (a2 - (3) - a~an + a~a3
+a~an - ~a2
(a~ + r!) (a2 - (3) - a n (a2 + (3) (a2 - (3)
+a2a3 (U2 - (3)

(u; + T;) (a2 - U3) - a n (U2 + U3) (a2 - (3) + U2 a 3 (a2 - U3)
= a~ (a2 - U3) - uª (al - (3) + a5 (al - (2)
(a; + T~) - a n (a2 + (3) + a2a3
-a~ (al - (3) + u~ (Ul - (2)
u 21 + _-=--'---_---.:~_...::....:c......:_--=~
(a2 - (3)
1.6. SOLUCiÓN GRÁFICA DE MOHR 21

El denominador se puede simplificar como sigue:

2 + +
-0"~0"1 0"~0"3 O"~O"I - 0"~0"2
0"1 +--~----~----~--~~
(0"2 - 0"3)
==>
2 (0"2 - 0"3) (-0"20"1 + 0"20"3 - 0"10"3) ( ) ( )
O" 1 + ( )
0"2 - 0"3
= 0"2 - 0"1 0"3 - 0"1

n2 = (O";t + T;J - O"n (0"2 + 0"3) + 0"20"3 (1.41)


x (0"2 - O"¡) (0"3 - 0"1)

De la misma manera se obtiene para n~ y n~:

(O"~ + T~) - O"n (0"1 + 0"3) + 0"30"1 (1.42)


(0"3 - 0"2) (0"1 - 0"2)

(O"; + T;) - O"n (0"1 + 0"2) + 0"10"2 (1.43)


(0"1 - 0"3) (0"2 - 0"3)

Procediendo por simple análisis de estas ecuaciones, podemos definir a la


región buscada. Así en la eco 1.41, observamos lo siguiente:
Dado que

concluimos que:

y como

concluimos que

ecuación que puede ser escrita como

Completando el cuadrado del primer miembro


22 1. ESTADO DE ESFUERZO

(J2+(J3)2 2 ((J2-(J3)2
(
\(Jn- 2 +7n~ 2 (1.44)

Obsérvese que con el signo de igualdad la ecuación anterior define a un círculo


de radio:

Rn" =
(J2 - (J3
2 Y centro en en" ((J2 +2(J3' O)
De esta manera concluimos que en lo referente a "nx", la región buscada
viene dada por los puntos del círculo definido y los exteriores a él.
Respecto a n~ obtenemos:

(J _ (JI + (J3)2 72 > ((JI - (J3)2


(1.45)
( n 2 + n - 2

Es decir, la región buscada queda definida por los puntos del círculo de radio:

R,'y =
(JI - (J3
2 Y centro en eny ((JI +2(J3) O) . .
, así como los puntos mtenores
a él.
Respecto a n~ obtenemos:

(JI +2 (J2)2 +7,,~


2
(1.46)
( (Jn-

Es decir, se trata de los puntos del círculo de radio:

R,'z =
(JI - (J2
2 Y centro en: en. ((JI +2(J2) O) .
,asÍ como los puntos ex tenores
a él.
Por lo tanto, la región buscada es la parte achurada que se muestra en la Fig.
1.12.
1.6. SOLUCiÓN GRÁFICA DE MOHR 23

't I

Rnz = 0"1 +
2
03

<J2+<J:3
2

Figura 1.12: Región de Mohr.

Ejemplo 1.4. Determine los esfuerzos (5n Y T n asociados al plano cuyos cosenos
directores están dados por n (~ , ~, ~). El tensor esfuerzo está dado por:

T ij = [180 g J 3
] [MPa]

Cálculo del vector esfuerzo tn.

_t = [10O 6O O]
n O
O O -3

Cálculo del esfuerzo normal (5 n'


U] =
- - 3-
5i + 3j - .J2k [MPa]
24 1. ESTADO DE ESFUERZO

Cálculo del esfuerzo cortante T".

Tu ±5.68 [M Pa]

1.7 Procedimiento gráfico de Mohr

Únicamente se enunciarán los pasos a seguir para obtener los esfuerzos a n y


Tn asociados a un plano cuya normal es n, estando fuera del alcance de estas
notas la demostración rigurosa del procedimiento.

Para obtener los esfuerzos a n Y T n mediante el procedimiento gráfico de Mohr,


se siguen los siguientes pasos:

1) Una vez definida la región de Mohr, por al se traza una recta perpendicular
al eje a y a partir de ella, pasando por al trazamos una recta que con la
primera forme un ángulo el (el = L cos n x ).

2) Por a3 se traza una recta perpendicular al eje a y a partir de ella,


pasando por a3 se traza una recta que con la primera forme un ángulo 3 e
(e 3 = L cos n::).
e
3) Haciendo centro en nx , se traza el arco de círculo que pasa por los puntos
e e
donde la recta A - A' intersect.a a los círculos cuyo centro es ny y nz •

4) Con centro en en::, se traza el arco de círculo que pasa por los puntos
donde la recta B - B' intersecta a los círculos cuyo centro es e,l,X y eny .

5) El punto buscado viene dado por la intersección de los dos arcos de


círculo anteriormente trazados; las coordenadas del punto Q representan a
los esfuerzos (an,T n ) actuantes en el plano cuya normal es n.

Solución del ejercicio 1.4, mediante el procedimiento gráfico de Mohr.


1.8. CASOS PARTICULARES DE ESTADOS DE ESFUERZO 25

T
I
A' I
",, ~

I '~
I 'L~\'1Xi'''--~

I I
Figura 1.13: Solución gráfica de Mohr.

Los ángulos (h y 03 se calculan como:

01 Lcosn x
03 Lcosnz

Siguiendo la construcción geométrica indicada en la Fig. 1.13 se obtiene:

(jn = 2.5 fMPaJ} Punto "Q"


{ T" = 5.7 MPa

1.8 Casos particulares de estados de esfuerzo


a) Estado de esfuerzo nulo; es aquel en el cual (jI = (T2 = (j3 = O.

T.,,= O O O]
'J [ OO OO OO
En el espacio de Lamé y de Mohr está representado por un punto, ubicado
en el origen del sistema de referencia.
26 1. ESTADO DE ESFUERZO

Fig'ura 1.14: Repre.:;entación gráfica de un estado de esfuerzo nulo.

b) Estado de esfuerzo uniaxialj es aquel en el cual O'} =1= Oj 0'2 = 0'3 = O.

~j =
O'}
O O O
O O]
[
O O O

Lamé

A-----~2

1
Figura 1.15: Estado de esfuerzo uniaxial, espacio de Larné.

Todos los vectores In, tienen la dirección del eje 1.

Mohr
"
Tn \
~------~ a
-4 al 1---
Figura 1.16: Estado de esfuerzo uniaxial, plano de Mohr.
- -- - - - -- - - -- - - -- -- - - -- - - --------------------

1.8. CASOS PARTICULARES DE ESTADOS ... 27

c) Estado de esfuerzo plano; es aquel en el cual uno de los esfuerzos principales


es nulo, por ejemplo, 0"1 i- O; 0"2 i- O; 0"3 = O.

[8 ~2 8]
1
Jij =

3
Lamé

1
Figura 1.17: Estado de esfuerzo plano, espacio de Lamé.

Mohr

v~"
~ \
IL......,_ _-----.l'--_ _ _ _---'-_~ cr
cr 2 cr 1
Figura 1.18: Estado de esfuerzo plano, región de Mohr, 0"3 = O.

Cuando dos círculos pasan por el origen se trata de un esfuerzo plano (Fig.
1.18).

d) Estado general de esfuerzos; es aquel en el cual 0"1 i- O; 0"2 i- O; 0"3 i- O.

El estado general de esfuerzo en el espacio de Lamé queda definido por un


elipsoide.

La representación en el plano de .Mohr quedará como en el caso anterior, con


la variante de que 0"3 i- O.
28 1. ESTADO DE ESFUERZO

e) Es(.ado ele esfuerzos hidrostático es aquel en el cual 0"1 = 0"2 = 0"3 = 0".

0"00]
Tij =
[
O O" O
O O O"

El e:;t.ado de esfuerzos hielrostático en el espacio de Lamé queda representado


por Hna esfera, mientras que en el plano de Mohr es un segmento de recta
sobre el eje 0".

Mohr

,----------~>~----~>7a
kIEé---- a --::>~I
Figura 1.19: Estado de esfuerzo hidrostático.

1.9 Descomposición de un estado general


de esfuerzos en sus componentes
volumétrica y desviadora
Supongamos que en un punto "P" de un medio continuo se establece el
siguiente tensor esfuerzo.

O"xx 7 yx 7 ZX ]
Tij = 7 xy O"yy 7 zy (1.47)
[ 7 xz 7 yz O"zz

Sea n la normal a un plano que pasa por un punto "P" de un medio continuo,
entonces:

T" [Tij] Inl


7 71 (O" xxnx + 7 yx ny + 7 zx n z ) I + (7 xynx + O"yyny + 7 zynz )] (1.48)
+ (7 xz nx + 7 yz ny + O" zz n z ) k

La ecuación 1.48 puede ser escrita en forma matricial como:


1.9. DESCOMPOSICION DE UN ESTADO GRAL... 29

tn [ u xx
T xy
T xz
-
~
3
U yy
T yx
_.!..L
T yz
3
u zz
T zx

r., .!..L 1
- 3
r nx ]
l~ +
[~O
O O
3
.!..L
O
O O .!..L
3 1[E~ ]
In = [Tollnl + [Tvllnl ( 1.49)

donde:
To = componente desviadora del tensor 7'¡ '.
Tv = componente volumétrica del tensor j , n
Obsérvese en la ec, 1.49 que la componente volumétrica representa un estado
de esfuerzos hidrostático, cuyos esfuerzos principales son iguales a --f, es
decir, el valor promedio de los esfuerzos normales que aparecen en 7'¡j.

La componente distorsional de un estado general de esfuerzos siempre


presenta la siguiente propiedad:

[Tol [U xx -
T xy
T xz
~
3
U yy
T yx
_.!..L
T yz
3
r..
T zy
U zz -
.!..L
3
1
T~x
r ]+ [ ~ T~y 1
O O
[Tol [O
T xy O]
O + [OO O O
T
O O O T xz O T yz O J

(1.50)

Obsérvese que:

_(u _1xx
1
3 ) - (u _1
zz 3)
1
-u xx - u zz
21
+31

Lo anterior indica que la componente distorsional o desviadora de un estado


general de esfuerzos resulta siempre ser igual a la suma de 5 estados de
esfuerzos plano,
30 1. ESTADO DE ESFUERZO

Algunos autores llaman a estos estados de esfuerzo plano "estados de cortante


puro" en los cuales el máximo esfuerzo cortante es igual al valor de
cualesquiera de los esfuezos principales diferentes de cero.
Es común considerar que en la componente desviadora se manifiesta la acción
del estado de cortante puro, mientras que en la componente volumétrica
se manifiesta la acción de los esfuerzos normales promedio en el punto en
estudio.
En general, en todos los estados de esfuerzo se acostumbra asociar los cambios
de volumen a la componente volumétrica, mientras que los cambios de forma
se asocian a la componente distorsional o desviadora.

1.10 Estado de esfuerzo plano


En la Fig. 1.20 se muestra un cuerpo deformable sometido a ciertas
solicitaciones en su plano. Si el elemento es de espesor pequeño, los esfuerzos
asociados a la dirección z se pueden despreciar, esto es:

O"zz = O T xz = O T yz =O

?f X
/'
/'
/'
/'
/'
/'
r----j--X
z

y
Figura 1.20: Esta.d9 de esfuerzo plano, en un sistema de referencia
cartesIano.

Por lo que el tensor esfuerzo T;j se reduce a:

(1.51)

Este tensor representa un estado de esfuerzo plano.


Las componentes normal y cortante del vector esfuerzo In, para un estado de
esfuerzo plano se pueden deducir como sigue:
1.10. ESTADO DE ESFUERZO PLANO 31

y n

/
/'
/
/
t n /

Figura 1.21: Derivación de los esfuerzos (Tn Y Tn para un estado de


esfuerzos plano.

Tomando en cuenta la Fig. 1.21 los vectores unitarios n y m, se pueden


expresar como:

n = (cos e) I + (sene)) (1.52)


m (sene)I - (cose) J (1.53)

En estas dos direcciones se pueden establecer las componentes normal y


tangencial del vector esfuerzo In, como:

[T.'..3 ] Inl = [(TXX


Txy ij;,x
YY
11
J
cos e
sene 1
(1.54)
(Tn In · n
((Txx cos 2 e + Tyxsene cos e) + (T xysenO cos e + (Tyysen 2 e) (1.55)

Si (T n es positivo, quiere decir que tiene la dirección de n.


Para el esfuerzo cortante, se tiene de manera similar:

Si Tn es positivo, esto quiere decir que tiene la dirección de m.


Haciendo simplificaciones, se tiene que :
32 1. ESTADO DE ESFUERZO

O" xx COS 2 O+ O"1/Y sen 2 O + 27 xy .senO cus O (1.57)


(O" xx - O"yy) senO cos 0+ 7 X1) (sen 2 0 - cos 2 O) (1.58)

Tomando en cuenta las siguientes identidades trigonométricas:

sen2fJ 2senOcos O (1.59)


1 - cos20
scn 2 0 (1.60)
2
1 + cos20
cos 2 O (1.61)
2

y sustituyendo estas últimas expresiones en las ecs 1.57 y 1.58, se llega a:

O"n O" xx (1 + ~os 20) + O"yy (1 - ~os 20) + 7 xysen20

7n (O"xx;O"YY)sen20+7 XY [(1-~OS20) _ (1+~OS20)]

Simplificando se tiene:

O"xx + O"yy 1
2 + "2 (O" xx - O" yy) cos 20 + 7 xy sen20 (1.62)
O"xx - O"yy) (1.63)
( 2 sen20 - 7 xy cos 20

Derivando O"n respecto a O, se tiene:

= - (O" x - O"y) sen20 + 27xy COS 20


-27" ( 1.64)

Este resultado muestra que 7 n es una medida del cambio de O"n con respecto, a
O, Y que cuando 7 n = O, se tendrá el máximo o mínimo valor de O"n (esfuerzos
principales) .
1.11. PLANOS Y EJES PRINCIPALES DE ESFUERZO 33

Si se acepta que los ejes de referencia coincidan con las direcciones en que
= 0, el tensor que representa al estado de esfuerzos en un punto del medio
T ..
continuo, para 0'2 = 0, quedará definido por:

T,"=[0'01 0
t] 0'3 ]
Para este tensor, las ecs 1.62 y 1.63 quedan como:

0'1 + 0'3 1
--2- + '2 (0'1 - 0'3) cos 2tJ (1.65)

T" == ( 0'1 ; 0'3 ) sen2tJ (1.66)

Estas últimas ecuaciones se pueden representar como sigue:

_0'1+0'3) 1
(
O'
n 2 '2 (0'1 - 0'3) cos 2tJ (1.67)

1
- (0'1 - 0'3) sen2tJ (1.68)
2
Elevando al cuadrado las ees 1.67 y 1.68 Y sumándolas, obtenemos:

0'1 +0'3)2 2 (0'1 - 0'3)2


( O'n- 2 +Tn = 2 (1.69)

La ecuación 1.69 nos muestra que el conjunto de vectores t n asociados a


planos cuyos vectores esfuerzo esten en el plano de esfuerzos deben tener
sus extremos sobre la circunferencia con centro en: (0'1; 0'3, O) y radio:

(0'1 ; 0'3 ).
Para definir a un vector esfuerzo dado, asociado al plano cuya normal forma
un ángulo tJ con el eje 1, es necesario girar un ángulo 2tJ, en dirección contraria
a () para definir al extremo del vector tn , cuyas proyecciones son O' n y T n'
Esta representación permite así definir a todos los vectores tn en planos
normales al plano de esfuerzos.

1.11 Planos y ejes principales de esfuerzo


Por definición los planos principales de esfuerzo son aquellos donde el esfuerzo
cortante T" = O.
34 1. ESTADO DE ESFUERZO

Por lo que tomando en cuenta la eco 1.63, obtenemos:

( (Jxx ; (Jyy) sen2(J _ T xy cos 2(J = O

Desarrollando se llega a :

(1. 70)

Las raíces de esta última ecuación representan las direcciones de los ejes
principales nI y n3.
Para encontar las raíces de la eco 1.70, hacemos uso de la Fig. 1.22.

2
± + 't xy

28
O"xx-ayy
2
Figura 1.22: Relación geométrica entre (J xx, (Jyy, T xy y (J.

(Jxx - (Jyy

cos 2(J 2 (1.71)


± ( (J xx ; (J yy ) 2 + T;y
T xy
sen2(J (1.72)
((JXX-(Jyy)2 2
± 2 +TXY

Sustituyendo las ecs. 1.71 y 1.72 en la 1.62, obtenemos a los esfuerzos


principales (JI y (J3:

2
(JXX - (Jyy) 2
( 2 + T xy (1. 73)
1.11. PLANOS Y EJES PRINCIPALES ... 35

Para el estado de esfuerzo plano, la ecuación característica se reduce a:

(1.74)
Donde:

O"xx + O"yy
det [ O"xx
T xy T yx ]
O"yy

Ejemplo 1.5. El estado de esfuerzo plano en un punto de un medio continuo


está dado por:

T ·=
IJ
[!1 O
-lO 8] [M
O
Pa]

. Usando los métodos analítico y gráfico del Mohr, determine:


a) Los esfuerzos normal y cortante asociados al plano de corte indicado en la
Fig. 1.23.
b) Los esfuerzos principales.
c) Las direcciones principales de esfuerzo, referidas a los ejes x, y en que se
define el tensor esfuerzo.
Los elementos del tensor esfuerzo son:
O"xx = 4 [MPa] ; O"yy = 2 [MPa] Txy = -1 [MPa]

Cara @
n

4 Cara@

Figura 1.23: Representación del estado de esfuerzo para el ejemplo 2.5.

El ángulo () que permite calcular los esfuerzos O"n Y T n, mediante las ecs. 1.62
y 1.63, vale: () = 135 0 •
Esfuerzo normal:

4+2 4-2
O"n -2- + -2- cos 270 + (-1) sen270°
0

O"n 3 + 1 = 4 [M Pa]
36 1. ESTADO DE ESFUERZO

Esfuerzo cortante:

4-2
-2-.sen270 o
- ()
-1 cos 270o
T'I

-1 [MPa]

LO!:i esfuerzos principales valen:

0"1,3
4+2
-2- y
± 1(4-2)2
-2- + (_1)2. = a ± Ka
0"1 3 + Ka = 4.41 [M Pa]
0"3 3- ffo = 1.60 [M Pa]

Cálculo de las direcciones principales:


-1
tan20~ (4;2) ~-1
Las raíces de la ecuación son:
el = -22.5°
e3 = -112.5°
Para comprobar la dirección asociada a los esfuerzos O"l y 0"3, se sustituyen
las raíces obtenidas en la eco 1.62.

4 +-2 + -
- 4 - -2 cos ( -45 0) + (-1 ) sen (-45 0)
2 2
3 + 0.7071 + 0.7071 = 4.41 [M Pa]

Por lo tanto el ángulo el = -22.5° esta asociado a la dirección del esfuerzo


principal mayor 0"], obviamente el ángulo e3 = -112.5° debe estar asociado
al esfuerzo principal menor.
Para la solución gráfica de Mohr haremos uso del procedimiento conocido
como método del "polo".
El polo es un punto en el círculo de Mohr el cual se puede ubicar como sigue
(Fig. 1.24).
a) Se traza primeramente el círculo de Mohr a partir de los puntos conocidos
A y B. Para que el signo del esfuerzo cortante Tu sea el correcto, es necesario
cambiar el signo del esfuerzo cortante T xy asociado a la cara A antes de trazar
el círculo de Mohr.
1.12. ESFUERZOS OCTAÉDRICOS 37

b) Por el punto A se traza un plano que sea paralelo a los esfuerzos asociados
a la cara A. Donde este plano corte el círculo, estará ubicado el polo.

c) Se hace lo mismo con el punto B para verificar la posición del polo.

1
/

---+--~~--~--~~/--~a

-1
-2

Figura 1.24: Solución gráfica de Mohr, ejemplo 2.5.

d) Una vez conocida la ubicación del polo, y a partir de éste, bastará con
trazar una recta que tenga la misma orientación del plano donde se desea
calcular los esfuerzos Un Y T TI '

e) El punto donde esta recta corte al círculo, determina los esfuerzos Un Y T n


correspondientes. En la Fig. 1.24 se muestra la aplicación del procedimiento
seguido para el problema planteado. Para este caso particular el punto
correspondiente al polo coincide con los esfuerzos Un Y T n buscados.

1.12 Esfuerzos octaédricos

Considérese un punto en el cual exista un estado general de esfuerzos. En ese


punto existen ocho planos octaédricos definidos por normales que dividen en
ángulos iguales a cada uno de los octantes del espacio.
38 1. ESTADO DE ESFUERZO

3
a=j3=y

°3 00ct

(a
~a/
a ~
A<-'' ':''--'---é;;'---;;¡'' 2
P
" V y.,plano '\ °2
1 ,é··
./
/c: 1
octaédrico

Figura 1.25: Definición de esfuerzos octaédricos.

Consideremos el plano octaédrico del primer octante, cuya normal es:

1- 1- 1-
noct = -i+ - j +-k (1.75)
v'3 v'3 v'3
Asociado a este plano existirá un vector esfuerzo total octaédrico, definido
por:

toct millnoctl (1.76)

toct [8 (f J,l [~ 1
1 o

1 - 1 - 1 -
toct = -O"Ii + -0"2j + -0"3 k (1.77)
v'3 v'3 v'3
El esfuerzo normal octaédrico es igual a:

- 1[- - I:l Ir:; - I:l


= t oct · noct = v'3 O"Ii + 0"2 j + 0"3 k J • v'3 Li + j + kJ
1
:3 [0"1 + 0"2 + 0"3]
JI
(1.78)
3

El esfuerzo cortante octaédrico puede calcularse como:


1.13. EQUILIBRIO DE PARTíCULAS EN UN MEDIO CONTINUO 39

De los resultados obtenidos se puede observar lo siguiente:


a) El esfuerzo normal octaédrico depende de la magnitud de la componente
volumétrica del tensor de esfuerzos.
b) El esfuerzo cortante octaédrico tiene como magnitud el valor medio
cuadrático de los elementos de la diagonal principal de la componente
distorsional.

1.13 Equilibrio de partículas en un medio


continuo
Consideremos que en cada punto del continuo es conocido el tensor T;.j, dado
por:

O'xx zx
T;.j = T xy T zy ]
[ T xz O'zz

Aceptemos que las componentes de T;.j son funciones continuas de (x, y, z, t).
Cuando el tensor es simétrico será necesario definir a seis funciones continuas
y derivables de variación suave, para conocer el estado de esfuerzo en todos
los puntos del continuo.
Para ello basta estudiar el equilibrio de partículas en el continuo, de ahí
surgen las condiciones de frontera y las fuerzas de cuerpo que producen a T;.j.
Consideremos en primer lugar a una partícula interior del continuo.
40 1. ESTADO DE ESFUERZO

F1
\
\

(~d'r--F,
\ 'v/--.... .. "

)"J-:d':~/
F3
Figura 1.26: Equilibrio estático de una partícula interior del continuo.

En la partícula elemental aparecen:

a) Fuerzas en las caras laterales, provocadas por la acción de partículas


vecinas.

b) Fuerza" de cuerpo generadas por la acción de cuerpos distantes.

dx ...J-------:7 y

dz
dV = dx dy dz

dy
Figura 1.27: Fuerzas de cuerpo y de superficie en una partícula elemental
de un medio continuo.
1.13. EQUILIBRIO DE PARTICULAS EN UN MEO... 41

En el elemento diferencial, que se muestra en la Fig. 1.27 no están incluidas


las fuerzas de cuerpo que pudieran existir.
Por equilibrio de fuerzas, se tiene:
Por L F.'t = O:

(u x + du x ) dydz - (uxdydz) + (T yx + dT yx ) dxdz - (Tyxdxdz)


+ (T zx + dT zx ) dxdy - (Tzxdxdy) + xdxdydz = O
duxdydz + dTyxdxdz + dTzxdxdy + xdxdydz = O
uUx éh yx OT zx
ux dxdydz + oy dxdydz + --¡¡;-dxdydz + xdxdydz = O

oUx dV OT yX dV OTzx dV f dV = O (1.81)


ox + oy + OZ + x
Si dV -+ O

ou x aT yX OTzx f - O (1.82)
ox + ay + az + x-
De manera similar; por L Fy = O, se obtiene:
OTxy alFil aTZY 1" _ O
ox + ay + az + JII- (1.83)

Por L Fz = O, se obtiene:

aT xz m yZ auz f - O (1.84)
ox + ay + az + z-

Las ecs. 1.82, 1.83 Y 1.84 pueden ser escritos en notación indicial como:

f1Iij + Ji = O
aXj

Esta última ecuación se conoce como ecuación de Cauchy.


Ejercicio 1.6. En un medio continuo se establece un tensor de esfuerzos
definido por:

T¡j
2xy
= A . [ c2 Oy2 g g1 donde:
él - y2 O
A, e = ctes.

Obtenga:
a) Las fuerzas de cuerpo.
b) Las fuerzas en las fronteras del continuo.
42 1. ESTADO DE ESFUERZO

Figura 1.28: Viga de sección prismática sometida a flexión simple.

Analizando cada una de las caras del continuo, se tiene:


En el plano y = Cj n = }
El vector esfuerzo In vale:

_
tn = [7ijllnl = A· 8 88] [b]
[ 2cx

In = O ~ Toda la cara está descargada.


En el plano y =- Cj n = -}

In = O ~ Cara descargada.
En el plano z = +~j n=k

In = O ~ Cara descargada.
En el plano z = -~j n =-k
In = O ~ Cara descargada.
En el plano x = O; n = -z
1.13. EQUILIBRIO DE PARTICULAS EN UN MEO... 43

1" nj (n) ~ A [ c' ~ y' e' t 8j[ Y1


in A (y2 - c2) :J
Un in . Ti = O = [A (y2 - c2 ):J] . [-7] = O

Tn = ±VJinJ2 -IO'nI 2
Tn = A (y2 - c2) = Ay2 - Ac2

Si y = O; Tn = - Ac2
En la Fig. 1.29 se muestra la distribución de esfuerzos cortantes
. correspondiente.

-----F----
I

r-1~
Figura 1.29: Distribución de esfuerzos cortantes en el plano x = O.

Integrando el volumen de esfuerzos cortantes para un ancho unitario, se


obtiene la fuerza cortante que actua en la cara x=O:

2
v :3 (_c 2 A) (2c)(1)
v 4
--e A
3
3

De donde la constante A vale:

3V
A=---
42
En el plano x = L; Ti = 7
44 1. ESTADO DE ESFUERZO

El esfuerzo normal y el esfuerzo cortante tienen por valor:

0"" = I". n = 2ALy


Tu = A (c 2 _ y2)

La distribuci6n de esfuerzos normales en x = L se muestra en la Fig. 1.:30.

2LcA

H~
T
e
2~e
+
e
L -H
)M

-2 LeA + 1 +
Figura 1.30: Distribución de esfuerzos normales en el plano x = L.

Dicha distribución genera un par interno que vale:

Siendo

H = 2LcAc = c2 LA
2
El brazo de palanca del par vale:
4
d= -e
3
Por lo tant.o el momento interno resulta igual a:

La distribuci6n de esfuerzos cortantes produce una fuerza cortante igual a:

Esta fuerza está en equilibrio con la que se produce en x = O.


1.14. ECUACIONES DE EQUILIBRIO DE MOMENTOS EN PARTíCULAS DE UN
CONTINUO 45

El momento de inercia centroidal de la sección se puede expresar como:

1 = ~C3 MiL
3 ,por lo que A = 21

Por equilibrio el momento interno que produce la distribución de esfuerzos


normales en x = L, debe ser igual al momento que produce la distribución
de esfuerzos cortantes en ambos extremos de la barra.
Dado que el esfuerzo normal que actua en la cara x = L, se puede expresar
como:
O"tt = 2ALy , Y sustituyendo el valor de A en esta última expresión, se llega a:

Ecuación conocida en los textos de mecánica de materiales como fórmula de


la escuadría.

1.14 Ecuaciones de equilibrio de momentos


en partículas de un continuo

dx
tiI )--,
x

dy (tzydxdy)
Figura 1.31: Ecuaciones de equilibrio de una partícula de un medio
continuo.

En la Fig. 1.31 se muestra el estado de esfuerzo a que se encuentra sometido


una partícula elemental de un medio continuo. Por facilidad únicamente se
indican las fuerzas en las caras que van a producir giros alrededor del eje x.
Por ¿Mx =O
46 1. ESTADO DE ESFUERZO

dy dy dz
(Tyz + dTyz) dxdz-;¡ + Tyzdxdz-;¡ - (T:;y + dTzy ) dxdy-;¡
dz
-Tzy dxdy- + 1nxdv = O
2
dT yz - dT zy + 2mx = O

aTYZd aT _ 2 - O
-- y - - ZYd
- "'+ m x - (1.85)
ay az

Dp. manera análoga; por L My = O, se obtiene:

aT-
- zxd z- - xzd x+ 2m = O
aT- (1.86)
az ax y

Finalmente, por L Fz = O, se llega a:


aT-
- xy dx - -
aT-
yx d y + 2m z = O (1.87)
ax ay

De las ecs. 1.85, 1.86 Y 1.87 se pueden definir los valores de m x , my y m z ,


debidos a la excentricidad de las acciones que se generan en la partícula
respecto al centroide.
Para el ejemplo 1.6, la aplicación de las ecuaciones de equilibrio de momentos
conduce a:
Por L Mx = O; O - O+ 2mx = O ==> mx = O

L My = O; O- O+ 2my = O ==> my = O
L Mz = O; O - (- 2y) Ady + 2m = O z

m z = -Aydy

Donde:

Este resultado muestra que existirán tendencias al giro de las partículas


alrededor del eje z, dependiendo de su posición en el continuo.

La existencia de estos momentos implica que las partículas no se mueven


como partículas rígidas, sino que deben presentarse distorsiones para
garantizar la existencia de continuidad.
1.14. ECUACIONES DE EQUILIBRIO DE MOMENTOS ... 47

Los pares de esfuerzo se autoequilibran con los correspondientes a las


partículas simétricas respecto al eje z, sin que exista la necesidad de aplicar
momentos externos para equilibrarlas. Estos momentos se producen en el
problema en estudio por la distribución no uniforme de esfuerzos en las caras
de las partículas.

Normalmente en campos de esfuerzos no uniformes existen estos pares de


esfuerzos que tratan de provocar rotaciones y distorsiones en las partículas.

Así, mediante las ecuaciones de equilibrio se pueden identificar la.')


7
características del vector fuerza de cuerpo y del vector momento m, dados
por:

7 xI +y] + zk
m mx2+my"]+mzk
48 1. ESTADO DE ESFUERZO
Capítulo 2
Estado de Deformación

Consideremos una partícula elemental de un medio continuo descargado,


y llamemos a este estado del cuerpo la configuración inicial (t=O) o no
deformada. Cada una de las partículas que forman el medio continuo en
esta configuración reciben el nombre de partículas materiales, cuya posición
se establecerá en el sistema de referencia cartesiano XI,X2,X3.

Al someter al cuerpo deformable a un sistema de cargas se generan estados de


esfuerzo los que a su vez provocarán estados de deformación. La forma que
adopta el medio para un tiempo t cualquiera se conoce como configuración
deformada y corresponde a la imagen que proporciona una fotografía del
cuerpo para el tiempo especificado. Cada una de las partículas que forman
el cuerpo en esta configuración reciben el nombre de partículas espaciales,
cuya posición se establecerá en el sistema cartesiano 311, X~, X;.

t=Q t=t
Configuración inicial x;

Configuración deformada

/
Xl
Figura 2.1: Definición de vector desplazamiento.

Para describir la cinemática del continuo se supondrá por comodidad que los
sistemas de referencia X¡,X2,X3 Y x;,x~,x;, están superpuestos (Fig. 2.1).

49
50 2. ESTADO DE DEFORMACiÓN

La posición de cada partícula del medio en la configuración inicial se puede


definir mediante un vector de posición O P cuyas componentes son:

(2.1)

Al transcurrir el tiempo de t = O a t = t , el vector de posición del punto que


define a la partícula quedará definido por el vector de posición:

(OP') = x;?: + x;} + x;k (2.2)


De la Fig. 2.1 se puede establecer que:

(2.3)
Siendo 8p el vector desplazamiento del punto P, el cual puede ser escrito
como:

8p = (x;~ + x;) + x;k) - (x]~ + X2J + X3k)


8p = (x; - x) ~ + (x; - y)} + (x; - z) k
Llamando a las componentes del vector desplazamiento:

x]I - X U (2.4)
x; -y v (2.5)
x~ - z w (2.6)
Tomando en cuenta estas últimas ecuaciones, el vector desplazamiento se
puede escribir como:

(2.7)
Evidentemente las componentes del vector desplazamiento u, v y w, dependen
de la posición del punto en el continuo y del tiempo transcurrido de t = O a
t = t.

Si el vector desplazamiento está asociado con cualesquier partícula en


la configuración no deformada, entonces sus componentes se pueden escribir
(notación indicial) como:
(2.8)
En el caso de que el vector desplazamiento esté asociado con cualquier
partícula en la configuración deformada, se tiene:
(2.9)
51

Así, las ecs. 2.8 y 2.9, continuas y derivables, permitirán describir el campo
vectorial de desplazamientos de todos los puntos de un continuo en cualquier
instante.
El lugar geométrico de los extremos de vectores desplazamiento 15 respecto
a la posición original describirá la trayectoria de cada punto material del
continuo.

El conocimiento de las trayectorias contribuye a describir la cinemática de las


partículas de un medio continuo, pero además se debe describir el movimiento
relativo de la partícula respecto a la trayectoria para conocer la cinemática
completa de las partículas.

Para ello se analizará una partícula elemental del continuo, tal como se
muestra en la Fig. 2.2, donde por facilidad de exposición se adoptará un
sistema de referencia cartesiano x, y, z.

dy _

~a'
z ~ Va .'

dz \
P P'
t =o t =t
y

/x
Figura 2.2: Definición de vector desplazamiento relativo.

Supongamos una partícula material del medio (t = O) y fijemos dos puntos


de dicha partícula, por ejemplo P y Q, separados una distancia diferencial.
La posición de dichos puntos está dada por las coordenadas:
P(x , y , z) ; Q(x+dx, y+dy , z +dz)

Al transcurrir el tiempo, la posición de la partícula en la configuración


deformada estará dada por las coordenadas de los puntos pi y Q', esto es:

pi (x + u, y + v , z + w)
Q'(x+dx+u+du, y+dy+v+dv , z +dz +w+dw)
Así, la nueva posición de la partícula quedará definida una vez que se
establezca el campo de desplazamientos 15, mediante funciones continuas y
derivables.
52 2. ESTADO DE DEFORMACiÓN

2.1 Características de rotación y cambio de


forma de la partícula

Con referencia a la Fig. 2.2, llamaremos al vector ó.J la diferencia vectorial


siguiente:

(2.10)

El vector ~J mide la rot.ación y el cambio de forma de la partícula y representa


el movimiento relativo del punto Q respecto al punto P.

Evaluando la ec.2.1O, el movimiento relativo de la partícula quedará descrito


por:

SJ = (u + du) z+ (v + dv) J + (w + dw) k - (uz + vJ + wk)


ó.J = duz + dVJ + dwk (2.11)

Desarrollando la eco 2.11, se tiene:

~J (8u dx + 8u dy + 8u dZ) z+ (8V dx + 8v dy + 8v dZ) J


8x 8y 8z 8x 8y 8z

+ (~:dX+ ~;dY+ ~:dZ) k (2.12)

La eco 2.12 se puede escribir de manera matricial como:

8u 8u 8u
- -
8x 8y 8z

~~ I
8v 8v 8v
ó.J = (2.13)
8x 8y
8w 8w
8z
8w
I
-
8x 8y 8z
2.1. CARACTERISTICAS DE ROTACION y CAMBIO... 53

De la Fig. 2.2 se puede observar que el vector que une los puntos P y Q
resulta igual a:

PQ = dxí + dyy + dzk (2.14)

Dividiendo la eco 2.13 entre el módulo de PQ, se obtiene:

8u 8u 8u dx
- - -
8x 8y 8z !PQ!
bo8 8v 8v 8v dy
= (2.15)
!PQ! 8x 8y 8z
8w 8w 8w
!PQI
. dz
-
8x 8y 8z IPQII
Obsérvese en esta última ecuación que los elementos del vector columna
representan los cosenos directores del vector que va de P a Q, el cual se
puede escribir como:

e = cos nI + cos {H + cos ,k (2.16)

Llamando a !~~I el vector deformación total, E, la ecuación 2.15 se puede


escribir como:

E = [Eij] lel (2.17)

Siendo

8u 8u 8u
-
8x 8y 8z
8v 8v 8v
Eij = (2.18)
8x 8y 8z
8w 8w 8w
-
8x 8y 8z

Eij recibe el nombre de tensor deformación total.


Este tensor se puede expresar como la suma de dos matrices, esto es:
54 2. ESTADO DE DEFORMACiÓN

on
OX ~(&+aw~
- ~ (OU + OV)
ox 2 oy 2 oz ox
Eij
r"
~12 owox'~ a,,~
-2 - +
oy
ou
-
ox OZ 2 oz ay
ov
-
ay
~ (OV+ ow)
1 ov
-2 -oz+ay
ow
OZ
ow
-

o ~ (OU _ ov) ~(au_aw~


2 ay ox. 2 OZ OX
1 OV OW
+ ~2 (OV _ Ou)
ox ay o -
2 oz
---
ay (2.19)

~(Ow_ou)
2 OX OZ ~(~;-~~) o

La primera matriz se conoce corno tensor deformación , E ij , mientras que la


segunda recibe el nombre de tensor rotacional, nij .
Así el tensor deformación total se puede expresar corno la suma de dos
tensores, esto es:

(2.20)

2.2 Interpretación física del tensor rotacional


Definiendo a los términos del tensor rotacional, nij, corno:

nyx ~ (ou _ Ov). n. = ~ (ou _ Ow)


2 oy OX ' ~x 2 OZ ox
n ~ (OV_ Ou) . n. = ~ (OV _ Ow)
XY
2 ox ay , 2 oz
~y ay
nxz ~ (OW _ ou). n z = ~ (OW _ Ov)
2 OX OZ ' y 2 oy OZ
Obsérvese que los elementos del tensor rotacional guardan la siguiente
relación:

Por lo que dicho tensor se puede representar como:


2.2. INTERPRETACION FISICA DEL TENSOR. .. 55

(2.21)

Para identificar los términos del tensor rotacional nij , recordemos la definición
de rotacional de un campo vectorial.

k
o -_ (OW
--- Ov) -;- (OU -Ow) -;- (OV
- - Ou)-
TOt"3 =
oz oy OZ z+ - - J+
oz ox ox -oy k
W

TOt"3 = 2 [nyz2 + n zx J+ nXykJ (2.22)

Calculando el producto vectorial !rot"3 x e, se tiene:

1 -
2roto x e -n dxf3 n:
cos a
ZY
cos
y
cos I
{n zx cos , - n xy cos f3} 1, - {-n zy cos , - nxy cos a} J
+ {-nzy cos f3 - n zx cos a} k
1 -
2roto x e {n,:x cos ,+ nyX cos f3} 1, + {n ZY COS, + nxy cos a} J
+ {n yZcos f3 + nxz cos a} k (2.23)

Por otra parte, el tensor rotacional aplicado a la dirección definida por el


vector e, resulta igual a:

{(flux) cos f3 + (n zx ) cos ,} I + {( n XY ) cos a + (n ZY ) cos ,} J


+ {(n xz ) cosa + (n yz ) cosf3} k (2.24)

Obsérvese que la ecs. 2.23 y 2.24 resultan ser iguales, por lo que:

(2.25)
56 2. ESTADO DE DEFORMACiÓN

1 -
La interpretación física de 2rotJ X e se muestra en la Fig. 2.3, donde se puede
observar que dicho término representa al conjunto de vectores ti generados
I
por una rotación alrededor del eje rot8, con un giro w de magnitud ~rot81·
El conjunto de vectores ti tienen magnitudes proporcionales a su distancia d
al vector rot8.

y
a'
-----E~
dv \
'a +
du/' I

aR dy
8
p+ dx + "" x
Fig·ura 2.4: Rotación pura de una barra rígida.

Para aclarar el significado físico del tensor rotacional, consideremos una barra
rigida sometida a un giro w z , tal como se muestra en la Fig.2.4. Consideremos
que en los extremos de la barra se localizan los puntos P y Q, separados una
distancia diferencial.

El vector desplazamiento relativo del punto Q respecto al punto P, es :

68 = dU2 + dVJ (2.26)


2.3. INTERPRETACiÓN FíSICA DEL TENSOR DEFORMACiÓN 57

El desplazamiento du se puede calcular como:

du dR cos( W z + O) - dR cos O
dR( cos W z cos O- ser/'wzsenO) - dR cos O
du
(cos W z - 1) cos O - senw z senO (2.27)
dR
dv dRsen(w z + O) - dRsenO
dR( senw:; cos O+ cos W z senO) - dRsenO
dv
senwz cosO + (cosw z -1)senO (2.28)
dR
Las ecs. 2.27 y 2.28 se pueden expresar en forma matricial como:

I ~ 1= [ co~;~~ 1 c~s~~z1] I ~~~~ I (2.29)

Si W z es pequeño, entonces senwz = Wz y cos W z = 1, por lo que la eco 2.29


se reduce a:

1 1~
áv
áR
I = [OW z -w
O
z ] I cos
senOO I (2.30)

Obsérvese que esta última ecuaci6n puede ser representada como:


e = [nij]lel ,con Eij = O , por lo tanto nyX = -W z y nxy = W z • Así los
términos del tensor rotacional nij representan los giros del elemento alrededor
del eje z.

2.3 Interpretación física del tensor


deformación
El tensor deformaci6n qued6 definido por:

r ~(8U+8V)
~~au+aw~
8u
-
8x 2 8y 8x 2 8z 8x

~2
8v 1 8v 8w
[Eij ] =
8x + 8y aUj -
8y
-2 - +-
8z 8y
(2.31)
1 8w
-2 - +
8u
- ~ (8v + 8W) 8w
8x 8z 2 8z 8y 8z
A continuaci6n se da la interpretaci6n física de cada uno de los elementos del
tensor deformaci6n.
58 2. ESTADO DE DEFORMACiÓN

El térmiJlo Du ~ 6u; 6u represent.a el desplazamiento del element.o dx en la


Dx 6x
8u
din-!cción x, mientras que 6x es la dimensión original del elemento. Así 8x
represent.a la deformación normal o longitudinal unitaria en la dirección x,
yen lo que sigue la designaremos como é xx . El primer subíndice representa
la dirección del desplazamiento del punto Q a Q/, en tanto que el segundo
represent.a la dirección de los puntos P a Q, dada por el vector unitario e.

··1ar (eslgnaremos
De manera slml l· a1 '·
os termmos 8v
8y = é yy Y 8w
8z = é zz ,
1os
cuales representan las deformaciones normales o longitudinales unitarias
en las direcciones y y z, respectivamente. Los valores de é serán positivos
cuando los lados de la partícula aumentan de longitud, en caso contrario
serán negativos.

, .
En 1o que respecta a 1 termmo (8U BV) ,a1 mterpretaclón
By + Bx . . f'·
lSlca se h ar á con
referencia a la Fig. 2.5. Cosidérese un elemento diferencial cuya configuración
no deformada está definida por los vértices O, P, Q, y R, fijo en O. Una vez
que el elemento ha sido deformado, los vértices en la configuración deformada
son PI, Q/, RJ.

y
ff dy !- R'

~dy
1
T Q rQ' -- -- ----1'
R 1
y
81 /
/ 1
/
1

-
/

dy / /
/
P' 1
/
1/ -e 1 av dx
ox
O "'- 2
P
----j 1-
X

~dX
·
F; 19ura 2.5: 1nterpretaclOn
. , fí· ' . 8u
Slca de ltermmo 8y 8v
+ Bx·

Con base en esta figura, se tiene:


2.3. INTERPRETACION FISICA DEL TENSOR D ... 59

BU
-y d
!!JL BH
tan (jI (2.32)
dy By
Bv
dx Bv
tan (j2 ~ (2.33)
dx Bx

Por otra parte si el ángulo (j es pequeño, entonces:

Bu + -Bv}
{ -By Bx
~ (jI + (j2 (2.34)

.
En 1o que slgue d '
eSlgnaremos al '
térmmo { Bu
By Bv }
+ Bx = 'i'yx,vemos que
éste representa el cambio angular entre los lados dx y dy Y se conoce como
deformación angular o deformación cortante.

Analizando las características de deformación en los planos xz y yz, se puede


establecer por analogía:

BW + BU}
{ Bx Bz

'i'zy =
BW + BV}
{ By Bz
siendo 'i'zx Y'i'zy los cambios angulares entre los lados dx y dz y los lados dy
y dz, respectivamente. Los valores de las deformaciones angulares 'i', serán
positivos cuando decrece el ángulo que forman los lados correspondientes, en
caso contrario serán negativos.

De esta manera el tensor deformación puede ser representado como:


1 1
Cxx 2'i' yx ~'i'zx
1
[Eij ] = Cyy (2.35)
~'i'xY 2'i'ZY
1
2'i'xz 2'i' yz Czz

Así, en el tensor Eij aparecen:

a) En su diagonal principal, las deformaciones longitudinales unitarias en


tres direcciones ortogonales.
60 2. ESTADO DE DEFORMACiÓN

h) SimHricament.e a la diagonal principal aparecen la mitad de las


deformaciones angulares entre tres lados originalment(~ ortogonales.
Ejemplo 2.1
En el coutinuo representado en la Fig. 2.6, se est.ablece el siguiente vector
desplazamient.o:

;5 = (-Oyz) i + (Oxz)J

/ -1 ;T Z
\...
~ V
/,/',
/'
/'
/'

X
X

Y
Figura 2.6: Barra cilíndrica sometida a un campo de desplazamientos.

Las componentes del vector desplazamiento son:

u = -Oyz; v = Oxz; w = O; O = cte

Los elementos del tensor deformación se calculan como:

8u = O 8v = Oz
8x 8x
81¿ = -Oz 8v = O
8y 8y
8u = -Oy 8v = Ox
8z 8z
810 810 = O
-=0
8x 8z

Por lo tanto:
2.3. INTERPRETACION FISICA DEL TENSOR D ... 61

k
Tot"8 a
aJ aJl E?ff
+k (avax _au)
ay
Tot"8 I (O - Ox) - J (O + Oy) + k (Oz + Oz)
Tot"8 I (-ex) - J (Oy) + k (20z)

Además:
Oy
O O
O O
oi
2
Oy Ox
O
2 2

e
a) Cinemática de la partícula (origen).
x = O, Y = O, z = O
Se tiene: "8 = O, Eij = O, roM = o.
De lo anterior se concluye que:
No hay traslación, ni rotación ni cambio de forma.
Analizando la cinemática de la partícula A, en z = O.
En este caso = O; V = O; w = O ~ "8 = o.
tL

rot"8 = - RoI ::::} la partícula gira alrededor del eje x en dirección negativa.
El tensor deformación vale:

Eij = [gO ORg ik]


O
::} hay cambio de forma

b) Manejando el mismo tensor, describa el movimiento de las partículas A,


B, e en las caras z = O Y z = L.
b.l) Establezca claramente la rotación de cada una de esas partículas.
b.2) Defina los cambios de forma y dimensiones en las partículas propuestas.
En z = O, se tiene:

"8 = (-Oyz)I + (Oxz)J


Por lo tanto, "8 = o. Este resultado implica que las partículas no experimentan
desplazamientos paralelos a ningún eje.
62 2. ESTADO DE DEFORMACiÓN

Para el rotacional se tiene:

1'Ot8 (-Ox) Z+ (-Oy)] + (20z) k


l'ot8 (-Ox)z-(Oy)]
El tensor deformación vale:

O O
Eij = [ O O
-y]
x
1
-O
-y x O 2

Si se compara el vector rotacional con el vector de posición de cada punto,


descrito por:

r = xi + yj
Se puede ver que ambos vectores tienen las mismas líneas de acción pero
difieren en dirección en 1800 • Además, las magnitudes de los vectores rot8
son proporcionales en O, a la magnitud de los vectores de posición.
Lo anterior indica que todos los vectores rot8 coinciden en el origen. Por lo
tanto,
[rot81.
u J ortgen = O ====} no hay giro
Analizando la sección z = 1, se tiene:

8 = (-Oy) Z+ (Ox)]
Ya que w = O, la secclon permanece plana, y sólo se presentan
desplazamientos en el plano de la sección.

y
Figura 2.7: Componentes del vector desplazamiento en la sección z = 1.

Dado que

r=xi+yj y 8= (-Oy)z+(Ox))
2.4. COMPONENTES NORMAL Y CORTANTE DEL VECTOR DEFORMACiÓN 63

Los desplazamientos son proporcionales a la magnitud del vector de posición


r y ortogonales a él, por lo que podemos escribir.

En consecuencia, la sección z = 1 permanece plana y gira rígidamente


alrededor del eje z, un ángulo O.
El ángulo O mide el giro respecto a z, de una sección distante z = 1, de la
sección que no se mueve.
Además,

[rotJ]Z=l = (-Ox)I - (Oy)) + (20) k


Este vector es igual al de la sección z = O, más un vector 20k, que implica
una rotación rígida alrededor de z = O.
Además de estos movimientos rígidos de las partículas, se cumple que

Finalmente, en z = L

La sección permanece plana, girando un ángulo OL respecto al eje z.

[rotJ]z=L = [rot8J z=o + (20L) k


Así, las rotaciones de las partículas son similares a las de z = O, pero son
una rotación rígida alrededor de z, con magnitud igual a OL.
Por lo tanto, 5 = (-Oyz) I+(Oxz)}, describe la torsión de una barra cilíndrica
cuando las secciones permanecen planas y la sección z = L gira un ángulo
OL respecto a la primera sección.

2.4 Componentes normal y cortante del


vector deformación
Previamente se estableció que el vector desplazamiento relativo de un punto
Q respecto a un punto P, se puede expresar como la diferencia vectorial
6.5 = 5Q - 5p. En la Fig. 2.8 se muestra gráficamente dicho vector. Se puede
observar en dicha figura que si P y Q son dos puntos en la configuración
no deformada del cuerpo, cuya dirección está dada por el vector unitario e,
entonces 6.5 es el vector que va del punto Q al punto Q', ubicado este último
en la configuración deformada del cuerpo. De esta manera, y con referencia
a la Fig. 2.8 se puede establecer la siguiente suma vectorial:
64 2. ESTADO DE DEFORMACiÓN

/
/
~~ ,/
/
/

e // "'). e
••_ _~~~_ _ _ _ _ _ _ .~
/ _ _- L_ _ _ _~~

p Q ~~n
Figura 2.8: Componentes normal y cortante del vector desplazamiento
relatIvo.

(2.36)

Dividiendo esta última ecuación entre el módulo del vector PQ, se obtiene:

1:18 l:1ó l:1óo


- - = - -n + - - (2.37)
IPQI IPQI IPQI
La eco 2.37 se puede escribir como:
(2.38)

siendo E'n el vector deformación normal o longitudinal y fQ el vector


deformación cortante o desviador. Obsérvese que la eco 2.38 tiene la misma
forma que la establecida para el vector esfuerzo In-

Dado que el vector deformación E' tiene la misma estructura matemática que
el vector esfuerzo In, todo lo establecido con anterioridad en el estudio del
estado de esfuerzo será valido para el estado de deformación.

2.5 Deformaciones principales y


direcciones principales de deformación
Para determinar las direcciones principales de deformación se hace un
planteamiento análogo al realizado para el caso del estado de esfuerzos.

Las magnitudes de las deformaciones principales quedan definidas por la


ecuación característica siguiente:
(2.39)
2.6. SIGNIFICADO FíSICO DEL INVARIANTE LINEAL JI 65

donde:
JI = Exx + Eyy + Ezz; Invariante lineal
J- 1 2 1 2 1 2 .
2 - ExxEyy + EyyEzz + EzzExx - ¡/'xy - ¡/'yz - ¡/'zx'
Invariante cuadrático
J 3 = I[Eij JI; Invariante cúbico

Por lo tanto el tensor deformación principal se puede expresar como:

(2.40)

Siendo El,E2, Y E3 las deformaciones principales mayor, intermedia y menor,


respectivamente.
Por analogía se puede llegar a establecer un elipsoide de Lamé que describe
el conjunto de vectores En.
Esto último implica la existencia de direcciones principales de deformación,
1, 2, 3; en dichas direcciones sólo se generan cambios de dimensiones lineales
pero no angulares.
Las direcciones principales de deformación se calculan con el modelo
matemático de valores y vectores característicos, de manera similar a lo que
se hizó para el estado de esfuerzo.
Diremos que dos estados de deformación son iguales si sus invariantes
correspondientes son iguales.

2.6 Significado físico del invariante lineal JI

Considérese un elemento de volumen de un medio continuo sometido a un


estado de esfuerzos tal que sus dimensiones se reducen pero el cuerpo
mantiene su forma (Fig.2.9). Por facilidad de presentación, supondremos
que las dimensiones originales del elemento son unitarias.
66 2. ESTADO DE DEFORMACiÓN

Figura 2.9: Cambio volumétrico de una partícula de un medio continuo.

El cambio de volumen b. V que experimenta el cuerpo resulta igual a:


b.V = (1 + él) (1 + é2) (1 + é3) - (1) (1) (1)

Desarrollando:

b. V = (él + é2 + é3) + (élé2 + é2é3 + é3él) + (é¡é2é 3)


Teniendo presente la definición de invariantes de deformación, resulta que:
(2.41)
La suma de los tres invariantes mide el cambio de volumen de una partícula
unitaria. Dado que en muchos problemas de ingeniería las deformaciones
longitudinales son menores de 1%, se acepta que:
(2.42)

2.7 Componentes volumétrica y desviadora


del tensor deformación

Supongamos que en un punto de un medio continuo se concoce el tensor


deformación,Eij, en un marco de referencia principal. Este tensor, al igual
que en el caso del tensor esfuerzo, puede ser descompuesto en su partes
volumétrica y distorsional o desviadora. Así,
(2.43)
2.8. ECUACIONES DE COMPATIBILIDAD DE DEFORMACIONES 67

donde:

E = JI [
u 3
6O OO1 ~ ] (2.44)

JI
101 - - O O
3
JI
Eo = O 102 - - O (2.45 )
3
JI
O O 103 - -
3
Eu mide el cambio de volumen de las partículas unitarias del cuerpo
deformable, mientras que Eo mide la distorsión de las partículas unitarias.

2.8 Ecuaciones de compatibilidad de


deformaciones
Supongamos que de alguna manera se establecen las dos ecuacIones
diferenciales siguientes, siendo la variable desconocida la componente
del vector desplazamiento u(x, y):

(]u
x+3y
Bx
Bu
(2.46)
By

El análisis de estas dos ecuaciones nos lleva a que no pueden ser resueltas
debido a que son inconsistentes. El cálculo de la derivada B~~y en ambas
ecuaciones pone en evidencia dicha inconsistencia. La derivada de la primera
ecuación da 3 mientras que la segunda da 2x.
Si las ecs. 2.46 se expresan como:

Bu
f (x, y)
Bx
Bu
9 (x , y) (2.47)
By

Entonces, para que estas ecuaciones puedan ser integradas se debe de cumplir
la siguiente condición:
68 2. ESTADO DE DEFORMACiÓN

DI = ()g
(2.48)
f)y ax
La eco 2.48 represent.a. la condición de int.cgTabilidad d<-! Riemann o ecuación
de compat.ibilidad.

Por lo expuesto, es claro que los elementos del tensor deformación deben de
gmu·dar cierta relaci6n ent.re sí para asegurar la exist.eucia de continuidad del
campo de desplazamientos (5 correspondient.e.

Para establecer la.<; ecuaciones de compat.ibilidad recordemos la definición de


cada. Hno de los element.os del tensor deformación Eij.

f)u
ex = (2.49)
f)x
f)v
ey (2.50)
f)y
f)w
e: (2.51)
f)z
f)u OV
Ixy -f)y +-
OX
(2.52)
ou ow
Ixz -oz + -ox (2.53)
ov ow
Iy: = -
oz
+oy
- (2.54)

Para relacionar a. los componentes de deformación hagamos las siguientes


operaciones.
Derivando la eco (2.49) respecto a y2 , se obtiene:

(2.55)

Haciendo lo mismo con la eco (2.50) pero ahora respecto a X2 , se llega a:

éPv
= (2.56)

Derivando la eco ( 2.51) respecto a x e y.


0 2 , xy [f3u 03 v
-- = + -::--.,-----::- (2.57)
axf)y oyaxoy axoyox
2.8. ECUACIONES DE COMPATIBILIDAD DE DEF. oo. 69

Si u y v son funciones continuas

(2.58)

Análogamente, si u y w son funciones continuas,

(2.59)

.Finalmente si v y w son funciones continuas,

(2.60)

Estas tres últimas ecuaciones forman el prImer grupo de ecuaciones de


compatibilidad.
Ahora, estableceremos las relaciones entre los valores de las deformaciones
angulares , que permitan definir explícitamente a las deformaciones
longitudinales €.
Derivando la eco (2.52) respecto a z,

O'Xy 02U 02V


-=--+-- (2.61)
OZ oyoz oxoz
La eco (2.53) respecto a y,

(2.62)

Finalmente, la eco (2.54) respecto a x,

(2.63)

Multiplicando la eco (2.63) por (-1) Y sumando miembro a miembro con las
ecS. (2.61) y (2.62), se llega a:

2 02u = o,xy _ o,YZ + o,zx (2.64)


oyoz oz ox oy
Derivando esta última ecuación respecto a x se obtiene:
70 2. ESTADO DE DEFORMACiÓN

(2.65)

De manera similar se puede obtener:

2
2 Dé y = ~ (D"/yz _ D,,/zx + D"/x y ) (2.66)
DxDz Dy Dx Dy Dz
2
2 Dé z = ~ (D"/yz + D,,/zx _ D"/x y ) (2.67)
DxDy Dz Dx Dy Dz

Las ecs. (2.65), (2.66) Y (2.67) forman el segtrndo grupo de ecuaciones de


compatibilidad. Si estas 6 ecuaciones son satisfechas se dirá que los elementos
de Eij son compatibles.
Cesaro (1905) demostró que las 6 ecuaciones previas garantizan que se
satisfagan las condiciones de Riemann para establecer las condiciones de
integrabilidad que permiten la aceptación de la existencia del campo (5 =
Ut + V] + wk, a partir de un tensor Eij compatible.
Ejemplo 2.2
Compruebe que el campo de desplazamientos (5 = (-Oyz) Z+ (Oxz)], siendo
O una constante, define un tensor esfuerzo compatible.
Los elementos del tensor esfuerzo son:

Eij
1 [OO
= -2 O O
O-y]
x
-y x O
Aplicando los dos grupos de ecuaciones de compatibilidad, se tiene:

D2,,/xy D2é x + D
__
2
__ éy
= ec.2.58
DxDy Dy2 DX2
O = 0+0
2
D ,,/xz D2é x D2é z
= ec.2.59
DxDz Dz 2 + DX2
O 0+0
Por lo tanto se satisface el primer grupo de ecuaciones.
2.8. ECUACIONES DE COMPATIBILIDAD DE DEF. .. , 71

82'Yy=
8y8z
O
Para el segundo grupo de ecuaciones se tiene:

2
28 éx ~ (_ 8'Y yz + 8'Yzx + 8'YXY) ec.2.65
8y8z 8x 8x 8y 8z
8
O 8x (- 2 - 2 + O) = O

Para las otras dos ecuaciones se obtiene:

2
28 éy ec.2.66
8x8z
O
2
28 éz
ec.2.67
8x8y
O

Se satisface el segundo grupo de ecuaciones, por lo que en consecuencia Eij


es compatible y existe.
72 2. ESTADO DE DEFORMACiÓN
Capítulo 3
Principios Generales de la
Mecánica

Estos principios constituyen las leyes fundamentales que rigen el


comportamiento mecánico de los medios continuos, siendo expresados muchas
de las veces como leyes de conservación de ciertas cantidades físicas. Tal es
el caso del principio de la conservación de la masa o ecuación de continuidad,
el de la conservación de la cantidad de movimiento, el de la conservación de
la energía (primera ley de la termodinámica) y el de aumento de entropía
(segunda ley de la termodinámica).

En algunos fenómenos mecánicos es frecuente ignorar algunos de estos


principios por considerar que sus efectos son despreciables, lo que permite
una formulación matemática del fenómeno más simple.

Las leyes de conservación mencionadas pueden ser aplicadas a un cierto


volumen de materia teniendo como frontera una superficie cerrada de forma
cualquiera. En el desarrollo del modelo de comportamiento es usual encontrar
que ciertas cantidades físicas aparecen como integrales de superficie y otras
como integrales de volumen. La transformación de una integral de volumen a
una de superficie y viceversa es una operación matematica requerida en esta
formulación. Esta transformación matemática es conocida como el teorema
de Green o de la divergencia, el cual establece que para una función espacial
f (x, y, z), continua y derivable, con primeras derivadas parciales también
continuas, se cumple que:

(3.1)

siendo ni los cosenos directores de la normal a una superficie cerrada A,


frontera de un volumen V , en el entorno a un punto definido por dA.

73
74 3. PRINCIPIOS GRALES. DE LA MEC . ...

3.1 Principio de la conservación de la masa o


ecuación de continuidad
Este principio establece que en el interior de un "volumen de control",
entendido éste como un elemento diferencial asociado a un sistema
de referencia fijo en el espacio, la masa no se crea ni se destruye. De esta
manera la existencia de cambios de masa en tal volumen de control tendrán
que estar asociados a un flujo de masa a través de la superficie de control.
Con referencia a la Fig.3.1, y suponiendo que la densidad, p, del medio llena
todo el volumen, V, la masa total, M, ocupada por dicho volumen en un
t.iempo t. , result.a:

M= ¡PdV (3.2)

n
1
/

Figura 3.1: Volumen de referencia.

Dado que la densidad del medio es una función de posición y del tiempo, ésta
se puede expresar como:

p=p(x,y,z,t) (3.3)
Por lo tanto, la rapidez de variación de la masa total respecto al tiempo, en
el volumen V, está dada por:

8M =
8t
r8p dV
Jv 8t (3.4)

Considerando que dentro del volumen V, la masa no se crea ni se destruye,


entonces la eco 3.4, es equivalente a la rapidez de variación del flujo de masa
hacia el interior del área A.
Por otra parte, el flujo de masa hacia el exterior del área dA en el entorno del
punto "P", es pvndA, siendo V n la componente normal del vector velocidad
V. Así, la rapidez de variación del flujo de masa total es:
3.1. PRINCIPIO DE LA CONSERVACION DE LA MASA ... 75

¡ (-pvn)dA = - ¡ p(V . n)dA (3.5)

En esta ecuación el signo menos obedece a que al entrar flujo el vector


velocidad va en sentido contrario de la dirección de la normal n a la superficie,
y puede ser expresada, de acuerdo con el t.eorema de Green, como

¡ (-pvn)dA =- ¡ dív(pv)dV (3.6)

Dado que la ecs. 3.4 y 3.6 representan el mismo fenómeno, se tiene:

aM
8t
-¡ dív(pv)dV

( ap dV (3.7)
lv at
Reordenando esta última ecuación se tiene:

¡ [~ + díV(PV)] dV = O (3.8)

La eco 3.8 se debe de satisfacer para cualquier volumen V, por lo que el


integrando necesariamente t.e ndrá que ser nulo, esto es:

ap + dív(pv) = O (3.9)
at
La eco 3.9 se puede representar en notación indicial como:

ap + a(pVi) = O (3.10)
at aXi
Desarrollando cada uno de los términos de la eco 3.10, se tiene:

(3.11)

Por otra parte:


dp ap ap
-=-+Vi- (3.12)
dt at aXi

Sustituyendo las ecs. 3.11 y 3.12 en la eco 3.10, se obtiene:

dp aVi _ O (3.13)
dt + p aXi -
76 3. PRINCIPIOS GRALES. DE LA MEC . ...

La eco a.13 se puede escribir en un sistema de referencia cartesiano como:

o
o (3.14)

Despejando la divv de esta últ.ima ecuación se tiene:

. _ 1dp
dwv = --- (3.15)
pdt

De la eco 3.15 se puede observar que si el medio es incompresible la divv = o.


Así en un sistema de referencia cartesiano se tiene que:

(3.16)

Esta ecuación diferencial representa el principio de conservación de la masa


y se conoce también como ecuación de continuidad.

3.2 Principio de conservación de la cantidad


de movimiento
La rapidez de variación con respecto al tiempo de la cantidad de movimiento
de un sistema de partículas es igual al vector fuerza resultante de todas las
fuerzas ext.ernas, act.uando sobre el conjunt.o de part.ículas, siempre y cuando
sea la t.ercera ley de Newton (acción y reacción) quien gobierne las fuerzas
int.ernas en el sist.ema.
En relación con la Fig.3.2, dicho principio queda expresado como:

f t"dA + f p!dV = ~ f pvdV (3.17)


lA lv dt lv

La eco 3.17 se puede expresar en notacion indicial como:

1 +¡tidA pfi dV = :t ¡ pVidV (3.18)


3.2. PRINCIPIO DE CONSERVACION DE LA CANTIDAD... 77

Figura 3.2: Conservación de la cantidad de movimiento.

Obsérvese que el término del lado derecho de la eco 3.18 se puede expresar
como:

(3.19)

Por otra parte, en el estudio del estado de esfuerzo se estableció que:


(3.20)

Sustituyendo la eco 3.20 en la 3.18 y aplicando el teorema de Green se obtiene:

Reordenando términos, se llega a:

Iv [:~ + PJi] dV Iv P~~i dV


= (3.21)

Finalmente, el principio de conservación de la cantidad de movimiento


conduce a:

(3.22)

Dado que la eco 3.22 se debe cumplir para todo volumen V, entonces el
integrando debe ser igual a cero, esto es:
87ij dVi
-+pJi-P- =0 (3.23)
8xj dt

Esta ecuación se conoce como la ecuación del balance de la cantidad de


movimiento o ecuación de Cauchy.
78 3. PRINCIPIOS GRALES. DE LA MEC . ...

En el caso de equilibrio estático la aceleración dV = 0, por lo que la eco 3.23


dt
se reduce a:

(3.24)

La eco 3.24 representa un sistema de tres ecuaciones diferenciales parciales


donde las incógnitas son los nueve elementos del tensor esfuerzo, que por
simet.ría del mismo, bast.ará con conocer seis elementos de dicho tensor.
Es obvio que el problema es estáticamente indeterminado, por lo que será
necesario incluir ecuaciones adicionales, por ejemplo, aquellas que relacionen
los esfuerzos con las deformaciones de un material en particular. Dichas
relaciones reciben el nombre de ecuaciones constitutivas, las cuales se
estudiarán en el capítulo 4, para el caso de los materiales elásticos lineales,
homogéneos e isótropos.

3.3 Princi pio de conservación de energía; ·


primera ley de la termodinámica
El principio de conservación de energía es una consecuencia de la primera
ley de la termodinámica, él establece que la energía no se crea ni se destruye
sólo se transforma.
Esta ecuación de energía involucra una incógnita adicional, la energía interna,
por lo que su utilidad radica en poder relacionar dicha energía interna con
alguna variable de estado.
En la mecánica de los medios continuos un sistema termodinámico se define
como una porción de materia continua, donde no existe intercambio de
mat.eria con cuerpos vecinos, lo que se ha dado en llamar un sistema cerrado.
Las superficies frontera del sistema se mueven en general con el flujo de
materia.
La rapidez de variación del trabajo realizado por las fuerzas de superficie y
de cuerpo sobre un sistema termodinámico, se puede expresar como:

w 1 t n ·vdA + ¡p]. vdV = 1 tiVi dA

+ ¡pjiVidV (3.25)

Sustituyendo en la integral de superficie el valor de ti = ~jnj, y aplicando el


teorema de Green, se tiene:
3.3. PRINCIPIO DE LA CONSERVACION DE ENERGIA... 79

(3.26)

Sustituyendo la eco 3.26 en la 3.25, se llega a:

(3.27)

Obsérvese que el término entre paréntesis, por la ecuación de movimiento,


resulta igual a:

(
8Tij
8x· + p,
J.) -_P 8t 8
V
i
(3.28)
J

De esta manera la rapidez de variación del trabajo, W, es igual a:


. = -d
W dt v -pv
2 ••
¡
1 ·v·dV + J T,. ·V·
'J "J
. dV (3.29)

Por otra parte:


(3.30)
siendo Dij el tensor rapidez de deformación y W i j el tensor vorticidad. Este
último es un tensor antisimétrico, esto es, W ij = - W j i
El producto tensorial que aparece en la segunda integral de la eco 3.29 es:
Tj i (D i j + W ij ) = Tj i Dij + T i jWi j (3.31)
En la eco 3.31 el término T ji W ij es igual a cero, por lo que:
TjiVi ,j = TijDi j (3.32)
De esta manera la eco 3.29 queda como:

W. d
= -dt 1 ¡
v -pv·v·dV
2 • • + J T, .D· ·dV
'J 'J
(3.33)

La primera integral de esta última ecuación representa la energía cinética


del sistema en tanto que la segunda representa la rapidez de variación de la
energía interna total, por lo tanto:

(3.34)

ESTA TESIS NO SALl


OE LA BIBIJOTECA
80 3. PRINCIPIOS GRALES. DE LA MEC . ...

El principio de la conservación de la energía establece que la variación de la


energía cinét.ica más la energía interna por unidad de tiempo es igual a la
~uma de la variación del trabajo mi:Í..c; cualquier otra energía suministrada o
ext.raída por unidad de tiempo en el sistema termomecánico.

Definiendo al vector q como fiujo de calor por unidad de área y tiempo en


el fenómeno de conducción calorífica y a r la constante de radiación de calor
por unidad de masa y tiempo, entonces la rapidez de aumento de la cantidad
de calor en el medio se puede expresar como:

(3.:35)

Para un medio continuo termomecánico es costumbre expresar la variación


de la energía interna total por unidad de tiempo como una función de la
energía específica interna u por unidad de tiempo, por lo que Ü, queda como:

Ü = :t [ pudV = [PüdV (3.36)

Aplicando el principio de la conservación de la energía, se tiene:

K+Ü=W+Q (3.37)
Sustituyendo las ecs. 3.33, 3.35 y 3.36 en 3.37, se obtiene:

d
-dt¡ (1
v -pvovodV
2 ••
\
+ pu ,o) dV = !!: r~pvo vodV + fr.o')'D')o°dV
dt lv 2 "

-1 q. TídA + [prdV (3.38)

Tomando en cuenta el teorema de Creen, el integrando q . TídA, se puede


expresar como:

r qjnjdA lvr 8xj


lA
8% dV
= (3.39)

De esta manera la eco 3.38 queda como:

[(;:: +P~~ -7ijD j +pr)dV=ü


i (3.40)

Para un volumen arbitrario V dentro del medio continuo, el integrando de la


eco 3.40 debe ser nulo, por lo que:

(3.41)
3.4. SEGUNDA LEY DE LA TERMODINÁMICA, DESIGUALDAD DE
CLAUSIUS-DUHEM 81

Esta última ecuación se conoce como la ecuación de la conservación de la


energía o primera ley de la t.ermodinámica .

3.4 Segunda ley de la termodinám.ica,


desigualdad de Clausius-Duhem

Esta ley establece que la variación con respecto al tiempo de la entropía total,
s, en un medio continuo de volumen V, siempre es mayor que la suma del
fiujo de entropía que entra a través de la superficie del medio, A, más la
entropía creada interiormente a causa del propio cuerpo. Esta ley se puede
expresar desde un punto de vista matemático en forma integral como:

:!:.- { psdV;;:: { p":..d\l - (:idA (3.42)


dt lv 1v e lA e
siendo e una función de estado denominada la temperatura absoluta.
En la eco 3.42 el signo " = " corresponde a procesos reversibles, el signo
" > " a procesos irreversibles y el signo " < )) indica que el proceso es no
factible.
82 3. PRINCIPIOS GRALES. DE LA MEC . ...
Capítulo 4
Elasticidad Lineal
Es la parte de la mecánica de los medios continuos que estudia el
comportamiento de sólidos cuyas propiedades del material son independientes
del tiempo.
Para establecer la relación entre las fuerzas que actuan en un continuo y los
desplazamientos que provocan, es necesario relacionar:

Fuerzas -t 1'.¡j ~ Eij ~ Desplazamientos


1
Relaciones constitutivas

4.1 Planteamiento matemático para


definir las relaciones constitutivas en
un continuo cualquiera

y t =t F2
Estado final
"B"
x
Figura 4.1 : Relaciones entre esfuerzos y deformaciones en un medio
continuo elástico lineal.

83
84 4. ELASTICIDAD LINEAL

Con~idérese que el estado de esfuerzos y deformaciones de un medio continuo,


corrt!spomle al estado A o inicial, dado por:

Una vez que el medio continuo se somete a un sistema de cargas se generan


e n el mismo un estado de f'-sfuerzos y deformaciones, que corresponden al
estado B o final, dado por:

Al aplicar las fuerzas {Pi} se establece un campo de desplazamientos (5 =


uI + vJ + wk . Este campo genera en cada punto del continuo a un tensor
[EijJn- Además, las fuerzas {Pi} generan en el continuo el tensor esfuerzo
[TijJn-
Para estudiar la dinámica del continuo, el volumen estará ocupado por un
material con relaciones constitutivas bien definidas, de manera que se pueden
relacionar los tensores Eij y T ij .
Sabemos que:

Ya que físicamente existe relación entre Eij y Tij, se puede escribir

(7' xx <pd Gxx , Gyy , Gzz, "fxy' "f xz ' "fyz )


{7'yy <P2 ( Gxx, Gyy, Gzz, "fxy' "f xz ' "fyz )
{7'zz <P3 ( G=, Gyy, Gzz, "fxy' "f xz ' "fyz )
Txy 'ljJd Gxx, Gyy, Gzz, "fxy' "f xz ' "fyz )
T xz = 'l/J2( Gxx, Gyy, Gzz, "fxy, "fxz, "fyz )
T yz 'l/J3( Gxx, Gyy , Gzz, "fxy' "f xz ' "fyz ) (4.1)

Todas las funciones <p y 'l/J deben ser funciones continuas.


Las funciones anteriores definen a las relaciones constitutivas de los medios
continuos.
Las relaciones constitutivas que se seleccionan en la mecánica del medio
continuo son relaciones probadas experimentalmente, que permiten describir
de manera razonable el comportamiento real de los materiales.
Se considera que el tiempo transcurrido entre el estado A y el estado B es
una cantidad pequeña, por lo que sería posible establecer, siguiendo las ideas
de continuidad que:
4.1. PLANTEAMIENTO MATEMÁTICO PARA DEFINIR LAS RELACS . ... 85

[o-xxlB [o- xx lA + do-xx


[O- yylB [o-yylA + do-yy
[o-yylB [o-yylA + do-yy
[Txyl B [Txyl A + dTxy
[Txzl B [Txzl A + dT xz
[Tyzl B [Tyzl A + dTyz (4.2)

Desarrollando las diferenciales totales de las ecs. 4.2 se tiene:


86 4. ELASTICIDAD LINEAL

Estas ecuaciones pueden ser escritas como:

[O"xxJA + Clléxx + C21 é yy + C31 é zz + C 4I'YXy


+G'sI'Yxz + C6 I'Yz y
[O"yyJ A + C 12 é xx + C 22 é yy + C 32 é zz + C421XY

+C521XZ + C621yZ
[O"zzJA + Cl3 éxx + C 23 é yy + C3Jézz + C431Xy
+C531xz + C 631yz (4.3)
[Tyz JA + C14é xx + C 24 é yy + C 34 é zz + C441XY
+G's4/xz + C641yz
[TyzJ A + C 15 é xx + C 25 é yy + C 35 é zz + C 451XY
+C45 I xz + C 65 I yz
[TyzJ A + C 16 é xx + C 26 é yy + C 36 é zz + C461XY

+C561xz + C 661yZ

Donde los coeficientes Cij representan operadores diferenciales respecto al


tiempo.
Obsérvese que en las expresiones planteadas aparecen 36 operadores
diferenciales Cij .
El desarrollo anterior indica que para definir a [TijJ B es necesario conocer los
tensores [7ijJ A , [Eij]B y los 36 operadores diferenciales Cij .

Estas relaciones indican que es posible seguir la siguiente secuencia:

(5 --t [EijJ B - t Relaciones constitutivas -t [TjkJ B - t {Fi}


Por lo tanto para relacionar a los desplazamientos que se generan en el
continuo, es necesario conocer a 36 operadores diferenciales respecto
al tiempo, lo que conduce a definir a las fuerzas {Fi} que provocaron la
aparición de (5 .
Para definir a esos 36 operadores es necesario analizar pruebas experimentales
en diversos materiales y observar las características de su respuesta, ésta
deberá compararse con formulaciones teóricas que permitan predecir dicha
respuesta con suficiente aproximación.
En Ingeniería es común establecer ciertas hipótesis que permitan simplificar
las formulaciones matemáticas.
Las hipótesis comunes son:
10 El material que ocupa el continuo es homogéneo (sólo existirán 6
relaciones) .
4.1. PLANTEAMIENTO MATEMÁTICO PARA DEFINIR LAS RELACS . ... 87

2° El material es isótropo.
3° Las direcciones principales de esfuerzos coinciden con las direcciones
principales de deformación.
Tomando en cuenta las hipótesis anteriores, supongamos que en un punto del
medio continuo se conocen los tensores esfuerzo y deformación en un sistema
de referencia principal:

Para este caso particular las relaciones constitutivas se reducen a:

(jI Cllél + C 2l é 2 + C 3l é 3
(j2 C 12 é l + C 22 é 2 + C 32é 3
(j3 = C 13 é l + C 23 é2 + C 33 é3 (4.4)

Obsérvese que la tercera hipótesis reduce a nueve el número de operadores


diferenciales necesarios.
Cambiando los ejes 2 y 3 por 3' y 2', al aplicar las ecuaciones constitutivas
anteriores, se tendrá:
(4.5)
Debido a la indiferencia respecto al marco de referencia, se debe tener que:
C 2l = C 3l (4.6)
Por lo anterior, la ecuación 4.5 queda como:

(jI Cllél + C2l (é2 + é3)


(jI (Cll - C 2l )él + C 2l (él +é2 +é3) (4.7)
llamando a
Cll - C2l = 2G (4.8)
y tomando en cuenta que JI = él + é2 + é3, entonces la eco 4.8 queda como:
(4.9)

siendo A Y G dos constantes elásticas conocidas como constantes de Lamé,


y representan a operadores diferenciales respecto al tiempo. Para (j2 y (j3 se
tiene:

2Gé2 + AJ1 (4.10)


2Gé3 + AJ1 (4.11 )
88 4. ELASTICIDAD LINEAL

Por las hipót.esis plant.eadas las ecs. 1.9 a 4.11, se pueden generalizar como:

(4.12)

De est.a manera el problema de la dinámica de continuos homogéneos,


isót.ropos y con coincidencia de direcciones principales de esfuerzos y
deformaciónes, se reduce a la búsqueda de dos operadores diferenciales >. y
G en lugar de 36.

4.2 Ecuaciones constitutivas de los materiales


elásticos lineales en un marco de
referencia principal
Considérese que en un punto de un medio continuo se establece el tensor
deformación:

(4.13)

El tensor esfuerzo correspondiente, tomando en cuenta la eco .4.13, resulta


igual a:

(4.14)

Estas expresiones pueden ser ligeramente transformadas, descomponiendo los


tensores esfuerzo y deformación en sus componentes volumétrica y desviadora,
esto es:

[E;j] [g c
O
er
8]
C3
= J¡ [ Ó O
3 O O 1
1 O]
O

Cl-- O O
3

+ O C2 - -
3
O
JI
O O C3 - -
3

Ahora las componentes del tensor ~j son:


4.3. ECUACIONES CONSTITUTIVAS DE LOS MATERIALES ELÁSTICOS LINEALES
HOMOGÉNEOS E ISÓTROPOS EN UN MARCO DE REFERENCIA CARTESIANO 89

Analizando las componentes de ambos tensores, y haciendo (3)' + 2G) = 3J{,


se concluye que:

(4.15)
(4.16)

Las ecs. 4.15 y 4.16 muestran una relación simple entre las componentes
volumétrica y desviadora de los tensores T;,j y Eij , a través de las constantes
J{ y G que reciben el nombre de módulo volumétrico y módulo de rigidez al
cortante, respectivamente.

4.3 Ecuaciones constitutivas de los materiales


elásticos lineales homogéneos e isótropos
en un marco de referencia cartesiano
Supongamos que en un punto del continuo se establecen los tensores Eij y
Jij en un marco de referencia cartesiano, esto es:

1 1
Exx "2'YX ~'zx
1 xx ;;:.x zx
Ejk = Tjk = [ U
T xy YY T zy ]
~'XY
Eyy
"2'ZY T xz T yz U zz
1
"2'xz "2'YZ Ezz

Tratemos de encontrar la forma que adquieren las relaciones constitutivas


para relacionar a esos tensores.
Para ello calculemos Un Y En, con el objeto de aplicar la relación constitutiva
previamente derivada para materiales isótropos, eco 4.12.
Por definición:
90 4. ELASTICIDAD LINEAL

Desarrollando:

(O" x:x + 7 yx COS (3 + 7 zx cos 7) ~


cos a
+ (7 xy cos a + O" yy cos (3 + 7zy cos 7 )]
+ (7 xz cos a + 7 yz COS (3 + O" zz cos 7) k (4.17)

El esfuerzo normal se calcula como:

-
t" . -n = O" xx cos a
2
+ O" yy cos (3 2 + O" zz cos 7 2 + 27xy COS o: cos (3
+27yz COS (3 cos 7 + 27zx COS 7 cos o: (4.18)

De manera análoga se puede establecer para la deformación normal én:

é xx cos a 2 + é yy COS (32 + é zz COS 7 2 + 7xy cos a cos (3


+7y z cos (3 cos 7 + 7zx cos 7 cos a (4.19)

Aplicando la eco 4.12, se tiene:

O"n )..)1 + 2Gén


O"n cos a 2 + O"yy COS (32 + O"zz COS 7 2 + 27xy COS o: COS (3
O" xx
+27 yz COS (3 cos 7 + 27zx COS 7 cos o:
= )..)1 (COS 0: + COS (32 + COS 7 2 )
2

+2G (éxx COS 0: + é yy COS (32 + é zz COS7 + 7xy cos o: COS (3


2 2

+7 yz cos (3 cos 7 + 7 zx cos 7 cos o: )

Para que exista igualdad se debe cumplir que:


O"xx = )..<b + 2Gé xx
O"yy = )..J1 + 2Géyy
O"zz = )..)1 + 2Gézz (4.20)
7xy = G7xy
7 yz = G7 y z
7 zx = G7zx

Estas ecuaciones reciben el nombre de ecuaciones constitutivas de los


materiales elásticos linelaes, homogéneos e isótropos.
Las ec.<¡. 4.20 se pueden expresar en forma matricial como:
[
u~~
u
U zz
T xy
T yz
1 _
-
l
4.3. ECUACIONES CONSTo MARCO DE REFERENCIA ...

A + 2G
O
O
A
A
A
O
O
A
A +A 2G A+2G
A

O
O
g g gj
G O O
O G O
[~~~l
YX;
Iyz
(4.21)
91

T zx O O O O O G Izx

Estas ecuaciones toman la forma:


[u] = [e] [E] (4.22)
siendo:
u]: matriz de esfuerzos
El: matriz de deformaciones
e]: matriz de constantes elásticas del medio
Las ecs. 4.21 pueden ser escritas en notación indicial como:
Uij = AÓijékk + 2Gé ij (4.23)
Para expresar las deformaciones en función de los esfuerzos, es necesario
invertir la matriz de constantes elásticas y operar matricialmente, esto es:

[E] = [Cj-l [u] (4.24)


En lugar de buscar invertir la matriz de constantes elásticas se seguirá el
siguiente procedimiento.
Las tres primeras expresiones de las relaciones constitutivas, eco 4.20 se
pueden escribir como:
u xx = (A + 2G)éxx + Aé yy + Aézz
U yy = Aéxx + (A + 2G) éyy + Aézz (4.25)
U zz = Aéxx + Aéyy + (A + 2G) ézz
Resolviendo el sistema para éxx , é yy Y ézz , empleando el método de Cramer,
se tiene:

(A + 2C) I (A +/G)
~ = A (A +A2G) A
A = A + 2G A
A A (A + 2G) (A + 2C)

+ (- A) I ~ (A ;\C) I+ A I ~ (A +/G) I
(A + 2C) [(A + 2G)2 - A2] - A [A (A + 2G) - A2]
+A [A 2 - ). (A + 2G)]
~ = (A+2C)(4GA+4CP)-A3_2A 2G+A3+ A3_ A3_2A2G
4CA 2 + 4C 2A + 8G 2A + 8G3 - A3 - 2A 2G + A3 + A3 - A3 - 2A 2C
12Q2 A + 8G3 = 4G2 [3A + 2G] (4.26)
92 4. ELASTICIDAD LINEAL

En consecuencia:

1 O" xx >. >.


Exx = t1 O"yy (>' + 2G) >.
0":;:; >. (>' + 2G)
Exx = ~ [O" xx { (>' + 2G)( >. + 2G) - >. 2 } - >. {O" yy (>' + 2G) - >.0"z:; }
+ >. {(O"l/Y>') - (>' + 2G) (O":;:;)}]
Exx = ~ [O"xx (4)'G + 4G 2) - O"yy>.2 - 20"yy>.G + >.20"zz + O"yy>.2
- O" zz>' 2 - 20":;z>.G]
1
Exx = t1 [O"xx (4G) (>' + G) - 2>'G (O"yy + o"z:;)]
(4G) (>' + G) [ 2 > ' G ]
Exx = 4G2 [3>' + 2G] O"xx - (4G)(>' + G) (O"yy + o":;z)

Haciendo:
1 (4G)(>' + G) 2>.G
E 4G2 [3>' + 2G] /J = (4G)(>' + G)

Simplificando:

1 (>' + G) (4.27)
E G[3>' + 2G] /J = 2 (>' + G)
(4.28)

De man(-!ra similar, para Eyy y E:;:;

1
E [O"yy - I/(O"xx + O"zz)] (4.29)
1
E [O":; z - 1/ (O"xx + O"yy)] (4.30)

Las ees 4.28,4.29 Y 4.30 son las leyes generalizadas de Hooke para relacionar
deformaciones unitarias con esfuerzos normales.
4.3. ECUACIONES CONSTo MARCO DE REFERENCIA ... 93

Para establecer la relación entre los esfuerzos cortantes T y las deformaciones


angulares "f, basta recurrir a lo anteriormente establecido.

T xy T xy
"fxy (4.31)
J-L G
T yz T yz
"fyz (4.32)
J-L G
T xz T xz
1xz (4.33)
J-L G

En notación indicial, las ecuaciones 4.28 a 4.33, quedan representadas como:

(4.34)

Ejemplo 4.1. Defina explícitamente a los valores de A y G en términos de E


y 1/.

De las expresiones:

1
C xx = E [CT xx - 1/ (CT yy + CT zz )]
1
Cyy E [CT yy - 1/ (CT xx + CT zz )]
1
Czz E [CT zz - 1/ (CT xx + CT yy )]

Resolviendo el sistema para CT xx , CT yy Y CT zz


94 4. ELASTICIDAD LINEAL

-v -// / 1
~ I }// 1 -v -
- / / -// 1 E
1
~ = [1(1-V 2)+V(-//-1/2)-V(V2 +V)] E3
[1 1
- //2 - // 2 - // 3 - //3 - 1/2] E3
1 1
.6. = - [1 - 31/2 - 21/3] = [(1 + v) (1 + v) (1 - 21/)]3
E3 . E
xx -v
.6. x ééyy 1 ~~ 1 = [exx (1 - v 2) + V (é yy + Vezz)
I
ez= -v
1
-v (-VeYY - e zz )1J __
E2

.6. x = [éxx (1 - v 2) + Véyy + V2ézz + V2éyy + Vézz ] ~2


= [éxx (1 - v 2) + éyy (v + v 2) + ezz (v + v 2)] ~2
.6. x = 2
[éxx (1 - v ) + (v + v 2) (éyy + é zz )] ~2

Eexx (1 - v) vE (éyy + ezz)


Cf xx =
(1 + v)(l - 2v) + (1 + v) (1 - 2v)
Eexx - Ee=v + v Eeyy + v Eezz + Eexxv - Eexxv
(1 + //) (1 - 2v)
Eé xx + vE (é xx + é yy + é zz ) - 2Eé xx V
(1 + v)(l - 2v)
Eé xx (1 - 2v) + vE (é= + éyy + é zz )
(1 + 1/) (1 - 21/)
Eéxx Ev)
Cf xx -:------:- + 1
(1 + //) (1 + 1/)(1 - 2v)
Cf xx )..)1 + 2Géxx

donde:

)..= Ev G= E
(1 + v) (1 - 2v) 2 (1 + v)
4.3. ECUACIONES CONSTo MARCO DE REFERENCIA ... 95

Para () yy y () zz se obtiene:

() yy >"J1 + 2Géyy
()z:: >"J1 + 2Gé zz
Apliquemos las ecuaciones 4.28 a 4.33 para buscar el significado físico de E
y 1/, Fig. 4.2.

Figura 4.2: Interpretación física de las constantes E y V.

En el continuo se establece el tensor:

T¡j = [()r g g]
El tensor de deformaciones unitarias quedará definido por:

Aplicando las relaciones constitutivas:

_v~~· l
o

[T -v-
O
E
()XX

Dado que:
1
é xx = E()xX
96 4. ELASTICIDAD LINEAL

y haciendo Cxx = x ===} y = Ex ; siendo la pendient.e m = E ; CT xx = Y


La constante E representa la pendiente del diagrama esfuerzo-deformación
de una barra prismática sometida a tensión uniaxial, y recibe el nombre de
módulo de Young.
El módulo E sólo puede ser aplicable cuando se defina la líne.-8, recta en el
diagrama esfuerzo-deformación y éste sea independiente del tiempo.
En caso contrario diremos que el material es inelástico y est.o implica que
las relaciones const.itutivas del continuo deben implicar operadores
diferenciales respecto al tiempo.
Analicemos ahora la interpretación física de 1/.

Del tensor Eij se tiene que:


CT
- v xx-
Cyy E = -v Cyy Czz
CT xx
v=--=--
cxx Cxx Cxx
E

1/resulta ser la relación entre las deformaciones unitarias transversales y la


longitudinal, en valor absoluto y recibe el nombre de relación de Poisson.
Si Cxx = x Cyy = y y= -vx
v resulta ser constante sólo cuando la relación entre Cxx y Cyy es una constante.

4.4 Solución de problemas elásticos haciendo


uso de funciones de esfuerzo
Las relaciones constitutivas de los materiales elásticos, lineales homogéneos
e isótropos se pueden expresar como:

(4.35)

donde J( y G son cantidades constantes independientes del tiempo.


Estas relaciones también pueden ser escritas como:

CTxx = >"J] + 2Gcxx T xy = 'YXyG


CT yy = >"J] + 2Gcyy T yz = 'YyzG (4.36)
CT zz = >"J) + 2Gczz T zx = 'Yzx G

De manera recíproca también se puede dar la siguiente forma a las relaciones


constitutivas:
4.5. FUNCiÓN DE ESFUERZOS DE AIRY 97

Cxx = f1 [O"xx - V (O"yy + o"::z)]


Cyy = f [O"yy - 1/ (O"xx + O"zz)]
( 4.37)
10:::: = E [O"z:: - 1/ (O"xx + O"yy)]
T xy T xz T yz
Ixy = G Ixz = G Iyz = G
Las relaciones anteriores, en cualquiera de sus formas, se pueden utilizar para
solucionar problemas elásticos empleando funciones de Airy.

4.5 Función de esfuerzos de Airy


Airy propusó el empleo de una función 1 = 1 (x, y), continua y derivable, que
permite definir a los elementos de un tensor ~j en el que no existan fuerzas
de cuerpo.
Mediante el uso de la Mecánica del Medio Continuo logra definir las
condiciones que debe satisfacer 1, para cumplir los requisitos de continuidad
y equilibrio. La secuencia a seguir es:
Estática ---t Dinámica ---t Cinemática
</J ---t Tij R C Eij

Para un estado de esfuerzo plano los elementos del tensor esfuerzo se pueden
calcular como:

Veamos que sucede desde el punto de vista de equilibrio al aceptar esta


~1 (4.38)

definición:

!x = O
98 4. ELASTICIDAD LINEAL

~F -O
D Z-

Eh xz 8T y Z 8{}zz f - O
8x + 8y + 8z + z-

fz = O
En consecuencia, al aceptar la definición de Airy, el equilibrio de medios
continuos sin fuerzas de campo, es automáticamente satisfecho. Así, si
<P (x, y) es continua y derivable, el equilibrio en el medio se cumple.

Ejemplo 4.2.
a) Sea <p = AX2
Por definición de función de Airy:

Fuerza resultante T
IF 1= 2Abt
v 1h

t
b
Figura 4.3: Barra prismática sometida a un estado de esfuerzo uniaxial.

Por lo tanto la función qJ = AX2 resuelve el problema de una barra sometida


a fuerzas colineales de magnitud 2Abt.
4.5. FUNCIÓN DE ESFUERZOS DE AIRY 99

b) Sea la función <p = By2; esta función represent.a la solución de un continuo


sometido a dos fuerzas horizontales colineales, de magnitud:

c) Sea la función <p = AX2 + By2¡ esta función resuelve el problema de una
barra prismática sometida a un estado de esfuerzos biaxial.

T
h

1 b
t
Figura 4.4: Barra prismática sometida a un estado de esfuerzos biaxial.

d) Sea la función <p = Cxy; esta función resuelve el estado de cortante puro.

e
Figura 4.5: Estado de cortante puro.

Para buscar los desplazamientos generados en el continuo por un conjunto


de fuerzas definidas por <p = <p(x, y) , apliquemos el procedimiento descrito
por la mecánica de medios continuos.

!
Relaciones constitutivas
100 4. ELASTICIDAD LINEAL

[(1 + 1/) 8Di?q,_1/ ( fhI 82 q,

[~~V'~l
2 2
1.
E
fPq, +8Di?q,)] -2C8x8y
Eij=
[
_~
2C8x8y ~ [(1 + v) ~-V\72<p]
O O i ]
(4.39)

Para afirmar la existencia del campo de desplazamientos J, debe comprobarse


la compatibilidad de E i j. La primera ecuación del primer grupo de
com patibilidad es:

Desarrollando:

De aquí puede escribirse:


4.5. FUNCiÓN DE ESFUERZOS DE AIRY 101

Esta última expresión puede escribirse como:

(4.40)

La 1a ecuación del primer grupo de ecuaciones de compatibilidad se satisface


si rp es una función biarmónica.
Verifiquemos ahora la segunda ecuación del primer grupo.

En consecuencia

(4.41)

Esta última ecuación se satisface si /J = ° o bien

~:~ es una función armón'ica cuando /J =1= 0, también se satisface la ecuación.


La tercera ecuación del primer grupo es:

La eco 4.42 es válida si /J = O, o bien verificando que ~:~ es una función


armónica cuando v =1= O.
1a ecuación del 2° grupo:
102 4. ELASTICIDAD LINEAL

2 (¡'12 Ex.",
oyoz
O
2a ecuación del 2° grupo:

202éyy ~ [O'Yyz _ o'Yxz + a'Yx y ]


oxaz ay ax ay oz
O O
:)n ecuación del 2° grupo:

(4.43)

La eco 4.43 es válida si /1 = O o bien si


a2 c/J
es una función armónica.
axay
En consecuencia, en medios elásticos en equilibrio, existe solución cuando es
posible definir a una función de Airy cp = cp (x, y) con los siguientes requisitos:
a) Si 1/ = O ; cp debe ser una función biarmónica.
b) Si 1/ =f. O ; cp debe tener derivadas segundas armónicas, es decir Tij debe
tener a sus elementos armónicos.

Ejemplo 4.3. Sea cp - C y3 una función continua, cuyo grado permite


afirmar que \J4qJ = O Y \J2T;j = O.

/'
/'
h~
b=1

Figura 4.6: Barra prismática sometida a flexión pura.


4.5. FUNCiÓN DE ESFUERZOS DE AIRY 103

Por lo tanto, la función propuesta es solución de algún problema elástico


lineal, el cual se identifica en lo que sigue:
Aplicando la definición de Airy se obtiene.

() xx=- 6Cy; ()yy =O ; Txy= O

Los elementos del tensor ~j resultan:

-6ey O O]
~j = [ 00 O O
O O

Ahora se analizarán cada una de las caras de la barra, para conocer si existe
un estado de esfuerzos.
En x = O; n =-?,

-6ey
O O O
O O]
[ O O O
[Tij llTiI = 6ey?'
In . Ti = -6ey ; Un = 6ey?'
In - Un = 6ey?' - 6ey?' = O

-3Ch
Figura 4.7: Volumen de esfuerzos normales en x = o.

El volumen de esfuerzos que se muestra en la Fig. 4.7, resulta ser un par M


alrededor del eje z.
En la cara x = L; n = ?,
104 4. ELASTICIDAD LINEAL

-6Cy O O]
[T·]
I.J
[ O O O
O O O
[7ij]In;1 = -6Cy'2
an In' n; = (-6Cy) '2. '2 = -6Cy'2
a" = 10',,1 Ti = -6CyI
Tn In - 0'" = -6CyI - (-6CyI)
T" O

El volumen de esfuerzos mostrado en la Fig. 4.7, es equivalente a un par -M


en la cara x = L.
h -
En la cara y = '2; Ti = j

-3Ch O O]
['lij] = [ 00 O O
O O
I" = ['lij]In;1 = O ~ la cara se encuentra descargada.
h -
En la cara y =- -'
2' Ti = -1'. .,

+3Ch O O]
['lij] = [ O O O
O O O
I" = ['lij]ITiI
= O ~ la cara se encuentra descargada
La función de esfuerzos propuesta resuelve el problema de una barra
prismática sometida a flexión para cualquier material elástico.
, M
Dado que e= 6Iz' entonces

Esta última ecuación se conoce como fórmula de la escuadría, la cual es de


suma importancia en mecánica de mat.eriales.

Ejemplo 4.4. Sea la función f/J = 3F xy -


4c
~'\ xy3
4c-
=
4c
y
3F {Xy - x :} , una
3c
función continua y derivable.

La función q) representa la solución de una ménsula sin fuerzas de cuerpo,


sometida a una carga en su extremo libre.
Ya que el tensor no tiene términos armónicos, la solución no es aplicable a
medios continuos con 1/ =1= O.
r
4.6. FUNCIONES DE ESFUERZOS EN COORDENADAS CILíNDRICAS 105

Obviamente V\p = O; es decir, resuelve el problema en medios elásticos en


los cuales 11 = O .

Figura 4.8: Ménsula sometida a una carga puntual en su extremo libre.

Queda como ejercicio al lector demostrar lo antes dicho.

4.6 Funciones de esfuerzos en coordenadas


cilíndricas
El tensor esfuerzo en coordenadas cilíndricas resulta igual a:

En coordenadas cilíndricas las ecuaciones de equilibrio se pueden expresar


como:

O (4.44)

O (4.45)

O (4.46)

Siendo fr> fu y fz, las componentes del vector fuerza de cuerpo en coordenadas
cilíndricas.
En coordenadas polares el tensor esfuerzo se reduce a

T,'J.. -- [O'rr
TrU
TUr]
O'U(J

Por lo tanto las ecuaciones de equilibrio que debe satisfacer la función de


Airy para un estado de esfuerzo plano y fuerzas de cuerpo nulas son:
106 4. ELASTICIDAD LINEAL

OCT'T 1 OTrO 1
- - + - - - + - (CT,.,. - CTOO)
01' l' oe T
o (4.47)
OTrO 10CTo 2
-01'- + - - + -T,·O o (4.48)
l' oe T

La función de Airy en términos de coordenadas polares resulta:

Haciendo el cambio de variables se puede establecer que:

02</J
CT(J(J (4.49)
01'2
10</J 1 02</J
CT rr - - +2- -2 (4.50)
l' 01' 1' 0e
Tr(J = -:1' (~~:) (4.51)

Para comprobar la compatibilidad de deformaciones bastará establecer la


definición de \72 en coordenadas cilíndricas.

Siendo

(4.52)

Con esta definición se puede checar la armonía y la biarmonía en problemas


específicos.

Ejemplo 4.5. Sea la función


p
</J = - - (1'esene)
7f
4.6. FUNCIONES DE ESFUERZOS oo. 107

continua y derivable donde P es una constante.

Cálculo de los elementos del tensor esfuerzo.

(J'rr Ter O1
T,..
I.J [ Tóe (J'dO 8J
DqJ P P
--f)senf) = --f)senf)O
01' 7r 7r
oc/J P P
- -1' senf) - - Tf) cos f)
of) 7r 7r
02c/J P Pr
of)2 = --1' cos f) - - (f)senf) (-1) + cos f))
7r 1f
Pr
- [f)senf) - 2 cos f)l
7r

( ~ Dc/J)
O
= = O
01' Df)
l'

P P 2P
(J'rr = --()senf) +- [()sen() - 2 cos f)l = - - cos ()
1fT 7rT 1fT
(J'ee O ; TrIJ = O

Por lo tanto los elementos del tensor esfuerzo son:

cos () 00 00] _ 2P
TrIJ = [ 00
O O 7rT

En cualquier medio este tensor genera un campo de esfuerzos en el cual


f x = fy = fz = o.

Ya que únicamente existe un sólo término (J'rn calculemos:


108 4. ELASTICIDAD LINEAL

2
(
82 +~ 8 +~ 8 ,2 ) U""
01'2 l' c')¡- 1'2 80
2P
U'T --cosO
'ir1'
8u rr 2P
-2 cosO
01' 'ir1'

8u,." 2P
- senO
80 'ir1'

rPu"" 4P
- - cosO
81'2 7lT:¡

8 2U .-r 2P
- cosO
O(p 'ir1'

- -
4P3 cos O+
2P
__"\ cos O+ -2 P)
3 cos O
( 'irT 'ir1~ 'ir1'

O \l2Tij = O

Por lo tanto el tensor propuesto es armónico; resuelve problemas elásticos


para cualquier medio con 1/ i= O.

La función propuesta representa al estado radial simple.

O O]
Tij = [ UÓr O O
O O O

7tr
Figura 4.9: Estado radial simple.
4.6. FUNCIONES DE ESFUERZOS ... 109

2P
(J'rr = - - cos()
7fr

La fuerza que actúa sobre un elemento dA = rd() es (J'rrrd() Y su componente


vertical es ((J'rrrd()) cos(). La resultante de tales componentes queda:

1
2 ~ (J' n ° cos ()rd() = - -4P¡~ cos 2 ()d()
O 7f O

_ 4P [~+ Sen2()] ~ = _P
7f 2 4 O

La función propuesta <p = -~ (r()sen()) genera el estado radial simple que


aplicado a un semiespacio elástico conduce a la solución de Flamant.
Variación de (J' rr con () oonstante y r variable.

2P
Haciendo - - = e r=x
7f
e
====> y = - ; la variación de esta última función se muestra en la Fig.4.lO.
x

IP
!
r------~
y
= cte = O

variación hiperbólica

X
Figura 4.10: Variación de (J'rr para () = cte y r variable.

La Fig.4.lO muestra que todos los puntos del semiespacio excepto los de la
frontera superior, contribuyen a soportar la carga P.
Ahora definiremos el lugar geométrico de los puntos del semiespacio con el
mismo valor del esfuerzo principal.
110 4. ELASTICIDAD LINEAL

le
~ r

1
1
2P
O"rr = - 7tO
Figura 4.11: Lugar geométrico del semiespacio con el mismo valor de U rr •

2P
(Trr =--
7rr

En relación con la Fig.4.11 el punto A que se encuentra sobre la circunferencia


tiene como radio vector:

r = DcosO

Por lo tanto el esfuerzo radial vale:

2P 2P cos O 2P
U'T = - - cosO = ---- = --
7rr 7r D cos () 7r D

En consecuencIa, ya que A se encuentra sobre la circunferencia, se debe


tener el mismo esfuerzo radial u rr - - ;~ en todos los puntos de la
circunferencia, exceptuando al punto de tangencia de la circunferencia con
la frontera superior.
Se llama isobara a la circunferencia de diámetro D, correspondiente al
2P
es fuerzo (T'T = - "D'
En la Fig. 4.12 se muestra la construcción de isobaras para diferentes valores
de (TrI" .
4.6. FUNCIONES DE ESFUERZOS ... 111

m,,,,, 1 """i~'~'0:¡¡¡\\"
/:f;,;¡J:2e
''líl r ¡ \\~""
iT TI
A k+Lt-----~J

i
\\
!

I (. \ \
I ¡ \

\ \. ." --0- . / /
\ \.... . . "'-"
\ \ \ \ 1/4D1/2D
I 1\ \

/ )1
4crrr ./ I
// /
. . . . . ._ --0._._--.... ///
\

/
\
;
11 1 3 4D
/
D

. " 2crrr _ / ; I
\ ---e--'- /

',."''''~:~~/ ~
crrr
Figura 4.12: Isobaras para el estado radial simple.

Ejemplo 4.6. Una roseta de deformación formada por tres "strain gages"
(SG), con el arreglo que se muestra en la FigA.13 se coloca en un punto de un
medio continuo, obteniéndose las mediciones que se reportan a continuación.

a) Determínense los elementos del tensor deformación é xx , é yy Y'Yxy'

Si los parámetros elásticos del medio continuo son: E = 2.1 X 105 M Pa y


/J = 0.2.

b) Calcúlense los elementos del tensor esfuerzo (}xx, (}yy Y Txy .

SG3

e1
"P" ~-E't===3-------::':'
SG1
Figura 4.13: Roseta de deformación en un punto de un cuerpo deformable.

La deformación longitudinal asociada a una dirección cualesquiera, está dada


por:
112 4. ELASTICIDAD LINEAL

C;cx + Cyy Exx - c1JY 1 ()


Ee = 2 + 2" cos 2() + 2')'Xysen 2
Aplicando est.a ecuación a las diferentes direcciones en que fueron colocadas
los strain gages, se tiene:

Para EC I = 3 X 10- 4

3 X 10 - 4 -
-
Cxx + Cyy ,1
,
E 3;x - Cyy
cos (00) + ~ "'.,ysen
. (00)
2 2 2 ,-
4
==? Exx = 3 X 10-
Para ce2 = 1.5 X 10- 4

1.5 X 10- 4 = Exx +


2
Cyy + Cxx - Eyy
2 cos
(900) + 1
2')'xysen
(900)

1.5 10- 4 =
Cxx + Cyy 1
X
2 + 2')'xy
o ==? Cyy = -')'xy

Para Ee3 =5X 10- 4

10- 4 C= + Cyy Cxx - Cyy 1


5 X 2 + . 2 cos 20 + 2')'xysen2()

5 X 10- 4 Cxx + Cyy + Cxx - Cyy cos (180°) + ~')'xysen (1800)


2 2 2
Cxx + Eyy Cxx - Cyy
5 X 10- 4
2 2
==? Cyy = 5 X 10- 4 ')' xy = -5 X 10- 4

Por lo tant.o el tensor deformación vale:


3 X 10- 4 -2.5 X 10- 4 ]
Eii = [ -2.5 X 10- 4 5 X 10- 4

Para calcular los elementos del tensor esfuerzo, haremos uso de las ecuaciones
constitutivas de los materiales elásticos lineales, homogéneos e is6tropos, esto
es,

{j xx 2GCX:l: + )..)1
{jyy 2GE yy + )..)1
T xy G')'XY
4.7. ENERGíA DE DEFORMACiÓN ELÁSTICA PARA UN ESTADO UNIAXIAL DE
ESFUERZOS 113

Cálculo de las constantes elásticas,

E 2.1 X 106 5
G = 0.875 x 10 M Pa
2 (1 + v ) 2 (1 + 0.2 )
6
l/E 0.2 x 2.1 x 10 = 0.583 X 105 M Pa
(1 + 1/) (1 - 21/) (1 + 0.2) (1 - 2 x 0.2)

El primer invariante del tensor deformación resulta,


JI = é xx + é yy = 8 X 10- 4
Por lo tanto los esfuerzos normales valen

CT xx 2 x 0.875 X 105 x 3 X 10- 4 + 0.583 X 105 x 8 X 10- 4


CT xx 52.50 + 46.64 = 99 .14 M Pa
CT yy 2 x 0.875 X 105 x 5 X 10- 4 + 0.583 X 105 x 8 X 10- 4
CT yy 87.5 + 46.64 = 134.14 M Pa
Txy = 0.875 X 105 x (-5 X 10- 4 ) = -43.75 MPa

Finalmente, el tensor esfuerzo vale:

4.7 Energía de deformación elástica para un


estado uniaxial de esfuerzos

8F+
+
F
+ L-~UL--~--------~~ 8i

+ 8 ++
d8
Figura 4.14: Definición de energía de deformación elástica.
114 4. ELASTICIDAD LINEAL

Supongamos que la curva esfuerzo-deformación para un material elástico


linf'..al, ·somet.ido a Ull estado de esfuerzos uniaxial es como se muestra en
la Fig. 4.14.
La energía que un cuerpo absorbe como resultado de su deformación bajo
ciert.a carga se llama energía de deformación.
Est.a se puede expresar como:

(4.53)

La energía de deformación por unidad de volumen o densidad de energia U,


se puede expresar como:

U= ~
V
t Fdo
Jo
(4.54)

Para una partícula elemental de un medio continuo de volumen dV = dxdydz


y asumiendo que la fuerza se aplica en la dirección x, se tiene:

Integrando se obtiene:

2
U = (J xxt: xx = (J xx
(4.55)
2 2E

El área bajo la parte lineal de la curva uniaxial (J - t: es una medida de


la capacidad del material para almacenar energía elástica. Esta medida se
llama módulo de Resiliencia (R) Y se puede calcular como:

(4.56)
4.8. ENERGíA DE DEFORMACiÓN ELÁSTICA PARA UN ESTADO TRIAXIAL DE
ESFUERZOS 115

cr

Figura 4.15: Definición de módulo de Resiliencia.

siendo eLE Y ('"LE la deformación longitudinal y el esfuerzo normal en el límite


elástico, respectivamente.

4.8 Energía de deformación elástica para un


estado triaxial de esfuerzos

dz

dx

dy
x""
Figura 4.16: Energía de deformación elástica para un estado general de
esfuerzos.
116 4. ELASTICIDAD LINEAL

Para un est.ado t.riaxial de esfuerzos el trahajo total realizado por los esfuerzos
0"1, 0"2 Y 0"3 será la suma de los trabajos efectuados por cada uno de ellos de
manera independiente. Fig. 4.16.

Por lo tant.o el trahajo realizado por al se calcula como:

1 1
dW - (aldydz) . (E¡dx) = -O"IEldV (4.57)
~ '-v-' 2
Fuerza media distancia

De esta forma la densidad de energía resulta:


dW 1
- = U = -O"lEl (4.58)
dV 2
Repitiendo el razonamiento para las demás caras de la partícula elemental,
se concluye que la densidad de energía de deformación elástica U almacenada
en el material, debido a un estado de esfuerzos al, 0"2 y 0"3, es

(4.59)

tomando en cuenta que

El

Sustituyendo estas últimas ecuaciones en la eco 4.59, se obtiene:

U
~ [0"1 (al - v (a2 + (3)) + a2(a2 - v (a1 + (3))
2E +0"3 (a3 - v (0"1 + 0"2))]
4.8. ENERGIA DE DEF. ELASTICA TRIAXIAL... 117

Por lo tanto, la densidad de energía se puede descomponer, como:

u= Uv + Uo (4.61)

Siendo Uv la densidad de energía volumétrica y Uo la densidad de energía


desviadora o distorsiona!.

A.<¡umiendo que existe la siguiente relación:


Uv ~ T v parte isotrópica (dilatación o compresión)
Uo ~ T o parte distorsional

La energía de deformación volumétrica se puede calcular como:

O"V 0"1 + 0"2 + 0"3


UV = 2 (el +e2 +e3) ; o"v =
3
E G=
o"v = Kev , K=
3 (1 - 211) . 2 (1 + 11)
E

O"v 2
O"v 12
1
Uv ev = 2K = 18K
2
Uv = _1_
1
2
Ir3 (1 - 211)
18K 18 E
(1 - 211) 1 2
Uv 6E 1
(4.62)

La energía de deformación desviadora se puede calcular como:

Uo U-Uv

Uo ~ (Ir _12) _ (1 - 211) 12


2 E G 6E 1

Desarrollando esta última ecuación se obtiene:

(4.63)
Las densidades de energía volumétrica y desviadora están relacionadas con los
esfuerzos normal y cortante octaédricos mediante las siguientes expresiones,
respectivamente:

2
O" oct 2KUv (4.64)
2 4
7 0ct -GUa (4.65)
3
118 4. ELASTICIDAD LINEAL
Capítulo 5
Teorías de Falla y Ruptura

Aceptemos que en una prueba de tensión el material deja de ser elástico


cuando se alcanza el esfuerzo límite SIen la probeta.

Si se aplican esfuerzos superiores a SI , el material puede fluir o se rompe,


por lo que SI representaría ellírnite de aplicabilidad de la teoría elástica para
el estado de esfuerzos unidimensionales que se presenta en la probeta. Este
límite puede ser establecido mediante pruebas de laboratorio.

Para un medio continuo sometido a un estado de esfuerzos principales, la


función que define la región donde el material tiene un comportamiento
elástico, se puede expresar como

(5.1)

En los párrafos siguientes se describirán algunas teorías de falla y ruptura


comúnmente utilizadas para establecer bajo qué condiciones se alcanza la
falla o ruptura de un material sometido a ciertas solicitaciones. Cabe señalar
que la mayoría de estas teorías tiene su aplicación principal en el
comportamiento de metales.

5.1 Teoría de Rankine

Esta teoría establece que en un material sometido a un estado de esfuerzos


principales se genera fluencia cuando cualesquiera de los esfuerzos principales
alcanza el valor límite SI , en una probeta representativa del mismo material
sometida a una prueba de tensión. Matemáticamente estas condiciones
quedan representadas por:

119
120 5. TEORíAS DE FALLA Y RUPTURA

181 1 (5.2)
18/1
(5.3)
ISII (5.4)

De manera gráfica, las condiciones anteriores se pueden representar en un


sist.ema de referencia donde los ejes representen a los esfuerzos O"}, 0"2 Y O"~.

De esta manera el límite de aplicabilidad de la teoría elástica queda definido


por seis superficies planas que conforman a un cubo de lado 2S1 (Fig. 5.1).

T
2Sf

1 (5\ / / /

q,S'
1--2Sf~
Figura 5.1: Volumen de fiuencia. Teoría de Rankine.

Si un punto P (0"1,0"2,0"3), que representa a un estado tridimensional de


esfuerzos, se ubica en el sistema de referencia establecido, se tendrá una
condición de fiuencia si P se encuentra fuera del cubo o en las caras.

Si P se localiza dentro del cubo, el material tiene un comportamiento elástico.

Para un estado de esfuerzo plano, con 0"3 = O, se obtendrá la curva límite


que se muestra en la Fig. 5.2, delimitada por las siguientes ecuaciones.
5.2. TEORíA DE COULOMB-TRESCA 121

f-<-Sf-+-Sf~
Figura 5.2: Superficie de fluencia para el estado de esfuerzo plano. Teoría
de Rankine.

(71 = ±S¡; (5.5)

Un punto en este plano representaría a un estado de esfuerzo plano. Si el


punto está dentro del cuadrado, el material tiene un comportamiento elástico,
mientras que si está fuera o en los bordes del cuadrado, el material habrá
dejado de ser elástico.

La teoría de Rankine es aplicable sobre todo a materiales frágiles.

5.2 Teoría de Coulomb-Tresca

El material deja de ser elástico cuando el esfuerzo cortante máximo generado


por un estado de esfuerws principales, en un punto cualesquiera del cuerpo,
iguala al esfuerzo cortante máximo que se engendra en una probeta sometida
a tensión.

Para un estado tridimensional de esfuerws el esfuerzo cortante máximo queda


definido por:

lil (5.6)

lil (5.7)

lil (5.8)
122 5. TEORíAS DE FALLA Y RUPTURA

Si cualquiera de las condiciones anteriores es satisfecha, se tendría la condición


límite para que el medio deje de ser elástico.

Para un estado de esfuerzos plano; <T3 = 0, las condiciones matemáticas que


indican iniciación de fluencia serían:

De la eco (5.6)
I<TI - <T21 = ISII ~ (5.9)
De la eco (5.7)
I<TII ISII ~
= <TI = +SI; <TI = -SI (5.10)
De la eco (5.8)
1<T21 = ISII ~ (5.11)

Representando estas condiciones gráficamente en el plano <TI> <T2,


se obtendrían seis líneas rectas límites, que definen la superficie de fluencia
(Fig. 5.3). .

Sr

Figura 5.3: SUI~erficie de fluencia para el estado de esfuerzos planó. Teoría


de Coulomb-Tresca.

De acuerdo con la teoría del esfuerzo cortante máximo, si se agregan esfuerzos


hidrostáticos de tensión o de compresión, no es posible predecir ningún
cambio en la respuesta del material. La suma de estos esfuerzos simplemente
desplaza el círculo de Mohr a lo largo del eje <T y T máx permanece constante.

Cuando los esfuerzos principales <TI y <T2 son del mismo signo, se tienen dos
condiciones para alcanzar la fluencia del material:
5.2. TEORíA DE COULOMB-TRESCA 123

Figura 5.4: Estado de esfuerzo plano; tensión.

(j2 Compresión

(j,

Figura 5.5: Estado de esfuerzo plano; compresión.

Para el estado de esfuerzos de compresión que se muestra en la Fig. 5.4, se


debe cumplir que:
(5.12)

Para el estado de esfuerzos de tensión, Fig. 5.5, se tiene:


(5.13)

Por lo tanto podemos concluir que cuando los esfuerzos principales 0"1 y 0"2
son del mismo signo, las teorías de Rankine y Coulomb-Tresca coinciden.

Si 0"1 Y 0"2 son de signo contrario, se tiene que:

(5.14)

Por lo tanto, cuando 0"1 y 0"2 son de signo contrario las teorías de Rankine y
Coulomb-Tresca difieren.
124 5. TEORíAS DE FALLA Y RUPTURA

5.3 Teoría de Saint Venant

8llnatcrial deja. de ::;er elástico cua.ndo una. de las deformaciones principales,


él, é2, {) é:¡
alcanza. el valor de la deformación principal, él, que se genera. en
Ulla probet.a. sometida a tensión.

Para. un estado uniaxial de esfuerzos, dicha deformación se puede expresar


como:

S, (5.15)
é/= -
E

Para un estado de esfuerzos tridimensional, las deformaciones principales


quedan definidas por:

léll I~ (UI - V (U2 + (3))1 = lé,1 (5.16)

lé21 = I~ (U2 - V(Ul +(3))1 = lé,1 (5.17)

IE 31 - I~ (U3 - V(Ul +(2))1 = lé,1 (5.18)

Para un estado de esfuerzo plano, con U3 = O , las condiciones de fluencia


resultarían ser:

De la eco (5.16)
IU l - l/u21 = ISII ==> Ul - VU2 = S,
UI - VU2 = -S, (5.19)
De la eco (5.17)
IU 2- vUll = IS,I ==> U2 - I/Ul = S,
U2 - VUI = -SI (5.20)
De la eco (5.18)
I-v (Ul + (2)1 = IS,I ==> v (Ul + (2) = S,
1/ (Ul + (2) = -SI (5.21)

Reordenando términos, las ecuaciones de las seis rectas límites resultan ser:
5.4. TEORíA DE NADAI 125

(TI _ S'I (TI SI


// //
(T2 = -
//
+-
//
(5.22)
+ S'I
//(Tl (T2 = //(Tl - SI (5.23)
S'I SI
-(TI +- (T2 = -(TI - - (5.24)
// //

La región que definen estas ecuaciones se muestra en la Fig. 5.6, en la que


se puede observar que la teoría de Saint Venant permite alcanzar niveles de
esfuerzos mayores que los que definen la teorías de Rankine y Tresca, antes
de alcanzar la falla.

T
SI Iv
--------t

SI Iv

I
.J...

1-- SI I V -1--- SI I v --l


Figura 5.6: Superficie de fluencia, criterio de Saint Venant.

5.4 Teoría de Nadai


Esta teoría es aplicable principalmente a materiales dúctiles, y establece que
se inicia la fluencia en una partícula de un medio continuo cuando se aplica
a ésta una energía de deformación igual a la energía de deformación que se
genera en una partícula de una probeta sometida a tensión.
126 5. TEORíAS DE FALLA Y RUPTURA

La f!lIergía de deformaciól! elá..,tica por unidad de volumen, o densidad de


energía, para l1U estado unixial de esfuerzos, tal como se definió en el capítulo
4, resl1lt.a igual a:

u= O"lél
2

Esta ecuación representa el área bajo la curva esfuerzo-deformación.

De acuerdo con esta teoría la probeta dejará de ser elástica cuando:

u _ SI él _ SI SI _ S} (5.25)
prob - 2 - 2E - 2E

Para un estado tridimensional de esfuerzos, la densidad de energía puede


expresarse como:

(5.26)

Por lo tanto, la fluencia se presentará cuando:

(5.27)

Por lo tanto la condición de fluencia queda representada como:

(5.28)

Para un estado de esfuerzo plano, con 0"3 = 0, la condición de fluencia queda


definida por:

(5.29)

La representación geométrica de la región que define la eco 5.29 se muestra


en la Fig. 5.7.
5.5. TEORíA DE VON MISES HENCKY(VMH) 127

Parav = 0.5

1
f--S

Figura 5.7: Superficie de fluenciaj teoría de Nadai.

5.5 Teoría de Von Mises Hencky(VMH)

Se alcanza la fluencia en una partícula de un medio continuo cuando la energía


de deformación distorsional en un estado de esfuerzos cualesquiera, igual a
la energía de deformación distorsional en una probeta sometida a tensión.

Para valuar Uo empleamos la siguiente relación:


Utotal = Uvol + Uo
De donde:
Uo = Utotal - Uvol

La energía de deformación total se puede expresar como, eco (4.60).


128 5. TEORíAS DE FALLA Y RUPTURA

I~l csfll~rzo volumétrico y la deformaci61l volumétrica son igl.t ales a:

De (-::;t.a man~ra la energía de deformación volumétrica resulta igual a:

Uvol = CT v
;E )

Dado que:
E
(Tv = ¡<Ev = 3 (1 - 2// )Ev

Despejando Ev.

Sustituyendo esta última expresión en Uvol , se tiene:

u _ 3 (1 - 2//) 2
vol - 2E O"v

Desarrollando:

3(1 - 2//) (0"1 + 0"2 +0"3)2


2E 9
1 - 2// 2
6E ((TI + 0"2 + 0"3)

Dado que:

La energía de deformación volumétrica queda como:

(5.30)
5.5. TEORíA DE VON MISES HENCKY(VMH) 129

La energía de deformación distorsional resulta:

Desarrollando se tiene:

(5.31)

En una probeta sometida a tensión: (j2 = (j3 = O

(1 + 1/) 2 (5.32)
Uprob = 3E (jI

Se alcanza la falla del material cuando: Uprob = Uo, por lo tanto:

De donde resulta
(5.33)

Condición de fluencia de Hencky, que coincide con Nadai si 1/ = ~.

Para el estado de esfuerzos plano.


2
(jI + (j22 - (jI(j2 = 2
SI (5.34)

Esta teoría resulta ser un caso particular de la de Nadai.

Experimentalmente se ha demostrado que la teoría de VMH, es la teoría que


más se apega a los valores experimentales si (jI y (j2 son positivos o si alguno
de los dos es positivo.

Si ambos esfuerzos principales son negativos, VMH, dá resultados


conservadores.

Esta teoría es aplicable sobre todo a metales.


130 5. TEORíAS DE FALLA Y RUPTURA

La teoría de VMH pre;enta la sig·uiente ventaja de tipo matemático.

Para un est.ado de esfuerzos principales se tiene que:

(0"1 + 0"2 + 0"3 ) (5.35)


(0"10"2 + 0"10"3 + 0"20"3) (5.36)

Elevando la eco 5.:35 al cuadrado

(5.37)

Multiplicando la eco 5.36 por 3.

-312 = -30"10"2 - 30"10"3 - 30"20"3 (5.38)

Sumando las ecs. 5.37 y 5.38 se tiene:

/ 12 - 312 = 2
0"1 + 0"22 + 0"32 - 0"10"2 - 0"10"3 - 0"20"3

lo cual indica que existe fiuencia en un material si

(5.39)

Esta últ.ima expresión puede ser aplicada para cualquier sist.ema de referencia.

Ejemplo 5.1. El tensor esfuerzo en el medio continuo que se muestra en la


Fig. 5.8, está dado por:
5.5. TEORíA DE VON MISES HENCKY(VMH) 131

La teoría del medio continuo será aplicable si


puntos del medio.
Ir - 31 2 < s/, en todos los

Si en algún punt.o Ir -312 = s /, se establecería el límite del tensor 'T.¡j hasta


donde sería aplicable la mecánica del medio continuo.

Si en algunas regiones Ir - 31
continuo no será aplicable.
2 > (T / , la teoría de la mecánica del medio

Si P se aplica al medio, se busca de definir la región del medio continuo en


el cual sea aplicable la teoría elástica.

Llamando:

Entonces las invariantes valen:

/1 2xyk
12 _(C2 _y2)2k2

Sustituyendo los valores de / 1 e 12 en la eco (5.39) existirá fluencia cuando:


4x2y 2k 2 + 3k 2 (c4 _ 2c2y2 + y4) = sJ

Desarrollando:

Esta última ecuación representa la condición de fluencia de VMH.


132 5. TEORíAS DE FALLA Y RUPTURA

La elal>t.icidad sería aplicable mient.ras 110 se plastifique algún punto. Los


punt.os más esforzados son a. ( L, -e) y b (L, +c).

Se iniciará la Huencia en el medie en el instante en que las coordenadas de


los punt.os él y b satisfagan la condición de VMH.

Sust.ituyendo las coordenadas de los puntos en la condición de Huencia, se


tiene:

Si la fuerza P del extremo, está comprendida entre los límites.


2
_ 2¿'af < P < 2c af
3L 3L

El medio es elástico.

2
Si IPI > 2c af el material deja de ser elástico.
- 3L
Casos límit.es de aplicación de la mecánica del medio continuo.

a) Pieza prismática sometida a flexión simple.

Me f - --+--C_ _-+I-t\-=--::.:M=---------.. y
L--------'--c_ _ __ --'I j
Figura 5.9: Viga sometida a flexión pura.
5.5. TEORíA DE VON MISES HENCKY(VMH) 133

El tensor esfuerzo está dado por:

existirá pla:stificación si Ir - 3/ 2 > sJ

Las rectas límites pasan por el borde del medio si c = I z


M,
s,
M, = s,-
Iz
c

Si

IM,I < s, I
e
z
todo el medio es elástico.

IM,I > Is,


e
z
; el medio deja de ser elástico y sólo una porción próxima al

eje y permanecerá elástica.

Esta distribución de esfuerzos aparecerá cuando la sección se ha plastificado


completamente, y será engendrada por un momento plástico total M p , tal
que:

Por estática se puede afirmar que el Mp es la resultante de la distribución de


esfuerzos en la sección transversal completamente plastificada.
134 5. TEORíAS DE FALLA Y RUPTURA

Por otra parle

M¡=~b(2e)38¡
12 e
2 2
M¡ = -be 8¡
3

Si el material que forma la pieza tiene el diagTama supuesto, se puede


incrementar el moment.o que provoca la la Huencia en un 50 % para alcanzar
el momento que provoca ftuencia en toda la sección transversal.

La relación que existe entre el momento de la Huencia y el momento de


plastificación total, Mp/ M¡, depende de la forma de la sección transversal y
puede oscilar entre 1 y 2.5.

1.07 < ~~ < 1.17 ~~ =1.5 ~~ > 1.5

Figura 5.10: Relación Mp/M¡ para diferentes secciones.

Ejemplo 5.2. Para el estado de esfuerzo:

~j = UTO]
[ TO UO uO siendo u y T constantes

Determine el esfuerzo de Huencia empleando los criterios de Coulomb-Tresca


y Von Mises Hencky.
5.5. TEORÍA DE VON MISES HENCKY(VMH) 135

Los esfuerzos principales se obtienen resolviendo la ecuación característica,


resultando:

Así el tensor esfuerzo queda como:

Condición de fluencia de VMH.

Los invariantes del tensor esfuerzo son:

/1 30"

/2 = 1(O"t 7 ) ~1+lo (O" -O 1+1 (0"0 7 ) 7)


O
(0"+7)
12 (O" + O" +O" (O" -
7) +(O" - 7)(0" +
7) 7)

12 0"2 +70" +0"2 - 70" +0"2 +70" - 70" - 7


2

12 30"2 - 72

Por lo que aplicando la condición de fluencia de VMH, se tiene:

(30")2 - 3 (3(72 -
2
7 ) = s}
9(72 - 9(72 + 37 2
= s}
2
SI = 37 2 ==} SI = ±V37
SI = ±1.737
136 5. TEORíAS DE FALLA Y RUPTURA

ApliC<LlJ(lo el crit.(~rio de Coulornh-Tresca.:

(al - (3) 8f
Tmáx= 2 2
(a+T)-(a-T) af
2 2
= } af = 2T

Obsérvese que este último criterio permite un mayor esfuerzo de ftuencia


comparado con la teoría de VMH.

Ejemplo 5.3. Una muestra cilíndrica de un material deformable está confinado


por un molde rígido que no le permite deformarse lateralmente, bajo una
presión constante "p". Aplicando el criterio de VMH diga si el material
alcanza la condición de ftuencia.

Para establecer el estado de esfuerzos y deformaciones suponga que el cuerpo


deformable es elástico lineal, homogéneo e isótropo. Suponga además que no
se produce esfuerzos cortantes en el contacto molde-muestra.

P molde
a) ~rígidO

~-Y

I tx :
b) I I
I

I
t
Figura 5.11: Muestra cilíndrica en un zmolde rígido, a) alzado, b) planta.

Datos:
5.5. TEORíA DE VON MISES HENCKY(VMH) 137

E 2.1 X 108 kPa


V 0.2
S, 4000 X 102 kPa
p 1000 X 102 kPa

Solución:

Las ecuaciones constitut.ivas de los materiales elásticos lineales, homogéneos


e isótropos son:

{7xx 2Géxx + >"J1


{7yy 2Géyy + >"J1
{7zz 2Gézz + >"J1

Cálculo de las constantes elásticas.

E 2.1 X 108 8
G 2 (1 + 0.2) = 0.875 x 10 kPa
2 (1 + v)
8
vE 0.2 x 2.1 x 10 = 0.583 X 108 kPa
(1 + 1/) (1 - 2v) (1 + 0.2) (1 - 2 x 0.2)

De los datos del problema se pueden establecer las siguientes condiciones.

{7xx = -p é xx =1- O
{7 yy =1- O é yy = O
{7 zz =1- O é zz =O

Para {7 yy se tiene:

Para {7 zz se tiene:
138 5. TEORíAS DE FALLA Y RUPTURA

Para u xx ~e tiene:

-1000 = 2Gé xx + >. (é xx + é yy + é zz )

-1000 X 102 2Gé xx + >.é xx = é xx (2G + >')


-1000 X 102 -1000 X 10 2
(2G + >') (2 x 0.875 x 108 + 0.583 X 108 )
-0.428 X 10-3 = JI
0.583 X 108 x -0.428 X 10- 3 = -249.5 X 102 kPa
u zz = -249.5 X 102 kPa

Por lo tanto el tensor esfuerzo resulta igual a:

'T¡j = [
-1000
O
O
-249.5
O]
O x 102 kPa
O O -249.5

Aplicando el criterio de VMH, se tiene:

Sj = 1i - 312

Cálculo del primer invariante 11 •

11 + U yy + U zz = -1000 x 102 -249.5 x 102 -249.5 x 102


u xx

-1499.0 X 102 kPa


12 = {(-1000 x -249.5)+(-249.5 x -249.5)+( -249.5 x -1000)} x 104
12 = 561250 X 104
sj 1499.02 -3 x 561250 x 104 =563251 x 10 4
SI = 750.5 x 104 kPa (esfuerzo de fluencia calculado)

Dado que 750.5 x 104 kPa < 4000 x 104 kPa no se presenta fluencia en el
material.
Apéndice A

Análisis Tensorial

El análisis tensorial se centra en el estudio de entes abstractos llamados


tensores, cuyas propiedades son independientes de los sistemas de referencia
empleados para determinarlos.

Un tensor está representado en un sistema de referencia particular mediante


un conjunto de funciones llamadas componentes.

El que un conjunto de funciones represente a un tensor depende de la ley de


transformación de estas funciones de un sistema coordenado a otro.

Los tensores se clasifican por su orden y el número de componentes en un


sistema de referencia cartesiano es 3N , siendo N el orden del tensor. En la
tabla A.1 se reportan algunos tensores, su representación en notación indicial
yel número de componentes.

Tabla A.1: Transformación de coordenadas.


Urden N Nombre del tensor ~lmº-olo No. de componentes
~
escalar a 1
vector bi 3
2 diádica Tij 9
3 triádica Eijk 27

139
140 A. ANÁLISIS TENSORIAL

X~ p
'\
\
~3
.-:.-~-----=:-~ X2
e2

!/
X;
Figura A.l: Coor~enadas de un punto "P" en un sistema de referencia
carteSIano.

La ubicación del punto "P" en el sistema de referencia XI, X2, X3, está dada
por el vector OP = x¡e¡ + X2e2 + X3e3, siendo el> e2 ye3 los vectores base
asociados al sistema de referencia X¡, X2, X3. La ubicación de dicho punto en
el sistema de referencia x'¡, x;, x~, es OP = x;e; + x;e; + ~e~ (Fig. A.l).

Un conjunto de ecuaciones

X: = Ji (Xl, X2, X3)


describe una transformación de las variables XI, X2, X3 a las variables x;, x;,
x~.
De la misma forma se puede plantear la transformación inversa como:
' , ,3)
Xi =gi (X¡'X2'X
Con el objeto de asegurar que la transformación inversa existe y es uno a
uno, en una región R, se debe cumplir que:

a) Las funciones Ji tienen valor único en el punto "P", son continuas así
como sus primeras derivadas.

b) El determinante del Jacobiano no se anula en ningún punto de la región


R, esto es,
OX; o'
x¡ ox;
oXl oX'f OX~
oX 2 oX2 oX 2
J =lox;
OXj
1= OX} fJx'f OX~
fO
oX3 oX3 oX3
oX¡ OX2 OX3
141

Si se cumple lo establecido en a) y b), se dice que la transformación es


admisible. Si el Jacobiano es positivo en cualquier punto de la región R,
se dice que la transformación es propia, en caso contrario se dice que la
transformación es impropia.

En el estudio de los medios continuos que se abordará en estas notas se acepta


que las transformaciones son admisibles y propias.

Consideremos un sistema de referencia cartesiano Xl, X2, X3 Y los vectores


base asociados a cada eje el, e2 y e3, tal corno se muestra en la figura A.2.

X3

X1
Figura A. 2: Sistema de referencia cartesiano.

Debido a que los vectores base son mutuamente ortogonales y tienen un


módulo unitario, se dice que forman una base ortonormal.

Cualquier vector v, puede ser expresado en esta base corno:


3

V = vlel + V2e 2 + V3 e 3 = L Viei (A.l)


i=l

siendo VI, V2 Y V3 los componentes del vector V.

Notación suma
Siempre que en una ecuación una letra aparezca corno subíndice o superíndice
repetido, se entenderá que se debe sumar sobre esta letra. Así si se está en
un espacio de n dimensiones, se sumará de 1 a n. Los índices repetidos son
mudos (o falsos), ya que el resultado final es independiente de la letra usada.
De esta manera la eco A.l se puede expresar corno:

(A.2)
142 A. ANÁLISIS TENSORIAL

Ejemplo (a.l): sin importar por el momento el significado de las siguientes


ecuaciones, representarlas en forma expandida, con variación de los subíndices
de 1 a 3.

Solución:

a) Desarrollando primero el subíndice i y luego el j

UiViWiCj = (1LlVl + U2V2 + U3V3) (wlel + W2e2 + W3 e 3)

b)

T1jVlej + T2jV2ej + T 3j V 3e j
Tllvlel + T 12 v le2 + Tl3vle3 + T2l v 2e l
+T22 v 2e 2 + T 23 v 2e 3 + r'U v 3 e l + T 32 v 3 e 2
+T33 v 3e 3
(TllVl + T2lV2 + T3lV3) el + (T12Vl + T22V2 + T 32 V3) e2
+ (Tl3 V l + T 23v2 + T33V3) e3

c)

(Tll + T 22 + T 33 ) vjej
(Tll + T22 + T 33 ) (Vlel + V2e 2 + V3 e 3)
Definiciones:
1) Suma de vectores en notación indicial.

w = TI + V w¡ei = (Ui + Vi) ei (A.3)


2) Multiplicación de un vector por un escalar.
(A.4)
Delta de Kronecker
Este símbolo se define como:

Ó¡j = { 6 82
82
"¿= J
i=l=j (A.5)
De acuerdo con la convención suma:
143

Este símbolo nos permitirá expresar el producto punto de dos vectores como:
Sean u y V, dos vectores cualesquiera cuyas componentes son:

u u¡e¡ + U2e2 + U3e3 = Uiei


V V1 e1 + V2e2 + V3e3 = Vjej

Por definición de producto punto, tenemos:

u· v = IUIIVI cosO (A.6)


Siendo () el ángulo que forman los dos vectores con origen común.
Así si () = 90°, u· v = O, lo cual significa que los vectores son ortogonales.
En notación indicialla eco A.6 se puede representar como:

u· v = (Uiei)· (vjej) = UiVjei . ej


= UI VI + U2V2 + U3 V 3
Por la definición del símbolo delta de Kronecker.
ei . ej =
Óij , por lo tanto

(A.7)
Símbolo de permutación

Este símbolo se define como:


1 permutación par: 1, 2, 3; 2, 3, 1; 3, 1, 2
éijl.: = { -1 permutación impar:, 1~. 3, 2; 3, ~~ 1; 1,3,2
O aparece al menos un md.1ce repetIdo
El símbolo de permutación tiene las siguientes propiedades

Definición:
El producto cruz de dos vectores u X v, es igual a:
u x V = lullvl senO (A.8)

uxv e1 U2e2 e3 ]
[ UI
V1 V2
U3
V3
e1 (U2 V 3 - U3V2) - e2 (UIV3 - U3VI) + e3 (UI V 2 - U2 V I)

(U2 V 3 -U3 V 2)eI + (U3VI -U1VJ)e2 + (UI V2 - U2 V I)e3

En notación indicial se tiene:


144 A. ANÁLISIS TENSORIAL

Obsérvese que:
(A.9)
Por lo t.ant.o:
uX V = 1tiVjEij/..;e/..; = Eij/..;lliVje/..; (A. 10)
Producto tensorial de dos vectores (diada)

uv = 1t;eiVjej = uivi~iej
Desarrollando esta ecuaci6n se tiene:

u v u¡vje¡ej + U2Vje2ej + U3Vje3ej


u¡v¡e¡e¡ + u¡v2e¡e2 + u¡v3e¡e3
+U2V¡ e2e¡ + u2v2e2e2 + u2 v 3e 2e 3
+1t3V¡ e3e¡ + u3v2e3e2 + u3 v 3e 3e 3

La suma de estos nueve términos se conoce como la forma "nonion" de la


diada ti v.
Frecuentemente se utiliza para el producto tensorial de dos vectores la
notación siguiente:
ti ® v= Uiei ® vjej = UiVjei ® ej
La suma de diadas, expresada por:
ti¡ v¡ + U2V2 + ... + V-NVN
se conoce como diádica.

Productos vector-diada

1) u· (v w) Uie; . (VjejWke/..;) = Uiv;Wk e /..;


2) (u v) . W (1tieiVjej) . w/..;ek = UiVjWjei
3) UX (vw) (1tiei x vjej) Wkek = Eij/";UiVjWkeqe/..;

4) (u v) x W 1t¡ei (vjej x Wkek) = Ej/..;qUiVjWkeieq

Producto diada-diada

(u v) . (w k) = Uiei (vjej . w/..;ek) sqeq

UiVjWjSqeieq
145

Producto vector-tensor

a) v· T = Viei . Tjkeiek = ViTjkOijek = ViTikek


b) T· v= Tijeiej . Vkek = ~jeiOjkVk = ~jVjei
c) Producto tensor - tensor
T· S = ~jeiej • Spqepeq = ~jSjqeíeq

Sunla de tensores
Los tensores cartesianos del mIsmo orden se pueden sumar o restar
componente a componente.

Multiplicación de cada componente de un tensor por un escalar

bi Aaí
Bij AAij

Producto externo de tensores


El producto externo de dos tensores de un orden arbitrario es un tensor
cuyas componentes se forman multiplicando cada componente de uno de
los tensores por todos los componentes del otro. Esta operación da como
resultado un tensor cuyo orden es la suma de los ordenes factores.
aibj = T ij ; Dij T km = <Pijkm
ViFjk = nijk ; éijk V m = 'l/Jíjkm
Contracción de tensores
La contracción de un tensor respecto a dos índices libres es la operación
que consiste en asignar a ambos índices una misma letra como subíndices,
cambiando de esta manera estos índices por seudoíndices.
La contracción produce un tensor que tiene un orden dos veces menor que el
original.
Contracción de T;j ~ T;i

aiVj ~ Uii

Eijak ~ Eijaj = bi

Eíjai Cj

Eiiak dk
146 A. ANÁLISIS TENSORIAL

EiiF,.m Gkm
EijFim H jm
EijFjm ~m

Producto interno
El producto int.erno de dos tensores es el resultado de una contracción, que
da lugar a un índice por cada tensor, realizando después el producto externo
de los do!; tensores.
Producto externo Cont.racción Producto interno
a·b · ---t '/, = J aibi
a:Elfr, ---t j = i a¡Eik = Jk
Eij I.-ro ---t k= j Eijf}m :. G im
EijEI.-ro ---t k=j Eijl!ljm - B,m

La derivada parcial de un campo tensorial respecto a una variable t se ·


simboliza como :t y sigue las mismas reglas del cálculo convencional.

Así, la derivada parcial respecto a la coordenada x q , se indica como !;Jo ,que


uX q
de manera compacta se puede expresar como Oq. La derivada parcial segunda
2
ox q0oX m ,se expresa de manera compacta como Oqm.
Frecuentemente la derivada parcial respecto a la variable Xi se representa en
notación indicial haciendo uso de la coma, esto es:
Oep
a) !:l = ep"
UXi
OVi
e) ~ = Vi,j
UXj
oTij
e) - - = 'nj,k
OXk

Algunos operadores usados en el cálculo vectorial se representan en notación


indicial como sigue:

Nabla; V'
147

Gradiente de un campo escalar:

Divergencia de un campo vectorial

Sea ti
div ti

Rotacional de un campo vectorial ti.

rotv
148 A. ANÁLISIS TENSORIAL
Apéndice B
Rotación de Ejes Coordenados

Las componentes del vector v en el sistema de referencia XI, X2, se pueden


expresar como

e
Si el sistema de referencia original se gira un ángulo respecto al eje Xl, el
nuevo sistema de referencia sera X~, x~. Ahora el vector v se puede expresar
en este nuevo sistema de referencia como

X'2
(i)
X'1

e2
( ~e---'--~x (j)
- 1 1
e2
Figura B.1: Rotación de ejes coordenados.

Primeramente se buscará la transformación de coordenadas

Los vectores base correspondientes al sistema de referencia x~, x~, se pueden


expresar como:

149
150 B. ROTACiÓN DE EJES COORDENADOS

cos oe l+ 8enf)e2 (B.l)


-senfJe l + cos (JC2 (B.2)
el allel + al2C2 = aljej
e; a21 el + a22e2 = aÚej
-f
e·,. (B.3)

El vector v puede ser expresado en las coordenadas sin primas como:


(B.4)

y en el sistema con primas como:

(B.5)

Sustituyendo la eco B.3 en la eco B.5 se obtiene:

v= v~a.¡iej (B.6)
Comparando la eco B.4 con la eco B.6 se obtiene

Vj = V~aij (B.7)
Ahora se busca la transformación inversa Xi = gi (x~, x;)
Los vectores base e) y e2 se pueden expresar como:

el cos fJe; - senBe;


e2 senBe; + cos Be;
alle) + a2)e 2 = ailei
-f -f -f

a12e) + a22e; = ai2 e ;

Generalizando se tiene:
(B.8)
El vector v puede ser expresado en el sistema coordenado con primas como:
(B.9)
y en el sistema coordenado sin primas como:
v= vjej (B. 10)
Sustituyendo la eco B.8 en la eco B.lO se tiene:

(B.11)
151

Comparando la eco B.9 con la eco 8.11 se obtiene:

(8.12)

Las ecs. 8.7 y B.12 representan la ley de transformación para un t.ensor de


orden (1).

Esta ley se puede generalizar para un tensor de orden 2, siendo su ley


de transformación

(8.13)

El término ajl representa la matriz de cosenos directores entre los ejes


coordenados x~ y x j, dada por [G.¡j 1= cos (x~, Xj), tal como se muestra en la
Fig. 8.2.

x;
Figura B.2: a) Transformación de un sistema de coordenadas x}, X2, X3 en
X;, x~ y x;.

Tabla 8.1: Tabla de deformación de coordenadas.


el e') e3
;¿. ;¿? ;¿3

el Xl all a12 a13


e') x? a21 an a23
e3 X3 a31 a32 a33
152 B. ROTACiÓN DE EJES COORDENADOS

,
Bibliografía

Eringen, A.C., 1967, Mechanics of Continua, First Edition, John Wiley and
Sonso .

Fung, Y.C., (1965), Foundations of Solid Mechanics, Printice-Hall.

Func, Y.G, 1994, A first Course in Continuum Mechanics, Third Edition,


Prentice-Hall.

Germain, P., Y Muller, P., 1980, Introduction á la mécanique des milieux


continus, Masson.

Levy, E., 1980, Elementos de Mecánica del Medio Continuo, primera edición,
Ed. Limusa.

Malvern, L.E., 1997, Introduction to the Mechanics of Continuous Medium,


Prentice-Hall.

Mase, G.E., 1978, Mecánica del Medio Continuo, Series Schaum.

Mase, G.T, Mase, G. T., and Mase, G.E., 1999, Continuum Mechanics for
Engineering, Second Edition, CRC Press.

Rodríguez GN., (1981), Apuntes de la clase de Mecánica del Medio Continuo,


División de Estudios de Posgrado, Facultad de Ingeniería, UNAM.

153

También podría gustarte