Revista de la Facultad de Ingeniería U.C.V., Vol. 24, N° 1, pp.

13–26, 2009

METODOLOGÍA DE MODELACIÓN MEDIANTE ANSYS DE LA HISTORIA TÉRMICA,

TENSIONES Y DEFORMACIONES DE SOLDADURA

Juan Pozo Morejón, Jorge García Jacomino, Félix Ramos Morales, Yenei García Rodríguez,
Amado Cruz Crespo, Eduardo Díaz Cedré, Alejandro Duffus Scott

Universidad Central “Marta Abreu” de Las Villas, Facultad de Ingeniería Mecánica, Centro de Investigación de
Soldadura (CIS). Santa Clara, Villa Clara, CP. 54830, Cuba. e-mail: jpozo@uclv.edu.cu

Recibido: junio de 2008 Recibido en forma final revisado: diciembre de 2008

RESUMEN

En este trabajo se presenta una metodología de modelación de la evolución de la historia térmica, tensiones y deformaciones
producidas por el proceso de soldadura, mediante el software de elementos finitos de uso general ANSYS Multiphysics
versión 9.0. Dicha metodología integra una información que se encuentra dispersa en la literatura, proponiendo variantes
que en ocasiones superan lo publicado en trabajos precedentes. La metodología es aplicada exitosamente en la modelación
tridimensional de una soldadura GMAW automática sobre chapas de aleación de aluminio al magnesio 5083-O, de amplio
empleo en la industria moderna, validándose los resultados que brinda el modelo contra datos experimentales publicados.
El estudio permitió establecer que las mejores respuestas, para este tipo de aleación, se obtienen empleando una temperatura
de corte igual a la temperatura de sólido del metal, utilizando los modelos de material bilineal sin endurecimiento, o con
endurecimiento cinemático, y activando los elementos de la soldadura en el instante en que sobre ellos se encuentra la
fuente de calor en ambas corridas térmica y mecánica.

Palabras clave: Soldadura, Modelación, ANSYS, Aluminio, Tensiones, Deformaciones.

MODELATION METHODOLOGY BY MEANS OF ANSYS OF THE WELDING

THERMAL HISTORY, STRESSES AND STRAINS

ABSTRACT

A modeling methodology of the evolution of thermal history, stresses and strains due to the welding process, using the
general purpose finite elements software ANSYS Multiphysics release 9.0, is exposed. This methodology integrates
existing information and new variants that, in occasions, provide better results than those previously published results. This
methodology was applied to the three-dimensional modeling of automatic GMAW on 5083-O aluminum - magnesium alloy
plates, which are frequently employed in the modern industry. Model results were validated against published experimental
dates. The study allowed to define that for this alloy type the best results are obtain using a cut-off temperature similar to
the metal solid temperature, employing the bilinear material models without hardening or using kinematics hardening, and
activating welding elements when thermal source is right over them, during both thermal and mechanic races.

Keywords: Welding, Modeling, ANSYS, Aluminum, Stress, Strains.

INTRODUCCIÓN transitoria con complicadas condiciones de frontera, trans-

formaciones de fases y estados tensionales y deformacio-
En décadas pasadas se han realizado múltiples investigacio- nales complejos, entre otros. Debido, por un lado, a que la
nes dirigidas a procedimientos analíticos que simulen con complejidad de los procesos físicos involucrados hace que
más precisión los procesos de soldadura. Como se conoce ninguna solución matemática simple logre dar respuestas
el proceso de soldeo es bastante complejo e involucra: el aceptables, mientras por el otro, que en muchos casos no es-
movimiento de una fuente de calor, depósito de material, tán disponibles técnicas experimentales que determinen los
propiedades termodependientes, transferencia de calor campos totales de tensiones y deformaciones en construc-
13
ciones soldadas, sumado al desarrollo cada vez mayor de comercial ANSYS Multiphysics versión 9.0, la que poste-
los medios de cómputo, es lo que ha hecho que desde hace riormente se trata de validar modelando un experimento de
varios años se vengan imponiendo las simulaciones reali- soldadura GMAW a tope de chapas de aleación de alumi-
zadas en computadoras, mediante métodos numéricos y en nio, realizado por Cañas et al. (1996).
especial mediante los métodos de elementos finitos (MEF).
ASPECTOS TEÓRICOS SOBRE LOS ANÁLISIS
A pesar de que las publicaciones sobre el tema de la mo- TÉRMICO Y ESTRUCTURAL
delación de soldaduras mediante MEF se han multiplica-
do en la actual década, sin dudas éste aún es un campo en El análisis mediante MEF del problema térmico en la sol-
el que existe mucho debate, se carece en muchos casos de dadura se basa en la ecuación (1) de transferencia del calor
criterios unificados, la información se encuentra dispersa y en un cuerpo (Zhu et al. 2002):
en el que se requiere seguir avanzando en aras de resolver
∂T  ∂  ∂T ∂  ∂T  ∂ ∂T 
importantes problemas de la industria. Dentro de los temas
r(T)c(T)
ñ(T)c(T) =q+
Kx(T)  + Ky(T)  + Kz(T) (1)
∂t ∂x  ∂x  ∂y  ∂y  ∂z  ∂z 
que actualmente se discuten se pueden mencionar: si es más
conveniente construir modelos bidimensionales (2D) o tri- donde:
dimensionales (3D), los modelos de material que se deben
usar, incluyendo si introducir propiedades termodependien- ρ(T) representa la densidad del material como propiedad
tes o constantes a temperatura ambiente, sobre la forma termodependiente, c(T) el calor específico, q el calor ge-
más conveniente de introducir el calor al modelo, sobre el nerado por unidad de volumen, Kx(T), Ky(T), Kz(T) son
momento más oportuno de activar los elementos cuando se los coeficientes de conductividad térmica en las tres direc-
emplea el método de activación-desactivación (birth-death) ciones (en materiales isotrópicos este valor es el mismo en
para simular el depósito de material, sobre como modelar cualquier dirección), T es la temperatura y t es el tiempo.
los cambios de fase, sobre si es conveniente o no el empleo Como se desprende, el problema térmico se rige por una
de una temperatura de corte, entre otros. ecuación diferencial no lineal, debido a que las propiedades
termofísicas del material son dependientes de la tempera-
En la literatura se reporta el empleo de dos tipos de códigos tura.
de elementos finitos en la modelación de tensiones y defor-
maciones de soldadura: códigos especializados en modelar Las pérdidas de calor al medio circundante por radiación y
este proceso como el SYSWELD o el WELDSIM, y códi- convección se consideran con el empleo de las ecuaciones
gos de uso general como ANSYS, ABAQUS, COSMOS M, (2) y (3) (Zhu et al. 2002):
COSMOS DESING STAR, entre otros.
qc= h(T-T00) (2)
ANSYS apareció en 1971, o sea, que es un código MEF que
qr= eB(T4-T004) (3)
posee más de treinta años de existencia y es capaz de rea-
lizar análisis estáticos, dinámicos, de transferencia de ca-
donde:
lor, de flujo de fluidos y de electromagnetismo, entre otros
(ANSYS Inc. 2004), que por más de veinticinco años ha
h representa el coeficiente de convección, T00 la temperatura
sido líder dentro de estos programas. ANSYS, a pesar de
ambiente, e la emisividad de la superficie del cuerpo y B la
no ser un código especializado, sin dudas constituye en la
constante de Stefan-Boltzmann.
actualidad uno de los más empleados en la modelación de
tensiones y deformaciones de soldadura, tanto en modelos El análisis estructural de la soldadura mediante MEF es aún
2D ó 3D y por las posibilidades que brinda constituye una más complejo que el anterior. Este análisis se basa en cal-
potente herramienta para tales efectos, cuyo empleo en los cular la deformación total que se produce producto del pro-
últimos años ha sido reportado por múltiples autores (Saha, ceso de soldeo, la que se expresa por la relación (Bezerra,
2005; Camilleri et al. 2005; Kim et al. 2006; Armentani et 2006; ANSYS Inc. 2004):
al. 2006; Bezerra, 2006; Camilleri et al. 2005; Hemanth,
2006; Mahapatra et al. 2006; Moraitis et al. 2007; Meo et eTotal = eterm + ee + ep (4)
al. 2007; He et al. 2007; Jerome et al. 2007; Pozo et al.
2008). donde:

