Cilindros Pared Delgada

UNIVERSIDAD MAYOR REAL Y PONTIFICIA DE SAN FRANCISCO XAVIER DE CHUQUISACA
MECANICA ESTRUCTURAL AVANZADA

ING. MOISES RTEAGA MIRANDA TRABAJO DE INVESTIGACIN
APLICACIN DEL MTODO DE ELEMENTO FINITO EN EL

CALCULO DE CILINDROS DE PARED DELGADA
1. INTRODUCCION.
El Mtodo de Elemento Finito (MEF en castellano o FEM en ingls) es un mtodo

numrico general para la aproximacin de soluciones de ecuaciones diferenciales
parciales muy utilizado en diversos problemas de ingeniera y fsica. Este mtodo est
pensado para ser usado en computadoras y permite resolver ecuaciones diferenciales
asociadas a un problema fsico sobre geometras complicadas. La variedad de
problemas a los que puede aplicarse ha crecido enormemente, siendo el requisito
bsico que las ecuaciones constitutivas y ecuaciones de evolucin temporal de
problema a considerar sean conocidas de antemano. Las aplicaciones van desde el
anlisis por deformacin y esfuerzos de estructuras hasta el anlisis en campos como
flujo de calor, de fluidos, magntico, filtraciones y otros.
Es as que en este trabajo de investigacin se usara este mtodo aplicado al clculo de
cilindros de pared delgada; y estos a su vez pueden ser calderas, depsitos (silos,
bunkers, etc.), tuberas, cpulas, etc. y para un mejor desarrollo de la investigacin se
indicaran adems los esfuerzos que trabajan sobre la pared delgada de un cilindro
como tambin se explicara como trabaja este mtodo de elemento finito, de qu
manera se proceden a obtener nodos y posteriormente mallas para as proceder al
clculo.
2. DESARROLLO DEL TEMA.
2.1. Mtodo de Elemento Finito.
El Mtodo del Elemento Finto (MEF), es una tcnica de anlisis matemtico,

especficamente matricial, que se emplea para obtener soluciones aproximadas a
problemas reales. Su aplicacin abarca diferentes campos de la ingeniera, tales como
los que se muestran en la Tabla
Tabla 1. Aplicaciones comunes en diferentes campos de ingeniera
MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

Debido a que al resolver un ejemplo aplicando este mtodo se forman muchos nodos y
a la vez a mallas es que en el mercado existe una gran cantidad de software basado en
el Mtodo de Elemento Finito, para una gran diversidad de aplicaciones. En la
siguiente tabla se muestran algunos de ellos.
- TIPOS DE ELEMENTOS FINITOS.
Los tipos de elementos finitos se pueden clasificar de la siguiente manera:
Elemento Finito Puntual.- Es un elemento que no tiene dimensin y su

definicin se reduce a un solo punto o nodo. Su aplicacin es, muy comnmente
para representar masas concentradas o representar resortes con rigidez donde
su longitud no es una caracterstica relevante.
Elemento Finito Lineal.- Es un elemento que tiene una sola dimensin,
generalmente es definido por dos nodos y su longitud si es relevante. Sin
embargo, aunque el elemento es lineal, tambin es posible definirle
propiedades para su seccin transversal
(Fig. 1) Elemento Finito Lineal
Elemento Finito Superficial Placa.- Los elementos planos, definidos por tres o
cuatro nodos (triangulares y cuadrilteros respectivamente), son las placas
delgadas que son cargadas en su propio plano (esfuerzo plano).

(Fig. 2)Elementos Finitos Superficiales Placa
Elementos Finitos Volumtricos Slidos.- Son la generalizacin tridimensional de

los elementos en esfuerzo plano: el tetraedro, el pentaedro, y el hexaedro son
las formas ms comunes de los elementos tridimensionales (Figura 5.3).
(Fig. 3) Elementos Finitos Volumtricos Solidos
Elementos Finitos Volumtricos Axisimtricos.- Utilizados para la resolucin de

una gran variedad de problemas de ingeniera incluyendo tanques y
contendores de acero y de concreto, recipientes de contenido nuclear, rotores,
mbolos y escapes de cohetes; en todos ellos se utilizan estos elementos. Su
caracterstica principal es que tanto la carga como la geometra, usualmente
son axisimetricos.
Placa Plana en Flexin.- Son elementos usados no solo entre s, si no tambin
junto con cascarones y miembros de pared delgada. Las formas geomtricas
son anlogas a las de los elementos en esfuerzo plano, con mayor nfasis
tambin en las formas triangulares, y cuadrilteras.
Cascaron Axisimtrico.- Tienen la misma importancia en aplicaciones prcticas
que los slidos axisimetricos, aunque aqu, las formulaciones se derivan de la
teora de la membrana. Dentro de esta formulacin se encuentra la diferencia
con respecto a los elementos placa en flexin y tensin y sirven para identificar
problemas clave.
Cascarn Curvo.- Cuando una estructura est curva, es preferible usar
elementos cascarn curvo para los modelos analticos. Dentro de las ventajas
esta la habilidad para describir de forma ms adecuada la geometra de una
superficie curva. Existe un gran nmero de alternativas para formular este tipo
de elementos.

2.2. Esfuerzo en Cilindros.
Con frecuencia los cilindros en sus diferentes usos pueden ser sujetos a presiones
internas y/o externas. Algunas aplicaciones comunes son cilindros de aire o
hidrulicos, depsitos de almacenamiento de fluidos y tuberas, as como caones
para armas. Algunos de estos dispositivos tienen uno o los dos extremos abiertos y
algunos estn cerrados. Si tienen un extremo abierto, aparecer un estado de
esfuerzos radiales r y tangenciales t (de anillo). Si estn cerrados, tambin estar
presente un esfuerzo en la tercera dimensin llamado longitudinal o axial z. Estos
tres esfuerzos aplicados son mutuamente ortogonales y son principales, ya que no hay
cortante debido a la presin uniformemente distribuida.
Se pueden estudiar cilindros de pared gruesa, cuyas ecuaciones son exactas y vlidas
para cualquier caso, y hay una teora ms especfica para cilindros de pared delgada,
en la que se usan ecuaciones aproximadas y simplificadas, vlidas a partir de un
determinado espesor.
2.3. Cilindros de pared delgada.
Los recipientes de pared delgada constituyen una aplicacin importante del anlisis
de esfuerzo plano. Como sus paredes oponen poca resistencia a la flexin, puede
suponerse que las fuerzas internas ejercidas sobre una parte de la pared son
tangentes a la superficie del recipiente.
Considerando recipiente cilndrico de radio interior r y espesor de pared t, que

