JP2021051467A - Analysis apparatus, analysis method and analysis program - Google Patents

Abstract

To shorten calculation time associated with readout of data from a storage unit for storing data.SOLUTION: An analysis apparatus that uses an FDTD method, comprises: a storage unit for storing a medium constant defined per each block for each of a plurality of models; a calculation unit for calculating an electromagnetic field component per each block using the medium constant; an area division unit for dividing a finite area to be analyzed into a common part, which consists of a plurality of blocks with the medium constant that is common among the plurality of models, and a non-common part, which consists of the plurality of blocks with different medium constants for each model; and a control unit for controlling the calculation unit to calculate to update the electromagnetic field component of the common part for the plurality of models by reading the medium constant for the common part from the storage unit once, and to calculate to update the electromagnetic field components for the non-common parts for each model by reading the medium constants for the non-common parts from the storage unit for each model.SELECTED DRAWING: Figure 6

Description

The present invention relates to an analyzer, an analysis method and an analysis program.

As the speed of wireless communication systems increases, there is an increasing need to properly analyze the electromagnetic fields generated by devices. For example, in the FDTD method (Finite-difference time-domain method), the analysis area is divided into a large number of blocks (cells), and the electric field and magnetic field are time-updated up to a predetermined number of repetitions for each block using Maxwell's equations. This is an electromagnetic field analysis method.

In the FDTD method, when the electric field and the magnetic field are time-updated, the electric field and the magnetic field in the previous time step, and the medium constants such as the conductivity and the magnetic permeability of each block are used. Therefore, the calculation for electromagnetic field analysis requires a memory for holding the electric field, the magnetic field, and the medium constant according to the number of blocks (see, for example, Patent Documents 1 to 4).

Japanese Unexamined Patent Publication No. 2012-003546 Japanese Unexamined Patent Publication No. 2014-006851 Japanese Unexamined Patent Publication No. 2015-022693 Japanese Unexamined Patent Publication No. 2017-011518

Conventionally, in the FDTD method, when there are a plurality of models in which the medium constant is set for each block, the medium constant must be read from a storage unit such as a memory for each model, and the time of the electromagnetic field must be updated again from the initial stage (0 time). There was a problem that it took a long time because it had to be done.

An object of the present invention is to provide an analysis device, an analysis method, and an analysis program capable of shortening the calculation time associated with reading data from a storage unit that stores data.

The analysis device according to one aspect of the present invention is an analysis device that performs electromagnetic field analysis using the FDTD method, and includes a storage unit that stores a medium constant determined for each block for each of a plurality of models, and the storage unit. A calculation unit that calculates the electromagnetic field component for each block using the medium constant stored in, and an intersection consisting of multiple blocks that have a common medium constant among multiple models for the finite region to be analyzed. By reading the area division unit that divides the area into non-common parts composed of a plurality of blocks having different medium constants for each model and the medium constants for the common parts once from the storage unit, the calculation unit makes the common parts. By performing a calculation to update the electromagnetic field component for a plurality of models and reading the medium constant for the non-intersection from the storage unit for each model, the calculation unit obtains the electromagnetic field component for the non-intersection. It is characterized by having a control unit that controls to perform a calculation to be updated for each model.

Further, the analysis method according to one aspect of the present invention is an analysis method in which electromagnetic field analysis is performed using the FDTD method, in which a finite region to be analyzed is divided into a plurality of blocks having a common medium constant among a plurality of models. A region division step of dividing the area into a common part and a non-common part consisting of a plurality of blocks having different medium constants for each model, and a storage unit for storing the medium constants determined for each block for each of the plurality of models. By reading the medium constant for the common portion once from, the electromagnetic field component for the common portion is calculated for a plurality of models, and the medium constant for the non-common portion is read from the storage unit for each model. It is characterized by including a calculation step of calculating the electromagnetic field component for the non-intersection for each model.

Further, the analysis program according to one aspect of the present invention causes a computer to function as each part of the above-mentioned analysis device.

According to the present invention, it is possible to shorten the calculation time associated with reading data from a storage unit that stores data.

