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

    Multiobjective optimization of friction welding of UNS S32205 duplex stainless steel

    2015-11-01 07:13:57AJITHBirendrKumrBARIKSATHIYAARAVINDAN
    Defence Technology 2015年2期
    關(guān)鍵詞:念書伯父熟人

    P.M.AJITH,Birendr Kumr BARIK,P.SATHIYA,*,S.ARAVINDAN

    aDepartment of Production Engineering,National Institute of Technology,Tiruchirappalli 620015,Tamilnadu,India

    bDepartment of Mechanical Engineering,Indian Institute of Technology Delhi,New Delhi 110016,India

    ?

    Multiobjective optimization of friction welding of UNS S32205 duplex stainless steel

    P.M.AJITHa,Birendra Kumar BARIKa,P.SATHIYAa,*,S.ARAVINDANb

    aDepartment of Production Engineering,National Institute of Technology,Tiruchirappalli 620015,Tamilnadu,India

    bDepartment of Mechanical Engineering,Indian Institute of Technology Delhi,New Delhi 110016,India

    The present study is to optimize the process parameters for friction welding of duplex stainless steel(DSS UNS S32205).Experiments were conducted according to central composite design.Process variables,as inputs of the neural network,included friction pressure,upsetting pressure,speed and burn-off length.Tensile strength and microhardness were selected as the outputs of the neural networks.The weld metals had higher hardness and tensile strength than the base material due to grain refinement which caused failures away from the joint interface during tensile testing.Due to shorter heating time,no secondary phase intermetallic precipitation was observed in the weld joint.A multi-layer perceptron neural network was established for modeling purpose.Five various training algorithms,belonging to three classes,namely gradient descent,genetic algorithm and Levenberg-Marquardt,were used to train artificial neural network.The optimization was carried out by using particle swarm optimization method.Confirmation test was carried out by setting the optimized parameters.In conformation test,maximum tensile strength and maximum hardness obtained are 822 MPa and 322 Hv,respectively.The metallurgical investigations revealed that base metal,partially deformed zone and weld zone maintain austenite/ferrite proportion of 50:50.

    Artificial neural network;Duplex stainless steel;Hardness;Tensile test;Friction welding;Particle swarm optimization

    1.Introduction

    Duplex stainless steel has equal phase balance of approximately equal amounts of ferrite and austenite.It has a mixed microstructure consisting of ferrite(bcc)and austenite(fcc)phases.The duplex stainless steel exhibits higher resistance to stress corrosion cracking and has higher strength than austenitic stainless steel.As a result of these positive factors,the duplex stainless steel is widely used in the oil and gas,petrochemical,pulp and paper,and pollution control industries.It is well known that the duplex stainless steel exhibits good weldability,but the melting and solidification associated with fusion welding processes destroy the favorableduplex microstructure of this stainless steel[1,2].A major concern with fusion welding of duplex stainless steel is the formation of detrimental intermetallic phases at elevated temperatures.Sigma and chi phases form in duplex stainless steels at elevated temperature and precipitate preferably in the ferrite.This will considerably affect the toughness of the welded joint[3].The formations of these phases are due to the high chromium and molybdenum contents.The problem mentioned above can be overcome by employing solid state welding process like friction welding.

    Friction welding is a solid state welding process.It makes use of frictional heat generated on the rubbing surfaces to raise the temperature at the interface,which is high enough to cause the two surfaces to be forged together at high pressure.Friction welding has significant economic and technical advantages.The present study utilized a continuous drive friction welding machine.In continuous drive friction welding,oneworkpiece is rotated at nominal constant speed and aligned with the second part at an applied pressure.The rotation and pressure are maintained for a specific period to ensure adequate thermal and mechanical conditioning of the interface region.Thereafter the rotation is stopped by forced braking,and at the same time the pressure is increased to forge the parts together.The application of an axial force maintains an intimate contact between the parts and causes the plastic deformation of the material near the weld interface.

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

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

    The friction welding finds widespread industrial use as a mass production process for the joining of materials.Friction welding process allows welding of several materials that are extremely difficult to fusion weld.The friction welding process parameters play a significant role in making good quality joints[4].To produce a good quality joint,it is important to set up proper welding process parameters.Therefore,identifying the suitable combinations of process input parameters to produce the desired output requires many experiments,making this process time-consuming and costly[5].

    So as to avoid this problem,various optimization methods can be applied to define the desired output variables by developing the mathematical models to specify the relationship among the input parameters and output variables. Generally,the quality of a weld joint is directly influenced by the welding input parameters during the welding process. Therefore,welding can be considered as a multi-input multioutput process.Though several studies have been made on weld quality by considering response variables separately,the report on simultaneous consideration of response variables in friction welding of DSS material is scarce.

    Artificial neural network(ANN)is a modelling technique which is inspired by the structure and functional aspects of biological neural networks.It has been widely used for modelling manufacturing related problems.ANN eliminates the limitations of the classical approaches by extracting the desired information using the input data.Applying ANN to a system needs sufficient input and output data instead of a mathematical equation[6,7].ANN is a good alternative to conventional empirical modeling based on polynomial and linear regressions[8].Employing the neural networks models would lead to time and cost saving by predicting the experimental results[9].Sathiya et al.[10]used the evolutionary computational techniques optimize the friction welding parameters.The suggested methods were used to determine the welding process parameters by which the desired tensile strength and minimized metal loss were obtained in friction welding.They described how to obtain near optimal welding conditions over a wide search space by conducting a smaller number of experiments.Paventhan et al.[11]have done the optimization of friction welding process parameters for joining carbon steel and stainless steel.They developed an empirical relationship to predict the tensile strengths of friction welded AISI 1040 grade medium carbon steel and AISI 304 austenitic stainless steel,incorporating the process parameters,namely friction force,forging force,friction time and forging time which have greater influence on strength of the joint.Response surface methodology was applied to optimize the friction weldingprocessparameterstoattainmaximumtensile strength of the joint.Koen et al.[12]developed a new welding method for fully automatic pipelines girth welding using a new friction welding machine.The proposed new welding procedure,called Friex,is a new variant of the well-known friction welding process.An intermediate ring is rotated between the pipes to be welded to generate the heat necessary to realize the weld.Luo et al.[13]designed a mixed-integrated approach to control the welding flashes in the continuous drive friction welding on small diameter tubes and then eliminate the problems of the inner friction welding flashes.Experimental results show that this mixed-integrated approach is used to reduce the sum of the inner welding flashes to control the crimping direction of the friction welding flashes,which promotes the forming of the outer friction welding flashes. Udayakumar et al.[14],carried out the experimental investigation on super duplex stainless steel and attempted to develop the mathematical models based on response surface methodology in order to predict the corrosion current and impact strength as a function of key input parameters in the friction welding process.They found that the friction force is a significant parameter for changing the impact strength.Friction force and burn-off length have a negative effect on impact strength.As the friction force increases,the impact strength decreases.

    In a very recent investigation,an attempt has been made to optimize the process parameters of activated tungsten inert gas(ATIG)welding process for ASTM/UNS S32205 DSS joints to obtain desirable aspect ratio and average ferrite number(FN). This investigation revealed that,the joints were produced by the optimized process parameters and their average ferrite number(FN)in the weld zone is 71.62,and the ferrite content is approximately 50.674%which is well within the acceptable range[15].

    From the above literatures,a very few investigations were carried out on friction welding and the parameter optimization of duplex stainless steel using evolutionary algorithms like artificial neural networks(multiobjective optimization).The aim of this work is to predict and optimize the friction welding parameters of UNS S32205 duplex stainless steel using ANN and PSO.The input(factors)parameters were friction pressure(FP),upsetting pressure(UP),rotational speed(N)and burnoff length(BOL),and the responses were tensile strength(TS)and hardness.The five training algorithms selected were batch back propagation(BBP),incremental back propagation(IBP),quick propagation(QP),Levenberg-Marquardt(LM),and genetic algorithm(GA).Among the five algorithms,the QP algorithm had a better performance.The QP algorithm was applied to the ANN network for modeling of friction welding parameters and also to study the direct effect of the individual parameters.

    2.Material and methods

    The base material chemical compositions of UNS S32205 duplex stainless steel were analyzed using an optical emission spectrometer,and their values are presented in Table 1.

    Table 1 Base material chemical composition(wt.%).

    The friction welding parameters were chosen based on the preliminary welding trials and their parameter levels were categorized as low,medium and high.The upper and lower limits were identified and the different levels of process parameters are presented in Table 2.

    The friction welding trials were conducted as per the central composite rotatable factorial design consisting of 30 sets of coded conditions[16].The microstructure samples were prepared as per the standard procedure.The samples were prepared by electrolytically etching them in 10%oxalic acid at 9 V for 30 s as per ASTM E3-11.The mechanical characteristics of friction welds were evaluated from tensile tests as per the ASTM E 8 standards.The tensile-tested samples are presented in Fig.1.Microhardness survey was carried out using a HMV-2000 Vickers microhardness tester at 500 g load for 10 s.

    3.Methodology

    Typical ANN model flow chart was used in this study,as shown in Fig.2.The basic steps considered for designing the neural network model are to collect the data required for training the network,designing the network architecture and training the network.

    Alldatasetswereobtainedfromtheexperiments mentioned above and were divided into three parts,i.e.21,5 and 4 data as training data,testing data and validation data,respectively,which are summarized in Table 4.The training data were used to compute the network parameters.The Neural Power professional version 2.5 software tool was employed in this study.The models were developed to establish the input-output correlations of the friction welding of duplex stainless steel(DSS UNS S32205)using the neural networks.ANNs were constructed with layers of units,and thus termed multilayer ANNs.Multilayer perceptron(MLP)neural network consisted of input,hidden and output units.A single hidden layer network was used in this study.The optimum number of neurons in hidden layer was determined by a series of topologies,in which the number of neurons was varied from 1 to 30.The root mean square error(RMSE)was used as the error function.

    Table 2 Upper and lower limits with different levels of the parameters.

    Fig.1.Tensile-tested weld samples.

    Experimental data were used to train the network.Scaled data were passed into the input layer and then were propagated from input layer to hidden layer and finally to the output layer of the network.Every node in hidden or output layer first acted as a summing junction which combined and modified the inputs from the previous layer using the following equation.

    where yiis the net input to node j in hidden or output layer;xiis the input to node j(or output of previous layer);wijis the weight representing the strength of the connection between the ith node and jth node;i is the number of nodes;and bjis the bias associated with node j.

    In order to perform a supervised training,ANN output error between the actual and predicted output results could be evaluated.The following equations were used to find R,DC,AAPD and RSME.DC reflects the degree of fit for themathematical model.The DC shows the level of model fitness. If value of DC is closer to 1,the model is considered as a better design and fits to the actual data.

    Fig.2.Typical ANN model flow chart.

    Table 3 ANN algorithm setting for five learning algorithms.

    where E=experimental value P=predicted value;E- is the mean value of Eiis the mean value of Pi.Coefficient of determination is

    Average absolute percentage deviation(AAPD)is

    So,we considered the ANN model with lowest RMSE,AAPD and highest DC,R as the best ANN design.

    3.1.Optimization procedure

    Particle swarm optimization(PSO)was used to optimize the friction welding parameters and its flow chart for PSO is presented in Fig.3.The trained data obtained from ANN should be imported into PSO as the basic input data.

    4.Results

    4.1.Selection of the best neural network model

    The results for various algorithms are summarized and presented in Table 3.As shown in Table 3,the QP algorithm has a better performance relative to BBP,IBP,LM and GA algorithms.Fig.4 shows the plots of ANN predicted response versus actual response with QP algorithms for the training and testing data.From Fig.4 it is seen that the calculated and observed hardness's and tensile strengths have almost the same values for QP algorithms.

    The percentages of contribution to the individual input parameters are presented in Fig.5.

    Fig.3.Flow chart of PSO.

    Fig.4.Training and testing data with experimental results for TS and mH using QP model.

    From Fig.5,it can be seen that the upsetting pressure is the most significant parameter followed by the friction pressure, speed of rotation and born-off length.The percentage contributions of individual parameter are upsetting pressure of 47.05%,friction pressure of 25.99%,speed of rotation of 21.96%and burn-off length of 4.96%.

    Fig.5.Percentages of contribution to the input parameters using QP model.

    4.2.Effects of process parameters

    The effects of upsetting pressure and friction pressure on tensile strength and microhardness are presented in Fig.6.

    From Fig.6(a)it is seen that the friction pressure varied from 45 MPa to 125 MPa and the upsetting pressure varied from 155 MPa to 185 MPa,and,the hardness and tensile strength increased at maximum friction pressure Maximum tensile strength obtained was 825 MPa at friction pressure of 185 MPa and upsetting pressure of 105 MPa.In Fig.6(b),the obtained maximum microhardness was 323 MPa at friction pressure of 85 MPa and upsetting pressure of 105 MPa.The predicted model using quick propagation algorithm was fitted so well to the actual values for both training and testing set. Therefore,it could be suggested that the model trained with QP algorithm is the most efficient model for this problem;hence this model was applied for further application.It was reported that the quick propagation learning algorithm could be adopted for the training of all the ANN models[17].Thepredicted values of the best model for training and testing set are presented in Table 4.

    Fig.6.Effects of friction pressure and upsetting pressure on TS and Microhardness(mH).

    Table 5 Results of confirmation test.

    From Table 4,it is clearly seen that the testing and training data are very closer to the experimental results.So QP model is the most appropriate method to predict the effect of the individual parameters of friction welding process.The input parameters were obtained from the QP model and fed in to the PSO algorithm to get the optimized values.The optimized values of friction pressure,upsetting pressure,rotational speed and burn-off length are 105 MPa,180 MPa,2000 rpm and 3.9 mm,respectively.The tensile strength of 827.17 MPa and the hardness of 325.61 Hv were obtained.

    4.3.Results of confirmation test

    因?yàn)槌潜緛?lái)是不大的,有許多熟人,也都是來(lái)看燈的都遇到了。其中我們本城里的在哈爾濱念書的幾個(gè)男學(xué)生,他們也來(lái)看燈了。哥哥都認(rèn)識(shí)他們。我也認(rèn)識(shí)他們,因?yàn)檫@時(shí)候我們到哈爾濱念書去了。所以一遇到了我們,他們就和我們?cè)谝黄?,他們出去看燈,看了一?huì),又回到我們的地方,和伯父談話,和哥哥談話。我曉得他們,因?yàn)槲覀兗冶容^有勢(shì)力,他們是很愿和我們講話的。

    In order to further validate the obtained results,the confirmation test was carried out to verify the PSO model results.The results of confirmation test are shown in Table 5.

    From Table 5,it is seen that the results obtained from test is very close to the results of PSO.The photograph and macrograph of optimized friction welded sample are presented in Fig.7(a)and(b).

    The different zones of microstructures for optimized parameters are presented in Fig.8.

    The joint interface is called as the weld zone(WZ)and the adjunct side is called as the partially deformed zone(PDZ)and also unaffected base metal(BM).The percentage of ferrite phase was measured using Fischer Feritscope MP 30 and average ferrite values are presented in Table 6.

    The microhardness test was carried out along the longitudinal direction of the joint interface with load of 500 g and dwell time of 10 s.The microhardness profile is shown in Fig.9.

    Fig.10 presents the SEM fracture surface of the confirmation tensile tested sample.The fracture of tensile weld sample takes place away from the joint interface.

    The XRD pattern of optimized parameter weld zone is presented in Fig.11.

    5.Discussions

    Fig.7.Confirmation test.

    Fig.8.SEM micrograph of different zones of friction weld joint..

    Table 6 Austenite-to-ferrite values in different weld zones.

    In the real time engineering applications of manufacturing industries,the selection of proper process parameters plays a crucial role in making quality products.In this study,the friction welding process parameters were optimized using PSO.The best result was obtained from QP algorithm with 4-4-2 topology that had minimum RMSE and AAP,maximum R and DC for both training and testing sets.From Fig.4,it is seen that a good agreement was made on predicted values with observed values of QP algorithms.Fig.5 shows that the joint performance of the friction welded joint mainly depends on the upsetting and friction pressures.Because of more refined grain formation during the high UP,a large amount of heat is generated.Fig.6(a)and(b)revels that more heat was generated at higher friction pressure and also more soften state was obtained.It was also observed that,when the upsetting pressure was high,less amount of soften material was retained in the weld zone.The grains were finer due to the faster cooling rate in the weld and HAZ regions.Coarse grains were observed at low upsetting pressure because of slow rate of cooling.The weld zone consisted of finer grains,while the adjunct side PDZ had the coarse grains.The refinement of grains in the weld zone was due to dynamic recrystallization and also due to higher temperature in the peripheral region[18].

    Fig.9.Microhardness profile for optimized parameter.

    During the friction welding,heat generated at both ends of the metal rod is high.Due to this,the ends of the rods are in red hot or plastic stage.Later the upsetting pressure is applied and welding is completed with immediate atmospheric air cooling.

    Fig.10.SEM fracture surface of tensile tested sample.

    Fig.11.XRD pattern of optimized parameter weld zone.

    So the simultaneous effect of short time cooling and high upsetting pressure will result in the formation of fine grains in PDZ and WZ regions.Another advantage is of the short timecooling prevents the phase changes.Therefore optimum amount of austenite was reformed.Table 6 reveals that the austenite-to-ferrite ratio is 50:50 in the three zones of the weld like BM,PDZ and WZ.It was found that hardness and tensile strength was mostly influenced by an interactive effect of upsetting pressure and heating pressure.

    From the confirmation experiment it is observed that PSO resultsareingoodagreementwiththeexperimentalresults.The deviationoftheexperimentalvalueswaslessthan2%compared with optimized PSO results of tensile strength and microhardnessvalues.Fig.8(a)and(b)revealthattheweldingcarriedwith the optimized values resulted in nowelding defects like cavities and cracks in the interface.Presence of secondary phases in duplex stainless steel microstructure can be very harmful for its corrosion resistance.No secondary intermetallic phases were found from Fig.8.From Fig.9,it is seen that PDZ has higher hardness values than that of WZ and base material.It is due to therefinementofgrainsoccurringinPDZ.Thiscanbeattributed to work hardening effect caused by heavy deformation in PDZ[19].The measured average grain sizes for BM,PDZ and WZ wereintheorderof28.7,10.6and22.3microns,respectively.In thetensiletestedsample,thefractureoccurredinbasemetaland was away from the joint zone.The tensile strength of the weld jointwasmuchhigherthanthatofthebasemetal(750Mpa).The maximum tensile strength obtained using optimized parameters was 827.17 Mpa and the experimental value was 822 Mpa.The strength and hardness were increased due to grain refinement in WZ and PDZ.It is also revealed that the hardness and tensile strength values are higher in both WZ and PDZ.Fig.10 reveals that the ductile fracturewas appeared and also the dimpleswere found.The dimples were elongated in the stress direction.Due to the high amount of ferrite content and more amounts of dimples present in the weld zone,the tensile strength of the joints is much higher than the strength of base material.Fig.11 reveals that the identified peaks represent the presence of only ferrite and austenite.No other intermetallic phases were identified in the XRD pattern.Less plastic state and slow cooling time restricted the weld zone phase transformation.

    6.Conclusions

    From this investigation the following conclusions are drawn:

    1)The percentage contributions of input parameters are: upsetting pressure=47.05%,friction pressure=25.99%,speed of rotation=21.96%,and burn-off length=4.96%.

    2)PDZ has the finer grain size compared to the weld zone and base material.

    3)Austenite and ferrite phases are present in the weld zone.

    4)The tensile strength of the friction joint(optimized)is higher than the base material strength and the fracture occurred is in ductile nature.

    5)Higher hardness values were obtained in PDZ compared to WZ and BM.This is due to the refinement of grain size in the weld region.

    6)The austenite-to-ferrite ratio is 50:50 for three zones of the weld like BM,PDZ and WZ.

    [1]TMR Stainless.Practical guidelines for the fabrication of duplex stainless steels.London:International Molybdenum Association;2009. p.1-69.

    [2]Jenney Cynthia L,Annette O'Brien.Welding handbookvol.2.Miami F L:AWS:American Welding Society;1991.p.523-30.

    [3]Olson David Leroy,Siewert Thomas A,Liu Stephen,Edwards Glen R. ASM handbook,vol.6.Materials Park:ASM International;1995. p.1130-52.

    [4]Dunkerton SB.Toughness properties of friction welds in steels.Weld J 1986;8:193-201.

    [5]Balamurugan Karupanan,Mishra Mahendra Kumar,Sathiya Paul,Sait Abdullah Naveen.Weldability studies and parameter optimization of AISI 904L super austenitic stainless steel using friction welding.Mater Res 2014;17(4):908-19.

    [6]Mandal S,Sivaprasad PV,Venugopal S,Murthy KPN.Artificial neural network modeling to evaluate and predict the deformation behavior of stainless steel type AISI 304L during hot torsion.Appl Soft Comput J 2009;9(1):237-44.

    [7]Ali Akcayol M,Cinar Can.Artificial neural network based modeling of heatcatalyticconverterperformance.ApplThermEng 2005;25(14-15):2341-50.

    [8]Erdogan Kose.Modelling of colour perception of different age groups usingartificialneuralnetworks.ExpertSystAppl 2008;34(3):2129-39.

    [9]Shabanzadeh Parvaneh,Norazak Senu,Shameli Kamyar,Ismail Fudziah,Maryam Mohagheghtabar.Application of artificial neural network(ANN)for the prediction of size of silver nano particles prepared by green method.Dig J Nanomater Biostructures 2013;8(2):541-9.

    [10]Sathiya P,Aravindan S,Noorul Haq A,Paneerselvam K.Optimization of friction welding parameters using evolutionary computational techniques.J Mater Process Technol 2009;209(5):2576-84.

    [11]Paventhan R,Lakshminarayanan PR,Balasubramanian V.Optimization of friction welding process parameters for joining carbon steel and stainless steel.J Iron Steel Res Int 2012;19(1):66-71.

    [12]Faes Koen,Dhooge Alfred,De Baets Patrick,Van Der Donckt Eric,De Waele Wim.Parameter optimisation for automatic pipeline girth welding using a new friction welding method.Mater Des 2009;30:581-9.

    [13]Luo J,Ye YH,Xu JJ,Luo JY,Chen SM,Wang XC.A new mixedintegrated approach to control welded flashes forming process of damping-tube-gland in continuous drive friction welding.Mater Des 2009;30:353-8.

    [14]Udayakumar T,Raja K,Afsal Husain TM,Sathiya P.Prediction and optimization of friction welding parameters for super duplex stainless steel(UNS S32760)joints.Mater Des 2014;53:226-35.

    [15]Magudeeswaran G,Nair Sreehari R,Sundar L,Harikannan N.Optimization of process parameters of the activated tungsten inert gas welding for aspect ratio of UNS S32205 duplex stainless steel welds.Def Technol 2014;10:251-60.

    [16]Montgomery DC.Design and analysis of experiments.7th ed.Wiley;2009.

    [17]Jain Sanjay K,Archana Sarkar,Vaibhav Garg.Impact of declining trend offlowonHarikeWetland.India.WaterResourManag 2008;22(4):409-21.

    [18]Satyanarayana VV,Madhusudhan Reddy G,Mohandas T.Dissimilar metal friction welding of austenitic-ferritic stainless steels.J Mater Process Technol 2005;160:128-37.

    [19]Satyanarayana VV,Madhusudhan Reddy G,Mohandas T.Continuous drive friction welding studies on AISI 304 austenitic stainless steel welds.MaterManufProcess2004;19:487-505.http://dx.doi.org/ 10.1081/AMP-120038657.

    4 August 2014;revised 16 February 2015;accepted 3 March 2015

    Available online 27 March 2015

    .Tel.:+91 431 2503510;fax:+91 431 2500133.

    E-mail address:psathiya@nitt.edu(P.SATHIYA).

    Peer review under responsibility of China Ordnance Society.

    Copyright?2015,China Ordnance Society.Production and hosting by Elsevier B.V.All rights reserved.

    猜你喜歡
    念書伯父熟人
    校園“老”熟人,我們的成長(zhǎng)大“師”
    理 想
    理想
    和熟人相處之道
    伯父的黃昏戀
    別忘記跟熟人打招呼
    北方人(2018年6期)2018-01-22 19:41:39
    一張油畫
    我們的路
    意林(2013年2期)2013-05-14 16:49:17
    借給伯父一百元
    文苑·感悟(2012年2期)2012-03-22 08:09:44
    叛徒
    国产亚洲5aaaaa淫片| 精品国产乱码久久久久久小说| 久久久欧美国产精品| 新久久久久国产一级毛片| 干丝袜人妻中文字幕| 亚洲人与动物交配视频| 亚洲国产色片| 99久久精品国产国产毛片| 久久精品久久精品一区二区三区| 国产精品一区二区在线不卡| 噜噜噜噜噜久久久久久91| 九色成人免费人妻av| 日本vs欧美在线观看视频 | 欧美变态另类bdsm刘玥| 精品亚洲成a人片在线观看 | 成人黄色视频免费在线看| 中文字幕av成人在线电影| 日韩欧美 国产精品| 免费观看av网站的网址| 久久人妻熟女aⅴ| 日本免费在线观看一区| 最近中文字幕2019免费版| 波野结衣二区三区在线| 97超碰精品成人国产| 男女无遮挡免费网站观看| 又黄又爽又刺激的免费视频.| 久久精品国产a三级三级三级| 国产精品伦人一区二区| 99热全是精品| 国产美女午夜福利| 亚洲欧美一区二区三区国产| 人妻少妇偷人精品九色| 国产精品一区www在线观看| 少妇猛男粗大的猛烈进出视频| 久久6这里有精品| 超碰av人人做人人爽久久| 国产男女内射视频| 一级a做视频免费观看| 大片电影免费在线观看免费| 欧美国产精品一级二级三级 | 一级片'在线观看视频| 国产男人的电影天堂91| 精品午夜福利在线看| 亚州av有码| 国产一区亚洲一区在线观看| 少妇的逼好多水| 18禁裸乳无遮挡动漫免费视频| 黄色欧美视频在线观看| 国产视频内射| 狂野欧美激情性bbbbbb| av女优亚洲男人天堂| 51国产日韩欧美| 欧美xxxx性猛交bbbb| 狂野欧美白嫩少妇大欣赏| 久久精品久久久久久久性| 在线免费观看不下载黄p国产| 在线 av 中文字幕| 午夜福利影视在线免费观看| 插阴视频在线观看视频| 午夜老司机福利剧场| 亚洲精品国产av成人精品| 亚洲精品国产av成人精品| av女优亚洲男人天堂| 少妇人妻精品综合一区二区| 免费看日本二区| 日日啪夜夜爽| 国产一区有黄有色的免费视频| 久久精品国产a三级三级三级| 最近最新中文字幕大全电影3| 秋霞伦理黄片| 日韩,欧美,国产一区二区三区| 国产精品一区二区三区四区免费观看| 老司机影院毛片| 高清黄色对白视频在线免费看 | 高清黄色对白视频在线免费看 | 深夜a级毛片| 看非洲黑人一级黄片| 亚洲国产av新网站| 午夜激情久久久久久久| 一边亲一边摸免费视频| 偷拍熟女少妇极品色| 久久97久久精品| 亚洲综合精品二区| 黄色视频在线播放观看不卡| 国产美女午夜福利| 国产精品蜜桃在线观看| 高清在线视频一区二区三区| 精品一品国产午夜福利视频| 国产精品.久久久| 成人亚洲精品一区在线观看 | 99久久精品热视频| 在线观看三级黄色| 亚洲国产欧美人成| 国产高清有码在线观看视频| 欧美xxxx黑人xx丫x性爽| 最近2019中文字幕mv第一页| 中文字幕亚洲精品专区| 国产大屁股一区二区在线视频| 成人一区二区视频在线观看| 成人综合一区亚洲| 欧美97在线视频| 成人美女网站在线观看视频| 中文字幕久久专区| 中文字幕人妻熟人妻熟丝袜美| 国产av一区二区精品久久 | 日韩制服骚丝袜av| xxx大片免费视频| 舔av片在线| 少妇猛男粗大的猛烈进出视频| 毛片女人毛片| 美女高潮的动态| 精品亚洲成a人片在线观看 | 国产精品三级大全| 一区二区三区四区激情视频| 一本久久精品| 黄色视频在线播放观看不卡| 超碰97精品在线观看| 中文字幕亚洲精品专区| 成人综合一区亚洲| 国产精品99久久久久久久久| 亚洲色图av天堂| 波野结衣二区三区在线| 国内少妇人妻偷人精品xxx网站| 男女边摸边吃奶| 国产一区二区在线观看日韩| 国产精品一及| av一本久久久久| 亚洲色图综合在线观看| 一级毛片我不卡| 亚洲国产色片| 午夜精品国产一区二区电影| 免费黄网站久久成人精品| 免费少妇av软件| 亚洲无线观看免费| 久久久久久久久久人人人人人人| 久久人妻熟女aⅴ| 亚洲欧美日韩无卡精品| 国产欧美另类精品又又久久亚洲欧美| 91久久精品国产一区二区三区| 精品少妇久久久久久888优播| 99久久精品国产国产毛片| 久久精品国产亚洲网站| 亚洲人成网站在线播| 亚洲三级黄色毛片| 色5月婷婷丁香| 国产男女超爽视频在线观看| 亚洲精品自拍成人| 欧美bdsm另类| 一本—道久久a久久精品蜜桃钙片| 丰满迷人的少妇在线观看| 久久久成人免费电影| 亚洲精品乱久久久久久| 国产精品爽爽va在线观看网站| 丰满少妇做爰视频| 久久精品夜色国产| 久久精品国产鲁丝片午夜精品| 国产亚洲av片在线观看秒播厂| 欧美老熟妇乱子伦牲交| 成人一区二区视频在线观看| 免费看光身美女| .国产精品久久| 只有这里有精品99| 亚洲天堂av无毛| 我的老师免费观看完整版| 精品熟女少妇av免费看| 熟女电影av网| 中文欧美无线码| 免费不卡的大黄色大毛片视频在线观看| 一级毛片久久久久久久久女| 王馨瑶露胸无遮挡在线观看| 麻豆精品久久久久久蜜桃| 在线观看人妻少妇| 麻豆成人av视频| 亚洲精品aⅴ在线观看| 最新中文字幕久久久久| 日产精品乱码卡一卡2卡三| 九九爱精品视频在线观看| 色婷婷久久久亚洲欧美| 精品亚洲成国产av| 身体一侧抽搐| 永久网站在线| 亚洲欧美日韩无卡精品| 亚洲精品乱久久久久久| 美女内射精品一级片tv| 亚洲精品第二区| 一级爰片在线观看| 免费av不卡在线播放| 街头女战士在线观看网站| 精品久久久久久久久av| 一本久久精品| 国产精品伦人一区二区| 99热这里只有是精品50| 伦理电影大哥的女人| 亚洲精品自拍成人| 99久久中文字幕三级久久日本| 亚洲,一卡二卡三卡| 久久人人爽人人片av| 国产精品麻豆人妻色哟哟久久| 日韩国内少妇激情av| 亚洲一区二区三区欧美精品| 日本wwww免费看| 观看美女的网站| 少妇 在线观看| 日日啪夜夜爽| 18+在线观看网站| 免费观看在线日韩| 99久久精品热视频| 免费大片18禁| 熟女人妻精品中文字幕| 亚洲图色成人| 99国产精品免费福利视频| 国产免费福利视频在线观看| 亚洲av男天堂| 久久国产亚洲av麻豆专区| 最近的中文字幕免费完整| 联通29元200g的流量卡| 日本免费在线观看一区| 精品国产三级普通话版| 亚洲精品国产色婷婷电影| 国产精品熟女久久久久浪| 国产乱人偷精品视频| 色综合色国产| 国产一区二区三区av在线| 日韩强制内射视频| 国产 精品1| 观看免费一级毛片| 国产成人精品福利久久| 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美在线一区| 国产成人免费无遮挡视频| 亚洲av综合色区一区| 国产成人午夜福利电影在线观看| 日本免费在线观看一区| 婷婷色综合大香蕉| 国产综合精华液| 国产成人aa在线观看| 少妇精品久久久久久久| 人体艺术视频欧美日本| 成人黄色视频免费在线看| 国产亚洲最大av| 久久国内精品自在自线图片| 欧美高清性xxxxhd video| 亚洲国产日韩一区二区| 日本av免费视频播放| 成人影院久久| 久久久色成人| 久久婷婷青草| 国产高潮美女av| 免费久久久久久久精品成人欧美视频 | 亚洲精华国产精华液的使用体验| 制服丝袜香蕉在线| 亚洲av.av天堂| 18禁裸乳无遮挡免费网站照片| 纯流量卡能插随身wifi吗| 亚洲久久久国产精品| av国产免费在线观看| 免费观看性生交大片5| 极品教师在线视频| 亚洲最大成人中文| 亚洲真实伦在线观看| 一级毛片电影观看| 色婷婷久久久亚洲欧美| 中文乱码字字幕精品一区二区三区| 精品一区二区三区视频在线| 日本欧美视频一区| 久久人妻熟女aⅴ| 日本猛色少妇xxxxx猛交久久| 日本欧美国产在线视频| 少妇精品久久久久久久| 午夜福利在线观看免费完整高清在| 国产成人91sexporn| 能在线免费看毛片的网站| 高清在线视频一区二区三区| 一级毛片黄色毛片免费观看视频| 精品国产三级普通话版| 天天躁日日操中文字幕| 亚洲精品中文字幕在线视频 | 久久久久性生活片| 80岁老熟妇乱子伦牲交| 亚洲精品视频女| 国产男女内射视频| 狂野欧美激情性bbbbbb| 亚洲精华国产精华液的使用体验| 亚洲精品乱码久久久久久按摩| 亚洲自偷自拍三级| 在线免费十八禁| 国产一区有黄有色的免费视频| 精品国产一区二区三区久久久樱花 | 欧美精品一区二区免费开放| 午夜免费鲁丝| 九九爱精品视频在线观看| 欧美日韩综合久久久久久| 18禁裸乳无遮挡动漫免费视频| 美女福利国产在线 | 国产av码专区亚洲av| 欧美bdsm另类| 91久久精品国产一区二区成人| 欧美高清性xxxxhd video| 国产高清不卡午夜福利| 99久久中文字幕三级久久日本| 国产精品久久久久成人av| 精品视频人人做人人爽| 日韩av在线免费看完整版不卡| 国产亚洲最大av| 久久精品久久久久久久性| 国产成人freesex在线| 免费播放大片免费观看视频在线观看| 久久99热这里只频精品6学生| 蜜桃亚洲精品一区二区三区| 亚洲成色77777| 我的女老师完整版在线观看| 97热精品久久久久久| 99久久人妻综合| 欧美三级亚洲精品| 国国产精品蜜臀av免费| 久久久久性生活片| 精品酒店卫生间| 一级毛片 在线播放| 女性被躁到高潮视频| 亚洲精品一二三| 最黄视频免费看| 一二三四中文在线观看免费高清| 国产精品久久久久久久久免| 亚洲色图av天堂| 菩萨蛮人人尽说江南好唐韦庄| 日韩欧美精品免费久久| tube8黄色片| 男的添女的下面高潮视频| 中国国产av一级| 亚洲精品国产av蜜桃| 国产男女超爽视频在线观看| 国产成人精品婷婷| 亚洲精品456在线播放app| 国产在线免费精品| 日韩中文字幕视频在线看片 | 老司机影院毛片| 久久人人爽人人爽人人片va| 国产老妇伦熟女老妇高清| 国产精品麻豆人妻色哟哟久久| 男女边摸边吃奶| 嫩草影院入口| 久久久久久久久久久丰满| 欧美xxⅹ黑人| 精品久久国产蜜桃| 一级黄片播放器| 黄色一级大片看看| 18禁裸乳无遮挡动漫免费视频| 国产 精品1| 伊人久久精品亚洲午夜| 久久久久久久久大av| 我要看黄色一级片免费的| 国产伦精品一区二区三区四那| 这个男人来自地球电影免费观看 | 欧美xxxx性猛交bbbb| 高清视频免费观看一区二区| 少妇人妻 视频| 男人爽女人下面视频在线观看| 欧美激情国产日韩精品一区| 久久精品久久久久久噜噜老黄| 性色avwww在线观看| 国产色爽女视频免费观看| 精品人妻偷拍中文字幕| 日韩成人av中文字幕在线观看| 高清毛片免费看| 国产女主播在线喷水免费视频网站| 搡老乐熟女国产| 中文在线观看免费www的网站| 午夜福利网站1000一区二区三区| 一级毛片久久久久久久久女| 又爽又黄a免费视频| 午夜福利在线在线| 插阴视频在线观看视频| 国产成人午夜福利电影在线观看| 免费观看a级毛片全部| 亚洲四区av| 伦精品一区二区三区| 一区二区av电影网| av女优亚洲男人天堂| 欧美+日韩+精品| 久久久久视频综合| 三级国产精品片| 欧美区成人在线视频| 舔av片在线| 国产久久久一区二区三区| 纵有疾风起免费观看全集完整版| 搡老乐熟女国产| 国产在线一区二区三区精| 22中文网久久字幕| 久久97久久精品| 国产在线男女| 日本wwww免费看| 99热全是精品| 国产精品久久久久久久电影| 久久影院123| 国产大屁股一区二区在线视频| 国产黄色免费在线视频| 深爱激情五月婷婷| 精品99又大又爽又粗少妇毛片| 精品人妻偷拍中文字幕| 啦啦啦视频在线资源免费观看| 国产日韩欧美在线精品| 久久国产乱子免费精品| 女性被躁到高潮视频| 亚洲最大成人中文| 亚洲内射少妇av| 九九爱精品视频在线观看| 久久精品人妻少妇| 国模一区二区三区四区视频| 免费看不卡的av| 91aial.com中文字幕在线观看| 青青草视频在线视频观看| 成年av动漫网址| 91久久精品电影网| 久久综合国产亚洲精品| 国产一区有黄有色的免费视频| 日本av免费视频播放| 日韩亚洲欧美综合| 不卡视频在线观看欧美| 亚洲内射少妇av| 熟妇人妻不卡中文字幕| 国产永久视频网站| 嫩草影院新地址| 亚洲丝袜综合中文字幕| 在线观看一区二区三区| 国产人妻一区二区三区在| 欧美激情极品国产一区二区三区 | 青春草亚洲视频在线观看| 嘟嘟电影网在线观看| 亚洲精品aⅴ在线观看| 亚洲欧美日韩卡通动漫| a级一级毛片免费在线观看| 久久99热这里只有精品18| 热re99久久精品国产66热6| 午夜福利网站1000一区二区三区| 亚洲成人手机| 久久av网站| 一级a做视频免费观看| 亚洲精品一区蜜桃| 日本与韩国留学比较| 五月天丁香电影| 国产一区二区三区综合在线观看 | 天美传媒精品一区二区| 欧美一级a爱片免费观看看| 日韩不卡一区二区三区视频在线| 免费观看的影片在线观看| 性高湖久久久久久久久免费观看| 久久久精品94久久精品| 国产在线一区二区三区精| 亚洲一区二区三区欧美精品| 欧美日韩国产mv在线观看视频 | 国产精品欧美亚洲77777| 成年av动漫网址| 一区二区三区精品91| 在线免费十八禁| 一个人看的www免费观看视频| 美女xxoo啪啪120秒动态图| av在线app专区| 在线观看一区二区三区| 国产乱来视频区| 国产免费一级a男人的天堂| 国产免费视频播放在线视频| 久久久精品免费免费高清| 亚洲怡红院男人天堂| 天堂俺去俺来也www色官网| 人妻夜夜爽99麻豆av| av卡一久久| 伊人久久精品亚洲午夜| 国产高清国产精品国产三级 | 久久久成人免费电影| 亚洲av综合色区一区| 如何舔出高潮| 一级爰片在线观看| 亚洲av成人精品一二三区| 少妇人妻精品综合一区二区| 全区人妻精品视频| 美女福利国产在线 | 久久久午夜欧美精品| 能在线免费看毛片的网站| 国产在线视频一区二区| 我要看日韩黄色一级片| 在线观看一区二区三区激情| 少妇丰满av| 久久99热6这里只有精品| 中文字幕av成人在线电影| 国产老妇伦熟女老妇高清| 夜夜看夜夜爽夜夜摸| 啦啦啦视频在线资源免费观看| 成人综合一区亚洲| 少妇丰满av| 亚洲精品成人av观看孕妇| 全区人妻精品视频| 国产精品久久久久久av不卡| 精品少妇黑人巨大在线播放| 激情 狠狠 欧美| 日韩欧美一区视频在线观看 | 午夜激情久久久久久久| 国产 一区 欧美 日韩| 欧美bdsm另类| 妹子高潮喷水视频| 国产精品国产三级专区第一集| 欧美xxxx性猛交bbbb| 国产女主播在线喷水免费视频网站| 免费观看的影片在线观看| 九九在线视频观看精品| 男男h啪啪无遮挡| 欧美激情极品国产一区二区三区 | 亚洲精品国产av成人精品| 黄色日韩在线| 国产亚洲午夜精品一区二区久久| 欧美97在线视频| 日日撸夜夜添| 精品人妻一区二区三区麻豆| 国产伦精品一区二区三区视频9| 国产欧美日韩一区二区三区在线 | 国产乱人偷精品视频| 精品国产三级普通话版| 国产欧美日韩一区二区三区在线 | 一个人免费看片子| av播播在线观看一区| 免费在线观看成人毛片| 天天躁夜夜躁狠狠久久av| 久久国内精品自在自线图片| 美女脱内裤让男人舔精品视频| 王馨瑶露胸无遮挡在线观看| 看免费成人av毛片| 99热这里只有是精品50| 免费播放大片免费观看视频在线观看| 晚上一个人看的免费电影| 亚洲精品日韩av片在线观看| 欧美3d第一页| 亚洲精品久久久久久婷婷小说| 人妻少妇偷人精品九色| 国产免费一级a男人的天堂| 国产精品嫩草影院av在线观看| 成年美女黄网站色视频大全免费 | 久久精品人妻少妇| 51国产日韩欧美| 婷婷色麻豆天堂久久| 亚洲精品成人av观看孕妇| 一级毛片黄色毛片免费观看视频| 亚洲欧美清纯卡通| 色网站视频免费| 性色avwww在线观看| 成人亚洲欧美一区二区av| av卡一久久| 老司机影院毛片| 丰满少妇做爰视频| 插逼视频在线观看| 午夜福利影视在线免费观看| 精品熟女少妇av免费看| 亚洲四区av| 国内少妇人妻偷人精品xxx网站| 欧美xxxx性猛交bbbb| 2022亚洲国产成人精品| 建设人人有责人人尽责人人享有的 | 色网站视频免费| 亚洲精品国产av成人精品| 午夜福利在线观看免费完整高清在| 日韩亚洲欧美综合| 亚洲av日韩在线播放| 亚洲四区av| 国产免费一级a男人的天堂| 麻豆成人午夜福利视频| 91久久精品电影网| 国产中年淑女户外野战色| 自拍偷自拍亚洲精品老妇| 91精品国产九色| 国产精品av视频在线免费观看| 亚洲不卡免费看| 一级毛片 在线播放| 亚洲国产毛片av蜜桃av| 欧美精品国产亚洲| 国产淫语在线视频| 亚洲欧洲日产国产| 蜜桃在线观看..| 久久精品久久精品一区二区三区| 欧美激情极品国产一区二区三区 | 久久精品夜色国产| 久久99蜜桃精品久久| 亚洲精品456在线播放app| 精品国产乱码久久久久久小说| 97在线人人人人妻| 一本色道久久久久久精品综合| 久久国产精品男人的天堂亚洲 | av在线播放精品| 人妻制服诱惑在线中文字幕| 日本爱情动作片www.在线观看| 大香蕉久久网| 王馨瑶露胸无遮挡在线观看| 99精国产麻豆久久婷婷| 大香蕉97超碰在线| 国产高潮美女av| 国产精品99久久99久久久不卡 | 国产精品欧美亚洲77777| 人人妻人人看人人澡| 国产精品欧美亚洲77777| 日日摸夜夜添夜夜添av毛片| 最近中文字幕高清免费大全6| 国产精品一区二区三区四区免费观看| 亚洲婷婷狠狠爱综合网| 免费观看av网站的网址| 简卡轻食公司| 国产亚洲5aaaaa淫片| 在线亚洲精品国产二区图片欧美 | 午夜免费观看性视频| 一级二级三级毛片免费看| 色综合色国产| 99久国产av精品国产电影| 这个男人来自地球电影免费观看 | 黑丝袜美女国产一区| 亚洲四区av| 91午夜精品亚洲一区二区三区| 成人二区视频| 少妇人妻久久综合中文| 精品一区二区三卡|