CN106777536A - Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation - Google Patents

Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation Download PDF

Info

Publication number
CN106777536A
CN106777536A CN201611056365.1A CN201611056365A CN106777536A CN 106777536 A CN106777536 A CN 106777536A CN 201611056365 A CN201611056365 A CN 201611056365A CN 106777536 A CN106777536 A CN 106777536A
Authority
CN
China
Prior art keywords
physical quantity
field
electro
dimensional
physical
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.)
Pending
Application number
CN201611056365.1A
Other languages
Chinese (zh)
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.)
Shandong Computer Science Center National Super Computing Center in Jinan
Shandong Computer Science Center
Original Assignee
Shandong Computer Science Center
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 Shandong Computer Science Center filed Critical Shandong Computer Science Center
Priority to CN201611056365.1A priority Critical patent/CN106777536A/en
Publication of CN106777536A publication Critical patent/CN106777536A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Abstract

Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation of the invention, study and select suitable data structure, display logic and visualization technique, after the calculating of fine Electromagnetic Simulation terminates, dimension and represented content according to Electromagnetic Simulation result data, 2-D data is visualized with geometric ways such as curve, broken lines under rectangular co-ordinate, polar coordinates, three-dimensional data is visualized in the way of color under rectangular coordinate system, convenient use person quickly and accurately understands the physical quantity result that fine Electromagnetic Simulation is calculated.By fine electromagnetism numerical simulation post processing, designer can be in the reasonability without processing test design in the case of in kind, it is possible to analyze its electromagnetic compatibility problem.Product design costses can be reduced, shortens the research and development of products cycle.Requirement to research staff's professional skill can be reduced simultaneously so that design simplerization of antenna and microwave system, popular.

