Electromagnetic response analysis method and device under broadband excitation signal and medium

Abstract

The invention discloses an electromagnetic response analysis method and device under a broadband excitation signal and a medium, and relates to the technical field of electromagnetism. Performing Fourier decomposition on the broadband time domain excitation signal to be analyzed to obtain a frequency spectrum signal of the broadband time domain excitation signal; calculating a time domain simulation duration and a frequency domain simulation duration of the frequency spectrum signal under each frequency, and calculating a critical frequency according to the time domain simulation duration and the frequency domain simulation duration under each frequency; acquiring a corresponding high-frequency signal component greater than or equal to the critical frequency from the frequency spectrum signal, and calculating the time domain electromagnetic response of the high-frequency signal component; acquiring a corresponding low-frequency signal component smaller than the critical frequency from the frequency spectrum signal, and calculating frequency domain electromagnetic response under each frequency in the low-frequency signal component; and synthesizing the time domain electromagnetic response and the frequency domain electromagnetic response to obtain the electromagnetic response of the broadband time domain excitation signal. According to the method, the calculation efficiency of the broadband electromagnetic response can be improved.

Description

Technical Field

[0001] The present invention relates to the field of electromagnetic technology, and in particular to an electromagnetic response analysis method, device and medium under a broadband excitation signal. Background Art

[0002] At present, the high-altitude electromagnetic pulse (HEMP) caused by the explosion of nuclear weapons above the atmosphere poses a serious security threat to basic energy facilities such as power plant electrical equipment due to its high field strength, wide spectrum range and long range. Modeling and analysis of the electromagnetic response of power plant electrical equipment under the action of HEMP is the key basis for predicting the interference / damage effect of HEMP on equipment, optimizing electromagnetic compatibility design, and taking HEMP protection measures.

[0003] In the prior art, the most intuitive means of solving the HEMP response of key electrical equipment in power plants is full-wave electromagnetic simulation, that is, accurately solving Maxwell's equations under HEMP excitation, calculating the precise field strength distribution inside the equipment, and further obtaining macroscopic response indicators such as induced voltage and current through post-processing. Although there are already abundant algorithm results and mature software in the field of electromagnetic field numerical calculation to complete this task, due to the unique transient waveform characteristics of HEMP and other similar ultra-wideband pulse excitations, the existing full-wave electromagnetic simulation algorithm cannot directly calculate the broadband electromagnetic response in practical applications, resulting in low computational efficiency in calculating broadband electromagnetic response in the prior art. Summary of the invention

[0004] Based on this, it is necessary to provide an electromagnetic response analysis method, device and medium under a broadband excitation signal to address the above technical problems. This method can improve the calculation efficiency of broadband electromagnetic response.

[0005] The present invention adopts the following technical solutions:

[0006] The present invention provides an electromagnetic response analysis method under a broadband excitation signal, comprising:

[0007] Perform Fourier decomposition on the broadband time-domain excitation signal to be analyzed to obtain a spectrum signal of the broadband time-domain excitation signal;

[0008] Calculate the time domain simulation time and the frequency domain simulation time of the spectrum signal at each frequency, and calculate the critical frequency according to the time domain simulation time and the frequency domain simulation time at each frequency;

[0009] Obtaining high-frequency signal components corresponding to frequencies greater than or equal to the critical frequency from the spectrum signal, and calculating the time-domain electromagnetic response of the high-frequency signal components;

[0010] Obtain low-frequency signal components corresponding to frequencies less than the critical frequency from the spectrum signal, and calculate the frequency domain electromagnetic response at each frequency in the low-frequency signal components;

[0011] The time domain electromagnetic response and the frequency domain electromagnetic response are synthesized to obtain the electromagnetic response of the broadband time domain excitation signal.

[0012] Preferably, calculating the time domain simulation duration of the spectrum signal at each frequency includes:

[0013] Determine the maximum frequency of the spectrum signal according to the number of truncated terms of Fourier decomposition and the time length of the broadband time domain excitation signal;

