ÍNDICE

RESUMEN…………………………….………………….….……………………………..3

NOTACIÓN……………....…..…………………………………………………………….4

1. ANTECEDENTES...…...…………………………………………………………..7
1.1 El problema…………………………………………………………………………7
1.2 Objetivos…………………………………………………………………………….8
1.3 Justificación…………………………………………………………………………8
1.4 Marco de referencia…………………………………………………………………9
2. ESTADO DEL ARTE…..………………………………………………………...11
2.1 Concepto de mecánica del medio continuo………...……………………………...11
2.2 Elemento finito triangular………………………………………………………….12
2.3 Ley generalizada de Hooke...……………………………………………………...24
2.4 Relación entre el módulo elástico transversal y el módulo elástico
longitudinal………………………………………...………………………………28
2.5 Matriz de tensiones [ T ] del sólido elástico…..……………………………………..31
2.6 Descripción del movimiento y la deformación…………………………………....37
2.7 Matriz de deformaciones [ ε ] del sólido elástico………………..………………….40
2.8 Modelo constitutivo de un material elástico lineal………………………………...45
2.9 Matriz constitutiva [ D ] del sólido elástico………………………………………...48
2.10 Condición plana de esfuerzos……………………………………………………...53
2.11 Condición plana de deformaciones………………………………………………...56
2.12 Formato matricial del elemento finito……………………………………………..58
2.13 Condiciones de carga………………………………………………………………61
2.14 Matriz de rigidez del elemento finito triangular lineal…………………………….67
2.15 Planteamiento general del problema elástico……………………………………...71
2.16 Formulación del problema elástico en desplazamientos…………………………..79
2.17 Formulación del problema elástico en tensiones………………………………......82
2.18 Método de los elementos finitos aplicado a la teoría de la elasticidad, partiendo de
la teoría de los trabajos virtuales…………………………………………………...84

1
3. DESCRIPCIÓN DEL PROGRAMA………...…...……………………………...90
3.1 Características técnicas del programa……………………………………………...90
3.2 Estructura y manejo del programa…….…………………………………………...91
3.3 Tipos de datos, formatos y matrices……………………………………………….94
3.4 Instrucciones de edición de matrices……………………………………………..106
3.5 Instrucciones de operaciones matriciales básicas………………………………...111
3.6 Instrucciones para crear la matriz de derivadas de funciones de forma………….114
3.7 Instrucciones para crear la matriz de rigidez……………………………………..115
3.8 Instrucciones para crear el vector de fuerzas……………………………………..116
3.9 Instrucciones para crear la matriz de constantes elásticas………………………..118
3.10 Instrucciones para numerar los grados de libertad……………………………….119
3.11 Instrucciones para ensamblar y extraer matrices…………………………………120
3.12 Instrucciones de organización de matrices……………………………………….121
4. EJERCICIOS DE APLICACIÓN……………………....……...….…………...123
4.1 Ejercicio N°1 (Viga)……………………………………………………………...123
4.2 Ejercicio N°2 (Suelo)…………………………………………………………….155
5. CONCLUSIONES Y RECOMENDACIONES…………..……………………181
5.1 Conclusiones……………………………………………………………………...181
5.2 Recomendaciones………………………………………………………………...182
6. BIBLIOGRAFÍA………………..………………………………………............184
7. ANEXOS………………………………...………………………………………185
7.1 Distribución de esfuerzos………………………………………………………...186

2
 RESUMEN

Este trabajo de investigación versa en la creación de un programa que permita hacer el

análisis estructural de cualquier tipo de solido, que goce de las condiciones de un material
lineal e isótropo en una condición plana de esfuerzos y de deformaciones.

Este programa es a código abierto, creado en el lenguaje de programación Visual Basic con
entorno Excel por medio del método de los elementos finitos, el elemento finito utilizado es
el elemento triangular lineal. Se ha elegido este tipo de elemento finito porque es el que se
acomoda geométricamente, el cual permite una mejor discretización de cualquier tipo de
solido a calcularse sin importar su geometría, tanto los esfuerzos como las deformaciones
son constante dentro del elemento finito, lo cual hace más simple su análisis.

La creación del programa a código abierto permite obtener el campo de los

desplazamientos, deformaciones y esfuerzos de un sólido sometido a fuerzas estáticas en
condición plana de esfuerzos o de deformaciones en los vértices del elemento finito,
considerando que las deformaciones son infinitesimales y que el material es elástico lineal
isótropo, para lo cual el programa cuenta con subrutinas capaces de solucionar para cada
problema en particular.
El método del elemento finito es simplemente una técnica numérica para obtener la
solución aproximada de un problema de campo, convierte el manejo de ecuaciones
diferenciales en un conjunto de ecuaciones algebraicas lineales, y su gran aceptación se
debe principalmente a la facilidad con que estas ecuaciones pueden reunirse y resolverse
por computadora.
Tiene por objetivo motivar, estimular y facilitar el aprendizaje del método de los elementos
finitos en el campo de la ingeniería civil, a través del uso de subrutinas preestablecidas.
Generar una herramienta didáctica para su uso en las materias teóricas de elementos finitos
y afines, ya que los métodos de los elementos finitos solo cobra real importancia cuando es
usado en un lenguaje de programación.

3
 NOTACIÓN

ϕ (x; y) Función de aproximación.

ϕi , ϕ j , ϕk Valores nodales del elemento finito triangular.

α 1 , α 2 ,α 3 Constantes a definir de la función de aproximación.

xi , x j , xk Coordenadas x del elemento finito triangular.

yi , y j , yk Coordenadas y del elemento finito triangular.

(e )
A=A Área del elemento finito triangular.

[ N ( e) ] Matriz de funciones de forma.

N (ie ) , N (je ) , N (ke ) Funciones de forma.

{ Φ (e ) } Vector de valores nodales.

∇ Operador nabla.

[ B( e ) ] Matriz de operadores diferenciales.

L 1 , L 2 , L3 Coordenadas de área.

a,b,c Potencias enteras.

σ xx , σ yy , σ zz Esfuerzos normales en dirección xx , yy , zz respectivamente.

ε xx , ε yy , ε zz Deformaciones unitarias longitudinales en dirección xx , yy , zz

respectivamente.

v Coeficiente de Poisson.

E Módulo de elasticidad longitudinal o módulo de Young.

G Módulo de elasticidad transversal o de corte.

[T ] Matriz de tensiones del sólido elástico.

4
τ xy , τ xz , τ yz Esfuerzos cortantes en dirección xy , xz , yz respectivamente.

un
⃗ Vector unitario normal a la superficie.

α ,β ,γ Cosenos directores del vector de esfuerzos.

[ε ] Matriz de deformaciones del sólido elástico.

⃗
U Vector desplazamiento.

⃗
U´ Cambio vectorial del desplazamiento.

γ xy , γ xz , γ yz Deformaciones angulares en dirección xy , xz , yz respectivamente.

[ D] Matriz constitutiva del sólido elástico.

{σ } Vector de esfuerzos.

{ε } Vector de deformaciones.

u, λ Constantes de Lamé.

{U ( e ) } Vector de desplazamientos.

{ a(e ) } Vector de desplazamientos nodales.

[ K ( e) ] Matriz de rigidez.

{ f (se ) } Vector de fuerzas de superficie.

{ f (be ) } Vector de fuerzas de cuerpo o de volumen.

{ f (e ) } Vector de fuerzas.

t Espesor.

⃗ρ Vector peso específico.

ρ x , ρ y , ρz Peso especifico del elemento en dirección x , y , z respectivamente.

5
u,v Componentes del vector de desplazamientos.

∇2 Operador laplaciano.

i^ , ^j, k^ Vectores unitarios direccionales cartesianos.

δWe Trabajo virtual externo.

δWi Trabajo virtual interno.

6
 1. ANTECEDENTES.

“CREACIÓN DE UN PROGRAMA DE ANÁLISIS ESTRUCTURAL”

(A código abierto en visual Basic con entorno Excel)

1.1 El problema.- En los años 60 después de establecidos los conceptos básicos del
análisis matricial, comienza el desarrollo de programas aislados de simulación
numérica con propósito especifico. Desde entonces se han desarrollado muchos
programas de computador con objetivos similares. En la actualidad se encuentran un
gran número de programas que permiten la simulación numérica de problemas
multifísicos con gran eficiencia y comodidad para el usuario. En este tipo de
herramientas un programa previamente compilado realiza el proceso de cálculo y
otro programa denominado post proceso, presenta los resultados, actuando en
conjunto como una “caja negra” para el usuario.
Debido a todo esto se hace necesaria una herramienta informática didáctica a código
abierto en el cual el usuario además de comprender el método, programe y
automatice el procedimiento de cálculo.
Planteamiento.- La inexistencia de una herramienta informática didáctica a código
abierto en el análisis estructural y debido a las limitaciones que se tiene en su
manejo, orienta a realizar el cálculo estructural bajo un programa que simplemente
realice la operación de cálculo y presente los resultados sin ninguna verificación o
en caso extremo hacerlo de forma manual el cual ha ocasionado sistemas
estructurales extremadamente complejos cuyas soluciones no representen
estrictamente los fenómenos más reales y en algunos casos cuyos procedimientos de
cálculo se han prácticamente irrealizables.
Formulación-. Entre la variedad de métodos propuestos y de mayor aceptación en
el estado vigente de la técnica bajo un lenguaje de programación adecuado se
presentan los siguientes métodos: análisis numérico, análisis matricial, diferencias
finitas y elementos finitos.
Sistematización.- De las cuatro alternativas de métodos propuestos se optara por el
método de los elementos finitos tipo triángulo bajo el lenguaje de programación
Visual Basic con entorno Excel

7
1.2 Objetivos:
General.- Crear un programa de análisis estructural a código abierto en Visual
Basic con entorno Excel de elementos finitos tipo triángulo.
Específicos.- Entre los objetivos específicos se establecen los siguientes:
 Determinar las deformaciones y esfuerzos de un sólido sometido a fuerzas
estáticas en condición plana de esfuerzos o de deformaciones en los vértices
del elemento finito tipo triángulo, considerando que las deformaciones son
infinitesimales y que el material es elástico lineal isótropo.
 Motivar, estimular y facilitar el aprendizaje del método de los elementos
finitos en el campo de la ingeniería civil, a través del uso de subrutinas
preestablecidas.
 Analizar cada procedimiento y desarrollar nuevos procesos de cálculo que le
permitan al usuario una transparencia y flexibilidad del programa.
 Generar una herramienta didáctica para su uso en las materias teóricas de
elementos finitos y afines, ya que los métodos de los elementos finitos solo
cobra real importancia cuando es usado en un lenguaje de programación.
 Generar una línea de investigación en análisis numérico utilizando una
metodología establecida.
1.3 Justificación.-
Teórica.- La investigación propuesta busca mediante la aplicación del lenguaje de
programación Visual Basic, enriquecer la base conceptual de programación y el
criterio para el análisis de elementos estructurales mediante el método de los
elementos finitos.
Metodológica.- Para logar el cumplimiento de los objetivos de estudio, se acude a
la creación de un programa informático o software didáctico a código abierto bajo el
método de los elementos finitos tipo triángulo, que expresen el comportamiento más
lógico y real de un elemento estructural.
Práctica.- Esta investigación propuesta ayudara a la mayor parte de los proyectistas
estructurales y usuarios a ser participes activamente del proceso de aprendizaje en el
tema, observando, entendiendo y desarrollando progresivamente cada una de las
etapas de cálculo en forma analítica y con el software, para la determinación de los

8
 esfuerzos y deformaciones en una condición plana de esfuerzos de manera eficiente
y rápida.
1.4 Marco de referencia.-
Teórico.- Los elementos finitos se utilizaron por primera vez en la década de los
años cincuenta en el diseño de un aeroplano, actualmente el método se enseña en
muchas universidades y colegios técnicos como una técnica numérica para el
análisis de esfuerzo. Sin embargo, en los últimos años se ha incrementado el uso de
este enfoque en otras ramas de la ingeniería. Este crecimiento en el área de
aplicación ha revelado lo que realmente es el método del elemento finito, un
procedimiento general para obtener soluciones aproximadas de ecuaciones
diferenciales parciales, y se ha convertido en un tema apropiado para incluirse en un
curso de matemáticas para ingeniería.
Este programa trata principalmente problema de campo, es decir, aquellos
problemas cuya solución implica la solución de ecuaciones diferenciales parciales
con las condiciones de frontera apropiadas en el campo de los desplazamientos y
deformaciones planas.
El método del elemento finito es simplemente una técnica numérica para obtener la
solución aproximada de un problema de campo, convierte el manejo de ecuaciones
diferenciales en un conjunto de ecuaciones algebraicas lineales, y su gran
aceptación se debe principalmente a la facilidad con que estas ecuaciones pueden
reunirse y resolverse por computadora.
El cálculo vectorial es el lenguaje natural para trabajar los problemas de campo, no
se requiere de amplios conocimientos pero si es necesario conocer los conceptos de
gradiente y divergencia.
El cálculo vectorial permite establecer los problemas de campo de manera
compacta, sin la necesidad de introducir un sistema de coordenadas especifico. Sin
embargo, en las soluciones numéricas los vectores se expresan usualmente en forma
de componentes. Las operaciones que se relacionan con dichos vectores están
convenientemente escritas utilizando la notación del álgebra de matrices.

9
El algebra de matrices es el lenguaje natural para la agrupación de ecuaciones
algebraicas lineales y se utiliza mucho en todos los libros sobre el método del
elemento finito.
La potencia y funcionalidad del método de los elementos finitos lo ha convertido en
una herramienta importante para la investigación teórica y experimental en
aplicaciones de la mecánica de sólidos. De igual forma, el método participa cada
vez más en la ingeniería de diseño y construcción en las áreas de Geotecnia,
Hidráulica y principalmente en la Ingeniería Estructural.
La intervención del computador como medio para la aplicación del método de los
elementos finitos en la simulación de problemas con muchos grados de libertad,
exige que el estudiante, además de comprender el método, programe y automatice el
procedimiento de cálculo.
Conceptual. En este esquema de investigación se deja claramente establecido los
siguientes conceptos que definen a un elemento elástico:
 Módulo de Elasticidad: Constante de proporcionalidad que define la relación
entre el esfuerzo y la deformación unitaria, es una medida de la rigidez del
elemento.
 Coeficiente de Poisson: Es la relación que existe entre las deformaciones
unitarias en dirección transversal y longitudinal, tiene un valor numérico que
es único para un material en particular que sea homogéneo e isotrópico.
Las variables de cálculo a determinar son las siguientes:
 Deformaciones: Es la relación que existe entre el desplazamiento y la
longitud inicial del elemento estructural.
 Esfuerzos: Son las fuerzas resultantes que se generan en una determinada
área debido a la acción de las cargas.
Espacial. El presente trabajo de investigación se realizara en la Universidad
Autónoma Juan Misael Saracho.

