JP2005308487A - Electromagnetic field analyzer, electromagnetic field analyzing method, computer program and computer readable recording medium - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To easily and certainly analyze the electromagnetic field generated in large-scaled equipment. <P>SOLUTION: The central part 20e1 of a magnetic body 20, the average magnetic flux density B<SB>ave</SB>generated in the equivalent element 31 occupied by the air on the lateral side of the magnetic body 20 and average magnetic field H<SB>ave</SB>are operated. At this time, so far as the value of magnetic flux density and the value in the direction of the easy magnetization axis X of a magnetic field are respectively constant on the line 22a in the equivalent element 31, the average magnetic flux density B<SB>ave</SB>, the average magnetic flux density B<SB>ave</SB>and the average magnetic field H<SB>ave</SB>become effective. Then, a B-H curves 90-93 showing the relation between the average magnetic field H<SB>ave</SB>and the average magnetic flux density B<SB>ave</SB>and the B-θ curves 100-103 showing the relation between the average magnetic flux density B<SB>ave</SB>and an angle θ<SB>BH</SB>are formed using angle θ<SB>B</SB>as a parameter. These formed curves are used to calculate electromagnetic field in a region 21 to be analyzed. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

  The present invention relates to an electromagnetic field analysis device, an electromagnetic field analysis method, a computer program, and a computer-readable recording medium, and in particular, used for analyzing an electromagnetic field generated in a region where a plurality of substances including a magnetic material exist. Is preferred.

  Generally, when an electric device is used, a magnetic field is generated around it. Therefore, a magnetic shield device is used to prevent the magnetic field generated from the electrical equipment from leaking to the surroundings.

And with the development of technology in recent years, facilities using large electric devices that use large currents are increasing. For such a large electric device, it is necessary to use a large magnetic shield device.
In addition, when a large current flows near the magnetic shield device or there is a magnetic material having a residual magnetic field, a magnetic field is generated therefrom. Such a magnetic field is a hindrance when trying to perform a precise measurement of a precision device within a magnetic shield device. It is also necessary to configure the magnetic shield device so that an external magnetic field does not enter.

  For example, 200 directional electrical steel sheets having a side length of 2 [m] and a thickness of 1 [mm] are prepared, and these 200 directional electrical steel sheets are arranged in a bowl shape at intervals. A large-sized magnetic shield device configured to have a size of 2 [m] square as a whole has been proposed (see, for example, Patent Document 1).

  By the way, in order to verify the usefulness of such a large magnetic shield device, it is desirable to numerically analyze the electromagnetic field generated in the large magnetic shield device. In such a case, conventionally, the finite element method is generally used. In this finite element method, an electromagnetic field is analyzed by dividing a region to be analyzed into a number of regions (cells) having a relatively simple shape.

  As described above, since the magnetic shield device is configured using a steel plate having magnetic anisotropy, when analyzing the electromagnetic field using the finite element method, it is necessary to make the region to be analyzed very fine. There is. More specifically, it is necessary to divide the region to be analyzed into a large number of cells having a size of about 0.2 [mm] square.

JP 2002-164686 A

However, when an electromagnetic field generated in a large-sized magnetic shield device as in the above-described example is analyzed by the conventional finite element method, a cubic region having a side of 2 [m] is converted into a region of 0.2 [mm] square. Must be divided into (cells). Therefore, the number of areas to be analyzed is about 1 trillion.
However, the capacity of the main storage device in the current personal computer is about 4 [GB] at the maximum. Accordingly, the maximum number of regions (cells) that can be analyzed is about one million.

  As described above, when an electromagnetic field in a large-scale magnetic shield device is analyzed using the conventional finite element method, the amount of data necessary for the analysis far exceeds the storage capacity of the personal computer. For this reason, there is a problem that it is extremely difficult to analyze an electromagnetic field generated in a large-scale facility such as the above-described magnetic shield device.

  The present invention has been made in view of the above-described problems, and an object thereof is to easily and reliably analyze an electromagnetic field generated in a large-scale facility.

  The electromagnetic field analysis apparatus of the present invention includes an average magnetic field calculation means for calculating an average magnetic flux density and an average magnetic field in an equivalent element occupied by a plurality of substances including a magnetic material, the average magnetic flux density, and the average magnetic field. And an electromagnetic field analysis means for analyzing an electromagnetic field generated in an analysis target area wider than the equivalent element, and the average magnetic field calculation means is an easy-magnetization-axis direction of magnetic characteristics at a predetermined position in the equivalent element. In this state, the average magnetic flux density and the average magnetic field in the equivalent element are calculated.

  The electromagnetic field analysis method of the present invention includes an average magnetic field calculation step for calculating an average magnetic flux density and an average magnetic field in an equivalent element occupied by a plurality of substances including a magnetic substance, the average magnetic flux density, and the average magnetic field. And an electromagnetic field analysis step for analyzing an electromagnetic field generated in an analysis target area wider than the equivalent element, and the average magnetic field calculation step includes a direction of an easy magnetization axis of magnetic characteristics at a predetermined position in the equivalent element. In this state, the average magnetic flux density and the average magnetic field in the equivalent element are calculated.

The computer program of the present invention includes an average magnetic field calculation step for calculating an average magnetic flux density and an average magnetic field in an equivalent element occupied by a plurality of substances including a magnetic material, the average magnetic flux density, and the average magnetic field. And an electromagnetic field analysis step for analyzing an electromagnetic field generated in an analysis target area wider than the equivalent element, and the average magnetic field calculation step includes an easy magnetization axis of magnetic characteristics at a predetermined position in the equivalent element. It is characterized in that the direction value becomes constant, and in this state, the average magnetic flux density and the average magnetic field in the equivalent element are calculated.
A computer-readable recording medium according to the present invention records the above-described computer program.

According to the present invention, an average magnetic flux density and an average magnetic field in an equivalent element occupied by a plurality of substances including a magnetic material are calculated, and the equivalent magnetic flux is calculated using the calculated average magnetic flux density and the average magnetic field. Since the electromagnetic field generated in the analysis target area wider than the element is obtained, the electromagnetic field generated in the analysis target area can be analyzed without dividing the analysis target area into a large number of elements as in the prior art. As a result, the storage capacity required for analyzing the electromagnetic field can be greatly reduced, and the electromagnetic field in a large-scale facility that has been difficult to analyze can be reliably analyzed.
In the present invention, when calculating the average magnetic flux density and the average magnetic field in the equivalent element as described above, the value of the magnetic property in the easy axis direction of the magnetic characteristic at a predetermined position in the equivalent element becomes constant. Since the state is set, it is possible to prevent the influence of the demagnetizing field and the edge generated inside the magnetic body in the equivalent element from appearing in the analysis result of the electromagnetic field generated in the analysis target region. Thereby, the electromagnetic field generated in the analysis target region can be accurately obtained.

(First embodiment)
Next, a first embodiment of the present invention will be described with reference to the drawings.
FIG. 1 is a block diagram showing an example of the configuration of the electromagnetic field analyzer in the present embodiment. In the present embodiment, as shown in FIG. 2, an electromagnetic field generated in an analysis target region 21 including five magnetic bodies (steel plates) 20 arranged in a bowl at regular intervals (24 [mm] intervals). The case of analyzing the above will be described as an example. The analysis target area 21 is a rectangular area having a size of 4000 [mm] in the vertical direction and 1155 [mm] in the horizontal direction. Each magnetic body (steel plate) 20 has a width of 1000 [mm] and a thickness of 1 [mm].

In FIG. 1, the electromagnetic field analysis device 1 includes an operation unit 2, a display unit 3, and a processing unit 4.
The operation unit 2 is a device configured by a keyboard, a mouse, and the like, and is a device for transmitting the contents executed by the user (analyzer) to the processing unit 4.

  The display unit 3 is a device configured with a display or the like, and is a device for displaying a processing result or the like executed by the processing unit 4. The user operates the operation unit 2 and inputs desired content while viewing the content displayed on the display unit 3.

  The processing unit 4 is a computer that includes a CPU, a ROM, a RAM, and the like. The processing unit 4 performs a processing operation in the electromagnetic field analysis device 1 by executing a program recorded in the ROM.

Specifically, the processing unit 4 includes an average magnetic field calculation unit 4a, a magnetic characteristic curve creation unit 4b, and a magnetic field distribution calculation unit 4c.
The average magnetic field calculation unit 4a regards a predetermined area in the analysis target area 21 shown in FIG. 2 as an equivalent element, and calculates an average magnetic flux density B _ave and an average magnetic field H _ave in the equivalent element. To do. In addition, the average magnetic field calculation unit 4a also calculates an angle θ _BH formed by the calculated average magnetic flux density B _ave and average magnetic field H _ave as the magnetic characteristics of equivalent elements. In the present invention, the above equivalent elements are referred to as equivalent elements. Further, in the above, bold characters indicate vectors.

  As shown in FIG. 2, in the analysis target region 21, an equivalent element 31 composed of a magnetic body (steel plate) 20 and a nonmagnetic body (air) is repeatedly present. That is, in FIG. 2, there are ten equivalent elements 31 in the X direction and five in the Y direction. Therefore, in the present embodiment, the magnetic body 20 in the analysis target region 21 in FIG. 2 is replaced with the equivalent element 31, and the electromagnetic field in the analysis target region 21 is calculated using the magnetic characteristics of the equivalent element 31. In this way, it is not necessary to divide the analysis target area 21 of FIG. 2 into a large number of areas as in the prior art. In the electromagnetic field analysis apparatus 1 of the present embodiment, the portions for obtaining the magnetic characteristics of the equivalent element 31 are the average magnetic field calculation section 4a and the magnetic characteristic curve creation section 4b of the processing section 4 in FIG. The part that analyzes the electromagnetic field in the analysis target region 21 with the number of divisions in the analysis target region 21 significantly reduced by using the characteristics is the magnetic field distribution calculation unit 4 c of the processing unit 4.

