Open Journal of Applied Sciences, 2012, 2, 47-53
doi:10.4236/ojapps.2012.21005 Published Online March 2012 (http://www.SciRP.org/journal/ojapps)
Numerical Simulation for Shaping Feature of Molten Pool
in Twin-Ar c Submerged Arc Welding
Kuanfang He*, Jun Chen, Siwen Xiao
Hunan Provincial Key Laboratory of Health Maintenance for Mechanical Equipment, Hunan University of
Science and Technology, Xiantan, China
Email: *hkf791113@163.com
Received January 4, 2012; revised February 6, 2012; accepted February 18, 2012
ABSTRACT
The notion of ratio of width to length is proposed to describe the shaping feature of molten pool of twin-Arc submerged
arc welding accurately, and analyze the law of molten pool variation and weld formation. The temperature field finite
element numerical simulation model of twin arc movement is established. The loading form of twin-arc with double
ellipsoid heat source is discussed. The molten pool temperature field of twin-arc submerged arc welding is calculated
and analyzed under different process parameters. The law of molten pool characteristics influenced by the welding
speed, current and voltage of twin-arc submerged arc welding parameters is analyzed. The relation between shaping
feature of molten pool and weld formation is discussed according to the ratio of width to length. The results manifested
that the width to length ratio of weld pool decreases with the improvement of welding speed, which result in generation
of weld defects. The width to length ratio of weld pool is increased by adjusting the proportion of the current, voltage
and the distance of the two arcs, which avoids the generation of weld defects.
Keywords: Twin-Arc Submerged Arc Welding; Shaping Feature; Ratio of Width to Length; Numerical Simulation
1. Introduction
Twin-arc submerged arc welding technology is very com-
plex, especially the heat characteristics is substantially
different, so it affects the shaping feature of molten pool
[1,2]. In the process of twin-arc submerged arc welding,
technological parameters are absolutely critical for the
molten pool forming, which directly affect the heat, mo-
mentum and quality of the molten pool, and has a further
decision on weld formation. As a result, it is significant for
the studies on the mechanism of weld formation to re-
search the molten pool morphological characters of twin-
arc submerged arc welding under different technological
parameters. With the development of computer technology,
numerical simulation is a powerful method to obtain quan-
titative understanding of welding process. The numerical
simulation of molten pool temperature field is an impor-
tant method to research the mechanism of weld formation
[3,4]. So it will provide essential data and theoretical
foundation for the optimization and digital control of
welding process.
Among the studies on loading form of arc with heat
source, first, an analytic model was most usually used as
point or line and face heat source was proposed [5]. The
point heat source was put forward initially to simulate
surface melt in welding process [6]. The units treated as
heat source in melting zone was used to design calcula-
tion program of welding heat transfer, difference method
was adopted in the program considering relationship be-
tween material physical property and temperature [7]. On
this basis, a Gaussian distribution mathematical model for
weld heat source in space was presented, then other mo-
dels based on Gaussion distribution such as hemispherical
distribution, ellipsoidal distribution, double ellipsoidal
distribution were used to calculate the temperature field
[8,9]. These heat sources mentioned above have been
applied in temperature field numerical simulation of vari-
ous welding situation. Recently, based on the double el-
lipsoid heat source model, a heat source model for twin-
wire welding was proposed [10], which was adopted to do
our work in this paper.
Recently, numerical simulation on welding had been
carried out. Dong et al. [11] have proposed the numerical
analysis of welding residual stress for the multi-pass
weldment, this paper investigates the two moving source
of ramp and double-ellipsoid in prediction of residual
stresses in a multi-pass-welded plate. In order to master
law of molten pool variation and weld formation, nu-
merical simulation of the molten pool temperature field
had been performed by finite element method.
A transient numerical model was established to provide
detailed insight about the nature of heat transfer and fluid
flow during laser hybrid stationary spot welding alumi-
*Corresponding autho
r
.
Copyright © 2012 SciRes. OJAppS
K. F. HE ET AL.
48
num alloy [12]. The dynamic development of weld pool
geometry during TIG welding was analyzed numerically,
and the effect of arc moving on the weld pool geometry
was discussed [13]. Deep penetration laser welding tem-
perature field of 5A06 aluminum alloy canister structure
was simulated using the surface body combination heat
source model by ANSYS, which was made up of Gauss
surface heat source model and Gauss revolved body heat
source model, radiation and conduction were all consi-
dered during the simulation process [14]. Aval et al. [15]
studied the theoretical and experimental of microstruc-
tures and weld pool geometry during GTAW of 304
stainless steel, temperature field and weld pool geometry
during gas tungsten arc welding of 304 stainless steel
were predicted by solving the governing equations of heat
transfer and fluid flow under quasi-steady state conditions.
Their research focus on the analysis of the shaping feature
of molten pool and the weld formation influenced by
process parameters. Taking into account the multi-field
coupling, the molten pool temperature field affected by
magnetic and other arc has been studied. Lin et al. [16] has
introduced the multicoupled analysis function of ANSYS
to analyze a moving gas tungsten arc weld pool with an
external longitudinal magnetic field applied, the distribu-
tions of current density and magnetic field, as well as fluid
flow and heat transfer in a moving weld pool, were sys-
tematically studied and investigated to understand and
reveal the effect of an external longitudinal magnetic field
on liquid metal in a moving GTA weld pool. Zhu et al. [17]
analyzed the projection welding on auto-body sheet metal
using a coupled finite element method. A comprehensive
finite element method employing a subroutine to link up
submodules of commercial code ANSYS was proposed to
perform analysis of projection welding in quantitative de-
tail. At present, based on analysis of submerged arc weld-
ing arc heat source model and droplet heat inputting uni-
form distribution, ANSYS parametric design language was
applied to develop sub-program for three-dimensional tem-
perature field numerical simulation of twin-arc high- speed
submerged arc welding process [18].
Based on the previous research, the notion of ratio of
width to length was proposed to describe the shaping fea-
ture of molten pool of twin-Arc submerged arc welding.
We analyzed the law of molten pool variation and weld
formation. The temperature field finite element numerical
simulation model of twin arc movement was established
in Section 2. Section 3 calculated the molten pool tem-
perature field of twin-arc submerged arc welding under
different technological parameters, and analyzed the law
of molten pool characteristics influenced by the welding
speed, current and voltage of twin-arc submerged arc
welding parameters. According to the ratio of width to
length, the relation between shaping feature of molten
pool and weld formation were further discussed in Sec-
tion 4. The conclusions were summarized in Section 5.
2. Numerical Simulation Modeling
Sketch of practical model was shown in Figure 1. Resur-
facing welding of medium plate surface was employed to
do simulation analysis, material of weldment was low-
carbon steel Q235, physical dimension of the workpiece
was 200 mm × 100 mm × 20 mm, and weld seam was
located on center line of x-y plane of the weldment. Arc
center moves along the y axis. Thermal physical perfor-
mance parameters were thermal coefficient (W/m·˚C),
coefficient of heat transfer (W/m2·˚C), density (Kg/m3),
specific heat (J/Kg·˚C), melting point (˚C) and initial
temperature of workpiece (˚C). In this paper, initial tem-
perature of workpiece was room temperature of 20˚C. The
material was isotropical and change with temperature
range, each thermal physical performance parameters was
calculated according to the following temperature func-
ion [19]: thermal coefficient: , t
2
54.3 0.000042 T
Welding Arc
Wire
Welding Direction
Torch
1
f
I
2f
I
1
f
U
2f
U
V
l
Weldpool
Fusion ZoneWorkpiece
o
y
z
x
x
y
Figure 1. Sketch of practic a l model.
Copyright © 2012 SciRes. OJAppS
K. F. HE ET AL. 49
specific heat: , density:
2
410 0.63T
0.2625T7800
, coefficient of heat transfer:
33.5
. Interpolation was conducted on thermal physi-
cal performance parameters in the unknown temperature
range [20].
2.1. Finite Element Model and Meshing
The three-dimensional eight-node element, SOLID70, was
employed as the finite element in the generation of finite
element model. Division density of elements were associ-
ated with the calculation, if the mesh was too coarse, the
results probably contain serious errors, the higher the ele-
ment density was, the more accurate the results was, but if
the mesh was too fine, it will cost too much computing
time and waste computer resources. As a result, fine mesh
was used in meshing of the welding seam to obtain a
higher element density and assure the computing accurate.
Coarse mesh was employed in the places far away from
the welding heat source, because the temperature gradient
was small. So it can reconcile both computing accurate
and speed. The meshing of finite element model shows as
Figure 2. In addition, because of the geometrical shape
and load distribution were symmetrical along the center of
welding line, half of the model was computed in order to
save computing time.
2.2. Boundary Conditions
Analysis of welding temperature field was a typical pro-
blem of non-transient heat conduction, the three-dimen-
sional temperature field controlling equation could be
given as:

  
 


 
   
 

TTT T
c
tx xy yzz

Q
(1)
where is internal heat source intensity in
solving domain V, T is function of temperature field dis-
tribution, λ is thermal coefficient, ρ is density and c is
specific heat of the material. The expression above is a
,,,Qxyzt
universal definite equation, in order to obtain certain so-
lution, conditions of determining solution were needed,
that were boundary conditions and initial conditions of
differential equation. The calculations of welding tem-
perature field often have these boundary conditions:
1) The first boundary condition, boundary tempera-
ture is given as:
,,,