[0014] The required number of time iteration steps is calculated according to the maximum time step length of the broadband time domain excitation signal that satisfies the time domain calculation stability condition and the time length of the broadband time domain excitation signal;

[0015] Obtaining the simulation time required for performing a preset number of iterations on the time length of the broadband time domain excitation signal by using the time domain discontinuous Galerkin algorithm;

[0016] Determine the time domain simulation duration at each frequency based on the maximum frequency, the required number of time iterations, and the simulation time.

[0017] Preferably, the time domain simulation duration at each frequency is calculated as follows:

[0018]

[0019] N = [T / Δt] + 1;

[0020] Among them, t i is the time domain simulation duration of the ith frequency, f max is the maximum frequency of the spectrum signal, f i represents the i-th frequency, Δt is the maximum time step for the broadband time-domain excitation signal to satisfy the stability condition of the time-domain calculation, T is the time length of the broadband time-domain excitation signal, N is the required number of time iteration steps, n is the number of truncation terms, and t0 is the simulation time.

[0021] Preferably, calculating the frequency domain simulation duration of the spectrum signal at each frequency includes:

[0022] According to the maximum frequency of the spectrum signal, determine the middle frequency of the spectrum signal;

[0023] The electromagnetic simulation of the spectrum signal at the intermediate frequency is performed using the boundary element-finite element solver to obtain the frequency domain calculation duration;

[0024] The frequency domain calculation time is determined as the frequency domain simulation time of the spectrum signal at each frequency.

[0025] Preferably, the critical frequency is calculated according to the time domain simulation duration and the frequency domain simulation duration at each frequency, including:

[0026] Construct a frequency optimization function with the goal of minimizing the frequency optimization function value. According to the time domain simulation time and frequency domain simulation time at each frequency, solve the frequency optimization function to obtain the critical frequency; the frequency optimization function is:

[0027]

[0028] Where n is the number of truncated terms in Fourier decomposition, t f0 is the frequency domain simulation duration, f i is the ith frequency, f max is the maximum frequency of the spectrum signal, t i is the time domain simulation duration at the i-th frequency.

[0029] Preferably, calculating the time domain electromagnetic response of the high-frequency signal component includes:

[0030] The time-domain electromagnetic response of the high-frequency signal components is calculated using a time-domain discontinuous Galerkin solver.

[0031] Preferably, calculating the frequency domain electromagnetic response at each frequency in the low-frequency signal component includes:

[0032] By sweeping the frequency, the boundary element-finite element solver is called to calculate the frequency domain electromagnetic response at each frequency in turn.

[0033] Preferably, synthesizing the time domain electromagnetic response and the frequency domain electromagnetic response to obtain the electromagnetic response of the broadband time domain excitation signal includes:

[0034] Performing inverse Fourier transform on the frequency domain electromagnetic response to obtain a transformed frequency domain electromagnetic response;

[0035] The transformed frequency domain electromagnetic response and time domain electromagnetic response are superimposed to obtain the electromagnetic response of the broadband time domain excitation signal.

[0036] The present invention provides an electromagnetic response analysis device under a broadband excitation signal, comprising:

[0037] A decomposition module, used for performing Fourier decomposition on the broadband time-domain excitation signal to be analyzed to obtain a spectrum signal of the broadband time-domain excitation signal;

[0038] The first calculation module is used to calculate the time domain simulation duration and the frequency domain simulation duration of the spectrum signal at each frequency, and calculate the critical frequency according to the time domain simulation duration and the frequency domain simulation duration at each frequency;

[0039] A second calculation module is used to obtain a high-frequency signal component greater than or equal to the critical frequency from the spectrum signal, and calculate the time-domain electromagnetic response of the high-frequency signal component;

[0040] A third calculation module is used to obtain low-frequency signal components corresponding to frequencies less than the critical frequency from the spectrum signal, and calculate the frequency domain electromagnetic response of the low-frequency signal components at each frequency;