Here, the equivalent element analysis region and equivalent elements in the present embodiment will be described in detail with reference to FIG. 3 and FIG.
3, the magnetic properties of the equivalent elements 31 (e.g., the average magnetic flux density B _ave relationships and the average magnetic field H _ave, the average magnetic flux density B _ave and the average magnetic field H _ave the angle theta _BH of the average magnetic flux density B _ave It is the figure which showed an example of the analysis model for calculating | requiring (relationship). As shown in FIG. 3, the equivalent element 31 is the same as the equivalent element 31 of FIG. 2, and is equivalent to the magnetic body (steel plate) 20 e that is one of the plurality of magnetic bodies (steel plates) 20 shown in FIG. 2. It is a two-dimensional region (plane) formed by the central portion 20e1 and air on its side, and has a size of 100 [mm] in the vertical direction and 25 [mm] in the horizontal direction (see FIG. 4). Here, the case where the equivalent element 31 is a part of the magnetic body 20e of FIG. 2 is described as an example, but the magnetic characteristics of the other magnetic bodies 20 of FIG. Therefore, in the analysis of FIG. 2 by this method, the same magnetic characteristic data calculated from FIG. 3 may be used.

  The equivalent element analysis region 30 is a two-dimensional region (plane) including the equivalent element 31 in the center, and has a size of 10000 [mm] in the vertical direction and 10000 [mm] in the horizontal direction. In the present embodiment, 41 magnetic bodies, which are the same as the magnetic bodies 20e, are arranged in a bowl shape at equal intervals (24 [mm] intervals) in the equivalent element analysis region 30 having such a size. The model is the analysis model. A two-dimensional region formed by the central portion of the magnetic body 20e located in the middle of these 41 magnetic bodies and the air on the side thereof is used as an equivalent element 31 (see FIG. 4). Thus, in the present embodiment, an electromagnetic field generated on a plane defined by the easy magnetization axis X and the hard magnetization axis Y of the magnetic body (steel plate) 20 is analyzed.

By the way, when an external magnetic field H _ext is applied to the equivalent element 31 in order to obtain the magnetic characteristics of the equivalent element 31, the magnetic body 20e1 in the equivalent element 31 is magnetized, but the shape of the magnetic body in the equivalent element analysis region 30 is obtained. Depending on the number, the influence of the demagnetizing field and the edge may appear, and predetermined magnetic characteristics may not be obtained. That is, the equivalent element analysis region 30, the number of magnetic bodies in the equivalent element analysis region 30, the length of the magnetic body in the equivalent element analysis region 30 in the easy axis direction, and the equivalent element analysis region 30 shown in FIG. Unless the distance from the end of the magnetic body inside to the end of the equivalent element analysis region is sufficiently increased, the influence of the demagnetizing field and the end generated inside the magnetic body (steel plate) 20e appears in the obtained electromagnetic field. . More specifically, the magnetic flux density and magnetic field in the equivalent element 31 are affected by the end of the magnetic body (steel plate) 20e, and the magnetic flux density and magnetic field distribution are not flat. For this reason, in the present embodiment, as shown in FIG. 3, a wide area in which the number and size of the magnetic materials are noted to such an extent that the influence of the demagnetizing field and the edge does not occur is defined as an equivalent element analysis area 30. The electromagnetic field in the element 31 is prevented from being affected as much as possible by the demagnetizing field and the edge generated inside the magnetic body (steel plate) 20e.
The details of the demagnetizing field are described in, for example, “Near-angle“ Physics of Ferromagnetic Material (above) ”裳 華 房 Physics Selections 4, 1978”, so only the outline is described here. Here, since the case where the equivalent element exists repeatedly is considered, if the repeatability is insufficient, an “end effect” appears. For this reason, the influence of the edge must be eliminated as much as possible to obtain the magnetic characteristics of the equivalent element.

Then, the average magnetic field calculation unit 4a according to the present embodiment calculates the magnetic flux density and the magnetic field in the equivalent element analysis region 30 set by the user in consideration of such a demagnetizing field and the influence of the edges.
In the present embodiment, whether or not the electromagnetic field in the equivalent element 31 is affected by the demagnetizing field and the edge generated inside the magnetic body (steel plate) 20e is determined as follows.
First, the average magnetic field calculation unit 4a passes through the center 22 of the magnetic body (steel plate) 20e in the equivalent element 31 (the center of the equivalent element 31 in the example of FIG. 3) 31a on the line 22 parallel to the hard axis Y. A value of the magnetic flux density in the easy magnetization axis X direction and a value of the magnetic field in the easy magnetization axis X direction are obtained. Then, the value in the equivalent element 31 is extracted from the value of the obtained magnetic flux density in the easy magnetization axis X direction and the value of the magnetic field in the easy magnetization axis X direction. When the value of the extracted magnetic flux density in the easy axis X direction is constant and the value of the extracted magnetic field in the easy axis X direction is constant, the equivalent element analysis region 30 set by the user is effective. It shall be

  On the other hand, if at least one of the value of the extracted magnetic flux density in the easy axis X direction and the value of the extracted magnetic field in the easy axis X direction is not constant, the equivalent element analysis region 30 is not appropriate. A message to the effect is displayed on the display unit 3. After this display, the user operates the operation unit 2 to reset the equivalent element analysis region 30. Then, the average magnetic field calculation unit 4a re-determines the electromagnetic field generated in the reset equivalent element analysis region 30 as described above, and determines again whether or not the equivalent element analysis region 30 is valid.

  Further, in the present embodiment, as shown in FIG. 5, the equivalent element analysis region 30 of FIG. 3 is divided into a plurality of square or rectangular regions so that the equivalent element analysis region 30 can be analyzed by the electromagnetic field analysis by the finite element method. The area 30 is divided. Each of the plurality of divided areas has an appropriate size that can be calculated by a finite element method (FEM). In the following description, this divided area is referred to as a divided area. In this embodiment, the finite element method is used as an electromagnetic field analysis method. However, even if it is used in other numerical analysis methods such as a difference method, Fourier transform, and Fourier series, it is the same as when using the finite element method. The electromagnetic field can be analyzed.

The average magnetic field calculation unit 4a obtains the magnetic flux density B and the magnetic field H generated in each of the plurality of divided regions when the external magnetic field H _ext is applied to the equivalent element analysis region 30 set as described above. The external magnetic field H _ext is expressed as (Formula 1) below.
H _ext = {H _extX , H _extY } (1 set)

Thus, the external magnetic field H _ext is a two-dimensional vector having a value H _EXTx easy axis direction, and a value H _ExtY the hard axis direction. In the present embodiment, the direction of the external magnetic field H _ext is specified by the angle φ formed by the external magnetic field H _ext and the easy axis X (see FIGS. 2 to 5).

When such an external magnetic field H _ext is applied, the magnetic flux density B (i, j) and magnetic field H (i, j) generated in the divided areas shown by the oblique lines in FIG. (Expression 3)

B (i, j) = {B _X (i, j), B _Y (i, j)} (Expression 2)
H (i, j) = {H _X (i, j), H _Y (i, j)} (Expression 3)

In the above, i and j are natural numbers for specifying the location of the divided region.
As described above, the magnetic flux density B (i, j) and the magnetic field H (i, j) generated in each divided region are also expressed as values B _X (i, j) and H _X (i, j) in the easy axis direction. It is a two-dimensional vector having values B _Y (i, j) and H _Y (i, j) in the magnetic field hard axis direction. The magnetic flux density B (i, j) and magnetic field H (i, j) generated in each divided region can be obtained by electromagnetic field analysis by a finite element method based on Maxwell's electromagnetic equations. Specifically, the equation regarding the static magnetic field shown in the following (formula 4) is used.

Here, [μ] ⁻¹ is the reciprocal of the magnetic permeability. A is a vector potential, and this vector potential A is defined as in the following (formula 5).

J ₀ is an applied current excited in the electric device. Based on this equation, using the interpolation function in the spatially discretized divided region as shown in FIG. 5, using the variational method or the Galerkin method, the constitutive equation is obtained, and the Gaussian elimination method or the ICCG method is used. Thus, the vector potential A in the divided region is obtained, and then the magnetic flux density B (i, j) and the magnetic field H (i, j) there are obtained. The details of the finite element method are described in, for example, “Nakada, Takahashi,“ Fine Element Method of Electrical Engineering, Second Edition ”, Morikita Publishing, 1982”, so only an outline is given here.

  The average magnetic field calculation unit 4a then generates the magnetic flux density generated in the equivalent element 31 from the magnetic flux density B (i, j) generated in each divided region and the magnetic field H (i, j) obtained as described above. B (i, j) and magnetic field H (i, j) are extracted.

Then, the average magnetic flux density B _ave and the average magnetic field H _ave in the equivalent element 31 are obtained from the extracted magnetic flux density B (i, j) and the magnetic field H (i, j). Specifically, it is obtained by the following (formula 6) to (formula 11).

In the above, B _X (i, j) is the value of the magnetic flux density B (i, j) in the easy axis direction. B _Y (i, j) is the value of the magnetic flux density B (i, j) in the hard axis direction.

H _X (i, j) is a value in the direction of the easy axis of the magnetic field H (i, j). H _Y (i, j) is the value of the magnetic field H (i, j) in the hard axis direction.
ΔS (i, j) is the size (area) of the divided region (see FIG. 5).

In this embodiment, the average magnetic flux density _Bave calculated as described above is regarded as the magnetic flux density generated in the equivalent element 31. The average magnetic field H _ave is regarded as a magnetic field generated in the equivalent element 31.

In the present embodiment, the direction of the average magnetic flux density B _ave is specified by the angle θ _B formed by the average magnetic flux density B _ave and the easy axis X. Further, the direction of the average magnetic field H _ave is specified by the angle θ _H formed by the average magnetic field H _ave and the easy magnetization axis X (see FIG. 4).

  Further, when the magnetic flux density B (i, j) and the magnetic field H (i, j) generated in each divided region are obtained as described above, a magnetic body (steel plate) measured in advance as shown in FIG. ) Use 20 BH curves. In FIG. 6, the BH curve 51 plotted with black circles is a BH curve in the easy magnetization axis direction of the magnetic bodies 20 to 24. A BH curve 52 plotted with a white square is a BH curve in the hard axis direction of the magnetic body (steel plate) 20.

  From the above, the magnetic characteristics in the equivalent element 31 are obtained by the following (Expression 12) to (Expression 14).

Then, the size B _ave average magnetic flux density, so that the angle theta _B becomes a desired value, the external magnetic field H _ext magnitude and direction (the external magnetic field H _ext, the angle φ between the axis of easy magnetization X) The calculation according to the above (Formula 2) to (Formula 14) is repeated.
More specifically, for example, the angle φ formed by the external magnetic field H _ext and the easy magnetization axis X is set to 0, 15, 30, 45, 60, 75, 90 [°], and the external magnetic field at each angle φ. An average magnetic flux density B _ave and an average magnetic field H _ave are obtained when the size of H _ext is 0, 100, 500, 1000, 2000, 5000, 10000 [Gauss].

As described above, the average magnetic field calculation unit 4a according to the present embodiment has the following four relations, respectively, _indicating the magnitude of the external magnetic field H _ext and the angle φ formed by the external magnetic field H _ext and the easy axis X. As a parameter.
1. External magnetic field H _{ext -angle} θ _B
2. External magnetic field H _{ext −average} magnetic flux density B _ave
3. External magnetic field H _{ext -average} magnetic field H _ave
4). External magnetic field H _{ext -angle} θ _BH
In this embodiment, the external magnetic field H _ext, the angle theta _B, the average magnetic flux density B _ave, the average magnetic field H _ave, and was to obtain the relationship between the angle theta _BH, and the external magnetic flux density B _ext, You may make it obtain _| require the relationship with angle (theta) _B , average magnetic flux density _Bave , average magnetic field _Have , and angle (theta) _BH .

The magnetic characteristic curve creation unit 4b creates a BH curve in the equivalent element 31 based on the average magnetic flux density B _ave and the average magnetic field H _ave calculated by the average magnetic field calculation unit 4a. That is, a curve representing the relationship between the average magnetic field H _ave and the average magnetic flux density B _ave is created. Then, the created BH curve is recorded on a recording medium. An example of a specific method for creating the BH curve in the equivalent element 31 will be described with reference to FIGS.

First, as shown in FIG. 7, curves 60 to 66 representing the relationship between the angle θ _B between the average magnetic flux density B _ave and the easy magnetization axis X and the external magnetic field H _ext are represented by the external magnetic field H _ext and the easy magnetization axis. An angle φ formed with X is created as a parameter.

From these curves 60 to 66, the average magnetic flux density B _ave and the average magnetic field H when the angle θ _B formed by the average magnetic flux density B _ave and the easy magnetization axis X is 0, 15, 30, 90 [°]. _Ask for _ave . Needless to say, the angle θ _B formed between the average magnetic flux density B _ave and the easy magnetization axis X used at this time is not limited to 0, 15, 30, 90 [°].

