CN108595848A - A kind of Penetration is hit by a bullet target action model modeling method and device - Google Patents

A kind of Penetration is hit by a bullet target action model modeling method and device Download PDF

Info

Publication number
CN108595848A
CN108595848A CN201810391492.XA CN201810391492A CN108595848A CN 108595848 A CN108595848 A CN 108595848A CN 201810391492 A CN201810391492 A CN 201810391492A CN 108595848 A CN108595848 A CN 108595848A
Authority
CN
China
Prior art keywords
warhead
axial vibration
rank
model
vibration
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810391492.XA
Other languages
Chinese (zh)
Other versions
CN108595848B (en
Inventor
程祥利
赵慧
李林川
叶海福
刘波
刘军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Electronic Engineering of CAEP
Original Assignee
Institute of Electronic Engineering of CAEP
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Electronic Engineering of CAEP filed Critical Institute of Electronic Engineering of CAEP
Priority to CN201810391492.XA priority Critical patent/CN108595848B/en
Publication of CN108595848A publication Critical patent/CN108595848A/en
Application granted granted Critical
Publication of CN108595848B publication Critical patent/CN108595848B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/18Manufacturability analysis or optimisation for manufacturability

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The embodiment of the present invention provides a kind of Penetration and is hit by a bullet target action model modeling method and device.The method includes:Obtain the parameter of warhead;The mass spring damper model of warhead is obtained based on axial vibration frequency, warhead quality and each rank axial vibration Damping calculating;The transmission function of warhead axial vibration is built according to mass spring damper model, and is calculated according to transmission function and obtained elastic vibration model;The rigid motion model of warhead is built, and is constituted with elastic vibration models coupling and plays target action model.Described device is for executing the method, the embodiment of the present invention is by building elastic vibration model, elastic vibration model and rigid motion model are combined composition and play target action model, influence of the warhead elastic vibration to warhead internal component is analyzed more fully hereinafter, to improve the accuracy for the stressing conditions for playing target action model reflection warhead internal component.

Description

A kind of Penetration is hit by a bullet target action model modeling method and device
Technical field
The present invention relates to data modeling technical field, built in particular to a kind of Penetration target action model of being hit by a bullet Mould method and device.
Background technology
In modern war, increasingly rely on ultrahigh speed Penetrating Warhead to fortification, surface ship, airfield runway, The high intensity hard goals such as face building, aircraft masking library are injured.And in Penetration, warhead (bullet) and high-strength rigid target The bullet target mechanism of action between (target) is extremely complex, is brought to the anti high overload optimization design of warhead internal component very big tired Difficulty, as mechanics inputs, indefinite, optimization direction is indefinite.
Currently, domestic and foreign scholars by establish rationally, effectively, simplify bullet target action model come disclose warhead with it is high-strength The bullet target mechanism of action between hard goal is spent, the anti high overload optimization design of warhead internal component is and guided, such as fuse, dress Medicine.
Existing bullet target action model modeling method only considers warhead rigid motion, i.e., warhead is not sent out in Penetration Raw any deformation.But violent flexible deformation can occur when warhead penetration high-strength rigid target, or even have to a certain degree Plastic deformation, directly by warhead be assumed to be rigid body cannot comprehensively, reasonably describe play target mechanism.In addition, weapon is electric Sub- system is more fragile under the conditions of sine impacts, and also illustrates that warhead elastic vibration is to influence the pass of internal component survivability Key factor.Therefore, only consider that the bullet target action model of warhead rigid motion cannot comprehensively, reasonably reflect group inside warhead The stressing conditions of part also really can not effectively instruct the anti high overload optimization design of warhead internal component.
Invention content
In view of this, a kind of Penetration of being designed to provide of the embodiment of the present invention is hit by a bullet target action model modeling method And device, to solve the above technical problems.
In a first aspect, be hit by a bullet target action model modeling method an embodiment of the present invention provides a kind of Penetration, including:
The parameter of warhead is obtained, the parameter includes:Case material elasticity modulus, case material density, fight minister Degree, warhead quality and each rank axial vibration damping;
According to the acquisition of the case material elasticity modulus, the case material density and the warhead length computation The axial vibration frequency of each rank of warhead;
Described in the axial vibration frequency, the warhead quality and each rank axial vibration Damping calculating acquisition The mass spring damper model of warhead;
The transmission function of warhead axial vibration is built according to the mass spring damper model, and according to the transmission letter Number, which calculates, obtains elastic vibration model;
The rigid motion model of the warhead is built, and is constituted with the elastic vibration models coupling and plays target effect mould Type.
Further, described according to the case material elasticity modulus, the case material density and the fight minister Degree calculates the axial vibration frequency for obtaining the warhead, including:
According to formulaCalculate the axial vibration frequency for obtaining the warhead, wherein fPlay iFor the axial vibration frequency of the i-th rank of the warhead, L is the warhead length, and E is the case material elasticity modulus, ρ is the case material density.
Further, the parameter further includes the dimensionless damping ratio of each rank axial vibration, correspondingly, described based on described Axial vibration frequency, the warhead quality and each rank axial vibration Damping calculating obtain the Tuned mass damper of the warhead Spring model, including:
The natural angular frequency that each rank axial vibration is calculated according to the axial vibration frequency, according to the warhead quality and The natural angular frequency calculates the rigidity of each rank axial vibration, according to the rigidity of each rank axial vibration, the fight Portion's quality and the dimensionless damping ratio calculate each rank axial vibration damping;
It is obtained according to the Rigidity Calculation of the warhead quality, each rank axial vibration damping and each rank axial vibration Obtain the mass spring damper model of the warhead.
Further, the natural angular frequency that each rank axial vibration is calculated according to the axial vibration frequency, including:
According to ωPlay i=2 π fPlay iCalculate the natural angular frequency of each rank axial vibration, wherein fPlay iFor the i-th rank of the warhead Axial vibration frequency;
The rigidity that each rank axial vibration is calculated according to the warhead quality and the natural angular frequency, including:
According toCalculate the rigidity of each rank axial vibration, wherein MBulletFor the warhead quality;
The rigidity, warhead quality and the dimensionless damping ratio meter according to each rank axial vibration Each rank axial vibration damping is calculated, including:
According toCalculate each rank axial vibration damping, wherein ξPlay iIt is damped for the dimensionless Than.
Further, the transmission function is:Wherein, Δ X (s) For the axial deflection, Fx(s) it is to play target active force, GPlay i(s) it is the transmission function of the i-th rank elastic vibration, andMBulletFor the warhead quality, s is multiple The independent variable of number field, CPlay iFor the damping of the i-th rank axial vibration, KPlay iFor the rigidity of the i-th rank axial vibration, ξPlay iIt is immeasurable for the i-th rank Guiding principle damping ratio.
Correspondingly, described calculated according to the transmission function obtains elastic vibration model, including:
The transmission function is subjected to partial differential calculating, the elastic vibration model of acquisition is:
Further, the rigid motion model of the structure warhead, including:
The bullet target active force in Penetration is calculated according to cavity-expansion theory and differential panel method;
Based on rigid motion theory, the rigid motion model is established using Newton's second law.
Further, described obtained according to the rigid motion model and the elastic vibration model plays target action model, Including:
The rigid motion model and the elastic vibration model are subjected to simultaneous, constitute the bullet target action model.
Second aspect is hit by a bullet target action model model building device an embodiment of the present invention provides a kind of Penetration, including:
Acquisition module, the parameter for obtaining warhead, the parameter include:Case material elasticity modulus, case material Density, warhead length, warhead quality and each rank axial vibration damping;
First computing module, for according to the case material elasticity modulus, the case material density and the fight Portion's length computation obtains the axial vibration frequency of each rank of the warhead;
Second computing module axially shakes for being based on the axial vibration frequency, the warhead quality and each rank Dynamic damping calculates the mass spring damper model for obtaining the warhead;
Third computing module, the transmission letter for building warhead axial vibration according to the mass spring damper model Number, and calculated according to the transmission function and obtain elastic vibration model;
Build module, the rigid motion model for building the warhead, and with the elastic vibration models coupling structure At bullet target action model.
Further, first computing module, is specifically used for:
According to formulaCalculate the axial vibration frequency for obtaining the warhead, wherein fPlay iFor the axial vibration frequency of the i-th rank of the warhead, L is the warhead length, and E is the case material elasticity modulus, ρ is the case material density.
Further, the parameter further includes the dimensionless damping ratio of each rank axial vibration, correspondingly, described second calculates Module is specifically used for:
The natural angular frequency that each rank axial vibration is calculated according to the axial vibration frequency, according to the warhead quality and The natural angular frequency calculates the rigidity of each rank axial vibration, according to the rigidity of each rank axial vibration, the fight Portion's quality and the dimensionless damping ratio calculate each rank axial vibration damping;
It is obtained according to the Rigidity Calculation of the warhead quality, each rank axial vibration damping and each rank axial vibration Obtain the mass spring damper model of the warhead.
A kind of Penetration provided in an embodiment of the present invention is hit by a bullet target action model modeling method and device, is passed through and is built bullet Property model of vibration, by elastic vibration model and rigid motion model be combined composition play target action model, analyze more fully hereinafter Influence of the warhead elastic vibration to warhead internal component plays group inside target action model reflection warhead to improve The accuracy of the stressing conditions of part.
Other features and advantages of the present invention will be illustrated in subsequent specification, also, partly be become from specification It is clear that by implementing understanding of the embodiment of the present invention.The purpose of the present invention and other advantages can be by saying what is write Specifically noted structure is realized and is obtained in bright book, claims and attached drawing.
Description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached Figure is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as pair The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this A little attached drawings obtain other relevant attached drawings.
Fig. 1 is that a kind of Penetration provided in an embodiment of the present invention is hit by a bullet target action model modeling method flow diagram;
Fig. 2 is mass spring damper model schematic diagram provided in an embodiment of the present invention;
Fig. 3 is the transmission function block diagram of warhead provided in an embodiment of the present invention vibration;
Fig. 4 is that a kind of Penetration provided in an embodiment of the present invention is hit by a bullet target action model model building device structural schematic diagram;
Fig. 5 is electronic devices structure block diagram provided by the embodiments of the present application.
Specific implementation mode
Below in conjunction with attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete Ground describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Usually exist The component of the embodiment of the present invention described and illustrated in attached drawing can be arranged and be designed with a variety of different configurations herein.Cause This, the detailed description of the embodiment of the present invention to providing in the accompanying drawings is not intended to limit claimed invention below Range, but it is merely representative of the selected embodiment of the present invention.Based on the embodiment of the present invention, those skilled in the art are not doing The every other embodiment obtained under the premise of going out creative work, shall fall within the protection scope of the present invention.
It should be noted that:Similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi It is defined, then it further need not be defined and explained in subsequent attached drawing in a attached drawing.Meanwhile the present invention's In description, term " first ", " second " etc. are only used for distinguishing description, are not understood to indicate or imply relative importance.
Fig. 1 is that a kind of Penetration provided in an embodiment of the present invention is hit by a bullet target action model modeling method flow diagram, As shown in Figure 1, the method, including:
Step 101:The parameter of warhead is obtained, the parameter includes:Case material elasticity modulus, case material density, Warhead length, warhead quality and each rank axial vibration damping;
In specific implementation process, since its different corresponding relevant parameter of the type of warhead is different, in structure It builds before the bullet target action model of warhead, to obtain the parameter of warhead, wherein parameter may include:Case material elasticity Modulus, case material density, warhead length, warhead quality and each rank axial vibration damping, it should be noted that warhead Parameter can also include other parameters, the embodiment of the present invention be not especially limited this.
Step 102:According to the case material elasticity modulus, the case material density and the warhead length computation Obtain the axial vibration frequency of each rank of the warhead;
In specific implementation process, when the shape of warhead is axial symmetry cylindrical bar, according to the shell material of warhead Expect that elasticity modulus, case material density and warhead length carry out model analysis, calculates the axial direction for obtaining each order of warhead Vibration frequency.
Step 103:Based on the axial vibration frequency, the warhead quality and each rank axial vibration Damping calculating Obtain the mass spring damper model of the warhead;
It, can be with base after calculating acquires the axial vibration frequency of each order of warhead in specific implementation process The mass spring damper model of warhead is established in axial vibration frequency, warhead quality and each rank axial vibration damping.It can be with Understand, warhead quality refers to warhead own wt.
Step 104:The transmission function of warhead axial vibration is built according to the mass spring damper model, and according to institute It states transmission function and calculates acquisition elastic vibration model;
In specific implementation process, after calculating acquisition mass spring damper model, Tuned mass damper spring can be utilized Model construction is using warhead active force as input, and transmission function of the axial deflection as the axial vibration of output is then right The transmission function being calculated carries out partial differential and calculates acquisition partial differential equation, wherein partial differential equation are elastic vibrating dynamic model Type.
Step 105:The rigid motion model of the warhead is built, and is constituted with the elastic vibration models coupling and plays target Action model.
In specific implementation process, the collision between warhead and high-strength rigid target is equivalent to axially to be applied to warhead Add the bullet target active force of a loading speed quickly, it is assumed that be Fx.On the one hand the bullet target active force can cause the deceleration of warhead Movement, i.e. rigid motion;On the other hand the multistage axial vibration of warhead, i.e. flexible deformation can have been excited, and has been transferred to fight Portion's internal component influences the survivability of internal component.Therefore, in the process of actual warhead high speed penetration high-strength rigid target In, existing warhead rigid motion, and have warhead elastic vibration, so, it can will play target action model and be divided into rigid motion Model and elastic vibration model two parts.Therefore, it is necessary to build the rigid motion model of warhead, finally it should be noted that rigid body is transported The structure of movable model is consistent with construction method in the prior art, i.e., the movement that Penetration is established by Newton's second law is micro- Equation group, the embodiment of the present invention is divided to repeat no more this.By the elastic vibrating dynamic model of the rigid motion model of structure and above-mentioned structure Type is combined composition and plays target action model.It is understood that rigid motion model and elastic vibration model can be joined Vertical be used as plays target action model.
Elastic vibration model and rigid motion model are combined structure by the embodiment of the present invention by building elastic vibration model At target action model is played, influence of the warhead elastic vibration to warhead internal component is analyzed more fully hereinafter, to improve Play the accuracy of the stressing conditions of target action model reflection warhead internal component.
On the basis of the above embodiments, it is described according to the case material elasticity modulus, the case material density and The warhead length computation obtains the axial vibration frequency of the warhead, including:
According to formulaCalculate the axial vibration frequency for obtaining the warhead, wherein fPlay iFor the axial vibration frequency of the i-th rank of the warhead, L is the warhead length, and E is the case material elasticity modulus, ρ is the case material density.
In specific implementation process, the inherent characteristic of warhead axial vibration can be obtained by model analysis, you can with According to formula (1) solve the axial vibration frequency for the warhead that one-dimensional wave equation obtains.Wherein, formula (1) is:
Wherein, fPlay iFor the axial vibration frequency of the i-th rank of warhead, i is positive integer, and optionally, i maximums can take 5, Other numerical value, the embodiment of the present invention can be taken to be not especially limited this, L is warhead length, and E is case material springform Amount, ρ are case material density, and therefore, the value of E and ρ are related to the material of warhead.
The embodiment of the present invention calculates the axial vibration frequency for obtaining each rank of warhead by formula, is elastic vibration model Structure provides data foundation, therefore, group inside warhead is combined as by elastic vibration model and rigid motion model The anti high overload optimization design of part provides more accurate mechanics design input and more specifically optimization design direction.
On the basis of the above embodiments, the parameter further includes the dimensionless damping ratio of each rank axial vibration, correspondingly, It is described that the fight is obtained based on the axial vibration frequency, the warhead quality and each rank axial vibration Damping calculating The mass spring damper model in portion, including:
The natural angular frequency that each rank axial vibration is calculated according to the axial vibration frequency, according to the warhead quality and The natural angular frequency calculates the rigidity of each rank axial vibration, according to the rigidity of each rank axial vibration, the fight Portion's quality and the dimensionless damping ratio calculate each rank axial vibration damping;
It is obtained according to the Rigidity Calculation of the warhead quality, each rank axial vibration damping and each rank axial vibration Obtain the mass spring damper model of the warhead.
In specific implementation process, it can be calculated according to the axial vibration frequency of each rank of warhead and obtain each rank of warhead Axial vibration natural angular frequency.After natural angular frequency is calculated, according to warhead quality and it is calculated intrinsic Angular frequency calculates the rigidity for obtaining each rank axial vibration of warhead.In addition, the corresponding parameter of warhead further comprises each rank axial direction The dimensionless damping ratio of vibration can be hindered according to the rigidity, warhead quality and dimensionless of each rank axial vibration being calculated Buddhist nun is than calculating each rank axial vibration damping.
Fig. 2 is mass spring damper model schematic diagram provided in an embodiment of the present invention, as shown in Fig. 2, according to warhead matter Measure MBullet, the corresponding each rank axial vibration of warhead damp CPlay i, each rank axial vibration stiffness KPlay iAnd axial deflection Δ X meters Calculate the mass spring damper model for the mass spring damper model composition warhead for obtaining warhead.It is understood that L is war Pan length, FxFor warhead active force.
The embodiment of the present invention is obtained by the Rigidity Calculation of warhead quality, each rank axial vibration damping and each rank axial vibration The mass spring damper model for obtaining warhead can accurately acquire the corresponding elastic vibration model of warhead, and then improve To warhead to the accuracy of the stressing conditions of internal component analysis in Penetration.
On the basis of the above embodiments, the intrinsic angle that each rank axial vibration is calculated according to the axial vibration frequency Frequency, including:
According to ωPlay i=2 π fPlay iCalculate the natural angular frequency of each rank axial vibration, wherein fPlay iFor the i-th rank of the warhead Axial vibration frequency.
It, can be according to public affairs to the calculating of the natural angular frequency of each rank axial vibration of warhead in specific implementation process Formula (2) obtains, wherein formula (2) is as follows:
ωPlay i=2 π fPlay i (2)
Wherein, ωPlay iFor the natural angular frequency of each rank axial vibration of warhead, fPlay iFor the axial vibration of the i-th rank of warhead Frequency, it is to be understood that i is positive integer, and maximum value can be set according to actual conditions.
The rigidity that each rank axial vibration is calculated according to the warhead quality and the natural angular frequency, including:
According toCalculate the rigidity of each rank axial vibration, wherein MBulletFor the warhead quality.
In specific implementation process, after the natural angular frequency that each rank axial vibration is calculated, formula (3) is substituted into In, the rigidity for obtaining each rank axial vibration of warhead can be calculated, formula (3) is as follows:
Therefore, according to warhead mass MBulletWith the natural angular frequency ω for each rank axial vibration of warhead being calculatedPlay iIt can To calculate the rigidity for obtaining each rank axial vibration.
The rigidity, warhead quality and the dimensionless damping ratio meter according to each rank axial vibration Each rank axial vibration damping is calculated, including:
According toCalculate each rank axial vibration damping, wherein ξPlay iIt is damped for the dimensionless Than.
In specific implementation process, after the rigidity for calculating each rank axial vibration for obtaining warhead, formula (4) is substituted into In, the corresponding each rank axial vibration damping of warhead can be calculated, formula (4) is as follows:
Wherein, ξPlay iFor the dimensionless damping ratio of warhead, MBulletFor warhead quality, KPlay iIt axially shakes for each rank of warhead Dynamic rigidity.
The embodiment of the present invention passes through warhead quality, the rigidity and axial direction of each rank axial vibration damping, each rank axial vibration Deformation calculation obtains the mass spring damper model of warhead, can accurately acquire the corresponding elastic vibrating dynamic model of warhead Type, and then improve the accuracy analyzed the stressing conditions of internal component in Penetration warhead.
On the basis of the above embodiments, the transmission function is: Wherein, Δ X (s) is the axial deflection, Fx(s) it is to play target active force, GPlay i(s) it is the transmission function of the i-th rank elastic vibration, AndMBulletFor the warhead quality, s is The independent variable of complex field, CPlay iFor the damping of the i-th rank axial vibration, KPlay iFor the rigidity of the i-th rank axial vibration, ξPlay iFor the i-th rank without Dimension damping ratio.
In specific implementation process, Fig. 3 is the transmission function block diagram of warhead provided in an embodiment of the present invention vibration, such as Shown in Fig. 3, from input Fx(s) shown in the transmission function such as formula (5) for arriving output Δ X (s):
Wherein, Δ X (s) is the axial deflection of warhead, Fx(s) it is warhead and high-strength rigid in Penetration Collision between target, to the bullet target active force that warhead applies, GPlay i(s) it is the transmission function of the i-th rank elastic vibration, and Optionally, the value of n can be 5,6 Or 7, it can also be other numerical value, can be set according to actual conditions, the embodiment of the present invention is not especially limited this, MBulletFor Warhead quality, s are the independent variable of complex field, CPlay iFor the damping of the i-th rank axial vibration, it can be calculated and be obtained according to formula (4), KPlay iFor the rigidity of the i-th rank axial vibration, it can be calculated and be obtained according to formula (3), ξPlay iFor the i-th rank dimensionless damping ratio.
Correspondingly, described calculated according to the transmission function obtains elastic vibration model, including:
The transmission function is subjected to partial differential calculating, the elastic vibration model of acquisition is:
In specific implementation process, formula (5) is subjected to partial differential calculating, shown in result of calculation such as formula (6):
Wherein, which is the partial differential equation of warhead axial vibration, which is elastic vibration Model.It is understood that consistent in the meaning and above-described embodiment of parameters in formula (6), details are not described herein again.
The embodiment of the present invention by from Solid Mechanics, automatically control, the angle of the Multidisciplinary Integrations such as mechanical oscillation, build Elastic vibration model has been found, and elastic vibration model and rigid motion model are combined to obtain and play target action model, it can More accurate mechanics design input is provided and is more specifically optimized for the anti high overload optimization design of warhead internal component Design direction.
On the basis of the above embodiments, the rigid motion model of the structure warhead, including:
The bullet target active force in Penetration is calculated according to cavity-expansion theory and differential panel method;
Based on rigid motion theory, the rigid motion model is established using Newton's second law.
In specific implementation process, for the structure of the rigid motion model of warhead, first, warhead penetration is established Process force analysis model, and using the bullet target active force in cavity-expansion theory and differential panel method solution Penetration;With Afterwards, it is based on rigid motion theory, terminal ballistics model, i.e. motion in Penetration are established using Newton's second law Equation group is finally input to play target active force, solves terminal ballistics model using the method for numerical integration, obtains displacement, speed Degree overloads the rule that changes with time, and and guides the anti high overload optimization design of warhead internal component.
The embodiment of the present invention establishes rigid motion model by using Newton's second law, in conjunction with the elastic vibration of structure Model constitutes the bullet target action model of warhead, can be reasonably to warhead inside in Penetration by playing target action model In force analysis.
It is described to be obtained according to the rigid motion model and the elastic vibration model on the basis of the various embodiments described above Target action model is played, including:
The rigid motion model and the elastic vibration model are subjected to simultaneous, constitute the bullet target action model.
Specifically, after building rigid motion model and elastic vibration model respectively, by rigid motion model and elasticity The equation that the corresponding equation of model of vibration carries out simultaneous composition is to play target action model.
Elastic vibration model and rigid motion model are combined structure by the embodiment of the present invention by building elastic vibration model At target action model is played, influence of the warhead elastic vibration to warhead internal component is analyzed more fully hereinafter, to improve Play the accuracy of the stressing conditions of target action model reflection warhead internal component.
Fig. 4 is that a kind of Penetration provided in an embodiment of the present invention is hit by a bullet target action model model building device structural schematic diagram, As shown in figure 4, described device, including:Acquisition module 401, the first computing module 402, the second computing module 403, third calculate Module 404 and structure module 405, wherein
Acquisition module 401 is used to obtain the parameter of warhead, and the parameter includes:Case material elasticity modulus, shell material Expect density, warhead length, warhead quality and each rank axial vibration damping;First computing module 402 is used for according to the shell The axial direction that body elasticity modulus of materials, the case material density and the warhead length computation obtain each rank of the warhead is shaken Dynamic frequency;Second computing module 403 based on the axial vibration frequency, the warhead quality and each rank for axially being shaken Dynamic damping calculates the mass spring damper model for obtaining the warhead;Third computing module 404 is used to be hindered according to the quality Buddhist nun's spring model builds the transmission function of warhead axial vibration, and is calculated according to the transmission function and obtain elastic vibrating dynamic model Type;Structure module 405 is used to build the rigid motion model of the warhead, and constitutes bullet with the elastic vibration models coupling Target action model.
In specific implementation process, before the bullet target action model of structure warhead, acquisition module 401 obtains fight The parameter in portion, wherein parameter may include:Case material elasticity modulus, case material density, warhead length, warhead matter Amount and each rank axial vibration damping, it should be noted that the parameter of warhead can also include other parameters, the embodiment of the present invention This is not especially limited.When the shape of warhead is axial symmetry cylindrical bar, the first computing module 402 is according to the shell of warhead Body elasticity modulus of materials, case material density and warhead length carry out model analysis, calculate and obtain each order of warhead Axial vibration frequency.After calculating acquires the axial vibration frequency of each order of warhead, the second computing module 403 can be with base The mass spring damper model of warhead is established in axial vibration frequency, warhead quality and each rank axial vibration damping.It can be with Understand, warhead quality refers to warhead own wt.After calculating acquisition mass spring damper model, third calculates mould Block 404 can be built using mass spring damper model using warhead active force as input, and axial deflection is as output Then the transmission function of axial vibration carries out partial differential to the transmission function being calculated and calculates acquisition partial differential equation, wherein Partial differential equation are elastic vibration model.It builds module 405 and builds rigid motion model, and by the rigid motion model of structure It is combined composition with the elastic vibration model of above-mentioned structure and plays target action model.It is understood that can be by rigid motion mould Type and elastic vibration model carry out simultaneous as bullet target action model.
Elastic vibration model and rigid motion model are combined structure by the embodiment of the present invention by building elastic vibration model At target action model is played, influence of the warhead elastic vibration to warhead internal component is analyzed more fully hereinafter, to improve Play the accuracy of the stressing conditions of target action model reflection warhead internal component.
On the basis of the above embodiments, first computing module, is specifically used for:
According to formulaCalculate the axial vibration frequency for obtaining the warhead, wherein fPlay iFor the axial vibration frequency of the i-th rank of the warhead, L is the warhead length, and E is the case material elasticity modulus, ρ is the case material density.
On the basis of the above embodiments, the parameter further includes the dimensionless damping ratio of each rank axial vibration, correspondingly, Second computing module, is specifically used for:
The natural angular frequency that each rank axial vibration is calculated according to the axial vibration frequency, according to the warhead quality and The natural angular frequency calculates the rigidity of each rank axial vibration, according to the rigidity of each rank axial vibration, the fight Portion's quality and the dimensionless damping ratio calculate each rank axial vibration damping;
According to the warhead quality, each rank axial vibration damping, the rigidity of each rank axial vibration and described Axial deflection calculates the mass spring damper model for obtaining the warhead.
On the basis of the above embodiments, the intrinsic angle that each rank axial vibration is calculated according to the axial vibration frequency Frequency, including:
According to ωPlay i=2 π fPlay iCalculate the natural angular frequency of each rank axial vibration, wherein fPlay iFor the i-th rank of the warhead Axial vibration frequency;
The rigidity that each rank axial vibration is calculated according to the warhead quality and the natural angular frequency, including:
According toCalculate the rigidity of each rank axial vibration, wherein MBulletFor the warhead quality;
The rigidity, warhead quality and the dimensionless damping ratio meter according to each rank axial vibration Each rank axial vibration damping is calculated, including:
According toCalculate each rank axial vibration damping, wherein ξPlay iIt is damped for the dimensionless Than.
On the basis of the above embodiments, the transmission function is: Wherein, Δ X (s) is the axial deflection, Fx(s) it is to play target active force, GPlay i(s) it is the transmission function of the i-th rank elastic vibration, AndMBulletFor the warhead quality, s is The independent variable of complex field, CPlay iFor the damping of the i-th rank axial vibration, KPlay iFor the rigidity of the i-th rank axial vibration, ξPlay iFor the i-th rank without Dimension damping ratio.
Correspondingly, the third computing module, is specifically used for:
The transmission function is subjected to partial differential calculating, the elastic vibration model of acquisition is:
On the basis of the above embodiments, the structure module, is specifically used for:
The bullet target active force in Penetration is calculated according to cavity-expansion theory and differential panel method;
Based on rigid motion theory, the rigid motion model is established using Newton's second law.
On the basis of the various embodiments described above, the structure module is specifically used for:
The rigid motion model and the elastic vibration model are subjected to simultaneous, constitute the bullet target action model.
It is apparent to those skilled in the art that for convenience and simplicity of description, the device of foregoing description Specific work process, can refer to preceding method in corresponding process, no longer excessively repeat herein.
In conclusion the embodiment of the present invention by build elastic vibration model, by elastic vibration model and rigid motion mould Type is combined composition and plays target action model, analyzes shadow of the warhead elastic vibration to warhead internal component more fully hereinafter It rings, to improve the accuracy for the stressing conditions for playing target action model reflection warhead internal component.
Fig. 5 is please referred to, Fig. 5 is electronic devices structure block diagram provided by the embodiments of the present application.Electronic equipment may include building Mold device 501, memory 502, storage control 503, processor 504, Peripheral Interface 505, input-output unit 506, audio Unit 507, display unit 508.
The memory 502, storage control 503, processor 504, Peripheral Interface 505, input-output unit 506, sound Frequency unit 507,508 each element of display unit are directly or indirectly electrically connected between each other, to realize the transmission or friendship of data Mutually.It is electrically connected for example, these elements can be realized between each other by one or more communication bus or signal wire.The modeling Device 501 can be stored in the memory 502 or be cured in the form of software or firmware (firmware) including at least one Software function module in the operating system (operating system, OS) of the model building device 501.The processor 504 are used to execute the executable module stored in memory 502, such as the software function module that the model building device 501 includes Or computer program.
Wherein, memory 502 may be, but not limited to, random access memory (Random Access Memory, RAM), read-only memory (Read Only Memory, ROM), programmable read only memory (Programmable Read-Only Memory, PROM), erasable read-only memory (Erasable Programmable Read-Only Memory, EPROM), Electricallyerasable ROM (EEROM) (Electric Erasable Programmable Read-Only Memory, EEPROM) etc.. Wherein, memory 502 is for storing program, and the processor 504 executes described program after receiving and executing instruction, aforementioned The method performed by server that the stream process that any embodiment of the embodiment of the present invention discloses defines can be applied to processor 504 In, or realized by processor 504.
Processor 504 can be a kind of IC chip, the processing capacity with signal.Above-mentioned processor 504 can To be general processor, including central processing unit (Central Processing Unit, abbreviation CPU), network processing unit (Network Processor, abbreviation NP) etc.;Can also be digital signal processor (DSP), application-specific integrated circuit (ASIC), Ready-made programmable gate array (FPGA) either other programmable logic device, discrete gate or transistor logic, discrete hard Part component.It may be implemented or execute disclosed each method, step and the logic diagram in the embodiment of the present invention.General processor Can be microprocessor or the processor 504 can also be any conventional processor etc..
The Peripheral Interface 505 couples various input/output devices to processor 504 and memory 502.At some In embodiment, Peripheral Interface 505, processor 504 and storage control 503 can be realized in one single chip.Other one In a little examples, they can be realized by independent chip respectively.
Input-output unit 506 is for being supplied to user input data to realize user and the server (or local terminal) Interaction.The input-output unit 506 may be, but not limited to, mouse and keyboard etc..
Audio unit 507 provides a user audio interface, may include that one or more microphones, one or more raises Sound device and voicefrequency circuit.
Display unit 508 provides an interactive interface (such as user interface) between the electronic equipment and user Or it is referred to user for display image data.In the present embodiment, the display unit 508 can be liquid crystal display or touch Control display.Can be that the capacitance type touch control screen or resistance-type of single-point and multi-point touch operation is supported to touch if touch control display Control screen etc..Single-point and multi-point touch operation is supported to refer to touch control display and can sense on the touch control display one or more The touch control operation generated simultaneously at a position, and transfer to processor 504 to be calculated and handled the touch control operation that this is sensed.
The Peripheral Interface 505 couples various input/output devices to processor 504 and memory 502.At some In embodiment, Peripheral Interface 505, processor 504 and storage control 503 can be realized in one single chip.Other one In a little examples, they can be realized by independent chip respectively.
The interaction that input-output unit 506 is used to that user input data to be supplied to realize user and processing terminal.It is described defeated Enter output unit 506 may be, but not limited to, mouse and keyboard etc..
It is appreciated that structure shown in fig. 5 is only to illustrate, the electronic equipment may also include it is more than shown in Fig. 5 or The less component of person, or with the configuration different from shown in Fig. 5.Each component shown in Fig. 5 may be used hardware, software or A combination thereof is realized.
In several embodiments provided herein, it should be understood that disclosed device and method can also pass through Other modes are realized.The apparatus embodiments described above are merely exemplary, for example, the flow chart in attached drawing and block diagram Show the device of multiple embodiments according to the present invention, the architectural framework in the cards of method and computer program product, Function and operation.In this regard, each box in flowchart or block diagram can represent the one of a module, section or code Part, a part for the module, section or code, which includes that one or more is for implementing the specified logical function, to be held Row instruction.It should also be noted that at some as in the realization method replaced, the function of being marked in box can also be to be different from The sequence marked in attached drawing occurs.For example, two continuous boxes can essentially be basically executed in parallel, they are sometimes It can execute in the opposite order, this is depended on the functions involved.It is also noted that every in block diagram and or flow chart The combination of box in a box and block diagram and or flow chart can use function or the dedicated base of action as defined in executing It realizes, or can be realized using a combination of dedicated hardware and computer instructions in the system of hardware.
In addition, each function module in each embodiment of the present invention can integrate to form an independent portion Point, can also be modules individualism, can also two or more modules be integrated to form an independent part.
It, can be with if the function is realized and when sold or used as an independent product in the form of software function module It is stored in a computer read/write memory medium.Based on this understanding, technical scheme of the present invention is substantially in other words The part of the part that contributes to existing technology or the technical solution can be expressed in the form of software products, the meter Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention. And storage medium above-mentioned includes:USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic disc or CD.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, any made by repair Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.It should be noted that:Similar label and letter exist Similar terms are indicated in following attached drawing, therefore, once being defined in a certain Xiang Yi attached drawing, are then not required in subsequent attached drawing It is further defined and is explained.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain Lid is within protection scope of the present invention.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
It should be noted that herein, relational terms such as first and second and the like are used merely to a reality Body or operation are distinguished with another entity or operation, are deposited without necessarily requiring or implying between these entities or operation In any actual relationship or order or sequence.Moreover, the terms "include", "comprise" or its any other variant are intended to Non-exclusive inclusion, so that the process, method, article or equipment including a series of elements is not only wanted including those Element, but also include other elements that are not explicitly listed, or further include for this process, method, article or equipment Intrinsic element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that There is also other identical elements in process, method, article or equipment including the element.

Claims (10)

  1. The target action model modeling method 1. a kind of Penetration is hit by a bullet, which is characterized in that including:
    The parameter of warhead is obtained, the parameter includes:Case material elasticity modulus, case material density, warhead length, war Pan quality and each rank axial vibration damping;
    The fight is obtained according to the case material elasticity modulus, the case material density and the warhead length computation The axial vibration frequency of each rank in portion;
    The fight is obtained based on the axial vibration frequency, the warhead quality and each rank axial vibration Damping calculating The mass spring damper model in portion;
    The transmission function of warhead axial vibration is built according to the mass spring damper model, and according to the transferometer It calculates and obtains elastic vibration model;
    The rigid motion model of the warhead is built, and is constituted with the elastic vibration models coupling and plays target action model.
  2. 2. according to the method described in claim 1, it is characterized in that, described according to the case material elasticity modulus, the shell Body density of material and the warhead length computation obtain the axial vibration frequency of the warhead, including:
    According to formulaCalculate the axial vibration frequency for obtaining the warhead, wherein fPlay iFor The axial vibration frequency of the i-th rank of the warhead, L are the warhead length, and E is the case material elasticity modulus, and ρ is institute State case material density.
  3. 3. according to the method described in claim 1, it is characterized in that, the parameter further includes the dimensionless resistance of each rank axial vibration Buddhist nun's ratio, correspondingly, described be based on the axial vibration frequency, the warhead quality and each rank axial vibration Damping calculating The mass spring damper model of the warhead is obtained, including:
    The natural angular frequency that each rank axial vibration is calculated according to the axial vibration frequency, according to the warhead quality and described Natural angular frequency calculates the rigidity of each rank axial vibration, according to the rigidity of each rank axial vibration, warhead matter Amount and the dimensionless damping ratio calculate each rank axial vibration damping;
    Institute is obtained according to the Rigidity Calculation of the warhead quality, each rank axial vibration damping and each rank axial vibration State the mass spring damper model of warhead.
  4. 4. according to the method described in claim 3, it is characterized in that, described calculate each rank axial direction according to the axial vibration frequency The natural angular frequency of vibration, including:
    According to ωPlay i=2 π fPlay iCalculate the natural angular frequency of each rank axial vibration, wherein fPlay iFor the axis of the i-th rank of the warhead To vibration frequency;
    The rigidity that each rank axial vibration is calculated according to the warhead quality and the natural angular frequency, including:
    According toCalculate the rigidity of each rank axial vibration, wherein MBulletFor the warhead quality;
    It is described that institute is calculated according to the rigidity, warhead quality and the dimensionless damping ratio of each rank axial vibration Each rank axial vibration damping is stated, including:
    According toCalculate each rank axial vibration damping, wherein ξPlay iFor the dimensionless damping ratio.
  5. 5. according to the method described in claim 1, it is characterized in that, the transmission function is:Wherein, Δ X (s) is the axial deflection, Fx(s) it is Play target active force, GPlay i(s) it is the transmission function of the i-th rank elastic vibration, andMBulletFor the warhead quality, s is multiple The independent variable of number field, CPlay iFor the damping of the i-th rank axial vibration, KPlay iFor the rigidity of the i-th rank axial vibration, ξPlay iIt is immeasurable for the i-th rank Guiding principle damping ratio;
    Correspondingly, described calculated according to the transmission function obtains elastic vibration model, including:
    The transmission function is subjected to partial differential calculating, the elastic vibration model of acquisition is:
  6. 6. according to the method described in claim 1, it is characterized in that, the rigid motion model of the structure warhead, packet It includes:
    The bullet target active force in Penetration is calculated according to cavity-expansion theory and differential panel method;
    Based on rigid motion theory, the rigid motion model is established using Newton's second law.
  7. 7. according to claim 1-6 any one of them methods, which is characterized in that described according to the rigid motion model and institute It states elastic vibration model and obtains bullet target action model, including:
    The rigid motion model and the elastic vibration model are subjected to simultaneous, constitute the bullet target action model.
  8. The target action model model building device 8. a kind of Penetration is hit by a bullet, which is characterized in that including:
    Acquisition module, the parameter for obtaining warhead, the parameter include:Case material elasticity modulus, case material density, Warhead length, warhead quality and each rank axial vibration damping;
    First computing module, for according to the case material elasticity modulus, the case material density and the fight minister Degree calculates the axial vibration frequency for obtaining each rank of warhead;
    Second computing module, for based on the axial vibration frequency, the warhead quality and each rank axial vibration resistance Buddhist nun calculates the mass spring damper model for obtaining the warhead;
    Third computing module, the transmission function for building warhead axial vibration according to the mass spring damper model, and It is calculated according to the transmission function and obtains elastic vibration model;
    Module, the rigid motion model for building the warhead are built, and bullet is constituted with the elastic vibration models coupling Target action model.
  9. 9. device according to claim 8, which is characterized in that first computing module is specifically used for:
    According to formulaCalculate the axial vibration frequency for obtaining the warhead, wherein fPlay iFor The axial vibration frequency of the i-th rank of the warhead, L are the warhead length, and E is the case material elasticity modulus, and ρ is institute State case material density.
  10. 10. device according to claim 8, which is characterized in that the parameter further includes the dimensionless of each rank axial vibration Damping ratio, correspondingly, second computing module, is specifically used for:
    The natural angular frequency that each rank axial vibration is calculated according to the axial vibration frequency, according to the warhead quality and described Natural angular frequency calculates the rigidity of each rank axial vibration, according to the rigidity of each rank axial vibration, warhead matter Amount and the dimensionless damping ratio calculate each rank axial vibration damping;
    Institute is obtained according to the Rigidity Calculation of the warhead quality, each rank axial vibration damping and each rank axial vibration State the mass spring damper model of warhead.
CN201810391492.XA 2018-04-27 2018-04-27 Method and device for modeling bullet target action model in penetration process Expired - Fee Related CN108595848B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810391492.XA CN108595848B (en) 2018-04-27 2018-04-27 Method and device for modeling bullet target action model in penetration process

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810391492.XA CN108595848B (en) 2018-04-27 2018-04-27 Method and device for modeling bullet target action model in penetration process

Publications (2)

Publication Number Publication Date
CN108595848A true CN108595848A (en) 2018-09-28
CN108595848B CN108595848B (en) 2022-02-01

Family

ID=63610707

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810391492.XA Expired - Fee Related CN108595848B (en) 2018-04-27 2018-04-27 Method and device for modeling bullet target action model in penetration process

Country Status (1)

Country Link
CN (1) CN108595848B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150123346A1 (en) * 2013-05-21 2015-05-07 Gregory T Mason Mason Target System
CN106214182A (en) * 2016-07-06 2016-12-14 西安交通大学 HIFU damaging shear ripple elastic characteristic method of estimation based on LK optical flow method
CN106446517A (en) * 2016-08-30 2017-02-22 江苏大学 Modeling and judgment method for laser shocking peening alloy plastic deformation depth

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150123346A1 (en) * 2013-05-21 2015-05-07 Gregory T Mason Mason Target System
CN106214182A (en) * 2016-07-06 2016-12-14 西安交通大学 HIFU damaging shear ripple elastic characteristic method of estimation based on LK optical flow method
CN106446517A (en) * 2016-08-30 2017-02-22 江苏大学 Modeling and judgment method for laser shocking peening alloy plastic deformation depth

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
X Z KONG ET AL: "Rigid and eroding projectile penetration into concrete targets based on an extended dynamic cavity expansion model", 《INTERNATIONAL JOURNAL OF IMPACT ENGINEERING》 *
李艳: "侵爆荷载作用下半无限岩土介质的力学性能及边坡稳定性研究", 《中国博士学位论文全文数据库电子期刊 工程科技II辑》 *
见雷雷: "钨合金结构侵彻弹体应用分析研究", 《中国优秀硕士学位论文全文数据库电子期刊 工程科技II辑》 *

Also Published As

Publication number Publication date
CN108595848B (en) 2022-02-01

Similar Documents

Publication Publication Date Title
CN102395409A (en) Systems and methods for an impact location and amplitude sensor
CN106779825A (en) A kind of item recommendation method, device and electronic equipment
Wiebe et al. On the experimental identification of unstable static equilibria
CN108827455B (en) A kind of jointed rock masses under blasting vibration attenuation parameter prediction technique and device
Qi et al. Experimental studying multi-impact phenomena exhibited during the collision of a sphere onto a steel beam
Das et al. Alignment of interstellar grains by mechanical torques: suprathermally rotating Gaussian random spheres
US10398980B2 (en) Method and apparatus for providing a computer implemented game
Zhao et al. Design method of dynamic parameters of a self-synchronization vibrating system with dual mass
Morais et al. Shock environment design for space equipment testing
CN108595848A (en) A kind of Penetration is hit by a bullet target action model modeling method and device
Troncossi et al. Identification of polymer concrete damping properties
Zhu et al. Three-dimensional dynamic modeling and analysis of moving elevator traveling cables
CN109740621A (en) A kind of video classification methods, device and equipment
US11951407B2 (en) Screen soccer system and method of providing screen soccer
CN108109052A (en) Article method for pushing, device and server
CN116628600A (en) Unbalanced data sampling method and device based on random forest
Lin et al. Design and implement a mobile badminton stroke classification system
EP3352123A1 (en) Server device, program, and method
CN105242854A (en) Method, apparatus and mobile terminal capable of realizing user interaction in collisions
Li et al. Parameter optimization and mechanism research of double-layer variable screen surface vibrating screen
Ganguly et al. A viscoelastic Timoshenko shaft finite element for dynamic analysis of rotors
Thompson et al. Multi-modal generative adversarial networks for synthesizing time-series structural impact responses
Kesilmiş et al. A geometric approach to beam type load cell response for fast weighing
Ni et al. Design and experimentation of piezoelectric crystal sensor array for grain cleaning loss
Hashemian et al. A comprehensive study of the effects of various design variables on the dynamic behavior and vibration characteristics of a front-loading washing machine

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220201

CF01 Termination of patent right due to non-payment of annual fee