(A) is a diagram illustrating a state in which an analysis area is divided into a large number of blocks when there is no person indoors. FIG. (B) is a diagram illustrating a state in which the analysis area when there is a person indoors is divided into a large number of blocks. It is a figure which shows typically the 1st example of the method of analyzing an electromagnetic field with respect to an analysis area using two models. It is a flowchart which shows the processing example which the analysis apparatus performs in order to carry out the method of analyzing the electromagnetic field shown in FIG. It is a figure which shows typically the 2nd example of the method of analyzing the electromagnetic field with respect to the analysis area using two models. It is a flowchart which shows the processing example which the analysis apparatus performs in order to carry out the method of analyzing the electromagnetic field shown in FIG. It is a block diagram which illustrates the function which the analysis apparatus which concerns on one Embodiment has. It is a figure which shows typically the method which the analyzer analyzes the electromagnetic field with respect to the analysis area using two models. It is a flowchart which shows the processing example which the analysis apparatus performs in order to carry out the method of analyzing the electromagnetic field shown in FIG. 7. It is a figure which shows the hardware configuration example of the analysis apparatus which concerns on one Embodiment.

An example of electromagnetic field analysis using the FDTD method will be described below. FIG. 1 is a diagram illustrating a state in which an analysis region when the FDTD method is used is divided into a large number of blocks in order to compare the difference in radio wave propagation characteristics depending on the presence or absence of a person (human body) indoors (indoors). FIG. 1A is a diagram illustrating a state in which an analysis area is divided into a large number of blocks when there is no person indoors. FIG. 1B is a diagram illustrating a state in which the analysis area when there is a person indoors is divided into a large number of blocks.

As shown in FIG. 1A, in the analysis region 1, for example, the indoor space is divided by a large number of blocks 2. In each block 2, the medium constant is set for each model. Further, the analysis region 1 is an analysis target in which absorbing boundary conditions are set around it, and the number of blocks 2 is finite.

It is assumed that the analysis region 1 is provided with a wave source 3 serving as a transmission point of an electromagnetic field, for example, at one upper corner. In the analysis region 1, when the medium constant is uniform, the radio wave is propagated in a spherical shape (radial in the electric field strength contour diagram) centered on the wave source 3. Further, it is assumed that the medium constant for an object fixed in the room does not change.

On the other hand, when there is a person in the analysis area 1, the electromagnetic field characteristics (shielding in radio wave propagation, etc.) in the analysis area 1 change depending on the position and posture of the person. This is because the medium constant of a person is different from the medium constant of the surroundings of a person.

Further, when a plurality of models are applied to the analysis region 1 in which a person exists in the room, the medium constants used for each of the plurality of models are almost the same in most of the blocks 2.

For example, as shown in FIG. 1 (b), in the analysis region 1, the intersection 4 in which the medium constants used in each of the plurality of models are common and the medium constants used in each of the plurality of models are different for each model. It is possible to divide the area into the common portion 5. Here, the common part 4 corresponds to an area where there are no people or the like. Further, the non-common portion 5 corresponds to an area occupied by a person or the like.

Next, a method of analyzing the electromagnetic field with respect to the analysis region 1 using a plurality of models will be described. Here, for example, a computer is used as an analysis device.

FIG. 2 is a diagram schematically showing a first example of a method of analyzing an electromagnetic field with respect to an analysis region 1 using two models. As shown in FIG. 2, first, the analyzer reads the medium constant for the first model determined for each block from the memory, and calculates the electric field and the magnetic field for each block using the first model. Then, the analysis device performs a calculation to update the electromagnetic field of each block at predetermined time intervals using the first model. At this time, the analyzer reads the medium constant from the memory every time the electromagnetic field is updated.

When the calculation of the first model is completed, the analyzer reads the medium constants for the second model defined for each block from the memory, and calculates the electric field and the magnetic field for each block using the second model. Then, the analysis device performs a calculation to update the electromagnetic field of each block at predetermined time intervals using the second model. At this time, the analyzer reads the medium constant from the memory every time the electromagnetic field is updated.

In this way, when the analysis device analyzes the electromagnetic field using a plurality of models, the calculation is completed for each single model, for example. Then, in each model, the analysis device reads the medium constant from the memory and updates the electric field / magnetic field.

In the example shown in FIG. 2, assuming that the number of models is M, the number of time steps until convergence of each model is T, and the number of blocks of each model is N, the number of times the memory is accessed to read the medium constant from the memory is , M × T × N (times).

FIG. 3 is a flowchart showing an example of processing performed by the analysis apparatus in order to carry out the method of analyzing the electromagnetic field shown in FIG.

First, the analysis device updates the calculation target model (S100) and inputs three-dimensional (i, j, k) absolute coordinate values to the model (S102).