More specifically, for example, as shown in FIG. 8, a line 70 having an angle θ _B of 30 [°] between the average magnetic flux density B _ave and the easy axis X, and a curve (for example, a curve) shown in FIG. 64) is obtained. Then, the average magnetic flux density B _ave and average magnetic field H _{ave at} the obtained intersection 71 are obtained.

At this time, performs linear interpolation using the average magnetic flux density B _ave before and after the point 72 and 73 of the intersection 71 and the average magnetic field H _ave, computing the the average magnetic flux density B _ave at the intersection 71 and the average magnetic field H _ave Can do.

Then, to create a B-H curve 80-83 representing the relationship between the average magnetic flux density B _ave and the average magnetic field H _ave obtained as described above (see Figure 9). Further, a leveling process is performed on the BH curves 80 to 83. In the present embodiment, the BH curves 90 to 93 that have been subjected to the leveling process are defined as BH curves in the equivalent element 31 (see FIG. 10).

The magnetic characteristic curve creation unit 4b creates a B-θ curve in the equivalent element 31 based on the average magnetic flux density B _ave and the average magnetic field H _ave calculated by the average magnetic field calculation unit 4a. That is, a curve representing the relationship between the angle θ _BH formed by the average magnetic flux density B _ave and the average magnetic field H _ave (see FIG. 4) and the average magnetic flux density B _ave is created. Then, the created B-θ curve is recorded on a recording medium.

More specifically, as shown in FIG. 11, the relationship between the average magnetic flux density B _ave obtained by using the above-described (Equation 12) and the angle θ _BH formed by the average magnetic field H _ave and the average magnetic flux density B _ave. B-θ curves 100 to 103 representing the angle θ _B formed by the average magnetic flux density B _ave and the easy magnetization axis X are created as parameters. In the example shown in FIG. 11, curves 100 to 103 are created when the angle θ _B formed between the average magnetic flux density B _ave and the easy axis X is 0, 15, 30, 90 [°]. The angle θ _B formed between the average magnetic flux density B _ave used as a parameter and the easy axis X is the same value as the angle θ _B used when creating the BH curves 90 to 93 shown in FIG. Use it.

  The magnetic field distribution calculation unit 4c generates an electromagnetic field (in the analysis target region 21) based on the BH curves 90 to 93 and the B-θ curves 100 to 103 created by the magnetic characteristic curve creation unit 4b as described above. Find the magnetic flux lines. Then, the magnetic flux density and magnetic field at a desired position in the analysis target region 21 are obtained from the obtained magnetic flux lines. Note that the electromagnetic field (magnetic flux lines) generated in the analysis target region 21 is obtained using a finite element method (FEM). The method of electromagnetic field analysis by the finite element method is as described above.

  Further, when using the finite element method, it is preferable to apply the side element finite element method. In this embodiment, a two-dimensional model is used. However, in a three-dimensional model described later, it is particularly important to perform an electromagnetic field analysis by applying side elements. This is because the amount of data used for the analysis is smaller when the edge element is applied than when the node element is applied. More specifically, when a node element is applied, the number (this embodiment) is obtained by multiplying the number of vertices (nodes) in the element (4 in this embodiment) by the number of dimensions (2 in this embodiment). Then, data based on 8) is required. On the other hand, when the side element is applied, there may be data based on the number of sides constituting the element (6 in the present embodiment).

  In addition, there is an advantage that the electromagnetic field (magnetic flux lines) can be obtained more accurately when the side element finite element method is applied than when the nodal element finite element method is applied. More specifically, for example, a nodal element finite element is applied to a rectangular region (element) at the boundary between different substances (in this embodiment, a magnetic body (steel plate) and air) on one of two opposing sides. When the method is applied, the physical quantity obtained from each node is a mixture of the physical quantities of the different substances. On the other hand, when the side element finite element method is applied, the physical quantity obtained from one of the two opposing sides is a mixture of physical quantities of the different substances, The physical quantity obtained from the other is the physical quantity of one substance. Therefore, in the rectangular region (element) as described above, the electromagnetic field (flux lines) can be obtained more accurately by applying the side element finite element method than by applying the nodal element finite element method.

Next, an example of the processing operation in the electromagnetic field analyzer 1 of the present embodiment will be described with reference to the flowchart of FIG.
First, in the first step S1, the average magnetic field calculation unit 4a sets the magnitude of the external magnetic field H _ext applied to the equivalent element analysis region 30 to an initial value (0 [Gauss]).

Next, in step S2, the average magnetic field calculation unit 4a sets an angle φ formed by the external magnetic field H _ext and the easy magnetization axis X to an initial value (0 [°]).
Thus, the direction of the external magnetic field H _ext is set by setting the angle φ formed by the external magnetic field H _ext and the easy magnetization axis X.

Next, in step S3, the average magnetic field calculation unit 4a sets the size and position of the equivalent element analysis region 30 according to the setting from the user. Further, the average magnetic field calculation unit 4a sets the boundary condition of the equivalent element analysis region 30 according to the external magnetic field _Hext . More specifically, in the present embodiment, when the external magnetic field H _ext is applied, for example, a value corresponding to the external magnetic field H _ext is inserted into the boundary element in the equivalent element analysis region 30 shown in FIG. To set boundary conditions.

Next, in step S4, the average magnetic field calculation unit 4a calculates the magnetic flux density B (i, j) and magnetic field H (i, j) generated in the equivalent element analysis region 30 in accordance with the setting contents of the external magnetic field _Hext .
More specifically, the above-described finite element method electromagnetic field analysis is performed on the equivalent element analysis region 30 with the external magnetic field H _ext as a boundary condition.

  Next, in step S5, the average magnetic field calculation unit 4a determines the magnetic flux density generated in the equivalent element 31 from the magnetic flux density B (i, j) and the magnetic field H (i, j) obtained by the processing in step S4. B (i, j) and magnetic field H (i, j) are extracted.

Next, in step S6, the average magnetic field calculation unit 4a uses the magnetic flux density B (i, j) and the magnetic field H (i, j) extracted in step S5 to determine the magnetic characteristics (in the equivalent element 31). Of the average magnetic flux density B _ave (i, j) and average magnetic field H _ave (i, j)) are obtained and stored in the recording medium. Further, an angle θ _BH formed by the obtained average magnetic flux density B _ave and average magnetic field H _ave is obtained and stored in the recording medium.

Next, in step S7, the average magnetic field calculation unit 4a determines whether or not the magnitude of the external magnetic field _Hext is a specified value.
In the present embodiment, the magnitude of the external magnetic field H _ext is increased in the order of 0, 100, 500, 1000, 2000, 5000, and 10000 [Gauss]. Therefore, 10000 [Gauss] is the specified value.

As a result of the determination in step S7, if the magnitude of the external magnetic field H _ext is not a specified value, the process proceeds to step S8 and the set value is changed. For example, when the set value of the magnitude of the external magnetic field H _ext is 0 [Gauss], the set value is changed to 100 [Gauss]. Then, the process returns to step S3.
On the other hand, when the set value of the magnitude of the external magnetic field H _ext is a specified value, the process proceeds to step S9 and the set value is reset to the initial value.

Then, the determination in step S10, the average magnetic field arithmetic section 4a, the set value of the direction of the external magnetic field H _ext (angle between the external magnetic field H _ext magnetization easy axis X phi) is, whether the specified value To do.
In the present embodiment, the angle φ formed by the external magnetic field H _ext and the easy axis X is increased in the order of 0, 15, 30, 45, 60, 75, 90 [°]. Therefore, 90 [°] is the specified value.

Result of judgment at step S10, the set value of the direction of the external magnetic field H _ext (angle between the external magnetic field H _ext magnetization easy axis X phi) is, if not specified value, the setting value proceeds to step S11 change. For example, when the set value of the angle φ formed by the external magnetic field H _ext and the easy axis X is 0 [°], the set value is changed to 15 [°]. Then, the process returns to step S3.

On the other hand, the set value of the direction of the external magnetic field H _ext (angle between the external magnetic field H _ext magnetization easy axis X phi) is, when it is the specified value is reset to the initial value settings proceeds to step S12 .

  Next, in step S13, the average magnetic field calculation unit 4a determines the magnetic flux on the line 22a parallel to the hard axis Y that passes through the center 31a of the equivalent element 31 from the electromagnetic field generated in the equivalent element analysis region 30 obtained in step S4. A value of the density in the easy magnetization axis X direction and a value of the magnetic field in the easy magnetization axis X direction are extracted. Then, it is determined whether or not the value of the extracted magnetic flux density in the easy axis X direction and the value of the magnetic field in the easy axis X direction are constant.

  If the result of this determination is not constant, the process proceeds to step S18, where the average magnetic field calculation unit 4a displays on the display unit 3 that the equivalent element analysis region 30 is not appropriate, and resets the equivalent element analysis region 30. Prompt the user and end the process.

On the other hand, if it is constant, the process proceeds to step S14, where the magnetic characteristic curve creation unit 4b determines from the average magnetic flux density B _ave (i, j) and the average magnetic field H _ave (i, j) obtained in step S6. The BH curves 90 to 93 in the equivalent element 31 are created.

Next, in step S15, the magnetic characteristic curve creation unit 4b determines from the average magnetic flux density B _ave (i, j) obtained in step S6 and the angle θ _BH formed by the average magnetic flux density B _ave and the average magnetic field H _ave. The B-θ curves 100 to 103 in the equivalent element 31 are created.

Next, in Step S16, the magnetic field distribution calculation unit 4c uses the BH curves 90 to 93 and the B-θ curves 100 to 103 in the equivalent element 31 created in Steps S13 and S14 to analyze the analysis target region. The electromagnetic field generated at 21 is obtained.
Finally, in step S17, the magnetic field distribution calculation unit 4c displays the electromagnetic field obtained in step S15 on the display unit 3 and ends the process.

Next, an example of the processing operation in the electromagnetic field analyzer 1 when creating the BH curve in step S14 in FIG. 12 will be described with reference to the flowchart in FIG.
First, in the first step S21, the magnetic characteristic curve creating unit 4b represents a curve representing the relationship between the external magnetic field H _ext and the angle θ _B formed by the average magnetic flux density B _ave obtained in step S6 and the easy axis X. 60 to 66 are created by performing the linear interpolation described above.

Next, in step S22, the magnetic characteristic curve generator unit 4b, from the curve 60 to 66 obtained in step S21, the angle theta _B between the average magnetic flux density B _ave magnetization easy axis X, 0, 15, 30, The average magnetic flux density B _ave and average magnetic field H _ave at 90 [°] are obtained.
Finally, in step S23, the magnetic characteristic curve generator unit 4b, the B-H curve 90 to 93 representing the relationship between the average magnetic flux density B _ave and the average magnetic field H _ave obtained in step S22, the leveling process described above Create by doing.

Next, an example of the processing operation in the electromagnetic field analyzer 1 when creating the B-θ curve in step S15 in FIG. 12 will be described with reference to the flowchart in FIG.
First, in the first step S31, the magnetic characteristic curve creation unit 4b determines that the angle θ _B formed by the average magnetic flux density _Bave obtained in step S6 and the easy magnetization axis X is 0, 15, 30, 90 [°]. The angle θ _BH formed by the average magnetic flux density B _ave and the average magnetic field H _ave is obtained.

