Large-scale simulation method and device for time domain field circuit

Abstract

The embodiment of the invention belongs to the field of electromagnetic field analysis, and particularly relates to a large-scale simulation method and device for a time domain field circuit, and the method comprises the steps: sending a magnetic field data obtaining request of a target circuit model to all first processes, and enabling the first processes to be responsible for different sub-regions of the target circuit model; determining an equivalent current of the target circuit model according to the magnetic field data returned by all the first processes; the equivalent current is sent to a second process, and the second process is used for determining a port voltage corresponding to a target circuit model according to the equivalent current and sending the port voltage to all the first processes, so that the first processes update magnetic field data of responsible sub-regions according to the port voltage; and iterating the steps until iteration conditions are met. According to the embodiment of the invention, the problem of low electromagnetic damage simulation precision of a large-size equipment electronic platform is solved.

Description

TECHNICAL FIELD

[0001] Embodiments of the present specification relate to the field of electromagnetic field analysis, and in particular, to a large-scale simulation method and device for time-domain field circuit. BACKGROUND

[0002] The time-domain field circuit simulation method can be used for electromagnetic damage research of electronic equipment. For example, by analyzing the voltage and / or current data of the circuit under nuclear electromagnetic pulse irradiation with time-domain field circuit simulation method, the damage trend of electronic equipment in extreme electromagnetic environment can be effectively predicted, so it is widely used in the military field. However, the existing time-domain field circuit simulation method is limited by the grid size, and the simulation accuracy of large circuits is not high, which is difficult to meet the needs of electromagnetic damage modeling and evaluation of large-size equipment electronic platforms. SUMMARY

[0003] Embodiments of the present specification provide a large-scale simulation method and device for time-domain field circuit, which is used to solve or at least partially solve the problem of low simulation accuracy of electromagnetic damage of large-size equipment electronic platform.

[0004] In order to solve the above technical problems, the first aspect of the embodiments of the present specification provides a large-scale simulation method for time-domain field circuit, the method comprising:

[0005] sending a magnetic field data acquisition request of a target circuit model to all first processes, each of the first processes being responsible for a different sub-region of the target circuit model;

[0006] determining an equivalent current of the target circuit model according to the magnetic field data returned by all first processes;

[0007] sending the equivalent current to a second process, the second process being configured to determine a port voltage corresponding to the target circuit model according to the equivalent current, and sending the port voltage to all first processes, so that the first processes update the magnetic field data of the responsible sub-region according to the port voltage;

[0008] iterating the above steps until the iteration condition is met.

[0009] Further, updating the magnetic field data of the responsible sub-region according to the port voltage comprises:

[0010] updating the electric field data of the responsible sub-region according to the port voltage and the topological information of the target circuit model in the responsible sub-region;

[0011] updating the magnetic field data of the responsible sub-region according to the updated electric field data of the responsible sub-region.

[0012] Further, updating the magnetic field data of the responsible sub-region according to the port voltage further comprises:

[0013] acquire new electromagnetic data corresponding to the target circuit model at the current time;

[0014] update the electric field data of the responsible sub-region according to the electric field data in the new electromagnetic data, the port voltage, and the topology information of the target circuit model in the responsible sub-region;

[0015] update the magnetic field data of the responsible sub-region according to the magnetic field data in the new electromagnetic data and the updated electric field data of the responsible sub-region.

[0016] Further, the large-scale time-domain field circuit simulation method further includes a plurality of third processes, each of which is responsible for a different irradiation sub-region outside the sub-region corresponding to the first process;

[0017] acquire new electromagnetic data corresponding to the target circuit model at the current time, including:

[0018] acquire electromagnetic information of the preset electromagnetic signal;

[0019] send the electromagnetic information to all first processes and all third processes, so that the first processes determine the new electromagnetic data corresponding to the target circuit model at the current time according to the electromagnetic information and the electromagnetic information returned by the third processes.

[0020] Further, the determination process of the electromagnetic information returned by the third process includes:

[0021] determine the passing information of the electromagnetic signal in the responsible sub-region of the process according to the received electromagnetic information;

[0022] determine the electromagnetic information returned to the next process according to the received electromagnetic information, the passing information, and the medium information of the responsible sub-region of the process.

[0023] Further, the responsible sub-region of the first process includes a plurality of secondary sub-regions, and the method further includes:

[0024] determine whether each secondary sub-region corresponds to the topology of the target circuit model;

[0025] If the determination is yes, update the electric field data of the secondary sub-region according to the port voltage and the topology information of the target circuit model in the secondary sub-region;

[0026] If the determination is no, acquire the magnetic field data of the adjacent region of the secondary sub-region, and update the electric field data of the secondary sub-region according to the magnetic field data of the secondary sub-region and the magnetic field data of the adjacent region.

[0027] Further, the responsible region of the first process includes a plurality of secondary sub-regions, and the method further includes:

[0028] sequentially determine whether each secondary sub-region corresponds to a topology of the target circuit model;

[0029] If the determination is yes, electric field data of the secondary sub-region is updated according to electric field data in the new electromagnetic data, the port voltage, and topology information of the target circuit model in the secondary sub-region;

[0030] If the determination is no, magnetic field data of a region adjacent to the secondary sub-region is obtained, and electric field data of the secondary sub-region is updated according to magnetic field data of the secondary sub-region, magnetic field data of the adjacent region, and electric field data in the new electromagnetic data.

[0031] Further, the iteration condition comprises:

[0032] The iteration of the above steps is performed for a preset number of times.

[0033] A second aspect of the embodiments of the present specification provides a large-scale simulation device of time-domain field circuit, the device comprising:

[0034] A first sending module is configured to send a magnetic field data acquisition request of a target circuit model to all first processes, each of the first processes being responsible for a different sub-region of the target circuit model;

[0035] A determination module is configured to determine an equivalent current of the target circuit model according to magnetic field data returned by all the first processes;

[0036] A second sending module is configured to send the equivalent current to a second process, the second process being configured to determine a port voltage corresponding to the target circuit model according to the equivalent current, and send the port voltage to all the first processes, so that the first processes update magnetic field data of the responsible sub-regions according to the port voltage;

[0037] An iteration module is configured to iterate the above steps until an iteration condition is met.

[0038] A third aspect of the embodiments of the present specification provides a computer device comprising a memory, a processor, and a computer program stored in the memory, the computer program being executed by the processor to execute instructions of the large-scale simulation method of time-domain field circuit of any of the preceding embodiments.

[0039] A fourth aspect of the embodiments of the present specification provides a computer storage medium having a computer program stored thereon, the computer program being executed by a processor of a computer device to execute instructions of the large-scale simulation method of time-domain field circuit of any of the preceding embodiments.

[0040] A fifth aspect of the embodiments of the present specification provides a computer program product, the computer program product comprising a computer program, when the computer program is executed by a processor of a computer device, instructions of a large-scale simulation method of time-domain field circuit are executed.

[0041] The large-scale simulation method and device of time-domain field circuit provided by the embodiments of the present specification decouple the electromagnetic field operation and the circuit operation during simulation and distribute them to different processes for processing, thereby accelerating the processing efficiency. Meanwhile, since the electromagnetic field operation is no longer coupled with the circuit operation, a basis is provided for further encryption of the electromagnetic field region, and high-precision simulation of the electromagnetic field region can be realized by using the multi-process parallel computing mode. Finally, through the interactive iteration of the process of the electromagnetic field region and the circuit process, the large-scale simulation of the time-domain field circuit is realized, and the problem of low simulation precision of electromagnetic damage simulation of large-size equipment electronic platforms is solved.

[0042] In order to make the above and other purposes, features and advantages of the embodiments of the present specification more obvious and easy to understand, the following preferred embodiments are described in detail below, and the accompanying drawings are described as follows. BRIEF DESCRIPTION OF DRAWINGS

[0043] In order to more clearly illustrate the technical solutions in the embodiments of the present specification or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or prior art description. Obviously, the drawings in the following description are only some embodiments of the present specification, and for those skilled in the art, other drawings can also be obtained without creative labor on the basis of these drawings.

[0044] Figure 1 A first flowchart of the large-scale simulation method of time-domain field circuit of the embodiments of the present specification is shown;

[0045] Figure 2 A first flowchart of the embodiments of the present specification for updating the magnetic field data of the responsible sub-region is shown;

[0046] Figure 3 A second flowchart of the embodiments of the present specification for updating the magnetic field data of the responsible sub-region is shown;

[0047] Figure 4 A flowchart of the embodiments of the present specification for obtaining new electromagnetic data is shown;

[0048] Figure 5 A flowchart of the determination process of the third process returning electromagnetic information of the embodiments of the present specification is shown;

[0049] Figure 6 A second flowchart of the large-scale simulation method of time-domain field circuit of the embodiments of the present specification is shown;

[0050] Figure 7 A third flow chart of a large-scale simulation method of a time-domain field path according to an embodiment of the present specification is shown;

[0051] Figure 8 A schematic diagram of a command vehicle model according to an embodiment of the present specification is shown;

[0052] Figure 9 A schematic diagram of an antenna radio frequency front-end circuit model according to an embodiment of the present specification is shown;

[0053] FIG. 10(a) shows a schematic diagram of the instantaneous electric field intensity distribution when the simulation time is 10 nanoseconds according to an embodiment of the present specification;

[0054] FIG. 10(b) shows a schematic diagram of the instantaneous electric field intensity distribution when the simulation time is 30 nanoseconds according to an embodiment of the present specification;

[0055] FIG. 11(a) shows a schematic diagram of the voltage change over time according to an embodiment of the present specification;

[0056] FIG. 11(b) shows a schematic diagram of the current change over time according to an embodiment of the present specification;

[0057] Figure 12 A structural diagram of a large-scale simulation device of a time-domain field path according to an embodiment of the present specification is shown;

[0058] Figure 13 A structural diagram of a computer device according to an embodiment of the present specification is shown.

[0059] Explanation of the drawing symbols:

[0060] 1210, first sending module;

[0061] 1220, determination module;

[0062] 1230, second sending module;

[0063] 1240, iteration module;

[0064] 1302, computer device;

[0065] 1304, processor;

[0066] 1306, memory;

[0067] 1308, driving mechanism;

[0068] 1310, input / output module;

[0069] 1312, input device;

[0070] 1314, output device;

[0071] 1316, presentation device;

[0072] 1318, graphical user interface;

[0073] 1320, network interface;

[0074] 1322, communication link;

[0075] 1324, communication bus. DETAILED DESCRIPTION

[0076] The technical solutions in the embodiments of the present specification will be described clearly and completely below in combination with the drawings in the embodiments of the present specification. Obviously, the described embodiments are only part of the embodiments of the present specification, rather than all the embodiments. Based on the embodiments in the present specification, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the embodiments of the present specification.

[0077] It should be noted that the terms "first", "second", and the like in the present specification and claims and the above-described drawings are used to distinguish similar objects, and do not necessarily have to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged under appropriate circumstances, so that the embodiments of the present specification described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion, for example, a process, method, device, product or equipment including a series of steps or units does not have to be limited to those steps or units clearly listed, but can include other steps or units not clearly listed or inherent to these processes, methods, products or equipment.