10
 2. ESTADO DEL ARTE

2.1 Concepto de Mecánica del Medio Continuo.

La mecánica del medio continuo o simplemente mecánica del continuo estudia el

comportamiento mecánico de los sólidos y los fluidos en una escala macroscópica,
ignorando la naturaleza discontinua del material, el cual se considera uniformemente
distribuido en la totalidad del cuerpo.

En la mecánica del continuo, una partícula de material se asocia a un punto del espacio
ocupado por el cuerpo, al cual se le pueden atribuir cantidades de campo como la densidad,
el desplazamiento y la velocidad entre otras. Dichas cantidades corresponden a funciones
continuas en términos de la posición.

La mecánica es la ciencia que describe la interacción entre la fuerza y el movimiento. En

consecuencia, las variables presentes en la mecánica del continuo son, por un lado, las
variables relacionadas con las acciones como las fuerzas por unidad de superficie o de
volumen, y por otro lado, las variables cinemáticas como el desplazamiento, la velocidad y
la aceleración.

Las ecuaciones de la mecánica del continuo se pueden clasificar en dos grupos. El primer
grupo corresponde a las ecuaciones aplicables a todos los materiales que describen las leyes
universales de la física, como por ejemplo, la conservación de la masa y la conservación de
la energía. Las expresiones del segundo grupo conocidas como ecuaciones constitutivas
describen el comportamiento mecánico de los materiales en situaciones interesantes para la
ingeniería o para la física.

Los problemas de la mecánica del continuo se observan dos etapas. En la primera se

establece la formulación de las ecuaciones constitutivas que describan adecuadamente el
comportamiento mecánico de los materiales. Mientras que en la segunda etapa se resuelve
la ecuación constitutiva en conjunto con las ecuaciones generales de la mecánica del
continuo, bajo las condiciones de borde apropiados.

11
 2.2 Elemento Finito Triangular Lineal

El método del elemento finito es simplemente una técnica numérica para obtener la
solución aproximada de un problema de campo, es decir aquellos problemas cuya solución
implica la solución de ecuaciones diferenciales parciales con las condiciones de borde o
frontera apropiadas. Convierte el manejo de ecuaciones diferenciales en un conjunto de
ecuaciones algebraicas lineales, y su gran aceptación se debe principalmente a la facilidad
con que estas ecuaciones pueden reunirse y resolverse por computadora.

El elemento finito triangular lineal es un elemento finito bidimensional de aproximación

lineal de 3 nudos, cuya función de aproximación o interpolación de todos los puntos
interiores del elemento finito triangular están definidas únicamente para el elemento y son
nulas en el exterior del elemento, la cual está definida para todo punto (x ; y) que pertenece
al dominio del elemento Ω (e) se puede escribir en términos de los valores nodales como:

ϕ ( x ; y )=α 1 +α 2 x +α 3 y

(e)
⩝ ΕΩ

ϕ i=ϕ ( xi ; y i )

ϕ i=α 1 +α 2 x i +α 3 y i

12
ϕ j=ϕ ( x j ; y j )

ϕ j=α 1+ α 2 x j +α 3 y j

ϕ k =ϕ ( x k ; y k )

ϕ k =α 1 + α 2 x k + α 3 y k

El cálculo vectorial es el lenguaje natural para tratar los problemas de campo, el álgebra de
matrices es el lenguaje natural para la agrupación de ecuaciones algebraicas lineales, por lo
tanto escribiendo en formato matricial los valores nodales se tiene:

{ } [ ]{ }
ϕi 1 xi y i α1
ϕj = 1 xj y j α2
ϕk 1 xk yk α3

Resolviendo el sistema de ecuaciones aplicando la regla de Cramer, se obtiene:

F1 F F
α 1= ; α 2= 2 y α 3= 3
Δ Δ Δ

Donde las variables F 1, F 2, F 3 y Δ usando la regla de Cramer pueden escribirse como:

| |
ϕi xi yi
F 1= ϕ j xj y j =( x j y k −x k y j ) ϕ i−( x i y k −x k y k ) ϕ j+ ( x i y j−x j y i ) ϕ k
ϕk xk yk

13
F 1=( x j y k −x k y j ) ϕ i + ( x k y i −xi y k ) ϕ j + ( x i y j −x j y i ) ϕ k

| |
1 ϕi yi
F 2= 1 ϕ j y j =−( y k − y j ) ϕ i+ ( y k − y i ) ϕ j−( y j− y i ) ϕ k
1 ϕk yk

F 2=( y j − y k ) ϕ i + ( y k − y i ) ϕ j+ ( y i− y j ) ϕ k

| |
1 xi ϕi
F 3= 1 x j ϕ j =( xk −x j ) ϕ i−( x k −x i ) ϕ j+ ( x j−x i ) ϕ k
1 xk ϕk

F 3=( x k −x j ) ϕi + ( x i−x k ) ϕ j+ ( x j−x i ) ϕ k

[ ]
1 xi yi
Δ= 1 x j y j =( x j y k −x k y j )−( xi y k −x k y i ) + ( x i y j−x j y i )
1 xk yk

Δ=( x j y k −x k y j ) −( x k y i−x i y k ) + ( xi y j −x j y i )

El área de un triangulo en función de las coordenadas de sus vértices es posible expresarlo

en un formato matricial, cuyo valor se obtiene mediante el determinante:

[ ]
1 xi yi
1 1
A= 1 x j y j = {( x j y k − xk y j ) −( x i y k −x k y i ) + ( x i y j −x j y i ) }
2 2
1 xk yk

1
A= {( x y −x y ) −( x k y i −xi y k ) +( x i y j−x j y i ) }
2 j k k j

Como consecuencia de esta ecuación el determinante Δ puede ser expresado en función del
área del triangulo, multiplicando y dividiendo por 2 la ecuación resulta:

Δ=2 {12 [ ( x y −x y )−( x y −x y )+ ( x y −x y )] }

j k k j k i i k i j j i

Δ=2 A

14
Definiendo notación:

a i=( x j y k −x k y j ) ; bi =( y j − y k ) ; ci =( x k −x j )

a j=( x k y i−x i y k ) ; b j =( y k − y i ) ; c j =( x i−x k )

a k =( x i y j−x j y i ) ; bk =( y i− y j ) ; ck =( x j−x i )

Entonces:

F 1 ai ϕ i+ a j ϕ j+ ak ϕ k
α 1= =
Δ 2A

F 2 bi ϕ i+ b j ϕ j+ bk ϕ k
α 2= =
Δ 2A

F3 c i ϕi + c j ϕ j +c k ϕ k
α 3= =
Δ 2A

ϕ =ϕ ( x ; y )=α 1+ α 2 x+ α 3 y

1
ϕ ( x ; y )= {( a ϕ +a ϕ +a ϕ ) + ( bi ϕi+ b j ϕ j+ bk ϕk ) + ( ci ϕ i +c j ϕ j + c k ϕk ) }
2A i i j j k k

1
ϕ ( x ; y )= {( a +b x +c i y ) ϕi +( a j +b j x+ c j y ) ϕ j + ( a k+ bk x+ c k y ) ϕk }
2A i i

Definiendo notación:

( e) ai +bi x+ ci y ( e) a j +b j x+ c j y (e ) a k + bk x+ c k y
Ni ( x ; y )= ; N j ( x ; y )= ; N k ( x ; y )=
2A 2A 2A

Por lo tanto puede escribirse como:

ϕ ( x ; y )=N (e) (e) (e)

i ϕ i+ N j ϕ j + N k ϕ k

En la ecuación anterior, A es el área del elemento, en el cual se puede observar que si el

elemento tiene área nula es decir dos nudos coincidentes, eso se manifiesta con A=0 y no
se podrían calcular las funciones de forma. Estas funciones son planos de valor 1 en uno de
los tres nodos y 0 en los otros dos nodos.

15
Al ser las tres funciones de interpolación planos el cual se aproxima mediante una
superficie poliédrica de facetas triangulares de la siguiente manera:

Función de interpolación para N i:

Función de interpolación para N j:

16
Función de interpolación para N k :

Donde:

N (e) (e) (e)

i , N j y N k → Son denominados “Funciones de forma”.

Matriz de funciones de forma:

[ N (e) ]=[ N (e)

i Nj
(e) (e)
Nk ]
Vector de valores nodales:

{}
ϕ (e)
i
(e) T
={ ϕ i }
(e) (e) (e)
ɸ = ϕ (e)
j
ϕj ϕk
(e)
ϕ k

Escribiendo las funciones de aproximación en formato matricial:

ϕ =[ N ] {ɸ(e) }
(e)

{}
(e)
ϕi
ϕ =[ N(e)
i N (e)
j k ] ϕj
N (e) (e)

(e)
ϕk

Derivando parcialmente respecto de las variables x e y la función de aproximación se

puede escribir como:

17
 (e) (e) (e)
∂ ϕ (x ; y ) ∂ N i (e) ∂ N j (e) ∂ N k (e)
= ϕ + ϕ + ϕ
∂x ∂x i ∂x j ∂x k

(e) (e) (e)

∂ ϕ (x ; y ) ∂ N i (e) ∂ N j (e) ∂ N k (e)
= ϕ + ϕ + ϕ
∂y ∂y i ∂y j ∂y k

En general, se indica con el símbolo ∇ al operador nabla definido como un vector operador
diferencial en términos de las derivadas parciales con respecto a la posición x e y de la
forma:

∂ ⃗ ∂ ⃗
∇= i+ j
∂x ∂ y

Expresado en formato matricial:

[]
∂
[ ∇ ]= ∂ x
∂
∂y

Por lo tanto el operador nabla afectado por la función de aproximación resulta:

[ ][ ]{ }
(e) (e) (e)
∂ϕ ∂ Ni ∂Nj ∂ Nk
ϕ(e)
i
∂x ∂x ∂x
[ ∇ ] ϕ= ∂ x = (e) (e) (e) ϕ(e)
∂ϕ ∂ Ni ∂Nj ∂ Nk j
(e)
∂y ϕk
∂y ∂y ∂y

Definiendo notación:

[ ] {}
∂ N (e) ∂ N (e) ∂ N (e) (e)
i j k
ϕi
∂x ∂x ∂x
[ B ]= (e)
e
(e)
{ (e) }= ϕ (e)j
(e) ; ɸ
∂Ni ∂N j ∂ Nk (e)
ϕk
∂y ∂y ∂y

[ Be ] → Matriz de operadores diferenciales

{ ɸ(e) } → Vector de valores nodales

18
Por lo tanto la matriz de operadores diferenciales puede escribirse como:

[]
∂
[ B(e)]= [ ∇ ] [ N (e) ]= ∂ x [ N(e)
i N (e)
j k ]
N (e)
∂
∂y

[]
∂
[ ∇ ] = ∂ x ; [ N (e) ]=[ N (e)
i Nj Nk ]
(e) (e)
∂
∂y

[ ∇ ] ϕ=[ ∇ ] [ N(e) ] { ɸ(e) }=[ B(e) ] { ɸ(e) }

O también:

T T
[ ∇ ] ϕ={ ɸ(e) } [ B(e) ]
T

Para el caso del elemento finito triangular lineal, operando las derivadas parciales de las
funciones de forma respecto de las variables x e y se obtiene:

∂ N (e)
∂x
i
=
∂x (
∂ ai +bi x+ c i y
2A
b
= i
2A )
( )
(e)
∂N j ∂ a j+ b j x+ c j y b
= = j
∂x ∂x 2A 2A

( )
(e)
∂Nk ∂ ak +b k x +c k y b
= = k
∂x ∂x 2A 2A

∂ N (e)
∂y
i
=
∂y (
∂ ai +b i x +c i y
2A
c
= i
2A )
∂ N (e)
∂y
j
=
∂y (
∂ a j+ b j x +c j y
2A
c
= j
2A )
∂ N (e)
∂y
k
=
∂y (
∂ a k + bk x + ck y
2A
c
= k
2A )
19
En conclusión la matriz de operadores diferenciales resulta:

[ ]
(e) (e) (e)
∂Ni ∂N j ∂ Nk
[ Be ]= ∂ x(e)
∂Ni
∂x
∂ N (e)
j
∂x
(e) =
∂ Nk [
1 bi b j bk
2 A ci c j ck ]
∂y ∂y ∂y

Coordenadas de área.

0 ≤ S i ≤ hi

Para convertir a una expresión adimensional dividimos a toda la inecuación por hi de lo

cual resulta:

0 Si hi
≤ ≤
hi hi hi

Si
0 ≤ Li ≤ 1; Li=
hi

S1 S2 S3
L1 = ; L2 = ; L3=
h1 h2 h3

20
Multiplicando y dividiendo por un medio de sus lados correspondientes a cada caso, se
tiene:

( )
1
S1 b1 2 A
L1 = = 1
h1 b 1 1 A
2

( )
1
S2 b2 2 A
L2 = = 2
h2 b2 1 A
2

( )
1
S3 b3 2 A
L3 = = 3
h3 b3 1 A
2

21
 | |
1 x y
1 1
A 1= 1 xj y j = {( x j y k − y k y i )− ( y k − y i ) x+ ( x k −x j ) y }
2 2
1 xk yk

1
A1= {( x y − y y ) +( y i− y k ) x + ( x k−x j ) y }
2 j k k i

Donde:

a i=( x j y k −x k y j ) ; bi =( y j − y k ) ; ci =( x k −x j )

22
Entonces:

1
A1= ( ai +bi x+ c i y )
2

A 1 ( ai +bi x+ c i y ) (e)
L 1= = =N i
A 2A

| |
1 x y
1 1
A2= 1 x k y k = {( x k y i−x i y k ) −( y i− y k ) x + ( x i−x k ) y }
2 2
1 xi yi

1
2{ k i i k
A2= ( x y −x y ) + ( y k − y i ) x + ( x i−x k ) y }

a j=( x k y i−x i y k ) ; b j =( y k − y i ) ; c j =( x i−x k )

1
A2= ( a j + b j x+ c j y )
2

A 2 ( a j+ b j x +c j y ) (e)
L 2= = =N j
A 2A

| |
1 x y
1 1
A 3 = 1 xi y i = {( x i y j−x j y i )−( y j− y i ) x+ ( x j−x i ) y }
2 2
1 xj yj

1
2{ i j j i
A3 = ( x y −x y ) + ( y i− y j ) x + ( x j−x i ) y }

a k =( x i y j−x j y i ) ;bk =( y i− y j ) ; ck =( x j−x i )

1
A3 = ( ak +b k x +c k y )
2

A 3 ( ak +b k x +c k y ) (e)
L 3= = =N k
A 2A

Por lo tanto se concluye:

23
 Las coordenadas Funciones
≣
de área de Forma

(e)
L1=N i

L2=N (e)
j

(e)
L3=N k

