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

    大型av网站在线播放| 久久午夜综合久久蜜桃| 国产精品熟女久久久久浪| 老司机靠b影院| 欧美午夜高清在线| 一夜夜www| 精品一区二区三区视频在线观看免费 | 国产熟女午夜一区二区三区| www.精华液| 老司机深夜福利视频在线观看| 久久婷婷成人综合色麻豆| 久久毛片免费看一区二区三区| 俄罗斯特黄特色一大片| 国产精品亚洲av一区麻豆| 久久人妻av系列| 国产精品98久久久久久宅男小说| 国产真人三级小视频在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲天堂av无毛| 最近最新中文字幕大全电影3 | 亚洲国产成人一精品久久久| 91精品三级在线观看| 纵有疾风起免费观看全集完整版| 法律面前人人平等表现在哪些方面| 不卡av一区二区三区| 精品人妻1区二区| 美女高潮到喷水免费观看| 成年人黄色毛片网站| 国产高清videossex| 亚洲午夜精品一区,二区,三区| 国产日韩欧美在线精品| av超薄肉色丝袜交足视频| 大码成人一级视频| 美女国产高潮福利片在线看| 少妇被粗大的猛进出69影院| 亚洲五月婷婷丁香| 看免费av毛片| 91字幕亚洲| 国产精品九九99| 欧美日韩福利视频一区二区| 亚洲成人手机| 国产精品麻豆人妻色哟哟久久| 久久中文字幕一级| 美女扒开内裤让男人捅视频| 岛国在线观看网站| 国产不卡av网站在线观看| 中亚洲国语对白在线视频| 法律面前人人平等表现在哪些方面| 亚洲一区中文字幕在线| 91麻豆精品激情在线观看国产 | 51午夜福利影视在线观看| 成人国产av品久久久| 90打野战视频偷拍视频| 成人国产av品久久久| 狠狠婷婷综合久久久久久88av| 日韩一卡2卡3卡4卡2021年| 美女福利国产在线| 色在线成人网| 18禁美女被吸乳视频| 色婷婷av一区二区三区视频| 一边摸一边抽搐一进一出视频| 亚洲国产看品久久| 久久精品成人免费网站| 精品少妇久久久久久888优播| 蜜桃国产av成人99| 视频区欧美日本亚洲| 国产国语露脸激情在线看| 久久国产亚洲av麻豆专区| 久久久久久人人人人人| 老司机影院毛片| 国产精品秋霞免费鲁丝片| 真人做人爱边吃奶动态| 波多野结衣av一区二区av| 免费观看人在逋| 中文字幕色久视频| 五月天丁香电影| 精品视频人人做人人爽| 欧美精品人与动牲交sv欧美| 精品少妇久久久久久888优播| 少妇粗大呻吟视频| 极品教师在线免费播放| 18禁裸乳无遮挡动漫免费视频| 成人永久免费在线观看视频 | 中文字幕最新亚洲高清| 麻豆乱淫一区二区| 日本黄色日本黄色录像| e午夜精品久久久久久久| 国产在视频线精品| 亚洲欧美色中文字幕在线| 满18在线观看网站| 黄色a级毛片大全视频| 婷婷成人精品国产| 国产精品自产拍在线观看55亚洲 | 天天躁日日躁夜夜躁夜夜| videos熟女内射| 免费观看人在逋| 国产亚洲欧美精品永久| 热re99久久精品国产66热6| 99re6热这里在线精品视频| 波多野结衣一区麻豆| 中文字幕av电影在线播放| 午夜精品国产一区二区电影| 亚洲精品中文字幕一二三四区 | 99re6热这里在线精品视频| 国产精品美女特级片免费视频播放器 | 人妻 亚洲 视频| 久久精品成人免费网站| 日日摸夜夜添夜夜添小说| 日韩制服丝袜自拍偷拍| 天堂中文最新版在线下载| 大片电影免费在线观看免费| 亚洲精品在线美女| 精品国产乱码久久久久久男人| 亚洲国产欧美在线一区| 下体分泌物呈黄色| 99在线人妻在线中文字幕 | 大片免费播放器 马上看| 新久久久久国产一级毛片| 天堂中文最新版在线下载| 多毛熟女@视频| 天堂俺去俺来也www色官网| 捣出白浆h1v1| 99久久精品国产亚洲精品| 超碰成人久久| 欧美黄色片欧美黄色片| 国产成人免费观看mmmm| 中文亚洲av片在线观看爽 | 视频区图区小说| 在线 av 中文字幕| 久久久水蜜桃国产精品网| 99久久国产精品久久久| 久久久久网色| 啦啦啦 在线观看视频| 亚洲国产av新网站| 亚洲午夜精品一区,二区,三区| 热99久久久久精品小说推荐| 香蕉国产在线看| 窝窝影院91人妻| av天堂在线播放| 久久人妻熟女aⅴ| 欧美国产精品一级二级三级| 三上悠亚av全集在线观看| 亚洲九九香蕉| av一本久久久久| 999久久久精品免费观看国产| 美女高潮到喷水免费观看| 一本色道久久久久久精品综合| 欧美激情久久久久久爽电影 | 亚洲午夜理论影院| 日本欧美视频一区| 免费观看a级毛片全部| 搡老熟女国产l中国老女人| 欧美在线黄色| 中文字幕人妻熟女乱码| 黄色视频在线播放观看不卡| 国产伦理片在线播放av一区| av又黄又爽大尺度在线免费看| 亚洲五月色婷婷综合| 精品人妻熟女毛片av久久网站| 日本wwww免费看| 老司机深夜福利视频在线观看| 国产不卡av网站在线观看| 免费在线观看黄色视频的| 国产在视频线精品| 亚洲久久久国产精品| 99热网站在线观看| 亚洲精品国产区一区二| 亚洲av日韩在线播放| 成人影院久久| 十八禁高潮呻吟视频| 国产福利在线免费观看视频| 国产av国产精品国产| 无遮挡黄片免费观看| 高清黄色对白视频在线免费看| 久久午夜综合久久蜜桃| 97在线人人人人妻| 免费人妻精品一区二区三区视频| 亚洲专区字幕在线| 十八禁网站网址无遮挡| 熟女少妇亚洲综合色aaa.| 999久久久精品免费观看国产| 国产av精品麻豆| 国产精品久久电影中文字幕 | 免费看a级黄色片| 在线看a的网站| 亚洲精品国产色婷婷电影| 好男人电影高清在线观看| 久久人妻熟女aⅴ| 国产一区二区 视频在线| 亚洲少妇的诱惑av| 久久久久国内视频| 啦啦啦视频在线资源免费观看| 精品亚洲成a人片在线观看| 在线看a的网站| 制服人妻中文乱码| 中文字幕制服av| 久久人妻熟女aⅴ| 久久人人爽av亚洲精品天堂| 日本欧美视频一区| 视频区图区小说| 精品福利观看| 成人亚洲精品一区在线观看| 国产精品一区二区在线不卡| 男女午夜视频在线观看| 亚洲熟女精品中文字幕| 成年版毛片免费区| 俄罗斯特黄特色一大片| 怎么达到女性高潮| 国产淫语在线视频| 日本一区二区免费在线视频| 黄色视频在线播放观看不卡| 91av网站免费观看| 一级毛片精品| 考比视频在线观看| 青草久久国产| 国产老妇伦熟女老妇高清| 一本一本久久a久久精品综合妖精| 国产麻豆69| 欧美精品亚洲一区二区| 视频区欧美日本亚洲| 欧美中文综合在线视频| 午夜免费成人在线视频| 亚洲av日韩在线播放| 看免费av毛片| 国产精品免费一区二区三区在线 | 黄色丝袜av网址大全| 乱人伦中国视频| 1024视频免费在线观看| 一本久久精品| 久久久水蜜桃国产精品网| 国产成人啪精品午夜网站| 大型黄色视频在线免费观看| tocl精华| 精品少妇一区二区三区视频日本电影| 日韩欧美免费精品| 丝袜美足系列| 99在线人妻在线中文字幕 | 黄色成人免费大全| 男女免费视频国产| 一级毛片精品| 天堂动漫精品| 国产精品一区二区在线不卡| 日韩 欧美 亚洲 中文字幕| 欧美日韩福利视频一区二区| 一区二区三区激情视频| 日韩人妻精品一区2区三区| netflix在线观看网站| 久久中文看片网| 久久久久国产一级毛片高清牌| 久久婷婷成人综合色麻豆| 国产精品久久久久久人妻精品电影 | 久久久欧美国产精品| 十分钟在线观看高清视频www| 妹子高潮喷水视频| 满18在线观看网站| 国产亚洲精品一区二区www | 国产麻豆69| 丰满饥渴人妻一区二区三| 国产成人精品久久二区二区91| 日日爽夜夜爽网站| 欧美黄色片欧美黄色片| 国产人伦9x9x在线观看| 人妻 亚洲 视频| 最新在线观看一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 91国产中文字幕| 在线十欧美十亚洲十日本专区| videosex国产| 中文字幕av电影在线播放| 欧美午夜高清在线| 一级毛片电影观看| 日韩 欧美 亚洲 中文字幕| 一进一出抽搐动态| 悠悠久久av| 一级,二级,三级黄色视频| 成人永久免费在线观看视频 | 女性被躁到高潮视频| videosex国产| 国产免费视频播放在线视频| www.999成人在线观看| 免费久久久久久久精品成人欧美视频| 亚洲欧美色中文字幕在线| 99久久国产精品久久久| 午夜精品国产一区二区电影| 中文字幕精品免费在线观看视频| 国产男女超爽视频在线观看| 怎么达到女性高潮| 99国产精品免费福利视频| 国产激情久久老熟女| 久久久国产欧美日韩av| 99久久国产精品久久久| 国产精品久久久久久精品古装| 精品福利永久在线观看| 久久精品亚洲熟妇少妇任你| 免费女性裸体啪啪无遮挡网站| 国产成人免费无遮挡视频| 一进一出好大好爽视频| 久久久久精品国产欧美久久久| 91成人精品电影| 美女主播在线视频| 婷婷成人精品国产| 欧美日韩视频精品一区| 国产无遮挡羞羞视频在线观看| 一级片免费观看大全| 一本大道久久a久久精品| 国产一区二区激情短视频| 黑人欧美特级aaaaaa片| 国产91精品成人一区二区三区 | 久久久水蜜桃国产精品网| 欧美大码av| 日韩 欧美 亚洲 中文字幕| 女人精品久久久久毛片| 精品亚洲成国产av| 国产精品亚洲av一区麻豆| 亚洲精品在线观看二区| 午夜老司机福利片| 老熟女久久久| 三级毛片av免费| 免费在线观看影片大全网站| 亚洲一卡2卡3卡4卡5卡精品中文| 老司机靠b影院| 在线观看免费日韩欧美大片| 在线观看66精品国产| 亚洲国产精品一区二区三区在线| 免费在线观看视频国产中文字幕亚洲| 大香蕉久久网| 日日爽夜夜爽网站| 成人免费观看视频高清| 国产免费视频播放在线视频| 视频区欧美日本亚洲| 午夜免费鲁丝| 黄色视频不卡| 精品国内亚洲2022精品成人 | 亚洲成av片中文字幕在线观看| 国产高清videossex| 久热这里只有精品99| 欧美日韩亚洲国产一区二区在线观看 | 精品少妇一区二区三区视频日本电影| 亚洲av国产av综合av卡| 日韩中文字幕欧美一区二区| 欧美乱妇无乱码| 一区二区三区激情视频| 如日韩欧美国产精品一区二区三区| 亚洲精品成人av观看孕妇| 色精品久久人妻99蜜桃| 久久国产精品影院| 亚洲五月婷婷丁香| 蜜桃国产av成人99| 中文字幕人妻熟女乱码| 欧美+亚洲+日韩+国产| 日韩有码中文字幕| 午夜视频精品福利| 在线观看免费午夜福利视频| 老熟女久久久| 在线av久久热| 日本欧美视频一区| 午夜久久久在线观看| 美女午夜性视频免费| 少妇的丰满在线观看| 新久久久久国产一级毛片| 人人妻人人澡人人看| 成人av一区二区三区在线看| 成人手机av| 黄网站色视频无遮挡免费观看| 一边摸一边抽搐一进一小说 | 国产欧美日韩一区二区三| 久久久精品94久久精品| 日韩欧美一区二区三区在线观看 | 久久 成人 亚洲| 久久久久精品人妻al黑| 国产免费视频播放在线视频| 老司机影院毛片| 亚洲欧美一区二区三区久久| 天堂中文最新版在线下载| 老司机深夜福利视频在线观看| 免费观看av网站的网址| 老熟妇乱子伦视频在线观看| 1024视频免费在线观看| 亚洲熟女毛片儿| 欧美日韩成人在线一区二区| 亚洲国产欧美日韩在线播放| 午夜老司机福利片| 丁香六月天网| 欧美一级毛片孕妇| 久久精品熟女亚洲av麻豆精品| 色综合欧美亚洲国产小说| 丁香欧美五月| 亚洲欧洲日产国产| 亚洲,欧美精品.| 淫妇啪啪啪对白视频| 国产一区二区三区综合在线观看| 国产成人精品久久二区二区免费| 日韩欧美免费精品| 手机成人av网站| 大型av网站在线播放| 国产精品一区二区精品视频观看| 免费人妻精品一区二区三区视频| 女人被躁到高潮嗷嗷叫费观| 午夜老司机福利片| 国产精品麻豆人妻色哟哟久久| 精品午夜福利视频在线观看一区 | 亚洲人成电影免费在线| 丝袜在线中文字幕| 最黄视频免费看| 91字幕亚洲| 性少妇av在线| 国产精品麻豆人妻色哟哟久久| 人妻久久中文字幕网| 一区二区日韩欧美中文字幕| 久久人人爽av亚洲精品天堂| 性高湖久久久久久久久免费观看| 一区二区三区乱码不卡18| 久久中文字幕一级| 老熟妇乱子伦视频在线观看| 久久精品亚洲av国产电影网| 欧美 日韩 精品 国产| 在线观看免费午夜福利视频| 在线观看人妻少妇| 人妻久久中文字幕网| 国产高清激情床上av| 成人影院久久| 欧美激情极品国产一区二区三区| 美国免费a级毛片| 中文字幕另类日韩欧美亚洲嫩草| 成年版毛片免费区| 成人手机av| 久久久久久久国产电影| 国产区一区二久久| 欧美黑人精品巨大| 50天的宝宝边吃奶边哭怎么回事| 人人妻人人澡人人爽人人夜夜| 成人永久免费在线观看视频 | 在线观看www视频免费| 亚洲天堂av无毛| 丝瓜视频免费看黄片| 国产区一区二久久| 精品人妻在线不人妻| 成人影院久久| 国产三级黄色录像| 超碰成人久久| 亚洲人成电影观看| 午夜视频精品福利| 国产精品亚洲av一区麻豆| 一区福利在线观看| 国产一区二区 视频在线| 欧美激情极品国产一区二区三区| 亚洲久久久国产精品| 少妇猛男粗大的猛烈进出视频| 久久中文字幕一级| 老司机亚洲免费影院| 在线看a的网站| 国产精品一区二区精品视频观看| 国产av又大| 黑人巨大精品欧美一区二区蜜桃| 久久久久视频综合| 国产91精品成人一区二区三区 | 国产免费现黄频在线看| 超碰成人久久| 麻豆av在线久日| 黄色 视频免费看| 午夜老司机福利片| 国产又色又爽无遮挡免费看| 淫妇啪啪啪对白视频| 五月开心婷婷网| 色综合婷婷激情| 狠狠狠狠99中文字幕| 在线观看免费视频日本深夜| 在线观看一区二区三区激情| 免费观看av网站的网址| 亚洲专区国产一区二区| 在线观看免费午夜福利视频| 国产高清视频在线播放一区| 久久久久视频综合| 欧美午夜高清在线| 中文字幕高清在线视频| 国产精品久久久久久人妻精品电影 | 亚洲av电影在线进入| 不卡av一区二区三区| 国产成人影院久久av| 蜜桃国产av成人99| 成人国语在线视频| 国产不卡av网站在线观看| 久久午夜综合久久蜜桃| 亚洲av电影在线进入| 女性被躁到高潮视频| 久久精品亚洲av国产电影网| 日韩欧美国产一区二区入口| 午夜久久久在线观看| 老汉色av国产亚洲站长工具| 亚洲午夜理论影院| a级毛片在线看网站| 老汉色∧v一级毛片| 在线av久久热| 国产不卡av网站在线观看| 99国产精品99久久久久| 国产福利在线免费观看视频| 亚洲精品在线美女| 在线av久久热| 1024香蕉在线观看| 欧美日韩国产mv在线观看视频| 电影成人av| 美女视频免费永久观看网站| 日韩大片免费观看网站| 蜜桃在线观看..| 大片免费播放器 马上看| 亚洲七黄色美女视频| 妹子高潮喷水视频| 狂野欧美激情性xxxx| kizo精华| 9191精品国产免费久久| 日韩人妻精品一区2区三区| 三上悠亚av全集在线观看| 狂野欧美激情性xxxx| 人妻一区二区av| 女人爽到高潮嗷嗷叫在线视频| 人成视频在线观看免费观看| 中文字幕制服av| 777米奇影视久久| 91精品三级在线观看| 午夜福利影视在线免费观看| 嫁个100分男人电影在线观看| 亚洲国产成人一精品久久久| 美女扒开内裤让男人捅视频| 国产精品九九99| 在线十欧美十亚洲十日本专区| 制服诱惑二区| 久久久国产欧美日韩av| 91老司机精品| 男女之事视频高清在线观看| 建设人人有责人人尽责人人享有的| 日韩一卡2卡3卡4卡2021年| 欧美黑人精品巨大| 亚洲成人手机| 国产黄色免费在线视频| 韩国精品一区二区三区| 亚洲欧洲日产国产| 五月开心婷婷网| 97人妻天天添夜夜摸| 最黄视频免费看| 日韩 欧美 亚洲 中文字幕| 国产99久久九九免费精品| 亚洲av第一区精品v没综合| 国产一区二区三区视频了| 一本色道久久久久久精品综合| 一本一本久久a久久精品综合妖精| 搡老熟女国产l中国老女人| 黄片大片在线免费观看| 午夜免费成人在线视频| 国产一区二区 视频在线| 亚洲九九香蕉| 国产精品久久久久久精品古装| 18禁裸乳无遮挡动漫免费视频| 男女床上黄色一级片免费看| 免费在线观看视频国产中文字幕亚洲| 亚洲熟女精品中文字幕| 国产又爽黄色视频| 啦啦啦在线免费观看视频4| 久久久久久久精品吃奶| 伊人久久大香线蕉亚洲五| 多毛熟女@视频| 老司机影院毛片| 亚洲精品成人av观看孕妇| 国产区一区二久久| 人妻一区二区av| 色婷婷久久久亚洲欧美| 亚洲人成电影免费在线| 午夜福利乱码中文字幕| 精品国产国语对白av| 久久精品91无色码中文字幕| 精品人妻1区二区| 黄色丝袜av网址大全| 亚洲九九香蕉| 亚洲成人免费电影在线观看| 老司机午夜福利在线观看视频 | 熟女少妇亚洲综合色aaa.| 欧美激情极品国产一区二区三区| 成年版毛片免费区| 菩萨蛮人人尽说江南好唐韦庄| 久久久精品免费免费高清| 欧美乱码精品一区二区三区| 国产黄频视频在线观看| 亚洲成人免费av在线播放| 久久午夜综合久久蜜桃| 51午夜福利影视在线观看| 纯流量卡能插随身wifi吗| 亚洲国产欧美在线一区| 19禁男女啪啪无遮挡网站| 777久久人妻少妇嫩草av网站| 日本vs欧美在线观看视频| 国产不卡av网站在线观看| 中文字幕精品免费在线观看视频| 亚洲熟女精品中文字幕| 一区二区日韩欧美中文字幕| 午夜福利视频精品| 国产黄色免费在线视频| 国产男靠女视频免费网站| 欧美久久黑人一区二区| 国内毛片毛片毛片毛片毛片| 国产成人欧美| 国产成人啪精品午夜网站| 精品国产一区二区三区久久久樱花| 在线观看www视频免费| 美女扒开内裤让男人捅视频| 国产成人精品在线电影| 亚洲精品在线美女| 亚洲国产中文字幕在线视频| 日韩熟女老妇一区二区性免费视频| 在线观看免费视频日本深夜| 亚洲成a人片在线一区二区| 十八禁网站免费在线| 久久久久网色| 99久久精品国产亚洲精品| 欧美激情 高清一区二区三区| 色老头精品视频在线观看| 视频区图区小说|