JP5236599B2 - Magnetic field analyzer and method of operating magnetic field analyzer - Google Patents

Description

  The present invention relates to a magnetic field analysis apparatus and a magnetic field analysis method.

  2. Description of the Related Art In recent years, electrical devices such as motors that use magnetic energy as a drive source have been required to have high precision and high efficiency, and at the same time, have been required to increase development speed.

  In order to satisfy such requirements, verification by magnetic field analysis of electrical equipment is indispensable.

  Conventionally, the finite element method has been widely used as a magnetic field analysis technique.

  For example, Patent Document 1 describes an analysis method using a finite element method as an example of magnetic field analysis when a magnetic core is present in air.

  In the magnetic field analysis, it is necessary to use a magnetization curve in order to obtain the magnetization of the magnetic material. Conventionally, as described in Patent Document 2, a curve ignoring hysteresis (initial magnetization curve) has been used.

Japanese Patent Laid-Open No. 10-188042 Japanese Patent Laid-Open No. 05-54161

  However, when magnetic field analysis is performed by the finite element method as in Patent Document 1, it is necessary to mesh-divide the entire space including the target electric device.

  Therefore, when analyzing large-scale, complex-shaped electrical equipment or dealing with high aspect ratio problems such as laminated iron cores, the analysis requires a tremendous amount of time and memory of the analytical equipment, and performs magnetic field analysis efficiently with sufficient accuracy. There was a problem that it was difficult.

  In addition, as in Patent Document 2, in the analysis using only the initial magnetization curve, the magnetization is obtained without considering the magnetic hysteresis inherent in the magnetic material, so that the behavior of the magnetic material in the actual machine cannot be accurately grasped. there were.

  The present invention has been made in view of such problems, and an object thereof is to provide a magnetic field analysis apparatus capable of performing magnetic field analysis with sufficient accuracy and efficiency for a system having a magnetic material. .

  In order to achieve the above-mentioned object, the present inventor, as a result of intensive studies, uses an information processing device to calculate the solution of the motion equation of the magnetic field of the particles constituting the magnetic material, and based on the solution of the motion equation of the magnetic field. It has been found that by calculating the magnetic field at any place in the space where the magnetic material is arranged, the magnetic field analysis can be efficiently performed with sufficient accuracy for the system having the magnetic material, and the present invention has been achieved. .

  In addition, the inventors have found that by using the Potter-Schmulian equation as the magnetization curve, it is possible to perform a magnetic field analysis in consideration of the magnetic hysteresis inherent in the magnetic material, and have led to the present invention.

  That is, the first invention is a magnetic field analysis device for analyzing a magnetic field at an arbitrary point in a space in which a magnetic material is arranged by using an information processing device, and represents the shape of the magnetic material and the magnetization curve of the magnetic material. From the analysis condition having the equation, the storage means for storing the motion equation of the magnetic field, and the analysis condition, the coefficient of the motion equation of the magnetic field of the particles constituting the magnetic body, the initial value of the variable and the constant are calculated, A calculation means for calculating a magnetic field motion equation and calculating a magnetic field at an arbitrary point in the space based on the solution, and an equation representing a magnetization curve of the magnetic material is a Potter-Schmulian equation It is a magnetic field analyzer characterized by being.

  The second invention is a program for causing a computer to function as the magnetic field analysis apparatus described in the first invention.

  A third invention is a magnetic field analysis method for analyzing a magnetic field at an arbitrary point in a space in which a magnetic material is arranged by using an information processing device, and includes an equation representing a shape of the magnetic material and a magnetization curve of the magnetic material. (A) storing the analysis conditions and the equation of motion of the magnetic field, and calculating the coefficients and the initial values and constants of the equations of motion of the magnetic field of the particles constituting the magnetic material from the analysis conditions, (B) calculating a solution of the equation of motion, and calculating a magnetic field at an arbitrary point in the space based on the solution of the equation of motion of the magnetic field, and an equation representing the magnetization curve of the magnetic substance is Potter -A magnetic field analysis method characterized by the Schmulian equation.

  According to the present invention, it is not necessary to divide the entire space including the analysis target into a mesh for a system having a magnetic material, and to perform magnetic field analysis efficiently with sufficient accuracy in consideration of the magnetic hysteresis inherent in the magnetic material. A magnetic field analysis apparatus that can be performed can be provided.

  In addition, since the solution of the equation of motion of the magnetic field derived from the Lagrangian of the particles constituting the magnetic material is calculated, the total energy of the system is preserved when the input does not change.