En 1973 Eisenberg y Malvern demostraron que:

a b c a!b! c!
∫ A L1 L2 L3 dA= 2A
( a+ b+c +2 ) !

a b c
L1 L2 L3= Coordenadas de área elevadas a potencias enteras.

2.3 Ley generalizada de Hooke.

La relación de Poisson permite generalizar la aplicación de la ley de Hooke al caso de

deformaciones por esfuerzos triaxiales. Por ejemplo, sea un paralelepípedo diferencial de
lados paralelos a los ejes x-y-z sometido exclusivamente a esfuerzos normales bajo las
siguientes condiciones:

a) σ xx ≠ 0; σ yy =0 y σ zz =0

24
b) σ xx =0; σ yy ≠ 0 y σ zz =0

c) σ xx =0; σ yy =0 y σ zz ≠ 0

25
En materiales elásticos lineales isótropos, la deformación lateral asociada a la contracción
transversal es proporcional a la deformación axial asociada al alargamiento longitudinal.

a) Sí :σ xx ≠ 0 ; σ yy =0 y σ zz =0

´ ´ σ xx
σ xx =E ε xx → ε xx=
E

El cociente entre la deformación en la dirección lateral y la deformación en dirección axial

se conoce como la relación de Poisson y se define como:

´ ´
ε yy ε zz ´ ´ ´ ´
´
= ´
=−ν → ε yy =−ν ε xx y ε zz =−ν ε xx
ε xx ε xx

Recordando que una deformación longitudinal positiva está relacionada con el alargamiento
del sólido, mientras que un valor negativo de deformación longitudinal está asociado al
acortamiento del mismo.

El signo menos en la expresión anterior indica que el alargamiento en la dirección axial está
acompañado de un acortamiento lateral y viceversa.

σ xx
Reemplazando ε ´xx= en ε ´yy y ε ´zz la ecuación queda:
E

26
 ´ σ xx ´ σ xx
ε yy=−ν y ε zz =−ν
E E

b) Sí :σ xx =0 ; σ yy ≠ 0 y σ zz =0

σ yy
σ yy=E ε ´´yy → ε ´´yy=
E

´´ ´´
ε xx ε zz
´´
= ´´
=−ν → ε ´xx´ =−ν ε ´´yy y ε ´´zz =−ν ε ´yy´
ε yy ε yy

σ yy
Reemplazando ε ´´yy = en ε ´´xx y ε ´´zz la ecuación queda:
E

´´ σ yy ´ ´ σ yy
ε xx =−ν y ε zz =−ν
E E

c) Sí :σ xx =0 ; σ yy =0 y σ zz ≠ 0

σ zz
σ zz =E ε ´zz´ ´ → ε ´zz´´ =
E

´´´ ´´ ´
ε xx ε yy
´´´
= ´´ ´
=−ν → ε ´xx´ ´ =−ν ε ´zz´ ´ y ε ´yy´´ =−ν ε ´zz´ ´
ε zz ε zz

Considerando que los esfuerzos cortantes solo producen deformaciones angulares, las
deformaciones longitudinales unitarias obtenidas de la acción de los esfuerzos normales son

La deformada unitaria total sobre el eje x-x:

´ ´´ ´´´
ε xx =ε xx+ ε xx+ ε xx

σ xx σ σ
ε xx = −ν yy −ν zz
E E E

1
ε xx = {σ −ν ( σ yy + σ zz ) }
E xx

La deformada unitaria sobre el eje y-y:

27
 ´ ´´ ´´´
ε yy =ε yy +ε yy +ε yy

σ xx σ yy σ
ε yy =−ν + −ν zz
E E E

1
ε yy = {σ −ν ( σ xx +σ zz ) }
E yy

La deformada unitaria sobre el eje z-z:

´ ´´ ´ ´´
ε zz =ε zz + ε zz + ε zz

σ xx σ yy σ zz
ε zz =−ν −ν +
E E E

1
E { zz
ε zz = σ −ν ( σ xx + σ yy ) }

28
 2.4 Relación entre el módulo elástico transversal y el módulo elástico longitudinal.

El estado de esfuerzos representado en la figura consiste en una tensión σ xx y una

compresión σ yy de la misma magnitud. Para esfuerzos en secciones inclinadas cuyo ángulo
con la vertical es 45° como el de la figura está sometido a esfuerzo cortante puro τ es
numéricamente igual a σ xx y a σ yy. Este esfuerzo cortante deforma al elemento abcd según
el contorno punteado a ´ b´ c ´ d´ de la figura. El ángulo recto en a y c ha disminuido a
(90° - γ), si γ es la distorsión o deformación angular. En el mismo instante, el ángulo recto
en b y d ha aumentado en γ, por lo que ahora su valor es (90° + γ).

Consideremos la deformación del triángulo rectángulo isósceles aob. Como:

σ xx =−σ yy =τ

Deduciendo de las deformaciones ε xx y ε yy en un plano x e y:

a) Sí: σ xx ≠ 0 y σ yy=0

´ ´ σ xx τ
σ xx =E ε xx → ε xx= =
E E

29
 ´
ε yy ´ ´ τ
´
=−ν → ε yy =−ν ε xx=−ν
ε xx
E

b) Sí: σ xx =0 y σ yy ≠ 0

´´ ´´ σ yy −τ
σ yy=E ε yy → ε yy= =
E E

´´
ε xx ´´ ´´ τ
´´
=−ν → ε xx =−ν ε yy =ν
ε yy
E

La deformada unitaria sobre el eje x – x:

´ ´´
ε xx =ε xx+ ε xx

τ τ τ (1+ ν )
ε xx = +ν =
E E E

La deformada unitaria sobre el eje y – y:

´ ´´
ε yy =ε yy +ε yy

τ τ −τ ( 1+ν )
ε yy =−ν − =
E E E

( )
´
γ ob
tagφ=tag 45− =
2 o a´

Del triángulo de deformaciones se tiene que:

b b ´ + o b´ =ob →o b´ =ob−b b ´

´ ´
a a + oa=o a

a a´
ε xx = → a a´ =oa ε xx
oa

´
−b b ´
ε yy = → b b =−ob ε yy
ob

30
Reemplazando estas ecuaciones en tagφ :

γ o b ob−b b ob +ob ε yy ob ( 1+ε yy )

( )
´ ´
tagφ=tag 45− = ´= ´ = =
2 o a a a +oa oa ε xx +oa oa ( 1+ε xx )

Pero oa=ob y reemplazando las deformaciones unitarias ε xx y ε yy :

τ ( 1+ν )
1−
(
γ
tag 45− =
2
1+
E
) =
E−τ ( 1+ ν )
τ ( 1+ ν ) E+τ ( 1+ν )
E

Por identidades trigonométricas:

γ
tag 45−tag
( γ
tag 45− =
2 ) 1+tag 45 tag
2
γ
2

Pero la tag 45=1 y para ángulos muy pequeños γ como son los de distorsión, la tangente
γ γ
coincide con el ángulo expresado en radianes, es decir tag = :
2 2

γ γ
tag 45−tag
1−
(γ
tag 45− =
2 )
1+tag 45 tag
2
γ
=
1+
2 2−γ
=
γ 2+ γ
2 2

Igualando estas dos ecuaciones:

γ
tag 45−tag
2 2−γ
=
γ 2+ γ
1+ tag 45 tag
2

2E-2 τ −2 τν+ γE−γτ −γτν =2E+2 τ+ 2 τν−γE−γτ−γτν

2 τν+ 2 τ=γE

τ E
=
γ 2 ( 1+ν )

31
El esfuerzo de corte puro es τ =Gγ despejando el modulo de elasticidad transversal o
τ
modulo de rigidez al corte G= la ecuación resulta:
γ

E
G=
2 ( 1+ ν )

2.5 Matriz de tensiones [ T ] del sólido elástico.

Para la obtención de la matriz de tensiones se analiza un elemento diferencial de volumen

de lados dx , dy y dz , el cual está sometido a esfuerzos normales como esfuerzos cortantes
en sus respectivas caras o superficies, para lo cual existen esfuerzos de igual magnitud pero
opuestos en dirección de tal manera que se garantice el equilibrio interno como externo en
el elemento diferencial.

Analizando los esfuerzos para cada cara en particular se tiene que:

Cara o superficie “ x ”

32
Cara o superficie “ y ”

Cara o superficie “ z ”

33
Las áreas respectivas en las caras o superficies son:

A z=dx dy ; A y =dx dz ; A x =dy dz

Haciendo un corte al diferencial de volumen se tiene:

34
Garantizando el equilibrio de rotación se tiene que:

∑ M x =0
τ zy ( A z ) ( dz )−τ yz ( A y ) ( dy )=0

τ zy ( dx dy ) ( dz )=τ yz ( dx dz ) ( dy )

τ zy =τ yz

∑ M y =0
τ zx ( A z ) ( dz )−τ xz ( A x ) ( dx )=0

τ zx ( dx dy ) ( dz )=τ xz ( dy dz )( dx )

τ zx =τ xz

∑ M z=0
τ xy ( A x ) ( dx )−τ yx ( A y ) ( dy )=0

τ xy ( dy dz )( dx )=τ yx ( dx dz )( dy )

35
 τ xy =τ yx

El diferencial de área se expresa en forma vectorial como:

⃗
dA=dA u^n

Donde el vector unitario normal a la superficie sobre el cual actúa puede escribirse como:

^ β ^j+ γ k^
u^n =α i+

Donde:

α , β , , γ → cosenos directores

Por lo tanto:

⃗ ^
dA . i=α dA

⃗
dA . ^j=β dA

⃗ ^
dA . k=γ dA

Entonces garantizando el equilibrio de traslación se tiene que:

. ∑ F x =0

σ x dA−σ xx ( α dA )−τ yx ( β dA )−τ zx ( γ dA ) =0

σ x =σ xx α + τ yx β+ τ zx γ

Como:

τ yx =τ xy ; τ zx=τ xz

Entonces:

σ x =σ xx α + τ xy β +τ xz γ

. ∑ F y =0

σ y dA−τ xy ( α dA ) −σ yy ( β dA ) −τ zy ( γ dA ) =0

36
 σ y =τ xy α + σ yy β+ τ zy γ

Como:

τ xy=τ yx ; τ zy =τ yz

Entonces:

σ y =τ yx α +σ yy β + τ yz γ

. ∑ F y =0

σ z dA−τ xz ( α dA ) −τ yz ( β dA )−σ zz ( γ dA )=0

σ z =τ xz α + τ yz β +τ zz γ

Como:

τ xz =τ zx ;τ yz=τ zy

Entonces:

σ z =τ zx α + τ zy β+ τ zz γ

En resumen se tiene que:

σ x =σ xx α + τ xy β +τ xz γ

σ y =τ xy α + σ yy β+ τ zy γ

σ z =τ zx α + τ zy β+ τ zz γ

Escribiendo en un formato matricial:

{ }[ ]{ }
σx σ xx τ xy τ xz α
σy = τ yx σ yy τ yz β
σz τ zx τ zy σ zz γ

37
Donde:

[ ]
σ xx τ xy τ xz
[ T ] = τ yx σ yy τ yz → Matriz de tensiones
τ zx τ zy σ zz

En conclusión la matriz de tensiones es aquella matriz que nos permite determinar todos los
posibles valores de esfuerzos en un punto cualquiera del sólido en función al plano de corte.
Es una matriz simétrica porque se cumple las leyes de la estática para cualquier punto del
solido elástico.

2.6 Descripción del movimiento y la deformación.

Conceptos básicos:

38
A continuación se presentan algunos conceptos básicos de la mecánica del medio continuo
indispensables para describir el movimiento y la deformación:

 El vector posición es aquel que determina la ubicación de un punto en el espacio y

en el tiempo. Como todo vector, las componentes de la posición están referidas a un
sistema coordenado específico.

 El tiempo es el sistema de medida de la secuencia de sucesos representado

mediante una cantidad escalar. En general el instante de tiempo inicial es igual a
cero y se denomina tiempo de referencia t=0 .

 Punto espacial se define como una posición fija en el espacio, es decir, invariante
con respecto del tiempo. En cambio, punto material o partícula es un lugar en el
espacio de posición relativa cambiante con el tiempo, relacionada con el
movimiento de un cuerpo.

 Cuerpo material es el conjunto de partículas que ocupan posiciones diferentes en

el espacio durante el movimiento a lo largo del tiempo.

 El lugar geométrico de las partículas de un cuerpo material con respecto a un

sistema coordenado fijo en un instante de tiempo se define como configuración del
cuerpo en ese instante. La configuración de un cuerpo en el instante de tiempo de
referencia t=0 se denomina configuración material o de referencia Ωo . En
cambio la configuración actual en t indicada como Ωt corresponde a la
configuración de un cuerpo en un instante de tiempo t=t .

 El volumen material se define como aquel volumen móvil conformado siempre por
las mismas partículas, el cual es útil cuando las partículas conservan cierta
proximidad durante el movimiento. En cambio, el volumen de control corresponde
a un volumen fijo en el espacio que es atravesado por las partículas durante su
movimiento y se utiliza cuando las partículas se dispersan con el movimiento.

Campo del desplazamiento y de la velocidad.

39
Se denomina desplazamiento de una partícula al cambio de posición de la misma entre la
configuración de referencia y la configuración actual. Al igual que la posición, el
desplazamiento es una cantidad vectorial, la cual está definida como la resta entre el vector
posición en la configuración actual y el vector posición en la configuración de referencia.

El campo del desplazamiento de un cuerpo es una función que define el vector

desplazamiento para cada partícula que lo conforma, en un instante de tiempo t. recordando
que las coordenadas materiales establecen la posición en la configuración de referencia, la
cual es independiente del tiempo, se puede considerar a la variación del desplazamiento con
respecto al tiempo como una definición alternativa de la velocidad.

Concepto y teorías de la deformación.

En general un cuerpo en movimiento cambia en su posición, orientación, tamaño y forma.

Si tal cuerpo cambia exclusivamente de posición y orientación se denomina cuerpo rígido,
en cambio, si el cuerpo modifica su tamaño o su forma se denomina cuerpo deformable.

La distorsión y el alargamiento o simplemente las distorsiones son los cambios de forma y

de tamaño de un cuerpo después de aplicado un desplazamiento presentado entre la
configuración de referencia o no deformada y la configuración actual o deformada. Las
distorsiones están medidas en unidades de longitud.

La deformación puede definirse como la medida geométrica de las distorsiones que

representa el desplazamiento relativo entre partículas de un cuerpo, es decir, una medida de
la cantidad de desplazamiento que difiere localmente del movimiento del cuerpo rígido. La
deformación es una cantidad adimensional la cual se expresa habitualmente en porcentaje o
fracciones decimales.

De acuerdo con la magnitud de la deformación, el análisis mecánico de sólidos se puede

clasificar según las siguientes teorías:

 En la teoría de deformación infinitesimal o teoría de las pequeñas