Por todo lo anterior, en el presente trabajo se exponen re- eterm representa la deformación térmica provocada por el ca-
sumidamente los aspectos generales de una metodología de lor del arco, ee la deformación elástica y ep la deformación
modelación de soldaduras, mediante el código MEF de uso plástica.
14
La deformación térmica del cuerpo se obtiene a través de METODOLOGÍA DE MODELACIÓN POR MEF DE
la conocida relación física de dilatación térmica (Bezerra, SOLDADURAS
2006; ANSYS Inc. 2004):
Con vistas a simplificar la simulación de soldaduras es efi-
eterm= a(T-Tref) (5) ciente computacionalmente ejecutar los análisis térmico y
estructural separadamente o de forma desacoplada. Físi-
donde: camente se asume que los cambios en el estado mecánico
no provocan cambios en el estado térmico; o sea, que las
α representa el coeficiente de dilatación lineal del material, deformaciones no provocan cambios de temperatura apre-
T es la temperatura instantánea y Tref es la temperatura de ciables, que merezca considerarse. Mientras que por otro
referencia para la que la deformación térmica es nula. lado se considera, que los cambios en el estado térmico
si causan cambios en el estado mecánico. La mayoría de
La componente elástica se determina, dentro del intervalo los autores han reportado el uso de esta forma de modelo
del límite elástico del material, directamente por la Ley de desacoplado (Chao et al. 1999; Francis, 2002; Zhu et al.
Hooke empleando el módulo de elasticidad termodepen- 2002; Saha, 2005; Kim et al. 2006; Armentani et al. 2006;
diente. La determinación de la componente plástica es la Bezerra, 2006; Moraitis et al. 2007; Meo et al. 2007; Pozo
más compleja, ya que primeramente se debe establecer el et al. 2008).
dominio elástico que normalmente se define en base al cri-
terio de fluencia de von Mises: El modelado de tensiones y deformaciones de soldadura
mediante MEF se realiza, según esquema mostrado, en dos
seq - se < 0 (6) etapas: primeramente una corrida térmica, y a continuación
la corrida estructural (figura 1). En la corrida térmica se in-
donde: troduce como carga al modelo, el calor proveniente del arco
y se obtiene como respuesta las temperaturas nodales para
se representa la tensión o límite de fluencia y seq la tensión cada paso de tiempo discretizado (historia térmica), hasta
equivalente de von Mises, definida por la conocida ecua- que se completa el enfriamiento. Seguidamente en la corri-
ción (7) (Bezerra, 2006; ANSYS Inc. 2004): da estructural se introducen como cargas estas temperaturas
nodales, para cada paso de tiempo, se calculan los despla-
1/2
seq= 1 [(sx - sy)2 + (sy - sz)2 + (sx - sz)2 + 6(txy+ tyz + txz)2] (7) zamientos y deformaciones elastoplásticas y finalmente se
2
obtiene como respuesta los desplazamientos y tensiones,
donde: transitorios y residuales.