xyzs
TTT
nnnTxy
xyz

zt
(2)
2) The second boundary condition, heat flux density
distribution was given as:
,,,



xyzx
TTT
nnnqxy
xyz

zt
(3)
3) The third boundary condition, heat exchange be-
tween boundary objects and ambient medium was given
as:

x
yz
TTT
nnnT
xyz
 


 s
T
(4)
In expressions (2), (3) and (4),
x
q is external heat in-
put per unit area, β is the surface heat transfer coefficient,
T
is the surface temperature,
s
T is ambient medium
temperature, ,,
yz
nn n is direction cosine of external
normal line respectively.
2.3. Heat Source Model
Gaussian heat flux distribution function was commonly
used in arc welding temperature field simulation currently,
which was applicable to situations of small specification
and arc momentum effect, such as manual arc welding
and argon tungsten-arc welding [6]. In submerged arc weld-
ing process, the specification was relatively large, it was a
welding method of large arc momentum effect, and the
double ellipsoidal heat source distribution function was
more suitable than the Gaussian heat flux distribution
unction in numerical simulation of temperature. Heat of f
x
y
V
o
z
Figure 2. Finite element griddling.
Copyright © 2012 SciRes. OJAppS
K. F. HE ET AL.
50
the model was mainly concentrated in the double ellipsoid
above the workpiece, the mathematical expression of dou-
ble ellipsoid heat source model of a moving single-arc heat
source could be given as:

222
222
1
1
222
222
2
2
63 333
exp 0
ππ
,,
63 333
exp 0
ππ

 


 
f
r
fQ xyzx
bac
ab c
qxyz
fQ xyzx
bac
ab c
(5)
where is heat flow density in double ellip-
soidal Gauss body of workpiece with its unit of
,,qxyz
2
Wm.
Q is effective power of arc provided by the welding
power source with unit of W
f
f
and r
f
are energy
distribution coefficients of first and latter half of the
model, which satisfy the relations: . 1, 2,
a and c are the intercept of heat source flow field on axis
of x, y and z, which represents shape of the ellipsoid.
2
fr
ff b b
Heat flow density of each node was calculated out ac-
cording to the given welding parameters and double el-
lipsoidal heat source distribution functions, and exerted
on the chosen node. In the twin-arc submerged arc weld-
ing, two arc heat sources were obtained according to the
given welding parameters and expression (5), then the
two heat sources were exerted following the spatial dis-
tribution of the actual location of the arc, repeatedly doing
this to the new nodes along the welding seam when arc
center moving, and realize exertion of the two moving
heat sources.
3. Modeling Results and Discussions
APDL (ANSYS parametric design language) of was em-
ployed to achieve pre-treatment, loading and boundary
conditions of the workpiece, linear search and automatic
time step were used to improve the computing accuracy
and treatment of element birth and death, and Full New-
ton-Raphson method was used to solution. The parame-
ters of twin-arc submerged arc welding were given in
Table 1, the computing results were shown in Figure 3.
Curve in Figure 3 were outline of weld pool that indi-
cates the melting temperature of steel Q235. It was simi-
lar to the pool surface size from of low carbon steel from
the practical testing.
In Figure 3, it can be found that outline of weld pool
presents a former-small and later-large double-oval shape
along the welding direction on the surface (X-Y) and sec-
tion (Y-Z), this was because the effection of the second
heat source to the first that leads to be wider at the end of
the first heat source, and have a great deepen penetration,
which indicates the twice heating process of base metal.
Geometric parameters presented the weld pool geome-
try mainly including: Maximum width of weld pool Bmax,
weld pool length L, weld pool depth H. Maximum width
of weld pool Bmax was defined as maximum distance be-
tween two weld pool boundary points perpendicular to
the welding direction. Weld pool length L was defined as
the distance between the head and tail of the weld pool.
Weld pool depth H was defined as maximum distance
between two weld pool boundary points parallel to the
welding direction. It can be seen from Figure 3 that the
max weld width and penetration were in different sec-
tions, the max penetration was behind of the max width,
the distance was supposed as X. The welding speed was
supposed as VB, the melting speed along the pool depth
direction was VH, arc heating effecting time in the max
penetration was t, then
H
DV t. Welding speed was
constant, the X was expressed by
Table 1. Calculation conditions.
Welding
current (A)
Welding
voltage (V)
Wire
spacing (mm)
Welding
speed (m/min)
Thermal
efficiency η
If1 If2 Uf1 Uf2
700 600363830 0.96 0.75
L
max
B
max
H
X
Figure 3. Calculated outline and practical surface of weld pool.
Copyright © 2012 SciRes. OJAppS
K. F. HE ET AL. 51