Here, i is an x-axis index (i = 1 ... Nx), j is a y-axis index (j = 1 ... Ny), and k is a z-axis index (k = 1 ... Nz). To do.

Next, the analysis device assigns the medium constants corresponding to each coordinate value as follows, for example (S104).
W (1,1,1) = ε ₀ (air)
_{W (1, 1,} 2) = ε 1 (Glass)
W (1,1,3) = ε ₀ (air)
W (1,1,4) = ε ₀ (air)
W (1,2,1) = ε ₂ (human body)
・
・
・
W (Nz, Ny, Nz) = ε ₀ (air)

Then, the analysis device updates the time step (S106). Next, the analysis device reads out the medium constant of the calculation target block (S108). At this time, the analysis device reads out the medium constants corresponding to the above-mentioned coordinate values.

Subsequently, the analysis device reads the electromagnetic field component of the previous time step (S110), calculates and updates the electromagnetic field component (S112), and then determines whether or not the calculation of all blocks is completed (S114). ). When the analyzer determines that the calculation of all blocks is completed (S114: Yes), the process proceeds to the process of S116, and when it is determined that the calculation of all blocks is not completed (S114: No), the analysis device proceeds to S108. Return to processing.

In the process of S116, the analyzer determines whether or not the electric field component calculation has converged. When the analyzer determines that the electric field component calculation has converged (S116: Yes), it proceeds to the processing of S118, and when it determines that the electric field component calculation has not converged (S116: No), it proceeds to the processing of S106. Return.

In the process of S118, the analysis device determines whether or not the calculation has been completed for all the models. When the analyzer determines that the calculation has been completed for all models (S118: Yes), the processing is completed, and when it is determined that the calculation has not been completed for all models (S118: No), the processing of S100 is completed. Return to.

FIG. 4 is a diagram schematically showing a second example of a method of analyzing an electromagnetic field with respect to an analysis region 1 using two models. As shown in FIG. 4, when analyzing an electromagnetic field using a plurality of models, the analysis device reads the medium constants of all models from the memory and sequentially updates the electric and magnetic fields of all models in each time step. Calculations may be performed. At this time, the analyzer reads the medium constants of all the models from the memory every time the electromagnetic field is updated.

In the example shown in FIG. 4, assuming that the number of models is M, the number of time steps until convergence of each model is T, and the number of blocks of each model is N, the number of times the memory is accessed to read the medium constant from the memory is , M × T × N (times).

FIG. 5 is a flowchart showing an example of processing performed by the analysis apparatus in order to carry out the method of analyzing the electromagnetic field shown in FIG.

First, the analysis device reads all the calculation models (S200), and inputs the three-dimensional (i, j, k) absolute coordinate values of each model to each model (S202).

Here, i is an x-axis index (i = 1 ... Nx), j is a y-axis index (j = 1 ... Ny), and k is a z-axis index (k = 1 ... Nz). To do.

Next, the analysis device assigns the medium constants corresponding to the coordinate values of each model as follows (S204).

For example, the analyzer assigns the following medium constants to the first model.
W ₁ (1,1,1) = ε ₀ (air)
W _{1 (1, 1,} 2) = ε ₁ (Glass)
W ₁ (1,1,3) = ε ₀ (air)
W ₁ (1,1,4) = ε ₀ (air)
W _{1 (1, 2,} 1) = ε 2 (human body)
・
・
・
W ₁ (Nz, Ny, Nz) = ε ₀ (air)

Further, the analyzer assigns the following medium constants to, for example, the second model.
W ₂ (1,1,1) = ε ₁ (Glass)
W _{2 (1, 1,} 2) = ε 1 (Glass)
W ₂ (1,1,3) = ε ₀ (air)
W ₂ (1,1,4) = ε ₀ (air)
W ₂ (1,2,1) = ε ₂ (human body)
・
・
・
W ₂ (Nz, Ny, Nz) = ε ₂ (human body)

When the number of models is M, the analysis device allocates, for example, WM (Nz, Ny, Nz) = ε ₁ (Glass).

Then, the analysis device updates the time step (S206). Next, the analysis device reads out the medium constant of the calculation target block (S208). At this time, the analysis device reads out the medium constant corresponding to each coordinate value of each model described above.

Subsequently, the analysis device reads the electromagnetic field component of the previous time step (S210), calculates and updates the electromagnetic field component (S212), and then determines whether or not the calculation of all blocks in each model is completed. (S214). When the analyzer determines that the calculation of all blocks in each model has been completed (S214: Yes), the process proceeds to S216, and when it is determined that the calculation of all blocks in each model has not been completed (S214: Yes). In No), the process returns to S208.