2 is a diagram illustrating a hardware configuration of the magnetic field analysis apparatus 1. FIG. 3 is a diagram showing a storage device 5. FIG. 2 is a diagram showing an outline of the configuration of an SPM motor 31. FIG. FIG. 4 is an enlarged perspective view around one of the stator teeth 43 of FIG. 3. 4 is a flowchart showing a procedure for analyzing a magnetic field using the magnetic field analysis apparatus 1. FIG. 6A is an enlarged perspective view of the coil 45, and FIG. 6B is a diagram showing an example in which the coil 45 is divided into local conductors. It is a figure which shows the relationship between the rectangular parallelepiped conductor 45a and a local coordinate. It is a figure which shows the relationship between the arc-shaped columnar conductor 45c and a local coordinate. It is explanatory drawing of the element used for the magnetic field calculation of a magnetic body. It is a flowchart which shows the specific calculation method of magnetization (magnetization vector) of a magnetic body.

  DESCRIPTION OF EMBODIMENTS Hereinafter, preferred embodiments of the present invention will be described in detail based on the drawings.

  First, the hardware configuration of the magnetic field analysis apparatus 1 according to the present embodiment will be described with reference to FIGS. 1 and 2.

  It should be noted that the hardware configuration in FIG. 1 is an example, and is not limited to this.

  As shown in FIG. 1, the magnetic field analysis apparatus 1 includes a control unit 3, a storage device 5, a media input / output unit 6, an input unit 7, a display unit 9, a printer port 11, and the like connected to each other via a bus 13.

  The control unit 3 includes a CPU (Central Processing Unit), a ROM (Read Only Memory), a RAM (Random Access Memory), and the like, and is connected via the bus 13 according to a program stored in the storage device 5 as a storage unit. Each of the devices is driven and controlled.

  As shown in FIG. 2, the storage device 5 stores a control program 15 for driving and controlling each component of the magnetic field analysis device 1, and a magnetic field analysis program 17 for carrying out the present invention.

  The magnetic field analysis program 17 includes input information 21 that is information having analysis conditions, and an arithmetic program 19 that calculates a magnetic field based on the motion equation of the magnetic field based on the input information 21.

  The media input / output unit 6 is a device for inputting / outputting information to / from a medium such as a floppy (registered trademark) disk, CD, or DVD.

  The input unit 7 is an input device such as a keyboard and a mouse, and the display unit 9 is a display device such as a display.

  A printer 12 as an output device is connected to the printer port 11.

  Next, the procedure of magnetic field analysis using the magnetic field analyzer 1 will be described with reference to FIGS.

  Here, a magnetic field analysis of an SPM motor 31 (Surface Permanent Magnet Motor), which is a kind of PM motor (Permanent Magnet Motor), will be described as an example.

  First, an outline of the configuration of the SPM motor 31 will be described with reference to FIGS. 3 and 4.

  As shown in FIG. 3, the SPM motor 31 has a rotor 33 as a rotor (moving element) and a stator 35 as a stator.

  The rotor 33 has a cylindrical rotor core 37 that is a magnetic body such as iron, and a permanent magnet 39 is provided on the surface of the rotor core 37.

  A rod-shaped rotor shaft 41 is provided at the axial center of the rotor core 37.

  The stator 35 includes a tooth-like stator teeth 43 that is a magnetic body, a core back 44 that is a cylindrical magnetic body provided outside the stator teeth 43, and a cylindrical frame 46 that is provided outside the core back 44. Has been.

  As shown in FIGS. 3 and 4, a coil 45 that is a conductor such as metal is wound around the stator teeth 43.

  In the actual SPM motor 31, the coil 45 is wound around the stator teeth 43 with the number of turns according to the specifications to form a bundle, but in this embodiment, as illustrated in FIGS. 3 and 4. The coil bundle is handled as one conductor without modeling each coil.

  In the SPM motor 31 having such a structure, the rotor 33 and the stator 35 are magnetized by the magnetic field of the permanent magnet 39 and the magnetic field generated by passing a current through the coil 45. The SPM motor 31 is driven by the deviation of the magnetic energy of the magnetic material.

  Therefore, in order to perform the magnetic field analysis of the SPM motor 31, it is necessary to calculate the magnetic field vector that the coil 45 and the permanent magnet 39 form on the magnetic body constituting the rotor 33 and the stator 35, and to analyze the magnetization phenomenon of these magnetic bodies. is there.

  Next, the analysis procedure will be described with reference to FIGS.

  In the following procedure, the SPM motor 31 is divided into a plurality of elements to form an aggregate of particles for each element, and is handled as a rigid model. However, the present invention is not limited to this, and particles are used. Instead, it may be handled as a rigid model for each element.

  In the following procedure, the magnetic material indicates the magnetic material constituting the rotor 33 and the stator 35 and the permanent magnet 39, and these are distinguished by a function representing a magnetization curve.

  Further, as is apparent from FIG. 5, in this embodiment, the analysis is performed by repeatedly executing a part of the steps of FIG. 5 until the end condition is satisfied. The step immediately before a certain step is referred to as “the most recent step”, and the step immediately before the most recent step is referred to as “one step before”.

  First, the control unit 3 of the magnetic field analysis apparatus 1 activates the magnetic field analysis program 17 and functions representing the three-dimensional structure (shape, coordinate point), mass density, and magnetization curve of the magnetic material as analysis conditions of the SPM motor 31 to be analyzed. The current density vector of the conductor is stored as input information 21 of the storage device 5 (step 101 in FIG. 5).

  These physical quantities may be read from a recording medium such as a CD-ROM via the media input / output unit 6, for example.

  The information of the three-dimensional structure of the SPM motor 31 is data such as CAD.

  Further, when the physical quantity is stored in advance as the input information 21, the above step is not necessary.

  The current density vector is necessary when calculating the magnetic field vector formed by the coil 45.

  Furthermore, the function representing the magnetization curve of the magnetic material is necessary when calculating the magnetization vector (in bold) M, and in this embodiment, the Potter-Schmulian equation is used as will be described later.

  The above is the details of step 101 in FIG.

  Next, the control unit 3 of the magnetic field analysis apparatus 1 determines that the magnetic material is N particles (a polyhedral element whose center is the position vector of the particle (cubic element in this embodiment) from the information of the three-dimensional structure included in the input information 21. ), And the particle position vector is calculated and stored in the storage device 5 (step 102 in FIG. 5).

  Here, N is an arbitrary integer, and the number N of particles and the position vector are calculated based on the three-dimensional structure of the input information 21 and the shape of the polyhedral element stored in the storage device 5 in advance.

  The above is the details of step 102 in FIG.

  Next, the control unit 3 determines that the conductor (coil 45 shown in FIG. 6A) is a local conductor (rectangular conductors 45a and 45b) as shown in FIG. And arc-shaped columnar conductors 45c and 45d), the coefficient for calculating the magnetic field vector created by each conductor is calculated and stored in the storage device 5 (step 103 in FIG. 5).

  Here, step 103 in FIG. 5 will be specifically described.

  First, the case where the local conductor is a rectangular parallelepiped conductor 45a will be described.

  In addition, since it is the same as that of the case where a local conductor is the rectangular parallelepiped conductor 45a when a local conductor is the rectangular parallelepiped conductor 45b, description is abbreviate | omitted.

When the local conductor is a rectangular parallelepiped conductor 45a, the control unit 3 applies a local coordinate system (x _s , y _s , z _s ) as shown in FIG.

What this in the local coordinate system to the center of gravity of the local conductor (rectangular conductor 45a) and the origin _{O s,} dimensions of rectangular conductor 45a is having a length of _{x s} direction 2a, _{y s} direction 2b, _{z s} directions 2c And

And shall particles is located at the origin O _s.

  The control unit 3 reads the three-dimensional structure of the conductor and the current density vector stored in the storage device 5 as the input information 21 from step 101, and a, b, c that are the dimensions of the local conductor (cuboid conductor 45a) and the particles The position vector is calculated and stored in the storage device 5.

  Next, a case where the local conductor is an arc-shaped columnar conductor 45c will be described.

  Note that the case where the local conductor is the arc-shaped columnar conductor 45d is the same as the case where the local conductor is the arc-shaped columnar conductor 45c, and a description thereof will be omitted.

When the local conductor is the arc-shaped columnar conductor 45c, the control unit 3 applies the local coordinate system (x _c , y _c , z _c ) as shown in FIG.

Origin O _c are present on an arc of the central axis, and that the arcuate columnar conductor 45c are symmetrical with respect to the height direction of the arcuate columnar conductor 45c (z _c direction in FIG. 8) in this local coordinate system Shall exist.

_{Moreover, x c, y c, z} c _, when viewed in x c -y _c plane, + _{x c-axis} as a base point, and so the arc of the arcuate columnar conductor 45c is + when viewed from the _{z c-axis} counterclockwise decide.

The average value of the inner diameter and the outer diameter of the arcuate columnar conductor 45c and _{R c, 2r} the thickness in the radial direction a, the _{z c} direction height and 2z _b.

