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

    Dependence on forming parameters of an integral panel during the electromagnetic incremental forming process

    2018-07-24 09:00:16JinqiangTANMeiZHANHongweiLI
    CHINESE JOURNAL OF AERONAUTICS 2018年7期

    Jinqiang TAN,Mei ZHAN,Hongwei LI

    State Key Laboratory of Solidification Process,School of Materials Science&Engineering,Northwestern Polytechnical University,Xi’an 710072,China

    KEYWORDS Electromagnetic incremental forming(EMIF);Forming uniformity;Integral panel;Parameters;Shape deviation

    Abstract Nowadays,more and more attentions are paid to electromagnetic incremental forming(EMIF),especially for a part with a large-scale size,e.g.,an integral panel with stiffened ribs.In this work,the bending of a panel into a double-curvature profile via EMIF is carried out experimentally and evaluated by comparing the formed profile with the desired profile.During the process,discharges at four positions along different discharge paths are designed.The effects of forming parameters on the die-fittingness of the workpiece are discussed,for which two evaluation indices are used to judge forming results.The results show that a discharge voltage in an incremental mode is helpful to improve thefittingness and avoid the collision rebound against the die at the same time.Discharging at the diagonal positions with the ‘X” discharge path exhibits the minimal shape deviation and the best forming uniformity.On the contrary,discharging at the parallel positions with the ‘Z”discharge path obtains the worst forming quality.Overlap of the coil at different positions should be given during EMIF;however,a lower overlap rate of the coil helps improve the forming quality.The results obtained in this work are useful for forming integral panels with stiffened ribs via the EMIF process.

    1.Introduction

    Integral panels,with high-stiffener having advantages of light weight,high stiffness,and high structural efficiency,have become important structural bearings and widely used in aircraft.However,because of their large sizes,net reinforced structures and structural stiffness,and requirements of good shape precision,the manufacture of integral panels has a high demand.1The conventional plastic forming methods for integral panels contain shot peen forming,2creep age forming,3press-bending,4roll forming,etc.Some successful examples and applications have been introduced in the ‘Integral Airframe Structures(IAS)”program.5However,due to the existence ofstiffened ribs,there are many difficulties in processing compared to flat plate forming,and limitations can also be enumerated.Garie′py et al.6pointed out that for shot peen forming,the process needs complex designs mainly based on experience,and sometimes additional prestress is necessary.According to Zhan et al.,7creep age forming is limited to certain materials or occasions,and excessive time is spent in the process time including subsequent springback calibration.Since punch forces are applied directly on ribs,mentioned by Yan et al.,8instability and springback phenomena become the main problems during the brake forming process.In a word,the problems in the traditional forming methods of integral panels prompt researchers to explore novel forming technologies.

    Electromagnetic forming(EMF)is an impulse or highspeed forming process using a pulsed magnetic field to apply Lorentz forces to workpieces preferably made of a highly electrically conductive material.Psyk et al.9reviewed that the formability limit can be increased and the production cycle can be shortened due to its high forming velocity.However,considering the large sizes of some parts,traditional EMF is difficult to apply directly because of the limited action scope and the depth of the magnetic field force from the coil with a fixed position.Therefore,some new methods have been proposed to overcome these difficulties.Kamal and Daehn10combined EMF and stamping via drawing to form a sheet with local features successfully.Lai et al.11utilized a large coil,of which the diameter is 860 mm,to form a large-scale sheet with a diameter of 1378 mm successfully.However,a large coil at a f i xed position was still adopted,and it is difficult to form parts with complex structures directly.

    In order to deal with the above challenges for large or complex parts,the idea of incremental forming was introduced by some researchers.Cui et al.12proposed an incremental forming method during EMF,called electromagnetic incremental forming(EMIF),to form large-scale sheets.In EMIF,a workpiece accumulates into the desired profile andfits well with the die gradually from one position to the adjacent via a small coil with small discharge energy moving and rotating along some certain paths in the space according to the requirement.In their study,comparing with traditional EMF,EMIF introduces some new forming parameters,such as discharge position,discharge pass,discharge path,coil overlap rate,etc.,which result in complex coupling effects.Therefore,EMIF is a more complicated process than traditional EMF.At present,only a few researchers have conducted studies on EMIF.Zhao et al.13analyzed the forming uniformity under different coil overlap rates and different loading paths during the electromagnetic incremental bulging process of tubes.The results of their study showed that the higher the coil overlap rate,the better the forming uniformity,whereas the lower forming ef ficiency.Cui et al.14analyzed the effects of coil moving path,discharge pass,and other discharge parameters on electromagnetic incremental deep drawing of sheet metals.They reported that sheet components with goodfittingness to the die were obtained by adopting reasonable process parameters.Xiong et al.15realized large-scale sheet forming via moving a coil along the axial direction.Tan et al.16investigated the EMIF process of several integral panels with different distributions of ribs.They established an FE model for the panels during the EMIF process,and discussed the distribution of electromagnetic force,the evolvement rule of velocity and height of the panels with different rib arrangements,as well as a comparison of forming quality between the panels obtained by brake forming and EMIF.

    The characteristics of large sizes,high ribs,high structural rigidity,and complex and staggered high ribs grid structures of integral panels lead to that the effects of forming parameters in EMIF are different from those of pipes and metal sheets.Therefore,it is vital to study the influences of parameters in the EMIF process for designing an optimized forming scheme to obtain the ideal profile of an integral panel.This paper studies the effects of forming parameters,such as discharge parameters,geometric parameter,and discharge path on thediefi ttingness of the workpiece to the die during EMIF,and the shape deviation and the forming uniformity are also evaluated.

    2.Experimental arrangement

    2.1.Experimental setup

    The aim of utilizing EMIF in this study is to realize the bending of a semi-finished integrally stiffened panel(by milling on a plate to a thin web with several ribs)to thefinial profile with double curvatures in two directions.The experimental installation for the integral panel is designed by Zhan et al.17and shown in Fig.1.In order to facilitate the movement of the coil in the EMIF process,it is fixed between two support trays and bolted to a platen,where the platen is pressed by a hydraulic machine and can move horizontally on the blank holders to realize the incremental movement.The distance between the coil and the blank is adjusted through changing the height of the pads set between the platen and the coil.Fig.2 shows the schematic of the forming process.

    In this study,the material used in the EMIF experiments of the integral panels(Fig.2)is 2A12-T4 aluminum alloy.18The web thickness and rib width of all the blanks are both 2 mm.A f l at spiral coil is adopted,in which a copper wire,with a section area of 4 mm×2 mm,is filled of resin with the roles of insulation and reinforcement.The wire is twined to three layers,and each layer has 10 winding turns wire.The inner radius is 7 mm and the gap between two adjacent layers is 0.4 mm.Two paths along theXandYaxes are defined to observe forming displacements,which will be discussed in Section 3.In the EMIF experiments,four discharge positions are designed,as shown in Fig.2.A die with double curvatures is used in the EMIF process,as shown in Fig.3.

    Fig.1 Installation schematic diagram of EMIF.

    Fig.2 Schematic of the EMIF process principle.

    Fig.3 Schematic diagram of the die.

    2.2.Experimental procedures

    In this study,the effects of process parameters,such as discharge voltage,discharge capacitance,ratio of rib-height to web-thickness,discharge path,and overlap rate of the coil,are obtained via the experimental setup described in the above section.

    Fig.4 Diagram of four discharge positions.

    Considering the action area of the electromagnetic forces and the sizes of the panel,four discharge positions with a rectangular distribution are designed,as shown in Fig.4,where 1,2,3,and 4 are discharge positions.Therefore,there are three different discharge paths,the ‘X” discharge path(the order of discharged at the diagonal position),the ‘U” discharge path(the order of discharged along the circumference)and the ‘Z”discharge path(the order of discharged at parallel positions),as shown in Fig.5.The discharge capacitance of the experiments is 160 μF,and the discharge paths and the position coordinates are shown in Tables 1 and 2,respectively.

    Generally,the curvature radii of the integral panels are very large,which means that the forming depths are shallow.Therefore,there is only one single discharge for each position in this study.

    2.3.Quantitative characterization methods

    To provide a good evaluation of the forming precision of the panel,the formed profiles along two paths(Fig.2)are observed,and two evaluation indices of feature points in the formed profile(intersection points of the marked lines on the surface without ribs),including the shape deviation(the diefi ttingness gap)to the expected profile and the forming uniformity,are proposed.

    The shape deviation of the panel is evaluated by the root mean square error(RMSE),i.e.,

    wherefmea,iis the measured value of the forming depth at feature pointi(i=1,2,...,n)along path X or path Y,nthe number of feature points,andfdes,ithe desired value of the forming depth at pointi,i.e.,the depth of the die corresponding to the feature point(Fig.6(a)).Fig.6(b)shows the measured data of these two indices,whereRmea,iandRavehave an identical curvature center with the desired profile(R).

    The forming uniformity degree of the panel is represented by the standard deviation of the shape deviation SD,i.e.,

    whereRmea,iis the measured distance between feature pointi(i=1,2,...,n)on the experimental panel and the curvature center of the die along path X or path Y,nthe number of feature points,andRavethe average value ofRmea,i.

    3.Results and discussion

    The effects of five parameters,two electromagnetic parameters(discharge voltage and discharge capacitance),one geometric parameter(ratio of rib-height to web-thickness),and two positional parameters(discharge path and overlap rate of the coil)on the forming displacement in theZdirection are studied.The die-fittingness of the panels to the die and the forming uniformity are evaluated.

    3.1.Effect of discharge voltage

    The discharge voltage is one of the most important parameters in the EMF process,and it is useful to improve the formability due to increased discharge energy(Psyk et al.9;Cui et al.12).The effect of discharge voltage for the workpiece on the diefi ttingness under the same capacitance(C=320 μF)and the same ratio of rib-height to web-thickness(5:1)is studied.The ‘X” discharge path is chosen for the four discharges,and the voltages at different positions are listed in Table 3.The discharge sequences and two paths defined to measure the profiles are shown in Fig.7(a).From Fig.7(b),the profile along path X indicates that the formed profile is uneven,and the largest forming depths occur in the region near the free edge along the negative direction of theXaxis,because the initial discharge voltage is lower than those of the subsequent three discharges.The incremental voltages in the second and third discharges enhance the formability.Since the workpiece has no contact with the die in the deformation region in the first set of experiments,a higher initial voltage is adopted in the subsequent experiments.It can be found in Fig.7(c)that the forming depths under higher initial voltages in the rest two sets of experiments get worse,especially in the third set of experiments.The reason is that a higher voltage makes the workpiece undergoes collision rebound,as shown in Fig.7(d).profile curves in Fig.7(b)and(c)show that the uniformity and forming depth are improved with a lower initial discharge.From Fig.7,the largest forming depth is mainly determined by the subsequent three discharges.The error analysis in Fig.8 shows that the shape deviation(RMSE)increases with the voltage groups,and the forming uniformity in the first group is the best,followed by the third,and the worst in the second group.

    Fig.5 Schematic diagrams of three discharge path modes.

    Table 1 Detailed discharge paths in EMIF experiments.

    Table 2 Discharge positions in EMIF experiments.

    Fig.6 Distribution graph of feature points and measured data in RMSE and SD for the formed profile.

    Table 3 Different discharge voltage groups along discharge path 1→4→3→2.

    Since the die used in this study is useful to expel the air,the problem occurred in the study of Cui et al.14can be neglected,where vent holes are necessary for a conical die with a flat bottom,and it is hard to improve the forming depth by increasing the voltage due to the resistance of high-speed compressed air.Thus,as shown in Fig.7,it is possible to obtain good forming quality via selecting suitable voltage groups for the four discharge positions.In addition,the uniformity of the forming depth in Fig.7(b)is lower than that in Fig.7(c),and the distance between two adjacent positions along path Y is larger than that along path X,which indicates that the discharge voltage must also be designed carefully.

    3.2.Effect of discharge capacitance

    The effect of the discharge capacitance with the same discharge energy at each position is studied,and the conditions are shown in Table 4.The same ratio of rib-height to webthickness(2:1)and discharge path(the ‘X” path)are used.Fig.9 shows the final profiles at different discharge capacitances along path X and path Y.It can be observed from Fig.9(a)that under the same discharge energy,the profile curves express different results.The workpiece shows a similar f i ttingness for capacitances equal to 80 μF and 160 μF,and a goodfit occurs in the region near the free end in the positive direction of theXaxis.However,the forming depth for a capacitance of 320 μF expresses an opposite phenomenon,where the largest depth is distributed in the region of the negative direction of theXaxis,and the die-fittingness ofC=320 μF is worse than those of the other two conditions.Fig.9(b)also shows that similar depth distributions and better f i t abilities for the two smaller capacitances.The shape deviation and forming uniformity of the surface profile in Fig.10 also indicate that largest errors come from the condition whenC=320 μF,and the obtained best shape deviation and forming uniformity are fromC=80 μF andC=160 μF,respectively.It is found that a lower capacitance and higher voltages are useful to improve the forming quality.

    The discharge capacitance is relative to the frequency of the discharge current and the skin depth,9and thus,it is necessary to select a reasonable capacitance in this study.Although the discharge energy is the same for each position of the three capacitances,the skin depth of the induced eddy is proportional to the capacitance;therefore,a larger capacitance increases the skin depth.In this study,the thickness of the web is 2 mm,and the skin depths for the three capacitances from small to large are about 1.73 mm,2.06 mm,and 2.45 mm,respectively.A skin depth near the thickness of the web means a higher energy efficiency.Therefore,from an economic standpoint,the capacitance should be set to 80 μF or 160 μF in this study.It can be found that the forming depths along the two paths are very close atC=80 μF andC=160 μF,and the worst forming depth is obtained atC=320 μF;therefore,the intermediate value of the capacitanceC=160 μF may be the most suitable option because the forming quality in this condition requires lower voltages.The study of the effect of discharge capacitance on forming of an integral panel is similar to that by Guo et al.,19where a panel was expected to be formed to a profile with a single curvature with three discharge positions.

    3.3.Effect of ratio of rib-height to web-thickness

    Obviously,the geometric stiffness of the workpiece is largely influenced by the geometric sizes of the panel,such as the height of the rib and the thickness of the web,especially the rib.In addition,the simulation research by Tan et al.16showed that there are electromagnetic forces on the rib and the web at the same time,and the largest forces on the rib are higher than those on the web.Therefore,it is important to discuss the effect of these sizes.Here,the ratio of rib-height to webthickness is selected as an observed variant.The consistency conditions for the experiments are the same capacitance(C=320 μF),discharge path(the ‘X”discharge path),the distance(6 mm)between the coil and the top of the rib,and discharge voltages of 4.5 kV,5 kV,6 kV,and 5.5 kV for positions 1,2,3,and 4,respectively.Considering an increased rib-height,the deformation for the workpiece becomes much more dif ficult.Therefore,the forming depth in Fig.11 shows inverse proportional to the rib-height generally for both of the two observed paths.Lei et al.20suggested that normal electromagnetic forces increase linearly on the ribs,and approximately remain unchanged on the web with an increase of the ribheight when discharging in the center of the integral panel once.That means the drive of deformation for the integral panel changes from the electromagnetic forces on the web at lower ribs to the electromagnetic forces on the ribs at higher ribs;however,the increased rib-height in turn results in greater deformation resistance and smaller electromagnetic forces on the web.

    The shape deviation and forming uniformity in Fig.12 both show that the forming quality gets worse at a high ratio of ribheight to web-thickness.The best results are obtained when the ratio equals to 3:1,not 2:1,which can be explained by that there occurs collision rebound when the ratio equals to 2:1 under the conditions in this study.

    Although the distance between the rib and the coil is identical among all the five ratios,the electromagnetic forces under the five conditions are different.The forming quality generally gets worse at high ratios,excepth/t=4:1,and this reduction of the forming depth is not linear.Therefore,varied suitable discharge conditions such as voltage groups also need to be devised according to the ratio.Considering the exception of a ratio ofh/t=4:1,it may be because the change of the geometric stiffness is nonlinear with an increase of the ratio.

    3.4.Effect of discharge path

    Fig.7 Results at different discharge voltage combinations.

    In the EMIF process,the discharge path has an important influence on the forming result.In order to design reasonable discharge paths for obtaining high forming precision,it is necessary to study the effect of the discharge path on the forming process.The forming depths of workpieces along path X and path Y under different discharge paths in Table 1 are shown in Fig.13,where the discharge voltages at the four discharge positions 1,2,3,and 4 are 4 kV,5 kV,6 kV,and 6.5 kV,respectively.The distances between the positions parallel to theYaxis and theXaxis are 70 mm and 50 mm,respectively.The capacitance is set to 320 μF,and the ratio of rib-height to webthickness is 2:1.From Fig.13(a),the discharge path mode has an effect on the forming depth.The overall forming depth under the ‘X” discharge path mode(discharge path-1 and 2)is the largest,followed by that of the ‘U” discharge path mode(discharge path-3 and 4),and the forming depth under the ‘Z”discharge path mode(discharge path-5 and 6)is the smallest.Except discharge path-1,the forming depth variation laws of the other five discharge paths are consistent:along the positive direction of theXaxis,the forming depth firstly increases rapidly,and then decreases slowly.For discharge path-1,the trend of the depth curve is opposite to those of the other paths.The difference between discharge path-1 and path-2,both of which under the same discharge path mode(the ‘X”discharge path),comes from the position selected in the third discharge.For discharge path-1,the position of the third discharge is located at position 3,while for discharge-2,at position 2.Therefore,the third discharge at position 2(X>0)results in a larger forming depth in the regionX>0 than that in the regionX<0,and so does the discharge at position 3.From Fig.13(a),the workpiece along discharge path-1 achieves the maximal forming depth and the bestfit with the die.The differences of the forming depth along path Y in Fig.13(b)are not obvious,and the largest depth is obtained when discharging along path-2.

    Fig.8 Error analysis of different voltage combinations during EMIF.

    Due to the voltage increasing gradually during the four discharges,the previous deformation region can be influenced by the subsequent deformation,as shown in Fig.13(a)where the forming quality varies on both sides of theXaxis,and the depth atX<0 is better than that atX>0.Meantime,the forming profiles in the same path mode,such as path-1 and path-2(both mode ‘X”),are also different from each other.That is because the panel is rectangular,and the distances between two adjacent discharge positions along the length direction(e.g.,distance between positions 1 and 2)and the width direction(e.g.,distance between positions 2 and 3)influence deformation results.Furthermore,for path mode ‘X”,the voltages at the subsequent positions atX>0 should be enhanced to improve the forming quality in this region when discharge path-1 is chosen.Conversely,the voltages should be reduced at the subsequent positions atX>0 in order to weaken the influence from the subsequent discharges on the previous ones when discharge path-2 is chosen.The effect of the discharge path on the shape of the panel indicates that for a large workpiece with a non-circular shape,it is necessary to design a suitable discharge path in the EMIF process.

    The shape deviations and the forming uniformity of the workpiece under different discharge paths are obtained,as shown in Fig.14.From Fig.14,it can be found that the workpiece under an ‘X” discharge path achieves the bestfit with the die.It is illustrated that the diagonal movement of the coil can reduce uncoordinated deformation,as well as the influence ofthe deformation of the subsequent passes on the formed area.The shape deviations of the workpiece under the ‘Z” discharge path are the worst comparing with those under the other two discharge path modes.Fig.14 also expresses that the general trends of the forming uniformity of these workpieces tend to be consistent with that of the shape deviation.

    Table 4 Different discharge capacitances under the same energy along discharge path 1→4→3→2.

    Fig.9 Panel profiles at different discharge capacitances along two paths.

    3.5.Effect of overlap rate of coil

    Fig.10 Error analysis at different discharge capacitances during EMIF.

    In the EMIF process,there may be overlap for the coil between two adjacent discharge positions,which means that the latter discharge has an effect on the former one.In this study,the concept of overlap rate represents the percentage of the overlap area for the coil at two positions along the width of the workpiece of the total coil area.Three overlap rates are designed,and the position coordinates of the coil with these three overlap rates are shown in Table 5.Fig.15 indicates that a smaller overlap rate is useful to improve the ability of die-fittingness along both paths.Therefore,a larger overlap results in a smaller effective action area,which can be explained that an increased overlap rate improves diefi ttingness in the region under the electromagnetic forces but weakens forming in the region far from the forces such as the free ends,as shown in Fig.15.In addition,insufficient deformation occurs if there is no overlap,and thus it can be concluded that there is an appropriate overlap rate which is the best for EMIF.From Fig.15,furthermore,inhomogeneous plastic deformation can also be observed,especially along path X in Fig.15(b).The error analysis in Fig.16 also shows that a lower overlap rate is suitable for obtaining good forming quality in EMIF.

    Cui et al.14discussed the effect of the number of discharge regions on final sheet shapes.In their study,three conditions were used,i.e.,the number was set to 3,4,and 5,corresponding to overlap rates of zero,smaller,and larger,respectively.Compared with the desired profile,the forming qualities with the three overlap rates from priority to inferiority were 4,5,and 3,and the forming results for the first two conditions were very close,except in the corner regions.It indicated that there was insufficient deformation at a larger overlap rate(number=5)due to decreased electromagnetic forces when the distance between the coil and the sheet was increased more obviously than that at a smaller overlap rate.In this study from Fig.15,however,the profiles at overlap rates of 20%and 40%are superior to that at 60%,which means a smaller overlap rate is useful to improve the forming quality,which differs from the results from Cui et al.14.That difference may be because the action area in this study is smaller than that in their study.Therefore,considering the forming efficiency,a smaller overlap rate is enough to obtain good forming quality in this study.

    Fig.11 Panel profiles at different rib-height/web-thickness ratios along the two paths.

    Fig.12 Error analysis at different ratios of rib-height to webthickness during EMIF.

    4.Conclusions

    Fig.13 Panel profiles at different discharge paths along the two paths.

    Fig.14 Error analysis at different discharge paths during EMIF.

    In this study,a double-curvature integral panel is formed with EMIF via discharging at four positions.Forming experiments are implemented to study the effects of forming parameters,including discharge voltage,capacitance,ratio of rib-height to web-thickness,discharge path,and overlap rate of the coil,on the die-fittingness,as well as to evaluate the forming quality of the panel workpiece.The obtained results show that an incremental voltage with a lower initial one is helpful to improve the die-fittingness and avoid the collision rebound of the workpiece against the die.Under the condition of the same discharge energy,the capacitance of a relatively lower value can achieve goodfittingness to the die and meanwhile enhance the energy efficiency due to the skin effect.Discharging at the diagonal positions along the ‘X” discharge path achieves the minimum shape deviation and the best forming uniformity of the workpiece,while discharging at the parallel positions with the ‘Z” discharge path gets the worst.A lower overlap rate of the coil is helpful to improve the forming quality.

    Table 5 Different overlap rates along discharge path 1→4→3→2.

    Fig.15 Panel profiles at different overlap rates along the two paths.

    Fig.16 Error analysis at different discharge overlap rates during EMIF.

    Acknowledgements

    This work was supported by the National Science Fund for Distinguished Young Scholars of China(No.51625505),the Key Program Project of the Joint Fund of Astronomy and NationalNaturalScience Foundation ofChina (No.U1537203),and the National Key Basic Research Program of China(No.2011CB012804).The authors would also like to thank the Wuhan National High Magnetic Field Center at Huazhong University of Science and Technology for the help on the experiments.

    Appendix A.Supplementary material

    Supplementary data associated with this article can be found,in the online version,at https://doi.org/10.1016/j.cja.2018.04.005.

    av一本久久久久| 91aial.com中文字幕在线观看| 亚洲美女黄色视频免费看| 亚洲国产av新网站| 精品一区二区三区四区五区乱码 | 国产一区二区在线观看av| 黑人欧美特级aaaaaa片| 水蜜桃什么品种好| 亚洲人成网站在线观看播放| 亚洲一码二码三码区别大吗| 永久免费av网站大全| 亚洲视频免费观看视频| 激情视频va一区二区三区| 一边亲一边摸免费视频| 街头女战士在线观看网站| 国产在线免费精品| 国产男女超爽视频在线观看| 久久婷婷青草| 香蕉精品网在线| 亚洲国产成人一精品久久久| 久久亚洲国产成人精品v| 边亲边吃奶的免费视频| 午夜久久久在线观看| 中文字幕精品免费在线观看视频| 欧美精品高潮呻吟av久久| 久久久久国产网址| 精品一区二区三区四区五区乱码 | 一区二区av电影网| 久久99热这里只频精品6学生| 日本av手机在线免费观看| 性色avwww在线观看| 人妻少妇偷人精品九色| 中文乱码字字幕精品一区二区三区| 国产片内射在线| 日韩av在线免费看完整版不卡| 日本免费在线观看一区| 大香蕉久久网| 黄色配什么色好看| 国产免费现黄频在线看| 亚洲精品一区蜜桃| av线在线观看网站| 日韩一区二区视频免费看| 男人添女人高潮全过程视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲三级黄色毛片| 欧美国产精品va在线观看不卡| 亚洲国产最新在线播放| 欧美日韩国产mv在线观看视频| 亚洲一级一片aⅴ在线观看| 七月丁香在线播放| 99久久精品国产国产毛片| 亚洲av综合色区一区| 老熟女久久久| 国产 精品1| 美女福利国产在线| 国产精品亚洲av一区麻豆 | 捣出白浆h1v1| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91精品国产国语对白视频| 建设人人有责人人尽责人人享有的| 国产精品蜜桃在线观看| 精品国产超薄肉色丝袜足j| 啦啦啦视频在线资源免费观看| 观看美女的网站| 一级毛片电影观看| 色视频在线一区二区三区| 青青草视频在线视频观看| 看免费av毛片| 亚洲国产精品成人久久小说| 日韩视频在线欧美| 国产精品成人在线| 亚洲内射少妇av| 久久久久网色| 伦精品一区二区三区| 国产 一区精品| 久久精品久久精品一区二区三区| 韩国精品一区二区三区| 国产免费福利视频在线观看| av又黄又爽大尺度在线免费看| 卡戴珊不雅视频在线播放| 超碰成人久久| 久久综合国产亚洲精品| 亚洲欧洲日产国产| 免费久久久久久久精品成人欧美视频| 国产亚洲午夜精品一区二区久久| 最近中文字幕高清免费大全6| 激情五月婷婷亚洲| 我的亚洲天堂| 一级黄片播放器| 久久亚洲国产成人精品v| 亚洲美女黄色视频免费看| 久久青草综合色| 亚洲精品国产av蜜桃| 宅男免费午夜| 成人亚洲精品一区在线观看| 女人久久www免费人成看片| 精品人妻一区二区三区麻豆| 国产在线一区二区三区精| 国产深夜福利视频在线观看| 亚洲精品在线美女| 亚洲国产欧美日韩在线播放| 久久久久久人妻| 精品亚洲乱码少妇综合久久| 国产深夜福利视频在线观看| 欧美激情 高清一区二区三区| 中文字幕人妻丝袜一区二区 | 午夜福利影视在线免费观看| 色视频在线一区二区三区| 国产精品一区二区在线不卡| 亚洲国产最新在线播放| 久久国产亚洲av麻豆专区| 男人爽女人下面视频在线观看| 午夜免费观看性视频| 最新的欧美精品一区二区| 大码成人一级视频| 久久97久久精品| 制服人妻中文乱码| 久久久久国产一级毛片高清牌| 色哟哟·www| 成人影院久久| 9191精品国产免费久久| 久久精品国产综合久久久| 免费观看无遮挡的男女| 久久久a久久爽久久v久久| 欧美精品一区二区大全| 久久99热这里只频精品6学生| 国产精品久久久久久久久免| 久久久久精品久久久久真实原创| 狠狠精品人妻久久久久久综合| 婷婷色av中文字幕| 亚洲精品国产av蜜桃| 久久国产精品男人的天堂亚洲| 欧美激情 高清一区二区三区| 成人毛片60女人毛片免费| 桃花免费在线播放| √禁漫天堂资源中文www| 国产精品一区二区在线不卡| 国产高清国产精品国产三级| 丰满乱子伦码专区| 欧美日韩视频高清一区二区三区二| 男男h啪啪无遮挡| 精品少妇一区二区三区视频日本电影 | a级毛片在线看网站| 看非洲黑人一级黄片| 国产欧美日韩综合在线一区二区| 精品国产乱码久久久久久小说| 99久久中文字幕三级久久日本| 97在线视频观看| 一区二区三区四区激情视频| 王馨瑶露胸无遮挡在线观看| 免费在线观看视频国产中文字幕亚洲 | 国产精品国产av在线观看| 国产成人精品无人区| av福利片在线| 99久久中文字幕三级久久日本| 一本—道久久a久久精品蜜桃钙片| 中文字幕色久视频| 美女国产视频在线观看| 成人国产av品久久久| 国产成人精品福利久久| 亚洲综合色惰| 一区二区三区激情视频| 狠狠婷婷综合久久久久久88av| 国产男人的电影天堂91| 欧美国产精品va在线观看不卡| 精品亚洲乱码少妇综合久久| 久久久精品免费免费高清| 成年av动漫网址| 午夜精品国产一区二区电影| 最近的中文字幕免费完整| 日韩一区二区三区影片| 亚洲av综合色区一区| 婷婷成人精品国产| 日韩精品免费视频一区二区三区| 国产精品一区二区在线观看99| 欧美日韩一级在线毛片| 黄色视频在线播放观看不卡| 欧美激情极品国产一区二区三区| 观看av在线不卡| av又黄又爽大尺度在线免费看| 美国免费a级毛片| 黄网站色视频无遮挡免费观看| 青春草国产在线视频| 国产精品久久久久久精品古装| 欧美成人精品欧美一级黄| 一区福利在线观看| 国产av精品麻豆| 日韩伦理黄色片| 黄色毛片三级朝国网站| 国产精品香港三级国产av潘金莲 | 国产日韩欧美亚洲二区| 久久狼人影院| 久久久久国产网址| 欧美激情高清一区二区三区 | 人人妻人人澡人人爽人人夜夜| 熟女电影av网| 69精品国产乱码久久久| 久久影院123| 国产精品99久久99久久久不卡 | 久久精品久久久久久久性| 亚洲av在线观看美女高潮| 欧美激情 高清一区二区三区| 人成视频在线观看免费观看| 在线观看免费日韩欧美大片| 久久99精品国语久久久| 国产男人的电影天堂91| av在线播放精品| 亚洲三区欧美一区| 日韩视频在线欧美| 国产成人精品久久久久久| 亚洲精品久久久久久婷婷小说| 啦啦啦在线免费观看视频4| 女人高潮潮喷娇喘18禁视频| 日本色播在线视频| 久久97久久精品| 国产精品久久久久久久久免| 久久亚洲国产成人精品v| 90打野战视频偷拍视频| 熟女电影av网| 男女啪啪激烈高潮av片| 成人国语在线视频| 亚洲国产最新在线播放| 免费黄频网站在线观看国产| 国产在线免费精品| 人人妻人人澡人人看| 国产免费一区二区三区四区乱码| 亚洲国产精品999| 少妇的丰满在线观看| 国产一级毛片在线| 国产极品天堂在线| 1024视频免费在线观看| 国产精品久久久久成人av| 久久久久精品人妻al黑| 亚洲国产欧美网| 亚洲情色 制服丝袜| av女优亚洲男人天堂| 免费黄频网站在线观看国产| 国产精品国产av在线观看| 国产精品国产三级专区第一集| 成人影院久久| 9191精品国产免费久久| 熟妇人妻不卡中文字幕| 午夜福利在线观看免费完整高清在| 综合色丁香网| 免费看不卡的av| 18禁裸乳无遮挡动漫免费视频| 欧美精品亚洲一区二区| 国产毛片在线视频| 日韩精品有码人妻一区| 午夜福利视频精品| 人妻系列 视频| 日韩电影二区| 成人毛片a级毛片在线播放| 免费播放大片免费观看视频在线观看| 好男人视频免费观看在线| 老熟女久久久| 黄频高清免费视频| 一级片'在线观看视频| 久久久久久久国产电影| 如日韩欧美国产精品一区二区三区| 中文乱码字字幕精品一区二区三区| 久久久久久久久久久免费av| 国产亚洲欧美精品永久| 老司机影院成人| 国产精品久久久av美女十八| 国产av精品麻豆| 国产欧美日韩综合在线一区二区| 欧美国产精品一级二级三级| 亚洲欧美中文字幕日韩二区| 九九爱精品视频在线观看| 叶爱在线成人免费视频播放| 少妇人妻 视频| 午夜福利影视在线免费观看| 美女午夜性视频免费| 亚洲av中文av极速乱| 中文字幕av电影在线播放| 国产乱人偷精品视频| 亚洲一码二码三码区别大吗| 综合色丁香网| 亚洲国产精品999| 日韩中文字幕视频在线看片| 少妇熟女欧美另类| 看免费成人av毛片| 美女主播在线视频| 男女高潮啪啪啪动态图| 在现免费观看毛片| 国产亚洲精品第一综合不卡| 成人亚洲欧美一区二区av| 色吧在线观看| 热99国产精品久久久久久7| 少妇被粗大的猛进出69影院| 亚洲精品自拍成人| 天堂8中文在线网| 另类精品久久| 国产成人精品久久二区二区91 | 久久99热这里只频精品6学生| 久久毛片免费看一区二区三区| av片东京热男人的天堂| 狂野欧美激情性bbbbbb| 免费大片黄手机在线观看| 久久久久精品久久久久真实原创| 黄片播放在线免费| 免费在线观看视频国产中文字幕亚洲 | 久久热在线av| 老汉色∧v一级毛片| 欧美激情 高清一区二区三区| 精品国产乱码久久久久久小说| 超碰成人久久| 中文字幕人妻丝袜制服| 99久久精品国产国产毛片| 美女高潮到喷水免费观看| 国产成人av激情在线播放| 亚洲视频免费观看视频| 中国国产av一级| 久久精品熟女亚洲av麻豆精品| 国产精品国产三级专区第一集| 亚洲婷婷狠狠爱综合网| 精品亚洲乱码少妇综合久久| 免费女性裸体啪啪无遮挡网站| 中文字幕人妻熟女乱码| 天美传媒精品一区二区| 国产一级毛片在线| av有码第一页| 七月丁香在线播放| 97人妻天天添夜夜摸| 国产伦理片在线播放av一区| 街头女战士在线观看网站| 大片免费播放器 马上看| 国产又爽黄色视频| 成人黄色视频免费在线看| 久久精品国产综合久久久| 亚洲少妇的诱惑av| 日韩av不卡免费在线播放| 亚洲综合色网址| 欧美人与性动交α欧美软件| 国产日韩欧美亚洲二区| 满18在线观看网站| 黄色怎么调成土黄色| 亚洲国产成人一精品久久久| 黄色配什么色好看| 2018国产大陆天天弄谢| 九九爱精品视频在线观看| 五月伊人婷婷丁香| 久久久久久人妻| 美女主播在线视频| 久久人人97超碰香蕉20202| 亚洲天堂av无毛| 另类精品久久| 中文字幕人妻熟女乱码| 自线自在国产av| 69精品国产乱码久久久| 国产爽快片一区二区三区| 看非洲黑人一级黄片| 亚洲国产欧美日韩在线播放| 欧美亚洲 丝袜 人妻 在线| 美国免费a级毛片| kizo精华| 亚洲欧洲国产日韩| 丝袜脚勾引网站| 久久久久久人妻| 欧美日本中文国产一区发布| 国产欧美日韩综合在线一区二区| 免费人妻精品一区二区三区视频| 大片免费播放器 马上看| 国产精品国产av在线观看| 免费播放大片免费观看视频在线观看| 国产精品久久久久成人av| 亚洲久久久国产精品| 大码成人一级视频| videos熟女内射| 国产成人91sexporn| 亚洲少妇的诱惑av| av又黄又爽大尺度在线免费看| 国产在线一区二区三区精| 亚洲激情五月婷婷啪啪| 国产成人免费无遮挡视频| 国产日韩欧美在线精品| 国产精品无大码| 国产成人av激情在线播放| 飞空精品影院首页| 久久鲁丝午夜福利片| 亚洲av男天堂| 亚洲五月色婷婷综合| 欧美激情 高清一区二区三区| 国产成人免费观看mmmm| 激情视频va一区二区三区| 成人毛片a级毛片在线播放| 丝袜脚勾引网站| 日韩av不卡免费在线播放| 少妇熟女欧美另类| 一区在线观看完整版| 三上悠亚av全集在线观看| 青草久久国产| 欧美日韩精品网址| 最近中文字幕2019免费版| 久久热在线av| 精品人妻一区二区三区麻豆| 黄片小视频在线播放| 国产精品偷伦视频观看了| 大片免费播放器 马上看| 国产精品香港三级国产av潘金莲 | 精品久久久精品久久久| 99re6热这里在线精品视频| 啦啦啦在线观看免费高清www| 日韩视频在线欧美| 亚洲欧洲精品一区二区精品久久久 | 亚洲少妇的诱惑av| 777米奇影视久久| 国产淫语在线视频| 色哟哟·www| 欧美在线黄色| 亚洲国产av新网站| 亚洲五月色婷婷综合| 国产免费现黄频在线看| 免费av中文字幕在线| 蜜桃国产av成人99| 亚洲欧美一区二区三区黑人 | av在线app专区| videosex国产| 国产亚洲欧美精品永久| 男人舔女人的私密视频| 国产精品无大码| 18禁动态无遮挡网站| 亚洲av成人精品一二三区| 国产在线一区二区三区精| 少妇的丰满在线观看| 久久精品熟女亚洲av麻豆精品| 伦理电影免费视频| 国产精品二区激情视频| 王馨瑶露胸无遮挡在线观看| 国产成人精品无人区| 在线免费观看不下载黄p国产| 两性夫妻黄色片| 国产爽快片一区二区三区| 一级毛片 在线播放| 青春草视频在线免费观看| 高清av免费在线| 久久国产精品男人的天堂亚洲| 性少妇av在线| 飞空精品影院首页| 黄片无遮挡物在线观看| 日产精品乱码卡一卡2卡三| 少妇人妻 视频| 国产亚洲av片在线观看秒播厂| 午夜激情久久久久久久| 亚洲第一青青草原| av国产精品久久久久影院| 久久久久久久久久人人人人人人| 亚洲国产精品成人久久小说| 国产精品国产三级专区第一集| 飞空精品影院首页| 色吧在线观看| 午夜免费男女啪啪视频观看| 高清不卡的av网站| 亚洲欧美精品自产自拍| 日韩一区二区三区影片| 精品国产一区二区久久| 久久精品熟女亚洲av麻豆精品| 制服人妻中文乱码| 亚洲经典国产精华液单| 曰老女人黄片| 亚洲色图 男人天堂 中文字幕| 久久精品国产a三级三级三级| 亚洲国产欧美在线一区| 欧美激情 高清一区二区三区| 国精品久久久久久国模美| 精品国产一区二区三区久久久樱花| a级毛片黄视频| 国产黄频视频在线观看| 国产乱人偷精品视频| 亚洲精品国产一区二区精华液| 国产毛片在线视频| 午夜老司机福利剧场| 精品福利永久在线观看| 国产精品久久久久久久久免| 爱豆传媒免费全集在线观看| 久久亚洲国产成人精品v| 美女福利国产在线| 交换朋友夫妻互换小说| 美女视频免费永久观看网站| 一本色道久久久久久精品综合| 久久毛片免费看一区二区三区| 亚洲 欧美一区二区三区| 欧美激情高清一区二区三区 | 伊人久久国产一区二区| 精品久久久久久电影网| 欧美成人午夜免费资源| 精品国产乱码久久久久久小说| 国产无遮挡羞羞视频在线观看| 最新的欧美精品一区二区| 亚洲欧美一区二区三区国产| 观看av在线不卡| 国产一区有黄有色的免费视频| 亚洲欧洲精品一区二区精品久久久 | 久久精品熟女亚洲av麻豆精品| 综合色丁香网| 免费观看在线日韩| 欧美国产精品va在线观看不卡| 免费观看在线日韩| 免费在线观看视频国产中文字幕亚洲 | 丁香六月天网| 日韩中字成人| 啦啦啦啦在线视频资源| 一级,二级,三级黄色视频| 老司机亚洲免费影院| 久久人妻熟女aⅴ| 色播在线永久视频| 五月开心婷婷网| 久久99热这里只频精品6学生| 伦精品一区二区三区| 欧美日韩国产mv在线观看视频| 一级,二级,三级黄色视频| 18禁观看日本| 国产精品欧美亚洲77777| 国产一区二区三区综合在线观看| 国产精品嫩草影院av在线观看| 香蕉国产在线看| 欧美精品高潮呻吟av久久| 久久综合国产亚洲精品| 国产av码专区亚洲av| 人妻系列 视频| 亚洲av.av天堂| 国产精品不卡视频一区二区| 国语对白做爰xxxⅹ性视频网站| 精品福利永久在线观看| 成年av动漫网址| 夫妻午夜视频| 99re6热这里在线精品视频| 久久国内精品自在自线图片| 人成视频在线观看免费观看| 亚洲精品成人av观看孕妇| 免费人妻精品一区二区三区视频| 日韩一区二区三区影片| 国产黄色视频一区二区在线观看| 日韩人妻精品一区2区三区| 久久久久网色| 最近最新中文字幕大全免费视频 | videosex国产| 国产av国产精品国产| 成年人免费黄色播放视频| 天堂8中文在线网| 少妇被粗大的猛进出69影院| 欧美少妇被猛烈插入视频| 亚洲精品久久成人aⅴ小说| 国产精品三级大全| 久久鲁丝午夜福利片| 亚洲一区中文字幕在线| 又大又黄又爽视频免费| 亚洲人成网站在线观看播放| 9191精品国产免费久久| 中文天堂在线官网| 午夜激情av网站| 免费高清在线观看视频在线观看| 亚洲精品成人av观看孕妇| 电影成人av| 男人爽女人下面视频在线观看| 你懂的网址亚洲精品在线观看| 最近手机中文字幕大全| 90打野战视频偷拍视频| 亚洲婷婷狠狠爱综合网| 久久久久久久久久人人人人人人| 卡戴珊不雅视频在线播放| 性高湖久久久久久久久免费观看| 久久97久久精品| 国产探花极品一区二区| 久久韩国三级中文字幕| 老司机影院毛片| 一区二区三区激情视频| 777米奇影视久久| 国产成人精品福利久久| 在线免费观看不下载黄p国产| 亚洲av.av天堂| 深夜精品福利| 国产精品久久久久久精品电影小说| 只有这里有精品99| 女的被弄到高潮叫床怎么办| 精品国产一区二区久久| 乱人伦中国视频| 日韩大片免费观看网站| 亚洲欧美色中文字幕在线| 亚洲国产av新网站| 久久人人97超碰香蕉20202| 日韩一卡2卡3卡4卡2021年| 久久久久久人妻| 99热网站在线观看| 亚洲av在线观看美女高潮| 一本色道久久久久久精品综合| 午夜影院在线不卡| 搡老乐熟女国产| 国产一区二区三区综合在线观看| a级毛片黄视频| 精品久久久久久电影网| 亚洲,欧美,日韩| 久久久久精品性色| 国产一区二区三区av在线| 新久久久久国产一级毛片| 国产精品免费大片| 尾随美女入室| 日日爽夜夜爽网站| 99国产综合亚洲精品| 亚洲国产欧美网| 亚洲成色77777| 丰满少妇做爰视频| 美女脱内裤让男人舔精品视频| 国语对白做爰xxxⅹ性视频网站| 久久99一区二区三区| 韩国av在线不卡| 欧美日韩精品成人综合77777| 日韩制服丝袜自拍偷拍| 爱豆传媒免费全集在线观看| 另类亚洲欧美激情| 亚洲av电影在线观看一区二区三区| 黄频高清免费视频| 亚洲美女搞黄在线观看| 青青草视频在线视频观看|