CN112733347A - Conductor-medium composite target and environment electromagnetic scattering rapid calculation method - Google Patents
Conductor-medium composite target and environment electromagnetic scattering rapid calculation method Download PDFInfo
- Publication number
- CN112733347A CN112733347A CN202011623614.7A CN202011623614A CN112733347A CN 112733347 A CN112733347 A CN 112733347A CN 202011623614 A CN202011623614 A CN 202011623614A CN 112733347 A CN112733347 A CN 112733347A
- Authority
- CN
- China
- Prior art keywords
- environment
- conductor
- target
- composite target
- medium composite
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a method for quickly calculating electromagnetic scattering of a conductor-medium composite target and an environment, which comprises the steps of subdividing the target and the environment to obtain an ultra-low-altitude conductor-medium composite target and environment composite scattering model, attributing the conductor-medium composite target to a JMCFIE-MSIE area, attributing the environment to a PO area, obtaining surface electricity and magnetic current of the PO area to obtain environmental surface induction electricity and magnetic field, and taking the environmental surface induction electricity and magnetic field and radar waves as a conductor-medium composite target irradiation source to obtain surface electricity and magnetic current to obtain surface induction electricity and magnetic field of the conductor-medium composite target area; taking the surface induction electricity, the magnetic field and the radar wave of the conductor-medium composite target area as an environment irradiation source to obtain environment surface electricity and magnetic current and obtain environment surface induction electricity and the magnetic field; and (5) performing repeated iteration, and stopping when the iteration error is smaller than a specific value to obtain a final field. The method can be used for quickly and accurately calculating the electromagnetic scattering property of the ultra-low-altitude conductor-dielectric composite target when the ultra-low-altitude conductor-dielectric composite target is positioned above the environment.
Description
Technical Field
The invention relates to the technical field of electromagnetic simulation, in particular to a method for quickly calculating electromagnetic scattering of a conductor-medium composite target and environment.
Background
In modern war, ultra-low altitude penetration has become the first choice of attack means for enemy air attack weapons. The enemy ultra-low altitude weapon mainly comprises: the air-raid weapons fly by depending on ultra-low altitude ground-attached and sea-swept, and avoid direct detection of a radar by utilizing the curvature of the earth surface and the fluctuation of the earth surface, and hide own echo signals by means of strong environmental clutter on the earth surface and multipath effects generated by multiple coupling of a target and the lower environment so as to avoid detection of air-defense weapons, and are one of main threats faced by an air-defense system. The actual ultra-low-altitude penetration weapon comprises various materials, wherein the conductor-dielectric material phase is compounded most typically, so that the ultra-low-altitude penetration weapon has extremely important military significance and practical value for the research of ultra-low-altitude conductor-dielectric compound target and environmental scattering modeling and compound scattering characteristics thereof.
With the highly informatization and intelligentization development of modern military technology, military weaponry is more dependent on the acquisition and control of all-round information of a battlefield. In the ultra-low altitude countermeasure, the target characteristic of the ultra-low altitude penetration weapon is the most important and basic resource of the air defense weapon system. Due to the existence of the underlying environment, the influence of the environment on the scattering characteristics of the ultra-low altitude target is not negligible. The composite electromagnetic scattering property of the ultra-low-altitude target and the environment is determined by various parameters, wherein the two most important items are the geometric construction and the electromagnetic parameters of the ultra-low-altitude target. The ultra-low altitude penetration weapon with excellent stealth performance firstly has unique appearance design to change the spatial distribution of radar wave scattering of a target so as to reduce the reflection and scattering of the ultra-low altitude target to radar waves in certain specific angle ranges; and then, by means of wave absorbing and wave transmitting materials, the electromagnetic parameters of the ultra-low altitude target are changed essentially, the radar scattering cross section is further reduced, and the observability of the target is reduced, so that the aim of stealth is fulfilled.
In past research, various scholars have conducted extensive research on the problem of electromagnetic scattering of targets, wherein the problem of scattering related to complex targets containing both conductors and media is not rare, and the scholars establish relevant electromagnetic scattering models and put forward or improve accurate and efficient electromagnetic calculation methods. However, when the environment is considered, the induced current generated on the surface of the conductor region of the composite target and the induced current, electricity and magnetic field generated on the surface of the medium region will act on the composite target as a secondary incident source of the environment, the induced current, electricity and magnetic field generated on the surface of the environment will act on each region of the composite target again as a new secondary incident source, the action is circulated until the electricity and magnetic current on the surface of all the regions tend to be stable, and the complexity of the scattering mechanism is far greater than that of the case of only considering the pure target. At present, research on the composite scattering characteristics of a conductor-medium composite target and an environment is rarely reported in domestic and foreign documents, few effective calculation methods are available, and the calculation efficiency is very low; in addition, the composite scattering characteristics of the composite target and the environment are obtained through an external field test, the cost is high, and the composite scattering characteristics are greatly limited, so that the research on the composite scattering characteristics of the conductor-medium composite target and the environment is innovative and challenging, and has strong military and civil application backgrounds and practical values.
Disclosure of Invention
The invention aims to provide a method for quickly calculating electromagnetic scattering of a conductor-medium composite target and an environment, which can quickly and accurately calculate the electromagnetic scattering property of the ultra-low-altitude conductor-medium composite target when the ultra-low-altitude conductor-medium composite target is positioned above the environment, thereby providing new theoretical and technical support for detecting the ultra-low-altitude conductor-medium composite target by a radar, and having important theoretical and application values in the fields of radar target detection and identification, ultra-low-altitude countermeasure, remote sensing and the like.
In order to achieve the purpose, the invention provides the following scheme: the invention provides a method for quickly calculating electromagnetic scattering of a conductor-medium composite target and environment, which comprises the following steps of:
s1, subdividing a target and an environment by adopting a triangular surface, performing geometric modeling on the ultra-low-altitude target and the environment to obtain an ultra-low-altitude conductor-medium composite target and environment composite scattering model, attributing the conductor-medium composite target to a JMCFIE-MSIE area, attributing the environment to a PO area, acquiring surface electricity and magnetic current of the PO area, and obtaining an induced electric field and an induced magnetic field of the surface of the environment;
s2, taking the induced electric field and the induced magnetic field on the surface of the environment and radar waves as conductor-medium composite target irradiation sources to obtain surface electricity and magnetic current of a conductor-medium composite target area and obtain surface induced electricity and magnetic field of the conductor-medium composite target area;
s3, taking induced electricity, magnetic field and radar wave of the surface of the conductor-medium composite target area as an environment irradiation source to obtain environment surface electricity and magnetic current, and obtain an induced electric field and an induced magnetic field of the environment surface;
and S4, repeating iteration of the step S2 and the step S3, stopping iteration when the iteration error is smaller than a specific value, and obtaining a final scattering electric field and a final scattering magnetic field through induced current and induced magnetic current of the composite target and the environment surface.
Preferably, the construction process of the composite electromagnetic scattering model of the ultra-low-altitude conductor-medium composite target and the environment is as follows: adopting a triangular surface element to subdivide a target and an environment, and constructing a geometric model of the target; generating an environment rough surface by adopting spectral functions with different statistical characteristics and a Monte Carlo method to obtain an environment geometric model; and placing the target geometric model at a specific distance above the environmental geometric model to obtain the ultra-low-altitude conductor-medium composite target and environmental composite electromagnetic scattering model.
Preferably, the subdivision density of the target and the environment by adopting the triangular surface element is 0.1 lambda-0.2 lambda, wherein lambda represents the radar working wavelength.
Preferably, the ultra-low-altitude-conductor medium composite target comprises a simple geometric body and a radar detection target;
the environment is a pure ideal conductor environment, and can also be a pure dielectric environment.
Preferably, the process of acquiring the induced electric field and the induced magnetic field of the environmental surface in step S1 is as follows: and supposing that the PO area surface current and the magnetic current are obtained by solving, and generating an induced electric field and an induced magnetic field according to the incident radar wave and the PO area surface current and the magnetic current.
Preferably, the iterative error of the current and the magnetic current of the composite target and the environmental surface is less than a specific value of 10 < -2 >.
The invention discloses the following technical effects:
the invention relates to the theoretical knowledge of electromagnetic scattering calculation, in the scattering calculation of an ultra-low-altitude conductor-medium composite target, a mixed integral equation with better iterative convergence is adopted, and the number of unknowns in a plurality of medium areas of the composite target is halved by combining a single integral equation method under the condition of a plurality of medium areas, so that the electromagnetic scattering of the composite target can be calculated more efficiently; in the calculation of the composite electromagnetic scattering, the induced current and the induced magnetic current on the surface of the environment are obtained by using a physical optical method, and an iterative algorithm structure of a numerical value-high frequency approximation algorithm is constructed, so that a mixed JMCFIE-MSIE-PO iterative algorithm is constructed.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a schematic flow chart of a method for rapidly calculating electromagnetic scattering between a conductor-medium composite target and an environment according to the present invention;
FIG. 2 is a schematic diagram of a conductor-dielectric composite target and an environmental geometric model according to the present invention;
FIG. 3 is a schematic diagram of a two-station RCS curve of a conductor-dielectric composite almond body and environment according to an embodiment of the present invention;
fig. 4 is a schematic diagram of a single station RCS curve of a conductor-dielectric composite helicopter and the environment according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
As shown in fig. 1, the present invention provides a method for fast calculating electromagnetic scattering between a conductor-medium composite target and an environment, comprising the following steps:
and S1, subdividing the target and the environment by adopting a triangular surface, and constructing a composite electromagnetic scattering model of the ultra-low-altitude conductor-medium composite target and the environment.
The ultra-low-altitude conductor-medium composite target specifically covers simple geometric bodies such as a sphere, a cube, an almond body and the like, and also comprises a general name of radar detection targets such as an ultra-low-altitude cruise missile, a fighter plane, a helicopter and the like; the environment refers to a rough surface generated by adopting spectral functions with different statistical characteristics through a Monte Carlo method, and refers to a specific rough surface geometric model, and the environment can be a pure ideal conductor environment or a pure medium environment.
The target subdivision adopts the existing CAD, 3Dmax and Rhino target geometric models; in the aspect of environment, spectrum functions with different statistical characteristics (Gaussian spectrum functions are mostly adopted on land, PM spectrum functions are mostly adopted on sea surface) are adopted, a rough surface is generated by a Monte Carlo method, finally, a geometric model of a target is placed above the geometric model of the environment for a certain distance, the target and the environment are divided by a triangular surface element, the division density is generally 0.1 lambda-0.2 lambda, wherein lambda represents the radar working wavelength. And finishing the construction of the composite electromagnetic scattering model of the target and the environment after subdivision is finished.
The geometric model of the ultra-low-altitude conductor-dielectric composite target and environment constructed in the embodiment is shown in fig. 2, wherein Ω0Representing free space; omega1Representing a conductor region; s1A surface representing a conductor region; omega2，Ω3，…ΩnRespectively representing 2,3, … n-th media areas, S2，S3，…SnSurface representing 2,3, … n media regions; sΓlpRepresenting the interface of the conductor region with the p-th dielectric region; ssurRepresenting an environmental surface; the inner surface of the p-th dielectric region is SdpIt is represented by SpAnd SΓ1pAre surrounded together to form; the value of the subscript p described above is taken to be 2,3, … n, where n represents the total number of media regions.
The ultra-low-altitude conductor-dielectric composite target has a certain distance from the surface of the environment, and the specific distance can be set according to the requirement.
S2, dividing the ultra-low-altitude conductor-medium composite target and the conductor-medium composite target in the environment composite scattering model into a JMCFIE-MSIE (the electric/magnetic current combined field integral equation, JMCFIE, mixed field integral equation, the single integral equation method for multiple dielectric objects in the composite target, MSIE, single integral equation method for the conductor and multiple dielectric composite target) area, and obtaining the induced electricity and the magnetic current of the conductor-medium composite target according to the induced electric field and the induced magnetic field of the environment surface, and obtaining the induced electricity and the magnetic field of the conductor-medium composite target.
According to FIG. 2, J1Denotes the conductor surface S1Induced current on; j. the design is a squarepAnd MpRespectively represent the p-th medium surface SpInduced current onThe JMCFIE equation for inducing magnetic flux, thereby establishing a conductor-dielectric composite target, is as follows:
wherein: e.g. of the typeinc＝n×EincX n and jinc＝n×HincRespectively representing the tangential electric field and the tangential magnetic field of the incident wave; α and β represent blending factors that satisfy the relationship: α is 0. ltoreq. α, β. ltoreq.1 and α + β ═ 1;
is the wave impedance in free space; l is0z(. a) and
respectively representing an electric field integral operator and a magnetic field integral operator in free space, and the expression is as follows:
wherein: g0(r, r') represents a Green function in free space, expressed as
Representing the wavenumber in free space; p.v. denotes the principal value integral.
Both the current and the magnetic current can be developed in the following way:
wherein: f. of1jAnd fpjRespectively represent the conductor surface S1With the p-th dielectric surface SpRWG basis functions on (p ═ 2.., n). The total number of edges on the surface of the conductor is N1The common edge of the conductor and the dielectric block is included; the total number of edges on the surface of the p-th medium is NpWherein no common edge is included.
Defining the expression of a complex operator as<u,v>Γ＝∫Γu · vd Γ. Selecting f by Galerkin method1jAnd fpjTo test the function, a matrix equation can be established by equations (1) and (2) as:
wherein:
the dimension of the matrix equation in equation (7) is (N)1+N2+...+Nn)×(N1+2N2+...+2Nn)，The relationship between the unknown number and the equation number cannot satisfy the final solution of the current and the magnetic current. Therefore, an additional set of equations must be established for matrix equation (7). Therefore, the embodiment adopts the single effective flow in the single integral equation method to represent the internal field of the pth block medium area
And
wherein
Showing the outer surface S of the dielectric blockpAnd the boundary surface SΓ1pAn enclosed dielectric region inner surface; they satisfy the following relationship:
and
here, the
Representing wave impedance in the pth dielectric region; in addition to this, an electric field integration operator in the p-th dielectric region
Magnetic field integral operator
And green function GpThe wave number in (r, r') is
It is to be noted that the internal fields are represented by significant flows, which must satisfy the boundary conditions, at the outer surface S of the p-th dielectric blockpAbove, there is the equation:
at the interface S between the conductor region and the p-th dielectric regionΓ1pThe boundary conditions are also satisfied:
interface S of conductor region and p dielectric regionΓ1pThe effective flow of upper stream disappears.
Thus, the effective flow of the inner surface of the pth dielectric block through the RWG basis function is developed in the form of:
in the formula (14), the total unknown edge number is Np+NΓ1pIn which N isΓ1pIs the interface S of the conductor region and the p-th dielectric regionΓ1pThe number of unknowns in (c). The expansion factor of equation (14) can be determined from the mean value x of the effective flow across the edgepjTo approximate the representation:
the simultaneous expression (11), the expression (12), and the expression (15), and the expansion coefficients { b } and { c } in the expression (7) satisfy the following matrix equation:
wherein: { xpDenotes a single significant flow expansion coefficient vector, and the expression of each element in equation (16) is as followsThe following:
here, when i is j, δ ij1 is ═ 1; when i ≠ j, δ ij0. Formula (16) is substituted for formula (7), and the new resultant matrix is as follows:
the expression of each element in the formula is as follows:
Ap1＝Qp1 (21)
in the formulas (21) to (26), the subscripts p and q take the values of: p is 1,2, n, q is 2, 3. By using the SIE method, the unknown quantity of the matrix equation (20) at this time is represented by N of the JMCFIE method1+N2+…+NnBecome into
And S3, obtaining a final scattering electric field and a final scattering magnetic field according to the induced electric field and the induced magnetic field of the ultra-low-altitude conductor-medium composite target and the induced electric field and the induced magnetic field of the environment surface.
It should be noted that, in order to take into account the edge diffraction effect of the environment when calculating the large scattering angle range RCS, the present embodiment is modified by the EEC method.
In this embodiment, a mixed JMCFIE-MSIE-PO iterative algorithm is obtained according to the JMCFIE-MSIE mixing algorithm of the ultra-low-altitude conductor-medium target, and the calculation of the composite scattering of the ultra-low-altitude conductor-medium target and the environment is completed.
In the iterative method, it is assumed that the induced electricity and magnetic flow of a PO (physical optics, PO) region are solved, and then the induced electricity and magnetic field excited by these surface flows can be used as a secondary radiation source and used as an irradiation source of a target together with an incident wave.
This example gives the surface S for the kth iteration1，S2，…，SnJMCFIE equation of (a):
wherein: when k is equal to 1, the first step is carried out,
and
respectively representing the tangential electric field and the tangential magnetic field of the environment surface of the k-1 st order;
and
the expression of (a) is:
on the surface of the medium environment, the surface electric field and the magnetic field are tangentially continuous, and the corresponding equation is as follows:
wherein: eincAnd HincRespectively representing an incident electric field and an incident magnetic field;
and
respectively representing incident electric field and incident magnetic field incident to the environment surface when the conductor-medium composite target is used as a secondary radiation source；
And
respectively representing the transmitted electric field and the transmitted magnetic field of the ambient surface.
Equations (31) and (32) can be further written as:
the overall iterative process starts with induced currents and induced magnetic currents of order 0 in the PO region. Based on the iterative model, the iterative process is circularly carried out until the induced currents and the induced magnetic currents on the surfaces of the PO area and the JMCFIE-MSIE area are stable. In the whole iteration process, the iteration error expression of the induction current and the induction magnetic current of the k-th order is as follows:
wherein: i | · | purple wind2Representing a 2 norm.
The iteration process of the JMCFIE-MSIE-PO mixing algorithm is to simultaneously reduce the iteration error of the induced current and the induced magnetic current to be less than 10-2And then stop. A multi-layer Fast Multipole Algorithm (MLFMA) can be applied to accelerate matrix multiplication in the computation process.
In order to verify the effectiveness of the present invention, the present embodiment performs simulation comparative analysis on the present invention and the MLFMA algorithm based on the traditional hybrid integral equation, specifically:
the simulation computing platform adopted by the embodiment is a 64-core 2.30GHz main frequency AMD processor, and the working frequency is set to be 1.0GHz, the environmental surface is generated by a Gaussian spectrum function with the size set to 40 lambda x 40 lambda and the root mean square height hrms0.1 lambda, and a correlation length l in the x and y axesx＝ly2.0 lambda, relative dielectric constant of environment is epsilonr(10.93, -j 6.48); the conductor-dielectric composite target is set as a composite almond body model with the three-dimensional size of 3.79m multiplied by 1.46m multiplied by 0.49m, the composite almond body model consists of an ideal conductor at the middle part and two different dielectric blocks at the two sides, and the height from the environment is set as 5 lambda. The incident angles are set as: pitch angle thetaiAt 45 deg. azimuth phi i0 °; the observation angle is set as: thetas＝-90°～90°，φs＝0°。
The two-station RCS curve of the conductor-medium composite almond body and the environment calculated by the JMCFIE-MSIE-PO hybrid iterative algorithm and the MLFMA algorithm obtained in this example is shown in fig. 3.
As can be seen from fig. 3, the number of environment unknowns is 388435, when the MLFMA algorithm is adopted, the number of conductor-medium composite almond bodies is 74045, and when the JMCFIE-MSIE-PO hybrid iterative algorithm is adopted, the number of unknowns of the composite target is reduced to 46943. Observing fig. 3, it can be obtained that the two curves are well matched at each scattering angle; both curves show a distinct peak around a scattering angle of-45. When the MLFMA algorithm is adopted for calculation, the memory requirement and the calculation time are 1.88GB and 1.95h respectively; when the JMCFIE-MSIE-PO hybrid iterative algorithm is adopted, the memory requirement and the calculation time are respectively 0.84GB and 0.52 h. Therefore, compared with the MLFMA algorithm, the JMCFIE-MSIE-PO hybrid iterative algorithm has the advantages that the memory requirement and the calculation speed are greatly improved.
In addition, in the embodiment, a JMCFIE-MSIE-PO mixed iterative algorithm is adopted for actual calculation, the working frequency is set to be 0.6GHz, the environment surface is generated by a Gaussian spectrum function, the size of the environment surface is set to be 100 lambda multiplied by 100 lambda, and the root-mean-square height is hrms0.1 lambda, and a correlation length l in the x and y axesx＝ly2.0 lambda, relative dielectric constant of environment is epsilonr(10.94, -j 7.18); the conductor-dielectric composite target is set to three-dimensional dimensions of 5.57m × 3.22m × 4.26The m composite helicopter model is composed of a body main body shown in figure 4 as an ideal conductor, a cartridge frame, wings and a tail wing which are all three media with different relative dielectric constants, and the height from the environment is set to be 5.0 m. The incident azimuth angle is fixed to phiiWhen the angle is equal to 0 degrees, the incident pitch angle is set as: thetai＝-90°～90°。
FIG. 4 is a single station RCS curve of the conductor-medium composite helicopter and the environment calculated by respectively adopting a JMCFIE-MSIE-PO hybrid iterative algorithm and an MLFMA algorithm. As can be seen from fig. 4, the composite scattering RCS curves of the three targets and the environment all show a distinct peak around the scattering angle of 0 °. If the parts 1,2 and 3 are respectively composed of dielectric constants epsilon as shown in the legend of FIG. 4, compared with the helicopter target made of pure ideal conductor material, the helicopter target is located above the environmentr1＝(20.0,-j10.0)、εr2═ 15.0, -j7.5 and ∈r3Medium (10.0, -j5.0) with a significant decrease in RCS values at other scattering angles than around 0 °; further the helicopter is integrally composed of a dielectric constant of epsilonr3The RCS value at other scattering angles than around a scattering angle of 0 ° is further reduced by replacing the medium (10.0, -j 5.0).
In addition, in this embodiment, the pure target situation is compared to show the target scattering property of the conductor-dielectric composite helicopter in detail, and compared with the helicopter made of pure ideal conductor material, when the dielectric constants of the parts 1,2 and 3 are respectively epsilonr1＝(20.0,-j10.0)、εr2═ 15.0, -j7.5 and ∈r3When the medium (10.0, -j5.0) is replaced, the target scattering properties will be significantly changed, especially at scattering angles of-75 °, -62 °, -45 °, -41 °, 34 °, 0 °, 6 °, 56 ° and 77 °; further the helicopter is integrally composed of a dielectric constant of epsilonr3The RCS values for the target at most scattering angles are further reduced by replacing the medium (10.0, -j 5.0).
Therefore, the conclusion that the multiple pieces of medium contained in the conductor-medium composite target cause the change of the scattering characteristics of the corresponding scattering angle range is also consistent with the principle of enhancing the stealth performance of the ultra-low altitude penetration prevention target.
The above-described embodiments are merely illustrative of the preferred embodiments of the present invention, and do not limit the scope of the present invention, and various modifications and improvements of the technical solutions of the present invention can be made by those skilled in the art without departing from the spirit of the present invention, and the technical solutions of the present invention are within the scope of the present invention defined by the claims.
Claims (6)
1. A conductor-medium composite target and environment electromagnetic scattering rapid calculation method is characterized by comprising the following steps:
s1, subdividing a target and an environment by adopting a triangular surface, performing geometric modeling on the ultra-low-altitude target and the environment to obtain an ultra-low-altitude conductor-medium composite target and environment composite scattering model, attributing the conductor-medium composite target to a JMCFIE-MSIE area, attributing the environment to a PO area, acquiring surface electricity and magnetic current of the PO area, and obtaining an induced electric field and an induced magnetic field of the surface of the environment;
s2, taking the induced electric field and the induced magnetic field on the surface of the environment and radar waves as conductor-medium composite target irradiation sources to obtain surface electricity and magnetic current of a conductor-medium composite target area and obtain surface induced electricity and magnetic field of the conductor-medium composite target area;
s3, taking induced electricity, magnetic field and radar wave of the surface of the conductor-medium composite target area as an environment irradiation source to obtain environment surface electricity and magnetic current, and obtain an induced electric field and an induced magnetic field of the environment surface;
and S4, repeating iteration of the step S2 and the step S3, stopping iteration when the iteration error is smaller than a specific value, and obtaining a final scattering electric field and a final scattering magnetic field through induced current and induced magnetic current of the composite target and the environment surface.
2. The conductor-medium composite target and environment electromagnetic scattering rapid calculation method according to claim 1, wherein the construction process of the composite electromagnetic scattering model of the ultra-low-altitude conductor-medium composite target and environment is as follows: adopting a triangular surface element to subdivide a target and an environment, and constructing a geometric model of the target; generating an environment rough surface by adopting spectral functions with different statistical characteristics and a Monte Carlo method to obtain an environment geometric model; and placing the target geometric model at a specific distance above the environmental geometric model to obtain the ultra-low-altitude conductor-medium composite target and environmental composite electromagnetic scattering model.
3. The method for rapidly calculating the electromagnetic scattering of the conductor-medium composite target and the environment according to claim 2, wherein the subdivision density of the target and the environment by adopting the triangular surface element is 0.1 λ -0.2 λ, wherein λ represents the radar working wavelength.
4. The conductor-medium composite target and environment electromagnetic scattering rapid calculation method according to claim 1, wherein the ultra-low-altitude-conductor-medium composite target comprises a simple geometric body and a radar detection target;
the environment is a pure ideal conductor environment, and can also be a pure dielectric environment.
5. The method for rapidly calculating the electromagnetic scattering between the conductor-medium composite target and the environment according to claim 1, wherein the step S1 is implemented by acquiring the induced electric field and the induced magnetic field on the surface of the environment through a process of: and supposing that the PO area surface current and the magnetic current are obtained by solving, and generating an induced electric field and an induced magnetic field according to the incident radar wave and the PO area surface current and the magnetic current.
6. The method for rapidly calculating conductor-medium composite target and environment electromagnetic scattering according to claim 1, wherein the iterative error of the current and the magnetic current of the composite target and the environment surface is less than a specific value of 10-2。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011623614.7A CN112733347B (en) | 2020-12-31 | 2020-12-31 | Conductor-medium composite target and environment electromagnetic scattering rapid calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011623614.7A CN112733347B (en) | 2020-12-31 | 2020-12-31 | Conductor-medium composite target and environment electromagnetic scattering rapid calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112733347A true CN112733347A (en) | 2021-04-30 |
| CN112733347B CN112733347B (en) | 2022-09-02 |
Family
ID=75609358
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011623614.7A Active CN112733347B (en) | 2020-12-31 | 2020-12-31 | Conductor-medium composite target and environment electromagnetic scattering rapid calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112733347B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113610988A (en) * | 2021-07-30 | 2021-11-05 | 上海无线电设备研究所 | Ray tracing-based electromagnetic scattering simulation method for complex structure target |
| CN114741880A (en) * | 2022-04-15 | 2022-07-12 | 北京理工大学 | Method and device for calculating electromagnetic characteristics of group targets |
| CN114842052A (en) * | 2022-05-06 | 2022-08-02 | 中国人民解放军空军工程大学 | Moving target identification method based on target and environment composite electromagnetic scattering field |
| CN117147977A (en) * | 2023-07-31 | 2023-12-01 | 西安电子科技大学 | Test quality estimation method for unmanned aerial vehicle external field composite scattering measurement |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030158715A1 (en) * | 2000-06-21 | 2003-08-21 | Ottusch John Jacob | Method and apparatus for modeling three-dimensional electromagnetic scattering from arbitrarily shaped three-dimensional objects |
| CN104778286A (en) * | 2014-01-09 | 2015-07-15 | 南京理工大学 | High-speed simulation method for electromagnetic scattering characteristics of sea skimmer |
| CN106295120A (en) * | 2016-07-26 | 2017-01-04 | 上海无线电设备研究所 | A kind of electromagnetic characteristic of scattering quick calculation method |
| CN106646402A (en) * | 2016-11-11 | 2017-05-10 | 上海无线电设备研究所 | Electromagnetic pulse radar echo simulation method in offshore area |
| CN110907901A (en) * | 2019-10-29 | 2020-03-24 | 西安电子科技大学 | Sea warfare passive interference calculation and evaluation method based on TDSBR |
| CN111948618A (en) * | 2020-06-30 | 2020-11-17 | 西安电子科技大学 | Forward scattering target detection method and system based on satellite external radiation source |
| CN112014815A (en) * | 2020-08-12 | 2020-12-01 | 网络通信与安全紫金山实验室 | Hybrid surface integral analysis method, device, equipment and storage medium |
-
2020
- 2020-12-31 CN CN202011623614.7A patent/CN112733347B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030158715A1 (en) * | 2000-06-21 | 2003-08-21 | Ottusch John Jacob | Method and apparatus for modeling three-dimensional electromagnetic scattering from arbitrarily shaped three-dimensional objects |
| CN104778286A (en) * | 2014-01-09 | 2015-07-15 | 南京理工大学 | High-speed simulation method for electromagnetic scattering characteristics of sea skimmer |
| CN106295120A (en) * | 2016-07-26 | 2017-01-04 | 上海无线电设备研究所 | A kind of electromagnetic characteristic of scattering quick calculation method |
| CN106646402A (en) * | 2016-11-11 | 2017-05-10 | 上海无线电设备研究所 | Electromagnetic pulse radar echo simulation method in offshore area |
| CN110907901A (en) * | 2019-10-29 | 2020-03-24 | 西安电子科技大学 | Sea warfare passive interference calculation and evaluation method based on TDSBR |
| CN111948618A (en) * | 2020-06-30 | 2020-11-17 | 西安电子科技大学 | Forward scattering target detection method and system based on satellite external radiation source |
| CN112014815A (en) * | 2020-08-12 | 2020-12-01 | 网络通信与安全紫金山实验室 | Hybrid surface integral analysis method, device, equipment and storage medium |
Non-Patent Citations (2)
| Title |
|---|
| 刘万萌 等: "不同类型土壤地面与其上方目标散射特性研究", 《微波学报》 * |
| 李宁 等: "船海复合场景面元化快速电磁建模方法及SAR成像仿真", 《国防科技大学学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113610988A (en) * | 2021-07-30 | 2021-11-05 | 上海无线电设备研究所 | Ray tracing-based electromagnetic scattering simulation method for complex structure target |
| CN114741880A (en) * | 2022-04-15 | 2022-07-12 | 北京理工大学 | Method and device for calculating electromagnetic characteristics of group targets |
| CN114842052A (en) * | 2022-05-06 | 2022-08-02 | 中国人民解放军空军工程大学 | Moving target identification method based on target and environment composite electromagnetic scattering field |
| CN114842052B (en) * | 2022-05-06 | 2024-04-05 | 中国人民解放军空军工程大学 | Moving target identification method based on target and environment composite electromagnetic scattering field |
| CN117147977A (en) * | 2023-07-31 | 2023-12-01 | 西安电子科技大学 | Test quality estimation method for unmanned aerial vehicle external field composite scattering measurement |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112733347B (en) | 2022-09-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112733347B (en) | Conductor-medium composite target and environment electromagnetic scattering rapid calculation method | |
| Jenn | Radar and laser cross section engineering | |
| Cools et al. | Accurate and conforming mixed discretization of the MFIE | |
| CN103870654A (en) | Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing | |
| Chatzigeorgiadis | Development of code for a physical optics radar cross section prediction and analysis application | |
| Sukharevsky et al. | Modeling of ultrawideband (UWB) impulse scattering by aerial and subsurface resonant objects based on integral equation solving | |
| Ananth et al. | Synthesis of radar absorption material for stealth application | |
| CN111460579B (en) | Three-dimensional target stealth optimization method | |
| Shin et al. | Analysis of radome cross section of an aircraft equipped with a FSS radome | |
| Chung et al. | Efficacy of an S-shaped air inlet on the reduction of front bistatic radar cross section of a fighter engine | |
| CN104732050A (en) | Method for estimating distribution of electromagnetism in flying targets made of carbon fiber materials under lightning pulses | |
| Chen et al. | Numerical simulation on the radar cross section of variable-sweep wing aircraft | |
| Chung et al. | Radar cross section analysis of stealth fighter design: Key factors and limitations of simulation | |
| CN107301301B (en) | Construction method of magnetic loss type shrinkage ratio wave-absorbing material | |
| Charris et al. | Analysis of radar cross section assessment methods and parameters affecting it for surface ships | |
| CN106156431B (en) | Target conductor electromagnetic scattering emulation mode based on nurbs surface modeling | |
| Zalevsky et al. | Calculation of scattering characteristics of aerial radar objects of resonant sizes based on iterative algorithm | |
| Sun | Influence evaluation of UAV inlet on electromagnetic scattering and time-frequency characteristics | |
| Chen et al. | A hybrid MoM-PO method combining ACA technique for electromagnetic scattering from target above a rough surface | |
| CN115169170A (en) | Composite target scattering semi-analytic rapid calculation method based on non-uniform grid model | |
| CN105303022A (en) | Gaussian beam method for quickly obtaining electromagnetic scattering property of target | |
| Huang et al. | Recognition of the major scattering sources on complex targets based on the high frequency radar cross section integrated calculation technique | |
| Heidar et al. | Simulation and reduction of radar cross section the unmanned aerial vehicle in X-band using Shaping technique | |
| Pelosi et al. | A preliminary assessment of specular radar cross section performance in the chengdu J-20 prototype | |
| Al Salameh et al. | COMSOL solutions for reducing radar cross section of drones |
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 |