Description

Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation
Technical field
The present invention relates to a kind of electro-magnetic far-field two, three-dimensional visual processing method, in particular, more particularly to a kind of base In the electro-magnetic far-field two, three-dimensional visual processing method of fine Electromagnetic Simulation.
Background technology
Electromagnetic field numerical simulation calculate refer to using numerical computation method frequency domain or time domain solve Maxwell equations and its Derivative other equations.With the development of computer technology and Numerical Calculation of Electromagnetic Fields method, electromagnetic field numerical simulation is extensive It is applied to the electromagnetism such as microwave and millimetre-wave attenuator, radar, precise guidance, electromagnetic compatibility, medical diagnosis, navigation and geological prospecting neck Domain, the modernization necessary means as installation electromagnetical specificity analysis with design, with huge practical value.
In recent years, the design of large-scale antenna array layout etc. " complicated, TV university " system proposes more smart to Electromagnetic Simulation calculating Really, more efficient active demand, generally requires the fine electromagnetism of the large-scale parallel of cores up to ten thousand or even hundreds of thousands core parallel scale The business softwares such as simulation calculation, traditional FEKO, HFSS is limited by parallel check figure cannot meet large-scale parallel emulation The demand of calculating.Fine Electromagnetic Simulation can carry out accurate simulation to electromagnetic problems such as antenna and microwave systems.Based on fine electricity Magnetic numerical method is emulated, and can analyze the indexs such as antenna pattern, gain, the S parameter of antenna designed by user.With reference to The modeling of CAD partial parameters can be optimized to antenna performance, can also calculate the RCS (RCS) of target, analysis emulation The design of the microwave devices such as the electromagnetic property of microwave device, assisted microwave synthesis wave filter, power splitter.
, it is necessary to the Simulation result data to producing is researched and analysed to know emulation after fine Electromagnetic Simulation is completed Series of results, the intuitively directional diagram in graphical representation far field, the electromagnetism distribution map near field and CURRENT DISTRIBUTION cloud atlas etc..Finely Electromagnetic Simulation engineering in general is calculated on supercomputer, can produce a large amount of, complicated data, if directly studying these Data, it is not only time-consuming, and easily produce deviation.With the raising of computing power, the scale of Solve problems is continued to increase, multiple Miscellaneous degree is improved constantly all to be increased, it is necessary to process and calculate the data volume for producing, and the analysis and explanation of result of calculation are increasingly stranded Difficulty, has had a strong impact on the efficiency and quality of visualization in scientific computing.
The physical significance in far field refers to the electromagnetic field of infinite point in theory, typically take in practice 10 wavelength or Electromagnetic field beyond 100 wavelength, is mainly used to characterize the radiation characteristic of electromagnetic field.It is three-dimensional in far field electromagnetic simulation calculation Directional diagram, two dimensional plot are related to numerous physical quantitys, and same physical quantity also has the component of different directions, each Component also includes real part, imaginary part as plural number.In addition, the normalized problem of component can be also related to.Therefore, construction display During figure, it is necessary to the accurate related data read in destination file.
The content of the invention
A kind of shortcoming in order to overcome above-mentioned technical problem of the invention, there is provided electro-magnetic far-field based on fine Electromagnetic Simulation 2nd, three-dimensional visual processing method.
Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation of the invention, it is characterised in that logical Following steps are crossed to realize:
A) obtains electro-magnetic far-field direction and value, because the data of electro-magnetic far-field are with spherical coordinates 's What form was given, it is first determinedThe span in directionAnd its value within the range, it is set toIt is then determined that span (the θ in θ directionsmin, θmax) and its value within the range, it is set to θ1、 θ2、...、θn2
B) obtains the value of physical quantity, works as directionTake successivelyDirection θ takes θ successively1、 θ2、...、θn2When, store every physical parameter of target to be analyzed;If what is analyzed is the antenna performance of target, store Every physical parameter beReal part,Imaginary part, E_ θ real parts, E_ θ imaginary parts, E absolute values, EdB values, gain G ain, side To coefficient;If what is analyzed is the RCS of target, the every physical parameter for storing isReal part, Imaginary part, E_ θ real parts, E_ θ imaginary parts, E absolute values, E dB values, RCS;
Wherein:E is electric-field intensity,E_ θ are respectively electric-field intensity edgeDirection, the numerical value in θ directions, EdB values are electricity The field ratio of gains, direction coefficient is the directivity factor of antenna, and RCS is the radar area of target;
C) normalized of physical quantitys, if treating normalized physical quantity set Gi={ Gi1、Gi2、...、Gij、...、 Gi(n1×n2) represent, GiMiddle i represents the species of physical quantity, GijRepresent j-th numerical value in i class physical quantitys, then it is right by formula (1) It carries out radius normalized:
Its kind of Gimax=max { Gi}、Gimin=min { Gi, rijIt is j-th normalized result of numerical value, 1 in i class physical quantitys ≤j≤n1×n2;
D) conversion from spherical coordinates to three-dimensional rectangular coordinate, if normalization after the i-th class physical quantity in j-th numerical value rijInstitute Deflection in corresponding spherical coordinates is respectivelyθl, corresponding three-dimensional rectangular coordinate is (xij,yij,zij), defining solid can The acquiescence maximum radius of regarding displaying is R;In order to realize displaying of the view in three-dimensional rectangular coordinate first quartile, using formula (2) by j-th numerical value r in the i-th class physical quantityijSpherical coordinates be converted into three-dimensional rectangular coordinate:
zij=R*rij cosθl
Wherein, 1≤k≤n1,1≤l≤n2;
E) two-dimensional curves are drawn, according to the section to be drawn physical quantity,A certain angle is selected in direction and θ directions Used as fixed value, and other direction by the institute of the respective physical amount to be drawn in the section a little, is carried out as independent variable The treatment of step c) and step d), you can obtain curve of the physical quantity to be shown in the section;
F) 3 dimensional drawings are drawn,, as independent variable, the selection physical quantity to be shown exists for direction and θ directionsAccording to It is secondary to takeθ is taken successively1、θ2、...、θn2When institute a little, carry out the place of step c) and step d) Reason, you can obtain the curved surface object figure of the physical quantity to be shown.
Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation of the invention, it is describedValue model Enclose forThe span of the θ is (0 °, 180 °), θ1、θ2、...、θn2=0,1,2 ..., 180 °.
The antenna parameters such as antenna pattern, gain, main lobe width, the minor level of aerial radiation problem Main Analysis antenna And S parameter.Need to provide several driving sources of conventional antenna, such as the delta sources of wire antenna, Waveguide slot day for radiation problem The ripple port of line, the feed of microstrip antenna, current source, magnetic current source.Scattering problems are primarily used to analyze the radar scattering of target Section (RCS), is used to judge the Stealth Fighter of the targets such as aircraft naval vessel.The driving source of scattering problems needs to support that plane wave swashs Encourage, any additional source forcing.RCS at distance objective any distance position of orientation can be calculated (as specified the near field on certain line RCS)。
The beneficial effects of the invention are as follows:Shortcoming in order to overcome above-mentioned technical problem of the invention, it is suitable to study and select Data structure, display logic and visualization technique, after the calculating of fine Electromagnetic Simulation terminates, according to Electromagnetic Simulation result data Dimension and represented content, 2-D data is carried out with geometric ways such as curve, broken lines under rectangular co-ordinate, polar coordinates Visualization, three-dimensional data is visualized in the way of color under rectangular coordinate system, and convenient use person is quickly and accurately Solve the physical quantity result that fine Electromagnetic Simulation is calculated.
By fine electromagnetism numerical simulation post processing, designer can check in the case of without processing material object The reasonability of design, it is possible to analyze its electromagnetic compatibility problem.Product design costses can be reduced, shortens the research and development of products cycle.Together When can reduce requirement to research staff's professional skill so that it is design simplerization of antenna and microwave system, popular.
Brief description of the drawings
Fig. 1 is in the present inventionIt is 0 °, the data storage format that θ changes in its span;
Fig. 2 is in the present inventionIt is 1 °, the data storage format that θ changes in its span;
Fig. 3 is in the present inventionIt is 360 °, the data storage format that θ changes in its span.
Specific embodiment
The invention will be further described with embodiment below in conjunction with the accompanying drawings.
Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation of the invention, it is characterised in that logical Following steps are crossed to realize:
A) obtains electro-magnetic far-field direction and value, because the data of electro-magnetic far-field are with spherical coordinates 's What form was given, it is first determinedThe span in directionAnd its value within the range, it is set toIt is then determined that span (the θ in θ directionsmin, θmax) and its value within the range, it is set to θ1、 θ2、...、θn2
Wherein,Span can be The span of θ can be (0 °, 180 °), θ1、θ2、...、θn2=0,1,2 ..., 180 °.
B) obtains the value of physical quantity, works as directionTake successivelyDirection θ takes θ successively1、 θ2、...、θn2When, store every physical parameter of target to be analyzed;If what is analyzed is the antenna performance of target, store Every physical parameter beReal part,Imaginary part, E_ θ real parts, E_ θ imaginary parts, E absolute values, E dB values, gain G ain, side To coefficient;If what is analyzed is the RCS of target, the every physical parameter for storing isReal part, Imaginary part, E_ θ real parts, E_ θ imaginary parts, E absolute values, E dB values, RCS;
Wherein:E is electric-field intensity,E_ θ are respectively electric-field intensity edgeDirection, the numerical value in θ directions, E dB values are The electric field ratio of gains, direction coefficient is the directivity factor of antenna, and RCS is the radar area of target;
As shown in figure 1, givingIt is 0 °, the data storage format that θ changes in its span, its kind 8000.000000 is frequency, and Pha is representedTheta represents θ;When Pha takes 0 °, Theta is by 0 ° to 180 ° change (due to a piece Width limitation only gives the changing value that Theta is by 0 ° to 29 °), the 3rd row to the 7th column data is followed successively byReal part,It is empty Portion, E_ θ real parts, E_ θ imaginary parts, E absolute values, the 8th row and the 9th row are corresponding RCS values, and last is classified as unnecessary row, does not do Treatment.
As shown in Fig. 2 in giving the present inventionIt it is 1 °, the data storage format that θ changes in its span, Fig. 3 gives Go out in the present inventionIt is 360 °, the data storage format that θ changes in its span, it can be seen that,At (0 °, 360 °) During scope, θ change in the range of (0 °, 180 °), the value that it includes is counted as 361 are multiplied by 181.
C) normalized of physical quantitys, if treating normalized physical quantity set Gi={ Gi1、Gi2、...、Gij、...、 Gi(n1×n2) represent, GiMiddle i represents the species of physical quantity, GijRepresent j-th numerical value in i class physical quantitys, then it is right by formula (1) It carries out radius normalized:
Its kind of Gimax=max { Gi}、Gimin=min { Gi, rijIt is j-th normalized result of numerical value, 1 in i class physical quantitys ≤j≤n1×n2;
D) conversion from spherical coordinates to three-dimensional rectangular coordinate, if normalization after the i-th class physical quantity in j-th numerical value rijInstitute Deflection in corresponding spherical coordinates is respectivelyθl, corresponding three-dimensional rectangular coordinate is (xij,yij,zij), defining solid can The acquiescence maximum radius of regarding displaying is R;In order to realize displaying of the view in three-dimensional rectangular coordinate first quartile, using formula (2) by j-th numerical value r in the i-th class physical quantityijSpherical coordinates be converted into three-dimensional rectangular coordinate:
zij=R*rij cosθl
Wherein, 1≤k≤n1,1≤l≤n2;
E) two-dimensional curves are drawn, according to the section to be drawn physical quantity,A certain angle is selected in direction and θ directions Used as fixed value, and other direction by the institute of the respective physical amount to be drawn in the section a little, is carried out as independent variable The treatment of step c) and step d), you can obtain curve of the physical quantity to be shown in the section;
F) 3 dimensional drawings are drawn,, as independent variable, the selection physical quantity to be shown exists for direction and θ directionsAccording to It is secondary to takeθ is taken successively1、θ2、...、θn2When institute a little, carry out the place of step c) and step d) Reason, you can obtain the curved surface object figure of the physical quantity to be shown.