sx, sy, sz,…, τxz son las tensiones normales y tangenciales El modo óptimo de modelación de soldaduras mediante AN-
que definen el estado tensional triaxial de un punto y seq la SYS es el modo combinado en el que se mezclan el modo
tensión normal uniaxial equivalente. interactivo y el modo “batch” o de programación. De esta
manera el usuario puede construir el modelo geométrico,
De modo que se considera que la deformación plástica se mallarlo, declarar el modelo de material, con las propieda-
inicia cuando se cumple la condición: des termodependientes, etc., por ejemplo mediante el modo
interactivo, y luego a través del modo “batch”, programar
seq - se = 0 (8) el movimiento de la fuente de calor en la corrida térmica, o
programar que ANSYS lea los resultados de temperaturas
Es importante destacar que la tensión equivalente nunca nodales para cada paso de tiempo y las aplique como cargas
puede exceder la tensión de fluencia del material, ya que en la corrida estructural, o simplemente chequear que no
si esto ocurriera se desarrollarían instantáneamente defor- se cometan errores durante la introducción de datos en el
maciones plásticas, que desplazarían nuevamente la tensión modo interactivo. Se debe aclarar, que lo anterior no es rí-
hacia dicho límite de fluencia (Bezerra, 2006; ANSYS Inc. gido, sino que cada usuario decidirá que partes del modelo
2004). La superficie representada por las ecuaciones (7) y le resulta más sencillo y rápido introducirlas en un modo o
(8) se conoce como superficie de fluencia del material. Esta en otro. A continuación se expone los aspectos generales de
superficie puede ampliarse, según el modelo de endureci- la metodología utilizada.
miento del material isotrópico, o desplazarse manteniendo
su forma y dimensiones según el modelo cinemático, am- La modelación mediante ANSYS de los análisis térmico y
bos implementados como modelos de material en el código estructural consta de tres etapas fundamentales: etapa de
ANSYS (ANSYS Inc. 2004). preprocesamiento, etapa de solución y etapa de postproce-
15
 Figura 1. Esquema general de la modelación de tensiones y deformaciones de soldadura mediante MEF.

samiento. calor, mediante el coeficiente de convección y la tem-

peratura del medio. Sólo se exceptúa el área de simetría
Metodología para el análisis térmico mediante ANSYS del modelo, si se modela la mitad de la unión aprove-
chando esta propiedad.
• Etapa de preprocesamiento:
Se debe destacar que las pérdidas de calor por radiación
1. Se elige un tipo de elemento, que posea la temperatura pueden ser de consideración a temperaturas sólo cercanas a
como grado de libertad. En el presente trabajo se emplea la de fusión, por lo que algunos autores únicamente le intro-
el elemento SOLID70 (ANSYS Inc. 2004), un elemento ducen al modelo las pérdidas por convección y consideran
sólido de ocho nodos, usado con éxito por otros auto- las pérdidas por radiación disminuyendo el coeficiente de
res para este propósito (Bezerra, 2006; Moraitis et al. eficiencia de la fuente de calor, o simplemente las despre-
2007). cian (Armentani et al. 2006; García et al. 2007; Ramos et
al. 2007). Si no se introduce ningún tipo de pérdida de calor
2. Se declaran como propiedades termodependientes del al modelo, luego de aplicado el calor de la soldadura, con
material, la conductividad térmica, el calor específico y el paso del tiempo la temperatura en dicho modelo tiende
la densidad, introduciendo el valor de la propiedad que a equilibrarse por conducción y nunca regresa a la tempe-
corresponde a cada valor de temperatura, entre la tem- ratura inicial o ambiente, independientemente del tiempo
peratura ambiente y la de fusión del metal. transcurrido, lo que es erróneo y no coincide con el fenó-
meno real.
3. Se construye el modelo geométrico tridimensional. El
modelo geométrico puede también ser importado a par- • Etapa de solución:
tir de un fichero con extensión IGES generado con otro
programa CAD (ANSYS Inc. 2004). 1. Se declara como tipo de análisis el transitorio en el tiem-
po. Se recomienda ejecutar la corrida térmica a través
4. Se malla el volumen, teniendo en cuenta refinar la malla de un análisis transitorio, como lo hacen Bezerra (2006)
en el lugar cercano a la soldadura, en lo que coinciden y Mahapatra (2006), que permite calcular como varían
todos los autores, por ser ésta la zona donde se produ- las temperaturas con el transcurso del tiempo, y abarca
cen los mayores gradientes térmicos. El usuario debe desde que se comienza a aplicar calor al modelo, hasta
cerciorarse de que el grado de refinamiento garantiza la que éste se termina de enfriar (Bezerra, 2006; Mahapa-
convergencia de resultados. tra et al. 2006).

5. Se declara como condición de frontera la pérdida de ca- 2. Se emplea el método iterativo de solución de ecuaciones
lor del modelo, por convección al medio, a través de no lineales Full Newton Raphson, que constituye el mé-
todas las superficies exteriores que físicamente emiten todo no lineal estándar de ANSYS (ANSYS Inc. 2004).
16
3. Se establece un valor de convergencia de la solución cuentra en la trayectoria del arco por donde se introduce el
para análisis no lineal de 0.1 ºC (Bezerra, 2006). calor (m3). Dicho flujo se introduce como carga al modelo,
para cada paso de tiempo y se retira en el siguiente paso.
4. Se selecciona el avance de tiempo automático (automa-
tic time stepping). Es importante destacar que el movimiento de la fuente se
programa a través de un pequeño y sencillo lazo *DO -
5. Se declara la temperatura inicial del metal base. *ENDDO, en lenguaje FORTRAN, dentro del que se in-
troduce el comando que aplica la carga térmica, durante el
tiempo en segundos calculado, permitiéndose que recorra
6. Se declara un tiempo de inicio de 1*10-6 s, suficiente-
paulatinamente todos los elementos por los que el arco
mente pequeño pero no nulo, como lo recomiendan pasa. El empleo de este lazo elimina la necesidad de cons-
Budgell (2007) y ANSYS Inc. (2004), ya que si se de- truir curvas de tiempo o funciones de carga, lo que reduce
clara un tiempo de inicio nulo en un análisis transitorio, el tiempo de preparación del modelo y supera una limita-
ANSYS asume el tiempo de inicio con valor igual a la ción de anteriores trabajos publicados (García et al. 2007;
unidad, lo que resulta erróneo. Ramos et al. 2007; Leal et al. 2007). La carga térmica se
aplica al modelo mediante la opción a salto STEPPED.
7. Si se modela una soldadura de recargue o de ranu-
ra de bordes rectos (sin preparación de bordes) no se Si se está empleando la técnica de activación-desactivación
debe emplear la técnica de activación – desactivación de elementos, es necesario, dentro del mismo lazo, activar
de elementos, por lo que este paso se obvia. Si por el los elementos finitos del cordón en el paso de tiempo en
contrario, se va a modelar el depósito de metal en una que se les va a aplicar la carga térmica a cada uno, o de
lo contrario los resultados serán erróneos. De esta manera
soldadura de ranura con preparación de bordes, o una
se garantiza que los elementos finitos debajo de la fuente
de filete, entonces es necesario desactivar todos los ele-
de calor retornen a su propiedad de conductividad térmica
mentos que forman parte del cordón, antes de aplicar original, con lo que físicamente se simula que haya metal
cualquier carga. Finalmente, se restringe la temperatura debajo de la fuente de calor, con su capacidad de conduc-
de los nodos que no están ligados a ningún elemento ción, y no aire.
activo, a un valor de temperatura igual a la ambiente, ya
que estos se encuentran en un lugar donde lo que hay es 9. Finalmente, en un último paso de tiempo se permite el
aire en el problema físico real. enfriamiento del modelo a través de la pérdida de calor
por convección hasta un tiempo estimado tal, en que se
8. Introducción de la carga térmica al modelo, a través de alcance el total enfriamiento.
la superficie superior o del volumen de los elementos
finitos que conforman la soldadura y modelación del 10. Solución mediante el comando respectivo, que se in-
desplazamiento de la fuente. De manera que la carga tercala al final de cada paso de tiempo para resolver el
térmica que se aplica al modelo se determina en base al problema.
régimen de soldadura utilizado, por la conocida ecua-
11. Una vez concluido el cálculo se guarda todo el modelo
ción (9):
elaborado y los resultados con la historia térmica en un
gran fichero con extensión DB.
Q = h*Is*Va (W) (9)
• Etapa de postprocesamiento:
donde:
Esta etapa abarca la obtención de los diferentes listados de
Q representa la potencia térmica que entra al modelo, h la resultados, su procesamiento, análisis y graficado.
eficiencia térmica de la fuente de calor (arco), Is la corriente
de soldeo empleada (A) y Va el voltaje de arco (V). Metodología para el análisis estructural