In the process of S216, the analyzer determines whether or not the electric field component calculation has converged. The analyzer ends the process when it is determined that the electric field component calculation has converged (S216: Yes), and returns to the process of S206 when it determines that the electric field component calculation has not converged (S216: No). ..

Next, an analysis device that shortens the calculation time involved in reading data from the storage unit that stores data and analyzes the electromagnetic field with respect to the analysis region 1 using a plurality of models will be described.

FIG. 6 is a block diagram illustrating the functions of the analysis device 10 according to the embodiment. As shown in FIG. 6, the analysis device 10 includes a storage unit 11, a parameter setting unit 12, an array initialization unit 13, an array data reading unit 14, a calculation unit 15, a region division unit 16, and a control unit 17, and includes an FDTD. Perform electromagnetic field analysis using the method. For example, the analysis device 10 analyzes the electromagnetic field with respect to the analysis region 1 using a plurality of models.

The storage unit 11 has an analytic space information storage unit 110 and a block model storage unit 112, and is composed of, for example, a memory capable of reading (writing) and reading data via the bus 100 and the bus 102. .. Further, the storage unit 11 is also designed to store the results processed by the calculation unit 15 and the area division unit 16.

The analysis space information storage unit 110 stores information about the analysis area 1. For example, the analytic space information storage unit 110 stores the medium constants of each of the above-mentioned blocks 2 for each model. The block model storage unit 112 stores a plurality of models such as the first model and the second model described above.

The parameter setting unit 12 sets the analysis space size, block (cell) size, discrete time interval definition, and analysis model (antenna, space conditions, etc.) for the analysis area 1 for analysis, and sets the analysis model (antenna, space conditions, etc.) via the bus 100. The set value is stored in the storage unit 11.

The array initialization unit 13 secures an array of a size required for performing the analysis for calculating the electric field and the magnetic field in the analysis space, and stores the array in the storage unit 11 via the bus 100.

The array data reading unit 14 reads the array data used for the calculation and stores it in the storage unit 11 via the bus 100.

The calculating unit 15 includes an electric field component calculating unit 150, an electric field absorbing boundary condition calculating unit 152, a magnetic field component calculating unit 154, and a magnetic field absorbing boundary condition calculating unit 156.

The electric field component calculation unit 150 reads the medium constant for each block from the analysis space information storage unit 110, calculates the electric field component for each block 2 for each model read from the block model storage unit 112, and calculates the calculation result via the bus 102. Is stored in the storage unit 11.

The electric field absorbing boundary condition calculation unit 152 calculates the absorbing boundary condition of the electric field component in the analysis region 1, performs the calculation of applying the absorbing boundary condition to the electric field component calculated by the electric field component calculating unit 150, and via the bus 102. The calculation result is stored in the storage unit 11.

The magnetic field component calculation unit 154 reads the medium constant for each block from the analytic space information storage unit 110, calculates the magnetic field component for each block 2 for each model read from the block model storage unit 112, and calculates the calculation result via the bus 102. Is stored in the storage unit 11.

The magnetic field absorbing boundary condition calculation unit 156 calculates the absorbing boundary condition of the magnetic field component in the analysis region 1, performs a calculation to apply the absorbing boundary condition to the magnetic field component calculated by the magnetic field component calculating unit 154, and performs a calculation via the bus 102. The calculation result is stored in the storage unit 11.

The region division unit 16 divides the analysis region 1, which is a finite region to be analyzed, into a common portion composed of a plurality of blocks 2 having a common medium constant among a plurality of models, and a plurality of blocks having a different medium constant for each model. The area is divided into a non-common part composed of 2, and the result of the area division via the bus 102 is stored in the storage unit 11.

The control unit 17 controls the storage unit 11, the calculation unit 15, and the area division unit 16 via the bus 102. For example, the control unit 17 controls the calculation unit 15 to perform a calculation for updating the electromagnetic field component for the common portion for a plurality of models by reading the medium constant for the common portion described above once from the storage unit 11. To do. Further, the control unit 17 reads the medium constant for the non-common portion described above from the storage unit 11 for each model, so that the calculation unit 15 performs a calculation for updating the electromagnetic field component for the non-common portion for each model. Control.