contiene un fluido a presin Se van a determinar los esfuerzos ejercidos sobre un
pequeo elemento de pared con lados respectivamente paralelos y perpendiculares al
eje del cilindro. Debido a la simetra axial del recipiente y de su contenido, no se
ejercen esfuerzos cortantes sobre el elemento.
El anlisis de pared delgada se refiere a los recipientes con la siguiente relacin:

3. ANALISIS MATEMATICO DEL TEMA.
El Mtodo de Elemento Finito viene dado por ecuaciones fundamentales y a partir de

ese punto se adecua a diferentes casos; a continuacin se ver estas ecuaciones:
- Funciones de los desplazamientos. En la figura 4 se muestra un elemento

triangular tpico considerado, con los nodos i, j, m numerados en sentido anti
horario.
(Fig. 4)
Los desplazamientos de un nodo tienen dos componentes:
Ec.1
Y las seis componentes de los desplazamientos del elemento se agrupan en un
vector
Ec.2
Los desplazamientos interiores a un elemento han de quedar definidos

unvocamente por esos seis valores. La representacin ms sencilla viene dada
evidentemente por dos polinomios de primer grado
Ec.3
Se pueden calcular fcilmente las seis constantes resolviendo los dos

sistemas de tres ecuaciones simultneas que se obtienen al sustituir las

coordenadas de los nodos e igualar las expresiones resultantes a los

desplazamientos correspondientes a los nodos: escribiendo, por ejemplo,
Ec.4
Podemos calcular fcilmente 1, 2 y 3 en funcin de los desplazamientos

nodales ui, uj y um para obtener finalmente
Ec.5
En la cual
Ec.6
Obtenindose los otros coeficientes mediante permutacin cclica de los

subndices i, j, m, y donde
Ec.7
Puesto que las ecuaciones para el desplazamiento vertical don similares, se

obtiene igualmente que
Ec.8
Aunque no sea estrictamente necesario en este momento, podemos expresar

las relaciones anteriores Ec.5 y Ec.7 en la siguiente forma general
Ec.9

Siendo I una matriz unidad de 2x2, y
Ec.10
La funcin de desplazamientos elegida garantiza automticamente la

continuidad de desplazamientos entre elementos adyacentes, debido a que los
desplazamientos varan linealmente a lo largo de cualquiera de los lados del
tringulo y, al imponer los mismos desplazamientos en los nodos para dos
elementos contiguos, evidentemente existir el mismo desplazamiento a lo
largo del contorno de separacin.
- Deformacin (total). La deformacin total en cualquier punto del elemento

puede definirse mediante sus tres componentes que contribuyen al trabajo
interno. As,
Ec.11
Sustituyendo la Ec.9 tenemos
Ec.12
Siendo Bi una matriz tpica, dada por
Ec.13
Con lo cual, la matriz B queda definida explcitamente. Se advertir que en este

caso la matriz B es independiente de la posicin del punto dentro del elemento
y, por consiguiente, las deformaciones son constantes en todo el mismo.

- Matriz de elasticidad. La matriz D de la relacin
Ec.14
Se puede establecer explcitamente para cualquier material (excluyendo ahora
o, que es simplemente un trmino aditivo).
Tensin plana de un material istropo. Para los estudios de tensin plana en

los materiales istropos se tiene, por definicin,
Ec.15
Despejando las tensiones en el sistema anterior, obtenemos la matriz D
Ec.16
En el cual E es el mdulo de elasticidad y el coeficiente de Poisson.
Deformacin plana de un material istropo. En este caso la tensin normal z

no es nula, debiendo aadirse a las otras tres componentes de la tensin, Para
el caso particular de dilatacin trmica istropa, tenemos
Ec.17
Pero, adems

Ec.18
Eliminando z y despejando las tres tensiones restantes de la expresin ya

citada de la deformacin inicial y comparando esta con la Ec.14, la matriz D
resulta
Ec.19
- Matriz de Rigidez. La matriz de rigidez de un elemento ijm viene definida a

partir de la siguiente relacin
Ec.20
En la que t es el espesor del elemento donde la integracin se efecta sobre la
superficie del tringulo. Si el espesor del elemento se supone constante,
hiptesis que tiende a ser tanto ms cierta a medida que disminuye el tamao
de los elementos, entonces, como ninguna de las matrices contiene a x o y,
tendremos sencillamente
Ec.21
Donde es el rea del tringulo (definido en la Ec.5). Esta expresin es ya
suficientemente explicita para el clculo dejando para el computador las
operaciones matriciales.
En nuestro trabajo de investigacin se analizara el Mtodo de Elementos Finitos

aplicado a un cilindro de pared delgada es por eso que entonces y a partir de lo
anterior se tendr lo siguiente:
Y cambiara nuestro campo de desplazamientos por lo siguiente:

Ec.22
Donde r y z representan respectivamente las coordenadas radial y axial de un punto y

u, v son los correspondientes desplazamientos radial y axial.
- Funciones de Interpolacin. Al ser la deformacin plana, el campo de

deformaciones en el interior del elemento se aproxima mediante la expresin
habitual:
...Ec.23
Siendo Ui las deformaciones en direccin radial y Vi las deformaciones en

direccin axial, que forman el vector de deformaciones nodales del elemento:
Ec.24
donde n es el nmero de nudos del elemento. La expresin tambien puede

ponerse en forma matricial:
Ec.25
La matriz de funciones de interpolacin N tiene dos filas y tantas columnas

como grados de libertad haya entre todos los nudos del elemento. La estructura
de esta matriz es la misma que en un problema de elasticidad plana:
Ec.26
- Deformaciones Unitarias. La deformacin unitaria en un problema con simetra
de revolucin tiene cuatro componentes, cuya expresin en funcin de los dos
desplazamientos del punto es:

Ec.27
En esta expresin es la deformacin unitaria circunferencial, producida por la

deformacin radial, al ser el dominio de revolucin. La expresin se puede
poner en la forma:
Ec.28
En esta expresin se identifica al operador matricial que pasa de las

deformaciones u a las deformaciones unitarias.
Deformaciones unitarias

Sustituyendo las deformaciones u en funcin de las deformaciones nodales, a