A partir de esta potencia térmica se determina el flujo de • Etapa de preprocesamiento:

calor volumétrico (FCv) que se introduce al modelo:
1. Se parte del modelo geométrico tridimensional mallado,
de la corrida térmica.
FCv = Q / Velem (W/m3) (10)
2. Se declaran como propiedades termodependientes del
donde: material, desde la temperatura ambiente hasta la tempe-
ratura de sólido (o la temperatura de corte), el módulo
Velem. representa el volumen del elemento finito que se en- de elasticidad, el coeficiente promedio (secante) de di-
17
 latación térmica lineal y el límite de fluencia del mate- nes de desplazamiento en los nodos y en el área de sime-
rial. Muchos autores, asumen el coeficiente de Poisson tría (en caso de que únicamente se modele la mitad de
con su valor constante a temperatura ambiente, aunque la pieza aprovechando la propiedad de simetría). Si se
de poseer datos confiables, se puede también declarar está modelando la soldadura de chapas en estado libre,
termodependiente (Cañas et al. 1996; Zhu et al. 2002; se debe garantizar que estas restricciones sólo eviten el
Bezerra, 2006). movimiento del cuerpo como un sólido rígido, permi-
tiendo su libre dilatación en los tres ejes coordenados.
El modelo de material implementado en ANSYS que De esta manera se evita el consiguiente error en la etapa
se debe emplear en estos casos es el de plasticidad con de solución, o que se alteren los campos de deformacio-
razón independiente que asume que las deformaciones nes y tensiones, motivo de estudio. Si, por otro lado, se
plásticas se desarrollan instantáneamente y son inde- está modelando la soldadura de piezas con otras restric-
pendientes del tiempo lo que se ajusta al fenómeno estu- ciones, como es el caso de una soldadura de reparación,
diado (ANSYS Inc. 2004; Bezerra, 2006). Luego dentro entonces éstas se declararán de manera que el modelo
de estos modelos se debe decidir cuál es más conve- se aproxime en su comportamiento al cuerpo físico real.
niente para describir el comportamiento elastoplástico
del metal, si el bilineal o el multilineal, y por otro lado, • Etapa de solución:
si se asume sin o con endurecimiento por deformación,
cinemático o isotrópico. El análisis estructural se realiza mediante una serie de
análisis de equilibrio estático. Se parte leyendo las tem-
Se considera que el modelo bilineal conduce a resulta- peraturas nodales de cada paso de tiempo, generadas
dos precisos y para definirlo sólo se necesita declarar la por la corrida térmica y se aplican como cargas nodales
variación del límite de fluencia, el módulo de elasticidad en la corrida estructural. De esta manera cada análisis
y el módulo tangente en relación con la temperatura. Si estructural estático parte del estado final de deforma-
se desprecia el endurecimiento entonces simplemente se ciones y tensiones del análisis previamente concluido
declara el módulo tangente como constante y nulo. Se hasta finalizar. En este análisis no es necesario efectuar
selecciona la superficie de fluencia de von Mises, como una integración en el tiempo debido a que fue estable-
criterio de comportamiento mecánico del material, en lo cido que las deformaciones plásticas son deformaciones
que coinciden prácticamente todos los autores. irreversibles independientes del tiempo (Bezerra, 2006;
ANSYS Inc. 2004).
Si se modela una soldadura de recargue o de ranura de
bordes rectos (sin preparación de bordes), como ya se 1. Se declara como tipo de análisis el estático.
analizó, no se debe emplear la técnica de activación -
desactivación de elementos y entonces se declara como 2. Se emplea el método iterativo de solución de ecuacio-
propiedad del material de toda la unión una temperatura nes no lineales Full Newton Raphson, que constituye el
de referencia igual a la inicial de la pieza. Por otro lado método no lineal estándar de ANSYS para problemas de
si se trata de una soldadura de ranura con preparación de plasticidad (ANSYS Inc. 2004).
bordes, o una de filete, se debe emplear dicha técnica de
activación-desactivación, y se declarará como propie- 3. Se establece un valor de convergencia de fuerza de la
dad de los elementos finitos de la soldadura una tempe- solución, para análisis no lineal de 1 N (Bezerra, 2006).
ratura de referencia igual a la temperatura de sólido de
la aleación, mientras que para el resto de los elementos 4. Se selecciona el avance de tiempo automático (automa-
que conforman el metal base se declara una temperatura tic time stepping).
de referencia igual a la inicial de la pieza.
5. Si se va a modelar el depósito de metal, entonces es ne-
3. Se declara como nuevo tipo de elemento uno con des- cesario desactivar todos los elementos que forman parte
plazamiento en los tres ejes coordenados como grados del cordón antes de aplicar cualquier carga. Si las piezas
de libertad, en este caso se prefiere el SOLID45, de 8 se sueldan con un punteado previo, entonces no se des-
nodos (ANSYS Inc. 2004), utilizado con éxito por otros activan los elementos que forman parte del punteado.
autores (Bezerra, 2006; Hemanth, 2006; Moraitis et al. En esta corrida no se deben restringir los nodos que no
2007). Luego se convierte el tipo de elemento del mo- están ligados a ningún elemento activo, para evitar que
delo mallado de SOLID70 a SOLID45. De esta manera el elemento se deforme en exceso al activarse (ANSYS
se garantiza que geométricamente el mallado no cambie Inc. 2004).
y que la posición y numeración de nodos permanezca
idéntica, lo que es imprescindible a la hora de aplicar las 6. Se declara que la aplicación de las cargas nodales de
cargas de temperatura. temperatura sea en forma de rampa RAMPED, de
modo que para la solución el programa asignará en cada
4. Se declaran, como condiciones de frontera, las restriccio- subpaso de carga un ascenso o descenso progresivo has-
18
 ta el valor final a aplicar, lo que facilita la convergencia mientos y deformaciones.