In step S32, the magnetic characteristic curve creating unit 4b determines that the angle θ _B formed by the average magnetic flux density _Bave and the easy axis X obtained in step 22 in FIG. 13 is 0, 15, 30, 90 [°. ] and the average magnetic flux density B _ave of time, determined in step S31, the angle theta _B between the average magnetic flux density B _ave easy magnetization axis X is the average magnetic flux when the 0,15,30,90 [°] B-θ curves 100 to 103 representing the relationship between an angle θ _BH formed by the density B _ave and the average magnetic field H _ave are created.

The above-described method (equivalent BH method) for obtaining an electromagnetic field (for example, magnetic flux lines) generated in the analysis target region 21 is summarized as shown in FIG.
The inventors of the present application use the method of the present embodiment as described above, even if the number of regions (number of elements) to be divided is smaller (in the order of two digits) than the conventional method. It was confirmed that substantially the same results as those obtained by the method were obtained.

Next, it will be shown how the value of the magnetic flux density in the easy magnetization axis X direction and the value of the magnetic field in the easy magnetization axis X direction change depending on the size of the equivalent element analysis region.
FIG. 16 shows an example of the relationship between the value of the magnetic flux density in the equivalent element analysis region in the easy magnetization axis X direction and the distance in the hard magnetization axis Y direction from the center 31a of the magnetic body (steel plate) 20e1 in the equivalent element 31. FIG. FIG. 17 shows the value in the easy magnetization axis X direction of the magnetic field (magnetic field strength) obtained by the magnetic field distribution calculation unit 4c and the distance in the hard magnetization axis Y direction from the center of the magnetic body (steel plate) 20e1 in the equivalent element 31. It is the figure which showed an example of the relationship.

A characteristic line 151 shown in FIG. 16A shows a relationship when the equivalent element analysis region 30 shown in FIG. 3 is obtained. In addition, a characteristic line 152 shown in FIG. 16B shows a relationship when the equivalent element analysis region 171 shown in FIG. Furthermore, a characteristic line 153 shown in FIG. 16C shows a relationship when the equivalent element analysis region 172 shown in FIG. 18B is obtained. The characteristic lines 151 to 153 indicate the characteristics when the average magnetic flux density B _ave in the equivalent element is 0.05 [T].

  As shown in FIG. 18A, the equivalent element analysis region 171 is a square region having a length of 10000 [mm] in the vertical and horizontal directions. Ten magnetic bodies (steel plates) arranged in a bowl shape at equal intervals (24 [mm] intervals) are arranged in the equivalent element analysis region 171. Each magnetic body (steel plate) has a width of 1000 [mm] and a thickness of 1 [mm]. Thus, the equivalent element analysis region 171 has the same size as the equivalent element analysis region 30 shown in FIG. However, the number of magnetic bodies (steel plates) disposed in the equivalent element analysis region 171 is smaller than the number of magnetic bodies (steel plates) disposed in the equivalent element analysis region 30 shown in FIG.

  As shown in FIG. 18B, the equivalent element analysis area 172 is a rectangular area having a length of 2000 [mm] and a width of 3000 [mm]. Ten magnetic bodies (steel plates) arranged in a bowl shape at equal intervals (24 [mm] intervals) are arranged in the equivalent element analysis region 172. Each magnetic body (steel plate) has a width of 1000 [mm] and a thickness of 1 [mm]. Thus, the equivalent element analysis region 172 is smaller in size than the equivalent element analysis region 30 shown in FIG. Further, the number of magnetic bodies (steel plates) disposed in the equivalent element analysis region 172 is smaller than the number of magnetic bodies (steel plates) disposed in the equivalent element analysis region 30 shown in FIG.

As shown in FIG. 16, when the equivalent element analysis region 30 is set as shown in FIG. 3, the value B _x of the magnetic flux density in the easy axis X direction is constant over a wide range, whereas FIG. When the equivalent element analysis regions 171 and 172 are set as shown in FIG. 5, it can be seen that the value B _x of the magnetic flux density in the easy axis X direction is not constant.

Similarly, as shown in FIG. 17, when the equivalent element analysis region 30 is set as shown in FIG. 3, the value H _{x in the} easy axis X direction of the magnetic field (magnetic field strength) is constant over a wide range. On the other hand, when the analysis target regions 171 and 172 are set as shown in FIG. 18, it can be seen that the value H _{x of the} magnetic field (magnetic field strength) in the easy magnetization axis X direction is not constant.

  As described above, when the equivalent element analysis region 30 is set as shown in FIG. 3, the influence of the demagnetizing field and the edge generated in the magnetic body (steel plate) 20 is not manifested in the obtained electromagnetic field. When the equivalent element analysis regions 171 and 172 are set as shown in FIG. 18, it can be seen that the demagnetizing field generated in the magnetic bodies (steel plates) 173 and 174 and the influence of the edges are manifested in the obtained electromagnetic field.

  From the results shown in FIG. 16 and FIG. 17, not only the size of the equivalent element analysis region but also the equivalent element so that the influence of the demagnetizing field and the edge generated in the magnetic body (steel plate) is not manifested in the obtained electromagnetic field. It can be seen that it is necessary to set the equivalent element analysis region in consideration of the number and size of the magnetic bodies (steel plates) included in the analysis region. The demagnetizing field generated in the magnetic body (steel plate) becomes smaller as the length of the magnetic body (steel plate) in the easy magnetization axis X direction becomes longer, and the magnetic body (steel plate) arranged around the magnetic body (steel plate). This is because the larger the number, the smaller, and the smaller the distance between the end of the equivalent element analysis region and the end of the magnetic body (steel plate), the smaller.

As described above, in this embodiment, a predetermined area occupied by the central portion 20e1 of the magnetic body (steel plate) 20 that repeatedly exists at a predetermined interval and the air on the side thereof is equivalent to one equivalent. The average magnetic flux density B _ave generated in the equivalent element 31, the average magnetic field H _ave, and the angle θ _BH formed by the average magnetic flux density B _ave and the average magnetic field H _ave are calculated, and the calculated average magnetic field H _ave the average magnetic flux density B _ave BH curve representing the relationship between the 90 to 93, and the average magnetic flux density B _ave and the angle theta _BH between the average magnetic field H _ave, B-theta representing the relationship between the average magnetic flux density B _ave The curves 100 to 103 are created using the angle θ _B as a parameter, and the electromagnetic fields (magnetic flux lines) in the analysis target region 21 are obtained using the created BH curves 90 to 93 and the B-θ curves 100 to 103. So do n’t divide it into many elements , It is possible to analyze the electromagnetic field generated in the analysis target area 21. As a result, the storage capacity required for analyzing the electromagnetic field can be greatly reduced, and the electromagnetic field in a large-scale facility that has been difficult to analyze can be reliably analyzed.
In this way, the electromagnetic field analysis method is called an equivalent BH method in that the electromagnetic field in the analysis target region 21 is analyzed using equivalent magnetic characteristics.

  In particular, in the present embodiment, a predetermined value on a line 22a in the equivalent element 31 that passes through the center 31a (the center of the equivalent element 31) 31a of the magnetic body (steel plate) 20e1 in the equivalent element 31 and is parallel to the hard axis Y. The equivalent element analysis region 30 is set so that the value of the magnetic flux density in the easy magnetization axis X direction and the value of the magnetic field (magnetic field strength) in the easy magnetization axis X direction are constant. . That is, since the size of the equivalent element analysis region is set in consideration of the number and size of the magnetic bodies (steel plates) included in the equivalent element analysis region, the demagnetizing field generated in the magnetic body (steel plate) and The influence of the edge can be prevented from appearing in the obtained electromagnetic field, and the electromagnetic field in the analysis target region 21 can be obtained accurately.

In this embodiment, in the magnetic field distribution calculation unit 4c, the value of the magnetic flux density in the equivalent element 31 in the easy axis X direction and the value of the magnetic field in the equivalent element 31 in the easy axis X direction are constant. Only in certain cases, the set equivalent element analysis region (the calculated average magnetic flux density B _ave and average magnetic field H _ave ) is made effective, but it is not always necessary to do so. For example, even when one of the value of the magnetic flux density in the equivalent element 31 in the easy axis X direction and the value of the magnetic field in the equivalent element 31 in the easy axis X direction is not constant, the set equivalent element is set. The analysis area may be effective. In this way, it is possible to inform the user how much the demagnetizing field generated in the magnetic body (steel plate) 20 affects the analysis result.

  In the present embodiment, when the value of the magnetic flux density in the equivalent element 31 in the easy axis X direction and the value of the magnetic field in the equivalent element 31 in the easy axis X direction are obtained, the value on the line 22 is obtained. The value on the line 22a is extracted from the above, but the value on the line 22a may be obtained directly.

In the present embodiment, the case where the electromagnetic field is analyzed using the average magnetic flux density B _ave as a reference (B center) has been described, but the average magnetic field H _ave may be used as a reference (H center). For example, in this embodiment, the angle theta _B between the average magnetic flux density B _ave easy axis of magnetization X as a parameter, has been to create a B-H curve 90 to 93, easy magnetization average magnetic field H _ave These may be created using the angle θ _H formed with the axis X as a parameter.

Moreover, you may make it obtain | require the iron loss which arises when the external magnetic field _Hext is given. Specifically, for example, the iron loss is obtained based on the BH curves 90 to 93 created by the magnetic characteristic curve creation unit 4b. Further, a curve representing the relationship between the obtained iron loss and the average magnetic flux density B _ave is created using an angle θ _B formed by the average magnetic flux density B _ave and the easy magnetization axis X as a parameter, and the created curve is used as a recording medium. It may be recorded.

  Furthermore, the equivalent element analysis region is not limited to those shown in FIGS. 2 and 3 as long as it is a region occupied by a plurality of substances including a magnetic material. That is, as long as it is a representative region in the analysis target region 21, the position and size of the equivalent elements are not limited to those shown in FIGS. 2 and 3, and the magnetic field distribution in the analysis target region 21 is appropriately set. It can be determined as appropriate within the range obtained.

For example, as shown in FIG. 19, when analyzing an electromagnetic field generated in a region to be analyzed 181 configured using a magnetic body (steel plate) 180 processed into a hollow rectangular parallelepiped shape, the electromagnetic field analysis apparatus of the present embodiment 1 can be applied. In this case, the equivalent element 182 is, for example, a two-dimensional region (plane) formed by a part of the magnetic body (steel plate) 180 and air on its side. Then, an external magnetic field H _ext is applied to the equivalent element 182 and the average magnetic flux density B _ave and the average magnetic field H _ave are calculated as described above. Further, the average magnetic flux density B _ave and the angle θ _BH formed by the average magnetic field H _ave are calculated. Then, magnetic flux lines generated in the analysis target area 181 are obtained.