[0078] The present specification provides method operation steps as described in the embodiments or flowcharts, but can include more or fewer operation steps based on routine or non-creative labor. The order of steps listed in the embodiments is only one of the many step execution orders, and does not represent the only execution order. When the system or device product is actually executed, it can be executed in sequence or in parallel according to the method order shown in the embodiments or drawings.

[0079] It should be noted that the acquisition, storage, use, processing, etc. of data in the technical solutions of the embodiments of the present specification comply with the relevant provisions of national laws and regulations.

[0080] It should be noted that in the embodiments of the present specification, some industry existing solutions of software, components, models, etc. may be mentioned, which should be considered as exemplary, and the purpose is only to illustrate the feasibility of the implementation of the technical solutions of the embodiments of the present specification, but does not mean that the applicant has or will necessarily use the solution.

[0081] In order to meet the needs of electromagnetic damage modeling and evaluation of dense electronic platforms, the electromagnetic environment modeling of the board-level circuit target of the vulnerable equipment of the dense electronic platform needs to be carried out. When facing the above target, the circuit elements such as resistor, capacitor, inductor, diode and triode need to be considered. According to the characteristics of the target electromagnetic pulse frequency band, the main energy electromagnetic wave corresponds to a wavelength much larger than the geometric size of the elements such as resistor, capacitor, inductor, diode and triode. At this time, the elements can be described by lumped circuit, so as to realize the electromagnetic simulation of the circuit elements.

[0082] At present, the time-domain field-circuit simulation of the target lumped circuit is realized by using the finite difference time domain (FDTD) method or the extended method based on the FDTD method, so as to obtain the voltage and / or current data of the circuit changing with time sequence according to the time-domain field-circuit simulation, and predict the damage trend of the electronic equipment in the extreme electromagnetic environment. However, the time-domain field-circuit simulation involves key steps such as electromagnetic field solving, circuit solving and electromagnetic field-circuit data exchange, and the calculation logic is complex. Therefore, the encryption of the FDTD grid is limited, and only the electromagnetic field-circuit cooperative simulation of the order of million to ten million grid size can be realized, which is difficult to meet the simulation precision of the electromagnetic damage of the large-size equipment electronic platform.

[0083] In order to solve the above problems, the time-domain field-circuit simulation of the target circuit is implemented based on the JASMIN framework independently developed by the embodiments of the present specification. In the application program based on the JASMIN framework, the numerical simulation process is completed by three levels of calling, including Main (program main body), grid layer (flow control) and grid sheet (specific implementation). Among them, Main is the program main body, responsible for starting the JASMIN framework, reading data files, creating grid layer and grid sheet classes, executing calculations, terminating the framework and the like; the grid layer is responsible for initialization, calculation time step, execution of calculation in one time step, data update and the like; the grid sheet is the specific implementation of each flow step on the grid layer, which is realized through the function interface reserved by the framework. Through the component, the connection between the grid layer and the grid sheet level is realized, and the same operation can be performed on all grid sheets by using the component on the grid layer level.

[0084] In an embodiment of the present specification, a large-scale simulation method of time-domain field-circuit is provided, which simulates a preset electromagnetic signal irradiating a target circuit model to solve the problem of low simulation precision of electromagnetic damage of large-size equipment electronic platform.

[0085] Specifically, as shown in Figure 1 The large-scale simulation method of time-domain field-circuit includes:

[0086] Step 110, sending a magnetic field data acquisition request of the target circuit model to all first processes, and each first process is responsible for a different sub-region of the target circuit model;

[0087] Step 120, determining the equivalent current of the target circuit model according to the magnetic field data returned by all first processes;

[0088] Step 130, sending the equivalent current to a second process, the second process being used to determine the port voltage corresponding to the target circuit model according to the equivalent current, and sending the port voltage to all first processes, so that the first processes update the magnetic field data of the responsible sub-regions according to the port voltage;

[0089] Step 140, iterating the above steps until the iteration condition is met.

[0090] The embodiment decouples the electromagnetic field operation in simulation from the circuit operation and distributes them to different processes for processing, thereby accelerating the processing efficiency; at the same time, since the electromagnetic field operation is no longer coupled with the circuit operation, a basis is provided for further encryption of the electromagnetic field region, and thus high-precision simulation of the electromagnetic field region can be realized by using the multi-process parallel computing mode; finally, through the interactive iteration of the processes of the electromagnetic field region and the circuit processes, large-scale simulation of the time-domain field circuit is realized, and the problem of low simulation precision of electromagnetic damage simulation of large-size equipment electronic platforms is solved.

[0091] In the embodiments of the present specification, single-core multi-process can be used to accelerate the simulation efficiency of the time-domain field circuit, or more computing resources can be invested to accelerate the simulation efficiency through multi-core multi-process; when multi-core multi-process is used, the upper limit of the grid encryption scale can be improved, and the grid can be encrypted to the order of hundreds of millions, thereby ensuring the time-domain field circuit simulation effect of large-size circuits, and thus solving the problem of low simulation precision of electromagnetic damage simulation of large-size equipment electronic platforms.

[0092] In the embodiments of the present specification, a plurality of first processes are responsible for different sub-regions of the target circuit model, and when a preset electromagnetic signal irradiates the target circuit model, each first process can determine the electromagnetic data (i.e., magnetic field data and electric field data) of the sub-region responsible by itself, so as to disperse the data processing pressure of a single process. When all first processes have electromagnetic data of the sub-regions responsible by themselves, the main process sends a magnetic field data acquisition request of the target circuit model to all first processes, all first processes return the magnetic field data to the main process in response to the request, and then the main process calculates the equivalent current of the region where the target circuit model is located according to the returned magnetic field data. The calculation of the equivalent current can refer to existing electromagnetic field related physical formulas, and the selection of the specific calculation mode is not limited in the present specification.

