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

    Fracture mechanism of a cylindrical shell cut by circumferential detonation collision

    2021-11-03 13:23:58HeyngXuWeiingLiWeninLiQingZhngjunWngXiowenHongYingYng
    Defence Technology 2021年5期

    He-yng Xu ,Wei-ing Li ,* ,Wen-in Li ,Qing Zhng ,Y-jun Wng ,Xio-wen Hong ,Ying Yng

    a Nanjing University of Science and Technology,Nanjing,210094,China

    b Jiangsu Yongfeng Mechanics Co.,Ltd,Nanjing,210094,China

    Keywords: Detonation wave Circumferential detonation collision Fracture mechanism Fragmentation control

    ABSTRACT The failure mechanism of a cylindrical shell cut into fragments by circumferential detonation collision was experimentally and numerically investigated.A self-designed detonation wave regulator was used to control the detonation and cut the shell.It was found that the self-designed regulator controlled the fragment shape.The macrostructure and micro-characteristics of fragments revealed that shear fracture was a prior mechanism,the shell fractured not only at the position of detonation collision,but the crack also penetrated the shell at the first contact position of the Chapmen-Jouguet (C-J) wave.The effects of groove number and outer layer thickness on the fracture behavior were tested by simulations.When the thickness of the outer layer was 5-18 mm,it has little effect on fragmentation of the shell,and shells all fractured at similar positions.The increase of the groove number reduced the fracture possibility of the first contact position of the C-J wave.When the groove number reached 7 with a 10 mm outer layer(1/4 model),the fracture only occurred at the position of detonation collision and the fragment width rebounded.

    1.Introduction

    The propagation of detonation waves is an important research topic in the field of explosive mechanics.With the increasing diversity of initiation methods,the collision of detonation waves has become common.The energy and the pressure converge and jump in a narrow area when two detonation collision.Due to the obvious flexibility and directionality,this accumulation effect of detonation collision is widely studied in basic research and military applications.

    Dunne et al.[1] first theoretically deduced the calculation method of the Mach rod pressure when two waves collided in a condensed explosive.Ivanov [2] conducted relevant experiments on the deformation of inert materials by detonation collision to verify the accumulation phenomenon.Müller F [3].studied the detonation interaction process in coaxial charges and captured the shape of the Mach disk by X-ray photography.However,Ivanov and Müller F.both did not measure the definite gain,although the enhancement effect of collision was experimentally proved.Later,with more advanced technology,researchers could accurately measure and capture the process of detonation interaction and energy changes.Tang Mingjun et al.[4]measured the pressure and velocity of detonation in each section with GSJ high-speed photography and found that the velocity of the Mach wave was 1.9 times the Chapman-Jouguet (C-J) detonation velocity and the pressure was 3 times the C-J detonation pressure.WEN Shang-gang et al.[5,6]and A V Shutov et al.[7,8]also examined the interactions between multiple detonations and confirmed that the pressure and velocity in the collision area were greater than those of the C-J wavefront.

    Detonation collision is widely studied in basic research and military applications due to the obvious flexibility and directionality.It can be applied to explosive shaped projectiles to achieve better flight stability [9-11] or can be applied to fragment warheads in a multi-point eccentric initiation method to achieve higher dispersion speed and better ejection angle [12-14].In addition,detonation collision can also be used to directly cut low ductility metals.W.Arnold et al.[15,16] proposed a method called high explosive pellet and applied the pressure enhancement effect at the collision area to cut shells or liners to generate multiple fragments;however,they mainly focused on the use of some devices to complete the exchange between different damage elements.Zipu Geng [17] simulated the fracture failure behavior of a liner resulting from the axial detonation collision by multi-point initiation and noticed that the detonation wave became superimposed on the liner surface and formed concentration zones with a strong stress,leading to fracture along these zones;however,no corresponding test was conducted to verify the result.Zhang K.F[18,19].conducted a similar test to further study the fracture mechanism of a steel plate cut by multi-detonation superposition.The fragment and simulation results were discussed and it was found that the steel plate shear fractured at local high-stress zones produced by the detonation collision.

    However,most of these studies have mainly focused on cutting metal plates by axial detonations and rarely studied the fracture mechanism of cylindrical shells under circumferential detonation collision.When the charge is detonated,a cylindrical shell becomes exposed to a large strain in the circumferential direction.The cylindrical shell then breaks under the combined effects of this large tensile strain and detonation collision,and this phenomenon does not exist in axial metal cutting.In the present paper,the fracture mechanism of a cylindrical shell under circumferential detonation collision is discussed.

    The soft recovery test was carried out for a cylindrical shell(with a self-designed detonation wave regulator) filled with explosives.The shape,quantity,and mass of fragments generated from the warhead with or without a regulator were compared.The fracture surfaces and shapes of fragments generated from the warhead with a regulator were analyzed by microscopic and macroscopic methods.The fracture mechanism of the cylindrical shell under circumferential detonation collision was subsequently studied.Further,the formation process of local high-pressure regions was analyzed by a finite element (FE) model.Finally,the influences of groove number and outer layer thickness on the fracture behavior of the cylindrical shell were investigated by the simulation model.

    2.Experimental procedure

    2.1.Self-designed regulator and action principle

    In order to achieve a circumferential detonation collision,it is necessary to generate independent circumferential multi-waves.There is a high demand for initiation synchronization to generate circumferential multi-waves;however,it has a big impact on the obtained results.Similarly,the multi-point initiation technology cannot produce independent circumferential multi-waves reliably.Therefore,a detonation wave regulator was designed to realize the collision and superposition of circumferential multi-waves by a simple,single detonator initiation at the end face.The structure of the regulator is displayed in Fig.1.The regulator was a ring-shaped cylinder with an outer diameter of 60 mm and a thickness of 10 mm.It was made of phenolic resin,which has a small density and good explosion-proof performance.The sidewall was processed with evenly distributed grooves that ran through the entire thickness.When the charge was detonated,the detonation propagated outward,then became blocked,and finally was divided by the regulator.A part of the wave propagated through the grooves to the outer layer,and new wavelets were formed.Adjacent wavelets collided with each other and generated localized high-pressure regions on the shell.As a result of that,the shell would experience different distribution of pressure and thus produced different effects.

    2.2.Experimental setup

    In the test,the warhead structure shown in Fig.2(a) was used.The shell was a 6 mm-thick AISI 1045 steel cylindrical tube (outer diameter=92 mm and height=80 mm),and its two ends were connected with a ring (height=10 mm) made of phenolic resin with the same diameter(the ring was used to reduce the impact of axial rarefaction waves on the shell).Different from the common structure,a self-designed regulator was placed coaxially in the shell to control the detonation wave,divided the explosive into inner charge and outer layer.The charge was made of Comp B with the shapes of the shell and the regulator,and the grooves were filled evenly to connect the explosives on both sides.A cover was used to fixing the detonator and the booster explosive at the top,the charge was point-initiated with a booster explosive from the top center.

    In order to retain more comprehensive fragment information and reduce the secondary fragmentation as much as possible,water was used as a soft medium to recover all fragments.Fig.2(b) presents the experimental setup that includes the specimen,an air chamber,and a recovery device.During the test,the specimen was placed in the middle of the air chamber and the water well,and a fine net was used to recover fragments at the bottom.Three warheads were used (Table 1):the treatment group had two repeated specimens with the structure shown in Fig.2(a),and the control group had a specimen without a regulator and filled with Comp B.

    Table 1 Experimental conditions.

    3.Experimental results and analysis

    3.1.Mass distribution of fragments

    All recovered fragments were cleaned,dried,and sorted,and as a result,Fig.3 displays partial typical fragments of the treatment group and the control group.It is evident that obvious differences in the fragment shape and size existed between these two groups.The fragment shape of the control group was normal and common,and the size was smaller than that of the treatment group(Fig.3(a)).In the treatment group,the fragment shape was more regular (Fig.3(b));however,two typical fragment shapes were observed:(i) large fragments with a regular shape (the overall shape tended to be rectangular),and fracture surfaces were oriented at~45°to the inner surface and (ii) long and elongated fragments;one side of these fragments was smooth,whereas the other side had a ridged surface.

    In order to understand the impact of the regulator comprehensively and also to prove that the above-mentioned large shape difference was not caused by accidental factors or artificial selection,all recovered fragments were weighed and counted (fragments of mass>1 g were considered).The mass distributions of the three cases are displayed in Fig.4.It is evident that when a regulator was in charge,it had a great impact on shell fragmentation,regardless of quantity or mass.

    The fragment mass distribution of case 3 (without a regulator,Fig.4(d)) is consistent with the known results of the shell:The fragmentation was relatively uniform,small fragments make up the majority,fragments were mainly concentrated in the range 0-4 g(accounting for~90% of the total number of fragments),the maximum fragment mass was less than 8 g,and the number of fragments decreased exponentially with the increase of the mass.In comparison to case 3,cases 1 and 2(with a regulator)manifested an obvious difference in mass distribution.The number of small mass fragments was much smaller than that without a regulator.Fragments of mass 1-2 g generated in the detonation cutting process or secondary breaking were not the main part.Moreover,an abnormal upward trend in the partition statistics was observed rather than a significant decrease in the fragment number with the increase of the fragment mass(Fig.4(a)).It is evident from Fig.4(b and c) that fragments were evenly distributed in the range of 5-20 g.These fragments were the main part from the shell cut and several individual fragments with a mass close to 30 g appeared as connected fragments.The regulator only processed grooves in the axial direction.Therefore,the shell became stretched under rarefaction waves (although two phenolic resin rings were connected with the shell to reduce rarefaction,it was not enough to completely offset the effect),and the fracture at different horizontal positions resulted in the fragment mass difference.

    Fig.1.Structure of the regulator:(a) 3D partial section view,(b) and (c) are the side view and the A-A section view,respectively.

    Fig.2.(a)Side view and the A-A section view of the specimen with a detonation wave regulator,(b)Experimental setup of water well(φ8 m×8 m)recovery,the centers of the well,the air chamber,and the specimen coincide and (c) Photos of the specimen with a regulator.

    Fig.3.Partial experimental results(a large grid is 10 mm):(a)fragments of the control group with a normal and common shape and(b)fragments of the treatment group with two typical shapes.

    Fig.4.Statistical results of fragment mass.

    By comparing the results of mass distribution and typical fragment shape,it is confirmed that the addition of the regulator had a significant impact on the fracture behavior of the cylindrical shell.Although the charge mass of the control group was larger than that of the treatment group,it could not impact the distribution trend and resulted in these two extremely different mass distributions and fragment morphologies.This happened because of the special detonation propagation path caused by the regulator rather than a larger charge mass.Therefore,it can be inferred that the shell could be cut into large and regular fragments by detonation after adding the regulator.

    The addition of the regulator completely changed the detonation propagation path,shown as Fig.5(a).After the detonation spread through regulator,formed wavelets,and collided at point I.With time,the high-pressure region moved to the shell at region M (M refers to the superposition position on the shell).Due to such high pressure,region M experienced a large deformation.In the subsequent expansion,the cylindrical shell underwent shear instability under a high strain.The shear band initiated at region M due to the large deformation and then grew into cracks through the shell,forming long strip fragments(shown as the red area in Fig.5(b)).The shell may also be fractured where the detonation first contacted position (region D) and formed regular large fragments (shown as the blue area in Fig.5(b));this will be discussed in the 3.2 section.

    3.2.Mechanism of waves cut the shell

    The macrostructure and microfracture characteristics of typical fragments in the treatment group were analyzed to understand the fracture mechanism of the cylindrical shell under circumferential detonation collision.Fig.6 presents the Scanning Electron Microscope (SEM) images of the fractured surface of a long-ridged fragment,and three different micro-characteristics were observed.Region I is the ridge part,which is the inner surface of the shell,and Fig.6(A)displays the enlarged microscopic feature of this part.This area presented a water ripple characteristic with obvious directionality.The direction of the microstructure was perpendicular to the direction of the fragment fracture,and similar characteristics were not observed on the inner surface of the shell with a normal charge.This happened because the pressure of the detonation collision was more than double the C-J pressure.Under such high pressure,the inner surface metal softened and flowed along the direction of detonation propagation.When the detonation collision area passed through the shell,the large pressure difference produced great shear stress and strain on both sides of region M.Severe plastic deformation and heat accumulation occurred in the metal,which reduced the shear strength of the material,slips occurred along a certain crystallographic plane and direction in the structure.Then micro-cracks appeared at region M and expanded rapidly to both sides,formed the elongated fragments with a ridged surface.During cracks propagated,fracture surfaces became squeezed and rubbed against each other,resulting in flat surfaces with obvious textures,consistent with Fig.6(B),which is indicated by a series of small bright surfaces macroscopically.However,the high pressure on region M only existed in a short period,and the detonation pressure gradually tended to be uniform during the subsequent expansion process.Therefore,region II did not cover the whole fracture surface,rather the uneven surface with elliptical dimples appeared in region III (Fig.6(C)).It is a typical shear fracture characteristic of the thin cylindrical shell without a regulator fractured at a high strain rate.Through the micro-characteristics,it can be known that a large compressive stress acted on ridged fragments,and embryos appeared after the passage of detonation collision.The initiation of embryos and the crack propagation occurred earlier than normal.These characteristics confirm that ridged fragments were formed in the position of detonation collision (Region M).

    The fracture in region D was not expected.The test results showed that the width of regular fragments(upper row in Fig.3(b))in the treatment group was smaller than expected.The width of all regular fragments ranged between 1.3 cm and 1.6 cm.Assuming that the initial radius and thickness of the shell before expansion area0and δ0,respectively,the radius and thickness of the shell after fracture areafand δf,respectively;in addition,εsatisfiesaf=εa0due to the volume of the shell being constant,and then [20].

    where,ν is shape factor,for cylindrical shell ν=1,ignoring highorder small quantities of δ0and δf,the thickness of the cylindrical shell at the moment of fracture is δf=δ0/ε.Ten fragments were selected randomly,and their thicknesses were measured.The average fragment thickness was calculated as 4.22 mm;thus ε=1.42.Further,the shell radius at fracture was calculated asaf=1.42× 4=5.68 cm.If the fracture only occurs at region M,the shell should be evenly broken into 12 parts and the width of each fragment is about 2.97 cm (without considering the loss of the shell).This value has a great deviation from the actual measured value (~2 times the actual value).

    Further,four connected fragments in the treatment group were measured and observed.All four pieces were located at the end of the shell because of the low pressure caused by the axial rarefaction wave.The inner surface and the section are displayed in Fig.7.It is clear that the section had a trapezoidal shape (the width of the outer surface was smaller than that of the inner surface),and the widths of these four fragments were 2.72 cm,2.81 cm,2.50 cm,and 2.73 cm.Those results are consistent with the calculated values and the fragment shape between adjacent detonation collision positions shown in Fig.5(b).What’s more,obvious indentations and incomplete cracks were noticed in the middle of the inner surface of these four connected fragments.Therefore,it is proved that the fracture also occurred at the midpoint between adjacent detonation collision positions (region D),We guess the reason of the indentation is due to the presence of an adjusted wavefront with a large curvature,and the indentation directly leads to the generation of crack embryos.This is a unique phenomenon that cannot be observed in the process of axial detonation wave cutting.

    Fig.5.(a) Detonation wave propagation and interaction process after the addition of the regulator and (b) Formation of two types of fragments in the treatment group.

    Fig.6.SEM images of the long-ridged fragment.

    Fig.7.Characters for four connected pieces.

    Fig.8(a) displays the fracture micro-characters of region D;it presents a serpentine sliding feature and featureless flat surfaces.Macroscopically,the fracture was relatively bright with obvious plastic deformation traces,and the direction was about 45°to the inner surface,thus revealing a typical shear fracture feature.This proves that region D was fractured after the shell had enough expansion,similar to the fragmentation of normal cylindrical shell under explosively loads,and later than region M.Fig.8(b)displays a metallograph of the thickness section of regular fragments.Cavities and micro-cracks can be observed,but these cavities were not connected to a through crack.This is due to the earlier initiation and propagation of cracks on both sides,namely regions M and D,and the released waves from both sides shield the crack growth.Hence,the width of regular fragments depends not only on the positions of region D and M,but also on the released waves speed from those positions.

    Therefore,the fracture process of the cylindrical shell with a regulator can be summarized as follows:the shell first deformed at regions M and D after the charge was detonated and generated a large number of important embryos.Embryos at region M developed into cracks earliest due to the high pressure.As the shell moved outward rapidly and deformed at a high strain rate,embryos at region D were gradually developed into cracks through the shell.Shear fracture is a prior mechanism for regions M and D.Microcracks in other locations were shielded by released waves and eventually could not penetrate the entire thickness.

    4.Numerical simulation

    4.1.Analysis of influencing factors

    It would be time-consuming to discuss different influencing factors for circumferential detonation shell cutting through individual experiments.Therefore,a finite element software was used to reconstruct the above fracture process and simulated the influences of different factors on the fracture behavior of the shell.

    When two waves propagate at an angle,the wave oblique collision occurs,and the pressure at the collision point of detonation increases.The entire process can be simplified as the detonation collides with the rigid wall at a certain angle,resulting in oblique shock waves after the collision.According to the conventional oblique collision theory and the unconventional oblique collision (Mach collision) theory,the pressure at the collision position [21] can be calculated as

    Fig.8.(a) SEM image of fracture surface at region D and (b) Metallograph of the section of a regular fragment.

    whereP2andPHare the collision position pressure and the explosion pressure,respectively,φis the incident angle,θis a parameter related to φ,φcis the critical angle that determines whether the collision belongs to the normal oblique collision or Mach collision,andkis the adiabatic index.

    When the explosive type is determined and the change ofkis ignored,P2/PHat the collision point is a univariate function of φ,P2/PHonly changes with φ.According to Fig.5(a),the incident angle φ depends on the distance between adjacent grooves and the thickness of the outer layer between the adjuster and the shell.Therefore,the groove number and the outer layer thickness can change the collision position pressure and affect the cutting effect of detonation.Hence,in the following sections,the influences of groove number and outer layer thickness on shell cutting by circumferential detonation collision are discussed.

    4.2.Finite element model

    The AUTODYN software was used to simulate the cylindrical shell with the regulator based on the multi-material fluid-solid coupling algorithm.Considering the large deformation of explosive and regulator elements,the 3D Multi-material Euler solver was used,and the Lagrange solver was used to the shell.The unit was gcm-us.The statistics of the regular fragment length revealed that the length distribution obeyed a step-like nature distribution,which signifies that the length could be quantized,and the average length was 4.11 cm.It could be assumed that no axial fracture occurred within 4 cm and the model had a symmetrical threedimensional geometry and load.Therefore,a 1/4 model with a length of 4 cm was established to reduce unnecessary calculation time,and the other sizes are same as the specimen of treatment group.The proposed finite element model,air not shown,is displayed in Fig.9.Mesh sizes play a significant role in computation efficiency and accuracy,the average mesh sizes for the specimen and the shell were 0.5 mm and 0.3 mm,respectively[22].

    The Johnson-Cook model [23] was applied to describe the deformation behavior of the specimen.This material model was sometimes used when the strain rate varied over a large range and the adiabatic temperature increased due to plastic heating,causing material softening.The equation of the material model can be expressed as

    where σyis the yield stress,is the effective plastic strain,ε*is the normalized effective plastic strain,T*=(T-Tr)(TM-Tr)is the homologous temperature,TMandTrare the melting temperature and the room temperature,respectively,andA0,B0,C,n,andmare constants.The Johnson-Cook failure model was used for the shell and the strain at failure can be expressed as

    where εfis the strain at failure,σ*is the ratio of pressure divided by effective stress,D1-D5are failure parameters.

    Fig.9.Finite element model.

    For the high explosive model,the Jones-Wilkins-Lee(JWL)[22]method was used to describe the detonation of Comp B.The pressure field was defined as a function of relative volume and internal energy per initial volume.

    wherePis the detonation pressure,A,B,R1,R2,and ω are constants,andE’is specific KE.

    The Von Mises model and the shock equation were applied to describe phenolic resin,and air used ideal gas equation with a density 0.001225 g/cm3,and an adiabatic exponent γ=1.4.The parameters used in simulations are listed in Table 2.

    4.3.Model validation

    The obtained simulation results are presented in Fig.10.Fig.10(a)displays the gauges and pressure distribution on the inner surface of the shell at 6.45 μs.The red and yellow areas indicate detonation collision areas,and they propagated along the entire length in the vertical direction.Five gauges were set on the shell at the same separation distance to further observe the pressure distribution,and Fig.10(b) displays the pressure history of these gauges.The purple line represents the interaction point of dual detonations.The pressure started to rise at 6.3 μs and then rapidly jumped to 65.79 GPa (about 2.5 times the C-J pressure) at 6.5 μs.Consequently,region M experienced a large deformation due to a huge pressure difference of~40 GPa.

    The simulation result of shell cutting at 60 μs is displayed in Fig.11.It is clear that numerical and test results were the same.The cylindrical shell fractured at regions M and D and was cut into six large fragments and three long-ridged fragments (1/4 model).To reduce measurement error,the inner width and outer width of fragments at three different sections were measured.Widths of six large fragments were 1.66 cm,1.44 cm,1.32 cm,1.34 cm,1.40 cm and 1.61 cm,the maximum error with the actual width(~1.4 cm)is 15.6%,and the average error is less than 10%.Therefore,it can be inferred that the model accurately reflected the fragmentation of the shell with a regulator.Subsequent simulations will be carried out base on this model.

    4.4.The effects of groove number

    The number of grooves on the regulator seriously affected the fragment number and fracture location of the cylindrical shell.The fracture mechanism of the shell was simulated as the groove number increased from 3 to 7 (1/4 model).Table 3 displays the fragmentation of the cylindrical shell with different groove numbers,red arrows indicate the positions of detonation collision and black arrows indicate the action positions of the C-J wave.

    Table 2 Material parameters for the simulation [24-27].

    Table 3 Simulation of the fragmentation of the cylindrical shell for different groove numbers.

    It is clear that the detonation superposition effectively cut the shell;however,the fracture at region D was affected by the groove number.When the number of grooves increases from 3 to 5,both regions M and D fractured and the fragment width decreased with the increase of the fracture position.When the groove number was 6,the partial region D generated connected fragments.When the groove number increased to 7,region D was not cut and the fragment width increased.The groove number did not affect the cutting of the shell at the detonation collision position;however,it weakened the effect of the C-J wave.With the increase of the groove number,the distance between adjacent detonation collision became smaller,and released waves spread to the entire distance in a shorter time,thus shielding the crack growth.Simultaneously,the smaller the curvature radius of the C-J wave also reduces the deformation at these regions.

    4.5.The effects of outer layer thickness

    The outer layer thickness affected the propagation of detonation and generated different pressures at the detonation collision position.When metal is cut axially,the distance between the initiation point and the metal plane seriously affects the cutting effect[19,28].In order to study the effect of outer layer thickness on the fracture behavior of the cylindrical shell,the fragmentation of the shell with the outer layer thickness of 5-18 mm was simulated.

    Fig.12 displays the peak pressures of detonation collision at the same position on the inner surface and the fracture of the shell with the increasing outer layer thickness from 5 mm to 18 mm.It is evident from Fig.12 that although the peak pressure of detonation collision changed greatly with the increase of the outer layerthickness,the fracture position of the shell did not change,even when the thickness was only 5 mm and the peak pressure was only 45.17 GPa(1.7 times the C-J pressure),the shell still fractured at the similar position.It is because that the shell bore a large circumferential strain when the explosive expanded violently,which is much larger than the axial strain.The detonation did not need to cut the shell completely.As long as defects are produced,they will develop into cracks under stress.We also simulated the cutting process of the 10-mm-thick outer layer with horizontally distributed grooves,and it was found that the shell could not be cut along high-pressure regions without a huge strain.

    Fig.10.Numerical results:(a) Distributions of five gauges and pressure at 6.45 μs and (b) Pressure history of gauges.

    Fig.11.Simulation result of cutting shell at 60 μs.

    Fig.12.Peak pressures of interaction on the shell and simulation results for outer layer thicknesses from 5 mm to 18 mm.

    5.Conclusion

    1.An annular detonation wave regulator with grooves was designed to realize circumferential detonation collision.Experiments and simulations were carried out to investigate the effect of the self-designed regulator on a cylindrical shell.The detonation wave interaction generated local high-pressure regions on the shell,and the shell fractured into large and regular fragments.

    2.Due to the circumferential detonation collision,two types of fragments were formed and the macrostructure and micro characteristics of typical fragments obtained.It is found that the interaction of dual waves and the first contact position of the C-J wave played a key role in the fracture mechanism of the shell.In the process of shell expansion,these two positions preferentially initiate embryos and evolve into cracks,eventually leading to shell fracture at these positions.

    3.The effects of outer layer thickness and groove number on the fracture behavior of the shell were analyzed by numerical simulations.When the thickness of the outer layer was 5-18 mm,the fracture intensity of the shell was very small,which was compensated by a large circumferential strain.When the outer layer thickness was 10 mm,the number of grooves on the regulator increased in the range of 3-6 and the fragment width became shorter with the increase of the fracture position.However,the groove number was not the main reason for the change of the fragment width.The increase of the groove number reduced the possibility of fracture at the first contact position of the C-J wave.When the groove number reached 7,the fracture only occurred at the position of detonation collision and the fragment width rebounded.

    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 work presented in this paper has been funded by the National Natural Science Foundation of China No.11972018,and the Defense Pre-Research Joint Foundation of Chinese Ordnance Industry No.6141B012858.

    欧美日本中文国产一区发布| 日本黄色视频三级网站网址| 一级毛片女人18水好多| 色尼玛亚洲综合影院| 国产成人av激情在线播放| 脱女人内裤的视频| 精品久久久精品久久久| 欧美人与性动交α欧美软件| 精品少妇一区二区三区视频日本电影| 超碰97精品在线观看| 午夜免费成人在线视频| 韩国精品一区二区三区| 少妇粗大呻吟视频| 久久欧美精品欧美久久欧美| 国产成人系列免费观看| 亚洲成人精品中文字幕电影 | 亚洲国产看品久久| 久久精品国产清高在天天线| 天堂√8在线中文| 欧美激情高清一区二区三区| 亚洲五月天丁香| 国产熟女午夜一区二区三区| 中出人妻视频一区二区| 丰满迷人的少妇在线观看| 久久 成人 亚洲| 国产无遮挡羞羞视频在线观看| 另类亚洲欧美激情| 久久性视频一级片| 久久青草综合色| 亚洲 欧美 日韩 在线 免费| 91精品三级在线观看| a在线观看视频网站| 多毛熟女@视频| 欧美日韩瑟瑟在线播放| 久久狼人影院| 日韩人妻精品一区2区三区| 少妇的丰满在线观看| 咕卡用的链子| 亚洲精品一区av在线观看| 国产三级黄色录像| 两人在一起打扑克的视频| 热re99久久国产66热| 国产欧美日韩一区二区三区在线| 国产一区二区激情短视频| 欧美色视频一区免费| 国产在线精品亚洲第一网站| 纯流量卡能插随身wifi吗| a级片在线免费高清观看视频| 啦啦啦免费观看视频1| 久久久久久久久久久久大奶| 香蕉丝袜av| 丝袜美腿诱惑在线| 欧美黄色片欧美黄色片| 国产一区二区三区综合在线观看| 99精国产麻豆久久婷婷| 国产精品久久视频播放| 国产精品美女特级片免费视频播放器 | 黑人欧美特级aaaaaa片| 一区二区三区国产精品乱码| 美女国产高潮福利片在线看| 丝袜在线中文字幕| 成人三级黄色视频| 国产精品国产av在线观看| 免费高清在线观看日韩| 桃红色精品国产亚洲av| 精品久久蜜臀av无| 黑人欧美特级aaaaaa片| 欧美日韩国产mv在线观看视频| 三上悠亚av全集在线观看| 国产三级黄色录像| 欧美一级毛片孕妇| 精品少妇一区二区三区视频日本电影| 黄色视频不卡| 97超级碰碰碰精品色视频在线观看| 午夜日韩欧美国产| 久久久国产欧美日韩av| 久久久久久久久久久久大奶| 国产日韩一区二区三区精品不卡| 久久人人97超碰香蕉20202| 国产一区二区激情短视频| 欧美日韩视频精品一区| 亚洲成国产人片在线观看| 精品欧美一区二区三区在线| 法律面前人人平等表现在哪些方面| www国产在线视频色| 国产av一区在线观看免费| 中文字幕av电影在线播放| 女性被躁到高潮视频| 色哟哟哟哟哟哟| 国产真人三级小视频在线观看| 在线观看免费视频网站a站| 亚洲国产中文字幕在线视频| 成人18禁高潮啪啪吃奶动态图| 最新美女视频免费是黄的| 50天的宝宝边吃奶边哭怎么回事| 性欧美人与动物交配| 久久天堂一区二区三区四区| 亚洲精品中文字幕在线视频| 91在线观看av| 亚洲一区中文字幕在线| 国产99白浆流出| 亚洲自拍偷在线| 欧美黑人欧美精品刺激| 国产亚洲精品久久久久5区| 国产熟女午夜一区二区三区| 精品久久久久久成人av| 91精品国产国语对白视频| av免费在线观看网站| 韩国精品一区二区三区| 99热只有精品国产| 欧美日韩亚洲综合一区二区三区_| 国产精品偷伦视频观看了| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕av电影在线播放| 久久香蕉精品热| 国产成人系列免费观看| 成人手机av| 亚洲精品av麻豆狂野| 嫩草影视91久久| 极品教师在线免费播放| 国产片内射在线| 黄网站色视频无遮挡免费观看| 亚洲精品美女久久av网站| 人人妻人人澡人人看| 中文字幕av电影在线播放| 国产男靠女视频免费网站| 高清黄色对白视频在线免费看| 国产成人啪精品午夜网站| 日韩中文字幕欧美一区二区| av有码第一页| 久久精品国产99精品国产亚洲性色 | 超色免费av| 亚洲精品久久午夜乱码| 精品人妻1区二区| av在线天堂中文字幕 | 免费av中文字幕在线| 性色av乱码一区二区三区2| 久久精品亚洲av国产电影网| 午夜老司机福利片| 亚洲欧美精品综合久久99| 久久影院123| 久久久久久久久中文| 亚洲专区国产一区二区| 亚洲精品在线观看二区| 99国产精品免费福利视频| 免费看十八禁软件| 中出人妻视频一区二区| 亚洲av片天天在线观看| 久久精品91蜜桃| 在线观看一区二区三区| 亚洲人成77777在线视频| 成人手机av| 国产片内射在线| 女性被躁到高潮视频| 亚洲精品在线观看二区| 视频区欧美日本亚洲| 真人一进一出gif抽搐免费| www.熟女人妻精品国产| 99久久久亚洲精品蜜臀av| 精品国产乱码久久久久久男人| 亚洲avbb在线观看| 日韩精品青青久久久久久| 日日干狠狠操夜夜爽| xxx96com| 新久久久久国产一级毛片| 亚洲aⅴ乱码一区二区在线播放 | 国产成人精品久久二区二区免费| 欧美国产精品va在线观看不卡| 一二三四在线观看免费中文在| 亚洲一区二区三区色噜噜 | 母亲3免费完整高清在线观看| e午夜精品久久久久久久| 欧美最黄视频在线播放免费 | 亚洲成av片中文字幕在线观看| 亚洲国产欧美一区二区综合| 老司机午夜福利在线观看视频| 女人爽到高潮嗷嗷叫在线视频| 大型黄色视频在线免费观看| 天堂√8在线中文| 桃红色精品国产亚洲av| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产精品合色在线| 国产一区在线观看成人免费| 国产亚洲精品久久久久5区| 久久久国产一区二区| 色婷婷久久久亚洲欧美| 亚洲av美国av| 精品久久久久久成人av| 日本五十路高清| 一级片免费观看大全| 国产激情欧美一区二区| 搡老岳熟女国产| 国产精品 国内视频| 欧美丝袜亚洲另类 | av在线播放免费不卡| 好看av亚洲va欧美ⅴa在| 亚洲精品成人av观看孕妇| 午夜a级毛片| 国产成人精品无人区| 9色porny在线观看| 老汉色∧v一级毛片| 亚洲av电影在线进入| 免费在线观看黄色视频的| 亚洲国产欧美网| 999精品在线视频| 亚洲精品美女久久av网站| 如日韩欧美国产精品一区二区三区| 精品一区二区三区av网在线观看| 精品一区二区三区视频在线观看免费 | 久久久久久人人人人人| 国产成人精品久久二区二区91| 日日干狠狠操夜夜爽| 夜夜夜夜夜久久久久| 日日夜夜操网爽| 久久久久亚洲av毛片大全| 亚洲av第一区精品v没综合| 成熟少妇高潮喷水视频| 亚洲专区国产一区二区| 亚洲精品美女久久久久99蜜臀| 国产成人一区二区三区免费视频网站| 脱女人内裤的视频| 久久久久亚洲av毛片大全| 亚洲av第一区精品v没综合| 十分钟在线观看高清视频www| 99精品在免费线老司机午夜| 亚洲精品一卡2卡三卡4卡5卡| 日本欧美视频一区| 国产精品九九99| 欧美人与性动交α欧美精品济南到| 脱女人内裤的视频| 国产三级在线视频| 国产一区二区激情短视频| 久久久国产成人免费| 久久午夜综合久久蜜桃| 欧美不卡视频在线免费观看 | 91麻豆精品激情在线观看国产 | 亚洲欧美激情在线| 18禁观看日本| 少妇裸体淫交视频免费看高清 | 国产91精品成人一区二区三区| 欧美在线一区亚洲| 老司机福利观看| 两性午夜刺激爽爽歪歪视频在线观看 | 少妇 在线观看| 久久亚洲精品不卡| 日韩av在线大香蕉| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲精品久久久久久毛片| 琪琪午夜伦伦电影理论片6080| 97超级碰碰碰精品色视频在线观看| 97人妻天天添夜夜摸| 久久久久亚洲av毛片大全| 国内毛片毛片毛片毛片毛片| 女人精品久久久久毛片| 老司机在亚洲福利影院| 在线观看一区二区三区| 午夜福利在线观看吧| 午夜成年电影在线免费观看| 亚洲三区欧美一区| 天堂√8在线中文| 操美女的视频在线观看| 老司机在亚洲福利影院| 99精品欧美一区二区三区四区| 欧美日韩福利视频一区二区| netflix在线观看网站| 亚洲美女黄片视频| 久久人妻熟女aⅴ| 国产单亲对白刺激| 欧美精品亚洲一区二区| 1024香蕉在线观看| 久久午夜亚洲精品久久| 美女福利国产在线| 久久人人精品亚洲av| 亚洲精品中文字幕一二三四区| 成人18禁在线播放| 亚洲一区二区三区色噜噜 | 中文字幕人妻熟女乱码| 成人精品一区二区免费| 亚洲国产精品sss在线观看 | 亚洲欧洲精品一区二区精品久久久| 精品人妻1区二区| 麻豆国产av国片精品| 国产精品98久久久久久宅男小说| 少妇裸体淫交视频免费看高清 | 亚洲久久久国产精品| 多毛熟女@视频| 日韩成人在线观看一区二区三区| 国产精华一区二区三区| 丝袜美足系列| 夜夜夜夜夜久久久久| 脱女人内裤的视频| 亚洲成人精品中文字幕电影 | 亚洲激情在线av| av网站免费在线观看视频| 级片在线观看| 身体一侧抽搐| 国产精品av久久久久免费| 国产精品久久电影中文字幕| 日韩中文字幕欧美一区二区| 国产午夜精品久久久久久| 欧美丝袜亚洲另类 | a级毛片在线看网站| 制服人妻中文乱码| 欧美 亚洲 国产 日韩一| 国产真人三级小视频在线观看| 欧美日韩乱码在线| 精品国产一区二区久久| 国产精品乱码一区二三区的特点 | 日韩人妻精品一区2区三区| 中文字幕最新亚洲高清| 久久久国产一区二区| 女人精品久久久久毛片| 老汉色∧v一级毛片| 国产伦人伦偷精品视频| 熟女少妇亚洲综合色aaa.| 麻豆久久精品国产亚洲av | 久久国产精品人妻蜜桃| 亚洲少妇的诱惑av| 免费看十八禁软件| 久久青草综合色| 免费少妇av软件| 中文欧美无线码| 咕卡用的链子| 亚洲片人在线观看| 在线观看免费午夜福利视频| 超碰成人久久| 免费高清视频大片| 欧美乱码精品一区二区三区| 亚洲av成人不卡在线观看播放网| 99riav亚洲国产免费| 免费日韩欧美在线观看| 欧美激情 高清一区二区三区| 久99久视频精品免费| 久久伊人香网站| 免费不卡黄色视频| 一级毛片女人18水好多| 三上悠亚av全集在线观看| 亚洲欧美一区二区三区黑人| 精品久久久久久久久久免费视频 | 两个人免费观看高清视频| 亚洲一区高清亚洲精品| 成人国产一区最新在线观看| 操出白浆在线播放| 亚洲精品国产区一区二| 国产精华一区二区三区| 黄网站色视频无遮挡免费观看| 国产黄色免费在线视频| 欧美精品亚洲一区二区| 精品福利永久在线观看| 美女国产高潮福利片在线看| 高清黄色对白视频在线免费看| 国产免费男女视频| 欧美日韩乱码在线| 国产伦一二天堂av在线观看| 欧美日韩亚洲综合一区二区三区_| 午夜91福利影院| 午夜福利在线免费观看网站| 久久久久久久久中文| 久久久久久大精品| 色综合婷婷激情| 一个人免费在线观看的高清视频| svipshipincom国产片| 欧美+亚洲+日韩+国产| 日韩 欧美 亚洲 中文字幕| 丰满的人妻完整版| 久9热在线精品视频| 无限看片的www在线观看| 两个人免费观看高清视频| 午夜精品久久久久久毛片777| 欧美 亚洲 国产 日韩一| 国产精品香港三级国产av潘金莲| a在线观看视频网站| 国产成+人综合+亚洲专区| 国产91精品成人一区二区三区| 美女 人体艺术 gogo| 午夜两性在线视频| 免费观看人在逋| 最新美女视频免费是黄的| 搡老岳熟女国产| 国产麻豆69| 亚洲午夜精品一区,二区,三区| 国产一区二区三区综合在线观看| 亚洲国产中文字幕在线视频| 悠悠久久av| 丝袜美腿诱惑在线| 伦理电影免费视频| 亚洲伊人色综图| a级毛片在线看网站| 9191精品国产免费久久| 久久精品国产清高在天天线| 最新在线观看一区二区三区| 超碰97精品在线观看| 亚洲aⅴ乱码一区二区在线播放 | 色播在线永久视频| 国产精品久久久人人做人人爽| 中文欧美无线码| 欧美人与性动交α欧美软件| 视频区欧美日本亚洲| 欧美成人免费av一区二区三区| 麻豆久久精品国产亚洲av | 国产99久久九九免费精品| 久久久久亚洲av毛片大全| 免费女性裸体啪啪无遮挡网站| 大码成人一级视频| 999久久久国产精品视频| 十八禁人妻一区二区| 亚洲一区中文字幕在线| 亚洲精品一二三| 色综合欧美亚洲国产小说| 国产高清视频在线播放一区| 日韩av在线大香蕉| 好看av亚洲va欧美ⅴa在| 19禁男女啪啪无遮挡网站| 在线播放国产精品三级| 国产真人三级小视频在线观看| 性少妇av在线| 日韩免费av在线播放| 亚洲五月婷婷丁香| 热99国产精品久久久久久7| 亚洲精品粉嫩美女一区| 久9热在线精品视频| av片东京热男人的天堂| 侵犯人妻中文字幕一二三四区| 黄色视频,在线免费观看| 午夜免费鲁丝| 久久中文看片网| 自拍欧美九色日韩亚洲蝌蚪91| 欧美一区二区精品小视频在线| av超薄肉色丝袜交足视频| 色综合婷婷激情| 亚洲国产精品sss在线观看 | 在线免费观看的www视频| 成人国语在线视频| 不卡一级毛片| 久久精品国产亚洲av高清一级| 可以在线观看毛片的网站| 成人免费观看视频高清| 国产野战对白在线观看| 操出白浆在线播放| 欧美日韩黄片免| 麻豆av在线久日| 欧美国产精品va在线观看不卡| 淫秽高清视频在线观看| 日韩大尺度精品在线看网址 | 涩涩av久久男人的天堂| 欧美日韩一级在线毛片| 女人高潮潮喷娇喘18禁视频| 狂野欧美激情性xxxx| 在线观看一区二区三区| 亚洲在线自拍视频| videosex国产| 欧美性长视频在线观看| 免费看十八禁软件| 老熟妇仑乱视频hdxx| e午夜精品久久久久久久| 久久狼人影院| 精品第一国产精品| 午夜福利免费观看在线| 国产精品乱码一区二三区的特点 | 男男h啪啪无遮挡| 好看av亚洲va欧美ⅴa在| 久久精品影院6| 91大片在线观看| av视频免费观看在线观看| 色综合欧美亚洲国产小说| 国产精品影院久久| 国产欧美日韩综合在线一区二区| 老汉色∧v一级毛片| 欧美成狂野欧美在线观看| 97碰自拍视频| 精品久久久久久成人av| 美女高潮喷水抽搐中文字幕| 在线国产一区二区在线| 久久久久久久久久久久大奶| 女人被躁到高潮嗷嗷叫费观| 女人被狂操c到高潮| 在线永久观看黄色视频| 丁香六月欧美| 日本a在线网址| 久久香蕉精品热| 黄色片一级片一级黄色片| 免费少妇av软件| 久久热在线av| 亚洲性夜色夜夜综合| 曰老女人黄片| 超色免费av| 国产精品影院久久| 免费av中文字幕在线| 成年人免费黄色播放视频| 欧美人与性动交α欧美软件| 一区二区三区精品91| 久久人人97超碰香蕉20202| 狂野欧美激情性xxxx| 午夜免费观看网址| 亚洲av片天天在线观看| 中亚洲国语对白在线视频| 国产精品综合久久久久久久免费 | 女人被狂操c到高潮| 91麻豆精品激情在线观看国产 | 免费人成视频x8x8入口观看| 亚洲一区二区三区色噜噜 | 大型黄色视频在线免费观看| 校园春色视频在线观看| 久久精品影院6| 超碰成人久久| 久久午夜亚洲精品久久| 日韩成人在线观看一区二区三区| xxx96com| 午夜亚洲福利在线播放| 夜夜看夜夜爽夜夜摸 | ponron亚洲| 人妻久久中文字幕网| 亚洲av成人av| 欧美日韩瑟瑟在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 国产xxxxx性猛交| 黄片播放在线免费| 国产伦人伦偷精品视频| 国产一区二区三区视频了| 女人被狂操c到高潮| 在线免费观看的www视频| 欧美乱码精品一区二区三区| 十八禁人妻一区二区| 精品人妻1区二区| 亚洲精品国产色婷婷电影| 一个人免费在线观看的高清视频| 侵犯人妻中文字幕一二三四区| 精品国产一区二区久久| 国产欧美日韩一区二区三| 9191精品国产免费久久| 99riav亚洲国产免费| 国产男靠女视频免费网站| 成人影院久久| 天天添夜夜摸| 日本黄色视频三级网站网址| 欧美中文日本在线观看视频| 91精品三级在线观看| 亚洲专区字幕在线| 亚洲色图 男人天堂 中文字幕| 亚洲精品在线观看二区| 最新美女视频免费是黄的| 亚洲九九香蕉| 在线观看www视频免费| 亚洲成国产人片在线观看| 男女床上黄色一级片免费看| 久久精品亚洲熟妇少妇任你| 欧美精品一区二区免费开放| 国产成人系列免费观看| 久久青草综合色| 欧美在线一区亚洲| www.自偷自拍.com| 久久午夜综合久久蜜桃| 欧美日韩乱码在线| 每晚都被弄得嗷嗷叫到高潮| 色综合欧美亚洲国产小说| 黄色视频不卡| 变态另类成人亚洲欧美熟女 | 美国免费a级毛片| 久久性视频一级片| 日本欧美视频一区| 国产99久久九九免费精品| 啪啪无遮挡十八禁网站| 天天躁夜夜躁狠狠躁躁| 天堂影院成人在线观看| 丝袜美腿诱惑在线| 亚洲精品国产区一区二| 美女扒开内裤让男人捅视频| 淫秽高清视频在线观看| 成人特级黄色片久久久久久久| 中文字幕色久视频| 91精品三级在线观看| 美女国产高潮福利片在线看| 午夜老司机福利片| 亚洲第一青青草原| 国产免费男女视频| 在线观看日韩欧美| 一区二区日韩欧美中文字幕| 欧美乱码精品一区二区三区| 男人操女人黄网站| 精品人妻在线不人妻| 国产1区2区3区精品| 99riav亚洲国产免费| 亚洲国产精品999在线| 国产精品野战在线观看 | 男人舔女人的私密视频| 免费av中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 亚洲精品一卡2卡三卡4卡5卡| 国产色视频综合| 国产乱人伦免费视频| 午夜a级毛片| 国产精品久久电影中文字幕| 欧美黑人精品巨大| 搡老乐熟女国产| 日韩三级视频一区二区三区| 亚洲自偷自拍图片 自拍| 久久天躁狠狠躁夜夜2o2o| 亚洲成人久久性| 99热国产这里只有精品6| 好男人电影高清在线观看| 夜夜夜夜夜久久久久| 亚洲片人在线观看| 精品国内亚洲2022精品成人| 久久久久久久午夜电影 | 久久亚洲精品不卡| 久久热在线av| 国产精品爽爽va在线观看网站 | 国产精品av久久久久免费| 手机成人av网站| 久久精品亚洲av国产电影网| 日韩免费av在线播放| 国产精品永久免费网站| 少妇裸体淫交视频免费看高清 | 欧美精品一区二区免费开放| 欧美黑人精品巨大|