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

    Numerical simulation and experimental validation of arc welding ofDMR-249A steel

    2016-04-18 10:14:20RishiPAMNANIVASUDEVANJAYAKUMARVASANTHARAJAGANESH
    Defence Technology 2016年4期

    Rishi PAMNANI*,M.VASUDEVAN,T.JAYAKUMAR,P.VASANTHARAJA,K.C.GANESH

    aPrototype Training Centre,Kalpakkam,Tamil Nadu 603 127,India

    bIGCAR,Tamil Nadu 603 102,India

    cMetallurgical and Materials Engg Dept,NIT,Warangal 506 004,India

    dDept of Mech Engg,NITT,Tiruchirappalli 620 015,India

    Numerical simulation and experimental validation of arc welding of
    DMR-249A steel

    Rishi PAMNANIa,*,M.VASUDEVANb,T.JAYAKUMARc,P.VASANTHARAJAb,K.C.GANESHd

    aPrototype Training Centre,Kalpakkam,Tamil Nadu 603 127,India

    bIGCAR,Tamil Nadu 603 102,India

    cMetallurgical and Materials Engg Dept,NIT,Warangal 506 004,India

    dDept of Mech Engg,NITT,Tiruchirappalli 620 015,India

    The thermo-mechanical attributes of DMR-249A steel weld joints manufactured by shielded metal arc welding (SMAW)and activated gas tungsten arc welding (A-GTAW)processes were studied using Finite Element Model (FEM)simulation.The thermal gradients and residual stresses were analyzed with SYSWELD software using double ellipsoidal heat source distribution model.The numerically estimated temperature distribution was validated with online temperature measurements using thermocouples.The predicted residual stresses prof i le across the weld joints was compared with the values experimentally measured using non-destructive techniques.The measured and predicted thermal cycles and residual stress prof i le was observed to be comparable.The residual stress developed in double sided A-GTAW joint were marginally higher in comparison to f i ve pass SMAW joint due to phase transformation associated with high heat input per weld pass forA-GTAW process.The present investigations suggest the applicability of numerical modeling as an effective approach for predicting the thermo-mechanical properties inf l uenced by welding techniques for DMR-249A steel weld joints.The tensile,impact and micro-hardness tests were carried to compare the welds. Considering benef i ts of high productivity and savings of labor and cost associated withA-GTAW compared to SMAW process,the minor variation in residual stress build up ofA-GTAW joint can be neglected to developA-GTAW as qualif i ed alternative welding technique for DMR-249A steel. ? 2016 China Ordnance Society.Production and hosting by Elsevier B.V.All rights reserved.

    DMR-249A;A-GTAW;SMAW;FEM;SYSWELD;Thermo-mechanical Simulation

    1.Introduction

    DMR-249A is a low carbon HSLA steel.The steel has micro-alloying additions of V,Nb and Ti for grain ref i nement resulting in increased strength and good impact toughness at sub-zero temperatures.The steel has microstructure of f i ne equiaxed ferrite with pearlite content of less than 10%by volume [1,2].The steel is being used in the construction of the hull of various vessels and ships.The conventional arc welding is extensively used during construction and fabrication of ships. The shielded metal arc welding (SMAW)is the most common process used while other arc welding processes also f i nd applications.Though gas tungsten arc welding (GTAW)processproduces high quality weld deposits,the limitations of shallow penetration and low productivity are associated with the process.A variant of GTAW process called activated f l ux gas tungsten arc welding (A-GTAW)has been reported to overcome the limitations of GTAW process.A signif i cant increase in penetration of up to 300%has been reported in with the use of activated f l ux in GTAW welding [3-5].

    The residual stresses (RS)in a component or structure are stresses caused by incompatible internal permanent strains. Welding is one of the signif i cant causes of RS.The measurement of residual stresses developed during welding of a ship’s hull helps in lowering risk to failure by predicting the inf l uence of residual stresses on fatigue,corrosion and other detrimental surface phenomena.Tensile residual stresses are harmful as they assist crack propagation and also contribute to fatigue failure and stress corrosion cracking [6,7].Compressive residual stresses increase wear and corrosion resistance and are benef i cial in preventing origination and propagation of fatigue cracks[8,9].The prediction of residual stress is important to optimize the welding process for reducing their detrimental effects.The estimation of RS prof i le across weld joints has been carried out using X-ray diffraction and ultrasonic techniques [10,11].

    Fig.1.Photos of weld joints (a)SMAW (b)A-GTAW.

    A lot of research has been undertaken to develop numerical models to understand the effect of heat input and welding process on the weld attributes [12-15].P.Duranton et al.[16]explained the simulation methodology for thermo-mechanical analysis of multi pass welding.P.Mollicone et al.[17]explained thermo-elasto-plastic stages of the welding process with f i nite-element(FE)modeling strategies.Srikas et al.[18]carried out research on thermo mechanical simulation of naval structure materials during arc welding.The SYSWELD software has been used to establish FEM analysis as a useful tool in thermo mechanical characterization of welding processes [5,12-18].

    The present work involves the f i nite element modeling of SMAW andA-GTAW processes for DMR-249A steel to predict thermal cycles and residual stresses on 10 mm thick plates and validation by experimental tools.The simulated thermal cycles were validated by thermocouple measurements at specif i ed locations during SMAW multipass welding and autogenous double-sided welding using A-GTAW process.The simulated residual stress prof i les were validated by non-destructive testing techniques.The objective of the research study was to establish the applicability of numerical simulation for residual stresses measurements of DMR-249A steel joints and ascertain the variation in residual stresses prof i le for conventional low heat input SMAW and high productivity A-GTAW technique.

    2.Materials and methods

    The material composition of the plates used in this study is given below in Table 1 and the welding parameters for weld joints are given in Table 2.

    For each weld joint, two plates of dimension 300 × 120 × 10 mm were used and welded to make a weld specimen plate of 300 × 240 × 10 mm with 300 mm weld length.The square butt joint was used for A-GTAW and 70° V-grove butt joint for SMAW (Fig.1).

    The GTAW using activated f l ux was carried out by using automatic welding machine (Fig.2).A DCEN power sourced tungsten electrode (2%thoriated)of 3.2 mm diameter with 60° tip angle and shielding gas (Argon)f l ow rate of 10 l/min was used.The mixture of SiO2,TiO2,NiO,CuO2and combinations of oxide powders specially developed for DMR-249A steel was used as activated f l ux [4].

    The metallographic samples of size 20 × 10 × 10 mm were cut to carry out microscopic studies of the weld joint.Samples were polished from 80 to 2400 grit SiC paper followed by alumina (sizes 5 μm and 1 μm)suspension to obtain mirror f i nish.The specimens were etched using 2%Nital solution,and optical microscopy was carried out to ascertain weld bead prof i le and changes in microstructure of base and weld metal. The microhardness,impact (55 × 10 × 10 mm3)and tensile(gauge length 260 mm,gauge diameter 4 mm)tests of the welded joint were conducted as per ASTM Standards E384,E23 and E8 respectively.

    Temperature measurements during the welding of butt joints were done with help of K-type thermocouples.The K-type thermocouples were spot welded on the plate before welding at a distance of 10 to 20 mm from the center line (Fig.2).For double sidedA-GTAW,both passes being symmetrical,thermocouples were used to measure temperature prof i les for second pass only.For multi pass SMAW,with limited storage capacity for thermocouple data,the time duration between changing of electrodes and removal of slag was not logged.The thermocouple values were used for the validation of thermal cycles obtained from welding model.

    Table 1Chemical composition (wt%)DMR-249A.

    Table 2Welding parameters for weld joints.

    Fig.2.Spot welded thermocouple (a)A-GTAW joint(b)SMAW joint.

    The measurement of residual stresses by Ultrasonic Testing involves selection of weld joint,determination of Acousto-Elastic Constant (AEC)i.e.change in speed or transit time of sound by applying varying loads on tensile specimen,ultrasonic velocity measurements in the weld joint and calculation of residual stresses using AEC.A 2 MHz probe with transmitter angle of 28°(f i rst critical angle for steels)was used to generate longitudinal critically refracted (LCR)waves.The waves have penetration depth of approximately 3 mm.The AEC for DMR-249A steel was calculated to be 0.069 ns/MPa.To measure transit time on the surface of the weld joint using the Ultrasonic probe,grid lines were marked parallel to the weld seam at regular interval of 5 mm starting from 15 mm away from the weld center line at each side of the weld seam.

    The XRD was carried out using Cr-Kα radiation X-ray tube operating with a target current of 7 mA at 30 kV (Model: Rigaku MSF 2M).The ψ angles were tilted in steps of 9°in the range of 0°to 45°.The residual stresses were estimated using the peak shift at ψ angles and d-spacing relationship of (211)plane.Scanning was done in the angular range of 150°to 162° in steps of 0.2°with a dwell time of 3 sec at each step to.The Young’s modulus of the DMR-249A is taken as 210 GPa to estimate the residual stress values.The calculation of residual stresses based on the shift in the peak position of diffracted X-rays of a selected set of planes is discussed in detail in available literature [19,20].

    Numerical model was developed using f i nite element module SYSWELD.The software is designed to carry out thermo-mechanical analysis of welding as sequentially coupled analysis.The sensitivity of mesh was analyzed by trial and error.The mesh was analyzed using variable mesh size with minimum size near the fusion zone and coarse mesh away from the fusion zone corresponding to the temperature gradient in the welded plates.Two different FE models were generated. The model was taken at the XY plane and the welding path parallel to the Y axis.For A-GTAW square butt joint,F(xiàn)E models for 300 × 120 × 10 mm3size plate was generated with 76,440 elements/82,830 nodes.The symmetric model was mirrored.For SMAW 70°V-Groove butt joint FE model for 300 × 240 × 10 mm3plates with included V-Groove at centerline was modeled with 103,350 elements/113,099 nodes. The FE models are shown in Fig.3.

    Fig.3.FE model of weld joint(a)A-GTAW (b)SMAW.

    The simulation consisting of thermo-metallurgical analysis and mechanical analysis requires temperature and phase dependent material properties.For thermal analysis;thermal conductivity,specif i c heat,coeff i cient of thermal expansion and density with respect to temperature is considered.For mechanical analysis;Young’s Modulus,ultimate tensile Strength,Poisson’s ratio with respect to temperature are considered.The available properties of equivalent HSLA steel with comparable chemical composition and mechanical properties were used to carry out simulations [18].

    During welding,most of the heat energy dispersed into the component by conduction mode heat transfer.In this model,the conduction heat transfer based welding simulation was undertaken by decoupling the welding arc from the welding component.The molten pool stirring was suppressed and the problem was considered as conduction heat transfer analysis.The thermal conductivity of the molten pool after melting point was artif i cially doubled to consider the molten pool stirring effect. The heat source model used represents the total volumetric heat f l ux generated by the molten weld pool irrespective of welding process.In FEA model,combined convection and radiation loss was considered as one of the boundary conditions.The model is considered as a “conduction model”,because the localized melting and solidif i cation phenomenon was not directly incorporated.In this FEA model,the liquid phase of weld pool was not formed,whereas the shape of molten pool and appropriate volumetric heat f l ux was represented using double ellipsoidal heat source.The major heat transfer took place in conduction heat transfer mode to the base metal.The convective heat transfer coeff i cient was taken as 25 W/m2.The elastic constraints with 1000 N/mm stiffness were applied for nodes/elements.

    The f i ller metal deposition was considered by activating each pass one by one during the simulation.The model was created with base plate and individual bead as an assembly model.In the solver when the simulation was set up,start and end time for each welding pass was appropriately given to maintain the inter-pass temperature.The remaining weld bead was deactivated when a welding pass was being carried out. Later the next bead was activated when the inter pass temperature reached a certain value.During the experiment the clamp was placed on the top of the plate such that it would not move in the vertical direction whereas the horizontal and transverse direction was allowed to expand.Similarly the model was assigned with the clamp that only restrains the plate in vertical direction.

    3.Theory/calculation

    3.1.XRD equations

    The residual strain present in the atomic planes is measured using XRD technique.The change in inter planer space “d”due to residual stress is measured by shift in peak intensity.The residual stress is calculated by using the relation

    where d0is the strain free inter planner spacing,ε is the calculated strain and angle Ψ is the angle between the surface normal and the strain measurement direction.

    3.2.Numerical simulation theory

    The assumption of welding heat source model for heat generation source in process is essential to simulate a welding process.One of the major characteristics of the heat source is its motion through time and space.Depending on the welding processes and weld bead prof i les,volumetric heat sources of double ellipsoidal [12-14,21],semi-spherical [22],3D conical[23]and combination of heat sources [5]have been widely used for the conduction mode welding simulation.The double ellipsoidal heat source can be used for shallow or high penetration welding process by adjusting its parameters [21].As the present experiments involve high penetration welding (A-GTAW)and shallow penetration/deposition welding (SMAW),a common heat source model(double ellipsoidal)to compare the integrity of both welding processes was used.The “double ellipsoidal”heat source model [21]is shown in Fig.4.

    The volumetric heat f l ux within two different ellipsoid (for any point x,y,z)is described by Eqs. (3)and (4).The semiellipsoids of different heat f l ux are combined to give the heat source.

    where a,b,cfand crare the geometric parameters of ellipsoidal heat source,shown in Fig.4.

    Q=Heat input inWatts,V=Voltage inVolts and I=Current in Amperes.

    Fig.4.Heat source model.

    The parameters rfand rrare proportion coeff i cients representing heat apportionment in front and back of the heat source respectively,where rf+rr=2.It is of great importance to note that the values of Q(x,y,z)given by Eqs. (3)and (4)must beequal at the x=0 plane to satisfy the condition of continuity of the overall volumetric heat source.This leads to def i nition of two constraints rf=2cf/(cf+cr)and rr=2cr/(cf+cr)such that rf/cf=rr/cr.However,the “double ellipsoidal”model is described by arc eff i ciency η,and four geometric parameters a,b,cfand cr.

    4.Results

    Experimentally measured macro bead prof i le and simulated bead prof i le for A-GTAW and SMAW weld joints are shown in Figs.5 and 6 respectively.The heat source was calibrated to achieve the bead prof i le as observed in macro cut section of the welded joints.The heat source parameters were measured from experimentally observed weld attributes and used in the calibration process.Dimensions of heat source were adjusted till it matched with the experimentally observed weld bead.The heat source f i tting parameters used for FEM simulation ofA-GTAW and SMAW joints are given in Table 3.Figs.5 and 6 show close agreement between simulated and experimental weld bead prof i les.

    The experimental validation of temperature distribution was carried out by physically measuring temperature on plate surface using thermocouples placed at 10 mm,15 mm and 20 mm away from the weld bead.Figs.7 and 8 shows the validation of predicted thermal cycles in A-GTAW and SMAW process respectively.

    Table 3Values of heat source calibrated function.

    Figs.7 and 8 indicate higher temperature gradients at areas closer to the welding heat input and lower temperatures away from the weld center.The highest heating and cooling rate is observed at the weld center point.The thermal cycle has a steep rise in the curve gradient and the slope of curve is lesser for cooling after the heat source passes,validating higher heating rate than the cooling rate.

    The heating rate for both experimental and simulated results was observed to be similar.Figs.7 and 8 exhibit that the experimental cooling rate was lesser than the simulated cooling rate.It was also observed that thermocouple readings are lower near to weld centerline as compared to FEM but higher as distance from centerline increases,showing higher cooling rate in FEM.The simulated thermal prof i le showed steeper slopeand higher cooling rate as the heat transfer mode was assumed to be primarily conduction with modif i ed higher thermal conductivity for molten weld pool.

    Fig.5.A-GTAW bead prof i le (a)Macro (b)First pass (c)Second pass.

    Fig.6.SMAW bead prof i les (a)Macro (b)5 passes meshing (c)-(g)Passes 1 to 5.

    Fig.7.A-GTAW Thermal cycle (a)transverse direction,(b)at 10 mm,(c)at 15 mm,(d)at 20 mm.

    The comparison of peak temperature values for each pass measured at the predetermined locations of thermocouple and as deduced from the FEM ofA-GTAW and SMAW processes is given in Table 4.

    The maximum peak temperature measured at molten weld cenetr in FEM for A-GTAW and SMAW depict the melting temperature and was found to be similar as mentioned in Table 4.The difference between maximum temperatures observed at different distances from weld center for A-GTAW and SMAW is attributed to the difference in heat input for the two welding techniques.The higher heat input of 5 kJ/mm for A-GTAW compared with 2 kJ/mm for SMAW process resulted in higher temperatures in the A-GTAW welded plate.Also,thecooling rate decreases at large heat inputs so that larger quantities of Widmanstatten ferrite are obtained with corresponding reduction in the amount of acicular ferrites.Widmanstatten structure is characterized by its low impact values so the A-GTAW weld joint is expected to have lower impact toughness value at sub zero temperatures as compared to SMAW weld joint.

    Table 4Comparison of peak temperatures measured by thermocouple and FEM.

    The residual stress prof i le in the A-GTAW and SMAW process are given in Fig.9(a)-(b).The residual stresses as measured experimentally using XRD,UT (Lcr)techniques and FEM simulation for plates using A-GTAW and SMAW processes are shown in Fig.10.

    The various studies on RS measurements reported existence of maximum tensile RS generally up to yield stress of the base material or even up to values matching the tensile strength.RS of 250 MPa for 2219 Al,300 MPa in 316LN steel,320 MPa in 316L and 750 MPA in Z8CD12 steels have been reported[24-27].In Fig.10(a)-(b),the tensile RS are shown to be present up to 20 mm to 40 mm on either side of the weld line.

    Fig.10 and Table 5 show that experimental stress measurements are comparable with the numerical prediction.The prof i les of the residual stresses measured experimentally and using FEM were observed to be similar.The minor variations in absolute values of stresses is attributable to the difference in residual stress gradients for different welding techniques and inherent volume of inspection of surface and bulk residualstresses of respective technique [5-14,19,20].The XRD measurements are sensitive to surface conditions with assumed depth of penetration in the order of 5 to 30 microns.For the LCRmeasurements,residual stress is the average value over effective penetration of 3 mm.The FEM simulated values are average stress values interpreted from the SYSWELD model.

    Fig.8.SMAW Thermal cycle (a)transverse direction,(b)at 10 mm,(c)at 15 mm,(d)at 20 mm.

    Fig.9.Simulated longitudinal stress using FEM (a)A-GTAW (b)SMAW.

    Fig.10.Residual stress comparison of FEM with XRD and LCR methods (a)A-GTAW (b)SMAW.

    Table 5Comparison of maximum residual stress (MPa).

    Table 6The volume fraction of grain boundary ferrite,Widmanstatten ferrite,polygonal ferrite and microphases using manual point counting.

    The microstructure studies of base metal DMR-249A and arc welded joints were undertaken using optical microscope. The microstructure of base metal and weld metal of different arc welded joints at 500× magnif i cation are given in Fig.11. The microstructure of base metal show predominantly f i ne grained equiaxed ferrite and some percentage of banded structured pearlite.For weld metals,the optical images showed arc welded joints with grain boundary ferrite,Widmanstatten ferrite with aligned second phase along with veins of ferrite,acicular ferrite,polygonal ferrite and microphases.The grain boundary ferrite has equiaxed form or thin veins delineating prior austenite grain boundaries.The sideplate Widmansttten ferrite is seen as the parallel ferrite laths emanating from prior austenite grain boundaries.

    The M-Shape prof i le for measured RS may be attributed to the phase transformation occurring during solid state transformation of equiaxial ferrite in base material to grain boundary ferrite,accicular ferrite,Widmanstatten ferrite,bainite and micro alloying phases,which causes volume changes in weld metal.The volume expansion due to phase transformation counters the volume contraction due to shrinkage,releasing the tensile RS.The RS induced by shrinkage of the molten region are usually tensile.The compressive residual stresses are formed in the transformed areas where the effect of phase transformations is dominant.The area fractions of various ferritic morphologies seen in optical images were ascertained with systematic manual point counting as per ASTM E 562.The SMAW exhibited higher percentage of microphases and polygonal ferrite and A-GTAW associated with high heat input was observed to contain about more of Widmanstatten ferrite. The volume fractions of grain boundary ferrite,Widmanstatten ferrite,polygonal ferrite and microphases for the SMAW and A-GTAW weld metals are given in Table 6.

    The transformation of austenite to these ferritic morphologies are distinguished by the atomic mechanism of transformation as reconstructive (grain boundary/polygonal ferrite)and displacive (Widmanstatten,acicular and bainitic ferrites). While reconstructive transformation is associated with volume change,the displacive transformation is accompanied by invariant plane strain (IPS)[28-31].The M-Shape prof i le for measured RS is attributed to the volume fraction changes due to transformation of equiaxial ferrite in base material to Grain Boundary Ferrite,Accicular Ferrite,Widmanstatten Ferrite,bainite and micro alloying phases in weld metal [32,33].The marginally higher residual stresses inA-GTAW as compared to SMAW is attributed to higher percentage of displacively transformedWidmanstatten ferrite and microphases (combination of acicular ferrite and bainite)in A-GTAW weld metal and also expected formation of higher percentage of bainite in HAZ due to higher heat input and slower cooling rates [34].The presence of higher percentage of reconstructively transformed grain boundary and polygonal ferrite in SMAW as compared to A-GTAW resulted in volume change in weld metal,leading to reduction in tensile residual stresses at weld center.A-GTAW being an automated process,the variations in experimentally measured residual stresses and derived from FEM were observed to be minimal(Fig.11(a)).Minor variation in residual stress values of experimental measurements and FEM for SMAW (Fig.11(b))is due to manual process being used.For SMAW,there is minor variation in the location of passes in experiment and FEM.Also,the energy absorbed during SMAW not being uniform might have impacted the f i eld under stress.

    Fig.11.Optical (500×)images of base metal and weld joints (a)DMR 249A (b)SMAW (c)A-GTAW.

    Fig.12.Stress vs Strain curve for SMAW and A-GTAW cross weld specimen.

    The impact toughness,with notch in weld metal,was carried out at room temperature for SMAW and A-GTAW joints.The impact toughness was found to be 150 J and 200 J for SMAW and A-GTAW joints respectively.The tensile test performed on cross weld specimen of both joints has shown comparable yield strength of about 480 MPa and percentage elongation of 20% but the ultimate tensile strength of 640 MPa in G-TAW specimen was higher as compared to about 600 MPa of SMAW specimen (Fig.12).The tensile fracture of the cross weld specimen of SMAW andA-GTAW occurred in the base metal region away from weld metal which conf i rmed the existence of adequate strength in weld metal of both processes.The comparison of mechanical tests of SMAW and A-GTAW joints are given in Table 7.

    The fractography studies of tensile specimen and V-notch Chapry impact specimen for both the arc weld joints revealed ductile fracture with characteristic cup and cone dimples structure.The fractography studies conf i rmed and corroborated thehigh values measured for UTS and impact tests as given in Table 7 above.The fractured tensile and impact specimens with scanning electron microscope images are shown in Figs.13 and 14 respectively.

    Table 7Comparison of mechanical tests of SMAW and A-GTAW joints.

    The micro hardness values measured across the weld joint at 2 mm and 5 mm below the surface for SMAW and A-GTAW joints is shown in Fig.15.In SMAW joint,hardness was found between 200 and 330 HV0.2whereas for A-GTAW joint hardness was within the scatter band of 200-280 HV0.2.The prof i le of hardness across the weld joint for both joints was observed to be similar.The hardness of weld metal was measured to be higher than the base metal.The micro-hardness values gradually increases from base metal to HAZ and a minor decrease in the hardness value was observed near the fusion zone followed by higher hardness values in weld metal.

    5.Conclusions

    The thermo-mechanical behavior of DMR-249A steel weld joints manufactured by SMAW and A-GTAW processes were studied using Finite Element Model (FEM)simulation.The conclusions based on the present investigations are summarized as follows:

    1)The thermo-mechanical characteristics of A-GTAW and

    Fig.13.Fractured tensile test specimen images (1)SMAW (2)A-GTAW.

    Fig.14.Fractured V-notch Charpy impact specimen images (1)SMAW (2)A-GTAW.

    SMAW for DMR-249A steel were investigated using numerical modeling and experimentally validated using non destructive testing.

    2)The higher heat input of 5 kJ/mm forA-GTAW compared with 2 kJ/mm for SMAW process resulted in higher temperatures in the A-GTAW welded plate.Also,the arcconstriction in A-GTAW inducesconcentrated temperature at the welding center line which attributes higher temperature.The measured and predicted thermal cycles for both the weld joints fabricated by SMAW and A-GTAW processes were observed to be comparable.

    3)The results presented by numerical models were consistent with the practical measurements by XRD and Ultrasonic Testing.The M-Shape prof i le for measured residual stress is attributed to the phase transformation occurring during transformation of equiaxed ferrite in base material to grain boundary ferrite,accicular ferrite,Widmanstatten ferrite,bainite and micro alloying phases which cause invariant plane strain and volume changes in weld metal.

    4)The study establishes reliability of FEM based thermomechanical analysis for interpretation of thermal cycles and residual stress in DMR-249A steel weld joints.

    5)Considering high productivity and savings of cost and labor associated with A-GTAW process,the minor variation of residual stresses buildup in A-GTAW joint can beignored to develop A-GTAW as qualif i ed alternative welding technique for DMR-249A.

    Fig.15.Micro-hardness prof i le across weld joints (a)SMAW (b)A-GTAW.

    Acknowledgement

    The authors express gratitude to the Indian Navy,Mazagaon Dockyard Limited,Mumbai and Cochin Shipyard Limited,Kochi for their constant encouragement and support.Authors would like to thank Dr GK Sharma,Mr S Mahadevan and Mr Arun Kumar for the help extended during residual stress measurements.Authors would also like to express sincere thanks to Dr.BPC Rao for extending the facilities at NDED,IGCAR to conduct the experiments and studies.

    [1]Mallik S,Minz BS,Mishra B.Production of DMR-249A steel at SAIL,bokaro steel plant.Mater Sci Forum 2012;710:149-54.

    [2]Malakondaiah G.Development of speciality low alloy steels,international conference on metals and alloys:past,present and future,December 7-10,2007,pp.17.

    [3]An international patent Tungsten Inert Gas Welding Flux for enhancing Penetration in Single Pass for austenitic stainless steels,Patent no. PCT/IN06/00442/dt.8/11/06.

    [4]Pamnani R,Vasudevan M,Vasantharaja P,Jayakumar T.Optimisation of A-GTAW welding parameters for naval steel (DMR-249A)by design of experiments approach.Proc Inst Mech Eng L J Mater Des Appl 2015;doi:10.1177/1464420715596455.

    [5]ChukkanJR,VasudevanM,MuthukumaranS,RaviKumarR,Chandrasekhar N.Simulation of laser butt welding ofAISI 316L stainless steel sheetusing various heat sources and experimental validation.J Mater Process Technol 2015;219:48-59.

    [6]Mochizuki M.Control of welding residual stress for ensuring integrity against fatigue and stress-corrosion cracking.Nucl Eng Des 2007;237:107-23.

    [7]Dong P,Brus FW.Welding residual stresses and effects on fracture in pressure vessel and piping components:a millennium review and beyond. J Press Vessel Technol 2000;122:329-38.

    [8]Lin Q,Chen J,Chen H.Possibility of inducing compressive residual stresses in welded joints of SS 400 steels.J Mater Sci Technol 2001;17(6):661-3.

    [9]Bate SK,Green D,Buttle D.A review of residual stress distribution in welded joints for defect assessment of offshore structures.HSE Books,AEA Tech PLc;1997.

    [10]Pamnani R,et al.Residual stress studies on arc welding joints of naval steel (DMR-249A).J Manuf Process 2015;doi:10.1016/j.jmapro .2015.09.004.

    [11]Vasantharaja P,Maduarimuthu V,Vasudevan M,Palanichamy P. Assessment of residual stresses and distortion in stainless steel weld joints.Mater Manuf Process 2012;27:1376-81.

    [12]Ravisankar A,Velaga SK,Rajput G,Venugopal S.Inf l uence of welding speed and power on residual stress during gas tungsten arc welding(GTAW)of thin sections with constant heat input:a study using numerical simulation and experimentalvalidation.JManufProcess2014;16(2):200-11.

    [13]Velaga Satish K,et al.Comparison of weld characteristics between longitudinal seam and circumferential butt weld joints of cylindrical components.J Manuf Process 2015;18:1-11.

    [14]Ganesh K,Vasudevan CM,Balasubramanian KR,Chandrasekhar N,Vasantharaja P.Thermo-mechanical analysis of TIG welding of AISI 316LN stainless steel.Mater Manuf Process 2014;29(8):903-9.

    [15]AlghamdiT, Liu S.Low transformation temperature welding consumables for residual stress management:a numerical model for the prediction of phase transformation-induced compressive residual stresses. Weld J 2014;93(12):458s-71s.

    [16]Duranton P,Devaux J,Robin V,Gilles P,Bergheau JM.3D modelling of multi pass welding ofA 316L stainless steel pipe.J Mater ProcessTechnol 2004;153-154:457-63.

    [17]Mollicone P, CamilleriD, Gray TGF, ComlekciT.Simple thermo-elastic-plastic models for welding distortion simulation.J Mater Process Technol 2006;176:77-86.

    [18]Srikas T,Papanikos P,Kermanidis T.Numerical simulation of laser welding process in butt joints joint specimens.J Mater Process Technol 2013;134(1):59-69.

    [19]Cullity BD,Stock SR.Elements of X-ray diffraction.1st ed.New Jersey: Prentice Hall;2001.

    [20]Raj B,Jaykumar T,Mahadevan S,Rai SK.X-ray diffraction based residual stress measurement for assessment of fatigue damage and rejuvenation process for undercarriages of aircraft.J Nondestruct Eval 2009;28(3-4):157-62.

    [21]Goldak JA,Chakravarti A,Bibby M.A f i nite element model for welding heat sources.Met Trans B 1984;15B:299-305.

    [22]Liu C,Zhang JX,Xue CB.Numerical investigation on residual stress distribution and evolution during multipass narrow gap welding of thick-walled stainless steel pipes.Fusion Eng Des 2011;86:288-95.

    [23]Heinze C,Schwenk C,Rethmeier M.Effect of heat source conf i guration on the result quality of numerical calculation of welding-induced distortion.Simul Model Pract Theory 2012;20:112-23.

    [24]Lu H X.,Liu S,Yang JG,Zhang SP,F(xiàn)ang HY.Ultrasonic stress evaluation on welded plates with LCRwave.STWJ 2008;13:70-4.

    [25]Gilles P,El-Ahmar W,Jullien J-F.Robustness analyses of numerical simulation of fusion welding NeT-TG1.Int J Pres Ves Pip 2009;86:3-12.

    [26]Pratihar S,Turski M,Edwards L,Bouchard PJ.Neutron diffraction residual stress measurements inA 316L stainless steel bead-on-plate weld specimen.Int J Pres Ves Pip 2009;86:13-19.

    [27]Belahcene F,Lu J.Determination of residual stress in Z8CDWV12 steel using critically refracted longitudinal waves.JSME Int J A Solid M 2000;43(4):367-73.

    [28]Withers PJ,Bhadeshia HKDH.Residual stress:part 1.J Mater SciTechnol 2001;17:355-65.

    [29]Withers PJ,Bhadeshia HKDH.Residual stress:part II.J Mater Sci Technol 2001;17:366-75.

    [30]DaiH, FrancisJA, StoneHJ, BhadeshiaHKDH,WithersPJ. Characterizing phase transformations and their effects on ferritic weld residual stresses with X-rays and neutrons.Metall Mater Trans A 2008;39A:3070-8.

    [31]Bhadeshia HKDH.Effect of Materials and Processing-ASM Handbook of Residual Stress and Deformation of Steel 2002.

    [32]Rossini NS,Dassisti M,Benyounis KY,Olabi AG.Methods of measuring residual stresses in components.Mater Des 2012;35:572-88.

    [33]Totten GE.Handbook on residual stress,vol.1.ASM International;2002. p.417.

    [34]Bose-Filho WW,Carvalho ALM,Strangwood M.Effect of alloying elements on the microstructure and inclustion formation in HSLA multipass welds.Mat Charac 2007;58:29-39.

    Received 24 October 2015;revised 11 December 2015;accepted 12 January 2016 Available online 2 March 2016

    Peer review under responsibility of China Ordnance Society.

    *Corresponding author.Tel.:+91 9445823932.

    E-mail address:rishi.rashmi.p@gmail.com (R.PAMNANI).

    http://dx.doi.org/10.1016/j.dt.2016.01.012

    2214-9147/? 2016 China Ordnance Society.Production and hosting by Elsevier B.V.All rights reserved.

    3wmmmm亚洲av在线观看| 亚洲人成77777在线视频| 亚洲情色 制服丝袜| 十八禁网站网址无遮挡| 亚洲精品av麻豆狂野| 国产精品不卡视频一区二区| 欧美日韩国产mv在线观看视频| 少妇的逼水好多| 观看美女的网站| 伦精品一区二区三区| 日韩成人av中文字幕在线观看| 精品人妻熟女av久视频| 一区在线观看完整版| 丝袜在线中文字幕| 免费不卡的大黄色大毛片视频在线观看| 日日啪夜夜爽| 一边亲一边摸免费视频| 免费黄网站久久成人精品| 欧美国产精品一级二级三级| 在线观看一区二区三区激情| 又大又黄又爽视频免费| 999精品在线视频| 国产精品秋霞免费鲁丝片| 亚洲精品久久成人aⅴ小说 | 国产精品蜜桃在线观看| 五月玫瑰六月丁香| 国产精品一区二区三区四区免费观看| 久久精品夜色国产| 欧美xxⅹ黑人| 亚洲精品成人av观看孕妇| 日韩av在线免费看完整版不卡| 九九爱精品视频在线观看| 国产成人av激情在线播放 | 成人18禁高潮啪啪吃奶动态图 | 秋霞在线观看毛片| 九九久久精品国产亚洲av麻豆| 免费看不卡的av| 曰老女人黄片| 亚州av有码| 伦精品一区二区三区| 免费黄色在线免费观看| 免费久久久久久久精品成人欧美视频 | 你懂的网址亚洲精品在线观看| 成年女人在线观看亚洲视频| 久久人人爽人人片av| 国产高清有码在线观看视频| 亚洲一级一片aⅴ在线观看| 国产女主播在线喷水免费视频网站| 亚洲av综合色区一区| 色婷婷av一区二区三区视频| 国产精品一区二区三区四区免费观看| 中文字幕亚洲精品专区| 视频中文字幕在线观看| 久久精品国产亚洲av涩爱| 大香蕉97超碰在线| 日韩制服骚丝袜av| 麻豆成人av视频| 最近的中文字幕免费完整| 18禁裸乳无遮挡动漫免费视频| 成人免费观看视频高清| 秋霞在线观看毛片| 有码 亚洲区| 不卡视频在线观看欧美| 亚洲精品第二区| 午夜久久久在线观看| 亚洲国产精品999| 蜜桃国产av成人99| 亚洲精品国产色婷婷电影| 欧美最新免费一区二区三区| 久久久久久久大尺度免费视频| 久久青草综合色| 欧美成人午夜免费资源| 精品一区在线观看国产| 成人亚洲欧美一区二区av| 日本黄色片子视频| 一级片'在线观看视频| 91精品三级在线观看| 五月天丁香电影| 寂寞人妻少妇视频99o| 日韩不卡一区二区三区视频在线| 卡戴珊不雅视频在线播放| 欧美三级亚洲精品| tube8黄色片| 黄色视频在线播放观看不卡| 亚洲精品色激情综合| 亚洲国产精品一区二区三区在线| 成人午夜精彩视频在线观看| 下体分泌物呈黄色| 黑人猛操日本美女一级片| 一级毛片黄色毛片免费观看视频| 老司机亚洲免费影院| 一本色道久久久久久精品综合| 久久鲁丝午夜福利片| 欧美xxxx性猛交bbbb| 亚洲情色 制服丝袜| 最近中文字幕2019免费版| av黄色大香蕉| 能在线免费看毛片的网站| 22中文网久久字幕| 丝瓜视频免费看黄片| 久久精品熟女亚洲av麻豆精品| 王馨瑶露胸无遮挡在线观看| 欧美日韩视频高清一区二区三区二| 午夜激情久久久久久久| 人妻夜夜爽99麻豆av| 伦理电影大哥的女人| 十八禁网站网址无遮挡| 少妇被粗大的猛进出69影院 | 欧美精品国产亚洲| 亚洲国产av新网站| 2022亚洲国产成人精品| 交换朋友夫妻互换小说| 国产成人精品在线电影| 女性生殖器流出的白浆| 午夜福利,免费看| 亚洲三级黄色毛片| 国产深夜福利视频在线观看| 日韩一区二区三区影片| 18禁在线无遮挡免费观看视频| 欧美精品国产亚洲| 亚洲精品乱久久久久久| 一区二区av电影网| 91国产中文字幕| 制服丝袜香蕉在线| 夜夜看夜夜爽夜夜摸| 考比视频在线观看| 日本欧美视频一区| 国产爽快片一区二区三区| 性色av一级| 黄色视频在线播放观看不卡| 欧美日韩av久久| 亚洲色图综合在线观看| 久久ye,这里只有精品| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品色激情综合| 少妇人妻久久综合中文| 少妇被粗大猛烈的视频| 日韩中字成人| 成人手机av| 亚洲精品乱码久久久久久按摩| 九九在线视频观看精品| 日韩精品免费视频一区二区三区 | 久久影院123| 国产成人精品婷婷| 久久99精品国语久久久| 男女边摸边吃奶| 日本免费在线观看一区| 高清午夜精品一区二区三区| 久久女婷五月综合色啪小说| 黑丝袜美女国产一区| 日韩熟女老妇一区二区性免费视频| 国产成人精品福利久久| 亚洲人成网站在线观看播放| 免费黄频网站在线观看国产| 黄色毛片三级朝国网站| 色吧在线观看| 久久免费观看电影| 99视频精品全部免费 在线| 不卡视频在线观看欧美| 国产毛片在线视频| 国产在线视频一区二区| 在线播放无遮挡| 亚洲精品一二三| 黄色怎么调成土黄色| 国产女主播在线喷水免费视频网站| 满18在线观看网站| 久久久久久久久久久免费av| 欧美日本中文国产一区发布| 国产色爽女视频免费观看| 建设人人有责人人尽责人人享有的| 美女国产高潮福利片在线看| 中文乱码字字幕精品一区二区三区| 搡老乐熟女国产| 日韩大片免费观看网站| 国产伦精品一区二区三区视频9| 免费人妻精品一区二区三区视频| 91久久精品电影网| 国产熟女欧美一区二区| 五月开心婷婷网| 一级爰片在线观看| 日本91视频免费播放| 91成人精品电影| av黄色大香蕉| 大又大粗又爽又黄少妇毛片口| 伊人久久国产一区二区| 亚洲美女搞黄在线观看| 乱人伦中国视频| 国产乱人偷精品视频| 精品99又大又爽又粗少妇毛片| 欧美亚洲日本最大视频资源| 国产乱来视频区| 午夜激情福利司机影院| 97精品久久久久久久久久精品| 精品久久久精品久久久| 中文字幕精品免费在线观看视频 | 亚洲伊人久久精品综合| 最新中文字幕久久久久| 国产伦精品一区二区三区视频9| 国产一区二区在线观看av| 午夜激情久久久久久久| 成人国语在线视频| 日韩制服骚丝袜av| 国产片内射在线| 国产亚洲最大av| 久久久久人妻精品一区果冻| 国产片内射在线| av一本久久久久| 国产精品一区二区在线不卡| 青春草亚洲视频在线观看| 免费大片18禁| 亚洲综合精品二区| 免费看av在线观看网站| 国产精品国产三级国产专区5o| 特大巨黑吊av在线直播| 精品亚洲乱码少妇综合久久| 国产色婷婷99| 蜜桃在线观看..| 九九在线视频观看精品| 国产精品国产三级专区第一集| 免费观看a级毛片全部| 免费播放大片免费观看视频在线观看| 久久久久精品久久久久真实原创| 人人妻人人澡人人看| 国产成人精品婷婷| 毛片一级片免费看久久久久| 建设人人有责人人尽责人人享有的| 精品视频人人做人人爽| 日韩av在线免费看完整版不卡| 欧美精品一区二区免费开放| 国产午夜精品久久久久久一区二区三区| 一区在线观看完整版| 亚洲av男天堂| 搡老乐熟女国产| 久久韩国三级中文字幕| 热99国产精品久久久久久7| 少妇熟女欧美另类| 伊人久久国产一区二区| 成人综合一区亚洲| 一级毛片aaaaaa免费看小| 亚洲国产精品999| 精品国产露脸久久av麻豆| 亚洲av欧美aⅴ国产| 26uuu在线亚洲综合色| 久久ye,这里只有精品| av国产久精品久网站免费入址| 成年人午夜在线观看视频| 少妇人妻精品综合一区二区| 亚洲不卡免费看| 免费看不卡的av| 欧美精品一区二区大全| 亚洲国产av影院在线观看| 亚洲色图综合在线观看| 18禁在线无遮挡免费观看视频| 热99国产精品久久久久久7| 国产精品一二三区在线看| 亚洲三级黄色毛片| 成年人免费黄色播放视频| 纵有疾风起免费观看全集完整版| 国产黄频视频在线观看| 女人精品久久久久毛片| 免费观看性生交大片5| 我要看黄色一级片免费的| 黄色毛片三级朝国网站| 麻豆乱淫一区二区| 精品一品国产午夜福利视频| 在线精品无人区一区二区三| 久久久久久久久久久久大奶| 国精品久久久久久国模美| 一本大道久久a久久精品| 狠狠精品人妻久久久久久综合| 成人午夜精彩视频在线观看| 你懂的网址亚洲精品在线观看| 草草在线视频免费看| 久久99热这里只频精品6学生| 九色亚洲精品在线播放| tube8黄色片| 色94色欧美一区二区| 18在线观看网站| av一本久久久久| 最近的中文字幕免费完整| 日韩精品免费视频一区二区三区 | 99久国产av精品国产电影| 欧美日韩亚洲高清精品| av视频免费观看在线观看| 午夜福利网站1000一区二区三区| 女人精品久久久久毛片| 国内精品宾馆在线| 国产精品麻豆人妻色哟哟久久| www.av在线官网国产| 日本色播在线视频| 色哟哟·www| 在线 av 中文字幕| 精品亚洲成国产av| 国产精品99久久久久久久久| 这个男人来自地球电影免费观看 | 国产精品欧美亚洲77777| 国产精品熟女久久久久浪| 伊人久久国产一区二区| 黄色一级大片看看| 视频区图区小说| 高清不卡的av网站| 女性被躁到高潮视频| 在线观看一区二区三区激情| 欧美最新免费一区二区三区| 制服诱惑二区| 青春草亚洲视频在线观看| 51国产日韩欧美| 美女中出高潮动态图| 欧美精品一区二区大全| 午夜久久久在线观看| 亚洲精品日本国产第一区| 久久久久久久久久久丰满| 国产在线免费精品| 如日韩欧美国产精品一区二区三区 | 久久免费观看电影| 国产精品国产三级国产av玫瑰| 亚洲婷婷狠狠爱综合网| 久久精品国产亚洲av天美| 精品酒店卫生间| 69精品国产乱码久久久| 制服人妻中文乱码| 午夜激情福利司机影院| 日韩一区二区视频免费看| 欧美精品人与动牲交sv欧美| 久久久国产一区二区| 日韩欧美精品免费久久| 久久久国产欧美日韩av| 国产av精品麻豆| 桃花免费在线播放| 美女福利国产在线| 蜜桃久久精品国产亚洲av| 日韩成人伦理影院| 亚洲av国产av综合av卡| 丰满迷人的少妇在线观看| 亚洲精品国产av蜜桃| 看十八女毛片水多多多| 少妇的逼好多水| 国产精品久久久久久久久免| 欧美精品人与动牲交sv欧美| 男男h啪啪无遮挡| 热99久久久久精品小说推荐| 99久久中文字幕三级久久日本| 成人二区视频| 亚洲激情五月婷婷啪啪| 婷婷色综合www| 亚洲在久久综合| 在线观看三级黄色| 久久久久国产精品人妻一区二区| 亚洲欧美成人精品一区二区| 国产永久视频网站| 国产高清国产精品国产三级| 久久99热6这里只有精品| 久久久久久久久久成人| 日韩中字成人| 毛片一级片免费看久久久久| 97在线人人人人妻| 成人亚洲精品一区在线观看| 精品久久久久久电影网| 精品久久久精品久久久| 免费黄频网站在线观看国产| 免费观看在线日韩| 亚洲内射少妇av| 黄色视频在线播放观看不卡| 国产精品国产av在线观看| 亚洲av欧美aⅴ国产| 精品一区在线观看国产| 99热这里只有是精品在线观看| 亚洲国产精品999| 亚洲国产av新网站| 夜夜骑夜夜射夜夜干| 久久久久久久亚洲中文字幕| 久久久久久久久久久免费av| 99久久中文字幕三级久久日本| 免费观看的影片在线观看| 亚洲欧美色中文字幕在线| 最近中文字幕2019免费版| 欧美精品国产亚洲| 欧美xxxx性猛交bbbb| 另类亚洲欧美激情| 亚洲欧美一区二区三区黑人 | 成人二区视频| 性高湖久久久久久久久免费观看| 日日撸夜夜添| 国内精品宾馆在线| 秋霞伦理黄片| 街头女战士在线观看网站| 啦啦啦视频在线资源免费观看| 亚洲少妇的诱惑av| 色94色欧美一区二区| 99热网站在线观看| 美女国产高潮福利片在线看| 18+在线观看网站| 日韩不卡一区二区三区视频在线| 国产精品熟女久久久久浪| 丝袜在线中文字幕| 久久精品人人爽人人爽视色| 日日撸夜夜添| 老司机影院成人| 热99久久久久精品小说推荐| 日本午夜av视频| 99九九在线精品视频| 亚洲丝袜综合中文字幕| 亚洲精品,欧美精品| 久久国产亚洲av麻豆专区| 国产在线一区二区三区精| 欧美日韩一区二区视频在线观看视频在线| 免费久久久久久久精品成人欧美视频 | 亚洲人成77777在线视频| 一边亲一边摸免费视频| 欧美丝袜亚洲另类| 久久女婷五月综合色啪小说| 男女边吃奶边做爰视频| 亚洲国产成人一精品久久久| 亚洲欧美一区二区三区黑人 | 99re6热这里在线精品视频| 国产伦精品一区二区三区视频9| 国产精品嫩草影院av在线观看| 亚洲人成网站在线观看播放| 观看av在线不卡| 久久国内精品自在自线图片| 午夜免费观看性视频| 国产国语露脸激情在线看| 亚洲国产色片| 人妻制服诱惑在线中文字幕| 大片电影免费在线观看免费| 国产精品久久久久久精品古装| 欧美日韩综合久久久久久| 亚洲精品国产色婷婷电影| 亚洲不卡免费看| a级毛片在线看网站| 亚洲,一卡二卡三卡| 日韩精品免费视频一区二区三区 | 中文欧美无线码| 国产免费视频播放在线视频| 欧美人与善性xxx| 亚洲国产欧美在线一区| 99热这里只有是精品在线观看| 日本黄色日本黄色录像| 日韩三级伦理在线观看| 性色av一级| 色网站视频免费| 成人影院久久| 亚洲精品一区蜜桃| 亚洲一区二区三区欧美精品| av国产精品久久久久影院| 国产高清三级在线| 街头女战士在线观看网站| 久久99一区二区三区| 精品一区二区免费观看| a 毛片基地| 国产成人精品福利久久| 亚洲一级一片aⅴ在线观看| 婷婷色av中文字幕| 日本黄大片高清| 最近中文字幕高清免费大全6| 日本欧美视频一区| 一区在线观看完整版| 国产精品熟女久久久久浪| 99国产精品免费福利视频| 日韩中文字幕视频在线看片| 午夜精品国产一区二区电影| 国产成人精品无人区| 亚洲精品av麻豆狂野| 午夜免费观看性视频| 国产精品秋霞免费鲁丝片| 在线免费观看不下载黄p国产| 国产男女超爽视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 汤姆久久久久久久影院中文字幕| 亚洲成人手机| 91精品一卡2卡3卡4卡| 黑人高潮一二区| 午夜精品国产一区二区电影| 久久热精品热| 久久午夜福利片| 熟女av电影| 女人精品久久久久毛片| 麻豆精品久久久久久蜜桃| 国产精品久久久久久久电影| 美女xxoo啪啪120秒动态图| 午夜影院在线不卡| 日韩伦理黄色片| 少妇人妻精品综合一区二区| 国产高清不卡午夜福利| 搡老乐熟女国产| 免费少妇av软件| 国产欧美日韩综合在线一区二区| 亚洲国产精品成人久久小说| 亚洲精品美女久久av网站| 九色成人免费人妻av| 精品亚洲乱码少妇综合久久| 超色免费av| 最近的中文字幕免费完整| 在线天堂最新版资源| 最近中文字幕高清免费大全6| 欧美一级a爱片免费观看看| 美女主播在线视频| 亚洲欧美一区二区三区国产| 欧美日韩综合久久久久久| 国产精品女同一区二区软件| 亚洲av成人精品一区久久| 久久精品国产自在天天线| 久久精品夜色国产| 天堂中文最新版在线下载| 亚洲精品久久午夜乱码| 少妇被粗大的猛进出69影院 | 99精国产麻豆久久婷婷| 色视频在线一区二区三区| 久久这里有精品视频免费| 十分钟在线观看高清视频www| 中文字幕制服av| 国产精品成人在线| 精品人妻熟女av久视频| 日本免费在线观看一区| 免费观看无遮挡的男女| 十八禁网站网址无遮挡| 欧美另类一区| 久久久欧美国产精品| 26uuu在线亚洲综合色| 亚洲精品第二区| 亚洲第一区二区三区不卡| av国产精品久久久久影院| 熟女av电影| 国产国拍精品亚洲av在线观看| 精品亚洲成a人片在线观看| 少妇被粗大猛烈的视频| 天堂俺去俺来也www色官网| 美女大奶头黄色视频| 中文字幕精品免费在线观看视频 | 91精品三级在线观看| 精品久久久久久久久av| 久久久a久久爽久久v久久| 亚洲欧美成人综合另类久久久| 亚洲成人一二三区av| av播播在线观看一区| 亚洲精品第二区| 成人毛片60女人毛片免费| 欧美日韩国产mv在线观看视频| 少妇高潮的动态图| 日本爱情动作片www.在线观看| 亚洲精品一二三| 91久久精品国产一区二区三区| 亚洲av成人精品一区久久| 日韩大片免费观看网站| 老熟女久久久| 少妇被粗大猛烈的视频| 韩国高清视频一区二区三区| 久久热精品热| 成年人午夜在线观看视频| 99视频精品全部免费 在线| 大片免费播放器 马上看| videosex国产| 精品少妇内射三级| 免费大片黄手机在线观看| 国产高清国产精品国产三级| 国产av一区二区精品久久| 久久人人爽人人爽人人片va| 欧美精品高潮呻吟av久久| 亚州av有码| 一个人免费看片子| 最近最新中文字幕免费大全7| 亚洲国产精品国产精品| 一区二区av电影网| 狠狠精品人妻久久久久久综合| 国产 一区精品| 成人二区视频| 日本午夜av视频| 免费看av在线观看网站| av播播在线观看一区| 国产视频内射| 亚洲精品,欧美精品| 国产亚洲一区二区精品| 在线观看免费日韩欧美大片 | 国产成人精品久久久久久| 一本大道久久a久久精品| 天堂中文最新版在线下载| 免费观看的影片在线观看| 免费观看av网站的网址| 大香蕉97超碰在线| 在线看a的网站| av女优亚洲男人天堂| 亚洲精品aⅴ在线观看| 人妻少妇偷人精品九色| 日本黄大片高清| 在线 av 中文字幕| 日韩av在线免费看完整版不卡| 秋霞在线观看毛片| 大香蕉久久成人网| 亚洲av综合色区一区| 制服诱惑二区| 香蕉精品网在线| 亚洲欧美成人综合另类久久久| 九九爱精品视频在线观看| 精品国产一区二区久久| 亚洲少妇的诱惑av| 成人国语在线视频| 黄色欧美视频在线观看| 亚洲欧洲日产国产| 另类亚洲欧美激情| 18禁在线播放成人免费| 91精品三级在线观看| 黄色毛片三级朝国网站| .国产精品久久| 一本久久精品| 在线看a的网站| 51国产日韩欧美| 欧美日韩成人在线一区二区| 亚洲av综合色区一区| 久久精品国产自在天天线| 国产在线一区二区三区精| 亚洲av男天堂| 一级爰片在线观看| 国产亚洲最大av| 欧美激情国产日韩精品一区| 日本爱情动作片www.在线观看|