CN113342049B - Multi-light source laser aircraft trajectory design method and system - Google Patents
Multi-light source laser aircraft trajectory design method and system Download PDFInfo
- Publication number
- CN113342049B CN113342049B CN202110709788.3A CN202110709788A CN113342049B CN 113342049 B CN113342049 B CN 113342049B CN 202110709788 A CN202110709788 A CN 202110709788A CN 113342049 B CN113342049 B CN 113342049B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- laser
- target laser
- light source
- target
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 81
- 230000001133 acceleration Effects 0.000 claims abstract description 54
- 230000008569 process Effects 0.000 claims abstract description 36
- 239000000126 substance Substances 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 23
- 238000005474 detonation Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 6
- 230000035939 shock Effects 0.000 claims description 4
- 239000007789 gas Substances 0.000 description 15
- 238000010586 diagram Methods 0.000 description 8
- 239000003380 propellant Substances 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 2
- 238000000608 laser ablation Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000004080 punching Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 108091092878 Microsatellite Proteins 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000001307 helium Substances 0.000 description 1
- 229910052734 helium Inorganic materials 0.000 description 1
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Lasers (AREA)
Abstract
The invention provides a method and a system for designing a trajectory of a multi-light-source laser aircraft, which are applied to a target laser aircraft, wherein the target laser aircraft provides laser beams for at least two light sources, and the model of the target laser aircraft is a Myrabo light ship structure; the method comprises the following steps: determining structural model parameters of the target laser aircraft; determining position information of a light source providing a laser beam for the target laser aircraft, and determining an initial pitch angle of the target laser aircraft; calculating dynamic model parameters of the target laser aircraft in an acceleration process and calculating pneumatic parameters of the target laser aircraft in an air suction mode based on the structural model parameters; and calculating the trajectory of the target laser aircraft based on the position information, the dynamic model parameters, the pneumatic parameters and the initial pitch angle of the light source. The invention solves the technical problem that the traditional chemical power carrier rocket trajectory design method in the prior art is not suitable for the trajectory design of the multi-light-source laser aircraft.
Description
Technical Field
The invention relates to the technical field of laser aircrafts, in particular to a method and a system for designing a trajectory of a multi-light-source laser aircraft.
Background
Laser propulsion is a novel aircraft propulsion technology, and can utilize laser beams emitted by a ground laser base station as a power source, so that the laser aircraft generates thrust to complete a satellite in-orbit emission task.
The laser aircraft launching scheme relates to the multidisciplinary design of structure, power, pneumatics, trajectory and the like. Related ballistic design methods of the domestic laser aircrafts are few, and are mostly designed for single-foundation laser stations, the existing design methods are mostly designed and integrated with various disciplines independently, coupling variables existing among the disciplines need to be iterated for multiple times, the time cost is large, and the overall design optimization is inconvenient to carry out.
Because the mode of generating thrust by laser propulsion is different from that of traditional chemical propulsion (a laser source provides a laser beam to act on a working medium to generate thrust), and the special requirement that the laser incidence angle (the included angle between the laser beam and the axis of the laser vehicle) is zero makes the traditional chemical carrier rocket ballistic design method not suitable for target laser vehicle ballistic design.
Disclosure of Invention
In view of the above, the present invention provides a method and a system for designing a trajectory of a multiple-light-source laser vehicle, so as to alleviate the technical problem that the conventional method for designing a trajectory of a chemical-powered launch vehicle in the prior art is not suitable for designing a trajectory of a multiple-light-source laser vehicle.
In a first aspect, an embodiment of the present invention provides a method for designing a trajectory of a multi-light source laser aircraft, where the method is applied to a target laser aircraft, the target laser aircraft is a laser aircraft that provides laser beams from at least two light sources, and the type of the target laser aircraft is a mylabo light ship configuration; the method comprises the following steps: determining structural model parameters of the target laser aircraft; determining position information of a light source providing a laser beam for the target laser vehicle, and determining an initial pitch angle of the target laser vehicle; calculating dynamic model parameters of the target laser aircraft in an acceleration process and calculating pneumatic parameters of the target laser aircraft in an air suction mode based on the structural model parameters; and calculating the trajectory of the target laser aircraft based on the position information of the light source, the dynamic model parameter, the pneumatic parameter and the initial pitch angle.
Further, a light source for providing a laser beam to the target laser vehicle includes: a first light source and a second light source; the first light source is a foundation laser station, and the second light source is a relay satellite; determining position information for a light source providing a laser beam for the target laser aerial vehicle, comprising: determining position information of the second light source based on: d α (t) ═ α0+ωαdt；
Wherein x is*、y*、z*For the position information of the relay satellite, x* max、y* max、z* maxA maximum coordinate value, alpha, projected on the geocentric equatorial coordinate system for the orbit of the relay satellitex、αy、αzRespectively the angle, alpha, traversed by the relay satellite when the maximum coordinate value is taken by projection0Angles, omega, between the line connecting the earth center and the line connecting the target laser vehicle and the earth center when the relay satellite starts to provide the laser beam for the target laser vehicleαIs the operating angular velocity, R, of the relay satelliteeThe radius of the earth, H is the height of the relay satellite from the ground, G is a universal gravitation constant, and M is the mass of the earth; d (t) is the connecting line between the relay satellite and the geocenter and the included angle of the connecting line between the laser aircraft and the geocenter, which changes along with time, when the relay satellite starts to provide the laser beam for the laser aircraft.
Further, the dynamic model parameters include impulse and thrust; calculating dynamic model parameters of the target laser aircraft in an acceleration process based on the structural model parameters, wherein the dynamic model parameters comprise: calculating the impulse to which the target laser aircraft is subjected in the air suction mode by the following formula: i ═
Calculating the thrust force of the target laser aircraft in the air suction mode according to the following formula: f ═ c (I)pp+Ism) PRF · cos δ; calculating the thrust force of the target laser aircraft in the rocket mode by the following formula: f ═ Isqmg; wherein, IppAnd IsmThe impulse components, t, of the laser pulse detonation and attenuation stages, respectively2DIs the time of propagation of the detonation wave to the momentum plate, RLSDIs the radius of the detonation wave, pLSDIs the impulse platePressure, w is the average perimeter of the skirt of the bare boat, w ═ pi (r)ciol+rceol)，rciolAnd rceolThe body radius of the skirt plate at the gas flow inlet and the body radius of the skirt plate at the gas flow outlet, tfIs the time of decay of the shock wave to ambient pressure, PaIs the ambient pressure; PRF is the maximum value of laser pulse frequency, delta is the inclination angle of the annular cover, and c is the laser frequency coefficient; q. q.smMass flow in rocket mode.
Further, the pneumatic parameters comprise air inlet channel ram resistance; calculating the aerodynamic parameters of the target laser aircraft in an air suction mode, wherein the calculation comprises the following steps: calculating the air inlet channel stamping resistance of the target laser aircraft in an air suction mode according to the following formula:
wherein D isRAMIn order to provide the air inlet channel with ram resistance,
ρin、uin、Ainrespectively the mass flow rate of the throat, the density of the throat, the flow velocity of the throat and the area of the throat, deltacone2Half cone angle of second-stage cone, rciolAnd rceolThe body radius of the apron at the gas flow inlet and the body radius of the apron at the gas flow outlet, respectively.
Further, performing ballistic computation on the target laser vehicle based on the position information of the light source, the dynamic model parameters, the pneumatic parameters, and the initial pitch angle, including: calculating a kinematic equation and a kinetic equation of the target laser aircraft in the launching process based on the dynamic model parameters, the pneumatic parameters and the initial pitch angle; and integrating the kinematic equation and the kinetic equation by adopting a fourth-order Runge Kutta method to obtain the speed information and the position coordinate of the target laser aircraft in the launching process.
Further, the acceleration process of the target laser vehicle includes: a first acceleration segment providing an acceleration segment of the laser beam to the first light source and a second acceleration segment providing an acceleration segment of the laser beam to the second light source; the method further comprises the following steps: determining the pitch angle change of the target laser aircraft in the transmitting process based on the position coordinates and the position information of the light source; wherein the change in pitch angle of the target laser vehicle in the first acceleration segment is determined by:
determining a change in pitch angle of the target laser vehicle over the second acceleration segment by:
and x, y and z are the position coordinates of the target laser aircraft.
In a second aspect, an embodiment of the present invention further provides a trajectory design system of a multi-light-source laser aircraft, which is applied to a target laser aircraft, where the target laser aircraft is a laser aircraft that provides laser beams from at least two light sources, and the model of the target laser aircraft is a mylabo light boat configuration; the method comprises the following steps: the device comprises a first determining module, a second determining module, a first calculating module and a second calculating module; the first determining module is used for determining the structural model parameters of the target laser aircraft; the second determination module is used for determining the position information of a light source providing a laser beam for the target laser aircraft and determining the initial pitch angle of the target laser aircraft; the first calculation module is used for calculating power model parameters of the target laser aircraft in an acceleration process and calculating pneumatic parameters of the target laser aircraft in an air suction mode based on the structural model parameters; the second calculation module is used for calculating the trajectory of the target laser aircraft based on the position information of the light source, the power model parameter, the pneumatic parameter and the initial pitch angle.
Further, the second calculation module is further configured to: calculating a kinematic equation and a kinetic equation of the target laser aircraft in the launching process based on the dynamic model parameters, the pneumatic parameters and the initial pitch angle; and integrating the kinematic equation and the kinetic equation by adopting a fourth-order Runge Kutta method to obtain the speed information and the position coordinate of the target laser aircraft in the launching process.
In a third aspect, an embodiment of the present invention further provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor implements the steps of the method according to the first aspect when executing the computer program.
In a fourth aspect, the present invention further provides a computer-readable medium having non-volatile program code executable by a processor, where the program code causes the processor to execute the method according to the first aspect.
The invention provides a method and a system for designing a trajectory of a multi-light-source laser aircraft, which are used for integrally designing the subjects such as structure, power, pneumatics and trajectories in the trajectory calculation of the launching process of the multi-light-source laser aircraft, so that the design iteration time is reduced, and the technical problem that the conventional method for designing the trajectory of the chemical power carrier rocket in the prior art is not suitable for designing the trajectory of the multi-light-source laser aircraft is solved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a flowchart of a method for designing a trajectory of a multi-light-source laser aircraft according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a Myrabo lightboat configuration provided by an embodiment of the present invention;
FIG. 3 is a schematic diagram of a laser aircraft launching process provided by an embodiment of the invention;
fig. 4 is a schematic diagram of a relay satellite position according to an embodiment of the present invention;
fig. 5 is a schematic diagram of a multiple light source laser aircraft trajectory design system according to an embodiment of the invention.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the accompanying drawings, and it should be understood that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The first embodiment is as follows:
fig. 1 is a flowchart of a method for designing a trajectory of a multi-light-source laser aircraft, which is applied to a target laser aircraft, where the target laser aircraft is a laser aircraft that provides laser beams from at least two light sources, and the target laser aircraft is a mylabo light ship configuration. As shown in fig. 1, the method specifically includes the following steps:
and S102, determining structural model parameters of the target laser aircraft.
Specifically, fig. 2 is a schematic diagram of a configuration of a Myrabo photoboat according to an embodiment of the present invention. As shown in fig. 2, the external structure includes: the thrust direction generated by the fairing 201, the annular sheath 202, the annular nozzle 203 and the parabolic reflector 204 is along the axis direction of the light boat. The internal structure includes: a payload (microsatellite) 205, a storage tank 206 filled with rocket mode working medium, a gas cylinder 207 filled with pressurized helium gas, and an electronic device, an attitude control device and other components 208.
Structural disciplines have a major impact on ballistic designMass properties and aerodynamic performance of a laser aircraft, the initial total mass of the laser aircraft comprising m0Containing the mass m of the payloadplMass m of propellantpStructural mass msAnd cowl mass mr. By the mass m of the payloadplAnd final on-track height requirement to determine propellant mass mpStructural mass msAnd cowl mass mr. Wherein m is0＝mpl+mp+ms+mr。
And step S104, determining the position information of a light source for providing the laser beam for the target laser aircraft, and determining the initial pitch angle of the target laser aircraft. In the embodiment of the invention, a double light source is taken as an example, and two light sources including a ground-based laser station and a relay satellite are used for providing laser beams for the laser aircraft.
And S106, calculating dynamic model parameters of the target laser aircraft in the acceleration process and calculating pneumatic parameters of the target laser aircraft in an air suction mode based on the structural model parameters.
And S108, performing ballistic computation on the target laser aircraft based on the position information of the light source, the dynamic model parameters, the pneumatic parameters and the initial pitch angle.
The invention provides a multi-light-source laser aircraft trajectory design method, which is used for integrally designing subjects such as structure, power, pneumatics and trajectories in the trajectory calculation of the launching process of a multi-light-source laser aircraft, so that the design iteration time is reduced, and the technical problem that the traditional chemical power carrier rocket trajectory design method in the prior art is not suitable for multi-light-source laser aircraft trajectory design is solved.
In the embodiment of the invention, the power system of the target laser aircraft adopts a pulse laser propulsion system, and the working modes are divided into an air suction mode and a rocket mode. In the air suction mode, the detonation wave generated by the laser breakdown gas is used for generating thrust. Specifically, the dynamic model parameters include impulse and thrust. The impulse of the target laser aircraft in the air suction mode is calculated by the following formula:
and calculating the thrust force of the target laser aircraft in the air suction mode by the following formula:
F＝c(Ipp+Ism)PRF·cosδ
the rocket mode is to generate thrust by using propellant working medium carried by a laser ablation laser aircraft. The thrust of the target laser aircraft in the rocket mode is calculated by the following formula:
F＝Is qm g
wherein, IppAnd IsmThe impulse components, t, of the laser pulse detonation and attenuation stages, respectively2DIs the time of propagation of the detonation wave to the momentum plate, RLSDIs the radius of the detonation wave, pLSDIs the pressure at the impulse plate, w is the average perimeter of the bare boat skirt, w ═ pi (r)ciol+rceol)，rciolAnd rceolThe body radius of the skirt plate at the gas flow inlet and the body radius of the skirt plate at the gas flow outlet, tfIs the time of decay of the shock wave to ambient pressure, PaIs the ambient pressure; PRF is the maximum value of laser pulse frequency, delta is the inclination angle of the annular cover, and c is the laser frequency coefficient; q. q.smMass flow in rocket mode.
Optionally, the aerodynamic parameter comprises air intake ram drag. In the embodiment of the invention, when the target laser aircraft works in the air suction mode, the resistance is mainly the air inlet channel stamping resistance. Specifically, the method provided by the embodiment of the invention calculates the air inlet channel ram resistance suffered by the target laser aircraft in the air suction mode through the following formula:
wherein D isRAMIn order to provide the punching resistance of the air inlet channel,
ρin、uin、Ainrespectively the mass flow rate of the throat, the density of the throat, the flow velocity of the throat and the area of the throat, deltacone2Half cone angle of second-stage cone, rciolAnd rceolThe body radius of the apron at the gas flow inlet and the body radius of the apron at the gas flow outlet, respectively.
Optionally, the light source for providing the laser beam to the target laser vehicle comprises: a first light source and a second light source; the first light source is a ground-based laser station and the second light source is a relay satellite.
In an embodiment of the invention, the ballistic calculation of the target laser vehicle includes a control calculation and a ballistic calculation of the launch process. The launching process of the target laser aircraft comprises four stages, specifically: an air-breathing mode acceleration section, a rocket mode acceleration section, a thrust-free gliding section and a horizontal acceleration section. Fig. 3 is a schematic diagram of a laser aircraft launching process provided by an embodiment of the invention. As shown in fig. 3, the transmission process includes: ground-based laser station 301, relay satellite 302, laser vehicle 303, laser beam 304, laser vehicle flight trajectory 305.
Optionally, the position of the second light source may also need to be calculated before ballistic calculations can be performed on the target laser vehicle. Specifically, step S104 includes:
determining position information of the second light source based on the following equation:
dα(t)＝α0+ωα dt
wherein x is*、y*、z*For relaying position information of satellites, x* max、y* max、z* maxMaximum coordinate value, alpha, projected on the geocentric equatorial coordinate system for the orbit of the relay satellitex、αy、αzRespectively the angle, alpha, traversed by the relay satellite when the projection takes the maximum coordinate value0The angle between the connecting line with the geocentric and the connecting line between the target laser aircraft and the geocentric when the relay satellite starts to provide the laser beam for the target laser aircraft, omegaaFor relaying the angular velocity of operation of the satellite, ReThe radius of the earth, H is the height of a relay satellite from the ground, G is a universal gravitation constant, and M is the mass of the earth; alpha (t) is an included angle between a connecting line of the relay satellite and the geocenter and a time-varying connecting line of the laser aircraft and the geocenter when the relay satellite starts to provide the laser beam for the laser aircraft. Fig. 4 is a schematic diagram of a relay satellite position provided according to an embodiment of the present invention. As shown in fig. 4, 401 is a relay satellite, and 402 is a target laser vehicle.
And then calculating the trajectory of the target laser aircraft. Specifically, step S108 further includes the following steps:
and step S1081, calculating a kinematic equation and a dynamic equation of the target laser aircraft in the launching process based on the dynamic model parameters, the pneumatic parameters and the initial pitch angle.
And step S1082, integrating the kinematic equation and the kinetic equation by adopting a fourth-order Runge Kutta method to obtain the speed information and the position coordinate of the target laser aircraft in the launching process.
Specifically, the kinematic equation of the centroid of the target laser vehicle can be determined by the following two equationsFormula (I) is determined, wherein xt、yt、ztIs the component, x, of the centre of mass position in the emitting coordinate system varying with time0＝0、y0＝3.65、z 00 is the initial position of the emission point, vx、vy、vzAs a component of velocity in the transmit coordinate system, ax、ay、azV is the component of the acceleration in the emission coordinate system and can be solved by the kinetic equationx0＝0、vy0＝0、vz0-0 is the initial speed.
Specifically, the centroid kinetic equation of the multi-light-source aircraft is shown as follows, wherein m is the total mass of the laser aircraft in the flight process, the mass is obtained from a flow-time curve, and the initial value is the takeoff mass m0(ii) a Left side of equation ax、ay、azIs the acceleration component of the laser aircraft in the inertial coordinate system of the launching point; px、Py、PzThe thrust under the inertial coordinate system of the launching point can be converted by coordinate transformation and the thrust under the elastic coordinate system without considering the eccentric deviation of the thrust. N is a radical ofx、Ny、NzThe laser aircraft is aerodynamic under an inertial coordinate system of a launching point, the aerodynamic force borne by the laser aircraft is mainly air resistance, is opposite to the flying speed direction in a speed coordinate system, and can be obtained through coordinate system transformation. g is a radical of formulax、gy、gzThe component of the acceleration of the earth gravity in an inertial coordinate system of the launching point is used; a isex、aey、aezProjecting the acceleration caused by the inertia force in an inertial coordinate system of a launching point, and calculating the acceleration according to the rotation angular velocity of the earth; in the formula acx、acy、aczThe acceleration caused by the Coriolis inertial force is projected in an inertial coordinate system of the launching point.
In an embodiment of the invention, the acceleration process of the target laser vehicle comprises: a first acceleration section providing an acceleration section of the laser beam for the first light source and a second acceleration section providing an acceleration section of the laser beam for the second light source; wherein the first acceleration section includes: an air-breathing mode acceleration section and a rocket mode acceleration section; the second acceleration section is a horizontal acceleration section. The pitch angle control equation of the target laser aircraft at each stage is as follows:
an air suction mode acceleration section: height H of emission point0Initial pitch angle at launch
The laser frequency coefficient c is determined by the optimization result, and the working mode of the laser aircraft is the air suction mode at the stage, the real-time thrust F is calculated according to the flight altitude and the Mach number1. The laser beam is directly supplied from a ground laser station, and the pitch angle is
Gradually reducing to ensure that the laser incidence angle is zero until the flying height reaches the height H of the working mode switching point1This stage is over; wherein the pitch angle is determined by the expression:
rocket mode acceleration section: to a height H1Then, the working mode of the aircraft is switched to a rocket mode to continue acceleration, and the thrust is F2And the fairing is thrown off at the appropriate height. The laser beam is directly provided by a foundation laser station, the pitch angle
The control is the same as the previous stage until the time T is set2Exhausted, and the phase ends.
A thrust-free gliding section: rate of change of pitch angle
And remain unchanged until the lower stage pitch angle condition is reduced. The laser station of the foundation no longer provides laser beams, the laser aircraft slides without thrust until the time TG of sliding without thrust is exhausted, and the stage is finished; the pitch angle change rate expression of the target laser aircraft at the stage is as follows:
a horizontal acceleration section: the ground-based laser station emits laser beams to the relay satellite, and the laser beams are shaped and repositioned by the relay lens system and then emitted to the laser aircraft to provide power for the laser aircraft to form a second light source with the thrust of F3And the pitch angle change ensures that the laser incidence angle is always zero, and the horizontal speed is increased until the laser enters the track. Specifically, the pitch angle change of the target laser aircraft in the horizontal acceleration section is determined by the following formula:
Therefore, the integrated design of the target laser aircraft trajectory can be completed, and the design variables mainly comprise the initial pitch angle
Laser frequency coefficient c and three-stage thrust F1、F2、F3Thrust-free sliding time TGHeight H of switching operation mode1And payload mass mp1。
In the above description, the relay satellite is used as the second light source, and the method provided by the embodiment of the present invention is described to realize the trajectory design of the multi-light-source laser aircraft, and for other multi-light-source schemes, for example, when the second light source is suspended on an airship at a certain height, and the second light source is on the ground at a certain distance from the first light source, etc., only the position coordinate (x) of the second light source needs to be changed*，y*，z*) According to the design method provided by the embodiment of the invention, the ballistic characteristics of the multi-light source laser aircraft can be obtained through integrated design. In the integrated design method provided by the embodiment of the invention, the laser aircraft, the first light source (foundation laser station) and the second light source (relay satellite) all have relative motion, and the method realizes real-time calculation of the positions of the laser aircraft, the first light source (foundation laser station) and the second light source (relay satellite).
The method for designing the trajectory of the multi-light-source laser aircraft provided by the embodiment of the invention has the following technical effects:
(1) the disciplines such as structure, power, pneumatics and trajectory are integrally designed, so that the design iteration time is reduced;
(2) aiming at the special characteristics of laser propulsion and multi-light source laser aircraft trajectory design, the pitch angle is controlled at each stage of the launching scheme, so that the laser incidence angle is always zero, the two structures of the laser aircraft are simplified, and the energy utilization rate is improved;
(3) the calculation of the position of the second light source (the relay satellite in the embodiment of the invention) ensures that the light source provides laser beams for the laser aircraft in real time;
(4) real-time calculation of the relative positions of the laser aircraft, the first light source (the ground-based laser station) and the second light source (the relay satellite in the embodiment of the invention) is realized.
(5) The method is suitable for the multi-light-source laser aircraft trajectory design which cannot be realized by the traditional chemical carrier rocket, and has strong applicability to the trajectory design of other multi-light-source laser aircraft except for the relay satellite serving as a second light source.
The second embodiment:
fig. 5 is a schematic diagram of a multi-light source laser aircraft trajectory design system according to an embodiment of the present invention, where the system is applied to a target laser aircraft, the target laser aircraft is a laser aircraft that provides laser beams from at least two light sources, and the target laser aircraft is a mylabo lightboat configuration. As shown in fig. 5, the system includes: a first determination module 10, a second determination module 20, a first calculation module 30 and a second calculation module 40.
Specifically, the first determination module 10 is configured to determine structural model parameters of the target laser aircraft.
And a second determination module 20 for determining position information of a light source providing a laser beam for the target laser flyer, and determining an initial pitch angle of the target laser flyer.
And the first calculation module 30 is used for calculating dynamic model parameters of the target laser aircraft in an acceleration process and calculating pneumatic parameters of the target laser aircraft in an air suction mode based on the structural model parameters.
And the second calculation module 40 is used for performing ballistic calculation on the target laser aircraft based on the position information of the light source, the power model parameters, the pneumatic parameters and the initial pitch angle.
The invention provides a multi-light-source laser aircraft trajectory design system, which is used for integrally designing subjects such as structure, power, pneumatics and trajectories in the trajectory calculation of the launching process of a multi-light-source laser aircraft, so that the design iteration time is reduced, and the technical problem that the traditional chemical power carrier rocket trajectory design method in the prior art is not suitable for multi-light-source laser aircraft trajectory design is solved.
Optionally, the first calculating module 30 is further configured to: calculating the impulse of the target laser aircraft in the air suction mode according to the following formula:
and calculating the thrust force of the target laser aircraft in the air suction mode by the following formula:
F＝c(Ipp+Ism)PRF·cosδ
the rocket mode is to use the propellant working medium carried by the laser ablation laser aircraft to generate thrust. The embodiment of the invention calculates the thrust of the target laser aircraft in the rocket mode by the following formula:
F＝Is qm g
wherein, IppAnd IsmThe impulse components, t, of the laser pulse detonation and attenuation stages, respectively2DIs the time of detonation wave propagation to the momentum plate, RLSDIs the radius of the detonation wave, pLSDIs the pressure at the impulse plate, w is the average perimeter of the bare boat skirt, w ═ pi (r)ciol+rceol)，rciolAnd rceolThe body radius of the skirt plate at the gas flow inlet and the body radius of the skirt plate at the gas flow outlet, tfIs the time of decay of the shock wave to ambient pressure, PaIs the ambient pressure; PRF is the maximum value of laser pulse frequency, delta is the inclination angle of the annular cover, and c is the laser frequency coefficient; q. q.smMass flow in rocket mode.
Optionally, the first calculating module 30 is further configured to: calculating the air inlet channel stamping resistance of the target laser aircraft in the air suction mode according to the following formula:
wherein D isRAMIn order to provide the punching resistance of the air inlet channel,
ρin、uin、Ainrespectively the mass flow rate of the throat, the density of the throat, the flow velocity of the throat and the area of the throat, deltacone2Half-cone angle of the second-order cone, rciolAnd rceolThe body radius of the skirt at the gas flow inlet and the body radius of the skirt at the gas flow outlet, respectively.
Optionally, the second calculating module 40 is further configured to: calculating a kinematic equation and a kinetic equation of the target laser aircraft in the launching process based on the dynamic model parameters, the pneumatic parameters and the initial pitch angle; and integrating the kinematic equation and the kinetic equation by adopting a fourth-order Runge Kutta method to obtain the speed information and the position coordinate of the target laser aircraft in the launching process.
An embodiment of the present invention further provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the steps of the method in the first embodiment are implemented.
Embodiments of the present invention further provide a computer-readable medium having a non-volatile program code executable by a processor, where the program code causes the processor to execute the method in the first embodiment.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (6)
1. A trajectory design method of a multi-light-source laser aircraft is applied to a target laser aircraft, wherein the target laser aircraft is a laser aircraft which provides laser beams for at least two light sources, and the model of the target laser aircraft is a Myrabo light boat configuration; it is characterized by comprising:
determining structural model parameters of the target laser aircraft;
determining position information of a light source providing a laser beam for the target laser vehicle, and determining an initial pitch angle of the target laser vehicle;
calculating dynamic model parameters of the target laser aircraft in an acceleration process and calculating pneumatic parameters of the target laser aircraft in an air suction mode based on the structural model parameters;
performing ballistic computation on the target laser vehicle based on the position information of the light source, the dynamic model parameters, the pneumatic parameters, and the initial pitch angle;
the light source for providing the target laser vehicle with a laser beam comprises: a first light source and a second light source; the first light source is a foundation laser station, and the second light source is a relay satellite; determining positional information for a light source providing a laser beam for the target laser aerial vehicle, comprising:
determining position information of the second light source based on:
wherein x is*、y*、z*For the position information of the relay satellite, x* max、y* max、z* maxA maximum coordinate value, alpha, projected on the geocentric equatorial coordinate system for the orbit of the relay satellitex、αy、αzRespectively the angle, alpha, traversed by the relay satellite when the maximum coordinate value is taken by projection0Angles, omega, between the line connecting the earth center and the line connecting the target laser vehicle and the earth center when the relay satellite starts to provide the laser beam for the target laser vehicleαIs the operating angular velocity, R, of the relay satelliteeThe radius of the earth, H is the height of the relay satellite from the ground, G is a universal gravitation constant, and M is the mass of the earth; alpha (t) is an included angle between a connecting line of the relay satellite and the geocenter and an included angle between the connecting line of the laser aircraft and the geocenter, which changes along with time, when the relay satellite starts to provide a laser beam for the laser aircraft;
performing ballistic computation on the target laser vehicle based on the position information of the light source, the dynamic model parameters, the pneumatic parameters, and the initial pitch angle, including:
calculating a kinematic equation and a kinetic equation of the target laser aircraft in the launching process based on the dynamic model parameters, the pneumatic parameters and the initial pitch angle;
integrating the kinematic equation and the kinetic equation by adopting a fourth-order Runge Kutta method to obtain speed information and position coordinates of the target laser aircraft in the launching process;
the acceleration process of the target laser aircraft comprises the following steps: a first acceleration section providing an acceleration section of the laser beam for the first light source and a second acceleration section providing an acceleration section of the laser beam for the second light source; the method further comprises the following steps:
determining the pitch angle change of the target laser aircraft in the transmitting process based on the position coordinates and the position information of the light source; wherein the content of the first and second substances,
determining a pitch angle of the target laser vehicle in the first acceleration segment by:
determining a pitch angle of the target laser vehicle in the second acceleration segment by:
phi is the pitch angle of the target laser aircraft, and x, y and z are the position coordinates of the target laser aircraft.
2. The method of claim 1, wherein the dynamic model parameters include impulse and thrust; calculating dynamic model parameters of the target laser aircraft in an acceleration process based on the structural model parameters, wherein the dynamic model parameters comprise:
calculating the impulse of the target laser aircraft in the air suction mode according to the following formula:
calculating the thrust force of the target laser aircraft in the air suction mode according to the following formula:
F=c(I pp +I sm )PRF·cosδ
calculating the thrust force of the target laser aircraft in the rocket mode by the following formula:
F=I s q m g
wherein, IppAnd IsmThe impulse components, t, of the laser pulse detonation and attenuation stages, respectively2DIs the time of propagation of the detonation wave to the momentum plate, RLSDIs the radius of the detonation wave, pLSDIs the pressure at the impulse plate, w is the average circumference of the smooth boat skirt, w = π (r)ciol+rceol），rciolAnd rceolThe body radius of the skirt plate at the gas flow inlet and the body radius of the skirt plate at the gas flow outlet, t f Is the time of decay of the shock wave to ambient pressure, PaIs the ambient pressure; PRF is the maximum value of laser pulse frequency, delta is the inclination angle of the annular cover, and c is the laser frequency coefficient; q. q ofmMass flow in rocket mode.
3. The method of claim 1, wherein the aerodynamic parameter comprises a port ram drag; calculating the aerodynamic parameters of the target laser aircraft in an air suction mode, wherein the calculation comprises the following steps:
calculating the air inlet channel stamping resistance of the target laser aircraft in the air suction mode according to the following formula:
wherein D isRAMIn order to provide the air inlet channel with ram resistance,
、ρin、uin、Ainletrespectively the mass flow rate of the throat, the density of the throat, the flow velocity of the throat and the area of the throat, deltacone2Half cone angle of second-stage cone, rciolAnd rceolThe body radius of the apron at the gas flow inlet and the body radius of the apron at the gas flow outlet, respectively.
4. A multi-light source laser aircraft trajectory design system is applied to a target laser aircraft, wherein the target laser aircraft is a laser aircraft which provides laser beams for at least two light sources, and the model of the target laser aircraft is a Myrabo light boat configuration; it is characterized by comprising: the device comprises a first determining module, a second determining module, a first calculating module and a second calculating module; wherein the content of the first and second substances,
the first determination module is used for determining the structural model parameters of the target laser aircraft;
the second determination module is used for determining the position information of a light source providing a laser beam for the target laser aircraft and determining the initial pitch angle of the target laser aircraft;
the first calculation module is used for calculating power model parameters of the target laser aircraft in an acceleration process and calculating pneumatic parameters of the target laser aircraft in an air suction mode based on the structural model parameters;
the second calculation module is used for performing ballistic calculation on the target laser aircraft based on the position information of the light source, the dynamic model parameter, the pneumatic parameter and the initial pitch angle;
the second computing module is further to:
calculating a kinematic equation and a dynamic equation of the target laser aircraft in the launching process based on the dynamic model parameters, the pneumatic parameters and the initial pitch angle;
integrating the kinematic equation and the kinetic equation by adopting a fourth-order Runge Kutta method to obtain speed information and position coordinates of the target laser aircraft in the launching process;
the light source for providing a laser beam to the target laser vehicle comprises: a first light source and a second light source; the first light source is a foundation laser station, and the second light source is a relay satellite; the second determining module is further configured to:
determining position information of the second light source based on:
wherein x is*、y*、z*For the position information of the relay satellite, x* max、y* max、z* maxA maximum coordinate value, alpha, projected on the geocentric equatorial coordinate system for the orbit of the relay satellitex、αy、αzRespectively the angle, alpha, traversed by the relay satellite when the maximum coordinate value is taken by projection0Angles, omega, between the line connecting the earth center and the line connecting the target laser vehicle and the earth center when the relay satellite starts to provide the laser beam for the target laser vehicleαIs the operating angular velocity, R, of the relay satelliteeThe radius of the earth, H is the height of the relay satellite from the ground, G is a universal gravitation constant, and M is the mass of the earth; alpha (t) is an included angle between a connecting line of the relay satellite and the geocenter and a connecting line of the laser aircraft and the geocenter, which changes along with time, when the relay satellite starts to provide a laser beam for the laser aircraft;
the acceleration process of the target laser aircraft comprises: a first acceleration section providing an acceleration section of the laser beam for the first light source and a second acceleration section providing an acceleration section of the laser beam for the second light source; the second computing module is further configured to:
determining the pitch angle change of the target laser aircraft in the transmitting process based on the position coordinates and the position information of the light source; wherein the content of the first and second substances,
determining the pitch angle of the target laser vehicle in the first acceleration section by the following formula:
determining a pitch angle of the target laser vehicle in the second acceleration segment by:
phi is the pitch angle of the target laser aircraft, and x, y and z are the position coordinates of the target laser aircraft.
5. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the steps of the method of any of the preceding claims 1 to 3 are implemented when the computer program is executed by the processor.
6. A computer-readable medium having non-volatile program code executable by a processor, wherein the program code causes the processor to perform the method of any of claims 1-3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110709788.3A CN113342049B (en) | 2021-06-25 | 2021-06-25 | Multi-light source laser aircraft trajectory design method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110709788.3A CN113342049B (en) | 2021-06-25 | 2021-06-25 | Multi-light source laser aircraft trajectory design method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113342049A CN113342049A (en) | 2021-09-03 |
| CN113342049B true CN113342049B (en) | 2022-07-19 |
Family
ID=77478741
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110709788.3A Active CN113342049B (en) | 2021-06-25 | 2021-06-25 | Multi-light source laser aircraft trajectory design method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113342049B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5796474A (en) * | 1996-06-21 | 1998-08-18 | Thermotrex Corporation | Projectile tracking system |
| CN106021628A (en) * | 2015-07-03 | 2016-10-12 | 中国运载火箭技术研究院 | Vertical returning trajectory design method for carrier rocket |
| CN109398762A (en) * | 2018-10-17 | 2019-03-01 | 湖北航天技术研究院总体设计所 | A kind of solid-rocket enters rail ballistic design method based on elliptical transfer orbit |
| CN109583041A (en) * | 2018-11-07 | 2019-04-05 | 湖北航天技术研究院总体设计所 | A kind of Craft Orbit design method |
| CN112464365A (en) * | 2020-11-19 | 2021-03-09 | 北京航空航天大学 | Parameter design method and device of aircraft |
-
2021
- 2021-06-25 CN CN202110709788.3A patent/CN113342049B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5796474A (en) * | 1996-06-21 | 1998-08-18 | Thermotrex Corporation | Projectile tracking system |
| CN106021628A (en) * | 2015-07-03 | 2016-10-12 | 中国运载火箭技术研究院 | Vertical returning trajectory design method for carrier rocket |
| CN109398762A (en) * | 2018-10-17 | 2019-03-01 | 湖北航天技术研究院总体设计所 | A kind of solid-rocket enters rail ballistic design method based on elliptical transfer orbit |
| CN109583041A (en) * | 2018-11-07 | 2019-04-05 | 湖北航天技术研究院总体设计所 | A kind of Craft Orbit design method |
| CN112464365A (en) * | 2020-11-19 | 2021-03-09 | 北京航空航天大学 | Parameter design method and device of aircraft |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113342049A (en) | 2021-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Smart et al. | Orbital delivery of small payloads using hypersonic airbreathing propulsion | |
| Josyula et al. | Applications of a counterflow drag reduction technique in high-speed systems | |
| Riddle et al. | Genetic-algorithm optimization of liquid-propellant missile systems | |
| CN112182772A (en) | Rocket propulsion control method, device and storage medium | |
| Tsuchiya et al. | Optimal conceptual design of two-stage reusable rocket vehicles including trajectory optimization | |
| Sarigul-Klijn et al. | Air-launching earth to orbit: Effects of launch conditions and vehicle aerodynamics | |
| Tsuchiya et al. | Optimal design of two-stage-to-orbit space planes with airbreathing engines | |
| CN113342049B (en) | Multi-light source laser aircraft trajectory design method and system | |
| Edquist et al. | Aerothermodynamic design of the Mars Science Laboratory backshell and parachute cone | |
| CN113428382B (en) | Laser aircraft launching system and method with relay satellite | |
| CN113741551A (en) | Overall process trajectory optimization method and device based on proxy model | |
| CN114279272A (en) | Combined mode laser propulsion system for tiny load emission | |
| DwyerCianciolo et al. | Overview of the nasa entry, descent and landing systems analysis exploration feed-forward study | |
| Johnston et al. | Studies of engine-airframe integration on hypersonic aircraft | |
| Lee | Mission and trajectory optimization of the air-launching rocket system using MDO techniques | |
| Lewis | Sharp leading edge hypersonic vehicles in the air and beyond | |
| Davis et al. | Review of laser lightcraft propulsion system | |
| Rodi et al. | High lift-to-drag ratio waveriders for missions in the Martian atmosphere | |
| Guo et al. | Optimization on interceptor trajectory for near space hypersonic target | |
| Boltz | Optimal ascent trajectory for efficient air launch into orbit | |
| Board | Why and Whither Hypersonics Research in the US Air Force | |
| Myrabo | Propulsion Systems Integration for a ‘Tractor Beam’Mercury Lightcraft: Liftoff Engine | |
| Ordoneza et al. | Constructal Law Conference (CLC 2023) Freedom, Design and Evolution Turin (Italy), 21–22 September | |
| Bonnefond et al. | Study of a generic SSTO vehicle using airbreathing propulsion | |
| KIRKHAM et al. | Joint USAF/NASA hypersonic research aircraft study |
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 |