CN108446471A - Mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method - Google Patents
Mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method Download PDFInfo
- Publication number
- CN108446471A CN108446471A CN201810193564.XA CN201810193564A CN108446471A CN 108446471 A CN108446471 A CN 108446471A CN 201810193564 A CN201810193564 A CN 201810193564A CN 108446471 A CN108446471 A CN 108446471A
- Authority
- CN
- China
- Prior art keywords
- massif
- field
- simulation model
- dimensional
- magnetic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention proposes a kind of mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method, under the premise of being meant to ensure that computational accuracy, computing resource consumption, and efficiently near field electromagnetic field distribution of the emulation mountain area under electromagnet source radiation event are reduced, realizes that step is：Establish mountain model；Obtain Gauss integration sample point coordinate；Build fritter massif three-dimensional simulation model；Calculate the scattering electric field and fringe magnetic field of the scattering electric field and fringe magnetic field and magnetic distribution region to be predicted site at massif D difference simulation model incidence field excitation；Obtain in-field of the massif D difference simulation model in TE waves and when TM waves；Calculate the scattering electric field and fringe magnetic field of massif D difference simulation model region site to be predicted in TE waves；Calculate the scattering electric field and fringe magnetic field of massif D difference simulation model region site to be predicted in TM waves；Mountain three-dimensional model is obtained in site field distribution in region to be predicted and Distribution of Magnetic Field.
Description
Technical field
The invention belongs to Electromagnetic Field and Microwave Technology fields, are related to a kind of mountain area near field electromagnetic field prediction method, specifically relate to
And a kind of mountain area near field electromagnetic field prediction method of the three-dimensional two-dimentional Fast Multiple Method algorithm of moment method mixing, it can be used for for shortwave day
The erection of line and base station provides reasonable proposal.
Background technology
China is the country on mountain more than one, and mountain area area accounts for national territory total area 69.1%.Mountain area is more subject to various
Disaster, disaster often lead to interruption of communication, and the commander of transmission and the rescue and relief work of information can only rely on communication system, grind
Study carefully distribution the preventing and reducing natural disasters with especially important meaning to mountain area of mountain area electromagnetic field.
The main method of existing detection mountain area near field electromagnetic field has a field survey and numerical simulation, field survey refer to
Short-wave antenna or base station predict the magnetic distribution situation in mountain area by way of measurement in the case of radiating, first this
Kind method cannot exclude the interference of survey crew and test probe to electromagnetic field, and secondly, this method needs to spend a large amount of
Manpower, material resources, financial resources cannot be satisfied growing demand.Numerical simulation refers to using Computer Simulation according to concrete model
Simulate mountain area electromagnetic field distribution, such as use full wave method --- three-dimensional moment method (MoM) combine High Performance Computing into
Row Digital calculation modelling, moment method are a kind of by the discrete method for turning to Algebraic Equation set of continuity equation, are a kind of electricity of classics
Field value simulation algorithm.It is by first establishing mountain three-dimensional model, then setting when in the prior art to mountain area electromagnetism field prediction
Simulation parameter is set, dense matrix equation is generated using three-dimensional moment method, it is thick then to solve this using dense matrix equation solver
Then close matrix equation solves electric field and the magnetic field of required point, needs to consume a large amount of time, memory and meter in this process
Resource is calculated, and when mountain three-dimensional model electric size becomes larger, emulation required time, memory and computing resource rapidly increase,
So that also cannot be satisfied growing demand using computer cluster.
Invention content
It is an object of the invention to overcome above-mentioned the shortcomings of the prior art, provide it is a kind of based on three-dimensional moment method and
The mountain area electromagnetic field prediction technique of two-dimentional Fast Multiple Method, it is intended under the premise of ensureing computational accuracy, computing resource consumption is reduced,
And efficiently near field electromagnetic field distribution of the emulation mountain area under electromagnet source radiation event.
The technical thought of the present invention is to obtain specific mountain model by field survey or GIS-Geographic Information System, radiation
Source with a certain range of fritter massif is integrally as new radiation source using centered on radiation source, using three-dimensional modeling and to new
Radiation source carry out simulation calculation using three-dimensional moment method, be more than this new radiation source massif used according to selected face it is two-dimentional
Model and two-dimentional Fast Multiple Method used to solve, in conjunction with key be to provide two-dimentional Fast Multiple Method emulation by three-dimensional moment method
Required in-field finally provides the scattered field that the scattered field on selected face is provided with two-dimentional Fast Multiple Method by three-dimensional moment method
The resultant field for being overlapped to the end.
According to above-mentioned technical thought, the technical solution that the present invention takes includes the following steps：
(1) mountain model is established：
(1a) obtains massif three-dimensional information and dielectric information using field survey or GIS-Geographic Information System, builds massif three
Dimension module, while building the fritter mountain three-dimensional model centered on radiation source for including in massif；
(1b) subtracts fritter mountain three-dimensional model from mountain three-dimensional model, obtains massif three-dimensional differential mode type, and by pre-
Plane where the electromagnetic field to be predicted first set cuts the massif three-dimensional differential mode type, obtains the massif two dimension differential mode for including height
Type；
(2) the Gauss integration sample point coordinate needed for massif D difference simulation model incidence field excitation is obtained：
The frequency of the radiation source of massif D difference model is arranged in (2a), and the radiation source is added to massif D difference model
In, then dielectric constant, conductivity and the magnetic conductivity for the massif D difference model for having added radiation source are set, obtain massif D difference
Simulation model；
(2b) emulates massif D difference simulation model in advance, obtains massif D difference simulation model incidence field excitation institute
The Gauss integration sample point coordinate needed；
(3) fritter massif three-dimensional simulation model is built：
The frequency of the radiation source of fritter mountain three-dimensional model is set, and the radiation source is added to fritter mountain three-dimensional model
In, then dielectric constant, conductivity and the magnetic conductivity for the fritter mountain three-dimensional model for having added radiation source are set, obtain fritter massif
Three-dimensional simulation model, wherein the frequency of the frequency of the radiation source of fritter mountain three-dimensional model and the radiation source of massif D difference model
Rate is identical, dielectric constant, the electricity of dielectric constant, conductivity and the magnetic conductivity and massif D difference model of fritter mountain three-dimensional model
Conductance is identical with magnetic conductivity；
(4) the scattering electric field E at massif D difference simulation model incidence field excitation is calculated using three-dimensional moment methodS-GWith dissipate
Penetrate magnetic field HS-GAnd the scattering electric field E of magnetic distribution region to be predicted siteS-IWith fringe magnetic field HS-I：
(4a) calculates electromagnetic current basic function coefficient vector X：
Using three-dimensional moment method, the impedance matrix A and voltage matrix B of fritter massif three-dimensional simulation model are calculated, and is utilized
Impedance matrix A and voltage matrix B build dense matrix equation AX=B, then solve the dense matrix equation, obtain electromagnetic current base letter
Number system number vector X, whereinInFor electric current basic function coefficient vector, MnFor magnetic current basic function coefficient vector；
(4b) calculates electric current J and magnetic current M on each curved surface of fritter massif three-dimensional simulation model：
According to electromagnetic current basic function coefficient vectorIt calculates on each curved surface of fritter massif three-dimensional simulation model
Electric current J and magnetic current M；
(4c) calculates the scattering electric field E in Gauss integration sample point coordinateS-GWith fringe magnetic field HS-GAnd fritter massif three
Tie up the scattering electric field E of magnetic distribution region to be predicted site in simulation modelS-IWith fringe magnetic field HS-I：
It is calculated high by the electric current J and magnetic current M on each curved surface in fritter massif three-dimensional simulation model using integral equation
Scattering electric field E on this integration sampling point coordinatesS-GWith fringe magnetic field HS-GAnd fritter massif three-dimensional simulation model is to be predicted
The scattering electric field E of magnetic distribution region siteS-IWith fringe magnetic field HS-I；
(5) in-field when obtaining massif D difference simulation model TE waves and in-field when TM waves：
(5a) using the faces xoy as the plane of reference, from scattering electric field ES-GExtract x durection components ExWith y durection components Ey, while from dissipate
Penetrate magnetic field HS-GExtract z durection components Hz, and by Ex、EyAnd HzCombination settings be massif D difference simulation model TE waves when enter
Penetrate field；
(5b) using the faces xoy as the plane of reference, from scattering electric field ES-GExtract z durection components Ez, while from fringe magnetic field HS-GIt carries
Take x durection components HxWith y durection components Hy, and by Ez、HxAnd HyCombination settings be massif D difference simulation model TM waves when
In-field；
(6) using two-dimentional Fast Multiple Method algorithm, massif D difference simulation model is calculated in region site to be predicted
Scattering electric field ES-TEWith fringe magnetic field HS-TE：
(6a) calculates electromagnetic current basis function vector when massif D difference simulation model in-field excitation types are TE waves
XFTE：
In-field when by TE waves that step (5a) obtains is input in two-dimentional Fast Multiple Method algorithm routine as at this time
The in-field of massif D difference simulation model, and using two-dimentional Fast Multiple Method algorithm, calculate massif D difference simulation model
Impedance matrix [AnearTE+AfarTE] and voltage matrix BFTE, recycle impedance matrix [AnearTE+AfarTE] and voltage matrix BFTEStructure
Build matrix equation [AnearTE+AfarTE]XFTE=BFTE, the matrix equation is then solved, electromagnetic current basic function coefficient vector is obtained
XFTE, wherein For electric current basic function coefficient vector,For magnetic current basic function coefficient vector；
(6b) calculates the electric current J on each line segment grid in massif D difference simulation modelFTEWith magnetic current MFTE：
According to electromagnetic current basic function coefficient vectorCalculate each line segment net in massif D difference simulation model
Electric current J on latticeFTEWith magnetic current MFTE；
Scattering electric field E of (6c) the calculating massif D difference simulation model in region site to be predictedS-TEAnd fringe magnetic field
HS-TE：
Using integral equation, pass through the electric current J on each line segment grid in massif D difference simulation modelFTEWith magnetic current MFTE,
Scattering electric field E of the calculating massif D difference simulation model in region site to be predictedS-TEWith fringe magnetic field HS-TE；
(7) using two-dimentional Fast Multiple Method algorithm, massif D difference simulation model is calculated in region site to be predicted
Scattering electric field ES-TMWith fringe magnetic field HS-TM：
(7a) calculates electromagnetic current basis function vector when massif D difference simulation model in-field excitation types are TM waves
XFTM：
In-field when by TM waves that step (5b) obtains is input in two-dimentional Fast Multiple Method algorithm routine as at this time
The in-field of massif D difference simulation model, and using two-dimentional Fast Multiple Method algorithm, calculate massif D difference simulation model
Impedance matrix [AnearTM+AfarTM] and voltage matrix BFTM, recycle impedance matrix [AnearTM+AfarTM] and voltage matrix BFTMStructure
Build matrix equation [AnearTM+AfarTM]XFTM=BFTM, the matrix equation is then solved, electromagnetic current basic function coefficient vector is obtained
XFTM, wherein For electric current basic function coefficient vector,For magnetic current basic function coefficient vector；
(7b) calculates the electric current J on each line segment grid in massif D difference simulation modelFTMWith magnetic current MFTM：
According to electric field basic function coefficient vectorCalculate each line segment grid in massif D difference simulation model
On electric current JFTMWith magnetic current MFTM；
Scattering electric field E of (7c) the calculating massif D difference simulation model in region site to be predictedS-TMAnd fringe magnetic field
HS-TM：
Using integral equation, pass through the electric current J on each line segment grid in massif D difference simulation modelFTMWith magnetic current MFTM,
Scattering electric field E of the calculating massif D difference simulation model in region site to be predictedS-TMWith fringe magnetic field HS-TM；
(8) mountain three-dimensional model is obtained in region site field distribution E to be predictedTOTWith Distribution of Magnetic Field HTOT：
By ES-I、ES-TEAnd ES-TMIt is overlapped, obtains mountain three-dimensional model in site field distribution in region to be predicted
ETOT：ETOT=ES-I+ES-TE+ES-TM, by HS-I、HS-TEAnd HS-TMIt is overlapped, obtains mountain three-dimensional model in regional field to be predicted
Distribution of Magnetic Field H at pointTOT：HTOT=HS-I+HS-TE+HS-TM。
Compared with prior art, the present invention having the following advantages that：
The present invention takes full advantage of three-dimensional moment method and is handled with two-dimentional Fast Multiple Method in the accuracy of processing complex target
The ability of super TV university target or environment；Compared with the prior art, the present invention is under the premise of ensureing computational accuracy, consumption calculations money
Source is less, more efficient.
Description of the drawings
Fig. 1 is the implementation process block diagram of the present invention；
Fig. 2 is the mountain three-dimensional model of the present invention；
Fig. 3 is the fritter mountain three-dimensional model of the present invention；
Fig. 4 is the geometry site of the mountain three-dimensional model and hill block threedimensional model of the present invention；
Fig. 5 is the massif D difference model comprising height of the present invention；
Fig. 6 is the electric field Comparative result of the mountain three-dimensional model global solution and combined method that are acquired using the prior art,
Middle Fig. 6 (a) is EyReal component Comparative result, Fig. 6 (b) is EyImaginary Comparative result；
Fig. 7 is the electric field Comparative result of the mountain three-dimensional model global solution and the present invention that are acquired using the prior art, wherein
Fig. 7 (a) is EzReal component Comparative result, Fig. 7 (b) is EzImaginary Comparative result.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, below in conjunction with attached drawing and specific implementation
Example, the present invention will be described in further detail.It should be appreciated that specific embodiment described herein is only used to explain this hair
It is bright, it is not intended to limit the present invention.
With reference to figure 1, a kind of mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method, including such as
Lower step：
Step 1) establishes mountain model：
Step 1a) using field survey or GIS-Geographic Information System massif three-dimensional information and dielectric information are obtained, use is soft
Part GID builds mountain three-dimensional model, as shown in Fig. 2, its position in reference axis is x-axis：1420m-2020m, y-axis：
1350m-1950m, z-axis：45m-200m, while using include in software GID structure massif being with (1720,1820,150)
The fritter mountain three-dimensional model of the heart, as shown in figure 3, its position in reference axis is x-axis：1570m-1870m, y-axis：1500m-
1800m, z-axis：The relationship of 71m-200m, mountain model and fritter mountain model is as shown in figure 4, fritter mountain model is in massif
In model；
Step 1b) using the Boolean calculation of software GID fritter mountain three-dimensional model, boolean are subtracted from mountain three-dimensional model
Operation is the logical deduction by reasoning method of numerical chracter, including joint, intersects and subtract each other, and is usually used in the place of figure in 3 d modeling software
Reason, obtains massif three-dimensional differential mode type, and selected x=1720m be electromagnetic field to be predicted where plane to cut the massif three-dimensional poor
Model obtains massif three-dimensional differential mode type in the sections yz of x=1720m, as shown in figure 5, being the massif D difference for including height
Model, it is magnetic distribution to be predicted to select and arrive point (1720,1820,270) from point (1720,1820,150) on the faces x=1720m
Region site, it is as shown in Figure 5 in the position of massif D difference model；
Step 2) obtains the Gauss integration sample point coordinate needed for massif D difference simulation model incidence field excitation：
Step 2a) radiation source type of massif D difference model is set for point source, coordinate is (1650,210), frequency
For 10MHz, and the radiation source is added in massif D difference model, then the massif D difference model for having added radiation source is set
Relative dielectric constant be 4, conductivity is 0.001 and relative permeability is 1, obtain massif D difference simulation model；
Step 2b) massif D difference simulation model is emulated in advance, it obtains massif D difference simulation model in-field and swashs
Encourage required Gauss integration sample point coordinate (yk,zk), k indicates k-th of sampled point；
Pre- emulation refers to running an existing simulation model, exports the value of needs, is referred here to the quick multipole of two dimension
Son runs a massif D difference simulation model；
Step 3) builds fritter massif three-dimensional simulation model：
Be arranged the radiation source of fritter mountain three-dimensional model type be half-wave doublet, centre coordinate be (1720,
1650,210), frequency 10MHz, and the radiation source is added in fritter mountain three-dimensional model, then be arranged and added radiation
The relative dielectric constant of the fritter mountain three-dimensional model in source is 4, conductivity is 0.001 and relative permeability is 1, obtains fritter mountain
Body three-dimensional simulation model；
Step 4) calculates the scattering electric field E at massif D difference simulation model incidence field excitation using three-dimensional moment methodS-GWith
Fringe magnetic field HS-GAnd the scattering electric field E of magnetic distribution region to be predicted siteS-IWith fringe magnetic field HS-I：
Step 4a) calculate electromagnetic current basic function coefficient vector X：
Using three-dimensional moment method, the impedance matrix A and voltage matrix B of fritter massif three-dimensional simulation model are calculated, and is utilized
Impedance matrix A and voltage matrix B build dense matrix equation AX=B, then solve the dense matrix equation, obtain electromagnetic current base letter
Number system number vector X；
Step 4b) calculate each curved surface of fritter massif three-dimensional simulation model on electric current J and magnetic current M：
According to electromagnetic current basic function coefficient vector
Use formula
Calculate the electric current J and magnetic current M on each curved surface, wherein fn(r) it is basic function, η0For the wave impedance in free space；
Step 4c) Gauss integration sample point coordinate that step 2) is obtained according to section completion is three-dimensional coordinate, because of step
The rapid Gauss integration sample point coordinate 2) obtained is two-dimensional coordinate, and it is three-dimensional sit that three-dimensional moment method, which calculates required coordinate type,
Mark, therefore need according to the face cut the three-dimensional coordinate that the two-dimensional coordinate completion is three-dimensional moment method needs, it is herein handle
Coordinate (yk,zk) completion be (1720, yk,zk),
Using integral equation
Wherein, J is the electric current on each curved surface, and M is the magnetic current on each curved surface, and j is imaginary unit, and k is medium space wave number,
η0For free space wave impedance, r ' is source point coordinate, and R is distance of the site to source point, and G (R) is free space Green's function；
By the electric current J and magnetic current M on each curved surface in fritter massif three-dimensional simulation model, calculates Gauss integration sampled point and sit
Mark (1720, yk,zk) on scattering electric field ES-GWith fringe magnetic field HS-GAnd fritter massif three-dimensional simulation model point (1720,
1820,150) the scattering electric field E at point (1720,1820,270) is arrivedS-IWith fringe magnetic field HS-I；
In-field when step 5) obtains massif D difference simulation model TE waves and in-field when TM waves：
Step 5a) using the faces xoy as the plane of reference, from scattering electric field ES-GExtract x durection components ExWith y durection components Ey, simultaneously
From fringe magnetic field HS-GExtract z durection components Hz, and by Ex、EyAnd HzCombination settings be massif D difference simulation model TE waves when
In-field；
Step 5b) using the faces xoy as the plane of reference, from scattering electric field ES-GExtract z durection components Ez, while from fringe magnetic field HS-G
Extract x durection components HxWith y durection components Hy, and by Ez、HxAnd HyCombination settings be massif D difference simulation model TM waves when
In-field；
Step 6) is calculated massif D difference simulation model and is arrived point in point (1820,150) using two-dimentional Fast Multiple Method algorithm
(1820,270) the scattering electric field E atS-TEWith fringe magnetic field HS-TE：
Step 6a) calculate electric field basis function vector when massif D difference simulation model in-field excitation types are TE waves
XFTE：
By step 5a) obtained TE waves when in-field be input in two-dimentional Fast Multiple Method algorithm routine as mountain at this time
The in-field of body D difference simulation model, and using two-dimentional Fast Multiple Method algorithm, calculate the resistance of massif D difference simulation model
Anti- matrix [AnearTE+AfarTE] and voltage matrix BFTE, recycle impedance matrix [AnearTE+AfarTE] and voltage matrix BFTEStructure
Matrix equation [AnearTE+AfarTE]XFTE=BFTE, the matrix equation is then solved, electromagnetic current basic function coefficient vector X is obtainedFTE；
Step 6b) calculate the electric current J on each line segment grid in massif D difference simulation modelFTEWith magnetic current MFTE：
According to electromagnetic current basic function coefficient vector
Use formula
Calculate the electric current J on each line segment gridFTEWith magnetic current MFTE,For basic function, η0For the wave in free space
Impedance；
Step 6c) calculate massif D difference simulation model region site to be predicted scattering electric field ES-TEWith scattering magnetic
Field HS-TE：
Using integral equation
Wherein,JFTEFor the electric current on each line segment grid, MFTEFor the magnetic on each line segment grid
Stream, j are imaginary unit, and k is free space wave number, η0For free space wave impedance, ρ ' is source point coordinate, and ρ is site coordinate,For the second class zeroth order Hankel function；
Pass through the electric current J on each line segment grid in massif D difference simulation modelFTEWith magnetic current MFTE, calculate massif D difference
Scattering electric field E of the simulation model at point (1820,150) to point (1820,270)S-TEWith fringe magnetic field HS-TE；
Step 7) is calculated massif D difference simulation model and is arrived point in point (1820,150) using two-dimentional Fast Multiple Method algorithm
(1820,270) the scattering electric field E atS-TMWith fringe magnetic field HS-TM：
Step 7a) calculate electromagnetic current basis function vector when massif D difference simulation model in-field excitation types are TM waves
XFTM：
In-field when by TM waves that step (5b) obtains is input in two-dimentional Fast Multiple Method algorithm routine as at this time
The in-field of massif D difference simulation model, and using two-dimentional Fast Multiple Method algorithm, calculate massif D difference simulation model
Impedance matrix [AnearTM+AfarTM] and voltage matrix BFTM, recycle impedance matrix [AnearTM+AfarTM] and voltage matrix BFTMStructure
Build matrix equation [AnearTM+AfarTM]XFTM=BFTM, the matrix equation is then solved, electromagnetic current basic function coefficient vector is obtained
XFTM；
Step 7b) calculate the electric current J on each line segment grid in massif D difference simulation modelFTMWith magnetic current MFTM：
According to electric field basic function coefficient vector
Use formula
Calculate the electric current J on each line segment gridFTMWith magnetic current MFTM, wherein For basic function,
η0For the wave impedance in free space；
Step 7c) calculate scattering electric field of the massif D difference simulation model at point (1820,150) to point (1820,270)
ES-TMWith fringe magnetic field HS-TM：
Using integral equation
Wherein,JFTMFor the electric current on each line segment grid, MFTMFor the magnetic on each line segment grid
Stream, j are imaginary unit, and k is free space wave number, η0For free space wave impedance, ρ ' is source point coordinate, and ρ is site coordinate,For the second class zeroth order Hankel function；
Pass through the electric current J on each line segment grid in massif D difference simulation modelFTMWith magnetic current MFTM, calculate massif D difference
Scattering electric field E of the simulation model at point (1820,150) to point (1820,270)S-TMWith fringe magnetic field HS-TM；
Step 8) obtains mountain three-dimensional model electric field minute at point (1720,1820,150) to point (1720,1820,270)
Cloth ETOTWith Distribution of Magnetic Field HTOT：
Script ES-TE、ES-TM、HS-TEAnd HS-TMMassif D difference simulation model point (1820,150) to point (1820,
270) scattering electric field and fringe magnetic field at place, because the plane cut is the faces yz of x=1720m, so scattered field is also to correspond to
Mountain three-dimensional model in point (1720,1820,150) to the scattered field at point (1720,1820,270)；By ES-I、ES-TEWith
ES-TMIt is overlapped, obtains mountain three-dimensional model field distribution at point (1720,1820,150) to point (1720,1820,270)
ETOT：ETOT=ES-I+ES-TE+ES-TM, by HS-I、HS-TEAnd HS-TMBe overlapped, obtain mountain three-dimensional model point (1720,
1820,150) Distribution of Magnetic Field H at point (1720,1820,270) is arrivedTOT：HTOT=HS-I+HS-TE+HS-TM。
Below in conjunction with emulation experiment, the technique effect of the present invention is further described：
1. simulated conditions and content
The prior art using computer cluster using three-dimensional moment method calculate massif three-dimensional simulation model point (1720,
1820,150) arrive point (1720,1820,270) electric field and Distribution of Magnetic Field, wherein mountain three-dimensional model simulation parameter setting and
The simulation parameter of fritter mountain three-dimensional model is consistent, as the contrast verification to combined method result；The present invention uses computer
Cluster using three-dimensional moment method calculate fritter massif three-dimensional simulation model point (1720,1820,150) to point (1720,1820,
270) electric field at Gauss integration sampled point needed for electric field and Distribution of Magnetic Field and D difference simulation model incidence field excitation
The magnetic field and；And D difference simulation model is calculated using two-dimentional Fast Multiple Method using computer and arrive point in point (1820,150)
(1820,270) electric field and Distribution of Magnetic Field；Each node configures the central processing unit of two 12 cores of Xeon E5-2692
64GB DDR3 memories, hard disk are several.
2. analysis of simulation result：
Following table is in specific embodiment, is that the prior art and the present invention use memory and CPU core calculation and calculating respectively
The comparison of time：
As can be seen from the above table data, this invention greatly reduces computing resources, and greatly improve emulation
Efficiency；
With reference to figure 6, Fig. 6 (a), Fig. 6 (b) are that the present invention and the prior art arrive point in point (1720,1820,150) respectively
(1720,1820,270) the electric field E onyReal and imaginary parts component；
With reference to figure 7, Fig. 7 (a), Fig. 7 (b) are that the present invention and the prior art arrive point in point (1720,1820,150) respectively
(1720,1820,270) the electric field E onzReal and imaginary parts component；
As can be seen from Figures 6 and 7, the present invention and prior art result coincide good, so this hair of the simulating, verifying
Bright accuracy.
Claims (8)
1. a kind of mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method, it is characterised in that including such as
Lower step：
(1) mountain model is established：
(1a) obtains massif three-dimensional information and dielectric information using field survey or GIS-Geographic Information System, builds massif three-dimensional mould
Type, while building the fritter mountain three-dimensional model centered on radiation source for including in massif；
(1b) subtracts fritter mountain three-dimensional model from mountain three-dimensional model, obtains massif three-dimensional differential mode type, and by setting in advance
Plane where fixed electromagnetic field to be predicted cuts the massif three-dimensional differential mode type, obtains the massif D difference model for including height；
(2) the Gauss integration sample point coordinate needed for massif D difference simulation model incidence field excitation is obtained：
The frequency of the radiation source of massif D difference model is arranged in (2a), and the radiation source is added in massif D difference model,
Dielectric constant, conductivity and the magnetic conductivity of the massif D difference model for having added radiation source are set again, it is imitative to obtain massif D difference
True mode；
(2b) emulates massif D difference simulation model in advance, obtains needed for massif D difference simulation model incidence field excitation
Gauss integration sample point coordinate；
(3) fritter massif three-dimensional simulation model is built：
The frequency of the radiation source of fritter mountain three-dimensional model is set, and the radiation source is added in fritter mountain three-dimensional model,
Dielectric constant, conductivity and the magnetic conductivity of the fritter mountain three-dimensional model for having added radiation source are set again, obtain fritter massif three
Tie up simulation model, wherein the frequency of the frequency of the radiation source of fritter mountain three-dimensional model and the radiation source of massif D difference model
It is identical, dielectric constant, conductivity and the magnetic conductivity of fritter mountain three-dimensional model and the dielectric constant of massif D difference model, conductance
Rate is identical with magnetic conductivity；
(4) the scattering electric field E at massif D difference simulation model incidence field excitation is calculated using three-dimensional moment methodS-GWith scattering magnetic
Field HS-GAnd the scattering electric field E of magnetic distribution region to be predicted siteS-IWith fringe magnetic field HS-I：
(4a) calculates electromagnetic current basic function coefficient vector X：
Using three-dimensional moment method, the impedance matrix A and voltage matrix B of fritter massif three-dimensional simulation model are calculated, and utilize impedance
Matrix A and voltage matrix B build dense matrix equation AX=B, then solve the dense matrix equation, obtain electromagnetic current basic function system
Number vector X, whereinInFor electric current basic function coefficient vector, MnFor magnetic current basic function coefficient vector；
(4b) calculates electric current J and magnetic current M on each curved surface of fritter massif three-dimensional simulation model：
According to electromagnetic current basic function coefficient vectorCalculate the electric current J on each curved surface of fritter massif three-dimensional simulation model
With magnetic current M；
(4c) calculates the scattering electric field E in Gauss integration sample point coordinateS-GWith fringe magnetic field HS-GAnd fritter massif three-dimensional is imitative
The scattering electric field E of magnetic distribution region to be predicted site in true modeS-IWith fringe magnetic field HS-I：
Gauss product is calculated by the electric current J and magnetic current M on each curved surface in fritter massif three-dimensional simulation model using integral equation
Divide the scattering electric field E in sample point coordinateS-GWith fringe magnetic field HS-GAnd fritter massif three-dimensional simulation model is in electromagnetism to be predicted
The scattering electric field E of field distribution region siteS-IWith fringe magnetic field HS-I；
(5) in-field when obtaining massif D difference simulation model TE waves and in-field when TM waves：
(5a) using the faces xoy as the plane of reference, from scattering electric field ES-GExtract x durection components ExWith y durection components Ey, while from scattering magnetic
Field HS-GExtract z durection components Hz, and by Ex、EyAnd HzCombination settings be massif D difference simulation model TE waves when incidence
；
(5b) using the faces xoy as the plane of reference, from scattering electric field ES-GExtract z durection components Ez, while from fringe magnetic field HS-GExtract the side x
To component HxWith y durection components Hy, and by Ez、HxAnd HyCombination settings be massif D difference simulation model TM waves when incidence
；
(6) using two-dimentional Fast Multiple Method algorithm, scattering of the massif D difference simulation model in region site to be predicted is calculated
Electric field ES-TEWith fringe magnetic field HS-TE：
(6a) calculates electromagnetic current basis function vector X when massif D difference simulation model in-field excitation types are TE wavesFTE：
In-field when by TE waves that step (5a) obtains is input in two-dimentional Fast Multiple Method algorithm routine as massif at this time
The in-field of D difference simulation model, and using two-dimentional Fast Multiple Method algorithm, calculate the impedance of massif D difference simulation model
Matrix [AnearTE+AfarTE] and voltage matrix BFTE, recycle impedance matrix [AnearTE+AfarTE] and voltage matrix BFTEBuild square
Battle array equation [AnearTE+AfarTE]·XFTE=BFTE, the matrix equation is then solved, electromagnetic current basic function coefficient vector X is obtainedFTE,
Wherein For electric current basic function coefficient vector,For magnetic current basic function coefficient vector；
(6b) calculates the electric current J on each line segment grid in massif D difference simulation modelFTEWith magnetic current MFTE：
According to electromagnetic current basic function coefficient vectorIt calculates in massif D difference simulation model on each line segment grid
Electric current JFTEWith magnetic current MFTE；
Scattering electric field E of (6c) the calculating massif D difference simulation model in region site to be predictedS-TEWith fringe magnetic field HS-TE：
Using integral equation, pass through the electric current J on each line segment grid in massif D difference simulation modelFTEWith magnetic current MFTE, calculate
Scattering electric field E of the massif D difference simulation model in region site to be predictedS-TEWith fringe magnetic field HS-TE；
(7) using two-dimentional Fast Multiple Method algorithm, scattering of the massif D difference simulation model in region site to be predicted is calculated
Electric field ES-TMWith fringe magnetic field HS-TM：
(7a) calculates electromagnetic current basis function vector X when massif D difference simulation model in-field excitation types are TM wavesFTM：
In-field when by TM waves that step (5b) obtains is input in two-dimentional Fast Multiple Method algorithm routine as massif at this time
The in-field of D difference simulation model, and using two-dimentional Fast Multiple Method algorithm, calculate the impedance of massif D difference simulation model
Matrix [AnearTM+AfarTM] and voltage matrix BFTM, recycle impedance matrix [AnearTM+AfarTM] and voltage matrix BFTMBuild square
Battle array equation [AnearTM+AfarTM]·XFTM=BFTM, the matrix equation is then solved, electromagnetic current basic function coefficient vector X is obtainedFTM,
Wherein For electric current basic function coefficient vector,For magnetic current basic function coefficient vector；
(7b) calculates the electric current J on each line segment grid in massif D difference simulation modelFTMWith magnetic current MFTM：
According to electric field basic function coefficient vectorIt calculates in massif D difference simulation model on each line segment grid
Electric current JFTMWith magnetic current MFTM；
Scattering electric field E of (7c) the calculating massif D difference simulation model in region site to be predictedS-TMWith fringe magnetic field HS-TM：
Using integral equation, pass through the electric current J on each line segment grid in massif D difference simulation modelFTMWith magnetic current MFTM, calculate
Scattering electric field E of the massif D difference simulation model in region site to be predictedS-TMWith fringe magnetic field HS-TM；
(8) mountain three-dimensional model is obtained in region site field distribution E to be predictedTOTWith Distribution of Magnetic Field HTOT：
By ES-I、ES-TEAnd ES-TMIt is overlapped, obtains mountain three-dimensional model in region site field distribution E to be predictedTOT：ETOT
=ES-I+ES-TE+ES-TM, by HS-I、HS-TEAnd HS-TMIt is overlapped, obtains mountain three-dimensional model in site magnetic in region to be predicted
Field distribution HTOT：HTOT=HS-I+HS-TE+HS-TM。
2. the mountain area electromagnetic field prediction technique according to claim 1 based on three-dimensional moment method and two-dimentional Fast Multiple Method,
It is characterized in that, the solution dense matrix equation described in step (4a), using dense matrix solver.
3. the mountain area electromagnetic field prediction technique according to claim 1 based on three-dimensional moment method and two-dimentional Fast Multiple Method,
It is characterized in that, calculating the electric current J and magnetic current M on each curved surface of fritter massif three-dimensional simulation model described in step (4b), calculate
Formula is respectively：
Wherein fn(r) it is basic function, η0For the wave impedance in free space, InFor electric current basic function coefficient vector, MnFor magnetic current base
Function coefficients vector.
4. the mountain area electromagnetic field prediction technique according to claim 1 based on three-dimensional moment method and two-dimentional Fast Multiple Method,
It is characterized in that, the scattering electric field E calculated in Gauss integration sample point coordinate described in step (4c)S-GWith fringe magnetic field HS-G
And the scattering electric field E of magnetic distribution region to be predicted siteS-IWith fringe magnetic field HS-I, calculation formula is respectively：
Wherein, J is the electric current on each curved surface, and M is the magnetic current on each curved surface, and j is imaginary unit, and k is medium space wave number, η0For
Free space wave impedance, r ' are source point coordinate, and R is distance of the site to source point, and G (R) is free space Green's function.
5. the mountain area electromagnetic field prediction technique according to claim 1 based on three-dimensional moment method and two-dimentional Fast Multiple Method,
It is characterized in that, the electric current J calculated in massif D difference simulation model on each line segment grid described in step (6b)FTEAnd magnetic
Flow MFTE, calculation formula is respectively：
Wherein,For basic function, η0For the wave impedance in free space,For electric current basic function coefficient vector,For
Magnetic current basic function coefficient vector.
6. the mountain area electromagnetic field prediction technique according to claim 1 based on three-dimensional moment method and two-dimentional Fast Multiple Method,
It is characterized in that, described in step (6c) calculate massif D difference simulation model region site to be predicted scattering electric field
ES-TEWith fringe magnetic field HS-TE, calculation formula is respectively：
Wherein,JFTEFor the electric current on each line segment grid, MFTEFor the magnetic current on each line segment grid, j is
Imaginary unit, k are free space wave number, η0For free space wave impedance, ρ ' is source point coordinate, and ρ is site coordinate,It is
Two class zeroth order Hankel functions.
7. the mountain area electromagnetic field prediction technique according to claim 1 based on three-dimensional moment method and two-dimentional Fast Multiple Method,
It is characterized in that, the electric current J calculated in massif D difference simulation model on each line segment grid described in step (7b)FTMAnd magnetic
Flow MFTM, calculation formula is respectively：
Wherein,For basic function, η0For the wave impedance in free space,For electric current basic function coefficient vector,For
Magnetic current basic function coefficient vector.
8. the mountain area electromagnetic field prediction technique according to claim 1 based on three-dimensional moment method and two-dimentional Fast Multiple Method,
It is characterized in that, described in step (7c) calculate massif D difference simulation model region site to be predicted scattering electric field
ES-TMWith fringe magnetic field HS-TM, calculation formula is respectively：
Wherein,JFTMFor the electric current on each line segment grid, MFTMFor the magnetic current on each line segment grid, j is
Imaginary unit, k are free space wave number, η0For free space wave impedance, ρ ' is source point coordinate, and ρ is site coordinate,It is
Two class zeroth order Hankel functions.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810193564.XA CN108446471B (en) | 2018-03-09 | 2018-03-09 | Mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810193564.XA CN108446471B (en) | 2018-03-09 | 2018-03-09 | Mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108446471A true CN108446471A (en) | 2018-08-24 |
| CN108446471B CN108446471B (en) | 2019-10-25 |
Family
ID=63193865
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810193564.XA Active CN108446471B (en) | 2018-03-09 | 2018-03-09 | Mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108446471B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109740265A (en) * | 2019-01-07 | 2019-05-10 | 西安电子科技大学 | Electromagnetic environment Tendency Prediction method outside city room based on MoM-UTD |
| CN110929375A (en) * | 2019-10-17 | 2020-03-27 | 中国科学院电子学研究所 | Efficient lens simulation and optimization method based on two-dimensional moment method and ray tracing method |
| CN113297763A (en) * | 2021-05-24 | 2021-08-24 | 北京航空航天大学 | Near-field data rapid lossless compression storage method suitable for moment method |
| CN113989440A (en) * | 2021-11-09 | 2022-01-28 | 北京航空航天大学 | Point fast matching method suitable for moment method |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6714008B1 (en) * | 2002-07-29 | 2004-03-30 | The United States Of America As Represented By The Secretary Of The Navy | Gradiometric measurement methodology for determining magnetic fields of large objects |
| CN1858937A (en) * | 2006-04-30 | 2006-11-08 | 西安电子科技大学 | Automatic lattice dividing method for analyzing multiple antenna and scatterer structure by method of moment |
| CN102156764A (en) * | 2010-12-27 | 2011-08-17 | 南京理工大学 | Multi-resolution precondition method for analyzing aerial radiation and electromagnetic scattering |
| US20140040836A1 (en) * | 2012-07-31 | 2014-02-06 | Taiwan Semiconductor Manufacturing Company Limited | Graded dummy insertion |
| CN105045033A (en) * | 2015-02-10 | 2015-11-11 | 北京理工大学 | Quick three-dimensional mask diffraction near-field calculation method based on sample library and data fitting |
| CN105354421A (en) * | 2015-11-06 | 2016-02-24 | 吉林大学 | Magnetotelluric meshless numerical simulation method for random conductive medium model |
| CN105354368A (en) * | 2015-10-12 | 2016-02-24 | 中国舰船研究设计中心 | Method for quickly predicting directional diagram of microwave antenna under occlusion of near-field complex barrier |
| CN106294284A (en) * | 2016-08-26 | 2017-01-04 | 河海大学 | Rainfall KPT Scatter computing accelerated method based on frequency angle two-dimensional mixing interpolation |
| CN107121624A (en) * | 2017-04-27 | 2017-09-01 | 国网电力科学研究院武汉南瑞有限责任公司 | Partial discharge detection multi-drop arrangement and localization method based on moment method |
-
2018
- 2018-03-09 CN CN201810193564.XA patent/CN108446471B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6714008B1 (en) * | 2002-07-29 | 2004-03-30 | The United States Of America As Represented By The Secretary Of The Navy | Gradiometric measurement methodology for determining magnetic fields of large objects |
| CN1858937A (en) * | 2006-04-30 | 2006-11-08 | 西安电子科技大学 | Automatic lattice dividing method for analyzing multiple antenna and scatterer structure by method of moment |
| CN102156764A (en) * | 2010-12-27 | 2011-08-17 | 南京理工大学 | Multi-resolution precondition method for analyzing aerial radiation and electromagnetic scattering |
| US20140040836A1 (en) * | 2012-07-31 | 2014-02-06 | Taiwan Semiconductor Manufacturing Company Limited | Graded dummy insertion |
| CN105045033A (en) * | 2015-02-10 | 2015-11-11 | 北京理工大学 | Quick three-dimensional mask diffraction near-field calculation method based on sample library and data fitting |
| CN105354368A (en) * | 2015-10-12 | 2016-02-24 | 中国舰船研究设计中心 | Method for quickly predicting directional diagram of microwave antenna under occlusion of near-field complex barrier |
| CN105354421A (en) * | 2015-11-06 | 2016-02-24 | 吉林大学 | Magnetotelluric meshless numerical simulation method for random conductive medium model |
| CN106294284A (en) * | 2016-08-26 | 2017-01-04 | 河海大学 | Rainfall KPT Scatter computing accelerated method based on frequency angle two-dimensional mixing interpolation |
| CN107121624A (en) * | 2017-04-27 | 2017-09-01 | 国网电力科学研究院武汉南瑞有限责任公司 | Partial discharge detection multi-drop arrangement and localization method based on moment method |
Non-Patent Citations (2)
| Title |
|---|
| ANDREW J. HESFORD ETC.: "《On Preconditioning and the Eigensystems of Electromagnetic Radiation Problems》", 《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 * |
| 朱湘琴等: "《垂直极化平行板有界波电磁脉冲模拟器辐射近场的快速估算方法》", 《强激光与粒子束》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109740265A (en) * | 2019-01-07 | 2019-05-10 | 西安电子科技大学 | Electromagnetic environment Tendency Prediction method outside city room based on MoM-UTD |
| CN109740265B (en) * | 2019-01-07 | 2023-01-06 | 西安电子科技大学 | Urban outdoor electromagnetic environment situation prediction method based on MoM-UTD |
| CN110929375A (en) * | 2019-10-17 | 2020-03-27 | 中国科学院电子学研究所 | Efficient lens simulation and optimization method based on two-dimensional moment method and ray tracing method |
| CN113297763A (en) * | 2021-05-24 | 2021-08-24 | 北京航空航天大学 | Near-field data rapid lossless compression storage method suitable for moment method |
| CN113297763B (en) * | 2021-05-24 | 2021-12-10 | 北京航空航天大学 | Near-field data rapid lossless compression storage method suitable for moment method |
| CN113989440A (en) * | 2021-11-09 | 2022-01-28 | 北京航空航天大学 | Point fast matching method suitable for moment method |
| CN113989440B (en) * | 2021-11-09 | 2022-05-13 | 北京航空航天大学 | Point fast matching method suitable for moment method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108446471B (en) | 2019-10-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108446471B (en) | Mountain area electromagnetic field prediction technique based on three-dimensional moment method and two-dimentional Fast Multiple Method | |
| Wang et al. | Multisolver domain decomposition method for modeling EMC effects of multiple antennas on a large air platform | |
| Botha et al. | On the variational formulation of hybrid finite element-boundary integral techniques for electromagnetic analysis | |
| Yang et al. | Novel extension of SBR–PO method for solving electrically large and complex electromagnetic scattering problem in half-space | |
| CN103278848B (en) | The seismic imaging forward modeling method of the fore condition iteration that walks abreast based on MPI | |
| Nakano et al. | Single precision in the dynamical core of a nonhydrostatic global atmospheric model: Evaluation using a baroclinic wave test case | |
| Li et al. | Dosimetry of infant exposure to power‐frequency magnetic fields: Variation of 99th percentile induced electric field value by posture and skin‐to‐skin contact | |
| Perrin et al. | Using a design-of-experiment technique to consider the wire harness load impedances in the FDTD model of an aircraft struck by lightning | |
| Liu et al. | Three-dimensional magnetotellurics modeling using edgebased finite-element unstructured meshes | |
| US10302684B2 (en) | Method for determining electric field and related applications | |
| CN109783829B (en) | Electromagnetic field prediction method for three-dimensional FEM (field emission modeling) and two-dimensional FMM (frequency modulation modeling) | |
| CN113033052B (en) | Electromagnetic rapid numerical modeling method for honeycomb wave-absorbing structure | |
| CN104778286B (en) | Sea skimming device Electromagnetic Scattering Characteristics rapid simulation method | |
| Torres et al. | Effects of the smoothness of partitions of unity on the quality of representation of singular enrichments for GFEM/XFEM stress approximations around brittle cracks | |
| CN105205299B (en) | The quick Dimension Reduction Analysis method of TV university electromagnetic characteristic of scattering | |
| Zhao et al. | Numerical approach for the sensitivity of a high-frequency magnetic induction tomography system based on boundary elements and perturbation method | |
| CN109212498B (en) | Rapid algorithm for radar scattering cross section of airplane formation | |
| Yang et al. | An enriched edge-based smoothed FEM for linear elastic fracture problems | |
| Sabariego et al. | Fast multipole acceleration of the hybrid finite-element/boundary-element analysis of 3-D eddy-current problems | |
| Alotto et al. | A boundary integral formulation on unstructured dual grids for eddy-current analysis in thin shields | |
| CN104091023B (en) | A kind of method for applying specification in electromagnetic field seamed edge elements method | |
| Chen et al. | Parallel DGTD method for transient electromagnetic problems | |
| Wang et al. | Edge diffraction in NURBS-UTD method | |
| JP3703812B2 (en) | Electromagnetic field analysis method using FDTD method, medium expression method in electromagnetic field analysis, simulation apparatus, and program | |
| Gómez-Sousa et al. | Junction modeling for piecewise non-homogeneous geometries involving arbitrary materials |
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 |