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
 238000004458 analytical method Methods 0.000 title claims abstract description 22
 230000005540 biological transmission Effects 0.000 claims abstract description 25
 238000005457 optimization Methods 0.000 claims abstract description 9
 230000005684 electric field Effects 0.000 claims description 21
 238000004364 calculation method Methods 0.000 claims description 9
 239000000463 material Substances 0.000 claims description 6
 239000003989 dielectric material Substances 0.000 claims description 4
 239000011159 matrix material Substances 0.000 claims description 3
 230000005611 electricity Effects 0.000 claims 2
 230000005672 electromagnetic field Effects 0.000 abstract 1
 238000010586 diagram Methods 0.000 description 5
 230000000694 effects Effects 0.000 description 4
 238000005516 engineering process Methods 0.000 description 3
 230000003993 interaction Effects 0.000 description 3
 230000003287 optical Effects 0.000 description 3
 230000001133 acceleration Effects 0.000 description 2
 241000208340 Araliaceae Species 0.000 description 1
 241000611421 Elia Species 0.000 description 1
 235000003140 Panax quinquefolius Nutrition 0.000 description 1
 238000006243 chemical reaction Methods 0.000 description 1
 239000011248 coating agent Substances 0.000 description 1
 238000000576 coating method Methods 0.000 description 1
 235000005035 ginseng Nutrition 0.000 description 1
 235000008434 ginseng Nutrition 0.000 description 1
 238000003780 insertion Methods 0.000 description 1
 239000002184 metal Substances 0.000 description 1
 238000000034 method Methods 0.000 description 1
 230000004048 modification Effects 0.000 description 1
 238000006011 modification reaction Methods 0.000 description 1
 238000006467 substitution reaction Methods 0.000 description 1
Classifications

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06F—ELECTRIC DIGITAL DATA PROCESSING
 G06F30/00—Computeraided design [CAD]
 G06F30/20—Design optimisation, verification or simulation

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06F—ELECTRIC DIGITAL DATA PROCESSING
 G06F30/00—Computeraided design [CAD]
 G06F30/10—Geometric CAD
 G06F30/17—Mechanical parametric or variational design