Claims (2)

1. a kind of electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation, it is characterised in that by following Step is realized:
A) obtains electro-magnetic far-field direction and value, because the data of electro-magnetic far-field are with spherical coordinates Form Be given, it is first determinedThe span in directionAnd its value within the range, it is set toIt is then determined that span (the θ in θ directionsmin, θmax) and its value within the range, it is set to θ1、 θ2、...、θn2
B) obtains the value of physical quantity, works as directionTake successivelyDirection θ takes θ successively1、θ2、...、θn2 When, store every physical parameter of target to be analyzed;If what is analyzed is the antenna performance of target, the every physics for storing Parameter isReal part,Imaginary part, E_ θ real parts, E_ θ imaginary parts, E absolute values, E dB values, gain G ain, direction coefficient;Such as What fruit to be analyzed is the RCS of target, then the every physical parameter for storing isReal part,Imaginary part, E_ θ realities Portion, E_ θ imaginary parts, E absolute values, E dB values, RCS;
Wherein:E is electric-field intensity,E_ θ are respectively electric-field intensity edgeDirection, the numerical value in θ directions, E dB values are electric field The ratio of gains, direction coefficient is the directivity factor of antenna, and RCS is the radar area of target;
C) normalized of physical quantitys, if treating normalized physical quantity set Gi={ Gi1、Gi2、...、Gij、...、 Gi(n1×n2)Represent, GiMiddle i represents the species of physical quantity, GijRepresent j-th numerical value in i class physical quantitys, then it is right by formula (1) It carries out radius normalized:
r i j = G i j - G i min G i m a x - G i min - - - ( 1 )
Its kind of Gimax=max { Gi}、Gimin=min { Gi, rijIt is j-th normalized result of numerical value, 1≤j in i class physical quantitys ≤n1×n2;
D) conversion from spherical coordinates to three-dimensional rectangular coordinate, if normalization after the i-th class physical quantity in j-th numerical value rijIt is corresponding Spherical coordinates in deflection be respectivelyθl, corresponding three-dimensional rectangular coordinate is (xij,yij,zij), define 3 d visualization The acquiescence maximum radius of displaying is R;In order to realize displaying of the view in three-dimensional rectangular coordinate first quartile, using formula (2) By j-th numerical value r in the i-th class physical quantityijSpherical coordinates be converted into three-dimensional rectangular coordinate:
Wherein, 1≤k≤n1,1≤l≤n2;
E) two-dimensional curves are drawn, according to the section to be drawn physical quantity,A certain angle conduct is selected in direction and θ directions Fixed value, and other direction by the institute of the respective physical amount to be drawn in the section a little, carries out step as independent variable C) with the treatment of step d), you can obtain curve of the physical quantity to be shown in the section;
F) 3 dimensional drawings are drawn,, as independent variable, the selection physical quantity to be shown exists for direction and θ directionsTake successivelyθ takes θ successively1、θ2、...、θn2When institute a little, carry out the treatment of step c) and step d), you can Obtain the curved surface object figure of the physical quantity to be shown.
2. the electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation according to claim 1, it is special Levy and be:It is describedSpan be (0 °, 360 °),The θ Span be (0 °, 180 °), θ1、θ2、...、θn2=0,1,2 ..., 180 °.
CN201611056365.1A 2016-11-26 2016-11-26 Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation Pending CN106777536A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611056365.1A CN106777536A (en) 2016-11-26 2016-11-26 Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611056365.1A CN106777536A (en) 2016-11-26 2016-11-26 Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation

