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

    Experimental and numerical study of the blast wave decrease using sandwich panel by granular materials core

    2021-11-03 13:24:04MsoudRhmniAlirezNddfOskoueiAminMoslemiPetrudi
    Defence Technology 2021年5期

    Msoud Rhmni ,Alirez Nddf Oskouei ,Amin Moslemi Petrudi

    a Department of Mechanical Engineering,IHU University,Tehran,Iran

    b Department of Mechanical Engineering,Tehran University,Tehran,Iran

    Keywords: Granular materials Shock tube Blast wave Numerical simulations

    ABSTRACT Among the intrinsic properties of some materials,e.g.,foams,porous materials,and granular materials,are their ability to mitigate shock waves.This paper investigated shock wave mitigation by a sandwich panel with a granular core.Numerical simulations and experimental tests were performed using Autodyn hydro-code software and a shock tube,respectively.The smoothed particle hydrodynamics (SPH)method was used to model granular materials.Sawdust and pumice,whose properties were determined by several compression tests,were used as granular materials in the sandwich panel core.These granular materials possess many mechanisms,including compacting (e.g.,sawdust) and crushing (e.g.,pumice)that mitigate shock/blast wave.The results indicated the ineffectiveness of using a core with low thickness,yet it was demonstrated to be effective with high thickness.Low-thickness pumice yielded better results for wave mitigation.The use of these materials with a core with appropriate core reduces up to 88%of the shock wave.The results of the experiments and numerical simulations were compared,suggesting a good agreement between the two.This indicates the accuracy of simulation and the ability of the SPH method to modeling granular material under shock loading.The effects of grain size and the coefficient of friction between grains have also been investigated using simulation,implying that increasing the grain size and coefficient of friction between grains both reduce overpressure.

    1.Introduction

    Shock waves in two or multi-phase environments include several phenomena,which have been considered by mathematicians,physicists,and mechanical and aerospace engineers.In some engineering designs,there is a need to reduce the power of waves.Wave dampers are required to build structures and armor to mitigate the blast wave in the military industry and the damping of the sound wave in the construction industry.The use of explosive wave dampers is not limited to these cases.Moreover,it is applicable in the explosive shaping industry,the construction of tunnels,and explosion test tanks.Foams,porous granular materials,and layered materials have been commonly used to reduce the impact and damping of the blast wave.These capabilities are due to their high energy absorption,among the energy-absorbing materials,relatively inexpensive grain materials,such as pumice,iron filings,and wood chips can be considered.Several studies have been carried out on the use of granular materials.Nesterenko[1]examined the shock wave damping of soft condensed materials,in which the effect of using sawdust in a metal tank was investigated in the explosion test and observed that the deformation of the tank was greatly reduced and the subsequent degradation was reduced.Serna et al.[2]investigated the reduction of the blast wave in grain flow.Nguyen [3] studied the shock wave dynamics in a chain of granular materials,in which this chain was a simplified specimen of the granular material used to study the dynamics of wave propagation.Gueders et al.[4]simulated the explosion-wave damping in granular materials of alumina ceramic by examining the pressure and impulse of the wave.Moreover,Levy[5]investigated the shock wave contact with granular materials.It is necessary to create the criteria for the properties and geometry of materials,which will ensure the mitigation for specific conditions of blast loading.The fact that such materials are appropriate energy absorbers does not guarantee their mitigation performance.For instance,the application of a porous layer with a small thickness equal to an effective radius of explosive charge did not result in the damping of strains in the explosive chamber [1].In the past two decades,granular materials have received significant attention due to their military and civil applications,such as impulse reduction,freeze plumbing,silos,and air purification [5].Liu et al.[6] presented a numerical simulation of the shock-powder interaction based on the particle kinetic theory and gave a close detection to the particle phase response behavior to the passage of shock wave within the dense particle powder.The simulations were carried out using a developed code compiled in FORTRAN.They used particle kinematic theory to model the powder and compared the results with the experimental results.According to the results,the propagation of waves in the particle powder showed a decreasing trend in the gas phase pressure.Burgoyne et al.[7] employed numerical simulation to investigate the two-and three-dimensional shock absorption among grain crystals.Their experimental and computational study showed that the behavior of three-dimensional grain crystals could be simulated.The rules of normal contact between components are well described,and it also turns out that the properties of a onedimensional granular chain can also be generalized and applied to two-and three-dimensional conditions.Luo et al.[8] investigated the shock wave contact with granular materials.In this study,a miniature shock tube was used.Furthermore,they used a highspeed camera to record wave propagation in the granular material and used metallic spheres as the geometry of granular material.Duan et al.[9] studied the shock wave in a one-dimensional granular chain by contacting Hertz between the particles.They found that the shock wave velocity increases with increasing piston velocity and initial overlap.Kandan et al.[10]studied the instability of the granular position surface for shock loading.In their study,the growth and instability of granular materials loaded by air shock waves were studied using the shock tube experiments and numerical simulation.The Eulerian-Lagrangian solver and the Druker-Prager model were used for granular materials.In a study,Fan et al.[11]examined the computational study on Peridynamics modeling and simulation of soil fragmentation under explosive loads.In their study,the SPH method was used.The results of the simulation were compared with experimental data,which are relatively close to each other.Moreover,the study showed that employing the Druker-Prager model,which describes steady-state and quasistatic loads behavior of granular soil material,may not be suitable for describing soil behavior at high strain rates.Mirova et al.[12]investigated the reduction of shock wave passing through a sand layer by changing the position of the sand layer installed at the end of the shack tube.Liang et al.[13]studied the shock waves in a onedimensional bead produced by a constant velocity impact in a short period using the numerical simulation.In their study,a cylinder and piston were used to generate the wave.It was found that the attenuation strength depends on the plasticity,the piston velocity,and the bead radius (see Tables 3 and 4).

    This paper is an experimental and numerical study on the blast wave decrease using a sandwich panel by granular materials core.In this study,mineral pumice and sawdust were used as granular materials.Several compression tests were performed on them to determine the behavior of these materials under pressure and to extract the required data of numerical simulation,to achieve the properties of these materials.The simulation results were validated using several experimental results of the shock tube.

    2.Numerical simulation

    The solid materials in regular loading often have elastic behavior,but in the shock loading,they are usually entered into the plastic phase.The equation of state is responsible for explaining the hydrodynamic behavior of the material.Hydrodynamic behavior in most materials is best described by high-rate dynamic loading of the material.For example,the response of gases and liquids that do not tolerate shear is hydrodynamic and the pressure in them is a function of density and internal energy.Also,the behavior of solids is influenced by high hydrodynamic strain rates,where the applied pressure to the solids is much higher than that of the yield stress,and the strength of the material against the applied stress can be neglected.A material model in general requires equations to relate stress to deformation and internal energy to temperature.The equation of state can be determined from knowledge of the thermodynamic properties of the material and ideally,should not require dynamic data to build up the relationship.However,in practice,the only practical way of obtaining data on the behavior of the material at high strain rates is to carry out well-characterized dynamic experiments [14].To define a material in Autodyn Hydro code software,it is necessary to define the Equation of state(EOS),the Strength model,and the fracture model,each of which is defined by the type of material and its behavior.Once the matter is properly defined,the problem solver must be solved,the appropriate solver for the problem is defined according to the type of matter and the physics of the problem.Eulerian solver is generally used for gases and the Lagrangian Solver is used for solids.But the granular material is out of these two states of matter.Some researchers have used discrete element method (DEM) codes for granular materials,a problem of the DEM method that has made it less developed in particle mechanics is that it requires a large number of coefficients to be able to use the coefficients of hardness the springs and estimate the viscosity between the particles and the friction.Some of these coefficients do not even have a physical meaning,and some of them can only be obtained from highprecision bi-axial tests.Also,the numerical computation of the DEM method requires a computer with high computational capacity and time-consuming.Arbitrary Lagrangian-Eulerian Equation(ALE)solvers also can be used in Autodyn software for granular materials,which have been used successfully to model these materials in some studies [15,16].In this study,using Smoothed Particle Hydrodynamics(SPH)has been investigated.This method is a mesh-independent method such as the DEM method,but it is also broader in scope and its algorithm and computer program are further developed.So far,very few articles have used this method,but the results of using this method can be promising.The equations of state,the Strength model,and the fracture model used for the materials in the simulation are presented in Table 1(see Table 2).

    Table 1 The equations of state,the strength model,and the fracture model used for materials.

    Table 2 Changes in density and speed of sound at different pressures for sawdust.

    Table 3 Changes in density and speed of sound at different pressures for the mineral pumice.

    Table 4 Sawdust and mineral pumice yield stress at different densities.

    Table 5 Shear Modulus of sawdust and mineral pumice at different densities.

    To perform numerical simulation,it is necessary to introduce the properties of matter to the software,Properties required for simulation of granular material in Autodyn software include changes in density to pressure,changes in yield stress to pressure,changes in shear modulus to pressure,and changes in the speed of sound in granular material in different densities.Using axial and non-axial compression tests other than the speed of sound at different densities,the remaining required parameters can be obtained[17].In Fig.1 the compression test on the specimen is shown.

    In this study,a press machine with a capacity of 30 tons and a compression rate of 2 mm/min as shown in Fig.1 has been used.The internal diameter of the matrix used in this experiment is 7 mm and the cavity depth was 100 mm.The matrix and punch used in the test are shown in Fig.1.The pumice grains are poured into the matrix cavity up to a height of 75 mm.The pressure to density changes can be obtained using an axial press.To obtain the most accurate results,three specimens of each material were tested and the data obtained from the three simulations were averaged using the data obtained from the three experiments.Two of the diagrams obtained from the compression test for sawdust and mineral pumice are shown in Figs.2 and 3.

    Using the mean graph data,the density-pressure changes can be introduced into software,which is needed to define the compression equations of state.Considering some simplification assumptions,other properties needed for numerical simulation can be obtained from the compression test.In compression test specimens only the volume changes with the length of the cylinder,to define the granular material model a nonlinear elastic-plastic model can be assumed,Properties are assumed isotropic (homogeneous).In the elastic portion,the bulk modulus k is a function of the material density and the Poisson ratio can also be constant and is approximately 0.29.Then the relation between the elastic constants is as Eq.(1).

    Fig.1.Compression test on granular material specimens.

    In this relation,G is the shear modulus.The speed of sound in the elastic state in a grain material is obtained by Eq.(2) [18].

    The bulk modulus is obtained by Eq.(3).

    Where P is the pressure and εvis the volume strain[18].Therefore,using Eqs.(1)-(3),the speed of sound and the shear modulus at different densities are calculated.To calculate the yield stress of granular materials,a simplified relation is used by Coulomb [19],which for non-adhesive granular materials,Eq.(4) is presented:

    Fig.2.The results of the sawdust compression test.

    Fig.3.The results of the mineral pumice compression test.

    Fig.4.The Mohr-Coulomb failure criterion,for the principal stresses and the angle of internal friction [19].

    Fig.5.Nonlinear loading pattern in compression equation of state.

    Fig.6.Simulation of the problem in the Autodyn software.

    The Coulomb yield criterion for the granular material can be written as Eq.(5) from the principal stresses[20]:

    In these relations p is the axial pressure,Cand φ represent the cohesion and internal friction angle of the granular material.The friction angle is about 43°for the sawdust and about 37°to 45°[21,22] for the mineral pumice.C is the coefficient of adhesion,which for non-adhesive granular materials such as pumice is a very small(C ≈0 kPa)[23,24]and for sawdust with a moisture content of less than 10%,it is about 7 kPa[24].The coefficient of friction for mineral pumice is between 0.3 and 0.7 and for sawdust is between 0.15 and 0.25.In this study,the average values,0.5 for mineral pumice,and 0.2 for sawdust according to Refs.[22,24] are considered.In both cases,the humidity is less than 10%.The principal stresses and internal friction angles are by Mohr-Coulomb yield as shown in Fig.4.

    Eq.(6) is provided for soil shear strength and can be used for other granular materials.

    In this relation σ is the normal stress[19].Using Eqs.(1)-(6)and laboratory results,the yield stress and shear modulus at different pressures and densities for the studied granular materials(soil and mineral pumice) are calculated.The calculation results are presented in Tables (2) to(5).

    For the granular materials,the compression equation of state is used in this study.Using this model,there is more control over the porous materials during the loading and unloading phase.In this equation of state,the speed of sound is considered as a function of density.This model is widely used in shock simulation in granular materials such as modeling of soil,sand,etc in Autodyn software.The compression model in the unloading state gives more accurate results for nonlinear behaviors.The pressure can be calculated using the bulk modulus and instantaneous density:

    Compression Equation of State loading and unloading is shown in Fig.5,at the point(Pn,ρn)the material reaches full compression.As can be seen in Fig.5,the unloading is done in a non-linear way,which is in the form of a dotted line,and their approximation is in the form of a solid line.The CTMDspeed of sound in the fully compressed state of matter [25].Because of the symmetry of the geometry and the loading of the problem,there is no need to simulate the problem in 3D and to reduce the cost of time,the problem is solved in two-dimensional axial symmetry,which shows in Fig.6 [25].

    3.Experiment results

    The explosive shock tube is used for experimental testing.Fig.7 shows a schematic diagram of the location of the specimens,explosive charge,and detonator in the test.In the first specimen,two sheets without core are adhesive together as shown in Fig.8.In subsequent experiments,the core of the specimens is 3,6,and 8 cm thick(see Table 5).

    Table 7 Specifications of experimental test samples.

    Fig.7.Schematic image of shock tube used in the experimental test.

    Fig.8.A specimen without core and measure specimen deformation.

    Fig.9.Specimens preparation.

    Table 6 Mechanical properties of AA3105 aluminum alloy[26].

    Fig.12.Explosion simulation for double layer aluminum foil.

    The dimensions of the used shock tube are shown in Fig.7.The inside of all samples has circularly meshed for ease of measurement.Fig.9 shows the prepare the specimens.Granular materials are put into cloth bags and then into cardboard cylinders.The aluminum plate used on both sides is AA3105,the mechanical properties of which are listed in Table 6 [26].The cardboard cylinder is adhesive to the aluminum sheets by adhesive.

    The explosive charge is of type C4 which is pasty-shaped and is shown in Fig.10.This charge is inserted into Teflon containers shown in Fig.11 left side and placed in the back of the shock tube as shown in Fig.11 in the middle.Teflon container back hole is for exit electric detonator wire ejection,and at the front is the location of the C4 explosive charge.In this study,an electric detonator with a score of 8(8 mm in diameter)is used,which has a power of 0.6 g of C4 explosive which should be considered in the simulation.The shape of the explosives is cylindrical.The diameter of the cylinder is twice its height and has been shaved inside the C4 and electric detonator holder due to the mass of the explosive charge.The sandwich panel specimens are fastened to the open end of the shock tube by 12 screws as shown in Fig.11 on the right side.The support is assumed to be completely clamped due to the high number of screws.After the test,the deformation of the outer sheet of the specimens at the mesh points(the mesh points on each line were 1 cm apart) is measured by a caliper as shown in Fig.8.The specifications of the experimental test samples are presented in Table 7.

    Fig.10.Type C4 explosive charge.

    Fig.11.Different views of the shock tube and explosive material holder.

    4.Results and discussion

    The simulation image for the non-core specimens with 4 and 5 g of the explosive charge is shown in Fig.12.

    Fig.13.Pressure wave before impacting the panel in a 1 mm gap.

    Fig.14.Explosion test of non-core specimens for 5 g of C4 explosive charge.

    The shock wave at 1 mm before impacting the sheet is visible for 4 and 5 g of C4 explosive in Fig.13.

    The result of the blast test by shock tube for the specimen with double layer aluminum foil without core is shown in Fig.14.

    Fig.15.Simulation result for the core specimen.

    In the second experiment,30 mm of sawdust was used as the core of the sandwich panel between the two sheets.The experiment was performed with 5 g of C4 explosive charge.In this case,the outer sheet (screwed to the open end of the shock tube) is completely torn and separated.The test result can be seen in Fig.21 on the left side.Fig.15 shows the simulation results of the Autodyn software for second specimen.It is evident that in the simulation,the specimen was also destroyed and the outer aluminum sheet was destroyed both from the stand and the center of the specimen.

    In the third simulation,a specimen with a core thickness of 30 mm mineral pumice and an explosive charge of 5 g was used,and the simulation results shown in Fig.16.In this case,the sheet is not torn and only deformed.The test result is shown in Fig.22 on the left side.

    Fig.16.Simulation result of a 30 mm thickness mineral pumice.

    Fig.17.Simulation result of an 80 mm thickness mineral pumice.

    Fig.18.Simulation result of a 60 mm-thickness sawdust.

    In the fourth test,80 mm of mineral pumice was used as the core of the sandwich panel and the experiment was performed with 5 g of the explosive charge.The simulation result and the amount of deformation of the center of the outer sheet are shown in Fig.17.According to Fig.17,the center of the sheet is about 15.5 mm deformed.The result of the empirical test is shown in Fig.22 on the right side.

    Fig.19.Is the result of simulating a specimen with an 80 mm thickness sawdust core.

    In the fifth test,60 mm of sawdust was used as the core of the sandwich panel and the experiment was performed with 3.5 g of the explosive charge.The result of the simulation is shown in Fig.18.According to Fig.18,the center of the outer sheet has displaced about 20 mm.The experimental test result is also shown in Fig.21 in the middle section.

    In the sixth simulation,80 mm of sawdust in the core and 3.5 g of explosive charge were used,according to the simulation result in Fig.19,the displacement of the center of the outer sheet is about 19 mm.Specimen of the empirical test shown in Fig.21 on the right side.

    The simulation result of the seventh sample of a sandwich panel with 60 mm of mineral pumice in the core and 5 g of the explosive charge is shown in Fig.20.According to Fig.20,the center of the sheet has been deformed about 18 mm.The seventh test specimen is shown in Fig.22.The results of the measured values from the empirical test and numerical simulation are presented and compared in Tables 8 and 9.The results of the experiments and numerical simulation were compared and there was a good agreement between these two modes of investigations,indicating the validity and accuracy of simulation assumptions.

    Fig.20.Simulation result of a 60 mm thickness core of mineral pumice.

    Fig.21.Status of specimens with sawdust core after explosion test.

    Fig.22.Status of specimens with mineral pumice core after the explosion test.

    The shape of the profile created by the experimental test and numerical simulation can also be compared.For this purpose,it is necessary to measure the displacement of all nodes of a line from a circular mesh in the samples.In the numerical simulation,instead of each of these nodes(mesh nodes),a gauge is placed in the same coordinates that report the displacement of that point on the sheet.Figs.23 and 24 show the comparison of sheet deformation profiles in simulation and test for different specimens with mineral pumice and sawdust,respectively.The main issue in the explosion problems is the pressure of the blast wave which hits the plate as shown in Fig.13 for 5 and 4 g of C4 charge.Each of these sandwich panel reduces the percentage of the wave.Overpressure and impulse are among the basic parameters that are considered in the issue of wave mitigation.Fig.25 shows the percentage of overpressure depreciation for the gauge on the center of the outer sheet in numerical simulation compared to the wave hitting the panel in shock tube (Fig.13).According to the diagrams,it is clear that increasing the core thickness of the specimens did not have much effect on increasing the damping of the blast wave overpressure;but it was observed in Tables 7 and 8 that the rate of deformation of the center of the plate decreased significantly with increasing core thickness(see Fig.4).

    Table 8 Comparison of Sheet Center Displacement in experimental testing and sawdust simulation (mm).

    Table 9 Comparison of Sheet Center Displacement in experimental testing and Mineral Pumice Simulation (mm).

    5.Investigating the effect of particle size and friction coefficient on shock wave depreciation

    After observing the results for experimental and numerical testing and comparing them,using numerical simulation,the effect of grain size parameters,and the coefficient of friction between grains were investigated.The particle size of granular material as a factor that is directly related to sample porosity can be one of the factors affecting the propagation speed of the wave.In this section,the collision of the blast wave with four samples with the properties of mineral pumice with a thickness of 60 mm and a grain diameter of 5,7,9,and 10 mm has been simulated.The pressure wave that hit the samples for 5 g of C4 is according to Fig.13,which was measured 1 mm before it hit the sample.The pressure diagrams compared to the time in the gauge on the back of the aluminum plate after passing the samples are observed for different grain sizes in Fig.26.The particle size of the granular material is directly related to porosity.By increasing the particle size,which means more porosity,the overpressure is further reduced.

    Fig.23.Comparison of deformation profile of numerical simulation and experimental panel testing with mineral pumice core.

    Fig.24.Comparison of the numerical simulation deformation profiles and the experimental test panel with sawdust core.

    Fig.25.Percent wave damping for different sandwich panel specimens.

    Fig.26.Effect of grain size on wave mitigation.

    The following are samples with mineral pumice properties and only with different coefficients of friction,the core thickness of the samples is 60 mm and the mass of C4 is 5 g.As mentioned earlier,the coefficient of friction for mineral pumice is between 0.3 and 0.7[20] Fig.27 shows changes in overpressure with a change in the coefficient of friction between the grains.

    According to Fig.27,it is clear that with increasing the friction coefficient,the amount of overpressure is slightly reduced,so that for friction coefficients of 0.3,0.5,and 0.7,the overpressure is equal to 1290 kPa,1480 kPa,and 1500 kPa,respectively is observed.

    Fig.27.The effect of the friction coefficient between grains in reducing shock wave.

    6.Conclusions

    This paper explored the ability of granular materials,such as sawdust and pumice to absorb blast energy was investigated.Pumice had a more significant effect at the same core thickness on blast energy absorption compared to sawdust for stronger waves,and its effect was noticeable even at lower thicknesses.Using these materials also had a great effect on reducing the sound of the blast.However,if the weight factor is also considered,sawdust had better performance,leading to more damped wave energy at a higher thickness,which can be highly effective.In the case of using a core with appropriate thickness,overpressure decreased to 88%,which indicates the ability of these materials (sawdust and pumice) to mitigate the blast wave.The deformation of the plate center varies significantly with the change in core thickness.Nevertheless,it is observed that the percentage of reduction in [shock/blast] wave overpressure has not changed significantly.The SPH method was used for simulation.The results revealed the appropriateness of this method in modeling granular materials under shock loading.An examination of the effects of grain size indicated a decrease in blast wave mitigation with a reduction in grain size.Investigating the coefficient of friction between particles exhibited that increasing the coefficient of friction between grains also led to a slight reduction in overpressure.

    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.

    久久亚洲真实| 777久久人妻少妇嫩草av网站| 国产精品1区2区在线观看. | 久久人人爽av亚洲精品天堂| 蜜桃在线观看..| 三上悠亚av全集在线观看| 热re99久久国产66热| 日韩中文字幕视频在线看片| 欧美黄色淫秽网站| 俄罗斯特黄特色一大片| 99精国产麻豆久久婷婷| 在线天堂中文资源库| 久久精品人人爽人人爽视色| 日韩视频一区二区在线观看| 中文字幕av电影在线播放| 精品第一国产精品| 十分钟在线观看高清视频www| 天堂动漫精品| 一本一本久久a久久精品综合妖精| www日本在线高清视频| a级毛片黄视频| 免费日韩欧美在线观看| 青青草视频在线视频观看| 国产一区二区在线观看av| 夜夜夜夜夜久久久久| 精品国产一区二区久久| 中文字幕人妻丝袜一区二区| 日韩欧美三级三区| 91麻豆精品激情在线观看国产 | 午夜成年电影在线免费观看| av天堂久久9| 99国产精品99久久久久| 欧美日韩国产mv在线观看视频| 天天躁日日躁夜夜躁夜夜| 国产亚洲精品一区二区www | 久久久精品94久久精品| 国产亚洲欧美精品永久| 丰满人妻熟妇乱又伦精品不卡| 久久久精品94久久精品| 国产精品电影一区二区三区 | 日本vs欧美在线观看视频| 国产单亲对白刺激| 亚洲综合色网址| 亚洲精品国产一区二区精华液| 久久久久久久久免费视频了| 精品卡一卡二卡四卡免费| 午夜日韩欧美国产| 99九九在线精品视频| 在线十欧美十亚洲十日本专区| a级毛片黄视频| 飞空精品影院首页| www.自偷自拍.com| 免费少妇av软件| 国产精品二区激情视频| 国产一区二区三区在线臀色熟女 | 两人在一起打扑克的视频| 一级片'在线观看视频| 久久热在线av| 国产男女内射视频| 国产日韩一区二区三区精品不卡| 日日夜夜操网爽| 免费日韩欧美在线观看| 侵犯人妻中文字幕一二三四区| 国产精品免费视频内射| 久久人人97超碰香蕉20202| 精品一区二区三卡| av视频免费观看在线观看| 91精品三级在线观看| 久久中文字幕一级| 69av精品久久久久久 | 久久影院123| 99国产精品99久久久久| 91精品国产国语对白视频| 脱女人内裤的视频| 一边摸一边抽搐一进一出视频| 建设人人有责人人尽责人人享有的| 丝瓜视频免费看黄片| 日韩欧美一区视频在线观看| 国产一区二区三区综合在线观看| 久久人妻福利社区极品人妻图片| 亚洲精品中文字幕一二三四区 | 成人18禁高潮啪啪吃奶动态图| 侵犯人妻中文字幕一二三四区| 美国免费a级毛片| 亚洲av成人一区二区三| 女性被躁到高潮视频| 999精品在线视频| www.自偷自拍.com| 亚洲国产毛片av蜜桃av| 人人妻人人爽人人添夜夜欢视频| 黑丝袜美女国产一区| 欧美成狂野欧美在线观看| 亚洲精品国产一区二区精华液| 99国产精品99久久久久| 一级a爱视频在线免费观看| 搡老熟女国产l中国老女人| 日日摸夜夜添夜夜添小说| 中文字幕精品免费在线观看视频| 精品一区二区三区av网在线观看 | 18在线观看网站| 人成视频在线观看免费观看| 国产精品成人在线| 国产区一区二久久| 首页视频小说图片口味搜索| 一级毛片电影观看| 久久热在线av| 亚洲专区字幕在线| 热re99久久精品国产66热6| 天天躁夜夜躁狠狠躁躁| 午夜91福利影院| 日韩免费高清中文字幕av| 搡老熟女国产l中国老女人| 色综合婷婷激情| 一边摸一边抽搐一进一出视频| 每晚都被弄得嗷嗷叫到高潮| 久久久精品免费免费高清| 亚洲七黄色美女视频| 亚洲成国产人片在线观看| 久久99热这里只频精品6学生| 曰老女人黄片| 国产日韩欧美亚洲二区| 国产精品免费一区二区三区在线 | 人人妻,人人澡人人爽秒播| 18禁国产床啪视频网站| 人成视频在线观看免费观看| 国产国语露脸激情在线看| 青青草视频在线视频观看| 亚洲精品成人av观看孕妇| 激情在线观看视频在线高清 | 美女福利国产在线| 久久99热这里只频精品6学生| 黄色怎么调成土黄色| 久久人人爽av亚洲精品天堂| 正在播放国产对白刺激| 国产男女内射视频| 自拍欧美九色日韩亚洲蝌蚪91| 一本—道久久a久久精品蜜桃钙片| 日韩一卡2卡3卡4卡2021年| 国产在线精品亚洲第一网站| a级毛片黄视频| 搡老岳熟女国产| 中文字幕人妻丝袜制服| 99久久精品国产亚洲精品| 国产精品欧美亚洲77777| 成人三级做爰电影| 国产精品亚洲av一区麻豆| 免费不卡黄色视频| 黄色片一级片一级黄色片| 免费观看a级毛片全部| 国产欧美日韩一区二区精品| 黄片播放在线免费| 黄色成人免费大全| 欧美变态另类bdsm刘玥| 这个男人来自地球电影免费观看| 男女高潮啪啪啪动态图| 免费一级毛片在线播放高清视频 | 黄色丝袜av网址大全| 免费在线观看视频国产中文字幕亚洲| 久久久久久免费高清国产稀缺| 涩涩av久久男人的天堂| 欧美日韩中文字幕国产精品一区二区三区 | 国产伦人伦偷精品视频| 中文欧美无线码| 色综合欧美亚洲国产小说| 中文字幕精品免费在线观看视频| 亚洲综合色网址| 亚洲av日韩精品久久久久久密| 国产伦人伦偷精品视频| 免费观看a级毛片全部| 国产精品亚洲av一区麻豆| 亚洲综合色网址| 成人永久免费在线观看视频 | 欧美午夜高清在线| 天天添夜夜摸| 99riav亚洲国产免费| 精品欧美一区二区三区在线| 黄色a级毛片大全视频| 欧美日韩精品网址| 色在线成人网| 久久精品熟女亚洲av麻豆精品| av在线播放免费不卡| 欧美变态另类bdsm刘玥| 日本精品一区二区三区蜜桃| 在线天堂中文资源库| 欧美激情高清一区二区三区| 国产精品 国内视频| 后天国语完整版免费观看| 免费不卡黄色视频| 搡老乐熟女国产| 亚洲自偷自拍图片 自拍| 一边摸一边抽搐一进一出视频| 午夜激情av网站| 亚洲专区国产一区二区| 亚洲中文字幕日韩| 每晚都被弄得嗷嗷叫到高潮| 国产精品影院久久| 久久精品aⅴ一区二区三区四区| 在线看a的网站| 精品欧美一区二区三区在线| 国精品久久久久久国模美| 大片免费播放器 马上看| 啦啦啦视频在线资源免费观看| 精品少妇内射三级| av免费在线观看网站| tube8黄色片| 人人妻人人澡人人爽人人夜夜| 国产黄色免费在线视频| 色尼玛亚洲综合影院| 国产老妇伦熟女老妇高清| 女人被躁到高潮嗷嗷叫费观| 欧美黑人精品巨大| 纵有疾风起免费观看全集完整版| 国产成人欧美在线观看 | 国产单亲对白刺激| 精品久久久精品久久久| 纯流量卡能插随身wifi吗| 另类亚洲欧美激情| 国产精品久久久久久人妻精品电影 | 美女扒开内裤让男人捅视频| 脱女人内裤的视频| 手机成人av网站| videos熟女内射| 一个人免费在线观看的高清视频| 亚洲熟女精品中文字幕| 丰满迷人的少妇在线观看| 成人特级黄色片久久久久久久 | 国产黄频视频在线观看| 久久影院123| 精品久久久精品久久久| 国产精品99久久99久久久不卡| 国产成人欧美| 搡老岳熟女国产| 久久av网站| 一个人免费在线观看的高清视频| 亚洲av欧美aⅴ国产| 国产成人免费观看mmmm| 亚洲avbb在线观看| 国产在线精品亚洲第一网站| 精品高清国产在线一区| 久久精品亚洲熟妇少妇任你| 午夜两性在线视频| 99久久人妻综合| 精品视频人人做人人爽| 亚洲久久久国产精品| 欧美激情高清一区二区三区| 汤姆久久久久久久影院中文字幕| 日韩一卡2卡3卡4卡2021年| 咕卡用的链子| 亚洲专区国产一区二区| 91字幕亚洲| 日韩欧美一区视频在线观看| 在线天堂中文资源库| 高清黄色对白视频在线免费看| 亚洲色图 男人天堂 中文字幕| e午夜精品久久久久久久| 午夜91福利影院| 在线观看人妻少妇| 露出奶头的视频| 1024视频免费在线观看| 国产成人一区二区三区免费视频网站| 欧美日韩亚洲国产一区二区在线观看 | 另类亚洲欧美激情| a在线观看视频网站| 正在播放国产对白刺激| 婷婷成人精品国产| 视频区欧美日本亚洲| 无人区码免费观看不卡 | 极品少妇高潮喷水抽搐| 亚洲av日韩精品久久久久久密| 女警被强在线播放| 色综合婷婷激情| 国产野战对白在线观看| 亚洲av日韩在线播放| 亚洲欧洲精品一区二区精品久久久| 一边摸一边做爽爽视频免费| 国产高清激情床上av| 久久久久国内视频| 国产亚洲精品第一综合不卡| 亚洲精品一二三| 亚洲av欧美aⅴ国产| 久久ye,这里只有精品| 欧美精品亚洲一区二区| 美女午夜性视频免费| 久久人妻福利社区极品人妻图片| 国产国语露脸激情在线看| 午夜福利视频精品| 韩国精品一区二区三区| 亚洲九九香蕉| 亚洲第一欧美日韩一区二区三区 | 成年版毛片免费区| 丁香欧美五月| 久久免费观看电影| 丰满迷人的少妇在线观看| 动漫黄色视频在线观看| 久久久精品免费免费高清| 1024视频免费在线观看| 久久99热这里只频精品6学生| 嫩草影视91久久| 国内毛片毛片毛片毛片毛片| 在线播放国产精品三级| 欧美成人午夜精品| 老汉色av国产亚洲站长工具| 成人国产av品久久久| 久久久久久久大尺度免费视频| 国产精品二区激情视频| 久久这里只有精品19| 亚洲欧美一区二区三区久久| 亚洲伊人色综图| 90打野战视频偷拍视频| 亚洲伊人久久精品综合| 一本大道久久a久久精品| 亚洲国产毛片av蜜桃av| 亚洲久久久国产精品| 国产精品 国内视频| 久久久久视频综合| 99国产精品99久久久久| 精品一区二区三区视频在线观看免费 | 午夜免费成人在线视频| 在线天堂中文资源库| 免费看十八禁软件| 777米奇影视久久| 王馨瑶露胸无遮挡在线观看| 老司机靠b影院| 性高湖久久久久久久久免费观看| e午夜精品久久久久久久| 国产av国产精品国产| 建设人人有责人人尽责人人享有的| 亚洲国产看品久久| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成a人片在线一区二区| 黄色毛片三级朝国网站| 国产精品久久久久久精品古装| 精品视频人人做人人爽| 男人操女人黄网站| 欧美国产精品va在线观看不卡| 午夜福利乱码中文字幕| 亚洲精品中文字幕一二三四区 | 日韩一区二区三区影片| 日韩视频在线欧美| 亚洲七黄色美女视频| 久久影院123| 午夜精品国产一区二区电影| 国产一区二区三区在线臀色熟女 | 久久久欧美国产精品| 久久精品成人免费网站| 法律面前人人平等表现在哪些方面| 一区二区三区激情视频| 午夜福利在线免费观看网站| 国产精品欧美亚洲77777| 在线观看免费视频日本深夜| 91成年电影在线观看| av网站在线播放免费| 成人特级黄色片久久久久久久 | 99riav亚洲国产免费| 91老司机精品| 午夜免费成人在线视频| 超色免费av| 正在播放国产对白刺激| 亚洲精品国产区一区二| 十八禁网站免费在线| 国产成人一区二区三区免费视频网站| 午夜福利影视在线免费观看| 午夜福利乱码中文字幕| 69精品国产乱码久久久| 不卡av一区二区三区| 美女高潮喷水抽搐中文字幕| 99香蕉大伊视频| 久久99热这里只频精品6学生| 中文字幕制服av| 啪啪无遮挡十八禁网站| 一级毛片精品| 国产野战对白在线观看| 丝袜人妻中文字幕| 热re99久久国产66热| 一区二区日韩欧美中文字幕| 黄色片一级片一级黄色片| 老熟妇乱子伦视频在线观看| 蜜桃在线观看..| 国产欧美日韩综合在线一区二区| 久久人妻av系列| 嫩草影视91久久| 亚洲精华国产精华精| 99国产精品99久久久久| 日韩视频在线欧美| 男女免费视频国产| xxxhd国产人妻xxx| 天天躁夜夜躁狠狠躁躁| 满18在线观看网站| 怎么达到女性高潮| 国产精品香港三级国产av潘金莲| 久久九九热精品免费| 精品少妇内射三级| av在线播放免费不卡| 久久香蕉激情| 757午夜福利合集在线观看| 亚洲第一青青草原| 国产成人免费观看mmmm| 精品福利观看| 女人久久www免费人成看片| 日韩中文字幕视频在线看片| 9热在线视频观看99| 久久精品人人爽人人爽视色| 大片电影免费在线观看免费| 另类精品久久| 满18在线观看网站| xxxhd国产人妻xxx| 国产成人系列免费观看| 精品高清国产在线一区| 精品一品国产午夜福利视频| e午夜精品久久久久久久| 亚洲国产欧美日韩在线播放| 国产一区二区三区综合在线观看| 大码成人一级视频| 99久久人妻综合| 最黄视频免费看| 欧美精品亚洲一区二区| 亚洲av日韩精品久久久久久密| 搡老乐熟女国产| 亚洲国产看品久久| 极品人妻少妇av视频| 每晚都被弄得嗷嗷叫到高潮| 久久精品成人免费网站| 不卡av一区二区三区| 日韩熟女老妇一区二区性免费视频| 欧美精品av麻豆av| 亚洲第一av免费看| 最黄视频免费看| 欧美精品亚洲一区二区| 一个人免费在线观看的高清视频| 色尼玛亚洲综合影院| 涩涩av久久男人的天堂| 亚洲色图 男人天堂 中文字幕| a级毛片黄视频| 午夜福利免费观看在线| 一本—道久久a久久精品蜜桃钙片| 动漫黄色视频在线观看| 一夜夜www| 亚洲精品久久午夜乱码| 在线观看免费视频网站a站| 久久九九热精品免费| 国产精品一区二区精品视频观看| 亚洲精品国产精品久久久不卡| 午夜久久久在线观看| 一本久久精品| 亚洲熟女精品中文字幕| 久久 成人 亚洲| 波多野结衣av一区二区av| av网站免费在线观看视频| 色视频在线一区二区三区| 十八禁高潮呻吟视频| 交换朋友夫妻互换小说| 亚洲情色 制服丝袜| av天堂在线播放| 视频在线观看一区二区三区| 精品高清国产在线一区| 国产不卡一卡二| 国产黄频视频在线观看| 精品久久久精品久久久| 91麻豆精品激情在线观看国产 | 日韩免费高清中文字幕av| 他把我摸到了高潮在线观看 | 丁香欧美五月| tocl精华| 日韩免费高清中文字幕av| 亚洲视频免费观看视频| 人人妻人人澡人人看| 一边摸一边抽搐一进一出视频| 男女床上黄色一级片免费看| 亚洲免费av在线视频| 亚洲成人免费电影在线观看| 18禁美女被吸乳视频| 制服人妻中文乱码| 嫩草影视91久久| 母亲3免费完整高清在线观看| 国产99久久九九免费精品| 国产xxxxx性猛交| 欧美日韩视频精品一区| 国产精品熟女久久久久浪| 亚洲久久久国产精品| 可以免费在线观看a视频的电影网站| 午夜两性在线视频| 国产成人精品久久二区二区91| 精品少妇内射三级| 国产亚洲欧美在线一区二区| 免费在线观看日本一区| 人妻一区二区av| 成人18禁在线播放| 91成人精品电影| 青草久久国产| 日韩免费高清中文字幕av| 美女高潮喷水抽搐中文字幕| 成年女人毛片免费观看观看9 | 高潮久久久久久久久久久不卡| 久久久精品免费免费高清| 麻豆国产av国片精品| 亚洲熟女毛片儿| 国产成人影院久久av| 久久精品亚洲熟妇少妇任你| 精品久久久久久久毛片微露脸| 色综合婷婷激情| 国产精品一区二区在线观看99| 久久99热这里只频精品6学生| 两个人免费观看高清视频| 欧美日韩亚洲国产一区二区在线观看 | 免费看十八禁软件| 99国产精品一区二区三区| 国产免费现黄频在线看| 国精品久久久久久国模美| 99在线人妻在线中文字幕 | 精品国产乱子伦一区二区三区| 一二三四社区在线视频社区8| 久久人人爽av亚洲精品天堂| 高清视频免费观看一区二区| 欧美久久黑人一区二区| 亚洲精品一二三| 午夜成年电影在线免费观看| 精品亚洲成a人片在线观看| av福利片在线| 欧美日韩亚洲高清精品| 视频在线观看一区二区三区| 手机成人av网站| 香蕉丝袜av| 香蕉久久夜色| 久久中文字幕人妻熟女| 一区二区三区国产精品乱码| 亚洲精品美女久久av网站| 久久久精品94久久精品| 国产又色又爽无遮挡免费看| 一进一出抽搐动态| 欧美日韩av久久| 天堂中文最新版在线下载| 精品卡一卡二卡四卡免费| 国产在线视频一区二区| 50天的宝宝边吃奶边哭怎么回事| kizo精华| 国产精品一区二区免费欧美| 18禁国产床啪视频网站| 可以免费在线观看a视频的电影网站| 日韩制服丝袜自拍偷拍| 高清毛片免费观看视频网站 | 久久久久久久国产电影| a在线观看视频网站| 亚洲午夜精品一区,二区,三区| 久久中文字幕人妻熟女| 日日摸夜夜添夜夜添小说| 欧美黄色片欧美黄色片| 黄色丝袜av网址大全| 国产淫语在线视频| 最新在线观看一区二区三区| 久久人人97超碰香蕉20202| 69精品国产乱码久久久| 天堂俺去俺来也www色官网| 亚洲五月婷婷丁香| 亚洲精品一卡2卡三卡4卡5卡| 18禁美女被吸乳视频| 久久精品国产综合久久久| 9热在线视频观看99| 久久久久久久精品吃奶| 一本大道久久a久久精品| 国产免费福利视频在线观看| 美女福利国产在线| 国产精品九九99| 看免费av毛片| 精品少妇黑人巨大在线播放| 久久国产精品人妻蜜桃| 少妇被粗大的猛进出69影院| 亚洲国产欧美在线一区| 99re在线观看精品视频| 欧美另类亚洲清纯唯美| 欧美在线黄色| 人妻久久中文字幕网| 国产免费av片在线观看野外av| 午夜免费鲁丝| 成人免费观看视频高清| 国产成+人综合+亚洲专区| 69av精品久久久久久 | 91成年电影在线观看| 久久青草综合色| 日韩熟女老妇一区二区性免费视频| 色视频在线一区二区三区| 啦啦啦 在线观看视频| 精品少妇内射三级| 欧美激情高清一区二区三区| 一区二区三区国产精品乱码| 肉色欧美久久久久久久蜜桃| 色婷婷久久久亚洲欧美| 最近最新免费中文字幕在线| 麻豆乱淫一区二区| 后天国语完整版免费观看| 成人国产一区最新在线观看| 成人国语在线视频| av欧美777| 国产在线精品亚洲第一网站| 国产午夜精品久久久久久| 搡老熟女国产l中国老女人| 一个人免费在线观看的高清视频| 欧美国产精品va在线观看不卡| 电影成人av| 飞空精品影院首页| 一本—道久久a久久精品蜜桃钙片| avwww免费| 两个人看的免费小视频| 捣出白浆h1v1| 一本久久精品| 亚洲一区二区三区欧美精品| 蜜桃国产av成人99| aaaaa片日本免费| 美女午夜性视频免费| 久久热在线av| 亚洲av日韩精品久久久久久密| av不卡在线播放| 这个男人来自地球电影免费观看| 亚洲专区中文字幕在线| av又黄又爽大尺度在线免费看| 精品一区二区三卡|