Here, the external magnetic field H _ext is given by, for example, an exciter 183 disposed at the center of a hollow region surrounded by a magnetic body (steel plate) 180 as shown in FIG. More specifically, an external magnetic field H _{ext is applied} to the equivalent element 182 by passing a current through a coil 183b wound around the core 183a. The magnitude and direction of the external magnetic field H _ext can be changed by changing the magnitude of the current in the coil 183b and the angle φ ′ of the core 183a.

  In addition, the actual shield device is reduced (for example, (1/16) times or (1/8) times) the shield device to be analyzed, and the magnetic flux lines generated in the shield device to be analyzed are You may make it obtain | require as mentioned above.

  Note that the above-described equivalent element analysis region may be composed of a plurality of types of equivalent elements, and in that case, there are a plurality of types of equivalent elements. In this case, since there are a plurality of types of equivalent elements, there are also a plurality of types of equivalent BH curves corresponding to FIG.

  In the present embodiment, the process of step S13 for determining whether or not the magnetic characteristics in the equivalent element 31 are constant is performed after step S12. However, it is not always necessary to perform the process at this timing. For example, the above processing may be performed after step S5. In this case, if the magnetic characteristics in the equivalent element 31 are constant, the process proceeds to step S6, and if not, the process may proceed to proceed to step S18.

  In the present embodiment, the two-dimensional region formed by the central portion 20e1 of the magnetic body 20e and the air on the side thereof is used as the equivalent element 31, but the substance included in the equivalent element is not limited to these. . That is, any region may be used as an equivalent element as long as it is a region occupied by a magnetic body and a non-magnetic body. Here, the magnetic body is not limited to a ferromagnetic body such as iron or cobalt, but may be a magnetic body other than a ferromagnetic body such as a paramagnetic body such as iron oxide. The non-magnetic material is not limited to air but refers to a non-magnetic material such as a copper coil, wood, and resin.

In this embodiment, the angle θ _BH is considered in order to improve the accuracy of the analysis result of the electromagnetic field. However, the angle θ _BH may be ignored depending on the simple formula and the analysis accuracy to be obtained.
In addition, the operation unit 2, the display unit 3, and the processing unit 4 in the electromagnetic field analysis device 1 may be in the same space, or may be coupled via a network.

(Second Embodiment)
Next, a second embodiment of the present invention will be described. In the first embodiment described above, when the electromagnetic field in the equivalent element is affected by the demagnetizing field and the edge generated in the magnetic body (steel plate), the user resets the equivalent element analysis region. . On the other hand, in the present embodiment, the electromagnetic field analysis device automatically resets the equivalent element analysis region. As described above, the present embodiment and the first embodiment described above differ only in the method of setting the equivalent element analysis region, and therefore, the same parts as those of the first embodiment described above. These are given the same reference numerals as those shown in FIGS.

The configuration of the electromagnetic field analyzer of this embodiment is the same as that shown in FIG.
However, as a result of obtaining the electromagnetic field in the equivalent element analysis region set by the user operating the operation unit 2, the average magnetic field calculation unit 4 a determines the magnetic flux density easy axis X direction on the line 22 a in the equivalent element 31. And the value of the magnetic field easy axis X direction are determined to be not constant, the equivalent element analysis region set by the user is changed.

  Specifically, the average magnetic field calculation unit 4a includes the number of magnetic bodies (steel plates) arranged in a bowl shape in the equivalent element analysis region, and the easy axis of magnetization of the magnetic bodies (steel plate) in the equivalent element analysis region. Based on the length in the X direction and the distance between the end of the magnetic element (steel plate) in the equivalent element analysis region and the end of the equivalent element analysis region, the magnetic element (steel plate) in the equivalent element analysis region is generated. A region having such a size that the influence of the demagnetizing field and the edge does not appear in the required electromagnetic field is reset as a new equivalent element analysis region.

  In particular, in the present embodiment, the average magnetic field calculation unit 4a resets a region as close as possible to the equivalent element analysis region set by the user as a new equivalent element analysis region. That is, the average magnetic field calculation unit 4a gradually increases the equivalent element analysis region set by the user stepwise until it determines that the obtained electromagnetic field is not affected by the demagnetizing field and the edges. The electromagnetic field when it is determined that the magnetic field and the edge are not affected is displayed on the display unit 3.

Next, an example of processing operation in the electromagnetic field analysis apparatus of the present embodiment will be described with reference to the flowchart of FIG.
The processing operations from step S1 to step S17 are the same as those shown in FIG. That is, a predetermined external magnetic field H _ext is given to the equivalent element 31, and an average magnetic flux density B _ave (i, j) and an average magnetic field H _ave (i, j) in the equivalent element 31 are obtained (steps S1 to S12). ).

Then, from the obtained average magnetic flux density B _ave (i, j) and average magnetic field H _ave (i, j), the magnetic flux density on the line 22 a parallel to the hard axis Y passing through the center 31 a of the equivalent element 31 is calculated. The value of the easy magnetization axis X direction and the value of the magnetic field easy magnetization axis X direction are extracted, and the value of the extracted magnetic flux density in the easy magnetization axis X direction and the value of the magnetic field easy magnetization axis X direction are constant. It is determined whether or not (step S13).

If the result of this determination is constant, the BH curves 90 to 93 in the equivalent element 31 are obtained from the obtained average magnetic flux density B _ave (i, j) and average magnetic field H _ave (i, j). Create (step S14). Further, B-θ curves 100 to 103 in the equivalent element 31 are created from the obtained average magnetic flux density B _ave (i, j) and the angle θ _BH formed by the average magnetic flux density B _ave and the average magnetic field H _ave (step). S15).
Then, using the created BH curves 90 to 93 and B-θ curves 100 to 103, an electromagnetic field generated in the analysis target region 21 is obtained, the obtained electromagnetic field is displayed on the display unit 3, and the process is terminated. (Steps S16 and S17).

  On the other hand, if at least one of the value of the magnetic flux density in the easy axis X direction and the value of the magnetic field in the easy axis X direction are not constant, the process proceeds to step S41 and the average magnetic field calculation unit 4a. Increases the size of the equivalent element analysis region 30 by a predetermined ratio and resets a new equivalent element analysis region 30. In step S13, the processing operations in steps S1 to S13 and S21 are repeated until the value of the magnetic flux density in the easy magnetization axis X direction and the value of the magnetic field in the easy magnetization axis X direction are constant. In this case, the process of step S3 is omitted.

  As described above, in the present embodiment, on the line 22a in the equivalent element 31, at least one of the value of the magnetic flux density in the easy axis X direction and the value of the magnetic field in the easy axis X direction is not constant. In this case, the equivalent element analysis region 30 is increased stepwise at a predetermined rate until these values become constant, and the electromagnetic field when these values become constant is displayed. In addition, an appropriate equivalent element analysis region 30 having no influence of the edges can be automatically set by the electromagnetic field analysis apparatus. Thereby, in addition to the effect in the first embodiment described above, the user does not have to reset the equivalent element analysis region 30 by operating the operation unit 2, and the electromagnetic field in the analysis target region 21 is analyzed. The burden on the user is reduced as much as possible.

  As in the first embodiment described above, in this embodiment, the process of step S13 for determining whether or not the magnetic characteristics in the equivalent element 31 are constant is performed after step S12. However, it is not necessarily performed at this timing. For example, the above processing may be performed after step S5. In this case, if the magnetic characteristics in the equivalent element 31 are constant, the process proceeds to step S6, and if not, the process may proceed to proceed to step S41.

(Third embodiment)
Next, a third embodiment of the present invention will be described. In the first and second embodiments described above, the case where a two-dimensional electromagnetic field is analyzed has been described. In the present embodiment, a case where a three-dimensional electromagnetic field is analyzed will be described. As described above, this embodiment is different from the first and second embodiments described above only in the dimension of the electromagnetic field to be analyzed, and thus is the same as the first and second embodiments described above. About the part, the detailed description is abbreviate | omitted by attaching | subjecting the code | symbol same as the code | symbol attached | subjected to FIGS.

  FIG. 21 is a diagram showing an example of the configuration of a magnetic shield device that performs analysis by the electromagnetic field analysis device of the present embodiment. In the present embodiment, an electromagnetic field analyzer for evaluating the shielding performance of this magnetic shield device will be described.

  In FIG. 21, the magnetic shield device 220 is configured by combining 32 magnetic bodies (steel plates) having a length of 300 mm, a width of 25 mm, and a thickness of 1 mm. More specifically, the magnetic shield device 220 is configured by stacking eight sets of four magnetic bodies (steel plates) arranged in a “mouth” shape at intervals of 30 mm. In this case, an equivalent element 222 and an equivalent element 223 can be considered as equivalent elements composed of a magnetic body and a non-magnetic body. In the present embodiment, an equivalent magnetic characteristic in the case of the equivalent element 222 will be described first. FIG. 21 corresponds to FIG. 2 in the first embodiment described above.

  FIG. 22 is a diagram showing an example of an analysis model for obtaining the magnetic characteristics of the equivalent element 222 shown in FIG. 21, and corresponds to FIG. 3 or FIG. 19 in the first embodiment described above. As shown in FIG. 22, an equivalent element 222 exists in the equivalent element analysis region 200, and the central part 221 a of the magnetic body 221 is in the equivalent element 222.

  In the present embodiment, the magnetic body (central portion) 221a in the equivalent element 222 has a width (a vertical axis (Z axis) direction perpendicular to the easy magnetization axis (X axis) and the hard magnetization axis (Y axis). ) Is 25 mm, the thickness (length in the hard axis (Y axis) direction) is 1 mm, and the length (length in the easy axis (X axis) direction) is 25 mm.

In addition, the equivalent element analysis region 200 in the present embodiment is a three-dimensional region (space) including the equivalent element 222 in the center.
Thus, in this embodiment, an electromagnetic field generated in a three-dimensional region (space) defined by the easy magnetization axis X, the hard magnetization axis Y, and the vertical axis Z of the magnetic body (steel plate) 221 is analyzed. I have to.

  Further, in the present embodiment, as shown in FIG. 23A, the equivalent element analysis region 200 is divided into a plurality of divided regions made of a cube or a rectangular parallelepiped. Each of the plurality of divided regions has an appropriate size that can be calculated by a finite element method (FEM).

When the external magnetic field H _ext is applied to the boundary region of the equivalent element analysis region 200, the average magnetic field calculation unit 4a performs a magnetic flux density B generated in each of a plurality of divided regions, a magnetic field H, and an electromagnetic field analysis by a finite element method. Ask for. The external magnetic field H _ext is expressed as the following (Equation 15).
H _ext = {H _extX , H _extY , H _extz } (Expression 15)

Thus, the external magnetic field H _ext is a three-dimensional vector having a direction of easy magnetization value H _EXTx, the values of H _ExtY the hard axis direction, and a value of the vertical axis value H _Extz .

When such an external magnetic field H _{ext is applied} to the equivalent element 200, the magnetic flux density B (i, j, k) generated in the divided area indicated by the oblique lines in FIG. 23A and the magnetic field H (i, j , K) are respectively expressed by the following (Expression 16) and (Expression 17).

B (i, j, k) = {B _X (i, j, k), B _Y (i, j, k), B _z (i, j, k)} (16)
H (i, j, k) = {H _X (i, j, k), H _Y (i, j, k), B _z (i, j, k)} (Expression 17)

In the above, i, j, and k are natural numbers for specifying the location of the divided region.
Thus, the magnetic flux density B (i, j, k) and magnetic field H (i, j, k) generated in each divided region are also values B _X (i, j, k), H _X ( i, j, k), values B _Y (i, j, k), H _Y (i, j, k) in the hard axis direction, and values B _z (i, j, k) in the vertical axis direction, A three-dimensional vector having H _z (i, j, k).