Current as a base point a + x _c-axis, + when viewed from the z _c the angle a θ of the counter-clockwise direction, current is assumed to flow in a uniform current density j in this direction.

It is assumed that the particles are located at (R _c , θ / 2, 0) in the cylindrical coordinate system.

The control unit 3 reads the three-dimensional structure of the conductor and the current density vector stored in the storage device 5 as the input information 21 from step 101, and the dimensions r _a and z _{b of the} local conductor (arc-shaped columnar conductor 45c). , Θ, R _c and the particle position vector are calculated and stored in the storage device 5.

  Further, when the physical quantity is stored in advance as the input information 21, the above step is not necessary.

  The above is the details of step 103 in FIG.

  Next, the control unit 3 of the magnetic field analysis device 1 stores the position vector of N particles constituting the magnetic body calculated up to this step and stored in the storage device 5 and stored in the storage device 5 in advance. The coefficient of the equation of motion of the magnetic field of the particles constituting the magnetic body is calculated according to the shape of the polyhedral element and stored in the storage device 5 (step 104 in FIG. 5).

  Here, step 104 in FIG. 5 will be specifically described.

  When the divided polyhedral elements (cubic elements in this embodiment) are displayed two-dimensionally, the shape shown in FIG. 9 is obtained.

Here, r _{g is} the particle position vector (bold), q is the midpoint of the element boundary surface of the polyhedral element centered on the particle position vector, and p is the midpoint between the particle position vector and q point. Define points.

  The control unit 3 reads the position vector of the particles (polyhedral elements (cubic elements in the present embodiment) having the center of gravity of the particles) constituting the magnetic body calculated in step 102 and stored in the storage device 5.

  If the particle position vector is updated and stored in the storage device 5 in step 113, the control unit 3 reads the value.

  From the particle position vector and the shape of the polyhedron element stored in the storage device 5 in advance, the control unit 3 calculates the points p and q in the polyhedron element whose center of gravity is the position vector of all particles constituting the magnetic body. The position vector and the normal vector (bold) n to the element boundary surface are calculated and stored in the storage device 5.

  Next, the control unit 3 includes a polyhedral element having the boundary area ΔS of the polyhedral element having the particle constituting the magnetic body as the center of gravity, and the vertex of the boundary surface of the polyhedral element having the particle position vector as the center of gravity and the particle position as the vertex. The volume ΔV is calculated and stored in the storage device 5.

  The above is the details of step 104 in FIG.

  Next, the control unit 3 of the magnetic field analysis apparatus 1 uses the size of the coil stored in the storage device 5 in step 103 and the current density vector stored in the storage device 5 up to this step to perform bio-savart. Based on the analytical solution obtained by integrating the above law, the energized coil 45 calculates the magnetic field vector formed on the particles constituting the magnetic material and stores it in the storage device 5 (step 105 in FIG. 5).

  Here, step 105 in FIG. 5 will be specifically described.

  Here, when it is assumed that an arbitrary position vector is a point P where the stator teeth 43 as a magnetic body in FIG. 4 are present, a procedure for calculating a magnetic field vector generated by the energized coil 45 on the position vector of the P point. Will be described as an example.

  First, the case where the local conductor is a rectangular parallelepiped conductor 45a will be described.

  In addition, since it is the same as that of the case where a local conductor is the rectangular parallelepiped conductor 45a when a local conductor is the rectangular parallelepiped conductor 45b, description is abbreviate | omitted.

When the local conductor is a rectangular parallelepiped conductor 45a, the local coordinate system (x _s , y _s , z _s ) is applied as shown in FIG.

  This local coordinate system is described in step 103.

Next, the position vector of the point P is converted into the local coordinate system (x _s , y _s , z _s ).

  Magnetic field vectors created at the point P by the energized rectangular parallelepiped conductor are described by the following equations (1) to (3).

Here, r _ps (bold) is the position vector of the point P in the local coordinate system (x _s , y _s , z _s ), and x _ps , y _ps , z _ps are x _s , y _s , z _s. Direction value.
π is the circumference ratio.

H _xs , H _ys , and H _zs are each component of the magnetic field vector in the local coordinate system (x _s , y _s , z _s ).
j is the current density.

X _i , y _j , and z _k represent the upper and lower limits of integration in the x _s , y _s , and z _s directions, and the relationship shown in Expression (4) is established.

  Here, a, b, and c are dimensions of a rectangular parallelepiped conductor.

  Next, a case where the local conductor is an arc-shaped columnar conductor 45c will be described.

  Note that the case where the local conductor is the arc-shaped columnar conductor 45d is the same as the case where the local conductor is the arc-shaped columnar conductor 45c, and a description thereof will be omitted.

When the local conductor is the arc-shaped columnar conductor 45c, the local coordinate system (x _c , y _c , z _c ) is applied as shown in FIG.

  This local coordinate system is described in step 103.

The position vector of the P point after being converted to the local coordinate system (x _c , y _c , z _c ) is expressed by the cylindrical coordinate system (bold) r _pc = (R _pc , φ _pc , Z _pc ).

  A magnetic field vector formed at the point P by the energized circular columnar conductor 45c is expressed by the following equations (6) to (11).

Here, H _rc , H _tc and H _zc are each component of the magnetic field vector in the cylindrical coordinate system.
j is the current density.

r _a , θ, and z _b are dimensions of the arc-shaped columnar conductor 45c, and R _c is an average value of the inner diameter and the outer diameter of the arc.

sgn is the sign of the _{Z k,} R _{j, Z k} is the upper limit of the integration represents the lower limit, the relationship shown in equation (12) holds.

  The above is the calculation procedure of the magnetic field vector created by the coil at an arbitrary point P.

  The control unit 3 of the magnetic field analysis apparatus 1 reads the position vector of the p point inside the polyhedron element having the center of gravity as the position vector of all particles constituting the magnetic body already calculated in step 104 and stored in the storage device 5. .

Next, the control unit 3 converts the position vector of the point p into the local coordinate system (x _s , y _s , z _s ), and calculates the dimensions of the rectangular parallelepiped conductor 45 a calculated in step 103 and stored in the storage device 5. From step 101, the current density vector stored in the storage device 5 as the input information 21 is read, and the magnetic field vector is calculated based on the equations (1) to (4).

  When the current density vector is updated in step 116 and stored in the storage device 5, the control unit 3 reads the value.

  Next, the control unit 3 converts the magnetic field vector calculated by the equations (1) to (4) into the global coordinate system (x, y, z), and stores the converted magnetic field vector in the storage device 5.

  Thereby, the magnetic field vector which the energized rectangular parallelepiped conductor 45a produces in p point is calculated | required.