[0093] The main process sends the equivalent current to a second process after calculating the equivalent current of the region where the target circuit model is located. The second process stores relevant information of the target circuit model, such as a topological structure, circuit elements, and the like. After receiving the equivalent current sent by the main process, the second process calculates the port voltage based on the target circuit model. The port voltage can be calculated by using SPICE simulation or by using other self-defined methods, which are not limited in the embodiments of the present application.

[0094] After the second process calculates the port voltage of the target circuit model, the second process sends the port voltage to all the first processes. The first processes update the magnetic field data of the sub-regions according to the port voltage. Then, step 140 is performed to complete the iterative update of the equivalent current and the port voltage. The iteration condition in step 140 includes that the iteration times of steps 110 to 130 reach a preset number of times. The preset number of times can be determined by a preset time length and a time step of each iteration, so as to study the influence of the circuit at each time node in a certain time period.

[0095] Through the cooperation of the main process, the second process and all the first processes, the key links of the time-domain field-circuit simulation, such as electromagnetic field solving, circuit solving and electromagnetic field-circuit data exchange, are completed. The circuit solving is performed by the second process, the electromagnetic field calculation is performed by the first processes, and the information collection, calculation and transmission are performed by the main process. The division of labor is clear, which provides a basis for realizing the time-domain field-circuit simulation of a large-scale grid.

[0096] In an embodiment of the present application, as described in Figure 2 The first process updates the magnetic field data of the sub-region according to the port voltage, including:

[0097] In step 210, the electric field data of the sub-region is updated according to the port voltage and the topological information of the target circuit model in the sub-region.

[0098] In step 220, the magnetic field data of the sub-region is updated according to the updated electric field data of the sub-region.

[0099] The first process updates the electric field data of the sub-region according to the received port voltage and the topological information of the target circuit model stored in the sub-region. Then, the magnetic field data of the sub-region is updated according to the updated electric field data of the sub-region, which provides a basis for the next iteration of the time-domain field-circuit simulation.

[0100] In some embodiments of the present specification, the region where the target circuit model is located is divided into multiple FDTD grids, and different first processes are responsible for different FDTD grid regions. When calculating the electric field data of the grid region they are responsible for, the first process calls the overall branch information of the internally stored target circuit model and the corresponding partial branch information of the grid region they are responsible for, obtains the electric field calculation weight of the grid region they are responsible for (which can be the ratio of the branch length of the grid region they are responsible for to the overall branch length of the target circuit model), and then updates the electric field data of the grid region they are responsible for based on the port voltage and the overall branch information of the target circuit model and the calculation weight. After obtaining the updated electric field data, the first process calls the conversion algorithm of the electric field and the magnetic field again to update the magnetic field data of the grid region they are responsible for.

[0101] Considering the continuous irradiation of the preset electromagnetic signal on the target circuit model, in addition to the continuous influence of the electromagnetic signal on the circuit at the previous time of irradiation (i.e., the influence of the inventory electromagnetic signal), there is also the influence caused by the newly added electromagnetic signal. Therefore, in some embodiments of the present specification, as shown in Figure 3 The first process updates the magnetic field data of the sub-region according to the port voltage, and further includes:

[0102] Step 310: Obtain the newly added electromagnetic data of the target circuit model corresponding to the current time;

[0103] Step 320: Update the electric field data of the sub-region according to the electric field data in the newly added electromagnetic data, the port voltage, and the topology information of the target circuit model in the sub-region;

[0104] Step 330: Update the magnetic field data of the sub-region according to the magnetic field data in the newly added electromagnetic data and the updated electric field data of the sub-region.

[0105] The present embodiment takes into account the continuous irradiation of the preset electromagnetic signal, and updates the electric field data and the magnetic field data in combination with the newly added electromagnetic data, further improving the real scene of the time-domain field-circuit simulation.

[0106] When updating the electric field data of the sub-region, if there is newly added electromagnetic data generated by the continuous irradiation of the preset electromagnetic signal at the current time, the first process will also add the electric field data in the newly added electromagnetic data to the calculation. After updating the electric field data, the first process calculates the magnetic field data after electromagnetic conversion based on the electric field data, and then adds the magnetic field data in the newly added electromagnetic data to the calculation to obtain the final updated magnetic field data.

[0107] In a real environment, various media (such as air, water, etc.) are usually present near the target circuit, and the electromagnetic signal will be affected by the media when propagating in the environment. Therefore, in some embodiments of the present specification, the first process further includes: Figure 8The command vehicle shown has four antennas on the top, and the bottom end of the antennas is connected to the target circuit. The target circuit is surrounded by media such as windows and vehicle shells, which can attenuate or enhance the gain of electromagnetic signal transmission. Therefore, in some embodiments of the present specification, the media near the target circuit are also taken into account, and a plurality of different third processes are responsible for the calculation of these media (i.e., the third process is responsible for a different irradiation sub-region outside the sub-region corresponding to the first process). Specifically, as shown in Figure 4 The electromagnetic data of the target circuit model corresponding to the current time is obtained, including:

[0108] Step 410, obtaining electromagnetic information of the preset electromagnetic signal;