The average magnetic field calculation unit 4a generates the magnetic flux density generated in the equivalent element 201 from the magnetic flux density B (i, j, k) and the magnetic field H (i, j, k) in each divided region obtained as described above. B (i, j, k) and magnetic field H (i, j, k) are extracted. From the extracted magnetic flux density B (i, j, k) and magnetic field H (i, j, k), the equivalent element 222 The average magnetic flux density B _ave and the average magnetic field H _ave are obtained.

As shown in FIG. 23, in this embodiment, the average magnetic flux density B _ave is _obtained by dividing it into a component B _{aveX-Y in the} easy axis-hard magnetization axis direction and a component B _{aveZ in} the vertical axis direction. In the following description, the component B _{aveX-Y in the} easy axis-hard magnetization axis direction and the component B _{aveZ in} the vertical axis direction of the average magnetic flux density B _ave are respectively _expressed as the second and third average magnetic flux densities B _{aveX. -Y} and _BaveZ .

The average magnetic field H _ave is also _obtained by dividing into a component H _{aveX-Y in the} easy axis-hard magnetization axis direction and a component H _{aveZ in} the vertical axis direction. In the following description, the component H _{aveX-Y in the} easy axis-hard magnetization axis direction and the component H _{aveZ in} the vertical axis direction of the average magnetic field H _ave are respectively _expressed as the second and third average magnetic fields H _{aveX-Y. And HaveZ} .

As for the easy magnetization axis-hard magnetization axis direction, the second average magnetic flux density magnitude B _aveX-Y and the angle θ _B are set to the desired values as in the first embodiment. The second average magnetic flux density B _aveX-Y , the second average magnetic field H _aveX-Y , the second average magnetic flux density B _aveX-Y and the second are changed by changing the magnitude and direction of the external magnetic field H _ext1 . the average magnetic field H _avex-Y and obtains the angle theta _BHX-Y of (see Figure 23).
In the above, the external magnetic field H _ext1 represents a component of the external magnetic field H _{ext in} the direction of the easy axis to the hard axis.

On the other hand, in the vertical axis direction, the magnitude of the external magnetic field H _ext2 is changed to change the third average magnetic flux density B _aveZ , the third average magnetic field H _aveZ , the third average magnetic flux density B _aveZ, and the third average magnetic flux density B _aveZ . An angle θ _BHZ formed with the average magnetic field H _aveZ is obtained.
In the above, the external magnetic field H _ext2 represents the vertical axis direction component of H _ext .

Also in the present embodiment, similarly to the first and second embodiments described above, by making the equivalent element analysis region 200 sufficiently large, the electromagnetic field in the equivalent element 222 is changed to the magnetic body (steel plate) 20. The demagnetizing field generated inside and the influence of the edges are avoided as much as possible.
Specifically, with reference to FIG. 26, the average magnetic field calculation unit 4a is difficult to magnetize through the center 201a of the magnetic body (steel plate) 221a in the equivalent element 222 (the center of the equivalent element 222 in the example of FIG. 26). At a predetermined position on the surface 2401 in the equivalent element 222 parallel to the axis Y and the vertical axis Z, the value of the magnetic flux density in the easy magnetization axis X direction and the value of the magnetic field (magnetic field strength) in the easy magnetization axis X direction However, when each is constant, the equivalent element analysis region 200 is regarded as effective, and the calculated magnetic characteristics (average magnetic flux density, average magnetic field, etc.) in the equivalent element 222 are output to the magnetic characteristic creation section 4b.

  On the other hand, when at least one of the value of the magnetic flux density in the easy magnetization axis X direction and the value of the magnetic field (magnetic field strength) in the easy magnetization axis X direction are not constant, the user operates the operation unit 2. Then, the equivalent element analysis area 200 is reset. Then, the average magnetic field calculation unit 4a obtains again the electromagnetic field in the reset analysis target region 200 as described above. At this time, as in the second embodiment, the electromagnetic field analysis device (average magnetic field calculation unit 4a) may automatically reset the equivalent element analysis region 200.

The magnetic characteristic curve creation unit 4b performs the second and third average magnetic flux densities B _aveX-Y and B _aveZ obtained by the average magnetic field calculation unit 4a as described above, and the second and third average magnetic fields H _{aveX. Based on -Y} and _HaveZ , BH curves 2201 and 2202 as shown in FIG. 24 are created. Further, the second and third average magnetic flux density B _avex-Y, based on the angle theta _BHX-Y, to create a B-theta curve 2301 as shown in FIG. 25.

That is, the BH curves 2201a to 2201n obtained from the second average magnetic flux density _BaveX-Y and the second average magnetic field _HaveX-Y are created in the same manner as in the first embodiment described above. Further, a BH curve 2202 obtained from the third average magnetic flux density B _aveZ and the third average magnetic field H _aveZ is created.
Further, a second average flux density B _avex-Y, the second average flux density B _avex-Y and B-theta curve 2301a determined from the angle theta _BHX-Y of the second average magnetic field H _avex-Y ˜2301n are created in the same manner as in the first embodiment described above.

  Then, the magnetic field distribution calculation unit 4c obtains an electromagnetic field in the analysis target area using the BH curves 2201 and 2202 and the B-θ curve 2301 obtained as described above.

As described above, in the present embodiment, the average magnetic flux density B _ave and the average magnetic field H _ave having three-dimensional values are obtained separately for the easy magnetization axis-hard magnetization axis direction component and the vertical direction component. As for the easy axis to hard axis direction component, a BH curve 2201 and a B-θ curve 2301 are created in the same manner as in the first embodiment described above, and the vertical axis direction is separately B. Since the −H curve 2202 is created, the BH curve and the B−θ curve in the three-dimensional equivalent element 201 can be created without performing complicated calculations.

  Also in the present embodiment, as in the first and second embodiments described above, the center of the magnetic body (steel plate) 221a in the equivalent element 222 (the center of the equivalent element 222 in the example of FIG. 26) 201a. At a predetermined position on the surface 2401 in the equivalent element 222 that is parallel to the hard axis Y and the vertical axis Z, and the value of the magnetic flux density in the easy axis X direction and the easy axis of the magnetic field (magnetic field strength). The equivalent element analysis region 200 is effective only when the values in the X direction are constant. That is, since the size of the three-dimensional analysis target region is set in consideration of the number and size of the magnetic bodies (steel plates) included in the three-dimensional equivalent element analysis region 200, the magnetic body (steel plate) It is possible to prevent the demagnetizing field and the influence of the edges generated inside from appearing in the obtained electromagnetic field, and the electromagnetic field can be obtained accurately even in the three-dimensional analysis target region.

In the present embodiment, the third average magnetic flux density B _avez and the third average magnetic field H _avez are calculated, but they may be prepared in advance without calculating them. That is, a BH curve prepared in advance may be used for the component in the vertical axis direction. In this case, the magnetic permeability (relative magnetic permeability) may be constant. For example, the value of relative permeability may be set to 1000.

Further, when the average magnetic flux density B _ave and the average magnetic field H _ave are obtained as described above, the three axes in the equivalent element 201 (the easy axis X, the hard axis Y, the vertical axis Z), and the equivalent element analysis region 200 If the three axes orthogonal to each other (vertical, horizontal, and height axes) do not coincide with each other, the equivalent element 201 is subjected to coordinate transformation or rotational transformation, etc. It is preferable to match the three axes at 200.

  Further, in the present embodiment, the equivalent element analysis region 200 and the equivalent element 222 are three-dimensional regions (spaces) determined from the easy magnetization axis X, the hard magnetization axis Y, and the vertical axis Z. Similarly to the first and second embodiments, the equivalent element analysis region and the equivalent element may be a two-dimensional region (plane) determined from the easy magnetization axis X and the hard magnetization axis Y.

  Further, the equivalent element 223 shown in FIG. 21 may be processed similarly to the case of obtaining the magnetic characteristics of the equivalent element 222. In this case, a three-dimensionally expanded model shown in FIG. 19 may be used as an analysis model for obtaining the magnetic characteristics of the equivalent element, and the corner portion of the analysis model may be used as the equivalent element.

(Fourth embodiment)
Next, a fourth embodiment of the present invention will be described. Note that the present embodiment and the third embodiment described above differ only in the calculation method of the average magnetic flux density B _ave and the average magnetic field H _ave , and therefore the first to third embodiments described above. About the same part, detailed description is abbreviate | omitted by attaching | subjecting the code | symbol same as the code | symbol attached | subjected to FIGS.

  As shown in FIG. 27, the equivalent element analysis region 200 and the equivalent element 222 in the present embodiment are the same as those in the third embodiment. Then, the equivalent element analysis region 200 is divided into a plurality of divided regions as shown in FIG.

When the external magnetic field H _ext is applied to the boundary region of the equivalent element analysis region 200, the average magnetic field calculation unit 4a performs a magnetic flux density B generated in each of a plurality of divided regions, a magnetic field H, and an electromagnetic field analysis by a finite element method. Ask for. The external magnetic field H _ext is a three-dimensional vector as in the third embodiment described above (see (Equation 15)).

Further, when the external magnetic field H _{ext is applied} to the equivalent element 222, the magnetic flux density B (i, j, k) and magnetic field H (i, j, k) generated in the divided areas shown by the oblique lines in FIG. Is a three-dimensional vector as in the third embodiment described above (see (16) and (17)).

The average magnetic field calculation unit 4a uses the magnetic flux density B (i, j, k) generated in the equivalent element 222 from the magnetic flux density B (i, j, k) generated in each divided region and the magnetic field H (i, j, k). ) And the magnetic field H (i, j, k), and the average magnetic flux density B _ave in the equivalent element 201 is calculated from the extracted magnetic flux density B (i, j, k) and the magnetic field H (i, j, k). An average magnetic field H _ave is obtained. Specifically, it is obtained by the following (formula 18) to (formula 25).

In the above, B _X (i, j, k) is a value in the easy axis direction of the magnetic flux density B (i, j, k). B _Y (i, j, k) is a value of the magnetic flux density B (i, j, k) in the hard axis direction. B _Z (i, j, k) is a value in the vertical axis direction of the magnetic flux density B (i, j, k).

H _X (i, j, k) is a value in the easy axis direction of the magnetic field H (i, j, k). H _Y (i, j, k) is a value in the hard axis direction of the magnetic field H (i, j, k). H _Z (i, j, k) is a value in the vertical axis direction of the magnetic field H (i, j, k).
ΔV (i, j, k) is the size (volume) of the divided area.

And as shown to Fig.28 (a), in this Embodiment, the direction of average magnetic flux density _Bave is pinpointed by angle (alpha) _B , (beta) _B. Specifically, the angle α _B is an angle formed between the reference line OP and the hard magnetization axis Y. Here, the reference line OP is a vector determined from the value B _{aveX in the} easy axis direction of the average magnetic flux density B _ave and the value B _aveY in the hard axis direction. The angle β _B is an angle formed by the reference line OP and the average magnetic flux density B _ave .

Here, the angle α _B is expressed by the following (Equation 26).
α _B = 90−θ _B (Expression 26)
In the above, θ _B is an angle formed by the reference line OP and the easy axis X, and corresponds to the angle θ _B described in the first to third embodiments.

Further, as shown in FIG. 28B, in the present embodiment, the direction of the average magnetic field H _ave is specified by the angles α _H and β _H. Specifically, the angle α _H is an angle formed by the reference line OQ and the hard magnetization axis Y. Here, the reference line OQ is a vector determined from the value H _aveX of the average magnetic field H _{ave in} the easy axis direction and the value H _aveY in the hard axis direction. The angle β _H is an angle formed by the reference line OQ and the average magnetic flux density H _ave .