[0041] The synthesis module is used to synthesize the time domain electromagnetic response and the frequency domain electromagnetic response to obtain the electromagnetic response of the broadband time domain excitation signal.

[0042] The present invention provides a computer-readable storage medium, wherein the storage medium stores a computer program, and when the computer program is executed by a processor, the electromagnetic response analysis method under the broadband excitation signal is implemented.

[0043] The present invention provides a computer device, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the electromagnetic response analysis method under the broadband excitation signal when executing the program.

[0044] At least one of the above technical solutions adopted by the present invention can achieve the following beneficial effects:

[0045] In the present invention, the critical frequency is calculated by the time domain simulation time and the frequency domain simulation time of the broadband time domain excitation signal at each frequency, and then the low-frequency signal component and the high-frequency signal component of the broadband time domain excitation signal are found by the critical frequency. In this way, the low-frequency signal component of the broadband time domain excitation signal is simulated in the frequency domain, and the high-frequency signal component of the broadband time domain excitation signal is sampled and simulated in the frequency domain, and then synthesized to obtain a complete electromagnetic response. This method combines the advantages of time domain and frequency domain algorithms, so that the broadband time domain excitation signal achieves the optimal computing efficiency, thereby improving the computing efficiency of the electromagnetic response of the broadband time domain excitation signal. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] The drawings described herein are used to provide a further understanding of the present invention and constitute a part of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0047] Figure 1 A schematic flow chart of an electromagnetic response analysis method under a broadband excitation signal provided by the present invention;

[0048] Figure 2 A schematic flow chart of another electromagnetic response analysis method under a broadband excitation signal provided by the present invention;

[0049] Figure 3A schematic diagram of an electromagnetic response analysis device under a broadband excitation signal provided by the present invention;

[0050] Figure 4 A schematic diagram of a computer device for realizing an electromagnetic response analysis method under a broadband excitation signal provided by the present invention. DETAILED DESCRIPTION

[0051] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the specific embodiments of the present invention and the corresponding drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0052] The spectrum range of typical ultra-wideband excitations such as HEMP spans 10kHz-100MHz. This wide spectrum makes the mainstream frequency domain simulation method computationally unbearable (a huge amount of frequency sweep operations are required to cover the HEMP spectrum range); in addition, for application scenarios involving nonlinear components (such as transformer cores), the assumption that the excitation source and the physical field distribution at any position in space in the frequency domain calculation conform to the sinusoidal law in the time domain is no longer valid, and a special equivalent modeling method is required. Therefore, the time domain method is a more widely applicable choice.

[0053] The time domain method, which is theoretically applicable to any broadband excitation, also has efficiency bottlenecks in practice. In order to simulate infinitely large open area problems on computers with limited memory, time domain electromagnetic simulation needs to truncate the computational domain with the help of virtual boundary surfaces or boundary layers, and use absorption boundary conditions, perfect matching layers and other technologies on the boundaries to ensure that the electromagnetic field in the computational domain can be absorbed by the boundaries without reflection. To ensure the accuracy of the calculation, it is usually required that the truncation boundary be no less than half a wavelength away from the research object. For low-frequency signals of several hundred kHz, the half-wavelength exceeds 300 meters, and the length of the computational domain in each coordinate axis direction is no less than 600 meters. The huge amount of meshing caused is unbearable in actual calculations. For low-frequency signals, frequency domain calculations can use boundary element discretization close to the research object as the open domain boundary condition, but this method requires the inversion of high-dimensional matrices in time domain calculations and cannot be directly applied.

[0054] In response to the aforementioned bottleneck, the present invention proposes an electromagnetic response analysis method under broadband excitation signals for ultra-wideband electromagnetic excitations such as HEMP, using a frequency domain simulation algorithm for the low-frequency components of the excitation signal and a time domain simulation algorithm for the high-frequency components, and then synthesizing them to obtain a complete broadband electromagnetic response. This method can combine the advantages of frequency domain and time domain algorithms to achieve optimal computational efficiency for ultra-wideband excitations.

