CN112325713B - Analysis method for pneumatic nonlinear angular motion characteristics of double spinning bombs - Google Patents

Analysis method for pneumatic nonlinear angular motion characteristics of double spinning bombs Download PDF

Info

Publication number
CN112325713B
CN112325713B CN202011173044.6A CN202011173044A CN112325713B CN 112325713 B CN112325713 B CN 112325713B CN 202011173044 A CN202011173044 A CN 202011173044A CN 112325713 B CN112325713 B CN 112325713B
Authority
CN
China
Prior art keywords
equation
spinning projectile
double
nonlinear
motion
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.)
Active
Application number
CN202011173044.6A
Other languages
Chinese (zh)
Other versions
CN112325713A (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Publication of CN112325713A publication Critical patent/CN112325713A/en
Application granted granted Critical
Publication of CN112325713B publication Critical patent/CN112325713B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42—AMMUNITION; BLASTING
    • F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
    • F42B35/00—Testing or checking of ammunition

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

The invention discloses a method for analyzing the characteristic of pneumatic nonlinear angular motion of a double spinning projectile, which comprises the following steps of: the method comprises the following steps: establishing the speed omega of the double spinning projectile relative to the attack angle alpha, the sideslip angle beta and the pitch anglezYaw rate ωyThe non-linear angular motion equation of (1); step two: determining a nonlinear complex attack angle equation of the spinning projectile according to the nonlinear angular motion equation established in the step one; step three: converting the stability analysis of the nonlinear complex attack angle equation solution into the analysis of the corresponding homogeneous equation solution; step four: determining an amplitude plane equation of the double spinning projectile under the action of the cubic static moment; step five: determining the necessary condition that the double spinning projectile generates stable conical motion under the action of cubic static moment; starting from pneumatic nonlinearity, the necessary condition that the double spinning projectile generates stable conical swing under the action of the cubic static moment is obtained, so that the reason that the generated attack angle is not attenuated and even the flight instability phenomenon is caused under the action of the nonlinear static moment is pointed out.

Description

Analysis method for pneumatic nonlinear angular motion characteristics of double spinning bombs
Technical Field
The invention relates to the technical field of analysis of angular motion characteristics of double spinning bullets, in particular to a method for analyzing the angular motion characteristics of pneumatic nonlinearity of the double spinning bullets.
Background
The rapid development of conventional projectile guidance technology is fueled by the tremendous demand in modern warfare for high-precision, low-cost, beyond-the-horizon, small collateral damage guided munitions. Aiming at the characteristic of high rotating speed of the rotation stabilization projectile, the duck rudder component is adopted to replace an original ammunition fuse to form a 'double-rotation' structure, and then the intelligent and flexible transformation of the rotation stabilization projectile is realized.
However, the double-spin projectile dynamics has the characteristics of strong coupling and strong nonlinearity, and meanwhile, the structural change and the introduction of a control system make a plurality of nonlinear motion phenomena of the projectile difficult to analyze and explain by using the theory of a linear system. The nonlinear characteristics severely restrict the application of the double-spin ammunition in the development of weaponry. Therefore, the method has important theoretical research value and engineering application significance for the research of the nonlinear dynamics analysis of the spinning projectile.
Disclosure of Invention
In view of the above, the invention provides a method for analyzing the pneumatic nonlinear lower angular motion characteristic of a double spinning projectile, and a necessary condition that the double spinning projectile generates stable conical oscillation under the action of cubic static moment is obtained by adopting an angular motion quasi-linear analysis method based on an averaging method.
The technical scheme of the invention is as follows: a method for analyzing the characteristic of the pneumatic nonlinear angular motion of a spinning projectile comprises the following steps:
the method comprises the following steps: establishing the speed omega of the double spinning projectile relative to the attack angle alpha, the sideslip angle beta and the pitch anglezYaw rate ωyThe non-linear angular motion equation of (1);
step two: determining a nonlinear complex attack angle equation of the spinning projectile according to the nonlinear angular motion equation established in the step one;
step three: converting the stability analysis of the nonlinear complex attack angle equation solution into the analysis of the corresponding homogeneous equation solution;
step four: determining an amplitude plane equation of the double spinning projectile under the action of the cubic static moment;
step five: determining the necessary condition of stable conical motion generated by the double spinning projectile under the action of the cubic static moment.
Preferably, in the first step:
selecting x ═ alpha beta omegaz ωy]TAccording to Newton's second law and momentum moment theorem, the equation of the nonlinear angular motion of the double spinning projectile is obtained as follows:
Figure GDA0002838041360000021
in the formula (I), the compound is shown in the specification,
Figure GDA0002838041360000022
Figure GDA0002838041360000023
wherein ρ, S, L, g,
Figure GDA00028380413600000210
m is air density, characteristic area of the double spinning projectile, characteristic length of the double spinning projectile, gravity acceleration of the double spinning projectile, pitch angle of the double spinning projectile, mass of the double spinning projectile, Jy＝Jz，JyAnd JzAre equatorial moments of inertia, JxIs the polar moment of inertia;
Figure GDA0002838041360000024
and
Figure GDA0002838041360000029
respectively representing the lift force, the static moment, the damping moment, the Magnus moment and the control moment generated by the attack angle.
Preferably, in the second step:
consider mΔNon-linearly induced conical motion, selecting
Figure GDA00028380413600000211
The complex variables delta beta + i alpha and delta are defined by a complex analysis method in Murphy stability theoryy+iδzEliminating ω from formula (1)z、ωyAnd the derivative thereof, the nonlinear complex attack angle equation of the spinning projectile is obtained as follows:
Figure GDA0002838041360000026
wherein H is cΔ-mΩDamping characteristics characterizing angular movement, M ═ MΔ+mΩcΔDetermines the frequency of angular movement, T ═ mω/P+cΔMiddle mωWhich is a magnus moment, affects the stability of the angular movement.
Preferably, in the third step:
m is expressed as M ═ M0+M2Δ2Δ ═ Δ |, the stability of the solution of equation (2) from the ordinary differential equation is determined by its corresponding solution of the homogeneous equation:
Figure GDA0002838041360000027
preferably, in step four:
let the two-circle motion form of the solution of the formula (3) be:
Figure GDA0002838041360000028
defining a damping factor lambdai：
Figure GDA0002838041360000031
Substituting the formula (4) and the formula (5) into the formula (3) to obtain:
Figure GDA0002838041360000032
on both sides of formula (6), is divided by
Figure GDA0002838041360000033
And order
Figure GDA0002838041360000034
Is amplified at
Figure GDA0002838041360000035
Averaged over one period of (a):
Figure GDA0002838041360000036
omitting its derivative
Figure GDA0002838041360000037
While omitting a small amount lambda1(λ1+ H) separating the real and imaginary parts of equation (7) to obtain:
Figure GDA0002838041360000038
Figure GDA0002838041360000039
resolving the frequency from the formula (8)
Figure GDA00028380413600000310
And
Figure GDA00028380413600000311
the substitution into formula (9) gives:
Figure GDA00028380413600000312
obtaining the following by the same method:
Figure GDA00028380413600000313
the stability problem of the fourth-order system is reduced to the stability problem of the second-order system, and the equations (10) and (11) are combined to construct
Figure GDA00028380413600000314
Figure GDA00028380413600000315
Amplitude plane equation for coordinate plane:
Figure GDA00028380413600000316
preferably, in the step five:
when in use
Figure GDA00028380413600000317
When it is used, order
Figure GDA00028380413600000318
Then, the equations (10) and (11) can be simplified to obtain:
Figure GDA00028380413600000319
Figure GDA0002838041360000041
here, λ1And λ2The expression is meaningful to satisfy:
Figure GDA0002838041360000042
at K2＜＜K1Under the condition (2), the necessary condition for the double-spin projectile to form a conical motion is that
Figure GDA00028380413600000413
With only one singular point on the axis
Figure GDA0002838041360000043
The singularity should be a zero damping factor curve lambda10 and
Figure GDA0002838041360000044
intersection of axes, let λ10 is gotTo:
Figure GDA0002838041360000045
for normal flying positive spin double spin ammunition, there is P>0，H>0, then singular point
Figure GDA0002838041360000046
The following constraints should be satisfied:
Figure GDA0002838041360000047
to make coordinate changes
Figure GDA0002838041360000048
Moving the origin of the coordinate system to an odd point, performing coordinate transformation by the equation (12), and substituting the coordinate transformation into the equations (13) and (14) to obtain a new amplitude plane equation:
Figure GDA0002838041360000049
(18) equation Jacobian matrix at equilibrium point (0,0)
Figure GDA00028380413600000410
Wherein b is1＝0，c1＝0，
Figure GDA00028380413600000411
Figure GDA00028380413600000412
The criterion of Lyapunov stability is known, when a1<0，d1<0, equilibrium point (0,0) is stable; the requirements for obtaining stable conic motion of double spin elasticity by combining the formulas (15) and (17) are as follows:
Figure GDA0002838041360000051
all in one
Figure GDA0002838041360000052
When deriving the
Figure GDA0002838041360000053
The requirements for stable conical motion of the double spin projectile are as follows:
Figure GDA0002838041360000054
wherein the content of the first and second substances,
Figure GDA0002838041360000055
Figure GDA0002838041360000056
has the advantages that:
(1) the invention starts from pneumatic nonlinearity, obtains the necessary condition that the bispin bomb generates stable conical swing under the action of cubic static moment, thereby indicating the reasons that the generated attack angle is not attenuated and even the flight instability phenomenon is caused under the action of nonlinear static moment, and having guiding significance for the structure and pneumatic design of the bispin bomb.
Detailed Description
The present invention will be described in detail below with reference to examples.
The embodiment provides a method for analyzing pneumatic nonlinear lower angular motion characteristics of a double-spinning projectile, and the method is based on an angular motion quasi-linear analysis method of an averaging method, so that the necessary condition that the double-spinning projectile generates stable conical swing under the action of cubic static moment is obtained.
The analysis method mainly comprises the following steps:
the method comprises the following steps: establishing the speed omega of the double spinning projectile relative to the attack angle alpha, the sideslip angle beta and the pitch anglezYaw rate ωyThe non-linear angular motion equation of (1);
to characterize angular motion, assume linear velocity V of the dual spinning projectile and rotational velocity ω of the rear bodyxSelecting x ═ alpha beta omega as slow variablez ωy]TIs the state quantity of the nonlinear angular motion of the double spinning projectile, and alpha and beta are relative to V and omega in the stable flying process of the double spinning projectilexThe fast variable is a small variable (when the attack angle and the sideslip angle are small, they can be approximated in the derivation process, that is, sin (α) ═ α, sin (β) ═ β, cos (α) ═ 1, cos (β) ═ 1), and the magnus force and the control surface control force are smaller than the aerodynamic force generated by the attack angle by two or more orders (i.e., the influence of the magnus force and the control surface control force on the angular motion of the double-spinning projectile can be ignored), so that the equation of the double-spinning projectile nonlinear angular motion can be obtained according to the second newton's law and the theorem of the moment of momentum as follows:
Figure GDA0002838041360000061
in the formula (I), the compound is shown in the specification,
Figure GDA0002838041360000062
Figure GDA0002838041360000063
wherein ρ, S, L, g,
Figure GDA00028380413600000610
m is air density, characteristic area of the double spinning projectile, characteristic length of the double spinning projectile, gravity acceleration of the double spinning projectile, pitch angle of the double spinning projectile and mass of the double spinning projectile respectively, JyAnd JzAre equatorial moments of inertia, J is due to the axial symmetry of the double spin projectiley＝Jz，JxIs the polar moment of inertia;
Figure GDA0002838041360000064
and
Figure GDA0002838041360000065
respectively representing a lift force, a static moment, a damping moment, a Magnus moment and a control moment generated by an attack angle;
step two: determining a nonlinear complex attack angle equation of the double spinning projectile;
for convenience of theoretical analysis, for the equation of the nonlinear angular motion of the double-spinning projectile described by the formula (21), the present embodiment mainly performs analysis on the phenomenon of the nonlinear motion of the double-spinning projectile from the pneumatic nonlinearity, and for the pneumatic nonlinearity, because the moment applied to the double-spinning projectile plays a dominant role in the angular motion, the cubic static moment applied to the double-spinning projectile is considered in a critical manner;
introducing horizontal ballistic assumptions, i.e. selection
Figure GDA0002838041360000069
The complex variables delta beta + i alpha and delta are defined by a complex analysis method in Murphy stability theoryy+iδzWherein, deltazFor pitching rudder deflection angle, deltayFor yaw rudder angle, the ω is eliminated from the equation (21)z、ωyAnd the derivative thereof, the nonlinear complex attack angle equation of the spinning projectile is obtained as follows:
Figure GDA0002838041360000066
in the formula (I), the compound is shown in the specification,
Figure GDA0002838041360000067
and
Figure GDA0002838041360000068
specifying the first and second derivatives of a complex-defining variable Δ, respectively, H ═ cΔ-mΩH characterizes the damping characteristic of angular motion, M ═ MΔ+mΩcΔM determinesHaving a frequency of angular movement, T ═ mω/P+cΔT is related to the Magnus moment (m)ωMagnus moments) that affect the stability of angular motion;
step three: converting the stability analysis of the nonlinear complex attack angle equation solution into the analysis of the corresponding homogeneous equation solution;
considering that the solution of the homogeneous equation is equivalent to the stability characterized by the solution of the nonlinear complex attack angle equation, and the solution of the homogeneous equation is easier to solve, meanwhile, the embodiment considers the influence of the cubic static moment on the angular motion, and M can be expressed as M ═ M0+M2Δ2Wherein M is0Is a linear term of moment, M2For a non-linear term of moment, Δ ═ Δ |, as can be derived from the knowledge of ordinary differential equations, the stability of the solution of equation (22) can be determined by its corresponding solution of the following homogeneous equation:
Figure GDA0002838041360000071
step four: determining an amplitude plane equation of the double spinning projectile under the action of the cubic static moment;
since the amplitude reflects the magnitude of the angular motion, let the form of the two-circle motion of the solution of equation (23) be:
Figure GDA0002838041360000072
wherein, K1And K2Is a die with two-circle motion,
Figure GDA0002838041360000073
and
Figure GDA0002838041360000074
the phase of the two-circle motion;
defining a damping factor lambdai：
Figure GDA0002838041360000075
Figure GDA0002838041360000076
Is KiA first derivative with respect to time;
substituting the formula (24) or (25) into the formula (23) to obtain:
Figure GDA0002838041360000077
wherein the content of the first and second substances,
Figure GDA0002838041360000078
and
Figure GDA0002838041360000079
are each lambdai，λi，
Figure GDA00028380413600000710
And
Figure GDA00028380413600000711
a first derivative with respect to time;
Figure GDA00028380413600000712
and
Figure GDA00028380413600000713
is composed of
Figure GDA00028380413600000714
And
Figure GDA00028380413600000715
a second derivative with respect to time;
on both sides of formula (26), divided by
Figure GDA00028380413600000716
And order
Figure GDA00028380413600000717
Is amplified at
Figure GDA00028380413600000718
Averaged over one period of (a):
Figure GDA00028380413600000719
due to lambda1Slowly changing, omitting its derivatives
Figure GDA00028380413600000720
While omitting a small amount lambda1(λ1+ H) separating the real and imaginary parts of equation (27) to obtain:
Figure GDA00028380413600000721
Figure GDA00028380413600000722
resolving the frequency from equation (28)
Figure GDA00028380413600000723
And
Figure GDA00028380413600000724
the substitution into the formula (29) gives:
Figure GDA00028380413600000725
the same can get:
Figure GDA0002838041360000081
the stability problem of the fourth-order system is reduced to that of the second-order system, and combined with the equations (30) and (31), the system can be constructed to
Figure GDA0002838041360000082
Amplitude plane equation for coordinate plane:
Figure GDA0002838041360000083
step five: determining the necessary condition that the double spinning projectile generates stable conical motion under the action of cubic static moment;
for the motion stability of the double spinning projectile, only the mode K of the two-circle motion is required1、K2If the double-spin projectile does not diverge, the motion is stable, and when the double-spin projectile makes two-circle motion, the mode K of two partial motions of the complex attack angle delta1、K2With K1＜＜K2Or K2＜＜K1Without setting K2＜＜K1(ii) a On the other hand, in the prior art, two angular movement frequencies
Figure GDA0002838041360000084
And
Figure GDA0002838041360000085
generally, the difference is large, and the present embodiment will be right
Figure GDA0002838041360000086
And
Figure GDA0002838041360000087
respectively discussing;
when in use
Figure GDA0002838041360000088
When it is used, order
Figure GDA0002838041360000089
Then, the equations (30) and (31) can be simplified to obtain:
Figure GDA00028380413600000810
Figure GDA00028380413600000811
here, λ1And λ2The expression is meaningful to satisfy:
Figure GDA00028380413600000812
at K2＜＜K1Under the condition (2), the necessary condition for the double-spin projectile to form a conical motion is that
Figure GDA00028380413600000813
With only one singular point on the axis
Figure GDA00028380413600000814
KeFor the final convergence of the conic motion of the spinning projectile, the singularity should be the zero damping factor curve λ10 and
Figure GDA00028380413600000815
intersection of axes, let λ10 can result in:
Figure GDA00028380413600000816
for normal flying positive spin double spin ammunition, there is P>0，H>0, then singular point
Figure GDA00028380413600000817
The following constraints should be satisfied:
Figure GDA00028380413600000818
for ease of discussion, coordinate changes are made
Figure GDA0002838041360000091
Moving the origin of the coordinate system to an odd point, performing coordinate transformation by the equation (12), and substituting the coordinate transformation into the equations (13) and (14) to obtain a new amplitude plane equation:
Figure GDA0002838041360000092
(38) equation Jacobian matrix at equilibrium point (0,0)
Figure GDA0002838041360000093
Wherein b is1＝0，c1＝0，
Figure GDA0002838041360000094
Figure GDA0002838041360000095
The criterion of Lyapunov stability is known, when a1<0，d1<0, equilibrium point (0,0) is stable; the requirements for obtaining stable conic motion of double spin can be obtained by combining the formulas (35) and (37):
Figure GDA0002838041360000096
all in one
Figure GDA0002838041360000097
When deriving the
Figure GDA0002838041360000098
The requirements for stable conical motion of the double spin projectile are as follows:
Figure GDA0002838041360000099
wherein the content of the first and second substances,
Figure GDA00028380413600000910
Figure GDA0002838041360000101
in summary, the equations (39) and (40) constitute the requirement for stabilizing the conical motion of the double spinning projectile under the action of the cubic static moment.
In summary, the above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (5)

