CN113704676A - Multi-physical-field-circuit cooperative computing method for electromagnetic rail gun - Google Patents
Multi-physical-field-circuit cooperative computing method for electromagnetic rail gun Download PDFInfo
- Publication number
- CN113704676A CN113704676A CN202111047061.XA CN202111047061A CN113704676A CN 113704676 A CN113704676 A CN 113704676A CN 202111047061 A CN202111047061 A CN 202111047061A CN 113704676 A CN113704676 A CN 113704676A
- Authority
- CN
- China
- Prior art keywords
- field
- current
- electromagnetic
- impedance
- physical
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Train Traffic Observation, Control, And Security (AREA)
Abstract
The invention belongs to the technical field of electronic information, and relates to a multi-physical-field-circuit collaborative computing method for an electromagnetic rail gun. The method comprises the steps of firstly, equally dividing the launching process of the whole electromagnetic orbit gun into N time slices based on set calculation time, calculating the end voltage V and the current I of impedance in one time slice, obtaining the results of the current, the magnetic field, the force field and the thermal field of the orbit gun and then calculating to obtain the impedance: a resistor R and an inductor L; and finally judging whether the transmitting process is finished. The invention has the advantages that: 1) the impedance of the rail-armature is determined by the fields φ and A and is updated at each time slice. 2) The multi-physical fields are determined in real time by the impedance terminal voltage and current. 3) Under the premise of coupling among multiple physical fields, the launching process of the electromagnetic rail gun can be described more accurately by performing cooperative calculation with a pulse power supply. 4) This method can be transplanted to the study of other electromagnetic emitting devices.
Description
Technical Field
The invention belongs to the technical field of electronic information, and relates to a multi-physical-field-circuit collaborative computing method.
Background
The electromagnetic rail gun is a new-concept weapon with wide application prospect and has incomparable advantages compared with the traditional weapon which takes gunpowder as the emission energy. The launching process of the electromagnetic rail gun is a moving process of the armature under the action of a pulse power supply, a phenomenon that various physical fields coexist can be generated in the process, the physical fields are mutually coupled, and meanwhile, the moving of the armature changes the impedance value of the armature-guide rail as a power supply load. This process is accompanied by interactions between the multiple physical fields and interactions with the pulsed power supply. The launching of the electromagnetic rail cannon involves extremely complex circuits (pulse power supply) and the combined action of various physical processes such as electromagnetism, force, movement and heat, and some complex physical phenomena are difficult to test.
The simulation calculation is an important and indispensable means for researching the launching process of the electromagnetic orbit cannon. At present, the research of numerical simulation at home and abroad almost relates to all aspects of an electromagnetic rail gun launching device, and mainly researches a pulse power supply, an electromagnetic field, a force field and a thermal field and the mutual influence process of the pulse power supply, the electromagnetic field, the force field and the thermal field. The pulse power supply is generally modeled by an equivalent circuit, and numerical simulation of multiple physical fields such as an electromagnetic field, a thermal field, a force field and the like mainly adopts a finite element algorithm or a method of combining finite elements with boundary elements and the like. In the conventional simulation method, the pulse power supply only acts as an excitation armature motion, and the influence of the multi-physical field accompanying the armature motion on the pulse power supply is not considered, so that certain deviation is caused to the result in the simulation calculation. In short, the track and the armature serve as a load of the pulse power supply, and the load value is changed along with the movement of the armature; therefore, the waveform of the current passing through the load also changes to some extent, and ultimately affects physical quantities such as the speed of the armature.
Disclosure of Invention
(I) technical problems to be solved by the invention
The invention aims to solve the technical problem of providing a multi-physical-field-circuit cooperative computing method, which solves the problem that the influence of the impedance values of an armature and a track changing along with time on a pulse power supply is neglected in the existing simulation method, and can more accurately realize the simulation of the transient process of an electromagnetic rail gun.
(II) clear and complete technical scheme of the invention
A multi-physical field-circuit cooperative computing method for an electromagnetic rail gun comprises the following steps:
firstly, time division is carried out;
equally dividing the launching process of the whole electromagnetic rail gun into N time slices based on the set calculation time;
secondly, calculating the end voltage V and the current I of the impedance in a time slice;
at the initial moment, solving the equivalent circuit of the electromagnetic rail gun by a Newton iteration method so as to obtain the end voltage V and the current I of the impedance;
acquiring the current, magnetic field, force field and thermal field multi-physical field results of the rail gun;
the current and magnetic field distribution of the rail gun can be respectively obtained through formulas (4) and (1);
wherein, B is magnetic induction intensity, mu is magnetic conductivity, and J is current density; e is the electric field intensity, A is the vector potential, and phi is the electric potential.
Based on the current and magnetic field distribution, a force field F can be obtained by equation (5);
based on the current and the force, a thermal field can be obtained by the formula (6);
wherein rho, Cp and k are respectively the specific heat of the material, the constant-pressure heat capacity and the heat conduction coefficient of the material, T is the temperature and sigma is the electrical conductivity; qfIs frictional heat;
fourthly, calculating to obtain impedance: a resistor R and an inductor L;
according to the joule heat quantity information, the resistance R can be obtained by the formulas (7) and (8);
Qe＝i2R (8)
wherein i is a current;
according to the magnetic energy information, the inductance L can be obtained by the formulas (12) and (14);
where H is the magnetic field strength.
Fifthly, judging whether the transmitting process is finished or not;
if the calculated armature movement displacement is larger than or equal to the length of the track, judging that the launching process is finished; if the time slice is less than the track length, returning to the step two, and calculating the next time slice; if the displacement of the armature movement has not reached the firing track length after cycling for N time slices, firing ends.
The invention has the advantages that: compared with the prior art, the invention has the advantages that: 1) the impedance of the rail-armature is determined by the fields φ and A and is updated at each time slice. 2) The multi-physical fields are determined in real time by the impedance terminal voltage and current. 3) Under the premise of coupling among multiple physical fields, the launching process of the electromagnetic rail gun can be described more accurately by performing cooperative calculation with a pulse power supply. 4) This method can be transplanted to the study of other electromagnetic emitting devices.
Drawings
FIG. 1 illustrates a multi-physics field-circuit cooperative computing. During the process of launching the described rail gun, data exchange between multiple physical fields and a circuit is realized, namely the movement of an armature needs the circuit to drive, and the load impedance of the circuit is changed when the armature moves.
Fig. 2 is an equivalent circuit structure of the electromagnetic orbital cannon. The impedance of the track-armature is described by a variable resistance and inductance.
FIG. 3 is a flow chart of multi-physics field-circuit cooperative computing.
Detailed Description
Referring to fig. 1, for an electromagnetic orbital cannon, it is generally divided into two parts: pulsed power supply and rail-armature. In simulation, it is called circuit analysis module and multi-physical field analysis module, and there is exchange between these two modules in the transmission process, so that such a process is field-circuit coupled. That is to say, the field-road cooperative computation includes two main processes, that is, within a time slice, the calculation from road to field and the calculation from field to road need to be completed. The fields here are mainly electric, magnetic, thermal and force fields.
Calculation of (one) way to field
As can be appreciated from the above description, the energy generated by the pulse circuit drives the armature to move, causing a change in the impedance of the switch-in, which value will be determined by the multiple physical fields. From a circuit point of view, an equivalent circuit structure thereof can be constructed. Therefore, the current I passing through the electromagnetic rail cannon and the corresponding end voltage V can be obtained by solving the circuit discharge equation by adopting a Newton iteration method.
And (4) carrying out mesh subdivision on the electromagnetic rail gun by adopting a finite element method, and obtaining corresponding field distribution through numerical calculation.
For magnetic fields, the maxwell's magnetic field rotation equation (here, the displacement current term is not considered) is considered for representation:
the distribution of the magnetic induction intensity B can be obtained by solving the expression, and then the distribution of the magnetic field intensity H is obtained.
Since the divergence of the magnetic induction B is constant 0, B can be expressed as the curl of a vector function a, i.e. as
The formula combines Faraday's law of electromagnetic induction to obtain
This means that the physical quantity in the above formula in parentheses is unrotated and a scalar function can be used
The final electric field strength can be expressed as:
wherein
The value of (b) is the front end pressure V.
The method for solving the electric and magnetic fields firstly solves the electric potential phi and the vector potential A and then obtains the electromagnetic distribution. This method is also called
The method is carried out.
Under the condition of magnetic field and current, the electromagnetic force can be obtained by Lorentz force formula, and the corresponding equation is as follows
The motion information of the armature, such as displacement, velocity and acceleration, is obtained by newton's second law.
Likewise, for the temperature distribution, the heat transfer equation is solved. The heat source includes Joule heat and friction heat, and the corresponding equation is as follows
Wherein QfThe friction heat is related to the pressure of the armature on the rail, the moving speed of the armature and the friction coefficient of a contact surface.
(II) field-to-road calculation
Not only does a large amount of joule heat be generated during the firing of an electromagnetic orbital cannon, but a certain amount of magnetic energy is stored, which describes the impedance of the orbital-armature as a combination of resistance and inductance. The calculation from the field to the circuit is therefore mainly to obtain the resistance R and the inductance L of the track-armature.
For the resistance, according to joule heat generated by the armature during the movement:
similarly, the calculation formula of joule heat in the circuit theory is as follows
Qe＝i2R (8)
Thus, it is possible to provide
Wherein the current density J can be represented by the field quantities (φ and A) as follows
The amount of the path (resistance R) is then represented by the field quantities (phi and a).
For the expression of the inductance, the idea of obtaining the expression of the resistance is the same. Except that the resistance is based on joule heat and the inductance is based on magnetic energy. In the electromagnetic module, magnetic field distribution can be obtained, and then magnetic energy can be determined, and the expression is as follows
The magnetic field strength H can be replaced by a magnetic induction B, i.e. H ═ B/μ.
The magnetic induction B can be expressed by a formula
The magnetic induction B is then characterized by a current density, which is in accordance with the above equivalent resistance.
Likewise, magnetic energy may in turn be expressed as
The inductance L can therefore be written as the following expression
The above formula current density J is replaced by formula (10), and the inductance L is also described by phi and a.
Compared with the prior art, the invention has the advantages that: 1) the impedance of the rail-armature is determined by the fields φ and A and is updated at each time slice. 2) The multi-physical fields are determined in real time by the impedance terminal voltage and current. 3) Under the premise of coupling among multiple physical fields, the launching process of the electromagnetic rail gun can be described more accurately by performing cooperative calculation with a pulse power supply. 4) This method can be transplanted to the study of other electromagnetic emitting devices.
Referring to fig. 3, the detailed calculation flow is as follows:
firstly, time division is carried out;
the launching process of the electromagnetic orbit gun is a transient process, time division is needed for carrying out numerical analysis on the process, the whole process is equally divided into N time slices, and field-path cooperative calculation is completed in each time slice.
Second, at the initial moment, the armature has not moved yet, and the track-armature impedances R and L at this moment are R0 and L0 (these two values are related to the track-armature parameters), so that the circuit shown in fig. 2 is solved by newton's iteration, so as to obtain the terminal voltage V and the current I of the impedance.
Thirdly, based on the obtained V and I, the current and the magnetic field distribution of the rail gun can be respectively obtained through formulas (4) and (1); based on the current and magnetic field distribution, a force field can be obtained by equation (5); based on the current and force, a thermal field can be obtained by equation (6).
Fourthly, according to the obtained joule heat quantity information, the resistance R can be obtained by the formulas (7) and (8), and the corresponding expression is the formula (9); from the obtained magnetic field distribution, the inductance L can be obtained from the expressions (12) and (14), and the corresponding expression is (15).
Fifthly, judging whether the transmitting process is finished or not;
based on the length of the track, as the length of the launching track is known, when the calculated armature movement displacement is more than or equal to the length of the launching track, the launching process is judged to be finished; if the length of the transmitting track is less than the length of the transmitting track, calculating the next time slice, and returning to the step two; if the displacement of the armature movement still does not reach the length of the firing track after cycling for N time slices, firing ends.
Claims (2)
1. A multi-physical field-circuit cooperative computing method for an electromagnetic rail gun is characterized by comprising the following steps:
firstly, time division is carried out;
equally dividing the launching process of the whole electromagnetic rail gun into N time slices based on the set calculation time;
secondly, calculating the end voltage V and the current I of the impedance in a time slice;
at the initial moment, solving an equivalent circuit of the electromagnetic rail gun to obtain end voltage V and current I of impedance;
acquiring the current, magnetic field, force field and thermal field multi-physical field results of the rail gun;
the current and magnetic field distribution of the rail gun can be respectively obtained through formulas (4) and (1);
wherein B is magnetic induction intensity, mu is magnetic conductivity, J is current density, E is electric field intensity, A is vector potential,
is an electric potential;
based on the current and magnetic field distribution, a force field F can be obtained by equation (5);
based on the current and the force, a thermal field can be obtained by the formula (6);
wherein rho, Cp and k are respectively the specific heat of the material, the constant-pressure heat capacity and the heat conduction coefficient of the material, T is the temperature and sigma is the electrical conductivity; qfIs frictional heat;
fourthly, calculating to obtain impedance: a resistor R and an inductor L;
according to the information of the Joule heat Qe, the resistance R can be obtained by the formulas (7) and (8);
Qe＝i2R (8)
wherein i is a current;
according to the magnetic energy Qm information, the inductance L can be obtained by the formulas (12) and (14);
wherein H is the magnetic field strength;
fifthly, judging whether the transmitting process is finished or not;
if the calculated armature movement displacement is more than or equal to the length of the launching track, judging that the launching process is finished; if the length of the launching track is less than the length of the launching track, returning to the step two, and calculating the next time slice; if the displacement of the armature movement has not reached the firing track length after cycling for N time slices, firing ends.
2. The method for the cooperative calculation of the multiple physical fields and the circuit of the electromagnetic orbital cannon according to claim 1, wherein in the second step, the equivalent circuit of the electromagnetic orbital cannon is solved by a Newton iteration method, so that the end voltage V and the current I of the impedance are obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111047061.XA CN113704676B (en) | 2021-09-06 | 2021-09-06 | Multi-physical-field-circuit cooperative computing method for electromagnetic rail gun |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111047061.XA CN113704676B (en) | 2021-09-06 | 2021-09-06 | Multi-physical-field-circuit cooperative computing method for electromagnetic rail gun |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113704676A true CN113704676A (en) | 2021-11-26 |
| CN113704676B CN113704676B (en) | 2022-11-01 |
Family
ID=78659102
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111047061.XA Active CN113704676B (en) | 2021-09-06 | 2021-09-06 | Multi-physical-field-circuit cooperative computing method for electromagnetic rail gun |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113704676B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105183993A (en) * | 2015-09-09 | 2015-12-23 | 哈尔滨工业大学 | Integrated simulation platform and method for electromagnetic rail gun |
-
2021
- 2021-09-06 CN CN202111047061.XA patent/CN113704676B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105183993A (en) * | 2015-09-09 | 2015-12-23 | 哈尔滨工业大学 | Integrated simulation platform and method for electromagnetic rail gun |
Non-Patent Citations (3)
| Title |
|---|
| 刘洋天: "数值解析轨道电磁炮发射速度与相关参数的关系", 《中国战略新兴产业》 * |
| 李刚等: "基于ANSYS的电磁轨道炮C型固体电枢有限元分析", 《现代机械》 * |
| 林庆华等: "电磁轨道炮瞬态磁场测量与数值模拟", 《兵工学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113704676B (en) | 2022-11-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | A microstrip circuit tool kit app with FDTD analysis including lumped elements | |
| CN107529612A (en) | A kind of electromagnetic railgun pulse power source control method and apparatus | |
| Barmada et al. | Analysis of the performance of a multi-stage pulsed linear induction launcher | |
| CN113704676B (en) | Multi-physical-field-circuit cooperative computing method for electromagnetic rail gun | |
| Deng et al. | Study on the influence of armature on the efficiency of reluctance accelerator | |
| Yang et al. | An improved genetic algorithm for multiobjective optimization of helical coil electromagnetic launchers | |
| Lu et al. | Advanced mathematical calculation model of single-stage RCG | |
| Bissal et al. | On the design of ultra-fast electromechanical actuators: A comprehensive multi-physical simulation model | |
| Niu et al. | Research on self-consistent control strategy of multistage synchronous induction coil launcher | |
| Fan et al. | Research on energy conversion efficiency of the reconfigurable reconnection electromagnetic launcher | |
| Preyer | Transient behavior of simple RC circuits | |
| Waindok et al. | Mutual inductances in a mathematical model of the three-stage reluctance accelerator | |
| Haghmaram et al. | Transient modeling of multiparallel tubular linear induction motors | |
| Ruddy | High force density linear permanent magnet motors:" electromagnetic muscle actuators" | |
| Cho et al. | Optimal design of electromagnet for Maglev vehicles using hybrid optimization algorithm | |
| Djuric et al. | The variable geometry inductors for the energy harvesting applications | |
| Keshtkar et al. | Inductance gradient variation with time and armature sliding along the rails | |
| Keshtkar et al. | Comparison of inductance gradient and electromagnetic force in two types of railguns with two projectiles by finite-element method | |
| CN104182589B (en) | A kind of modeling method of electromagnetic rail launching device | |
| Patel et al. | Effect of design parameter variation on force acting on projectile of reluctance coilgun type electromagnetic launcher | |
| CN117077473B (en) | Two-dimensional axisymmetric electromagnetic thermal multi-field coupling calculation method based on transmission line method | |
| Ruddy | High Force Density Linear Permanent Magnet Motors | |
| CN116362098A (en) | Close coupling implementation method of 'path-field-particle' collaborative simulation | |
| Yaxiu et al. | A Novel Algorithm Based on Two-step CN-FDTD | |
| Rajput et al. | Modeling and Simulation of High Energy Capacitor Bank based Electromagnetic Railgun |
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 |