[0055] The device for executing the solution of the present invention may be a server arranged on a service platform, or a device such as a desktop computer, a notebook computer, etc. that can execute the solution of the present invention.

[0056] The technical solutions provided by various embodiments of the present invention are described in detail below in conjunction with the accompanying drawings.

[0057] Figure 1 The present invention is a schematic flow chart of an electromagnetic response analysis method under a broadband excitation signal, which specifically includes the following steps:

[0058] S101, performing Fourier decomposition on the broadband time-domain excitation signal to be analyzed to obtain a frequency spectrum signal of the broadband time-domain excitation signal.

[0059] The transient waveform expression of broadband time-domain excitation is G(t), which can be expanded into Fourier series form as shown in formula (1), i.e., the spectrum signal.

[0060] G(t)≈g(t)=a0+a1 cos(ωt)+b1 sin(ωt)+Ka n cos(nωt)+b n sin(nωt) (1)

[0061] Among them, the basic angular frequency ω=1 / 2πT, T is the time length of the broadband time domain excitation signal, and n is the number of truncation terms. The selection of the number of truncation terms n can satisfy:

[0062]

[0063] Among them, the second norm of any broadband time-domain excitation signal waveform x(t) can be calculated by the following formula:

[0064]

[0065] S102, calculating the time domain simulation duration and the frequency domain simulation duration of the spectrum signal at each frequency, and calculating the critical frequency according to the time domain simulation duration and the frequency domain simulation duration at each frequency.

[0066] Optionally, the time domain simulation duration of the spectrum signal at each frequency is calculated, including: determining the maximum frequency of the spectrum signal according to the number of truncated terms of Fourier decomposition and the time length of the broadband time domain excitation signal; calculating the required number of time iteration steps according to the maximum time step that satisfies the time domain calculation stability condition of the broadband time domain excitation signal and the time length of the broadband time domain excitation signal; obtaining the simulation time required for iterating the time length of the broadband time domain excitation signal a preset number of times through the time domain discontinuous Galerkin algorithm; determining the time domain simulation duration at each frequency according to the maximum frequency, the required number of time iteration steps and the simulation time.

[0067] Specifically, the characteristic scale of the problem at a certain frequency f is defined as

[0068]

[0069] Where Δl1 is the minimum operating wavelength, Δl2 is the minimum penetration depth of the system when the research object is conductive, and v C is the speed of light, μ and σ are the magnetic permeability and electrical conductivity respectively.

[0070] Assume that the frequency of the spectrum signal is the maximum frequency, that is, take According to the research object (equipment) distance from the boundary surface (launch point) at least d = v c / 2f max The computational domain is selected according to the principle of , and the geometric partitioning of the model is completed, requiring the grid feature size not to exceed Δl / 10. Calculate the maximum time step Δt and the time length T of the broadband time domain excitation signal that meet the stability conditions of the time domain calculation, and calculate the required number of time iteration steps N = [T / Δt] + 1. Call the Discontinuous Galerkin Time-Domain (DGTD) algorithm for 10 time iterations, with the simulation time being t0, and then for the frequency in the spectrum signal not less than f i The time estimation of electromagnetic response of all frequency components (1≤i≤n) obtained by time domain simulation is:

[0071]

[0072] Therefore, the time domain simulation duration of the spectrum signal at each frequency can be calculated according to formula (5), t i is the time domain simulation duration of the ith frequency.

[0073] It should be noted that the Discontinuous Galerkin Time-Domain (DGTD, DGTD) algorithm is an improvement on the classical time-domain finite element method. Its essential advantage is that it does not need to invert the high-dimensional matrix in each time iteration like the time-domain finite element method, allowing purely explicit time-domain iteration with high efficiency and low memory consumption. Therefore, the present invention selects the DGTD algorithm as the solver for time-domain electromagnetic simulation.