deformaciones se considera que los desplazamientos, las rotaciones y las
deformaciones de las partículas de un cuerpo son pequeñas. En este caso, la
configuración actual es muy parecida a la configuración de referencia. Las

40
 deformaciones de los materiales más comunes en la construcción de obras civiles
tales como el concreto y el acero se pueden considerar infinitesimales bajo
condiciones de carga habituales.

 La teoría de deformación finita o teoría de grandes deformaciones establece que

los desplazamientos, las rotaciones y las deformaciones son arbitrariamente
grandes. En tal caso, la configuración actual de un cuerpo es sustancialmente
diferente con respecto a su configuración de referencia. Algunos materiales
descritos mediante la teoría de deformación finita son las gomas, los elastómeros y
los tejidos biológicos suaves.

 La teoría de grandes desplazamientos o teoría de grandes rotaciones considera

que los desplazamientos y las rotaciones pueden ser grandes mientras las
deformaciones son pequeñas.

2.7 Matriz de deformaciones [ ε ] del sólido elástico

41
Sean A y B dos partículas ubicadas inicialmente antes de la acción de las cargas en un
entorno diferencial dentro de un cuerpo material, para la cual el vector definido por la línea
AB es ⃗
dr , y sea A´ y B´ ´ las nuevas posiciones de dichas partículas luego de la acción de
las cargas para la cual el vector definido por la línea A´ B ´ es ⃗
dr .
´

Escribiendo el vector ⃗
dr , resulta:

⃗ ^ d y ^j+ d z k^
dr =d x i+

Escribiendo el vector de desplazamiento ⃗

U en forma vectorial:

⃗ ^ v ( x ; y ; z ) ^j + w ( x ; y ; z ) k^
U =u ( x ; y ; z ) i+

El diferencial ⃗
dU puede expresarse como:

⃗
dU =
∂u
∂x(dx+
∂u
∂y
∂u
dy + dz i+
∂z ∂x ∂y ) (
^ ∂ v dx+ ∂ v dy + ∂ v dz ^j
∂z )
+ ( ∂∂wx dx+ ∂∂wy dy + ∂∂wz ) k^

El cambio vectorial puede escribirse como:

42
 ⃗ ⃗ +⃗
U ´ =U dU

Escribiendo en un formato matricial el cambio vectorial se tiene que:

{⃗
U }= { ⃗
´
U }+[ A ] {⃗
dr }

[ ]
∂u ∂u ∂u

{ }{}
∂x ∂y ∂z

{}
u´ u dx
´ ∂v ∂v ∂v
v = v + dy
∂x ∂y ∂z
w
´
w dz
∂w ∂w ∂w
∂x ∂y ∂z

Escribiendo la matriz [ A ] de otra forma valiéndose de artificios se tiene que:

T T
[A] [A] [A] [A]
[ A ]= + + −
2 2 2 2

[ A ] + [ A ]T [ A ] −[ A ]T
[ A ]= +
2 2

Donde:

T
[ A ] +[ A ]
[ ε ]= → Matriz simétrica
2

[ A ] −[ A ] T
[ H ]= → Matriz antisimétrica
2

Entonces el vector {⃗
U } queda definido como:
´

{⃗
U }= { ⃗
´
U }+[ A ] {⃗
dr }

{⃗
U ´ }= { ⃗
U }+( [ ε ]+[ H ] ) {⃗
dr }

La resultante vectorial puede escribirse como:

{ ⃗R }={⃗
U ´ }+ { ⃗
dr }

43
Ó también:

{⃗ U } + {⃗
R }={ ⃗ dr }
´

Igualando estas dos últimas expresiones se tiene:

{⃗
U }+ { ⃗
´
⃗ }+ {⃗
dr }= {U dr }
´

Reemplazando {⃗
U } en esta última expresión:
´