1. A method for analyzing the pneumatic nonlinear angular motion characteristics of a spinning projectile is characterized by comprising the following steps:
the method comprises the following steps: establishing the speed omega of the double spinning projectile relative to the attack angle alpha, the sideslip angle beta and the pitch anglezYaw rate ωyThe non-linear angular motion equation of (1);
step two: determining a nonlinear complex attack angle equation of the spinning projectile according to the nonlinear angular motion equation established in the step one;
step three: converting the stability analysis of the nonlinear complex attack angle equation solution into the analysis of the corresponding homogeneous equation solution;
step four: determining an amplitude plane equation of the double spinning projectile under the action of the cubic static moment;
step five: determining the necessary condition that the double spinning projectile generates stable conical motion under the action of cubic static moment;
in the first step:
selecting x ═ alpha beta omegaz ωy]TAccording to Newton's second law and momentum moment theorem, the equation of the nonlinear angular motion of the double spinning projectile is obtained as follows:
Figure FDA0003182481350000011
in the formula (I), the compound is shown in the specification,
Figure FDA0003182481350000012
Figure FDA0003182481350000013
wherein V is the linear velocity of the double spinning projectile, omegaxRho, S, L, g, theta and m are respectively the air density, the characteristic area of the double spinning projectile, the characteristic length of the double spinning projectile, the gravity acceleration of the double spinning projectile, the pitch angle of the double spinning projectile, the mass of the double spinning projectile, and Jy＝Jz，JyAnd JzAre equatorial moments of inertia, JxIs the polar moment of inertia;
Figure FDA0003182481350000014
and
Figure FDA0003182481350000015
respectively representing the lift force, the static moment, the damping moment, the Magnus moment and the control moment generated by the attack angle.
2. The method for analyzing the characteristic of the pneumatic nonlinear lower angular motion of the spinning projectile as set forth in claim 1, wherein in the second step:
consider mΔSelecting theta as 0 degree, and defining complex variables delta as beta + i alpha and delta as delta by using a complex analysis method in Murphy stability theoryy+iδzWherein, deltayFor yaw rudder angle, deltazIs the pitch rudder deflection angle; elimination of omega from formula (1)z、ωyAnd the derivative thereof, the nonlinear complex attack angle equation of the spinning projectile is obtained as follows:
Figure FDA0003182481350000021
wherein H is cΔ-mΩDamping characteristics characterizing angular movement, M ═ MΔ+mΩcΔDetermines the frequency of angular movement, T ═ mω/P+cΔMiddle mωWhich is a magnus moment, affects the stability of the angular movement.
3. The method for analyzing the characteristic of the pneumatic nonlinear lower angular motion of the spinning projectile as claimed in claim 2, wherein in the third step:
m is expressed as M ═ M0+M2Δ2Δ ═ Δ |, the stability of the solution of equation (2) from the ordinary differential equation is determined by its corresponding solution of the homogeneous equation:
Figure FDA0003182481350000022
4. the method for analyzing the characteristics of the pneumatic nonlinear lower angular motion of the spinning projectile as recited in claim 3, wherein in the fourth step:
let the two-circle motion form of the solution of the formula (3) be:
Figure FDA0003182481350000023
wherein, K1And K2Is a die with two-circle motion,
Figure FDA0003182481350000024
and
Figure FDA0003182481350000025
the phase of the two-circle motion;
defining a damping factor lambdai：
Figure FDA0003182481350000026
Substituting the formula (4) and the formula (5) into the formula (3) to obtain:
Figure FDA0003182481350000027
on both sides of formula (6), is divided by
Figure FDA0003182481350000028
And order
Figure FDA0003182481350000029
Is amplified at
Figure FDA00031824813500000210
Averaged over one period of (a):
Figure FDA00031824813500000211
omitting its derivative
Figure FDA00031824813500000212
While omitting a small amount lambda1(λ1+ H) separating the real and imaginary parts of equation (7) to obtain:
Figure FDA00031824813500000213
Figure FDA00031824813500000214
resolving the frequency from the formula (8)
Figure FDA00031824813500000215
And
Figure FDA00031824813500000216
the substitution into formula (9) gives:
Figure FDA0003182481350000031
obtaining the following by the same method:
Figure FDA0003182481350000032
the stability problem of the fourth-order system is reduced to the stability problem of the second-order system, and the equations (10) and (11) are combined to construct
Figure FDA00031824813500000316
Figure FDA00031824813500000315
Amplitude plane equation for coordinate plane:
Figure FDA0003182481350000033
5. the method for analyzing the characteristic of the pneumatic nonlinear lower angular motion of the spinning projectile as set forth in claim 4, wherein in the fifth step:
when in use
Figure FDA0003182481350000034
When it is used, order
Figure FDA0003182481350000035
Then, the equations (10) and (11) can be simplified to obtain:
Figure FDA0003182481350000036
Figure FDA0003182481350000037
here, λ1And λ2The expression is meaningful to satisfy:
Figure FDA0003182481350000038
at K2＜＜K1Under the condition (2), the necessary condition for the double-spin projectile to form a conical motion is that
Figure FDA0003182481350000039
With only one singular point on the axis
Figure FDA00031824813500000310
The singularity should be a zero damping factor curve lambda10 and
Figure FDA00031824813500000311
intersection of axes, let λ1Obtaining when the yield is 0:
Figure FDA00031824813500000312
for normal flying positive spin double spin ammunition, there is P>0，H>0, then singular point
Figure FDA00031824813500000313
The following constraints should be satisfied:
Figure FDA00031824813500000314
to make coordinate changes
Figure FDA0003182481350000041
Moving the origin of the coordinate system to an odd point, performing coordinate transformation by the equation (12), and substituting the coordinate transformation into the equations (13) and (14) to obtain a new amplitude plane equation:
Figure FDA0003182481350000042
(18) equation Jacobian matrix at equilibrium point (0,0)
Figure FDA0003182481350000043
Wherein b is1＝0，c1＝0，
Figure FDA0003182481350000044
Figure FDA0003182481350000045
The criterion of Lyapunov stability is known, when a1<0，d1<0, equilibrium point (0,0) is stable; the requirements for obtaining stable conic motion of double spin elasticity by combining the formulas (15) and (17) are as follows:
Figure FDA0003182481350000046
all in one
Figure FDA0003182481350000047
When deriving the
Figure FDA0003182481350000048
The stable conical motion of the double-spinning projectile must be madeThe following conditions are required:
Figure FDA0003182481350000049
wherein the content of the first and second substances,
Figure FDA00031824813500000410
Figure FDA00031824813500000411
CN202011173044.6A 2019-12-24 2020-10-28 Analysis method for pneumatic nonlinear angular motion characteristics of double spinning bombs Active CN112325713B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201911349090 2019-12-24
CN2019113490904 2019-12-24