Here, the angle α _H is expressed by the following (Expression 27).
α _H = 90−θ _H (Expression 27)
In the above, θ _H is an angle formed by the reference line OQ and the easy axis X, and corresponds to the angle θ _H described in the first to third embodiments.

Then, the average magnetic field calculation unit 4a changes the magnitude and direction of the external magnetic field H _ext so that the average magnetic flux density magnitude B _ave and the angles α _B and β _B have desired values, and the above (18 Calculations according to (Expression) to (25) are performed to obtain the average magnetic flux density B _ave and the average magnetic field H _ave .

In the present embodiment, the direction of the external magnetic field H _ext is specified by the angles φ ₁ and φ ₂ as shown in FIG. Here, the angle φ ₁ is an angle between the easy axis X and the reference line OR determined from the value H _{extx in the} easy axis direction of the external magnetic field H _ext and the value H _extY in the hard axis direction. Further, the angle φ ₂ is an angle formed by the reference line OR and the external magnetic field H _ext . In the following description, the case where the direction of the external magnetic field H _ext is specified by using the angles φ ₁ and φ ₂ will be described as an example. However, if the direction of the external magnetic field H _ext can be specified, the angle φ ₁ , The direction of the external magnetic field H _ext need not be specified using φ ₂ .

Further, the average magnetic field calculation unit 4a obtains an angle α _BH formed by the reference line OP and the reference line OQ, as shown in FIG. The angle α _BH corresponds to the angle θ _BH described in the first to third embodiments.

As described above, the average magnetic field calculation unit 4a of the present embodiment obtains the following six relationships using the magnitude of the external magnetic field H _ext and the angles φ ₁ and φ ₂ as parameters, respectively.
1. External magnetic field H _{ext −average} magnetic flux density B _ave
2. External magnetic field H _{ext -average} magnetic field H _ave
3. External magnetic field H _{ext -angle} α _B
4). External magnetic field H _{ext -angle} β _B
5). External magnetic field H _{ext -angle} β _H
6). External magnetic field H _{ext -angle} α _H
In the above, the relationship between the external magnetic field H _ext and the angle β _H (external magnetic field H _{ext −angle} β _H ) or the relationship between the external magnetic field H _ext and the angle α _H (external magnetic field H _{ext −angle} α _H ). Instead, the relationship between the external magnetic field H _ext and the angle α _BH (external magnetic field H _{ext −angle} α _BH ) may be obtained.

  Also in the present embodiment, similarly to the above-described third embodiment, by making the equivalent element analysis region 200 sufficiently large, the electromagnetic field in the equivalent element 222 is inside the magnetic body (steel plate) 20. The demagnetizing field generated and the influence of the edges are avoided as much as possible.

  That is, as shown in FIG. 26, the average magnetic field calculation unit 4a, at a predetermined position on the surface 2401 in the equivalent element 222, the value of the magnetic flux density in the easy axis X direction and the magnetization of the magnetic field (magnetic field strength). When the values in the easy axis X direction are constant, the equivalent element analysis region 200 is regarded as effective, and the calculated magnetic characteristics (average magnetic flux density, average magnetic field, etc.) in the equivalent element 222 are calculated as a magnetic characteristic creation unit. Output to 4b.

  On the other hand, when at least one of these values is not constant, the equivalent element analysis region 200 is reset, and the electromagnetic field in the reset equivalent element analysis region 200 is obtained again as described above. At this time, as in the second embodiment described above, the electromagnetic field analysis device (average magnetic field calculation unit 4a) may automatically reset the equivalent element analysis region 200.

Then, as shown in FIG. 29, the magnetic characteristic curve creation unit 4b uses the angles α _B and β _B as parameters, and the BH curves 2701a to 2701n in the equivalent element 222, that is, the average magnetic flux density B _ave and the average magnetic field H. Create a curve representing the relationship with _ave .

Further, as shown in FIG. 30, the magnetic characteristic curve creation unit 4b uses the angles α _B and β _B as parameters, and the first B-θ curves 2801a to 2801n in the equivalent element 222, that is, the average magnetic flux density B _ave Then, a curve representing the relationship between the angle α _BH formed by the reference line OP and the reference line OQ is created.
Further, as shown in FIG. 31, the magnetic characteristic curve creation unit 4b uses the angles α _B and β _B as parameters, and the second B-θ curves 2901a to 2901n in the equivalent element 222, that is, the average magnetic flux density B _ave Then, a curve representing the relationship between the reference line OQ and the angle β _H formed by the average magnetic field H _ave is created.

  Then, the magnetic field distribution calculation unit 4c uses the BH curve 2701, the first B-θ curve 2801, and the second B-θ curve 2901 obtained as described above in the analysis target region. Find the electromagnetic field.

As described above, in the present embodiment, the relationship between the average magnetic flux density B _ave having a three-dimensional value and the average magnetic field H _ave , the relationship between the average magnetic flux density B _ave and the angle α _BH , and the average magnetic flux density B _ave And the angle β _H are obtained using the angles α _B and β _B as parameters, the average magnetic flux density B _ave and average magnetic field H _ave having three-dimensional values can be obtained more accurately.

  Also in the present embodiment, similarly to the third embodiment described above, the value of the magnetic flux density in the easy axis X direction and the magnetic field (magnetic field) at a predetermined position on the surface 2401 in the equivalent element 222. Since the equivalent element analysis region 200 is effective only when the value of the strength) in the easy magnetization axis X direction is constant, the demagnetizing field generated in the magnetic body (steel plate) and the influence of the edges are It is possible not to reveal the obtained electromagnetic field, and it is possible to accurately obtain the electromagnetic field in the three-dimensional analysis target region.

  In the present embodiment, the magnetic characteristic curve creation unit 4b creates the BH curves 2701a to 2701n and the first and second B-θ curves 2801a to 2801n and 2901a to 2901n. Needless to say, if the magnetic field distribution calculation unit 4c can obtain the electromagnetic field in the analysis target region, the curve created by the magnetic characteristic curve creation unit 4b is not limited to these.

(Other embodiments of the present invention)
The control operation by the electromagnetic field analysis device and the overall model analysis device in each of the above-described embodiments can be realized by using a computer system as shown in FIG.
FIG. 32 is a block diagram showing an example of the configuration of a computer system arranged in the electromagnetic field analysis apparatus 1.
32, a computer system 3100 includes a CPU 3101, a ROM 3102, a RAM 3103, a keyboard controller (KBC) 3105 of a keyboard (KB) 3104, a CRT controller (CRTC) 3107 of a CRT display (CRT) 3106 as a display unit, The disk controller (DKC) 3110 of the hard disk (HD) 3108 and flexible disk (FD) 3109 and the network interface controller (NIC) 3112 for connection to the network 3111 can communicate with each other via the system bus 3113. Connected configuration.

The CPU 3101 comprehensively controls each component connected to the system bus 3103 by executing software stored in the ROM 3102 or the HD 3108 or software supplied from the FD 3109.
That is, the CPU 3101 performs a control for realizing an operation to be described later by reading a processing program according to a predetermined processing sequence from the ROM 3102, the HD 3108, or the FD 3109 and executing it.

The RAM 3103 functions as a main memory or work area for the CPU 3101.
The KBC 3105 controls an instruction input from the KB 3104 or a pointing device (not shown).

A CRTC 3107 controls the display of the CRT 3106.
The DKC 3110 controls access to the HD 3108 and the FD 3109 that store a boot program, various applications, edit files, user files, a network management program, a predetermined processing program in the present embodiment, and the like.
The NIC 3112 exchanges data bidirectionally with devices or systems on the network 3111.

  In order to operate various devices so as to realize the functions of the above-described embodiments, the functions of the above-described embodiments can be realized for an apparatus connected to the various devices or a computer in the system. Implementations by supplying software program codes and operating the various devices in accordance with programs stored in a computer (CPU or MPU) of the system or apparatus are also included in the scope of the present invention.

  In this case, the program code itself of the software realizes the functions of the above-described embodiments, and the program code itself and means for supplying the program code to the computer, for example, the program code are stored. The recorded medium constitutes the present invention. As a recording medium for storing the program code, for example, a flexible disk, a hard disk, an optical disk, a magneto-optical disk, a CD-ROM, a magnetic tape, a nonvolatile memory card, a ROM, or the like can be used.

  Further, by executing the program code supplied by the computer, not only the functions of the above-described embodiments are realized, but also the OS (operating system) or other application software in which the program code is running on the computer. Such program code is also included in the embodiment of the present invention when the functions of the above-described embodiment are realized in cooperation with the above.

  Further, after the supplied program code is stored in the memory provided in the function expansion board of the computer or the function expansion unit connected to the computer, the CPU provided in the function expansion board or function expansion unit based on the instruction of the program code The present invention also includes a case where the functions of the above-described embodiment are realized by performing part or all of the actual processing.

It is the block diagram which showed the 1st Embodiment of this invention and showed an example of the structure of the electromagnetic field analyzer. It is the figure which showed the 1st Embodiment of this invention and showed an example of the analysis object area | region analyzed with an electromagnetic field analyzer. It is the figure which showed the 1st Embodiment of this invention and showed an example of the analysis model for calculating | requiring the magnetic characteristic of an equivalent element. It is the figure which showed the 1st Embodiment of this invention and expanded and showed an example of the equivalent element. It is the figure which showed the 1st Embodiment of this invention and showed an example of the division area obtained by dividing | segmenting an equivalent element analysis area. It is the figure which showed the 1st Embodiment of this invention and showed an example of the BH curve of a magnetic body (steel plate). It is the figure which showed the 1st Embodiment of this invention and showed an example of the curve showing the relationship between the angle which an average magnetic flux density and an easy magnetization axis make, and an external magnetic field. It is the figure which showed the 1st Embodiment of this invention and expanded and showed the curve showing the relationship between the angle which an average magnetic flux density and an easy magnetization axis make, and an external magnetic field. It is the figure which showed the 1st Embodiment of this invention and showed an example of the BH curve in an equivalent element. It is the figure which showed the 1st Embodiment of this invention and showed an example of the BH curve in the equivalent element which performed the leveling process. It is the figure which showed an example of the B-theta curve in an equivalent element. It is a flowchart which shows the 1st Embodiment of this invention and demonstrates an example of the processing operation in an electromagnetic field analyzer. It is a flowchart which shows the 1st Embodiment of this invention and demonstrates an example of the specific process operation at the time of creating the BH curve in an electromagnetic field analyzer. It is a flowchart which shows the 1st Embodiment of this invention and demonstrates an example of the specific process operation at the time of creating the B-theta curve in an electromagnetic field analyzer. It is the figure which showed the 1st Embodiment of this invention and showed collectively the method of calculating | requiring the electromagnetic field which arises in an analysis object area | region. 1 shows the first embodiment of the present invention, the value in the easy axis direction of the magnetic flux density obtained by the magnetic field distribution calculation unit, and the distance in the hard axis direction from the center of the magnetic body (steel plate) in the equivalent element It is the figure which showed an example of the relationship. The 1st Embodiment of this invention is shown, the value of the magnetization easy axis direction of the magnetic field (magnetic field strength) calculated | required by the magnetic field distribution calculating part, and the magnetization difficult axis from the center of the magnetic body (steel plate) in an equivalent element It is the figure which showed an example of the relationship with the distance in a direction. It is the figure which showed the 1st Embodiment of this invention and showed the other example of the analysis object area | region. It is the figure which showed the 1st Embodiment of this invention and showed the other example of the equivalent element. It is a flowchart which shows the 2nd Embodiment of this invention and demonstrates an example of the processing operation in an electromagnetic field analyzer. It is the figure which showed the 3rd Embodiment of this invention and showed an example of the structure of the magnetic shielding apparatus. It is the figure which showed the 3rd Embodiment of this invention and showed an example of an equivalent element analysis area | region and an equivalent element. It is the figure which showed the 3rd Embodiment of this invention and showed an example of the division area obtained by dividing | segmenting an equivalent element analysis area. It is the figure which showed the 3rd Embodiment of this invention and showed an example of the BH curve in an equivalent element. It is the figure which showed the 3rd Embodiment of this invention and showed an example of the B-theta curve in an equivalent element. It is the figure which showed the 3rd Embodiment of this invention and expanded and showed the equivalent element. It is the figure which showed the 4th Embodiment of this invention and showed an example of an equivalent element analysis area | region and an equivalent element. It is a figure which shows the 4th Embodiment of this invention and demonstrates an example of the method of specifying the direction of an average magnetic flux density and the direction of an average magnetic field. It is the figure which showed the 4th Embodiment of this invention and showed an example of the BH curve in an equivalent element. It is the figure which showed the 4th Embodiment of this invention and showed an example of the 1st B-theta curve in an equivalent element. It is the figure which showed the 4th Embodiment of this invention and showed an example of the 2nd B-theta curve in an equivalent element. It is the block diagram which showed other embodiment of this invention and showed an example of the structure of the computer system arrange | positioned at the electromagnetic field analyzer.