U }+ ( [ ε ] + [ H ] ) {⃗
( {⃗ dr } ) + { ⃗ U } + {⃗
dr } ={ ⃗ dr ´ }

{⃗
dr }=( [ ε ] + [ H ] + [ I ] ) { ⃗
´
dr }

{⃗
dr }=[ ε ] { ⃗
´
dr } + ( [ H ] + [ I ] ) { ⃗
dr }

{⃗
dr }=[ ε ] { ⃗
´
dr } + [ I ] { ⃗
dr }+ [ H ] { ⃗
dr }

Como:

[ I ] {⃗
dr } → es un vector en la posición inicial de las partículas sin deformarse, por lo
tanto no afecta en la distorsión que experimentan las partículas.

[ H ] {⃗
dr } → es una matriz ortogonal que solo cambia la dirección del vector ⃗
dr y no así
su magnitud, en consecuencia no afecta la distorsión que experimentan las partículas.

44
Por lo tanto según las dos hipótesis anteriores haciendo las respectivas simplificaciones
resulta que la deformación es igual a:

{⃗
dr }=[ ε ] { ⃗
´
dr }

Escribiendo la matriz [ ε ] como:

[ A ] +[ A ]T
[ ε ]=
2

Donde:

[ ]
∂u ∂u ∂u
∂x ∂y ∂z
[ A ]= ∂ v ∂v ∂v
∂x ∂y ∂z
∂w ∂w ∂w
∂x ∂y ∂z

[ ]
∂u ∂v ∂w
∂x ∂x ∂x
T ∂u ∂v ∂w
[A] =
∂y ∂y ∂y
∂u ∂v ∂w
∂z ∂v ∂z

Por lo tanto:

[ ][ ]
∂u ∂u ∂u ∂u ∂v ∂w
∂x ∂y ∂z ∂x ∂x ∂x
[ ε ]= 1 ∂ v ∂v ∂ v 1 ∂u
+
∂v ∂w
2 ∂x ∂y ∂z 2 ∂ y ∂y ∂y
∂w ∂w ∂w ∂u ∂v ∂w
∂x ∂y ∂z ∂z ∂v ∂z

45
 [ ]
( ∂∂ux + ∂∂ ux ) ( ∂∂ uy + ∂∂ vx ) ( ∂∂uz + ∂∂ wx )
[ ε ]= 1 ( ∂ v + ∂ u ) ( ∂ v + ∂ v ) ( ∂ v + ∂ w )
2 ∂x ∂ y ∂y ∂y ∂z ∂ y

( ∂∂wx + ∂u∂z ) ( ∂∂wy + ∂∂ vz ) ( ∂∂wz + ∂∂wz )

[ ]
∂u
∂x (
1 ∂u ∂ v
+
2 ∂ y ∂x ) 12 ( ∂∂ uz + ∂∂wx )
[ ε ]= 1
2 ( ∂v ∂u
+
∂x ∂ y ) ∂v
∂y 2(∂ z ∂ y )
1 ∂v ∂w
+

1
2 ( ∂ w ∂u
+
∂x ∂ z ) ( 1 ∂ w ∂v
+
2 ∂ y ∂z ) ∂∂wz
Como:

∂u ∂v ∂w
ε xx = ε yy= ε zz =
∂x ∂y ∂z

γ xy = ( ∂∂ uy + ∂∂ vx ) γ =( ∂∂uz + ∂∂ wx ) γ =( ∂∂ vz + ∂∂wy )
xz yz

Entonces:

[ ]
1 1
ε xx γ γ
2 xy 2 xz
1 1
[ ε ]= γ ε yy γ → Matriz de deformaciones
2 yx 2 yz
1 1
γ γ ε zz
2 zx 2 zy

En conclusión la matriz de deformaciones es aquella matriz que relaciona el

desplazamiento presentado entre la configuración de referencia no deformada y la
configuración actual o deformada. Es una matriz simétrica porque al ser compatible con los
esfuerzos deben reflejar igual simetría.

46
 2.8 Modelo constitutivo de un material elástico lineal

Modelos constitutivos:

El comportamiento mecánico real de los materiales es bastante complejo, sin embargo es

posible establecer la relación entre el esfuerzo y la deformación durante el proceso de carga
mediante modelos matemáticos denominados modelos constitutivos. La aplicación de
dichos modelos depende del orden de magnitud de la deformación y de las características
mecánicas del material, entre otras.

Los modelos constitutivos describen el comportamiento de un material bajo diferentes

acciones, principalmente mecánicas. La respuesta mecánica de los materiales está
expresada en términos del esfuerzo, la deformación y las variables internas. Estas últimas
describen el efecto de la historia de carga sobre las propiedades actuales.

En general los modelos constitutivos están formados por la ecuación constitutiva, el criterio
de fluencia o de fractura, las leyes de evolución y las condiciones de carga y descarga.

La ecuación constitutiva es la expresión que define al esfuerzo (o el incremento del

esfuerzo) en función de la deformación (o del incremento de la deformación), en un
instante de tiempo y para todo punto material.

El criterio de fluencia o de fractura establece la condición límite para el cual el estado de

esfuerzos o deformaciones de una partícula describe un comportamiento elástico. Para los
materiales dúctiles tal condición se denomina criterio de fluencia y para materiales frágiles
se llama criterio de fractura.

Las leyes de evolución son expresiones que definen la historia en el tiempo de las variables
internas del modelo constitutivo. Como su nombre lo indica, las condiciones de carga y
descarga establecen cuando un punto material se encuentra en régimen de carga o descarga.

A pesar de la diversidad en la representación del comportamiento mecánico de un material,

a continuación se indican algunas características indispensables de todo modelo
constitutivo:

47
  Homogeneidad dimensional. La dimensión de cada uno de los términos de la
relación constitutiva debe ser la misma.

 Sistema coordenado indiferente. Aunque se pueda expresar con respecto a un

sistema coordenado especifico, las relaciones constitutivas son independientes del
sistema coordenado escogido.

 Objetividad material. Las relaciones constitutivas no cambian con respecto a la

traslación y rotación del marco de referencia. Por ejemplo, el movimiento del sólido
rígido no afecta su estado de esfuerzos.

Los tipos de formulación matemática utilizados para obtener las relaciones constitutivas
son los siguientes:

 El formato algebraico relaciona los esfuerzos y de deformaciones mediante la

relación constitutiva de un modelo elástico lineal que está dada de la forma σ =c :ε
donde C la matriz constitutiva elástica y depende de las constantes elásticas del
material y ε el vector de deformaciones.

 El formato integral es una representación funcional donde las componentes de la

matriz de esfuerzos son el resultado de las derivadas parciales de un potencial W ( ε ij )
∂W
con respecto a las componentes de deformación, es decir σ ij=
∂ ε ij

 El formato diferencial establece una relación constitutiva tangente de la forma

dσ =C : dε , definida a partir de elementos diferenciales de esfuerzo dσ y de
tan

deformación dε . Este formato permite almacenar las variables internas las cuales
memorizan los cambios inelásticos de las propiedades materiales.

Material elástico:

Elasticidad es la propiedad física de un sólido de regresar a su configuración de referencia

después de retiradas las acciones impuestas sobre él, en otras palabras, no hay deformación
en el sólido después de la descarga.

48
Un material se denomina elástico si cumple con tal propiedad y además el tensor de
esfuerzos es una función única del tensor de deformaciones. En cambio un material se
define como inelástico si después de la descarga existe una deformación permanente.

Material elástico lineal:

Para un material elástico lineal se establece que el campo tensorial de esfuerzos es una
función lineal del campo de deformaciones. Asimismo, considera que los gradientes de
desplazamiento son pequeños con respecto a la unidad.

Para muchos problemas de ingeniería el modelo elástico lineal ofrece resultados

satisfactorios, tal es el caso de estructuras bajo cargas de servicio con materiales como el
concreto y el acero.

49
 2.9 Matriz constitutiva [ D ] del sólido elástico

Escribiendo en formato matricial la matriz [ D ] que relaciona los esfuerzos con las
deformaciones:

{ σ }=[ D ] { ε }

Donde:

{ σ } → Vector de esfuerzos

[ D ] → Matriz constitutiva

{ ε } → Vector de deformaciones

De la ley generalizada de Hooke:

σ xx −ν σ yy −ν σ zz =E ε xx

−ν σ xx +σ yy −ν σ zz =E ε yy

−ν σ xx −ν σ yy + σ zz =E ε zz

Escribiendo en un formato matricial:

[ ]{ } { }
1 −ν −ν σ xx E ε xx
−ν 1 −ν σ yy = E ε yy
−ν −ν 1 σ zz E ε zz

Resolviendo el sistema aplicando la regla de Cramer:

F1 F2 F3
σ xx = ; σ yy = y σ zz=
Δ Δ Δ

Donde:

| |
E ε xx −ν −ν
F 1= E ε yy 1 −ν =( 1−ν 2 ) E ε xx−(−ν−ν 2) E ε yy + ( ν+ ν 2 ) E ε zz
E ε zz −ν 1

50
F 1=( 1−ν 2 ) E ε xx + ( ν +ν 2 ) E ε yy + ( ν + ν 2) E ε zz

F 1=( 1+ ν ) ( 1−ν ) E ε xx + ( 1+ν ) νE ε yy + ( 1+ ν ) νE ε zz

| |
1 E ε xx −ν
E ε yy −ν =−( −ν−ν ) E ε xx + ( 1−ν ) E ε yy −(−ν−ν ) E ε zz
2 2 2
F 2= −ν
−ν E ε zz 1

F 2=( ν+ ν ) E ε xx + ( 1−ν ) E ε yy ( ν +ν ) E ε zz
2 2 2

F 2=( 1+ν ) ν E ε xx+ ( 1+ ν ) ( 1−ν ) E ε yy+ ( 1+ ν ) νE ε zz

| |
1 −ν E ε xx
F 3= −ν 1 E ε yy =( ν +ν ) E ε xx−(−ν−ν ) E ε yy + ( 1−ν ) E ε zz
2 2 2

−ν −ν E ε zz

F 3=( ν+ ν 2 ) E ε xx+ ( ν + ν2 ) E ε yy + ( 1−ν 2 ) E ε zz

F 3=( 1+ν ) νE ε xx + ( 1+ ν ) νE ε yy + ( 1+ ν ) (1−ν ) E ε zz

| |
1 −ν −ν
Δ= −ν 1 −ν =( 1−ν ) + ( −ν−ν ) ν−( ν + ν ) ν
2 2 2

−ν −ν 1

Δ=( 1−ν )− ( ν + ν ) ν−( ν+ ν ) ν=( 1−ν ) −2 ( ν + ν ) ν=1−3 ν −2 ν

2 2 2 2 2 2 3

Δ=( 1+ ν )( 1−2 ν )( 1+ ν )

Entonces:

F 1 ( 1+ν ) ( 1−ν ) E ε xx + ( 1+ ν ) νE ε yy + ( 1+ ν ) νE ε zz
σ xx = =
Δ ( 1+ ν ) ( 1−2 ν ) ( 1+ ν )

( 1−ν ) E νE νE
σ xx = ε xx + ε yy + ε
( 1+ν ) ( 1−2 ν ) ( 1+ v )( 1−2 ν ) ( 1+ν ) ( 1−2 ν ) zz

E
σ xx =
( 1+ν ) ( 1−2 ν )
{( 1−ν ) ε xx + ν ε yy + ν ε zz }

51
 F2 ( 1+ ν ) ν E ε xx+ ( 1+ ν ) (1−ν ) E ε yy + (1+ ν ) νE ε zz
σ yy = =
Δ ( 1+ ν )( 1−2 ν )( 1+ ν )

νE ( 1−ν ) E νE
σ yy = ε xx + ε yy + ε
( 1+ν ) ( 1−2 ν ) ( 1+v ) ( 1−2 ν ) ( 1+ ν ) ( 1−2 ν ) zz

E
σ yy =
( 1+ν ) ( 1−2 ν )
{ν ε xx+ ( 1−ν ) ε yy +ν ε zz }

F3 ( 1+ν ) νE ε xx+ (1+ ν ) νE ε yy + ( 1+ ν ) ( 1−ν ) E ε zz

σ zz = =
Δ ( 1+ ν ) ( 1−2 ν ) ( 1+ ν )

νE νE ( 1−ν ) E
σ zz = ε xx + ε yy + ε
( 1+ν ) ( 1−2 ν ) ( 1+v ) ( 1−2 ν ) ( 1+ν ) ( 1−2 ν ) zz

E
σ zz =
( 1+ν ) ( 1−2 ν )
{ν ε xx +ν ε yy+ ( 1−ν ) ε zz }

Por otro lado, la deformación angular es proporcional al esfuerzo cortante o tangencial en

dicho plano de la forma:

τ xy τ xz τ yz
γ xy= ; γ xz = ; γ yz =
G G G

Despejando las componentes de esfuerzo cortante se tiene que:

τ xy =G γ xy =
E
2 ( 1+ ν )
γ xy =
E ( 1−2 ν )
2 ( 1+ν ) ( 1−2 ν )
γ xy=
E
( 1+ ν )( 1−2 ν ) 2{
( 1−2 ν )
γ xy }
τ xz =G γ xz =
E
2 ( 1+ν )
γ xz =
E ( 1−2 ν )
2 ( 1+ν ) ( 1−2 ν )
γ xz =
E
( 1+ ν ) (1−2 ν ) 2{
( 1−2 ν )
γ xz }
τ yz=G γ yz =
E
2 ( 1+ ν )
γ yz =
E ( 1−2 ν )
2 ( 1+ν ) ( 1−2 ν )
γ yz =
E
( 1+ ν )( 1−2 ν ) {
( 1−2 ν )
2
γ yz }
Escribiendo las ecuaciones resultantes de los esfuerzos:

E
σ xx =
( 1+ν ) ( 1−2 ν )
{( 1−ν ) ε xx + ν ε yy + ν ε zz }

52
 E
σ yy=
( 1+ν ) ( 1−2 ν )
{ν ε xx+ ( 1−ν ) ε yy +ν ε zz }

E
σ zz =
( 1+ν ) ( 1−2 ν )
{ν ε xx +ν ε yy+ ( 1−ν ) ε zz }

τ xy=
E
( 1+ν ) ( 1−2 ν ){( 1−2 ν )
2
γ xy }
τ xz =
E
( 1+ ν )( 1−2 ν ){( 1−2 ν )
2
γ xz }
τ yz=
E
( 1+ ν )( 1−2 ν ){( 1−2 ν )
2
γ yz }
Escribiendo en un formato matricial:

{} [ ]{ }
σ xx D11 D12 D13 D14 D 15 D16 ε xx
σ yy D21 D22 D 23 D24 D 25 D26 ε yy
σ zz E D31 D32 D33 D34 D 35 D36 ε zz
=
τ xy ( 1+ν ) ( 1−2 ν ) D41 D42 D 43 D44 D 45 D46 γ xy
τ xz D51 D52 D53 D54 D 55 D56 γ xz
τ yz D61 D62 D 63 D64 D 65 D66 γ yz

Por lo tanto la matriz constitutiva [ D ] es igual a la siguiente expresión:

[ ]
1−v v v 0 0 0
v 1−v v 0 0 0
E v v 1−v 0 0 0
[ D ]=
( 1+ ν ) (1−2 ν ) 0 0 0 ( 1−2 v ) /2 0 0
0 0 0 0 (1−2 v ) /2 0
0 0 0 0 0 ( 1−2 v ) /2

La matriz constitutiva de un material lineal isótropo está definida por las constantes
elásticas y las componentes de esfuerzo y de deformación las cuales también se pueden
expresar en función a las constantes de Lamé que son las siguientes:

E νE
μ=G= , λ=
2 ( 1+ ν ) ( 1+ ν )( 1−2 ν )

53
Escribiendo los esfuerzos en función a las constantes de Lamé resulta:

σ xx = ( λ+2 μ ) ε xx + λ ε yy + λ ε zz

σ yy =λ ε xx + ( λ+2 μ ) ε yy + λ ε zz

σ zz = λ ε xx+ λ ε yy + ( λ+2 μ ) ε zz

τ xy=μ γ xy

τ xz =μ γ xz

τ yz=μ γ yz

Escribiendo en notación de Voigt:

{ }[ ]{ }
σ xx ( λ+2 μ ) λ λ 0 0 0 ε xx
σ yy λ ( λ +2 μ ) λ 0 0 0 ε yy
σ zz λ λ ( λ+2 μ ) 0 0 0 ε zz
=
τ xy 0 0 0 μ 0 0 γ xy
τ xz 0 0 0 0 μ 0 γ xz
τ yz 0 0 0 0 0 μ γ yz

En conclusión la matriz constitutiva, es aquella matriz que relaciona los esfuerzos con las
deformaciones, para un material elástico lineal e isótropo solo depende de dos constantes
que son: el coeficiente de Poisson y el Módulo de elasticidad longitudinal o Módulo de
Young. Es simétrica debido a que las matrices de esfuerzos y deformaciones reflejan igual
simetría.

54
 2.10 Condición Plana de Esfuerzos

El estado plano o condición plana de esfuerzos es una simplificación del estado

tridimensional, en el cual las componentes de esfuerzo contenidas en un plano coordenado
especifico pueden ser diferentes de cero, mientras que las componentes de esfuerzo que no
están contenidas en un plano especifico o en planos paralelos son nulas.

Por lo tanto, si el sólido describe una condición plana de esfuerzos en el plano xy los
esfuerzos y deformaciones resultantes son las siguientes:

σ xx ≠ 0→ ε xx ≠ 0

σ yy ≠ 0 → ε yy ≠ 0

σ zz =0→ ε zz ≠ 0

τ xy ≠ 0→ γ xy ≠ 0

τ xz =0 → γ xz =0

55
 τ yz=0→ γ yz =0

Se observa que a pesar que las componentes de esfuerzo fuera del plano xy son nulas, la
componente de deformación longitudinal ε zz en la dirección z es diferente de cero y
depende de las otras dos componentes de deformación longitudinal.

Entonces las ecuaciones de las deformaciones longitudinales resultan:

1
ε xx= {σ −ν σ yy }
E xx

1
ε yy = {σ −ν σ xx }
E yy

Escribiendo en un formato matricial:

[−ν1 −ν1 ] {σσ }={EE εε }

xx

yy
xx

yy

Resolviendo el sistema aplicando la regla de Cramer:

F1
σ xx =
Δ

F2
σ yy =
Δ

Donde:

Δ= |−ν1 −ν1 |=1−ν 2

F 1=
| E ε xx −ν
E ε yy 1 |
=E ε xx + νE ε yy

56
F 2=
| 1
−ν
E ε xx
E ε yy|=E ε yy +νE ε xx

Entonces:

F 1 E ε xx +νE ε yy E νE E
2 ( xx
σ xx = = 2
= ε +
2 xx
ε =
2 yy
ε +ν ε yy )
Δ 1−ν 1−ν 1−ν 1−ν

F2 E ε yy + νE ε xx νE E E
2(
σ yy= = 2
= ε +
2 xx
ε =
2 yy
νε xx +ε yy )
Δ 1−ν 1−ν 1−ν 1−ν

τ xy =G γ xy =
E
2 ( 1+ ν )
γ xy( )
1−ν
=
E 1−ν
γ
1−ν 1−ν 2 2 xy ( )
Escribiendo en un formato matricial:

{} [ ]{ }
σ xx 1 ν 0 ε xx
E ν 1 0
σ yy = ε yy
1−ν 0 0 1−ν
2
σ zz 2 γ xy

La deformada en la dirección z es:

−v
ε zz = ( σ +σ )
E xx yy

Pero:

E
σ zz =
( 1+ν ) ( 1−2 ν )
{ν ε xx +ν ε yy+ ( 1−ν ) ε zz }=0

ν
ε zz = (ε +ε )
ν −1 xx yy

Luego:

−v ν
ε zz =
E
( σ xx +σ yy )= ( ε +ε )
ν −1 xx yy

57
 2.11 Condición Plana de Deformaciones

El estado plano o condición plana de deformaciones es una simplificación del estado

tridimensional, en el cual las componentes de deformación contenidas en un plano
coordenado especifico pueden ser diferentes de cero y en aquel en el cual las componentes
de deformación que no están contenidas en un plano específico (o en planos paralelos) son
nulas.

Por lo tanto, si el sólido describe una condición plana de deformaciones en el plano xy las
deformaciones y esfuerzos resultantes son los siguientes:

ε xx ≠ 0→ σ xx ≠ 0

ε yy ≠ 0 → σ yy ≠ 0

ε zz =0→ σ zz ≠ 0

58
 γ xy ≠ 0 → τ xy ≠ 0

γ xz=0→ τ xz =0

γ yz =0 → τ yz =0

Se observa que a pesar que las componentes de deformación fuera del plano xy son nulas,
la componente de esfuerzo normal en la dirección z es diferente de cero y depende de las
otras dos componentes de esfuerzo normal.

Entonces los esfuerzos resultantes son:

( 1−ν ) E νE
σ xx = ε xx + ε
( 1+ν ) ( 1−2 ν ) ( 1+ v )( 1−2 ν ) yy

E
σ xx =
( 1+ν ) ( 1−2 ν )
{( 1−ν ) ε xx + ν ε yy }

νE ( 1−ν ) E
σ yy= ε xx + ε
( 1+ν ) ( 1−2 ν ) ( 1+v ) ( 1−2 ν ) yy

E
σ yy=
( 1+ν ) ( 1−2 ν )
{ν ε xx+ ( 1−ν ) ε yy }

E E ( 1−2 ν )
τ xy =G γ xy = γ xy = γ
2 ( 1+ ν ) 2 ( 1+ν ) ( 1−2 ν ) xy

τ xy=
E
( 1+ν ) ( 1−2 ν ) {
( 1−2 ν )
2
γ xy }
Escribiendo en un formato matricial:

[ ]{ }
( 1−ν )

{}
ν 0
σ xx ε xx
E ν ( 1−ν ) 0
σ yy = ε
( 1+ ν )( 1−2 ν ) ( 1−2 ν ) yy
τ xy 0 0 γ xy
2

El esfuerzo en la dirección z es:

59
 νE
σ zz = {ε + ε }
( 1+ν ) ( 1−2 ν ) xx yy

Pero:

1
E { zz
ε zz = σ −ν ( σ xx + σ yy ) }=0

σ zz =ν ( σ xx + σ yy )

Luego:

νE
σ zz = {ε + ε }=ν ( σ xx + σ yy )
( 1+ν ) ( 1−2 ν ) xx yy

2.12 Formato matricial del elemento finito

El vector de desplazamientos dentro del elemento finito triangular puede expresarse como:

{U (e) }=[ N (e) ] {a(e) }

Donde:

{U (e) }→ Vector de desplazamientos

[ N (e) ] → Matriz de funciones de forma

{ a(e) } →Vector de desplazamientos nodales

Por lo tanto:

{ }
(e)
{U (e) }= u(e)
v

60
La matriz de funciones de forma en un formato matricial puede expresarse como:

[
[ N (e) ]= N (e)
1

0
0
N(e)
1
N (e)
2
2

0 N (e) 0
0 N (e)
3 0
N(e)
3
]

El vector de desplazamientos nodales en un formato matricial puede escribirse como:

{}
(e)
u1
v(e)
1
(e)
{ a }= u(e)
(e) 2

v2
u(e)
3

v(e)
3

Entonces:

{U (e) }=[ N (e) ] {a(e) }

61
 {}
u(e)
1

v(e)
1

{ }[
u(e)
v (e) =
N(e)
1 0 N (e)
(e)
2 0 N (e)
(e)
3

0 N 1 0 N 2 0 N 3 v2
0 u(e)
(e)
2
(e)
] (e)
u3
(e)
v3

Expresando las matrices y vectores para el elemento finito triangular:

[ ]
(e)
Ni 0
[ N ]=
(e)
i (e)
0 Ni

[ ]
bi 0
1
[ B ]= (e) 0 c i
(e)
i
2A
ci bi

{ a }=
(e)
i
{ } u(e)
i

v (e)
i

[ N (e) ]=[ [ N (e)

1 ] [ N2 ] [ N3 ]]
(e) (e)

[ B(e) ]= [ [ B(e)
1 ] [ B(e)2 ] [ B(e)3 ] ]

[{ }]
(e)
a1
{ a(e) }= a(e)
2
(e)
a3

El vector de deformaciones también puede escribirse como:

{ ε(e) }=[ B(e) ] {a(e) }

Donde la matriz de rigidez puede escribirse como:

T
[ K (e) ]=∫ [ B(e) ] [ D(e) ][ B(e) ] dV
(e)
V

62
La fuerza en el elemento finito triangular queda expresado como:

T T
{ f (e) }=∫ [ N(e) ] { P(e) } dS+∫ [ N (e) ] {b(e) } dV
Γ V

Donde:

T
{ f (e)s }=∫ [ N (e) ] {P(e)} dS → Vector de fuerzas de superficie
Γ

(e) T
b }=∫ [ N
{ f (e) ] {b(e) }dV → Vector de fuerzas de cuerpo
V

Entonces:

T T
{ f (e) }=∫ [ N(e) ] { P(e) } dS+∫ [ N (e) ] {b(e) } dV ={ f (e)
s }+ { f b }
(e)

Γ V

2.13 Condiciones de carga

Las condiciones de carga que intervienen en el elemento finito triangular lineal en una
condición plana de esfuerzos y deformaciones son las fuerzas de superficie y las fuerzas de
cuerpo o de volumen, las cuales dependen del tipo de problema, la geometría y la función
de aproximación del elemento.

Fuerzas de superficie:

Las fuerzas de superficie son acciones distribuidas en una superficie específica en el

contorno del sólido o dentro de él, que actúan por contacto directo con otro cuerpo.

El vector de fuerzas de superficie puede escribirse como:

63
 T
{ f (e)s }=∫ [ N (e) ] {P(e) } dS
Γ

Donde:

[ ]
(e)
N1 0
(e)
0 N1
(e)
= N2 0
T
[ N (e) ] (e)
0 N2
N (e)
3 0
(e)
0 N3

{ }
(e)
Px
{ P }=(e)
(e)
Py

El diferencial de superficie puede expresarse como:

dS=t dL

Entonces el vector de fuerzas de superficie en las tres caras puede escribirse de la siguiente
forma:

T T T
{ f (e)s }=∫ [ N (e) ] { P(e)1−2 } dS+∫ [ N (e)] { P(e)
2−3 } dS+∫ [ N ] { P1−3 } dS
(e) (e)

Γ Γ Γ

64
 [ ] [ ]
(e) (e)
N1 0 N1 0
(e) (e)
0 N 0 N1

{ } { }
1
(e) (e)
0 P 0 P x tdl
(e) (e)
{ f }= ∫
(e) N2 x1−2
tdL + ∫ N2 2−3
s (e) (e) (e) (e)
(e)
L
1−2
0 N2 P y 1−2 L
( e)
2−3
0 N2 Py 2−3

N (e)
3 0 N (e)
3 0
(e) (e)
0 N3 0 N3

[ ]
(e)
N1 0
(e)
0 N1

{ }
(e)
0 P x tdL
(e)
+∫ N2 1−3

(e) (e)
L
(e)
1−3
0 N2 Py 1−3
(e)
N3 0
(e)
0 N3

Se va trabajar con presión uniforme es decir:

(e) (e)
P x y P y =ctte
2−3 2−3

Por lo que:

P(e) (e)
x L2−3
f s =f s = 2−3

x2 x3
2

(e) (e)
P y L2−3
f s =f s = 2−3

y2 y3
2

65
Cualquier elemento finito triangular nunca va a tener una distribución de esfuerzos en las 3
caras, asumiendo como máximo en 2 caras, para demostrar las anteriores 2 ecuaciones se
tiene que:

a!b! c!
∫ La1 Lb2 L3c dA= ∫ N a1 N 2b N c3 = ( a+ b+c +2 ) ! 2 A
(e) ( e)
A A

Como las coordenadas de área son iguales a las funciones de forma se tiene que:

N i=l i=coordenadas de área

Entonces:

a! b !
∫ Lai Lbj dL=∫ N ai N bj dL= ( a+b +1 ) ! L(e)ij
(e) (e)
Lij Lij

1! 0! 1
∫ N (e)
i dL=∫ N j dL=
(e)
(e) (e)
L = L
(e)

( 1+0+1 ) ij 2 ij
(e)
L ij L ij

Reemplazando esta última expresión en { f s } e integrando se tiene:

(e)

{} {} {}
(e)
Px 1−2
0 (e)
Px 1−3

P(e) 0 P (e)
y 1−2 (e) y1−3
Px
{ f (e)s }= 12 L(e)1−2 t 1 1
(e)
P + L(e)
x 1−2 t (e)
2−3
+ L(e) t 0
2 2−3 P y
(e) 2 1−3 0
P y 1−2
2−3

0 P(e)
x 2−3
P(e)
x 1−3

(e) (e)
0 Py 2−3
Py 1−3

Fuerzas de cuerpo o de volumen:

Las fuerzas de cuerpo son acciones distribuidas en el volumen del sólido las cuales actúan a
distancia como por ejemplo las fuerzas gravitacionales, inerciales o magnéticas.

El vector de fuerzas de cuerpo puede escribirse como:

66
 T
{ f (e)b }=∫ [ N (e) ] {b(e) }dV
V

Donde:

[ ]
(e)
N1 0
(e)
0 N1
(e)
= N2 0
T
[ N (e) ] (e)
0 N2
(e)
N3 0
(e)
0 N3

{ }
(e)
bx
{b }=(e)
(e)
by

El diferencial de volumen puede expresarse como:

dV =t dA

Entonces el vector de fuerzas de cuerpo en el elemento finito, puede escribirse como:

T
{ f (e)b }=∫ [ N (e) ] {b(e) }dV
V

67
 [ ]
(e)
N1 0
(e)
0 N1

{ }
(e)
0 b x tdA
(e)
{ f }= ∫
(e) N2
b (e) (e)
A
( e)
0 N2 by
N(e)
3 0
(e)
0 N3

Utilizando sub-matrices:

()
{f } (e)
b1

{ f }= { f }
(e)
b
(e)
b2

{f } (e)
b3

(e) T
{ f }= ∫ [ N ] {b } tdA para i=1 ,2 , 3
(e)
bi
(e)
i
(e)

Sustituyendo las funciones de forma:

1
N (e)
i = ( e) ( i
a + bi x +c i y )
( x; y )
2A

Integrando la siguiente expresión se tiene que:

1! 0 !0 ! 1
∫ N (e)
i dA= ∫ ( N i ) ( N j ) ( N k ) dA=
(e) (e) 1
(e) (e)0 (e) 0
2A = A
A
(e)
A
(e) ( 1+ 0+0+2 ) ! 3

68
Entonces el vector de fuerzas de cuerpo resulta:

{}
b(e)
x

b(e)
y

1 (e) b(e)
{ f (e)
b }=
3
A t x
(e)
by
(e)
bx
(e)
by

Esta expresión demuestra que las fuerzas de cuerpo se obtienen multiplicando el peso
específico por el espesor t y el área del elemento finito triangular y se lo divide entre 3, lo
cual significa que se reparte a cada nodo equitativamente, sin importar la forma del
triangulo, en consecuencia se reparte en magnitudes iguales a los 3 nodos, lo cual es un
caso interesante ya que demostraría que ningún vértice de un triángulo es privilegiado ya
que:

x 1+ x 2 + x 3
x=
3

y1+ y2 + y3
y=
3

2.14 Matriz de Rigidez del Elemento Finito Triangular Lineal

Del planteamiento de los trabajos virtuales aplicable a la mecánica del medio continuo y el
análisis estructural, el cual establece que un cuerpo está en equilibrio si y solo si, el trabajo

69
virtual realizado por todas las acciones internas y externas es nulo, se define como la matriz
de rigidez del elemento como:

T T
[ K (e) ]=∫ [ B(e) ] [ D(e) ][ B(e) ] dV = ∫ [ B(e) ] [ D(e) ] [ B(e) ] t dA
(e) (e)
V A

Como [ B(e) ] y [ D (e) ] son matrices constantes dentro del elemento finito la ecuación anterior
T

queda como.

T
[ K ( e ) ]=[ B (e ) ] [ D(e ) ][ B(e ) ] t ∫ dA
(e)
A

T
[ K ( e ) ]=[ B (e ) ] [ D(e ) ][ B(e ) ] t A(e)

Donde:

Escribiendo la matriz de operadores diferenciales en un formato matricial:

[ ]
b1 0 b 2 0 b 3 0
1
[ B ]= (e) 0 c 1 0 c 2 0 c3
(e)

2A
c 1 b1 c 2 b2 c 3 b 3

Escribiendo la matriz de constantes elásticas en un formato matricial:

70
 [ ]
D(e)
11 D(e)
12 0
[ D ]=
(e)
D (e)
12 D(e)
22 0
0 0 D (e)
33

Escribiendo la matriz transpuesta de la matriz de operadores diferenciales se tiene:

[ ]
b1 0 c1
0 c1 b1
T 1 b2 0 c2
[ B(e) ] = (e)
2A 0 c2 b2
b3 0 c3
0 c3 b3

Entonces:

T
[ K ( e ) ]=[ B (e ) ] [ D(e ) ][ B(e ) ] t A(e)

Escribiendo el producto matricial entre matrices se tiene que:

[ ]
b1 0 c1

[ ]
0 c1 b1 (e) (e)
D11 D12 0
T 1 b2 0 c2
[ B( e ) ] [ D ( e ) ] = (e)
(e) (e)
D12 D22 0
2A 0 c2 b2 (e)
b3 0 c3 0 0 D33
0 c3 b3

[ ]
b 1 D(e)
11 b1 D(e)
12 c 1 D(e)
33

c1 D(e)
12 c 1 D(e)
22 b1 D(e)
33
(e) (e) (e)
T 1 b 2 D11 b2 D12 c 2 D33
[ B( e ) ] [ D ( e ) ] = (e) (e) (e) (e)
2 A c2 D12 c 2 D22 b2 D33
(e) (e) (e)
b3 D11 b3 D12 c 3 D33
(e) (e) (e)
c 3 D12 c 3 D 22 b 3 D33

71
 [ ]
b1 D(e)
11 b1 D(e)
12 c 1 D(e)
33

c 1 D(e) c 1 D(e) b 1 D(e)

[ ]
12 22 33
(e) (e) (e) b1 0 b 2 0 b 3 0
1 b2 D11 b2 D12 c 2 D33 1
[ K ]= (e)
( e)
(e) (e) (e) (e) 0 c1 0 c 2 0 c 3 t A
(e)

2 A c 2 D12 c 2 D22 b 2 D33 2 A

(e) (e) (e)
c 1 b 1 c 2 b2 c3 b3
b3 D11 b3 D12 c 3 D33
(e) (e) (e)
c 3 D12 c 3 D22 b 3 D33

En conclusión la matriz de rigidez [ K (e) ] operando el producto matricial entre la matriz de

operadores diferenciales con la matriz constitutiva se tiene que:

[ ]
K 11 K 12 K 13 K 14 K 15 K 16
K 21 K 22 K 23 K 24 K 25 K 26
T 1 K 31 K 32 K 33 K 34 K 35 K 36
[ K (e) ]=[ B (e ) ] [ D(e ) ][ B (e ) ] t A(e) =
4 A(e) K 41 K 42 K 43 K 44 K 45 K 46
K 51 K 52 K 53 K 54 K 55 K 56
K 61 K 62 K 63 K 64 K 65 K 66

Donde:

K 11=b 21 D(e) 2 (e) (e) (e) (e) (e)

11 + c1 D 12 K 12 =b1 c 1 D 12 + b1 c1 D 33 K 13 =b1 b2 D 11 +c 1 c 2 D 33

(e) (e) (e) (e) (e) (e)

K 14 =b1 c2 D12 +b 1 c 2 D 33 K 15=b1 b3 D11 +c 1 c3 D 33 K 16=b 1 c 3 D12 +b3 c 1 D33

(e) (e) 2 (e) 2 (e) (e) (e)

K 21 =b 1 c 1 D 12 +b1 c 1 D 33 K 22=c 1 D22 +b1 D33 K 23=b 2 c 1 D 12 +b1 c 2 D33

(e )
K 24 =c 1 c 2 D(e) (e) (e) (e) (e)
22 +b1 b2 D 33 K 25 =b 3 c1 D 12 +b1 c 3 D 33 K 26 =c 1 c 3 D 22 +b1 b3 D 33

(e) (e) (e) (e) 2 (e) 2 (e)

K 31 =b 1 b 2 D11 + c1 c 2 D 33 K 32=b 2 c 1 D 12 +b1 c 2 D33 K 33 =b2 D11 + c 2 D33

(e) (e) (e) (e) (e) (e)

K 34 =b2 c2 D12 +b 2 c 2 D33 K 35=b2 b3 D11 +c 2 c 3 D33 K 36=b 2 c 3 D12 +b 3 c 2 D33

K 41 =b1 c2 D(e) (e) (e) (e) (e) (e)

12 +b 2 c 1 D 33 K 42 =c 1 c 2 D 22 +b1 b2 D 33 K 43 =b 2 c 2 D 12 +b2 c 2 D 33

2 (e) 2 (e) (e) (e) (e) (e)

K 44 =c 2 D22 + b2 D 33 K 45=b3 c 2 D12 + b2 c 3 D33 K 46 =c 2 c 3 D22 +b 2 b3 D33

(e) (e) (e) (e) (e) (e)

K 51 =b 1 b3 D11 + c1 c 3 D33 K 52=b 3 c1 D12 +b 1 c 3 D 33 K 53=b2 b3 D11 +c 2 c3 D33

72
 (e) (e) 2 (e) 2 (e) (e) (e)
K 54 =b3 c 2 D12 +b 2 c 3 D33 K 55=b3 D11 + c3 D33 K 56 =b 3 c 3 D12 +b3 c 3 D 33

(e) (e) (e) (e) (e) (e)

K 61 =b 1 c 3 D12 +b3 c 1 D33 K 62=c1 c 3 D22 +b1 b3 D 33 K 63=b2 c 3 D12 + b3 c 2 D33

K 64 =c 2 c 3 D(e) (e) (e) (e) 2 (e) 2 (e)

22 +b2 b3 D 33 K 65 =b3 c 3 D 12 +b3 c 3 D 33 K 66 =c 3 D 22 + b3 D 33

En conclusión la matriz de rigidez, es aquella que relaciona los elementos de la matriz de

operadores diferenciales junto con los elementos de la matriz constitutiva, es una matriz
simétrica, debido a que su demostración a sido por métodos energéticos lo cual significa
que el trabajo virtual externo es igual al trabajo virtual interno.

Recordar que:

1
N (e)
1 = (e) ( 1
a +b 1 x +c 1 y )
2A

1
N (e)
2 = (e) ( 2
a +b 2 x +c 2 y )
2A

(e) 1
N3 = ( a3 +b 3 x +c 3 y )
2 A(e)

a 1=x 2 y 3−x 3 y 2 b 1= y 2− y 3 c 1=x 3−x 2

a 2=x 3 y 1−x 1 y3 b 2= y 3− y 1 c 2=x 1−x3

a 3=x 1 y 2−x 2 y 1 b3 = y 1− y 2 c3 =x 2−x 1

( e) (e ) (e ) ( e) ( e) (e ) (e )
U ( x ; y )=N 1 U 1 + N 2 U 2 + N 3 U 3

( e) (e ) (e ) (e ) (e ) ( e) ( e)
V ( x ; y )=N 1 V 1 + N 2 V 2 + N 3 V 3

73
 2.15 Planteamiento General del Problema Elástico

El problema elástico lineal consiste en obtener, para todo instante de tiempo, el

desplazamiento, la deformación infinitesimal y el esfuerzo de cada una de las partículas que
conforma un sólido de material elástico lineal isótropo, el cual está sometido a diferentes
fuerzas externas.

Ecuaciones de equilibrio interno o ecuación de Cauchy.

Se considera un cuerpo en equilibrio si la fuerza y el momento resultante alrededor de un

punto específico son iguales a cero. Asimismo, una parte cualquiera del cuerpo de volumen
y de superficie de contorno también estará en equilibrio, si todo el sólido cumple con tal
condición.

Analizando los esfuerzos en las caras correspondientes en un diferencial de volumen se

tiene que:

Superficie o cara “ x ”

74
Superficie o cara “ y ”

Superficie o cara “ z ”

Las fuerzas que actúan sobre una parte del cuerpo corresponden a la suma de las fuerzas
como el peso propio, expresando en forma vectorial:

75
 ⃗ ^ dF ^j+ dF k^
dF =dF X i+ y Z

⃗ ^ ρ dV ^j+ ρ dV k^
dF =ρx dV i+ y z

La primera condición de equilibrio estático indica que la fuerza resultante es igual a cero y
se puede expresar mediante las siguientes ecuaciones:

∑ F x=0
∂ τ yx ∂τ ∂σ
τ yx+ dy A y −τ yx +τ zx + zx dz A z−τ zx +σ xx + xx dx A x −σ xx + ρ x dV =0
∂y ∂z ∂x

∂ τ yx ∂τ ∂σ
dy ( dxdz ) + zx dz ( dydx ) + xx dx ( dydz )+ ρx ( dxdydz )=0
∂y ∂z ∂x

∂ σ xx ∂ τ xy ∂ τ xz
ρ x+ + + =0
∂ x ∂ y ∂z

∑ F y =0
∂ σ yy ∂ τ zy ∂ τ xy
∂ yy + dy A y −σ yy + τ zy + dz A z −τ zy + τ xy + dx A x −τ xy + ρ y dV =0
∂y ∂z ∂x

∂ σ yy ∂ τ zy ∂ τ xy
dy ( dxdz ) + dz ( dydx ) + dx ( dydz ) + ρ y ( dxdydz )=0
∂y ∂z ∂x

∂ τ yx ∂ σ yy ∂ τ yz
ρy+ + + =0
∂x ∂y ∂z

∑ F z=0
∂ τ yz ∂σ ∂τ
τ yz + dy A y −τ yz + σ zz + zz dz A z −σ zz +τ xz + xz dx A x −τ xz + ρ z dV =0
∂y ∂z ∂x

∂ τ yz ∂σ ∂τ
dy ( dxdz )+ zz dz ( dydx ) + xz dx ( dydz ) + ρ z ( dxdydz )=0
∂y ∂z ∂x

∂ τ zx ∂ τ zy ∂ σ zz
ρz+ + + =0
∂x ∂ y ∂z

76
 Por lo tanto las tres ecuaciones de equilibrio interno o ecuación de Cauchy son las
siguientes:
∂ σ xx ∂ τ xy ∂ τ xz
ρ x+ + + =0
∂ x ∂ y ∂z

∂ τ yx ∂ σ yy ∂ τ yz
ρy+ + + =0
∂x ∂y ∂z

∂ τ zx ∂ τ zy ∂ σ zz
ρz+ + + =0
∂x ∂ y ∂z

Componentes de la matriz de deformación.

La matriz de componentes de deformaciones, cuyas componentes están en función de los

gradientes del desplazamiento:

[ ]
∂u
∂x (
1 ∂u ∂ v
+
2 ∂ y ∂x ) 12 ( ∂∂ uz + ∂∂wx )
[ ε ]= 1
2 ( ∂v ∂u
+
∂x ∂ y ) ∂v
∂y 2(∂ z ∂ y )
1 ∂v ∂w
+

1
2 ( ∂ w ∂u
+
∂x ∂ z ) ( 1 ∂ w ∂v
+
2 ∂ y ∂z ) ∂∂wz
Por lo tanto las componentes de deformación infinitesimal y de desplazamiento
corresponden a un sistema de seis ecuaciones diferenciales parciales de la forma:

∂u ∂v ∂w
ε xx = ε yy= ε zz =
∂x ∂y ∂z

γ xy = ( ∂∂ uy + ∂∂ vx ) γ =( ∂∂uz + ∂∂ wx ) γ =( ∂∂ vz + ∂∂wy )
xz yz

Siendo u , v , w tres funciones continuas y derivables para todo punto material en el dominio
del medio continuo.

Ecuaciones de compatibilidad de deformaciones infinitesimales.

77
Como se indico anteriormente, el campo de las deformaciones infinitesimales es función
del gradiente del campo del desplazamiento, siendo este último un grupo de funciones
continuas, al igual que sus derivadas, que depende de la posición del punto material y
tiempo.

La relación obtenida entre las componentes de deformación infinitesimal y de

desplazamiento corresponde a un sistema de seis ecuaciones diferenciales parciales de la
forma:

∂u ∂v ∂w
ε xx = ε yy= ε zz =
∂x ∂y ∂z

( ∂∂ uy + ∂∂ vx ) γ =( ∂∂uz + ∂∂ wx ) γ =( ∂∂ vz + ∂∂wy )
γ xy = xz yz

Siendo u , v , w tres funciones continuas y derivables para todo punto material en el dominio
del medio continuo.

Si las funciones de desplazamiento u , v , w en términos de la posición ( x , y , z ) son

conocidas, las seis componentes de deformación se obtienen evaluando las derivadas de
u , v , w con respecto a x , y , z. En cambio, cuando se desean calcular las tres componentes
de desplazamiento, siendo conocidas las seis componentes de deformación, se conforma un
sistema de seis ecuaciones con tres incógnitas, para el cual podría o no existir una solución,
es decir, un valor de las funciones de desplazamiento obtenido de un grupo de componentes
de deformación escogidas de forma arbitraria.

Si se deriva dos veces la componente ε xx con respecto a y, al igual que se deriva dos veces
la componente ε yy con respecto a x se tiene que:

∂2 ε xx 2

∂ y2
= ( )
∂2 ∂u ∂ ε yy ∂2 ∂ v
; =
∂ y 2 ∂ x ∂ x2 ∂ x 2 ∂ y ( )
Derivando la componente de deformación γ xy con respecto a x y a y se tiene que:

∂2 γ xy
( ) ( ) ( )
2 2 2
∂ ∂u ∂ v ∂ ∂u ∂ ∂v
= + = 2 + 2
∂x ∂ y ∂x ∂ y ∂ y ∂ x ∂ y ∂ x ∂x ∂ y

78
 2 2
∂ ε xx ∂ ε yy
Sustituyendo 2
y 2
en la expresión anterior se obtiene la primera ecuación de
∂y ∂x
compatibilidad de deformaciones infinitesimales de la forma:

∂2 ε xx ∂2 ε yy ∂2 γ xy
+ =
∂ y2 ∂ x2 ∂ x ∂ y

De la misma manera se obtienen otras dos ecuaciones de compatibilidad relacionadas con

las componentes de deformación γ xz y γ yz indicadas a continuación:

2 2 2
∂ ε xx ∂ ε zz ∂ γ xz
2
+ 2
=
∂z ∂x ∂ x∂ z

2 2 2
∂ ε yy ∂ ε zz ∂ γ yz
2
+ 2
=
∂z ∂y ∂ y∂ z

Ahora, derivando la componente ε xx con respecto a y y a z se tiene que:

∂2 ε xx
=
∂2 ∂u
∂ y∂ z ∂ y ∂ z ∂ x( )
Y la segunda derivada de las componentes de deformación angular γ yz , γ xz y γ xy con respecto
a x, y y z respectivamente son las siguientes:

−∂2 γ yz
∂ x2
= ( +
∂ x2 ∂ z ∂ y
= )
−∂2 ∂ v ∂ w −∂2 ∂ v
− ( )
∂2 ∂ w
( )
∂ x2 ∂ z ∂ x2 ∂ y

∂2 γ xz
=
∂2 ∂u ∂w
+ (∂2 ∂ w
= 2 +
∂2
) ∂u
( )
∂x ∂ y ∂x ∂ y ∂ z ∂ x ∂ x ∂ y ∂ y ∂ z ∂ x ( )
∂2 γ xy
( ) ( ) ( )
2 2 2
∂ ∂u ∂ v ∂ ∂v ∂ ∂u
= + = 2 +
∂x ∂ z ∂ x ∂z ∂ y ∂x ∂x ∂ z ∂ y ∂ z ∂ x

De la suma de las tres ecuaciones anteriores se obtiene:

−∂2 γ yz ∂2 γ xz ∂2 γ xy −∂ 2 ∂ v
( ) ( ) ( ) ( )
2 2 2
∂ ∂w ∂ ∂w ∂ ∂u
+ + = 2 − 2 + 2 +
∂ x2 ∂ x ∂ y ∂x ∂ z ∂ x ∂ z ∂x ∂ y ∂x ∂ y ∂ y ∂z ∂x

79
 ( ) ( )
2 2
+∂ ∂ v ∂ ∂u
+
∂x ∂ z ∂ y∂ z ∂x
2

−∂2 γ yz ∂2 γ xz ∂2 γ xy ∂2 ε xx
( )
2
∂ ∂u
+ + =2 =2
∂ x2 ∂ x ∂ y ∂x ∂ z ∂ y∂z ∂x ∂ y∂z

Por lo tanto se concluye que:

∂2 ε xx −∂ 2 γ yz ∂2 γ xz ∂2 γ xy
2 = + +
∂ y∂z ∂ x2 ∂ x∂ y ∂ x ∂z

De la misma manera se obtienen otras dos ecuaciones de compatibilidad relacionadas con

las componentes de deformación ε yy y ε zz indicadas a continuación:

2 2 2 2
∂ ε yy ∂ γ yz ∂ γ xz ∂ γ xy
2 = − +
∂ x ∂ z ∂ x ∂ y ∂ y2 ∂ y ∂ z

2 2 2 2
∂ ε zz ∂ γ yz ∂ γ xz ∂ γ xy
2 = + −
∂ x ∂ y ∂ x ∂ z ∂ y ∂ z ∂ z2

En resumen las seis ecuaciones de compatibilidad de deformaciones infinitesimales son las

siguientes:

∂2 γ xy ∂2 ε yy ∂2 ε xx
= +
∂ x ∂ y ∂ x2 ∂ y2

∂2 γ xz ∂2 ε zz ∂2 ε xx
= +
∂ x ∂ z ∂ x2 ∂ z 2

∂2 γ yz ∂2 ε zz ∂2 ε yy
= +
∂ y ∂ z ∂ y2 ∂ z2

2 2 2 2
∂ ε xx −∂ γ yz ∂ γ xz ∂ γ xy
2 = + +
∂ y∂z ∂x
2
∂ x∂ y ∂ x ∂z

2 2 2 2
∂ ε yy ∂ γ yz ∂ γ xz ∂ γ xy
2 = − +
∂ x ∂ z ∂ x ∂ y ∂ y2 ∂ y ∂ z

80
 2 2 2 2
∂ ε zz ∂ γ yz ∂ γ xz ∂ γ xy
2 = + −
∂ x ∂ y ∂ x ∂ z ∂ y ∂ z ∂ z2

Por lo tanto se puede concluir en lo siguiente:

Número de incógnitas:

 6 componentes de la matriz [ T ]

 6 componentes de la matriz [ ε ]

 3 componentes del vector desplazamiento ( ⃗ ^ v ^j+w k^ )

U =u i+

Número de ecuaciones:

 3 ecuaciones de equilibrio interno.

 6 ecuaciones de las componentes de la matriz de deformación.

 6 ecuaciones de compatibilidad de deformaciones infinitesimales.

Por lo tanto para resolver el problema elástico de la forma general, es necesario resolver un
sistema de 15 ecuaciones diferenciales en derivadas parciales con 15 incógnitas.

81
 2.16 Formulación del problema elástico en desplazamientos

Ecuaciones de Navier: El problema elástico lineal se puede plantear en función del campo
de los desplazamientos, de acuerdo a las ecuaciones de gobierno y de las condiciones de
borde apropiadas, por lo tanto escribiendo las ecuaciones de Lamé:

σ xx = ( λ+2 G ) ε xx + λ ε yy + λ ε zz =( ε xx +ε yy+ ε zz ) λ +2G ε xx

σ yy=λ ε xx + ( λ+2 G ) ε yy + λ ε zz =( ε xx+ ε yy + ε zz ) λ+ 2G ε yy

σ zz = λ ε xx+ λ ε yy + ( λ+2 G ) ε zz = ( ε xx + ε yy + ε zz ) λ +2G ε zz

τ xy =G γ xy

τ xz =G γ xz

τ yz=G γ yz

Pero:

∂u ∂v ∂w
ε xx = ε = ε =
∂ x yy ∂ y zz ∂ z

γ xy = ( ∂∂ uy + ∂∂ vx ) γ =( ∂∂uz + ∂∂ wx ) γ =( ∂∂ vz + ∂∂wy )
xz yz

Por otra parte:

⃗ ^ v ^j+ w k^
U =u i+

⃗ ∂ ^ ∂ ^ ∂ ^
∇= i+ j+ k
∂x ∂y ∂z

⃗ ⃗= ∂ u + ∂ v + ∂ w
∇ .U
∂x ∂ y ∂z

Entonces:

82
 σ xx =( ε xx+ ε yy + ε zz ) λ+ 2G ε xx= ( ∂∂ ux + ∂∂ vy + ∂∂wz ) λ+2 G ∂∂ux =(⃗∇ .U⃗ ) λ+2 G ∂∂ ux
De las ecuaciones de equilibrio interno:

∂ σ xx ∂ τ xy ∂ τ xz
ρ x+ + + =0
∂ x ∂ y ∂z

∂ τ yx ∂ σ yy ∂ τ yz
ρy+ + + =0
∂x ∂y ∂z

∂ τ zx ∂ τ zy ∂ σ zz
ρz+ + + =0
∂x ∂ y ∂z

Partiendo de la primera ecuación de equilibrio interno:

ρ x+
∂ ⃗ ⃗
∂x { ∂x ∂ y } {(
∂ y ∂x ∂z ∂ z ∂x )} { (
( ∇ . U ) λ+ 2G ∂ u + ∂ G ∂u + ∂ v + ∂ G ∂ u + ∂ w =0 )}
2 2 2 2 2
∂ ⃗ ⃗
ρ x+ λ ( ∇ . U ) +2 G ∂ u2 +G ∂ u2 +G ∂ v +G ∂ u2 +G ∂ w =0
∂x ∂x ∂y ∂x ∂ y ∂z ∂x∂y

( ) ( )
2 2 2
∂ (⃗ ⃗) ∂ ∂u ∂ v ∂ w ∂ u ∂ u ∂u
ρ x+ λ ∇ . U +G + + +G + + =0
∂x ∂x ∂x ∂y ∂z ∂ x2 ∂ y 2 ∂ z 2

Por otra parte:

2 2 2
2 ∂ ^ ∂ ^ ∂ ^
∇= 2
i+ 2 j+ 2 k
∂x ∂y ∂z

2 2 2
∂u ∂ u ∂ u
∇ 2 . u= 2
+ 2+ 2
∂x ∂ y ∂ z

Entonces:

∂ ⃗ ⃗
ρ x+ λ ( ∇ . U ) +G ∂ ( ⃗ ⃗ ) +G ( ∇2 . u )=0
∇ .U
∂x ∂x

∂ ⃗ ⃗
ρ x + ( λ+G ) ( ∇ . U ) +G ( ∇ 2 . u )=0
∂x

83
De la misma manera se obtienen las otras dos ecuaciones relacionando las ecuaciones de
equilibrio con las deformaciones infinitesimales indicadas a continuación:

∂ ⃗ ⃗
ρ y + ( λ +G ) ( ∇ . U ) +G ( ∇ 2 . u )=0
∂y

∂ ⃗ ⃗
ρ z + ( λ+G ) ( ∇ . U ) +G ( ∇ 2 .u ) =0
∂z

Si multiplicamos estas tres ecuaciones por los vectores unitarios i^ , ^j, k^ respectivamente y
sumamos miembro a miembro, tenemos:

∂ ⃗ ⃗ ^
ρ x i^ + ρ y ^j + ρ z k^ + ( λ+G ) ( ∇ . U ) i + ( λ+G ) ∂ ( ⃗ U ) ^j + ( λ +G )
∇.⃗
∂ ⃗ ⃗ ^
(∇ .U ) k
∂x ∂y ∂z

^
+G ( ∇ . u ) i+G
2
( ∇2 . u ) ^j+G ( ∇ 2 . u ) k^ =0

Por otra parte:

^ ρ y ^j+ ρz k^
⃗ρ =ρx i+

⃗
0=0 ^ 0 ^j+0 k^
i+

∂ ⃗ ⃗ ^ ∂ ⃗ ⃗ ^ ∂ ⃗ ⃗ ^
∇ .( ⃗ U )=
∇.⃗ ( ∇ . U ) i+ ( ∇ . U ) j + ( ∇ . U ) k
∂x ∂y ∂z

∇2 [ U
⃗ ] =∇ 2 ( u i^ +v ^j+ w k^ )

Finalmente:

( ρx i^ + ρ y ^j+ ρz k^ ) + ( λ+G ) {∂∂x (⃗∇ .U⃗ )i^ + ∂∂y (⃗∇ .U⃗ ) ^j+ ∂∂z ( ⃗∇ . ⃗U ) k^ }
+G {( ∇ .u ) i^ + ( ∇ .u ) ^j + ( ∇ . u ) k^ }=0
2 2 2

⃗ρ + ( λ +G ) ⃗
∇ . (⃗ U ) +G ∇ 2 ( U
∇.⃗ ⃗ )=0⃗

Ecuación fundamental de la elasticidad

84
La ecuación fundamental de la elasticidad es la única a la que debe satisfacer el vector
desplazamiento, en todos los puntos del sólido elástico.

2.17 Formulación del problema elástico en tensiones

Otra forma de resolver el problema elástico lineal estático consiste en plantear un grupo de
ecuaciones diferenciales parciales en función del estado de esfuerzos.

Ecuaciones de Michell: Estas ecuaciones resultan del procedimiento de sustituir la

ecuación constitutiva en la ecuación de compatibilidad de deformación infinitesimal para
luego reemplazar en la ecuación de equilibrio, lo cual resulta:

∇ 2 σ xx +
1 ∂2
1+ν ∂ x 2 ( xx
σ +σ yy +σ zz )=−
ν ⃗
1−ν ( )
∇ . ⃗ρ−2
∂δx
∂x

2
∇ σ yy +
1 ∂2
1+ ν ∂ y 2
( σ xx +σ yy +σ zz )=− ( )
ν ⃗
1−ν
∇ . ⃗ρ−2
∂ δy
∂y

2
∇ σ zz +
1 ∂2
( σ +σ + σ )=−
1+ν ∂ z 2 xx yy zz
ν ⃗
1−ν
∇ . ⃗ρ−2( )
∂ δz
∂z

2
∇ τ xy +
1 ∂2
1+ν ∂ x ∂ y
( σ xx+ σ yy + σ zz) =− [
∂ δx ∂ δy
+
∂ y ∂x ]
[ ]
2
1 ∂ ∂ δx ∂ δz
∇ 2 τ xz + ( σ xx +σ yy + σ zz )=− +
1+ ν ∂ x ∂ z ∂ z ∂x

2
∇ τ yz +
1 ∂2
1+ν ∂ y ∂ z
( σ xx +σ yy + σ zz )=− [
∂ δy ∂ δz
+
∂z ∂ y ]
Si ρ es constante estas ecuaciones se llaman Ecuaciones de Beltrami las cuales se pueden
escribir como:

2
( 1+ν ) ∇ 2 σ xx + ∂ 2 ( σ xx +σ yy+ σ zz) =0
∂x

85
 2
( 1+ν ) ∇ 2 σ yy + ∂ 2 ( σ xx +σ yy+ σ zz )=0
∂y

2
( 1+ν ) ∇ 2 σ zz + ∂ 2 ( σ xx + σ yy + σ zz) =0
∂z

2
( 1+ν ) ∇ 2 τ xy + ∂
( σ + σ + σ ) =0
∂ x ∂ y xx yy zz

2
( 1+ν ) ∇ 2 τ xz + ∂ ( σ xx + σ yy + σ zz ) =0
∂ x ∂z

2
( 1+ν ) ∇ 2 τ yz + ∂ ( σ xx+ σ yy + σ zz ) =0
∂ y ∂z

86
 2.18 Método de los elementos finitos aplicado a la teoría de la elasticidad
partiendo de la teoría de los trabajos virtuales

El principio de los trabajos virtuales es el primer principio variacional aplicable a la

mecánica del medio continuo y el análisis estructural.

En cuerpos deformables, el trabajo virtual hecho por las fuerzas reales se puede dividir en
dos partes: el trabajo virtual hecho por las fuerzas internas δ W i denominado trabajo virtual
interno y el trabajo virtual realizado por las fuerzas externas δ W e llamado trabajo virtual
externo.

El trabajo virtual externo, es decir, aquel realizado por las fuerzas reales externas mientras
se presenta un desplazamiento virtual δU es igual a:

δ W e =∫ { δU } { P } dS+∫ { δU } { b } dV + ∑ {δU ( x ) } {f ( x ) }
T T T
n n
Γ V

El trabajo virtual interno por unidad de volumen es igual al trabajo realizado por el esfuerzo
T
asociado a desplazamientos virtuales por unidad de volumen, es decir [ [ ∇ ] { δU } ] { σ }. Por lo
tanto el trabajo virtual interno almacenado en el cuerpo está definido como:

δ W i=∫ [ [ ∇ ] { δU } ] { σ } dV
T

El principio de los trabajos virtuales establece que un cuerpo está en equilibrio si y solo si,
el trabajo virtual realizado por todas las acciones internas y externas es nulo, es decir:

δW =−δ W e +δ W i=0

Por lo tanto se concluye que:

δ W i=δ W e

87
 ∫ [ [ ∇ ] { δU } ] {σ } dV =∫ { δU } { P } dS+∫ { δU } {b } dV + ∑ { δU ( x ) } { f ( x ) }
T T T T
n n
V Γ V

Donde:

[ ]
∂/∂ x 0 0
0 ∂/∂ y 0
[ ∇ ]= 0 0 ∂ /∂ z
∂ /∂ y ∂/∂ x 0
∂ /∂ z 0 ∂/∂ x
0 ∂/∂ z ∂ /∂ y

{}
δu
{ δU }= δv → Vector de desplazamientos virtuales
δw

{}
σ xx
σ yy
σ
{ σ }= zz →Vector de esfuerzos
σ xy
σ xz
σ yz

{}
Px
{ P }= P y → Vector de fuerzas de superficie
Pz

{}
bx
{ b }= b y → Vector de fuerzas de cuerpo
bz

{}
fx
{ ( x ) } f y →Vector de fuerzas puntuales (escalares)
f n
=
fz

El vector de deformaciones puede escribirse como:

88
 { ε }= [ ∇ ] { U }

[ ∇ ] { δ U (e) }=[ B(e) ] { δ a(e) }

Elevando a ambos miembros por la transpuesta a esta última expresión, se tiene que:

T
( [ ∇ ] {δ U (e) }) =( [ B ] {δ a })
T (e) (e)

T T
( [ ∇ ] {δ U (e) }) = {δ a } [ B ]
T (e) (e)

Expresando el vector de esfuerzos en función a la matriz constitutiva y al vector de

deformación:

{ σ(e) }=[ D(e) ] {ε (e) }

{ σ(e) }=[ D(e) ] ( [ B(e) ] {a(e) })

La suma de los trabajos virtuales de las cargas puntuales es igual a un formato matricial
para toda la estructura igual a:

∑ {δU ( x ) }T {f ( x ) }=¿ { δa}T {f (n) } ¿

n n

{ δa } { f (n) } → en el tamaño de toda la estructura

T

{ δU (e) }=[ N (e) ] {δa (e) }

Elevando a ambos miembros por la transpuesta a esta última expresión, se tiene que:

T T
( {δU (e) }) =( [ N (e) ] {δa(e) })
T T T
( {δU (e) }) = {δa(e) } [ N (e) ]
n n n
T T T
∑ ∫ ([ ∇ ] {δU (e) }) {σ (e) } dV =∑ ∫ {δU (e) } { P(e) }dS +∑ ∫ {δU (e) } {b(e) } dV
e=1 V ( e) e=1 Γ (e) e=1 V (e)

+ { δa } { f }
T (n )

89
 n n
T T T T
∑ ∫ { δa(e) } [ B(e) ] [ D(e) ][ B(e) ] {a(e) } dV =∑ ∫ {δa(e) } [ N (e) ] {P(e) } dS
e=1 V ( e) e=1 Γ (e)

n
T T
+ ∑ ∫ { δa(e) } [ N (e) ] { b(e) } dV + { δa } { f (n) }
T

e=1 V ( e)

Como δa (e) son valores virtuales que no dependen de dV y el vector de desplazamiento a(e)
son valores fijos salen fuera del signo de la integral, entonces:

[ ] [ ]
n n
T T T T
∑ {δa(e) } {a(e) } ∫ [ B(e) ] [ D(e) ][ B(e) ] dV =∑ {δa(e) } ∫ [ N (e) ] { P(e) } dS
e=1 V
(e) e=1 Γ
( e)

[ ]
n
T T
+ ∑ { δa(e) } ∫ [ N (e) ] {b(e) } dV + { δa } { f (n) }
T

e=1 V
(e)

Recordar que:

T
[ K (e) ]=∫ [ B(e) ] [ D(e) ][ B(e) ] dV
(e)
V

Entonces:

[ ]
n n

∑ [ {δa(e) } [ K (e)] {a(e) }]=∑ {δa(e) } ∫ [ N (e) ] { P(e) } dS

T T T

e=1 e=1 Γ
(e)

[ ]
n
T T
+ ∑ { δa(e) } ∫ [ N (e) ] {b(e) } dV + { δa } { f (n) }
T

e=1 V
(e)

Como:

(e) T
s }= ∫ [ N ] { P } dS
{ f (e) (e)
(e)
Γ

T
{ f (e)b }=∫ [ N (e)] {b(e) } dV
(e)
V

Por lo tanto resulta:

90
 n n n

∑ [ {δa } [ K ] {a }]=∑ [ {δa } {f }] +∑ [ {δa } { f }]+¿ { δa }

(e) T (e) (e) (e) T (e)
s
(e) T (e)
b
T
{ f (n ) }¿
e=1 e=1 e=1

n n

∑ [ {δa (e) } [ K (e) ] {a(e) }]=∑ [ {δa (e) } ( { f (e)s }+{ f (e)b }) ] +{ δa }T { f (n) }
T T

e=1 e=1

Como la suma de las fuerzas de superficie y de volumen puede expresarse como la fuerza
total en el elemento, resulta:

{ f (e) }={ f (e)

s }+ {f b }
(e)

Luego:

n n

∑ [ {δa } [ K ] {a }]=∑ {δa

(e) T (e) (e) (e) T
} { f (e) }+ { δa }T {f (n) }
e=1 e=1

n n

∑ [ {δa } [ K ] {a }]−∑ {δa

(e) T (e) (e) (e) T
} { f (e) }+ { δa }T {f (n) }= {0 }
e=1 e=1

Recordar que:

n
T
∑ {δa(e) } ={ δa }T
e=1

∑ [ K (e)]=[ K ]
e=1

∑ {a(e) }={ a }
e=1

Como la suma representa el ensamblaje de toda la estructura, la ecuación queda:

{ δa } [ K ] { a }− {δa } { f (e) } + { δa } { f (n) }= { 0 }

T T T

{ δa } [ K ] { a }− {δa } ( { f (e) } + { f (n) } ) = {0 }

T T

91
Como:

n
{ f }=∑ { f (e) }+ ¿ { f (n) } → Fuerza resultante en toda laestructura ¿
e=1

Resulta:

T T
{ δa } [ K ] { a }− {δa } { f }= {0 }

T
{ δa } ( [ K ] { a }− { f } )= { 0 }

Como { δa }T son desplazamientos virtuales independientes no pueden ser cero es decir

T
{ δa } ≠ 0 por lo tanto la segunda condición es:

[ K ] { a }− { f }= {0 }

[ K ] { a }= { f }

Se puede concluir que existe una relación entre la matriz de rigidez [ K ] , el vector de
desplazamientos { a } y el vector de fuerzas { f }, en la cual son conocidos la matriz de rigidez
y el vector de fuerzas lo que permite encontrar el vector de desplazamientos.

92
 3. DECRIPCIÓN DEL PROGRAMA

3.1 Características técnicas del programa

Las principales características técnicas del programa de elementos finitos a código abierto
son las siguientes:
 Funcionalidad: El programa tiene más de 100 instrucciones escritas en lenguaje
Visual Basic, que permiten resolver problemas mediante el método de los elementos
finitos.
 Estabilidad= Las subrutinas del programa cuentan con controles y trampas de
errores que garantizan la estabilidad del programa durante su ejecución. Sin
embargo, la posibilidad de modificar las subrutinas de los programas a código
abierto como este, reducen su robustez.
 Usabilidad: La aplicación Microsoft Excel proporciona una interfaz de entrada y
salida de datos ágil y amigable para el usuario.
 Eficiencia: Debido a sus fines didácticos, el programa guarda en la memoria RAM
algunas matrices características de las etapas intermedias del procedimiento, sin
embargo los usuarios avanzados podrán evitar que dichas matrices sean guardadas.
El tiempo de ejecución depende del número de grados de libertad que tenga el
problema y de la cantidad de resultados intermedios que se desean observar.
 Extendibilidad: Al ser un programa a código abierto es fácil adicionar nuevas
subrutinas que permitan la solución de otros problemas mediante el método de los
elementos finitos.
 Portabilidad: Tanto el programa como los datos de un problema están contenidos
en un libro de Microsoft Excel, esto facilita su portabilidad a computadores que
cuenten con un sistema operativo Microsoft Windows y con el programa Microsoft
Excel, lo cual es común en el medio.

93