Further, the control unit 3 converts the position vector of the point p into the local coordinate system (x _c , y _c , z _c ), and further, the cylindrical coordinate system (bold) r _pc = (R _pc , φ _pc , Z _pc ). Convert to

Next, the control unit 3 reads the dimensions r _a , θ, z _b and R _c of the arc-shaped columnar conductor 45 _c calculated in step 103 and stored in the storage device 5, and the control unit 3 inputs in step 101. A current density vector stored in the storage device 5 as information 21 is read.

  When the current density vector is updated in step 116 and stored in the storage device 5, the control unit 3 reads the value.

  The controller 3 calculates a magnetic field vector based on the equations (6) to (12).

  Next, the calculated magnetic field vector is converted into an orthogonal coordinate system, further converted into a global coordinate system (x, y, z), and stored in the storage device 5.

  As a result, a magnetic field vector created at the point p by the energized circular columnar conductor 45c is obtained.

  The above is the details of step 105 in FIG.

  Next, the control unit 3 of the magnetic field analysis apparatus 1 uses the calculation program 19 to calculate the magnetic field vector after the virtual time step δt without considering the constraint from the motion equation of the magnetic field of the particles constituting the magnetic body, and stores it. It memorize | stores in the apparatus 5 (step 106 of FIG. 5).

  Here, step 106 in FIG. 5 will be specifically described.

  The Lagrangian of the N particles constituting the magnetic material is assumed to have a shape represented by the equations (13) to (15).

Here, in Expression (13), α is a virtual mass, bold r is a position vector, bold H is a magnetic field vector, bold bold H is a time derivative of the magnetic field vector, bold M is a magnetization vector, and bold n Is the normal vector of the polyhedral element boundary surface with the particle position vector as the center of gravity, bold Hex is the externally applied magnetic field vector, ΔS is the area of the element boundary surface of the polyhedral element with the particle position vector as the center of gravity, and ΔV is The vertex of the polyhedral element boundary surface with the particle position vector as the center of gravity and the polyhedral element volume with the particle position as the vertex, χ is the magnetic susceptibility, μ ₀ is the permeability of vacuum, λ is the Lagrange's undetermined constant, and π is the circumference The rate, s, is the number of element boundary faces of a polyhedral element whose center of gravity is the particle position vector.

  The subscript ip of each physical quantity is the physical quantity of the p point in the polyhedron element whose center is the position vector of the i-th particle, and the subscript jq is the polyhedron element whose center is the position vector of the j-th particle. The physical quantity of q point is shown.

  The points p and q are described in step 104.

In equation (13), m is the mass of the particle, v is the velocity, and φ (r _i −r _j ) is the interaction potential energy between the i-th particle and the j-th particle. Subscripts i and j indicate physical quantities of the i and j-th particles, respectively.

Next, the canonical variables (in _{bold) H ip,} and (bold with _{neighbors) H ip,} and substituting Lagrangian of formula (13) to (15) to the Lagrange equation of motion, the motion equation of the magnetic field Can be described as in equation (16).

  Here, the second term on the right side of Equation (16) imposes a constraint (the divergence of the magnetization vector is 0) through Lagrange's undetermined constant.

  The third term on the right side of equation (16) is an attenuation term for quickly following when the externally applied magnetic field vector changes, and γ is an attenuation constant.

  In solving the equation of motion including the constraint shown in the second term on the right side of Expression (16), the present embodiment employs the SHAKE method, which is a generalized constraint introduction method.

  When equation (16) is discretized by the jumping method without considering the constraint, the following equations (17), (18), and (19) are obtained.

  Here, δt is a virtual time step used for calculating the convergence of the magnetization phenomenon.

  The subscript n is an arbitrary integer, a physical quantity at nδt, n−1 / 2 is a physical quantity at (n−1 / 2) δt, n + 1/2 is a physical quantity at (n + 1/2) δt, and n + 1 is (n + 1). This corresponds to the physical quantity at δt.

  The control unit 3 of the magnetic field analysis apparatus 1 reads the position vectors of the p point and the q point in the polyhedron element whose center of gravity is the particle position vector already calculated in step 104 and stored in the storage device 5.

  Next, the control unit 3 generates a magnetic field vector generated at the point p inside the polyhedron element whose center of gravity is the particle constituting the magnetic body, which is already calculated in step 105 and stored in the storage device 5. Read as a vector.

  Further, the control unit 3 calculates the polyhedral element boundary surface having the centroid of the element boundary area, normal vector, and particle position vector of the particles constituting the magnetic body, which are already calculated in step 104 and stored in the storage device 5. The volume of a polyhedron element whose vertex is the particle position is read.

  Further, the control unit 3 reads the attenuation constant, virtual mass, and virtual time increment stored in advance in the storage device 5.

  Next, the control unit 3 reads the magnetic field vector stored in advance in the storage device 5 and the initial value of the time differentiation of the magnetic field vector.

  In addition, the control part 3 reads the value, when the magnetic field vector and the time differentiation of the magnetic field vector are updated in step 107 mentioned later.

  The magnetization vector is calculated by substituting the magnetic field vector into a function (Potter-Schmulian equation) representing a magnetization curve stored in the storage device 5 as the input information 21 from step 101.

  Here, a specific calculation method of the magnetization vector using the Potter-Schmulian equation will be specifically described with reference to FIG.

  First, the Potter-Schmulian equation will be described.

  The Potter-Schmulian equation is a representation of the magnetization M (magnetization vector) as a function of the magnetic field H (magnetic field vector), as expressed by the following equations (a) and (b). Magnetic hysteresis can be expressed by setting three parameters. That is, a function representing a magnetization curve considering magnetic hysteresis can be obtained.

  Specifically, by setting three parameters Hc, Ms, and Mr from the equation (a), a function representing a magnetization curve corresponding to the parameters can be obtained.

  In addition, a change in magnetic hysteresis can be expressed by changing α.

  Next, a specific calculation method of the magnetization vector will be described.

  First, the control unit 3 of the magnetic field analysis apparatus 1 reads the initial conditions (Hc, Ms, Mr, α) stored in the storage device 5 using the calculation program 19 and substitutes them into the equations (a) and (b). (Step 201).

  Step 201 is executed only when step 106 is executed for the first time, and is not executed after the second time.

  Next, the control unit 3 of the magnetic field analysis apparatus 1 uses the calculation program 19 to calculate the calculation result of the most recent step calculated by the magnetic field analysis apparatus 1, specifically, the magnetic flux density B, the magnetic field H, and the like (A) , (A) are substituted (step 202).

  When the magnetic field analysis apparatus 1 executes step 106 for the first time, step 202 is not executed, and the initial conditions of step 201 are used for the magnetic flux density B, the magnetic field H, and the like.

  Next, the control unit 3 of the magnetic field analysis device 1 uses the calculation program 19 and the most recent magnetic flux density B (newly obtained magnetic flux density, first magnetic flux density) calculated by the magnetic field analysis device 1 and one step before. The difference dB = B−B ′ of the magnetic flux density B ′ (preliminarily calculated magnetic flux density or the magnetic flux density of the initial condition, the second magnetic flux density) obtained in the above is obtained (step 203).

  Next, the control unit 3 of the magnetic field analysis apparatus 1 obtains the α change determination value J using the calculation program 19 (step 204). Specifically, J is obtained from the following equation (c).