[0074] Optionally, the frequency domain simulation time of the spectrum signal at each frequency is calculated, including: determining the intermediate frequency of the spectrum signal according to the maximum frequency of the spectrum signal; performing electromagnetic simulation on the spectrum signal at the intermediate frequency using a boundary element-finite element solver to obtain the frequency domain calculation time; and determining the frequency domain calculation time as the frequency domain simulation time of the spectrum signal at each frequency.

[0075] Specifically, take f=f max / 2=n / 4πT, the frequency domain finite element and boundary element methods are used to calculate the electromagnetic simulation at this frequency, and the calculation time is t f0 , can be used as an approximate estimate of the frequency domain simulation calculation time of each frequency component, that is, t f0 As the frequency domain simulation duration of the spectrum signal at each frequency.

[0076] Based on the time domain simulation time and frequency domain simulation time at each frequency, the critical frequency is calculated: First, a frequency optimization function is constructed to minimize the value of the frequency optimization function. According to the time domain simulation time and frequency domain simulation time at each frequency, the frequency optimization function is solved to obtain the critical frequency; the frequency optimization function is:

[0077]

[0078] Where n is the number of truncated terms in Fourier decomposition, t f0 is the frequency domain simulation duration, f i is the ith frequency, f max is the maximum frequency of the spectrum signal, t i is the time domain simulation duration at the i-th frequency.

[0079] The independent variable of the above optimization problem is the critical frequency number i. Since the number is a positive integer between 0 and n, formula (6) defines an integer optimization problem. In practice, the typical value of n can be within 1000. It is easy to obtain the number i and the corresponding optimal critical frequency f by exhaustive enumeration. i .

[0080] The critical frequency f in the present invention i The definition is: The frequency in formula (1) is not less than f i The components are calculated in the time domain, and the frequency is lower than f i The electromagnetic response of the component is obtained by multiple frequency domain simulations. The optimal critical frequency makes the total calculation time the shortest, that is, it satisfies the optimization problem.

[0081] S103, obtaining a high-frequency signal component corresponding to a frequency greater than or equal to a critical frequency from the frequency spectrum signal, and calculating a time-domain electromagnetic response of the high-frequency signal component.

[0082] According to the distance between the research object and the boundary surface, at least d = v c / 2f i , the computational domain is selected and the mesh is divided according to the principle that the characteristic size of the grid unit does not exceed Δl / 10, and the critical frequency is f i For example, the high-frequency signal component is defined as:

[0083] R(t)=G(t)-[a0+a1cos(ωt)+b1sin(ωt)+Kai-1 cos((i-1)ωt)+b i-1 sin((i-1)ωt)](7)

[0084] The time domain discontinuous Galerkin solver can be called to calculate the time domain electromagnetic response of the high-frequency signal component, where the high-frequency signal component is the high-frequency part of the broadband time domain excitation signal.

[0085] S104, obtaining low-frequency signal components corresponding to frequencies less than the critical frequency from the frequency spectrum signal, and calculating frequency domain electromagnetic responses at each frequency in the low-frequency signal components.

[0086] By sweeping the frequency, the boundary element-finite element solver is called to calculate the frequency domain electromagnetic response at each frequency in turn. For example, taking the critical frequency as f i For example, for frequencies below f i The low-frequency signal component is calculated by sweeping the frequency domain electromagnetic response of each frequency point. The boundary element-finite element solver is called to calculate the frequencies f0, f1, L, f i-1 The frequency domain electromagnetic response under (the DC component is regarded as a special case of f0=0) is calculated, and each group of calculation results is stored. Among them, the low-frequency signal component is the low-frequency part of the broadband time domain excitation signal.

[0087] S105, synthesizing the time domain electromagnetic response and the frequency domain electromagnetic response to obtain an electromagnetic response of a broadband time domain excitation signal.