de la solución.
• Etapa de postprocesamiento:
7. Dentro de un lazo *DO - *ENDDO, en caso de empleo
de la técnica de activación – desactivación de elemen- Es similar a la corrida térmica.
tos, activar los elementos de la soldadura en el momento
oportuno, sobre el que hay discrepancias en la literatura. PROPIEDADES FÍSICAS Y MECÁNICAS DE LA
Se considera que el momento es cuando la fuente de ca- ALEACIÓN 5083 O
lor se encuentra sobre los elementos y la temperatura de
estos sobrepasa la temperatura de sólido del metal. Es- Una gran problemática durante los trabajos de modelación
tos elementos se activan con una deformación térmica de soldadura consiste en establecer, con la mayor precisión
nula, ya que previamente se declaró como temperatura posible, las propiedades mecánicas y físicas termodepen-
de referencia de estos la temperatura de sólido, lo que dientes del material en cuestión. En muchas ocasiones esta
se corresponde con el fenómeno físico real de solidifi- información no está disponible de forma completa, con las
cación del metal de aporte que se deposita a una tempe- propiedades termodependientes desde la temperatura am-
ratura cercana a la de fusión. biente hasta la de fusión del material.

8. Dentro del mismo lazo, mediante el comando LDREAD, La propiedad clave en la corrida térmica es la conductivi-
se programa que ANSYS lea los valores nodales de tem- dad térmica y en este trabajo se tomaron los valores ofreci-
peratura de cada paso de tiempo de la corrida térmica, dos por el código ASME para la aleación 5083, extendidos
almacenados en un fichero con extensión RTH y los según línea de tendencia hasta la temperatura de 580 ºC
aplique como cargas en esta corrida estructural. (temperatura de sólido), lo que concuerda con los valores
brindados por el Eurocódigo para las aleaciones de alu-
9. En el interior del mismo lazo se seleccionan todos los minio al magnesio (ASME, 2002 y Comité Técnico AEN/
nodos cuyo valor de temperatura sobrepasa la tempera- CTN, 2000). La conductividad se duplicó para cualquier
tura de sólido (o de corte) y a los mismos se les aplica temperatura por encima de la temperatura de líquido de la
como carga dicha temperatura, que sustituye la anterior aleación para simular el efecto físico de la transferencia de
almacenada en el fichero RTH, garantizándose que nin- calor por convección que ocurre en el baño de soldadura.
gún nodo se someta a temperaturas superiores a este
valor de corte. Es reconocido que la propiedad clave en la corrida estruc-
tural es el límite de fluencia y los valores empleados se ob-
De esta manera se evita la realización de cálculos iterati- tuvieron a partir de los datos que brinda el ASM Handbook
vos para valores por encima de esta temperatura y cuan- (1990), para la aleación 5083 O (figura 2), extendiéndola
tificar falsas deformaciones plásticas en los elementos según Zhu (2002). El valor mínimo se situó en 6.00E+06
de la soldadura, se hace innecesario también el modelar Pa, y no en un valor nulo, para evitar problemas desde el
el cambio de fase sólido-líquido y líquido-sólido y por punto de vista numérico.
tanto la necesidad de estudiar y declarar la variación de
las propiedades termodependientes del metal a tempera- El resto de las propiedades se asumieron básicamente en
turas superiores a la temperatura de sólido, lo que puede base a las propuestas por Zhu (2002) para la aleación 5052
ser complejo, eliminándose posibles fuentes de error. H32, aleación similar de aluminio al magnesio. La tabla 1
resume las propiedades físico mecánicas empleadas en la
10. En caso de empleo del método de activación – desacti- modelación.
vación de elementos, también dentro del lazo se elimi-
nan las cargas de temperatura aplicadas por el comando
anterior a los nodos que no pertenecen a ningún elemen-
to activo, ya que en estos lugares en la unión soldada
física real o no hay material, o el metal se encuentra
fundido por lo que no ocurren termodeformaciones que
sea de interés cuantificar.

11. Solución mediante el comando respectivo que se in-

tercala al final de cada paso de carga para resolver el
problema.

12. Una vez concluido el cálculo se salva todo el modelo

Figura 2. Límite de fluencia termodependiente
elaborado y los resultados con la historia de desplaza-
de la aleación 5083 O.
19
 Tabla 1. Propiedades asumidas en la modelación de la aleación de aluminio 5083 O.

Propiedades Físicas Propiedades Mecánicas