[0109] Step 420, sending the electromagnetic information to all first processes and all third processes, so that the first processes determine the new electromagnetic data of the target circuit model corresponding to the current time according to the electromagnetic information and the electromagnetic information returned by the third processes.

[0110] In this embodiment, the media near the target circuit are added to the calculation by a plurality of third processes, which restores the changes of the preset electromagnetic signal caused by passing through the media during transmission, and improves the accuracy of the new electromagnetic data corresponding to the target circuit model when the preset electromagnetic signal reaches the target circuit model.

[0111] For example Figure 8 As shown in the command vehicle, when the preset electromagnetic signal (i.e., the source term) irradiates the area directly above the front of the command vehicle, the electromagnetic signal may directly reach the area where the target circuit (i.e., the part connected to the bottom of the antenna) is located. At this time, all first processes can directly obtain the new electromagnetic data corresponding to the responsible sub-region. However, part of the electromagnetic signal may pass through other media near the target circuit and be refracted into the target circuit. Therefore, all third processes calculate the changes of the preset electromagnetic signal when it passes through the sub-region responsible for by themselves according to the electromagnetic information of the preset electromagnetic signal, and finally return the changed electromagnetic signal to the first process. The electromagnetic information includes incident signal and incident intensity of the electromagnetic signal.

[0112] In an embodiment of the present specification, as shown in Figure 5 The determination process of the electromagnetic information returned by the third process is given, including:

[0113] Step 510, determining the passing information of the electromagnetic signal in the process responsible sub-region according to the received electromagnetic information;

[0114] Step 520, determining the electromagnetic information returned to the next process according to the received electromagnetic information, the passing information, and the medium information of the process responsible sub-region.

[0115] The third process realizes processing of the electromagnetic signal. According to the received electromagnetic information, the passing information of the electromagnetic signal is determined, and the updated electromagnetic information is obtained according to the medium in the responsible sub-region.

[0116] Because different materials have different properties, and different properties have different influences on electromagnetic signals, the third process stores material information corresponding to the medium in the responsible sub-region. When updating the electromagnetic information, the third process adds the material information of the medium to the calculation to obtain the intensity information of the weakened or enhanced electromagnetic signal and the next propagation direction. Then, the third process returns the updated electromagnetic information to the next process. If the next process is still the third process, the next process continues to perform the same calculation and returns the updated electromagnetic information to the process after it. This process is repeated until the next process is the first process. In this way, the first process obtains new electromagnetic data, which is obtained by directly irradiating the preset electromagnetic signal and by transmitting through the medium. The real scene of the influence of the electromagnetic signal irradiation on the target circuit is restored.

[0117] In some embodiments of the present application, the sub-region responsible for the first process also includes a plurality of secondary sub-regions, that is, a plurality of secondary sub-grids. At this time, some of the secondary sub-grids do not correspond to the target circuit model (that is, there is no topology of the target circuit model on the edges of the secondary sub-grid). Therefore, when updating the electric field data of the secondary sub-region, as shown in Figure 6 The large-scale simulation of the time-domain field path also includes:

[0118] Step 610: Determine whether each secondary sub-region corresponds to the topology of the target circuit model in sequence. If the determination is yes, perform step 620. If the determination is no, perform step 630.

[0119] Step 620: Update the electric field data of the secondary sub-region according to the port voltage and the topology information of the target circuit model in the secondary sub-region.

[0120] Step 630: Obtain the magnetic field data of the adjacent region of the secondary sub-region, and update the electric field data of the secondary sub-region according to the magnetic field data of the secondary sub-region and the magnetic field data of the adjacent region.

[0121] The present embodiment relates to the case where the first process is responsible for a sub-region with a plurality of secondary sub-regions. At this time, the secondary sub-regions need to be judged. If the secondary sub-region corresponds to part of the branch of the target circuit model, the electric field data is updated and calculated based on the aforementioned port voltage. If the grid edge of the secondary sub-region does not correspond to part of the branch of the target circuit model, the electric field data is updated and calculated according to the magnetic field data of the secondary sub-region and the surrounding adjacent region. In this way, the electric field data update of the first process responsible for the multi-grid region is realized.

[0122] Similarly, considering the continuous irradiation of the preset electromagnetic signal to the target circuit model, in addition to the continuous influence of the electromagnetic signal on the circuit at the previous time of irradiation (i.e., the influence of the inventory electromagnetic signal), there is also the influence caused by the newly added electromagnetic signal, and therefore, in some embodiments of the present specification, as shown in Figure 7 As shown in the first process, when the responsible area includes a plurality of secondary sub-areas, the large-scale simulation of the time domain field path further includes:

[0123] Step 710, in turn, determine whether each secondary sub-area corresponds to the topology of the target circuit model, if the judgment is yes, execute step 720, if the judgment is no, execute step 730;

[0124] Step 720, according to the electric field data in the new electromagnetic data, the port voltage, and the topology information of the target circuit model in the secondary sub-area, update the electric field data of the secondary sub-area;

[0125] Step 730, get the magnetic field data of the adjacent area of the secondary sub-area, and update the electric field data of the secondary sub-area according to the magnetic field data of the secondary sub-area, the magnetic field data of the adjacent area, and the electric field data in the new electromagnetic data.

[0126] The present embodiment considers the continuous irradiation of the preset electromagnetic signal when the first process responsible sub-area has a plurality of secondary sub-areas, and realizes the electric field data update of the first process responsible multi-grid area by combining the branch corresponding condition of the secondary sub-area and the target circuit model and the new electromagnetic data.

[0127] In an embodiment of the present specification, a software implementation process of time domain field path simulation based on self-developed JASMIN framework is provided:

[0128] Step S1. Initialize parameters (start JASMIN framework, read in data files, create grid layer and grid piece strategy class, discrete circuit path, create circuit solver, etc.);

[0129] Step S2. The first process and the third process call numerical component 1 to realize the magnetic field update of each sub-area in space;

[0130] Step S3. The first process and the third process call numerical component 2 to realize the electric field update of each sub-area in space;

[0131] Step S4. The main process calls reduction component 1 to collect the magnetic field data of each first process, and realizes the update of the current item (i.e., the equivalent current) of the circuit port;

[0132] Step S5. The second process calls the circuit solver to realize the update of the port voltage of the circuit, and broadcasts and distributes the port voltage to all second and third processes;

[0133] Step S6. The first process calls the numerical component 3 to determine the spatial electric field of the sub-region by using the circuit port voltage.

[0134] Step S7. Determine whether the time step loop is ended, if the result is true, terminate the calculation, if the result is false, repeat steps S2-S6.

[0135] In another embodiment of the present specification, the above steps can also be adjusted as follows, and the software implementation process after adjustment includes:

[0136] Step S1. Initialize parameters (start JASMIN framework, read data files, create grid layer and grid slice strategy class, discrete circuit path, create circuit solver, etc.);

[0137] Step S2. The first process and the third process call the numerical component 1 to realize the magnetic field update of each sub-region in space;

[0138] Step S3. The main process calls the reduction component 1 to collect the magnetic field data of each first process, and realizes the update of the current item (i.e. equivalent current) of the circuit port;

[0139] Step S4. The second process calls the circuit solver to realize the update of the port voltage of the circuit, and broadcasts and distributes the port voltage to all second and third processes;

[0140] Step S5. The first process and the third process call the numerical component 2 to realize the electric field update of each sub-region in space;

[0141] Step S6. The first process calls the numerical component 3 to determine the spatial electric field of the sub-region by using the circuit port voltage.

[0142] Step S7. Determine whether the time step loop is ended, if the result is true, terminate the calculation, if the result is false, repeat steps S2-S6.

[0143] Based on the above steps, the large-scale parallel calculation of time-domain field-circuit cooperation can be realized.

[0144] In an embodiment of the present specification, a specific experiment based on a communication command vehicle is provided, as shown in Figure 8 The vehicle is 6.5 meters long, 3.2 meters wide, and 4.2 meters high, and includes a metal body, glass windows, rubber tires and other structures. Four short wave antennas are arranged at the top corners of the vehicle, and a target circuit is arranged at the antenna port. The radio frequency front-end model of the target circuit is as shown in Figure 9As shown, the band-pass filter and the amplifier are included. The standard high-altitude nuclear electromagnetic pulse (peak field 50 kV / m) plane wave is used for oblique incidence, with an elevation angle of 135°, an azimuth angle of 180°, and a polarization angle of 0°. The time-domain electromagnetic field distribution and the SPICE circuit port voltage / current are calculated. The number of calculation grids is 1300x640x840, and the number of grids is about 700 million. Parallel computing is used with 140 CPU cores.

[0145] During the experiment, the time-domain electric field distribution of the communication command vehicle coupling is obtained as shown in FIG. 10(a) and FIG. 10(b), wherein FIG. 10(a) shows the instantaneous electric field intensity distribution at the simulation time of 10 nanoseconds, and FIG. 10(b) shows the instantaneous electric field intensity distribution at the simulation time of 10 nanoseconds (from the horizontal angle, FIG. 10(a) and FIG. 10(b) show the electric field distribution of the leftmost antenna on the top of the communication command vehicle, and the deeper the color, the higher the electric field intensity); as shown in FIG. 10(a), the incident wave front is in the area far away from the vehicle, and since there is no electromagnetic coupling with the vehicle structure, the incident wave front is perpendicular to the incident direction, which conforms to the characteristics of free space plane wave transmission; as shown in FIG. 10(b), when the incident wave is close to the vehicle and the antenna structure, the reflection of the incident wave at the vehicle metal surface and the antenna structure, the transmission at the vehicle window glass medium structure, etc. can be observed, which conforms to the electromagnetic coupling characteristics of plane wave and metal, dielectric and other materials. The electric field intensity distribution conforms to the physical law.

[0146] After the time-domain field simulation of the experimental duration of 200 ns is completed, the voltage curve of the target circuit model changing with time is obtained as shown in FIG. 11(a) (in the figure, Voltage represents the voltage change per unit-meter, and time represents the affected duration of the target circuit model in ns), and the current curve of the target circuit model changing with time is obtained as shown in FIG. 11(b) (in the figure, CURRENT represents the current size, the unit is A, and time represents the affected duration of the target circuit model in ns); in the figure, S1 and S2 are the calculation results of the two antenna ports at the front of the vehicle, and S3 and S4 are the calculation results of the two antenna ports at the rear of the vehicle; in FIG. 11(a), the maximum value (absolute value) of the voltage coupled by S1 and S2 is about 120 volts (V), and the maximum value (absolute value) of the voltage coupled by S3 and S4 is about 130 volts (V); in FIG. 11(b), the maximum value (absolute value) of the current coupled by S1 and S2 is about 80 amperes (A), and the maximum value (absolute value) of the current coupled by S3 and S4 is about 130 amperes (A).

[0147] In the above experiment, the circuit models of S1, S2, S3 and S4 corresponding to the bottom connection of the antenna are consistent. At the same time, the geometric structure of the top of the communication command vehicle and the high-altitude nuclear electromagnetic pulse incidence conditions are basically kept symmetrical on the left and right. Observing Figures 11(a) and 11(b), it can be seen that the calculated results of S1 and S2 are basically consistent, and the calculated results of S3 and S4 are basically consistent, showing a bilaterally symmetrical distribution feature. Therefore, the calculation results of the coupled voltage / current are consistent with the physical laws, verifying the accuracy of the experimental calculation.

