During the 10th Chinese Arctic scientific expedition carried out in the summer of 2019,the surface current in the high-latitude areas of the Arctic Ocean was observed using a self-developed surface drifting buoy,which...During the 10th Chinese Arctic scientific expedition carried out in the summer of 2019,the surface current in the high-latitude areas of the Arctic Ocean was observed using a self-developed surface drifting buoy,which was initially deployed in the Chukchi Sea.The buoy traversed the Chukchi Sea,Chukchi Abyssal Plain,Mendeleev Ridge,Makarov Basin,and Canada Basin over a period of 632 d.After returning to the Mendeleev Ridge,it continued to drift toward the pole.Overall,the track of the buoy reflected the characteristics of the transpolar drift and Chukchi Slope Current,as well as the inertial flow,cross-ridge surface flow,and even the surface disorganized flow for some time intervals.The results showed that:(1)the transpolar drift mainly occurs in the Chukchi Abyssal Plain,Mendeleev Ridge,and western Canada Basin to the east of the ridge where sea ice concentration is high,and the average northward flow velocity in the region between 79.41°N and 86.32°N was 5.1 cm/s;(2)the average surface velocity of the Chukchi Slope Current was 13.5 cm/s,and while this current moves westward along the continental slope,it also extends northwestward across the continental slope and flows to the deep sea;and(3)when sea ice concentration was less than 50%,the inertial flow was more significant(the maximum observed inertial flow was 26 cm/s,and the radius of the inertia circle was 3.6 km).展开更多
Mesoscale eddies are common oceanic phenomena.Although many related studies have been conducted,the ability for mesoscale eddies to carry real particles remains poorly addressed.We considered the drifters as real part...Mesoscale eddies are common oceanic phenomena.Although many related studies have been conducted,the ability for mesoscale eddies to carry real particles remains poorly addressed.We considered the drifters as real particles to characterize the capability of mesoscale eddies to carry particles.Firstly,mesoscale eddies in the northwest Pacific(99°E-180°E,0°-66°N)were identified using sea level anomaly(SLA)data from 1993 to 2015.Secondly,three important parameters(the carrying days,the number of circles the drifter revolving around the eddy center,and the carrying distances)were calculated by colocalizing eddy data with drifters.Finally,statistical analysis and composite analysis were conducted,reflecting the capability of mesoscale eddies to carry particles.The mechanisms on the carrying capability of eddies were also discussed.Results show that(1)the motion of carried drifters reflects the upper limit of rotational speed of eddies that the drifters revolve around the eddy center by≤90°for one day in most cases;(2)the drifters tend to be carried for a longer time when their minimal distances to the eddy center measured with normalized distance are small;(3)there are two types of eddies(cyclonic and anticyclonic eddies)in different subregions of northwest Pacific,and each has a different carrying capability(on average,similar in the tropical ocean and Subtropical Countercurrent,cyclonic eddies tend to have stronger carrying capability in Southern Kuroshio Extension,and anticyclonic eddies tend to have stronger carrying capability in Northern Kuroshio Extension and Subarctic Gyre);(4)on average,the carried drifters tend to travel for a longer time around the normalized eddy radii ranging from 0.41 to 0.76;(5)the carrying days are related to the Rossby number of the eddy(in general when the Rossby number is smaller,the carrying days are longer).展开更多
Because of the low convergence efficiency of the typical Vicsek model,a Vicsek with static summoning points(VSSP)algorithm based on the Vicsek model considering static summoning points is proposed.Firstly,the mathemat...Because of the low convergence efficiency of the typical Vicsek model,a Vicsek with static summoning points(VSSP)algorithm based on the Vicsek model considering static summoning points is proposed.Firstly,the mathematical model of the individual movement total cost on each summoning point is established.Then the individual classification rule is designed according to the initial state of the cluster to obtain the subclusters guided by each summoning point.Finally,the summoning factor is introduced to modify the course angle updating formula of the Vicsek model.To verify the effectiveness of the proposed algorithm and study the effect of the cluster summoning factor on the convergence rate,three groups of simulation experiments under different summoning factors are designed in this paper.To verify the superiority of the VSSP algorithm,the performance of the VSSP algorithm is compared with the classic algorithm by designing the algorithm performance comparison verification experiment.The results show that the algorithm proposed in this paper has good convergence and course angle consistency.The summoning factor is the sensitive factor of cluster convergence.This algorithm can provide a reference for efficient cluster segmentation movement.展开更多
An improved least mean square/fourth direct adaptive equalizer(LMS/F-DAE)is proposed in this paper for underwater acoustic communication in the Arctic.It is able to process complex-valued baseband signals and has bett...An improved least mean square/fourth direct adaptive equalizer(LMS/F-DAE)is proposed in this paper for underwater acoustic communication in the Arctic.It is able to process complex-valued baseband signals and has better equalization performance than LMS.Considering the sparsity feature of equalizer tap coefficients,an adaptive norm(AN)is incorporated into the cost function which is utilized as a sparse regularization.The norm constraint changes adaptively according to the amplitude of each coefficient.For small-scale coefficients,the sparse constraint exists to accelerate the convergence speed.For large-scale coefficients,it disappears to ensure smaller equalization error.The performance of the proposed AN-LMS/F-DAE is verified by the experimental data from the 9th Chinese National Arctic Research Expedition.The results show that compared with the standard LMS/F-DAE,AN-LMS/F-DAE can promote the sparse level of the equalizer and achieve better performance.展开更多
Compression tests on twenty unidirectional(UD) carbon fibre reinforced plastic(CFRP) specimens are conducted, the statistics on the measured compressive strength is calculated, and the fracture surface is characterize...Compression tests on twenty unidirectional(UD) carbon fibre reinforced plastic(CFRP) specimens are conducted, the statistics on the measured compressive strength is calculated, and the fracture surface is characterized. Two types of different fracture surface are experimentally observed, and they are corresponding to very different values on the compressive strength. A finite element(FE) analysis is conducted to investigate the influence of random fibre packing on the compressive strength. And a riks method(provided in ABAQUS software) is applied in FE model to analyze fibre buckling behaviour in the vicinity of compressive failure. The FE analysis agrees well with the experimental observation on the two types of buckling modes and also the partition of compressive strength. It is clearly shown that the random fibre packing lays a significant influence on the random variability of compressive strength of CFRP.展开更多
On a narrow warship platform,the coordinated use of shipborne weapon systems may cause firepower conflicts,which seriously endangers the ship safety.Meanwhile,with directed-energy weapons mounted on ships,firepower co...On a narrow warship platform,the coordinated use of shipborne weapon systems may cause firepower conflicts,which seriously endangers the ship safety.Meanwhile,with directed-energy weapons mounted on ships,firepower conflicts between weapons become a“high probability event”.Aiming at the problem of firepower safety control,based on the research about the collision probability model of air crafts and space targets and according to the cone of fire model of conventional weapons and directed-energy weapons,this paper solved the firepower conflict probabilities between conventional weapons as well as between conventional weapons and directed-energy weapons respectively using the methods of probability theory,and established the firepower safety control model.Then the calculation of firepower conflict probability was carried out using the dimensionality reduction method based on the equivalent conversion of polar coordinates and the power series method based on Laplace transform.The simulation results revealed that the proposed model and calculation methods are effective and reliable,which can provide theoretical basis and technical support for resolution of firepower conflicts between weapons.展开更多
Resolving the contradiction between Marine economic development and Marine ecological environment protection has become an unavoidable and sharp problem.The uncontrolled use of Marine antifouling technology will bring...Resolving the contradiction between Marine economic development and Marine ecological environment protection has become an unavoidable and sharp problem.The uncontrolled use of Marine antifouling technology will bring uncontrollable and even irreversible damage to the Marine biosphere,which will lead to ecological disaster and threaten the survival of human beings.Therefore,it is an urgent task to find antifouling technology with lower environmental toxicity under the premise of considering economy.More attention should be paid to the long-term impact of mature and new technologies on the Marine ecological environment.This paper introduces the development status of antifouling technology,its influence on Marine ecological environment and puts forward the design strategy of comprehensive biological fouling prevention and control technology.展开更多
A metamaterial vibration isolator,termed as wave-insulating isolator,is proposed,which preserves enough load-bearing capability and offers ultra-low and broad bandgaps for greatly enhanced wave insulation.It consists ...A metamaterial vibration isolator,termed as wave-insulating isolator,is proposed,which preserves enough load-bearing capability and offers ultra-low and broad bandgaps for greatly enhanced wave insulation.It consists of plate-shaped metacells,whose symmetric and antisymmetric local resonant modes offer several low and broad mode bandgaps although the complete bandgap remains high and narrow.The bandgap mechanisms,vibration isolation properties,effects of key parameters,and robustness to complex conditions are clarified.As experimentally demonstrated,the wave-insulating isolator can improve the vibration insulation in the ranges of[50 Hz,180 Hz]and[260 Hz,400 Hz]by 15 dB and 25 dB,respectively,in contrast to the conventional isolator with the same first resonant frequency.展开更多
The steering vectors estimated by the existing beamforming methods generally depart from their ground-truth values when they are inffuenced by both steering deviation and random error. To solve this problem, a steerin...The steering vectors estimated by the existing beamforming methods generally depart from their ground-truth values when they are inffuenced by both steering deviation and random error. To solve this problem, a steering vector random error modification robust beamforming method is proposed. The steering vector is divided into both inside and outside parts of the pre-defined sector feature subspace, and then the robust beamforming for steering deviation is employed to estimate the component of the steering vector inside the pre-defined sector feature subspace. After that, the uncertainty set constraint is further added to estimate the component of the steering vector outside the pre-defined observation subspace, which then brings high accuracy in steering vector estimation. Simulation results indicate that the proposed method possesses a more stable signal-to-interference-plus-noise ratio and higher steering vector estimation accuracy than the existing robust beamforming methods under both steering deviation and random error. Additionally, the peak of the steering vector beam response is an unbiased estimation of the target direction. Furthermore, sea–trial results illustrate the outstanding robustness of the proposed method in low signal-to-noise-ratio conditions.展开更多
As a part of underwater acoustics, the study and development of Arctic Acoustics is later than underwater acoustics about several hundred years. After World War II, the study of Arctic Acoustics attracts many interest...As a part of underwater acoustics, the study and development of Arctic Acoustics is later than underwater acoustics about several hundred years. After World War II, the study of Arctic Acoustics attracts many interests by some developed countries, especially U.S. and pre-Soviet Union. The research works obviously have some kinds of cold war brand. After the Cold War, with the gradual warming trend of Arctic area, the ocean & acoustic environment in Arctic and its adjacent area have been considerably concerned. The 8 countries of Arctic rim organized exclusive "Arctic Council" in 1996. The white paper of "China Arctic Policy"was published in January 26, 2018, the Chinese government declared that China is a close Arctic country, and is the responsible stakeholder of Arctic interests. The new advances in Arctic research works axe introduced in this paper, including the results of scientific survey and studies about Arctic underwater acoustics of Chinese researchers. It is showed that the Arctic underwater acoustical research area is not only parallely copy the topics what traditional underwater acoustics covered, e.g. environment noise,reverberation,and propagation, but also the topics which is specifically based on the Arctic environment. Some of these research fields cannot be included in the traditional shallow water, deep water acoustics, e.g. the ice-water interface feature, ice covered semi-acoustic channel effect, the communication and navigation of UUV in the condition of under ice, and the adaptation of sonar technique, equipment in the ice-covered environment, etc.展开更多
Objective:Fragment injury is a type of blast injury that is becoming more and more common in military campaigns and terrorist attacks.Numerical simulation methods investigating the formation of natural fragments and i...Objective:Fragment injury is a type of blast injury that is becoming more and more common in military campaigns and terrorist attacks.Numerical simulation methods investigating the formation of natural fragments and injuries to biological targets are expected to be developed.Methods:A cylindrical warhead model was established and the formation process of natural fragments was simulated using the approach of tied nodes with failure through the explicit finite element(FE)software of LS-DYNA.The interaction between the detonation product and the warhead shell was simulated using the fluidestructure interaction algorithm.A method to simulate the injury of natural fragments to a biological target was presented by transforming Lagrange elements into smooth particle hydrodynamics(SPH)particles after the natural fragments were successfully formed.A computational model of the human thorax was established to simulate the injury induced by natural fragments by the node-to-surface contact algorithm with erosion.Results:The discontinuous velocities of the warhead shell at different locations resulted in the formation of natural fragments with different sizes.The velocities of natural fragments increased rapidly at the initial stage and slowly after the warhead shell fractured.The initial velocities of natural fragments at the central part of the warhead shell were the largest,whereas those at both ends of the warhead shell were the smallest.The natural fragments resulted in bullet holes that were of the same shape as that of the fragments but slightly larger in size than the fragments in the human thorax after they penetrated through.Stress waves propagated in the ribs and enhanced the injury to soft tissues;additionally,ballistic pressure waves ahead of the natural fragments were also an injury factor to the soft tissues.Conclusion:The proposed method is effective in simulating the formation of natural fragments and their injury to biological targets.Moreover,this method will be beneficial for simulating the combined injuries of natural fragments and shock waves to biological targets.展开更多
基金The Fundamental Research Fund Project of the First Institute of OceanographyMinistry of Natural Resources+1 种基金under contract No.GY022Y07the National Natural Science Foundation of China under contract No.42106232。
文摘During the 10th Chinese Arctic scientific expedition carried out in the summer of 2019,the surface current in the high-latitude areas of the Arctic Ocean was observed using a self-developed surface drifting buoy,which was initially deployed in the Chukchi Sea.The buoy traversed the Chukchi Sea,Chukchi Abyssal Plain,Mendeleev Ridge,Makarov Basin,and Canada Basin over a period of 632 d.After returning to the Mendeleev Ridge,it continued to drift toward the pole.Overall,the track of the buoy reflected the characteristics of the transpolar drift and Chukchi Slope Current,as well as the inertial flow,cross-ridge surface flow,and even the surface disorganized flow for some time intervals.The results showed that:(1)the transpolar drift mainly occurs in the Chukchi Abyssal Plain,Mendeleev Ridge,and western Canada Basin to the east of the ridge where sea ice concentration is high,and the average northward flow velocity in the region between 79.41°N and 86.32°N was 5.1 cm/s;(2)the average surface velocity of the Chukchi Slope Current was 13.5 cm/s,and while this current moves westward along the continental slope,it also extends northwestward across the continental slope and flows to the deep sea;and(3)when sea ice concentration was less than 50%,the inertial flow was more significant(the maximum observed inertial flow was 26 cm/s,and the radius of the inertia circle was 3.6 km).
基金Supported by the National Key R&D Program of China(No.2018YFC1406202)the National Natural Science Foundation of China(No.41976188)。
文摘Mesoscale eddies are common oceanic phenomena.Although many related studies have been conducted,the ability for mesoscale eddies to carry real particles remains poorly addressed.We considered the drifters as real particles to characterize the capability of mesoscale eddies to carry particles.Firstly,mesoscale eddies in the northwest Pacific(99°E-180°E,0°-66°N)were identified using sea level anomaly(SLA)data from 1993 to 2015.Secondly,three important parameters(the carrying days,the number of circles the drifter revolving around the eddy center,and the carrying distances)were calculated by colocalizing eddy data with drifters.Finally,statistical analysis and composite analysis were conducted,reflecting the capability of mesoscale eddies to carry particles.The mechanisms on the carrying capability of eddies were also discussed.Results show that(1)the motion of carried drifters reflects the upper limit of rotational speed of eddies that the drifters revolve around the eddy center by≤90°for one day in most cases;(2)the drifters tend to be carried for a longer time when their minimal distances to the eddy center measured with normalized distance are small;(3)there are two types of eddies(cyclonic and anticyclonic eddies)in different subregions of northwest Pacific,and each has a different carrying capability(on average,similar in the tropical ocean and Subtropical Countercurrent,cyclonic eddies tend to have stronger carrying capability in Southern Kuroshio Extension,and anticyclonic eddies tend to have stronger carrying capability in Northern Kuroshio Extension and Subarctic Gyre);(4)on average,the carried drifters tend to travel for a longer time around the normalized eddy radii ranging from 0.41 to 0.76;(5)the carrying days are related to the Rossby number of the eddy(in general when the Rossby number is smaller,the carrying days are longer).
基金supported by the National Natural Science Foundation of China(51979193)the China Scholarship Council(201506290080)+1 种基金the China Postdoctoral Science Foundation(2019M653652)the Natural Science Basic Research Plan in Shaanxi Province of China(2019JQ-607).
文摘Because of the low convergence efficiency of the typical Vicsek model,a Vicsek with static summoning points(VSSP)algorithm based on the Vicsek model considering static summoning points is proposed.Firstly,the mathematical model of the individual movement total cost on each summoning point is established.Then the individual classification rule is designed according to the initial state of the cluster to obtain the subclusters guided by each summoning point.Finally,the summoning factor is introduced to modify the course angle updating formula of the Vicsek model.To verify the effectiveness of the proposed algorithm and study the effect of the cluster summoning factor on the convergence rate,three groups of simulation experiments under different summoning factors are designed in this paper.To verify the superiority of the VSSP algorithm,the performance of the VSSP algorithm is compared with the classic algorithm by designing the algorithm performance comparison verification experiment.The results show that the algorithm proposed in this paper has good convergence and course angle consistency.The summoning factor is the sensitive factor of cluster convergence.This algorithm can provide a reference for efficient cluster segmentation movement.
基金The National Natural Science Foundation of China under contract Nos 61631008 and 61901136the National Key Research and Development Program of China under contract No.2018YFC1405904+3 种基金the Fok Ying-Tong Education Foundation under contract No.151007the Heilongjiang Province Outstanding Youth Science Fund under contract No.JC2017017the Opening Funding of Science and Technology on Sonar Laboratory under contract No.6142109KF201802the Innovation Special Zone of National Defense Science and Technology.
文摘An improved least mean square/fourth direct adaptive equalizer(LMS/F-DAE)is proposed in this paper for underwater acoustic communication in the Arctic.It is able to process complex-valued baseband signals and has better equalization performance than LMS.Considering the sparsity feature of equalizer tap coefficients,an adaptive norm(AN)is incorporated into the cost function which is utilized as a sparse regularization.The norm constraint changes adaptively according to the amplitude of each coefficient.For small-scale coefficients,the sparse constraint exists to accelerate the convergence speed.For large-scale coefficients,it disappears to ensure smaller equalization error.The performance of the proposed AN-LMS/F-DAE is verified by the experimental data from the 9th Chinese National Arctic Research Expedition.The results show that compared with the standard LMS/F-DAE,AN-LMS/F-DAE can promote the sparse level of the equalizer and achieve better performance.
文摘Compression tests on twenty unidirectional(UD) carbon fibre reinforced plastic(CFRP) specimens are conducted, the statistics on the measured compressive strength is calculated, and the fracture surface is characterized. Two types of different fracture surface are experimentally observed, and they are corresponding to very different values on the compressive strength. A finite element(FE) analysis is conducted to investigate the influence of random fibre packing on the compressive strength. And a riks method(provided in ABAQUS software) is applied in FE model to analyze fibre buckling behaviour in the vicinity of compressive failure. The FE analysis agrees well with the experimental observation on the two types of buckling modes and also the partition of compressive strength. It is clearly shown that the random fibre packing lays a significant influence on the random variability of compressive strength of CFRP.
文摘On a narrow warship platform,the coordinated use of shipborne weapon systems may cause firepower conflicts,which seriously endangers the ship safety.Meanwhile,with directed-energy weapons mounted on ships,firepower conflicts between weapons become a“high probability event”.Aiming at the problem of firepower safety control,based on the research about the collision probability model of air crafts and space targets and according to the cone of fire model of conventional weapons and directed-energy weapons,this paper solved the firepower conflict probabilities between conventional weapons as well as between conventional weapons and directed-energy weapons respectively using the methods of probability theory,and established the firepower safety control model.Then the calculation of firepower conflict probability was carried out using the dimensionality reduction method based on the equivalent conversion of polar coordinates and the power series method based on Laplace transform.The simulation results revealed that the proposed model and calculation methods are effective and reliable,which can provide theoretical basis and technical support for resolution of firepower conflicts between weapons.
文摘Resolving the contradiction between Marine economic development and Marine ecological environment protection has become an unavoidable and sharp problem.The uncontrolled use of Marine antifouling technology will bring uncontrollable and even irreversible damage to the Marine biosphere,which will lead to ecological disaster and threaten the survival of human beings.Therefore,it is an urgent task to find antifouling technology with lower environmental toxicity under the premise of considering economy.More attention should be paid to the long-term impact of mature and new technologies on the Marine ecological environment.This paper introduces the development status of antifouling technology,its influence on Marine ecological environment and puts forward the design strategy of comprehensive biological fouling prevention and control technology.
基金supported by the National Natural Science Foundation of China(Nos.52241103 and 52322505)the Natural Science Fund for Distinguished Young Scholars of Hunan Province of China(No.2023JJ10055)。
文摘A metamaterial vibration isolator,termed as wave-insulating isolator,is proposed,which preserves enough load-bearing capability and offers ultra-low and broad bandgaps for greatly enhanced wave insulation.It consists of plate-shaped metacells,whose symmetric and antisymmetric local resonant modes offer several low and broad mode bandgaps although the complete bandgap remains high and narrow.The bandgap mechanisms,vibration isolation properties,effects of key parameters,and robustness to complex conditions are clarified.As experimentally demonstrated,the wave-insulating isolator can improve the vibration insulation in the ranges of[50 Hz,180 Hz]and[260 Hz,400 Hz]by 15 dB and 25 dB,respectively,in contrast to the conventional isolator with the same first resonant frequency.
文摘The steering vectors estimated by the existing beamforming methods generally depart from their ground-truth values when they are inffuenced by both steering deviation and random error. To solve this problem, a steering vector random error modification robust beamforming method is proposed. The steering vector is divided into both inside and outside parts of the pre-defined sector feature subspace, and then the robust beamforming for steering deviation is employed to estimate the component of the steering vector inside the pre-defined sector feature subspace. After that, the uncertainty set constraint is further added to estimate the component of the steering vector outside the pre-defined observation subspace, which then brings high accuracy in steering vector estimation. Simulation results indicate that the proposed method possesses a more stable signal-to-interference-plus-noise ratio and higher steering vector estimation accuracy than the existing robust beamforming methods under both steering deviation and random error. Additionally, the peak of the steering vector beam response is an unbiased estimation of the target direction. Furthermore, sea–trial results illustrate the outstanding robustness of the proposed method in low signal-to-noise-ratio conditions.
文摘As a part of underwater acoustics, the study and development of Arctic Acoustics is later than underwater acoustics about several hundred years. After World War II, the study of Arctic Acoustics attracts many interests by some developed countries, especially U.S. and pre-Soviet Union. The research works obviously have some kinds of cold war brand. After the Cold War, with the gradual warming trend of Arctic area, the ocean & acoustic environment in Arctic and its adjacent area have been considerably concerned. The 8 countries of Arctic rim organized exclusive "Arctic Council" in 1996. The white paper of "China Arctic Policy"was published in January 26, 2018, the Chinese government declared that China is a close Arctic country, and is the responsible stakeholder of Arctic interests. The new advances in Arctic research works axe introduced in this paper, including the results of scientific survey and studies about Arctic underwater acoustics of Chinese researchers. It is showed that the Arctic underwater acoustical research area is not only parallely copy the topics what traditional underwater acoustics covered, e.g. environment noise,reverberation,and propagation, but also the topics which is specifically based on the Arctic environment. Some of these research fields cannot be included in the traditional shallow water, deep water acoustics, e.g. the ice-water interface feature, ice covered semi-acoustic channel effect, the communication and navigation of UUV in the condition of under ice, and the adaptation of sonar technique, equipment in the ice-covered environment, etc.
基金The work was funded by the National Science Foundation for Young Scientists of China(11902356)China Postdoctoral Science Foundation(2018M633715)+1 种基金Innovation and Cultivation Fund of the Sixth Medical Center of PLA General Hospital(No.CXPY201825)the Army Scientific Research(LB20182D040012).
文摘Objective:Fragment injury is a type of blast injury that is becoming more and more common in military campaigns and terrorist attacks.Numerical simulation methods investigating the formation of natural fragments and injuries to biological targets are expected to be developed.Methods:A cylindrical warhead model was established and the formation process of natural fragments was simulated using the approach of tied nodes with failure through the explicit finite element(FE)software of LS-DYNA.The interaction between the detonation product and the warhead shell was simulated using the fluidestructure interaction algorithm.A method to simulate the injury of natural fragments to a biological target was presented by transforming Lagrange elements into smooth particle hydrodynamics(SPH)particles after the natural fragments were successfully formed.A computational model of the human thorax was established to simulate the injury induced by natural fragments by the node-to-surface contact algorithm with erosion.Results:The discontinuous velocities of the warhead shell at different locations resulted in the formation of natural fragments with different sizes.The velocities of natural fragments increased rapidly at the initial stage and slowly after the warhead shell fractured.The initial velocities of natural fragments at the central part of the warhead shell were the largest,whereas those at both ends of the warhead shell were the smallest.The natural fragments resulted in bullet holes that were of the same shape as that of the fragments but slightly larger in size than the fragments in the human thorax after they penetrated through.Stress waves propagated in the ribs and enhanced the injury to soft tissues;additionally,ballistic pressure waves ahead of the natural fragments were also an injury factor to the soft tissues.Conclusion:The proposed method is effective in simulating the formation of natural fragments and their injury to biological targets.Moreover,this method will be beneficial for simulating the combined injuries of natural fragments and shock waves to biological targets.