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

    A Stroke-Limitation AMD Control System with Variable Gain and Limited Area for High-Rise Buildings

    2024-02-19 12:03:06ZuoHuaLiQingGuiWuJunTengandChaoJunChen

    Zuo-Hua Li,Qing-Gui Wu,★,Jun Teng,★and Chao-Jun Chen,2

    1School of Civil and Environmental Engineering,Harbin Institute of Technology(Shenzhen),Shenzhen,518055,China

    2Department of Civil and Environmental Engineering,The Hong Kong Polytechnic University,Hong Kong,100872,China

    ABSTRACT

    Collisions between a moving mass and an anti-collision device increase structural responses and threaten structural safety.An active mass damper(AMD)with stroke limitations is often used to avoid collisions.However,a strokelimited AMD control system with a fixed limited area shortens the available AMD stroke and leads to significant control power.To solve this problem,the design approach with variable gain and limited area(VGLA)is proposed in this study.First,the boundary of variable-limited areas is calculated based on the real-time status of the moving mass.The variable gain(VG)expression at the variable limited area is deduced by considering the saturation of AMD stroke.Then,numerical simulations of a stroke-limited AMD control system with VGLA are conducted on a high-rise building structure.These numerical simulations show that the proposed approach has superior strokelimitation performance compared with a stroke-limited AMD control system with a fixed limited area.Finally,the proposed approach is validated through experiments on a four-story steel frame.

    KEYWORDS

    High-rise buildings;active control;stroke limitations;variable gain;variable limited area

    1 Introduction

    Tuned mass dampers(TMDs)[1–9],semi-actively tuned mass dampers(SATMDs)[10–15],active tuned mass dampers (ATMDs),and active mass dampers (AMDs) are often used for engineering vibration control[16–22].In theory,the control effect of AMD is better than that of the other control systems.However,AMD’s moving mass can experience high speeds during operations,resulting in excessive strokes.In particular,when the moving mass moves rapidly and runs close to the maximum permitted stroke,it may collide with the anti-collision device,thus increasing the structural response and even endangering the safety of the AMD control system and structure.Therefore,it is important to limit moving mass stroke and improve the reliability of AMD control systems.

    A vibration control system should exhibit robust performance.For continuous systems,a robust H∞control method based on linear matrix inequalities(LIM)was proposed[23].By combining H∞and H2control,the active control system of a single-story frame structure is established by considering the interaction between the irregularly controlled structure and foundation.The simulation results illustrate that the control system exhibits robust performance [24].A pole assignment control algorithm is applied to the vibration control of high-rise structures to ensure the robust performance of the control system [25,26].Because the gain can determine the control effects and parameters of a control system,some scholars have proposed the saturating control and variable gain(VG)method to improve the control performance effectively.For example,applying a saturating control algorithm based on linear quadratics to an AMD control system can limit the moving mass stroke [27].To improve the control performance of the control system for high-rise buildings,an optimal design considering the saturation of actuators and stroke limitations was proposed [28].A control system with a VG was established using linear control theory,and the optimal gain was selected by solving inequality problems to achieve the optimal control force [29].An active vibration control method with variable gain(VG)control was proposed for high-rise buildings under wind load or earthquake action[30].Variable gain-based structural control was used to control civil engineering structures in the case of earthquakes[31].The AMD control method of VG was proposed,and the gain matrix was calculated using a linear quadratic regulator (LQR) control algorithm.The relative optimal weight matrix within the specified range was selected based on the magnitude of external excitation to obtain the optimal feedback gain matrix [32].An algorithm of VG for the stroke-limitation control system was designed,and the design was verified by conducting an experiment [33].The algorithm of VG with a fixed and limited area requires a large external energy supply.Its limit mechanism still needs improvement to enhance the stroke-limited performance and further avoid the collisions of AMD,ensuring the safe operation of AMD.

    To reduce the influence of stroke limitations on the performance of AMD control systems,a stroke-limited AMD control system with variable gain and limited area (VGLA) was proposed in this study.The running states of the VGLA control system are divided into unlimited and variablelimited areas,and the control gain is calculated using an algorithm with H∞/H2and pole assignment constraints.The deceleration mechanism of the moving mass is determined based on the control force and maximum stroke.Then,based on the real-time states of the moving mass,the gain expression in the variable limited area was deduced considering the saturation of the AMD stroke.Furthermore,numerical simulations were conducted using a stroke-limited AMD control system.Finally,the effectiveness of the proposed methodology was verified by conducting an experiment with a four-story steel frame structure.

    2 The Stroke-Limitation Principle of AMD with VGLA

    The main objective of structural control is to reduce the acceleration response of buildings to a comfortable level.An AMD device is commonly installed on the structure,as shown in Fig.1.Under control force u,the building structure vibrations caused by external forces W can be significantly reduced.

    The running states of AMD are divided into eight states,as shown in Fig.2;the AMD can run freely during running states 1,4,5,and 8.When the moving mass is in a limited area and moves toward the equilibrium position(running states 3 and 7,as illustrated in Fig.2),there is no need to limit the stroke.However,the moving mass may collide with the anti-collision device during running states 2 and 6,thus reducing the control performance and even endangering the safety of AMD.Therefore,to ensure AMD’s safety and high control performance,it is critical to optimize the stroke-limitation strategy during running states 2 and 6.

    Figure 1 :Building structure equipped with AMD

    Figure 2 :The AMD operation diagram

    The force equilibrium of an AMD model(Fig.1)with uncertain mass and stiffness is

    whereM,C,andKrepresent the nominal mass,damping,and stiffness matrices of the system,respectively;ΔMandΔKrepresent the uncertain quantities of the mass and stiffness matrices,respectively;X(t) represents the displacement vector;BsandBwrepresent the position matrices of the control forces and excitations,respectively;andu(t)andw(t)represent the control force vector and external excitation vector,respectively.

    MatricesΔK,ΔM,Bs,andBware described as follows:

    MatricesM,C,K,andX(t)are described as follows:

    Assuming that the state vectors areZ(t)=X(t)(t)T,then the AMD control system with uncertain parameters,such as quality and stiffness,can be expressed as follows:

    MatricesC1,C2,D11,D12,D21,andD22are described as follows:

    Based on reference[34],the AMD control force is obtained using the following equation:

    whereui(t)denotes the component of the force andGdenotes the feedback gain of the AMD.

    It is assumed that the matrix of uncertain parameters of structural parameters is

    Substituting Eqs.(6)and(7)into Eq.(4),the AMD control system is expressed as follows:

    Based on reference[35],the H2control law with state feedback can be expressed as follows:

    Based on reference[36],the H∞control law with state feedback can be expressed as follows:

    Based on reference [37],the performance requirements for regional pole placement can be expressed as follows:

    This problem can be solved by finding a common Lyapunov matrix,denoted as X,that satisfies a set of linear matrix inequalities(9)–(11):

    To solve this multi-objective control problem,the following convex optimization problem is required:

    whereα,β,γ0andη0denote the given parameters.

    The state feedback controller design method has H2/H∞performance requirements and closedloop region pole constraints,which are solved using LMI techniques.LetW′andX′denote the optimal solutions to the above problem.Then,the state feedback control law of the system is

    The gain of the AMD control system is expressed as follows:

    wherega,gn,andgn+1denote the gains based on the displacement of the floor in which the AMD is installed,the top floor,and AMD,respectively.

    When AMD runs into a variable-limited area,as shown in Fig.3,the displacement of the floor with AMD and moving mass must be calculated in real time to determine the gains and control forces.Because the building structure has many degrees of freedom,the displacement of the other floors is not considered when calculating the control force,which reduces the computation.Then,the control forces are expressed as follows:

    wherezadenotes the displacement of the floor where the AMD is installed andzn+1denotes the displacement of the AMD.

    Figure 3 :The curve of the stroke-limitation mechanism

    Assuming thatga=-gn+1,then

    wherex(t)is the stroke of the moving mass andx(t)>0.

    Based on Eq.(17),we have

    From Eq.(17),ifgn+1>0,the control force at the variable-limited area has the opposite sign tox(t).When the AMD runs into a variable limited-area,the moving mass is subjected to braking forceu′(t)reversing tox(t),sogn+1>0.

    Based on Eqs.(15)–(18),the gain of the AMD control system is expressed as follows:

    When the AMD runs into the variable limited area (the shadow area shown in Fig.3),the gain at the variable-limited area is calculated by considering the maximum allowable stroke L and the realtime velocity of moving massVd.If AMD is not in a variable-limited area,the gain of the AMD control system is calculated using the regional pole-assignment algorithm.

    In the variable-limited area,the velocity of the moving mass decreases.When the stroke of the moving mass reachesL,its velocity is zero(Fig.3).The moving mass’s velocity in a variable limited area is

    wherex(t),L,andFmaxdenote the stroke of the moving mass,the maximum allowable stroke,and the control force of the AMD control system,respectively.Vddenotes the real-time velocity of the moving mass at the boundaries of a variable-limited area.The accelerations can be determined by differentiating the velocity of the moving mass as follows:

    The expected control forces are

    Substituting Eq.(19)–(22)into Eq.(18),the gain is based on the displacement of the floor where AMD is installed and the moving mass is

    3 Numerical Simulation of an AMD Control System Using VGLA

    3.1 Modeling of the AMD Control System Using VGLA

    The high-rise building named KingKey Financial Center in Shenzhen,shown in Fig.4a,has 100 floors above the ground and four floors underground.The height above the ground is 441.8 m,and their structural modal parameters are listed in Table 1.The AMD control system depicted in Figs.4b and 4c includes two AMD control systems.The AMD’s control parameters are listed in Table 2.When the moving mass is in an unlimited area.The control forces can be represented by Eq.(14).The VGLA control system with stroke limitations is illustrated using a Simulink block diagram(Fig.5).

    Table 1 : The periods and frequencies of KK100 along the minor axis

    Table 2 : Key parameters of the AMD control system

    Figure 4 : The high-rise building and AMD control system.(a) The high-rise building; (b) The installation position of AMD;(c)AMD control system

    The Simulink block diagram shown in Fig.5 can be used to depict a VGLA system with stroke limitations.The 10-year return period fluctuating wind speed is based on the Davenport spectrum,and a mixed auto-regressive moving average model is used to simulate its stochastic process[38].The fluctuating wind load of theith floor is calculated using the following equation[33]:

    Figure 5 :Simulink module of the stroke-limitation control system of VGLA

    To verify and compare the control performance of the VGLA,the method of stroke-limited AMD control system with VG in reference[33]is analyzed.The VG velocities of the moving mass in a fixed limited area are calculated using Eq.(25)(Fig.6).The VG gain of the moving mass in a fixed limited area is calculated using Eq.(26).Its fixed boundary of fixed limited area,dis 1.5 m.

    Figure 6 :The curve of the stroke-limitation mechanism of the VG control system

    3.2 Simulation Results for VG and VGLA

    Fig.7 and Table 3 show that AMD without stroke limitations can effectively suppress the dynamic response of the structure.The VG,VGLA,and AMD control systems without stroke limitations have the same displacement and acceleration control effects.When the stroke of the moving mass exceeds 1.5 m,the control force and power of the VG increase rapidly(Figs.7e,7f,7i,and 7j).In particular,as demonstrated in Fig.7i,when the moving mass of the VG control system enters a fixed limited area,it moves rapidly.It is subjected to a large braking force,resulting in a sudden increase in power to 1459 kW.

    Table 3 : Control effects and parameters of VG,VGLA,and AMD control systems without stroke limitations

    Figure 7 : (Continued)

    Figure 7 :The top floor’s responses and the AMD parameters

    Furthermore,the peak stroke of the VG is 1.73 m,and the peak stroke of VGLA is 1.65 m.The peak stroke of VG is 5.98% lower than that of the AMD system without stroke limitations.The peak stroke of VGLA is 10.33% lower than that of AMD without stroke limitation.However,the peak control force and power of the VGLA are 12.6% and 37.27% lower than those of the VG,respectively.The mean-square control force and power of the VGLA are 3.20% and 16.47% lower than those of the VG,respectively.Therefore,the control performance of VGLA is significantly better than that of VG.

    4 Experimental Verification

    As shown in Fig.8,the experimental system comprises a single-span four-story steel frame structure and its control system.The AMD control system mainly comprises a servo driver,control motor,DSPACE of board model DS1103,Ethernet for Control Automation Technology(EtherCAT)bus,and microcomputer.The system is equipped with an EtherCAT bus system composed of an LMCF210201 linear control motor produced by Shenzhen Dazu Laser Company,Guangdong Province,EL3008 input terminal,EL4034 output terminal,and EK1100 coupling terminal; an AX5000 servo driver produced by Beckhoff Company,Germany.

    Figure 8 : (Continued)

    Figure 8 :Picture of the AMD experimental system and its schematic diagram

    The weak axial displacement and acceleration of the four-floor steel frame structure can be measured by applying the measurement system.The force balance acceleration sensor is used in the system,and the micro-epsilon series laser displacement sensor is used in the system.The structural mass matrix can be obtained using the concentrated mass method.That is,the mass of each floor of the structure is obtained by calculating the geometric size,density,and other parameters of the controlled structure as well as the related parameters of the equipment; the mass matrix is finally combined.The structural stiffness matrix is obtained by utilizing FEM software.The aforementioned components of the AMD control system are shown in Fig.8.

    Table 4 shows that the first two modes of the structure mainly contribute to the dynamic response of the structure because the mass participation coefficient of the first two modes is 0.946,which is approximately 1.The loading frequency of the loading system is 1 Hz.The corresponding peak value of the excitation force is 45.89 N.Under this excitation,the maximum stroke is 0.4 m.Under the action of excitation load,the time-history comparison of the measured displacement and acceleration responses of each floor of the structure is shown in Fig.9,and the control effects are shown in Table 5.

    Table 4 : Mode mass participation coefficient and natural frequency

    Table 5 : Control effects and AMD parameters of the test

    Figure 9 : (Continued)

    Figure 9 :Experimental responses and AMD parameters

    Table 5 and Fig.9 show that the structural response does not completely obey the sinusoidal variation law because of the interaction between the AMD control system and the structure and the coupling of the horizontal and vertical vibrations of the structure.The control effects of VG and VGLA are higher than that of an AMD without stroke limitations.The reason for this phenomenon is that the stroke of AMD without stroke limitation exceeds the permitted stroke of the control system,and the collision with the anti-collision device reduces its control effects.

    The peak stroke of VGLA is 19.72% less than that of VG,but the peak control force and power of VGLA are 23.05% and 24.11% lower than those of VG,respectively.The mean-square stroke of VGLA is 3.15% less than that of VG.However,the mean-square control power of VGLA is 10.40% lower than that of VG.The experimental results agree with the numerical simulation results(Section 3).Therefore,the control performance of VGLA is better than that of the traditional stroke-limitation control method VG.The proposed VGLA method exhibits a good control effect and can significantly decrease the stroke to guarantee the safe operation of the AMD control system.

    5 Conclusions

    In order to limit the stroke of the moving mass in the AMD control system,the relation between the stroke and the corresponding feedback gain was established,and the design method of a strokelimitation system with variable limited area and gain was proposed.Finally,the numerical model of an actual building structure and the experimental model of the four-story frame were taken to verify the necessity and effectiveness of the stroke-limitation system with VGLA.The main conclusions are drawn as follows:

    (1) The stroke-limited AMD control system of the VGLA exhibits a good control effect.In a situation where no collision of the control system without stroke limitation occurred,the control effect of the VGLA is approximately equal to that of the control system without stroke limitation.In a situation where a collision of the control system without stroke limitation occurred,the control effects of the VGLA are better than those of the control system without stroke limitation.

    (2)The maximum force and stroke of the AMD are considered in the stroke-limitation principle,which can limit the control force and power of the VGLA to permissible ranges.

    (3) Considering the variable limited area,the stroke limit effect of VGLA is better than that of VG;however,the control force and power of VGLA are lower than those of VG.

    (4) The numerical simulation results agree with the experimental results.The stroke-limitation system with VGLA proposed in this study can effectively limit the stroke to a reasonable range,thus ensuring the control effects and safety of the AMD system.

    Based on the variable gain feedback control method,the VGLA method proposed in this paper achieves a stroke limitation control method for building structures with variable-limit areas.The application results show that the proposed method has a short stroke and requires low control power,indicating its potential for effectively controlling structural vibrations in high-rise buildings.

    Acknowledgement: The authors would like to thank Houbing Xing for advice on experimental design.

    Funding Statement: This research was founded by the Funds for Creative Research Groups of National Natural Science Foundation of China(Grant No.51921006),the National Natural Science Foundations of China (Grant No.51978224),the National Major Scientific Research Instrument Development Program of China (Grant No.51827811),the National Natural Science Foundation of China,(Grant No.52008141) and the Shenzhen Technology Innovation Program (Grant Nos.JCYJ20170811160003571,JCYJ20180508152238111 and JCYJ20200109112803851).

    Author Contributions:Zuohua Li and Qinggui Wu contributed equally to this work.Qinggui Wu wrote the paper and summarized the results.Zuohua Li and Jun Teng participated in the data analysis and conceived the study.Chaojun Chen reviewed the study plan.

    Availability of Data and Materials:All data generated and analyzed during this study are included in this published article.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    欧美日韩黄片免| 亚洲国产欧美人成| 在线观看美女被高潮喷水网站 | 最近在线观看免费完整版| 日本在线视频免费播放| 欧美激情久久久久久爽电影| 最好的美女福利视频网| 最近视频中文字幕2019在线8| 午夜福利视频1000在线观看| 成人18禁在线播放| 日本成人三级电影网站| 黄色女人牲交| 女人十人毛片免费观看3o分钟| 久9热在线精品视频| 精品熟女少妇八av免费久了| 精品久久久久久久人妻蜜臀av| 色哟哟哟哟哟哟| 日本 av在线| 69人妻影院| 国产视频内射| 国产精品亚洲美女久久久| 此物有八面人人有两片| 国产黄a三级三级三级人| 嫁个100分男人电影在线观看| 国内久久婷婷六月综合欲色啪| 在线十欧美十亚洲十日本专区| 国产爱豆传媒在线观看| 欧美性猛交黑人性爽| 国产欧美日韩一区二区精品| 国产激情偷乱视频一区二区| 久久精品人妻少妇| 黄色视频,在线免费观看| 狂野欧美白嫩少妇大欣赏| 怎么达到女性高潮| 久久精品人妻少妇| 我要搜黄色片| 欧美极品一区二区三区四区| 精品福利观看| 91麻豆精品激情在线观看国产| 天天躁日日操中文字幕| 欧美中文综合在线视频| 精品福利观看| 美女高潮喷水抽搐中文字幕| 欧美日韩福利视频一区二区| 欧美高清成人免费视频www| 伊人久久大香线蕉亚洲五| 19禁男女啪啪无遮挡网站| 久久午夜亚洲精品久久| 啦啦啦观看免费观看视频高清| 国产午夜福利久久久久久| 国产精品,欧美在线| 免费看十八禁软件| svipshipincom国产片| 老司机在亚洲福利影院| 国产免费av片在线观看野外av| 色精品久久人妻99蜜桃| 精品久久久久久久久久免费视频| 国产午夜福利久久久久久| 深夜精品福利| 在线免费观看不下载黄p国产 | 黑人欧美特级aaaaaa片| 久久精品91蜜桃| 18禁国产床啪视频网站| 国产精品嫩草影院av在线观看 | 亚洲国产色片| 国产亚洲精品综合一区在线观看| av天堂在线播放| 舔av片在线| 欧美色欧美亚洲另类二区| 免费av观看视频| 国产av在哪里看| 国内毛片毛片毛片毛片毛片| 国产精品久久电影中文字幕| 老司机在亚洲福利影院| 99久久精品热视频| a在线观看视频网站| 最近最新中文字幕大全电影3| 国产视频一区二区在线看| 在线免费观看的www视频| 99riav亚洲国产免费| 亚洲男人的天堂狠狠| 老司机在亚洲福利影院| 十八禁网站免费在线| 日韩亚洲欧美综合| 9191精品国产免费久久| 亚洲精品一区av在线观看| 国产精品女同一区二区软件 | 一进一出抽搐动态| 国产精品99久久久久久久久| 国产不卡一卡二| 两个人的视频大全免费| 麻豆国产97在线/欧美| 少妇人妻精品综合一区二区 | 黄色日韩在线| 久久中文看片网| 久久九九热精品免费| 淫妇啪啪啪对白视频| 不卡一级毛片| 久久香蕉国产精品| 成年女人毛片免费观看观看9| 日韩精品中文字幕看吧| 久久99热这里只有精品18| 日韩高清综合在线| 天堂网av新在线| 一进一出抽搐gif免费好疼| 女生性感内裤真人,穿戴方法视频| 无人区码免费观看不卡| 日本熟妇午夜| 99精品欧美一区二区三区四区| 日韩亚洲欧美综合| 无限看片的www在线观看| 欧美在线黄色| 午夜免费激情av| 啦啦啦韩国在线观看视频| 久久精品人妻少妇| 91在线精品国自产拍蜜月 | 免费高清视频大片| 亚洲一区高清亚洲精品| 激情在线观看视频在线高清| 亚洲av电影不卡..在线观看| 精品一区二区三区视频在线观看免费| 国产精品自产拍在线观看55亚洲| 久久久国产成人精品二区| 亚洲最大成人中文| 美女高潮喷水抽搐中文字幕| 桃色一区二区三区在线观看| 黄色成人免费大全| 色综合站精品国产| 日本在线视频免费播放| 成人国产一区最新在线观看| 欧美成人a在线观看| 久久人妻av系列| 日韩欧美 国产精品| 国产精品爽爽va在线观看网站| 首页视频小说图片口味搜索| 亚洲成人免费电影在线观看| 九色成人免费人妻av| 亚洲人成伊人成综合网2020| 日韩精品中文字幕看吧| 十八禁网站免费在线| 一级黄片播放器| 亚洲最大成人中文| 高清日韩中文字幕在线| 51国产日韩欧美| 无限看片的www在线观看| 搡女人真爽免费视频火全软件 | 亚洲欧美日韩东京热| 久久精品国产亚洲av香蕉五月| 男人的好看免费观看在线视频| 成人一区二区视频在线观看| 国产不卡一卡二| 久久性视频一级片| 亚洲中文日韩欧美视频| 国产一区二区在线观看日韩 | 国产精品国产高清国产av| 91麻豆精品激情在线观看国产| 欧美日韩福利视频一区二区| h日本视频在线播放| 久久久久久久久大av| 精品一区二区三区视频在线 | 午夜福利在线观看吧| 夜夜躁狠狠躁天天躁| 亚洲精品亚洲一区二区| 亚洲五月天丁香| 久久久精品大字幕| 国产成人系列免费观看| 桃红色精品国产亚洲av| 免费一级毛片在线播放高清视频| avwww免费| 男人和女人高潮做爰伦理| 在线观看舔阴道视频| 国内揄拍国产精品人妻在线| 国产成人av教育| 757午夜福利合集在线观看| 香蕉av资源在线| bbb黄色大片| 又黄又爽又免费观看的视频| www.999成人在线观看| 亚洲欧美激情综合另类| 国产视频一区二区在线看| 久久久久久久午夜电影| 蜜桃久久精品国产亚洲av| 1000部很黄的大片| 日本 欧美在线| 免费在线观看日本一区| 午夜福利欧美成人| 非洲黑人性xxxx精品又粗又长| 成年版毛片免费区| 欧美日韩中文字幕国产精品一区二区三区| 欧美成人免费av一区二区三区| 不卡一级毛片| 一个人观看的视频www高清免费观看| 一级毛片女人18水好多| 亚洲成人中文字幕在线播放| 国产麻豆成人av免费视频| 老熟妇乱子伦视频在线观看| 中文字幕熟女人妻在线| 色尼玛亚洲综合影院| 国产精品综合久久久久久久免费| 岛国在线免费视频观看| 免费人成在线观看视频色| 91av网一区二区| 亚洲av免费高清在线观看| 欧美+日韩+精品| 99热只有精品国产| 成熟少妇高潮喷水视频| 亚洲精品亚洲一区二区| 亚洲av免费高清在线观看| 国产伦在线观看视频一区| 高清日韩中文字幕在线| 人妻久久中文字幕网| 麻豆久久精品国产亚洲av| 性色av乱码一区二区三区2| 观看免费一级毛片| 欧美bdsm另类| 夜夜夜夜夜久久久久| 在线天堂最新版资源| 嫩草影院精品99| 国产av一区在线观看免费| 神马国产精品三级电影在线观看| 久久久色成人| 亚洲精品一卡2卡三卡4卡5卡| 国产精品爽爽va在线观看网站| 黄色丝袜av网址大全| 神马国产精品三级电影在线观看| 最后的刺客免费高清国语| 午夜日韩欧美国产| 十八禁网站免费在线| 色av中文字幕| 日韩精品青青久久久久久| 国产精品 国内视频| 在线国产一区二区在线| 18禁裸乳无遮挡免费网站照片| 亚洲国产欧美网| 18禁国产床啪视频网站| 99久国产av精品| 99久久久亚洲精品蜜臀av| av视频在线观看入口| 国产欧美日韩精品一区二区| 在线视频色国产色| 91久久精品国产一区二区成人 | 精品人妻1区二区| 国产成人影院久久av| 一本综合久久免费| 琪琪午夜伦伦电影理论片6080| 亚洲第一电影网av| 无遮挡黄片免费观看| 一本一本综合久久| 久久香蕉精品热| 亚洲av一区综合| 人妻丰满熟妇av一区二区三区| 国产一区在线观看成人免费| 欧美黄色淫秽网站| 中文字幕av成人在线电影| 亚洲精品久久国产高清桃花| 老司机在亚洲福利影院| 18禁美女被吸乳视频| 高潮久久久久久久久久久不卡| 欧美av亚洲av综合av国产av| 99久久精品国产亚洲精品| 在线天堂最新版资源| 亚洲七黄色美女视频| 国产欧美日韩精品一区二区| 亚洲熟妇中文字幕五十中出| 两个人看的免费小视频| 美女免费视频网站| 麻豆国产av国片精品| 网址你懂的国产日韩在线| 久久久久久国产a免费观看| 成人18禁在线播放| 日韩成人在线观看一区二区三区| 很黄的视频免费| 成年女人看的毛片在线观看| 老鸭窝网址在线观看| netflix在线观看网站| 岛国视频午夜一区免费看| 亚洲av不卡在线观看| 嫩草影视91久久| 国产亚洲精品久久久com| 欧美精品啪啪一区二区三区| 午夜亚洲福利在线播放| 香蕉久久夜色| 久久精品影院6| 内地一区二区视频在线| 3wmmmm亚洲av在线观看| 成年女人毛片免费观看观看9| 伊人久久大香线蕉亚洲五| 男插女下体视频免费在线播放| 国内精品一区二区在线观看| 久久久久久大精品| 99热精品在线国产| 18禁黄网站禁片免费观看直播| 国产成人影院久久av| 免费人成在线观看视频色| 亚洲精品一区av在线观看| 久久久久国内视频| 手机成人av网站| 噜噜噜噜噜久久久久久91| 欧美一级a爱片免费观看看| 欧美bdsm另类| 亚洲av电影不卡..在线观看| 中文字幕人成人乱码亚洲影| 国产亚洲欧美98| 大型黄色视频在线免费观看| 国产美女午夜福利| 麻豆久久精品国产亚洲av| 久久久久国内视频| 女同久久另类99精品国产91| 中亚洲国语对白在线视频| 欧美+日韩+精品| 免费av毛片视频| 免费在线观看亚洲国产| 此物有八面人人有两片| 美女 人体艺术 gogo| 欧美日韩黄片免| a在线观看视频网站| 亚洲色图av天堂| 91在线精品国自产拍蜜月 | 亚洲精华国产精华精| 成人高潮视频无遮挡免费网站| 中文字幕熟女人妻在线| 欧美日本视频| 噜噜噜噜噜久久久久久91| 国产av麻豆久久久久久久| 美女高潮喷水抽搐中文字幕| 男女那种视频在线观看| 岛国视频午夜一区免费看| 尤物成人国产欧美一区二区三区| 国产不卡一卡二| 成年女人毛片免费观看观看9| 亚洲精品一区av在线观看| 无人区码免费观看不卡| 国内少妇人妻偷人精品xxx网站| 日本一二三区视频观看| 久久精品91蜜桃| 日韩高清综合在线| 亚洲国产高清在线一区二区三| 欧美xxxx黑人xx丫x性爽| 亚洲色图av天堂| 岛国在线免费视频观看| 手机成人av网站| 亚洲欧美一区二区三区黑人| 婷婷亚洲欧美| 久久天躁狠狠躁夜夜2o2o| 亚洲色图av天堂| 99国产精品一区二区蜜桃av| 一个人看的www免费观看视频| 性色av乱码一区二区三区2| 亚洲在线自拍视频| 欧美bdsm另类| 国产成人啪精品午夜网站| 色噜噜av男人的天堂激情| 中出人妻视频一区二区| 搡老岳熟女国产| 免费看日本二区| 亚洲自拍偷在线| 亚洲人成伊人成综合网2020| 日本与韩国留学比较| 亚洲人成伊人成综合网2020| 久久人人精品亚洲av| 免费在线观看成人毛片| 午夜免费成人在线视频| 性欧美人与动物交配| 亚洲成人精品中文字幕电影| 蜜桃亚洲精品一区二区三区| 欧美bdsm另类| www.熟女人妻精品国产| 高清日韩中文字幕在线| 夜夜躁狠狠躁天天躁| 波多野结衣高清作品| 亚洲在线自拍视频| 麻豆国产97在线/欧美| 成人欧美大片| 露出奶头的视频| 亚洲中文日韩欧美视频| 国内毛片毛片毛片毛片毛片| 天天一区二区日本电影三级| 久久亚洲真实| 国产精品久久久久久久电影 | 女人高潮潮喷娇喘18禁视频| 久久草成人影院| 国产一区二区在线观看日韩 | 色在线成人网| 欧美极品一区二区三区四区| www.色视频.com| 日韩欧美国产在线观看| 国产精品一区二区三区四区久久| 免费人成在线观看视频色| 有码 亚洲区| 免费人成在线观看视频色| 有码 亚洲区| 精品久久久久久久人妻蜜臀av| 成人特级av手机在线观看| 精品久久久久久久人妻蜜臀av| 成人特级av手机在线观看| 成人18禁在线播放| 搡老岳熟女国产| 成熟少妇高潮喷水视频| 久久久久久久久大av| 国产色爽女视频免费观看| 久久久久久久久久黄片| 大型黄色视频在线免费观看| 久久6这里有精品| 亚洲国产精品合色在线| 啪啪无遮挡十八禁网站| 在线观看66精品国产| www.www免费av| 午夜两性在线视频| 国产真实乱freesex| 国产亚洲精品久久久久久毛片| 最后的刺客免费高清国语| 国产亚洲精品久久久久久毛片| 久久九九热精品免费| 99久久久亚洲精品蜜臀av| 午夜激情欧美在线| 91av网一区二区| 欧美成人免费av一区二区三区| 午夜福利欧美成人| 精品久久久久久久久久久久久| 免费电影在线观看免费观看| 午夜精品一区二区三区免费看| 亚洲18禁久久av| 国产精品野战在线观看| 美女高潮喷水抽搐中文字幕| 欧美一级a爱片免费观看看| 老汉色av国产亚洲站长工具| 18+在线观看网站| 在线观看美女被高潮喷水网站 | 黄色成人免费大全| 男人和女人高潮做爰伦理| 两个人的视频大全免费| 国产中年淑女户外野战色| 大型黄色视频在线免费观看| 久久精品国产99精品国产亚洲性色| 午夜影院日韩av| 成人亚洲精品av一区二区| 黄色女人牲交| 欧美最黄视频在线播放免费| 欧美性猛交╳xxx乱大交人| 欧美日韩瑟瑟在线播放| 一边摸一边抽搐一进一小说| 每晚都被弄得嗷嗷叫到高潮| 欧美色欧美亚洲另类二区| 亚洲精品成人久久久久久| www.熟女人妻精品国产| 亚洲内射少妇av| 久久久久久久久中文| 搡老妇女老女人老熟妇| 麻豆成人av在线观看| 亚洲最大成人手机在线| 成人国产一区最新在线观看| 色综合欧美亚洲国产小说| 一夜夜www| 99久久精品一区二区三区| 99国产综合亚洲精品| 欧美日韩瑟瑟在线播放| 免费观看的影片在线观看| 一进一出抽搐gif免费好疼| xxxwww97欧美| 最新中文字幕久久久久| 午夜精品一区二区三区免费看| 欧美性猛交黑人性爽| 久久精品国产亚洲av涩爱 | 欧美极品一区二区三区四区| 久久久久久久久久黄片| 午夜视频国产福利| 亚洲欧美日韩卡通动漫| 国产极品精品免费视频能看的| 51国产日韩欧美| 搡老岳熟女国产| av黄色大香蕉| 欧美性猛交黑人性爽| 亚洲欧美日韩东京热| 18禁裸乳无遮挡免费网站照片| 国产高清视频在线播放一区| 久久草成人影院| 精品国产美女av久久久久小说| 高潮久久久久久久久久久不卡| 中文字幕熟女人妻在线| 亚洲av二区三区四区| 精品一区二区三区视频在线 | 国产aⅴ精品一区二区三区波| 国产三级黄色录像| 天天一区二区日本电影三级| 毛片女人毛片| 国产av在哪里看| 欧美中文综合在线视频| 久久精品影院6| 欧美乱色亚洲激情| www.www免费av| 国产乱人视频| 色av中文字幕| 国产精品久久久久久精品电影| 久久人妻av系列| 亚洲熟妇中文字幕五十中出| 一级毛片女人18水好多| 亚洲无线在线观看| 婷婷亚洲欧美| 一本一本综合久久| 欧美av亚洲av综合av国产av| 两个人看的免费小视频| 叶爱在线成人免费视频播放| 99视频精品全部免费 在线| 亚洲人与动物交配视频| 最新中文字幕久久久久| 欧美最黄视频在线播放免费| 国产高清有码在线观看视频| 身体一侧抽搐| 久久精品国产99精品国产亚洲性色| 国产91精品成人一区二区三区| 看片在线看免费视频| а√天堂www在线а√下载| 成年女人看的毛片在线观看| 成人鲁丝片一二三区免费| 18禁国产床啪视频网站| 91九色精品人成在线观看| 一边摸一边抽搐一进一小说| 人人妻,人人澡人人爽秒播| 国产精品永久免费网站| 午夜福利免费观看在线| 黄色日韩在线| 午夜激情福利司机影院| 国产精品av视频在线免费观看| 欧美最黄视频在线播放免费| 久久久久久久精品吃奶| 三级国产精品欧美在线观看| 老司机午夜福利在线观看视频| 国产私拍福利视频在线观看| 18+在线观看网站| 欧美黑人巨大hd| 欧美一区二区精品小视频在线| 中出人妻视频一区二区| 老汉色av国产亚洲站长工具| 一个人看视频在线观看www免费 | 日本与韩国留学比较| 国产av一区在线观看免费| 99在线人妻在线中文字幕| 国产探花极品一区二区| 一本久久中文字幕| 精品人妻1区二区| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲欧美在线一区二区| 手机成人av网站| av女优亚洲男人天堂| 在线视频色国产色| 99热这里只有是精品50| 国产精品女同一区二区软件 | 激情在线观看视频在线高清| 精品欧美国产一区二区三| 国产真实伦视频高清在线观看 | 国产精品久久电影中文字幕| 亚洲精品国产精品久久久不卡| 久久久久国产精品人妻aⅴ院| 欧美色视频一区免费| 无遮挡黄片免费观看| 国产色爽女视频免费观看| 欧美性感艳星| 亚洲中文字幕一区二区三区有码在线看| 美女高潮的动态| 在线免费观看的www视频| 香蕉久久夜色| 国产精品嫩草影院av在线观看 | 在线十欧美十亚洲十日本专区| 亚洲av免费高清在线观看| 97人妻精品一区二区三区麻豆| 夜夜躁狠狠躁天天躁| 欧美zozozo另类| 老熟妇仑乱视频hdxx| 欧美国产日韩亚洲一区| 极品教师在线免费播放| 国产一级毛片七仙女欲春2| 久久草成人影院| 国产一区二区在线av高清观看| 婷婷精品国产亚洲av| 可以在线观看的亚洲视频| 国产蜜桃级精品一区二区三区| 韩国av一区二区三区四区| 欧美成狂野欧美在线观看| www日本黄色视频网| 美女被艹到高潮喷水动态| 亚洲五月天丁香| 人人妻人人澡欧美一区二区| 欧美一级a爱片免费观看看| 69av精品久久久久久| 日韩精品青青久久久久久| 国产 一区 欧美 日韩| 免费观看人在逋| 成人av一区二区三区在线看| 亚洲国产中文字幕在线视频| 丰满乱子伦码专区| 麻豆成人av在线观看| 老司机午夜十八禁免费视频| 久久久久久久久大av| 国产精品亚洲美女久久久| 亚洲国产色片| 婷婷六月久久综合丁香| 午夜免费男女啪啪视频观看 | 日韩欧美 国产精品| 内地一区二区视频在线| 噜噜噜噜噜久久久久久91| 国产极品精品免费视频能看的| www国产在线视频色| 国产高清视频在线观看网站| 无限看片的www在线观看| 亚洲最大成人手机在线| 97人妻精品一区二区三区麻豆| 国产激情欧美一区二区| 国产精品久久视频播放| 国产精品久久久久久人妻精品电影| 日韩欧美 国产精品| 精品人妻一区二区三区麻豆 | 中文字幕久久专区| 中亚洲国语对白在线视频| 国产真实伦视频高清在线观看 | 国产成人av激情在线播放|