Publications (2)

Publication Number Publication Date
CN112325713A CN112325713A (en) 2021-02-05
CN112325713B true CN112325713B (en) 2021-09-21

Family

ID=74297154

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011173044.6A Active CN112325713B (en) 2019-12-24 2020-10-28 Analysis method for pneumatic nonlinear angular motion characteristics of double spinning bombs

Country Status (1)

Country Link
CN (1) CN112325713B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113821996B (en) * 2021-07-14 2024-01-30 南京理工大学 Novel rapid calculation method for high-speed entry trajectory of projectile

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013022507A2 (en) * 2011-05-13 2013-02-14 Leigh Aerosystems Corporation Ground-projectile guidance system
CN104267733A (en) * 2014-10-25 2015-01-07 哈尔滨工业大学 Attitude control type direct lateral force and aerodynamic force composite missile attitude control method based on mixed forecasting control
CN106886625A (en) * 2017-01-05 2017-06-23 北京航天自动控制研究所 A kind of Design of Aerodynamic Configuration method of the dual-spin stabilization bullet based on fixed-wing duck rudder
CN107726927A (en) * 2017-10-19 2018-02-23 西安交通大学 A kind of amendment torque control system based on fixed-wing duck rudder bispin bullet
CN109614633A (en) * 2018-10-25 2019-04-12 南京航空航天大学 A kind of composite rotor craft non-linear modeling method and Calculate Ways

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9702674B2 (en) * 2014-08-07 2017-07-11 The United States Of America As Represented By The Secretary Of The Army Method and apparatus for GPS-denied navigation of spin-stabilized projectiles

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013022507A2 (en) * 2011-05-13 2013-02-14 Leigh Aerosystems Corporation Ground-projectile guidance system
CN104267733A (en) * 2014-10-25 2015-01-07 哈尔滨工业大学 Attitude control type direct lateral force and aerodynamic force composite missile attitude control method based on mixed forecasting control
CN106886625A (en) * 2017-01-05 2017-06-23 北京航天自动控制研究所 A kind of Design of Aerodynamic Configuration method of the dual-spin stabilization bullet based on fixed-wing duck rudder
CN107726927A (en) * 2017-10-19 2018-02-23 西安交通大学 A kind of amendment torque control system based on fixed-wing duck rudder bispin bullet
CN109614633A (en) * 2018-10-25 2019-04-12 南京航空航天大学 A kind of composite rotor craft non-linear modeling method and Calculate Ways

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
非线性力矩作用下气动偏心弹丸强迫圆锥运动稳定性条件;舒敬荣等;《兵工学报》;20180531;第39卷(第5期);第875-882页 *

