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

    Inf l uence of conf i ning prestress on the transition from interface defeat to penetration in ceramic targets

    2016-04-18 09:13:00PatrikLUNDBERGRenRENSTROlofANDERSSON
    Defence Technology 2016年3期

    Patrik LUNDBERG*,René RENSTR?M,Olof ANDERSSON

    FOI,Swedish Defence Research Agency,SE-164 90 Stockholm,Sweden

    Inf l uence of conf i ning prestress on the transition from interface defeat to penetration in ceramic targets

    Patrik LUNDBERG*,René RENSTR?M,Olof ANDERSSON

    FOI,Swedish Defence Research Agency,SE-164 90 Stockholm,Sweden

    Replica scaled impact experiments with unconf i ned ceramic targets have shown that the transition velocity,i.e.,the impact velocity at which interface defeat ceases and ceramic penetration occurs,decreased as the length scale increased.A possible explanation of how this scale effect is related to the formation of a cone crack in the ceramic has been presented by the authors in an earlier paper.Here,the inf l uence of conf i nement and prestress on cone cracking and transition velocity is investigated.The hypothesis is that prestress will suppress the formation and growth of the cone crack by lowering the driving stress.A set of impact experiments has been performed in which the transition velocity for four different levels of prestress has been determined.The transition velocities as a function of the level of conf i ning prestress is compared to an analytical model for the inf l uence of prestress on the formation and extension of the cone crack in the ceramic material.Both experiments and model indicate that prestress has a strong inf l uence on the transition from interface defeat to penetration,although the model underestimates the inf l uence of prestress.

    Impact;Ceramic;Armour;Interface defeat;Dwell;Conf i nement;Prestress

    1.Introduction

    The high strength of armour ceramics [1-3]makes it possible to partially or totally defeat high velocity projectiles directly at the surface of the ceramic material.This phenomenon is called interface defeat or dwell[4-17]and is an important defeat mechanism in,e.g.,light armour applications.

    One limitation when applying this in heavier armour designs is that it appears to be length scale dependent.Replica scaled impact experiments with unconf i ned ceramic targets show that the transition velocity,i.e.,the velocity at which interface defeat ceased and ceramic penetration occurred,decreased as the length scale increased [11].A probable explanation of the observed scale effect is that although maximum shear strength determines the upper bound for the transition from interface defeat to penetration,it is usually limited by the formation and growth of macroscopic cracks.Since the crack resistance of ceramic materials decreases with increasing length scale,in contrast to the otherwise scale-invariant stress f i eld,the extension of a crack to a critical size will occur at a lower impactvelocity in a larger target.An analytical model in [11]for the inf l uence of length scale on the growth of a cone shaped modus I crack in thick unconf i ned ceramic targets gave reasonable results compared to the replica scaled impact experiments. The model showed that the projectile pressure at transition,i.e.,the impact velocity at which the contact pressure exceeds the strength of the ceramic material and penetration initiates,is proportional to one over the square root of the length scale.

    A possible way to suppress the formation and growth of macroscopiccracksistoprestresstheceramicmaterial.Theinf l uence of prestress and the related failure modes of impacted ceramics have been studied by several authors.The papers [18-20]report experimental data on small calibre projectiles impacting thin prestressed ceramics (i.e.,the thickness is of the same order as the diameter of the projectile).These studies show that prestress reducesdamageintheformoffewermacroscopiccracksandthat the trajectory of possible cone cracks becomes shallower.An increase in protective performance was also observed.The papers [8,9]report experimental data on model scale long rod projectiles impacting thick ceramic targets (i.e.,the thickness is much larger than the diameter of the projectile).The experiments in [8]with large and heavily conf i ned and prestressed targets showed similar interface defeat velocities as small,unconf i ned targetsin [9].Thisindicatestheneedofprestressinlargertargets.Holmquist and Johnson [21]and later Runqiang et al.[22]conducted a computational study on the responses of a small scale thick prestressed ceramic target tested by Lundberg et al. [7]. Various levels of prestress and stress states were simulated.Their studies showed that prestress enhanced the performance and that the velocity at which ceramic penetration occurred,i.e.,the transition velocity,could be increased by prestress.

    This paper explores the inf l uence of a radial conf i ning prestress on the transition from interface defeat to penetration for a thick ceramic target.Although the physical background of the inf l uence of prestress on the transition velocity in ceramic targets is not fully explained,impact experiments as well as modelling indicate that it is intimately linked to ceramic fracture.A hypothesis proposed in [11]is that the centre part of the ceramic suddenly loses radial support as a result of the cone cracking.A conf i ning prestress will suppress the growth of the cone crack by lowering the stress intensity over the crack tip.In order to overcome this virtual toughening of the ceramic,the projectile pressure on the surface of the target must be increased relative to that for an unconf i ned target in order to initiate critical fracture.A set of impact experiments have been performed in which the transition velocities for four different levels of prestress were determined.Two grades of silicon carbide ceramics with slightly different mechanical properties were used.The experimental technique used is presented in the paper together with the determined transition velocities versus radial conf i ning prestress.The experimental data are compared to an extended version of the model presented in [11].

    2.Model of cone crack under conf i ning prestress

    The inf l uence of a radial conf i ning prestress on the formation and extension of a cone shaped crack to a critical size is approximated in the model by the assumption that the crack extension occurs along a surface of principal stress.The normal stressonthissurfaceiscalculatedforthecaseofanaxi-symmetric contact pressure from a projectile in a state of interface defeat and for a radial conf i ning prestress,respectively.The critical normal stress for propagating the crack is determined through:(i)a stress intensity factor at the tip of the crack and (ii)a function of the inf l uence of external load and geometry on the path of the crack.

    The detailed description of the present model is divided into four sections:2.1 Projectile contact pressure and stresses in the target,2.2 Principal surfaces and stresses,2.3 Crack initiation and propagation under conf i ning prestress and 2.4 Crack opening under conf i ning prestress.

    2.1.Projectile contact pressure and stresses in the target

    A long cylindrical projectile is assumed to fl ow axisymmetrically on the fl at and friction-free surface of an otherwise unbounded,but radially prestressed target,see Fig.1.The fl ow and the loading on the surface are steady,i.e.,the initial transient part of the impact process is not considered.The target material is linearly elastic with Young's modulus E and Poisson's ratio ν,respectively.The material of the projectile is linearly elastic and perfectly plastic with bulk modulus Kp,yield strengthσypand densityρp.

    Fig.1.Projectile shape (solid line),axi-symmetric pressure of projectile (dotted line)and crack trajectory in target(gray)during a state of interface defeat.

    With the assumption that the effects of yield strength and compressibility are small relative to that of inertia,the axisymmetric contact pressure of projectile load can be approximated [7]by

    whereq( r)is the radial pressure distribution corresponding to inertia,

    and

    Herev0is the impact velocity of the projectile andqpis the stagnation pressure of an ideal fl uid with densityρp.The dimensionless parameters α >> 1 and β <<1 relate elastic and plastic effects to the effect of inertia.The in fl uence of β in Eq.(1)is evaluated from simulations in [17]and the radial distribution ofq( r)is taken from a low-velocity water jet [7].

    The stress fi eld σij(r, z)of projectile pressure p in the semiin fi nite elastic target half-space,is expressed by a Boussinesq's potential as

    whereIijis an in fl uence-function (negative in compression)for a point load of the amount of p? d ? d υat a surface point (υ, ?),with p according to Eq.(1).Indicesi, jare the generic spatial variables r,φ,and z and A is the radius?of the circular limit of the distribution of projectile pressure.The radial stress componentσrris affected by the con fi ning prestresspcfso that the stress components from projectile and con fi nement together are expressed as

    2.2.Principal surfaces and stresses

    The assumed crack path along the principal surface z1( r)is determined as

    whereθ(r)is the angle betweenz1( r)and the target surface z=0.The target surface is a principal surface since friction free,thus θ= π2forz=0.Rewriting Eq.(6),by the use of tan 2θ= 2 tan θ(1-tan2θ),Eq. (5) and the well-known expressions for stress-geometry on principal surfaces gives

    where the sign is chosen so that the resulting surface crosses the target surface perpendicularly.The crack arc lengthc( r)along the surfacez1( r ),the principal stressesσ1,σ2andσ3(normal and tangential to this surface)and the strainε1,are expressed as

    and

    The principal surfacez1( r )is chosen in such a way that the stressσ1is tensile on z1( r ).

    2.3.Crack initiation and propagation under con fi ning prestress

    The maximum tensile stress at radiusr =rion the surface of the target initiates fracture for

    and propagates the crack alongz1( r )as long asσ1exceeds a critical stressσc.The critical stress will be determined through a stress intensity factor for small crack extension c,as

    where h is a function of external load and geometry.Although not explicitly mentioned in [11],h was found to correspond to unit value for pcf=0.Since the geometry and loading conditions here are similar to the one used in [11],except for the radial con fi ning prestress 0 < pcf,h is sought as a function ofpcfand p0in such a way that the criterion of crack propagation and the critical stressσccan be written in comparison to that for pcf=0as

    whereK1cis the fracture toughness of the ceramic material. Eqs.(7)-(12)result in an equation for critical crack extensions rcas roots of the equation

    A more explicit expression of parameter in fl uence in Eq.(13)is achieved by the use of dimensionless variables.The principal stressσ1,con fi ning prestresspcf,crack arc lengthc and spatial coordinatesr,z are expressed in units of projectile pressure p0and radiusa through,pcf=p0c =,Eq.(8)andr =,z =,respectively.The overbar sign designates dimensionless.

    Introducing dimensionless variables into Eq. (13)gives a relation between the projectile radiusaand a critical projectile pressure.This critical projectile pressure corresponds to the development of a cone crack with radial extensionrcand is named the transition pressure where the superscript*is used to denote critical/threshold quantities.The transition pressure corresponds to the transition velocityvia Eqs. (1)-(3)and is here assumed to be the lowest bound for transition from interface defeat to penetration for different levels of con fi ning prestresspcfaccording to

    where the right hand member of Eq. (14)is valid for roots rc=of Eq. (13).

    The transition pressure for a con fi ned target relative to that of zero con fi nement is approximately expressed as

    2.4.Crack opening under conf i ning prestress

    Consider a sub-critical crack and its expansion to a stable,open crack and that the projectile force act to propagate thecrack and the force of con fi ning prestress is counteracting the crack propagation.The function of external load and geometry h in Eq. (11)is approximated through the opening process of the crack,that is,the changeΔU of the accessible partUof the internal elastic energy of the ceramic material.The energy term ΔUis the elastic energy released in crack-opening,so that the derivative ofΔUwith respect to crack propagation will be the relevant energy to drive the crack.Therefore,ΔUis proportional to the square of the stress intensity and thereby the square of h,Eq.(11).

    With the force of con fi ning prestress in place and with a projectile load close to maximum,but before cracking,the accessible internal energy may be expressed as

    where F,P and ur( F ),uz(F ),ur( P)are the forces from projectile,con fi nement pressure and corresponding boundary displacements of the ceramic cylinder,respectively.The internal energy is expressed as if the ceramic body is fi rst squeezed at its cylindrical surface by a constant pressurepcfand then at its cross-end surfaces,each of area A0,by constant pressureFA0.The displacements in Eq. (16)may then be expressed as

    whereRandL are radius and length of the ceramic cylinder and A0=πR2its cross-end surface.The differenceΔUcan be expressed as

    whereU′and u′denote accessible internal energy and displacement as if the open stable crack is already present when F is applied.As a fi rst order approximation of Eq. (18),it is noted that(i)the volume of the open crack is small and(ii)at the moment of its instability,the critical geometry and volume at the crack are unaffected by P.It is here assumed thatΔuican be expressed as

    where ηiare assumed to be small constants.Eq.(18)can now be expressed as

    where,according to inferred assumptions and Eq. [11],

    so that

    Eq.(22)gives the effect of con fi ning prestress for a constant P and an approximate expression ofh.

    Using Eq.(17),the projectile radiusa,the cross-end area of the ceramicA0=πR2and the expressions of force,viz.,F(xiàn)=p02πa2and P = pcfπ2RL ,the terms in the third of Eqs.(22)are expressed as

    Finally,Eqs. (14)and (15)are substituted byEq. (23)to express the parameters of the transition pressurefor a confi ning prestress 0 ≤relative to that of zero con fi nement.

    3.Experiments

    The impact experiments were performed using a reverse impact technique.The stationary projectiles were suspended in blocks of Divinycell material (density 45 kg/m3)and mounted in front of the muzzle of the gun,see Fig.2.Two different qualities of silicon carbide materials have been used:SiC-B and SiC-X1,both materials are from CoorsTek (former BAE Systems Advanced Ceramics Division and Cercom Inc,Vista,CA).The SIC-B material was initially delivered as large cylinders with diameter 50 mm and length 50 mm.From these,smaller cylinders with diameter 20 mm and length 20 mm were produced.The SiC-X1 material was delivered as cylinders with diameter 20 mm and length 20 mm.

    Data on microstructure,Young's modulus and fracture toughness of the SiC-B used has been published byWereszczak et al. [23]and properties for SiC-X1 has been provided by CoorsTek.The properties are given in Table 1.

    Fig.2.Experimental set-up showing the stationary projectile suspended in a block of Divinycell material and mounted in front of the muzzle of the gun. Images plates and normal intensif i er screens and f i lm for the f l ash X-ray pictures are placed in the U-shaped cassette below the projectile.

    Table 1Material data for SiC-B and SiC-X1.

    Table 2Hardness and fracture toughness for SiC-B and SiC-X1.

    The Vickers hardnessH and the fracture toughnessKIcwere estimated by means of a Wholpert macro hardness indenter equipped with video-system and imaging software. SiC-B and SIC-X1 samples were carefully polished using a semi-automatic polishing machine,and each material was indented 9 times.The fracture toughness was determined according to Anstis et al.[24].The estimated Vickers hardness and fracture toughness of SiC-B and SiC-X1,normalised to these quantities for SiC-B,are given in Table 2.

    The projectiles were fl at-ended cylinders made of a sintered tungsten alloy,Y925 from Kennametal HertelAG [25].Material data for the projectile are given in Table 3.

    The target con fi nement consisted of a steel tube (B?hler W725, Poisson'sratio νs=0. 3andYoung'smodulus Es=186GPa ).The external diameter of the steel tube was originally 30 mm and its internal diameter was 0.07 mm smaller than the diameter of the ceramic cylinder.Shrink fi t was achieved by heating the steel tube to about 500 °C before inserting the ceramic cylinder.After cooling,the con fi nement was machined to fi nal shape with tube wall thicknessest=1,2 and 4 mm,respectively.The shrink fi t resulted in a con fi ning prestress pcfon the ceramic cylinder estimated by

    whereδis the difference between the radius of the ceramic cylinder and the interior radius of the con fi ning steel tube.Both the ceramic cylinders and the con fi ning steel tubes were ground to fi nal dimensions.The variation in ceramic and steel tube diameters gave a maximal variation inδof the order of 0.005 mm.

    The front surface of the ceramic target was protected with a circular copper cover.The copper cover is expanded in its central part to a cylinder.It was glued onto the front surfacealong the rim.The geometries of the four target types used are shown in Fig.3.The different target types are labelled (a)-(d),corresponding to uncon fi ned,1,2 and 4 mm thick con fi nement,respectively.The different test series are summarised in Table 4.

    Table 3Material data for the projectile material Y925.

    The impact velocity and target response were evaluated from fl ash X-ray pictures.A 150 kV X-ray system was used together with image plates for the velocity pictures taken before the interaction.For the penetration pictures,a 450 kV X-ray system was used together with both the image plates and normal intensi fi er screens and X-ray fi lm.The 450 kV fl ashes were positioned at the same distance from the barrel and radially separated by 30°.The image plate picture was digitised using a laser scanner.The X-ray fi lms were digitised using a fl at-bed scanner.Enhancement was achieved by image-processing (contrasts,edges etc.)before evaluation.The uncertainty in the impact velocity evaluation,due to limited X-ray image resolution,measurement errors etc.,was within ±5 m/s.

    The transition velocity was estimated by systematically varying the impact velocityv0.Ideally,the transition occurs at a well-de fi ned impact velocity,but in practice it is only possible to determine a velocity interval in which the transition occurs. The lower and upper limits of this interval correspond to the highest impact velocity observed without penetration and the lowest impact velocity with penetration.The transition velocity was estimated as the centre-point of this interval.The velocity interval was indicated by adding ± half the length of the interval to the estimated transition velocity.Typically 4-8 impact tests were needed in order to determine the transition velocitywithin ± 50 m/s,viz.,a total number of 44 impact experiments has been performed.

    Fig.3.Target types: (a)unconf i ned,(b)conf i ned t=1 mm,(c)conf i ned t=2 mm,(d)conf i ned t=4 mm.Dimensions are in mm.

    Table 4Targets and projectiles used in the impact experiments.

    4.Results

    Examples of X-ray pictures from the impact tests are shown in Figs.4 and 5.Fig.4 shows an uncon fi ned SiC-X1 target under a state of interface defeat while Fig.5 illustrate the phenomena of interface defeat for four different levels of prestress:pcf=0,56,101 and 16?8 MPa.

    Fig.4.X-ray pictures of Y925 projectile and uncon fi ned SiC-X1 target at four different times after impact.The time interval between the pictures is 10 μs.The? i mpact velocity v0 = 967m sis slightly below the transition velocityv0 = 982 ±15ms.

    Fig.5.X-ray pictures of Y925 projectiles and SiC-X1 targets with four levels of prestress under a state of interface defeat: (a)t=0 mm,v0 = 967m s,(b)t=1mm,v0 = 1347m s,(c)t=2 mm,v0 = 1508m s and (d)t=4 mm,v0 = 1485m s.

    The transition velocitiesv0determined from the impact experiments and corresponding projectile pressures at transitioncalculated from Eqs. (1)-(3)are shown in Table 5 together with estimated levels of con fi ning prestresspcfaccording to Eq. (24).

    The principal surfacez1a used as crack-path in the model and the principal stressσ1p0along the crack are shown in Fig.6 versus radial crack extensionra for different levels of con fi ning prestresspcfp0,Eqs.(4),(5)and (7)-(9).The critical crack extensionrca ,determined as a root of Eq.(14),is illustrated in Fig.7(a)for an uncon fi ned target and Fig.7(b)show the function1h versus the con fi ning prestress pcfp0.

    The in fl uence of the fracture toughness KIcof the ceramic material on the transition pressure is illustrated in Fig.8. Experimentally determined as well as estimated transition velocities and pressures versus prestress are shown in Fig.9. Unless otherwise stated,the material data for SiC-B in Table 1 have been used in Figs.5-8.

    5.Discussion

    Fig.4 shows the development of interface defeat at an impact velocity just below the transition velocity in an unconfi ned SiC-X1 target.The copper cover in front of the ceramic reduces the initial effect of impact by establishing erosion of the projectile before the latter reaches the ceramic surface.Other studies [10,26]have shown that after radial fl ow had been established,it continued steadily for a long time,i.e.,the projectile load can be seen as quasi-static.Although the transition velocity increases with prestress,the velocity interval from interface defeat to penetration remains narrow and is not affected by the prestress and the fl ow onto of the ceramic surface look similar for the different targets,see Fig.5.

    The effect of prestress is clearly seen in the experiments;the transition from interface defeat to penetration is moved to signi fi cantly higher levels as the con fi ning pressure increases.The velocity at transition in uncon fi ned SiC-X1 target was found to be 982 m/s whereas a con fi ning pressure of pcf≈56MPa increased this transition velocity to 1367 m/s.This relatively low level of prestress almost doubles the projectile pressure at transition.The two leftmost open circles in Fig.9 correspond to these values.Further increase of the con fi ning pressure did not show the same strong in fl uence and a transition velocity of≈ 1500 ms was found unaffected in spite of increased confi ning pressure from 100 to 168 MPa.The SiC-B seams to behave in a similar way as the SiC-X1,though no data is available in-between uncon fi ned and 168 MPa of prestress.Although the number of experiments is not suff i cient for a statistical analysis,the SiC-B seems to perform slightly better than SiC-X1.This could be a result of the slightly higher fracture toughness for SiC-B or due to natural variations in the properties of the projectile which determines the transition pressure.Use of Eqs. (1)-(3)and typical uncertainties for the material data in Table 3 give an indication of an maximum error in transition pressure within ±1.7 GPa,i.e.,the observed differences between SiC-B and SiC-X1 may simply be within statistical f l uctuations.

    Table 5Transition velocities and corresponding estimates of the projectile pressure.The superscript*is used to denote critical/threshold quantities.

    The SiC-X1 experiments indicate a shift in behaviour for a conf i ning prestress around 100 MPa.This is illustrated with two grey sectors in Fig.9(a)and (b).The shift may indicate that another fracture mode,e.g.,modus-II cracking,has been activated.A higher transition velocity for a similar combination of ceramic and projectile materials has been reported [8]but thenwith strong axial and radial conf i nement.Such a high overall conf i ning pressure could affect both modus-I and II fracture.

    Fig.6. (a)Crack path z a and (b)principal tensile stress along crack path σ1 p0versus radial crack extension r a for different levels of con fi ning prestress.

    Fig.7. (a)Principal tensile stressσ1 p0(solid curve)and critical stress σc/p0(dashed curve)along the crack (r/a)for an estimated stress state at the point of transition in an uncon fi ned targetpcf p0 =0.(b)The function1 h (solid curve)versus con fi ning prestresspcf p0.

    Fig.8.Estimate of maximum projectile pressure at transition p?0versus confi ning prestress pcffor three different levels of ceramic fracture toughness: KIc=3(dashed curve),4 (solid curve)and 5 MPam?(dotted curve),respectively.

    The model for the projectile pressure in Eqs. (1)-(3)gives similar results as numerical simulations in [17]where more detailed material models for both the tungsten and the SiC material were used.The projectile pressure distribution used here is not identical to the one in [17]but gives only minor changes in the overall stress distribution and does not appreciably change the results.

    The inf l uence of the conf i ning prestress in the model consists of two parts:(i)the change in tensile stress f i eld due to the prestress and (ii)the inf l uence of prestress on the opening process of the crack itself.The model for the tensile stress over the crack assumes that the stress f i eld will be unaffected by the crack.This is a simplif i cation but does not signif i cantly change the general behaviour of the model.Fig.6 shows that the principal tensile stress over the crack decreases and changes direction towards the impact surface as a result of the con fi ning prestress.This will gradually reduce the effect of a radial prestress on the crack propagation.

    A valid solution rca of Eq. (13)is assumed to exist for aprojectile pressure suf fi cient to satisfy Eq. (12)and so that is the largest of the three-root solution.The critical crack length is therefore a constant for each grade of prestress since

    in Eq. (13)is a fi x function.In agreement with this assumption,F(xiàn)ig.7(a)shows the estimated critical radial crack extensionrca =6. 4for an uncon fi ned target.

    The in fl uence of prestress on the opening of the crack is approximated by an energy model,that is,the changeΔU of the accessible partU of the internal elastic energy of the ceramic material.The effect is described by the function of external load and crack geometry h,see Fig.7(b).

    Eq.(14)and Fig.8 show that the model in fl uence of fracture toughness on the transition from interface defeat to penetration is strong;the transition pressure at transition is directly proportional to the fracture toughness of the ceramic material.This is in line with the experimental fi ndings in [10]where the transition velocity was determined for four different silicon carbide materials with slightly different fracture toughness.

    The model shows an approximately linear relation between con fi ning prestress and the transition pressure and that this corresponds to the experimental data for prestress 0 ≤pcf≤100MPa .But the model underestimates the in fl uence of prestress with a factor of two relative to the experiments in this interval.The main hypothesis for the deviation between model and experiments is that the real con fi ning prestress is higher than estimated by Eq. (24).The reason for this is that dynamic effects,e.g.,the inertia of the con fi nement,have not been taken into account,neither in the model nor in Eq.(24).A dynamic effect would probably shift the experimental data towards higher prestress values and at the same time increase the gradient of the function h.A more detailed analysis of a possible dynamic effect is not within the scope of this paper and will probably require some type of continuum mechanic code.

    Fig.9. (a)Estimated maximum projectile pressure at transitionand (b)transition velocityversus con fi ning prestresspcf.Experiments with SiC-B are marked with fi lled circles and SiC-X1 with open circles,respectively.Error bars indicate maximum and minimum values for projectile velocity,projectile pressure and con fi ning prestress,respectively.Solid curve:model estimate,grey sectors:graphical illustration of a possible two mode behaviour,respectively.

    6.Conclusions

    The main conclusions of this study can be summarised as follows: (i)An analytical model for the relation between projectile pressure and propagation of a cone crack under a state of interface defeat in a ceramic target has been formulated.The model connects the effects of fracture toughness and con fi ning prestress of the ceramic to the transition from interface defeat to penetration. (ii)The model shows a strong in fl uence of radial prestress;the projectile pressure at transition increases linearly with the level of prestress. (iii)Impact experiments with four different levels of prestress show that prestress has a stronger in fl uence than predicted by the model;the maximum possible contact pressure at transition is more than doubled if prestress is increased from zero to around 100 MPa.(iv)Further increase of the con fi ning pressure did not show the same in fl uence and a transition velocity of ≈ 1500 ms was found unaffected in spite of increased conf i ning pressure. (v)The model shows a similar behaviour compared to the lower interval of prestress,but underestimates the inf l uence of prestress with a factor of two.(vi)The main hypothesis for the deviation between model and experiments is due to dynamic effects which will increase the conf i ning prestress relative to the one estimated by the quasistatic model.

    Acknowledgements

    This research was funded by the SwedishArmed Forces and by the Army Research Laboratory through US Naval Regional Contracting Centre,Contract No.W911NF0810271.

    [1]Vogler TJ,Chhabildas LC.Strength behaviour of materials at high pressure.Int J Impact Eng 2006;33:812-25.

    [2]Kanel GI,Zaretsky EB,Rajendran AM,Razorenov SV,Savinykh AS,Paris V.Search for conditions of compressive fracture of hard brittle ceramics at impact loading.Int J Plasticity 2009;25:649-70.

    [3]Paris V,F(xiàn)rage N,Dariel MP,Zaretsky E.Divergent impact study of the compressive failure threshold in SiC and B4C.Int J Impact Eng 2011;38:228-37.

    [4]HauverGE,NetherwoodPH,BenckRF,KecskesLJ.Ballistic performance of ceramic targets.Army Symposium on Solid Mechanics,USA;1993.

    [5]Hauver GE,Netherwood PH,Benck RF,Kecskes LJ.Enhanced ballistic performance of ceramic targets.19thArmy Science Conference.USA;1994.

    [6]Rapacki EJ,Hauver GE,Netherwood PH,Benck RF.Ceramics for armours-a material system perspective.7thAnnual TARDEC Ground Vehicle Survivability Symposium.USA;1996.

    [7]Lundberg P,Renstr?m R,Lundberg B.Impact of metallic projectiles on ceramic targets:transition between interface defeat and penetration.Int J Impact Eng 2000;24:259-75.

    [8]Hauver GE,Rapacki EJ,Netherwood PH,Benck RF.Interface defeat of long rod projectiles by ceramic armor.ARL-TR-3590;2005.

    [9]Holmquist TJ,Anderson CE Jr,Behner T.Design,analysis and testing of an unconf i ned ceramic target to induce dwell.22ndInternational Symposium on Ballistics,Vancouver,BC,Canada,14-18 November;2005.

    [10]Lundberg P,Lundberg B.Transition between interface defeat and penetration for tungsten projectiles and four silicon carbide materials.Int J Impact Eng 2005;31:781-92.

    [11]Lundberg P,Renstr?m R,Anderson O.Inf l uence of length scale on the transition from interface defeat to penetration in unconf i ned ceramic targets.J Appl Mech 2013;80(3):031804-1-9.

    [12]Behner T,Anderson CE Jr,Holmquist TJ,Wickert M,Templeton DW. Interface defeat for unconf i ned SiC ceramics.Proc 24thInt Symp Ballistics 2008;1:35-42.

    [13]Anderson CE Jr,Behner T,Orphal DL,Nicholls AE,Holmquist TJ,Wickert M.Long rod penetration into intact and pre-damaged SiC ceramic.Proc 24thInt Symp Ballistics 2008;2:822-9.

    [14]Anderson CE Jr,Behner T,Holmquist TJ,Orphal DL,Wickert M Dwell,Interface Defeat,and Penetration of Long Rods Impacting Silicon Carbide,Southwest Research Institute Technical Report 18.12544/008.

    [15]LaSalviaJC,HorwathEJ,RapackiEJ,ShihCJ,MeyersMA. Microstructuraland micromechanicalaspectsofceramic/long-rod projectiles interactions:dwell/penetration transitions.In:Staudhammer KP,Murr KP,Meyers MA,editors.Fundamental issues and applications of shock-wave and high-strain-rate phenomena.New York:Elsevier Science.2001 p.437-46.

    [16]Renstr?m R,Lundberg P,Lundberg B.Stationary contact between a cylindrical metallic projectile and a f l at target under conditions of dwell. Int J Impact Eng 2004;30:1265-82.

    [17]Lundberg P,Renstr?m R,Lundberg B.Impact of conical tungsten projectiles on f l at silicon carbide targets:transition from interface defeat to penetration.Int J Impact Eng 2006;32:1842-56.

    [18]Sherman D,Ben-Shushan T.Quasi-static impact damage in conf i ned ceramic tiles.Int J Impact Eng 1998;21:245-65.

    [19]Sherman D.Impact failure mechanisms in alumina tiles on f i nite thickness support and the effect of conf i nement.Int J Impact Eng 2000;24:313-28.

    [20]Savio SG,Ramanjaneyulu K,Madhu V,Bhat TB.An experimental study on ballistic performance of boron carbide tiles.Int J Impact Eng 2011;38:535-41.

    [21]Holmquist TJ,Johnson GR.Modelling prestressed ceramic and its effect on ballistic performance.Int J Impact Eng 2005;31:113-27.

    [22]Runqiang C,Ahmad S,Idapalapati S,Tan G.Pre-stress effect on conf i ned ceramic armor ballistic performance.Int J Impact Eng 2015;84:159-70.

    [23]Wereszczak AA,Johanns KE,Jadaan OM.Hertzian ring crack initiation in hot-pressed silicon carbides.J Am Ceram Soc 2009;92(8):1788-95.

    [24]Anstis GR,Chantikul P,Lawn BR,Marshall DB.A critical evaluation of indentation techniques for measuring fracture toughness:I direct crack measurements.J Am Ceram Soc 1981;64(9):533-8.

    [25]Skoglund P.Constitutive modelling of a tungsten heavy metal alloy. J Phys IV 2003;110:207-12.

    [26]Lundberg P,Renstr?m R,Holmberg L.An experimental investigation of interface defeat at extended interaction times.19thInt Symp Ballistics,Switzerland 2001;3:1463-9.

    Received 15 September 2015;revised 23 January 2016;accepted 15 February 2016 Available online 2 March 2016

    Peer review under responsibility of China Ordnance Society.

    *Corresponding author.Tel.:+46 709277154.

    E-mail address:patrik.lundberg@foi.se (P.LUNDBERG).

    http://dx.doi.org/10.1016/j.dt.2016.02.002

    2214-9147/? 2016 China Ordnance Society.Production and hosting by Elsevier B.V.All rights reserved.

    ? 2016 China Ordnance Society.Production and hosting by Elsevier B.V.All rights reserved.

    91成年电影在线观看| 亚洲中文日韩欧美视频| 精品电影一区二区在线| 欧美黄色片欧美黄色片| 制服诱惑二区| 国产在线精品亚洲第一网站| 亚洲精品成人av观看孕妇| 一a级毛片在线观看| 欧美日韩亚洲综合一区二区三区_| 欧美色视频一区免费| 国产aⅴ精品一区二区三区波| 欧美日韩精品网址| 久久天堂一区二区三区四区| 久久精品亚洲av国产电影网| 精品国产乱子伦一区二区三区| 亚洲欧美日韩另类电影网站| 午夜精品国产一区二区电影| 女性生殖器流出的白浆| aaaaa片日本免费| 久久久久国产一级毛片高清牌| 国产激情久久老熟女| 欧美久久黑人一区二区| 欧美日韩瑟瑟在线播放| 久久久精品免费免费高清| 天堂中文最新版在线下载| 老司机午夜十八禁免费视频| 成人精品一区二区免费| 免费看a级黄色片| 午夜福利免费观看在线| 女人被狂操c到高潮| 91九色精品人成在线观看| 91字幕亚洲| 啦啦啦 在线观看视频| 香蕉国产在线看| 欧美大码av| av免费在线观看网站| 亚洲成人免费av在线播放| 国产精品av久久久久免费| 男人操女人黄网站| 成人精品一区二区免费| 国产男靠女视频免费网站| 午夜福利免费观看在线| 亚洲人成伊人成综合网2020| 两性夫妻黄色片| 搡老乐熟女国产| 亚洲五月色婷婷综合| 日韩欧美一区二区三区在线观看 | www.熟女人妻精品国产| 一级黄色大片毛片| 另类亚洲欧美激情| 成人手机av| 黑人欧美特级aaaaaa片| 天天躁夜夜躁狠狠躁躁| 欧美激情高清一区二区三区| e午夜精品久久久久久久| 99久久99久久久精品蜜桃| 大型黄色视频在线免费观看| 99在线人妻在线中文字幕 | 三上悠亚av全集在线观看| 成年人免费黄色播放视频| 午夜精品久久久久久毛片777| 50天的宝宝边吃奶边哭怎么回事| 91字幕亚洲| 久久久国产精品麻豆| 午夜老司机福利片| 99久久99久久久精品蜜桃| 亚洲精品在线观看二区| 亚洲性夜色夜夜综合| 性少妇av在线| 王馨瑶露胸无遮挡在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 一边摸一边抽搐一进一小说 | 成人免费观看视频高清| 狂野欧美激情性xxxx| 欧美日韩瑟瑟在线播放| 免费在线观看完整版高清| 如日韩欧美国产精品一区二区三区| 亚洲在线自拍视频| 韩国精品一区二区三区| 国产男女超爽视频在线观看| 亚洲avbb在线观看| 亚洲熟女毛片儿| 国产一卡二卡三卡精品| 一级作爱视频免费观看| 免费一级毛片在线播放高清视频 | 麻豆乱淫一区二区| 国产精品.久久久| 精品国产一区二区久久| 国产在视频线精品| 国产蜜桃级精品一区二区三区 | 亚洲专区中文字幕在线| 一夜夜www| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲va日本ⅴa欧美va伊人久久| 国产又爽黄色视频| 青草久久国产| 热re99久久国产66热| 精品亚洲成a人片在线观看| 免费观看人在逋| 久久久久视频综合| 久久香蕉精品热| 人妻一区二区av| 波多野结衣一区麻豆| 国产成人免费无遮挡视频| 人人妻,人人澡人人爽秒播| 大陆偷拍与自拍| 成年女人毛片免费观看观看9 | 日韩熟女老妇一区二区性免费视频| 精品国产一区二区久久| 欧美日韩精品网址| 日韩三级视频一区二区三区| 国产单亲对白刺激| 麻豆成人av在线观看| 一本一本久久a久久精品综合妖精| 免费观看a级毛片全部| 日日摸夜夜添夜夜添小说| 好男人电影高清在线观看| 国产av精品麻豆| 婷婷丁香在线五月| 老汉色∧v一级毛片| 国产一区二区激情短视频| 性少妇av在线| 在线观看免费视频日本深夜| 国产精品.久久久| 99久久精品国产亚洲精品| 欧美性长视频在线观看| 日韩欧美三级三区| 国产亚洲欧美98| 老司机亚洲免费影院| 精品人妻在线不人妻| 久久性视频一级片| 欧美乱色亚洲激情| 99精品欧美一区二区三区四区| 中文字幕人妻熟女乱码| 久久久久国产精品人妻aⅴ院 | 99精品在免费线老司机午夜| 99re在线观看精品视频| 又大又爽又粗| 在线观看舔阴道视频| 国产又色又爽无遮挡免费看| 一区在线观看完整版| 亚洲国产欧美网| 免费少妇av软件| 久久亚洲真实| xxxhd国产人妻xxx| 国产一区在线观看成人免费| 麻豆国产av国片精品| 亚洲午夜理论影院| 在线观看免费视频网站a站| 在线观看免费日韩欧美大片| 国产成人免费无遮挡视频| 亚洲美女黄片视频| 最新的欧美精品一区二区| av国产精品久久久久影院| 涩涩av久久男人的天堂| 成人手机av| 亚洲av片天天在线观看| 欧美午夜高清在线| 操美女的视频在线观看| 97人妻天天添夜夜摸| 欧美黑人欧美精品刺激| 亚洲一区高清亚洲精品| 999久久久精品免费观看国产| 91成年电影在线观看| 久久精品国产99精品国产亚洲性色 | 日本欧美视频一区| 欧美日韩视频精品一区| av电影中文网址| 黄色女人牲交| 三级毛片av免费| 免费日韩欧美在线观看| 国产欧美亚洲国产| svipshipincom国产片| 12—13女人毛片做爰片一| 久久久久精品人妻al黑| 国产主播在线观看一区二区| 女性生殖器流出的白浆| 亚洲精品中文字幕一二三四区| 国产极品粉嫩免费观看在线| 亚洲在线自拍视频| 久久国产精品影院| 视频区欧美日本亚洲| 无人区码免费观看不卡| 免费日韩欧美在线观看| 18禁国产床啪视频网站| 99香蕉大伊视频| 亚洲欧美一区二区三区久久| 国产91精品成人一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | tube8黄色片| 99国产精品一区二区蜜桃av | 女性生殖器流出的白浆| 精品久久蜜臀av无| 一边摸一边做爽爽视频免费| 久久人妻福利社区极品人妻图片| av天堂久久9| 中文字幕制服av| 少妇裸体淫交视频免费看高清 | 国产99久久九九免费精品| av线在线观看网站| 亚洲av熟女| 国产精品一区二区免费欧美| 丁香欧美五月| av福利片在线| 国产乱人伦免费视频| 亚洲熟妇中文字幕五十中出 | 91九色精品人成在线观看| 国产淫语在线视频| 三级毛片av免费| 色老头精品视频在线观看| 日韩欧美三级三区| 中出人妻视频一区二区| 国产高清视频在线播放一区| 日韩人妻精品一区2区三区| 国产99久久九九免费精品| 一级a爱片免费观看的视频| 男人的好看免费观看在线视频 | 色在线成人网| 欧美日韩视频精品一区| 成年人午夜在线观看视频| 欧美黑人欧美精品刺激| 成熟少妇高潮喷水视频| 午夜免费成人在线视频| 两性夫妻黄色片| 欧美乱码精品一区二区三区| 国产精品免费一区二区三区在线 | 国产成人欧美| 久久精品91无色码中文字幕| 国产亚洲精品一区二区www | 两性午夜刺激爽爽歪歪视频在线观看 | 国产视频一区二区在线看| 亚洲熟妇熟女久久| 国产精品久久电影中文字幕 | 十八禁人妻一区二区| 久久 成人 亚洲| av中文乱码字幕在线| 电影成人av| 真人做人爱边吃奶动态| 久久中文看片网| 日本vs欧美在线观看视频| 麻豆成人av在线观看| 久久久国产精品麻豆| 午夜成年电影在线免费观看| 涩涩av久久男人的天堂| 午夜精品久久久久久毛片777| 欧美精品啪啪一区二区三区| 人人妻人人澡人人爽人人夜夜| 免费av中文字幕在线| 亚洲伊人色综图| 亚洲精品自拍成人| 欧美性长视频在线观看| 成人亚洲精品一区在线观看| 亚洲情色 制服丝袜| 国产又爽黄色视频| 亚洲人成伊人成综合网2020| 亚洲精品国产区一区二| 少妇裸体淫交视频免费看高清 | 国产精品久久久久成人av| 日本精品一区二区三区蜜桃| www.精华液| 嫩草影视91久久| 一级,二级,三级黄色视频| 色播在线永久视频| 天堂√8在线中文| 校园春色视频在线观看| videos熟女内射| 成人黄色视频免费在线看| www日本在线高清视频| 99精品欧美一区二区三区四区| 可以免费在线观看a视频的电影网站| 精品国产亚洲在线| 亚洲国产欧美一区二区综合| 国产亚洲精品一区二区www | 欧美日韩国产mv在线观看视频| 宅男免费午夜| 人妻久久中文字幕网| 悠悠久久av| 国产伦人伦偷精品视频| 日日摸夜夜添夜夜添小说| 日本欧美视频一区| 国产精品秋霞免费鲁丝片| 国产亚洲精品久久久久久毛片 | 一区二区三区激情视频| 亚洲国产欧美一区二区综合| 亚洲一区中文字幕在线| 精品国产乱码久久久久久男人| 两个人看的免费小视频| 人妻 亚洲 视频| 日韩一卡2卡3卡4卡2021年| 激情视频va一区二区三区| 国产成人精品无人区| 99国产精品99久久久久| 亚洲欧美精品综合一区二区三区| 黄频高清免费视频| 亚洲欧美一区二区三区久久| 91成年电影在线观看| 国产日韩一区二区三区精品不卡| 国产主播在线观看一区二区| 一级片'在线观看视频| 日韩欧美国产一区二区入口| 国产极品粉嫩免费观看在线| 国产高清视频在线播放一区| 成人黄色视频免费在线看| 久热爱精品视频在线9| 日韩欧美免费精品| 亚洲国产欧美网| 中文字幕人妻丝袜一区二区| 国产精品一区二区在线不卡| 日韩欧美一区视频在线观看| 欧美日韩福利视频一区二区| 黄片播放在线免费| 999久久久精品免费观看国产| 亚洲第一av免费看| 99在线人妻在线中文字幕 | 999精品在线视频| videos熟女内射| 久久久久精品人妻al黑| av网站免费在线观看视频| 两人在一起打扑克的视频| 亚洲成人手机| 极品少妇高潮喷水抽搐| 午夜福利免费观看在线| 三级毛片av免费| 成人永久免费在线观看视频| 丁香欧美五月| 老鸭窝网址在线观看| 正在播放国产对白刺激| 国产99白浆流出| 十八禁网站免费在线| 一区在线观看完整版| 久久国产亚洲av麻豆专区| 亚洲国产精品合色在线| 一级a爱片免费观看的视频| 亚洲人成77777在线视频| 一级,二级,三级黄色视频| 99精品在免费线老司机午夜| 男女床上黄色一级片免费看| 最近最新中文字幕大全电影3 | 亚洲免费av在线视频| 国产深夜福利视频在线观看| 99久久国产精品久久久| 国产成人欧美| 免费在线观看黄色视频的| 国产精品久久久久久精品古装| 在线观看午夜福利视频| 亚洲精品成人av观看孕妇| 国产精品久久久久久人妻精品电影| 无限看片的www在线观看| 90打野战视频偷拍视频| 色综合婷婷激情| 午夜福利免费观看在线| 精品久久久精品久久久| 免费女性裸体啪啪无遮挡网站| netflix在线观看网站| 亚洲av电影在线进入| 精品少妇久久久久久888优播| 美女高潮喷水抽搐中文字幕| 精品第一国产精品| 日韩欧美一区视频在线观看| 国产人伦9x9x在线观看| 一本一本久久a久久精品综合妖精| 亚洲第一欧美日韩一区二区三区| 国产野战对白在线观看| 首页视频小说图片口味搜索| 精品电影一区二区在线| 免费在线观看亚洲国产| 国产日韩欧美亚洲二区| 国产精品亚洲av一区麻豆| 精品人妻1区二区| 久久久精品国产亚洲av高清涩受| 老熟女久久久| 欧美黄色淫秽网站| 天天躁狠狠躁夜夜躁狠狠躁| 中出人妻视频一区二区| 久久精品aⅴ一区二区三区四区| 老司机福利观看| 天天影视国产精品| av网站在线播放免费| 日本欧美视频一区| 熟女少妇亚洲综合色aaa.| 99re在线观看精品视频| 高清在线国产一区| 亚洲成人国产一区在线观看| 国产精品一区二区免费欧美| 老汉色av国产亚洲站长工具| 精品高清国产在线一区| 国产精品自产拍在线观看55亚洲 | 亚洲国产欧美一区二区综合| 国产亚洲精品久久久久久毛片 | 极品少妇高潮喷水抽搐| 精品欧美一区二区三区在线| avwww免费| 黑人猛操日本美女一级片| 日韩欧美一区视频在线观看| 两个人看的免费小视频| 国产99白浆流出| 欧美亚洲 丝袜 人妻 在线| 国产亚洲av高清不卡| 男人的好看免费观看在线视频 | 国产免费av片在线观看野外av| 一进一出抽搐动态| 国产精品一区二区精品视频观看| 国产精品av久久久久免费| 精品乱码久久久久久99久播| a在线观看视频网站| 99国产精品99久久久久| 两个人免费观看高清视频| 亚洲成人免费av在线播放| 亚洲av日韩在线播放| 人人妻人人添人人爽欧美一区卜| 天天影视国产精品| 亚洲欧美激情综合另类| 一二三四在线观看免费中文在| 亚洲av成人av| 精品午夜福利视频在线观看一区| 一进一出抽搐gif免费好疼 | 极品教师在线免费播放| 久久久久久久久免费视频了| avwww免费| 99热网站在线观看| 欧美日韩黄片免| 啪啪无遮挡十八禁网站| 精品久久久精品久久久| 美女高潮到喷水免费观看| 黄色丝袜av网址大全| 丁香欧美五月| 亚洲欧美激情综合另类| 老司机靠b影院| 在线观看午夜福利视频| 国产精品国产高清国产av | 高清欧美精品videossex| 亚洲av第一区精品v没综合| 777久久人妻少妇嫩草av网站| 国产成人影院久久av| 激情视频va一区二区三区| 中文字幕人妻熟女乱码| 极品少妇高潮喷水抽搐| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人国产一区在线观看| 看片在线看免费视频| 国产精品av久久久久免费| 亚洲一区中文字幕在线| 久久性视频一级片| 国产精品影院久久| 欧美一级毛片孕妇| 亚洲成人国产一区在线观看| 真人做人爱边吃奶动态| 日本a在线网址| 水蜜桃什么品种好| 中文欧美无线码| 三上悠亚av全集在线观看| x7x7x7水蜜桃| 国产xxxxx性猛交| 成人18禁在线播放| www.999成人在线观看| 成人特级黄色片久久久久久久| 777米奇影视久久| 天天影视国产精品| 大片电影免费在线观看免费| 午夜福利,免费看| 男男h啪啪无遮挡| 日韩欧美免费精品| 黑人巨大精品欧美一区二区mp4| 中文欧美无线码| 黑人巨大精品欧美一区二区蜜桃| 在线av久久热| 国产精品久久久人人做人人爽| 黑人欧美特级aaaaaa片| 久久精品人人爽人人爽视色| 成年人午夜在线观看视频| 男人的好看免费观看在线视频 | 99精国产麻豆久久婷婷| 97人妻天天添夜夜摸| 国产精品偷伦视频观看了| 午夜视频精品福利| 在线av久久热| 1024香蕉在线观看| 真人做人爱边吃奶动态| 大型av网站在线播放| 无限看片的www在线观看| 午夜福利在线观看吧| 婷婷成人精品国产| 国产成人免费观看mmmm| 一级a爱片免费观看的视频| 成人18禁在线播放| 成人免费观看视频高清| 女性被躁到高潮视频| 啪啪无遮挡十八禁网站| 激情视频va一区二区三区| 午夜精品在线福利| av中文乱码字幕在线| 国产精品 欧美亚洲| 91精品三级在线观看| 一区二区三区国产精品乱码| 满18在线观看网站| 丝袜在线中文字幕| 久久久久国产一级毛片高清牌| 欧美午夜高清在线| 成人精品一区二区免费| 国产成人精品无人区| 不卡一级毛片| 久久久久久亚洲精品国产蜜桃av| 午夜福利一区二区在线看| 人人澡人人妻人| 老司机影院毛片| 黄色成人免费大全| 国产色视频综合| videos熟女内射| 欧美精品av麻豆av| 久久国产精品大桥未久av| 曰老女人黄片| 久久久国产成人免费| 丁香欧美五月| 精品少妇久久久久久888优播| 90打野战视频偷拍视频| 久热爱精品视频在线9| 波多野结衣一区麻豆| 亚洲国产精品sss在线观看 | 国产免费男女视频| 超色免费av| 精品人妻1区二区| 黄色丝袜av网址大全| 91成年电影在线观看| 九色亚洲精品在线播放| 午夜老司机福利片| 午夜福利乱码中文字幕| 美女福利国产在线| 夜夜爽天天搞| ponron亚洲| 亚洲精品国产一区二区精华液| 十八禁网站免费在线| 欧美性长视频在线观看| 国产欧美日韩一区二区精品| 99re在线观看精品视频| 中文字幕av电影在线播放| 狠狠狠狠99中文字幕| 丝袜在线中文字幕| 国产亚洲欧美精品永久| 大香蕉久久网| 下体分泌物呈黄色| 国产成人一区二区三区免费视频网站| 黄色女人牲交| 咕卡用的链子| a在线观看视频网站| 男女午夜视频在线观看| 国产精品一区二区在线不卡| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美日韩另类电影网站| 两人在一起打扑克的视频| 在线免费观看的www视频| 亚洲欧美一区二区三区久久| 亚洲一区二区三区欧美精品| videos熟女内射| 黄片播放在线免费| 亚洲精品中文字幕一二三四区| 飞空精品影院首页| 国产成人啪精品午夜网站| 久久久久久久午夜电影 | 国产av一区二区精品久久| 少妇猛男粗大的猛烈进出视频| 啪啪无遮挡十八禁网站| 99久久人妻综合| 91大片在线观看| 精品国产一区二区三区久久久樱花| 国产视频一区二区在线看| 最新在线观看一区二区三区| 91精品三级在线观看| 国产精品国产高清国产av | 亚洲国产看品久久| 少妇猛男粗大的猛烈进出视频| 最近最新免费中文字幕在线| 亚洲少妇的诱惑av| 国产日韩一区二区三区精品不卡| 日本a在线网址| 久久久久久久国产电影| av一本久久久久| 天天添夜夜摸| 一区二区三区激情视频| 777米奇影视久久| 女人久久www免费人成看片| www日本在线高清视频| 女警被强在线播放| 国产一卡二卡三卡精品| 午夜福利影视在线免费观看| 国产熟女午夜一区二区三区| 日韩欧美一区视频在线观看| 咕卡用的链子| 一区在线观看完整版| 久热这里只有精品99| 色播在线永久视频| 久久九九热精品免费| 一a级毛片在线观看| 色综合婷婷激情| 欧美精品一区二区免费开放| 啦啦啦免费观看视频1| 男男h啪啪无遮挡| 国产精品偷伦视频观看了| 午夜福利欧美成人| 看片在线看免费视频| 一区二区三区精品91| 亚洲国产毛片av蜜桃av| 久久国产亚洲av麻豆专区| 欧美日韩成人在线一区二区| 国产精品久久久av美女十八| 成人精品一区二区免费| 91国产中文字幕| 巨乳人妻的诱惑在线观看| 黄色毛片三级朝国网站| 黄片大片在线免费观看| 一a级毛片在线观看| 91在线观看av| 69精品国产乱码久久久| 久久久精品免费免费高清| 又黄又粗又硬又大视频| 人人妻人人澡人人看|