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

    Dynamic crushing behaviors and enhanced energy absorption of bioinspired hierarchical honeycombs with different topologies

    2023-05-06 01:59:26XinchunZhngNnnnLiuChochoAnHexingWuLiKemingHo
    Defence Technology 2023年4期

    Xin-chun Zhng ,Nn-nn Liu ,Cho-cho An ,He-xing Wu ,N Li ,Ke-ming Ho

    a Department of Mechanical Engineering,North China Electric Power University,Baoding,071003,China

    b Hebei Key Laboratory of Electric Machinery Health Maintenance & Failure Prevention,North China Electric Power University,Baoding,071003,China

    c School of Civil Engineering,Northeast Forestry University,Harbin,150040,China

    d Department of Intelligent Engineering,Hebei Software Institute,Baoding,071000,China

    Keywords:Bio-inspired honeycombs Crushing load uniformity Energy absorption Hierarchical topologies Finite element analysis

    ABSTRACT In order to pursue good crushing load uniformity and enchance energy absorption efficiency of conventional honeycombs,a kind of bio-inspired hierarchical honeycomb model is proposed by mimicking the arched crab shell structures.Three bio-inspired hierarchical honeycombs (BHHs) with different topologies are designed by replacing each vertex of square honeycombs with smaller arc-shaped structures.The effects of hierarchical topologies and multi-material layout on in-plane dynamic crushings and absorbed-energy capacities of the BHHs are explored based on the explicit finite element (FE) analysis.Different deformation modes can be observed from the BHHs,which mainly depend upon hierarchical topologies and impact velocities.According to energy efficiency method and one-dimensional(1D)shock theory,calculation formulas of densification strains and plateau stresses for the BHHs are derived to characterize the dynamic bearing capacity,which is consistent well with FE results.Compared with conventional honeycombs,the crushing load efficiency and energy absorption capacity of the BHHs can be improved by changing the proper hierarchical topology and multi-material layout.These researches will provide theoretical guidance for innovative design and dynamic response performance controllability of honeycombs.

    1.Introduction

    Cellular solids have been widely used in many engineering applications of aerospace,automobile,packaging engineering,defence technology and so on due to their light weight and excellent performance ratio[1—5].With the further research on dynamic crushings and energy absorption capacities of cellular solids,the dynamic topological and hierarchical design have attracted great attention of the scholars at home and abroad[6—9].Therefore,how to explore the optimally hierarchical topology and then realize their controllable design of dynamic crushing performance according to engineering demands are one of important goals of sandwich structures and materials.

    In recent years,hierarchical honeycombs have been developed based on conventional honeycombs with different topologies like triangle,square and hexagon [10—12].Many novel hierarchical structures with enhanced in-plane impact resistance and energy absorption capacities have been proposed by lots of researchers[11,13—15].And they indicated that introducing the hierarchy into conventional structures can effectively reduce the whole deformation of honeycombs.For example,based on the concept of hierarchical structures,An et al.[11]introduced self-similar hexagonal hierarchy at the vertex of cell structure to improve the dynamic strengths of honeycombs.Nian et al.[13]proposed an innovative self-similar graded honeycomb-filled thin-walled structure to systematically study the bending behaviors.They found that this kind of structure can present better energy absorption characteristics and crashworthiness performance than the corresponding uniform honeycomb-filled circular tube.Zhang et al.[14]obtained self-similar hierarchical circular tubes with higher average impact force and specific energy absorption by choosing suitable configurations.Lu et al.[15]proposed a novel hierarchical structure to effectively enhance energy absorption capacity of chiral structures.These studies reveal that mechanical behaviors of hierarchical honeycombs are affected by the low-order hierarchical structure and the order of hierarchical cell [16—18].Inspired by multi-level biological structure in nature,Xu et al.[16]constructed a self-similar layered hexagonal column with multi-order characteristics to improve the crashworthiness.Hou et al.[17]and Li et al.[18]developed self-similar multi-order hierarchical models to predict the dynamic response characteristics.Carrying out FE analysis on the dynamic collapse of self-similar hierarchical corrugated sandwich structures,Zhang et al.[19]found that the main deformation modes of collapse mechanism diagram and Euler buckling region of small compression bar will obviously change with increasing impact velocity.Du et al.[20]designed a novel hierarchical thermoplastic composite honeycomb cylindrical structures and experimentally investigated the mechanical response and energy absorption characteristics.Rahman et al.[21]combined hierarchical structure with cellular metamaterials to construct 2-order and 3-order positive Poisson's ratio hierarchy to explore the influence of defects on hierarchical structures.Kandasamy et al.[22]predicted the transient vibration behavior of hierarchical honeycomb plates with surface bonded actuators by controlling gain and hierarchical level.Although extensive researches on the mechanical properties and energy absorption capacities of hierarchical honeycombs have been carried out,most of them are directed toward the small deformation.Further research is needed to contribute to maintain good crushing load uniformity while enhancing in-plane energy absorption capacities of hierarchical honeycombs under large deformation.However,it should be pointed out that the maximum stress is mainly concentrated at the apex of conventional honeycomb where plastic hinge generates under in-plane impacts,which is not conducive to obtain good crushing load uniformity.Hence,how to select and design hierarchical topology is extremely important for the dynamic performance controllability of hierarchical honeycombs.

    When subjected to impact loadings,the bio-inspired honeycombs based on structural bionics can absorb more energy due to their excellent impact resistance,which has the potential engineering application for the large-scale manufacturing of energyabsorbing structure protection and attracts a great deal of attention.With the development of bionic technology,many novel biomorphic honeycombs have been designed based on structural bionics aiming to improve the mechanical properties and enhance energy absorption capacities of conventional honeycombs[23—28].Xiang et al.[23]designed a new bio-inspired multi-layered structure which can effectively enhance energy absorption and crashworthiness by imitating the structure characteristics of human skeleton.Inspired by the shape of tortoise shell,Zhang et al.[24]proposed a bio-inspired re-entrant arc-shaped honeycomb to improve the absorbed-energy abilities of honeycombs.These results show that good crushing load uniformity is given in bioinspired honeycombs compared with traditional concave ones under in-plane impacts.On the basis of the concept of bionic design,Yang et al.[25]devised a series of bidirectional corrugated plates and revealed that impact resistance of corrugated structure can be effectively improved by reasonably drafting wave amplitude and wavelength.Ma et al.[26]constructed a biological materials with horseshoe microstructures to tailor its mechanical responses.Ha et al.[27]developed a bio-inspired sandwich wavy panel by imitating a woodpecker's beak,which expressed superior specific energy absorption than conventional sandwich panel.Ufodike et al.[28]proposed a bamboo biomorphic structure for in-plane energy absorption improvement and they indicated that the bio-inspired design showed significantly higher densification strain.Although bionic structure has good crashworthiness,the researches on the dynamic dominant mechanism and crushing performance adjustable design of the macro-/micro-structure with bio-inspired cellular materials under large deformation are far from enough.

    Based on the inspiration of arched crab shell structures with surpassing impact resistance,three BHHs with different topologies are designed by replacing each vertex of regular square honeycombs(RSHs)with smaller arc-shaped structures in this paper.The effects of hierarchical topologies and multi-material configuration on in-plane dynamic responses and energy absorption capacities of the BHHs are numerically discussed.A comparative study between the BHHs and RSHs are carried out.Empirical formulas of densification strains and plateau stresses to characterize the dynamic properties are derived.Lastly,energy absorption efficiencies of the BHHs with different topologies and material configuration parameters are also explored.

    2.Methodologies and models

    2.1.Geometric configuration

    In the nature,natural organisms can provide many excellent mechanical properties and energy absorption capacities after a long period of natural selection,such as turtle shell [24],horseshoe pattern [26],the beaks of woodpeckers [27]and the crab chela structures[29].It is interesting that all these configurations have a common feature,that is,they have the similar arched structures.The arched structure of crab chela can sustain strong impact loadings while cutting objects,which can provide excellent impact resistance for crab chela.Inspired by the arched crab shell structures,a kind of BHH model is proposed in this paper.Fig.1 shows the detail design strategy and representative volume element for a bio-inspired honeycomb.As shown in Fig.2,three BHHs are designed by replacing the nodes of the RSHs (Fig.2(a)) with different arc-shaped hierarchical topologies: (a) the BHH with circle cell (BHH—C,Fig.2(b)): the nodes of RSHs are replaced by circular structures;(b)the BHH with diamond cell(BHH-D,Fig.2(c)):the nodes of the RSHs are replaced by diamond-shaped structures which arc-edge of the circle is folded along the side of inscribed square with the circle to form a concave quadrilateral and the apex of the quadrilateral is connected with the ligament;(c) the BHH with re-entrant structure (BHH-R,Fig.2(d)): rotating the concave quadrilateral 45°in Fig.2(c) and connecting the concave arched edge with the ligament.With the evolvement of metal 3D printing rapid prototyping technology,the preparation difficulty of the BHHs will be greatly reduced and the manufacture quality can be effectively controlled.

    Unit cell geometries of regular and bio-inspired honeycombs are displayed in Fig.3,wherelis the half length of unit cell with a RSH,erefers to the ligament length,randtare respectively the radius of the arched cell wall and the ledge thickness.Both the ligament and the introduced hierarchical cell thickness remain the same.In the filling,the half cell length of hierarchical honeycomb is a constant withl=2.4 mm.The relationship betweenr,eandlcan be given as

    Relative density is a key indicator to describe the mechanical behavior of honeycombs.According to the definition of cellular material theory(CMT)[30],the relative densities of the BHHs with different topologies can be formulated as

    Fig.1.Design strategy and representative element for a BHH.

    Fig.2.Cross-sectional geometries of regular and bio-inspired honeycombs.

    Fig.3.Unit cell geometries of regular and bio-inspired honeycombs.

    where ΔρBHH-C,ΔρBHH-Dand ΔρBHH-Rare respectively the relative density of the BHH—C,BHH-D and BHH-R specimens,ρ* presents the density of the BHHs,and ρsindicates the density of matrix material.In particular,when the ligament length of cell wall with the BHH—C is 0,its relative density can be specified as

    which is the same as the relative density of circular honeycomb[31].

    2.2.FE models and the validation

    In order to disclose the dynamic behaviors and energy absorption efficiencies of the BHHs,a FE model was established using ABAQUS/EXPLICIT.Fig.4 illustrates the diagrammatic sketches of the calculating model for the BHHs under in-plane crushing.The specimens are located between two rigid plates.The fixed rigid plate is on the bottom,and the impact rigid plate is crushing a specimen alongy-direction with a constant velocity [32,33].The matrix material are assumed to be elastic-perfectly plastic.The strain rate sensitivity is very low and ignored [12,24,32—34].To explore the effects of hierarchical topologies on dynamic crushing of the BHHs,the matrix material remains the same with typical aluminum alloy (Material-1) [34]and the corresponding material parameters are listed in Table 1.As for the multi-material layout models,the specimens are made of aluminum alloy [34]and aluminum (Al) (Material-2) [35],the material properties are also listed in Table 1.The materials of impact and fixed plates are assumed to be rigid.The width of the specimen isLxand the length isLy.Researches have shown that size effect on the dynamic response characteristics can be ignored when the number of filled cells is more than 5 [31—35].Here,the cell number is 13 in thexdirection and 14 in they-direction,which could effectively capture the dynamic mechanical behaviors.

    Fig.4.Diagrammatic sketches of the calculating model for bio-inspired honeycombs.

    Some researches have shown that shell element could ensure sufficiently reliable results[32—35].In this study,the specimens are also selected as shell element with S4R.The impact and fixed plates are employed by the discrete element with R3D4.For the consideration of convergence and calculation precision,five integration points are used along the thickness direction of each shell element.Mesh convergence analysises on element size are conducted firstly to test the validation of FE models.The best choice of meshing length is 0.21 mm for the BHHs according to the mesh convergence and computing efficiency conservatively.The out-of-plane displacement of FE model is limited during the crushing to prevent the out-of-plane buckling,and the specimen boundaries on the left and right edges are free.In this simulation,the fracture failure of matrix material was ignored.Moreover,automatic surface-to-surface contact was adopted to define the contact between the specimens and two rigid plates,and the friction coefficients were assumed to 0.2 [35].And general contact was defined for each surface of the specimen.The out-of-plane thickness of the BHHs was taken as unit thickness.To investigate the effects of hierarchical topologies,the BHHs can be obtained by changinge/l.Table 2 provides more details on the geometries and dimensions of the specimens.

    Table 1 Matrix material parameters of the BHHs.

    Table 2 Dimensions of the BHHs (unit: mm).

    The validation of FE model was provided,the same model as in Ref.[31]was established to compare the dynamic response curve of a regular circular honeycomb.The nominal stress and strain of the specimens are respectively defined as

    whereF(ΔLy) is the reaction force on impact rigid plate when the compression displacement alongy-direction is ΔLy,A0is the original cross-section area of the specimen before compression,bis the specimen thickness in out-of-plane displacement (in this paper,b=1 mm),LxandLyare respectively the width and the length of the specimen.At present,lots of researches on the dynamic responses of cellular structures have been explored based on FE simulation[10—12,24,31—35].It is proved that the FE modeling method is validated to be effective and accurate.Sun et al.[31]studied inplane crushing behavior and energy absorption performance of multi-layer regularly arranged circular honeycombs by using ANSYS/LS-DYNA.For the BHH—C specimen in this paper,whene=0,the model is the same with multi-layer regularly arranged circular honeycomb in Ref.[31].Fig.5 shows the dynamic response curve of a regular circular honeycomb atv=50 m/s.When matrix material,model size,impact velocity and boundary conditions are all the same,the nominal stress-strain curve of FE model in this paper is in good agreement with that of Ref.[31].But there is the difference between two curves at the beginning of impact.On the one hand,two different simulation products are exployed to investigate the dynamic responses.Two softwares will cause the error difference at the beginning of impact.On the other hand,more output datas are captured in our simulation to achieve higher calculation accuracy.It should be noted that both the trends and amplitudes between two curves at the beginning of impact are the same.It is seen that FE model established in this paper is reliable and the same method will be used to explore the dynamicbehaviors and absorbed-energy capacities of the BHHs with different topologies for subsequent study.

    Fig.5.The nominal stress-strain curves for the BHH—C with e/l=0 at v=50 m/s.

    2.3.Basic properties and evaluation indexes

    Many researches indicated that impact velocity was an important factor affecting the dynamic response characteristics of cellular solids[32—38].Different stress waves are generated at the proximal end under different impact velocities.When impact velocity is higher than wave trapping velocity,local deformation of cellular solids will begin to occur.According to the ‘wave trapping’ theory[33],wave trapping velocity can be specified as the first critical velocity,which can be written as

    where εcris initial strain of cellular solid when the crushing stress reaches initial peak stress,and σ′(ε) refers to elastic modulus in linear stage.For Δρ=0.115,the first critical velocity for the BHH—C,BHH-D and BHH-R specimens are estimated as 9.4,12.7 and 8.5 m/s,respectively.With the increase of impact velocity,the deformation gradually tends to be localized.Local collapse zone begins to propagate from the impact end to the fixed end in the form of shock waves.The critical velocity at which shock waves propagate in cellular solid is considered as the second critical velocity[32,37,38].

    where σ0represents the quasi-static plateau stress and εdrefers to the densification strain.When Δρ=0.115,the second critical velocity for the BHH—C,BHH-D and BHH-R are 62.6,65.1 and 53.9 m/s,respectively.In order to explore the effect of impact velocity on dynamic behaviors,impact velocity ranges from 1 to 120 m/s.Based on the above definition of critical velocity,the range of impact velocity is divided into three sections.Whenv

    Fig.6.Nominal stress-strain curve(black line),energy efficiency-strain curve(red line)and the corresponding specific energy absorption curve for a BHH.

    Fig.6 shows the typical dynamic response curve(black line)of a BHH,which contains three stages: the elastic stage with the maximum peak stress σp,the plateau stage with the plateau stress σmand the compressed densification stage.σmis a significant index to describe energy absorption capacity of honeycombs,which can be obtained by

    where εcris the initial strain corresponding to the peak stress.Our researches show that εcris greatly small during the crushing,a valure 0.02 is assumed for initial strain to ensure the higher calculation accuracy.εdis the densification strain.When the strain is over the densification strain,the stress of honeycombs will increase sharply.It is obvious that plateau stress is directly correlated to εd.The perfect densification strain is the porosity of a honeycomb,that is εd=1-Δρ.However,it should be mentioned that εddepends upon impact velocities and topological structures[32—35].In order to avoid the randomness of artificial selection,εdcan be obtained by using the energy efficiency method,which can be given as

    whereEis energy absorption efficiency of a honeycomb,which can be defined by the ratio of the strain energy to corresponding nominal stress as

    When honeycombs are used for energy-absorbing components,good crushing load uniformity should be maintained during energy dissipation.The maximum peak force should be lower than the critical one of destroying the protected components.Furthermore,in order to prevent the dynamic damage to the protected structure,the reaction force of honeycombs should be kept as constant as possible.There are many indexes to describe the crushing load uniformity of honeycombs [12,24],which can be given as

    where CLE is the crushing load efficiency,σmis the plateau stress,Δσ refers to the average stress fluctuation and σpis the maximum peak stress in the dynamic response curve.For an ideal energy absorber,a smaller σpand a higher σmare needed.Due to σm≤σp,that is,0

    Moreover,a non-dimensional indicator (NLC) is used to characterize the effective load-carrying capacity of honeycombs [39],which can be obtained by

    whereFmis the average reaction force on the impact plate,M0is the fully plastic bending moment per unit length,which can be formulated as

    in which σysis yield stress of matrix materials andtrepresents the edge thickness of honeycombs.Moreover,undulation of loadcarrying capacity (ULC) is another dimensionless indicator to reflect the undulation degree of impacting stress around plateau stress during the plateau stage,which can be defined as

    where σmis plateau stress,εdis densification strain.

    Due to the demands of lightweight structure design,the sensitivity of material weight should be considered in many energyabsorbing components.A specific energy absorption (SEA) is a crucial index to evaluate the energy absorption capacity of sandwich structures[24,33,35],which can be defined as

    whereEVis the densification strain energy which can be determined by,Δρ is the relative density of the specimen and ρsis the density of matrix material.Furthermore,is also adopted to characterize the influence of impact velocity on energy absorption capacity [40].It normalizes the plastic energyU0by the perfect static plastic energy absorption,which can be defined as [41].

    where σ0is the static plateau stress,A0andLyare the originally sectional area and length of the specimen,respectively.

    3.Numerical results and discussion

    3.1.Deformation characteristics

    Many researches have shown that different collapsed modes are observed in honeycombs with different topologies under impacts.And the deformation mechanisms for hierarchical honeycombs depend upon cell micro-structures [11,12,15,24].In order to study the influences of hierarchical topologies and impact velocities on the deformation mechanisms of the BHHs,the deformation modes for the BHH—C,BHH-D and BHH-R with different impact velocities are discussed in this section.

    Fig.7 displays the macro-scopic deformation of three BHHs at ε=0.35.For low-velocity crushing (v=5 m/s

    3.2.Densification strain and plateau stress

    Plateau stress (σm) and densification strain (εd) are two key indexes in describing the dynamic behaviors of honeycombs.Fig.10 shows the effects of impact velocities on the densification strains for the BHHs with different topologies.For the same hierarchical honeycombs,εdwill increase with impact velocity [35].It can be predicted by the following equation with impact velocity [12].

    Fig.7.Deformation modes of the BHHs at ε=0.35 under different impact velocities.

    where bothkandware constants,which depend on hierarchical topologies.The fitting values forkandwof the BHHs are obtained from energy efficiency method and FE results.Therefore,the empirical formulas of εdfor the BHHs with different topologies can be predicted as

    According to the theory of 1D shock waves,Reid and Peng [42]proposed a simple formula of plateau stress for cellular solid,which is related to relative density and impact velocity

    whereAis a parameter which can be decided by cell topologies of cellular solids.Gibson and Ashby indicated that Eq.(20)can also be applied to other cellular solids [29].Using the least square fitting method and introducing Eq.(19)into Eq.(20),the dynamic plateau stress at the impact end for the BHHs with different topologies are obtained,which can be written as

    Fig.8.Nominal stress-strain curves of the BHHs with different velocities for e/l=0.583.

    Fig.9.Nominal stress-strain curves of a BHH-R with different e/l at v=20 m/s.

    Fig.10.Variation of densification strains for BHHs with respect to impact velocities.

    Fig.11.Variation of plateau stresses for the BHHs with respect to impact velocities.

    Fig.12.Δσ and CLE for the BHHs with different e/l at v=20 m/s.

    3.3.Crushing load uniformity

    As an excellent energy absorbing structure,the higher crushing load efficiency(CLE),and the lower average stress fluctuation(Δσ)and peak stress(σp)are required.Fig.12 illustrates the Δσ and CLE for the BHHs with differente/latv=20 m/s.For comparison,the corresponding results of the RSHs are also given.Obviously,the CLEs of the BHHs are higher than that of the RSHs.They first increase to the maximum value ate/l=0.417,and then follow a decline with the increasing ofe/l.The CLE of the BHH-R is the highest,which is 71.89%,164.76% and 521.26% higher than the BHH—C,BHH-D and RSH specimens,respectively.As shown in Figs.12 and 13,the mean stress fluctuation and maximum peak stresses of the BHHs could achieve the minimum value ate/l=0.417.For the BHH-R,σpis larger than that of other hierarchical honeycombs (as shown in Fig.13).The σpof the BHH—C is the lowest compared with other specimens for differente/l.This indicates that the crushing load uniformity of the BHHs can be improved by changing the propere/l.The BHH—C specimens are more suitable for energy absorbing components than other honeycombs when considering the crushing load uniformity under large deformation.

    Fig.13.Variation of σp for the BHHs with different e/l at v=20 m/s.

    Fig.14.CLE for the BHHs with multi-material configuration.

    In order to explore the influence of multi-material configuration on crushing load uniformity of the BHHs,two kinds of matrix materials (as shown in Table 1) are considered to establish four layout models,which are the H—H,S—H,H—S and S—S,respectively.Take H—S model for example,S (Material-2) refers to matrix material of the introduced joint hierarchy,and H (Material-1) is the corresponding ligament material.Similarly,S—H indicates the material of the introduced joint hierarchy as Material-1.H—H and S—S define the same matrix material of the whole model as Material-1 and Material-2,respectively.Additionally,in order to compare multi-material configuration between the BHHs and RSHs,two materials are also considered in the RSHs,that is the RSH(H—H) and the RSH(S—S).

    Whenv=20 m/s ande/l=0.334,the CLE,Δσ and σpfor the BHHs with different material layouts are shown in Figs.14—16.In contrast,the corresponding results of the RSHs are also given.It is clear that the BHHs with H—S configuration can obtain the higher CLE (as shown in Fig.14).Especially,the BHH-R with H—S has higher CLE than other hierarchical models,which is 48.07% and 63.91%higher than BHH—C and BHH-D,respectively.The CLE of the BHHs with multi-material layout is higher than that of the RSHs,as shown in Fig.14.However,the Δσ and σpof the BHH-R are always greater than those of BHH—C and BHH-D(Figs.15 and 16),which is mainly related to the deformation mechanism.Moreover,the BHH—C with multi-material configuration has the low stress fluctuation,which is even lower than the RSH with S—S.It can be observed from Fig.16 that σpof the BHHs with H—S are more advantageous for energy absorption than that of other layouts due to the lower peak stresses.Therefore,the BHHs with H—S perform well in all aspects for the crushing load uniformity.The reason is that due to the introduced joint hierarchy of low yield stress,the plastic deformation is first observed in introduced hierarchical structures and then the ligaments,which causes the whole deformation of the specimens.The BHH-R specimens exhibit excellent performance in CLE under large deformation,but they exhibit the higher Δσ and σp.It can be seen from the above analysis that compared with the RSHs,the crushing load uniformity of the BHHs can be improved if the proper multi-material layouts are chosen.

    As an excellent energy absorbing component,the higher dynamic load-carrying capacity and smaller undulation of loadcarrying should be considered.Fig.17 shows the variation of the NLC for the BHHs with different topologies atv=20 m/s.The NLC of the BHH-R is obviously higher than that of the BHH—C,BHH-D and RSH specimens,which is consistent with Fig.12.However,for the BHH—C and BHH-D,the NLCs are relatively stable (they will fluctuate around 30)due to their insensitivity to the variation ofe/l.The NLC of the BHH-D is slightly higher than that of the BHH—C.It can be seen from Fig.17 that the NLCs of the BHHs are obviously higher than that of the RSHs.Research results show that the NLC of the BHHs can be improved and effectively controlled through the proper choice ofe/land the introduced hierarchical topologies.

    Fig.15.Δσ for the BHHs with multi-material configuration.

    Fig.16.σp for the BHHs with multi-material configuration.

    Fig.18 shows the variation of ULCs for the BHHs with respect toe/latv=20 m/s.Along with the increase ofe/l,the ULCs of the BHHs will decrease first to the minimum ate/l=0.417,and then followed by an upward trend.The ULC of the BHH-R is higher than the RSH.It is seen that all BHHs can obtain better ULC ate/l=0.417,and the BHH—C has lower ULC than other specimens,which is 117.17%,108.29% and 16.59% lower than the RSH,BHH-R and BHH-D,respectively.Fig.19 shows the effects of multi-material layout on the ULC of the BHHs with different topologies.For the BHHs with H—S,the introduced joint hierarchy first deforms due to less yield stress,which prompts the whole deformation of the BHHs.It is the reason that the ULC of the BHHs with H—S can obtain the optimal value.Research results show that selecting the appropriatee/land using multi-material layout can decrease the undulation degree of impact stress for the BHHs.

    Fig.17.Variation of NLC for BHHs with respect to e/l at v=20 m/s.

    Fig.18.Variation of ULC for the BHHs with respect to e/l at v=20 m/s.

    3.4.Energy absorption characteristics

    Based on Eq.(17),the variation offor the BHHs with respect to impact velocity is shown in Fig.20.A higher value ofindicates that inertial effects have more significant enhancement on energy absorption of honeycombs.increases with the increase of impact velocity.The BHH-R has the higherthan the BHH—C and BHH-D.Under low-velocity impacting(v

    Fig.19.ULC for the BHHs with multi-material configuration for e/l=0.334 at v=20 m/s.

    Fig.20.Variation of for the BHHs with respect to impact velocity.

    Based on Eq.(16) and energy absorption efficiency method,Fig.21 shows the effects of impact velocity on the SEA of the BHHs with the same cell-wall thickness.As a comparison,the SEA of the RSH is also given.With the increase of impact velocity,the SEAs of the BHHs are obviously improved.For the same impact velocity,the SEA of the BHH-R is higher than that of other hierarchical topologies and the RSH,which is consistent with dynamic load-carrying capacity(as shown in Fig.17).Fig.22 shows the variation of the SEA for the BHHs with respect toe/latv=20 m/s.Compared with the RSH,due to the introducing of joint hierarchical topologies,the BHHs have higher SEA than the RSHs.It is noted that the SEAs of the BHH—C and BHH-D are almost equal for the samee/l,and it has little sensitivity one/l.That is because these two models have the same relative density for the samee/l.However,with the increase ofe/l,the SEA of the BHH-R displays a manifest decrease,which will receive a higher value ate/l=0.334.In addition,the SEA of the BHH-R is higher than that of the BHH—C and BHH-D for the same material layout ande/l,as shown in Figs.22 and 23.The SEAs of the BHHs with H—H and S—H exhibit significantly higher value than ones with H—S and S—S,but they have the lower CLE and higher σp,as shown in Figs.14 and 16.

    The energy absorption efficiency curves(that is,SEA-σout)of the BHHs with Δρ=0.115 under different impact velocities are shown in Fig.24.In the figure,σoutis the output stresses of the BHHs during the crushing,which can be given by the ratio of the contact force on fixed rigid plate to the original cross-section area of the specimen.The most effective absorbing energy structure is expected to be located in the upper-left corner of the chart.For lowvelocity impacting (v=5 m/s

    Fig.21.SEAs for the BHHs under different impact velocities.

    Fig.22.Effects of e/l on SEA for the BHHs with different topologies at v=20 m/s.

    Fig.23.Effects of multi-material configuration on SEA for the BHHs at v=20 m/s.

    Fig.24.Energy absorption efficiency for the BHHs with different topologies (⊿ρ=0.115).

    4.Conclusions

    In this paper,by imitating the arched shell structure of crab chela and replacing the structural joints of the RSHs with smaller arched topologies,three BHHs were designed to improve the crashworthiness and energy absorbing efficiency.The effects of hierarchical topologies and multi-material configuration on inplane dynamic behaviors of the BHHs were numerically explored.The main conclusions were drawn as follows:

    (1) For the BHHs,the crushing deformation modes mainly depends upon hierarchical topology and impact velocity.Introducing of a hierarchical arched structure into the joints of the RSHs will result in changing the macro-/micro-structure deformation characteristics.For low-velocity crushing(v=5 m/s

    (2) The BHHs have higher CLE than the RSHs,which is mainly related toe/l.Whene/l=0.417,the BHHs show the higher CLE,and lower σp,Δσ and ULC.The introduced hierarchical cell will deform first for the BHHs with H—S,which causes more stable collapse deformation and good crushing load uniformity.The CLE of the BHHs can be effectively controlled through the properly choice of hierarchical topolgies and multi-material layout.Moreover,empirical formulas of densification strains and plateau stresses are deduced to characterize the dynamic behaviors of the BHHs,which are in keep well with FE results.

    (3) When the cell-wall thickness is same,the SEA of the BHH-R is higher than that of other hierarchical topologies under the same impact velocity,which can be improved by reducinge/l.However,the SEAs of the BHH—C and BHH-D have little sensitivity one/l.The SEA of the BHHs can be controlled bye/land multi-material layout.The SEA of the BHH-R is higher than that of the BHH—C and BHH-D for the same material layout ande/l.It is interesting that the BHH-D will have the higher energy absorbing efficiency than the BHH—C and BHH-R for the same relative density.

    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.

    Acknowledgments

    The authors are grateful for the financial support provided by the Natural Science Foundation of Hebei Province of China [No.A2020502005],the Fundamental Research Funds for the Central Universities [No.2020MS113]and Science &Technology Program of Baoding [No.1911ZG019].

    日本-黄色视频高清免费观看| 国产国拍精品亚洲av在线观看| 中国国产av一级| 中文字幕精品免费在线观看视频 | 国产av精品麻豆| 只有这里有精品99| 亚洲精品国产色婷婷电影| 能在线免费看毛片的网站| 如何舔出高潮| 大香蕉久久网| 精品人妻视频免费看| 亚洲精品国产成人久久av| 一本一本综合久久| 日韩电影二区| 永久网站在线| 十分钟在线观看高清视频www | 夜夜骑夜夜射夜夜干| 亚洲精品乱码久久久v下载方式| 国产伦在线观看视频一区| 91午夜精品亚洲一区二区三区| 国产精品人妻久久久影院| 欧美xxxx黑人xx丫x性爽| 在线亚洲精品国产二区图片欧美 | 精品一品国产午夜福利视频| 22中文网久久字幕| 激情五月婷婷亚洲| 三级国产精品片| 久久久久网色| 亚洲丝袜综合中文字幕| 国产免费一级a男人的天堂| 最近最新中文字幕免费大全7| 女人十人毛片免费观看3o分钟| 美女高潮的动态| 18+在线观看网站| 22中文网久久字幕| 免费大片18禁| 中文字幕亚洲精品专区| 国产欧美亚洲国产| 欧美bdsm另类| 日日摸夜夜添夜夜爱| 婷婷色av中文字幕| 少妇高潮的动态图| 国产一区二区在线观看日韩| 精品久久久久久久末码| 久久久a久久爽久久v久久| 久久国产精品大桥未久av | 好男人视频免费观看在线| 精品99又大又爽又粗少妇毛片| 99国产精品免费福利视频| freevideosex欧美| 噜噜噜噜噜久久久久久91| 国产成人aa在线观看| 日韩电影二区| av播播在线观看一区| 久久久久人妻精品一区果冻| 成人二区视频| 亚洲av二区三区四区| 男女下面进入的视频免费午夜| 国产免费视频播放在线视频| 久久久久久久国产电影| 插逼视频在线观看| 亚洲内射少妇av| 国产精品久久久久久精品电影小说 | 日韩av在线免费看完整版不卡| 一边亲一边摸免费视频| 水蜜桃什么品种好| 91精品一卡2卡3卡4卡| 国产黄色免费在线视频| 丰满迷人的少妇在线观看| 国产乱来视频区| h视频一区二区三区| 日韩欧美 国产精品| 日韩av免费高清视频| 老司机影院成人| 中文字幕久久专区| 欧美亚洲 丝袜 人妻 在线| 国产视频首页在线观看| 男人狂女人下面高潮的视频| 亚洲成色77777| www.色视频.com| 能在线免费看毛片的网站| 丝袜喷水一区| 大话2 男鬼变身卡| 日本猛色少妇xxxxx猛交久久| 久久久久精品性色| 伦精品一区二区三区| 国产精品成人在线| 91精品国产国语对白视频| 国产欧美日韩一区二区三区在线 | 另类亚洲欧美激情| 国产午夜精品久久久久久一区二区三区| 99九九线精品视频在线观看视频| 国产91av在线免费观看| 免费观看av网站的网址| 97在线视频观看| 国产精品一区www在线观看| 精品人妻一区二区三区麻豆| kizo精华| 美女中出高潮动态图| 免费观看的影片在线观看| 伊人久久国产一区二区| 日韩国内少妇激情av| 男女啪啪激烈高潮av片| 日韩大片免费观看网站| 国产高清有码在线观看视频| 国产视频内射| 国产成人午夜福利电影在线观看| 国产av码专区亚洲av| 国产黄频视频在线观看| 99热网站在线观看| 日本av手机在线免费观看| 日韩 亚洲 欧美在线| 免费人成在线观看视频色| av在线蜜桃| 最黄视频免费看| 中文精品一卡2卡3卡4更新| 亚洲欧洲日产国产| 成人影院久久| 亚洲精品国产色婷婷电影| 一本一本综合久久| 内地一区二区视频在线| 免费看光身美女| 亚洲成人中文字幕在线播放| 成人18禁高潮啪啪吃奶动态图 | 啦啦啦啦在线视频资源| 久久精品国产鲁丝片午夜精品| 国产成人午夜福利电影在线观看| 波野结衣二区三区在线| 久久久亚洲精品成人影院| 少妇丰满av| 男女边吃奶边做爰视频| 一二三四中文在线观看免费高清| 毛片女人毛片| 国产亚洲一区二区精品| 99久久人妻综合| 人妻制服诱惑在线中文字幕| 国产乱人视频| 3wmmmm亚洲av在线观看| 深爱激情五月婷婷| 免费黄色在线免费观看| 国产精品久久久久久精品古装| 乱码一卡2卡4卡精品| 女人久久www免费人成看片| 午夜福利网站1000一区二区三区| 午夜日本视频在线| 少妇人妻久久综合中文| 亚洲欧洲日产国产| 日韩亚洲欧美综合| 99久久精品国产国产毛片| 97精品久久久久久久久久精品| 97在线人人人人妻| 高清视频免费观看一区二区| 九草在线视频观看| 久久热精品热| 搡女人真爽免费视频火全软件| av在线播放精品| 热99国产精品久久久久久7| av女优亚洲男人天堂| 99热这里只有精品一区| 亚洲国产毛片av蜜桃av| 人妻一区二区av| 亚洲国产高清在线一区二区三| 国产精品一区二区性色av| 国产精品福利在线免费观看| 国产亚洲91精品色在线| 精品久久久久久久末码| 欧美高清成人免费视频www| 日本午夜av视频| 日韩欧美精品免费久久| 久热这里只有精品99| 只有这里有精品99| 夜夜骑夜夜射夜夜干| 亚洲精华国产精华液的使用体验| 91午夜精品亚洲一区二区三区| 免费av中文字幕在线| 久久综合国产亚洲精品| 午夜福利影视在线免费观看| 成年人午夜在线观看视频| 午夜精品国产一区二区电影| 久热这里只有精品99| 插逼视频在线观看| 日产精品乱码卡一卡2卡三| 欧美bdsm另类| 男女边吃奶边做爰视频| 一本久久精品| 嘟嘟电影网在线观看| 成人一区二区视频在线观看| 中文资源天堂在线| 免费av不卡在线播放| 亚洲精品自拍成人| 国产成人91sexporn| 亚洲色图综合在线观看| 久久婷婷青草| 久久99热这里只有精品18| 国产精品久久久久久精品电影小说 | 亚洲欧洲国产日韩| 丝袜脚勾引网站| 日本黄色日本黄色录像| 亚洲国产欧美人成| 国产91av在线免费观看| 亚洲真实伦在线观看| 久久韩国三级中文字幕| 久久久精品免费免费高清| 婷婷色综合大香蕉| 免费观看的影片在线观看| av视频免费观看在线观看| 国产精品无大码| 午夜老司机福利剧场| 少妇的逼水好多| 又大又黄又爽视频免费| 日韩一区二区三区影片| 最近中文字幕2019免费版| 国产亚洲午夜精品一区二区久久| 国产乱人偷精品视频| 蜜臀久久99精品久久宅男| 精品久久久久久久末码| 日本黄色日本黄色录像| 激情 狠狠 欧美| 亚洲国产精品国产精品| 日韩av不卡免费在线播放| 久久国产亚洲av麻豆专区| 大又大粗又爽又黄少妇毛片口| 亚洲熟女精品中文字幕| 中文乱码字字幕精品一区二区三区| 人妻制服诱惑在线中文字幕| 久久影院123| videos熟女内射| 亚洲欧美日韩无卡精品| 午夜激情福利司机影院| 制服丝袜香蕉在线| 最近中文字幕高清免费大全6| 成人美女网站在线观看视频| 美女主播在线视频| 狂野欧美白嫩少妇大欣赏| 国产黄片美女视频| 啦啦啦中文免费视频观看日本| 国产色婷婷99| 99热网站在线观看| 久久久久久人妻| 久久久亚洲精品成人影院| 欧美xxxx黑人xx丫x性爽| 尾随美女入室| 夫妻午夜视频| 久久久久久久亚洲中文字幕| 亚洲综合色惰| 国产欧美日韩一区二区三区在线 | 免费观看在线日韩| 精品久久久久久久末码| www.av在线官网国产| 国产精品免费大片| 亚洲第一av免费看| 色吧在线观看| 精品亚洲乱码少妇综合久久| 日本午夜av视频| 亚洲欧美一区二区三区国产| 最近中文字幕2019免费版| 综合色丁香网| 欧美老熟妇乱子伦牲交| 国产一区亚洲一区在线观看| 国产乱来视频区| 亚洲无线观看免费| 九九爱精品视频在线观看| 国产精品不卡视频一区二区| 国产成人91sexporn| 国产人妻一区二区三区在| 国产精品人妻久久久久久| 婷婷色综合大香蕉| 春色校园在线视频观看| 成人午夜精彩视频在线观看| 26uuu在线亚洲综合色| 国产亚洲精品久久久com| 成人综合一区亚洲| 久久人妻熟女aⅴ| 亚州av有码| 久久国产亚洲av麻豆专区| 国产精品一二三区在线看| 亚洲三级黄色毛片| 欧美三级亚洲精品| 亚洲美女视频黄频| 超碰av人人做人人爽久久| a级毛色黄片| 中国美白少妇内射xxxbb| 男的添女的下面高潮视频| 日韩av免费高清视频| 久久久久国产精品人妻一区二区| 中文资源天堂在线| 麻豆国产97在线/欧美| 观看免费一级毛片| 精品久久久久久电影网| 又黄又爽又刺激的免费视频.| 午夜福利在线在线| 国内精品宾馆在线| 波野结衣二区三区在线| 欧美成人a在线观看| 尾随美女入室| 亚洲国产精品一区三区| 久久亚洲国产成人精品v| 日本猛色少妇xxxxx猛交久久| 国产亚洲av片在线观看秒播厂| 少妇丰满av| 欧美3d第一页| 只有这里有精品99| 国产视频内射| 日本与韩国留学比较| av女优亚洲男人天堂| 精品国产露脸久久av麻豆| 91精品国产国语对白视频| 99国产精品免费福利视频| 成人国产av品久久久| 99久国产av精品国产电影| 国产一级毛片在线| 国产在线男女| 欧美亚洲 丝袜 人妻 在线| 午夜老司机福利剧场| 青春草视频在线免费观看| 亚洲欧美日韩无卡精品| 男人和女人高潮做爰伦理| 99久久精品一区二区三区| 国产精品一区二区在线观看99| 国产精品av视频在线免费观看| 99久久精品国产国产毛片| 人妻 亚洲 视频| 各种免费的搞黄视频| 亚洲美女搞黄在线观看| 少妇被粗大猛烈的视频| 建设人人有责人人尽责人人享有的 | 搡女人真爽免费视频火全软件| 美女视频免费永久观看网站| 一区在线观看完整版| 成年免费大片在线观看| 只有这里有精品99| 超碰97精品在线观看| 婷婷色麻豆天堂久久| 欧美精品亚洲一区二区| 91午夜精品亚洲一区二区三区| 久久久久网色| 国内少妇人妻偷人精品xxx网站| 亚洲精品日本国产第一区| 精品国产露脸久久av麻豆| 欧美成人a在线观看| 国产成人精品福利久久| 国产精品久久久久久av不卡| 网址你懂的国产日韩在线| 嫩草影院新地址| 亚洲精华国产精华液的使用体验| 成人18禁高潮啪啪吃奶动态图 | 国产视频内射| 国产一区二区三区av在线| 少妇熟女欧美另类| 国产精品国产三级专区第一集| 国产精品久久久久久精品电影小说 | 一二三四中文在线观看免费高清| 伦理电影大哥的女人| 高清不卡的av网站| 伦理电影大哥的女人| 国产成人免费无遮挡视频| 国产精品一区二区三区四区免费观看| 啦啦啦视频在线资源免费观看| 亚洲久久久国产精品| 成人18禁高潮啪啪吃奶动态图 | 啦啦啦啦在线视频资源| 一区在线观看完整版| 男人狂女人下面高潮的视频| 少妇被粗大猛烈的视频| av线在线观看网站| 久久99热这里只有精品18| 亚洲图色成人| 中国三级夫妇交换| 99热这里只有是精品50| 中文精品一卡2卡3卡4更新| 黄色欧美视频在线观看| 国产精品欧美亚洲77777| 成人午夜精彩视频在线观看| 久久影院123| 中国美白少妇内射xxxbb| 成人漫画全彩无遮挡| 亚洲无线观看免费| 久久97久久精品| 精品酒店卫生间| 久久韩国三级中文字幕| 插逼视频在线观看| 中文乱码字字幕精品一区二区三区| 大话2 男鬼变身卡| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 51国产日韩欧美| 亚洲精品乱码久久久久久按摩| 精品人妻熟女av久视频| 久久99精品国语久久久| 特大巨黑吊av在线直播| 精品一区二区免费观看| 亚洲国产毛片av蜜桃av| 国产亚洲午夜精品一区二区久久| 在线播放无遮挡| 美女xxoo啪啪120秒动态图| 亚洲性久久影院| 日本爱情动作片www.在线观看| 黄色配什么色好看| 18禁在线无遮挡免费观看视频| 插逼视频在线观看| av女优亚洲男人天堂| 日韩av不卡免费在线播放| 99热全是精品| 欧美日韩精品成人综合77777| 成年av动漫网址| 少妇人妻久久综合中文| 黄色欧美视频在线观看| 亚洲精品国产成人久久av| 中国美白少妇内射xxxbb| 午夜免费鲁丝| av福利片在线观看| 五月天丁香电影| 日韩一本色道免费dvd| 国产精品国产三级国产专区5o| 日本av免费视频播放| 蜜桃亚洲精品一区二区三区| 日韩不卡一区二区三区视频在线| 老师上课跳d突然被开到最大视频| 十分钟在线观看高清视频www | 麻豆成人av视频| 黄色怎么调成土黄色| 女人久久www免费人成看片| 欧美日韩国产mv在线观看视频 | 欧美丝袜亚洲另类| 日韩伦理黄色片| 丝瓜视频免费看黄片| 99热网站在线观看| 免费人成在线观看视频色| 国产视频首页在线观看| 九草在线视频观看| 最近手机中文字幕大全| 在线观看av片永久免费下载| 国产精品久久久久久精品古装| 久久久久精品久久久久真实原创| 又黄又爽又刺激的免费视频.| 日本猛色少妇xxxxx猛交久久| 亚洲精华国产精华液的使用体验| 女人十人毛片免费观看3o分钟| 久久久久久久久久人人人人人人| av不卡在线播放| 男女边摸边吃奶| 成年免费大片在线观看| 国产精品av视频在线免费观看| 又黄又爽又刺激的免费视频.| 国产精品免费大片| 久久久久久久久久成人| 97热精品久久久久久| 国产精品一区二区性色av| 亚洲精品第二区| 亚洲成人中文字幕在线播放| 新久久久久国产一级毛片| 大片电影免费在线观看免费| 熟女av电影| 亚洲aⅴ乱码一区二区在线播放| 99视频精品全部免费 在线| 亚洲欧美中文字幕日韩二区| 国产一区二区三区av在线| 五月开心婷婷网| 啦啦啦视频在线资源免费观看| 美女脱内裤让男人舔精品视频| 色婷婷久久久亚洲欧美| 国产精品福利在线免费观看| 国产深夜福利视频在线观看| 少妇的逼好多水| 尾随美女入室| 国产精品久久久久久av不卡| 午夜福利网站1000一区二区三区| 妹子高潮喷水视频| 日韩制服骚丝袜av| 久热这里只有精品99| 日本av手机在线免费观看| 大陆偷拍与自拍| 色婷婷久久久亚洲欧美| 国产午夜精品久久久久久一区二区三区| 欧美高清性xxxxhd video| 亚洲一级一片aⅴ在线观看| 成人一区二区视频在线观看| 黄色怎么调成土黄色| 精品国产露脸久久av麻豆| 亚洲欧美一区二区三区国产| 久久精品国产自在天天线| tube8黄色片| 成人影院久久| av福利片在线观看| 亚洲欧美一区二区三区国产| 99热国产这里只有精品6| 免费黄频网站在线观看国产| 伦理电影免费视频| 人人妻人人澡人人爽人人夜夜| 中文乱码字字幕精品一区二区三区| 国产黄频视频在线观看| 国产日韩欧美亚洲二区| 日韩强制内射视频| 国产精品.久久久| 中文乱码字字幕精品一区二区三区| 午夜福利在线在线| 在线亚洲精品国产二区图片欧美 | 欧美精品一区二区免费开放| 欧美丝袜亚洲另类| 蜜臀久久99精品久久宅男| 亚洲国产精品成人久久小说| 日韩欧美 国产精品| 99久国产av精品国产电影| 欧美最新免费一区二区三区| 亚洲综合精品二区| 三级国产精品欧美在线观看| 春色校园在线视频观看| 亚洲国产高清在线一区二区三| 日韩亚洲欧美综合| 日本黄大片高清| av卡一久久| 大香蕉久久网| av国产精品久久久久影院| 日本免费在线观看一区| 99热全是精品| 国产成人免费观看mmmm| 亚洲精华国产精华液的使用体验| videossex国产| 亚洲激情五月婷婷啪啪| 亚洲图色成人| 亚洲精品国产色婷婷电影| 国产精品一区二区三区四区免费观看| 麻豆成人午夜福利视频| 欧美bdsm另类| 午夜免费观看性视频| 高清不卡的av网站| 国产爱豆传媒在线观看| 纯流量卡能插随身wifi吗| 大码成人一级视频| 99re6热这里在线精品视频| 亚洲va在线va天堂va国产| 欧美亚洲 丝袜 人妻 在线| 尤物成人国产欧美一区二区三区| 国产精品一区二区在线不卡| 欧美最新免费一区二区三区| 久久久久性生活片| 丰满人妻一区二区三区视频av| 91狼人影院| 免费观看的影片在线观看| 欧美日韩在线观看h| 免费黄频网站在线观看国产| 在线观看美女被高潮喷水网站| 精品熟女少妇av免费看| 啦啦啦啦在线视频资源| 亚洲自偷自拍三级| 亚洲国产日韩一区二区| 简卡轻食公司| 草草在线视频免费看| 久久久欧美国产精品| 男男h啪啪无遮挡| 国产爱豆传媒在线观看| 久久精品国产亚洲av天美| 国产精品欧美亚洲77777| 黄色日韩在线| 国产 精品1| 免费播放大片免费观看视频在线观看| 午夜视频国产福利| 日韩欧美精品免费久久| 街头女战士在线观看网站| 在线观看一区二区三区激情| 午夜免费鲁丝| 交换朋友夫妻互换小说| 欧美高清性xxxxhd video| 亚洲精品国产成人久久av| 欧美区成人在线视频| 亚洲成人一二三区av| 国产欧美另类精品又又久久亚洲欧美| 免费看光身美女| 丝袜脚勾引网站| 九九久久精品国产亚洲av麻豆| 久久毛片免费看一区二区三区| 午夜老司机福利剧场| 国产精品无大码| 日韩电影二区| 18禁在线无遮挡免费观看视频| 久久99精品国语久久久| 成人美女网站在线观看视频| 啦啦啦视频在线资源免费观看| 久久精品人妻少妇| 老女人水多毛片| 欧美日韩视频高清一区二区三区二| 久久人妻熟女aⅴ| 国产欧美日韩一区二区三区在线 | 最黄视频免费看| 赤兔流量卡办理| 午夜老司机福利剧场| 国产 一区精品| 久久这里有精品视频免费| 欧美成人一区二区免费高清观看| 亚洲欧美精品专区久久| 日韩电影二区| 国模一区二区三区四区视频| 国产一区二区三区av在线| 成人一区二区视频在线观看| 欧美xxxx黑人xx丫x性爽| 国产精品无大码| 久久这里有精品视频免费| 午夜老司机福利剧场| 一区在线观看完整版| h视频一区二区三区| 一级毛片我不卡| 深夜a级毛片| 日韩欧美一区视频在线观看 | 中文在线观看免费www的网站| 日本爱情动作片www.在线观看| 日韩人妻高清精品专区| 极品少妇高潮喷水抽搐| 久久青草综合色| 22中文网久久字幕| 国产 一区 欧美 日韩| 最近中文字幕高清免费大全6| 日韩av在线免费看完整版不卡| 国产一区有黄有色的免费视频| 在线观看免费高清a一片| 国产精品福利在线免费观看|