Next, a specific operation example of the analyzer 10 will be described with reference to FIGS. 7 and 8. FIG. 7 is a diagram schematically showing a method in which the analysis device 10 analyzes an electromagnetic field with respect to the analysis region 1 shown in FIG. 1 (b) using two models.

As shown in FIG. 7, when analyzing an electromagnetic field using a plurality of models, the analysis device 10 reads the medium constants of all the models from the memory and sequentially updates the electric and magnetic fields of all the models in each time step. Do the calculation. At this time, the analysis device 10 reads the medium constant for each block 2 of the common portion 4 in all models once from the storage unit 11 made of memory, using the wave source 3 (see FIG. 1B) as a reference point, for each model. The medium constant for each block 2 of the non-intersection 5 is read from the storage unit 11. Further, the analysis device 10 reads the medium constant from the memory every time the electromagnetic field is updated.

That is, even if the number of models is plural, the analysis device 10 reads the medium constant for the common portion 4 from the memory once, reads only the medium constant for the non-common portion 5 from the memory for each model, and calculates the electromagnetic field. Do.

At this time, the control unit 17 (FIG. 6) controls the common portion 4 and the non-common portion 5 so that the boundary conditions at the boundary surfaces of each other are common. Further, the calculation unit 15 calculates the electromagnetic field component in parallel for each model.

FIG. 8 is a flowchart showing an example of processing performed by the analyzer 10 in order to carry out the method of analyzing the electromagnetic field shown in FIG. 7.

First, the analysis device 10 reads all the calculation models (S300), and inputs the three-dimensional (i, j, k) absolute coordinate values of each model to each model (S302).

Here, i is an x-axis index (i = 1 ... Nx), j is a y-axis index (j = 1 ... Ny), and k is a z-axis index (k = 1 ... Nz). To do.

Next, the analysis device 10 assigns the medium constants corresponding to the coordinate values of each model as follows (S304).

For example, the analysis device 10 assigns the following medium constants to the intersection 4 (FIG. 1 (b)) based on the coordinate values.
_{W (1, 1,} 2) = ε 1 (Glass)
W (1,1,3) = ε ₀ (air)
W (1,1,4) = ε ₀ (air)
・
・
・

Further, the analysis device 10 assigns the following medium constants to the non-intersection 5 based on the coordinate values. For example, the analyzer 10 assigns the following medium constants to the first model of the non-intersection 5.
W ₁ (1,1,1) = ε ₀ (air)
・
・
・
W ₁ (Nz, Ny, Nz) = ε ₀ (air)

Further, the analysis device assigns the following medium constants to the second model of the non-intersection 5 based on the coordinate values.
W ₂ (1,1,1) = ε ₁ (Glass)
・
・
・
W ₂ (Nz, Ny, Nz) = ε ₂ (human body)

When the number of models is M, the analysis device 10 allocates, for example, WM (Nz, Ny, Nz) = ε ₁ (Glass).

Next, in the analysis device 10, the region dividing unit 16 divides the common portion 4 and the non-common portion 5 into regions based on the medium constant (S306).

Then, the analysis device 10 updates the time step (S308). Next, the analyzer 10 reads the electromagnetic field component of the previous time step (S310), calculates and updates the electromagnetic field component for the common part 4 (S312), and then sets the common part 4 and the non-common part 5 together. The boundary conditions are shared (intersectioned) based on the evaluation frequency and the medium constant (S314).

Next, the analysis device 10 calculates and updates the electromagnetic field component for the non-intersection 5 (S316), and then shares the boundary conditions between the common part 4 and the non-common part 5 based on the evaluation frequency and the medium constant. (Common) (S318).

After that, in the process of S320, the analysis device 10 determines whether or not the electric field component calculation has converged. The analysis device 10 ends the process when it is determined that the electric field component calculation has converged (S320: Yes), and proceeds to the process of S308 when it determines that the electric field component calculation has not converged (S320: No). Return.

In this way, the analysis device 10 collectively reads the medium constants from the storage unit 11 for the intersection 4 once in a plurality of models, calculates the time update of the electromagnetic field in parallel, and determines the number of accesses to the storage unit 11. It is reduced and the calculation time associated with reading data from the storage unit 11 is shortened.

When the medium of the non-intersection 5 is a perfect conductor, or the evaluation frequency band is the 5th generation mobile communication (5G) candidate frequency band, the medium of the non-intersection 5 has a complex dielectric constant of about that of a human body. In the case, the analyzer 10 may consider only the boundary condition between the common portion 4 and the non-common portion 5, and omit the calculation of the electromagnetic field inside the non-common portion 5.