Publications (1)

Publication Number Publication Date
CN106777536A true CN106777536A (en) 2017-05-31

Family

ID=58910748

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611056365.1A Pending CN106777536A (en) 2016-11-26 2016-11-26 Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation

Country Status (1)

Country Link
CN (1) CN106777536A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109918858A (en) * 2019-04-22 2019-06-21 西北工业大学 A kind of CST phantom antenna data visualization method
CN110502864A (en) * 2019-08-29 2019-11-26 中国航空工业集团公司沈阳飞机设计研究所 A kind of airframe and its electromagnet shield effect appraisal procedure
CN111208171A (en) * 2020-01-12 2020-05-29 中国人民解放军国防科技大学 Method for judging transmittance performance based on electromagnetic periodic structure impedance characteristic
CN112149271A (en) * 2019-06-28 2020-12-29 Aptiv技术有限公司 Method for simulating antenna

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104808204A (en) * 2015-04-13 2015-07-29 电子科技大学 Moving-target detecting method and imaging method of stationary transmitter bistatic forward-looking synthetic aperture radar (SAR)
CN105244636A (en) * 2015-10-13 2016-01-13 北京宇航系统工程研究所 Calculation method for conversion coefficient of electric field measurement antenna in spacecraft cabin

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104808204A (en) * 2015-04-13 2015-07-29 电子科技大学 Moving-target detecting method and imaging method of stationary transmitter bistatic forward-looking synthetic aperture radar (SAR)
CN105244636A (en) * 2015-10-13 2016-01-13 北京宇航系统工程研究所 Calculation method for conversion coefficient of electric field measurement antenna in spacecraft cabin

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孔龙星等: "一种新的电磁场数据三维绘制算法", 《现代电子技术》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109918858A (en) * 2019-04-22 2019-06-21 西北工业大学 A kind of CST phantom antenna data visualization method
CN112149271A (en) * 2019-06-28 2020-12-29 Aptiv技术有限公司 Method for simulating antenna
CN112149271B (en) * 2019-06-28 2024-02-13 Aptiv技术有限公司 Method for simulating antenna
CN110502864A (en) * 2019-08-29 2019-11-26 中国航空工业集团公司沈阳飞机设计研究所 A kind of airframe and its electromagnet shield effect appraisal procedure
CN111208171A (en) * 2020-01-12 2020-05-29 中国人民解放军国防科技大学 Method for judging transmittance performance based on electromagnetic periodic structure impedance characteristic