[0148] Based on the same inventive concept, the embodiments of this specification also provide a large-scale simulation device for time-domain field circuits. This device, as described in the following embodiments, addresses the issue of low accuracy in electromagnetic damage simulation for large-scale electronic platforms by simulating a preset electromagnetic signal to illuminate a target circuit model. Because the principles underlying the large-scale simulation device for time-domain field circuits are similar to those of the large-scale simulation method for time-domain field circuits, the implementation of the large-scale simulation device for time-domain field circuits can be referenced to the large-scale simulation method for time-domain field circuits, and any repetitions will not be repeated.

[0149] Specifically, such as Figure 12 As shown in the figure, the large-scale simulation device of the time domain field circuit includes:

[0150] A first sending module 1210 is configured to send a magnetic field data acquisition request of a target circuit model to all first processes, where each first process is responsible for a different sub-region of the target circuit model;

[0151] a determination module 1220, configured to determine an equivalent current of a target circuit model based on the magnetic field data returned by all first processes;

[0152] a second sending module 1230 configured to send the equivalent current to a second process, the second process being configured to determine a port voltage corresponding to a target circuit model based on the equivalent current, and to send the port voltage to all first processes, so that the first processes update magnetic field data of a responsible sub-region based on the port voltage;

[0153] The iteration module 1240 is used to iterate the above steps until an iteration condition is met.

[0154] The large-scale simulation method and device of time domain field circuits provided in the embodiments of this specification decouple the electromagnetic field operations from the circuit operations during simulation and assign them to different processes for processing, thereby improving processing efficiency; at the same time, since the electromagnetic field operations are no longer coupled with the circuit operations, it provides a basis for further encryption of the electromagnetic field area, and thus can use multi-process parallel computing to achieve high-precision simulation of the electromagnetic field area; finally, through the interactive iteration of the process of the electromagnetic field area and the circuit process, large-scale simulation of time domain field circuits is achieved, solving the problem of low accuracy of electromagnetic damage simulation of large-scale equipment electronic platforms.

[0155] An embodiment of the present specification also provides a computer device for implementing the method described in any of the above embodiments. Figure 13 As shown in FIG. 13, a structure schematic diagram of a computer device in an embodiment of the present specification is shown, the computer device 1302 can include one or more processors 1304, such as one or more central processing units (CPUs), each of which can implement one or more hardware threads. The computer device 1302 can also include any memory 1306 for storing any kind of information, such as code, settings, data, etc. Without limitation, for example, the memory 1306 can include any one or a combination of the following: any type of RAM, any type of ROM, a flash memory device, a hard disk, an optical disk, etc. More generally, any memory can store information using any technology. Further, any memory can provide volatile or non-volatile retention of information. Further, any memory can represent a fixed or removable component of the computer device 1302. In one case, the computer device 1302 can perform any operation of the associated instructions when the processor 1304 executes the associated instructions stored in any memory or combination of memories. The computer device 1302 also includes one or more drive mechanisms 1308, such as a hard disk drive mechanism, an optical disk drive mechanism, etc., for interacting with any memory.

[0156] The computer device 1302 can also include an input / output module 1310 (I / O) for receiving various inputs (via input devices 1312) and for providing various outputs (via output devices 1314). One particular output mechanism can include a presentation device 1316 and an associated graphical user interface (GUI) 1318. In other embodiments, the input / output module 1310 (I / O), the input devices 1312, and the output devices 1314 can also not be included, just as a computer device in a network. The computer device 1302 can also include one or more network interfaces 1320 for exchanging data with other devices via one or more communication links 1322. One or more communication buses 1324 couple the above-described components together.

[0157] The communication links 1322 can be implemented in any manner, for example, through a local area network, a wide area network (e.g., the Internet), a point-to-point connection, etc., or any combination thereof. The communication links 1322 can include any combination of hardwired links, wireless links, routers, gateway functionality, name servers, etc., governed by any protocol or combination of protocols.

[0158] Corresponding to Figures 1 to 7The method in this specification also provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the steps of the above method are executed.

[0159] The embodiment of this specification also provides a computer-readable instruction, wherein when the processor executes the instruction, the program therein causes the processor to execute the following Figures 1 to 7 The method shown.

[0160] It should be understood that in the various embodiments of this specification, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this specification.

[0161] It should also be understood that in the embodiments of this specification, the term "and / or" is merely a description of the relationship between associated objects, indicating that three possible relationships exist. For example, "A and / or B" can represent three situations: A exists alone, A and B exist simultaneously, and B exists alone. Furthermore, the character " / " in this specification generally indicates that the associated objects are in an "or" relationship.

[0162] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed in the embodiments of this specification can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the embodiments of this specification.

[0163] Those skilled in the art will clearly understand that, for the sake of convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0164] In several embodiments provided in the specification, it should be understood that the disclosed system, device and method can be implemented in other ways. For example, the above-described device embodiments are only illustrative, and for example, the division of the units is only a logical function division, and actual implementation can have another division manner, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the units shown or discussed can be indirect coupling or communication connection through some interfaces, devices or units, and can also be electrical, mechanical or other forms of connection.

[0165] The units described as separate components can or can not be physically separated, and the components shown as units can or can not be physical units, i.e., can be located in one place or can be distributed to multiple network units. Part or all of the units can be selected according to actual needs to achieve the purpose of the embodiments of the specification.