Each function of the analysis device 10 may be partially or wholly configured by hardware, or may be configured as a program executed by a processor such as a CPU.

That is, the analysis device 10 according to the present invention can be realized by using a computer and a program, and the program can be recorded on a storage medium or provided through a network.

FIG. 9 is a diagram showing a hardware configuration example of the analysis device 10 according to the embodiment. As shown in FIG. 9, the analysis device 10 has, for example, an input unit 20, an output unit 21, a communication unit 22, a CPU 23, a memory 24, and an HDD 25 connected via a bus 26, and has a function as a computer. Further, the analysis device 10 is capable of inputting / outputting data to / from the storage medium 27.

The input unit 20 is, for example, a keyboard, a mouse, or the like. The output unit 21 is a display device such as a display. The communication unit 22 is, for example, a wired and wireless network interface.

The CPU 23 controls each part constituting the analysis device 10 and performs the above-mentioned calculation and the like. The memory 24 and the HDD 25 form a storage unit 11 for storing data. In particular, the memory 24 stores each data used in the above-mentioned calculation. The storage medium 27 is made capable of storing an analysis program or the like for executing the function of the analysis device 10. The architecture constituting the analysis device 10 is not limited to the example shown in FIG.

The embodiments described above are exemplary and not limited to the embodiments of the present invention, and the present invention can also be implemented in various other modifications and modifications.

1 ... Analysis area, 2 ... Block, 3 ... Wave source, 4 ... Common part, 5 ... Non-common part, 10 ... Analysis device, 11 ... Storage unit, 12. .. Parameter setting unit, 13 ... Array initialization unit, 14 ... Sequence data reading unit, 15 ... Calculation unit, 16 ... Area division unit, 17 ... Control unit, 20 ... Input unit, 21 ... Output unit, 22 ... Communication unit, 23 ... CPU, 24 ... Memory, 25 ... HDD, 26 ... Bus, 27 ... Storage medium, 110.・ ・ Analysis space information storage unit, 112 ・ ・ ・ block model storage unit, 150 ・ ・ ・ electric field component calculation unit, 152 ・ ・ ・ electric field absorption boundary condition calculation unit, 154 ・ ・ ・ magnetic field component calculation unit, 156 ・ ・ ・Magnetic field absorption boundary condition calculation unit

Claims (7)

In an analyzer that performs electromagnetic field analysis using the FDTD method,
A storage unit that stores medium constants defined for each block for each of multiple models,
A calculation unit that calculates the electromagnetic field component for each block using the medium constants stored in the storage unit,
Region division that divides the finite region to be analyzed into a common part consisting of a plurality of blocks having a common medium constant among a plurality of models and a non-common part consisting of a plurality of blocks having a different medium constant for each model. Department and
By reading the medium constant for the common portion once from the storage unit, the calculation unit performs a calculation to update the electromagnetic field component for the common portion with respect to a plurality of models, and obtains the medium constant for the non-common portion. An analyzer characterized by having a control unit that controls the calculation unit to perform a calculation for updating the electromagnetic field component for the non-intersection for each model by reading from the storage unit for each model. The control unit
The analysis apparatus according to claim 1, wherein the common portion and the non-common portion are controlled so that the boundary conditions at the interface between the common portions and the non-common portions are controlled to be common. The calculation unit
The analysis apparatus according to claim 1 or 2, wherein the electromagnetic field components are calculated in parallel for each model. In an analysis method that performs electromagnetic field analysis using the FDTD method,
Region division that divides the finite region to be analyzed into a common part consisting of a plurality of blocks having a common medium constant among a plurality of models and a non-common part consisting of a plurality of blocks having a different medium constant for each model. Process and
By reading the medium constant for the common part once from the storage unit that stores the medium constant determined for each block for each of the plurality of models, the electromagnetic field component for the common part is calculated for the plurality of models. An analysis method including a calculation step of calculating the electromagnetic field component for the non-intersection for each model by reading the medium constant for the non-intersection for each model from the storage unit. The calculation process is
The analysis method according to claim 4, wherein the electromagnetic field component is calculated for the common portion and the non-common portion by making the boundary conditions at the boundary surfaces common to each other. The calculation process is
The analysis method according to claim 4 or 5, wherein the electromagnetic field components are calculated in parallel for each model. An analysis program for operating a computer as each part of the analysis device according to any one of claims 1 to 3.

Images