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 PDF

Info

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
Application number
CN201710712133.5A
Other languages
Chinese (zh)
Other versions
CN107391880B (en
Inventor
顾昊
樊康
李国平
眭韵
冯红全
余兴
玄晓波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Radio Equipment Research Institute
Original Assignee
Shanghai Radio Equipment Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Radio Equipment Research Institute filed Critical Shanghai Radio Equipment Research Institute
Priority to CN201710712133.5A priority Critical patent/CN107391880B/en
Publication of CN107391880A publication Critical patent/CN107391880A/en
Application granted granted Critical
Publication of CN107391880B publication Critical patent/CN107391880B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical 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 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

A kind of accurate efficient radome electrical property energy analysis method
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>&amp;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>&amp;perp;</mo> </msubsup> <mo>)</mo> <mo>+</mo> <msubsup> <mi>Z</mi> <mn>0</mn> <mo>&amp;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>&amp;perp;</mo> </msubsup> <mo>)</mo> </mrow> </mfrac> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>jksin&amp;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&amp;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>&amp;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>&amp;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>&amp;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>&amp;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.
CN201710712133.5A 2017-08-18 2017-08-18 Accurate and efficient antenna housing electrical property analysis method Active CN107391880B (en)

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 (3)

* Cited by examiner, † Cited by third party
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
CN113625062A (en) * 2021-07-29 2021-11-09 西安电子科技大学 Antenna housing electrical property estimation method based on Taylor expansion method

Citations (8)

* Cited by examiner, † Cited by third party
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

Patent Citations (8)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
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

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
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
CN105468909B (en) Time-lag 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 structure-electromagnetic hybrid unit
CN112989680B (en) FVFD far-field integration boundary condition calculation method for reducing grid usage
Chen et al. Improved scattering-matrix 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 surface-and leaky-wave 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 2-D isotropic-dispersion finite difference time domain (ID-FDTD) 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 higher-order 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 CFIE-IEFIE
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 1-d 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 wide-band 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