[0166] In addition, each functional unit in each embodiment of the specification can be integrated in one processing unit, or each unit can be physically present separately, or two or more units can be integrated in one unit. The integrated unit can be realized in the form of hardware or in the form of a software functional unit.

[0167] The integrated unit, if realized in the form of a software functional unit and sold or used as an independent product, can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the embodiments of the specification essentially or say the part of the prior art that contributes, or all or part of the technical solutions can be embodied in the form of a software product, which is stored in a storage medium and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the method described in each embodiment of the specification. The foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various program code storage media.

[0168] The principles and implementation manners of the specification are described in the specific embodiments in the specification, and the above embodiment descriptions are only used to help understand the method and core idea of the embodiments of the specification; at the same time, for those skilled in the art, according to the idea of the embodiments of the specification, the specific implementation manner and application range will be changed, and the above-mentioned content of the specification should not be understood as a limitation of the embodiments of the specification.

Claims

1. A large-scale simulation method for time domain field circuit, characterized in that: Simulating a preset electromagnetic signal to irradiate a target circuit model, the method comprising: Sending a magnetic field data acquisition request of the target circuit model to all first processes, where each first process is responsible for a different sub-area of ​​the target circuit model; determining an equivalent current of a target circuit model based on all magnetic field data returned by the first process; Sending the equivalent current to a second process, where the second process is configured to determine a port voltage corresponding to a target circuit model based on the equivalent current, and sending the port voltage to all first processes, so that the first processes update magnetic field data of their responsible sub-regions based on the port voltage; Iterate the above steps until the iteration condition is met.

2. The method according to claim 1, wherein Updating magnetic field data of the responsible sub-region according to the port voltage includes: updating the electric field data of the responsible sub-region according to the port voltage and the topology information of the target circuit model in the responsible sub-region; The magnetic field data of the responsible sub-region is updated according to the updated electric field data of the responsible sub-region.

3. The method according to claim 1, wherein Updating the magnetic field data of the responsible sub-region according to the port voltage also includes: Obtain the newly added electromagnetic data corresponding to the target circuit model at the current moment; updating the electric field data of the responsible sub-region according to the electric field data in the newly added electromagnetic data, the port voltage, and the topology information of the target circuit model in the responsible sub-region; The magnetic field data of the responsible sub-region is updated according to the magnetic field data in the newly added electromagnetic data and the updated electric field data of the responsible sub-region.

4. The method according to claim 3, wherein The method further includes a plurality of third processes, each of the third processes being responsible for irradiating a different sub-region other than the sub-region corresponding to the first process; Get the newly added electromagnetic data corresponding to the target circuit model at the current moment, including: Acquiring electromagnetic information of a preset electromagnetic signal; The electromagnetic information is sent to all first processes and all third processes, so that the first processes determine the newly added electromagnetic data corresponding to the target circuit model at the current moment according to the electromagnetic information and the electromagnetic information returned by the third processes.

5. The method according to claim 4, wherein The third process returns the determination process of electromagnetic information, including: According to the received electromagnetic information, determine the passage information of the electromagnetic signal in the sub-area responsible for this process; The electromagnetic information returned to the next process is determined based on the received electromagnetic information, the travel information and the medium information of the sub-area that the current process is responsible for.

6. The method according to claim 2, wherein The responsible sub-area of ​​the first process includes a plurality of secondary sub-areas, and the method further includes: Determine in turn whether each secondary sub-region corresponds to the topology of the target circuit model; If the judgment is yes, updating the electric field data of the secondary sub-region according to the port voltage and the topology information of the target circuit model in the secondary sub-region; If the judgment is no, the magnetic field data of the area adjacent to the secondary sub-area is obtained, and the electric field data of the secondary sub-area is updated according to the magnetic field data of the secondary sub-area and the magnetic field data of the adjacent area.

7. The method according to claim 3, wherein The responsible area of ​​the first process includes multiple secondary sub-areas, and the method further includes: Determine in turn whether each secondary sub-region corresponds to the topology of the target circuit model; If the judgment is yes, updating the electric field data of the secondary sub-region according to the electric field data in the newly added electromagnetic data, the port voltage, and the topology information of the target circuit model in the secondary sub-region; If the judgment is no, obtain the magnetic field data of the area adjacent to the secondary sub-area, and update the electric field data of the secondary sub-area based on the magnetic field data of the secondary sub-area, the magnetic field data of the adjacent area, and the electric field data in the newly added electromagnetic data.

8. The method according to claim 1, wherein The iteration conditions include: The number of iterations of the above steps reaches the preset number.

9. A large-scale simulation device for time domain field circuit, characterized in that: The device simulates a preset electromagnetic signal irradiating a target circuit model, the device comprising: A first sending module is configured to send a magnetic field data acquisition request of a target circuit model to all first processes, where each first process is responsible for a different sub-area of ​​the target circuit model; a determination module, configured to determine an equivalent current of a target circuit model according to the magnetic field data returned by all the first processes; a second sending module, configured to send the equivalent current to a second process, the second process being configured to determine a port voltage corresponding to a target circuit model based on the equivalent current, and to send the port voltage to all first processes, so that the first processes update magnetic field data of a responsible sub-region based on the port voltage; The iteration module is used to iterate the above steps until the iteration condition is met.

10. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the method according to any one of claims 1 to 8 is implemented.

11. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor of a computer device, the method according to any one of claims 1 to 8 is implemented.

12. A computer program product, comprising a computer program, characterized in that: When the computer program is executed by a processor of a computer device, the method according to any one of claims 1 to 8 is implemented.

Images