Temp. k c r Temp. E a sf Et
(ºC) (W/mºC) (J/kgºC) (kg/m3) (ºC) (Pa) (1/ºC) (Pa) (Pa)
25 116.2 930 2680 25 7.00E+10 2.178E-05 1.45E+08 1.800E+09
38 118.2 100 6.60E+10 2.345E-05 1.45E+08 1.615E+09
65 122.4 200 6.00E+10 2.448E-05 1.15E+08 1.358E+09
80 958 2669 260 5.40E+10 2.502E-05 7.50E+07 1.220E+09
93 125.9 430 3.15E+10 2.635E-05 6.00E+06 8.000E+08
121 129.0 580 7.00E+08 2.750E-05 6.00E+06 0
149 131.8
177 134.5
180 1006 2651
204 137.0
240 141.9
280 146.5 1053 2634
330 152.2
380 157.9 1101 2616
420 1120
430 163.6
480 169.3 1140 2598
530 175.0
540 1160
580 180.8 1240 2580
640 360
k- conductividad térmica; c- calor específico; ρ- densidad; E- módulo de elasticidad;
a- coeficiente de dilatación térmica secante; sf- tensión de fluencia; Et- módulo tangente.

MODELACIÓN DE SOLDADURA GMAW SOBRE

ALEACIÓN DE ALUMINIO

Con vistas a validar la anterior metodología propuesta, se

modela mediante ANSYS un experimento desarrollado por
Cañas (1996), y retomado por Chao (1999), de soldadura
GMAW de dos chapas a tope, con preparación de bordes Figura 3. Geometría de la unión soldada
en V y mediante pasada única (figura 3). Las chapas son de del experimento de Cañas.
aleación de aluminio al magnesio 5083-O. El régimen de
soldeo empleado es: Is = 170 A, Va= 23,4 V y vs= 11 mm/s. tos), una intermedia con tamaño de 5 * 2,5 * 3,3 mm (2193
Se emplea como eficiencia de la fuente η = 0,64, similar a nodos, 1300 elementos), y una fina con 5 * 1,25 * 2,2 mm
la usada por ambos autores. Cañas determinó experimen- (4080 nodos, 2800 elementos) (figura 4).
talmente las tensiones residuales en la sección media de la
unión, las que se consideran ideales para validar los resul- En la corrida térmica se declara 25 ºC como temperatura
tados del modelo. inicial de la placa y una pérdida de calor por convección
a través de todas las áreas, con un coeficiente de película
En este caso se modela sólo la mitad de la unión aprove- constante de 15 W/(m2 * ºC), similar al empleado por Ar-
chando la simetría de la misma. Se emplearon tres tipos de mentani (2006) y Mahapatra (2006), y una temperatura del
mallas, una gruesa de partida con un tamaño de elemento en medio de 25 ºC; sólo se exceptúa el área de simetría del mo-
la soldadura de 5 * 5 * 6,6 mm (2244 nodos, 1050 elemen- delo que se comporta como adiabática. Por tratarse de una
20
 Figura 4. Modelos geométricos: a) malla gruesa, b) malla intermedia, c) malla fina.

soldadura de ranura con preparación de bordes se decidió

introducir el calor como flujo volumétrico, con un valor de
12858181818 W/m3 por cada elemento de la misma, el que
actúa durante un tiempo de 0,4545 s. Luego de aplicado el
calor al modelo durante un tiempo de 22.72 s, que represen-
ta el tiempo de desplazamiento de la fuente, se permite su
total enfriamiento hasta los 800 s.

Como propiedades del material se emplearon en ambas co- Figura 5. Restricciones de movimiento
rridas térmica y mecánica las propiedades termodependien- impuestas al modelo mecánico.
tes mostradas en la tabla 1 y en la figura 2.
RESULTADOS Y DISCUSIÓN

Como condición de frontera en la corrida mecánica se de- En la figura 6 se muestra a modo de ejemplo el campo tér-
claró una restricción de simetría en el área de simetría del mico en la superficie superior de la unión, obtenido para
modelo y restricciones adicionales de movimiento en dos un tiempo de 17,7 s desde el inicio de la soldadura, con
nodos (figura 5), de manera que se simulara la soldadura movimiento de la fuente de derecha a izquierda y mediante
de la unión en estado libre. Con un fin comparativo se em- el modelo de malla intermedia. Se observan las isotermas
plearon las distintas variantes posibles de comportamiento alargadas características del proceso de soldadura, con un
elastoplástico bilineal del material. pequeño vacío delante donde los elementos del cordón es-
tán desactivados, por no existir allí metal de aporte en el
proceso físico real.
21
 Figura 6. Campo de temperaturas para una tiempo t = 17,7 s en modelo con malla intermedia.

En todos los modelos estructurales, cuyos resultados de y se ajustan a los valores experimentales de Cañas (1996),
tensiones residuales longitudinales y transversales a la sol- por lo que se considera que en ambos tipos de modelos la
dadura se muestran en las figuras 7, 8 y 9, se activaron los solución prácticamente converge y por tanto es innecesario
elementos de la soldadura en el paso de tiempo en que sobre un refinamiento mayor, que sólo aumentaría los tiempos de
ellos se encuentra la fuente de calor, donde sus nodos pre- cálculo; el error relativo entre ambos modelos en la tensión
sentan temperaturas superiores a la temperatura de fusión. longitudinal sobre el eje de la soldadura es inferior al 10%.

Los modelos geométricos, empleando los tres tipos de ma- Aunque los tres modelos de material (endurecimiento ci-
llas, entregaron resultados bastante semejantes para cada nemático, isotrópico y sin endurecimiento) se ajustan a los
caso de modelo de material (endurecimiento cinemática, datos experimentales medidos por Cañas (1996), se con-
isotrópico y sin endurecimiento). Sólo la malla gruesa falla sidera que los modelos con endurecimiento cinemático y
por el salto que presentan los gráficos de tensión residual sin endurecimiento son los que más fielmente describen el
longitudinal y transversal en la zona de la soldadura, lo que fenómeno, tomando en consideración que el otro modelo
es lógico porque en esta malla la mitad de la sección trans- se aleja algo más de los datos experimentales de la tensión
versal de la soldadura está representada por un único ele- longitudinal en la zona de compresión. Por tanto para esta
mento finito (figura 4-a). Las mallas intermedia y fina si co- soldadura GMAW de pasada única sobre aleación de alu-
inciden prácticamente a todo lo ancho de la unión soldada minio – magnesio se considera que ambos modelos son
adecuados.