Similar Documents

Publication Publication Date Title
CN104992001B (en) The accurate quick calculation method of extensive MIMO array Antenna Far Field radiation field
CN106777536A (en) Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation
CN104993251B (en) A kind of large planar array Antenna measuring table cascades optimization method
CN102385655B (en) Simulation method for radiation coupling electromagnetic susceptibility of electronic equipment
CN103246781A (en) Array antenna radar cross section reduction method based on space mapping
CN108957450A (en) A kind of millimetre-wave radar GPU real time three-dimensional imaging method
CN104063544A (en) Near-field effect error analysis method
CN102156764A (en) Multi-resolution precondition method for analyzing aerial radiation and electromagnetic scattering
CN102054094A (en) Fast directional multilevel simulation method for planar microstrip circuit
He et al. A vector meshless parabolic equation method for three-dimensional electromagnetic scatterings
CN109765538A (en) The scattered field of inhomogeneous dielectric object body determines method
CN104485513A (en) Real-time integrating method for wideband antenna array
CN113158485B (en) Electromagnetic scattering simulation method for electrically large-size target under near-field condition
Cakir et al. An FDTD-based parallel virtual tool for RCS calculations of complex targets
CN103065007A (en) Antenna modeling method based on test data of directional diagram and input impedance
CN106250657A (en) A kind of wireless charging emulation mode combining MATLAB and HFSS software
CN110362877A (en) The environment electromagnetics analysis of scattering method of uncertain factor
CN104915465A (en) Metal target transient electromagnetic scattering analysis method based on delay laguerre polynomials
CN104573257A (en) High-speed algorithm for judging electromagnetic shelter based on surface element space diversity
CN106709150A (en) Fine simulation-based current distribution and near-field electromagnetic distribution three-dimensional visual method
CN106503349B (en) One type periodic structure electromagnetic characteristic of scattering quick calculation method
Liu et al. Application of Chaos Mind Evolutionary Algorithm in Antenna Arrays Synthesis.
CN105205299B (en) The quick Dimension Reduction Analysis method of TV university electromagnetic characteristic of scattering
CN107305536A (en) Mix the discontinuous Jia Lvejin methods of rank time domain
CN106446377A (en) Method and system for calculating pilot frequency multi-antenna electromagnetic field distribution

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20170531

RJ01 Rejection of invention patent application after publication