B
BB
HH
V
D
X
VtV D
VV (6)
It can be seen from the expression that X was propor-
tional to the welding speed VB, and inversely proportional
to the melting speed VH along the pool depth direction.
The larger the welding speed VB was, the greater the dis-
tance X was. In high speed welding process, X have a
significant impact on weld pool behavior when V
B was
larger, which makes it difficult to guarantee welding qua-
lity, and also leads to phenomenon of undercut and hump
weld seam.
In order to describe the influence rules of the weld pool
behavior of high-speed welding on weld seam forming,
Ratio of width to length of weld pool was introduced, and
expressed as
, which was the ratio of the max weld
width and weld pool length, the max weld width was ex-
pressed as B, the weld pool length was expressed as L,
then the
can be expressed by
BL
(7)
The purpose of definition of these dimension parameters
was to further analyze the influencing factors of weld pool
behavior of twin-arc high speed submerged arc welding.
4. Discussions
4.1. Effect of Welding Speed on Weld Pool Shape
States of weld pool shape in single-arc and twin-arc weld-
ing when welding speed change were computed according
to the parameters in Table 2 and Table 3, the results
shows in Table 4.
Table 2. Calculation conditions.
Welding
current
(A)
Welding
voltage
(V)
Wire
diameter
(mm)
Plate
thickness
(mm)
Plate
width
(mm)
Welding
speed
(m/min)
600 36 4 12 100 0.48/0.6/0.72
Table 3. Calculation conditions.
Welding
current
(A)
Welding
voltage
(V)
Wire
diameter
(mm)
Plate
thickness
(mm)
Wire
spacing
(mm)
Welding
speed
(m/min)
If1 If2 Uf1 Uf2
600 600 36 38 4 12 30 0.96/1.2/1.44
Table 4. Calculation results.
Arc state Single arc Double arc
Welding speed
(m/min) 0.48 0.6 0.72 0.96 1.2 1.44
Weld width (mm) 16.70 14.90 13.70 17.20 15.6014.20
Weld pool
length (mm) 44.10 45.20 46.00 44.80 45.60 46.40
Ratio of width
to length 0.38 0.33 0.30 0.39 0.340.31
It can be known from Table 4 that the weld width de-
creased while weld pool length increased and the ratio of
width to length of the weld pool decreased gradually with
the increase of welding speed. The ratio of width to
length the weld pool in twin-arc welding was approxi-
mately equal to that of the single arc welding. As a result,
the maximum Ratio of width to length was employed in
this paper to analyze the effect of weld pool stability on
weld seam forming in high speed welding process.
According to the above calculation and analysis, the
twin-arc was arranged along the direction of welding, a
weld pool was formed, detention time of the arc was
lengthened, the critical welding speed of poor weld pro-
ducing was effectively improved, tendency of production
of poor weld seam was reduced. Meanwhile, owing to
the mutual thermal effect of the twin-arc, heat energy
efficiency and filling capacity of weld wire were increased,
and also had a retarding effect on generation of hump in
high speed submerged arc welding.
4.2. Effect of Voltage on Weld Pool
According to the parameters in Table 5, temperature dis-
tribution on the weld pool surface of twin-arc submerged
arc welding in different arc voltage were computed, the
results were shown in Table 5 and Figure 4. It can be
seen from the results that a former small and later large
superimposed double-oval shape was presented with its
axis along the welding direction on the Y-Z surface, be-
cause the twice heating effect of the heat source lead to a
deeper penetration at the beginning of the second heat
source, it indicated the twice heating process of base
metal. On the X-Y surface of weld pool, when the voltage
of the front arc was larger than that of the later arc, the
center of weld pool surface was tightened, this kind of
weld shape characteristic in high speed welding easily
Table 5. Calculation condition sand results.
sequence
number If1 (A) If2 (A) Uf1 (V) Uf2 (V)
Welding
speed
(m/min)
Breadth
length
ratio
(a) 60060038 36 0.6 0.36
(b) 60060036 38 0.6 0.38
X-Y section X-Y section
Z-Y section Z-Y section
(a) (b)
Figure 4. Weld pool contour of different voltage combination.
Copyright © 2012 SciRes. OJAppS
K. F. HE ET AL.
52
lead to generation of weld defects such biting edge and
the hump. When the voltage of the front arc was smaller
than that of the later arc, a former-small and later-large
superimposed double-oval shape was presented with its
axis along the welding direction on the X-Y surface, which
indicated the twice heating process of base metal, and the
calculated results of the ratio of width to length was lager
than the former. Conclusions can be obtained that weld
defects such as hump can be partly resisted when the
voltage of the front arc was slightly less than that of the
later arc along the welding direction, and perfect weld
seam performance can be obtained.
4.3. Effect of Current on Weld Pool
The size and collocation of the twin-arc current were
most important to the temperature field distribution of
weld pool. The front and later arc voltage were selected
as 36 V and 38 V respectively, the total inputting current
was 1300 A, temperature field distribution of weld pool
was calculated according to the parameters in Table 6,
the calculated results were shown in Figure 5.
It can be seen from the calculated results, when the
current ratio of the front and later arc If1:If2 changed from
1:1 to 7:6, the ratio of width to length was relative large.
Conclusions can be obtained that weld defects such as
hump can be partly resisted when the heat inputting of
the front arc was slightly larger than that of the later arc
along the welding direction, and perfect weld seam per-
formance can be obtained.
Table 6. Calculation conditions.
Welding
speed
(m/min)
Wire
diameter
(mm)
Plate
thickness
(mm)
Plate
width
(mm)
Current ratio
If1: If2
0.96 4 12 100 2:8/3:7/4:6/1:1/6:4/7:3/8:2
2:8 3:7 6:7 1:1 7:6 7:3 8:2
0.3 0
0.3 2
0.3 4
0.3 6
0.3 8
0.4 0
0.4 2
熔池宽长比( B /L )
两电弧电流比( I
f1
/I
f2
)
The
r
atio of wi
d
th to length
(B/L)
Figure 5. The ratio of width and length of weld pool in dif-
ferent current of two arc.
4.4. Effect of Different Wire Space on Weld Pool
The space of twin-arc effect on the shape of weld pool
was obvious. If the space between the two wires was zero,
its effect was equivalent to the work of a single arc, this
was an ideal state. If distance was too large, two weld
pools will be formed and independently of each other.
According to the parameters in Table 7, the calculated
results are shown in Tab le 8. It can be seen that the ratio
of width to length of weld pool reduced gradually with
the increasing of space between the two wires.
In order to get the ratio of width to length of the weld
pool and ensure the weld seam performance, space be-
tween the two wires should be closer. If the space be-
tween the two wires was too closer, because the interac-
tion was existed in the twin-arc, it made large interfere-
ence and reduced the arc stability. As a result of this, the
arc stability and effect factors on the weld shaping should
be considered in practice.
5. Conclusions
Temperature field finite element numerical simulation
model of twin arc movement was established. The loading
form of twin-arc with double ellipsoid heat source was
realized by the APDL of ANSYS. The three-dimensional
dynamic simulation of temperature field of plate sub-
merged arc resurfacing welding was conducted. The geo-
metric parameters of weld pool were obtained through
calculation data results. The effect rules of welding pro-
cess on the morphological character of weld pool were
analyzed by introducing the ratio of width to length pro-
posed in this paper. Numerical simulation of twin-arc
submerged arc welding indicated the law of molten pool
variation and weld formation as following. With increase-
ing of the welding speed, the welding width decreased
and the pool length increased gradually, the ratio of width
to length of weld pool reduced gradually, it would raise
the generation tendency of poor weld shaping. When the
two heat sources was arrayed in tandem along the welding
direction, the ratio of width to length of weld pool could
Table 7. Calculation conditions.
Welding
current (A)
Welding
voltage (V)
Wire
diameter
(mm)
Welding
speed
(m/min)
Wire
Space
(mm)
If1 If2 Uf1 Uf2
60050036 38 4 0.96
30/40/
50/60
Table 8. Calculation results.
Wire spacing l (mm) 30 40 50 60
Weld width (mm) 17.20 16.23 15.9114.44
Pool length (mm) 43.70 54.8 65.7176.82
Breadth length ratio 0.400 0.300 0.2420.188
Copyright © 2012 SciRes. OJAppS
K. F. HE ET AL.
Copyright © 2012 SciRes. OJAppS
53
be effectively increased, poor weld shaping was restrained
and higher speed welding was realized.
6. Acknowledgements
Project supported by Hunan Provincial Natural Science
Foundation of China (11JJ2027), National Natural Science
Foundation of China (51005073), Project of Hunan Pro-
vincial Research Scheme (2011GK3052), Scientific Re-
search Fund of Hunan Provincial Education Department
(10C0682), Ph. D Start Fund (E51088), CEEUSRO spe-
cial plan of Hunan province (2010XK6066), Industrial
Cultivation Program of Scientific and Technological
Achievements in Higher Educational Institutions of Hu-
nan Province (10CY008), also from Aid program for
Science and Technology Innovative Research Team in
Higher Educational Institutions of Hunan Province, are
gratefully acknowledged.
REFERENCES
[1] D. J. Tusek, “Mathematical Modeling of Melting Rate in
Twin-Wire Welding,” Journal of Materials Processing
Technology, Vol. 100, No. 1-3, 2000, pp. 250-256.
doi:10.1016/S0924-0136(99)00485-9
[2] H. Q. Li, X. Wang and X. P. Cai, “Developing Status of
SAW at Home and Abroad,” Electric Welding Machine,
Vol. 36, No. 4, 2006, pp. 1-6.
[3] J. Dong, C. L. Mo and Q. C. Shan, “The Simulation of the
Automatic Welding Temperature Field in the Hidden Arc
Butt Joint Welding of 20 Steel Tube,” Transactions of
Shenyang Ligong University, Vol. 26, No. 2, 2007, pp.
38-41.
[4] X. F. Tan and H. Zhang, “Analysis of Temperature and
Stress Fields in Welding Process by Thermo-Elastoplastic
Finite Element Method,” Journal of Plasticity Engineer-
ing, Vol. 11, No. 5, 2004, pp. 71-74.
[5] C. L. Mo, B. N. Qian and X. M. Guo, “The Development
of Models about Welding Heat Sources Calculation,”
Transactions of the China Welding Institution, Vol. 22,
No. 3, 2001, pp. 93-96.
[6] W. F. Zhu, Z. Q. Lin and X. M. Lai, “Numerical Analysis
of Projection Welding on Auto-Body Sheet Metal Using a
Coupled Finite Element Method,” The International Jour-
nal of Advanced Manufacturing Technology, Vol. 28, No.
1-2, 2006, pp. 45-52. doi:10.1007/s00170-004-2336-8
[7] Z, Paley and P. D. Nbbrt, “Computation of Temperatures
in Actual Weld Designs,” Welding Journal, Vol. 54 , No.
11, 1975, pp. 385-392.
[8] T. W. Eager and N. S .Tsai, “Temperature Fields Pro-
duced by Traveling Distributed Heat Sources,” Welding
Journal, Vol. 62, No. 12, 1983, pp. 346-355.
[9] G. A. John, “A New Finite Model for Welding Hear
Source,” Metallurgical Transactions, Vol. 15, No. 2,
1984, pp. 299-305. doi:10.1007/BF02667333
[10] Q. G. Meng, “Heat Source Model For Twin Wire Weld-
ing,” Chinese Journal of Mechanical Engineering, Vol.
41, No. 4, 2005, pp. 110-113.
doi:10.3901/JME.2005.04.110
[11] H. B. Dong, H. K. Chul and Y. C. Seon, “Numerical
Analysis of Welding Residual Stress Using Heat Source
Models for the Multi-Pass Weldment,” Korean Society of
Mechanical Engineers International Journal, Vol. 16, No.
9, 2002, pp. 1054-1064.
[12] Z. G. Gao, Y. X. Wu and J. Huang, “Analysis of Weld
Pool Dynamic during Stationary Laser-MIG Hybrid
Welding,” The International Journal of Advanced Manu-
facturing Technology, Vol. 44, No. 9-10, 2009, pp. 870-
879. doi:10.1007/s00170-008-1896-4
[13] M. Zhao and R. Y. Li, “Numerical Analysis of Dynamic
Variation of Weld Pool Geometry in Fully-Penetrated
TIG Welding,” China Welding, Vol. 17, No. 2, 2008, pp.
47-53.
[14] S. L Chen, Z. S. Wu and C. R. Liu, “Numerical Simula-
tion of Temperature Field in Deep Penetration Laser
Welding of 5A06 Aluminum Cylinder,” China Welding,
Vol. 17, No. 4, 2008, pp. 74-78.
[15] H. J. Aval, A. Farzadi and S. Serajzadeh, “Theoretical
and Experimental Study of Microstructures and Weld
Pool Geometry during GTAW of 304 Stainless Steel,”
The International Journal of Advanced Manufacturing
Technology, Vol. 42, No. 11-12, 2009, pp. 1043-1051.
doi:10.1007/s00170-008-1663-6
[16] Z. Q. Lin, Y. B. Li and Y. S. Wang, “Numerical Analysis
of a Moving Gas Tungsten Arc Weld Pool with an Exter-
nal Longitudinal Magnetic Field Applied,” The Interna-
tional Journal of Advanced Manufacturing Technology,
Vol. 27, No. 3-4, 2005, pp. 288-295.
doi:10.1007/s00170-004-2194-4
[17] W. F. Zhu, Z. Q. Lin and X. M. Lai, “Numerical Analysis
of Projection Welding on Auto-Body Sheet Metal Using a
Coupled Finite Element Method,” The International Jour-
nal of Advanced Manufacturing Technology, Vol. 28, No.
1-2, 2006, pp. 45-52. doi:10.1007/s00170-004-2336-8
[18] K. F. He, X. J. Li and J. G. Wu, “Three-Dimensional
Temperature Field Numerical Simulation of Twin-Arc,”
Advanced Materials Research, Vol. 216, 2011, pp. 188-
193. doi:10.4028/www.scientific.net/AMR.216.188
[19] S. L. Huang, Y Xue and J. X. Zhang, “3D Finite Element
Simulation on Stress-Strain Distributions in All Position
Root Butt Welding of Pipe,” Transactions of the China
Welding Institution, Vol. 27, No. 4, 2006, pp. 73-75.
[20] H. N. Cai and M. Y. Tang, “FEM Analysis on TIG Weld-
ing Temperature Field,” Chinese Journal of Mechanical
Engineering, Vol. 32, No. 2, 1996, pp. 34-39.