• <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麻豆精品| 国产老妇伦熟女老妇高清| 日韩中字成人| 大香蕉久久成人网| 久久久午夜欧美精品| 久久久久久久久久久久大奶| 精品少妇黑人巨大在线播放| 赤兔流量卡办理| 日日啪夜夜爽| 另类精品久久| 国产精品99久久99久久久不卡 | 精品少妇久久久久久888优播| 最新中文字幕久久久久| 欧美日韩精品成人综合77777| 永久免费av网站大全| 精品熟女少妇av免费看| 国产高清不卡午夜福利| 色网站视频免费| 性高湖久久久久久久久免费观看| 亚洲国产日韩一区二区| 久久av网站| 久久av网站| 国产精品.久久久| 日韩熟女老妇一区二区性免费视频| 亚洲性久久影院| 亚洲精品美女久久av网站| 国产探花极品一区二区| 亚洲精品一二三| 一级爰片在线观看| 天天影视国产精品| 精品亚洲乱码少妇综合久久| 国产亚洲精品第一综合不卡 | 亚洲精品成人av观看孕妇| 日韩亚洲欧美综合| 国产男人的电影天堂91| 亚洲av中文av极速乱| 国产极品天堂在线| 日韩,欧美,国产一区二区三区| 久久久久久久亚洲中文字幕| 亚洲综合精品二区| 国产成人a∨麻豆精品| 91精品国产国语对白视频| 亚洲高清免费不卡视频| av国产久精品久网站免费入址| 中文字幕av电影在线播放| 91精品国产国语对白视频| 国产老妇伦熟女老妇高清| 青青草视频在线视频观看| 综合色丁香网| 久久国内精品自在自线图片| 一区二区三区免费毛片| 在线观看三级黄色| 国产精品国产三级专区第一集| 欧美人与性动交α欧美精品济南到 | 永久网站在线| 99国产综合亚洲精品| 成人亚洲精品一区在线观看| 精品人妻在线不人妻| www.av在线官网国产| 免费观看在线日韩| av免费观看日本| 色视频在线一区二区三区| 一本久久精品| 少妇精品久久久久久久| 九九在线视频观看精品| 99久久中文字幕三级久久日本| 好男人视频免费观看在线| 高清在线视频一区二区三区| 在线看a的网站| 一边亲一边摸免费视频| 99热全是精品| 十八禁网站网址无遮挡| 看非洲黑人一级黄片| 高清不卡的av网站| 晚上一个人看的免费电影| 国产一区二区三区综合在线观看 | 亚洲精品久久午夜乱码| 成年av动漫网址| 性高湖久久久久久久久免费观看| 成人国语在线视频| av女优亚洲男人天堂| 男女无遮挡免费网站观看| 午夜福利网站1000一区二区三区| 男女高潮啪啪啪动态图| 欧美3d第一页| 亚洲国产日韩一区二区| 女人精品久久久久毛片| av女优亚洲男人天堂| 久久久精品免费免费高清| 亚洲av中文av极速乱| 欧美性感艳星| 久久精品国产亚洲av天美| 国产精品一区二区在线观看99| 精品国产一区二区久久| 亚洲精品日本国产第一区| 精品视频人人做人人爽| 免费高清在线观看视频在线观看| 天美传媒精品一区二区| 久久精品国产a三级三级三级| 插逼视频在线观看| 国产亚洲午夜精品一区二区久久| 国产黄色视频一区二区在线观看| 国产成人一区二区在线| 精品一区二区免费观看| a 毛片基地| 国产午夜精品一二区理论片| 69精品国产乱码久久久| 午夜91福利影院| 在线观看人妻少妇| 男的添女的下面高潮视频| a级毛片免费高清观看在线播放| 99热6这里只有精品| 国产伦理片在线播放av一区| 视频区图区小说| 国产免费一区二区三区四区乱码| 青春草视频在线免费观看| av网站免费在线观看视频| 97在线视频观看| 麻豆精品久久久久久蜜桃| 亚洲久久久国产精品| 最黄视频免费看| 亚洲婷婷狠狠爱综合网| 亚洲三级黄色毛片| 国产午夜精品久久久久久一区二区三区| 成人影院久久| 18禁在线无遮挡免费观看视频| 国产精品 国内视频| 亚洲精品日韩av片在线观看| 亚洲av国产av综合av卡| 飞空精品影院首页| 九色亚洲精品在线播放| 少妇人妻久久综合中文| 精品久久久久久久久亚洲| 人人妻人人澡人人看| 女人精品久久久久毛片| 亚洲欧洲日产国产| 欧美性感艳星| 久久久国产一区二区| 亚洲,一卡二卡三卡| 美女cb高潮喷水在线观看| 高清不卡的av网站| 国产精品久久久久久av不卡| 综合色丁香网| 黄色怎么调成土黄色| 插阴视频在线观看视频| 午夜激情av网站| av线在线观看网站| 精品久久久噜噜| 久久狼人影院| 成人免费观看视频高清| 一级,二级,三级黄色视频| 免费播放大片免费观看视频在线观看| 亚洲国产精品一区三区| 国产在视频线精品| 久久人人爽人人片av| 国产精品麻豆人妻色哟哟久久| 国产在线免费精品| 欧美精品一区二区免费开放| 97在线视频观看| 丝袜脚勾引网站| 日韩视频在线欧美| 亚洲第一区二区三区不卡| 国产精品偷伦视频观看了| 国产午夜精品久久久久久一区二区三区| 国产日韩欧美亚洲二区| 午夜免费男女啪啪视频观看| 国产免费福利视频在线观看| 草草在线视频免费看| 免费观看性生交大片5| 三级国产精品片| 成年女人在线观看亚洲视频| 大陆偷拍与自拍| 最新中文字幕久久久久| 9色porny在线观看| 久久精品夜色国产| 久久这里有精品视频免费| 国产精品久久久久久久电影| 国产精品秋霞免费鲁丝片| 精品少妇久久久久久888优播| 青春草国产在线视频| 亚洲三级黄色毛片| 天堂中文最新版在线下载| 日韩中文字幕视频在线看片| 美女cb高潮喷水在线观看| 国产精品麻豆人妻色哟哟久久| 伦精品一区二区三区| 亚洲欧美一区二区三区黑人 | 热99国产精品久久久久久7| 久久精品国产a三级三级三级| 看免费成人av毛片| 久久99蜜桃精品久久| 欧美另类一区| 亚洲国产色片| 高清在线视频一区二区三区| 久久久久国产网址| 中文欧美无线码| 久久综合国产亚洲精品| 青春草国产在线视频| 考比视频在线观看| 高清黄色对白视频在线免费看| 一本色道久久久久久精品综合| 国产精品久久久久久av不卡| 欧美成人精品欧美一级黄| 日韩在线高清观看一区二区三区| 母亲3免费完整高清在线观看 | 如何舔出高潮| 国产成人一区二区在线| 日本与韩国留学比较| 制服人妻中文乱码| 久久精品夜色国产| 日韩欧美一区视频在线观看| 免费观看av网站的网址| 精品国产一区二区久久| 午夜福利视频精品| 国产深夜福利视频在线观看| 精品人妻熟女毛片av久久网站| 在线观看免费日韩欧美大片 | 亚洲欧美成人精品一区二区| 日产精品乱码卡一卡2卡三| 久久久国产一区二区| 老司机影院成人| 成人综合一区亚洲| 国产永久视频网站| 九色成人免费人妻av| 日韩亚洲欧美综合| 丝袜喷水一区| 欧美日韩综合久久久久久| 国产午夜精品久久久久久一区二区三区| 视频中文字幕在线观看| 日本黄色片子视频| kizo精华| 精品熟女少妇av免费看| 国产精品秋霞免费鲁丝片| 91精品伊人久久大香线蕉| 国产69精品久久久久777片| 日本欧美视频一区| 国产有黄有色有爽视频| 国产精品国产三级国产av玫瑰| 精品卡一卡二卡四卡免费| 婷婷色av中文字幕| 一级毛片aaaaaa免费看小| 91精品一卡2卡3卡4卡| 国产永久视频网站| 精品久久国产蜜桃| 欧美精品亚洲一区二区| 成年女人在线观看亚洲视频| 欧美亚洲 丝袜 人妻 在线| 国产成人av激情在线播放 | 精品久久久噜噜| 国产成人一区二区在线| 亚洲精品乱码久久久久久按摩| 视频区图区小说| 夫妻午夜视频| 青青草视频在线视频观看| 亚洲高清免费不卡视频| 成年av动漫网址| 高清视频免费观看一区二区| 99视频精品全部免费 在线| 午夜av观看不卡| 成人国语在线视频| 在线观看免费高清a一片| 亚洲成人一二三区av| 久久毛片免费看一区二区三区| av视频免费观看在线观看| av在线老鸭窝| 久久久久久久久久成人| 日本免费在线观看一区| 九色成人免费人妻av| 最近2019中文字幕mv第一页| 国产探花极品一区二区| 亚洲国产av新网站| 我要看黄色一级片免费的| 丝袜喷水一区| 天天躁夜夜躁狠狠久久av| 午夜久久久在线观看| 欧美日韩av久久| 80岁老熟妇乱子伦牲交| 欧美精品一区二区免费开放| 美女内射精品一级片tv| 免费人成在线观看视频色| 在线看a的网站| 日韩av免费高清视频| 亚洲四区av| 国产av码专区亚洲av| 九九久久精品国产亚洲av麻豆| 久久久久人妻精品一区果冻| 亚洲av不卡在线观看| 精品少妇黑人巨大在线播放| 日本vs欧美在线观看视频| 午夜91福利影院| 亚洲国产精品国产精品| 寂寞人妻少妇视频99o| 高清av免费在线| 亚洲欧美一区二区三区国产| 日韩在线高清观看一区二区三区| 国产亚洲欧美精品永久| 精品人妻一区二区三区麻豆| 免费播放大片免费观看视频在线观看| 国产精品国产三级国产专区5o| 午夜福利在线观看免费完整高清在| 国产黄片视频在线免费观看| a级片在线免费高清观看视频| 国产成人免费观看mmmm| 丰满饥渴人妻一区二区三| 交换朋友夫妻互换小说| 免费观看的影片在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久狼人影院| freevideosex欧美| 激情五月婷婷亚洲| 欧美3d第一页| 日韩强制内射视频| 看非洲黑人一级黄片| 在线观看免费视频网站a站| 在线天堂最新版资源| 国产精品不卡视频一区二区| 22中文网久久字幕| 日韩三级伦理在线观看| 精品久久久精品久久久| av在线app专区| 五月开心婷婷网| 国产精品国产三级专区第一集| 桃花免费在线播放| 99九九在线精品视频| 18禁观看日本| 九草在线视频观看| 26uuu在线亚洲综合色| 中文字幕最新亚洲高清| 久久精品夜色国产| 丁香六月天网| 满18在线观看网站| 国产精品免费大片| 国产精品一二三区在线看| 亚洲少妇的诱惑av| 交换朋友夫妻互换小说| av天堂久久9| 久久久a久久爽久久v久久| 天天躁夜夜躁狠狠久久av| av国产久精品久网站免费入址| 人人妻人人爽人人添夜夜欢视频| 国产男女内射视频| 亚洲第一av免费看| www.色视频.com| 亚洲美女搞黄在线观看| 亚洲美女视频黄频| 亚洲精品国产av蜜桃| 91aial.com中文字幕在线观看| 春色校园在线视频观看| 成人亚洲欧美一区二区av| 一区二区三区四区激情视频| 伊人久久精品亚洲午夜| 国产综合精华液| h视频一区二区三区| 免费播放大片免费观看视频在线观看| 亚洲av不卡在线观看| 久久久午夜欧美精品| 亚洲欧美中文字幕日韩二区| 成人毛片60女人毛片免费| 美女cb高潮喷水在线观看| 精品国产一区二区三区久久久樱花| 久久 成人 亚洲| 91精品三级在线观看| 黄片播放在线免费| 尾随美女入室| 天堂8中文在线网| 久久亚洲国产成人精品v| 中文字幕最新亚洲高清| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久精品国产亚洲av天美| 免费观看a级毛片全部| 狂野欧美激情性bbbbbb| 伊人亚洲综合成人网| 老司机影院成人| 亚洲国产精品一区二区三区在线| 亚洲欧美中文字幕日韩二区| 亚洲精品视频女| 丝袜脚勾引网站| 欧美精品人与动牲交sv欧美| 在线免费观看不下载黄p国产| 中文天堂在线官网| 少妇 在线观看| 伊人亚洲综合成人网| 黄色视频在线播放观看不卡| 久久免费观看电影| 少妇 在线观看| 日韩一本色道免费dvd| 欧美日韩国产mv在线观看视频| 欧美97在线视频| 高清av免费在线| 久久这里有精品视频免费| 精品久久国产蜜桃| 少妇的逼水好多| 丰满乱子伦码专区| 80岁老熟妇乱子伦牲交| 99九九线精品视频在线观看视频| 日本黄色日本黄色录像| 精品少妇内射三级| 高清午夜精品一区二区三区| 精品一区二区三区视频在线| 亚洲av国产av综合av卡| 亚洲精品456在线播放app| 美女主播在线视频| 飞空精品影院首页| 国产伦精品一区二区三区视频9| 国产欧美日韩一区二区三区在线 | av黄色大香蕉| 十分钟在线观看高清视频www| 精品亚洲成a人片在线观看| 97超视频在线观看视频| 午夜免费观看性视频| 精品一区在线观看国产| 色5月婷婷丁香| 欧美日韩视频高清一区二区三区二| 日韩,欧美,国产一区二区三区| 天美传媒精品一区二区| 国产精品一区二区三区四区免费观看| 国产精品久久久久久av不卡| av.在线天堂| 免费观看的影片在线观看| 满18在线观看网站| 韩国av在线不卡| 亚洲av中文av极速乱| 最近手机中文字幕大全| 久久鲁丝午夜福利片| 久久97久久精品| 久久女婷五月综合色啪小说| 亚洲av免费高清在线观看| 美女大奶头黄色视频| 少妇人妻精品综合一区二区| 哪个播放器可以免费观看大片| 亚洲国产欧美日韩在线播放| 日韩伦理黄色片| 伊人久久国产一区二区| 国产亚洲欧美精品永久| 日韩一区二区视频免费看| 精品久久国产蜜桃| 高清在线视频一区二区三区| 亚洲美女黄色视频免费看| 人妻少妇偷人精品九色| 午夜福利在线观看免费完整高清在| 丝袜美足系列| 大香蕉久久网| xxxhd国产人妻xxx| 伦精品一区二区三区| 久久午夜综合久久蜜桃| 精品一品国产午夜福利视频| 又大又黄又爽视频免费| 一二三四中文在线观看免费高清| 久久久欧美国产精品| 亚洲国产精品专区欧美| 国产黄色视频一区二区在线观看| 色94色欧美一区二区| 黑人欧美特级aaaaaa片| 亚洲少妇的诱惑av| 亚洲熟女精品中文字幕| 国内精品宾馆在线| 久久久久国产网址| 国产精品三级大全| 亚洲国产毛片av蜜桃av| 精品国产一区二区久久| 女的被弄到高潮叫床怎么办| 日本av免费视频播放| 日本黄色片子视频| 亚洲精品日韩在线中文字幕| 亚洲欧美中文字幕日韩二区| 亚洲人成网站在线观看播放| 国产精品蜜桃在线观看| 日产精品乱码卡一卡2卡三| 久久久久国产精品人妻一区二区| 男女免费视频国产| 免费观看无遮挡的男女| 十八禁网站网址无遮挡| 青青草视频在线视频观看| 亚洲性久久影院| 国产综合精华液| 一个人看视频在线观看www免费| 免费大片18禁| 亚洲性久久影院| 美女xxoo啪啪120秒动态图| 午夜视频国产福利| 大片电影免费在线观看免费| 又黄又爽又刺激的免费视频.| 女人精品久久久久毛片| 老熟女久久久| 成人国产麻豆网| 如何舔出高潮| 91aial.com中文字幕在线观看| 人妻系列 视频| 一级,二级,三级黄色视频| 黄色视频在线播放观看不卡| 亚洲第一区二区三区不卡| 亚洲国产最新在线播放| 免费高清在线观看日韩| 热re99久久国产66热| 欧美日本中文国产一区发布| 色视频在线一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美性感艳星| 18+在线观看网站| av又黄又爽大尺度在线免费看| 中文字幕人妻熟人妻熟丝袜美| 国产成人91sexporn| 老熟女久久久| 色婷婷av一区二区三区视频| 丝瓜视频免费看黄片| 亚洲精品日韩av片在线观看| 韩国高清视频一区二区三区| 国产黄色免费在线视频| 一边摸一边做爽爽视频免费| 亚洲精品自拍成人| 51国产日韩欧美| 人成视频在线观看免费观看| 大香蕉97超碰在线| 国产成人精品久久久久久| 在现免费观看毛片| 我的老师免费观看完整版| 国产精品嫩草影院av在线观看| 亚洲第一区二区三区不卡| 久久国内精品自在自线图片| 高清午夜精品一区二区三区| 最近手机中文字幕大全| 全区人妻精品视频| 2021少妇久久久久久久久久久| 免费观看在线日韩| 99久久精品国产国产毛片| 日本-黄色视频高清免费观看| 一区二区三区乱码不卡18| 黄色一级大片看看| 热re99久久国产66热| 黑丝袜美女国产一区| 亚洲精品国产av成人精品| 你懂的网址亚洲精品在线观看| 中文字幕av电影在线播放| 看免费成人av毛片| 国产av国产精品国产| 中国美白少妇内射xxxbb| 一级二级三级毛片免费看| 中国美白少妇内射xxxbb| 国产精品女同一区二区软件| 国产伦理片在线播放av一区| 观看美女的网站| 在线精品无人区一区二区三| 国产高清有码在线观看视频| 久久久久久久久久久免费av| 精品久久久精品久久久| 久久99一区二区三区| 午夜福利影视在线免费观看| 亚洲第一区二区三区不卡| 国内精品宾馆在线| 边亲边吃奶的免费视频| 免费黄色在线免费观看| 欧美日韩成人在线一区二区| 秋霞在线观看毛片| 亚洲精品视频女| 国产乱来视频区| 国产熟女午夜一区二区三区 | 久久热精品热| 国产乱来视频区| 少妇被粗大猛烈的视频| 国产深夜福利视频在线观看| 色吧在线观看| 久久狼人影院| 人人澡人人妻人| 老司机影院毛片| 免费日韩欧美在线观看| freevideosex欧美| 18禁动态无遮挡网站| 精品一区二区免费观看| 日韩视频在线欧美| 成人影院久久| 欧美激情极品国产一区二区三区 | 热re99久久国产66热| 五月开心婷婷网| 亚洲国产最新在线播放| 永久网站在线| 国产 精品1| 国产精品成人在线| 晚上一个人看的免费电影| 国产精品秋霞免费鲁丝片| 又黄又爽又刺激的免费视频.| 多毛熟女@视频| 国产伦理片在线播放av一区| 国精品久久久久久国模美| 最近中文字幕2019免费版| 国产精品久久久久成人av| 久久热精品热| 一级黄片播放器| av网站免费在线观看视频| 免费观看的影片在线观看| 国产女主播在线喷水免费视频网站| 久久人人爽人人片av| 日韩中字成人| 91成人精品电影| 久久人人爽人人爽人人片va| 国产亚洲av片在线观看秒播厂| av免费观看日本| 日韩视频在线欧美| 91精品伊人久久大香线蕉| 成人免费观看视频高清| 最近中文字幕高清免费大全6| 亚洲国产日韩一区二区| 一区二区三区免费毛片| 永久网站在线| 久久精品夜色国产| av免费观看日本| 国产亚洲av片在线观看秒播厂| 免费久久久久久久精品成人欧美视频 | 免费大片18禁| 久久精品人人爽人人爽视色| 制服诱惑二区| 97在线人人人人妻| 九色成人免费人妻av| 国产毛片在线视频| 少妇猛男粗大的猛烈进出视频| 精品熟女少妇av免费看|