Also Published As

Publication number Publication date
CN112325713A (en) 2021-02-05

Similar Documents

Publication Publication Date Title
Fusato et al. Flight dynamics of an articulated rotor helicopter with an external slung load
CN112578805B (en) Attitude control method of rotor craft
CN112325713B (en) Analysis method for pneumatic nonlinear angular motion characteristics of double spinning bombs
CN111176315B (en) Variant cruise missile control method based on L1 adaptive control
CN110209192A (en) Fighter plane course augmentation control design method
CN107977009A (en) A kind of airbreather attitude control law design method for considering coupling
CN109460055B (en) Aircraft control capability determining method and device and electronic equipment
CN111240204B (en) Model reference sliding mode variable structure control-based flying projectile patrol control method
Wernert et al. Modelling and stability analysis for a class of 155 mm spin-stabilized projectiles with course correction fuse (CCF)
CN115930698A (en) Method for converting rudder deflection angle of control channel into X-shaped rudder deflection angle under non-zero rolling angle condition
CN113885543B (en) Design method of navigation body auto-disturbance rejection roll controller under complex sea condition disturbance
CN106406333A (en) Integral-type-terminal-sliding-form-based method for tracking pitching angle of stratospheric airship
Yuanchuan et al. Observer-based adaptive sliding mode backstepping output-feedback DSC for spin-stabilized canard-controlled projectiles
CN109634110B (en) Tailstock type unmanned aerial vehicle transition section robust control method based on nonlinear disturbance observer
He et al. Coning motion stability of spinning missiles with strapdown seekers
Mao et al. Coning motion stability of wrap around fin rockets
Stefański et al. Analysis of the guiding of bombs on ground targets using a gyroscope system
Yao et al. Trajectory tracking controller based on PID-NLADRC
CN112417650B (en) Analysis method for geometrical nonlinear angular motion characteristics of double spinning bombs
Ye et al. The stability of rolling motion of hypersonic vehicles with slender configuration under pitching maneuvering
Yu et al. Flight Trajectory Correction of Target Missile and Design of Tracking Guidance
Özkan et al. Comparison of the notable acceleration-and angle-based guidance laws for a short-range air-to-surface missile
Li et al. Attitude coupling control and lateral stability analysis of hypersonic vehicle
Sun Research on Thrust Vector/Aerodynamics Compound Control Method Based on Linear Quadratic Regulator Control for Solid Rocket
Tan et al. Analysis of aerodynamic characteristics for a flying-wing UAV with asymmetric wing damage

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