Figura 7. Tensiones residuales para modelo de malla gruesa.

22
 Figura 8. Tensiones residuales para modelo de malla intermedia.

Figura 9. Tensiones residuales para modelo de malla fina.

Todos los modelos elaborados produjeron resultados simi- Las isobaras mostradas en las figuras 10 y 11 confirman
lares a los reportados en la literatura para la tensión residual el carácter completamente tridimensional de los campos de
longitudinal, con los valores máximos de tracción próximos tensiones residuales propios del proceso de soldadura. En la
al límite de fluencia del material en la zona cercana a la figura 10 se observa que las tensiones longitudinales mayo-
soldadura altamente calentada y de compresión en la zona res en valor absoluto son de tracción y actúan en una amplia
más alejada de ésta (Frolov, 1988). zona media de la longitud de la soldadura, con valor aproxi-
mado de 166 MPa, mientras que en la figura 11 las mayores
En la figura 9 se muestran los resultados de las tensiones tensiones transversales son de compresión y se encuentran
residuales que produjo un modelo estructural en el que los en las reducidas zonas de inicio y fin de la soldadura, con
elementos de la soldadura se activaron de manera que siem- un valor de casi 200 MPa por lo que son las más críticas, lo
pre existieran dos elementos activos (10 mm) por delante que se corresponde con lo reportado en la literatura (Chao
del elemento sobre el que se encuentra la fuente de calor. et al. 1999). Se considera, por tanto, que ningún modelo
En este caso, y en otros estudiados que aquí no se muestran, bidimensional es capaz de registrar estas variaciones de
se produjo un sobredimensionado de la tensión transversal los niveles de tensiones en las secciones transversales de
respecto a los valores experimentales por lo que se desechó la unión soldada.
esta variante. La causa de esto es que los elementos activa-
dos delante de la fuente de calor provocan un mayor grado
de restricción a la libre dilatación y contracción del metal
en el sentido transversal, lo que se refleja en estos valores
excedidos.

23
 Figura 10. Campo de tensiones residuales longitudinales obtenido mediante tamaño de malla intermedio
con modelo de material sin endurecimiento (valores de tensiones en Pa).

Figura 11. Campo de tensiones residuales transversales obtenido mediante tamaño de malla intermedio
con modelo de material sin endurecimiento (valores de tensiones en Pa).

CONCLUSIONES Activar los elementos de la soldadura, tanto en la corrida

térmica como en la estructural, en el paso de tiempo en el
Es eficiente computacionalmente y conduce a resultados que la fuente de calor se encuentra sobre ellos, conjugado
precisos el emplear modelos geométricos con un mallado con el empleo de una temperatura de referencia de los ele-
fino únicamente en la zona altamente calentada, donde los mentos de la soldadura igual a la temperatura de fusión del
gradientes térmicos y de tensiones son altos, y que se hace metal, conduce a resultados precisos en el cómputo de la
mucho más basto en la medida que se aleja de ésta. Siempre historia térmica, tensiones y deformaciones de soldadura,
el mínimo posible de nodos y elementos garantiza un tiem- evitándose falsas deformaciones elasto-plásticas en los ele-
po de cálculo más eficiente. mentos que conforman la soldadura. Por otro lado activar
los elementos de la soldadura delante de la fuente de ca-
24
lor provoca un cierto sobredimensionamiento de la tensión Armentani, E., Esposito, R., Sepe, R. (2006). Valutazione
transversal. numerica delle tensioni residue in processi di saldatura
GMAW. Associazione Italiana per l’Analisi delle Solle-
El empleo del modelo de material con comportamiento citazioni (AIAS). XXXV Convegno Nazionale.13-16
elastoplástico del tipo bilineal conduce a resultados preci- Settembre 2006, Università Politecnica delle Marche.
sos, y para definirlo sólo se necesita declarar la variación
del límite de fluencia, el módulo de elasticidad y el módulo Bezerra, A.C. (2006). Simulação Numérica da Soldagem
tangente en relación con la temperatura. Aunque los tres com Aplicação à Caracterização do Comportamento
modelos de material (endurecimiento cinemático, isotró- Dinâmico de Estruturas Soldadas. Tesis de Doctorado
pico y sin endurecimiento) se ajustan a los datos experi- no publicada, Universidad Federal de Uberlândia, Uber-
mentales, se considera que los modelos con endurecimiento lândia. Brasil. p. 124.
cinemático y sin endurecimiento son los que más fielmente
describen, en este caso, el fenómeno termodeformacional Budgell, P. (2007). Tips and ANSYS Tricks. Recuperado
de la soldadura GMAW de aleaciones aluminio – magnesio.
en Junio de 2007. http://www3.sympatico.ca/peter_bud-
gell/ANSYS_tips.html.
El empleo de una temperatura de corte en la corrida estruc-
tural igual a la temperatura de fusión de la aleación evita
Camilleri, D., Gray, G.F. (2005). Computationally effi-
cuantificar falsas deformaciones elastoplásticas a tempera-
cient welding distortion simulation techniques. Mode-
turas superiores, cuando en el proceso físico real el metal
lling Simul. Mater. Sci. Eng. (13); pp. 1365–1382.
se encuentra en fase líquida o sólida - líquida, evitándose
también la necesidad de ser exquisitos en ajustar las propie-
Cañas, J., Picón, R., París, F.A. (1996). A simplified nume-
dades a temperaturas superiores, lo que pudiera representar
alguna dificultad. rical analysis of residual stresses in aluminum welded
plates. Computers & Structures. Vol. 58. No. I; pp. 59-
Se considera que solamente los modelos tridimensionales 69.
son capaces de reproducir con precisión el comportamiento
tensional de las uniones soldadas, y permiten determinar Chao, Y. & Qi, X. (1999). Three dimensional modeling of
el valor máximo de tensión transversal de compresión que gas metal arc welding process. Transaction of NAMRI/
se produce en las zonas al inicio y fin de la soldadura, lo SME. Volumen XXVII. pp. 117-122.
que es prácticamente inalcanzable para los modelos en dos
dimensiones. Comité Técnico AEN/CTN 140 Eurocódigos Estructu-
rales (1999). Eurocódigo 9: Proyecto de estructuras de

