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

    Design of crashworthy attenuator structures as a part of vehicle safety against impact: Application of waste aluminum can-based material

    2021-07-30 13:34:56LksmnWidiPrsetyAdityRioProwoUidillhUidillhIwnIstntoNurAzmhBintiNordin

    Lksmn Widi Prsety , Adity Rio Pr owo , ?, Uidillh Uidillh , Iwn Istnto ,Nur Azmh Binti Nordin

    a Department of Mechanical Engineering, Universitas Sebelas Maret, Surakarta 57126, Indonesia

    b Agency for the Assessment and Application of Technology (BPPT), Jakarta 10340, Indonesia

    c Malaysia-Japan International Institute of Technology - Universiti Teknologi Malaysia, Kuala Lumpur 54100, Malaysia

    Keywords: Attenuator structure Waste can Vehicle safety Crashworthiness analysis Impact phenomenon

    A B S T R A C T The impact attenuator is an essential system in both race cars and urban vehicles. The structure of an impact attenuator serves as a safety barrier between the impacted surface and the driver in an accident.Attenuator materials tend to have a high price; thus, alternative materials were explored in the current work, i.e., used cans from food and beverage containers. The study deployed a nonlinear finite element algorithm to calculate a series of impacts on the attenuator structures. The thickness of the cans and velocity of the impact were considered as the main parameters. Analysis results concluded that the attenuator’s average energy was 160 0 0 J for a can thickness of 1 mm. This value is more than two times the 0.5 mm thick used cans. The attenuator’s new design was then matched with an attenuator regulation,and the results surpassed the standard value of 7350 J.

    The Student Formula is a competition formed by the Society of Automotive Engineers, which aims to challenge students from each university to design and conceptualize small-model formula racing cars. In making this car, protection issues are a significant demand that is very well regulated in the Formula Society of Automotive Engineers [1-3] . Recently, there has been an increasing interest in vehicle safety, which has led to increased research and analysis on the vehicle accident response from experimental trials and simulations [4] . The development of active and passive safety systems that work for passenger safety is also being conducted in the industrial world [5] .

    The competition needs every participant to the concept and builds a racecar similar to in the same company. In creating a vehicle, remember to concentrate on the protective factors regulated within the race rules (Formula Society of Automotive Engineers -FSAE Rules). In these regulations, each team must make a safety device called an impact attenuator. Based on these rules, the impact attenuator is declared fit for use to absorb at least 7350 J of energy. The lightweight design will improve the vehicle acceleration performance and fuel economy [6] . An attenuator must be specially made to absorb the collision’s energy—namely the vehicle’s kinetic energy—and reduce the driver’s deceleration acting [7] .The technical requirements of an attenuator are that:

    1.The attenuator is set up at the front of a system vehicle with a whole mass of 300 kg, the average and highest deceleration have to be no more than 20 g and 40 g, respectively, while the automobile is traveling at a constant velocity of 7 mm/ms (gstands for gravity acceleration 9.8 m/s2).

    2.The entire attenuator absorbed ought to be, as a minimum,7350 J. There also are geometric limitations to the design. The minimum length is 200 mm with a minimal height of 100 mm,200 mm width and maximum distance to the front bulkhead is 200 mm [1] .

    3.In impact/collision situations, the attenuator is disallowed to penetrate the front bulkhead while it cannot be a part of bodywork structure. It ought to be securely connected to the front bulkhead [1] .

    4.When off-center and off-axis impact/collisions occur, there ought to be sufficient load paths for transverse and vertical loads [8,9] . Furthermore, the policies suggest that the antiintrusion plates has to be inserted to all racing motors and automobile. These may be manufactured from 5 mm robust steel or 4 mm solid aluminum.

    5.Vital components, e.g., the critical hoop, battery, and hydraulic reservoir has to be placed behind the bulkhead as they are not allowed on the impact attenuator sector [1] .

    In designing an attenuator, the elements that may be considered for fabric dedication include the cost, reliability, weight, and availability [10] . Attenuators are normally manufactured from honeycomb, aluminum, Kevlar, Nomex, aluminum foil, or a combination of the mentioned materials [11-13] . Apart from those materials, different materials, including used cans, are typically manufactured from aluminum and different metal alloys.

    Cans are one form of steel waste that is regularly reused and,barring that, may pollute nature. If rust is exposed to water after entering the soil, it will interfere with the soil fertility [14] .Another trouble is the developing use of aluminum for meals and gentle drink packaging [15] . Ordinary aluminum is mixed for added mechanical properties and strength. Cans are composed approximately 92–99% of aluminum, and other minors, e.g., copper,magnesium, silica, manganese [16] .

    The purpose of this study is to design an impact attenuator as development effort of vehicle safety. Consideration to reduce the pollution caused by used cans, the waste can-based material is used as material for the proposed design. The design will be analyzed to quantify its crashworthiness against impact/collision. Nonlinear finite element method is used to solve several impact cases in this study.

    An impact attenuator that functions as a safety system capable of blocking the surface that is impacted by a driver is a fundamental structure in a formula racing car [11] . This safety system is designed to absorb kinetic energy at the time of impact,which is converted into a deformation that maintains a reasonably low force level and uniform deformation [17,18] . As the deformation occurs unevenly, it will cause a driver to experience a G force.Energy absorption occurs in the form of extensively destroying the structure [19] . Transverse and vertical load paths are provided for the impact attenuator if the impact occurs outside the center axis[20] . The specific energy absorbed and the total energy absorbed can be used to assess the viability of the impact attenuator. The rate of energy dissipation concentrated in a narrow zone is a characteristic of impact attenuators, whereas the structural reminders experience rigid motion [21,22] .

    The FSAE rules state that there should be an Anti-Intrusion Plate and Impact Attenuator in the front of the Front Bulkhead,with the Anti-Intrusion Plate in the front of the Impact Attenuator. The Anti-Intrusion plate should be 1.5 mm (0.060 in.) thick of a strong metal material, and 4.0 mm (0.157 in.) thick for aluminum.The Anti-Intrusion Plate should be connected to the Front Seal securely, and for the Anti-Intrusion Plate’s outer profile for the weld joint, the face should be at the centerline of the Front Baffle tube on all sides. The shape should match the Front Bulkhead’s outer dimensions around the rims for bolted joints [1] .

    There also are provisions for effect absorbers in which the scale of the Impact Absorbers are at the least 200 mm (7.8 in.) long,with the duration oriented alongside the front/rear axle of the frame. A minimal peak of 100 mm (3.9 in.) and width of 200 mm(7.8 in.) are required for a minimum distance of 200 mm (7.8 in.)in the front of the Front Bulkhead. Finally, they should be attached securely to the Anti-Intrusion Plates or to the Front Seals. It is likewise assumed that the surprise absorber is established on the front of the car with a total mass of 300 kg (661 lbs.), and hitting a concrete effect barrier with a sufficient velocity of 7.0 m/s (23.0 ft/s),should gradually slow down the car to a pace not exceeding a mean of 20 g and a maximum of 40 g. The strength absorbed on this occasion should meet or exceed 7350 J [1] .

    Dynamic testing (sledge, pendulum, drop tower, etc.) of the Impact Reducer can be completed in a committed testing facility.While these centers can be a part of a university, they should be supervised by expert personnel or University faculty. Teams are not allowed to construct their own dynamic check kits. Whereas quasistatic testing may be completed through companies using their University centers/equipment, the groups are recommended to exercise warning whilst undertaking testing [1] .

    To maintain fuel and enhance the acceleration performance,the Impact Attenuator should be designed with light-weight substances [6] . Typically aluminum, honeycomb, Nomex, carbon composites, or a mixture of those substances are the substances used to make Impact Attenuators to provide the most safely for drivers as they are robust and tend to be light. Tube and plate types, honeycomb systems, trapezoidal truncated shapes, and sandwich systems are examples of effect-soaking systems [11-13] .

    The first impact and failure phenomenon concerns occupant protection and particularly pedestrian protection, which has become increasingly important over the years. The reason is inadequate pedestrian safety against severe injuries during a collision with vehicles. To reduce injuries and their severity, design changes in the front structure are required to absorb energy and minimize the force and pedestrians’ acceleration. Therefore, a growing number of energy absorption structures are often designed with different materials, such as highly porous and energy-absorbing foams with plastic deformation behavior, such as egg box structures. The combination of these parts interact under accident loads, for example, the bumper system. For the layout of the new protection system, the deformation behavior and force levels of these "chain interactions" must be predicted [23] .

    We concluded that, although the individual modules of materials of different stiffness were excellent to simulate, the adjustment of the parameters of the selected material, as well as the adjustment of the parameters for the calculation of contacts, proved necessary for the interaction—a chain with two or four components to produce a good fit between the simulation and test. The examination approach just described shows its respective limitations for simulation capability, which implies a tremendous importance for interpreting the calculation results. In comparison, there are currently opportunities for more detailed component modeling and new or enlarged material models; testing to verify the simulation results is still necessary to achieve improved prediction accuracy for complex and interacting systems.

    The restitution coefficient is a variable that considerably impacts the output fee received to use item collision evaluations.Typically, mathematical–mechanical fashions count on cognition of the numerical values of the restitution coefficients. In exercise, it is predicted that the restitution coefficient’s numerical fee is decided primarily based on experience. This exercise is mainly used within the case of inhomogeneous bodies of complicated shape. One realistic instance is automobile effect evaluation [24] .

    The methodological elimination has to be suitable for practical implementation. The reason for this technique is to lessen mistakes within the output of the automobile effect evaluation system.The benefit of this technique is its clean and rapid application. The damping ratio utilized in passenger cars is generally 40% to 60%better the issue of leap damping. Motorcycle dampers are characterized with the aid of using more significant full-size asymmetry.Asymmetry may also rely on drag working conditions (compression fee or extension).

    Occasionally, there may be an expressive asymmetry within the dampers positioned at the rear axle of a passenger automobile. This phenomenon is uncommon in racing cars and motorbikes. One viable clarification is that higher effect absorption even with decreasing consolation presents a motive force with a more elevated "trip feel." The research accomplished up to now has not yielded a reasonable cause of the presence of asymmetrical damper. The results of theoretical studies which state a requirement that the dampers have to be symmetrical, now no longer offers a legitimate reason. This changed with the discovery that damper asymmetry changed as an instantaneous result of its construction, i.e., it no longer considered due to outside influences[25] .

    The combined effects of this consider that every person is almost the same. Still, while choosing the most valuable part, it is essential to remember that the automobile’s favored conduct is ensured inside a wide variety of velocities and street traits (bumps,potholes, etc.).

    The distinction in amplitude from vertical acceleration and the frequency reached while an automobile travels through a street with minimal unevenness is invisible. To recognize the phenomenon of asymmetry, it is vital to look at the development of the damper itself and the impact of choosing the damping traits on the improvement and manufacturing charges of the drag. In practice, there may be a semi-lively suspension machine that can alternate the damper’s traits to acquire an excessive stage of consolation even maintaining the stableness of the automobile with unique riding conditions. The damper traits may be modified by adjusting the magnetorheological or electrorheological fluid viscosity within the damper, which offers a better or decreasing damping ratio value.

    Impact attenuator collision was considered and classified into a nonlinear dynamic phenomenon that desires an appropriate calculation approach to remedy its almost-infinite scenarios [26-29] .These assumptions are described in the shape of quick finite detail set of rules as provided in Eqs. (1) to (5) . Nonlinear is embedded within the evaluation as the essential assumptions of linear evaluation are violated because the fabric and shape exceed their yield limit. In this state, each component revels in the massive volume of deformation. On the other hand, the dynamic function is considered, as the calculation outcomes will extrude gradually depending on the time. To remedy the phenomena with those criteria, the specific approach is favored because numerous assumptions are advantageous:

    1.Inversion at the stiffened matrix is not always required as all nonlinear parameters (which includes contact) are protected within the inner pressure vector.

    2.A small time step is required to hold the balance limit, which is in shape with a small length of time procedure of the effect attenuator collision (commonly much less than 1 s.).

    3.Simple inversion is desired because the challenge is a lumped mass matrix.

    where {at} is the acceleration; {Fext} is the applied external and body force vector; {Fint} is the internal force vector;Fhgis the hourglass resistance force;Fcontis the contact force;BTis the form identical to the linear discrete strain–displacement matrix;Ωis the solid volume;σnis the internal stress; {vt} is the velocity at timet; {ut} is the displacement at timet; {x0} is the initial geometry;{xt} is the updated geometry at timet; and Δttis the difference in time of at timetcompared to the initial/selected condition.

    Based at the 2019 FSAE rules, the attenuator cannot be a part of the bodywork shape; however, the alternative is to be mounted securely and immediately into the front bulkhead so that once a collision occurs, the tool needs to no longer penetrate the front bulkhead.

    The impact attenuator includes three parts designed to use the thin-walled concept. The first is used cans, then rigid plates, and finally anti-intrusion plates. The used cans are set up in two levels. The early-stage includes six cans, and the higher step includes five cans. The used cans have a diameter of 72.4 mm with a thickness of two mm and a peak of 113 mm. A rigid plate has a 250 ×250 mm length with a thickness of 2 mm, which features a holder for used cans on the pinnacle stage. The closing is the anti-intrusion plate, which has a length of 325.4 ×305.4 mm. The general peak of the effect attenuator is 232 mm, as shown in Fig. 1 .

    The applied material on the used cans is aluminum alloy 3004.The mechanical properties of these materials are shown in Table 1 .In the finite element calculation, it is idealized as piecewise linear plasticity with failure. This material is assumed to be a deformable component. The correlation of the strain with increasing stress is presented in Fig. 2 .

    Table 1 Mechanical properties of waste food cans [30] .

    Fig. 1. Geometric model of the impact attenuator (a) isometric view; (b) left side view.

    Fig. 2. Correlation of the stress–strain for the deformable component.

    On the other hand, Table 2 presents the mechanical properties used in rigid parts, the first is the rigid wall as a hammer, the rigid plate, and the last is the anti-intrusion plate as a buffer.

    Table 2 Mechanical properties of rigid parts [31] .

    Documentation of an impact attenuator ready for compression testing using a universal testing machine (UTM) apparatus is presented in Fig. 3 . In this experimental test, the Impact Attenuator was made with the first step, namely the front bulkhead. The front bulkhead was made of steel pipe, which was cut and then welded.After that, construction was continued by cutting the basic aluminum plate with a thickness of 4 mm as the Anti-Intrusion Plate,the second aluminum plate with a thickness of 2 mm was for separating the upper and lower stack cans, and the third was an aluminum plate with a thickness of 2 mm. Then, the cans were arranged with five cans at the top and six cans at the bottom. The cans were glued so that they did not slide when pressed. Finally,the press test was carried out using a UTM tool.

    Fig. 3. Impact attenuator for the experimental tests.

    Fig. 4. Impact attenuator boundary condition.

    Fig. 5. Finite model of the impact attenuator (a) isometric view; (b) left side view.

    Fig. 6. (a) Crushing force A1 and A2; and (b) Crushing force A3 and A4.

    The set up in the compression test was applied and adjusted to the simulation test. For the first, the impact attenuator was assumed to be a shell element type consisting of a can and a rigid plate. For the boundary condition ( Fig. 4 ), the pounder plate and the separation plate between the top and bottom cans array were fully motion-restricted (ALL DOF) unless the displacement in z (vertical) direction was specified as ?200 mm. The anti-intrusion plate components were considered as non-translating and fully rotating (ALL DOF).

    Fig. 7. (a) Internal Energy A1 and A2; and (b) Internal Energy A3 and A4.

    For the used cans ( Fig. 5 ), all quadrilateral factors are modelled with a length of 5 mm. Determining the form of fabric for meal cans makes use of piecewise linear plasticity and all inflexible elements utilizing the type of severe material. This fabric choice was so that, after the collision occurs, the can will alternate in shape.The hard partitions will push the can under with a compressive pressure transmitted from the can at the top. In the input material properties for the used cans, the effective plastic strain curve shown in Fig. 2 is inputted into the define curve feature by loading the.csv data. The curve can be used as the material type featuring piecewise linear plasticity.

    Fig. 8. (a) Kinetic energy A1 and A2; and (b) Kinetic energy A3 and A4.

    Fig. 9. (a) Total energy A1 and A2; and (b) Total energy A3 and A4.

    To adjust the simulation test so that it has the same conditions as the quasi-static test, a prescribed motion rigid feature was used by defining a displacement curve, where this curve will be applied to the colliding part. This displacement curve was adjusted so that they-axis andx-axis were equal to 200 to obtain a speed of 1 mm/ms, while, for a speed of 7 mm/ms, they-axis was set at 28.5. The termination time used also adjusted to the speed applied.In the numerical modelling of the impact attenuator, contacts were used, mainly automatic node to surface and automatic single surface. Automatic node to surface contact was used to contact the rigid plate and among cans with a rigid plate and anti-intrusion plate as support.

    Based on the preceding collision analysis, this kind of contact is appropriate for contacts concerning inflexible objects. The type of automatic single surface contacts was used on the wall of the can itself. This contact kind is suitable for crashworthiness applications.Based on the body interaction, the contact among the plate and the attenuator was defined as a friction phenomenon with the values of static and dynamic friction coefficients at 0.4 and 0.3, consecutively. In this research, we used two different methods to compare the final results. The first was a quasi-static test with a velocity of 1 mm/ms, and the second one had velocity according to the rules,namely 7 mm/ms. We also varied the thickness of the used cans,with the first thickness being 0.5 mm (designated as A1 and A3)and the second being 1 mm (designated as A2 and A4). Table 3 presents the variations to simulate the various conditions.

    Table 3 Impact attenuator variations.

    The first result is the crushing force curve as shown in Fig. 6 .The crushing force curve was derived from the total amount of force received by the can. Based on the presented results in this curve, the crushing forces A1 and A2 had velocities of 1 mm/ms.The force that works at the beginning experiences a fluctuating force value so that the curve goes up and down. With deformation around 175 mm, the curve tends to go up with a maximum force value of 1792.75 kN for A1 and 693.228 kN for A2.

    For the crushing forces A3 and A4, the curve has a similar trend where, initially, the force received was fluctuating, and, when the displacement reached about 100 mm, the force value tended to increase. The force value tended to decrease but still appeared to fluctuate. In both graphs, the curve appears to have increased in the displacement area of about 100 mm. This is due to the holding force on a plate that separates the two levels of the cans so that the force’s value increases so that the bottom can begin to deform.

    The internal energy ( Fig. 7 ) that occurred in A1 and A2 appears to have a similar curve trend. The same is true for the curves for A3 and A4. The A2 and A4 curves, which were cans with a thickness variation of 1 mm, had higher values of 16782.6 J for A2 and 16380.3 J for A4. From this, it can be seen that, with the same thickness variation, the results were almost the same even though the velocity of deviation was different. A1 and A3 were variations with a thickness of 0.5 mm. The results above show that the two curves were under the curve with a thickness variation of 1 mm.For A1, the internal energy reached a value of 5974.88 J, and for A3, it obtained a value of 7529.11 J.

    Fig. 10. Progressive contours of the von-Mises stress: (a) t = 50 ms; (b) t = 100 ms; (c) t = 150 ms; and (d) t = 200 ms.

    Fig. 11. Progressive contours of the z-stress: (a) t = 50 ms; (b) t = 100 ms; (c) t = 150 ms; and (d) t = 200 ms.

    For kinetic energy ( Fig. 8 ), the curves from the compression test with a velocity of 1 mm/ms are for A1 and A2, and the impact test curves with a velocity of 7 mm/ms are for A3 and A4. It can be seen from the figure that the value of the kinetic energy with a velocity of 7 mm/ms had a stronger value with a large difference when compared to the kinetic energy curve with a speed of 1 mm/ms. This happens because kinetic energy is affected by the rate and, the higher the collision velocity, the greater the kinetic energy. It can be seen from the two graphs that they have almost similar curve trends.

    The total energy ( Fig. 9 ) is the total amount of energy generated at the impact attenuator collision. From this curve, it can be seen that the highest total energy value was the total energy A4—namely the impact attenuator with a thickness of 1 mm and a velocity of 7 mm/ms with a value reaching 26833.4 J.

    In addition to displaying the values of the energy and crushing force curves, the results also show the deformation display of the impact attenuator. The values that will be displayed in this research are the von-Misses stress, z-stress, z-displacement, and plastic strain.

    Stress behavior can only be seen in certain circumstances. The red fringe levels in Fig. 10 show the maximum stress values in GPa.Therefore, in both display contours, the red color at the edge level represents the most vulnerable part of the simulated outline, while the blue color represents the most durable part. In contrast to the von-Mises stress, the stress on the z-axis has a negative sign on stress value, which indicates that the direction of stress is opposite to the direction of the axis. The red contour represents the maximum stress value along the axis and vice versa in blue, as shown in Fig. 11 .

    Fig. 12. Progressive contours of the z-displacement: (a) t = 50 ms; (b) t = 100 ms; (c) t = 150 ms; and (d) t = 200 ms.

    Fig. 13. Progressive contours of the plastic strain: (a) t = 50 ms; (b) t = 100 ms; (c) t = 150 ms; and (d) t = 200 ms.

    The results of the structural displacement in Fig. 12 show the displacement result on the z-axis. A positive sign on the displacement value indicates that the element is moving in the direction of the axis and vice versa for the negative sign. The maximum values are shown in red for the fringe levels and blue for the minimum values.

    Based on the contour of the effective plastic strain, we can observe that the impact attenuator did not experience heavy strain because almost all impact attenuators are blue, indicating that the element experienced a minimum strain. As shown in Fig. 13 , the maximum strain in the fringe scale is 0.8234 at the 50 ms step,and the maximum strain is 1.103 at 100 ms with the can at the top already deforming all over the surface. At the time of 150 ms,the maximum strain increased compared to the time of 100 ms,which was 1.677. At about 150 ms, the arrangement of the lower cans appeared as deformed.

    From the internal energy curve, the highest energy value is found on the A2 curve with a value of 16782.6 J. Based on the terms and conditions that apply to the 2019 FSAE rules, the energy must be absorbed as at least 7350 J. Based on the above statement, we concluded that the A2 impact attenuator met the requirements of the 2019 FSAE regulations. The A4 impact attenuator also completed the requirements, while A1 and A3 did not meet the requirements because the value did not reach 7350 J. From this study, the collector’s velocity did not have a significant influence on the internal energy because the same thickness variation resulted in an internal energy value with a low difference and a similar curve trend.

    From the results of the deformation, when viewed in the zstress section, the deformation of the impact attenuator at each step tended to increase. For the displacement that occurred on thez-axis, because the load moves against thez-axis, the displacement results displayed a high value. As for the strain, part of the impact attenuator experienced heavy strain at the end of the collision.This statement can be concluded based on contours where almost all the elements are blue, which indicates that the deformable element was flattened by the intrusion plate (a rigid component in the structure).

    The impact attenuator in this study produced two variations that have values that exceeded the requirements of the FSAE rules where the passed variation had a can thickness of 1 mm. For impact attenuators with used cans, when compared to other materials that tend to be expensive and difficult to make, this canned food waste can be an alternative material to reduce costs. The availability of goods is relatively easy to obtain, so it can be said that this material is cheap because it can take advantage of used cans in relatively good condition. It is also possible that this research can contribute to protecting the environment by recycling the waste can into fundamental material in manufacturing the impact attenuator.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    听说在线观看完整版免费高清| 国产精品99久久99久久久不卡 | 国产女主播在线喷水免费视频网站| 在线观看av片永久免费下载| 亚洲av福利一区| 王馨瑶露胸无遮挡在线观看| 特大巨黑吊av在线直播| 国产毛片a区久久久久| 2021天堂中文幕一二区在线观| 一级a做视频免费观看| 美女cb高潮喷水在线观看| 中文字幕免费在线视频6| 亚洲av中文字字幕乱码综合| 亚洲精品视频女| 中国三级夫妇交换| 黄色怎么调成土黄色| 乱系列少妇在线播放| 中文字幕人妻熟人妻熟丝袜美| 人妻制服诱惑在线中文字幕| 日日啪夜夜撸| 中国国产av一级| 成年人午夜在线观看视频| 精品人妻一区二区三区麻豆| 欧美人与善性xxx| 国产黄色免费在线视频| 国产亚洲最大av| 最后的刺客免费高清国语| 亚洲精品中文字幕在线视频 | 插阴视频在线观看视频| 卡戴珊不雅视频在线播放| 一区二区三区免费毛片| 日韩免费高清中文字幕av| 久久精品夜色国产| 国产亚洲最大av| 男女边吃奶边做爰视频| 午夜激情久久久久久久| 久久99热这里只有精品18| 国国产精品蜜臀av免费| 久久韩国三级中文字幕| 天堂俺去俺来也www色官网| 欧美精品一区二区大全| 亚洲天堂av无毛| 国产精品久久久久久久久免| 人妻少妇偷人精品九色| 在线观看美女被高潮喷水网站| xxx大片免费视频| 欧美一区二区亚洲| 午夜爱爱视频在线播放| 狠狠精品人妻久久久久久综合| 国产成人午夜福利电影在线观看| 久久热精品热| 亚洲自拍偷在线| 久久精品久久久久久噜噜老黄| av播播在线观看一区| 精品国产三级普通话版| 性插视频无遮挡在线免费观看| 免费看av在线观看网站| 欧美+日韩+精品| 天天躁夜夜躁狠狠久久av| 亚洲一级一片aⅴ在线观看| 一二三四中文在线观看免费高清| 少妇 在线观看| 男女那种视频在线观看| 国产有黄有色有爽视频| 99久久中文字幕三级久久日本| 久久精品久久精品一区二区三区| 特大巨黑吊av在线直播| 人妻制服诱惑在线中文字幕| 亚洲自拍偷在线| 18禁在线播放成人免费| 国产精品久久久久久精品古装| 久久久欧美国产精品| 久久久久久九九精品二区国产| 亚洲国产欧美在线一区| 国产老妇伦熟女老妇高清| 免费黄频网站在线观看国产| 日本黄大片高清| 久久久久久伊人网av| 午夜福利视频精品| 日日啪夜夜爽| 精品熟女少妇av免费看| 2018国产大陆天天弄谢| 乱系列少妇在线播放| 青春草视频在线免费观看| 亚洲成人av在线免费| 少妇人妻 视频| 日本av手机在线免费观看| 青春草国产在线视频| 日韩欧美一区视频在线观看 | 国产永久视频网站| 亚洲精品国产色婷婷电影| 乱系列少妇在线播放| 免费人成在线观看视频色| 国产精品一区二区性色av| 爱豆传媒免费全集在线观看| av女优亚洲男人天堂| 在线观看一区二区三区激情| av网站免费在线观看视频| 国产亚洲91精品色在线| 精品人妻熟女av久视频| 日日摸夜夜添夜夜添av毛片| 五月伊人婷婷丁香| 18禁裸乳无遮挡免费网站照片| 伊人久久国产一区二区| 寂寞人妻少妇视频99o| 免费不卡的大黄色大毛片视频在线观看| 97热精品久久久久久| av线在线观看网站| 夜夜看夜夜爽夜夜摸| 日韩三级伦理在线观看| 亚洲最大成人手机在线| 久久久精品免费免费高清| 午夜福利高清视频| 久久ye,这里只有精品| 亚洲精品成人av观看孕妇| 国产精品国产av在线观看| 男女啪啪激烈高潮av片| 在线观看美女被高潮喷水网站| 蜜桃亚洲精品一区二区三区| 免费av不卡在线播放| 一级毛片 在线播放| av卡一久久| 亚洲av在线观看美女高潮| .国产精品久久| 免费av毛片视频| 国语对白做爰xxxⅹ性视频网站| 又爽又黄a免费视频| 国产av码专区亚洲av| 最近最新中文字幕免费大全7| 插阴视频在线观看视频| 久久精品熟女亚洲av麻豆精品| 欧美xxⅹ黑人| 激情 狠狠 欧美| 日韩一区二区三区影片| 伦理电影大哥的女人| 草草在线视频免费看| 欧美日韩国产mv在线观看视频 | 免费av不卡在线播放| 日韩大片免费观看网站| 最近的中文字幕免费完整| 免费av观看视频| 2021天堂中文幕一二区在线观| 日本av手机在线免费观看| 免费黄网站久久成人精品| 波野结衣二区三区在线| 自拍欧美九色日韩亚洲蝌蚪91 | 国产亚洲av嫩草精品影院| 天堂网av新在线| 久久久色成人| 国产又色又爽无遮挡免| 国产一区二区三区av在线| av国产精品久久久久影院| 亚洲av不卡在线观看| 久久久久久久久久久免费av| 又粗又硬又长又爽又黄的视频| 国产精品一区www在线观看| 国产又色又爽无遮挡免| www.色视频.com| 禁无遮挡网站| 国产伦理片在线播放av一区| 亚洲伊人久久精品综合| 国产美女午夜福利| 精品酒店卫生间| 亚洲av男天堂| 一级黄片播放器| 特级一级黄色大片| 最近最新中文字幕免费大全7| 国产淫片久久久久久久久| 亚洲人成网站在线观看播放| 又大又黄又爽视频免费| av女优亚洲男人天堂| 中文资源天堂在线| 日日撸夜夜添| 国产人妻一区二区三区在| 午夜免费男女啪啪视频观看| 搡女人真爽免费视频火全软件| 能在线免费看毛片的网站| 成人一区二区视频在线观看| 精品久久久久久久久亚洲| 欧美xxxx黑人xx丫x性爽| 亚洲av成人精品一区久久| 日韩人妻高清精品专区| 91久久精品国产一区二区成人| 新久久久久国产一级毛片| 1000部很黄的大片| 日本av手机在线免费观看| 国产日韩欧美亚洲二区| 亚洲精品日韩在线中文字幕| 王馨瑶露胸无遮挡在线观看| 日韩在线高清观看一区二区三区| 七月丁香在线播放| 韩国av在线不卡| 黄色一级大片看看| 亚洲内射少妇av| 国产一级毛片在线| 国产高潮美女av| 黄色怎么调成土黄色| 国产日韩欧美亚洲二区| 欧美日韩精品成人综合77777| 亚洲精品国产av成人精品| 国产又色又爽无遮挡免| 国产一区亚洲一区在线观看| 特大巨黑吊av在线直播| 亚洲精品日韩在线中文字幕| 欧美xxxx黑人xx丫x性爽| 看黄色毛片网站| 免费播放大片免费观看视频在线观看| 久久久久网色| 国产成人freesex在线| 免费av毛片视频| 少妇熟女欧美另类| 大陆偷拍与自拍| 欧美另类一区| 高清毛片免费看| 男女啪啪激烈高潮av片| 亚洲婷婷狠狠爱综合网| 欧美高清成人免费视频www| 亚洲国产精品专区欧美| 女人被狂操c到高潮| 肉色欧美久久久久久久蜜桃 | av福利片在线观看| 深爱激情五月婷婷| 人妻 亚洲 视频| 日韩中字成人| 国产亚洲av嫩草精品影院| 99re6热这里在线精品视频| 大码成人一级视频| 亚洲不卡免费看| 欧美精品人与动牲交sv欧美| 日日摸夜夜添夜夜爱| 亚洲在线观看片| 哪个播放器可以免费观看大片| 99久国产av精品国产电影| 国产伦在线观看视频一区| 我的女老师完整版在线观看| 热re99久久精品国产66热6| 国产精品久久久久久久电影| 亚洲aⅴ乱码一区二区在线播放| 老女人水多毛片| 亚洲精品,欧美精品| 国内精品宾馆在线| 三级经典国产精品| 免费大片18禁| 国产伦精品一区二区三区四那| 男女下面进入的视频免费午夜| 国产老妇伦熟女老妇高清| www.av在线官网国产| 97在线视频观看| 99久久人妻综合| 自拍偷自拍亚洲精品老妇| 国产精品熟女久久久久浪| 国产成人freesex在线| 国产一区二区亚洲精品在线观看| 久久久久久久久久久丰满| 男人舔奶头视频| 亚洲国产最新在线播放| 国产成人freesex在线| 国产亚洲午夜精品一区二区久久 | 久久久精品94久久精品| av黄色大香蕉| 亚洲天堂av无毛| 亚洲欧美一区二区三区国产| 我的老师免费观看完整版| 一二三四中文在线观看免费高清| 精品少妇黑人巨大在线播放| 99九九线精品视频在线观看视频| 成年女人在线观看亚洲视频 | 亚洲丝袜综合中文字幕| 亚洲国产av新网站| 18禁动态无遮挡网站| 亚洲av国产av综合av卡| 最近2019中文字幕mv第一页| 亚洲第一区二区三区不卡| av国产精品久久久久影院| 午夜免费男女啪啪视频观看| 亚洲欧洲日产国产| 天堂网av新在线| 久久久成人免费电影| 国产亚洲一区二区精品| 国产精品国产三级专区第一集| 男人和女人高潮做爰伦理| 国产毛片a区久久久久| av一本久久久久| 亚洲av成人精品一区久久| 免费播放大片免费观看视频在线观看| 晚上一个人看的免费电影| 日日啪夜夜撸| 日韩 亚洲 欧美在线| av.在线天堂| 男男h啪啪无遮挡| 亚洲,欧美,日韩| 一边亲一边摸免费视频| 青春草视频在线免费观看| 亚洲国产精品专区欧美| 久久久久久久久久久丰满| 日本午夜av视频| 久久精品夜色国产| 男女边摸边吃奶| 男插女下体视频免费在线播放| 国产精品.久久久| 欧美日韩视频高清一区二区三区二| 日韩电影二区| xxx大片免费视频| 日日啪夜夜爽| 国产成人91sexporn| 在线亚洲精品国产二区图片欧美 | 国产美女午夜福利| 精品酒店卫生间| 精品少妇黑人巨大在线播放| 午夜福利网站1000一区二区三区| 激情五月婷婷亚洲| 18禁动态无遮挡网站| 白带黄色成豆腐渣| 成人国产av品久久久| 免费人成在线观看视频色| 亚洲第一区二区三区不卡| 97超视频在线观看视频| 最近2019中文字幕mv第一页| 九草在线视频观看| 18禁裸乳无遮挡动漫免费视频 | 国产 一区精品| 美女cb高潮喷水在线观看| 亚洲av中文字字幕乱码综合| 久久99蜜桃精品久久| 久久影院123| 欧美xxxx性猛交bbbb| 一个人看视频在线观看www免费| 色吧在线观看| 99久久人妻综合| 久久精品国产自在天天线| 亚洲成人一二三区av| av在线观看视频网站免费| 久久久久国产精品人妻一区二区| av网站免费在线观看视频| 久久久色成人| 中文天堂在线官网| 一级毛片我不卡| av在线app专区| 天天躁夜夜躁狠狠久久av| 99久久人妻综合| 伦理电影大哥的女人| 亚洲自偷自拍三级| 中文乱码字字幕精品一区二区三区| 内射极品少妇av片p| 亚洲精品色激情综合| 精品一区二区三区视频在线| 成人午夜精彩视频在线观看| 亚洲欧美精品专区久久| 91久久精品国产一区二区成人| 免费黄频网站在线观看国产| 日韩av在线免费看完整版不卡| 国产高清有码在线观看视频| 日日啪夜夜爽| 成人毛片60女人毛片免费| 日韩强制内射视频| 女的被弄到高潮叫床怎么办| 国产有黄有色有爽视频| 久久ye,这里只有精品| 亚洲图色成人| 22中文网久久字幕| 欧美性感艳星| 久久精品国产亚洲av天美| 91久久精品电影网| 久久97久久精品| 国产一区有黄有色的免费视频| 午夜福利在线在线| av在线天堂中文字幕| 午夜精品国产一区二区电影 | 精品国产乱码久久久久久小说| 成人无遮挡网站| 狠狠精品人妻久久久久久综合| 亚洲国产高清在线一区二区三| 久久久久久国产a免费观看| 建设人人有责人人尽责人人享有的 | 亚洲精品日韩在线中文字幕| 天天躁日日操中文字幕| tube8黄色片| 国产精品偷伦视频观看了| 性色avwww在线观看| 亚洲无线观看免费| 亚洲一级一片aⅴ在线观看| a级毛色黄片| 在线a可以看的网站| 丰满乱子伦码专区| 亚洲不卡免费看| 边亲边吃奶的免费视频| 一二三四中文在线观看免费高清| 国产在线男女| 久久久久久伊人网av| 人妻系列 视频| 人妻制服诱惑在线中文字幕| 亚洲av男天堂| 日本黄色片子视频| 国产欧美亚洲国产| 国产成人免费观看mmmm| 一本久久精品| 18+在线观看网站| 亚洲成人中文字幕在线播放| 久久久色成人| 久久久久国产网址| 菩萨蛮人人尽说江南好唐韦庄| 丝袜脚勾引网站| 看十八女毛片水多多多| 麻豆国产97在线/欧美| 精品亚洲乱码少妇综合久久| 亚洲真实伦在线观看| videossex国产| 欧美老熟妇乱子伦牲交| av免费在线看不卡| 国产亚洲一区二区精品| 日本-黄色视频高清免费观看| 久久99热这里只频精品6学生| av专区在线播放| 国内揄拍国产精品人妻在线| 免费黄色在线免费观看| 你懂的网址亚洲精品在线观看| 日日摸夜夜添夜夜爱| 国产精品不卡视频一区二区| 久久99热这里只频精品6学生| 国产精品.久久久| 最后的刺客免费高清国语| 高清午夜精品一区二区三区| 男人爽女人下面视频在线观看| 亚洲精品乱码久久久久久按摩| 久久人人爽人人片av| 少妇人妻一区二区三区视频| 欧美成人午夜免费资源| 精品午夜福利在线看| 人妻 亚洲 视频| 99九九线精品视频在线观看视频| 女人久久www免费人成看片| 能在线免费看毛片的网站| 制服丝袜香蕉在线| 韩国av在线不卡| av免费观看日本| 国产在线一区二区三区精| 亚洲av免费在线观看| 青春草视频在线免费观看| 亚洲精品国产色婷婷电影| 久久精品国产亚洲av涩爱| 女人十人毛片免费观看3o分钟| 久久久久久久久久人人人人人人| 亚洲国产色片| 街头女战士在线观看网站| 亚洲人成网站在线播| 日韩大片免费观看网站| 少妇的逼水好多| 欧美一区二区亚洲| .国产精品久久| 99久国产av精品国产电影| 日韩av在线免费看完整版不卡| 一二三四中文在线观看免费高清| 久久精品久久久久久噜噜老黄| 建设人人有责人人尽责人人享有的 | 日日摸夜夜添夜夜爱| 欧美高清成人免费视频www| 精品久久久久久久久av| 涩涩av久久男人的天堂| av国产免费在线观看| 久久人人爽人人片av| 男女边摸边吃奶| 人妻少妇偷人精品九色| 免费观看av网站的网址| 国语对白做爰xxxⅹ性视频网站| 三级男女做爰猛烈吃奶摸视频| 亚洲精品久久午夜乱码| 免费看光身美女| 欧美日韩国产mv在线观看视频 | 亚洲欧洲日产国产| 欧美+日韩+精品| 好男人视频免费观看在线| 两个人的视频大全免费| 性色av一级| 精品国产三级普通话版| 国产淫语在线视频| 久久人人爽av亚洲精品天堂 | 日日撸夜夜添| 看黄色毛片网站| 欧美激情在线99| 精品熟女少妇av免费看| 久久99热6这里只有精品| 少妇人妻 视频| 联通29元200g的流量卡| 日本av手机在线免费观看| 人体艺术视频欧美日本| 亚洲欧美一区二区三区国产| 精品国产乱码久久久久久小说| 又黄又爽又刺激的免费视频.| 国产亚洲最大av| 99久久九九国产精品国产免费| 国内精品宾馆在线| 亚洲自拍偷在线| 中文字幕人妻熟人妻熟丝袜美| 欧美少妇被猛烈插入视频| 搞女人的毛片| 一本久久精品| 一区二区三区精品91| 久久精品国产鲁丝片午夜精品| 色哟哟·www| 内射极品少妇av片p| 日韩免费高清中文字幕av| 亚洲精品乱码久久久久久按摩| 欧美极品一区二区三区四区| 欧美精品国产亚洲| 精品少妇黑人巨大在线播放| 免费观看的影片在线观看| 久久99精品国语久久久| 99久久人妻综合| 久久精品综合一区二区三区| 久久久久久久久久久免费av| 国产精品99久久99久久久不卡 | 精品视频人人做人人爽| av国产精品久久久久影院| 日本-黄色视频高清免费观看| 99热这里只有是精品在线观看| 日韩制服骚丝袜av| 午夜免费男女啪啪视频观看| 好男人视频免费观看在线| 日韩av不卡免费在线播放| 亚洲人成网站在线播| 亚洲内射少妇av| 欧美激情久久久久久爽电影| 丝袜脚勾引网站| 日韩伦理黄色片| 成年女人在线观看亚洲视频 | 国内精品美女久久久久久| 婷婷色av中文字幕| xxx大片免费视频| 内射极品少妇av片p| 亚洲最大成人手机在线| 成人黄色视频免费在线看| 在线观看人妻少妇| 一级毛片久久久久久久久女| 日本-黄色视频高清免费观看| 超碰97精品在线观看| 各种免费的搞黄视频| 99re6热这里在线精品视频| 大香蕉97超碰在线| 免费观看a级毛片全部| 成人美女网站在线观看视频| av国产精品久久久久影院| 国产女主播在线喷水免费视频网站| 蜜桃亚洲精品一区二区三区| 2018国产大陆天天弄谢| 白带黄色成豆腐渣| 狂野欧美白嫩少妇大欣赏| 久久99热这里只有精品18| 亚洲av中文字字幕乱码综合| videossex国产| 欧美亚洲 丝袜 人妻 在线| 日日啪夜夜爽| 欧美一区二区亚洲| 亚洲欧美精品自产自拍| 中文字幕制服av| 久久久久久国产a免费观看| 亚洲av二区三区四区| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品专区欧美| av在线亚洲专区| 热99国产精品久久久久久7| 亚洲欧美日韩无卡精品| 免费少妇av软件| 国精品久久久久久国模美| 国产女主播在线喷水免费视频网站| 亚洲人成网站在线播| 亚洲精品国产成人久久av| 人体艺术视频欧美日本| 国产精品久久久久久精品古装| 国产成人精品婷婷| 久久久久久久久久久免费av| 免费在线观看成人毛片| av免费在线看不卡| 国产探花极品一区二区| 午夜爱爱视频在线播放| 国产成人91sexporn| 麻豆成人av视频| 国产片特级美女逼逼视频| 搡女人真爽免费视频火全软件| 亚洲精品一二三| 2021少妇久久久久久久久久久| 亚洲自拍偷在线| 成人一区二区视频在线观看| 国产免费一区二区三区四区乱码| 国产黄频视频在线观看| 99热这里只有是精品50| 丝袜美腿在线中文| 国产精品国产三级专区第一集| 一级爰片在线观看| 国产成人免费无遮挡视频| 汤姆久久久久久久影院中文字幕| 一级爰片在线观看| 在线观看免费高清a一片| 日韩av在线免费看完整版不卡| 国产 一区 欧美 日韩| 国产综合精华液| 我的老师免费观看完整版| 别揉我奶头 嗯啊视频| 久久人人爽av亚洲精品天堂 | 97在线视频观看| 在线天堂最新版资源| 18禁动态无遮挡网站| 亚洲精品国产色婷婷电影| 午夜福利在线观看免费完整高清在| 香蕉精品网在线| 亚洲精品国产色婷婷电影| 日本-黄色视频高清免费观看| 99热这里只有精品一区| 99九九线精品视频在线观看视频| 两个人的视频大全免费| 尤物成人国产欧美一区二区三区| 深爱激情五月婷婷| 男的添女的下面高潮视频| 国产乱人偷精品视频| 国产黄片视频在线免费观看| 亚洲精品视频女|