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

    汽車磁流變阻尼器多目標(biāo)優(yōu)化分析

    2015-11-26 09:28:42鄧國(guó)紅楊鄂川
    機(jī)床與液壓 2015年24期
    關(guān)鍵詞:工程學(xué)院阻尼器理工大學(xué)

    鄧國(guó)紅,李 飛*,楊鄂川,歐 健,張 勇

    1.重慶理工大學(xué)車輛工程學(xué)院,重慶 400054 2.重慶理工大學(xué)機(jī)械工程學(xué)院,重慶 400054

    1 Introduction

    With the rapid development of automobile industry and the rapid increasing of car ownership,car collision accidents happened under various traffic conditions also show a trend of raising,among them,the low-speed collision is one of the highest frequency of traffic accidents.But automobile parts of energy absorption play an important role during a low-speed collision,so,improving the energy absorption characteristics of components in low-speed collision can effectively improve the passive safety performance of car.Compared with energy absorption mechanism and responsiveness of traditional bumpers,magnetorheological damper can better play its cushioning effect with fast response time,continuous adjustable damping force.However,at present the design of magnetorheological dampers mainly bases on the quasi static models[1-3],which are only applicable to low speed environment,when magnetorheological dampers are applied to high speed impact occasions(such as car crash),it is necessary to adopt a new design theory and method.So,on the basis of Bingham-plastic nonlinear flow modified model(BPM)[4],regarding the maximum damping force and the damping dynamic range as the optimization goals,the NSGA-II algorithm of modeFRONTIER was used to design the structural parameters of the magnetorheological damper for multiobjective optimization,and the optimization results were analyzed.

    2 The working p rincip le of m agnetorheological dam per

    Magnetorheological dampers make mainly use of magnetorheological fluid flowing properties to control the adjustable damping force.Under the action of magnetic field,MRF can transform quickly and reversibly Newtonian fluid with good liquidity into Bingham plastic solid with high yield strength,low liquidity.When the current size of excitation coil changes,applied magnetic field changes gradually,magnetorheological fluid viscosity also changes gradually.Magnetorheological fluid’s shear yield stress changes,as magnetic induction intensity of applied magnetic field changes,so as to adjust the damping force of damper continuously and reversible.According to the different forms of stress and flow of MRF applied to dampers magnetorheological damper mainly consists of valve type,shear,extrusion and shear valve type four basic patterns.According to the background requirements of magnetorheological damper used in the front energy absorption structure in car,based on the single rod cylinder monocular lever shear valve type magnetorheological damper as the research object,the piston structure diagram is shown in Fig.1.

    Fig.1 The piston structure diagram

    Where,nr means cylinder diameter,d1 means piston rod diameter,d2 means piston diameter,la means piston effective length,lc means coil width,dc means the clearance of the coil and damping channel,wc means the depth of coil,dd means damping clearance.

    3 Bingham-p lastic non linear flow m odified model(BPM)

    The mechanical model of the current research of magnetorheological damper have glue Bingham plastic model,Bouc-Wen model,the phenomenon model of Spencer,Herschel-Bulkley model.Now,the most common used constitutive model is Bingham plasticity model[5],referred to as BP model,BP model regards the yielded viscosity of MRF as a fixed value,so the relationship between shear stress and shear strain

    When speed is low ,the flow of MRF can be approximate to laminar flow,the design of the damper adopts BP model can meet the engineering requirements,but as more and more wide engineering application of magnetorheological damper,its running speed has been improved,gradually develops to the impact,collision,drop and other fields of the high speed,at this moment,magnetorheological fluid flow of activation regional of the magnetorheological damper transforms layer flow into turbulent flow[6],BP model applied to low speed condition becomes no longer accurate under high speed,it will produce great error when design,so,in this paper,adopting Bingham-plastic nonlinear flow modified model(BPM)[7-9]based on Bingham plastic model.According to the fluid mechanics knowledge,the maximum damping force can be obtained when magnetorheological fluid flows through the damping channel.rate is the type(1),where,τ0means the shear yield stress;H means the magnetic field strength;˙γmeans the shear strain rate;ηmeans the liquid viscosity of yielded.

    Among them:

    Where,F(xiàn)maxmeans maximum damping force produced by piston structure;Foffmeans the viscous force of piston compression stroke without magnetic field;Fsmeans the damping force of the shear mode without magnetic field;Fτmeans the coulomb damping force with magnetic field;D means the damping dynamic range;Poffmeans pressure loss caused by the viscous damping;Apmeans the piston cross-sectional area;η means MRF viscosity;V1means MRF average flow velocity in work gap;V2means MRF average flow velocity in coil and work cylinder clearance;τymeans MRF shear yield stress with magnetic field;Pηmeans the pressure loss of MRF flows through the piston outside diameter and the inner wall of the work cylinder;Pmlmeans the expansion and compression pressure loss of MRF flows through the coil and the inner wall of the work cylinder;Pcoilmeans the pressure loss of MRF flows through the coil and the inner wall of the work cylinder;Peemeans the expansion and compression pressure loss of MRF flows through the piston outside diameter and working cylinder inner diameter;ρ means the density of magnetorheological fluid;f1and f2mean respectively Darcy friction coefficient related to the Reynolds number of the piston outside diameter and working cylinder inner diameter,the coil and working cylinder inner diameter;Kentryand Kexitmean respectively inlet pressure loss coefficient and the outlet pressure loss coefficient,respectively 0.5 and 1;Kscand Ksemean respectively inlet pressure loss coefficient and outlet pressure loss coefficient,Admeans the sectional area of the piston outside diameter and working cylinder inner diameter;Acmeans the sectional area of the coil and working cylinder inner diameter.

    4 M u lti-ob jective op tim ization of magnetorheological dam per

    4.1 Damper design requirements

    In order to design reasonable magnetorheological damper,it is necessary to estimate the range of control damper damping force and stroke,it is assumed that the vehicle collides with fixed rigid wall at a certain speed,the collision process does not occur secondary collision,all the kinetic energy is dissipated by two magnetorheological dampers,ignoring the friction between the tires and the ground,car is 1 421 kg,the initial velocity of the collision is 3 m/s(10.8 km/h),collision time is approximate56 ms,the required average damping force of magnetorheological damper is 38 kN calculated by the impulse theorem in head-on collision,other dimensions are estimated according to the limiting conditions of damper replaces crash box:working cylinder outside diameter is 60 mm,inner diameter is 54 mm,piston length is 50 mm,the damper’s biggest length is320 mm.The proposed design of this paper is that the magnetorheological damper’s biggest damping force is greater than 38.5 kN,the damping dynamic range is greater than 1.2 when the maximum current is 2 A at 3 m/s.

    Considering the requirements of the magnetorheological damper applied to frontal crash.Setting the maximum speed of the piston as 3 m/s.The adopted magnetorheological fluid is SG-MRF2035 developed by a company,the fitting relationship of the shear yield stressτyand the magnetic flux density B is that:

    The materials of piston and cylinder are industrial pure iron,the piston rod is 45#steel,the wire is enamelled copper withφ0.8 mm,the relative permeability is 1.

    4.2 Optimization design variables,constraints and objective functions

    When the performance parameters of magnetorheological fluid are confirmed,the size of the magnetorheological damper damping channel and the size of the excitation magnetic field will have a decisive effect on the magnetorheological damper performance.The size of excitation magnetic field is related to coil loading current and the size of piston,loading current can be controlled through constant current power,but the piston sizes need to be optimized.So,setting la,lc,dc,wc,dd as optimization variables.Optimization variables and scope is shown Table 1.

    Table 1 Optimization variables and scope

    The initial value for the initial variables are:dd=1 mm,dc=2 mm,wc=5 mm,la=30 mm,lc=20 mm.

    Other parameters values are wr=60 mm,nr=54 mm,d1=24 mm.

    The objective functions:to make Fmaxand D the largest.

    4.3 Establishing optimization model

    This optimization model is established using multiobjective optimization analysis software modeFRONTIER,first simulating the magnetic flux density produced by piston and coil by using magnetic field analysis module of ANSYS and optimization variables[10],inputting the magnetic flux density into the Matlab,and then Matlab brings the optimization variables and the magnetic flux density into the mechanical model(BPM)for magnetorheological damper to calculate the maximum damping force and the damping dynamic range.According to the need to select 30 groups optimization variables,setting optimization variables according to the scope of table 1,setting optimization initial conditions combined with the constraints.Choosing NSGA-II algorithm for the multi-objective optimization,the maximum generation is 200,crossover rate is 0.9,mutation rate is 0.1,multi-objective optimization model is set up as shown in Fig.2.

    Fig.2 The optimization model

    5 Op tim ization resu lts

    Seeing optimized Pareto frontier is very clear(as shown in Fig.3),it’s impossible for two optimization objectives D and Fmaxto achieve the optimal at the same time,the relationship of them is inversely proportional,so only to choose eligible optimized solutions in many frontier solutions.

    Through the correlation coefficient matrix analysis of five optimization variables and two optimization targets(as shown in Fig.4),seeing that the two main factors of influencing the maximum damping force are the clearance of damping channel dd and the piston effective length la,the correlation coefficient are respectively-0.909 and 0.128;The two main factors of influencing the damping dynamic range are the clearance of damping channel dd and the clearance of coil and damping channel dc,the correlation coefficient are respectively 0.630 and-0.433.

    Fig.3 Pareto solutions for the chart

    Fig.4 The correlation coefficient matrix

    When analyzing the relationship of multiple independent variables and dependent variables,it’s more intuitive to observe the effect of each variable on the dependent variables through the response surface analysis.Fig.5 is the response surface analysis of the clearance of damping channel dd and the piston effective length la about the maximum damping force.Fig.6 is the response surface analysis of the clearance of damping channel dd and the clearance of coil and damping channel dc about the damping dynamic range.

    Fig.5 The influence of dd,la to the F max

    Fig.6 The influence of dd,dc to D

    As can be seen from the Fig.5:the change of the maximum damping force along the clearance of damping channel dd direction is faster than the piston effective length la direction,i.e,the influence of the clearance of damping channel on the maximum damping force is greater than the piston effective length,the damping clearance is more smaller and the maximum damping force is more bigger,but the piston effective length makes the maximum damping force decrease within a certain scope.

    From Fig.6 can be obtained:the change of the damping dynamic range along the clearance of damping channel dd direction is faster than the clearance of coil and damping channel dc direction,i.e,the influence of the clearance of damping channel on the damping dynamic range is greater than the clearance of coil and damping channel,the clearance of damping channel is more bigger and the damping dynamic range is more bigger,but the clearance of coil and damping channel is more bigger,the damping dynamic range is more smaller.

    Considering the application background of automotive magnetorheological damper,table 2 lists the five groups typical Pareto solutions which can satisfy the maximum damping force and the damping dynamic range,making the effective length of piston la as big as possible when it’s possible for the two optimization goals to obtain the maximum value,at the same time,according to the design requirements of vehicle under low-speed collision(it is high speed for magnetorheological damper),choosing the fourth group solution as the final optimization solution.

    Simulating the magnetic induction intensity produced by piston coil for before optimization and optimized structure parameters of magnetorheological damper in magnetic field analysis module of ANSYS,as shown in Figs.7,8,the simulation results showed that the optimized magnetic field distribution is more concentrated and reasonable.

    Table 2 Typical Pareto solutions

    Fig.7 Magnetic induction intensity before optimization

    Fig.8 The optimized magnetic induction intensity

    6 Conclusions

    The method of using the non dominated sorting genetic algorithm(NSGA II)to optimize the structure parameters of the magnetorheological damper for multiobjective is feasible.Optimization results showed that it’s impossible for two optimization goals to achieve the optimal at the same time,only to choose the optimized solutions that meet the conditions from all the Pareto solutions,the magnetic field distribution of optimized magnetorheological damper is more concentrated and reasonable.The NSGA-II algorithm can well solve the magnetorheological damper structure parameters optimization problem.

    Applying the optimized magnetorheological damper to the automobile front energy-absorbing structure to study its energy absorption characteristics will be the work that next step need to be done.

    Acknow ledgem ent

    This paper is supported by Graduate Student Innovation fund of Chongqing University of Technology(YCX2014201).

    [1]Zang Chuanxiang.The mechanical properties of magnetorheological damper and its application in artillery Reverse recoil[D].Nanjing:Nanjing University of Science and Technology,2011:33-38.

    [2]Guan Xinchun,Guo Pengfei,Ou Jinping.The multi-objective optimization design and analysis of the magnetorheological damper[J].Engineering Mechanics,2009,26(9):30-35.

    [3]Metered H.Application of Nonparametric Magnetorheological Damper Model in Vehicle Semi-active Suspension System[J].SAE Int.J.Passeng.Cars-Mech.Syst.2012,5(1):715-726.doi:10.4271/2012-01-0977.

    [4]Min Mao,Wei Hu,Norman M W.A nonlinear analtical model for magnetorheological energy absorbers under impact conditions[C]//Proceedings of ASME 2009 Conference on Smart Materials,Adaptive Structures and Intelligent Systems.USA:2009,1447:393-404.

    [5]Shames I H,Cozzarelli F A.Elastic and inelastic stress A-nalysis[M]. Prentice Hall, Englewood Cliffs, New Jersey,1992.

    [6]Yang Guangqiang,SPENCER B F,CARLSON JD,et a1.Large-scale MR fluid dampers:Modeling and dynamic performance considerations[J].Engineering Structures,2002,24(3):309-323.

    [7]Mao Min.Adaptive Magnetorheologocial Sliding Seat System for Ground Vehicles [D]. USA:University of Maryland,2011.

    [8]Gu Xiaoke.Magnetorheological damper multi-objective optimization design and semi-active suspension simulation research[D].Chongqing:Chongqing University,2014:13-17.

    [9]Guan Zhi.The multi-objective optimization design and experimental study of automobile magnetorheological collision buffer[D].Chongqing:Chongqing University,2012.

    [10]Wang Lingling,Song Aiguo,Kuang Yinghui,et al.The design of new magnetorheological damper and magnetic field finite element analysis[J].Machinery Design and Manufacture,2011(4):10-12.

    猜你喜歡
    工程學(xué)院阻尼器理工大學(xué)
    核電廠機(jī)械式阻尼器故障分析及處理
    福建工程學(xué)院
    昆明理工大學(xué)
    福建工程學(xué)院
    幾種軟鋼阻尼器對(duì)建筑安全的影響
    昆明理工大學(xué)
    福建工程學(xué)院
    昆明理工大學(xué)
    浙江理工大學(xué)
    連梁阻尼器的初步設(shè)計(jì)方法
    四川建筑(2018年4期)2018-09-14 00:16:52
    久久亚洲国产成人精品v| 欧美日韩精品网址| 手机成人av网站| 亚洲精品国产av蜜桃| 亚洲精品国产av蜜桃| 国产免费av片在线观看野外av| www.999成人在线观看| 99国产综合亚洲精品| 欧美亚洲 丝袜 人妻 在线| 色婷婷av一区二区三区视频| 亚洲 欧美一区二区三区| 男人爽女人下面视频在线观看| 国产精品免费视频内射| 男女无遮挡免费网站观看| 亚洲精品国产精品久久久不卡| 亚洲av欧美aⅴ国产| 精品国产一区二区久久| 色播在线永久视频| 亚洲伊人色综图| 日韩,欧美,国产一区二区三区| 国产伦理片在线播放av一区| 亚洲国产看品久久| 亚洲 国产 在线| 亚洲少妇的诱惑av| 国产亚洲av片在线观看秒播厂| 欧美日韩av久久| 真人做人爱边吃奶动态| 色综合欧美亚洲国产小说| 精品国产乱码久久久久久小说| 日韩中文字幕视频在线看片| 国产精品一区二区免费欧美 | h视频一区二区三区| 99国产综合亚洲精品| 热99re8久久精品国产| 搡老乐熟女国产| 不卡av一区二区三区| 99国产极品粉嫩在线观看| 黄色视频不卡| 国产欧美日韩一区二区精品| 黄色怎么调成土黄色| 亚洲午夜精品一区,二区,三区| 韩国精品一区二区三区| 91成年电影在线观看| av福利片在线| 美女福利国产在线| 少妇粗大呻吟视频| 日本黄色日本黄色录像| 精品国产乱码久久久久久小说| 欧美精品av麻豆av| 天堂俺去俺来也www色官网| 日韩视频在线欧美| 久久精品国产亚洲av高清一级| 三级毛片av免费| 啦啦啦在线免费观看视频4| 啦啦啦啦在线视频资源| 精品人妻熟女毛片av久久网站| 色精品久久人妻99蜜桃| 国产精品一区二区在线不卡| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲人成电影观看| 青青草视频在线视频观看| 777久久人妻少妇嫩草av网站| 女人高潮潮喷娇喘18禁视频| 99精国产麻豆久久婷婷| 99香蕉大伊视频| av在线播放精品| 亚洲精品中文字幕在线视频| 日韩三级视频一区二区三区| 日韩制服骚丝袜av| 黑人操中国人逼视频| 午夜精品国产一区二区电影| 老汉色av国产亚洲站长工具| 99久久国产精品久久久| 国产亚洲精品一区二区www | 秋霞在线观看毛片| 18禁国产床啪视频网站| 亚洲色图 男人天堂 中文字幕| 久久ye,这里只有精品| 久久久久网色| 国产精品香港三级国产av潘金莲| 最近最新中文字幕大全免费视频| 久久精品久久久久久噜噜老黄| 纵有疾风起免费观看全集完整版| 深夜精品福利| 国产在线一区二区三区精| 丝袜在线中文字幕| 欧美日韩黄片免| 老司机靠b影院| 亚洲精品一区蜜桃| 国产日韩欧美视频二区| 国产精品免费视频内射| 免费在线观看视频国产中文字幕亚洲 | 男女免费视频国产| 美女大奶头黄色视频| 两性夫妻黄色片| 国产男女内射视频| 别揉我奶头~嗯~啊~动态视频 | 91九色精品人成在线观看| 久久精品国产综合久久久| 99国产精品一区二区三区| 777久久人妻少妇嫩草av网站| 成年人午夜在线观看视频| 一本大道久久a久久精品| 一级a爱视频在线免费观看| 婷婷成人精品国产| 女人高潮潮喷娇喘18禁视频| 亚洲精品久久成人aⅴ小说| 高清视频免费观看一区二区| 成年人免费黄色播放视频| 国产一区二区三区综合在线观看| 欧美在线黄色| 99久久人妻综合| 日韩欧美国产一区二区入口| 国产主播在线观看一区二区| 日韩三级视频一区二区三区| 青草久久国产| 日本wwww免费看| 美女扒开内裤让男人捅视频| 91精品国产国语对白视频| 亚洲精品中文字幕在线视频| 精品国产一区二区三区四区第35| 在线观看www视频免费| 老司机在亚洲福利影院| 色婷婷久久久亚洲欧美| 亚洲五月色婷婷综合| 12—13女人毛片做爰片一| 国产精品久久久av美女十八| 欧美午夜高清在线| 黄色片一级片一级黄色片| 一区二区三区乱码不卡18| 涩涩av久久男人的天堂| 亚洲一码二码三码区别大吗| av天堂久久9| 久久久久久免费高清国产稀缺| 欧美激情 高清一区二区三区| 久久久久国内视频| 美女福利国产在线| 欧美黑人精品巨大| 欧美日韩一级在线毛片| 国产老妇伦熟女老妇高清| 日韩电影二区| 亚洲国产欧美一区二区综合| 91大片在线观看| 国产无遮挡羞羞视频在线观看| 久久久久精品人妻al黑| 91麻豆av在线| 欧美 日韩 精品 国产| 一区二区日韩欧美中文字幕| 老司机影院毛片| 少妇裸体淫交视频免费看高清 | 亚洲av国产av综合av卡| 1024视频免费在线观看| 日韩大码丰满熟妇| 18在线观看网站| 一区二区av电影网| av片东京热男人的天堂| 99久久精品国产亚洲精品| 精品久久蜜臀av无| 国产成人啪精品午夜网站| 丰满人妻熟妇乱又伦精品不卡| 国产97色在线日韩免费| 人人妻人人添人人爽欧美一区卜| 欧美精品啪啪一区二区三区 | 国产深夜福利视频在线观看| 国产真人三级小视频在线观看| 国产欧美亚洲国产| 啦啦啦 在线观看视频| 99国产综合亚洲精品| 久久免费观看电影| 日韩一卡2卡3卡4卡2021年| 久久综合国产亚洲精品| 欧美 亚洲 国产 日韩一| 国产主播在线观看一区二区| 中文欧美无线码| 50天的宝宝边吃奶边哭怎么回事| 悠悠久久av| 亚洲性夜色夜夜综合| 欧美中文综合在线视频| 国产在线一区二区三区精| 精品免费久久久久久久清纯 | 一级,二级,三级黄色视频| av天堂久久9| 自线自在国产av| 亚洲七黄色美女视频| 久久午夜综合久久蜜桃| 老熟女久久久| 久久久久久久久久久久大奶| 色播在线永久视频| 18禁国产床啪视频网站| 国产精品秋霞免费鲁丝片| 美国免费a级毛片| 久久亚洲精品不卡| 久久这里只有精品19| 黄色视频在线播放观看不卡| 亚洲精品美女久久av网站| 在线av久久热| 纯流量卡能插随身wifi吗| 国产亚洲精品久久久久5区| 久久ye,这里只有精品| 欧美 日韩 精品 国产| 黄色视频,在线免费观看| 夫妻午夜视频| 狠狠精品人妻久久久久久综合| 俄罗斯特黄特色一大片| 久久久久久久久免费视频了| 午夜免费成人在线视频| 久久九九热精品免费| 精品久久久精品久久久| 欧美老熟妇乱子伦牲交| 亚洲伊人久久精品综合| 欧美成人午夜精品| 欧美xxⅹ黑人| 国产av国产精品国产| 极品少妇高潮喷水抽搐| 在线观看免费日韩欧美大片| 欧美日韩国产mv在线观看视频| 国产成人系列免费观看| 99精品久久久久人妻精品| 色婷婷久久久亚洲欧美| 人人妻人人爽人人添夜夜欢视频| 国产97色在线日韩免费| 国产精品免费视频内射| 十分钟在线观看高清视频www| 亚洲av日韩精品久久久久久密| 欧美日韩成人在线一区二区| 久久久精品94久久精品| 国产精品1区2区在线观看. | 亚洲精品日韩在线中文字幕| 色综合欧美亚洲国产小说| 亚洲国产精品成人久久小说| 国产精品久久久av美女十八| 中文字幕人妻丝袜制服| 正在播放国产对白刺激| 日韩视频在线欧美| 免费在线观看黄色视频的| 亚洲午夜精品一区,二区,三区| 深夜精品福利| 欧美另类一区| 欧美+亚洲+日韩+国产| 交换朋友夫妻互换小说| 精品卡一卡二卡四卡免费| 在线永久观看黄色视频| 爱豆传媒免费全集在线观看| 国产欧美日韩一区二区三区在线| 午夜福利在线观看吧| 91精品伊人久久大香线蕉| 欧美97在线视频| 久久久久久久久久久久大奶| 国产精品久久久久成人av| 国产欧美日韩精品亚洲av| 国产区一区二久久| 国产高清视频在线播放一区 | 女性被躁到高潮视频| 欧美日本中文国产一区发布| 制服人妻中文乱码| 亚洲精品久久成人aⅴ小说| 久久久精品国产亚洲av高清涩受| 日韩中文字幕欧美一区二区| 伊人亚洲综合成人网| 99久久99久久久精品蜜桃| 欧美老熟妇乱子伦牲交| 亚洲成国产人片在线观看| 欧美亚洲日本最大视频资源| 汤姆久久久久久久影院中文字幕| 在线观看舔阴道视频| 在线精品无人区一区二区三| www.熟女人妻精品国产| 久久久国产欧美日韩av| 国产成人一区二区三区免费视频网站| 国产亚洲一区二区精品| 91av网站免费观看| 欧美日韩国产mv在线观看视频| 国产91精品成人一区二区三区 | 国产成人精品久久二区二区免费| 午夜影院在线不卡| 色精品久久人妻99蜜桃| 黑人猛操日本美女一级片| 欧美精品啪啪一区二区三区 | 美女午夜性视频免费| 日韩三级视频一区二区三区| 后天国语完整版免费观看| 国产精品一区二区免费欧美 | 国产精品自产拍在线观看55亚洲 | 免费久久久久久久精品成人欧美视频| av天堂久久9| av不卡在线播放| 精品国产一区二区三区四区第35| 国产免费福利视频在线观看| 免费在线观看黄色视频的| 国产精品一区二区在线不卡| 在线十欧美十亚洲十日本专区| 一个人免费看片子| 男女免费视频国产| 丰满人妻熟妇乱又伦精品不卡| 九色亚洲精品在线播放| 欧美日本中文国产一区发布| 老司机午夜十八禁免费视频| 欧美精品啪啪一区二区三区 | 亚洲精品粉嫩美女一区| 亚洲国产av影院在线观看| 一区二区三区乱码不卡18| 男人舔女人的私密视频| 啦啦啦 在线观看视频| 国产黄频视频在线观看| 欧美另类亚洲清纯唯美| 久久国产精品影院| av欧美777| 精品福利永久在线观看| 欧美精品av麻豆av| 亚洲精品中文字幕在线视频| 久久精品人人爽人人爽视色| 自拍欧美九色日韩亚洲蝌蚪91| 精品福利观看| 免费在线观看完整版高清| 人人妻人人澡人人看| 亚洲欧美清纯卡通| 午夜免费观看性视频| 亚洲激情五月婷婷啪啪| 淫妇啪啪啪对白视频 | 他把我摸到了高潮在线观看 | 国产欧美日韩一区二区精品| 欧美人与性动交α欧美软件| 午夜福利视频在线观看免费| 成人黄色视频免费在线看| 欧美国产精品va在线观看不卡| 久久中文看片网| 国产一区二区在线观看av| 欧美日韩黄片免| 久久久久国内视频| 国产精品国产三级国产专区5o| 日本a在线网址| 美国免费a级毛片| 人人妻人人澡人人爽人人夜夜| 亚洲av电影在线观看一区二区三区| 久久国产精品男人的天堂亚洲| 欧美日韩一级在线毛片| 午夜精品国产一区二区电影| 国产国语露脸激情在线看| 欧美亚洲日本最大视频资源| 欧美人与性动交α欧美精品济南到| 水蜜桃什么品种好| 99国产精品一区二区蜜桃av | 日韩欧美一区视频在线观看| 亚洲av男天堂| 日韩制服丝袜自拍偷拍| 最黄视频免费看| 丰满饥渴人妻一区二区三| 电影成人av| 亚洲一区中文字幕在线| 美女大奶头黄色视频| 色老头精品视频在线观看| 亚洲精品中文字幕一二三四区 | 国产欧美日韩一区二区三区在线| 欧美亚洲 丝袜 人妻 在线| netflix在线观看网站| 男人添女人高潮全过程视频| 天堂中文最新版在线下载| 成人国产av品久久久| 精品国产一区二区久久| 老司机深夜福利视频在线观看 | 国产黄色免费在线视频| 黑人猛操日本美女一级片| 在线观看www视频免费| 久久国产精品大桥未久av| 免费久久久久久久精品成人欧美视频| 欧美激情久久久久久爽电影 | 精品国产乱子伦一区二区三区 | 久久九九热精品免费| av电影中文网址| 另类亚洲欧美激情| 久久精品久久久久久噜噜老黄| 欧美 亚洲 国产 日韩一| 男女床上黄色一级片免费看| 国产成人免费观看mmmm| 国产精品一区二区精品视频观看| 久久久精品免费免费高清| 欧美日韩福利视频一区二区| 手机成人av网站| 老汉色av国产亚洲站长工具| 欧美中文综合在线视频| 一个人免费看片子| 久9热在线精品视频| 老熟妇乱子伦视频在线观看 | 中文字幕av电影在线播放| 亚洲精品美女久久久久99蜜臀| 亚洲av电影在线观看一区二区三区| 一区二区日韩欧美中文字幕| 搡老乐熟女国产| 人妻 亚洲 视频| 亚洲avbb在线观看| 久久毛片免费看一区二区三区| 大片电影免费在线观看免费| 男女无遮挡免费网站观看| 久久久久国产一级毛片高清牌| 精品亚洲乱码少妇综合久久| 精品一区二区三区四区五区乱码| 伊人久久大香线蕉亚洲五| 久久久国产成人免费| 老熟妇仑乱视频hdxx| 国产欧美日韩综合在线一区二区| 午夜精品国产一区二区电影| 色婷婷av一区二区三区视频| 巨乳人妻的诱惑在线观看| av又黄又爽大尺度在线免费看| 亚洲少妇的诱惑av| 成人三级做爰电影| 1024香蕉在线观看| 亚洲国产精品一区三区| www.av在线官网国产| 啦啦啦中文免费视频观看日本| 亚洲精品一二三| 777米奇影视久久| 午夜影院在线不卡| 欧美另类一区| 亚洲第一青青草原| 国产成人免费无遮挡视频| 亚洲精品国产精品久久久不卡| 一级毛片女人18水好多| 国产三级黄色录像| 真人做人爱边吃奶动态| 欧美黄色片欧美黄色片| 可以免费在线观看a视频的电影网站| 男人舔女人的私密视频| 免费看十八禁软件| 热re99久久国产66热| 美女高潮喷水抽搐中文字幕| 国产精品一区二区精品视频观看| 黄色 视频免费看| 亚洲精品粉嫩美女一区| e午夜精品久久久久久久| 国产老妇伦熟女老妇高清| 午夜免费观看性视频| 99精品久久久久人妻精品| 少妇 在线观看| 精品乱码久久久久久99久播| 黄色片一级片一级黄色片| 亚洲一区中文字幕在线| www.999成人在线观看| e午夜精品久久久久久久| 丰满饥渴人妻一区二区三| 亚洲国产欧美一区二区综合| 免费在线观看视频国产中文字幕亚洲 | av有码第一页| √禁漫天堂资源中文www| 亚洲av成人不卡在线观看播放网 | 老司机影院成人| 欧美激情高清一区二区三区| 黑人操中国人逼视频| 国产精品自产拍在线观看55亚洲 | 免费高清在线观看日韩| 亚洲av电影在线观看一区二区三区| 亚洲 国产 在线| 啪啪无遮挡十八禁网站| 亚洲av成人一区二区三| 国产真人三级小视频在线观看| 欧美精品一区二区免费开放| 下体分泌物呈黄色| 亚洲av美国av| 成人国产av品久久久| 黄色怎么调成土黄色| 999久久久精品免费观看国产| 中文字幕人妻丝袜制服| 91精品国产国语对白视频| 国产主播在线观看一区二区| 日韩一卡2卡3卡4卡2021年| 黄色 视频免费看| 亚洲精品自拍成人| 亚洲国产精品一区三区| avwww免费| 老司机午夜十八禁免费视频| 亚洲精品在线美女| 国产免费现黄频在线看| 精品国产乱子伦一区二区三区 | 高清av免费在线| 伊人久久大香线蕉亚洲五| 亚洲美女黄色视频免费看| 亚洲精品一卡2卡三卡4卡5卡 | 日韩欧美国产一区二区入口| 国产成人影院久久av| 妹子高潮喷水视频| 操美女的视频在线观看| 两个人免费观看高清视频| 精品国产乱子伦一区二区三区 | 亚洲午夜精品一区,二区,三区| 中亚洲国语对白在线视频| 蜜桃国产av成人99| 亚洲免费av在线视频| 久热这里只有精品99| 国产日韩一区二区三区精品不卡| 国产在线视频一区二区| 超碰成人久久| 亚洲美女黄色视频免费看| av天堂在线播放| 免费久久久久久久精品成人欧美视频| 精品一品国产午夜福利视频| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品在线电影| 亚洲精品国产av成人精品| 深夜精品福利| 青春草亚洲视频在线观看| 美女中出高潮动态图| 久久久国产一区二区| 国产淫语在线视频| 老司机福利观看| 久久精品aⅴ一区二区三区四区| 少妇的丰满在线观看| 999精品在线视频| 日韩大码丰满熟妇| 国产亚洲av高清不卡| 亚洲精品日韩在线中文字幕| 在线十欧美十亚洲十日本专区| 99精国产麻豆久久婷婷| 亚洲中文日韩欧美视频| 国产亚洲欧美在线一区二区| 美女福利国产在线| 亚洲av电影在线进入| 国产成人a∨麻豆精品| 国产免费现黄频在线看| 色播在线永久视频| 天天添夜夜摸| 一本久久精品| 亚洲国产日韩一区二区| 黑人欧美特级aaaaaa片| 久久av网站| 我的亚洲天堂| 男女床上黄色一级片免费看| 欧美乱码精品一区二区三区| 亚洲精品自拍成人| 中文字幕制服av| 国产成人免费观看mmmm| 国产av一区二区精品久久| 91精品国产国语对白视频| 亚洲黑人精品在线| 天天躁夜夜躁狠狠躁躁| 国产精品 欧美亚洲| 成人亚洲精品一区在线观看| 一级a爱视频在线免费观看| 精品熟女少妇八av免费久了| 欧美97在线视频| 国产精品一区二区免费欧美 | 水蜜桃什么品种好| 最黄视频免费看| 日韩精品免费视频一区二区三区| 色视频在线一区二区三区| 国产免费现黄频在线看| 少妇 在线观看| 国产极品粉嫩免费观看在线| 久热这里只有精品99| 窝窝影院91人妻| 在线看a的网站| 黄网站色视频无遮挡免费观看| 日本wwww免费看| 欧美黑人精品巨大| 91精品国产国语对白视频| 亚洲国产欧美一区二区综合| 99久久人妻综合| 夫妻午夜视频| 国产精品.久久久| 老司机亚洲免费影院| 国产熟女午夜一区二区三区| 日韩电影二区| 精品国内亚洲2022精品成人 | 在线观看人妻少妇| 91大片在线观看| 香蕉丝袜av| 蜜桃国产av成人99| 美女扒开内裤让男人捅视频| 成人亚洲精品一区在线观看| 日韩一区二区三区影片| 少妇人妻久久综合中文| 国产欧美日韩一区二区三 | 国产又色又爽无遮挡免| 真人做人爱边吃奶动态| 1024视频免费在线观看| 丁香六月天网| xxxhd国产人妻xxx| 一级片免费观看大全| 97在线人人人人妻| 亚洲欧洲日产国产| 在线天堂中文资源库| av免费在线观看网站| 国产一区二区在线观看av| 亚洲,欧美精品.| 老司机亚洲免费影院| 18禁黄网站禁片午夜丰满| 久热这里只有精品99| 美女福利国产在线| 十八禁人妻一区二区| 搡老岳熟女国产| 国产免费视频播放在线视频| 亚洲国产欧美一区二区综合| 精品亚洲成a人片在线观看| 久久热在线av| 亚洲欧美一区二区三区黑人| 亚洲全国av大片| 成人影院久久| 免费观看人在逋| 亚洲国产欧美日韩在线播放| 日本猛色少妇xxxxx猛交久久| 五月开心婷婷网| 日韩 亚洲 欧美在线| 最新在线观看一区二区三区| 人人妻人人爽人人添夜夜欢视频| 老司机深夜福利视频在线观看 | 波多野结衣av一区二区av| 亚洲色图综合在线观看| 伊人亚洲综合成人网| videos熟女内射| 青青草视频在线视频观看| 亚洲七黄色美女视频| 日韩制服丝袜自拍偷拍| 久久久久久久久免费视频了| 永久免费av网站大全| 免费在线观看影片大全网站| 男女国产视频网站|