CN117057098B - Frequency-selective radome analysis method based on physical optical method and equivalent medium model - Google Patents
Frequency-selective radome analysis method based on physical optical method and equivalent medium model Download PDFInfo
- Publication number
- CN117057098B CN117057098B CN202310862769.3A CN202310862769A CN117057098B CN 117057098 B CN117057098 B CN 117057098B CN 202310862769 A CN202310862769 A CN 202310862769A CN 117057098 B CN117057098 B CN 117057098B
- Authority
- CN
- China
- Prior art keywords
- electric field
- magnetic field
- antenna aperture
- magnetic
- antenna
- 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 46
- 238000004458 analytical method Methods 0.000 title claims abstract description 23
- 230000003287 optical effect Effects 0.000 title claims abstract description 20
- 230000005684 electric field Effects 0.000 claims abstract description 92
- 230000005540 biological transmission Effects 0.000 claims description 42
- 230000010287 polarization Effects 0.000 claims description 25
- 230000035699 permeability Effects 0.000 claims description 19
- 238000010586 diagram Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 3
- 239000002184 metal Substances 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aerials With Secondary Devices (AREA)
Abstract
The invention provides a frequency-selective radome analysis method based on a physical optical method and an equivalent medium model, wherein the frequency-selective radome comprises a radome and a frequency-selective surface, and the method comprises the following steps: and the frequency-selective radome is equivalent to a medium model, the medium model is analyzed to obtain an outer surface electric field and an outer surface magnetic field of the medium model, and then the outer surface electric field and the outer surface magnetic field are subjected to integral operation to obtain an electric field and a magnetic field at any point of a space, so that the electric performance analysis of the frequency-selective radome is realized, and the modeling and calculating time of the frequency-selective radome is greatly shortened while the electric performance analysis precision is ensured.
Description
Technical Field
The invention relates to the technical field of radomes, in particular to a frequency-selective radome analysis method based on a physical optical method and an equivalent medium model.
Background
The frequency selective surface is usually an infinite array formed by periodically arranging identical metal patch units or open pore units on a metal screen along one-dimensional or two-dimensional directions, and has selectivity to the frequency, polarization mode, incidence angle and the like of the incident electromagnetic wave, which is equivalent to an electromagnetic wave filter in an open space.
The radome is an electromagnetic window positioned at the front end of the aircraft, has good aerodynamic characteristics, bears the erosion of aerodynamic load, aerodynamic heat, rain, snow and the like in the flying process, ensures the good electromagnetic working environment of an internal antenna system, and furthest reduces the influence on the communication electromagnetic wave of the radar system. The radome loaded with the frequency selective surface is called a frequency selective surface radome, and is called a frequency selective radome for short.
The radome must maintain the streamline shape and satisfy the electrical performance, but the streamline shape and the electromagnetic performance are a pair of contradictors. Therefore, after designing the radome, it is necessary to analyze the electrical performance thereof. The analysis method for the electrical performance of the common radome includes a physical optical method, a ray tracing method, a full wave method and the like. However, the current analysis method for the electrical performance of the frequency-selective radome only comprises a full-wave method, and although the full-wave method can realize accurate electrical performance analysis, the time is long.
In view of the above, the invention provides a frequency-selective radome analysis method based on a physical optical method and an equivalent medium model, which can greatly shorten the modeling and calculation time of the frequency-selective radome while ensuring the accuracy of the electrical performance analysis of the frequency-selective radome.
Disclosure of Invention
According to the technical problem that the existing frequency-selective radome analysis consumes longer time, the frequency-selective radome analysis method based on the physical optical method and the equivalent medium model is provided. According to the invention, the frequency-selective radome is equivalent to a medium model, and the medium model is used for analyzing the electrical properties of the frequency-selective radome, so that the method is simple, and the modeling and calculating time of the frequency-selective radome is greatly shortened.
The invention adopts the following technical means:
the invention provides a frequency-selective radome analysis method based on a physical optical method and an equivalent medium model, wherein the frequency-selective radome comprises a radome and a frequency-selective surface, and the method comprises the following steps:
The frequency-selective radome is equivalent to a medium model, the medium model has the same size and shape as the frequency-selective radome, the medium model comprises a first curved surface plate, a second curved surface plate and a third curved surface plate which are laminated, and the thickness of the second curved surface plate is the same as that of the frequency-selective surface;
Acquiring an antenna aperture surface electric field and an antenna aperture surface magnetic field;
Obtaining antenna aperture surface magnetic current according to the antenna aperture surface electric field, and obtaining antenna aperture surface current according to the antenna aperture surface magnetic field;
Obtaining an inner surface electric field and an inner surface magnetic field of the medium model according to the antenna aperture surface magnetic current and the antenna aperture surface current;
Obtaining an outer surface electric field of the medium model according to the inner surface electric field, and obtaining an outer surface magnetic field of the medium model according to the inner surface magnetic field;
And carrying out integral operation on the outer surface electric field and the outer surface magnetic field to obtain an electric field and a magnetic field at any point in space.
Further, the antenna aperture surface magnetic current is obtained according to the antenna aperture surface electric field, and is calculated according to the following mode:
M＝E×n；
wherein M is the magnetic current of the aperture surface of the antenna, E is the electric field of the aperture surface of the antenna, and n is the normal vector of the aperture surface of the antenna.
Further, the antenna aperture plane current is obtained according to the antenna aperture plane magnetic field, and is calculated according to the following mode:
J＝n×H；
Wherein J is the current of the aperture plane of the antenna, n is the normal vector of the aperture plane of the antenna, and H is the magnetic field of the aperture plane of the antenna.
Further, the obtaining the internal surface electric field and the internal surface magnetic field of the dielectric model according to the antenna aperture surface magnetic current and the antenna aperture surface current sum comprises the following steps:
obtaining a magnetic vector according to the current of the antenna aperture surface, and obtaining an electric vector according to the magnetic current of the antenna aperture surface;
and obtaining the inner surface electric field and the inner surface magnetic field according to the magnetic sagittal position and the electric sagittal position.
Further, the magnetic vector position is obtained according to the antenna aperture plane current, and is calculated according to the following mode:
Wherein A is the magnetic vector, mu is magnetic permeability, S 1 is an antenna aperture plane, J is the antenna aperture plane current, J is an imaginary number, k is a free space wave vector, and R is the distance from a point P to the antenna aperture plane.
Further, the electric vector position is obtained according to the magnetic current of the antenna aperture surface, and is calculated according to the following mode:
Wherein F is the electrical vector, epsilon is the dielectric constant, S 1 is the antenna aperture plane, M is the magnetic flow of the antenna aperture plane, j is the imaginary number, k is the free space wave vector, and R is the distance from the point P to the antenna aperture plane.
Further, the inner surface electric field and the inner surface magnetic field are obtained according to the magnetic sagittal position and the electric sagittal position, and are calculated according to the following modes:
Wherein E (P) is the inner surface electric field, A is the magnetic sagittal plane, F is the electric sagittal plane, ε is the dielectric constant, j is the imaginary number, w is the angular frequency, μ is the magnetic permeability, To represent the sign of rotation, H (P) is the internal surface magnetic field.
Further, obtaining an outer surface electric field of the dielectric model according to the inner surface electric field comprises the following steps:
the inner surface electric field is decomposed into an inner surface electric field horizontal component and an inner surface electric field vertical component, calculated as follows:
E’＝E⊥×T⊥+E∥×T∥；
Wherein E' is the outer surface electric field, E ⊥ is the inner surface electric field vertical component, T ⊥ is the vertical polarization transmission coefficient, E ∥ is the inner surface electric field horizontal component, and T ∥ is the horizontal polarization transmission coefficient.
Further, the obtaining the external surface magnetic field of the medium model according to the internal surface magnetic field includes:
decomposing the internal surface magnetic field into an internal surface magnetic field horizontal component and an internal surface magnetic field vertical component, calculated as follows:
H’＝H⊥×T⊥+H∥×T∥；
Wherein H' is the external surface magnetic field, H ⊥ is the internal surface magnetic field vertical component, T ⊥ is the vertical polarization transmission coefficient, H ∥ is the internal surface magnetic field horizontal component, and T ∥ is the horizontal polarization transmission coefficient.
Further, the vertical polarization transmission coefficient and the horizontal polarization transmission coefficient are obtained as follows:
obtaining a normal incidence transmission coefficient and a normal incidence reflection coefficient of electromagnetic waves perpendicularly incident on the second curved plate;
obtaining the wave impedance and complex refractive index of the second curved plate according to the normal incidence transmission coefficient and the normal incidence reflection coefficient;
Obtaining the magnetic permeability and the dielectric constant of the second curved plate according to the wave impedance and the complex refractive index;
determining an incident angle of electromagnetic waves according to the inner surface electric field and the inner surface magnetic field;
And calculating the vertical polarization transmission coefficient and the horizontal polarization transmission coefficient by adopting a transmission line theory according to the incidence angle, the magnetic permeability and the dielectric constant.
Compared with the prior art, the invention has the following advantages:
The invention provides a frequency-selective radome analysis method based on a physical optical method and an equivalent medium model, wherein the frequency-selective radome comprises a radome and a frequency-selective surface, and the method comprises the following steps: the frequency-selective radome is equivalent to a medium model, the medium model has the same size and shape as the frequency-selective radome, the medium model comprises a first curved surface plate, a second curved surface plate and a third curved surface plate which are laminated, and the thickness of the second curved surface plate is the same as that of the frequency-selective surface; acquiring an antenna aperture surface electric field and an antenna aperture surface magnetic field; obtaining antenna aperture surface magnetic current according to the antenna aperture surface electric field and obtaining antenna aperture surface current according to the antenna aperture surface magnetic field; obtaining an inner surface electric field and an inner surface magnetic field of the medium model according to the magnetic current of the antenna caliber surface and the current sum of the antenna caliber surface; obtaining an outer surface electric field of the medium model according to the inner surface electric field, and obtaining an outer surface magnetic field of the medium model according to the inner surface magnetic field; and carrying out integral operation on the outer surface electric field and the outer surface magnetic field to obtain the electric field and the magnetic field at any point in space. The frequency-selective radome is equivalent to a medium model, the medium model is analyzed to obtain an outer surface electric field and an outer surface magnetic field of the medium model, and then integral operation is carried out on the outer surface electric field and the outer surface magnetic field to obtain an electric field and a magnetic field at any point of a space, so that the electric performance analysis of the frequency-selective radome is realized, the accuracy of the electric performance analysis is ensured, and meanwhile, the modeling and calculating time of the frequency-selective radome is greatly shortened.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings may be obtained according to the drawings without inventive effort to a person skilled in the art.
Fig. 1 is a flowchart of a frequency selective radome analysis method based on a physical optical method and an equivalent medium model.
Fig. 2 is a schematic structural diagram of a dielectric model according to the present invention.
Fig. 3 is a block diagram of a frequency selective surface provided by the present invention.
Fig. 4 is an enlarged view at a in fig. 3.
Fig. 5 is an antenna far field pattern.
In the figure: 0. a dielectric model 1, a first curved plate; 2. a second curved plate member; 3. and a third curved plate.
Detailed Description
In order that those skilled in the art will better understand the present invention, a technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in which it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present invention without making any inventive effort, shall fall within the scope of the present invention.
It should be noted that the terms "first," "second," and the like in the description and the claims of the present invention and the above figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that the embodiments of the invention described herein may be implemented in sequences other than those illustrated or otherwise described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
Referring to fig. 1, fig. 2, fig. 3, and fig. 4, fig. 1 is a flowchart of a method for analyzing a frequency selective radome based on a physical optical method and an equivalent medium model according to the present invention, fig. 2 is a schematic structural diagram of a medium model according to the present invention, fig. 3 is a structural diagram of a frequency selective surface according to the present invention, and fig. 4 is an enlarged view at a in fig. 3, to illustrate a specific embodiment of a method for analyzing a frequency selective radome based on a physical optical method and an equivalent medium model according to the present invention, where the frequency selective radome includes a radome and a frequency selective surface, and includes:
S1: the frequency-selective radome is equivalent to a medium model 0, the medium model 0 has the same size and shape as the frequency-selective radome, the medium model 0 comprises a first curved surface plate 1, a second curved surface plate 2 and a third curved surface plate 3 which are laminated, and the thickness of the second curved surface plate 2 is the same as that of the frequency-selective surface;
S2: acquiring an antenna aperture surface electric field and an antenna aperture surface magnetic field;
s3: obtaining antenna aperture surface magnetic current according to the antenna aperture surface electric field and obtaining antenna aperture surface current according to the antenna aperture surface magnetic field;
S4: obtaining an inner surface electric field and an inner surface magnetic field of the medium model 0 according to the magnetic current of the antenna caliber surface and the current sum of the antenna caliber surface;
s5: obtaining an outer surface electric field of the medium model 0 according to the inner surface electric field, and obtaining an outer surface magnetic field of the medium model 0 according to the inner surface magnetic field;
S6: and carrying out integral operation on the outer surface electric field and the outer surface magnetic field to obtain the electric field and the magnetic field at any point in space.
It can be understood that the first curved plate 1, the second curved plate 2 and the third curved plate 3 are curved plates with constant thickness and uniform material, the first curved plate 1 and the third curved plate 3 are made of the same material, the first curved plate 1 and the second curved plate 2 are made of different materials, the second curved plate 2 is an equivalent medium surface of a frequency selective surface, referring to fig. 3 and 4, fig. 4 only illustrates a grid square annular frequency selective surface topology structure, the grid square annular frequency selective surface topology structure shown in fig. 4 is arranged in an array to obtain a frequency selective surface as shown in fig. 3, and the side length p of the grid square annular frequency selective surface topology structure shown in fig. 4 is 9 mm, so the side length p of the frequency selective surface shown in fig. 4 is 3p-2w 1, i.e. 26.6 mm. When the frequency selective surface is equivalent to the second curved plate 2, the outer contour of the second curved plate 2 is made to be the same as the outer contour of the frequency selective surface, and then the thickness of the second curved plate 2 is made to be the same as the thickness of the frequency selective surface. And the frequency-selective radome is equivalent to the medium model 0, the medium model 0 is analyzed to obtain an outer surface electric field and an outer surface magnetic field of the medium model 0, and then the outer surface electric field and the outer surface magnetic field are subjected to integral operation to obtain an electric field and a magnetic field at any point of a space, so that the electrical performance analysis of the frequency-selective radome is realized, the electrical performance analysis precision is ensured, and the modeling and calculation time of the frequency-selective radome is greatly shortened.
In some alternative embodiments, with continued reference to FIG. 1, the antenna aperture face magnetic flow is derived from the antenna aperture face electric field, calculated as follows:
M＝E×n；
Wherein M is the magnetic current of the antenna aperture surface, E is the electric field of the antenna aperture surface, and n is the normal vector of the antenna aperture surface.
It can be understood that the electric field of the aperture surface of the antenna is a readily available known parameter, and the magnetic flow of the aperture surface of the antenna is calculated through the known parameter, so that the calculation process is simple.
In some alternative embodiments, with continued reference to fig. 1, the antenna aperture plane current is obtained from the antenna aperture plane magnetic field, calculated as follows:
J＝n×H；
wherein J is the current of the aperture plane of the antenna, n is the normal vector of the aperture plane of the antenna, and H is the magnetic field of the aperture plane of the antenna.
It will be appreciated that the antenna aperture plane magnetic field is also a readily available known parameter by which the antenna aperture plane current is calculated, and the calculation process is relatively simple.
In some alternative embodiments, with continued reference to FIG. 1, deriving the internal surface electric field and the internal surface magnetic field of the dielectric model 0 from the antenna aperture face magnetic current and the antenna aperture face current sum includes:
obtaining a magnetic vector according to the current of the antenna aperture surface, and obtaining an electric vector according to the magnetic current of the antenna aperture surface;
the internal surface electric field and the internal surface magnetic field are obtained according to the magnetic sagittal position and the electric sagittal position.
Specifically, the magnetic vector position is obtained according to the antenna aperture plane current, and is calculated according to the following mode:
Wherein A is magnetic vector, mu is magnetic permeability, S 1 is antenna aperture plane, J is antenna aperture plane current, J is imaginary number, k is free space wave vector, and R is distance from point P to antenna aperture plane.
Specifically, an electric vector position is obtained according to the magnetic current of the antenna aperture surface, and is calculated according to the following mode:
Wherein F is the electric vector position, epsilon is the dielectric constant, S 1 is the antenna aperture plane, M is the magnetic current of the antenna aperture plane, j is the imaginary number, k is the free space wave vector, and R is the distance from the point P to the antenna aperture plane.
It will be appreciated that the permeability and permittivity can be obtained by analysing the equivalent dielectric surface of the frequency selective surface, i.e. the second curved plate member 2, comprising:
Obtaining a normal incidence transmission coefficient and a normal incidence reflection coefficient of electromagnetic waves perpendicularly entering the second curved plate 2;
obtaining the wave impedance and complex refractive index of the second curved plate 2 according to the normal incidence transmission coefficient and the normal incidence reflection coefficient;
obtaining the magnetic permeability and the dielectric constant of the second curved plate 2 according to the wave impedance and the complex refractive index;
wherein the normal incidence transmission coefficient and the normal incidence reflection coefficient are both known parameters which are readily available, and the wave impedance and the complex refractive index of the second curved plate member 2 are calculated from the normal incidence transmission coefficient and the normal incidence reflection coefficient in the following manner:
Where z is the wave impedance, which takes its real non-negative value, S 11 'is the normal incidence reflection coefficient, and S 21' is the normal incidence transmission coefficient. n' is the complex refractive index, the complex impedance takes its imaginary part non-negative value, k 0 is the wave vector of free space, and d is the thickness of the dielectric model 0.
And obtaining the magnetic permeability and the dielectric constant of the second curved plate 2 according to the wave impedance and the complex refractive index, wherein the magnetic permeability and the dielectric constant are calculated according to the following modes:
μ＝n′z；
Wherein μ is permeability, ε is permittivity.
In some alternative embodiments, with continued reference to FIG. 1, the internal surface electric field and the internal surface magnetic field are derived from the magnetic and electrical sagittal positions, calculated as follows:
Wherein E (P) is an internal surface electric field, A is a magnetic sagittal position, F is an electric sagittal position, ε is a dielectric constant, j is an imaginary number, w is an angular frequency, μ is magnetic permeability, In order to represent the rotation, H (P) is the internal surface magnetic field.
In some alternative embodiments, with continued reference to FIG. 1, deriving the outer surface electric field of media model 0 from the inner surface electric field includes:
the inner surface electric field is decomposed into an inner surface electric field horizontal component and an inner surface electric field vertical component, calculated as follows:
E’＝E⊥×T⊥+E∥×T∥；
Wherein E' is the external surface electric field, E ⊥ is the internal surface electric field vertical component, T ⊥ is the vertical polarization transmission coefficient, E ∥ is the internal surface electric field horizontal component, and T ∥ is the horizontal polarization transmission coefficient.
In some alternative embodiments, with continued reference to FIG. 1, deriving the external surface magnetic field of the media model 0 from the internal surface magnetic field includes:
the inner surface magnetic field is decomposed into an inner surface magnetic field horizontal component and an inner surface magnetic field vertical component, calculated as follows:
H’＝H⊥×T⊥+H∥×T∥；
Wherein H' is the external surface magnetic field, H ⊥ is the internal surface magnetic field vertical component, T ⊥ is the vertical polarization transmission coefficient, H ∥ is the internal surface magnetic field horizontal component, and T ∥ is the horizontal polarization transmission coefficient.
The foregoing embodiment numbers of the present invention are merely for the purpose of description, and do not represent the advantages or disadvantages of the embodiments.
In the foregoing embodiments of the present invention, the descriptions of the embodiments are emphasized, and for a portion of this disclosure that is not described in detail in this embodiment, reference is made to the related descriptions of other embodiments.
In some alternative embodiments, the vertically polarized transmission coefficient and the horizontally polarized transmission coefficient are obtained as follows:
s11: obtaining a normal incidence transmission coefficient and a normal incidence reflection coefficient of electromagnetic waves perpendicularly entering the second curved plate 2;
s12: obtaining the wave impedance and complex refractive index of the second curved plate 2 according to the normal incidence transmission coefficient and the normal incidence reflection coefficient;
S13: obtaining the magnetic permeability and the dielectric constant of the second curved plate 2 according to the wave impedance and the complex refractive index;
s14: determining the incident angle of the electromagnetic wave according to the inner surface electric field and the inner surface magnetic field;
S15: and calculating to obtain the vertical polarization transmission coefficient and the horizontal polarization transmission coefficient by adopting a transmission line theory according to the incidence angle, the magnetic permeability and the dielectric constant.
It is to be understood that the specific methods of steps S11 to S13 have been described in detail in the above embodiments, and are not described in detail in this embodiment. The inner surface electric field and the inner surface magnetic field can determine the incidence angle of electromagnetic waves, and when the outer surface electric field and the outer surface magnetic field are solved according to the inner surface electric field and the inner surface magnetic field, the vertical polarization transmission coefficient and the horizontal polarization transmission coefficient under the corresponding incidence angles are needed. Under the condition that the normal incidence transmission coefficient and the normal incidence reflection coefficient are known, the solved magnetic permeability and dielectric constant are calculated by adopting a transmission line theory according to the incidence angle, the magnetic permeability and the dielectric constant, and the normal polarization transmission coefficient and the horizontal polarization transmission coefficient are obtained under the incidence angle. After the outer surface electric field and the outer surface magnetic field are solved according to the inner surface electric field and the inner surface magnetic field, the electric field and the magnetic field at any point of the space are obtained by performing integral operation on the outer surface electric field and the outer surface magnetic field in the same way as the calculation mode of obtaining the inner surface electric field and the inner surface magnetic field according to the magnetic sagittal position and the electric sagittal position in the embodiment, that is, the electric field at any point of the solving space needs to be calculated by using the outer surface electric field and the outer surface magnetic field, and the magnetic field at any point of the solving space needs to be calculated by using the outer surface electric field and the outer surface magnetic field.
In some alternative embodiments, referring to fig. 5, fig. 5 is an antenna far-field pattern, where the antenna far-field pattern is obtained by the frequency-selective radome analysis method based on the physical optical method and the equivalent medium model provided by the present invention.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention.
Claims (8)
1. A frequency-selective radome analysis method based on a physical optical method and an equivalent medium model is characterized by comprising the following steps of:
The frequency-selective radome is equivalent to a medium model, the medium model has the same size and shape as the frequency-selective radome, the medium model comprises a first curved surface plate, a second curved surface plate and a third curved surface plate which are laminated, and the thickness of the second curved surface plate is the same as that of the frequency-selective surface;
Acquiring an antenna aperture surface electric field and an antenna aperture surface magnetic field;
Obtaining antenna aperture surface magnetic current according to the antenna aperture surface electric field, and obtaining antenna aperture surface current according to the antenna aperture surface magnetic field;
Obtaining an inner surface electric field and an inner surface magnetic field of the medium model according to the antenna aperture surface magnetic current and the antenna aperture surface current;
Obtaining an outer surface electric field of the medium model according to the inner surface electric field, wherein the method comprises the following steps:
the inner surface electric field is decomposed into an inner surface electric field horizontal component and an inner surface electric field vertical component, calculated as follows:
E’＝E⊥×T⊥+E∥×T∥；
wherein E' is the outer surface electric field, E ⊥ is the inner surface electric field vertical component, T ⊥ is the vertical polarization transmission coefficient, E ∥ is the inner surface electric field horizontal component, and T ∥ is the horizontal polarization transmission coefficient;
obtaining an external surface magnetic field of the medium model according to the internal surface magnetic field, wherein the external surface magnetic field comprises the following steps:
decomposing the internal surface magnetic field into an internal surface magnetic field horizontal component and an internal surface magnetic field vertical component, calculated as follows:
H’＝H⊥×T⊥+H∥×T∥；
wherein H' is the outer surface magnetic field, H ⊥ is the inner surface magnetic field vertical component, and H ∥ is the inner surface magnetic field horizontal component;
And carrying out integral operation on the outer surface electric field and the outer surface magnetic field to obtain an electric field and a magnetic field at any point in space.
2. The method for analyzing the frequency selective radome based on the physical optical method and the equivalent medium model according to claim 1, wherein the antenna aperture plane magnetic flow is obtained according to the antenna aperture plane electric field, and is calculated according to the following mode:
M＝E×n；
wherein M is the magnetic current of the aperture surface of the antenna, E is the electric field of the aperture surface of the antenna, and n is the normal vector of the aperture surface of the antenna.
3. The method for analyzing the frequency selective radome based on the physical optical method and the equivalent medium model according to claim 1, wherein the antenna aperture plane current is obtained according to the antenna aperture plane magnetic field, and is calculated according to the following mode:
J＝n×H；
Wherein J is the current of the aperture plane of the antenna, n is the normal vector of the aperture plane of the antenna, and H is the magnetic field of the aperture plane of the antenna.
4. The method for analyzing a frequency selective radome based on a physical optical method and an equivalent medium model according to claim 1, wherein the obtaining an inner surface electric field and an inner surface magnetic field of the medium model according to the antenna aperture surface magnetic current and the antenna aperture surface current sum comprises:
obtaining a magnetic vector according to the current of the antenna aperture surface, and obtaining an electric vector according to the magnetic current of the antenna aperture surface;
and obtaining the inner surface electric field and the inner surface magnetic field according to the magnetic sagittal position and the electric sagittal position.
5. The method for analyzing the frequency selective radome based on the physical optical method and the equivalent medium model according to claim 4, wherein the magnetic vector position is obtained according to the antenna aperture plane current, and is calculated according to the following mode:
Wherein A is the magnetic vector, mu is magnetic permeability, S 1 is an antenna aperture plane, J is the antenna aperture plane current, J is an imaginary number, k is a free space wave vector, and R is the distance from a point P to the antenna aperture plane.
6. The method for analyzing the frequency selective radome based on the physical optical method and the equivalent medium model according to claim 4, wherein the electric vector position is obtained according to the magnetic current of the antenna aperture surface, and is calculated according to the following mode:
Wherein F is the electrical vector, epsilon is the dielectric constant, S 1 is the antenna aperture plane, M is the magnetic flow of the antenna aperture plane, j is the imaginary number, k is the free space wave vector, and R is the distance from the point P to the antenna aperture plane.
7. The method for frequency selective radome analysis based on physical optics and equivalent medium models according to claim 4, wherein the internal surface electric field and the internal surface magnetic field are obtained according to the magnetic sagittal position and the electric sagittal position, and are calculated according to the following modes:
Wherein E (P) is the inner surface electric field, A is the magnetic sagittal plane, F is the electric sagittal plane, ε is the dielectric constant, j is the imaginary number, w is the angular frequency, μ is the magnetic permeability, To represent the sign of rotation, H (P) is the internal surface magnetic field.
8. The method for analyzing a frequency selective radome based on a physical optical method and an equivalent medium model according to claim 1, wherein the vertical polarization transmission coefficient and the horizontal polarization transmission coefficient are obtained in the following manner:
obtaining a normal incidence transmission coefficient and a normal incidence reflection coefficient of electromagnetic waves perpendicularly incident on the second curved plate;
obtaining the wave impedance and complex refractive index of the second curved plate according to the normal incidence transmission coefficient and the normal incidence reflection coefficient;
Obtaining the magnetic permeability and the dielectric constant of the second curved plate according to the wave impedance and the complex refractive index;
determining an incident angle of electromagnetic waves according to the inner surface electric field and the inner surface magnetic field;
And calculating the vertical polarization transmission coefficient and the horizontal polarization transmission coefficient by adopting a transmission line theory according to the incidence angle, the magnetic permeability and the dielectric constant.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310862769.3A CN117057098B (en) | 2023-07-13 | 2023-07-13 | Frequency-selective radome analysis method based on physical optical method and equivalent medium model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310862769.3A CN117057098B (en) | 2023-07-13 | 2023-07-13 | Frequency-selective radome analysis method based on physical optical method and equivalent medium model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117057098A CN117057098A (en) | 2023-11-14 |
| CN117057098B true CN117057098B (en) | 2024-04-26 |
Family
ID=88668196
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310862769.3A Active CN117057098B (en) | 2023-07-13 | 2023-07-13 | Frequency-selective radome analysis method based on physical optical method and equivalent medium model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117057098B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107391880A (en) * | 2017-08-18 | 2017-11-24 | 上海无线电设备研究所 | A kind of accurate efficient radome electrical property energy analysis method |
| CN108268677A (en) * | 2016-12-30 | 2018-07-10 | 深圳光启高等理工研究院 | The electrical performance evaluation method and device of antenna house |
| CN110059422A (en) * | 2019-04-23 | 2019-07-26 | 北京环境特性研究所 | The Electromagnetic Scattering Characteristics emulation mode of frequency-selective surfaces curve antenna cover |
| KR20200056134A (en) * | 2018-11-14 | 2020-05-22 | 아주대학교산학협력단 | Method and apparatus for analyzing radar cross section (RCS) of complex structure |
| CN113609646A (en) * | 2021-07-08 | 2021-11-05 | 中国人民解放军32215部队 | Modeling and simulation method for coupling electromagnetic scattering characteristics of complex land environment and equipment |
-
2023
- 2023-07-13 CN CN202310862769.3A patent/CN117057098B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108268677A (en) * | 2016-12-30 | 2018-07-10 | 深圳光启高等理工研究院 | The electrical performance evaluation method and device of antenna house |
| CN107391880A (en) * | 2017-08-18 | 2017-11-24 | 上海无线电设备研究所 | A kind of accurate efficient radome electrical property energy analysis method |
| KR20200056134A (en) * | 2018-11-14 | 2020-05-22 | 아주대학교산학협력단 | Method and apparatus for analyzing radar cross section (RCS) of complex structure |
| CN110059422A (en) * | 2019-04-23 | 2019-07-26 | 北京环境特性研究所 | The Electromagnetic Scattering Characteristics emulation mode of frequency-selective surfaces curve antenna cover |
| CN113609646A (en) * | 2021-07-08 | 2021-11-05 | 中国人民解放军32215部队 | Modeling and simulation method for coupling electromagnetic scattering characteristics of complex land environment and equipment |
Non-Patent Citations (2)
| Title |
|---|
| AI-SI法分析X波段机载雷达天线罩;史维光;王建;张云祥;罗天光;郑俊;;现代雷达;20140215(第02期);72-75 * |
| 基于等效电路分析方法的微型频率选择表面设计;陈丙根;;火控雷达技术;20180925(第03期);87-90+95 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117057098A (en) | 2023-11-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109374985A (en) | Electromagnetic environment monitor method, system and storage medium | |
| CN106252873B (en) | A kind of Novel Interval Methods of conformal bearer Antenna Far Field power radiation pattern | |
| CN107527035A (en) | Earthquake damage to building information extracting method and device | |
| CN117057098B (en) | Frequency-selective radome analysis method based on physical optical method and equivalent medium model | |
| CN112818762B (en) | Large-size composite material and rapid nondestructive testing method for sandwich structure thereof | |
| Li et al. | Research on an improved metal surface defect detection sensor based on a 3D RFID tag antenna | |
| CN103852758B (en) | The on-the-spot method for rapidly judging of invisbile plane scattering properties and device | |
| Jeong et al. | Analysis of RCS of low observable aircraft in VHF band | |
| Fang et al. | Visualization and quantitative evaluation of delamination defects in GFRPs via sparse millimeter-wave imaging and image processing | |
| Michalski | Complex image method analysis of a plane wave-excited subwavelength circular aperture in a planar screen | |
| CN112858794B (en) | System and method for measuring scattering parameters of ultrathin reflection-type electromagnetic super surface | |
| Tian et al. | The analysis and measurement of FSS radome for antenna RCS reduction | |
| Monorchio et al. | A hybrid mode-matching finite-elements approach to the analysis of thick dichroic screens with arbitrarily shaped apertures | |
| CN113177307A (en) | Short-wave ground wave simulation method and system applied to complex terrain | |
| Datta et al. | Subwavelength Microwave Imaging System using a Negative Index Metamaterial Lens | |
| CN116400320B (en) | Sea fog effective particle radius estimation method based on laser and W-band radar | |
| CN112462153B (en) | Fresnel coefficient-based method for estimating radiation characteristics of antennas in dielectric plate | |
| Zeyde et al. | A free-space samples material parameters validation technique | |
| CN110907912A (en) | Pole characteristic solving method in target narrow-band mode of resonance region | |
| Din et al. | Rotation included 3-axis scanning free-space measurement and curvature compensation for electromagnetic evaluation of leading-edge and curved stealth structures | |
| CN114076924B (en) | Method for realizing large-scale complex target classification recognition based on multi-frequency echo data | |
| CN109117594A (en) | A kind of submarine mast searching probability acquisition methods based on radar wave | |
| CN116092072B (en) | Spacecraft target detection method, spacecraft target detection system, storage medium and electronic equipment | |
| CN110210177B (en) | Method and system for determining frequency selective surface structure scattered field | |
| CN116299239B (en) | Radar scattering sectional area error optimization method for near-far field transformation |
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 |