J = α × dB Formula (c)
Next, the control unit 3 of the magnetic field analysis apparatus 1 uses the calculation program 19 to determine whether the α change determination value J is 0 or more. If it is 0 or more, the process proceeds to step 208; Proceed to 206 (step 205).

  If J is smaller than 0, α is changed (step 206). This is because that α is smaller than 0 (a negative value) means that the current magnetization M (magnetization vector) reaches the corner of the magnetization curve (magnetic hysteresis).

  When α is changed, it is necessary to redo the calculation of magnetization, so the process returns to the analysis one step before the latest calculation calculated by the magnetic field analyzer 1 (step 207).

  If J is 0 or more, the magnetization M (magnetization vector) is calculated from the equation (a) without changing α (step 208).

  Note that the magnetic susceptibility χ can also be calculated as a function of the magnetic field H as shown in the following equation (d).

  Thus, in the present embodiment, the Potter-Schmulian equation is used as the magnetization curve.

  Therefore, it is possible to perform a magnetic field analysis considering the magnetic hysteresis inherent in the magnetic material.

  The above is the specific calculation method of the magnetization vector.

  The control unit 3 calculates the equations (17), (18), and (19) with respect to all p points in the polyhedron element having the position vector of the particles constituting the magnetic body as the center of gravity. The magnetic field vector is stored in the storage device 5.

  The details of step 106 in FIG. 5 have been described above.

  Next, the control unit 3 of the magnetic field analysis device 1 is bound to the magnetic field vector at the point p in the polyhedron element having the center of gravity as the position vector of the particles constituting the magnetic material, which is calculated in step 106 and stored in the storage device 5. A force is applied, and the calculated magnetic field vector is stored in the storage device 5 (step 107 in FIG. 5).

  Here, step 107 in FIG. 5 will be specifically described.

  A binding force is applied to the magnetic field vector at the point p in the polyhedron element with the position vector of the particle constituting the magnetic body as the center of gravity according to the following equations (20) and (21).

  Here, bold H is a magnetic field vector, α is a virtual mass, δt is a virtual time step, γ is an attenuation constant, N is the number of particles constituting the magnetic body, s is a position vector of the particles constituting the magnetic body, and the center of gravity. The subscript i is the physical quantity on the position vector of the i-th particle among the particles constituting the magnetic body, and ip is the i-th particle among the particles constituting the magnetic body. The physical quantity at point p in the polyhedron element with the position vector at the center of gravity, the subscript n is an arbitrary integer, the physical quantity at nδt, and n + 1 corresponds to the physical quantity at (n + 1) δt.

  The control unit 3 of the magnetic field analysis apparatus 1 uses the calculation program 19 to calculate the magnetic field vector at the point p in the polyhedron element with the particle constituting the magnetic body calculated in step 106 and stored in the storage device 5 as the center of gravity. Read.

  The control unit 3 reads the virtual mass, virtual time step, and attenuation constant stored in advance in the storage device 5.

  The control unit 3 performs a calculation based on Expression (20) and Expression (21) for each particle constituting the magnetic body, and stores the calculated magnetic field vector in the storage device 5.

  The above is the detail of step 107 in FIG.

  Next, the control unit 3 of the magnetic field analysis apparatus 1 determines whether or not the magnetic field vector obtained in step 107 satisfies the constraint condition, and if satisfied, proceeds to the next step, and if not satisfied, returns to step 107 (FIG. 5 step 108).

  Specifically, the control unit 3 performs the calculation based on the equation (22) using the magnetization vector calculated from the magnetic field vector of each of the particles constituting the magnetic body calculated in step 107 and stored in the storage device 5. .

  Here, the subscript n + 1 corresponds to the physical quantity at (n + 1) δt.

err _i is an error value for the binding condition of the i-th particle among the particles constituting the magnetic material.

  The magnetization vector (bold) M is calculated in step 107 and the magnetic field vector stored in the storage device 5 is read by the control unit 3 and substituted in the function representing the magnetization curve stored in the input information 21 from step 101. Calculated by The specific calculation method is the same as the method described in step 106 (see step 201 to step 208).

  For the normal vector (bold) n, the control unit 3 reads and substitutes the one already calculated in step 104 and stored in the storage device 5.

  If the error value does not satisfy the equation (23) for all particles, the control unit 3 returns to Step 107.

In Expression (23), A ₁ is an arbitrary error determination value, and an arbitrary value is stored in the storage device 5 in advance.

The above is the details of step 108 in FIG.

Next, the control unit 3 of the magnetic field analysis apparatus 1 determines whether the magnetization phenomenon of the magnetic material has reached a steady state, and proceeds to the next step if the condition is satisfied. (Step 109 in FIG. 5)
Here, step 109 in FIG. 5 will be specifically described.

Whether the particles constituting the magnetic material have reached a steady state is determined by the following equation (24).

Here, the bold letter H is the magnetic field vector, μ ₀ is the magnetic permeability of the vacuum, N is the number of particles constituting the magnetic body, and s is the number of faces of the polyhedral element whose center of gravity is the position vector of the particles constituting the magnetic body. The subscript ip is a physical quantity at the point p in the polyhedron element whose center of gravity is the position vector of the i-th particle among the particles constituting the magnetic body.

  Further, the subscript n is an arbitrary integer, which corresponds to a physical quantity at nδt, n−1 corresponds to (n−1) δt, and n + 1 corresponds to a physical quantity at (n + 1) δt.

A ₂ is an arbitrary error determination value for determining whether the magnetic field vector of the magnetic material has reached a steady state.

The control unit 3 of the magnetic field analysis apparatus 1 reads A ₂ and μ ₀ stored in advance in the storage device 5.

  Further, the control unit 3 reads the magnetic field vector at the point p in the polyhedron element having the center of gravity of the particles constituting the magnetic body, which is calculated in step 107 and stored in the storage device 5.

  Next, the control unit 3 calculates Expression (24), and proceeds to the next step if Expression (24) is satisfied, and returns to Step 106 if not satisfied.