[0088] In one embodiment, the time domain electromagnetic response and the frequency domain electromagnetic response are synthesized to obtain the electromagnetic response of the broadband time domain excitation signal, including: performing an inverse Fourier transform on the frequency domain electromagnetic response to obtain a transformed frequency domain electromagnetic response; and superimposing the transformed frequency domain electromagnetic response and the time domain electromagnetic response to obtain the electromagnetic response of the broadband time domain excitation signal.

[0089] Specifically, the frequency domain electromagnetic response of each frequency point obtained in S104 is extracted through inverse Fourier transform to obtain its corresponding time domain waveform, and then synthesized with the time domain electromagnetic response (high frequency part) obtained in S103 based on the principle of linear superposition to obtain a complete electromagnetic response under a broadband excitation signal.

[0090] For example, when studying the electric field strength at a certain spatial point, assuming that the high-frequency time domain response calculated in S103 is H(t), the frequency responses obtained in S104 are l0(f), l1(f), …, l i-1 (f), the actual electric field strength time domain waveform (electromagnetic response of broadband time domain excitation signal) is:

[0091]

[0092] Among them, the operator represents the inverse Fourier transform, and l k The inverse transformation of (f) is L k (t)(0≤k≤i-1).

[0093] In an exemplary embodiment, the present invention also provides an electromagnetic response analysis method under a broadband excitation signal, such as Figure 2 As shown, this embodiment includes the following steps:

[0094] S201, performing Fourier decomposition on the broadband time-domain excitation signal to obtain a frequency spectrum signal of the broadband time-domain excitation signal.

[0095] S202, estimating the computational complexity of the lower frequency domain and time domain simulation of the broadband time domain excitation signal by means of typical frequency point trial calculation and the like.

[0096] S203, calculating the optimal critical frequency by solving the integer optimization problem defined by formula (6).

[0097] S204, for the part of the spectrum that is not lower than the optimal critical frequency, call the time domain discontinuous Galerkin solver once to obtain the corresponding time domain electromagnetic response; for the part of the spectrum that is lower than the optimal critical frequency, call the boundary element-finite element solver by frequency sweeping to calculate and store the frequency domain electromagnetic response corresponding to each frequency point in turn.

[0098] S205, the frequency domain electromagnetic response of each frequency point in the low frequency part is obtained by inverse Fourier transform to obtain the corresponding time domain waveform, which is then synthesized with the time domain electromagnetic response in the high frequency part based on the principle of linear superposition to obtain a complete electromagnetic response under broadband excitation such as HEMP.

[0099] In the present invention, the present invention is oriented to ultra-wideband electromagnetic excitations such as HEMP, performs frequency domain simulation on the low-frequency components of the broadband excitation, uses time domain simulation algorithms on the high-frequency components, and then synthesizes to obtain a complete broadband electromagnetic response. This method fully combines the advantages of frequency domain and time domain algorithms, achieves optimal computational efficiency for ultra-wideband excitation, and has the following advantages:

[0100] First, the most efficient solution method is customized for each band of the broadband excitation spectrum: for low-frequency components, a finite element-boundary element hybrid solver is selected (the boundary element discrete domain close to the research object is used as the open domain boundary condition of the finite element algorithm), avoiding the bottleneck of setting up a huge calculation area for applying open domain boundary conditions in low-frequency simulation; for high-frequency components, a time-domain discontinuous Galerkin (DGTD) solver is selected, which eliminates the huge burden of frequency domain simulation sweep calculations while also overcoming the problem of the need to invert large-scale matrices in each time iteration of the classical time-domain finite element method.

[0101] Second, the optimal critical frequency for frequency-domain and time-domain simulation calculations is determined through a standard process of trial calculations at typical frequency points and then solving the optimization model, avoiding the burden and efficiency bottleneck of relying on manual experience to select this parameter.

[0102] When applying the electromagnetic response analysis method under the broadband excitation signal provided by the present invention, it is not necessary to Figure 1 The steps are executed in the order shown. The specific execution order of the steps can be determined according to needs, and the present invention does not limit this.

