• <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.

    亚洲精品在线观看二区| 满18在线观看网站| 丝袜人妻中文字幕| 人妻久久中文字幕网| 免费在线观看日本一区| 亚洲精品av麻豆狂野| 下体分泌物呈黄色| 无遮挡黄片免费观看| 91精品三级在线观看| 国产免费福利视频在线观看| 国产一区二区在线观看av| 国产aⅴ精品一区二区三区波| 两性午夜刺激爽爽歪歪视频在线观看 | 老汉色av国产亚洲站长工具| 国产色视频综合| 法律面前人人平等表现在哪些方面| 欧美乱码精品一区二区三区| 亚洲成人免费av在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 欧美国产精品一级二级三级| 国产精品99久久99久久久不卡| tube8黄色片| 久久久国产精品麻豆| 黄片大片在线免费观看| 久久午夜综合久久蜜桃| 精品国产亚洲在线| 国产精品二区激情视频| 嫁个100分男人电影在线观看| 美女扒开内裤让男人捅视频| av电影中文网址| 少妇粗大呻吟视频| 老汉色av国产亚洲站长工具| 三级毛片av免费| 啦啦啦中文免费视频观看日本| 建设人人有责人人尽责人人享有的| 欧美黑人精品巨大| 亚洲一区二区三区欧美精品| 丰满迷人的少妇在线观看| 777久久人妻少妇嫩草av网站| 国产日韩欧美在线精品| 欧美精品一区二区大全| 亚洲欧洲精品一区二区精品久久久| 在线 av 中文字幕| 人妻 亚洲 视频| 人成视频在线观看免费观看| 伦理电影免费视频| 国产精品久久久久久精品古装| 亚洲av美国av| 国产极品粉嫩免费观看在线| 真人做人爱边吃奶动态| 91成人精品电影| 成在线人永久免费视频| 亚洲国产欧美在线一区| 变态另类成人亚洲欧美熟女 | 亚洲 欧美一区二区三区| 国产无遮挡羞羞视频在线观看| 国产伦理片在线播放av一区| 在线观看免费日韩欧美大片| 黄色a级毛片大全视频| 麻豆乱淫一区二区| 老汉色av国产亚洲站长工具| 日本一区二区免费在线视频| 91精品国产国语对白视频| 12—13女人毛片做爰片一| 99re6热这里在线精品视频| 在线观看免费日韩欧美大片| 日日爽夜夜爽网站| 亚洲精品在线观看二区| 在线看a的网站| av又黄又爽大尺度在线免费看| 少妇猛男粗大的猛烈进出视频| 久久婷婷成人综合色麻豆| 国产精品国产av在线观看| 最近最新免费中文字幕在线| 男女高潮啪啪啪动态图| 精品一区二区三区视频在线观看免费 | 国产人伦9x9x在线观看| 麻豆成人av在线观看| 一进一出好大好爽视频| 如日韩欧美国产精品一区二区三区| 黑人欧美特级aaaaaa片| 久久久久久久大尺度免费视频| 成年人免费黄色播放视频| 一夜夜www| 国产单亲对白刺激| 91精品国产国语对白视频| 色综合欧美亚洲国产小说| 我要看黄色一级片免费的| 欧美精品高潮呻吟av久久| 香蕉丝袜av| 三级毛片av免费| 热re99久久国产66热| 国产精品影院久久| cao死你这个sao货| 亚洲欧洲精品一区二区精品久久久| 国产亚洲一区二区精品| 深夜精品福利| 后天国语完整版免费观看| 国产成人av激情在线播放| 欧美亚洲日本最大视频资源| 久久狼人影院| 啦啦啦 在线观看视频| 1024香蕉在线观看| 国产精品免费视频内射| 亚洲全国av大片| 日韩欧美一区视频在线观看| 国产三级黄色录像| 电影成人av| 国产精品免费一区二区三区在线 | 蜜桃在线观看..| 精品国产乱码久久久久久男人| 精品国产乱码久久久久久小说| 免费av中文字幕在线| 亚洲一区中文字幕在线| 丁香欧美五月| 一个人免费在线观看的高清视频| 日本av免费视频播放| 日韩中文字幕视频在线看片| 亚洲 欧美一区二区三区| 水蜜桃什么品种好| 韩国精品一区二区三区| 少妇被粗大的猛进出69影院| av有码第一页| 麻豆av在线久日| 中文欧美无线码| 在线观看免费午夜福利视频| 成年女人毛片免费观看观看9 | 国产极品粉嫩免费观看在线| 人妻久久中文字幕网| 欧美精品av麻豆av| 大片免费播放器 马上看| 侵犯人妻中文字幕一二三四区| 999精品在线视频| av在线播放免费不卡| 国产三级黄色录像| 亚洲黑人精品在线| 免费不卡黄色视频| 精品卡一卡二卡四卡免费| 国产深夜福利视频在线观看| 久久人妻熟女aⅴ| 操美女的视频在线观看| 黄色视频在线播放观看不卡| 亚洲精品美女久久久久99蜜臀| 久久这里只有精品19| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜激情av网站| 国产伦理片在线播放av一区| 在线av久久热| 一边摸一边做爽爽视频免费| 免费高清在线观看日韩| www.精华液| 丰满饥渴人妻一区二区三| 黄色视频不卡| 一本色道久久久久久精品综合| 国产精品国产高清国产av | 亚洲色图 男人天堂 中文字幕| 黄色毛片三级朝国网站| 国产亚洲精品第一综合不卡| 高清毛片免费观看视频网站 | 一进一出抽搐动态| 欧美日韩精品网址| 久久久水蜜桃国产精品网| 一区二区av电影网| 午夜日韩欧美国产| 久久久久久久国产电影| 国产精品免费一区二区三区在线 | 中文字幕最新亚洲高清| 色老头精品视频在线观看| 久久人人97超碰香蕉20202| 国产精品.久久久| 操出白浆在线播放| 啦啦啦中文免费视频观看日本| 亚洲精品国产区一区二| 亚洲成av片中文字幕在线观看| 久久亚洲真实| √禁漫天堂资源中文www| 国产精品久久久久成人av| 女人被躁到高潮嗷嗷叫费观| 宅男免费午夜| 黄色视频,在线免费观看| 天天操日日干夜夜撸| 99热国产这里只有精品6| 欧美久久黑人一区二区| 丝袜喷水一区| 十八禁人妻一区二区| 免费一级毛片在线播放高清视频 | 国产亚洲欧美精品永久| 18禁观看日本| 麻豆av在线久日| 亚洲 欧美一区二区三区| 亚洲成人免费av在线播放| 深夜精品福利| 日韩人妻精品一区2区三区| 欧美精品一区二区大全| 欧美一级毛片孕妇| av不卡在线播放| 亚洲美女黄片视频| 日韩欧美免费精品| 别揉我奶头~嗯~啊~动态视频| 成人18禁在线播放| 曰老女人黄片| 一区二区三区激情视频| 97在线人人人人妻| 99国产综合亚洲精品| 电影成人av| 色尼玛亚洲综合影院| 在线 av 中文字幕| 亚洲第一av免费看| 欧美在线一区亚洲| 高清av免费在线| www日本在线高清视频| 久久久久久人人人人人| 欧美黑人精品巨大| 亚洲一码二码三码区别大吗| 首页视频小说图片口味搜索| 大片免费播放器 马上看| 精品一区二区三区视频在线观看免费 | 九色亚洲精品在线播放| 亚洲国产欧美日韩在线播放| 乱人伦中国视频| 欧美黑人欧美精品刺激| 亚洲一区二区三区欧美精品| 欧美日韩亚洲国产一区二区在线观看 | 高清黄色对白视频在线免费看| 无人区码免费观看不卡 | 国产欧美日韩一区二区三| 纵有疾风起免费观看全集完整版| 国产成人av教育| av国产精品久久久久影院| 女同久久另类99精品国产91| 一区二区日韩欧美中文字幕| 中文字幕人妻丝袜制服| 精品熟女少妇八av免费久了| 脱女人内裤的视频| 久久精品亚洲熟妇少妇任你| 午夜久久久在线观看| 久久国产精品影院| 91麻豆精品激情在线观看国产 | 999精品在线视频| 久久中文字幕一级| 久久久久精品国产欧美久久久| 亚洲精品中文字幕在线视频| 亚洲成人国产一区在线观看| 最新的欧美精品一区二区| 两个人看的免费小视频| 欧美 日韩 精品 国产| 成人亚洲精品一区在线观看| 亚洲免费av在线视频| 久久精品亚洲精品国产色婷小说| 国产在线观看jvid| 亚洲av第一区精品v没综合| 法律面前人人平等表现在哪些方面| 婷婷成人精品国产| 丝袜美腿诱惑在线| 亚洲avbb在线观看| 另类精品久久| 精品国产一区二区三区久久久樱花| 美女扒开内裤让男人捅视频| 亚洲熟女毛片儿| 亚洲精品粉嫩美女一区| 一本—道久久a久久精品蜜桃钙片| 在线十欧美十亚洲十日本专区| 欧美午夜高清在线| 99riav亚洲国产免费| 精品高清国产在线一区| 亚洲欧美色中文字幕在线| 男女边摸边吃奶| 日日夜夜操网爽| 国产成人系列免费观看| 久久天躁狠狠躁夜夜2o2o| 久久久久久久久免费视频了| 新久久久久国产一级毛片| 人成视频在线观看免费观看| 人人澡人人妻人| 精品人妻1区二区| 亚洲七黄色美女视频| 久久影院123| 天天躁夜夜躁狠狠躁躁| 免费女性裸体啪啪无遮挡网站| 国产精品久久电影中文字幕 | 看免费av毛片| 亚洲精品在线观看二区| 18在线观看网站| 国产免费视频播放在线视频| 丁香欧美五月| 亚洲天堂av无毛| 大片免费播放器 马上看| 精品一区二区三卡| 超色免费av| 91九色精品人成在线观看| 欧美变态另类bdsm刘玥| 中文字幕人妻熟女乱码| www日本在线高清视频| 亚洲男人天堂网一区| 精品久久久久久电影网| 不卡av一区二区三区| 考比视频在线观看| 日韩有码中文字幕| 欧美国产精品va在线观看不卡| 搡老岳熟女国产| 国产精品 国内视频| 亚洲av成人一区二区三| 91成人精品电影| 最新的欧美精品一区二区| 搡老熟女国产l中国老女人| 国产在线观看jvid| 色精品久久人妻99蜜桃| 午夜成年电影在线免费观看| 一二三四在线观看免费中文在| 成年动漫av网址| 一边摸一边抽搐一进一小说 | 伦理电影免费视频| 国产高清视频在线播放一区| 欧美激情极品国产一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久国产亚洲av麻豆专区| 天天躁夜夜躁狠狠躁躁| 久久人妻熟女aⅴ| 女警被强在线播放| 日韩中文字幕视频在线看片| 色婷婷久久久亚洲欧美| 极品少妇高潮喷水抽搐| 亚洲成av片中文字幕在线观看| 欧美在线一区亚洲| 美国免费a级毛片| 国产成人精品无人区| 亚洲黑人精品在线| 久久天躁狠狠躁夜夜2o2o| 69av精品久久久久久 | 亚洲欧美日韩高清在线视频 | 国产成人精品无人区| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美在线一区亚洲| 午夜精品国产一区二区电影| 男人操女人黄网站| 美女视频免费永久观看网站| 极品人妻少妇av视频| 少妇精品久久久久久久| 中文字幕色久视频| 亚洲美女黄片视频| 看免费av毛片| 日本一区二区免费在线视频| 国产精品国产av在线观看| 久久久久网色| 99久久精品国产亚洲精品| 男女床上黄色一级片免费看| 欧美国产精品一级二级三级| 久久中文字幕一级| 三上悠亚av全集在线观看| av网站在线播放免费| 国产av精品麻豆| 午夜福利视频精品| 久久久久精品人妻al黑| 免费日韩欧美在线观看| 亚洲国产欧美一区二区综合| 色在线成人网| 欧美大码av| 18在线观看网站| 国产精品美女特级片免费视频播放器 | 精品少妇黑人巨大在线播放| 久久久久精品人妻al黑| 亚洲,欧美精品.| 三上悠亚av全集在线观看| 视频区欧美日本亚洲| 视频区图区小说| 亚洲av成人不卡在线观看播放网| 亚洲人成电影观看| 久久精品亚洲精品国产色婷小说| videosex国产| www.999成人在线观看| 青青草视频在线视频观看| 十八禁高潮呻吟视频| 亚洲人成伊人成综合网2020| 欧美日韩亚洲综合一区二区三区_| 亚洲第一av免费看| 婷婷成人精品国产| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕最新亚洲高清| 91精品国产国语对白视频| 日本vs欧美在线观看视频| 亚洲精品国产一区二区精华液| 在线播放国产精品三级| 久久精品亚洲av国产电影网| 国产一区二区 视频在线| 亚洲午夜理论影院| 老熟妇仑乱视频hdxx| 1024香蕉在线观看| 91成人精品电影| 国产欧美日韩一区二区三区在线| 99久久人妻综合| 少妇被粗大的猛进出69影院| 久久久国产精品麻豆| 又大又爽又粗| 国产精品久久电影中文字幕 | 1024视频免费在线观看| 中文欧美无线码| 无人区码免费观看不卡 | 国产精品九九99| 纵有疾风起免费观看全集完整版| 欧美精品一区二区免费开放| 夫妻午夜视频| 最近最新中文字幕大全电影3 | 亚洲第一欧美日韩一区二区三区 | 亚洲av成人不卡在线观看播放网| 黄片播放在线免费| 视频在线观看一区二区三区| 啦啦啦 在线观看视频| 亚洲国产看品久久| av视频免费观看在线观看| 久久av网站| 亚洲一区二区三区欧美精品| av一本久久久久| 亚洲午夜理论影院| 变态另类成人亚洲欧美熟女 | 91字幕亚洲| 久久久精品免费免费高清| 日韩欧美一区二区三区在线观看 | 999精品在线视频| 国产91精品成人一区二区三区 | 日韩欧美免费精品| 男女边摸边吃奶| 成人国产av品久久久| 久久精品成人免费网站| cao死你这个sao货| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产欧美一区二区综合| 男女之事视频高清在线观看| 一区二区日韩欧美中文字幕| 在线观看免费视频日本深夜| 香蕉久久夜色| 黑人欧美特级aaaaaa片| 夜夜夜夜夜久久久久| 丝袜人妻中文字幕| 男女免费视频国产| 90打野战视频偷拍视频| 在线观看免费视频网站a站| 久久精品亚洲熟妇少妇任你| 欧美一级毛片孕妇| 99久久99久久久精品蜜桃| 精品亚洲成a人片在线观看| av有码第一页| 十八禁人妻一区二区| 久久人妻熟女aⅴ| 欧美黑人欧美精品刺激| 精品国产一区二区三区久久久樱花| 精品国产乱码久久久久久男人| 亚洲国产av影院在线观看| 黑人欧美特级aaaaaa片| 18在线观看网站| 色视频在线一区二区三区| 一级黄色大片毛片| 免费黄频网站在线观看国产| 黑丝袜美女国产一区| 亚洲国产av新网站| 啦啦啦视频在线资源免费观看| 亚洲第一av免费看| 亚洲欧美日韩另类电影网站| 国产在线一区二区三区精| 久久精品人人爽人人爽视色| 久久免费观看电影| 黄色a级毛片大全视频| 日本欧美视频一区| 人成视频在线观看免费观看| 午夜成年电影在线免费观看| 成人特级黄色片久久久久久久 | 久久久久久久久免费视频了| 可以免费在线观看a视频的电影网站| 自拍欧美九色日韩亚洲蝌蚪91| 大码成人一级视频| 午夜日韩欧美国产| 亚洲伊人色综图| 欧美日韩亚洲综合一区二区三区_| 国产成人影院久久av| xxxhd国产人妻xxx| 欧美日韩黄片免| 一进一出好大好爽视频| 黄色视频,在线免费观看| 亚洲天堂av无毛| 久久久国产一区二区| 一本久久精品| 超碰97精品在线观看| 91麻豆精品激情在线观看国产 | 日本wwww免费看| 久久精品亚洲精品国产色婷小说| 欧美精品高潮呻吟av久久| 久久国产精品男人的天堂亚洲| 精品福利永久在线观看| 制服人妻中文乱码| 91字幕亚洲| 老汉色av国产亚洲站长工具| 天天操日日干夜夜撸| av福利片在线| 亚洲色图综合在线观看| 亚洲免费av在线视频| 日韩免费高清中文字幕av| 大型黄色视频在线免费观看| www.自偷自拍.com| www.999成人在线观看| 肉色欧美久久久久久久蜜桃| 欧美日韩成人在线一区二区| 精品少妇久久久久久888优播| 黄色a级毛片大全视频| 久久人人爽av亚洲精品天堂| 老司机福利观看| 国产精品亚洲一级av第二区| 大香蕉久久网| 老鸭窝网址在线观看| 国产免费现黄频在线看| 国产一区二区三区在线臀色熟女 | 18禁黄网站禁片午夜丰满| 亚洲精品国产色婷婷电影| 亚洲自偷自拍图片 自拍| 国产野战对白在线观看| 飞空精品影院首页| 狠狠狠狠99中文字幕| 国产色视频综合| 国产男女超爽视频在线观看| 亚洲第一欧美日韩一区二区三区 | 超碰97精品在线观看| 精品人妻熟女毛片av久久网站| 亚洲性夜色夜夜综合| 多毛熟女@视频| 国产精品免费视频内射| 久久久欧美国产精品| 最新的欧美精品一区二区| 一级,二级,三级黄色视频| 日韩视频一区二区在线观看| 久久香蕉激情| 国产aⅴ精品一区二区三区波| 一本综合久久免费| 国产xxxxx性猛交| 精品国产一区二区三区久久久樱花| 狂野欧美激情性xxxx| 飞空精品影院首页| 精品国产一区二区三区四区第35| 麻豆国产av国片精品| 成人18禁高潮啪啪吃奶动态图| 黄色成人免费大全| a级毛片黄视频| 美女高潮到喷水免费观看| 国产精品欧美亚洲77777| 精品人妻1区二区| 国产亚洲精品久久久久5区| 国产激情久久老熟女| 999久久久精品免费观看国产| 女性生殖器流出的白浆| 一级毛片电影观看| 香蕉久久夜色| 国产深夜福利视频在线观看| 久久婷婷成人综合色麻豆| 黄色视频不卡| 肉色欧美久久久久久久蜜桃| 自拍欧美九色日韩亚洲蝌蚪91| 人人妻人人添人人爽欧美一区卜| 在线观看舔阴道视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲男人天堂网一区| 99国产精品一区二区蜜桃av | 欧美日韩福利视频一区二区| 亚洲av欧美aⅴ国产| 亚洲九九香蕉| videosex国产| 色尼玛亚洲综合影院| 丝袜在线中文字幕| 亚洲视频免费观看视频| 无人区码免费观看不卡 | 亚洲黑人精品在线| 黄片播放在线免费| 老司机亚洲免费影院| 色94色欧美一区二区| 19禁男女啪啪无遮挡网站| 视频在线观看一区二区三区| 天堂8中文在线网| 成年人黄色毛片网站| 丝袜美足系列| 脱女人内裤的视频| 国产精品影院久久| 欧美日韩亚洲综合一区二区三区_| 国产男女内射视频| 老司机福利观看| 国产亚洲欧美精品永久| 天堂俺去俺来也www色官网| 国产又爽黄色视频| 国产在线视频一区二区| 日韩有码中文字幕| 亚洲精品自拍成人| 18禁黄网站禁片午夜丰满| 啦啦啦视频在线资源免费观看| 97人妻天天添夜夜摸| 欧美日韩黄片免| av视频免费观看在线观看| 如日韩欧美国产精品一区二区三区| 啦啦啦在线免费观看视频4| 日韩欧美国产一区二区入口| 成年人免费黄色播放视频| 国产有黄有色有爽视频| 69av精品久久久久久 | 日韩中文字幕欧美一区二区| 91成年电影在线观看| 亚洲人成伊人成综合网2020| 国产亚洲精品一区二区www | 国产av精品麻豆| 亚洲成av片中文字幕在线观看| 精品国产乱码久久久久久男人| 日本wwww免费看| 国产高清videossex| 黄色成人免费大全| 国产精品久久久久久精品电影小说| 久久婷婷成人综合色麻豆| 亚洲av第一区精品v没综合| 国产片内射在线| 99久久99久久久精品蜜桃| 精品少妇内射三级| 一级a爱视频在线免费观看| 日韩三级视频一区二区三区|