travs de las funciones de interpolacin se obtiene:
Ec.29
La expresin de la matriz B es:
Ec.30
Esta matriz se puede poner en la forma:
Ec.31
siendo cada una de las matrices:
Ec.32
- Estado de Tensiones. Ecuacin constitutiva. El estado de tensiones
correspondiente a un problema con simetra de revolucin, en coordenadas
cilndricas, tiene cuatro componentes:
Ec.33
La ecuacin constitutiva, en ausencia de temperaturas, es:

Ec.34
donde se identifica la matriz elstica D.
- Formulacin Isoparametrica. La interpolacin de coordenadas se emplea de la

misma forma que en los problemas de elasticidad bidimensional, pero aplicada a la
definicin de las coordenadas radial y axial:
Ec.35
El vector xe agrupa a las coordenadas (r,z) de todos los nudos del elemento.
Interpolacin de coordenadas
De la misma forma, existen elementos de tipo sub y superparamtrico, en

funcin del nmero de nudos usados para definir la forma del elemento.

Pueden emplearse asimismo elementos triangulares, tal y como se indica en el

estudio de la elasticidad bidimensional.
El clculo de las propiedades del elemento requiere conocer las derivadas de
las funciones de interpolacin respecto a las coordenadas generales. Siguiendo
las reglas de la derivacin se puede poner la siguiente relacin matricial entre
las derivadas en ambos sistemas de coordenadas:
La matriz J es la matriz Jacobiana de la transformacin de coordenadas (r,z) a

(). Despejando se obtiene:
Ec.36
El vector de la derecha se puede calcular sin ms que derivar las Ni respecto a ,.

Conociendo J se pueden obtener todas las derivadas que forman la matriz Bi. El
clculo de J se hace apoyndose en la interpolacin de coordenadas:
Ec.37
Esta expresin puede ser evaluada fcilmente, ya que las funciones N son
conocidas en funcin de y ri,zi son las coordenadas de los nudos que
definen la forma del elemento.
- Matriz de Rigidez. La matriz de rigidez K del elemento es:
Ec.38
El dominio de integracin expresado en las coordenadas locales es:

Ec.39
El valor del radio se aproxima segn la ley de interpolacin de coordenadas,
con lo que el diferencial de volumen queda:
Ec.40
Considerando la estructura de B, la matriz de rigidez K es:
Ec.41
La matriz K se puede dividir en nxn submatrices, que relacionan a los n nudos entre s:
Ec.42
Cada una de ellas tiene la expresin:
Ec.43
Sobre la naturaleza de estas integrales, pueden hacerse las mismas consideraciones

que en el caso de elasticidad bidimensional. nicamente debe hacerse notar que el
orden de los polinomios es mayor que en aquel caso, debido a la presencia del trmino
N i r i en el integrando.
Ya conociendo de manera general las ecuaciones que dominan este Mtodo de

Elemento Finito a continuacin se obtendrn las ecuaciones ms generales para los
cilindros sometidos a presin. A partir de ese resultado general, la solucin se
simplificar dependiendo las distintas hiptesis, como puede ser el caso de los
cilindros de pared delgada.

(Fig. 4) Esquema General de un Cilindro con Presin Interna y Externa
Aplicando la ley de Hook tenemos la siguiente ecuacin para la deformacin a lo largo

del cilindro (z), que, como dijimos en la introduccin, consideraremos constante, por
lo que tenemos:
A continuacin, se estudia un elemento diferencial del cilindro (Fig. 4) y, a partir del

equilibrio del mismo, se determinan las diferentes tensiones a las que est sometido
(Fig. 5). Se inicia as el desarrollo de la demostracin hasta llegar a las ecuaciones de
Lam.
(Fig. 5) Elemento Diferencial dentro de un Cilindro

(Fig. 6) Detalle de las Tensiones de un Elemento Estructural
El elemento diferencial tendr un volumen infinitesimal que viene dado por la

siguiente expresin diferencial:
Como el sumatorio de fuerzas tiene que ser cero (elemento en equilibrio) se tiene lo
siguiente:
Integrando la ltima ecuacin obtenemos:
A continuacin, se introducen las condiciones y se determinan los valores de c1 y c2:

De donde:
Conociendo estos dos valores podemos deducir expresiones para los esfuerzos
tangenciales y radiales en funcin de los radios y las presiones (interna y externa)
para el caso ms general posible, a estas ecuaciones se las conoce como ecuaciones de
Lam:
Con estas dos ecuaciones se pueden obtener ambas tensiones en cualquier cilindro,
independientemente de las presiones y del espesor, son vlidas para cualquier caso, A
partir de aqu analizaremos los casos ms particulares suponiendo varias hiptesis.
Antes hay que sealar que estas tensiones son principales. La distribucin de ambas
se muestra en la figura 2.6 cuando no hay presin externa.
(Fig. 7) Esfuerzos principales

(Fig. 8) Distribucin de las tensiones tangenciales y radiales
Por otro lado, tenemos que hablar de los esfuerzos longitudinales (eje Z) cuando el
cilindro est cerrado en sus extremos y po= 0. Son constantes (Fig.9), y vienen dados
por la siguiente expresin:
(Fig. 9) Esfuerzos longitudinales

Cilindros de Pared Delgada.

En caso de que los extremos del cilindro estn cerrados, el esfuerzo axial sobre los
mismos es:
Se puede observar en esta ecuacin la ausencia de r, ya que el esfuerzo axial es

uniforme en todo el espesor de la pared de cierre.
Si la presin externa es nula, p0 = 0, entonces las ecuaciones se reducen a:
Y si estn cerrado en sus extremos:
Otro caso que se puede suponer, es un cilndrico macizo sometido a una presin p0
sobre su superficie lateral y a una presin q sobre sus bases, con esta hiptesis
aparece un estado de esfuerzos uniforme en el que sobre su superficie lateral o
cualquier superfie paralela de radio r actuarn los esfuerzos = -Po, y sobre las
superficies paralelas a las bases los esfuerzos = -q.
Haciendo en las ecuaciones de Lam:
Cilindro macizo: ri = 0
Sin presin interna: Pi = 0
Resulta:
t = -Po r = -Po
Y, adems, del equilibrio sobre cualquier seccin paralela a las bases:
Se considera que un cilindro es de pared delgada cuando tiene un espesor t pequeo

tal que:
El esfuerzo radial producido por la presin del contenido del cilindro es muy pequeo
en comparacin con el esfuerzo tangencial.
Suponemos Po = 0 (no hay presin exterior):

Valor medio:
Esfuerzo radial despreciable en comparacin con el esfuerzo tangencial. Despus de