The above is the details of step 109 in FIG.

  Next, the control unit 3 of the magnetic field analysis apparatus 1 performs the magnetic field vector, the magnetization vector, the magnetic flux density vector on the position vector of all particles constituting the magnetic body, and the magnetic field on the position vector of all particles constituting the coil. The vector and the magnetic flux density vector are calculated and stored in the storage device 5 (step 110 in FIG. 5).

  Step 110 will be specifically described below.

(Particles constituting magnetic material)
The magnetic field vector on the position vector of the particles constituting the magnetic body is expressed by the following equation (25).

Here, (bold) e _x = (1, 0, 0), (bold) e _y = (0, 1, 0), (bold) e _z = (0, 0, 1).

  The bold letter H is a magnetic field vector, the bold letter n is a normal vector, the suffix i represents the physical quantity on the position vector of the i-th particle of the particles constituting the magnetic body, and the suffix ip is the magnetic body. Is a physical quantity of p point in the polyhedron element whose center of gravity is the position vector of the i-th particle.

  The magnetization vector of the particles constituting the magnetic material can be obtained by substituting the magnetic field vector into a function representing the magnetization curve. The specific calculation method is the same as the method described in step 106 (see step 201 to step 208).

  The magnetic flux density vector on the position vector of the particles constituting the magnetic body is expressed by Expression (26).

  Here, the bold letter B is the magnetic flux density vector, the bold letter H is the magnetic field vector, the bold letter M is the magnetization vector, and the suffix i is the position vector of the i-th particle among the particles constituting the magnetic body. Indicates the physical quantity.

μ ₀ is the vacuum permeability.

  The control unit 3 of the magnetic field analysis apparatus 1 reads the magnetic field vector at the point p in the polyhedron element with the particle constituting the magnetic body as the center of gravity calculated up to step 110 from the storage device 5.

  Next, the control unit 3 reads the normal vector already calculated in step 104 and stored in the storage device 5, and calculates the magnetic field vector on the position vector of the particles constituting the magnetic body based on the equation (25). Calculate and store in the storage device 5.

  Further, the control unit 3 substitutes the magnetic field vector calculated in this step into a function representing the magnetization curve stored in the storage device 5 as the input information 21 from step 101, and the magnetization vector of the particles constituting the magnetic body. Is calculated and stored in the storage device 5. The specific calculation method is the same as the method described in step 106 (see step 201 to step 208).

  Next, the control unit 3 substitutes the magnetic field vector and the magnetization vector calculated in this step into Expression (26), calculates the magnetic flux density vector on the position vector of the particles constituting the magnetic body, and stores it in the storage device 5. .

  The control unit 3 reads the vacuum permeability stored in the storage device 5 in advance.

(Particulate particles)
The magnetic field vector on the position vector of the particles constituting the coil is described by Expression (27) and Expression (28).

  Here, the bold letter r is a position vector, the bold letter H is a magnetic field vector, N is the number of particles constituting the magnetic body, s and ΔS are the number of boundary faces and boundaries of the polyhedron element having the particle constituting the magnetic body as the center of gravity. Area, bold M is a magnetization vector, bold n is a normal vector, subscript i is a physical quantity on the position vector of the i-th particle among the particles constituting the coil, and j is a magnetic material Jq is the physical quantity on the position vector of the j-th particle, and jq is the physical quantity at the point q of the polyhedron element whose center of gravity is the position vector of the j-th particle among the particles constituting the magnetic body. Point to.

  The first term on the right side of Equation (27) is a term describing a magnetic field vector created on the position vector of the particles constituting the coil by the particles constituting the magnetic material, and is represented by Equation (28).

  The second term on the right side of the equation (27) is a magnetic field vector formed on the position vector of the particles constituting the coil by the particles constituting the coil, and is represented by equations (1) to (12).

  The control unit 3 of the magnetic field analysis apparatus 1 reads the position vector of the particles constituting the coil calculated in step 103 and stored in the storage device 5 from the storage device 5.

  In addition, when the position vector of the particle | grains which comprise a coil is updated in step 113, the control part 3 reads the value.

  Next, the control unit 3 calculates the number of particles constituting the magnetic body and the number of faces of the polyhedron element calculated from step 102 and stored in the storage device 5, and q points in the polyhedron element having the particle position vector as the center of gravity. The position vector of is read.

  Further, the control unit 3 reads the magnetization vector of the particles constituting the magnetic body, which is calculated in this step and stored in the storage device 5.

  Further, the control unit 3 reads the coil size calculated in step 103 and stored in the storage device 5 and the current density vector stored in the storage device 5 as the input information 21 in step 101.

  If the current density vector is updated in step 116, the control unit 3 reads the value.

  Thereafter, the control unit 3 calculates the magnetic field vector on the position vector of the particles constituting the coil in accordance with the equations (27), (28), and (1) to (12), and stores them in the storage device 5.

The control unit 3 converts the magnetic flux density vector (bold) B _i on the position vector of the particles constituting the coil into the magnetic field vector on the position vector of the particles constituting the coil that has already been calculated, and the vacuum permeability. Multiply and store in the storage device 5.

  The magnetic permeability of the vacuum is stored in the storage device 5 in advance.

  The details of step 110 have been described above.

  Thus, the control unit 3 analyzes the magnetization phenomenon of the magnetic material by solving the motion equation of the magnetic field.

  Therefore, it is not necessary to divide the entire space into meshes unlike the finite element method, and the magnetic field analysis can be efficiently performed with sufficient accuracy for a system having a magnetic body.

  In addition, since the equation of motion is solved, a matrix is not handled and the amount of memory necessary for the calculation is proportional to the number of particles.

  Next, the control unit 3 of the magnetic field analysis apparatus 1 calculates a force vector acting on all the particles constituting the SPM motor 31 and stores it in the storage device 5 (step 111 in FIG. 5).

  The step 111 in FIG. 5 will be specifically described below.

(Force vector acting on particles constituting magnetic material)
Let the canonical variables be bold r _i and bold r _i with side points, and the Lagrangian of equation (13) is substituted into the Lagrangian equation of motion.

  Then, the force acting on the particles constituting the magnetic material is described as shown in Expression (29).

Here, μ ₀ is the vacuum permeability, bold H is the magnetic field vector, bold M is the magnetization vector, bold r is the position vector, ΔV is the volume, and N is the number of particles constituting the magnetic body. I and j indicate the physical quantities of the i-th and j-th particles among the particles constituting the magnetic material, and iq represents the q in the polyhedral element having the ith particle as the center of gravity among the particles constituting the magnetic material. Indicates the physical quantity of a point.

