CN107391880A - A kind of accurate efficient radome electrical property energy analysis method - Google Patents
A kind of accurate efficient radome electrical property energy analysis method Download PDFInfo
- Publication number
- CN107391880A CN107391880A CN201710712133.5A CN201710712133A CN107391880A CN 107391880 A CN107391880 A CN 107391880A CN 201710712133 A CN201710712133 A CN 201710712133A CN 107391880 A CN107391880 A CN 107391880A
- Authority
- CN
- China
- Prior art keywords
- mrow
- antenna house
- electrical property
- msup
- msub
- 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
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/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Details Of Aerials (AREA)
Abstract
The present invention relates to a kind of accurate efficient radome electrical property energy analysis method, comprising：S1, the geometry physical three-dimensional model for establishing antenna house, and carry out mesh generation；S2, according to antenna near-field Aperture distribution, calculate the incident electromagnetic field for inciding antenna house inwall；S3, the reflectance factor and transmission coefficient for calculating antenna house surface；S4, the initial electromagnetic stream for calculating antenna house surface；S5, accelerated using iterative physical optics formula and GPU, iterate and calculate the electromagnetic current on antenna house surface, until obtaining stable surface electromagnetic current；S6, analysis calculate the electrical property of antenna house；S7, when unit for electrical property parameters is unsatisfactory for design requirement, by adjusting the parameters such as antenna house dielectric constant, antenna house wall thickness, repeatedly perform S1~S6 optimize.The present invention solves the radiation problem of antenna house cover body by the iterative physical optics of amendment, can take into account the speed of numerical method and the precision of high-frequency methods, can be designed for radome electrical property and save substantial amounts of optimization time early stage.
Description
Technical field
The present invention relates to radome electrical property energy analysis method, in particular to a kind of radome electrical property accurately efficiently to emulate
Computing technique, belong to antenna house emulation technology.
Background technology
Antenna house is the important device for protecting radar antenna system, and it is located at system front end.Antenna house not only needs full
The demand of pedal system aerodynamic configuration, while also need to the requirement for reaching the electrical property that system is proposed.The electrical property of antenna house point
Analysis belongs to electrically large sizes problem, and requires high.Existing emulated computation method can not meet wanting for accuracy and speed simultaneously
Ask.
Nanjing Aero-Space University is in patent《It is a kind of to weigh method of the antenna house to the lateral performance impact of antenna permutation》It is (public
The number of opening：CN105388449A illustrated in), each antenna electric field magnitude caused by far field outside the cover is calculated using physical optical method
And phase；To any antenna element in aerial array, the insertion phase delay of this antenna element is calculated；Then aerial array system is calculated
The phase error average value of any two antenna in system；The Sensor gain and phase perturbations of antenna house are finally calculated.
Xian Electronics Science and Technology University is in patent《A kind of Thickness Design Method of antenna house》(publication number：CN106469850A)
Middle elaboration, the thickness of antenna house is subjected to discretization according to the short transverse of cover body；Determine the span of thickness value, and to from
The thickness of scatterplot assigns initial value；The transmission coefficient of cover body is calculated with transmission line theory；Establish mathematical optimization models.And and then in patent
《A kind of quick Thickness Design Method of aircraft antenna house》(publication number：CN106654566A illustrated in), with the addition of determination day
Gain loss caused by irdome and collimating fault, judge gained radome design scheme electrical performance indexes and thickness distribution whether
Meet preset requirement.
Above three patent is all mainly that the electrical property of antenna house, physics are obtained by physioptial high frequency method
Optical method is with the obvious advantage in calculating speed, but is slightly short of in precision, it is difficult to meets high-precision demand.
In current existing open source literature and periodical, radome electrical property energy analysis method relate generally to accurate numerical methods and
High frequency Asymptotical Method.Wherein, university of Dalian University of Technology Zhang Yunshang is in Master's thesis《Antenna-radowi Systems electrical property numerical analysis
With optimizing research》Middle elaboration, the main algorithm of numerical analysis is used as using plane wave spectrum-surface integration method.Ugo d ' Elia exist
The paper delivered on PIER periodicals《A Physical Optics Approach to the Analysis of Large
Frequency Selective Radomes》Middle elaboration, having carried out radome electrical property using physical optical method can analyze.But
The calculating speed of accurate numerical methods, and the computational accuracy of high frequency Asymptotical Method can not all meet actual demand.
Therefore, the Computational Electromagnetic Methods for proposing a kind of radome electrical property energy are needed badly at present, radome electrical property can be applied to
The quick calculating and required precision of energy, it is to save substantial amounts of design time early stage to the practical conversion of achievement.
The content of the invention
It is an object of the invention to provide a kind of accurate efficient radome electrical property energy analysis method, pass through the iteration thing of amendment
The radiation problem for solving antenna house cover body is learned by Ricoh, can take into account the speed of numerical method and the precision of high-frequency methods, be antenna
Cover electrical property design saves substantial amounts of optimization time early stage.
To achieve the above object, the present invention provides a kind of accurate efficient radome electrical property energy analysis method, comprising following
Step：
S1, according to design requirement, the geometry physical three-dimensional model of antenna house is established in optimization, and carries out mesh generation to it；
S2, according to antenna near-field Aperture distribution, calculate the incident electric fields and incident magnetic for inciding antenna house inwall；
S3, the reflectance factor and transmission coefficient for calculating antenna house surface diverse location；
S4, according to incident electric fields, incident magnetic, reflectance factor and transmission coefficient, calculate the initial current on antenna house surface
With initial magnetic current, initial value is provided for successive iterations physical optics；
S5, the initial current according to antenna house surface and initial magnetic current, accelerated using iterative physical optics formula and GPU,
Iterate and calculate the electromagnetism and electric current on antenna house surface, until obtaining stable surface electromagnetism and electric current；
S6, the electrical property for analyzing antenna house：Radiation characteristic is obtained according to stable surface electromagnetism and electric current, calculates antenna house
Unit for electrical property parameters；
S7, judge whether the unit for electrical property parameters of antenna house meets design requirement；It is if not, normal by adjusting antenna house dielectric
The parameters such as number, antenna house wall thickness, S1~S6 is performed repeatedly, the unit for electrical property parameters of antenna house is optimized, until its satisfaction is set
Meter requires.
In described S1, required according to the subdivision of iterative physical optics, the several of interval subdivision antenna house are used as using λ/3~λ 4
What physical three-dimensional model, λ represent the wavelength of incident electromagnetic wave.
In described S2, the incident electric fields being calculated are：
The incident magnetic being calculated is：
Wherein, S be antenna near field distribution surface, ES、HSThe respectively electric field on the near field distribution surface of antenna and magnetic field,
For normal vector, G0For Green's function, ω, μ and ε represent the dielectric constant and magnetic conductivity of electromagnetism angular frequency, antenna house respectively,
R and r ' represents site position and source location respectively,Expression acts on the Hamiltonian of source coordinate system.
In described S3, the calculation formula of transmission coefficient is expressed as：
Wherein,For the transmission matrix of multilayered medium material；Symbol ⊥ and // represent vertical polarization and level respectively
Polarization；θ0For incidence angle；T=d1tanθ1+d2tanθ2+…+dN tanθN, d is the thickness of every layer of dielectric material；Z is wave impedance；
K is wave number.
In described S4, the initial current on antenna house surface and initial magnetic current are respectively：
Wherein, J1Represent the initial current of antenna house inner wall surface；M1Represent the initial magnetic current of antenna house inner wall surface；J2
Represent the initial current of antenna house outer wall surface；M2Represent the initial magnetic current of antenna house outer wall surface；EincRepresent antenna house surface
Incident electric fields；ErRepresent the reflected field on antenna house surface；EtRepresent the transmitted electric fields on antenna house surface；HincRepresent antenna house
The incident magnetic on surface；HrRepresent the reflection magnetic field on antenna house surface；HtRepresent the transmission magnetic field on antenna house surface.
In described S5, the iterative physical optics formula of use is expressed as：
J (r ')=J (Ηinc,r′)+J(Η(J,M),r′)；
M (r ')=M (Einc,r′)+M(E(J,M),r′)；
Wherein, J () and M () represents electric current, the action function of magnetic current respectively, and Η (J, M) and E (J, M) represent upper one respectively
Contribution of the magnetic field and electric field that secondary electromagnetic current excites to electric current and magnetic current.
In summary, accurate efficient radome electrical property energy analysis method provided by the present invention, have advantages below and
Beneficial effect：1) electrical property of calculation medium material antenna house is extended to by introducing magnetic current by iterative physical optics；2) propose
More accurate transmission coefficient calculation formula, the phase offset introduced in traditional calculating formula is compensate for, what Accelerated iteration calculated
Speed；3) interaction in GPU speed-up computation iterative physical optics between each bin is used；4) the amendment iteration that GPU accelerates
Physical optics reaches λ/3~λ/4 in subdivision size and can meet to calculate requirement, calculates and the dual acceleration of grid and iteration thing
The characteristic for being themselves based on magnetic field integral equation is learned by Ricoh, ensures the fast and accurate of radome electrical property energy optimization design.
Brief description of the drawings
Fig. 1 is the flow chart of accurate efficient radome electrical property energy analysis method in the present invention；
Fig. 2 is the amendment schematic diagram of multilayered medium material transmission phase in the present invention；
Fig. 3 is the computation model schematic diagram of antenna house surface electromagnetic current in the present invention；
Fig. 4 is the schematic diagram of near field Aperture distribution specifically used in the present invention；
Fig. 5 is the initial distribution of antenna house surface current and the contrast schematic diagram of steady-state distribution in the present invention.
Embodiment
Below in conjunction with Fig. 1~Fig. 5, a preferred embodiment of the present invention is described in detail.
As shown in figure 1, for accurate efficient radome electrical property energy analysis method provided by the present invention, following step is included
Suddenly：
S1, according to design requirement, the geometry physical three-dimensional model of antenna house is established in optimization, and to the geometry physical three-dimensional mould
Type carries out mesh generation；
S2, according to antenna near-field Aperture distribution, calculate the incident electric fields and incident magnetic for inciding antenna house inwall；
S3, Fresnel (Fresnel) reflection Rs and transmission coefficient t for calculating antenna house surface diverse location；
S4, according to incident electric fields, incident magnetic, reflectance factor and transmission coefficient, calculate the initial current on antenna house surface
With initial magnetic current, initial value is provided for successive iterations physical optics；
S5, the initial current according to antenna house surface and initial magnetic current, accelerated using iterative physical optics formula and GPU
(being carried out using graphics processing unit GPU hardware-accelerated), is iterated and calculates the electromagnetism and electric current on antenna house surface, until table
Face electromagnetism and electric current reach surface electromagnetism and electric current stable, and that record finally stable；
S6, the electrical property for analyzing antenna house：Radiation characteristic is obtained according to stable surface electromagnetism and electric current, calculates antenna house
Unit for electrical property parameters；
S7, judge whether the unit for electrical property parameters of antenna house meets design requirement；It is if not, normal by adjusting antenna house dielectric
The parameters such as number, antenna house wall thickness, S1~S6 is performed repeatedly, the unit for electrical property parameters of antenna house is optimized, until its satisfaction is set
Meter requires.
In described S1, required according to the subdivision of iterative physical optics, be typically used as interval subdivision antenna using λ/3~λ/4
The geometry physical three-dimensional model of cover, so as to meet the calculating requirement of antenna house, while reach the purpose for reducing amount of calculation；Wherein, λ
Represent the wavelength of incident electromagnetic wave.
As shown in figure 4, in described S2, the incident electric fields being calculated are：
The incident magnetic being calculated is：
Wherein, S be antenna near field distribution surface, ES、HSThe respectively electric field on the near field distribution surface of antenna and magnetic field,For normal vector, G0For Green's function, ω, μ and ε represent electromagnetism angular frequency, the dielectric constant of antenna house and magnetic conductance respectively
Rate, r and r ' represent site position and source location respectively,Expression acts on the Hamiltonian of source coordinate system.
In described S3, in order to more be precisely calculated transmission coefficient t, it is proposed that a kind of more accurate transmission coefficient meter
Formula is calculated, enables to the calculating of surface electromagnetic current more accurate；The calculation formula of transmission coefficient t is expressed as：
Wherein,For the transmission matrix of multilayered medium material；Symbol ⊥ and // represent vertical polarization and level respectively
Polarization；θ0For incidence angle；T=d1tanθ1+d2tanθ2+…+dN tanθN, d is the thickness of every layer of dielectric material；Z is wave impedance；
K is wave number, specific as shown in Figure 1.
As shown in figure 3, in described S4, the initial current on antenna house surface and initial magnetic current are respectively：
Wherein, J1Represent the initial current of antenna house inner wall surface；M1Represent the initial magnetic current of antenna house inner wall surface；J2
Represent the initial current of antenna house outer wall surface；M2Represent the initial magnetic current of antenna house outer wall surface；EincRepresent antenna house surface
Incident electric fields；ErRepresent the reflected field on antenna house surface；EtRepresent the transmitted electric fields on antenna house surface；HincRepresent antenna house
The incident magnetic on surface；HrRepresent the reflection magnetic field on antenna house surface；HtRepresent the transmission magnetic field on antenna house surface.
In described S5, the iterative physical optics formula of use is expressed as：
J (r ')=J (Ηinc,r′)+J(Η(J,M),r′)；
M (r ')=M (Einc,r′)+M(E(J,M),r′)；
Wherein, J () and M () represents electric current, the action function of magnetic current respectively, and Η (J, M) and E (J, M) represent upper one respectively
Contribution of the magnetic field and electric field that secondary electromagnetic current excites to electric current and magnetic current.
As shown in figure 5, being the contrast schematic diagram of antenna house Surface current distribution, wherein initial distribution only embodies direct photograph
The effect penetrated, and steady-state distribution by successive ignition due to just having embodied mirror effect.
In summary, a kind of accurate efficient radome electrical property energy analysis method proposed by the present invention, is one specifically
The computational methods of kind fast optimal design radome electrical property energy parameter, the double requirements of accuracy and speed can be kept.This method
The accurate initial electromagnetic stream in antenna house cover body surface is obtained using accurate transmission coefficient first, initial ginseng is provided for successive iterations
Number；Then traditional iterative physical optics are extended to the calculating of dielectric material antenna house, cover body surface is only existed by original
Current expansion to surface current and all existing situation of surface magnetic current；Accelerate to iterate finally by GPU and calculate cover body table
The electromagnetic current in face, obtain the surface electromagnetic current that cover body surface is finally stablized；And then the unit for electrical property parameters of antenna house is obtained, and can lead to
The parameters such as adjustment dielectric constant, wall thickness are crossed to continue to optimize until meeting design requirement.
The inventive method solves the radiation problem of antenna house cover body by the iterative physical optics of amendment, can take into account numerical value
The speed of method and the precision of high-frequency methods, substantial amounts of optimization time early stage is saved so as to be designed for radome electrical property.Specifically
For, there is advantages below and beneficial effect：
1) use range of iterative physical optics has been promoted.Iterative physical optics can only originally calculate metal, coating problem,
The electrical property of calculation medium material antenna house is extended to by introducing magnetic current by iterative physical optics；
2) in order to more be precisely calculated transmission coefficient, it is proposed that a kind of more accurate transmission coefficient calculation formula, more
The phase offset introduced in traditional calculating formula, the speed that Accelerated iteration calculates are mended；
3) interaction in GPU speed-up computation iterative physical optics between each bin is used.Due to phase between each bin
The calculation formula that interaction uses is identical, few logic judgments, is applicable very much GPU parallel computations；
4) the amendment iterative physical optics that GPU accelerates reach λ/3~λ/4 in subdivision size and can meet to calculate requirement,
The characteristic that magnetic field integral equation is themselves based on the dual acceleration of grid and iterative physical optics is calculated, ensures radome electrical property energy
Optimization design it is fast and accurate.
Although present disclosure is discussed in detail by above preferred embodiment, but it should be appreciated that above-mentioned
Description is not considered as limitation of the present invention.After those skilled in the art have read the above, for the present invention's
A variety of modifications and substitutions all will be apparent.Therefore, protection scope of the present invention should be limited to the appended claims.
Claims (6)
1. a kind of accurate efficient radome electrical property energy analysis method, it is characterised in that comprise the steps of：
S1, according to design requirement, the geometry physical three-dimensional model of antenna house is established in optimization, and carries out mesh generation to it；
S2, according to antenna near-field Aperture distribution, calculate the incident electric fields and incident magnetic for inciding antenna house inwall；
S3, the reflectance factor and transmission coefficient for calculating antenna house surface diverse location；
S4, according to incident electric fields, incident magnetic, reflectance factor and transmission coefficient, calculate the initial current and just on antenna house surface
Beginning magnetic current, initial value is provided for successive iterations physical optics；
S5, the initial current according to antenna house surface and initial magnetic current, accelerated using iterative physical optics formula and GPU, repeatedly
The electromagnetism and electric current on antenna house surface are iterated to calculate, until obtaining stable surface electromagnetism and electric current；
S6, the electrical property for analyzing antenna house：Radiation characteristic is obtained according to stable surface electromagnetism and electric current, calculates the electricity of antenna house
Performance parameter；
S7, judge whether the unit for electrical property parameters of antenna house meets design requirement；If not, by adjust antenna house dielectric constant,
The parameters such as antenna house wall thickness, S1~S6 is performed repeatedly, the unit for electrical property parameters of antenna house is optimized, until its satisfaction design will
Ask.
2. accurate efficient radome electrical property energy analysis method as claimed in claim 1, it is characterised in that in described S1,
Required according to the subdivision of iterative physical optics, the geometry physical three-dimensional model of interval subdivision antenna house, λ generations are used as using λ/3~λ/4
The wavelength of table incident electromagnetic wave.
3. accurate efficient radome electrical property energy analysis method as claimed in claim 1, it is characterised in that in described S2,
The incident electric fields being calculated are：
The incident magnetic being calculated is：
Wherein, S be antenna near field distribution surface, ES、HSThe respectively electric field on the near field distribution surface of antenna and magnetic field,For method
To vector, G0For Green's function, ω, μ and ε represent the dielectric constant and magnetic conductivity of electromagnetism angular frequency, antenna house respectively, r and
R ' represents site position and source location respectively,Expression acts on the Hamiltonian of source coordinate system.
4. accurate efficient radome electrical property energy analysis method as claimed in claim 3, it is characterised in that in described S3,
The calculation formula of transmission coefficient is expressed as：
<mrow>
<msub>
<mi>T</mi>
<mo>&perp;</mo>
</msub>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>+</mo>
<mi>B</mi>
<mo>/</mo>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>&perp;</mo>
</msubsup>
<mo>)</mo>
<mo>+</mo>
<msubsup>
<mi>Z</mi>
<mn>0</mn>
<mo>&perp;</mo>
</msubsup>
<mo>(</mo>
<mi>C</mi>
<mo>+</mo>
<mi>D</mi>
<mo>/</mo>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>&perp;</mo>
</msubsup>
<mo>)</mo>
</mrow>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>jksin&theta;</mi>
<mn>0</mn>
</msub>
<mi>t</mi>
</mrow>
</msup>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mo>/</mo>
<mo>/</mo>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>+</mo>
<mi>B</mi>
<mo>/</mo>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>/</mo>
<mo>/</mo>
</mrow>
</msubsup>
<mo>)</mo>
<mo>+</mo>
<msubsup>
<mi>Z</mi>
<mn>0</mn>
<mrow>
<mo>/</mo>
<mo>/</mo>
</mrow>
</msubsup>
<mo>(</mo>
<mi>C</mi>
<mo>+</mo>
<mi>D</mi>
<mo>/</mo>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>/</mo>
<mo>/</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>jksin&theta;</mi>
<mn>0</mn>
</msub>
<mi>t</mi>
</mrow>
</msup>
<mo>;</mo>
</mrow>
Wherein,For the transmission matrix of multilayered medium material；Symbol ⊥ and // represent vertical polarization and horizontal polarization respectively；
θ0For incidence angle；T=d1tanθ1+d2tanθ2+…+dNtanθN, d is the thickness of every layer of dielectric material；Z is wave impedance；K is ripple
Number.
5. accurate efficient radome electrical property energy analysis method as claimed in claim 4, it is characterised in that in described S4,
The initial current on antenna house surface and initial magnetic current are respectively：
<mrow>
<msub>
<mi>J</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>E</mi>
<mi>r</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
1
<mrow>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>H</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>H</mi>
<mi>r</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&times;</mo>
<msup>
<mi>E</mi>
<mi>t</mi>
</msup>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mover>
<mi>n</mi>
<mo>^</mo>
</mover>
<mo>&times;</mo>
<msup>
<mi>H</mi>
<mi>t</mi>
</msup>
<mo>;</mo>
</mrow>
Wherein, J1Represent the initial current of antenna house inner wall surface；M1Represent the initial magnetic current of antenna house inner wall surface；J2Represent
The initial current of antenna house outer wall surface；M2Represent the initial magnetic current of antenna house outer wall surface；EincRepresent entering for antenna house surface
Radio field；ErRepresent the reflected field on antenna house surface；EtRepresent the transmitted electric fields on antenna house surface；HincRepresent antenna house surface
Incident magnetic；HrRepresent the reflection magnetic field on antenna house surface；HtRepresent the transmission magnetic field on antenna house surface.
6. accurate efficient radome electrical property energy analysis method as claimed in claim 5, it is characterised in that in described S5,
The iterative physical optics formula of use is expressed as：
J (r ')=J (Ηinc,r′)+J(Η(J,M),r′)；
M (r ')=M (Einc,r′)+M(E(J,M),r′)；
Wherein, J () and M () represents electric current, the action function of magnetic current respectively, and Η (J, M) and E (J, M) represent last electricity respectively
Contribution of the magnetic field and electric field that magnetic current excites to electric current and magnetic current.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710712133.5A CN107391880B (en) | 2017-08-18 | 2017-08-18 | Accurate and efficient antenna housing electrical property analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710712133.5A CN107391880B (en) | 2017-08-18 | 2017-08-18 | Accurate and efficient antenna housing electrical property analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107391880A true CN107391880A (en) | 2017-11-24 |
| CN107391880B CN107391880B (en) | 2020-07-31 |
Family
ID=60352813
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710712133.5A Active CN107391880B (en) | 2017-08-18 | 2017-08-18 | Accurate and efficient antenna housing electrical property analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107391880B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108920831A (en) * | 2018-06-29 | 2018-11-30 | 西安电子科技大学 | The quick calculation method that high-speed aircraft cover high temperature ablation influences antenna electric performance |
| CN110380208A (en) * | 2019-07-03 | 2019-10-25 | 惠州市德赛西威智能交通技术研究院有限公司 | A kind of double arc millimetre-wave radar antenna houses of Varying-thickness and design method |
| CN110738002A (en) * | 2019-10-18 | 2020-01-31 | 上海无线电设备研究所 | Comprehensive performance design method for microwave/infrared composite antenna housing |
| CN113625062A (en) * | 2021-07-29 | 2021-11-09 | 西安电子科技大学 | Antenna housing electrical property estimation method based on Taylor expansion method |
| CN116484522A (en) * | 2023-04-13 | 2023-07-25 | 大连理工大学 | Radome electrical performance compensation processing amount determining method based on cascade network |
| CN117057098A (en) * | 2023-07-13 | 2023-11-14 | 大连理工大学 | Frequency-selective radome analysis method based on physical optical method and equivalent medium model |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100066595A1 (en) * | 2008-09-18 | 2010-03-18 | Lee Chul J | Electromagnetic (em) solver using a shooting bouncing ray (sbr) technique |
| CN102590656A (en) * | 2012-01-03 | 2012-07-18 | 西安电子科技大学 | Antenna cover electric property forecasting method based on distant field |
| CN103870654A (en) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing |
| CN104750960A (en) * | 2013-12-30 | 2015-07-01 | 南京理工大学 | Method for rapidly extracting electrical property parameter of metal truss-type radome |
| CN104794262A (en) * | 2015-03-30 | 2015-07-22 | 西安电子科技大学 | Electromechanical-integration-based shaped beam design method for satellite-borne mesh antenna |
| CN105388449A (en) * | 2015-11-27 | 2016-03-09 | 南京航空航天大学 | Method of measuring influence of radome on antenna array direction-finding performance |
| CN106772301A (en) * | 2016-12-02 | 2017-05-31 | 上海无线电设备研究所 | A kind of multilayer non-parallel interfaces medium electromagnetic scattering simulation algorithm based on medium ray tracing |
| CN106934097A (en) * | 2017-02-09 | 2017-07-07 | 西安电子科技大学 | Towards the spatial networks antenna key dynamics mode choosing method of electrical property |
-
2017
- 2017-08-18 CN CN201710712133.5A patent/CN107391880B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100066595A1 (en) * | 2008-09-18 | 2010-03-18 | Lee Chul J | Electromagnetic (em) solver using a shooting bouncing ray (sbr) technique |
| CN102590656A (en) * | 2012-01-03 | 2012-07-18 | 西安电子科技大学 | Antenna cover electric property forecasting method based on distant field |
| CN104750960A (en) * | 2013-12-30 | 2015-07-01 | 南京理工大学 | Method for rapidly extracting electrical property parameter of metal truss-type radome |
| CN103870654A (en) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing |
| CN104794262A (en) * | 2015-03-30 | 2015-07-22 | 西安电子科技大学 | Electromechanical-integration-based shaped beam design method for satellite-borne mesh antenna |
| CN105388449A (en) * | 2015-11-27 | 2016-03-09 | 南京航空航天大学 | Method of measuring influence of radome on antenna array direction-finding performance |
| CN106772301A (en) * | 2016-12-02 | 2017-05-31 | 上海无线电设备研究所 | A kind of multilayer non-parallel interfaces medium electromagnetic scattering simulation algorithm based on medium ray tracing |
| CN106934097A (en) * | 2017-02-09 | 2017-07-07 | 西安电子科技大学 | Towards the spatial networks antenna key dynamics mode choosing method of electrical property |
Non-Patent Citations (4)
| Title |
|---|
| JI HYUNG KIM ET AL: ""Analysis of FSS Radomes Based on Physical Optics Method and Ray Tracing Technique"", 《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》 * |
| UGO D"ELIA ET AL: ""A PHYSICAL OPTICS APPROACH TO THE ANALYSIS OF LARGE FREQUENCY SELECTIVE RADOMES"", 《PROGRESS IN ELECTROMAGNETICS RESEARCH》 * |
| 隋淼 等: ""基于曲面片迭代物理光学法的电磁耦合计算"", 《电波科学学报》 * |
| 顾长青 等: ""一种改进的物理光学迭代法"", 《电子与信息学报》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108920831A (en) * | 2018-06-29 | 2018-11-30 | 西安电子科技大学 | The quick calculation method that high-speed aircraft cover high temperature ablation influences antenna electric performance |
| CN110380208A (en) * | 2019-07-03 | 2019-10-25 | 惠州市德赛西威智能交通技术研究院有限公司 | A kind of double arc millimetre-wave radar antenna houses of Varying-thickness and design method |
| CN110380208B (en) * | 2019-07-03 | 2021-02-05 | 惠州市德赛西威智能交通技术研究院有限公司 | Variable-thickness double-arc millimeter wave radome and design method |
| CN110738002A (en) * | 2019-10-18 | 2020-01-31 | 上海无线电设备研究所 | Comprehensive performance design method for microwave/infrared composite antenna housing |
| CN113625062A (en) * | 2021-07-29 | 2021-11-09 | 西安电子科技大学 | Antenna housing electrical property estimation method based on Taylor expansion method |
| CN113625062B (en) * | 2021-07-29 | 2022-05-17 | 西安电子科技大学 | Antenna housing electrical property estimation method based on Taylor expansion method |
| CN116484522A (en) * | 2023-04-13 | 2023-07-25 | 大连理工大学 | Radome electrical performance compensation processing amount determining method based on cascade network |
| CN116484522B (en) * | 2023-04-13 | 2024-04-19 | 大连理工大学 | Radome electrical performance compensation processing amount determining method based on cascade network |
| CN117057098A (en) * | 2023-07-13 | 2023-11-14 | 大连理工大学 | Frequency-selective radome analysis method based on physical optical method and equivalent medium model |
| CN117057098B (en) * | 2023-07-13 | 2024-04-26 | 大连理工大学 | Frequency-selective radome analysis method based on physical optical method and equivalent medium model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107391880B (en) | 2020-07-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107391880A (en) | A kind of accurate efficient radome electrical property energy analysis method | |
| CN112989680B (en) | FVFD far-field integration boundary condition calculation method for reducing grid usage | |
| Huang et al. | Electromagnetic scattering of randomly rough soil surfaces based on numerical solutions of Maxwell equations in three-dimensional simulations using a hybrid UV/PBTG/SMCG method | |
| CN105160115B (en) | Based on the approximate reflector antenna electromechanics Integrated Optimal Design method with sensitivity analysis | |
| CN113158527B (en) | Method for calculating frequency domain electromagnetic field based on implicit FVFD | |
| Chen et al. | Improved scattering-matrix method and its application to analysis of electromagnetic wave reflected by reentry plasma sheath | |
| CN105468909A (en) | Time delay power system electromechanical oscillation mode computing method based on SOD-PS-R R algorithm | |
| CN105302962A (en) | Electromechanical integration optimization design method of reflector antenna on the basis of structure-electromagnetic hybrid unit | |
| CN105786765A (en) | Method for generating incentive irrelevant characteristic basis function rapidly in self-adaption mode | |
| Choudhari et al. | Computations of crossflow instability in hypersonic boundary layers | |
| CN104992035A (en) | Quick calculation method for terahertz band surface rough target electromagnetic scattering | |
| Song et al. | Accurate evaluation of Green's functions for a lossy layered medium by fast extraction of surface-and leaky-wave modes | |
| CN106528954A (en) | Electromechanical integration optimization design method for space mesh antenna under action of temperature load | |
| Zhou et al. | Efficient EDM-PO method for the scattering from electrically large objects with the high-order impedance boundary condition | |
| CN106485071B (en) | A kind of quick far to field transformation method of multilayer packet configuration | |
| CN105549031B (en) | A kind of time domain numerical computation method of the ionospheric propagation time delay of satellite-signal | |
| Dyrud et al. | Plasma and electromagnetic wave simulations of meteors | |
| Koh et al. | New scaling factors of 2-D isotropic-dispersion finite difference time domain (ID-FDTD) algorithm for lossy media | |
| Fan et al. | The efficient analysis of electromagnetic scattering from composite structures using hybrid CFIE-IEFIE | |
| Franek et al. | Application of numerical dispersion compensation of the Yee-FDTD algorithm on elongated domains | |
| Yin et al. | Study on calculation and verification of radiowave propagation using parabolic equation for the antenna near the ground | |
| Butler et al. | Towards a fast and precise forward model for air shower radio simulation | |
| Sun et al. | Analysis of scattering from composite conductor and dielectric objects using single integral equation method and MLFMA based on JMCFIE | |
| Ding et al. | Analysis of electromagnetic scattering of cavities with nonuniform plasma coating | |
| Lipatov et al. | The light (H+, H2+, He+) and heavy (Na+) pickup ion dynamics in the lunar-like plasma environment: 3D hybrid kinetic modeling |
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 |