todas estas consideraciones se trabaja, cuando el cilindro es de pared delgada, con las
siguientes ecuaciones:
- En todo el espesor:
- Cilindro cerrado en los extremos:

- Deformaciones:
4. IDEALIZACIONES DE CLCULO.
Para entender los fundamentos bsicos de este mtodo se har un breve desarrollo de
un algoritmo de elementos finitos para resolver un problema definido mediante
ecuaciones diferenciales y condiciones de contorno requiere en general cuatro etapas:
- El problema debe reformularse en forma variacional.
- El dominio de variables independientes (usualmente un dominio espacial para

problemas dependientes del tiempo) debe dividirse mediante una particin en
subdominios o elementos finitos. Asociadas a la particin anterior se construye
un espacio vectorial de dimensin finita, llamado espacio de elementos finitos.
Siendo la solucin numrica aproximada obtenida por elementos finitos una
combinacin lineal en dicho espacio vectorial.
- Se obtiene la proyeccin del problema variacional original sobre el espacio de

elementos finitos obtenido de la particin. Esto da lugar a un sistema con un
nmero de ecuaciones finito, aunque en general con un nmero elevado de
ecuaciones incgnitas. El nmero de incgnitas ser igual a la dimensin del
espacio vectorial de elementos finitos obtenido y, en general, cuanto mayor
ser dicha dimensin tanto mejor ser la aproximacin numrica obtenida.
- El ltimo paso es el clculo numrico de la solucin del sistema de ecuaciones.
El Mtodo de Elemento Finito permite obtener una solucin numrica aproximada

sobre un cuerpo, estructura o dominio (medio continuo) sobre el que estn definidas
ciertas ecuaciones diferenciales en forma dbil o integral que caracterizan el
comportamiento fsico del problema, dividindolo en un nmero elevado de
subdominios, no intersectantes entre s, denominados elementos finitos. El conjunto
de elementos finitos forma una particin del dominio tambin denominada
discretizacin. Dentro de cada elemento se distinguen una serie de puntos
representativos llamados nodos. Dos nodos son adyacentes si pertenecen al mismo

elemento finito, adems, un nodo sobre la frontera de un elemento finito puede

pertenecer a varios elementos. El conjunto de nodos considerando sus relaciones de
adyacencia se llama malla.
(Fig. 10) Discretizacion de nodos y mallas en la superficie
Los clculos se realizan sobre una malla de puntos o nodos, que sirven a su vez de
base para la discretizacin del dominio en elementos finitos. La generacin de la malla
se realiza usualmente con programas especiales llamados generadores de mallas, en
una etapa previa a los clculos que se denomina pre-proceso. De acuerdo con estas
relaciones de adyacencia o conectividad se relaciona el valor de un conjunto de
variables incgnitas definidas en cada nodo y denominadas grados de libertad. El
conjunto de relaciones entre el valor de una determinada variable entre los nodos se
puede escribir en forma de sistema de ecuaciones lineales. La matriz de dicho sistema
de ecuaciones se llama matriz de rigidez del sistema. El nmero de ecuaciones de dicho
sistema es proporcional al nmero de nodos.
(Fig. 11) Discretizacion de dominios en elementos finitos
Los resultados pueden fcilmente graficarse y as obtener una mejor representacin

visual del comportamiento de la estructura. Esto es lo que hace particularmente tan
atractivo al Mtodo de Elemento Finito al contar con herramientas visuales para
representar una simulacin del comportamiento de la estructura.

Si la malla del modelo contiene un nmero elevado de nodos, mayor ser la

complejidad para la resolucin matemtica. Si se tiene una estructura con un nmero
demasiado grande de nodos, resolver las ecuaciones de este elemento de forma
manual sera una tarea prcticamente imposible. Por ellos es necesaria la utilizacin
de un software basado en el Mtodo de Elemento Finito.
Los softwares de elemento finito tienen su funcionalidad bsicamente en tres partes:
Funciones del preprocesador. Consiste en la definicin de geometra, generacin

de la malla, las condiciones de contorno y asignacin de propiedades a los
materiales y otras propiedades. En ocasiones existen operaciones cosmticas
de regularizacin de la malla y preacondicionamiento para garantizar una
mejor aproximacin o una mejor convergencia del clculo.
Funciones del procesador. Abarca la solucin matemtica de las ecuaciones de

equilibrio y de compatibilidad establecidas en la matriz de rigidez. Adems se
definen los parmetros y criterios de resolucin.
Funciones del postprocesador. Contiene las herramientas, rutinas y macros para

poder visualizar a travs de grficos los resultados de la solucin matemtica.
(Fig. 12) Procedimiento de resolucin de problemas mediante software
A continuacin se describir ms detalladamente todo lo indicado anteriormente:

- Preproceso y generacin de la malla.
La malla se genera y sta en general consta de miles (e incluso de centenares de

miles) de puntos. La informacin sobre las propiedades del material y otras
caractersticas del problema se almacena junto con la informacin que describe
la malla. Por otro lado, las fuerzas, los flujos trmicos o las temperaturas se
reasignan a los puntos de la malla. A los nodos de la malla se les asigna una
densidad por todo el material dependiendo del nivel de la tensin mecnica u
otra propiedad. Las regiones que recibirn gran cantidad de tensin tienen
normalmente una mayor densidad de nodos (densidad de malla) que aquellos
que experimentan poco o ninguno. Puntos de inters consisten en: puntos de
fractura previamente probados del material, entrantes, esquinas, detalles
complejos, y reas de elevada tensin. La malla acta como la red de una araa
en la que desde cada nodo se extiende un elemento de malla a cada nodo
adyacente. Este tipo de red vectorial es la que lleva las propiedades del
material al objeto, creando varios elementos.
Las tareas asignadas al preproceso son:
El conjunto se divide, mediante lneas o superficies imaginarias en un

