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

    Local blast wave interaction with tire structure

    2020-06-28 03:02:14PaweBaranowskiJerzyMaachowskiukaszMazurkiewicz
    Defence Technology 2020年3期

    Pawe? Baranowski, Jerzy Ma?achowski, ?ukasz Mazurkiewicz

    Faculty of Mechanical Engineering, Military University of Technology, Department of Mechanics and Applied Computer Science, Gen. Sylwester Kaliski Street, 00-908 Warsaw, Poland

    Keywords:Blast wave Tire structure Finite element analysis Optimization

    ABSTRACT The paper presents modelling and simulation of a local blast wave interaction with a tire of a logistic truck supporting military operations. In the military industry, it is desired to improve effectiveness and strength of vehicle components and simultaneously to minimize the risk of crew injuries.Therefore,the main goal of this paper is to present an attempt to improve blast resistance of a tire. Based on the developed, verified and validated finite element model an optimization procedure was conducted to minimize the damage of a tire subjected to a blast wave.The main issue in the performed computations was to estimate an influence of the cords angle in each layer.For this purpose,a pre-processor script was developed in order to easily modify the finite element model so that the generation process was performed automatically on the basis of optimization variables.Dynamic response of a tire subjected to blast wave in different cases (cords configurations) was analyzed. It was shown that the optimized cords angles configuration reduces tire local damage and increases its blast resistance.

    1. Introduction

    The paper relates to the events where high explosives are commonly used. In many cases the blast wave can cause serious injuries of people and unfortunately these injures are very often deadly.It shows that there is a need to increase the blast resistance of structures which are exposed to explosion. Therefore, an application of various protective shields was found to be very effective.The main goal of them is to minimize the blast energy as much as possible [1-9]. In the military industry, it is desired to improve effectiveness and strength of vehicle components and simultaneously to minimize the risk of crew injuries.The explosions under vehicle can destroy wheels or even the whole suspension system resulting in an immobilization of the vehicle(Fig.1).Therefore,the presented study is focused on an attempt to improve blast resistance of a tire by modifying its interior properties.At first glance a tire does not resemble a composite. However, it is a very complex structure consisting of rubber parts (tread, sidewall, etc.), which can be considered as the matrix. Additionally, it is reinforced with cords, which, depending on tire application, can be made of different materials: nylon, polyester or steel and can have various configurations.Immersed with different orientations in the rubber parts,they form a special ring-like laminate [10-15].

    Development of a pneumatic tire with a high operating standard is associated with a series of experimental studies and numerical simulations to determine stability and reliability of its implementation. It also involves large financial costs and time constraints. Scientists have conducted tire (wheel) testing under different load conditions. In Ref. [14] the measurement of force acting on the tire surface is shown which is a result of low-profile tire tread slipping. In Ref. [16] the numerical and experimental study of high explosive (HE) explosion under an armored vehicle(LAV) is presented, the mesh topology influence on the obtained results is also discussed.

    Finite element analysis (FEA) is also adopted in tire industry which enables estimation of deformation and stress state analysis in a tire and to perform necessary design modifications,even before the production stage. Thus, for effective and correct analyses, a precise numerical model of tire should be developed[10,11,13,14,17,18]. It was shown that the final form and characteristics of a developed finite element(FE)model of a tire is affected by such factors as knowledge and availability of the input and material data.

    Optimization analyses or parametric studies of tire were also carried out by several authors [19-23]. The neural network together with the finite element method (FEM) was used in Ref.[19]to optimize the cornering property of the motorcycle tire.A new tire design procedure with FEM, analytical method and optimization technique using genetic algorithm was developed and verified in drum tests and field tests in paper [22]. Moreover, the authors of the paper[20]presented the possibility of FEA utilization to optimize the tire operating characteristics,where such factors as material property or spacing and angles of the cords were analyzed.Some attempts to find an optimum belt construction for tire weight and rolling resistance can be found in Ref. [21].

    The above-mentioned optimization studies were applied to optimize basic tires operating characteristics. Strongly dynamic load conditions such as blast loading impacts were not introduced in the modelling.Therefore,in this paper a method of optimization of tire subjected to blast loading based on computational mechanics is presented. The objective is to increase strength and resistance of a tire loaded with blast wave.Based on the developed,verified and validated finite element (FE) model an optimization procedure was conducted to minimize the damage of a tire and the blast energy. The influence of the cords angle in each layer on obtained results was analyzed.For this purpose,a pre-processor script was developed to easily modify the finite element model so that the generation process could be performed automatically on the basis of optimization variables.Dynamic response of the tire subjected to the blast wave in different cases was analyzed. Obtained results clearly showed that application of the optimized cords angles configuration reduces the tire damage and increases the blast resistance. The study is a continuation of the previous authors'investigations related to tires modelling and simulation[14,17,24-26].

    2. Analyzed problem

    In the previous studies a novel concept of a tire modelling was validated [17,26]. For this purpose, radial deflection test, bounce test and blast loading test were considered.The last one,in fact, is used by the authors in the presented optimization procedure study.Therefore,only experimental and numerical blast loading tests will be briefly discussed.

    Fig. 2. Experimental blast loading tests: a) test 1: HE located 95 mm above tire surface, b) test 2: HE directly on the tire surface [26].

    Fig. 3. Tread damage experimental and FEA comparison, a) test 1, b) test 2 [26].

    Fig. 4. CT scan of a) tread and b) sidewall with visible cords layers arrangement (additional photo of cords size in microscope scale) [26].

    The explosive charge (HE) with cylindrical shape (50 mm diameter and 159 mm height) and a mass equivalent to 0.5 kg of TNT was used in two different experimental tests. Such HE was selected in order to easily test a tire under strongly dynamic conditions and to analyze rubber behavior under high-strain rates. In test 1 HE was 95 mm above the tire surface with tire inflation pressure of 0.60 MPa (Fig. 2(a)). Test 2 considered HE placed directly on the tire surface(close range explosion with a larger local damage) without internal tire pressure (Fig. 2(b)). Numerical simulations with correspondent initial boundary conditions were performed using LS-DYNA code with an implementation of Multi-Material Arbitrary Lagrangian-Eulerian (MM-ALE) formulation with Jones-Wilkins-Lee (JWL) equation for defining HE [25,27,28].

    During the experimental tests as well as numerical simulations a significant damage of the tire was obtained. Using the 3D laser scanner system,the geometry of the actual tire after experimental tests was measured.It was compared with discrete representation of the tire (Fig. 3). It was shown that the damage of the tread is similar and discrepancies between the results are relatively small.

    Fig. 5. Cords layers width and placement within the tread and sidewall [17,26].

    The results indicated that there is a room for tire strength improvement.Thus,the authors took an advantage of optimization algorithms coupled with computational mechanics in analyses of the vehicle tire.

    Fig. 6. Tire partly embedded in Eulerian air domain.

    Fig. 7. Mechanical characteristics obtained from experimental tests.

    Table 1 Steel parameters for the simplified JC model used in analyses [26,36].

    Table 2 TNT parameters for the HE material and JWL equation of state[36].

    Table 3 Optimization variables used in the study.

    3. Model description

    3.1. Tire numerical modelling

    In the paper,the developed,verified and validated FE model of a Michelin 315/70 R22.5 tire was used. A detailed description of adopted modelling techniques, material data and a validation procedure can be found in authors'previous papers[17,26].One of the crucial elements of the finite modelling was to verify tire cords arrangement and their size.By using a CT scanning technique and a microscope it was observed that a single cord consists of 9 smaller wires, with a total average diameter of 1.20 mm (Fig. 4). The steel cords are arranged radially inside the tire sidewall (Fig. 4(a)),whereas within the tread area,4 layers are placed circumferentially under different angles (Fig. 4(b)).

    Brick elements were adopted (constant stress solid elements)for modelling of the tire,whereas cable(truss)elements were used for representing the steel cords [27]. Their nodes were separated from the brick elements nodes which was possible by using the special feature[29]based on the penalty method.Such a method is original and different than those found in literature[13,14,30,31].In Fig. 5 the FE model of the tire is presented where different width and angle of each cord layer can be distinguished.The airbag model was adopted for representing the pressure inside [29].

    3.2. Blast loading modelling

    In order to carry out an optimization procedure,it is required to analyze a very large number of FE models with different design variables.Results obtained from the analyses using these FE models can be used for developing a metamodel(response surface).Such a model describes an influence of design variables on objective function and can be used for selecting the best variable set[2].The large number of FE models is also required in applying the Direct Genetic Algorithm optimization, which in fact was used in the presented study.This requires many simulations.Therefore,the FE model needs to be relatively simple to minimize computational time and,on the other hand,accurate enough to represent a trend of objective function change.Based on the authors'studies[17,26]it was decided to modify the adopted “global” FE model of the tire with Eulerian air domain, where mainly local effects were observed. To simplify the model a slice of the tire was considered,the Eulerian domain was reduced, and the following components were removed: fabric surface for airbag, support stand elements,rigid wall and rim. An additional closing surface for the airbag modelled with the use of the Lagrangian Belytschko-Tsay(BT)shell elements was added. The ambient domain was replaced with a non-reflecting boundary condition to prevent artificial stress wave reflections generated at the model boundaries. For the optimization,only test 1 was considered(described briefly in section 2)with the cylinder-shaped explosive charge(HE)made of 0.5 kg TNT with diameter of 50 mm and height of 159 mm. The HE was placed 95 mm above the tire surface with tire inflation pressure of 0.60 MPa.

    Fig. 8. Cords layers angles in respect to layer of radial cords taken as variables.

    Fig. 9. Computational strategy of the optimization procedure.

    Fig.10. Graphical presentation of three optimization procedures stages.

    The procedure of blast loading simulations was the same as for the “global” FE model presented in Refs. [17,26]. The boundary conditions applied to the tire at the bottom airbag surface and displacement in the perpendicular axis to the side airbag surfaces were fixed. The tire was partially embedded in the Eulerian air domain(Fig.6)and the above-mentioned non-reflecting boundary condition was implemented. The air was considered as vacuum.The interaction between the Eulerian blast wave and Lagrangian parts was simulated using the penalty coupling method utilizing the Van Leer's second-order advection theory [29].

    Fig.11. Sampling points results; a) cords layer no.1: α1, b) cords layer no. 2: α2, c) cords layer no. 3: α3, d) cords layer no. 4: α4.

    Fig.12. Objective function vs. time histories of all cases (sampling points - all considered cords configurations).

    Fig.13. Internal energy history of tire components before and after optimization.

    In the presented study the main attention was to conduct very detail studies oriented on the local effects (deformation, failure,blast wave reflection, Fluid Structure Interaction (FSI), etc.).Therefore, within the area of direct interaction between the pressure wave and the tire structure, a very fine mesh was applied(2-3 mm elements in the central area of the model). Based on previous investigations and author's experience the following parameters to prevent leakage in FSI interface were introduced.User defined interface stiffness as relation between pressure and penetration was applied.This linear relation was based on element size(maximum allowable penetration should be lower than edge length of structure elements). Moreover, reflected wave pressure from preliminary analysis (defined maximum interface pressure)should be larger than maximum pressure obtained from analysis.Second very important parameter was minimum volume fraction of detonation products(PD)in Eulerian elements to activate coupling(FRCMIN). It was significantly reduced to 10e-2which drastically minimizes the leakage problem, especially in problems with high velocity fluids impacts. Third parameter taken into consideration was a number of coupling points distributed over each coupled Lagrangian surface segment (NQUAD). It was increased to 3×3 points in order to ensure that minimum two coupling points per each Eulerian element will be generated. This guaranteed that no leakage was observed on the pressure distribution on the tire surface. This also gave the possibility to evaluate more effectively the influence of the cords angle on the tire strength.

    3.3. Material description

    Simplified Rubber Material (SRM) model was used for the rubber parts [17,26,32,33], where strain rate sensitivity is directly implemented using a stress-strain curves table for different strain rates [33,34]. The mechanical characteristics were obtained from the quasi-static and dynamic tests,which are described in Ref.[35](Fig. 7). The erosion in the rubber parts was simulated using a failure criterion based on the fracture strain value taken from the experimental tension test.

    For the cords elements the simplified Johnson-Cook (JC) model with isotropic hardening and strain rate effects was applied with the material parameters taken from Refs.[26,36](Table 1).The HE material as well as the detonation process was simulated using an“explosive burn” material model with (JWL) equation of state.The parameters for the TNT were obtained from analytical calculations supplemented with literature data [36] (Table 2). For a detailed description of material data please see authors other papers[17,26,35].

    4. Optimization procedure

    4.1. Procedure definition

    The main goal of the study is to increase strength and resistance of the tire loaded with blast wave. In order to achieve this goal a method of tire structure optimization subjected to blast loading based on computational mechanics was implemented. The optimization procedure was based on genetic algorithms[37]with 50 populations and 10 generations.The optimization objective was to minimize the tire damage, which from the FE analyses point of view is represented by the eroding internal energy. In each numerical simulation a value of the eroding internal energy at t=0.3 ms was taken as optimization objective.It was observed that after this time the tire damage did not increase. Therefore, it was taken as termination time in each simulation. The problem can be described by the following equation:

    where:mindobjectiveis tire tread eroding internal energy;α1,α2,α3,α4is angle of corresponded cords layers (Table 3).

    At the first stage,the angle of the circumferential cords(Table 3)in each layer with respect to wheel axis were considered as the variable (Fig. 8). For this purpose, a pre-processor script was developed in order to easily modify the 2D FE model which could be automatically generated depending on the angles used.Despite different angles of the cords in each case,the width of layers as well as the distance between them was the same.The 2D FE model was developed using a cross section of the tire which was dragged in one direction.Thanks to such an approach,generation procedure of each cords layer with different angle was easier than in case of the 3D FE model.In the second stage,generated 2D FE model with a set of applied variables was transformed into the 3D FE model using Fortran script.Therefore,the 2D FE model was curved to represent a slice of the tire. At the last stage, Eulerian domain was added to the curved 3D FE model and the simulation of tire inflation followed by the blast wave interaction with the tire was carried out.The computational strategy is presented in Fig. 9, whereas the multistage optimization procedure is graphically presented in Fig.10.

    5. Optimization results and discussion

    The optimization procedure optimum case was found after analyzing 5 generation and 150 sampling points. In Fig. 11 all sampling points for each variable are presented. The optimum angles are also highlighted for which the minimum objective value(eroding internal energy) was obtained at t=0.3 ms. It can be noticed that the optimum variables are as follows:α1=45°,α2=90°,α3= 105°,α4= 150°. This is reflected in Fig. 12 where historical data of eroding internal energy of the tire for all sampling points (all considered cords configurations) is presented. Additionally, an area within the termination time t=0.3 ms was enlarged for a better presentation of the results. The final point of each curve corresponds to the objective value that was checked during optimization.The minimum value of the objective function reached 0.806 kJ, whereas in the basic model this value was 0.829 kJ.

    Fig.14. Blast wave with tire interaction comparison in selected time frames.

    The influence of cords arrangement had a direct impact on the amount of the energy that the tire absorbs.In Fig.13 internal energy absorbed by the tire components is presented for the basic tire model and the optimized model. Internal energy is calculated incrementally (for each element) based on the six components of stress and strain tensorial values. Bigger tire deformation and consequently greater strains in the basic model, resulted in a greater number of eroded elements in the model. In this case the cords absorbed smaller amount of energy. An inverse relationship can be noticed for the optimized model: internal energy of tread and sidewall was smaller, whereas for the cords it was larger.

    For a better presentation of the described phenomena,the blast wave distribution within the Eulerian domain for both the basic and the optimized tire model is shown in Fig. 14 for three time frames.It can be noticed that at the initial stage the characteristic of tire damage is different and larger damage was obtained in the basic tire model. At t2=0.2 ms blast, pressure flows throughgenerated hole, thus the greater flow is also clearly visible in the model before optimization.Due to an increased strength of the tire a significantly reduced flow is evident in the optimized model at t3=0.3 ms. Comparison of the optimized and basic models is presented in Table 4 where objective function is compared as well as internal energy values of tread,sidewall and cords at t=0.3 ms are listed.Despite a relatively small influence of the angles of cords on the objective function(~3%) the optimum configuration decreased the amount of consumed energy by tread and sidewall up to approximately 10%. Moreover, maximum differences between the objective function values are 6.4%(see Fig.12).It shows a tendency of tire modification which can ultimately lead to a military tire which will be resistant to local damage and consequently will increase safety of a vehicle.

    Table 4 Optimization results comparison.

    6. Conclusions

    A tire is a part which is directly exposed to the blast wave,which comes from a detonation of HE under a vehicle.After the blast wave interaction with the wheel(tire),the reflected wave propagates and starts to interact with other LAV mechanical parts. Therefore, to precisely predict effects of IED explosion under armored vehicles a detailed modelling of the tire is very important. In this paper the authors demonstrate capabilities of numerical optimization based on computational mechanics to increase strength and resistance of tire structure loaded with blast wave. The optimization was based on the simplified FE model of the tire with Eulerian air. The FE model which was used in the study was based on the verified and validated tire, which was investigated in the previous authors'studies [17,26].

    Based on the above-mentioned papers but also on the presented study it was found that the tire FE model with SRM constitutive model can be used and is adequate for simulating the strongly dynamic problems.Additionally,the implemented failure criterion based on the strain value allowed simulating of erosion in the rubber parts with overall reasonable reliability. Moreover, the method used for representing the cords inside tire tread and sidewall is novel and original. The main advantage of such an approach is that the cords nodes can be separated from the brick elements of tread and sidewall.This enables accurate modelling of tire without the need to coincide the nodes of cords and tire elements. Based on this approach it was possible to analyze the influence of the cords angle in each layer and to implement a preprocessor script to easily modify and automatically generate the FE model depending on the optimization variables.

    As a results of the optimization investigations dynamic response of the tire structure subjected to the blast wave in different cases was analyzed. From the optimization procedure, 5 generation and 150 sampling points were obtained within which the optimum case was found. The tire after optimization contains cords layers with the following cords angles: α1= 45°,α2= 90°,α3= 105°,α4=150°. Obtained results showed that application of the optimized cords angles configuration reduces the tire damage and increases the blast resistance.

    The presented optimization methodology has proved its effectiveness and that it can be applied for modifying and analyzing the tire structure. Proposed FE tire model representation with cords coupled with tread and sidewall, gives the possibility to apply the pre-processor script, which greatly facilitates the generation of many models needed for optimization. Authors are aware of the fact, that in the optimum model the observed damage (erosion of the finite elements)was approximately 3%higher than in the basic model.However,in the study only cords angles were modified,and the authors' intention was to present the possibilities of optimization procedure implementation for tire structure. Other studies,i.e.[38]show that there is still a room for improving the strength of tire which is locally damaged by blast wave or gunshot puncturing.The presented paper shows a possible direction of tire modification which consequently will increase safety of a vehicle. In future studies the investigations will be continued using cords configurations with different materials, thickness, etc. Hopefully, an optimum cords arrangement and properties will be obtained which will result in significant reduction of the damage or in complete prevention of tire rupture. However, the new cords configuration needs to be verified in tire operating conditions and only after obtaining positive results it can be considered as ready for implementation in a military vehicle.

    Acknowledgments

    The research was carried out under a research grant no. RMN 723 and the support of the Interdisciplinary Centre for Mathematical and Computational Modelling (ICM) University of Warsaw under grant no GB73-19. This support is gratefully acknowledged.

    少妇的逼好多水| 亚洲欧美一区二区三区国产| 六月丁香七月| 久久6这里有精品| 有码 亚洲区| 九九久久精品国产亚洲av麻豆| 中文天堂在线官网| 最近中文字幕2019免费版| 大陆偷拍与自拍| h日本视频在线播放| 久久久欧美国产精品| 天堂av国产一区二区熟女人妻| 欧美一级a爱片免费观看看| 国产淫语在线视频| 啦啦啦啦在线视频资源| 国产成人aa在线观看| 亚洲人与动物交配视频| 亚洲精品第二区| 国产乱来视频区| 青春草国产在线视频| 美女高潮的动态| 久久草成人影院| 蜜桃亚洲精品一区二区三区| 成人二区视频| 中文字幕制服av| 久久久久久久国产电影| 亚洲精品乱码久久久v下载方式| 麻豆乱淫一区二区| 国产精品一二三区在线看| 99热6这里只有精品| 三级毛片av免费| 熟女人妻精品中文字幕| 成人av在线播放网站| 免费在线观看成人毛片| av女优亚洲男人天堂| 国产片特级美女逼逼视频| 2021少妇久久久久久久久久久| 亚洲最大成人手机在线| 国产精品一区www在线观看| 三级毛片av免费| 建设人人有责人人尽责人人享有的 | 伦理电影大哥的女人| 国产黄片美女视频| 插逼视频在线观看| 日本免费在线观看一区| 全区人妻精品视频| h日本视频在线播放| 久久精品国产亚洲av涩爱| 内射极品少妇av片p| 亚洲av二区三区四区| 久久久久久久久久人人人人人人| 精品不卡国产一区二区三区| 国内精品美女久久久久久| 天堂中文最新版在线下载 | 久久久久久久亚洲中文字幕| 成人一区二区视频在线观看| 中文字幕久久专区| 亚洲精品一区蜜桃| 国产大屁股一区二区在线视频| 一级爰片在线观看| 精品久久久噜噜| 赤兔流量卡办理| 波多野结衣巨乳人妻| 成人午夜精彩视频在线观看| 人人妻人人看人人澡| 高清视频免费观看一区二区 | 国产一区亚洲一区在线观看| 精品久久国产蜜桃| 国产黄片视频在线免费观看| 中国国产av一级| 一级二级三级毛片免费看| 一夜夜www| 国模一区二区三区四区视频| 亚洲精品日本国产第一区| 简卡轻食公司| 国产成人aa在线观看| 99久国产av精品国产电影| 大香蕉97超碰在线| 亚洲va在线va天堂va国产| 五月天丁香电影| 51国产日韩欧美| 色视频www国产| 久热久热在线精品观看| 亚洲不卡免费看| 欧美变态另类bdsm刘玥| 日韩国内少妇激情av| 人妻制服诱惑在线中文字幕| 国产有黄有色有爽视频| 久久热精品热| 日韩欧美精品v在线| 少妇被粗大猛烈的视频| 美女国产视频在线观看| 午夜激情久久久久久久| 久久精品综合一区二区三区| 一夜夜www| 91久久精品国产一区二区成人| 超碰av人人做人人爽久久| 天天躁日日操中文字幕| 成人无遮挡网站| 在线免费十八禁| 精品一区二区三区视频在线| 26uuu在线亚洲综合色| 舔av片在线| 色吧在线观看| 亚洲成色77777| 在线观看人妻少妇| 久久久久久久大尺度免费视频| 日日啪夜夜爽| 精品久久久久久久久av| 成人漫画全彩无遮挡| 毛片女人毛片| 久久精品久久久久久久性| 国产久久久一区二区三区| 亚洲三级黄色毛片| 色5月婷婷丁香| 插逼视频在线观看| 91午夜精品亚洲一区二区三区| 黄片wwwwww| 国产成年人精品一区二区| 在线观看美女被高潮喷水网站| 18禁在线无遮挡免费观看视频| 51国产日韩欧美| 亚洲国产成人一精品久久久| 狠狠精品人妻久久久久久综合| 一级毛片aaaaaa免费看小| 国产午夜精品论理片| 欧美日韩一区二区视频在线观看视频在线 | 亚洲怡红院男人天堂| 日韩欧美三级三区| 亚洲精品成人久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲精品久久久com| 女人被狂操c到高潮| 国产成人精品一,二区| 亚洲最大成人中文| av福利片在线观看| 亚洲av日韩在线播放| 亚洲av国产av综合av卡| 日韩精品有码人妻一区| 日韩三级伦理在线观看| 国产黄片视频在线免费观看| 大香蕉97超碰在线| 亚洲人成网站在线播| 18禁动态无遮挡网站| 少妇熟女欧美另类| 免费看日本二区| 欧美97在线视频| 国内精品一区二区在线观看| 神马国产精品三级电影在线观看| 丰满人妻一区二区三区视频av| 在线观看美女被高潮喷水网站| 成年人午夜在线观看视频 | 国产成人一区二区在线| 国产精品1区2区在线观看.| 国产av在哪里看| 床上黄色一级片| 大片免费播放器 马上看| 国产69精品久久久久777片| 亚洲丝袜综合中文字幕| 国产男女超爽视频在线观看| 国产一区二区三区av在线| 国产av不卡久久| 一级毛片我不卡| 青春草国产在线视频| 国产黄片视频在线免费观看| 国产黄频视频在线观看| 高清日韩中文字幕在线| a级毛片免费高清观看在线播放| 美女内射精品一级片tv| 啦啦啦啦在线视频资源| 久久精品熟女亚洲av麻豆精品 | 亚洲精品亚洲一区二区| 久久精品久久久久久噜噜老黄| 噜噜噜噜噜久久久久久91| 麻豆成人午夜福利视频| 国产大屁股一区二区在线视频| 日韩国内少妇激情av| 午夜福利高清视频| 亚洲av免费在线观看| 嫩草影院入口| 特大巨黑吊av在线直播| 国产精品人妻久久久久久| 国产高潮美女av| 久久这里只有精品中国| 18禁在线播放成人免费| 国产精品熟女久久久久浪| 欧美日韩国产mv在线观看视频 | 亚洲精品,欧美精品| 午夜激情久久久久久久| av国产久精品久网站免费入址| 成人二区视频| 国产在线男女| 成人毛片60女人毛片免费| 欧美bdsm另类| 免费看美女性在线毛片视频| 国产精品嫩草影院av在线观看| 特大巨黑吊av在线直播| 在线播放无遮挡| 国产免费福利视频在线观看| 国产真实伦视频高清在线观看| 国产精品三级大全| 能在线免费看毛片的网站| 国产69精品久久久久777片| 成人综合一区亚洲| 精品国产一区二区三区久久久樱花 | 大话2 男鬼变身卡| 亚洲欧美成人综合另类久久久| 最近视频中文字幕2019在线8| 日韩欧美精品免费久久| 中文字幕免费在线视频6| 中文字幕久久专区| 美女脱内裤让男人舔精品视频| 国产一区二区在线观看日韩| 成人漫画全彩无遮挡| 成年人午夜在线观看视频 | 九草在线视频观看| 亚洲欧美日韩无卡精品| 床上黄色一级片| 国产精品美女特级片免费视频播放器| 免费黄频网站在线观看国产| 日本熟妇午夜| 亚洲乱码一区二区免费版| 夫妻午夜视频| 最近手机中文字幕大全| 乱码一卡2卡4卡精品| 女人久久www免费人成看片| 深夜a级毛片| 校园人妻丝袜中文字幕| 欧美丝袜亚洲另类| 国产成人a区在线观看| 久久久久久久久久人人人人人人| 国产真实伦视频高清在线观看| 18禁在线播放成人免费| 少妇熟女欧美另类| 99久久精品一区二区三区| 日本免费在线观看一区| 最新中文字幕久久久久| 日韩强制内射视频| 成人性生交大片免费视频hd| 国产乱人偷精品视频| 成年版毛片免费区| 3wmmmm亚洲av在线观看| 91在线精品国自产拍蜜月| 免费看av在线观看网站| 亚洲性久久影院| 久久精品人妻少妇| 一级毛片 在线播放| 69人妻影院| 成人鲁丝片一二三区免费| 国产成人午夜福利电影在线观看| 免费无遮挡裸体视频| 26uuu在线亚洲综合色| 久久精品久久精品一区二区三区| 秋霞在线观看毛片| 日韩电影二区| 精品一区二区三区人妻视频| 777米奇影视久久| 男人爽女人下面视频在线观看| 国产在线一区二区三区精| 免费av不卡在线播放| 亚洲欧洲国产日韩| 日日摸夜夜添夜夜爱| 久久久久久久国产电影| 国内精品一区二区在线观看| 亚洲精品成人久久久久久| 超碰av人人做人人爽久久| 男人狂女人下面高潮的视频| 亚洲无线观看免费| 国产精品国产三级专区第一集| 成人特级av手机在线观看| 非洲黑人性xxxx精品又粗又长| 欧美潮喷喷水| 午夜福利在线在线| 国产亚洲最大av| 国产精品国产三级国产专区5o| 国产探花极品一区二区| 美女主播在线视频| 一区二区三区乱码不卡18| 日韩视频在线欧美| av网站免费在线观看视频 | 男人爽女人下面视频在线观看| 午夜免费观看性视频| 亚洲激情五月婷婷啪啪| 欧美另类一区| 久久人人爽人人片av| 色5月婷婷丁香| 看非洲黑人一级黄片| 91精品伊人久久大香线蕉| 高清毛片免费看| 美女脱内裤让男人舔精品视频| 亚洲av福利一区| 国产片特级美女逼逼视频| 国产欧美日韩精品一区二区| 亚洲精品自拍成人| 床上黄色一级片| 嘟嘟电影网在线观看| 久久鲁丝午夜福利片| 精品久久久久久久久亚洲| 亚洲精品影视一区二区三区av| 老女人水多毛片| 午夜福利网站1000一区二区三区| 床上黄色一级片| 少妇熟女欧美另类| 亚洲最大成人中文| 久久精品夜夜夜夜夜久久蜜豆| videossex国产| 成人综合一区亚洲| 精品人妻熟女av久视频| 国产精品蜜桃在线观看| 日韩国内少妇激情av| 久久久久精品久久久久真实原创| 久久精品久久精品一区二区三区| 午夜精品国产一区二区电影 | 别揉我奶头 嗯啊视频| 夜夜看夜夜爽夜夜摸| 少妇的逼好多水| 美女高潮的动态| 人妻制服诱惑在线中文字幕| 搡女人真爽免费视频火全软件| 18禁在线无遮挡免费观看视频| 老女人水多毛片| 午夜福利高清视频| 免费大片18禁| 亚洲精品乱码久久久v下载方式| 夫妻性生交免费视频一级片| 久久久成人免费电影| 亚洲精品国产av成人精品| 成人综合一区亚洲| 婷婷色麻豆天堂久久| 久久精品久久精品一区二区三区| 又爽又黄a免费视频| 男的添女的下面高潮视频| 国产亚洲最大av| 欧美成人精品欧美一级黄| 在线a可以看的网站| 国产成人免费观看mmmm| 国产亚洲av嫩草精品影院| 久久久久久久久中文| 亚洲欧美精品专区久久| 国产黄片视频在线免费观看| 高清在线视频一区二区三区| 三级国产精品欧美在线观看| 91久久精品国产一区二区三区| 天堂俺去俺来也www色官网 | 国产极品天堂在线| 亚洲aⅴ乱码一区二区在线播放| videos熟女内射| 亚洲精品自拍成人| 卡戴珊不雅视频在线播放| 天堂影院成人在线观看| 人妻夜夜爽99麻豆av| 九色成人免费人妻av| 欧美丝袜亚洲另类| 免费av观看视频| 精华霜和精华液先用哪个| 一本久久精品| 久久精品熟女亚洲av麻豆精品 | 亚洲国产欧美在线一区| 51国产日韩欧美| 伊人久久精品亚洲午夜| 日韩大片免费观看网站| 熟女人妻精品中文字幕| 80岁老熟妇乱子伦牲交| 免费观看a级毛片全部| 男女视频在线观看网站免费| 国产亚洲av片在线观看秒播厂 | 国产欧美另类精品又又久久亚洲欧美| av女优亚洲男人天堂| 国产亚洲精品av在线| 国产精品av视频在线免费观看| 日本一二三区视频观看| 18+在线观看网站| 亚洲欧洲国产日韩| 久久久久网色| 国产一区二区三区综合在线观看 | 成人漫画全彩无遮挡| 99久国产av精品| 国产av码专区亚洲av| 国产精品爽爽va在线观看网站| 中文字幕制服av| 国产免费一级a男人的天堂| 综合色av麻豆| 女的被弄到高潮叫床怎么办| 欧美三级亚洲精品| 欧美性猛交╳xxx乱大交人| 丝瓜视频免费看黄片| 日韩欧美 国产精品| 国产免费视频播放在线视频 | 久久精品国产自在天天线| 最近手机中文字幕大全| 久久久久久久亚洲中文字幕| 亚洲精品自拍成人| 国产免费视频播放在线视频 | xxx大片免费视频| 熟女人妻精品中文字幕| 日韩,欧美,国产一区二区三区| 婷婷色综合大香蕉| 99久久精品国产国产毛片| 日韩成人伦理影院| 国产欧美另类精品又又久久亚洲欧美| 美女高潮的动态| 亚洲国产高清在线一区二区三| 99九九线精品视频在线观看视频| 国产 一区 欧美 日韩| 亚洲自拍偷在线| 2021少妇久久久久久久久久久| 欧美+日韩+精品| 国产成人精品婷婷| 成人鲁丝片一二三区免费| 国产熟女欧美一区二区| 欧美性感艳星| 天堂网av新在线| 国产黄频视频在线观看| 插阴视频在线观看视频| av专区在线播放| 热99在线观看视频| 麻豆成人午夜福利视频| 久久久久国产网址| 亚洲av成人av| a级毛色黄片| 1000部很黄的大片| 嫩草影院新地址| 日韩电影二区| 少妇被粗大猛烈的视频| 久久久久精品性色| 在线a可以看的网站| 中国国产av一级| 夜夜爽夜夜爽视频| 搡女人真爽免费视频火全软件| 男人爽女人下面视频在线观看| 岛国毛片在线播放| 国产亚洲最大av| 国产成人a区在线观看| 国产一区二区三区av在线| 国产精品美女特级片免费视频播放器| 女人久久www免费人成看片| 最近视频中文字幕2019在线8| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美精品专区久久| 在现免费观看毛片| 国产毛片a区久久久久| 国产亚洲午夜精品一区二区久久 | 国产一级毛片在线| 男女边摸边吃奶| 97热精品久久久久久| 国产成人一区二区在线| 人人妻人人澡欧美一区二区| 99久国产av精品国产电影| 三级经典国产精品| av在线天堂中文字幕| 日日啪夜夜爽| 日韩欧美精品v在线| 三级毛片av免费| 身体一侧抽搐| 91久久精品国产一区二区成人| 欧美日本视频| 免费不卡的大黄色大毛片视频在线观看 | 日本一本二区三区精品| 黄片wwwwww| 青春草国产在线视频| 国内精品美女久久久久久| 午夜福利视频精品| 久久精品国产亚洲av涩爱| 国产精品三级大全| 精品久久久噜噜| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人二区视频| 日韩国内少妇激情av| 久久人人爽人人爽人人片va| 午夜久久久久精精品| 五月伊人婷婷丁香| 最近的中文字幕免费完整| 永久网站在线| 一区二区三区乱码不卡18| 十八禁网站网址无遮挡 | 偷拍熟女少妇极品色| 久久精品熟女亚洲av麻豆精品 | videossex国产| 国产精品.久久久| 69人妻影院| 99久久精品一区二区三区| 国产成人91sexporn| 亚洲精品456在线播放app| ponron亚洲| 国产精品久久久久久av不卡| 国产精品99久久久久久久久| 夫妻午夜视频| 亚洲在线自拍视频| 色综合亚洲欧美另类图片| 18禁动态无遮挡网站| 亚洲真实伦在线观看| 国产成人a∨麻豆精品| 久久精品国产亚洲av涩爱| 舔av片在线| 午夜免费男女啪啪视频观看| 有码 亚洲区| 国产探花在线观看一区二区| 免费不卡的大黄色大毛片视频在线观看 | 欧美高清成人免费视频www| 九色成人免费人妻av| 日本一本二区三区精品| 亚洲怡红院男人天堂| 亚洲av电影不卡..在线观看| 97人妻精品一区二区三区麻豆| 麻豆av噜噜一区二区三区| 中文天堂在线官网| 美女cb高潮喷水在线观看| 国产白丝娇喘喷水9色精品| 嫩草影院新地址| av在线亚洲专区| 久久久精品欧美日韩精品| 久久久欧美国产精品| 日本猛色少妇xxxxx猛交久久| 欧美成人a在线观看| 人人妻人人看人人澡| 久久久久久久午夜电影| 久久这里有精品视频免费| 精品午夜福利在线看| 观看免费一级毛片| 99久久中文字幕三级久久日本| 国产大屁股一区二区在线视频| 亚洲欧美成人综合另类久久久| 亚洲国产成人一精品久久久| 97人妻精品一区二区三区麻豆| 国产中年淑女户外野战色| 国内精品美女久久久久久| 亚洲四区av| 免费观看av网站的网址| 男插女下体视频免费在线播放| 欧美最新免费一区二区三区| 毛片一级片免费看久久久久| 岛国毛片在线播放| 91久久精品国产一区二区成人| 国产精品.久久久| 六月丁香七月| 人人妻人人澡人人爽人人夜夜 | 如何舔出高潮| 亚洲熟女精品中文字幕| 亚洲精品456在线播放app| 欧美日韩在线观看h| 国产中年淑女户外野战色| 午夜日本视频在线| 午夜激情福利司机影院| freevideosex欧美| 国产黄色视频一区二区在线观看| 久久精品人妻少妇| 少妇人妻精品综合一区二区| 国产精品一区二区性色av| 麻豆成人午夜福利视频| 亚洲精品aⅴ在线观看| 老司机影院成人| 边亲边吃奶的免费视频| 久久久久精品性色| 久久热精品热| 最近的中文字幕免费完整| 久久久午夜欧美精品| 日韩一区二区三区影片| 美女大奶头视频| 99热这里只有精品一区| 精品久久久久久电影网| 久久久a久久爽久久v久久| 亚洲精品色激情综合| 国产欧美日韩精品一区二区| 久久久精品94久久精品| 久久久久免费精品人妻一区二区| 偷拍熟女少妇极品色| 亚洲最大成人中文| 欧美激情在线99| 欧美最新免费一区二区三区| 最近手机中文字幕大全| 寂寞人妻少妇视频99o| 国产成人aa在线观看| 亚洲精品日韩在线中文字幕| 亚洲成人一二三区av| 最新中文字幕久久久久| 亚洲国产精品国产精品| 久久久久久久久久黄片| 亚洲在线观看片| 天堂√8在线中文| 亚洲精华国产精华液的使用体验| 青春草亚洲视频在线观看| 十八禁网站网址无遮挡 | 男人爽女人下面视频在线观看| 国产爱豆传媒在线观看| 青春草国产在线视频| 亚洲精品成人久久久久久| 直男gayav资源| 亚洲av福利一区| 午夜免费男女啪啪视频观看| av在线天堂中文字幕| 亚洲在线观看片| av一本久久久久| 不卡视频在线观看欧美| 91精品伊人久久大香线蕉| 免费观看av网站的网址| 国产精品一区www在线观看| 日韩av不卡免费在线播放| 久久久久久久久久成人| 国产黄色小视频在线观看| 国产毛片a区久久久久| 欧美成人a在线观看| 国产精品久久久久久久电影| 国产精品一区www在线观看| av在线观看视频网站免费| 亚洲av电影不卡..在线观看| 99热网站在线观看| 国产淫语在线视频| 国产精品日韩av在线免费观看| 久久亚洲国产成人精品v| 国产免费视频播放在线视频 | 欧美成人a在线观看| 身体一侧抽搐| 青春草国产在线视频| 又大又黄又爽视频免费| 最近最新中文字幕大全电影3|