La metodología de modelación aquí expuesta integra una aluminio. Parte 1-2: Reglas generales. Proyecto de es-
información que se encuentra dispersa en la literatura, pro- tructuras sometidas a fuego. ARNOR 2000. pp. 25-26.
poniendo variantes que en ocasiones superan lo publicado
en trabajos precedentes, como la modelación del desplaza- Francis, J.D. (2002). Welding Simulations of Aluminum
miento de la fuente de calor mediante un sencillo lazo *do- Alloy Joints by Finite Element Analysis. Thesis presen-
*enddo dentro del que se introduce el comando que aplica ted in option to the degree of Master of Science in Ae-
la carga térmica, en lugar del empleo de curvas de tiempo rospace Engineering. Faculty of the Virginia Polytech-
o funciones de carga, o el momento de activación de los nic Institute and State University. p. 178. Recuperado en
elementos de la soldadura, constituyendo un aporte para el enero de 2007 http://scholar.lib.vt.edu/theses/available/
desarrollo de esta temática. etd-05062002- 00056/unrestricted/JFThesis4_12_02.
pdf.
Con este trabajo se confirma que el Software ANSYS Mul-
tiphysics versión 9.0 de uso general puede ser empleado Frolov, V.V. (1988). Teoriya svarochnij protsessov. Vissha-
satisfactoriamente en la modelación de la historia térmica, ya Shkola. Moscú; pp. 425-429.
tensiones y deformaciones de soldaduras, brindando resul-
tados tan fiables como cualquier código de elementos fini-
García, Y., Burgos, J. (2007). Procedimiento para la obten-
tos especializado en este proceso.
ción mediante MEF de campos térmicos, deformacio-
nes y tensiones residuales en uniones soldadas. Revista
REFERENCIAS
Técnica de Ingeniería, Universidad de Zulia, Venezuela,
Vol. 30, Nº 1 pp. 13-22.
ANSYS INC. Tutorial de ANSYS Multiphysics versión 9.0.
25
He, Y., Shouju, L., Yingxi, L. (2007). Finite Element Simu- Pozo, J. A., Ramos, F., García, Y., Cruz, A., Díaz, E.,
lation for Residual Stresses in Welding Process. Engi- Duffus, A. (2008). Modelación Mediante MEF de Ten-
neering Materials Vols. 353-358; pp. 1915-1918. Recu- siones y Deformaciones en Soldadura GMAW sobre
perado en enero de 2008. http://www.scientific.net. Aleación de Aluminio 5052 H32. Soldagem & Insp. São
Paulo, Vol. 13, No. 4, pp. 329-337. Recuperado en enero
Hemanth, K. (2006). Three dimensional thermal modeling de 2008. http://www.abs-soldagem.org.br/s&i/.
of friction stir processing. Tesis para el grado de Master
of Science. Florida State University. College of Engi- Shah, M. (2005). Structural integrity and fatigue crack pro-
neering. p. 59. pagation life assessment of welded and weld-repaired
structures. A Dissertation Submitted to the Graduate
Jerome, S., Kumaresh, S.P., Ravisankar, B., Pandey, K. S. Faculty of the Louisiana State University and Agricul-
(2007). Prediction of welding residual stresses in butt tural and Mechanical College in partial fulfillment of the
welded plate by finite element analysis. International requirements for the degree of Doctor of Philosophy in
Conference on Advanced Materials and Composites The Department of Mechanical Engineering, pp. 164.
(ICAMC-2007). Oct 24-26. pp. 553-558.
The American Society of Mechanical Engineers-ASME
Kim, K. S., Choi, S. B., Lee, J. H., Park, S.M. (2006). A Stu- (2002). ASME Boiler and pressure Vessel Code. ASME.
dy on Measurement of Welding Residual Stress Using Sección II, Parte D. New York, N.Y, 2002. p. 668.
ESPI System. Engineering Materials Vols. 324-325;
pp. 859-862. Recuperado en enero de 2008 http://www. The American Society of Metals-ASM (1990). Properties
scientific.net. and selection: nonferrous alloys and special purpose
materials. Metals Handbook. Vol 2. 10ma Edición. 1990
Leal, M., Pozo, J. A., García, Y. (2007). Modelación por USA International.
MEF de los Campos Térmicos y de las Tensiones Resi-
duales de Soldadura en Aleaciones de Aluminio. Univer- Zhu, X.K., Chao, Y.J. (2002). Effects of temperature-
sidad Central “Marta Abreu” de Las Villas, Facultad de dependent material properties on welding simulation.
Ingeniería Mecánica, Centro de Investigación de la Sol- Computers and Structures 80 (2002); pp. 967-976.
dadura (CIS). Trabajo de Diploma no publicado. p. 98.

Mahapatra, M.M., Datta, G.L., Pradhan, B., Mandalb,

N.R. (2006). Three-dimensional finite element analysis
to predict the effects of SAW process parameters on tem-
perature distribution and angular distortions in single-
pass butt joints with top and bottom reinforcements.
International Journal of Pressure Vessels and Piping 83;
pp. 721–729.

Meo, M., Vignjevic, R. (2007). Welding simulation using

FEA. Recuperado en diciembre de 2007 http://www.
eng.upm.edu.my/~kaa/WEC/WECPapers-FinalVersion/
AE16.doc.

Moraitis, G.A., Labeas, G.N. (2007). Residual stress and

distortion calculation of laser beam welding for alumi-
num lap joints. J. Mater. Process. Tech. Recuperado en
enero de 2008. http://www.elsevier.com/locate/jmatpro-
tec.

Ramos, F., Duffus, A., Rodríguez, M., Díaz, E., Pozo, J. A.,
Miguel, J.V. (2007). Modelo de elementos finitos para
determinar el campo de temperatura en una costura de
filete. Ingeniería, Revista Académica de la FI-UADY,
11-3, ISSN: 1665-529X; pp. 33-42
26