nmero de elementos finitos. Esta parte del proceso se desarrolla
habitualmente mediante algoritmos incorporados a programas
informticos de mallado durante la etapa de preproceso.
Se supone que los elementos estn conectados entre s mediante un
nmero discreto de puntos o nodos, situados en sus contornos. Los
desplazamientos de estos nodos sern las incgnitas fundamentales del
problema tal y como ocurre en el anlisis simple de estructuras o
elementos mecnicos por el mtodo matricial.
Se toma un conjunto de funciones que definan de manera nica el
campo de desplazamientos dentro de cada elemento finito, en funcin
de los desplazamientos nodales de dicho elemento. Por ejemplo, el
campo de desplazamientos dentro de un elemento lineal de dos nodos
podra venir definido por: u = N1 u1 + N2 u2, siendo N1 y N2 las
funciones comentadas (funciones de forma) y u1 u u2 los
desplazamientos en el nodo 1 y en el nodo 2.
Estas funciones de desplazamientos definirn entonces de manera nica
el estado de deformacin del elemento en funcin de los
desplazamientos nodales. Estas deformaciones, junto con las
propiedades constitutivas del material, definirn a su vez el estado de
tensiones en todo el elemento y, por consiguiente, en sus contornos.
Se determina un sistema de fuerzas concentradas en los nodos, tal que
equilibre las tensiones en el contorno y cualesquiera carga repartidas,
resultando as una relacin entre fuerzas F y desplazamientos u de la
forma F = k u , que es similar a la del clculo matricial.

- Clculo y resolucin de sistema de ecuaciones.
En un problema mecnico lineal no dependiente del tiempo, como un problema

de anlisis estructural esttico o un problema elstico, el clculo generalmente
se reduce a obtener los desplazamientos en los nodos y con ellos definir de
manera aproximada el campo de desplazamientos en el elemento finito.
Cuando el problema es no lineal, en general, la aplicacin de las fuerzas
requiere la aplicacin incremental de las mismas y considerar incrementos
numricos, y calcular en cada incremento algunas magnitudes referidas a los
nodos. Algo similar sucede con los problemas dependientes del tiempo, para
los que se considera una sucesin de instantes, en general bastante cercanos al
tiempo, y se considera el equilibrio instantneo en cada instante. En general
estos dos ltimos tipos de problemas requieren un tiempo de clculo
substancialmente ms elevado que en un problema estacionario y lineal.
- Postproceso.
Actualmente, el MEF es usado para calcular problemas tan complejos que los
ficheros de resultados que generan los programas de clculo MEF tienen tal
cantidad de datos que resulta conveniente procesarlos de alguna manera
adicional para hacerlos ms comprensibles e ilustrar diferentes aspectos del
problema. En la etapa de postproceso los resultados obtenidos de la resolucin
del sistema son tratados para obtener representaciones grficas y obtener
magnitudes derivadas, que permiten extraer conclusiones del problema.
El postproceso del MEF generalmente requiere software adicional para
organizar los datos de salida, de tal manera que sea ms fcilmente
comprensible el resultado y permita decidir si ciertas consecuencias del
problema son o no aceptables. En el clculo de estructuras por ejemplo, el
postproceso puede incluir comprobaciones adicionales de si una estructura
cumple los requisitos de las normas pertinentes, calculando si se sobrepasan
tensiones admisibles o existe la posibilidad del pandeo en la estructura
Con objeto de clarificar las ideas del apartado anterior aplicaremos los conceptos all
expuestos a la resolucin de un caso. Se trata de obtener las ecuaciones (matriz de
rigidez y vectores de cargas y desplazamientos) para resolver el problema elstico.
En primer lugar trataremos de obtener las ecuaciones que rigen el comportamiento de

un elemento triangular como el de la figura inferior.

Las funciones de interpolacin de los desplazamientos dentro del elemento se

consideran lineales. Es decir
donde u y v son los desplazamientos horizontal y vertical respectivamente. La

ecuacin anterior puede ser escrita en forma matricial
Particularizando las coordenadas y los desplazamientos para cada nodo obtenemos la

expresin matricial

Esta expresin nos permite obtener los parmetros de las funciones de interpolacin
en funcin de los desplazamientos nodales sin ms que invertir una matriz.
Reordenando los distintos trminos podemos escribir
donde
Luego ya conocemos la matriz [N] que nos relaciona el campo de desplazamientos en

el elemento con los desplazamientos en los nodos
Realizando la inversa de la matriz A, podemos reescribir la matriz N en funcin de las

caractersticas geomtricas del elemento
donde los valores de Ni vienen dados por
La matriz [D] que relaciona deformaciones y desplazamientos es

Y podemos construir la matriz [B] que nos relaciona las deformaciones con los
desplazamientos en los nodos.
Sustituyendo los valores de las distintas matrices
dado que los valores de las funciones Ni son conocidos en funcin de las coordenadas
nodales, es posible determinar la matriz [B] a partir de datos puramente geomtricos
siendo el rea del tringulo.

La matriz de comportamiento [C] que relaciona deformaciones y tensiones, para el
caso de tensin plana, viene dada por la relacin
donde E es el mdulo de elasticidad y m el coeficiente de Poisson. Con las matrices

definidas o calculadas hasta el momento ya es posible determinar las matrices de
rigidez y el vector de cargas de cada uno de los dos elementos del modelo que hemos
realizado
y dado que las matrices estn en funcin de las coordenadas nodales es posible
escribir

Siendo t el espesor de la placa.
5. DESARROLLO DE EJEMPLOS DE APLICACIN.
Cilindros a Presin.
Para la resolucin de este problema se va a utilizar el programa Abaqus 6.8. Se

resolver por el mtodo de los elementos finitos (MEF) problemas sobre los cilindros
a presin (es el caso, por ejemplo de depsitos o recipientes) (Fig.12).
(Fig.12)Esquema general de un cilindro a presin

En este problema se plantean varias cuestiones. Se trata de mostrar, por ejemplo,
como los resultados de la teora aproximada de cilindros de pared delgada se va
ajustando a los resultados de los clculos MEF segn se va reduciendo el espesor de la
pared del cilindro. Calcularemos los errores que se producen entre la teora exacta, la
aproximacin de pared delgada y los resultados obtenidos con el programa de clculo
MEF. Tambin se mostrar la distribucin de tensiones para cilindros estudiados
(cuatro) en la seccin central de todos ellos, y como vara segn nos movemos en
direccin radial por la pared de cada cilindro.
Se ha considerado que la presin externa es cero (p0 = 0) en todos los casos, lo que
reduce las ecuaciones para las tensiones tangenciales y radiales de la teora exacta a
las siguientes expresiones:
Donde pi es la presin interior y r0 y ri el radio exterior e interior del cilindro.

Las frmulas equivalentes segn la teora aproximada para cilindros de pared delgada
son:
Donde p es la presin, r el radio medio y t el espesor.

MODELADO DEL PROBLEMA CON ELEMENTOS FINITOS

