<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">WJET</journal-id><journal-title-group><journal-title>World Journal of Engineering and Technology</journal-title></journal-title-group><issn pub-type="epub">2331-4222</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/wjet.2018.61010</article-id><article-id pub-id-type="publisher-id">WJET-82678</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Chemistry&amp;Materials Science</subject><subject> Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Myung</surname><given-names>Su Yi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Chung</surname><given-names>Min Hyun</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jeom</surname><given-names>Kee Paik</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Central Research Institute, Samsung Heavy Industries Co., Ltd., Geoje, Korea</addr-line></aff><aff id="aff1"><addr-line>Department of Naval Architecture and Ocean Engineering, Pusan National University, Busan, Korea</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>jeompaik@pusan.ac.kr(JKP)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>25</day><month>12</month><year>2017</year></pub-date><volume>06</volume><issue>01</issue><fpage>176</fpage><lpage>200</lpage><history><date date-type="received"><day>18,</day>	<month>December</month>	<year>2017</year></date><date date-type="rev-recd"><day>24,</day>	<month>February</month>	<year>2018</year>	</date><date date-type="accepted"><day>27,</day>	<month>February</month>	<year>2018</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, t
  hree-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.
 
</p></abstract><kwd-group><kwd>Steel Stiffened-Plate Structures</kwd><kwd> Weld-Induced Initial Distortion</kwd><kwd> Weld-Induced Residual Stress</kwd><kwd> Nonlinear Finite Element Method</kwd><kwd> Three-Dimensional Ther-mo-Elastic-Plastic Finite Element Analysis</kwd><kwd> Full Scale Measurements</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Ships and offshore structures are fabricated by welding and thus weld-induced initial imperfections in the form of initial distortions and residual stresses are unavoidable. It is required to predict the weld-induced initial imperfections in magnitude and pattern in association with structural design and construction.</p><p>A large number of studies are found in the literature. In recent years, these studies are associated with weld-induced initial distortions [<xref ref-type="bibr" rid="scirp.82678-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref5">5</xref>] and welding-induced residual stresses [<xref ref-type="bibr" rid="scirp.82678-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref10">10</xref>] . <xref ref-type="fig" rid="fig1">Figure 1</xref> shows typical modes of weld-induced initial distortions in plated structures [<xref ref-type="bibr" rid="scirp.82678-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref12">12</xref>] .</p><p>Simplified methods to predict weld-induced initial distortions are available in the literature [<xref ref-type="bibr" rid="scirp.82678-ref13">13</xref>] - [<xref ref-type="bibr" rid="scirp.82678-ref18">18</xref>] . However, weld-induced initial imperfections are affected by various parameters such as weld joint details, plate thickness, thickness transition (when the joint consists of plates with different thicknesses), stiffener spacing, numbers of attachments, corrugated constructions, mechanical restraint conditions, assembly sequences and overall construction planning [<xref ref-type="bibr" rid="scirp.82678-ref19">19</xref>] . In this regard, more refined method is still required to accurately predict them.</p><p>Nonlinear finite element method is a powerful tool to predict the weld-induced initial imperfections as it can deal with the distribution of heat transfer and general deformation due to welding. Ueda and Murakawa [<xref ref-type="bibr" rid="scirp.82678-ref20">20</xref>] undertook a pioneering work using two-dimensional finite element method to predict weld-induced residual stresses. Lindgrean and Karlsson [<xref ref-type="bibr" rid="scirp.82678-ref21">21</xref>] are the first to use three-dimensional finite element method for predicting weld-induced residual stresses. Chen and Guedes Soares [<xref ref-type="bibr" rid="scirp.82678-ref22">22</xref>] also used the three-dimensional finite element method to predict the weld-induced residual stresses with a validation by a comparison with test database on small scale models.</p><p>Studies on predicting the weld-induced initial distortions using finite element method are also found in the literature [<xref ref-type="bibr" rid="scirp.82678-ref23">23</xref>] - [<xref ref-type="bibr" rid="scirp.82678-ref28">28</xref>] . It was realized from the studies that the computing times are very expensive, and thus simplified methods using inherent strain have been suggested [<xref ref-type="bibr" rid="scirp.82678-ref27">27</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref29">29</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref30">30</xref>] . Michaleris and DeBiccari [<xref ref-type="bibr" rid="scirp.82678-ref31">31</xref>] developed a useful method to predict weld-induced distortions of large and complex structures, where two-dimensional thermal analysis was performed while three-dimensional structural analysis was undertaken. Regardless of the useful contributions in the literature, more efforts are still required for predicting the weld-induced initial imperfections of ships and offshore structures in association with welding and repair methods as well as robust design [<xref ref-type="bibr" rid="scirp.82678-ref32">32</xref>] .</p><p>And, studies on repair welding using finite element method and experiment are also found in the literature [<xref ref-type="bibr" rid="scirp.82678-ref33">33</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref34">34</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref35">35</xref>] .</p><p>In this paper, a three-dimensional thermo-elastic-plastic finite element method model is developed to predict the weld-induced distortions and residual stresses of thin-walled structures while achieving the goals in resulting accuracy and computational cost. The validity of the method is confirmed by a comparison with test database separately obtained from the full scale measurements by the authors.</p></sec><sec id="s2"><title>2. Finite Element Modeling for Predicting the Weld-Induced Residual Stresses</title><p>In this section, the finite element model is presented to predict the weld-induced residual stresses. The three-dimensional thermo-elastic-plastic finite element analysis method is applied.</p><p>Welding residual stress analysis is a complex problem requiring a lot of knowledge such as solid mechanics, heat transfer, materials science and so on. An important feature of this problem is the nonlinearity of the material. That is, material properties such as thermal conductivity, heat capacity, elastic modulus, yield stress, and Poisson’s ratio vary sensitively with temperature. Particularly in the case of welding, there is sudden change of the temperature from the atmospheric temperature to the melting temperature of the metal. In addition, during the thermal transition period, the material experiences state changes from the solidus to the liquidus, or vice versa. These state changes produce changes in material properties and specific volume. This is generally explained by the coefficient of thermal expansion. For this reason, it is very difficult to obtain an analytical solution to the temperature and stress field. However, recent advances in finite element analysis using computer have led to many improvements in the residual stress analysis of welding. In general, the finite element method provides very detailed information on transition temperature, displacement, strain and stress. This FE method provides an accurate solution to the welding problem. This is possible because the temperature-dependent changes in the material properties of the material can be taken into account and properties such as the state of the weld and the phase change can be sufficiently modeled. As described above, the features of the welding problem requires consideration of many physical phenomena. <xref ref-type="fig" rid="fig2">Figure 2</xref> shows the form of the connection between thermal, mechanical and metallurgical processes for the phenomena occurring during the welding process. This method is used to analyze the weld-induced residual stress.</p><p>The validation of the finite element model is made by a comparison with test database obtained by Paik and Yi [<xref ref-type="bibr" rid="scirp.82678-ref36">36</xref>] . The MSC/MARC computer program [<xref ref-type="bibr" rid="scirp.82678-ref37">37</xref>] is used for computations in the present study.</p><sec id="s2_1"><title>2.1. Modeling of Heat Source</title><p>As heat source moves associated with welding, the heat input can change accordingly. Rosenthal [<xref ref-type="bibr" rid="scirp.82678-ref38">38</xref>] suggested a model to simulate the moving heat source although the distribution of heat input was not considered due to the change of heat source. Goldak et al. [<xref ref-type="bibr" rid="scirp.82678-ref39">39</xref>] suggested a model to consider thermal distribution through moving heat source using Gaussian surface or double-ellipsoidal distributions. For welding with low depth of penetration, the Gaussian surface distribution model is well more suitable to apply. However, as the weld penetration becomes deep like in shipbuilding industry practice, heat transfer is of significance in the thickness direction. In this case, the double-ellipsoidal distribution model is more suitable.</p><p><xref ref-type="fig" rid="fig3">Figure 3</xref> presents the configuration of the heat source distribution with a double ellipsoidal shape. The heat flux for the first half-ellipsoidal internal area located at the front of the welding arc can then be defined as follows:</p><p>q f ( x , y , z ) = 6 3 f f Q a b c f π π exp ( − 3 x 2 a 2 ) exp ( − 3 y 2 b 2 ) exp ( − 3 z 2 c f 2 ) (1)</p><p>at z ≥ 0</p><p>where q f ( x , y , z ) is the heat flux for the first half-ellipsoidal internal area located at the front of the welding arc, f f is the heat input proportion in the front part, Q is the heat flux of arc which is taken as Q = η I U , a , b and c f are geometric parameters, and η is the arc efficiency, I is the current and U is</p><p>the voltage. x, y and z are the coordinate where the origin is located on the weld surface below the heat source as shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p>The arc efficiency η is defined as the ratio of the heat transferred from the arc to the base material and the total heat generated by the arc, which is a constant simply indicative of the complex heat loss between the substrates in the arc. Since the arc efficiency is influenced by various factors such as welding conditions, melting depth and protective gas, it is very challenging to predict it theoretically and it is usually determined experimentally. This is because the arc efficiencies given in the literature are somewhat different from each other. For flux cored arc welding considered in the present method, η = 0.75 is usually taken.</p><p>The heat flux at any point ( x , y , z ) in the second semi-ellipsoid covering the rear part of the arc can be defined as follows:</p><p>q r ( x , y , z ) = 6 3 f r Q a b c r π π exp ( − 3 x 2 a 2 ) exp ( − 3 y 2 b 2 ) exp ( − 3 z 2 c r 2 ) (2)</p><p>at z &lt; 0</p><p>where q r ( x , y , z ) is the heat flux at any point ( x , y , z ) in the second semi-ellipsoid covering the rear part of the arc, c r is the geometric parameter, and f r is the heat input proportion in the rear part. It is realized that f f + f r = 2 is approximately taken as far as the following conditions are satisfied:</p><p>f f = 2 ( 1 + c r / c f ) , f r = 2 ( 1 + c f / c r ) (3)</p></sec><sec id="s2_2"><title>2.2. Geometric Properties</title><p><xref ref-type="fig" rid="fig4">Figure 4</xref> shows a steel stiffened-plate structure which is extracted from living quarter structures of a very large crude oil carrier (VLCC) class floating, production, storage and offloading (FPSO) unit which is used for the production of oil in deep water. Typically, the longitudinal stiffener spacing (plate breadth) is in the range of 800 - 1000 mm, the transverse stiffener spacing (plate length) is in the range of 1800 - 3600 mm, and the plate thickness is in the range of 6 - 13 mm.</p><p>Paik and Yi [<xref ref-type="bibr" rid="scirp.82678-ref36">36</xref>] measured the weld-induced residual stresses of the structure shown in <xref ref-type="fig" rid="fig4">Figure 4</xref> where angle type longitudinal stiffener and T type transverse stiffener as shown in <xref ref-type="fig" rid="fig5">Figure 5</xref> are used. <xref ref-type="table" rid="table1">Table 1</xref> indicates the overall dimensions of the test structures. The weld-induced residual stresses for a total of four structures were measured with varying the plate thickness. It was noted that 1/2 + 1 + 1/2 bay structure model could minimize the interacting effect of adjacent plate panels. While details of the measurements are referred to in Paik and Yi [<xref ref-type="bibr" rid="scirp.82678-ref36">36</xref>] , the finite element model presented in this paper is validated by a comparison with the measurements. A quarter of the structure is taken as the extent of the present finite element analysis for the geometric symmetry as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>. <xref ref-type="table" rid="table1">Table 1</xref> indicates the geometric properties of the structures where the</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Dimensions of the test structures considered [<xref ref-type="bibr" rid="scirp.82678-ref33">33</xref>] </title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Type</th><th align="center" valign="middle" >a [mm]</th><th align="center" valign="middle" >b [mm]</th><th align="center" valign="middle" >t<sub>p</sub> [mm]</th><th align="center" valign="middle" >L<sub>w</sub> [mm]</th><th align="center" valign="middle" >Stiffener Type</th><th align="center" valign="middle" >h<sub>w</sub> [mm]</th><th align="center" valign="middle" >t<sub>w</sub> [mm]</th><th align="center" valign="middle" >b<sub>f</sub> [mm]</th><th align="center" valign="middle" >t<sub>f</sub> [mm]</th></tr></thead><tr><td align="center" valign="middle"  rowspan="4"  >Deck</td><td align="center" valign="middle"  rowspan="4"  >3200</td><td align="center" valign="middle"  rowspan="4"  >800</td><td align="center" valign="middle" >6</td><td align="center" valign="middle"  rowspan="4"  >6.5</td><td align="center" valign="middle"  rowspan="2"  >Longi.</td><td align="center" valign="middle"  rowspan="2"  >125</td><td align="center" valign="middle"  rowspan="2"  >7</td><td align="center" valign="middle"  rowspan="2"  >75</td><td align="center" valign="middle"  rowspan="2"  >7</td></tr><tr><td align="center" valign="middle" >8</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle"  rowspan="2"  >Trans.</td><td align="center" valign="middle"  rowspan="2"  >350</td><td align="center" valign="middle"  rowspan="2"  >12</td><td align="center" valign="middle"  rowspan="2"  >100</td><td align="center" valign="middle"  rowspan="2"  >17</td></tr><tr><td align="center" valign="middle" >12</td></tr></tbody></table></table-wrap><p>Note: L<sub>w</sub> indicates the weld bead length (leg length of welding).</p><p>plate thickness is varied as 6, 8, 10 and 12 mm.</p></sec><sec id="s2_3"><title>2.3. Material Properties</title><p><xref ref-type="table" rid="table2">Table 2</xref> indicates the mechanical properties of material at room temperature. During welding, the material properties greatly change with elevated temperatures as described in <xref ref-type="fig" rid="fig6">Figure 6</xref> [<xref ref-type="bibr" rid="scirp.82678-ref40">40</xref>] . It is noted that the minimum requirement of the yield strength specified by classification societies is 235 MPa for mild steel which is used for construction of the living quarter structures, but the actual yield strength of 309 MPa is used for the finite element analysis in the present study.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Mechanical properties of steel at room temperature used for the living quarter structures</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Kind of material</th><th align="center" valign="middle" >E [GPa]</th><th align="center" valign="middle" >σ<sub>Y</sub> [MPa]</th><th align="center" valign="middle" >σ<sub>u</sub> [MPa]</th><th align="center" valign="middle" >n [−]</th><th align="center" valign="middle" >Elongation [mm/mm]</th></tr></thead><tr><td align="center" valign="middle" >Mild A</td><td align="center" valign="middle" >203</td><td align="center" valign="middle" >309</td><td align="center" valign="middle" >458</td><td align="center" valign="middle" >0.3</td><td align="center" valign="middle" >0.23</td></tr></tbody></table></table-wrap><p>Note: E is the elastic modulus, σ<sub>Y</sub> is the yield strength, σ<sub>u</sub> is the ultimate tensile strength, and n is the Poisson’s ratio.</p></sec><sec id="s2_4"><title>2.4. Thermal Loading Load and Boundary Conditions</title><p>As mentioned earlier, a quarter model of the target structure is applied as the extent of the finite element analysis as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref> in association with the symmetric condition. <xref ref-type="fig" rid="fig7">Figure 7</xref> shows the sample model of the thermo-elastic-plastic finite element analysis together with loading and boundary conditions where the welding sequence is divided into 8 zones in the same order as the actual welding sequence. The structure is modelled using eight-node 3D brick thermal-mechanical finite elements or Solid 7 (full integration type) elements which treat plasticity, hyper-elasticity, stress stiffening, creep, large deflection, and large strain capabilities [<xref ref-type="bibr" rid="scirp.82678-ref37">37</xref>] .</p><p>Welding heat input is applied using the dual heat source model where thermal loading and boundary conditions with convection and radiation are allocated. With convective conditions, the related constants depending on elevated temperature are defined from the database obtained from experiments and computational fluid dynamics simulations [<xref ref-type="bibr" rid="scirp.82678-ref41">41</xref>] . <xref ref-type="fig" rid="fig8">Figure 8</xref> shows the variation of the temperature-dependent convection coefficient values [<xref ref-type="bibr" rid="scirp.82678-ref41">41</xref>] .</p><p><xref ref-type="table" rid="table3">Table 3</xref> indicates the heat input defined from welding parameters associated with weld bead length applied in the present study. <xref ref-type="table" rid="table4">Table 4</xref> gives the values of double ellipsoidal moving source model [<xref ref-type="bibr" rid="scirp.82678-ref22">22</xref>] .</p></sec><sec id="s2_5"><title>2.5. Computed Results and Discussion</title><p><xref ref-type="fig" rid="fig9">Figure 9</xref>(a) shows the process of the moving heat source and resulting temperature distributions when welding is performed at the second weld zone. <xref ref-type="fig" rid="fig9">Figure 9</xref>(b) shows the equivalent von Mises stress distribution for the structure with a plate thickness of 6 mm after welding and cooling are completed.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>0 shows the residual stress perpendicular to the weld line in relation to the section surface in the direction of each stiffener for a 6 mm thickness model.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Welding parameter for leg length</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Leg length [mm]</th><th align="center" valign="middle" >Current [A]</th><th align="center" valign="middle" >Voltage [V]</th><th align="center" valign="middle" >Speed [CPM]</th><th align="center" valign="middle" >Speed [mm/s]</th><th align="center" valign="middle" >Heat input [KJ/mm]</th></tr></thead><tr><td align="center" valign="middle" >6.5</td><td align="center" valign="middle" >320</td><td align="center" valign="middle" >32</td><td align="center" valign="middle" >36</td><td align="center" valign="middle" >6.0</td><td align="center" valign="middle" >1.71</td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Values of the double ellipsoidal heat source model [<xref ref-type="bibr" rid="scirp.82678-ref22">22</xref>] </title></caption><table><tbody><thead><tr><th align="center" valign="middle" >a [mm]</th><th align="center" valign="middle" >b [mm]</th><th align="center" valign="middle" >c<sub>f</sub> [mm]</th><th align="center" valign="middle" >c<sub>f</sub> [mm]</th><th align="center" valign="middle" >f<sub>f</sub></th><th align="center" valign="middle" >f<sub>r</sub></th></tr></thead><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >16</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >1.6</td></tr></tbody></table></table-wrap><p>The observed residual stress distribution is typically characterized by a high tensile stress near the weld bead and a distribution of compressive stress as it moves away from the weld.</p></sec><sec id="s2_6"><title>2.6. Comparison with Test Data</title><p>Paik and Yi [<xref ref-type="bibr" rid="scirp.82678-ref36">36</xref>] performed the experiment to obtained test database of the welding-induced residual stresses in a full scale prototype structure. <xref ref-type="fig" rid="fig1">Figure 1</xref>1 compares the experiments and the present computations. It is obvious that the welding-induced residual stresses are affected by plate thickness. With increase in the plate thickness, the compressive residual stresses decrease. It is also found that the magnitude of tensile residual stresses developed in the heat-affected zone reaches the material yield stress. The breadth of the heat-affected zone tends to decrease as the plate thickness increases for the same heat input. This may be due to the fact that the thermal load due to heat transfer is increased in the plate thickness direction and subsequently the magnitude of compressive residual stresses tends to decrease with increase in the plate thickness.</p><p>It is concluded that the present finite element analysis model gives very good agreements with test database.</p></sec></sec><sec id="s3"><title>3. Finite Element Modelling for Predicting the Weld-Induced Distortions</title><p>In this section, the finite element model is presented to predict the weld-induced distortions. The thermal buckling analysis method is applied. The validation of the finite element model is made by a comparison with test database obtained from Paik and Yi [<xref ref-type="bibr" rid="scirp.82678-ref36">36</xref>] . The MSC/MARC computer program [<xref ref-type="bibr" rid="scirp.82678-ref37">37</xref>] is used for computations in the present study.</p><sec id="s3_1"><title>3.1. Procedure of the Analysis</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref>2 presents the procedure for the nonlinear finite element analysis to predict the weld-induced distortions comprising 3 steps [<xref ref-type="bibr" rid="scirp.82678-ref42">42</xref>] . In step I, the welding</p><p>analysis is performed using strain-as-direct-boundary (SDB) method [<xref ref-type="bibr" rid="scirp.82678-ref43">43</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref44">44</xref>] . In step II, shrinkage forces are calculated in association with the weld-induced residual stresses of welding. In step III, thermal buckling analysis is performed to determine the weld-induced distortions.</p></sec><sec id="s3_2"><title>3.2. Material Properties</title><p>The mechanical properties of material in <xref ref-type="table" rid="table2">Table 2</xref> are used. The yield strength of both plating and stiffeners is 309 MPa. The elastic modulus and Poisson’s ratio are 203 GPa and 0.3, respectively. The elastic-perfectly material model without strain-hardening effect is applied.</p></sec><sec id="s3_3"><title>3.3. Thermal Loading and Boundary Conditions</title><p>For welding analysis to obtain initial irregularity in thermal buckling analysis, contact analysis is carried out using self-weight and spring elements connected on four vertexes. The spring node constrains all six degrees of freedom to prevent rigid body motion. In order to simulate the actual welding process situation, the planar state is analyzed to the plane and the contact analysis with the structural element is carried out. The contact friction uses a stick-slip model with excellent convergence. A two-bay (1/2 + 1 + 1/2) model both in the longitudinal and transverse directions is taken as the extent for the analysis as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>3. The symmetry condition at all edges is then applied in association with structural continuity.</p><p>The thermal loads obtained from steps I and II indicated in <xref ref-type="table" rid="table5">Table 5</xref> are applied to individual finite elements. <xref ref-type="fig" rid="fig1">Figure 1</xref>3(a) shows the boundary conditions. The imaginary welding temperature (T<sub>top</sub> and T<sub>bot</sub>) loads indicated in <xref ref-type="table" rid="table5">Table 5</xref></p><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Shrinkage forces depending on the plate thickness [<xref ref-type="bibr" rid="scirp.82678-ref45">45</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref46">46</xref>] </title></caption><table><tbody><thead><tr><th align="center" valign="middle" >t<sub>p</sub><sub> </sub> [mm]</th><th align="center" valign="middle" >L<sub>w</sub> [mm]</th><th align="center" valign="middle" >Shrinkage force [N]</th><th align="center" valign="middle" >Imaginary temperature at top, T<sub>top</sub> [˚C]</th><th align="center" valign="middle" >Imaginary temperature at bottom, T<sub>bot</sub> [˚C]</th></tr></thead><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle"  rowspan="4"  >6.5</td><td align="center" valign="middle" >45,466.49</td><td align="center" valign="middle"  rowspan="4"  >0.174</td><td align="center" valign="middle"  rowspan="4"  >−0.174</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >36,963.12</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >31,861.09</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >28,459.75</td></tr></tbody></table></table-wrap><p>[<xref ref-type="bibr" rid="scirp.82678-ref45">45</xref>] [<xref ref-type="bibr" rid="scirp.82678-ref46">46</xref>] are applied at the top and bottom of the plates along the welding lines. After the thermal buckling analysis is completed, thermal shrinkage forces are applied to the target structure with modified boundary conditions as described in <xref ref-type="fig" rid="fig1">Figure 1</xref>3(b).</p></sec><sec id="s3_4"><title>3.4. Computed Results and Discussion</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref>4 shows the computed results of the weld-induced distortions. <xref ref-type="fig" rid="fig1">Figure 1</xref>5</p><p>presents the computed results of the weld-induced distortions after welding at the plate center in the longitudinal direction. According to the analysis procedure shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>2, at first welding analysis for initial imperfection was performed. The initial imperfection is assumed to be the deformation that occurs due to welding, and then thermal buckling analysis is performed by applying the shrinkage force as the boundary condition to this model. It is observed that the deflection at the welding line, i.e., on the transverse stiffeners is zero and has the maximum values at the middle of plating between stiffeners. The magnitude of the initial deflection in a thick plate is smaller than that in a thin plate.</p></sec><sec id="s3_5"><title>3.5. Comparison with Test Data</title><p>Paik and Yi [<xref ref-type="bibr" rid="scirp.82678-ref36">36</xref>] performed the experiment to obtained test database of the welding-induced distortions in a full scale prototype structure. <xref ref-type="fig" rid="fig1">Figure 1</xref>6 compares the weld-induced distortions obtained from experiment as well as the computations. It is concluded that the computations are in very good agreement with the experimental results.</p></sec></sec><sec id="s4"><title>4. Conclusions</title><p>The objectives of the present study were to develop finite element modelling techniques to predict weld-induced residual stresses and weld-induced distortions in a stiffened plate structures. The developed techniques were validated using test database obtained from direct measurements from a full scale prototype stiffened plate structures. Based on the results obtained from this study, the following conclusions are drawn.</p><p>1) In terms of predicting the weld-induced residual stresses, the modelling of heat source is important. In the present study, a double ellipsoidal shape model was applied to describe the heat source distribution during welding. It was confirmed that the numerical computations applying the proposed model are in good agreement with test data.</p><p>2) The procedure for predicting the weld-induced distortions comprises three steps, namely the weld temperature analysis, the shrinkage force computations and the thermal buckling analysis. The proposed procedure was implemented into a commercial computer code. It was confirmed that the numerical computations are in good agreement with test data.</p></sec><sec id="s5"><title>Acknowledgements</title><p>The present study was undertaken in the Korea Ship and Offshore Research Institute at Pusan National University which has been a Lloyd’s Register Foundation Research Centre of Excellence. The support of Samsung Heavy Industries is greatly acknowledged.</p></sec><sec id="s6"><title>Cite this paper</title><p>Yi, M.S., Hyun, C.M. and Paik, J.K. (2018) Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures. World Journal of Engineering and Technology, 6, 176-200. https://doi.org/10.4236/wjet.2018.61010</p></sec></body><back><ref-list><title>References</title><ref id="scirp.82678-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Bruno, G., Yordan, G. and Soares, C.G. (2011) Effect of Weld Shape Imperfections on the Structural Hot-Spot Stress Distribution. Ships and Offshore Structures, 6, 145-159. https://doi.org/10.1080/17445302.2010.497052</mixed-citation></ref><ref id="scirp.82678-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Teresa, M. and Craig, F. (2013) Effect of Weld-Induced Imperfections on the Ultimate Strength of an Aluminium Patrol Boat Determined by the ISFEM Rapid Assessment Method. Ships and Offshore Structures, 9, 218-235.</mixed-citation></ref><ref id="scirp.82678-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Luís, R.M., Soares, C.G. and Nikolov, P.I. (2009) Collapse Strength of Longitudinal Plate Assemblies with Dimple Imperfections. Ships and Offshore Structures, 3, 359-370. https://doi.org/10.1080/17445300802370024</mixed-citation></ref><ref id="scirp.82678-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Vhanmane, S. and Bhattacharya, B. (2008) Estimation of Ultimate Hull Girder Strength with Initial Imperfections. Ships and Offshore Structures, 3, 149-158.  
https://doi.org/10.1080/17445300802204389</mixed-citation></ref><ref id="scirp.82678-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Khedmati, M.R., Pedram, M. and Rigo, P. (2012) The Effects of Geometrical Imperfections on the Ultimate Strength of Aluminium Stiffened Plates Subject to Combined Uniaxial Compression and Lateral Pressure. Ships and Offshore Structures, 9, 88-109. https://doi.org/10.1080/17445302.2012.726761</mixed-citation></ref><ref id="scirp.82678-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Farajkhah, V., Liu, Y. and Gannon, L. (2016) Finite Element Study of 3D Simulated Welding Effect in Aluminium Plates. Ships and Offshore Structures, 12, 196-208.  
https://doi.org/10.1080/17445302.2015.1123865</mixed-citation></ref><ref id="scirp.82678-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Gannon, L.G., Liu, Y., Pegg, N.G. and Smith, M.J. (2012) Effect of Three-Dimensional Welding-Induced Residual Stress and Distortion Fields on Strength and Behaviour of Flat-Bar Stiffened Panels. Ships and Offshore Structures, 8, 565-578.  
https://doi.org/10.1080/17445302.2012.707386</mixed-citation></ref><ref id="scirp.82678-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Gannon, L., Liu, Y., Pegg, N. and Smith, M.J. (2015) Nonlinear Collapse Analysis of Stiffened Plates Considering Welding-Induced Residual Stress and Distortion. Ships and Offshore Structures, 11, 228-244.  
https://doi.org/10.1080/17445302.2014.985428</mixed-citation></ref><ref id="scirp.82678-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Khan, I. and Zhang, S. (2011) Effects of Welding-Induced Residual Stress on Ultimate Strength of Plates and Stiffened Panels. Ships and Offshore Structures, 6, 297-309. https://doi.org/10.1080/17445301003776209</mixed-citation></ref><ref id="scirp.82678-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Gannon, L.G., Pegg, N.G., Smith, M.J. and Liu, Y. (2013) Effect of Residual Stress Shakedown on Stiffened Plate Strength and Behaviour. Ships and Offshore Structures, 8, 638-652. https://doi.org/10.1080/17445302.2012.664429</mixed-citation></ref><ref id="scirp.82678-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Ueda, Y. (1999) Computational Welding Mechanics (Selected from Papers Published in English). Japan Osaka: Joining and Welding Research Institute, Osaka University, Osaka.</mixed-citation></ref><ref id="scirp.82678-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Paik, J.K. and Thayamballi, A.K. (2003) Ultimate Limit State Design of Steel-Plated Structures. Wiley, Chichester.</mixed-citation></ref><ref id="scirp.82678-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Masubuchi, K. (1980) Analysis of Welded Structures: Residual Stresses, Distortion and Their Consequences. Pergamon Press, Oxford.</mixed-citation></ref><ref id="scirp.82678-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Radaj, D. (1992) Heat Effects of Welding: Temperature Field, Residual Stress and Distortion. Springer Verlag Publishing, Berlin Heidelberg.  
https://doi.org/10.1007/978-3-642-48640-1</mixed-citation></ref><ref id="scirp.82678-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Lindgren, L.E. (2007) Computational Welding Mechanics: Thermo Mechanical and Microstructural Simulations. Woodhead Publishing Limited, Cambridge.</mixed-citation></ref><ref id="scirp.82678-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Goldak, J.A. and Akhlaghi, M. (2005) Computational Welding Mechanics. Springer Science &amp; Business Media, New York.</mixed-citation></ref><ref id="scirp.82678-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Feng, Z. (2005) Processes and Mechanisms of Welding Residual Stress and Distortion. Woodhead Publishing Limited, Cambridge.</mixed-citation></ref><ref id="scirp.82678-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Michaleris, P. (2011) Minimization of Welding Distortion and Buckling: Modelling and Implementation. Woodhead Publishing Limited, Cambridge.  
https://doi.org/10.1533/9780857092908</mixed-citation></ref><ref id="scirp.82678-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Tsai, C.L., Park, S.C. and Chueng, W.T. (1999) Welding Distortion of a Thin-Plate Panel Structure. Welding Research Supplement, 78, 157-165.</mixed-citation></ref><ref id="scirp.82678-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Ueda, Y. and Yamakawa, T. (1971) Analysis of Thermal Elastic-Plastic Stress and Strain during Welding by Finite Element Method. Transactions of Japan Welding Research Institute, 2, 90-100.</mixed-citation></ref><ref id="scirp.82678-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Lindgren, L.E. and Karlsson, L. (1988) Deformations and Stresses in Welding of Shell Structures. International Journal for Numerical Methods in Engineering, 25, 635-655. https://doi.org/10.1002/nme.1620250223</mixed-citation></ref><ref id="scirp.82678-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Chen, B.Q., Hashemzadeh, M., Garbatov, Y. and Guedes Soares, C. (2015) Numerical and Parametric Modeling and Analysis of Weld-Induced Residual Stresses. International Journal of Mechanical and Material in Design, 11, 439-453.  
https://doi.org/10.1007/s10999-014-9269-7</mixed-citation></ref><ref id="scirp.82678-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Guangming, F., Marcelo, I.L., Menglan, D. and Segen, F.E. (2016) Influence of the Welding Sequence on Residual Stress and Distortion of Fillet Welded Structures. Marine Structures, 46, 30-55. https://doi.org/10.1016/j.marstruc.2015.12.001</mixed-citation></ref><ref id="scirp.82678-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Wang, J., Yuan, H., Ma, N. and Murakawa, H. (2016) Recent Research on Welding Distortion Prediction in Thin Plate Fabrication by Means of Elastic FE Computation. Marine Structures, 47, 42-59. https://doi.org/10.1016/j.marstruc.2016.02.004</mixed-citation></ref><ref id="scirp.82678-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Gray, T., Camilleri, D. and McPherson, N. (2014) Control of Welding Distortion in Thin-Plate Fabrication: Design Support Exploiting Computational Simulation. Woodhead Publishing Limited, Cambridge.</mixed-citation></ref><ref id="scirp.82678-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Deng, D., Murakawa, H. and Liang, W. (2007) Numerical Simulation of Welding Distortion in Large Structures. Computer Methods in Applied Mechanics and Engineering, 196, 4613-4627. https://doi.org/10.1016/j.cma.2007.05.023</mixed-citation></ref><ref id="scirp.82678-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Luo, Y., Deng, D., Xie, L. and Murakawa, H. (2004) Prediction of Deformation for Large Welded Structures Based on Inherent Strain. Transactions of JWRI, 33, 65-70.</mixed-citation></ref><ref id="scirp.82678-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Wang, R., Zhang, J.X., Serizawa, H. and Murakawa, H. (2009) Study of Welding Inherent Deformations in Thin Plates Based on Finite Analysis using Interactive Substructure Method. Materials &amp; Design, 30, 3474-3481.  
https://doi.org/10.1016/j.matdes.2009.03.015</mixed-citation></ref><ref id="scirp.82678-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Daniewicz, S.R., McAninch, M.D., Mc-Farland, B. and Knoll, D. (1993) Application of Distortion Control Technology during Fabrication of Large Offshore Structures. Proceedings of AWS/ORNL International Conference on Modeling and Control of Joining Processes, Orlando, 8-10 December 1993, 239-246.</mixed-citation></ref><ref id="scirp.82678-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Murakawa, H., Luo, Y. and Ueda, Y. (1996) Prediction of Welding Deformation and Residual Stress by Elastic FEM Based on Inherent Strain (First Report)—Mechanism of Inherent Strain Production. Journal of the Society of Naval Architects of Japan, 180, 739-751. (In Japanese) https://doi.org/10.2534/jjasnaoe1968.1996.180_739</mixed-citation></ref><ref id="scirp.82678-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Michaleris, P. and DeBiccari, A. (1997) Prediction of Welding Distortion. Welding Journal, 76, 172-179.</mixed-citation></ref><ref id="scirp.82678-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Yi, M.S., Kim, S.H., Kim, S.J., Seo, J.K. and Paik, J.K. (2017) Development of an Empirical Formula for Predicting Weld-Induced Initial Deflection in Steel Stiffened-Plate Structures: Numerical and Experimental Studies.</mixed-citation></ref><ref id="scirp.82678-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Jiang, W., Luo, Y., Zhang, G. and Tu, S.-T. (2013) Experimental to Study the Effect of Multiple Weld-Repairs on Microstructure, Hardness and Residual Stress for a Stainless Steel Clad Plate. Materials &amp; Design, 51, 1052-1059.  
https://doi.org/10.1016/j.matdes.2013.05.027</mixed-citation></ref><ref id="scirp.82678-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Jiang, W., Xu, X.P., Gong, J.M. and Tu, S.-T. (2012) Influence of Repair Length on Residual Stress in the Repair Weld of a Clad Plate. Nuclear Engineering and Design, 246, 211-219. https://doi.org/10.1016/j.nucengdes.2012.01.021</mixed-citation></ref><ref id="scirp.82678-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Jiang, W., Liu, Z., Gong, J.M. and Tu, S.-T. (2010) Numerical Simulation to Study the Effect of Repair Width on Residual Stresses of a Stainless Steel Clad Plate. International Journal of Pressure Vessels and Piping, 87, 457-463.  
https://doi.org/10.1016/j.ijpvp.2010.06.003</mixed-citation></ref><ref id="scirp.82678-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Paik, J.K. and Yi, M.S. (2016) Experimental and Numerical Investigations of Welding-Induced Distortions and Stresses in Steel Stiffened Plate Structures. The Korea Ship and Offshore Research Institute, Pusan National University, Busan.</mixed-citation></ref><ref id="scirp.82678-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">MSC (2013) Online User’s Guide. MSC Software Corporation, Newport Beach.</mixed-citation></ref><ref id="scirp.82678-ref38"><label>38</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Rosenthal</surname><given-names> D. </given-names></name>,<etal>et al</etal>. (<year>1941</year>)<article-title>Mathmatical Theory of Heat Distribution during Welding and Cutting</article-title><source> Welding Journal</source><volume> 20</volume>,<fpage> 220</fpage>-<lpage>234</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.82678-ref39"><label>39</label><mixed-citation publication-type="other" xlink:type="simple">Goldak, J., Chakravarti, A. and Bibby, M. (1984) A New Finite-Element Model for Welding Heat-Sources. Metallurgical and Materials Transactions B, 15, 299-305.  
https://doi.org/10.1007/BF02667333</mixed-citation></ref><ref id="scirp.82678-ref40"><label>40</label><mixed-citation publication-type="other" xlink:type="simple">Yi, M.S. (2006) The Effect of Geometry Non-Linearity in Welding Deformation. Ms.D. Thesis, Pusan National University, Busan. (In Korean)</mixed-citation></ref><ref id="scirp.82678-ref41"><label>41</label><mixed-citation publication-type="other" xlink:type="simple">Yi, M.S., Hyun, C.M., Cho, S.H. and Jang, T.W. (2009) A Study of Correction Effect in T-Bar According to Cooling Method. Proceeding of the 2009 Spring Annual Meeting of Korea Welding Society, Changwon, 21-22 April 2009, 60. (In Korean)</mixed-citation></ref><ref id="scirp.82678-ref42"><label>42</label><mixed-citation publication-type="other" xlink:type="simple">Yi, M.S., Hyun, C.M. and Yoo, J.S. (2016) Development of New Thermal Buckling Analysis Technique Caused by Welding. International Welding/Joining Conference, Gyeongju, 11-14 April 2017, 21.</mixed-citation></ref><ref id="scirp.82678-ref43"><label>43</label><mixed-citation publication-type="other" xlink:type="simple">Ha, Y.S., Cho, S.H. and Jang, T.W. (2008) Development of Welding Distortion Analysis Method using Residual Strain as Boundary Condition. Material Science Forum, 580-582, 649-654.  
https://doi.org/10.4028/www.scientific.net/MSF.580-582.649</mixed-citation></ref><ref id="scirp.82678-ref44"><label>44</label><mixed-citation publication-type="other" xlink:type="simple">Ha, Y. (2011) A Study on Weldment Boundary Condition for Elasto-Plastic Thermal Distortion Analysis of Large Welded Structures. Journal of Welding and Joining, 29, 410-415. https://doi.org/10.5781/KWJS.2011.29.4.410</mixed-citation></ref><ref id="scirp.82678-ref45"><label>45</label><mixed-citation publication-type="other" xlink:type="simple">Shin, S.B., Kim, K.K. and You, J.G. (2006) A Study on Control of Buckling Distortion at the Thin Panel during Welding. International Offshore and Polar Engineering Conference, San Francisco, CA, 28 May-2 June 2006, 236-242.</mixed-citation></ref><ref id="scirp.82678-ref46"><label>46</label><mixed-citation publication-type="other" xlink:type="simple">Lee, D.J. and Shin, S.B. (2003) A Study on the Prediction of Shrinkage during the Manufacturing of a Deckhouse of RIG. International Offshore and Polar Engineering Conference, Honolulu, 25-30 May 2003, 166-171.</mixed-citation></ref></ref-list></back></article>