Φ (r _i −r _j ) is an interaction potential energy between the i-th particle and the j-th particle.

  In this embodiment, the SPM motor is a rigid body (relative distance between particles is unchanged).

  Assuming that the model is a rigid model, the force acting on the particles constituting the magnetic material is described as in equation (30).

  The control unit 3 of the magnetic field analysis apparatus 1 reads the magnetic field vector at the point p in the polyhedron element with the particle constituting the magnetic body calculated up to step 110 and stored in the storage device 5 as the center of gravity.

  The control unit 3 calculates the magnetization vector at the p point by substituting the magnetic field vector at the p point into a function representing the magnetization curve stored in the storage device 5 as the input information 21 from step 101 to calculate the storage device 5. To remember. The specific calculation method is the same as the method described in step 106 (see step 201 to step 208).

  Further, the control unit 3 reads the volume of the polyhedron element having the vertex of the polyhedral element boundary surface having the particle center and the particle constituting the magnetic body calculated in step 104 and stored in the storage device 5 as the center of gravity.

  Next, the control unit 3 calculates the partial differentiation of the magnetic field vector at the point p in the polyhedral element with the particle constituting the magnetic body as the center of gravity, and calculates the force vector from the equation (30) based on the value. Store in the storage device 5.

(Force vector acting on the particles constituting the coil)
Of the particles constituting the coil, the position vector of the i-th particle is bold r _i , the current unit vector at that position vector is bold t _i , the magnetic flux density vector is bold B _i , and the current flow direction When the length of the coil and L _i, of the particles constituting the coil, the force vector acting on the i-th particle is described by the following equation (31).

  The control unit 3 of the magnetic field analysis apparatus 1 reads the position vector of the particles constituting the coil calculated in step 103 and stored in the storage device 5.

  If the position vector of the particles constituting the coil has been updated in step 113, the control unit 3 reads the value.

  Further, the control unit 3 reads the dimensions of the coil calculated in step 103 and stored in the storage device 5.

  Next, the control unit 3 reads the current density vector of the coil stored in the storage device 5 as the input information 21 from step 101.

  Note that if the current density vector of the coil is updated in step 116, the control unit 3 uses the value.

  The control unit 3 calculates a force vector acting on the particles constituting the coil based on the equation (31) and stores the force vector in the storage device 5.

  The above is the details of step 111 in FIG.

  Next, the control unit 3 of the magnetic field analyzer 1 calculates the rotation amount, angular velocity, translation amount, translation of the movable unit (moving element) from the force vector acting on the particles constituting the moving part (moving element) calculated in step 111. The speed is calculated and stored in the storage device 5 (step 112 in FIG. 5).

  The step 112 in FIG. 5 will be specifically described below.

  The equation of motion of the translational motion of the rigid body is expressed by equation (32).

Here, m _tot constitutes the mass of the movable portion (movable element), the position vector of the representative point arbitrarily determined in bold r _c is the movable portion (movable element), bold F movable portion (slider) The force vector acting on the particles, N, is the number of particles constituting the movable part (moving element).

  Next, the equation of motion of the rotational motion around an arbitrary fixed point is expressed by Equation (33-1), Equation (33-2), and Equation (33-3).

Here, I ₁ , I ₂ , and I ₃ are the main inertia moment of the rigid body, ω ₁ , ω ₂ , and ω ₃ are the angular velocities in the inertial principal axis coordinate system, and N ₁ , N ₂ , and N ₃ are the inertial principal axis coordinate systems. Torque.

  The rotor 31 of the SPM motor 31 of this embodiment is fixed to a certain rotating shaft through a rotor shaft 41.

  If the fixed rotation axis is the z axis, the motion of the rotor 31 is calculated as follows.

  Since the rotor 31 is fixed to the z-axis through the rotor shaft 41, the motion of the rotor 31 is calculated by solving the motion equation of the rotational motion around the z-axis.

  The equation of motion of the rotational motion of the rotor 31 around the z-axis is expressed by equation (34).

Here, I _z is the moment of inertia of the rotor 31, ω _z is the angular velocity in the z direction, and N _z is the torque in the z direction.

  When the equation (34) is discretized by the flying method, the following equations (35) and (36) are obtained.

