CN108959834A - A kind of time-domain radiation rapid simulation method of the electromagnetic pulse simulator of coaxial feeding - Google Patents
A kind of time-domain radiation rapid simulation method of the electromagnetic pulse simulator of coaxial feeding Download PDFInfo
- Publication number
- CN108959834A CN108959834A CN201811133253.0A CN201811133253A CN108959834A CN 108959834 A CN108959834 A CN 108959834A CN 201811133253 A CN201811133253 A CN 201811133253A CN 108959834 A CN108959834 A CN 108959834A
- Authority
- CN
- China
- Prior art keywords
- space
- pmc
- electric field
- face
- magnetic field
- 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/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
The present invention proposes a kind of time-domain radiation analogy method of the electromagnetic pulse simulator of coaxial feeding, and it is the face PMC that this method, which selects simulator to calculate the plane of symmetry of the overall model in space, and according to the plane of symmetry, establishes a half model of the overall model；A, the medium of the setting face PMC and UPML layers of intersection is perfect magnetic conductor；B, it at non-UPML, when tangential electric field iteration on the face PMC, if using the magnetic field to the space other than the half model established, is substituted with the negative value in the magnetic field for symmetrically calculating space with the face PMC；The normal electric field indirect assignment in space is 0 other than the half model established；C, it when magnetic field iteration at coaxial line excited source, if using the electric field to the space other than the half model established, is substituted with the electric field for symmetrically calculating space with the face PMC；At non-actuated source and at non-UPML, the tangential magnetic field indirect assignment in space is 0 other than the half model established.
Description
Technical field
The present invention relates to a kind of method for numerical simulation of electromagnetic impulse radiation wave simulation device time domain field computation.
Background technique
Due to being capable of providing electromagnetic pulse (EMP) radiation environment, the large-scale vertical polarization Radiating-wave simulator of coaxial feed
(Baum C E is published in report " Sensor and Simulation " nineteen eighty-two Note 277, " Review of Hybrid
And Equivalent-Electric-Dipole EMP Simulators ") in foreign countries it is widely used to large complicated electronics
(Singaraju etc. is published in report " ATHAMAS in the effect experiments such as electromagnetic compatibility, electromagnetic interference and the electromagnetic failure of system
Memos " 1976 years Memo 12, " Use of ATHAMAS II (Larger VPD) for Testing Aircraft in
Flight").It is therefore desirable to carry out sunykatuib analysis to large-scale vertical polarization Radiating-wave simulator.
(Zhu Xiangqin et al. is published in and " calculates physics " network edition page 2 of periodical 2018, and " discrete resistance adds for existing document
The parallel FDTD of the large-scale vertical polarization EMP Radiating-wave simulator of load is simulated ", network address: http://kns.cnki.net/kcms/
Detail/11.2011.O4.20180503.1308.008.html it) mentions, it is limited using the three-dimensional parallel time domain of coaxial feed
Difference (FDTD) method analyzes the time-domain radiation field of such large-scale simulator, and when simulation assumes the metal that ground is laid with
Net is infinity, therefore the metal mesh on ground directly is truncated using uniaxial completely permutation (UPML) absorbing boundary.But it is real
When the engineering construction of border, the metal mesh that ground is laid with is limited size, therefore when board design before engineering construction, using parallel
What UPML absorbing boundary was truncated is the earth around metal mesh when FDTD is calculated.The situation of infinite metal net is laid in ground
It compares, the increase of zoning can be caused, to increase calculation amount and calculate the time.
In view of when the metal mesh of laying finite region, the structure of entire model has symmetry below simulator, and
Tangential electric field is parallel with the plane of symmetry, therefore, calculates time and calculating space to save, can set the plane of symmetry to ideal
The magnetic conductor surface (PMC) is simulated.(Zhu Xiangqin et al. is published in the tenth national High-Power Microwave in 2015 to existing document
Scientific seminar article number HPM-2015-D17, " PMC image method is calculated in vertical polarization parallel-plate Bounded-wave Simulator FDTD
In application ") mention, the PMC image method based on image theory is used for the FDTD of vertical polarization parallel-plate Bounded-wave Simulator
In simulation；Also have document (Lu Feng et al. is published in " light laser and the particle beams " periodical the 3rd phase page 259 of volume 15 in 2003,
" electromagnetic pulse emitter near region Numerical Simulation ") it mentions, the PMC image method based on image theory is used for portable nuclear electromagnetic
The simulation of pulse simulator field.
But the driving source that simulator uses in above-mentioned document is plane source feed, therefore is not related to coaxial feed
When PMC image method processing, have not yet to see a kind of side of the vertical polarization simulator radiation field of Fast simulation coaxial feed
Case.
Summary of the invention
The large-scale vertical polarization radiated wave mould that it is an object of the invention to solve to be laid with limited regional metal net on the ground
Zoning is excessive when quasi- device three-dimensional FDTD simulation, calculates the problem of overlong time, proposes a kind of electromagnetism arteries and veins of coaxial feeding
The time-domain radiation analogy method of pulse simulator.
Solution of the invention mainly comprises the steps that
1) it selects the plane of symmetry of entire model for the face PMC, and according to the plane of symmetry, only establishes the half of this overall model
Model；
2) space electric field, magnetic field interative computation
A, the medium of the setting face PMC and UPML layers of intersection is perfect magnetic conductor, and the electric and magnetic fields for carrying out this region change
Generation；
B, the electric field iteration in other regions of space
At non-UPML, when tangential electric field iteration on the face PMC, if using to the space other than the half model established
Magnetic field, just substituted with the negative value in the magnetic field for symmetrically calculating space with the face PMC；At non-UPML, the half model established
With the normal electric field indirect assignment of the external space for 0；
C, the magnetic field iteration in other regions of space
When magnetic field iteration at coaxial line excited source, if using the electric field to the space other than the half model established,
Just substituted with the electric field for symmetrically calculating space with the face PMC；At non-actuated source and at non-UPML, the half model established with
The tangential magnetic field indirect assignment of the external space is 0；
3) space electric field iteration and space magnetic field iteration terminate to get the time-domain radiation analog result of test point is arrived.
The beneficial effects of the present invention are:
(1) the present invention is based on finite region metal mesh floor file, coaxial feed large-scale vertical polarization Radiating-wave simulators
The symmetry of radiation space field, it is only necessary to establish symmetry model half therein, the large size simulator radiation field can be completed
Numerical simulation, breach existing three-dimensional simulation algorithm and calculate that space is excessive, calculates the limitation of overlong time.
(2) present invention can be used in the situation on infinite metal net laying ground, can be used for symmetrical configuration and tangential electric field
The time domain of other types of simulator (such as EMP Radiating-wave simulator bipyramid cage modle antenna) radiation field of symmetrical coaxial feed
Radiomimesis provides effectively for the radiation field of Fast simulation analysis mode device, for the Numeric simulation design before Practical Project implementation
Technological means.
Detailed description of the invention
Fig. 1 is the method flow diagram of the embodiment of the present invention.
Fig. 2 is the laying work area of metal mesh on ground in present example below simulator, is R by radiusdSemicircle
Arc, radius RwCircular arc and two parallel lines constitute.In addition, R is the radius that simulator top projects to ground in figure, survey
Point B is the position of a measuring point on ground.
Fig. 3 and Fig. 4 is respectively integrally to simulate in the embodiment of the present invention and intend side using the PMC mirror image digital-to-analogue of coaxial feeding
The Comparative result of the electric field component of the position test point B and B ' that method obtains perpendicular to the ground.
Specific embodiment
In order to which the radiation field of the large-scale vertical polarization radiated wave EMP simulator to coaxial feed is simulated, the application couple
The field of coaxial feed part carries out specially treated, realizes the Fast simulation of the vertical polarization simulator radiation field of coaxial feed.Under
Face is in conjunction with attached drawing, and by embodiment, the present invention is further described.
In the present embodiment, the height that simulator is arranged is 32m；Simulator top projects to ground when metal mesh is laid in Fig. 2
The radius R in face is 20m；Constitute the radius R of the circular arc of ground metal netdFor 30m, RwFor 40m；Select measuring point B and origin O apart
30m.In addition selection measuring point B ' is located at the position of 5m right above measuring point B.Calculating total time step is 5000 steps.The embodiment of the present invention
Techniqueflow chart is as shown in Figure 1.It is to establish model and mesh generation, followed by main program reads in model parameter and according to sharp first
The simulation that source carries out electric and magnetic fields is encouraged, including at UPML absorbing boundary iteration.After final time step terminates, it will calculate
To result compared with the result integrally simulated, the two meets preferably, and can be saved greatly using technology of the invention
The calculating space of amount and calculating time.
Specific link is described below:
(1) it selects the plane of symmetry of entire model for the face PMC, and according to the plane of symmetry, establishes comprising large-scale vertical polarization spoke
The half of metal mesh and the earth overall model that ejected wave simulator, below ground are laid with.
(2) model parameter is read in, it is 0 that setting initial time, which walks, is divided according to the process sum for participating in parallel computation parallel
The subdomain of calculating.
(3) the discrete obtained one-dimensional FDTD equation of equation for transmission line met according to TEM wave in coaxial line carries out coaxial line
The iteration of middle voltage and current.
(4) after adjacent subdomain exchange magnetic field data, all subdomains is set and complete entire space electric field iteration together, including
At UPML absorbing boundary, the face PMC and its lateral area.One when the iterative manner of the electric field at UPML absorbing boundary and whole simulation
Sample；At non-UPML, when tangential electric field iteration on the face PMC, if using to the space other than the half model established
Magnetic field is just substituted with the negative value in the magnetic field for symmetrically calculating space with the face PMC；At non-UPML, the half model established with
The normal electric field indirect assignment of the external space is 0；
(5) after adjacent subdomain exchange electric field data, all subdomains complete entire space magnetic field iteration, including coaxial line together
The load of driving source, UPML absorbing boundary, the face PMC and its lateral area when feed.It changes for the magnetic field at coaxial line excited source
Dai Shi, if as soon as the electric field to the space other than the half model established is used, with the electric field for symmetrically calculating space with the face PMC
Substitution；Unit permeance and the unit permeance phase in vacuum herein is arranged in the face PMC and UPML layer intersection, when UPML layers of magnetic field iteration
Together；At non-actuated source and at non-UPML, the tangential magnetic field indirect assignment in space is 0 other than the half model established；It is other
The magnetic field in region is consistent with iterative manner when whole simulation.
(6) when calculating time step is walked less than 5000 total times, step (3)~(5) are repeated, until calculating time step and always
It calculates and terminates when time step is equal.The electric field component of the position test point B and B ' finally simulated perpendicular to the ground is at any time
Variation is as shown in Figures 2 and 3 respectively.As a comparison, giving the result that entire model is integrally simulated in figure.From Fig. 2 and Fig. 3
As can be seen that being met very well using the time domain waveform that the present invention simulates resulting measuring point with what is integrally simulated.And according to simulation
It is found that the relative mistake 0.52% of measuring point B peak, the relative mistake of rising edge is 2.40%, and halfwidth is that relative mistake is 0.68%；It surveys
The relative mistake 0.52% of point B ' peak value, the relative mistake of rising edge are 2.61%, and halfwidth is that relative mistake is 0.35%.To sum up may be used
Know, the present invention simulates obtained result and meets very well with the result integrally simulated.In addition, compared with integrally simulation, the present invention
Zoning reduce about 47.7%, it is entire to calculate the saving of time time about 43.4%.
The PMC mirror image number analogy method of this coaxial feeding of the present embodiment is hung down mainly for the large size of coaxial feed
The time-domain radiation simulation of straight polarized radiation wave electromagnetic pulse simulator (is laid with the limited big and infinitely great feelings of metal mesh including ground
Shape).The present invention can also be for other type simulation device (such as EMP with symmetrical configuration and the symmetrical coaxial feed of tangential electric field
Radiating-wave simulator bipyramid cage modle antenna etc.) radiation field time-domain radiation simulation.
Claims (1)
1. a kind of time-domain radiation rapid simulation method of the electromagnetic pulse simulator of coaxial feeding, which is characterized in that including with
Lower link:
1) selecting simulator to calculate the plane of symmetry of the overall model in space is the face PMC, and according to the plane of symmetry, establishes the entirety
One half model of model；
2) space electric field, magnetic field interative computation
A, the medium of the setting face PMC and UPML layers of intersection is perfect magnetic conductor, carries out the electric and magnetic fields iteration in this region；
B, the electric field iteration in other regions of space
At non-UPML, when tangential electric field iteration on the face PMC, if using the magnetic to the space other than the half model established
, then it is substituted with the negative value in the magnetic field for symmetrically calculating space with the face PMC；
At non-UPML, the normal electric field indirect assignment in space is 0 other than the half model established；
C, the magnetic field iteration in other regions of space
When magnetic field iteration at coaxial line excited source, if using the electric field to the space other than the half model established, use
The electric field for symmetrically calculating space with the face PMC substitutes；
At non-actuated source and at non-UPML, the tangential magnetic field indirect assignment in space is 0 other than the half model established；
3) space electric field iteration and space magnetic field iteration terminate to get the time-domain radiation analog result of test point is arrived.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811133253.0A CN108959834B (en) | 2018-09-27 | 2018-09-27 | Time domain radiation rapid simulation method of electromagnetic pulse simulator with coaxial line feed |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811133253.0A CN108959834B (en) | 2018-09-27 | 2018-09-27 | Time domain radiation rapid simulation method of electromagnetic pulse simulator with coaxial line feed |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108959834A true CN108959834A (en) | 2018-12-07 |
| CN108959834B CN108959834B (en) | 2023-04-07 |
Family
ID=64471105
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811133253.0A Active CN108959834B (en) | 2018-09-27 | 2018-09-27 | Time domain radiation rapid simulation method of electromagnetic pulse simulator with coaxial line feed |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108959834B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005274233A (en) * | 2004-03-23 | 2005-10-06 | Fujitsu Ltd | Electromagnetic field simulator, electromagnetic field analyzer, electromagnetic field simulation program, and electromagnetic field analysis program |
| CN102207987A (en) * | 2011-05-31 | 2011-10-05 | 中国航天标准化研究所 | Method for accelerating three-dimensional finite-difference time-domain electromagnetic field simulation by using graphic processing unit (GPU) based on Open computer language (OpenCL) |
| CN105891563A (en) * | 2014-12-16 | 2016-08-24 | 中国人民解放军63973部队 | High-altitude nuclear explosion electromagnetic pulse standard signal analog device |
-
2018
- 2018-09-27 CN CN201811133253.0A patent/CN108959834B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005274233A (en) * | 2004-03-23 | 2005-10-06 | Fujitsu Ltd | Electromagnetic field simulator, electromagnetic field analyzer, electromagnetic field simulation program, and electromagnetic field analysis program |
| CN102207987A (en) * | 2011-05-31 | 2011-10-05 | 中国航天标准化研究所 | Method for accelerating three-dimensional finite-difference time-domain electromagnetic field simulation by using graphic processing unit (GPU) based on Open computer language (OpenCL) |
| CN105891563A (en) * | 2014-12-16 | 2016-08-24 | 中国人民解放军63973部队 | High-altitude nuclear explosion electromagnetic pulse standard signal analog device |
Non-Patent Citations (1)
| Title |
|---|
| 朱湘琴等: "辐射波电磁脉冲模拟器笼形天线辐射特性的并行计算", 《强激光与粒子束》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108959834B (en) | 2023-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shanker et al. | Analysis of transient electromagnetic scattering phenomena using a two-level plane wave time-domain algorithm | |
| Bagci et al. | A fast stroud-based collocation method for statistically characterizing EMI/EMC phenomena on complex platforms | |
| Aygun et al. | A two-level plane wave time-domain algorithm for fast analysis of EMC/EMI problems | |
| Liu et al. | Simulations of GPR in dispersive media using a frequency-dependent PSTD algorithm | |
| CN104462713A (en) | Electromagnetic compatibility modeling method and system for rail transit vehicle | |
| Sevgi | Electromagnetic modeling and simulation: challenges in validation, verification, and calibration | |
| Godinho et al. | Efficient analysis of sound propagation in sonic crystals using an ACA–MFS approach | |
| Bettini et al. | Fast Fourier transform-volume integral: a smart approach for the electromagnetic design of complex systems in large fusion devices | |
| CN110119586B (en) | Axial conductivity anisotropy transient electromagnetic three-component three-dimensional FDTD forward modeling method | |
| CN109783829B (en) | Electromagnetic field prediction method for three-dimensional FEM (field emission modeling) and two-dimensional FMM (frequency modulation modeling) | |
| CN105844019A (en) | Ground point selection method based on vehicle ground current distribution regularities | |
| Mota et al. | Analysis of tower surge impedance using the finite element method | |
| EP2287757A1 (en) | Multilevel-Multigrid simulation techniques | |
| CN108959834A (en) | A kind of time-domain radiation rapid simulation method of the electromagnetic pulse simulator of coaxial feeding | |
| Garza et al. | High-order Chebyshev-based Nyström Methods for electromagnetics | |
| Arzag et al. | Lightning electric and magnetic fields computation using the 3D-FDTD method and electromagnetic models in presence of different ground configurations | |
| Wu et al. | Modeling of bandpass GPR problem by HIE procedure with enhanced absorption | |
| Yilmaz et al. | A single-boundary implicit and FFT-accelerated time-domain finite element-boundary integral solver | |
| Asadpourahmadchali et al. | Hybrid continuous circuit‐trapezoidal integration method analysis of multi‐cross structure of grounding system | |
| JP3703812B2 (en) | Electromagnetic field analysis method using FDTD method, medium expression method in electromagnetic field analysis, simulation apparatus, and program | |
| Teixeira et al. | Perfectly matched layer and piecewise-linear recursive convolution for the FDTD solution of the 3D dispersive half-space problem | |
| Xu et al. | A new efficient algorithm for the 2D WLP-FDTD method based on domain decomposition technique | |
| Hafner | The multiple multipole program (MMP) and the generalized multipole technique (GMT) | |
| Xie et al. | Simulation and analysis of radiated electromagnetic environment from cable in cabin | |
| CN110569526B (en) | Method for analyzing sound radiation characteristics of power transformer in multiple operating states |
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 |