• <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.

    熟女人妻精品中文字幕| 一本久久中文字幕| 听说在线观看完整版免费高清| 精品乱码久久久久久99久播| 在线免费观看的www视频| 国产黄a三级三级三级人| 哪里可以看免费的av片| 97超视频在线观看视频| 国产精品三级大全| 免费高清视频大片| 高清日韩中文字幕在线| 午夜a级毛片| 国产真实伦视频高清在线观看| www日本黄色视频网| 亚洲乱码一区二区免费版| 国产三级在线视频| 你懂的网址亚洲精品在线观看 | 久久午夜亚洲精品久久| 久久综合国产亚洲精品| 日本a在线网址| 美女内射精品一级片tv| 亚洲欧美清纯卡通| 国产精品久久久久久av不卡| or卡值多少钱| 国产亚洲91精品色在线| 国产私拍福利视频在线观看| 国产伦精品一区二区三区视频9| av国产免费在线观看| 日本 av在线| 日韩制服骚丝袜av| 全区人妻精品视频| 真人做人爱边吃奶动态| 老女人水多毛片| 麻豆国产av国片精品| 亚洲五月天丁香| 亚洲精品色激情综合| 丰满的人妻完整版| 婷婷六月久久综合丁香| 日本精品一区二区三区蜜桃| 哪里可以看免费的av片| 最近视频中文字幕2019在线8| 女人十人毛片免费观看3o分钟| 国产成人a∨麻豆精品| 欧美成人a在线观看| 日韩亚洲欧美综合| 国产真实乱freesex| 99热全是精品| 日本三级黄在线观看| 最近最新中文字幕大全电影3| 欧美高清成人免费视频www| 99久久九九国产精品国产免费| 欧美日韩一区二区视频在线观看视频在线 | 国产黄片美女视频| 国产乱人视频| 少妇被粗大猛烈的视频| 丰满人妻一区二区三区视频av| 中国国产av一级| 亚洲av一区综合| 毛片一级片免费看久久久久| 一进一出抽搐动态| 干丝袜人妻中文字幕| 99热这里只有是精品50| 国产亚洲精品av在线| 日韩国内少妇激情av| 亚洲人成网站在线播放欧美日韩| 此物有八面人人有两片| 色综合亚洲欧美另类图片| 成人特级av手机在线观看| 卡戴珊不雅视频在线播放| 久久天躁狠狠躁夜夜2o2o| 欧美色视频一区免费| 美女高潮的动态| 别揉我奶头 嗯啊视频| 久久久精品大字幕| 中国美白少妇内射xxxbb| 欧美激情久久久久久爽电影| 亚洲中文字幕一区二区三区有码在线看| 亚洲成人中文字幕在线播放| 国产成人91sexporn| a级毛片a级免费在线| 一个人观看的视频www高清免费观看| 麻豆成人午夜福利视频| 精品久久久久久久末码| 午夜爱爱视频在线播放| 国内久久婷婷六月综合欲色啪| 中文字幕人妻熟人妻熟丝袜美| 国产欧美日韩精品一区二区| 三级经典国产精品| 亚洲精品在线观看二区| 久久久成人免费电影| 欧美一区二区国产精品久久精品| 麻豆国产97在线/欧美| 成人午夜高清在线视频| 久久午夜亚洲精品久久| 99热这里只有精品一区| 国产精品电影一区二区三区| 国产午夜精品论理片| 中文字幕av在线有码专区| 国产精品1区2区在线观看.| 中文字幕av成人在线电影| 国产不卡一卡二| 亚洲国产欧美人成| 免费看美女性在线毛片视频| 男插女下体视频免费在线播放| 又爽又黄无遮挡网站| 国产探花在线观看一区二区| 黄色一级大片看看| 欧美三级亚洲精品| 一进一出抽搐gif免费好疼| 午夜a级毛片| 国产欧美日韩精品亚洲av| 丰满乱子伦码专区| 夜夜爽天天搞| 美女 人体艺术 gogo| 内地一区二区视频在线| 国产亚洲精品综合一区在线观看| 欧美不卡视频在线免费观看| 亚洲第一区二区三区不卡| av.在线天堂| 人人妻人人澡欧美一区二区| 黑人高潮一二区| 观看美女的网站| 中国美女看黄片| 久久99热6这里只有精品| 男女那种视频在线观看| 91在线观看av| 一个人看的www免费观看视频| 亚洲精品亚洲一区二区| 日韩中字成人| 亚洲熟妇熟女久久| 免费观看人在逋| 久99久视频精品免费| 久久久久久久久久久丰满| 97超碰精品成人国产| 欧美日韩一区二区视频在线观看视频在线 | 热99在线观看视频| 久久精品综合一区二区三区| 99国产精品一区二区蜜桃av| 亚洲va在线va天堂va国产| 99热全是精品| 国产又黄又爽又无遮挡在线| 一个人看视频在线观看www免费| 不卡一级毛片| 男女下面进入的视频免费午夜| 久久久久久大精品| 毛片女人毛片| 久久久久久久久久久丰满| 成人永久免费在线观看视频| 寂寞人妻少妇视频99o| 搡女人真爽免费视频火全软件 | 桃色一区二区三区在线观看| 国产 一区 欧美 日韩| 国产激情偷乱视频一区二区| 变态另类成人亚洲欧美熟女| 无遮挡黄片免费观看| 国产精品三级大全| 亚洲最大成人手机在线| 亚洲不卡免费看| 日日摸夜夜添夜夜添av毛片| 麻豆久久精品国产亚洲av| 亚洲综合色惰| 91在线精品国自产拍蜜月| 国产麻豆成人av免费视频| 天堂网av新在线| 亚洲最大成人中文| 91午夜精品亚洲一区二区三区| 99久久中文字幕三级久久日本| 夜夜爽天天搞| 成熟少妇高潮喷水视频| 国产伦在线观看视频一区| 国产亚洲精品久久久久久毛片| eeuss影院久久| 国产精品女同一区二区软件| 最近的中文字幕免费完整| 国产大屁股一区二区在线视频| 欧美激情久久久久久爽电影| 午夜影院日韩av| 午夜影院日韩av| 亚洲成人精品中文字幕电影| 日韩一本色道免费dvd| 久久久久久九九精品二区国产| 色哟哟哟哟哟哟| 久久久久国产网址| 网址你懂的国产日韩在线| 国产伦一二天堂av在线观看| 久久久久免费精品人妻一区二区| a级一级毛片免费在线观看| 联通29元200g的流量卡| 精品久久久噜噜| 联通29元200g的流量卡| 在线国产一区二区在线| 欧美成人a在线观看| 一进一出好大好爽视频| 欧美精品国产亚洲| 毛片女人毛片| 午夜精品一区二区三区免费看| 高清午夜精品一区二区三区 | 深夜精品福利| 久久久久国内视频| 欧洲精品卡2卡3卡4卡5卡区| 大香蕉久久网| 亚洲美女黄片视频| a级毛片a级免费在线| 99热6这里只有精品| 国产在线男女| 国产爱豆传媒在线观看| 亚洲av.av天堂| 毛片一级片免费看久久久久| 国产精品久久电影中文字幕| 搞女人的毛片| 久久欧美精品欧美久久欧美| 精品久久久久久久久亚洲| av在线播放精品| 国产精品一二三区在线看| 看片在线看免费视频| 国内少妇人妻偷人精品xxx网站| 久久中文看片网| 淫秽高清视频在线观看| 久久韩国三级中文字幕| 一级黄片播放器| 亚洲在线观看片| 国产精品国产高清国产av| 亚洲精品日韩av片在线观看| 亚洲18禁久久av| 综合色丁香网| 欧美丝袜亚洲另类| 99久久成人亚洲精品观看| 少妇熟女aⅴ在线视频| 欧美+亚洲+日韩+国产| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久午夜电影| 国内精品宾馆在线| 亚洲在线自拍视频| 精品少妇黑人巨大在线播放 | 我要看日韩黄色一级片| 欧美成人一区二区免费高清观看| 国产精品一区www在线观看| 免费看av在线观看网站| 啦啦啦韩国在线观看视频| 日韩精品中文字幕看吧| 日日啪夜夜撸| 亚洲av免费在线观看| 日日干狠狠操夜夜爽| 精品人妻熟女av久视频| 亚洲国产精品国产精品| 看片在线看免费视频| 日韩精品中文字幕看吧| 国产午夜福利久久久久久| 一区二区三区高清视频在线| 少妇高潮的动态图| 亚洲国产高清在线一区二区三| 嫩草影视91久久| 国产白丝娇喘喷水9色精品| 国产白丝娇喘喷水9色精品| 国产伦精品一区二区三区四那| 国产成人freesex在线 | 日日干狠狠操夜夜爽| 久久午夜亚洲精品久久| 最近2019中文字幕mv第一页| 最新在线观看一区二区三区| 色尼玛亚洲综合影院| 国产激情偷乱视频一区二区| 亚洲色图av天堂| 高清午夜精品一区二区三区 | 3wmmmm亚洲av在线观看| 我的女老师完整版在线观看| 日日干狠狠操夜夜爽| 九九爱精品视频在线观看| 免费无遮挡裸体视频| 日韩制服骚丝袜av| 国产精品爽爽va在线观看网站| 美女 人体艺术 gogo| 久久久久免费精品人妻一区二区| 美女免费视频网站| 欧美xxxx黑人xx丫x性爽| 国产又黄又爽又无遮挡在线| 日本成人三级电影网站| 麻豆国产97在线/欧美| a级毛色黄片| 大香蕉久久网| 久久中文看片网| 狂野欧美白嫩少妇大欣赏| 国产精品不卡视频一区二区| 国产一区二区在线观看日韩| 最后的刺客免费高清国语| 可以在线观看毛片的网站| 此物有八面人人有两片| 在线免费十八禁| 国产精品亚洲一级av第二区| 日本爱情动作片www.在线观看 | 亚洲七黄色美女视频| 久久久久久久久中文| 精品乱码久久久久久99久播| 乱码一卡2卡4卡精品| 亚洲精品成人久久久久久| 久久精品夜色国产| 老熟妇仑乱视频hdxx| 男插女下体视频免费在线播放| 欧美激情久久久久久爽电影| 岛国在线免费视频观看| 亚洲欧美日韩东京热| 淫秽高清视频在线观看| 亚洲中文日韩欧美视频| 女生性感内裤真人,穿戴方法视频| 老女人水多毛片| 国产一区二区三区av在线 | 91在线观看av| 亚洲最大成人av| 午夜精品一区二区三区免费看| 久久久成人免费电影| 亚洲欧美精品综合久久99| 人人妻人人澡欧美一区二区| 日韩 亚洲 欧美在线| 日本黄色视频三级网站网址| 国产精品久久久久久亚洲av鲁大| 国产精品美女特级片免费视频播放器| 麻豆久久精品国产亚洲av| 亚洲va在线va天堂va国产| 国产白丝娇喘喷水9色精品| 精品免费久久久久久久清纯| 久久久久免费精品人妻一区二区| 免费av毛片视频| 少妇人妻精品综合一区二区 | 欧美人与善性xxx| 少妇人妻一区二区三区视频| 亚洲精品影视一区二区三区av| 中文字幕久久专区| 村上凉子中文字幕在线| 亚洲高清免费不卡视频| 中文字幕av在线有码专区| avwww免费| 性欧美人与动物交配| 日本五十路高清| 国产精品美女特级片免费视频播放器| 一区二区三区高清视频在线| 亚洲成人久久爱视频| 我的老师免费观看完整版| 日韩高清综合在线| 免费观看在线日韩| 亚洲精品国产av成人精品 | 麻豆成人午夜福利视频| 欧美一区二区精品小视频在线| 在现免费观看毛片| 麻豆久久精品国产亚洲av| 给我免费播放毛片高清在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久久久久精品电影| 久99久视频精品免费| 成人欧美大片| 成人亚洲精品av一区二区| 三级经典国产精品| 欧美一区二区国产精品久久精品| 亚洲精品456在线播放app| 国产老妇女一区| 亚洲最大成人av| 久久久久精品国产欧美久久久| 日本-黄色视频高清免费观看| 国产欧美日韩精品亚洲av| 干丝袜人妻中文字幕| 成人特级黄色片久久久久久久| 成人漫画全彩无遮挡| 国产在视频线在精品| 嫩草影视91久久| 国产精品久久视频播放| 中文字幕av在线有码专区| 久久这里只有精品中国| 国产久久久一区二区三区| 久久精品国产自在天天线| 特大巨黑吊av在线直播| 国产伦精品一区二区三区视频9| 九九爱精品视频在线观看| 久久亚洲精品不卡| 国产aⅴ精品一区二区三区波| 国产精品一及| 国产在线精品亚洲第一网站| 亚洲熟妇中文字幕五十中出| 99久久精品一区二区三区| 哪里可以看免费的av片| 在线国产一区二区在线| 综合色丁香网| 免费无遮挡裸体视频| 六月丁香七月| 日本黄色视频三级网站网址| 成人av在线播放网站| 国产爱豆传媒在线观看| 麻豆精品久久久久久蜜桃| 色在线成人网| 亚洲久久久久久中文字幕| 国产高清激情床上av| 欧美高清性xxxxhd video| 99久久精品热视频| 十八禁国产超污无遮挡网站| av免费在线看不卡| 欧美最黄视频在线播放免费| 色综合站精品国产| 欧美人与善性xxx| 狠狠狠狠99中文字幕| 国内精品一区二区在线观看| 欧美成人一区二区免费高清观看| 国产成人精品久久久久久| 精品一区二区三区视频在线| 波多野结衣高清无吗| 国产视频内射| 3wmmmm亚洲av在线观看| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美三级三区| 国产激情偷乱视频一区二区| 免费看美女性在线毛片视频| 伦理电影大哥的女人| 毛片一级片免费看久久久久| 欧美三级亚洲精品| 国产黄色视频一区二区在线观看 | 搡老熟女国产l中国老女人| 少妇裸体淫交视频免费看高清| 日本熟妇午夜| 亚洲精品久久国产高清桃花| 1000部很黄的大片| 99热这里只有是精品在线观看| 美女大奶头视频| 精品久久久久久成人av| 久久久久久久久久黄片| 成人欧美大片| 少妇裸体淫交视频免费看高清| 欧美激情久久久久久爽电影| 午夜免费激情av| 国产精品一区二区性色av| 亚洲人与动物交配视频| 欧美日韩在线观看h| 中文字幕免费在线视频6| 久久婷婷人人爽人人干人人爱| 国产av一区在线观看免费| 精品乱码久久久久久99久播| 99久久精品一区二区三区| 国产一区二区三区在线臀色熟女| 欧美人与善性xxx| 久久6这里有精品| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美清纯卡通| 国产女主播在线喷水免费视频网站 | 给我免费播放毛片高清在线观看| 网址你懂的国产日韩在线| 久久久久国内视频| 一进一出抽搐gif免费好疼| 在线看三级毛片| 婷婷色综合大香蕉| 精品久久久久久久久久久久久| 日本精品一区二区三区蜜桃| 午夜日韩欧美国产| 亚洲欧美日韩卡通动漫| av.在线天堂| 欧美xxxx性猛交bbbb| 亚洲欧美日韩无卡精品| 在线观看免费视频日本深夜| 熟女人妻精品中文字幕| 国产真实伦视频高清在线观看| 亚洲精品乱码久久久v下载方式| 高清日韩中文字幕在线| 天堂av国产一区二区熟女人妻| 免费看美女性在线毛片视频| .国产精品久久| 好男人在线观看高清免费视频| 久久热精品热| 亚洲欧美精品综合久久99| 91麻豆精品激情在线观看国产| 精品国内亚洲2022精品成人| 长腿黑丝高跟| 免费观看精品视频网站| 淫秽高清视频在线观看| 成年版毛片免费区| 亚洲三级黄色毛片| 我要看日韩黄色一级片| 亚洲精品一卡2卡三卡4卡5卡| 1024手机看黄色片| 亚洲成人精品中文字幕电影| 天堂影院成人在线观看| 九色成人免费人妻av| 啦啦啦观看免费观看视频高清| 日本精品一区二区三区蜜桃| 欧美色欧美亚洲另类二区| 久久久久久久久久久丰满| 国产男靠女视频免费网站| 欧美色视频一区免费| 在线观看免费视频日本深夜| 男女视频在线观看网站免费| 成人综合一区亚洲| 欧美人与善性xxx| 男女边吃奶边做爰视频| av黄色大香蕉| 久久婷婷人人爽人人干人人爱| 校园春色视频在线观看| 亚洲精品成人久久久久久| 中出人妻视频一区二区| 久久国产乱子免费精品| 亚洲一区高清亚洲精品| 亚洲自拍偷在线| 成人鲁丝片一二三区免费| 欧美日本亚洲视频在线播放| 少妇熟女欧美另类| 午夜福利在线在线| 国产视频内射| 久久久午夜欧美精品| 99热全是精品| 最好的美女福利视频网| 国产中年淑女户外野战色| 日产精品乱码卡一卡2卡三| 欧美色欧美亚洲另类二区| 午夜福利18| 久久草成人影院| 国产不卡一卡二| 蜜桃久久精品国产亚洲av| 国产精品久久久久久精品电影| 日韩欧美国产在线观看| 在线免费观看的www视频| 精品久久久久久久久久免费视频| 深夜精品福利| 精品国内亚洲2022精品成人| 床上黄色一级片| 人人妻人人看人人澡| 最近视频中文字幕2019在线8| 国产免费男女视频| 99九九线精品视频在线观看视频| 国产亚洲精品av在线| 看十八女毛片水多多多| 麻豆国产97在线/欧美| 中文字幕久久专区| 国产高清视频在线观看网站| 乱码一卡2卡4卡精品| 亚洲一级一片aⅴ在线观看| 成人美女网站在线观看视频| 欧美精品国产亚洲| 国产av麻豆久久久久久久| 欧美性猛交黑人性爽| 国产精品不卡视频一区二区| 深爱激情五月婷婷| 欧美+日韩+精品| 91在线精品国自产拍蜜月| 看十八女毛片水多多多| videossex国产| 你懂的网址亚洲精品在线观看 | 国产成人freesex在线 | 亚洲国产精品合色在线| 嫩草影院入口| 在线观看免费视频日本深夜| 国产精华一区二区三区| 99久久久亚洲精品蜜臀av| 欧美xxxx性猛交bbbb| 一a级毛片在线观看| 99久久精品一区二区三区| av卡一久久| 久久精品人妻少妇| av黄色大香蕉| 我要搜黄色片| 色哟哟哟哟哟哟| 欧美极品一区二区三区四区| 国国产精品蜜臀av免费| 婷婷精品国产亚洲av| 欧美xxxx性猛交bbbb| 亚洲综合色惰| 欧美激情国产日韩精品一区| 午夜精品在线福利| 国产一区二区在线观看日韩| 国产在线男女| 性色avwww在线观看| 欧美日韩在线观看h| 啦啦啦观看免费观看视频高清| 欧美在线一区亚洲| 中文字幕av成人在线电影| 精品久久久久久久久久久久久| 18禁在线播放成人免费| 色哟哟哟哟哟哟| 日本黄色片子视频| 亚洲美女视频黄频| 亚洲欧美日韩高清在线视频| 国产精品女同一区二区软件| 麻豆久久精品国产亚洲av| 亚洲久久久久久中文字幕| 亚洲最大成人av| 成年女人看的毛片在线观看| 干丝袜人妻中文字幕| 久久精品国产亚洲av香蕉五月| 亚洲在线自拍视频| 美女xxoo啪啪120秒动态图| 亚洲一区二区三区色噜噜| 晚上一个人看的免费电影| 国产伦精品一区二区三区视频9| 精品一区二区三区av网在线观看| 免费大片18禁| 一进一出抽搐gif免费好疼| 亚洲第一电影网av| 我要搜黄色片| ponron亚洲| 精品人妻一区二区三区麻豆 | 午夜免费男女啪啪视频观看 | 国产av在哪里看| 国产高清三级在线| 欧美在线一区亚洲| av天堂在线播放| 天堂网av新在线| 性插视频无遮挡在线免费观看| 人人妻人人澡欧美一区二区| 成人综合一区亚洲| 久久天躁狠狠躁夜夜2o2o| 欧美日韩在线观看h| 日日摸夜夜添夜夜添小说| 69av精品久久久久久| 亚洲第一区二区三区不卡| 一进一出抽搐动态| 欧美bdsm另类| 51国产日韩欧美| 色综合色国产| 亚洲国产精品成人综合色| 天天一区二区日本电影三级| 国产成人影院久久av| 麻豆久久精品国产亚洲av| 91精品国产九色|