A continuacin, se pasa a describir el mtodo utilizado para modelar el problema con
el programa de clculos avanzados por elementos finitos ABAQUS 6.8. Se describirn
las operaciones a realizar as como los comandos necesarios para llevar a cabo el
modelado. Todo clculo con elementos finitos se puede dividir en tres operaciones un
preproceso, proceso y postproceso. El proceso consiste fundamentalmente en los
clculos internos que efecta el programa utilizado, por lo que aqu se explicarn las
otras dos partes.
- Preproceso.
Geometra.
En este problema que se ha decidido trabajar en milmetros y milmetros cuadrados,

se ha credo que son unidades ms convenientes para analizar de mejor manera la
variacin de todos los factores que son estudio de este problema.
Lo primero es definir la geometra del modelo a estudiar. Aunque el problema real sea
en 3 dimensiones, se puede estudiar una seccin longitudinal del cilindro reduciendo
el problema a uno bidimensional. Dada la simetra cilndrica, se ha optado por trabajar
en coordenadas cilndricas siendo el eje radial el horizontal, el tangencial el
perpendicular a la pantalla y el eje Z el del sentido longitudinal del cilindro (Vertical
en este caso). Bastar con trazar un nico rectngulo con la altura del cilindro y el
ancho correspondiente al espesor (seccin longitudinal). Se crear con la opcin de
axisimetra para que as nuestro rectngulo se interprete como una superficie que
tiene que girar alrededor de un eje para formar un cilindro.
La longitud l de los cuatro cilindros considerados ser de 120 mm.
Geometras estudiadas:
El caso 1 tiene las siguientes dimensiones: un dimetro interior di de 236 mm y un

dimetro exterior d0 de 240 mm, por lo que tiene un espesor t de 2 mm.
El caso 2 tiene las siguientes: un dimetro interior di de 232 mm y un dimetro

exterior d0 de 240 mm, por lo que tiene un espesor t de 4 mm.
Ambos casos cumplen lo que es un Cilindro de pared delgada ya que cumplen la

siguiente relacin:
Para el caso 1: 2(mm) < 11.8 (mm) OK
Para el caso 2: 4(mm) < 11.6 (mm) OK

Materiales.
En cuanto a sus dimensiones, ya estn todos los cilindros definidos. A continuacin

sea describir el material. Ser comn a todos, lo que es lgico, ya que se desea
comparar la variacin que sufren las tensiones y deformaciones simplemente
modificando el espesor. Se trabajar con cilindros de acero estndar con las siguientes
caractersticas:
Mdulo de elasticidad E: 210000 N/mm2

Coeficiente de Poisson : 0.3.
Siempre se trabajar en la zona elstica del material, no se considerar ningn efecto

plstico. Tampoco interviene en ninguna parte del problema otro tipo de efectos como
dilataciones trmicas.
Condiciones de contorno.
Otra propiedad comn para los cuatro casos antes expuestos es que todos ellos son cilindros
sin tapas, es decir, toda la presin la soportar la pared del cilindro.
Hay que introducir tambin las condiciones de contorno y las cargas a las que est
sometido el cilindro. Aplicaremos una carga distribuida a lo largo de la lnea interior
del rectngulo que representa el cilindro, el programa la reinterpreta como una carga
distribuida o presin aplicada a toda la superficie interior del cilindro. En nuestro
ejercicio se ha empleado una presin de 2.4 MPa o, lo que es lo mismo, 2.4 N/mm2,
igual para los cuatro casos. Por ltimo, en esta parte hay que restringir los grados de
libertad o los desplazamientos y rotaciones necesarias de los puntos adecuados. En
este caso hemos impedido el desplazamiento vertical (U2 = 0) de todos los nodos de la
seccin central transversal del cilindro (Fig.13).
(Fig.13)Esquema de condiciones de contorno

Mallado.
En esta parte es donde se define el tipo de elemento con el que se va a mallar y el

nmero de elementos y de nodos con los que se va a discretizar el slido. Gracias al
grado de finura de esta discretizacin, el programa, en principio, puede obtener
resultados tan exactos como se desee, esta es la base del procedimiento de clculo por
elementos finitos.
En este caso hemos utilizado elementos rectangulares y de tipo estructural, definido
como cax4r. Aunque el tipo de elemento del mallado es igual en los cuatro cilindros,
no es igual el nmero de elementos y, por lo tanto, el nmero de nodos de cada caso,
debido a la geometra. Incluso no es el mismo nmero de nodos en la seccin central
pues vara segn el espesor. A continuacin, se indican los datos concretos de cada
caso:
Caso 1: tiene 828 nodos y 618 elementos. La seccin central contiene 4 nodos.
Caso 2: tiene 966 nodos y 800 elementos. La seccin central contiene 6 nodos.
Hay que destacar que la versin del programa disponible est restringida a 1000
nodos.
(Fig.14) Detalle de mallado del caso 2

Despus de introducir toda esta informacin en el programa, se lanza el clculo (proceso) y

se obtienen los resultados.
- Postproceso.
En esta parte se interpretan los resultados del clculo del programa, ya que simplemente
nos proporciona la solucin numrica del problema pero, en todo caso, cmo interpretarlos
y cmo utilizarlos es tarea del ingeniero.
A partir de aqu se analizan los resultados obtenidos y si utilizando MEF, hemos llegado a
lo que la teora predice.
Se van a estudiar y comparar los resultados obtenidos en cada uno de los 2 casos. Slo se
va a considerar la seccin transversal central del cilindro, ya que las tensiones son
independientes de la variable longitudinal Z.
Caso 1.
Este caso, como se dijo en la introduccin, corresponde al cilindro de menor espesor de

todos (2 mm.) y con 4 nodos en la seccin central. Con los datos de la presin interior y la
geometra del cilindro, tambin conocida, se calculan las tensiones por la teora de cilindros
a presin y segn la aproximacin para pared delgada. Por otro lado, Abaqus aporta
resultados para estas tensiones una vez llevadas a cabo todas las operaciones sealadas en
el apartado del preproceso. En la tabla siguiente se pueden ver todos los resultados
conjuntamente:
Aqu vemos las tensiones tangenciales y radiales para los nodos centrales obtenidas
por los tres mtodos: la teora exacta, la aproximada y por MEF. A primera vista se
pueden extraer dos conclusiones, primera, que la teora aproximada de paredes
delgadas es bastante exacta, y, segunda, que los tericos se parecen mucho a los
resultados MEF obtenidos va Abaqus. Para verificar todo esto, a continuacin se
calculan los errores relativos entre las tensiones calculadas de diferente forma.