[0103] The above is an electromagnetic response analysis method under a broadband excitation signal provided by one or more embodiments of the present invention. Based on the same idea, the present invention also provides a corresponding electromagnetic response analysis device under a broadband excitation signal, such as Figure 3 shown.

[0104] Figure 3 A schematic diagram of an electromagnetic response analysis device under a broadband excitation signal provided by the present invention, the device 300 comprises:

[0105] A decomposition module 301 is used to perform Fourier decomposition on the broadband time-domain excitation signal to be analyzed to obtain a spectrum signal of the broadband time-domain excitation signal;

[0106] The first calculation module 302 is used to calculate the time domain simulation duration and the frequency domain simulation duration of the spectrum signal at each frequency, and calculate the critical frequency according to the time domain simulation duration and the frequency domain simulation duration at each frequency;

[0107] The second calculation module 303 is used to obtain a high-frequency signal component greater than or equal to the critical frequency from the spectrum signal, and calculate the time-domain electromagnetic response of the high-frequency signal component;

[0108] The third calculation module 304 is used to obtain the low-frequency signal component corresponding to the critical frequency from the spectrum signal, and calculate the frequency domain electromagnetic response at each frequency in the low-frequency signal component;

[0109] The synthesis module 305 is used to synthesize the time domain electromagnetic response and the frequency domain electromagnetic response to obtain the electromagnetic response of the broadband time domain excitation signal.

[0110] For the specific definition of the electromagnetic response analysis device under the broadband excitation signal, please refer to the definition of the electromagnetic response analysis method under the broadband excitation signal above, which will not be repeated here. Each module in the above-mentioned electromagnetic response analysis device under the broadband excitation signal can be implemented in whole or in part by software, hardware and a combination thereof. The above-mentioned modules can be embedded in or independent of the processor in the computer device in the form of hardware, or can be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.

[0111] The present invention also provides a computer-readable storage medium, which stores a computer program, which can be used to execute the above Figure 1 The electromagnetic response analysis method under broadband excitation signal is provided.

[0112] The present invention also provides Figure 4 The structural diagram of the computer device shown in FIG. Figure 4 As shown in the figure, at the hardware level, the computer device includes a processor, an internal bus, a network interface, a memory, and a non-volatile memory, and may also include other hardware required for the business. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs it to achieve the above Figure 1 The electromagnetic response analysis method under broadband excitation signal is provided.

[0113] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided by the present invention can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory or optical memory, etc. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0114] The technical features of the above embodiments may be arbitrarily combined. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of the present invention.

Claims

1. A method for analyzing electromagnetic response under a broadband excitation signal, characterized in that: include: Performing Fourier decomposition on the broadband time-domain excitation signal to be analyzed to obtain a frequency spectrum signal of the broadband time-domain excitation signal; Calculating the time domain simulation duration and the frequency domain simulation duration of the spectrum signal at each frequency, and calculating the critical frequency according to the time domain simulation duration and the frequency domain simulation duration at each frequency; Acquire a high-frequency signal component greater than or equal to the critical frequency from the frequency spectrum signal, and calculate a time-domain electromagnetic response of the high-frequency signal component; Acquire low-frequency signal components corresponding to frequencies less than the critical frequency from the frequency spectrum signal, and calculate frequency-domain electromagnetic responses at each frequency in the low-frequency signal components; The time domain electromagnetic response and the frequency domain electromagnetic response are synthesized to obtain the electromagnetic response of the broadband time domain excitation signal.

2. The method according to claim 1, characterized in that Calculating the time domain simulation duration of the spectrum signal at each frequency includes: Determining the maximum frequency of the spectrum signal according to the number of truncated terms of the Fourier decomposition and the time length of the broadband time-domain excitation signal; Calculate the required number of time iteration steps according to the maximum time step length of the broadband time domain excitation signal that satisfies the time domain calculation stability condition and the time length of the broadband time domain excitation signal; Acquire the simulation time required for performing a preset number of iterations on the time length of the broadband time-domain excitation signal by using a time-domain discontinuous Galerkin algorithm; The time domain simulation duration at each frequency is determined according to the maximum frequency, the required number of time iteration steps and the simulation time.