Here, θ _z is a rotation amount around the z axis, and dt is a minute time step.

  The subscripts n−1 / 2, n, n + 1/2, and n + 1 correspond to physical quantities in (n−1 / 2) dt, ndt, (n + 1/2) dt, and (n + 1) dt, respectively. Yes.

  The control unit 3 reads the time step width in which an arbitrary value is stored in the storage device 5 in advance.

  Initial values of the rotation amount and the angular velocity are stored in the storage device 5 in advance in the storage device 5.

  If the rotation amount and angular velocity of the rotor 31 have already been updated in step 112, those values are used.

  The controller 3 calculates the moment of inertia from the three-dimensional structure and mass density of the rotor 31 stored in the storage device 5 as the input information 21 from step 101.

  Further, the control unit 3 calculates the equation (35) based on the position vector of the particles constituting the rotor 31 and the force vector acting on the particles constituting the rotor 31 calculated in step 111 and stored in the storage device 5. Equation (36) is calculated, and the rotation amount and angular velocity of the rotor 31 after dt seconds are calculated and stored in the storage device 5.

  The above is the details of step 112 in FIG.

  Next, the control unit 3 of the magnetic field analysis apparatus 1 calculates and stores the position vector of the particles constituting the moving element based on the rotation amount after a minute time and the movement amount of the center of gravity calculated in step 112 ( Step 113 in FIG.

  Next, the control unit 3 of the magnetic field analysis apparatus 1 sends the position vector, force vector, magnetic field vector, magnetic flux density vector, magnetization vector, and magnetic susceptibility χ of the particles constituting the SPM motor 31 via the printer port 11 as necessary. Output from the printer 12 (step 114 in FIG. 5).

  Next, the control unit 3 of the magnetic field analysis apparatus 1 determines whether or not a predetermined end condition (time, amount of movement, etc.) is satisfied. If it is satisfied, the analysis ends, and if not, the process proceeds to step 116. (Step 115 in FIG. 5).

  Next, the control unit 3 of the magnetic field analysis apparatus 1 updates the current density vector flowing in the coil 45 as necessary, and returns to Step 104 (Step 116 in FIG. 5).

  The above is the magnetic field analysis procedure according to the present embodiment.

  As described above, according to the present embodiment, the magnetic field analysis apparatus 1 calculates the magnetization phenomenon of the magnetic material by solving the equation of motion derived from the Lagrangian of the particle.

  Therefore, it is not necessary to divide the entire space into meshes unlike the finite element method, and the magnetic field analysis can be efficiently performed with sufficient accuracy for a system having a magnetic body.

  Further, since a matrix is not handled, the memory required for calculation is proportional to the number of particles (elements). Therefore, the memory can be greatly reduced as compared with the calculation that handles the asymmetric dense matrix.

  Furthermore, according to the present embodiment, the Potter-Schmulian equation is used as the magnetization curve.

  Therefore, it is possible to perform a magnetic field analysis considering the magnetic hysteresis inherent in the magnetic material.

  Further, according to the present embodiment, since the solution of the motion equation of the magnetic field derived from the Lagrangian of the particles constituting the magnetic body is calculated, the total energy of the system is preserved when the external input does not change. The

  In the above-described embodiment, the case where the magnetization phenomenon of the magnetic body is calculated by solving the equation of motion derived from the Lagrangian of the particles constituting the magnetic body has been described, but the present invention is not limited to this. Instead, any equation other than that derived from a Lagrangian may be used as long as it is an equation of motion.

  In the above-described embodiment, the present invention is used for the magnetic field analysis of the SPM motor 31. However, the present invention is not limited to this, and the present invention analyzes an arbitrary magnetic field in the space where the magnetic body is arranged. If so, for example, a motor other than the SPM motor 31 such as a linear motor or a magnetic shield surrounding the space with a magnetic material may be used for magnetic field analysis in the space shielded by the magnetic material.

DESCRIPTION OF SYMBOLS 1 ………… Magnetic field analyzer 3 ………… Control unit 5 ………… Storage device 7 ………… Input unit 9 ………… Display unit 11 ………… Printer port 12 ……… Printer 13 …… … Bus 31 ……… SPM motor 33 ……… Rotor 35 ……… Stator 37 ……… Rotor core 39 ……… Permanent magnet 41 ……… Rotor shaft 43 ……… Stator teeth 44 ……… Core back 45 …… ... coil 46 ... …… frame

Claims (19)

A magnetic field analysis device for analyzing a magnetic field at an arbitrary point in a space where a magnetic material is arranged using an information processing device,
A storage means for storing a shape of the magnetic body, an analysis condition having an expression representing a magnetization curve of the magnetic body, and a motion equation of the magnetic field;
From the analysis conditions, the coefficient of the motion equation of the magnetic field of the particles constituting the magnetic body, the initial value and the constant of the variable are calculated, the solution of the motion equation of the magnetic field is calculated, and based on the solution, A computing means for computing the magnetic field at an arbitrary point;
Have
The magnetic field analyzer according to claim 1, wherein an equation representing a magnetization curve of the magnetic material is a Potter-Schmulian equation. The magnetic field analyzer according to claim 1, wherein the Potter-Schmulian equation is expressed by the following equations (a) and (b).

  The magnetic field analysis apparatus according to claim 2, further comprising a calculation unit that obtains a magnetization curve using the Potter-Schmulian equation. The calculating means includes
Based on the magnetic flux density difference dB which is the difference between the first magnetic flux density B newly calculated by the computing means and the second magnetic flux density B ′ calculated in advance or in advance, α in the equation (A) is calculated. The magnetic field analysis apparatus according to claim 3, further comprising a determination unit that determines. The magnetic field analysis apparatus according to claim 4, wherein the determination unit changes α when an α change determination value J obtained from the following equation (c) is smaller than 0.
J = α × dB (C)
dB = B−B ′ (magnetic flux density difference)   The magnetic field analysis apparatus according to any one of claims 1 to 5, wherein the variable of the equation of motion of the magnetic field is a magnetic field vector of particles constituting the magnetic body. The coefficient of the equation of motion of the magnetic field is
A position vector of particles constituting the magnetic body;
A function representing a magnetization curve of the magnetic material;
The magnetic field analysis device according to claim 1, wherein the magnetic field analysis device is calculated from the following. The equation of motion of the magnetic field is
7. The particle Lagrangian including a magnetic field term is obtained by substituting the Lagrangian equation of motion with a magnetic field vector and a time derivative of the magnetic field vector as a canonical variable. The magnetic field analysis apparatus described. A conductor is further arranged in the space,
The magnetic field analysis apparatus according to claim 1, wherein the analysis condition includes a shape of a conductor and a current.   The program for functioning a computer as a magnetic field analyzer as described in any one of Claims 1-9. A storage means for storing a shape of the magnetic body, an analysis condition having an expression representing a magnetization curve of the magnetic body, and a motion equation of the magnetic field;
Based on the analysis conditions, the coefficient of the motion equation of the magnetic field of the particles constituting the magnetic body, the initial value of the variable and the constant are calculated, the solution of the motion equation of the magnetic field is calculated, and based on the solution, the magnetic body Computing means for computing the magnetic field at an arbitrary point in the space where
Control means for controlling the storage means and the computing means;
And an operation method of a magnetic field analysis apparatus for analyzing a magnetic field at an arbitrary point in a space where a magnetic body is arranged, wherein the control means includes:
Wherein the calculating means, from the analysis conditions, said initial value of the coefficients and variables of the magnetic field equations of motion of the particles constituting the magnetic body and the constants calculated, by calculating the solution of the equation of motion of the magnetic field,
Based on the solution of the equation of motion of the magnetic field, and step (b) of Ru was calculated magnetic field at an arbitrary point in the space,
Have
A method for operating a magnetic field analysis apparatus, wherein an equation representing a magnetization curve of the magnetic material is a Potter-Schmulian equation. 12. The operation method of the magnetic field analysis apparatus according to claim 11, wherein the Potter-Schmulian equation is expressed by the following equations (a) and (b).

13. The method of operating a magnetic field analysis apparatus according to claim 12 , wherein the control means includes a step (c) of obtaining the magnetization of the magnetic body using the Potter-Schmulian equation. In the step (c), the control means
A second magnetic flux density B determined anew in the step (b), the initial value or the first of the formula based on the magnetic flux density difference dB which is the difference between the magnetic flux density B 'calculated in advance (A) The method of operating a magnetic field analysis apparatus according to claim 13, further comprising a step (d) of determining α. 15. The method of operating a magnetic field analysis apparatus according to claim 14, wherein in the step (d), α is changed when an α change determination value J obtained from the following equation (c) is smaller than 0.
J = α × dB (C)
dB = B−B ′ (magnetic flux density difference) 16. The method of operating a magnetic field analysis apparatus according to claim 11, wherein the variable of the equation of motion of the magnetic field is a magnetic field vector of particles constituting the magnetic body. The coefficient of the equation of motion of the magnetic field is
A position vector of particles constituting the magnetic body;
A function representing a magnetization curve of the magnetic material;
The operation method of the magnetic field analysis apparatus according to claim 11, wherein the operation method is calculated from: The equation of motion is
18. The method according to claim 11, wherein the Lagrangian of a particle including a magnetic field term is obtained by substituting a Lagrangian equation of motion with a magnetic field vector and a time derivative of the magnetic field vector as a canonical variable. An operation method of the magnetic field analysis apparatus described. A conductor is further arranged in the space,
The operation method of the magnetic field analysis apparatus according to claim 11, wherein the analysis condition includes a shape of a conductor and a current.

Images