El error relativo depende del valor de r donde se calcula segn:
Tangenciales:
Errores entre teora exacta y teora aproximada para pared delgada.
Errores entre teora exacta y MEF.
Errores entre teora aproximada para pared delgada y MEF.

Radiales:
Las tensiones radiales en la teora aproximada para cilindros de pared delgada tienen
un valor de cero. En la tabla 3.1 se puede ver que las tensiones radiales en las
columnas pared gruesa y clculo MEF son despreciables frente a las tangenciales,
como predice la teora cilindros de pared delgada.
Todo esto quedar a continuacin perfectamente claro en los grficos comparativos

de los tres procedimientos de clculo. En este apartado, para las tensiones radiales
solamente se muestran los errores relativos entre la teora exacta y MEF.
Como puede observarse, los resultados obtenidos por el mtodo de elementos finitos
con el programa Abaqus son bastante fieles a la teora exacta. Por otro lado, vemos
como al ser un cilindro de pequeo espesor, puede aplicarse la teora aproximada
para cilindros de pared delgada sin cometer un error apreciable. En la Fig.15 se puede
apreciar la distribucin de tensiones tangenciales t dibujadas segn los tres mtodos
de clculo en funcin del radio. Se puede observar el gran acuerdo entre la teora
exacta y el MEF siendo coincidentes en los valores de r interiores al cilindro y
difiriendo muy poco en los puntos de las paredes (esto se puede deber al tipo de
mallado).
Por otro lado, se aprecia la lnea de tensin constante que se obtiene con la teora
aproximada de pared delgada. La diferencia entre una y otras es ms bien aparente,
pues la mxima diferencia est en torno a 3 N/mm2 (en r = r0) que, frente a valores de
la tensin en torno a 140 N/mm2, resulta ser irrelevante.

(Fig.15) Distribucin de tensiones tangenciales para el caso 1 (t = 2mm).
En la figura 16 se ha dibujado la distribucin de tensiones radiales para el mismo caso.

Al igual que en la grfica anterior, vienen representadas las tensiones calculadas
segn los tres mtodos utilizados en el problema, y vienen dadas en funcin del radio
r. Aqu destaca, en un primer lugar, el pequeo valor absoluto (son tensiones
negativas de compresin) que toman las tensiones en comparacin con las tensiones
tangenciales, lo que indica cules van a ser las tensiones ms importantes a la hora del
clculo y diseo de piezas de este tipo.
Por otro lado, se aprecia tambin la gran coincidencia entre la teora exacta y MEF,
excepto una pequea distorsin en los nodos sobre las paredes del cilindro. La teora
aproximada considera que la tensin radial es de valor nulo y, frente a este, el valor
mximo en valor absoluto de r, segn los otros dos procedimientos de clculo, no
sobrepasa los 2.5 N/mm2, casi despreciable. Por ello, r puede aproximarse por el
valor nulo que proporciona la teora de cilindros de pared delgada.

(Fig.16) Distribucin de tensiones radiales para el caso 1 (t = 2mm)
Otro consideracin a tener en cuenta que se puede apreciar en las grficas anteriores,
es la variacin de las tensiones en funcin del radio. Como vemos, en la grfica 1
(tensiones t) la tensin es mayor en el interior del cilindro y segn nos movemos
hacia la parte exterior, es decir, conforme r crece, va disminuyendo. (Fig.17). Si nos
fijamos en la grfica de las tensiones radiales podemos ver que segn nos acercamos
al exterior desde el interior del cilindro,r tiende a cero (Fig.18). En realidad, las
tensiones radiales crecen con r al ser negativas aunque, obviamente, disminuyen si se
considera la variacin de su valor absoluto.
(Fig.17) detalle de la distribucin de las tensiones tangenciales en el caso 1.

(Fig.18) detalle de la distribucin de las tensiones radiales en el caso 1
En estas dos figuras, 17 y 18, se puede observar la variacin cromtica de ambas

tensiones en el espesor segn los resultados obtenidos con Abaqus. Estas imgenes
permiten visualizar grficamente los valores numricos de las distribuciones de
tensiones y, as, saber cul es la zona ms esforzada del cilindro.
Caso 2.
Este caso, como se explic al principio de este ejercicio, corresponde a un espesor

mayor en relacin con el anterior (caso 1). Ahora tenemos un espesor de 4 mm. con 6
nodos en la seccin central. Se muestra en la siguiente tabla los resultados de las
tensiones segn los tres mtodos de clculo: teora exacta, teora aproximada y MEF.
En muchos aspectos los resultados son parecidos al caso 1.

A primera vista se puede observar que los resultados son bastante parecidos en los
tres casos. Se han calculado los errores relativos para ver el grado de acuerdo entre
los tres mtodos de clculo.
Tangenciales.
Errores entre teora exacta y teora pared delgada.
Errores entre teora pared delgada y MEF.

Radiales:
En este caso se puede observar como, al aumentar el espesor de la pared, el error

relativo entre la teora exacta y la aproximacin para cilindros de pared delgada
tambin lo ha hecho, aunque todava hay un gran acuerdo entre ellas al ser todava el
espesor relativamente pequeo. Por otro lado, los resultados obtenidos con Abaqus
son casi coincidentes con los calculados por la teora de cilindros a presin.
A continuacin, se muestra la figura 19, donde se ha dibujado la distribucin de
tensiones tangenciales t para este caso en funcin del radio. Volvemos a apreciar el
gran acuerdo entre la teora exacta y MEF con la pequea discrepancia en el primer y
ltimo nodo (nodos en r = ri y r = ro). En cuanto a la teora aproximada, volvemos a ver
como hay una diferencia mxima de aproximadamente 2.5 N/mm2 con respecto a los
otros resultados (igual que ocurra en el caso 1). Sin embargo, dado que un aumento
del espesor reduce el valor de las tensiones, han pasado del orden 140 N/mm2 (caso
1) a estar en torno a70 N/mm2, una misma diferencia conduce un mayor error
relativo.
(Fig.19) Distribucin de tensiones tangenciales para el caso 2 (t = 4mm).

La figura 20 es muy similar y muestra unos valores muy parecidos a los de la figura 16
los valores mximos de las tensiones r calculadas va MEF o con la teora exacta son
muy pequeos (no llegan, en valor absoluto a 2.5N/mm2). Aunque tambin en este
caso podran aproximarse por el valor nulo que predice la teora de cilindros de pared
delgada, frente a las tensiones tangenciales se han hecho ms relevantes al reducirse
significativamente el valor mximo de stas como se comentado anteriormente.
(Fig.20) Distribucin de tensiones radiales para el caso 2 (t = 4mm).
Por ltimo, se muestra las distribuciones de tensiones cromticas que nos