Abstract
The present invention relates to a kind of accurate efficient radome electrical property energy analysis method, comprising：S1, the geometry physical threedimensional model for establishing antenna house, and carry out mesh generation；S2, according to antenna nearfield 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 highfrequency 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 AeroSpace 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 highprecision 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《Antennaradowi Systems electrical property numerical analysis
With optimizing research》Middle elaboration, the main algorithm of numerical analysis is used as using plane wave spectrumsurface 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 highfrequency 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 threedimensional model of antenna house is established in optimization, and carries out mesh generation to it；
S2, according to antenna nearfield 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 threedimensional 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, E_{S}、H_{S}The respectively electric field on the near field distribution surface of antenna and magnetic field,
For normal vector, G_{0}For 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；θ_{0}For incidence angle；T=d_{1}tanθ_{1}+d_{2}tanθ_{2}+…+d_{N} 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, J_{1}Represent the initial current of antenna house inner wall surface；M_{1}Represent the initial magnetic current of antenna house inner wall surface；J_{2}
Represent the initial current of antenna house outer wall surface；M_{2}Represent the initial magnetic current of antenna house outer wall surface；E^{inc}Represent antenna house surface
Incident electric fields；E^{r}Represent the reflected field on antenna house surface；E^{t}Represent the transmitted electric fields on antenna house surface；H^{inc}Represent antenna house
The incident magnetic on surface；H^{r}Represent the reflection magnetic field on antenna house surface；H^{t}Represent 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 (E^{inc},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 speedup 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 steadystate 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 threedimensional model of antenna house is established in optimization, and to the geometry physical threedimensional mould
Type carries out mesh generation；
S2, according to antenna nearfield 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 hardwareaccelerated), 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 threedimensional 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, E_{S}、H_{S}The respectively electric field on the near field distribution surface of antenna and magnetic field,For normal vector, G_{0}For 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；θ_{0}For incidence angle；T=d_{1}tanθ_{1}+d_{2}tanθ_{2}+…+d_{N} 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, J_{1}Represent the initial current of antenna house inner wall surface；M_{1}Represent the initial magnetic current of antenna house inner wall surface；J_{2}
Represent the initial current of antenna house outer wall surface；M_{2}Represent the initial magnetic current of antenna house outer wall surface；E^{inc}Represent antenna house surface
Incident electric fields；E^{r}Represent the reflected field on antenna house surface；E^{t}Represent the transmitted electric fields on antenna house surface；H^{inc}Represent antenna house
The incident magnetic on surface；H^{r}Represent the reflection magnetic field on antenna house surface；H^{t}Represent 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 (E^{inc},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 steadystate 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 highfrequency 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 speedup 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 abovementioned
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 threedimensional model of antenna house is established in optimization, and carries out mesh generation to it；
S2, according to antenna nearfield 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 threedimensional 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, E_{S}、H_{S}The respectively electric field on the near field distribution surface of antenna and magnetic field,For method
To vector, G_{0}For 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；
θ_{0}For incidence angle；T=d_{1}tanθ_{1}+d_{2}tanθ_{2}+…+d_{N}tanθ_{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, J_{1}Represent the initial current of antenna house inner wall surface；M_{1}Represent the initial magnetic current of antenna house inner wall surface；J_{2}Represent
The initial current of antenna house outer wall surface；M_{2}Represent the initial magnetic current of antenna house outer wall surface；E^{inc}Represent entering for antenna house surface
Radio field；E^{r}Represent the reflected field on antenna house surface；E^{t}Represent the transmitted electric fields on antenna house surface；H^{inc}Represent antenna house surface
Incident magnetic；H^{r}Represent the reflection magnetic field on antenna house surface；H^{t}Represent 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 (E^{inc},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)  20170818  20170818  Accurate and efficient antenna housing electrical property analysis method 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201710712133.5A CN107391880B (en)  20170818  20170818  Accurate and efficient antenna housing electrical property analysis method 
Publications (2)
Publication Number  Publication Date 

CN107391880A true CN107391880A (en)  20171124 
CN107391880B CN107391880B (en)  20200731 
Family
ID=60352813
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201710712133.5A Active CN107391880B (en)  20170818  20170818  Accurate and efficient antenna housing electrical property analysis method 
Country Status (1)
Country  Link 

CN (1)  CN107391880B (en) 
Cited By (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN108920831A (en) *  20180629  20181130  西安电子科技大学  The quick calculation method that highspeed aircraft cover high temperature ablation influences antenna electric performance 
CN110380208A (en) *  20190703  20191025  惠州市德赛西威智能交通技术研究院有限公司  A kind of double arc millimetrewave radar antenna houses of Varyingthickness and design method 
CN113625062A (en) *  20210729  20211109  西安电子科技大学  Antenna housing electrical property estimation method based on Taylor expansion method 
Citations (8)
Publication number  Priority date  Publication date  Assignee  Title 

US20100066595A1 (en) *  20080918  20100318  Lee Chul J  Electromagnetic (em) solver using a shooting bouncing ray (sbr) technique 
CN102590656A (en) *  20120103  20120718  西安电子科技大学  Antenna cover electric property forecasting method based on distant field 
CN103870654A (en) *  20140326  20140618  西安电子科技大学  Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing 
CN104750960A (en) *  20131230  20150701  南京理工大学  Method for rapidly extracting electrical property parameter of metal trusstype radome 
CN104794262A (en) *  20150330  20150722  西安电子科技大学  Electromechanicalintegrationbased shaped beam design method for satelliteborne mesh antenna 
CN105388449A (en) *  20151127  20160309  南京航空航天大学  Method of measuring influence of radome on antenna array directionfinding performance 
CN106772301A (en) *  20161202  20170531  上海无线电设备研究所  A kind of multilayer nonparallel interfaces medium electromagnetic scattering simulation algorithm based on medium ray tracing 
CN106934097A (en) *  20170209  20170707  西安电子科技大学  Towards the spatial networks antenna key dynamics mode choosing method of electrical property 

2017
 20170818 CN CN201710712133.5A patent/CN107391880B/en active Active
Patent Citations (8)
Publication number  Priority date  Publication date  Assignee  Title 

US20100066595A1 (en) *  20080918  20100318  Lee Chul J  Electromagnetic (em) solver using a shooting bouncing ray (sbr) technique 
CN102590656A (en) *  20120103  20120718  西安电子科技大学  Antenna cover electric property forecasting method based on distant field 
CN104750960A (en) *  20131230  20150701  南京理工大学  Method for rapidly extracting electrical property parameter of metal trusstype radome 
CN103870654A (en) *  20140326  20140618  西安电子科技大学  Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing 
CN104794262A (en) *  20150330  20150722  西安电子科技大学  Electromechanicalintegrationbased shaped beam design method for satelliteborne mesh antenna 
CN105388449A (en) *  20151127  20160309  南京航空航天大学  Method of measuring influence of radome on antenna array directionfinding performance 
CN106772301A (en) *  20161202  20170531  上海无线电设备研究所  A kind of multilayer nonparallel interfaces medium electromagnetic scattering simulation algorithm based on medium ray tracing 
CN106934097A (en) *  20170209  20170707  西安电子科技大学  Towards the spatial networks antenna key dynamics mode choosing method of electrical property 
NonPatent 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 (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN108920831A (en) *  20180629  20181130  西安电子科技大学  The quick calculation method that highspeed aircraft cover high temperature ablation influences antenna electric performance 
CN110380208A (en) *  20190703  20191025  惠州市德赛西威智能交通技术研究院有限公司  A kind of double arc millimetrewave radar antenna houses of Varyingthickness and design method 
CN110380208B (en) *  20190703  20210205  惠州市德赛西威智能交通技术研究院有限公司  Variablethickness doublearc millimeter wave radome and design method 
CN113625062A (en) *  20210729  20211109  西安电子科技大学  Antenna housing electrical property estimation method based on Taylor expansion method 
CN113625062B (en) *  20210729  20220517  西安电子科技大学  Antenna housing electrical property estimation method based on Taylor expansion method 
Also Published As
Publication number  Publication date 

CN107391880B (en)  20200731 
Similar Documents
Publication  Publication Date  Title 

CN107391880A (en)  A kind of accurate efficient radome electrical property energy analysis method  
Huang et al.  Electromagnetic scattering of randomly rough soil surfaces based on numerical solutions of Maxwell equations in threedimensional simulations using a hybrid UV/PBTG/SMCG method  
CN105468909B (en)  Timelag power system electromechanic oscillation mode computational methods based on SOD PS R algorithms  
CN105302962A (en)  Electromechanical integration optimization design method of reflector antenna on the basis of structureelectromagnetic hybrid unit  
CN112989680B (en)  FVFD farfield integration boundary condition calculation method for reducing grid usage  
Chen et al.  Improved scatteringmatrix method and its application to analysis of electromagnetic wave reflected by reentry plasma sheath  
Wu et al.  The modern high frequency methods for solving electromagnetic scattering problems  
Song et al.  Accurate evaluation of Green's functions for a lossy layered medium by fast extraction of surfaceand leakywave modes  
CN105160115B (en)  Based on the approximate reflector antenna electromechanics Integrated Optimal Design method with sensitivity analysis  
CN104915465A (en)  Metal target transient electromagnetic scattering analysis method based on delay laguerre polynomials  
CN104992035A (en)  Quick calculation method for terahertz band surface rough target electromagnetic scattering  
CN113158527B (en)  Method for calculating frequency domain electromagnetic field based on implicit FVFD  
CN105630740B (en)  Emi analysis method based on matrix Taylor series expansion  
Koh et al.  New scaling factors of 2D isotropicdispersion finite difference time domain (IDFDTD) algorithm for lossy media  
CN107305536A (en)  Mix the discontinuous Jia Lvejin methods of rank time domain  
CN106528954A (en)  Electromechanical integration optimization design method for space mesh antenna under action of temperature load  
Lai et al.  Hybrid method of higherorder MoM and nyström disretization PO for 3D PEC problems  
Yin et al.  Study on calculation and verification of radiowave propagation using parabolic equation for the antenna near the ground  
Oğuzer et al.  E‐polarized beam scattering by an open cylindrical PEC strip having an arbitrary “conical‐section” profile  
Fan et al.  The efficient analysis of electromagnetic scattering from composite structures using hybrid CFIEIEFIE  
Sun et al.  Analysis of scattering from composite conductor and dielectric objects using single integral equation method and MLFMA based on JMCFIE  
Wang et al.  Fast simulations of electromagnetic scattering from 1d rough surface through adaptively modified characteristic basis function method  
CN106649197A (en)  Calculation method of internal scattering characteristic of complicated cavity  
CN108268677A (en)  The electrical performance evaluation method and device of antenna house  
Nie et al.  Fast solutions of wideband RCS of objects using general characteristic basis function method with the model based parameter estimation 
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 