Explanation of symbols

DESCRIPTION OF SYMBOLS 1 Electromagnetic field analyzer 2 Operation part 3 Display part 4 Processing part 4a Average magnetic field calculating part 4b Magnetic characteristic curve preparation part 4c Magnetic field distribution calculating part 20, 221 Magnetic body (steel plate)
21 Analysis target region 30, 200 Equivalent element analysis region 31, 222, 223 Equivalent element 90-93, 2201, 2202, 2701 BH curve 100-103, 2301 B-θ curve 220 Magnetic shield device 2801 First B- θ curve 2901 Second B-θ curve

Claims (24)

An average magnetic field calculating means for calculating an average magnetic flux density and an average magnetic field in an equivalent element occupied by a plurality of substances including a magnetic material;
Electromagnetic field analysis means for analyzing an electromagnetic field generated in an analysis target area wider than the equivalent element using the average magnetic flux density and the average magnetic field,
The average magnetic field calculation means sets the value of the magnetic characteristic in the easy axis direction of the magnetic property at a predetermined position in the equivalent element to be constant, and in that state, calculates the average magnetic flux density and the average magnetic field in the equivalent element. An electromagnetic field analysis device characterized by calculating. The average magnetic field calculation means calculates an average magnetic flux density and an average magnetic field in an equivalent element occupied by a plurality of substances including a magnetic material, and calculates an angle formed by the average magnetic flux density and the average magnetic field. And
The electromagnetic field analysis means analyzes an electromagnetic field generated in an analysis target area wider than the equivalent element using the average magnetic flux density, the average magnetic field, and the angle formed by the average magnetic flux density and the average magnetic field. The electromagnetic field analysis apparatus according to claim 1.   The average magnetic field calculation means calculates a magnetic characteristic in an equivalent element analysis region wider than an equivalent element occupied by a plurality of substances including a magnetic substance, and uses the calculated magnetic characteristic to calculate an average in the equivalent element. The electromagnetic field analysis apparatus according to claim 1, wherein a magnetic flux density and an average magnetic field are obtained.   Based on the size of the equivalent element analysis area, the average magnetic field calculation means makes the value of the magnetic characteristic in the easy axis direction of the magnetic property at a predetermined position in the equivalent element constant, and in that state, The electromagnetic field analysis apparatus according to claim 3, wherein an average magnetic flux density and an average magnetic field in the element are obtained.   The average magnetic field calculation means includes the number of magnetic bodies in the equivalent element analysis region, the length in the easy axis direction of the magnetic body in the equivalent element analysis region, and the magnetic body in the equivalent element analysis region. Based on at least one of the distances from the end of the equivalent element analysis region to the end of the equivalent element analysis region, the value of the magnetic property in the easy magnetization direction in a predetermined position in the equivalent element is made constant, 5. The electromagnetic field analysis apparatus according to claim 3, wherein an average magnetic flux density and an average magnetic field in the equivalent element are obtained in that state.   The electromagnetic field analysis apparatus according to claim 1, wherein the electromagnetic field analysis unit analyzes an electromagnetic field generated in a region where a plurality of substances including the magnetic substance are repeatedly present.   The average magnetic field calculation means makes the value of the easy magnetization direction of the magnetic property at a predetermined position on a line parallel to the hard magnetization axis of the magnetic body located in the equivalent element constant, and in that state, The electromagnetic field analysis apparatus according to claim 6, wherein an average magnetic flux density and an average magnetic field in the equivalent element are obtained.   The average magnetic field calculation means is configured to provide a magnetic field at a predetermined position on a plane parallel to the hard axis of magnetization of the magnetic body located within the equivalent element and the easy axis of the magnetic body and the vertical axis perpendicular to the hard axis of magnetization. 7. The electromagnetic field analysis apparatus according to claim 6, wherein a value of the characteristic easy axis direction is made constant and an average magnetic flux density and an average magnetic field in the equivalent element are obtained in that state. The plurality of substances are a magnetic material and a non-magnetic material,
The average magnetic field calculating means calculates an average magnetic flux density and an average magnetic field in an equivalent element occupied by a partial region of the magnetic material and a nonmagnetic material around the magnetic material. Item 9. The electromagnetic field analysis device according to any one of Items 1 to 8.   The electromagnetic field analysis apparatus according to claim 1, wherein the magnetic characteristics include a magnetic flux density and a magnetic field. A BH curve representing the relationship between the average magnetic flux density calculated by the average magnetic field calculation means and the average magnetic field, the average magnetic flux density, and the relationship between the average magnetic flux density and the angle formed by the average magnetic field. A magnetic characteristic curve creating means for creating a -θ curve;
The electromagnetic field analysis means analyzes an electromagnetic field generated in an analysis target area wider than the equivalent element based on the BH curve and the B-θ curve created by the magnetic characteristic curve creation means. The electromagnetic field analysis apparatus according to claim 10. An average magnetic field calculation step for calculating an average magnetic flux density in an equivalent element occupied by a plurality of substances including a magnetic material and an average magnetic field;
Using the average magnetic flux density and the average magnetic field, an electromagnetic field analysis step for analyzing an electromagnetic field generated in an analysis target area wider than the equivalent element,
In the average magnetic field calculation step, the value of the easy magnetization direction of the magnetic characteristics at a predetermined position in the equivalent element is made constant, and in that state, the average magnetic flux density and the average magnetic field in the equivalent element are calculated. An electromagnetic field analysis method characterized by calculating. The average magnetic field calculation step calculates an average magnetic flux density and an average magnetic field in an equivalent element occupied by a plurality of substances including a magnetic material, and calculates an angle formed by the average magnetic flux density and the average magnetic field. And
The electromagnetic field analysis step analyzes an electromagnetic field generated in an analysis target area wider than the equivalent element by using the average magnetic flux density, the average magnetic field, and the angle formed by the average magnetic flux density and the average magnetic field. The electromagnetic field analysis method according to claim 12.   The average magnetic field calculation step calculates a magnetic characteristic in an equivalent element analysis region wider than an equivalent element occupied by a plurality of substances including a magnetic substance, and uses the calculated magnetic characteristic to calculate an average in the equivalent element. The electromagnetic field analysis method according to claim 12 or 13, wherein a magnetic flux density and an average magnetic field are obtained.   In the average magnetic field calculation step, based on the size of the equivalent element analysis region, the value of the magnetic characteristic in the easy axis direction of the magnetic characteristic at a predetermined position in the equivalent element is made constant, and the equivalent The electromagnetic field analysis method according to claim 14, wherein an average magnetic flux density in the element and an average magnetic field are obtained.   The average magnetic field calculation step includes the number of magnetic bodies in the equivalent element analysis region, the length of the magnetic body in the equivalent element analysis region in the easy axis direction, and the magnetic body in the equivalent element analysis region. Based on at least one of the distances from the end of the equivalent element analysis region to the end of the equivalent element analysis region, the value of the magnetic property in the easy magnetization direction in a predetermined position in the equivalent element is made constant, The electromagnetic field analysis method according to claim 14 or 15, wherein an average magnetic flux density and an average magnetic field in the equivalent element are obtained in that state.   The electromagnetic field analysis method according to any one of claims 12 to 16, wherein the electromagnetic field analysis step analyzes an electromagnetic field generated in a region where a plurality of substances including the magnetic substance are repeatedly present.   In the average magnetic field calculation step, the value of the magnetic property in the easy axis direction of the magnetic property at a predetermined position on a line parallel to the hard axis of magnetization of the magnetic substance located in the equivalent element is made constant, and in that state, The electromagnetic field analysis method according to claim 17, wherein an average magnetic flux density and an average magnetic field in the equivalent element are obtained.   The step of calculating the average magnetic field includes the magnetic field at a predetermined position on a plane parallel to the hard axis of magnetization of the magnetic body located within the equivalent element and the easy axis of the magnetic body and the vertical axis perpendicular to the hard axis of magnetization. 18. The electromagnetic field analysis method according to claim 17, wherein the value of the characteristic easy axis direction is made constant, and the average magnetic flux density and the average magnetic field in the equivalent element are obtained in that state. The plurality of substances are a magnetic material and a non-magnetic material,
The average magnetic field calculating step calculates an average magnetic flux density and an average magnetic field in an equivalent element occupied by a partial region of the magnetic material and a nonmagnetic material surrounding the region. Item 20. The electromagnetic field analysis method according to any one of Items 12 to 19.   The electromagnetic field analysis method according to claim 12, wherein the magnetic characteristics include a magnetic flux density and a magnetic field. A BH curve representing the relationship between the average magnetic flux density calculated in the average magnetic field calculating step and the average magnetic field, the average magnetic flux density, and the relationship between the average magnetic flux density and the angle formed by the average magnetic field. A magnetic characteristic curve creating step for creating a -θ curve,
In the electromagnetic field analysis step, an electromagnetic field generated in an analysis target area wider than the equivalent element is analyzed based on the BH curve and the B-θ curve created in the magnetic characteristic curve creation step. The electromagnetic field analysis method according to claim 21. An average magnetic field calculation step for calculating an average magnetic flux density in an equivalent element occupied by a plurality of substances including a magnetic material and an average magnetic field;
Using the average magnetic flux density and the average magnetic field, the computer executes an electromagnetic field analysis step for analyzing an electromagnetic field generated in an analysis target area wider than the equivalent element,
In the average magnetic field calculation step, the value of the easy magnetization direction of the magnetic characteristics at a predetermined position in the equivalent element is made constant, and in that state, the average magnetic flux density and the average magnetic field in the equivalent element are calculated. A computer program characterized by computing.   A computer-readable recording medium on which the computer program according to claim 23 is recorded.

Images