proporciona Abaqus. En primer lugar se muestra la distribucin de tensiones
tangenciales (Fig.21) y luego las radiales (Fig.22). Las figuras son similares a las del
caso 1.
(Fig.21) Detalle de la distribucin de las tensiones tangenciales para el caso 2.

(Fig.22) Detalle de la distribucin de las tensiones radiales para el caso 2.
A modo de comparacin se va a introducir la siguiente tabla en forma de resumen con

los errores relativos promedios entre Abaqus y la teora exacta; y entre la teora
exacta y la aproximada, para poder comparar de un solo vistazo los resultados
obtenidos para los dos casos. Se har solamente paras las tensiones tangenciales t
por alcanzar valores mayores.
Esta pequea tabla resumen nos da una muy buena visin de cmo la teora
aproximada para cilindros de pared delgada no se puede aplicar a partir de un cierto
valor de la relacin entre espesor y radio (la bibliografa propone una relacin lmite
de1/10). Resumiendo, la teora de cilindros de pared delgada podra aplicarse sin
cometer un error significativo en ejemplos 1 y 2.

6. CONCLUSIONES.
En este trabajo de investigacin se vio que el Mtodo de Elementos Finitos se puede

aplicar a diferentes tipos de estructuras, modelos geomtricos, etc. en si se acomoda
de acuerdo a lo que necesitemos calcular; y para eso simplemente debemos conocer
bien las ecuaciones que rigen el problema que se analizara.
Para nuestro caso se pudo resolver un cilindro de pared delgada en el cual se calcul
las presiones que actuaban sobre ella, tambin se dio a conocer que para que un
cilindro sea considerado de pared delgada debe cumplir una relacin en la que el
espesor es menor a el radio dividido entre diez.
De los resultados alcanzados se concluye que se han encontrado niveles de

aproximacin muy buenos con respecto a los que proporcionan las correspondientes
teoras. En la mayora de los casos no se han superado errores relativos del 5%. Para
la visualizacin de todos los resultados se han dibujado multitud de grficas y
comparando curvas tericas con calculadas. Se han presentado tambin multitud de
tablas con errores relativos y resultados numricos.
Tambin como recomendacin podemos decir que es aconsejable el uso de software

por que incluso en una estructura pequea o placas simples se tendrn varios nudos y
a la vez una infinidad de mallas, las cuales nos darn como resultados muchas
matrices para calcular y hacerlas manualmente seria trabajo de das.

Cilindros Pared Delgada

Cargado por

Copyright:

Formatos disponibles

Cilindros Pared Delgada

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Cilindros Pared Delgada

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD MAYOR REAL Y PONTIFICIA DE SAN FRANCISCO XAVIER DE CHUQUISACA

MECANICA ESTRUCTURAL AVANZADA

APLICACIN DEL MTODO DE ELEMENTO FINITO EN EL

El Mtodo de Elemento Finito (MEF en castellano o FEM en ingls) es un mtodo

2.1. Mtodo de Elemento Finito.

El Mtodo del Elemento Finto (MEF), es una tcnica de anlisis matemtico,

Tabla 1. Aplicaciones comunes en diferentes campos de ingeniera

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

- TIPOS DE ELEMENTOS FINITOS.

Los tipos de elementos finitos se pueden clasificar de la siguiente manera:

Elemento Finito Puntual.- Es un elemento que no tiene dimensin y su

(Fig. 1) Elemento Finito Lineal

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

(Fig. 2)Elementos Finitos Superficiales Placa

Elementos Finitos Volumtricos Slidos.- Son la generalizacin tridimensional de

(Fig. 3) Elementos Finitos Volumtricos Solidos

Elementos Finitos Volumtricos Axisimtricos.- Utilizados para la resolucin de

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

2.2. Esfuerzo en Cilindros.

2.3. Cilindros de pared delgada.

Considerando recipiente cilndrico de radio interior r y espesor de pared t, que

El anlisis de pared delgada se refiere a los recipientes con la siguiente relacin:

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

3. ANALISIS MATEMATICO DEL TEMA.

El Mtodo de Elemento Finito viene dado por ecuaciones fundamentales y a partir de

- Funciones de los desplazamientos. En la figura 4 se muestra un elemento

Los desplazamientos interiores a un elemento han de quedar definidos

Se pueden calcular fcilmente las seis constantes resolviendo los dos

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

coordenadas de los nodos e igualar las expresiones resultantes a los

Podemos calcular fcilmente 1, 2 y 3 en funcin de los desplazamientos

Obtenindose los otros coeficientes mediante permutacin cclica de los

Puesto que las ecuaciones para el desplazamiento vertical don similares, se

Aunque no sea estrictamente necesario en este momento, podemos expresar

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

Siendo I una matriz unidad de 2x2, y

La funcin de desplazamientos elegida garantiza automticamente la

- Deformacin (total). La deformacin total en cualquier punto del elemento

Con lo cual, la matriz B queda definida explcitamente. Se advertir que en este

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

- Matriz de elasticidad. La matriz D de la relacin

Tensin plana de un material istropo. Para los estudios de tensin plana en

Despejando las tensiones en el sistema anterior, obtenemos la matriz D

En el cual E es el mdulo de elasticidad y el coeficiente de Poisson.

Deformacin plana de un material istropo. En este caso la tensin normal z

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

Eliminando z y despejando las tres tensiones restantes de la expresin ya

- Matriz de Rigidez. La matriz de rigidez de un elemento ijm viene definida a

En nuestro trabajo de investigacin se analizara el Mtodo de Elementos Finitos

Y cambiara nuestro campo de desplazamientos por lo siguiente:

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

Donde r y z representan respectivamente las coordenadas radial y axial de un punto y

- Funciones de Interpolacin. Al ser la deformacin plana, el campo de

Siendo Ui las deformaciones en direccin radial y Vi las deformaciones en

donde n es el nmero de nudos del elemento. La expresin tambien puede

La matriz de funciones de interpolacin N tiene dos filas y tantas columnas

MAESTRIA EN INGENIERA ESTRUCTURAL GRUPO 2

En esta expresin es la deformacin unitaria circunferencial, producida por la

En esta expresin se identifica al operador matricial que pasa de las