• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Transient out-of-plane distortion of multi-pass fillet welded tube to pipe T-joints

    2017-06-01 11:35:04VetriselvnDevkumrnPSthiyRvichndrn
    Defence Technology 2017年2期

    R.Vetriselvn*,K.DevkumrnP.Sthiy,G.Rvichndrn

    aWelding Research Institute,Tiruchirappalli 620014,India

    bNational Institute of Technology,Tiruchirappalli 620015,India

    Transient out-of-plane distortion of multi-pass fillet welded tube to pipe T-joints

    R.Vetriselvana,*,K.Devakumarana,P.Sathiyab,G.Ravichandrana

    aWelding Research Institute,Tiruchirappalli 620014,India

    bNational Institute of Technology,Tiruchirappalli 620015,India

    A R T I C L E I N F O

    Article history:

    Received 5 February 2016

    Received in revised form

    30 May 2016

    Accepted 12 June 2016

    Available online 19 July 2016

    Multi-pass welds

    FEM

    Temperature field

    Out-of-plane distortion

    tube to pipe T-joints

    Circumferential welding

    The out-of-plane distortion induced in a multi-pass circumferential fillet welding of tube to pipe under different weld sequences and directions was studied using Finite Element Method(FEM)based Sysweld software and veri fied experimentally.The FEM analyses consisted of thermal and mechanical analyses. Thermal analysis was validated with experimental transient temperature measurements.In the mechanical analysis,three different weld sequences and directions were considered to understand the mechanism of out-of-plane distortion in the tube to pipe T-joints.It was learnt that the welding direction plays a major role in minimizing the out-of-plane distortion.Further,during circumferential fillet welding of the tube to pipe component,the out-of-plane distortion generated in the x direction was primarily influenced by heat input due to the start and stop points,whereas the distortion in the z direction was influenced by time lag and welding direction.The FEM predicted distortion was compared with experimental measurements and the mechanism of out-of-plane distortion was con firmed.

    ?2016 The Authors.Published by Elsevier Ltd.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    1.Introduction

    Welded T-joints are one of the commonly used joint types.Often pipes and tubes are welded in perpendicular orientation constituting a T-joint.This is possible with circumferential welding.These circumferential fillet welded socket joints are widely used in various industries such as line pipes,automobiles,bridges,defense components,pressure vessels and structural applications.Owing to the relatively large wall thickness in certain tube to pipe systems, the fillet weld is often constructed with multiple weld passes. During the fabrication of these components welding induced distortion occurs.

    T-joints experience out-of-plane distortion in both longitudinal and circumferential types of welding.Many authors have studied the distortion experienced due to longitudinal welding on a plate [1-6].The temperature field of an arc and its mechanical effect on fillet welding of plates using FEM are also reported[6]and it was found that the arc start and stop points enhance the heat input and accordingly increase distortion in the case of longitudinal welding.

    The out-of-plane distortion due to circumferential welding of thin-walled butt-welded pipes using FEM has been reported with the effect of weld thermal cycle on out-of-plane distortion generated under the different segment of the circumferential welding [7,8].The buckling distortion is also another type of out-of-plane distortion experienced in case of thin cylindrical shells during circumferential welding[9,10].Generally,out-of-plane distortion resulted in both longitudinal and circumferential fillet welding and this leads to misalignment,stress concentration,buckling and improper fit-up etc.However very few articles reported the out-ofplane distortion experienced during the circumferential welding of t-joints.The start and stop points have signi ficant effect on out-ofplane distortion in the case of half-circumferential fillet welding in a T-joint con figuration,and it has been reported that the stop points have more in fluence on out-of-plane distortion than the start points[11,12].In the case of full circumferential welding,both the start and stop points merge with each other and the effect is not studied so far.It is inferred from the literature survey that transient analysis of the full circumferential multi-pass welding on a tube to pipe component in the t-joint orientation is not available either experimentally or by other modes such as FEM simulation.These data,when made available,are of vital importance for various industrial and defense applications.

    Thus,the present work describes the out-of-plane distortiongenerated in multi-pass full circumferential fillet welding of tube to pipe under different weld sequences and directions.The studies were carried out using FEM based Sysweld[13]software.Firstly,the thermal analysis of the FEM model for the tube to pipe welding was carried out followed by mechanical analysis.In the mechanical analysis,three different weld sequences and directions were considered to understand the mechanism of out-of-plane distortion in the tube to pipe welding.Further,the predicted FEM model was veri fied experimentally.Based on the FEM and experimental results,the suitable welding procedure was explored for the tube to pipe fabrication.

    2.Experimentation

    The welding arrangement considered for the out-of-plane distortion analysis included the thick walled pipe consisting of lots of small tubes welded with the pipe.The schematic diagram and photographic view of the pipe and tube are shown in Fig.1(a) and(b)respectively.The thickness and outer diameter of the pipe is 20 mm and 130 mm.The thickness and outer diameter of the tube are 5 mm and 44.5 mm.The material used for pipe and tube are low carbon steel of grade EN S355J2G3.The chemical composition of the material is given in Table 1.The multiple numbers of tubes are fillet welded using Shielded Metal Arc Welding(SMAW)process and multi-pass welding procedure used for the fillet welding of the tube to pipe component.Prior to welding,holes are drilled in the pipe and the tubes are inserted into it.The joint type is set on type and is shown schematically in Fig.2(a).

    2.1.Welding

    The tubes are inserted into the holes in the pipe and it is tack welded.After that,the tubes are aligned with the pipe’s central axis using a standard jig.Three weld sequences in reference to the start and stop points of welding direction were considered for the tube to pipe fabrication.The welding parameters used in the experimentation are given in Table 2.The welding parameters such as welding current and arc voltage were recorded using digital meter fitted in the welding power source during welding.

    2.2.Measurement of temperature during welding

    The temperature was measured during the tube to pipe welding using a thermocouple.The temperature measurement requires thermocouples and a data acquisition card.Two thermocouples are positioned in the weld region as shown in schematic representation in Fig.2(b).The K-type thermocouple was used to measure the transient temperature in the weld vicinity.One thermocouple was positioned near the tube(TC1,5 mm from weld)and another was positioned nearer to the pipe(TC2,15 mm from weld metal).The data logging frequency for the thermocouple was set to 600 Hz.A computer-controlled data logger was used to capture the transient temperature distribution during welding of the tube with pipe.The process of data logging was shielded against electrical interference. The temperature sensors and data logger instrument were calibrated before the experimental trial.Data logging was started before welding and the measurements were taken during first, second and third weld passes and until the distortion reached steady state values.Photographic view of the thermocouple on the tube to pipe set up with the data logger is shown in Fig.3(a)and(b).

    Fig.1.The schematic diagram and the photographic view of the tube to pipe arrangement.

    Fig.2.The schematic diagram of set on type joint used in the tube to pipe welding.

    Table 1The chemical composition of the carbon steel used for the study.

    Table 2The welding parameters used for tube to pipe fabrication.

    2.3.Measurement of out-of-plane distortion

    The transient measurement of distortion of the tube during welding could not be measured as the length of the tube was only 110 mm and positioning displacement sensors at the end of the tube caused accessibility issues to the welder during all around circumferential welding.Hence,the out-of-plane distortionmeasurement was carried out using Vernier caliper after completion of every weld pass.Prior to welding of tube to pipe,the standard jig was mounted on the pipe to monitor the movement of the tube during circumferential weld deposition as shown in Fig.4.The jig allows free movement of the tube in all the directions and enables the measurement of distortion from the initial position of the tube.The movement of the tube from its original position is referred to as out plane distortion in the tube to pipe fabrication. Under different weld sequences and weld directions,the out-ofplane distortion generated after welding(reached to room temperature)was measured using a Vernier-caliper having the least count of 0.02 mm.

    Fig.3.(a)The photographic view of thermocouple on the tube to pipe component.(b)Data logger used during the study.

    3.Out-of-plane distortion analysis

    3.1.Finite element model

    The out-of-plane distortion generated during multi-pass fillet welding in the tube to pipe fabrication under different weld sequences was simulated.The material properties considered for the investigation are temperature dependent.The sizes of the pipe and the tube considered for FEM simulation are shown in Fig.1(a).The tube to pipe model was meshed with solid elements.The meshed model consists of 25,000 nodes and 20,000 solid elements including tetrahedral and brick elements as shown in Fig.5.The three weld passes were considered for the FEM analysis as shown in Fig.6.FEM analysis of multipass weld deposition was carried out using element activation and deactivation technique[14].The technique implies that the element is active until the temperature is dissipated to the surrounding through conduction,convection, and radiation when the heat is applied to the element.The death of the element will occur when the element reaches ambient temperature.The size of the finite element considered in the arc zone is about 2 mm and it increases away from the weld zone.

    The isotropic elastic plastic approach was used for the analyses. The thermal analysis is based on the heat conduction formulation using temperature dependent thermophysical properties (Fig.7(a)).The temperature dependent material properties considered for the mechanical investigation are shown in Fig.7(b) [12].The Goldak’s non-axisymmetric 3D double ellipsoid heat source was considered for the transient thermal analysis of multipass weld deposition of the tube to pipe welding[15].The arc heat source moves through the de fined one-dimensional circumferential trajectory at a given welding speed.The present FEM model also includes the heat losses due to convection and radiation during welding.The ambient temperature of 25°C was considered.The heat input is estimated as follows

    where V is the arc voltage(volt),I is the welding current(ampere),S is the welding speed(mm·s-1)and?is process efficiency 65(%). The welding parameter used for the present investigation is given in Table 2.Based on the above considerations,the out-of-plane distortion generated on the tube to pipe fabrication during multipass weld deposition was simulated under three weld sequences in reference to the start and stop points as well as welding direction as shown schematically in Fig.8.In addition,the clamping conditions contributed to the plastic stresses induced in the weldments. Hence,for the simulation,minimum clamping conditions were assumed to avoid the rigid body motion of the component.Three nodal points were rigidly clamped in the x,y,z directions.The clamping conditions considered for the simulation are shown in Fig.9.

    Fig.4.The measurement of distortion using standard jig during tube to pipe welding.

    The start and stop points of all three-weld passes are the same and follow the same welding directions(0°-90°-180°-270°-0°, Fig.8(a))in weld sequence-1.In the case of weld Sequence-2,the start point of weld pass 1 is at 0°and weld passes 2 and 3 are at 180°.It means that,in the case of weld pass 1,the direction of arc movement is 0°-90°-180°-270°-0°,whereas weld passes 2 and 3 are 180°-90°-0°-270°-180°(Fig.8(b)).In Sequence-3(Fig.8(c)),the start and stop points of weld passes 1 and 2 are at 0°and 180°, respectively,and the corresponding arc directions of respective weld passes are 0°-90°-180°-270°-0°and 180°-90°-0°-270°-180°.In the case of weld pass 3,the start and stop points of the arc is the same as the weld pass 2(180°),but the direction of welding changes opposite to weld pass 2 (180°-270°-0°-90°-180°).Firstly,the computation of the temperature history during welding and subsequent cooling was completed and then,the temperature field was applied to the mechanical model as a body force to perform the out-of-plane distortion analysis.

    Fig.5.The finite element model of a pipe with tube using solid elements.

    Fig.6.The cross sectional view of three weld passes between tube and pipe.

    4.Results and discussion

    4.1.Multi-pass weld thermal analysis of tube to pipe welding

    In order to study the effect of weld thermal cycle on distortion generated during multipass weld deposition of the tube to pipe welding,two different locations,TC1 and TC2,were considered as shown in Fig.2.At a given heat input and welding direction,the simulated and experimental weld thermal cycles under two different locations are shown in Figs.10 and 11 respectively.In Fig.10,it is observed that there are three peaks in the weld thermal cycle,con firming three weld passes involved in the tube to pipe component.Also,the maximum temperatures simulated at each weld pass in locations TC1 and TC2 differ from each other and are given in Table 3.This is primarily due to the fact that the location of TC1 is 5 mm away from the weld region and the location of TC2 is 15 mm away from the weld region.

    The variation in peak temperatures in TC1,between three weld passes,is attributed to the fact that weld pass 1 was carried out between the tube and the pipe with equal heat input distribution, whereas weld pass 2 was oriented between the pipe and the first weld pass.Hence,the tube receives less heat input comparatively and so peak 2 shows a lower temperature than peak 1.However, Weld pass 3 is oriented between the first weld pass and the tube and so peak 3 shows a higher temperature compared to peak 2.The maximum temperature recorded at TC1 is 1120°C as the thermocouple is located in the tube nearer to weld.

    In Fig.11,it is observed that each weld pass shows two peak temperatures.It is due to the fact that the TC2 is heated twice as start and stop points lie at the same location.These peaks are the characteristics of full circumferential welding.However,it is further observed that the peak 2 temperature is higher than the peak 1 temperature in all the three weld passes.During circumferential welding,the start point is once again reheated by the welding arc at the end of weld completion and therefore shows a higher temperature.

    The cooling rate of each weld pass can be estimated in the region from its peak to 500°C because this region is considered as a critical cooling condition in case of EN S355J2G3[12].It is observed that due to lower peak temperature,location TC2 shows lesser cooling rate than location TC1 irrespective of change in weld passes. Based on the results discussed above,the temperature contour plot obtained in the two different locations of TC1 and TC2 duringmultipass welding of tube to pipe at a given heat input and welding direction is shown in Fig.12.

    Fig.7.The variation of properties with temperature.(a)Thermal properties and(b)mechanical properties.

    The peak temperatures obtained in the thermal analyses for experimental trials are also given in Table 3.It is found that both the FEM simulated and experimental curves follow a similar trend and the maximum temperatures reached during the welding are in agreement.It is due to the assumption of reasonable(material) thermal properties such as convection,heat transfer coefficients, co-efficient of thermal expansion and contraction,and thermal conductivity during the formulation of the thermal model.Similar material properties assumed for carbon steel in the study yielded good results elsewhere[14].

    The deviation of maximum temperature reached by the FEM model and experimental result is 10%.It can also be noted that there is an offset between FEM and experimental curves in time(xaxis).It is attributed to two reasons majorly.

    Fig.8.The schematic diagram showing the three welding sequences considered for out-of-plane distortion analysis on multipass weld deposition of the tube to pipe component.

    Fig.9.The clamping condition considered to avoid the rigid body motions of the tube to pipe component.

    Fig.10.Thermal cycle at a node in TC1.

    1)The position of the thermocouple in the FEM model is at the exact location,whereas identifying the same location in the experimental trial is arbitrary and has a slight offset.

    2)The welding process used for the investigation is SMAW.The manual welding time between different weld passes is not uniform in the case of the experimental trial.This causes a slight time offset in the x-axis,whereas the FEM simulations follow a uniform time gap between the weld passes.Due to this time shift in the x axis,the FEM and experimental curves are slightly shifted.

    Table 3The simulated and experimental peak temperatures at locations TC1 and TC2.

    Fig.11.The thermal cycle at a node in TC2.

    Fig.12.The temperature contour of a tube to pipe set-on-type joint.

    4.2.Out-of-plane distortion analysis of tube to pipe welding

    The transient out-of-plane distortion simulated during multipass weld deposition of the tube to pipe welding at a given weld sequence-1 is shown in Fig.13.From the figure,it is observed that an increase in the number of weld passes increases the out-ofplane distortion in both directions(x and z)of the tube.This is primarily because an increase of quantum of weld metal deposition results in the accumulation of heat input giving rise to weld metal contraction.It is also observed that out-of-plane distortion in the x direction is higher than that in the z direction because the overall heat input in x direction is higher than the z direction.The excess heat input in the x direction is contributed by the weld start and stop points of all the three weld passes,whereas in the z direction there is no such start and stop points in Sequence-1.Similarly,the out-of-plane distortion in the z direction is the effect of time lag associated with the welding direction between different weld passes.Therefore,it is inferred that out-of-plane distortion in the z direction is primarily influenced by the welding direction,whereas distortion in the x direction is influenced by heat input due to the arc start and stop points in the same location.

    Fig.13.The effect of weld Sequence-1 on out-of-plane distortion under multipass weld deposition of tube to pipe welding.

    In order to establish a clearer understanding about the effect of welding sequence on out-of-plane distortion induced during tube to pipe welding,the trends of tube movements around the circumferential welding of each weld pass were analyzed.In Fig.11, it is observed that there are two different peaks in each weld pass. These peaks indicate the out-of-plane distortion behavior during heating and cooling cycles of circumferential fillet welding of the tube to pipe.Points A,E and I are the start points of each weld pass.Curves A-B,E-F and I-J represent the heating cycle of weld deposition from 0°to 180°of weld passes 1,2 and 3,respectively, and corresponding cooling cycles are B-C,F-G and J-K respectively.During circumferential welding,the curve(A-B,E-F and I-J) moving to a positive direction indicates that expansion takes place during the heating cycle.The curve(B-C,F-G and J-K)moving to the negative direction indicates the contraction that takes place during the cooling cycle.After that,the arc moves to the opposite direction(180°-0°)during the circumferential welding.Curves C-D,G-H and K-L and D-E,H-I and L-M are the heating and cooling cycles from 180°to 0°of weld passes 1,2 and 3 respectively. As in the case of weld deposition from 0°to 180°,it is observed that the expansion and contraction take place during heating and cooling cycles respectively.However,the magnitude of distortion generated in the direction 180°-0°is higher in comparison to that of the direction 0°-180°.This is due to the effect of the time lag involved in the welding direction of weld sequence-1.

    The effect of weld sequence-2 on out-of-plane distortion induced during multipass weld deposition of the tube to pipe fabrication is shown in Fig.14.It is observed that the increase of weld passes increases the distortion in both directions.In welding sequence-2,weld pass 2 showed a negative distortion trend during heating E-F when the weld deposition was from 0°to 180°.This is in contrast to Sequence-1 of weld pass 2.This behavior primarily occurred due to the shifting of the starting point of the arc to 180°. In Fig.14,it is further observed that a little shift in distortion trend is observed in points G and K because of the time lag effects due to the change in welding direction.It can also be noted that the out-ofplane distortion induced in the x direction is lower than the z direction especially at weld pass 3,which is in contrast to sequence-1 (Fig.13)due to the effect of the same weld direction involved in weld passes 2 and 3 than the effect of the arc start and stop points. This behavior resulted in the reduction of out-of-plane distortion in both x and z directions compared to Sequence-1.

    The effect of weld sequence-3 on out-of-plane distortion induced during multi-pass weld deposition of the tube to pipe fabrication is shown in Fig.15.Among the three weld sequences, sequence-3 induces the lowest distortion in both directions.This is primarily due to the balancing effect of heat input due to the start and stop points and welding direction.In sequence 3,the welding direction of each weld pass is entirely different from each other and accordingly balances the heat input during circumferential welding of the tube results in the reduction of out-of-plane distortion.Based on the above discussions,it is well understood that the welding sequence plays a major role in minimizing the out-of-plane distortion induced during circumferential welding of the tube to pipe component.

    Fig.14.The effect of weld Sequence-2 on out-of-plane distortion under multipass weld deposition of tube to pipe welding.

    Fig.15.The effect of weld Sequence-3 on out-of-plane distortion under multipass weld deposition of tube to pipe welding.

    Fig.16.The FEM-simulated displacement contour for three weld sequences.

    After weld metal deposition at the 1500th second(at this time, the weld metal deposition reached to room temperature ensures the complete contraction takes place during the cooling cycle),the steady state contours of different weld sequences 1,2 and 3 are shown in Fig.16(a)-(c).This contour shows the out-of-planedistortion occurred on the tube to pipe component after completion of all weld passes.In Fig.16,it is observed that the tube movements due to circumferential welding under different weld sequences are entirely different from each other because of the change of welding direction and weld start/stop points as explained earlier.

    The simulated and experimentally measured out-of-plane distortion data are presented in Table 4.The steady state out-ofplane distortion magnitudes of FEM and experimental analyses are validated and are shown in Figs.17 and 18.The FEM predicted and experimentally measured values of out-of-plane distortion during tube to pipe welding under different weld sequences and directions are in agreement with the deviation of around 15%. Compared to thermal analysis,mechanical analysis has 5%more deviation.This could be attributed to the manual error involved in the distortion measurement.Based on the FEM and experimental results,the suitable welding sequences and directions were established for the tube to pipe component.

    Table 4The experimental validation of FEM simulated out-of-plane distortion.

    5.Conclusion

    The FEM model was used to study the out-of-plane distortion during multi-pass circumferential fillet welding of the tube to pipe component under three different welding sequences and directions.The FEM model was validated through experiments.Based on the study,the control of out-of-plane distortion was explored and the following conclusions were made.

    1)The thermal analysis was carried out using thermal cycles captured at two locations around the tube and it was observed that the distortion trend was influenced by thermal cycle.The change in welding direction causes the signi ficant change in thermal cycle thereby out-of-plane distortion trend is changed in case of circumferential fillet welding.

    2)The out-of-plane distortion increases when all the three weld passes follow the same circumferential direction around the tube.This is due to the combined effect of heat input due to the same start and stop points and welding direction of the arc. However,the out-of-plane distortion decreases when all the three weld passes follow different circumferential directions to each other around the tube.This is due to the balancing of heat effects of weld passes due to change in welding direction.

    3)The multi-pass circumferential fillet welding of the tube to pipe induces out-of-plane distortion of the tube.The heat input due to weld start and stop points in fluences the out-of-plane distortion in the x direction.The welding direction and associated time lag effects between different weld passes in fluence the out-of-plane distortion in the z direction.

    4)Experimental results con firm the mechanism of out-of-plane distortion explored by FEM analysis,hence the suitable welding procedure was established for the tube to pipe component.

    Fig.17.The steady state distortion magnitudes of FEM and experimental results in x direction.

    Fig.18.The steady state distortion magnitudes of FEM and experimental results in z direction.

    Acknowledgment

    The authors wish to express their sincere thanks to M/s Bharat Heavy Electricals Ltd,Tiruchirappalli,India for providing permission to carry out this research work at their works.

    [1]Chul KY.Production mechanism of out-of-plane deformation in fillet welding. Trans JWRI 1998;27(2).

    [2]Jung GH,Tsai CL.Fundamental studies on the effect of distortion control plans on angular distortion in fillet welded T-joints.Weld J 2004:213s-23s.

    [3]Rui W,Sherif R,Hisashi S,Hidekazu M,Zhang J.Numerical and experimental investigations on welding deformation.Trans JWRI 2008:79-90.

    [4]Wang J,Rashed S,Murakawa H,Luo Y.Numerical production and mitigation of out-of-plane welding distortion in ship panel structure by elastic FE analysis. Mar Struct 2013;34:135-55.

    [5]Vetriselvan R,Raju N,Ravichandran G,Suresh S.Analysis of angular distortion in stainless steel T joint using FEM.WRI J 2011;32:29-38.

    [6]Deng D,Murakawa H,Liang W.Numerical simulation of welding distortion in large structures.Comput Methods Appl Mech Eng 2007:4613-27.

    [7]Ravichandran G,Raghupathy VP,Ganesan N,Krishnakumar R.Prediction of axis shift distortion during circumferential welding of thin pipes using the finite element method.Weld Res Suppl 1997:39s-55s.

    [8]Dar NU,Qureshi EM,Hammouda MMI.Analysis of weld induced residual stresses and distortions in thin walled cylinders.J Mech Sci Technol 2009;23: 1118-31.

    [9]Yu C-L,Chen Z,Wang J,Yan S,Yang L-C.Effect of weld reinforcement on axial plastic buckling of welded steel cylindrical shells.J Zhejiang Univ Sci A 2012;2012(2).ISSN 1862-1775(online).

    [10]Sagalevich VM,Shvetsov VA,Baumana NE.Distortion in welding circumferential seams in thin walled shells.In:Welding production.Cambridge:The Welding Institute;1970.

    [11]Vetriselvan R,Sudharsanam V,Raju N,Ravichandran G.Analysis of angular distortion in header to nipple welding using FEM.In:International symposium on joining of materials.Tiruchirappalli:IWS;2012.

    [12]Vetriselvan R,Sathiya P,Ravichandran G.Characterization of transient out-ofplane distortion of nipple welding with header component.J Manuf Process 2015:67-72.

    [13]Sysweld.Manual.France:ESI;2014.

    [14]Dialami N,Chiumenti M,Agelet de Scracibar C,Cervera M.A strategy for the numerical simulation of metal deposition and multipass welding process.In: International conference on computational plasticity.CIMNE;2009. Barcelona.

    [15]Goldak J,Chakravarti A,Bibby M.A double ellipsoidal finite element model for welding heat sources.IIW Doc 212-603-85.1985.

    *Corresponding author.

    E-mail addresses:vetri@bheltry.co.in,vetriselvan2003@gmail.com (R.Vetriselvan).

    Peer review under responsibility of China Ordnance Society.

    建设人人有责人人尽责人人享有的| 久久鲁丝午夜福利片| 飞空精品影院首页| 日韩免费高清中文字幕av| 午夜久久久在线观看| 一级毛片 在线播放| 亚洲激情五月婷婷啪啪| 日本免费在线观看一区| 成人免费观看视频高清| 久久毛片免费看一区二区三区| 欧美精品高潮呻吟av久久| 久久久久国产精品人妻一区二区| 午夜日韩欧美国产| 91久久精品国产一区二区三区| 亚洲国产av新网站| 亚洲第一区二区三区不卡| 亚洲国产日韩一区二区| 欧美变态另类bdsm刘玥| 日韩一区二区三区影片| 午夜免费鲁丝| 啦啦啦在线观看免费高清www| freevideosex欧美| 亚洲熟女精品中文字幕| 国产男女内射视频| 国产男女内射视频| 日韩三级伦理在线观看| 又粗又硬又长又爽又黄的视频| 亚洲欧洲日产国产| 校园人妻丝袜中文字幕| 卡戴珊不雅视频在线播放| av不卡在线播放| 亚洲av免费高清在线观看| 秋霞在线观看毛片| videos熟女内射| 成人午夜精彩视频在线观看| 波野结衣二区三区在线| 1024视频免费在线观看| 免费不卡的大黄色大毛片视频在线观看| 校园人妻丝袜中文字幕| 精品久久久精品久久久| 人妻系列 视频| 久久精品国产a三级三级三级| 在线观看一区二区三区激情| 国产又色又爽无遮挡免| 亚洲经典国产精华液单| 男女边摸边吃奶| 最新中文字幕久久久久| 国产亚洲欧美精品永久| 国精品久久久久久国模美| 曰老女人黄片| 久热久热在线精品观看| 成人午夜精彩视频在线观看| 国产成人精品福利久久| 午夜福利网站1000一区二区三区| 99香蕉大伊视频| 熟女少妇亚洲综合色aaa.| 最近中文字幕高清免费大全6| 亚洲四区av| 超碰97精品在线观看| 狠狠精品人妻久久久久久综合| 亚洲国产看品久久| 一区福利在线观看| 国产精品一区二区在线观看99| 亚洲人成电影观看| 国产国语露脸激情在线看| 大香蕉久久成人网| 亚洲一区二区三区欧美精品| 如日韩欧美国产精品一区二区三区| 1024视频免费在线观看| 一级,二级,三级黄色视频| av福利片在线| 国产欧美日韩综合在线一区二区| 中文字幕人妻熟女乱码| 中国三级夫妇交换| 亚洲国产欧美在线一区| 日韩成人av中文字幕在线观看| 国产在视频线精品| 老司机影院成人| freevideosex欧美| 美女脱内裤让男人舔精品视频| 欧美精品人与动牲交sv欧美| 18在线观看网站| 国产成人91sexporn| 亚洲av成人精品一二三区| 高清av免费在线| 免费观看无遮挡的男女| 水蜜桃什么品种好| 女的被弄到高潮叫床怎么办| 国产淫语在线视频| 日韩欧美一区视频在线观看| 亚洲在久久综合| videosex国产| 亚洲欧美一区二区三区黑人 | 亚洲色图 男人天堂 中文字幕| 99热网站在线观看| 青春草亚洲视频在线观看| 丝袜美足系列| 在线 av 中文字幕| 国产日韩欧美视频二区| 久久精品国产亚洲av天美| h视频一区二区三区| 欧美激情极品国产一区二区三区| 人成视频在线观看免费观看| 人人妻人人澡人人看| 成人18禁高潮啪啪吃奶动态图| 午夜日韩欧美国产| 97精品久久久久久久久久精品| 女人被躁到高潮嗷嗷叫费观| 国产色婷婷99| 精品人妻熟女毛片av久久网站| 国产av精品麻豆| 久久精品亚洲av国产电影网| 成人国产麻豆网| 欧美精品av麻豆av| 国产免费一区二区三区四区乱码| 午夜老司机福利剧场| 免费看av在线观看网站| av在线老鸭窝| 国产毛片在线视频| 成年美女黄网站色视频大全免费| 大香蕉久久成人网| 精品午夜福利在线看| 国产av国产精品国产| 下体分泌物呈黄色| 女人被躁到高潮嗷嗷叫费观| 两个人免费观看高清视频| 国产精品二区激情视频| 99国产综合亚洲精品| 亚洲美女黄色视频免费看| 日韩三级伦理在线观看| 午夜久久久在线观看| 在线观看免费高清a一片| 亚洲成人一二三区av| 亚洲第一av免费看| 99re6热这里在线精品视频| 国产成人精品久久二区二区91 | 国产白丝娇喘喷水9色精品| 交换朋友夫妻互换小说| 看非洲黑人一级黄片| 男女边吃奶边做爰视频| 在线观看免费高清a一片| 久久人人爽人人片av| 纯流量卡能插随身wifi吗| 一区福利在线观看| 久久女婷五月综合色啪小说| 超色免费av| 大码成人一级视频| 青青草视频在线视频观看| 99久国产av精品国产电影| 亚洲av综合色区一区| 另类精品久久| 男人添女人高潮全过程视频| 欧美最新免费一区二区三区| 街头女战士在线观看网站| 99久久中文字幕三级久久日本| 黄片播放在线免费| 美女国产高潮福利片在线看| 久久久久久人人人人人| 午夜福利视频在线观看免费| 久久精品国产亚洲av涩爱| 国产成人精品无人区| 人人妻人人澡人人看| 色视频在线一区二区三区| a 毛片基地| 亚洲av在线观看美女高潮| 亚洲人成网站在线观看播放| 亚洲成人av在线免费| 汤姆久久久久久久影院中文字幕| 中文字幕色久视频| 国产精品女同一区二区软件| 久久久久久免费高清国产稀缺| 国产不卡av网站在线观看| 精品福利永久在线观看| 成人国产av品久久久| 国产男女超爽视频在线观看| 少妇被粗大猛烈的视频| 综合色丁香网| 国产无遮挡羞羞视频在线观看| 午夜福利,免费看| 丝袜人妻中文字幕| 激情五月婷婷亚洲| 日韩 亚洲 欧美在线| 成人午夜精彩视频在线观看| 高清黄色对白视频在线免费看| 亚洲一级一片aⅴ在线观看| 精品久久久精品久久久| 精品亚洲乱码少妇综合久久| 九草在线视频观看| 一区二区av电影网| 一二三四中文在线观看免费高清| 性色avwww在线观看| 日日啪夜夜爽| 黄色配什么色好看| av免费观看日本| 久久久久久久久免费视频了| 青春草国产在线视频| 午夜福利,免费看| 少妇被粗大猛烈的视频| av在线播放精品| 在线 av 中文字幕| 天天影视国产精品| 亚洲国产精品一区二区三区在线| 国产成人91sexporn| 色婷婷久久久亚洲欧美| 亚洲第一av免费看| 日韩av不卡免费在线播放| 国产日韩一区二区三区精品不卡| 男人爽女人下面视频在线观看| 日韩视频在线欧美| 午夜福利网站1000一区二区三区| 亚洲国产欧美在线一区| 女性被躁到高潮视频| 久久精品熟女亚洲av麻豆精品| 交换朋友夫妻互换小说| 91久久精品国产一区二区三区| 少妇熟女欧美另类| 欧美激情 高清一区二区三区| 男女啪啪激烈高潮av片| 国产亚洲午夜精品一区二区久久| 一级a爱视频在线免费观看| 欧美老熟妇乱子伦牲交| 婷婷色麻豆天堂久久| 男女高潮啪啪啪动态图| 精品少妇一区二区三区视频日本电影 | 蜜桃在线观看..| 国产成人精品无人区| 免费不卡的大黄色大毛片视频在线观看| av国产久精品久网站免费入址| 婷婷成人精品国产| 美女中出高潮动态图| 亚洲国产精品一区二区三区在线| 日韩大片免费观看网站| 性色avwww在线观看| 亚洲国产精品国产精品| 国产免费现黄频在线看| 成人手机av| 亚洲精品国产色婷婷电影| 最近2019中文字幕mv第一页| 黄片小视频在线播放| 久久99精品国语久久久| 男女高潮啪啪啪动态图| 黄网站色视频无遮挡免费观看| 日韩在线高清观看一区二区三区| 男人操女人黄网站| av国产精品久久久久影院| 精品99又大又爽又粗少妇毛片| 人体艺术视频欧美日本| 亚洲av.av天堂| 国产精品一区二区在线不卡| 99久国产av精品国产电影| av在线app专区| 精品酒店卫生间| av福利片在线| 久久韩国三级中文字幕| 亚洲国产欧美日韩在线播放| 成年人免费黄色播放视频| 日韩av在线免费看完整版不卡| 高清av免费在线| 最近中文字幕高清免费大全6| 黑人巨大精品欧美一区二区蜜桃| 18+在线观看网站| 亚洲人成电影观看| 黄色怎么调成土黄色| 美女视频免费永久观看网站| 久久久久国产网址| 精品国产乱码久久久久久小说| 午夜日韩欧美国产| 精品国产一区二区三区四区第35| 亚洲激情五月婷婷啪啪| 9色porny在线观看| 香蕉精品网在线| 亚洲成色77777| 日韩大片免费观看网站| 水蜜桃什么品种好| 免费在线观看完整版高清| 日本午夜av视频| 精品国产国语对白av| 国产野战对白在线观看| 欧美激情高清一区二区三区 | 99热网站在线观看| av国产久精品久网站免费入址| 巨乳人妻的诱惑在线观看| 欧美xxⅹ黑人| 国产亚洲午夜精品一区二区久久| 午夜激情av网站| 亚洲第一青青草原| 免费观看a级毛片全部| 丝袜美足系列| 久久狼人影院| 成人手机av| 毛片一级片免费看久久久久| 美女午夜性视频免费| 日韩av不卡免费在线播放| 在线观看国产h片| av不卡在线播放| 国产淫语在线视频| 日韩人妻精品一区2区三区| 成人漫画全彩无遮挡| 男人舔女人的私密视频| 在线天堂最新版资源| 一本—道久久a久久精品蜜桃钙片| 最近2019中文字幕mv第一页| 深夜精品福利| 亚洲激情五月婷婷啪啪| 一级片'在线观看视频| 水蜜桃什么品种好| a级片在线免费高清观看视频| 婷婷色麻豆天堂久久| 久久青草综合色| 精品亚洲乱码少妇综合久久| 一级片免费观看大全| 9色porny在线观看| 伊人久久国产一区二区| 国产成人精品久久久久久| 国产精品欧美亚洲77777| 国产老妇伦熟女老妇高清| 美女午夜性视频免费| 国产极品天堂在线| 久久99热这里只频精品6学生| av在线app专区| 精品少妇久久久久久888优播| 久久精品久久久久久久性| 人人妻人人澡人人爽人人夜夜| 777米奇影视久久| 韩国av在线不卡| 国产在线视频一区二区| 亚洲欧洲国产日韩| 国产片特级美女逼逼视频| 另类亚洲欧美激情| 搡女人真爽免费视频火全软件| 五月天丁香电影| 亚洲国产日韩一区二区| 精品久久蜜臀av无| 亚洲五月色婷婷综合| 十八禁网站网址无遮挡| 99久久中文字幕三级久久日本| 亚洲欧美精品综合一区二区三区 | a级毛片在线看网站| 精品少妇久久久久久888优播| 亚洲av在线观看美女高潮| 成人毛片60女人毛片免费| 综合色丁香网| 伊人久久国产一区二区| 纵有疾风起免费观看全集完整版| 国产一区二区 视频在线| 啦啦啦啦在线视频资源| 久久鲁丝午夜福利片| 国产97色在线日韩免费| 美女国产高潮福利片在线看| 国产精品二区激情视频| 爱豆传媒免费全集在线观看| 亚洲精品第二区| 波野结衣二区三区在线| 777米奇影视久久| 久久久久久免费高清国产稀缺| 午夜福利在线免费观看网站| 免费久久久久久久精品成人欧美视频| 久久久精品区二区三区| 免费高清在线观看日韩| 亚洲av福利一区| 国产亚洲一区二区精品| 免费高清在线观看日韩| 国产亚洲最大av| 国产欧美日韩一区二区三区在线| 黄频高清免费视频| 一级毛片电影观看| www.自偷自拍.com| 天天躁夜夜躁狠狠久久av| 啦啦啦中文免费视频观看日本| 免费黄频网站在线观看国产| 两性夫妻黄色片| 一区二区三区乱码不卡18| 国产av一区二区精品久久| 日韩伦理黄色片| 欧美日韩av久久| 国产白丝娇喘喷水9色精品| 久久久精品94久久精品| 黄色 视频免费看| 日本-黄色视频高清免费观看| 制服人妻中文乱码| 女的被弄到高潮叫床怎么办| 欧美日韩综合久久久久久| 多毛熟女@视频| 狠狠精品人妻久久久久久综合| 日韩一本色道免费dvd| 亚洲av男天堂| 热re99久久精品国产66热6| 高清黄色对白视频在线免费看| 国产av国产精品国产| 18禁国产床啪视频网站| 亚洲内射少妇av| 日本91视频免费播放| 免费播放大片免费观看视频在线观看| 天天躁日日躁夜夜躁夜夜| 亚洲精华国产精华液的使用体验| 久久久久久伊人网av| 狂野欧美激情性bbbbbb| 女性生殖器流出的白浆| 国产亚洲午夜精品一区二区久久| 纯流量卡能插随身wifi吗| 亚洲精品视频女| 夜夜骑夜夜射夜夜干| 中文字幕人妻丝袜一区二区 | 一级片'在线观看视频| 国产欧美日韩一区二区三区在线| 亚洲av欧美aⅴ国产| 久久午夜福利片| 中文字幕人妻丝袜一区二区 | 亚洲欧洲日产国产| 国产精品蜜桃在线观看| 成人毛片60女人毛片免费| 美女福利国产在线| 国产成人一区二区在线| 又粗又硬又长又爽又黄的视频| 狠狠精品人妻久久久久久综合| 男的添女的下面高潮视频| 青青草视频在线视频观看| 国产淫语在线视频| 亚洲av成人精品一二三区| 久久人人爽人人片av| 国产成人精品久久久久久| 嫩草影院入口| 亚洲国产日韩一区二区| 欧美中文综合在线视频| 欧美成人午夜免费资源| 欧美日韩视频精品一区| 99久久精品国产国产毛片| av有码第一页| 国产亚洲最大av| 国产免费又黄又爽又色| 日韩不卡一区二区三区视频在线| 国产精品久久久久久精品古装| 最近中文字幕2019免费版| 在线观看免费高清a一片| 男女无遮挡免费网站观看| 国产成人精品婷婷| 久久精品国产鲁丝片午夜精品| 最新中文字幕久久久久| 午夜精品国产一区二区电影| 亚洲色图 男人天堂 中文字幕| 久久久精品区二区三区| 国产精品偷伦视频观看了| 久久精品夜色国产| 成人毛片60女人毛片免费| 咕卡用的链子| 国产一区二区激情短视频 | 又黄又粗又硬又大视频| 国产高清不卡午夜福利| 免费观看无遮挡的男女| 夫妻性生交免费视频一级片| 国产av一区二区精品久久| 少妇人妻久久综合中文| 精品少妇一区二区三区视频日本电影 | 黄色毛片三级朝国网站| 亚洲精华国产精华液的使用体验| 日韩中文字幕视频在线看片| 999久久久国产精品视频| 最新的欧美精品一区二区| 人妻少妇偷人精品九色| 亚洲国产精品一区二区三区在线| 好男人视频免费观看在线| 亚洲男人天堂网一区| 大码成人一级视频| 亚洲少妇的诱惑av| 男人操女人黄网站| 最近2019中文字幕mv第一页| 寂寞人妻少妇视频99o| 国产成人a∨麻豆精品| 中文字幕人妻丝袜一区二区 | 久久久久久久久免费视频了| 久久精品熟女亚洲av麻豆精品| 五月伊人婷婷丁香| 色婷婷久久久亚洲欧美| 免费高清在线观看日韩| 校园人妻丝袜中文字幕| 赤兔流量卡办理| 色哟哟·www| 国产日韩欧美在线精品| 久久久久久伊人网av| 国产爽快片一区二区三区| 国产精品成人在线| 国产精品久久久久久精品电影小说| 成人亚洲精品一区在线观看| 天美传媒精品一区二区| 成人午夜精彩视频在线观看| 777米奇影视久久| 精品亚洲成a人片在线观看| 精品国产乱码久久久久久小说| 亚洲精品在线美女| 日产精品乱码卡一卡2卡三| 亚洲四区av| 久久久久久久久久人人人人人人| 久久精品久久精品一区二区三区| a级片在线免费高清观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品国产鲁丝片午夜精品| 国产熟女午夜一区二区三区| 建设人人有责人人尽责人人享有的| av网站免费在线观看视频| 春色校园在线视频观看| 看十八女毛片水多多多| 久久久久久久精品精品| 在线精品无人区一区二区三| 久久午夜综合久久蜜桃| 日韩电影二区| 在线观看三级黄色| 极品少妇高潮喷水抽搐| av福利片在线| 婷婷成人精品国产| 制服诱惑二区| 黄网站色视频无遮挡免费观看| 亚洲av国产av综合av卡| 日本欧美国产在线视频| 黄色一级大片看看| 国产精品欧美亚洲77777| 秋霞伦理黄片| 高清av免费在线| 久久国产精品男人的天堂亚洲| 亚洲伊人久久精品综合| 高清在线视频一区二区三区| 亚洲精品,欧美精品| 啦啦啦啦在线视频资源| 免费人妻精品一区二区三区视频| 最近2019中文字幕mv第一页| 男女下面插进去视频免费观看| 春色校园在线视频观看| 丝袜美腿诱惑在线| 美国免费a级毛片| 叶爱在线成人免费视频播放| 精品国产国语对白av| 国产白丝娇喘喷水9色精品| 久久青草综合色| 国产欧美亚洲国产| 精品少妇久久久久久888优播| 黄片小视频在线播放| 国产精品一区二区在线不卡| 免费看不卡的av| 夫妻午夜视频| 成年女人在线观看亚洲视频| 国产精品亚洲av一区麻豆 | 香蕉精品网在线| 亚洲色图 男人天堂 中文字幕| 国产一级毛片在线| 街头女战士在线观看网站| 建设人人有责人人尽责人人享有的| 午夜av观看不卡| 蜜桃国产av成人99| 欧美日韩精品成人综合77777| 女人被躁到高潮嗷嗷叫费观| 精品国产乱码久久久久久小说| 狠狠精品人妻久久久久久综合| 久久久久精品久久久久真实原创| 午夜精品国产一区二区电影| 哪个播放器可以免费观看大片| 丝袜在线中文字幕| 26uuu在线亚洲综合色| 久久久久久久国产电影| 亚洲欧美成人综合另类久久久| 伊人久久大香线蕉亚洲五| 日韩一卡2卡3卡4卡2021年| 国产亚洲精品第一综合不卡| 国产一区亚洲一区在线观看| 久久久久久久亚洲中文字幕| 国产精品人妻久久久影院| 男女无遮挡免费网站观看| 亚洲欧洲日产国产| 1024香蕉在线观看| 黄片播放在线免费| www.av在线官网国产| 日韩成人av中文字幕在线观看| 不卡视频在线观看欧美| 啦啦啦在线免费观看视频4| 成人亚洲欧美一区二区av| 日本欧美国产在线视频| 边亲边吃奶的免费视频| 亚洲欧美一区二区三区国产| 免费人妻精品一区二区三区视频| 婷婷成人精品国产| 亚洲欧美成人综合另类久久久| 十八禁高潮呻吟视频| 久久久久视频综合| 精品一品国产午夜福利视频| 久久久久人妻精品一区果冻| 久久精品国产亚洲av涩爱| 男男h啪啪无遮挡| 你懂的网址亚洲精品在线观看| www日本在线高清视频| 日韩一本色道免费dvd| 久久久精品区二区三区| 哪个播放器可以免费观看大片| 999久久久国产精品视频| 七月丁香在线播放| 国产色婷婷99| 黑人猛操日本美女一级片| 亚洲,欧美精品.| 日日摸夜夜添夜夜爱| 夫妻性生交免费视频一级片| 夜夜骑夜夜射夜夜干| 美女大奶头黄色视频| 国产精品一二三区在线看| av免费观看日本| 亚洲国产精品国产精品| 中文精品一卡2卡3卡4更新| 一区二区日韩欧美中文字幕| 亚洲情色 制服丝袜| 久久久久久久久久久免费av| 精品卡一卡二卡四卡免费| 在线看a的网站| 欧美成人午夜精品| 国产欧美日韩综合在线一区二区| 丁香六月天网| 男女边吃奶边做爰视频| 深夜精品福利| 91精品伊人久久大香线蕉| 免费观看在线日韩|