3. The method according to claim 2, characterized in that The calculation method of the time domain simulation duration at each frequency is: N = [T / Δt] + 1; Among them, t i is the time domain simulation duration of the ith frequency, f max is the maximum frequency of the spectrum signal, f i represents the i-th frequency, Δt is the maximum time step for the broadband time-domain excitation signal to satisfy the stability condition of the time-domain calculation, T is the time length of the broadband time-domain excitation signal, N is the required number of time iteration steps, n is the number of truncation terms, and t0 is the simulation time.

4. The method according to claim 1, characterized in that Calculating the frequency domain simulation time of the spectrum signal at each frequency, including: Determining the middle frequency of the spectrum signal according to the maximum frequency of the spectrum signal; Using a boundary element-finite element solver to perform electromagnetic simulation on the spectrum signal at the intermediate frequency to obtain a frequency domain calculation duration; The frequency domain calculation duration is determined as the frequency domain simulation duration of the spectrum signal at each of the frequencies.

5. The method according to claim 1, characterized in that The calculating of the critical frequency according to the time domain simulation duration and the frequency domain simulation duration at each of the frequencies comprises: A frequency optimization function is constructed with the goal of minimizing the frequency optimization function value. The frequency optimization function is solved according to the time domain simulation duration and the frequency domain simulation duration at each frequency to obtain the critical frequency. The frequency optimization function is: Where n is the number of truncated terms in Fourier decomposition, t f0 is the frequency domain simulation duration, f i is the ith frequency, f max is the maximum frequency of the spectrum signal, t i is the time domain simulation duration at the i-th frequency.

6. The method according to claim 1, characterized in that The calculating the time domain electromagnetic response of the high frequency signal component comprises: The time-domain electromagnetic response of the high-frequency signal component is calculated by a time-domain discontinuous Galerkin solver.

7. The method according to claim 1, characterized in that The calculating the frequency domain electromagnetic response at each frequency in the low frequency signal component comprises: By means of frequency sweeping, the boundary element-finite element solver is called to calculate the frequency domain electromagnetic response at each of the frequencies in turn.

8. The method according to claim 1, characterized in that The synthesizing the time domain electromagnetic response and the frequency domain electromagnetic response to obtain the electromagnetic response of the broadband time domain excitation signal includes: Performing an inverse Fourier transform on the frequency domain electromagnetic response to obtain a transformed frequency domain electromagnetic response; The transformed frequency domain electromagnetic response and the time domain electromagnetic response are superimposed to obtain the electromagnetic response of the broadband time domain excitation signal.

9. An electromagnetic response analysis device under a broadband excitation signal, characterized in that: include: A decomposition module, used for performing Fourier decomposition on the broadband time-domain excitation signal to be analyzed to obtain a spectrum signal of the broadband time-domain excitation signal; A first calculation module is used to calculate the time domain simulation duration and the frequency domain simulation duration of the spectrum signal at each frequency, and calculate the critical frequency according to the time domain simulation duration and the frequency domain simulation duration at each frequency; A second calculation module, used for obtaining a high-frequency signal component greater than or equal to the critical frequency from the spectrum signal, and calculating a time-domain electromagnetic response of the high-frequency signal component; A third calculation module is used to obtain a low-frequency signal component corresponding to a frequency less than the critical frequency from the spectrum signal, and calculate the frequency domain electromagnetic response of each frequency in the low-frequency signal component; The synthesis module is used to synthesize the time domain electromagnetic response and the frequency domain electromagnetic response to obtain the electromagnetic response of the broadband time domain excitation signal.

10. A computer-readable storage medium, characterized in that: The storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 8 is implemented.

Images