JP2007280210A - Simulator, simulation method and program - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide a simulator for simulating a circuit at a high speed, while accurately calculating a temporal change in inductance of an inductance element, in the circuit including the inductance element for generating electromagnetic induction. <P>SOLUTION: This simulator simulates motion of a circuit including a coil and a movable part magnetically joined to the coil, and determines whether or not an electric current value to the coil between a predetermined time step and a previous time step, an electric current variation and a position variation as differences of a position of the movable part, are larger than a predetermined threshold value, and calculates the inductance L of the coil on the basis of the electric current value in the predetermined time step and the position of the movable part, when the electric current variation or the position variation is larger than the threshold value, and stores the calculated L, and outputs the stored L. When the electric current variation and the position variation are the threshold value or less, the L is not calculated, and the L at the previous time step is outputted. <P>COPYRIGHT: (C)2008,JPO&INPIT

Description

  The present invention relates to a simulator, a simulation method, and a program. In particular, the present invention relates to a simulator and a simulation method for simulating a circuit including an inductance element, and a program for the simulator.

Analyzing the control state of an electromagnetic field analysis unit that analyzes the electromagnetic field state of the electromagnetic mechanism unit of the electromagnetic drive device, a circuit analysis unit that analyzes the circuit state of the circuit configuration unit of the electromagnetic drive device, and a control unit of the electromagnetic drive device A simulation system including an analysis device having a control analysis unit is known (see, for example, Patent Document 1). In this simulation system, an electromagnetic field analysis unit, a circuit analysis unit, and a control analysis unit have independent configurations. The values representing current, inductance, induced voltage, force, and displacement as time functions are mutually coupled as parameters.
JP-A-8-0776814

  For example, when simulating the operation of a motor, the inductance of the circuit through which the drive current flows changes with time in accordance with the drive current passed through the motor and the rotation of the rotor. By calculating such a temporal change in inductance with higher accuracy and simulating the circuit, it is possible to simulate the movement of the motor with higher accuracy. However, if it is attempted to calculate the time variation of the inductance with higher accuracy, the amount of calculation required for the magnetic field analysis increases, and as a result, the time required for the entire simulation increases. Patent Document 1 does not describe a method of simulating a circuit at high speed while calculating such a change in inductance value over time with higher accuracy.

  Then, an object of this invention is to provide the simulator, simulation method, and program which can solve said subject. This object is achieved by a combination of features described in the independent claims. The dependent claims define further advantageous specific examples of the present invention.

  According to the first aspect of the present invention, there is provided a simulator for simulating a circuit including an inductance element, a current calculation unit for calculating a current value flowing through the inductance element for each time step, and a predetermined value calculated by the current calculation unit A current change amount determination unit for determining whether or not a current change amount indicating a difference between the current value in the time step and the current value in the time step before the predetermined time step is larger than a predetermined current change amount; When the current change amount determination unit determines that the current change amount is larger than the predetermined current change amount, the current change amount determination unit determines whether the current change amount is calculated based on the current value at the predetermined time step calculated by the current calculation unit. The magnetic field analysis unit that calculates the inductance of the inductance element, and the inductance at a predetermined time step calculated by the magnetic field analysis unit are recorded. Output parameter storage unit, the parameter output unit that outputs the inductance at a predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step, and the current change amount determination unit, the current change amount is When it is determined that the current change amount is equal to or less than a predetermined current change amount, the magnetic field analysis unit does not calculate the inductance at a predetermined time step of the inductance element, and the inductance at the time step before the output parameter storage unit stores Is output from the parameter output unit to the current calculation unit.

  A time change amount determination unit that determines whether or not a time change amount indicating a difference between a time indicating a predetermined time step and a time indicating a previous time step is larger than a predetermined time change amount; The magnetic field analysis unit determines that the time change amount determination unit determines that the time change amount is greater than a predetermined time change amount, and the current change amount determination unit determines that the current change amount is a predetermined current change amount. When it is determined that the value is larger, the inductance of the inductance element in the predetermined time step is calculated based on the current value in the predetermined time step, and the output parameter storage unit is in the predetermined time step calculated by the magnetic field analysis unit. The inductance is stored, and the parameter output unit outputs the inductance stored in the output parameter storage unit to the current calculation unit at each time step, The analysis control unit causes the magnetic field analysis unit to calculate the inductance of the inductance element at a predetermined time step when the time change amount determination unit determines that the time change amount is equal to or less than a predetermined time change amount. Instead, the inductance at the previous time step stored in the output parameter storage unit may be output from the parameter output unit to the current calculation unit.

  The current calculation unit calculates a current value flowing through each of the plurality of inductance elements magnetically coupled to each other for each time step, and the current change amount determination unit calculates the maximum of the current change amount in each of the plurality of inductance elements. The magnetic field analysis unit determines that the maximum value of the current change amount is larger than the predetermined current change amount. In this case, the inductance of each of the plurality of inductance elements in the predetermined time step is calculated based on the current value of each of the plurality of inductance elements in the predetermined time step, and the output parameter storage unit is calculated by the magnetic field analysis unit. Memorize each inductance of multiple inductance elements at a given time step Outputs the inductance of each of the plurality of inductance elements stored in the output parameter storage unit to the current calculation unit for each time step, and the magnetic field analysis control unit determines that the current change amount determination unit has a maximum current change amount. If it is determined that the value is equal to or less than the current change amount, the magnetic field analysis unit does not calculate the inductance of each of the plurality of inductance elements at a predetermined time step, and the time before the output parameter storage unit stores The inductance of each of the plurality of inductance elements in the step may be output from the parameter output unit to the current calculation unit.

  The current calculation unit calculates a current value to be supplied to each of the plurality of inductance elements magnetically coupled to each other for each time step, and the current change amount determination unit has a current to be supplied to each of the plurality of inductance elements in a predetermined time step. The total current change amount indicating the difference between the total value of the values and the current value in the previous time step is greater than a predetermined total current change amount, and the magnetic field analysis unit When the change amount determination unit determines that the total current change amount is larger than a predetermined current change amount, based on the current values of the plurality of inductance elements in the predetermined time step, in the predetermined time step The inductance of each of the plurality of inductance elements is calculated, and the output parameter storage unit is a predetermined time step calculated by the magnetic field analysis unit. The inductance output unit stores the inductances of the plurality of inductance elements, and the parameter output unit stores the inductances of the plurality of inductance elements in the predetermined time step stored in the output parameter storage unit in the current calculation unit for each time step. When the current change amount determination unit determines that the total current change amount is equal to or less than a predetermined total current change amount, the magnetic field analysis control unit outputs a plurality of values in a predetermined time step to the magnetic field analysis unit. The inductance of each of the plurality of inductance elements in the previous time step stored in the output parameter storage unit may be output from the parameter output unit to the current calculation unit without calculating the inductance of each of the inductance elements.

  A position calculation unit that calculates the position of the movable part that is magnetically coupled to the inductance element and is driven by the magnetic field generated by the inductance element for each time step, and the position of the movable part at a predetermined time step calculated by the position calculation unit And a position change amount determination unit that determines whether or not a position change amount indicating a difference from the position of the movable part in the previous time step is larger than a predetermined position change amount, and the magnetic field analysis unit includes: When the position change amount determination unit determines that the position change amount is larger than the predetermined position change amount, or the current change amount determination unit indicates that the current change amount is larger than the predetermined current change amount. Is determined based on the current value at the predetermined time step calculated by the current calculation unit and the position at the predetermined time step calculated by the position calculation unit. The output parameter storage unit stores the inductance at a predetermined time step calculated by the magnetic field analysis unit, and the parameter output unit stores the inductance stored in the output parameter storage unit. Output to the current calculation unit at each time step, and the magnetic field analysis control unit determines that the position change amount determination unit determines that the position change amount is equal to or less than a predetermined position change amount, and the current change amount determination unit. However, when it is determined that the current change amount is equal to or less than the predetermined current change amount, the output parameter storage unit does not cause the magnetic field analysis unit to calculate the inductance of the inductance element at a predetermined time step. The inductance in the time step before the current may be output from the parameter output unit to the current calculation unit.

  A position calculating unit that calculates the position of a movable unit that is magnetically coupled to the inductance element and is driven by a magnetic field generated by the inductance element is provided for each time step. When it is determined that the amount of change is greater than a predetermined amount of time change, a predetermined value is determined based on the current value calculated in the predetermined time step by the current calculation unit and the position in the predetermined time step calculated by the position calculation unit. The induced electromotive force due to the movable part in the inductance element at the time step is calculated, the output parameter storage unit stores the induced electromotive force due to the movable part at the predetermined time step calculated by the magnetic field analysis unit, and the magnetic field analysis control unit is When the time variation determination unit determines that the time variation is equal to or less than a predetermined time variation, the magnetic field solution The parameter output unit does not calculate the induced electromotive force due to the movable part stored in the output parameter storage unit for each time step. May be output.

  The magnetic field analysis unit is configured such that when the time change amount is greater than a predetermined time change amount, the time change amount determination unit calculates the current value at the predetermined time step calculated by the current calculation unit and the predetermined value calculated by the position calculation unit. Based on the position in the time step, the electromagnetic force of the movable part in the predetermined time step is calculated, and the output parameter storage unit stores the electromagnetic force of the movable part in the predetermined time step calculated by the magnetic field analysis unit, and the magnetic field analysis When the time change amount determination unit determines that the time change amount is equal to or less than a predetermined time change amount, the control unit causes the magnetic field analysis unit to calculate the electromagnetic force of the movable unit at a predetermined time step. Instead, the parameter output unit may output the electromagnetic force of the movable unit stored in the output parameter storage unit to the position calculation unit for each time step.

  The magnetic field analysis unit, based on the current value of the inductance element at a predetermined time step calculated by the current calculation unit, calculates the flux linkage of the inductance element by the current indicated by the current value, and the current calculation The interlinkage flux of the inductance element is calculated by the current indicated by the current value changed by a minute current value from the current value at the predetermined time step calculated by the unit, and the interlinkage flux calculated by the interlinkage flux calculation unit is calculated. An inductance calculating unit that calculates the inductance of the inductance element based on a value obtained by dividing a difference from the magnetic flux by a minute current value may be included.

  The current calculation unit calculates a current value flowing through each of the plurality of inductance elements magnetically coupled to each other for each time step, and the flux linkage calculation unit calculates the first inductance in the predetermined time step calculated by the current calculation unit. The interlinkage magnetic flux of each of the first inductance element and the second inductance element is calculated from the respective current values of the element and the second inductance element, and the inductance calculation unit is configured to calculate the first time step in a predetermined time step calculated by the current calculation unit. A first interlinkage magnetic flux calculation unit for calculating the interlinkage magnetic flux of each of the first inductance element and the second inductance element based on the current indicated by the current value changed by a minute current value from the current value of the inductance element; The flux linkage of the first inductance element calculated by the magnetic flux calculator and the first flux calculated by the flux linkage calculator A first self-inductance calculation unit that calculates the self-inductance of the first inductance element based on a value obtained by dividing the difference from the flux linkage of the conductance element by a minute current value, and a first self-inductance calculation unit calculated by the first linkage magnetic flux calculation unit. Based on the value obtained by dividing the difference between the linkage flux of the two inductance elements and the linkage flux of the second inductance element calculated by the linkage flux calculating unit by the minute current value, the second inductance between the first inductance element and the second inductance element is calculated. A first mutual inductance calculation unit for calculating a mutual inductance of the element, and a first inductance based on a current indicated by a current value changed by a minute current value from the current value of the second inductance element at a predetermined time step calculated by the current calculation unit. A second interlinkage magnetic flux calculator that calculates the interlinkage magnetic flux of each of the element and the second inductance element, and a second interlinkage magnetic flux calculator Based on the value obtained by dividing the difference between the flux linkage of the second inductance element and the flux linkage of the second inductance element calculated by the flux linkage calculator by the minute current value, the self inductance of the second inductance element is calculated. A difference between the linkage magnetic flux of the first inductance element calculated by the second self-inductance calculation unit calculated by the second flux linkage calculation unit and the linkage flux of the first inductance element calculated by the linkage flux calculation unit is minute. A second mutual inductance calculation unit that calculates the mutual inductance of the first inductance element between the second inductance element and the second inductance element based on the value divided by the current value may be included.

  The inductance calculation unit includes a self-inductance of the first inductance element in the time step before the first self-inductance calculation unit calculates, a mutual inductance of the second inductance element in the time step before the first mutual inductance calculation unit, and An inductance storage unit that stores the self-inductance of the second inductance element in the previous time step calculated by the second self-inductance calculation unit and the mutual inductance of the first inductance element in the previous time step calculated by the second mutual inductance calculation unit The first inductance in the time step before the self-inductance of the first inductance element and the inductance storage unit stored in the predetermined time step calculated by the first self-inductance calculation unit The difference between the inductance element and the self-inductance of the conductance element is smaller than a predetermined first value, and the mutual inductance of the second inductance element and the inductance storage section at a predetermined time step calculated by the first mutual inductance calculation section are stored in the memory. When the difference from the mutual inductance of the second inductance element in the previous time step is smaller than a predetermined second value, the flux linkage of the first inductance element is relative to the current value of the first inductance element. When determining that the linear region determining unit and the linear region determining unit are not in the linear region, the second self-inductance calculating unit and the second mutual inductance calculating unit include Calculate the mutual inductance of the inductance and the first inductance element When the linear region determination unit determines that the linear region is in the linear region, the second self-inductance calculation unit and the second mutual inductance calculation unit calculate the self-inductance of the second inductance element and the mutual inductance of the first inductance element, respectively. In addition, an inductance calculation control unit that outputs the self-inductance of the second inductance element and the mutual inductance of the first inductance element in the previous time step stored in the inductance storage unit may be further included.

  A time change amount determination unit for determining whether a time change amount indicating a difference between a time indicating a predetermined time step and a time indicating a previous time step is larger than a predetermined time change amount; and an inductance element A position calculation unit that calculates the position of the movable unit that is magnetically coupled and driven by the magnetic field generated by the inductance element for each time step, and the position of the movable unit at the predetermined time step calculated by the position calculation unit, A position change amount determination unit for determining whether or not a position change amount indicating a difference from the position of the movable part in the previous time step is larger than a predetermined position change amount; The position change indicates that the position change amount is greater than the predetermined position change amount, provided that the time change amount determination unit determines that the change amount is greater than the predetermined time change amount. The current value and position at the predetermined time step calculated by the current calculation unit when the determination unit determines, or when the current change amount determination unit determines that the current change amount is larger than a predetermined current change amount. Based on the position at the predetermined time step calculated by the calculation unit, the inductance of the inductance element at the predetermined time step is calculated, and the output parameter storage unit stores the inductance at the predetermined time step calculated by the magnetic field analysis unit, The parameter output unit outputs the inductance stored in the output parameter storage unit to the current calculation unit for each time step, and the magnetic field analysis control unit indicates that the time change amount is equal to or less than a predetermined time change amount. When the time change amount determination unit makes a determination, or the position change amount is equal to or less than a predetermined position change amount. If the current change amount determination unit determines that the current change amount is equal to or less than a predetermined current change amount, the magnetic field analysis unit determines whether the current change amount is equal to or less than a predetermined current step. The inductance in the previous time step stored in the output parameter storage unit may be output from the parameter output unit to the current calculation unit without calculating the inductance.

  A position calculating unit that calculates the position of the movable unit that is magnetically coupled to the inductance element and is driven by the magnetic field generated by the inductance element, for each time step; and the magnetic field analysis unit includes a predetermined value calculated by the current calculation unit. And further including an induced electromotive force calculation unit for calculating an induced electromotive force by the movable part in the inductance element in the predetermined time step based on the current value in the time step and the position in the predetermined time step calculated by the position calculation unit, The magnetic flux calculation unit indicates the current value based on the current value flowing through the plurality of inductance elements at the predetermined time step calculated by the current calculation unit and the position of the movable unit at the predetermined time step calculated by the position calculation unit. Calculate the flux linkage of the inductance element by the current and the moving part, the inductance calculation part is Based on the current value when only a minute current value is changed from the current value calculated at a predetermined time step calculated by the current calculation unit and the position of the movable unit at the predetermined time step calculated by the position calculation unit, the current value is changed by a small current value. Based on the value obtained by dividing the difference between the calculated flux linkage and the linkage flux calculated by the linkage flux calculation unit by the minute current value. The inductance of the inductance element is calculated, and the induced electromotive force calculation unit calculates the magnetic flux change value obtained by subtracting the linkage magnetic flux in the previous time step from the linkage magnetic flux in the predetermined time step calculated by the linkage magnetic flux calculation unit. From the inductance calculated by the unit at a predetermined time step and the current value at the predetermined time step calculated by the current calculation unit. Current change value obtained by subtracting the current value of flops, and based on the difference between the time indicated by the time the previous time step indicated by a predetermined time step may calculate the induced electromotive force by the movable portion.

  According to a second aspect of the present invention, there is provided a simulation method for simulating a circuit including an inductance element, wherein a current value flowing through the inductance element is calculated at each time step, and is calculated at a current calculation stage. Current change amount determination for determining whether or not a current change amount indicating a difference between a current value in a predetermined time step and a current value in a time step before the predetermined time step is larger than a predetermined current change amount And when the current change amount is determined to be larger than a predetermined current change amount in the current change amount determination step, based on the current value in the predetermined time step calculated in the current calculation step. In the magnetic field analysis stage to calculate the inductance of the inductance element in the time step of An output parameter storage stage for storing the output inductance in the predetermined time step, a parameter output stage for outputting the inductance in the predetermined time step stored in the output parameter storage stage for each time step in the current calculation stage; In the current change amount determination step, when it is determined that the current change amount is equal to or less than a predetermined current change amount, the output parameter storage step is performed without causing the magnetic field analysis step to calculate the inductance at a predetermined time step of the inductance element. And a magnetic field analysis control stage for outputting the inductance in the previous time step stored in the current calculation stage in the parameter output stage.

  According to the third aspect of the present invention, there is provided a program for a simulator for simulating a circuit including an inductance element, wherein the simulator calculates a current value flowing through the inductance element for each time step, a current calculation unit A current for determining whether or not a current change amount indicating a difference between a current value at a predetermined time step calculated by the current value and a current value at a time step before the predetermined time step is larger than a predetermined current change amount. When the change amount determination unit and the current change amount determination unit determine that the current change amount is larger than a predetermined current change amount, a predetermined value is determined based on the current value at a predetermined time step calculated by the current calculation unit. The magnetic field analysis unit that calculates the inductance of the inductance element at the time step of Output parameter storage unit for storing the inductance, a parameter output unit for outputting the inductance at a predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step, and a current change amount determination unit for changing the current When it is determined that the amount is equal to or less than a predetermined current change amount, the magnetic field analysis unit does not calculate the inductance at a predetermined time step of the inductance element, and the time step before the output parameter storage unit stores Is caused to function as a magnetic field analysis control unit that outputs the inductance from the parameter output unit to the current calculation unit.

  According to a fourth aspect of the present invention, there is provided a simulator for simulating a circuit including an inductance element, wherein the position of a movable part magnetically coupled to the inductance element and driven by a magnetic field generated by the inductance element is time-stepped. A position calculation unit that is calculated every time, a current calculation unit that calculates a current value flowing through the inductance element for each time step, a position of the movable unit at a predetermined time step calculated by the position calculation unit, and a position before the time step. A position change amount determination unit that determines whether or not a position change amount indicating a difference from the position of the movable part in the time step is larger than a predetermined position change amount, and a position change amount determination unit, the position change amount is The current value and position calculation at a predetermined time step calculated by the current calculation unit when it is determined that the position change amount is greater than a predetermined amount. A magnetic field analysis unit that calculates the inductance of the inductance element at the predetermined time step based on the position at the predetermined time step calculated by the unit, and an output parameter storage unit that stores the inductance at the predetermined time step calculated by the magnetic field analysis unit A parameter output unit that outputs the inductance at a predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step, and a position change amount determination unit that determines the position change amount in advance. When it is determined that the change amount is less than or equal to the amount of change, the magnetic field analysis unit does not calculate the inductance of the inductance element at the predetermined time step, and the inductance at the previous time step stored in the output parameter storage unit is used as the parameter output unit. Analysis of the magnetic field output from the current to the current calculator And a control unit.

  According to a fifth aspect of the present invention, there is provided a simulation method for simulating a circuit including an inductance element, wherein the position of a movable part that is magnetically coupled to the inductance element and is driven by a magnetic field generated by the inductance element is determined in time. A position calculation stage to be calculated for each step; a current calculation stage to calculate a current value to be passed through the inductance element for each time step; a position of the movable part at a predetermined time step calculated in the position calculation stage; In the position change amount determination stage for determining whether or not the position change amount indicating the difference from the position of the movable part in the previous time step is larger than a predetermined position change amount, the position change amount in the position change amount determination stage When it is determined that the amount is larger than a predetermined position change amount, it is calculated in the current calculation stage. A magnetic field analysis stage for calculating the inductance of the inductance element at the predetermined time step based on the current value at the predetermined time step and the position at the predetermined time step calculated at the position calculation stage, and the predetermined value calculated at the magnetic field analysis stage. An output parameter storage stage for storing the inductance at each time step, a parameter output stage for outputting the inductance at a predetermined time step stored at the output parameter storage stage at each time step to the current calculation stage, and position change determination In the stage, when it is determined that the position change amount is equal to or less than a predetermined position change amount, the inductance at a predetermined time step of the inductance element is not calculated in the magnetic field analysis stage, and the output parameter storage stage. The inductance in the time step before being stored Te, and a magnetic field analysis control step of outputting the current calculation step from the parameter output stage.

  According to a sixth aspect of the present invention, there is provided a program for a simulator for simulating a circuit including an inductance element, the simulator being magnetically coupled to the inductance element and being driven by a magnetic field generated by the inductance element. A position calculation unit that calculates the position of the unit at each time step, a current calculation unit that calculates a current value flowing through the inductance element at each time step, a position of the movable unit at a predetermined time step calculated by the position calculation unit, and the time The position change amount determination unit and the position change amount determination unit for determining whether or not the position change amount indicating the difference from the position of the movable part in the time step before the step is larger than a predetermined position change amount A predetermined time step calculated by the current calculation unit when it is determined that the change amount is larger than a predetermined position change amount. A magnetic field analysis unit that calculates the inductance of the inductance element at a predetermined time step based on the current value and the position at the predetermined time step calculated by the position calculation unit, and stores the inductance at the predetermined time step calculated by the magnetic field analysis unit Output parameter storage unit, a parameter output unit that outputs the inductance at a predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step, and a position change amount determination unit. When it is determined that the position change amount is equal to or less than the position change amount, the magnetic field analysis unit does not calculate the inductance at a predetermined time step of the inductance element, and the inductance at the previous time step stored in the output parameter storage unit is From parameter output section to current calculation section To function as a magnetic field to force the analysis control unit.

  According to the seventh aspect of the present invention, there is provided a simulator for simulating a circuit including an inductance element, and a current value at a predetermined time step calculated by a current calculation unit that calculates a current value flowing through the inductance element for each time step. And a current change amount determination unit for determining whether or not a current change amount indicating a difference between the current value in the time step before the predetermined time step is larger than a predetermined current change amount, and a current change amount determination When the unit determines that the current change amount is larger than the predetermined current change amount, the inductance of the inductance element at the predetermined time step is calculated based on the current value at the predetermined time step calculated by the current calculation unit. A magnetic field analysis unit to be calculated, and an output parameter for storing the inductance at a predetermined time step calculated by the magnetic field analysis unit. Data storage unit, a parameter output unit that outputs the inductance at a predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step, and a current change amount determination unit, the current change amount When it is determined that the current change amount is equal to or less than the predetermined current change amount, the magnetic field analysis unit does not calculate the inductance of the inductance element at the predetermined time step, and the inductance at the previous time step stored in the output parameter storage unit is calculated. And a magnetic field analysis control unit that outputs the parameter output unit to the current calculation unit.

  According to an eighth aspect of the present invention, there is provided a simulation method for simulating a circuit including an inductance element, wherein the current at a predetermined time step is calculated by a current calculation unit that calculates a current value flowing through the inductance element for each time step. A current change amount determining step for determining whether or not a current change amount indicating a difference between the value and a current value in a time step before a predetermined time step is greater than a predetermined current change amount; and a current change amount In the determination stage, when it is determined that the current change amount is larger than the predetermined current change amount, the inductance element at the predetermined time step is determined based on the current value at the predetermined time step calculated by the current calculation unit. The magnetic field analysis stage for calculating the inductance of the current and the predetermined time step calculated in the magnetic field analysis stage An output parameter storage stage for storing the inductance, a parameter output stage for outputting the inductance at a predetermined time step stored in the output parameter storage stage to the current calculation unit for each time step, and a current change amount determination stage When it is determined that the amount of change is equal to or less than a predetermined amount of current change, the inductance before the inductance element is calculated at a predetermined time step in the magnetic field analysis stage, and the time before being stored in the output parameter storage stage A magnetic field analysis control step of causing the current calculation unit to output the inductance in the step in the parameter output step.

  According to the ninth aspect of the present invention, there is provided a program for a simulator for simulating a circuit including an inductance element, wherein the simulator calculates a current value flowing through the inductance element for each time step. Current change amount determination unit for determining whether or not a current change amount indicating a difference between a current value in a time step of the current step and a current value in a time step before a predetermined time step is greater than a predetermined current change amount When the current change amount determination unit determines that the current change amount is larger than the predetermined current change amount, the current change amount determination unit determines whether the current change amount is calculated based on the current value at the predetermined time step calculated by the current calculation unit. Magnetic field analysis unit for calculating the inductance of the inductance element, the inductor at a predetermined time step calculated by the magnetic field analysis unit An output parameter storage unit that stores the output, an inductance at a predetermined time step stored in the output parameter storage unit, a parameter output unit that outputs the current calculation unit for each time step, and a current change amount determination unit Is determined to be equal to or less than a predetermined current change amount, the magnetic field analysis unit does not calculate the inductance of the inductance element at a predetermined time step, and the output parameter storage unit stores the inductance at the previous time step. The inductance is caused to function as a magnetic field analysis control unit that outputs the parameter output unit to the current calculation unit.

  Note that the above summary of the invention does not enumerate all the necessary features of the present invention, and sub-combinations of these feature groups can also be the invention.

  According to the present invention, in a circuit including an inductance element in which electromagnetic induction occurs, a simulation capable of reducing the time required to simulate the circuit while calculating the temporal change in inductance of the inductance element with high accuracy. An environment can be provided.

  Hereinafter, the present invention will be described through embodiments of the invention. However, the following embodiments do not limit the claimed invention, and all combinations of features described in the embodiments are inventions. It is not always essential to the solution.

  FIG. 1 shows an example of a control target of the simulator 100. Simulator 100 simulates the operation of electromagnetic drive system 150 including driver 160, inverter 170, and motor 180. Motor 180 is driven by a magnetic field generated by the current supplied from inverter 170. The driver 160 controls the operation of the motor 180 by controlling the voltage applied to the motor 180 by the inverter 170. The simulator 100 includes a circuit / control simulator 110 and an electromagnetic field simulator 120. The circuit / control simulator 110 simulates a circuit through which a current for driving the motor 180 flows. Specifically, the circuit / control simulator 110 simulates an equivalent circuit 190 of the motor 180. The electromagnetic field simulator 120 analyzes the magnetic field using a method such as a finite element method based on the current supplied to the motor 180, the material of the constituent material constituting the motor 180, and the physical arrangement thereof. Simulate the characteristics of the motor 180 in space.

Taking the case where the motor 180 is a three-phase synchronous motor as an example, the equivalent circuit 190 includes an inductance L, a resistance R of each of the three armature circuits through which the armature current flows, and an induced electromotive force generated by the rotor of the motor 180. V _i is a circuit parameter. The circuit / control simulator 110 calculates the angle θ of the rotor of the motor 180 based on the equation of motion of the electromagnetic drive system 150 including the load on the motor 180 at one time step in the simulation, and sets the control target value as the control target value. The voltage V applied to each armature circuit by the inverter 170 is determined according to the control target value. Then, the circuit / control simulator 110 determines the current I flowing through each of the armature circuits based on the applied voltage V and the circuit parameters of the equivalent circuit 190. Θ and I determined by the circuit / control simulator 110 are passed to the electromagnetic field simulator 120 together with a time t indicating a time step. Then, the electromagnetic field simulator 120 calculates the magnetic field in the motor 180 based on the current I and the angle θ, and analyzes the calculated magnetic field and the change in the magnetic field from the previous time step, thereby obtaining the L of the armature circuit, R and V _i are calculated. The electromagnetic field simulator 120 calculates the torque T of the motor 180 by analyzing the calculated magnetic field. Then, the electromagnetic field simulator 120 passes the calculated L, R, V _i , and T to the circuit / control simulator 110. Thereby, the calculation in one time step is completed. In the next time step, the circuit / control simulator 110 calculates the angle θ based on T passed from the electromagnetic field simulator 120, determines the new control target value, and passes it from the electromagnetic field simulator 120. A current is calculated based on L, R, and V _i . The circuit / control simulator 110 and the electromagnetic field simulator 120 simulate the operation of the electromagnetic drive system 150 by calculating the above calculation steps at appropriate time steps. As described above, the simulator 100 performs a coupled simulation by the circuit / control simulator 110 and the electromagnetic field simulator 120.

  The electromagnetic field simulator 120 calculates L on the condition that the change amount of the current value I from the current value in the previous time step when the circuit parameter was calculated last time is larger than a predetermined threshold value. The electromagnetic field simulator 120 may calculate the inductance L even when the amount of change in the angle θ from the rotor angle in the previous time step is larger than a predetermined threshold. Further, the electromagnetic field simulator 120 calculates L on the condition that the amount of change in time from the time indicated by the previous time step is larger than a predetermined threshold. When the circuit parameter is calculated, the electromagnetic field simulator 120 stores the value. When the circuit parameter is not calculated, the electromagnetic field simulator 120 passes the stored value of each circuit parameter to the circuit / control simulator 110. In this way, the simulator 100 can reduce the amount of calculation in the coupled simulation.

  The electromagnetic field simulator 120 calculates a value obtained by dividing the change in the linkage flux of the armature circuit by the minute change when the value of the current I is minutely changed as the inductance. Thereby, the electromagnetic field simulator 120 can include the nonlinear characteristic in the armature circuit in the inductance, and as a result, the circuit / control simulator 110 can perform a highly accurate simulation. On the other hand, the electromagnetic field simulator 120 needs to analyze the magnetic fields when the armature currents slightly change. However, since the simulator 100 calculates the inductance L on condition that the current I or the angle θ exceeds a predetermined threshold, the calculation amount of the inductance L can be reduced while maintaining the calculation accuracy of the inductance L.

  As described above, according to the simulator 100, the time required for the simulation can be reduced while accurately simulating the motor 180. Thereby, the development time of the motor 180 can be reduced. Further, even when the simulator 100 performs a simulation of the motor 180 for the purpose of checking the operation state of the motor 180 that is actually operated, the real-time property can be further improved.

  In addition to the three-phase synchronous motor, the motor 180 may be a rotating machine having a rotating movable part that is magnetically coupled to a circuit, such as a hysteresis motor, a stepping motor (pulse motor), or an induction motor. The simulator 100 can also control a generator, an electromagnetic brake, and the like that operate on substantially the same principle as these rotating machines. The simulator 100 can also control an electromagnetically driven device having a movable part magnetically coupled to a circuit including a solenoid valve (electromagnetic valve). Further, the simulator 100 can be a control target as long as it is an electromagnetic induction device regardless of the presence or absence of a movable part magnetically coupled to the circuit. For example, the simulator 100 can control a transformer.

  When the circuit / control simulator 110 simulates a system having n coils, a voltage equation to be solved by the circuit / control simulator 110 is expressed by the following formula (1).

  However, [V], [R], [I], and [Ψ] are the voltage applied to both ends of the coil, the resistance of the coil, the current flowing through the coil, and the interlinkage magnetic flux of the coil, respectively. 2). In the following formulas, parameters corresponding to the coils 1 to n are identified by subscripts 1 to n.

Here, the flux linkage [Ψ] of the coil is expressed as the sum of the contribution [Ψ _cur ] due to the current flowing through the coil and the contribution [Ψ _mag ] due to the magnet (including the electromagnet). That is, [Ψ] is represented by Expression (3).

  Therefore, the voltage equation to be solved by the circuit / control simulator 110 is expressed by Equation (4).

  Here, the inductance [L] of the coil is defined by Equation (5) in consideration of the non-linear effect of the material.

  Here, the inductance [L] is expressed by Equation (6). In Expression (6), it is assumed that the change in inductance [L] is negligible for a minute current change.

  Therefore, from Equations (1) and (5), the voltage equation to be solved by the circuit / control simulator 110 can be expressed by Equation (7).

The electromagnetic field simulator 120 calculates circuit parameters [L] and [R] and an induced electromotive force [V _i ] = [dΨ _mag / dt] for solving the voltage equation represented by Expression (7), It is passed to the circuit / control simulator 110. When the circuit / control simulator 110 simulates the motion of a movable part such as a motor or a valve using electromagnetic induction, the electromagnetic field simulator 120 is an electromagnetic force (parameter for solving the motion equation of the movable part). For example, torque or the like can be calculated and provided to the circuit / control simulator 110.

  FIG. 2 shows an example of a block configuration of the simulator 100. The circuit / control simulator 110 includes a current calculation unit 200, a time step determination unit 210, and a position calculation unit 220. The electromagnetic field simulator 120 includes a current change amount determination unit 230, a time change amount determination unit 240, a position change amount determination unit 250, a magnetic field analysis control unit 260, an output parameter storage unit 270, a parameter output unit 280, and a magnetic field analysis unit 290. Have.

  The simulator 100 simulates a circuit including an inductance element. The time step determining unit 210 determines a time step of simulation. And the current calculation part 200 calculates the electric current value sent through an inductance element for every time step. Specifically, the current calculation unit 200 calculates a current value for each time step based on Expression (7). The time step determination unit 210 may determine the time step of the simulation by determining the time width from the time indicated by the previous time step according to the accuracy of the simulation.

  Then, the current change amount determination unit 230 determines in advance a current change amount indicating a difference between the current value at the predetermined time step calculated by the current calculation unit 200 and the current value at the time step before the predetermined time step. It is determined whether or not the current change amount is larger. The magnetic field analysis unit 290 is based on the current value at a predetermined time step calculated by the current calculation unit 200 when the current change amount determination unit 230 determines that the current change amount is larger than a predetermined current change amount. Thus, the inductance of the inductance element at a predetermined time step is calculated. The output parameter storage unit 270 stores the inductance at a predetermined time step calculated by the magnetic field analysis unit 290. Then, the parameter output unit 280 outputs the inductance at a predetermined time step stored in the output parameter storage unit 270 to the current calculation unit 200 for each time step. Then, the current calculation unit 200 uses the inductance output from the parameter output unit 280 to calculate a current value that flows through the inductance element in a time step subsequent to a predetermined time step.

  Further, when the current change amount determination unit 230 determines that the current change amount is equal to or less than a predetermined current change amount, the magnetic field analysis control unit 260 causes the magnetic field analysis unit 290 to perform a predetermined time step of the inductance element. The inductance at the previous time step stored in the output parameter storage unit 270 is output from the parameter output unit 280 to the current calculation unit 200 without calculating the inductance at. For this reason, the simulator 100 can reduce the time required for the simulation. The previous time step may be the previous time step when the magnetic field analysis unit 290 calculates the inductance.

  The time change amount determination unit 240 determines whether or not the time change amount indicating the difference between the time indicating the predetermined time step and the time indicating the previous time step is larger than a predetermined time change amount. To do. The magnetic field analysis unit 290 determines that the time change amount determination unit 240 determines that the time change amount is greater than the predetermined time change amount, and the current change amount determination unit 230 determines the current change amount in advance. When it is determined that the current change amount is larger than the calculated current change amount, the inductance of the inductance element at the predetermined time step is calculated based on the current value at the predetermined time step. The output parameter storage unit 270 stores the inductance at a predetermined time step calculated by the magnetic field analysis unit 290, and the parameter output unit 280 stores the inductance stored in the output parameter storage unit 270 in the current calculation unit 200. Output at every time step.

  When the time variation determination unit 240 determines that the time variation is equal to or less than a predetermined time variation, the magnetic field analysis control unit 260 causes the magnetic field analysis unit 290 to perform a predetermined time of the inductance element. The inductance at the previous time step stored in the output parameter storage unit 270 is output from the parameter output unit 280 to the current calculation unit 200 without calculating the inductance in the step.

  When the simulator 100 simulates a circuit including a plurality of inductance elements that are magnetically coupled to each other, the current calculation unit 200 calculates a current value that flows through each of the plurality of inductance elements that are magnetically coupled to each other. Calculate every time step. Then, the current change amount determination unit 230 determines whether or not the maximum value among the current change amounts in each of the plurality of inductance elements is larger than a predetermined current change amount. Then, when the current change amount determination unit 230 determines that the maximum value of the current change amount is larger than the predetermined current change amount, the magnetic field analysis unit 290 determines each of the plurality of inductance elements in a predetermined time step. Based on the current value, the inductance of each of the plurality of inductance elements in a predetermined time step is calculated. Then, the output parameter storage unit 270 stores the inductances of the plurality of inductance elements in the predetermined time step calculated by the magnetic field analysis unit 290.

  Then, the parameter output unit 280 outputs the inductances of the plurality of inductance elements stored in the output parameter storage unit 270 to the current calculation unit 200 for each time step. Then, when the current change amount determination unit 230 determines that the maximum value of the current change amount is equal to or less than the current change amount, the magnetic field analysis control unit 260 causes the magnetic field analysis unit 290 to perform a plurality of times at a predetermined time step. The inductances of the plurality of inductance elements in the previous time step stored in the output parameter storage unit 270 are output from the parameter output unit 280 to the current calculation unit 200 without calculating the inductances of the inductance elements.

  Note that the current change amount determination unit 230 has a total current change amount indicating a difference between the total value of the current values passed through each of the plurality of inductance elements in a predetermined time step and the total value of the current values in the previous time step. It may be determined whether or not the current change amount is larger than a predetermined amount. Then, when the current change amount determination unit 230 determines that the total current change amount is larger than the predetermined total current change amount, the magnetic field analysis unit 290 determines each of the plurality of inductance elements in a predetermined time step. Based on the current value, the inductance of each of the plurality of inductance elements in a predetermined time step is calculated. The output parameter storage unit 270 stores the inductances of the plurality of inductance elements in the predetermined time step calculated by the magnetic field analysis unit 290, and the parameter output unit 280 stores the predetermined parameter stored in the output parameter storage unit 270. The inductances of the plurality of inductance elements in the time step are output to the current calculation unit 200 for each time step.

  When the current change amount determination unit 230 determines that the total current change amount is equal to or less than the predetermined total current change amount, the magnetic field analysis control unit 260 causes the magnetic field analysis unit 290 to perform a predetermined time step. The respective inductances of the plurality of inductance elements in the previous time step stored in the output parameter storage unit 270 are output from the parameter output unit 280 to the current calculation unit 200 without calculating the inductances of the plurality of inductance elements in FIG. Let For example, when a three-phase armature circuit is arranged so as to share one magnetic circuit, for example, as in a three-phase permanent magnet synchronous motor, the current value of a specific one phase has changed considerably. However, the interlinkage magnetic flux by all three phases may not change so much. In such a case, the electromagnetic field simulator 120 can more appropriately determine whether or not to calculate the inductance by determining whether or not the total value of the current values flowing through the three phases has changed.

  The position calculating unit 220 calculates the position of the movable unit that is magnetically coupled to the inductance element and is driven by the magnetic field generated by the inductance element for each time step. Specifically, the position calculation unit 220 determines the angle of the rotor by solving an equation of motion related to the movable part. Then, the position change amount determination unit 250 determines in advance a position change amount indicating a difference between the position of the movable unit at the predetermined time step calculated by the position calculation unit 220 and the position of the movable unit at the previous time step. It is determined whether or not the position change amount is larger. Then, the magnetic field analysis unit 290 determines that the position change amount determination unit 250 determines that the position change amount is larger than a predetermined position change amount, or the current change amount determination unit 230 determines that the current change amount is in advance. When it is determined that the current change amount is larger than the predetermined current change amount, a predetermined time is determined based on the current value calculated by the current calculation unit 200 in the predetermined time step and the position calculated by the position calculation unit 220 in the predetermined time step. The inductance of the inductance element in the step is calculated. The output parameter storage unit 270 stores the inductance at a predetermined time step calculated by the magnetic field analysis unit 290, and the parameter output unit 280 stores the inductance stored in the output parameter storage unit 270 in the current calculation unit 200. Output at every time step.

  Then, the magnetic field analysis control unit 260 determines that the position change amount determination unit 250 is less than or equal to a predetermined position change amount, and the current change amount determination unit 230 determines that the current change amount is When it is determined that the current change amount is equal to or less than the predetermined current change amount, the magnetic field analysis unit 290 does not calculate the inductance at a predetermined time step of the inductance element, and the time before the output parameter storage unit 270 stores The inductance in the step is output from the parameter output unit 280 to the current calculation unit 200. As described above, the electromagnetic field simulator 120 can calculate the inductance as necessary according to the position of the rotor of the motor 180.

  The magnetic field analysis unit 290 determines the current value at the predetermined time step calculated by the current calculation unit 200 when the time change amount determination unit 240 determines that the time change amount is larger than a predetermined time change amount. And based on the position in the predetermined time step calculated by the position calculation unit 220, the induced electromotive force by the movable part in the inductance element in the predetermined time step is calculated. The output parameter storage unit 270 stores the induced electromotive force generated by the movable unit at a predetermined time step calculated by the magnetic field analysis unit 290. When the time variation determination unit 240 determines that the time variation is equal to or less than a predetermined time variation, the magnetic field analysis control unit 260 causes the magnetic field analysis unit 290 to perform a predetermined time of the inductance element. The parameter output unit 280 outputs the induced electromotive force generated by the movable unit stored in the output parameter storage unit 270 to the current calculation unit 200 for each time step without calculating the induced electromotive force generated by the movable unit in the step. Then, the current calculation unit 200 calculates a current value to be passed through the inductance element in the next time step based on the induced electromotive force and the inductance output from the parameter output unit 280.

  The magnetic field analysis unit 290 determines that the time change amount determination unit 240 determines that the time change amount is greater than a predetermined time change amount, or the position change amount determination unit 250 determines that the position change amount is When it is determined that the position change amount is larger than the predetermined position change amount, a predetermined time is determined based on the current value calculated by the current calculation unit 200 in the predetermined time step and the position calculated by the position calculation unit 220 in the predetermined time step. You may calculate the induced electromotive force by the movable part in the inductance element in a step. Thereby, the electromagnetic field simulator 120 can appropriately determine whether or not to calculate the induced electromotive force according to the position of the rotor. The details of the operation of calculating the induced electromotive force by the magnetic field analysis unit 290 will be described with reference to FIG.

  In the magnetic field analysis unit 290, the current value and position calculation unit 220 in the predetermined time step calculated by the current calculation unit 200 when the time change amount determination unit 240 is larger than the predetermined time change amount. Based on the calculated position at the predetermined time step, the electromagnetic force of the movable part at the predetermined time step is calculated. Then, the output parameter storage unit 270 stores the electromagnetic force of the movable unit at a predetermined time step calculated by the magnetic field analysis unit 290. The magnetic field analysis control unit 260 allows the magnetic field analysis unit 290 to move at a predetermined time step when the time change amount determination unit 240 determines that the time change amount is equal to or less than a predetermined time change amount. The parameter output unit 280 outputs the electromagnetic force of the movable part stored in the output parameter storage unit 270 to the position calculation unit 220 for each time step without calculating the electromagnetic force of the unit. The position calculation unit 220 calculates the position of the movable unit in the next time step based on the electromagnetic force output by the parameter output unit 280.

  As described above, the simulator 100 can appropriately determine whether to calculate the inductance, the induced electromotive force, and the torque according to the armature current of the motor 180, the angle of the rotor, and the time. . Thereby, the simulator 100 can reduce the simulation time while maintaining the accuracy of the simulation.

  FIG. 3 shows an example of a block configuration of the magnetic field analysis unit 290. The magnetic field analysis unit 290 includes an interlinkage magnetic flux calculation unit 300, an inductance calculation unit 302, a linear region determination unit 340, an inductance storage unit 350, an inductance calculation control unit 360, an induced electromotive force calculation unit 380, and an electromagnetic force calculation unit 390. . The inductance calculation unit 302 includes a first interlinkage magnetic flux calculation unit 311, a first self-inductance calculation unit 321, a first mutual inductance calculation unit 331, a second interlinkage magnetic flux calculation unit 312, a second self-inductance calculation unit 322, and 2 mutual inductance calculation part 332 is included.

  Based on the current value of the inductance element in the predetermined time step calculated by the current calculation unit 200, the flux linkage calculation unit 300 calculates the flux linkage of the inductance element based on the current indicated by the current value. Then, the inductance calculation unit 302 calculates the interlinkage magnetic flux of the inductance element based on the current indicated by the current value changed by a minute current value from the current value at the predetermined time step calculated by the current calculation unit 200, and the calculated interlinkage magnetic flux The inductance of the inductance element is calculated based on a value obtained by dividing the difference between the flux linkage calculated by the flux linkage calculation unit 300 by the minute current value. Specifically, the inductance calculation unit 302 calculates the inductance based on Expression (5). As a result, the electromagnetic field simulator 120 can include non-linearity with respect to the current value in the inductance, so that the circuit / control simulator 110 can accurately calculate the current flowing through the circuit.

  The induced electromotive force calculation unit 380 is based on the current value calculated in the predetermined time step by the current calculation unit 200 and the position in the predetermined time step calculated by the position calculation unit 220, and the movable part in the inductance element in the predetermined time step. Calculate the induced electromotive force. Specifically, the flux linkage calculation unit 300 is configured to determine the value of the current flowing through the plurality of inductance elements at the predetermined time step calculated by the current calculation unit 200 and the position of the movable unit at the predetermined time step calculated by the position calculation unit 220. Based on the above, the current indicated by the current value and the flux linkage of the inductance element by the movable part are calculated.

  Then, the inductance calculation unit 302 is a current value when the current value is changed by a minute current value from the current value at the predetermined time step calculated by the current calculation unit 200 and the position of the movable unit at the predetermined time step calculated by the position calculation unit 220. Based on the current indicated by the current value changed by a minute current value and the linkage flux of the inductance element by the movable part, and the difference between the calculated linkage flux and the linkage flux calculated by the linkage flux calculation unit 300 The inductance of the inductance element is calculated based on the value obtained by dividing the value by the minute current value. The induced electromotive force calculation unit 380 has a magnetic flux change value obtained by subtracting the interlinkage magnetic flux in the previous time step from the interlinkage magnetic flux in the predetermined time step calculated by the interlinkage magnetic flux calculation unit 300, and the inductance calculation unit 302 has the predetermined value. The inductance calculated in the time step, the current change value obtained by subtracting the current value in the previous time step from the current value in the predetermined time step calculated by the current calculation unit 200, and the time indicated by the predetermined time step and the previous time step Based on the difference from the indicated time, an induced electromotive force by the movable part is calculated. As described above, the electromagnetic field simulator 120 calculates the nonlinearity with respect to the current value by including the inductance in the inductance, and subtracts the change in the linkage flux due to the inductance from the change in the total linkage magnetic flux, and calculates the remainder of the linkage by the movable part. Calculated as the amount of change in magnetic flux. Since the electromagnetic field simulator 120 calculates the induced electromotive force by the movable part from the change in the interlinkage magnetic flux by the movable part, the circuit / control simulator 110 can simulate the circuit with high accuracy.

  The case where the current calculation unit 200 calculates the current value flowing through each of the plurality of inductance elements magnetically coupled to each other for each time step will be described. The linkage flux calculation unit 300 is calculated by the current calculation unit 200. The interlinkage magnetic flux of each of the first inductance element and the second inductance element is calculated from the respective current values of the first inductance element and the second inductance element at the predetermined time step. Then, the first flux linkage calculation unit 311 includes the first inductance element and the first inductance element based on the current indicated by the current value changed by a minute current value from the current value of the first inductance element in the predetermined time step calculated by the current calculation unit 200. The interlinkage magnetic flux of each second inductance element is calculated. The first self-inductance calculation unit 321 calculates the linkage flux between the first inductance element calculated by the first linkage flux calculation unit 311 and the linkage flux of the first inductance element calculated by the linkage flux calculation unit 300. Based on the value obtained by dividing the difference by the minute current value, the self-inductance of the first inductance element is calculated. Then, the first mutual inductance calculation unit 331 calculates the relationship between the linkage flux of the second inductance element calculated by the first linkage flux calculation unit 311 and the linkage flux of the second inductance element calculated by the linkage flux calculation unit 300. Based on the value obtained by dividing the difference by the minute current value, the mutual inductance of the second inductance element with the first inductance element is calculated.

  In addition, the second flux linkage calculating unit 312 includes the first inductance element and the first inductance element based on the current indicated by the current value changed by a minute current value from the current value of the second inductance element in the predetermined time step calculated by the current calculation unit 200. The interlinkage magnetic flux of each second inductance element is calculated. Then, the second self-inductance calculating unit 322 calculates the link between the flux linkage of the second inductance element calculated by the second flux linkage calculating unit 312 and the flux linkage of the second inductance element calculated by the flux linkage calculating unit 300. Based on the value obtained by dividing the difference by the minute current value, the self-inductance of the second inductance element is calculated. Then, the second mutual inductance calculation unit 332 calculates the relationship between the linkage flux of the first inductance element calculated by the second linkage flux calculation unit 312 and the linkage flux of the first inductance element calculated by the linkage flux calculation unit 300. Based on the value obtained by dividing the difference by the minute current value, the mutual inductance of the first inductance element with the second inductance element is calculated. Specifically, the first self-inductance calculator 321, the first mutual inductance calculator 331, the second self-inductance calculator 322, and the second mutual inductance calculator 332 calculate the inductance based on the equation (5). To do. Thereby, the inductance calculation part 302 can calculate the inductance including the nonlinearity with respect to the current value. Moreover, the inductance calculation part 302 can calculate appropriately the inductance of each coil in the apparatus in which several coils exist.

  The inductance storage unit 350 has a mutual relationship between the self-inductance of the first inductance element in the time step before the time calculated by the first self-inductance calculation unit 321 and the second inductance element in the time step before the time calculated by the first mutual inductance calculation unit 331. The inductance, the self-inductance of the second inductance element in the previous time step calculated by the second self-inductance calculating unit 322, and the mutual inductance of the first inductance element in the previous time step calculated by the second mutual inductance calculating unit 332 are calculated. Remember. Then, the linear region determination unit 340 includes the first inductance element at the time step before the self-inductance of the first inductance element at the predetermined time step calculated by the first self-inductance calculation unit 321 and the inductance storage unit 350 stores the first inductance element. And the inductance storage unit 350 stores the mutual inductance of the second inductance element in the predetermined time step calculated by the first mutual inductance calculation unit 331 and the difference from the self-inductance of the second inductance element is smaller than a predetermined first value. When the difference from the mutual inductance of the second inductance element in the previous time step is smaller than a predetermined second value, the flux linkage of the first inductance element is relative to the current value of the first inductance element. Determine that it is in the linear region

  Then, when the linear region determination unit 340 determines that the linear region determination unit 340 is not in the linear region, the inductance calculation control unit 360 transmits the self inductance of the second inductance element to the second self inductance calculation unit 322 and the second mutual inductance calculation unit 332. When the mutual inductance of the first inductance element is calculated and the linear region determining unit 340 determines that the first inductance element is in the linear region, the second self-inductance calculating unit 322 and the second mutual inductance calculating unit 332 have the second inductance element The self-inductance and the mutual inductance of the first inductance element are not calculated, and the self-inductance of the second inductance element and the mutual inductance of the first inductance element in the time step before the inductance storage unit 350 stores them. And outputs the Nsu. As described above, when the inductance in the linear region is calculated for one armature circuit, the electromagnetic field simulator 120 calculates the inductance calculated last time without calculating the inductance of the other armature circuit. To provide. Thereby, the electromagnetic field simulator 120 can reduce simulation time.

  The magnetic field analysis unit 260 determines that the position change amount is greater than the predetermined position change amount on the condition that the time change amount determination unit 240 determines that the time change amount is larger than the predetermined time change amount. When the position change amount determination unit 250 determines that the current change amount is large, or when the current change amount determination unit 230 determines that the current change amount is larger than a predetermined current change amount, the current calculation unit 200 calculates the current change amount. Based on the current value at the predetermined time step and the position at the predetermined time step calculated by the position calculation unit 220, the inductance of the inductance element at the predetermined time step is calculated. The magnetic field analysis control unit 260 then determines that the time change amount determination unit 240 determines that the time change amount is equal to or less than the predetermined time change amount, or the position change amount is equal to or less than the predetermined position change amount. If the position change amount determination unit 250 determines that the current change amount is equal to or less than the predetermined current change amount, the magnetic field analysis unit 260 The inductance at a predetermined time step of the inductance element is not calculated, and the inductance at the previous time step stored in the output parameter storage unit 270 is output from the parameter output unit 280 to the current calculation unit 200.

  In the above description, the operation in which the electromagnetic field simulator 120 calculates the inductance of the two inductance elements that are magnetically coupled to each other has been described. However, the electromagnetic field simulator 120 uses the same method to calculate the inductance of two or more inductance elements. Needless to say, can be calculated. The inductance element in the present invention means a component that generates electromagnetic induction. Therefore, the coil is an example of the inductance element in the present invention, and the circuit itself in which electromagnetic induction occurs can be the inductance element in the present invention.

FIG. 4 shows a case where the circuit parameters (inductance [L] and induced electromotive force [V _i ]) and torque T are calculated by the magnetic field analysis unit 290 when the simulator 100 simulates the operation of the n-phase permanent magnet synchronous motor. An example of the operation flow is shown. The following description is an example of an operation flow when the circuit parameters and torque are calculated at time step N. It is assumed that the magnetic field analysis unit 290 calculates the circuit parameters at time step M, which is a time step before time step N. Note that the times indicated by the time steps N and M are t _N and t _M , respectively.

At time step N, the magnetic field analysis unit 290 outputs the current value I _m (t _N ) (where m = 1 to n) and the angle of the rotor from the circuit / control simulator 110 to the time t _N. Obtain θ (t _N ). Then, the time change amount determination unit 240 calculates a difference Δt = t _N −t _M between the time t _M and the time t _N indicated by the time step M when the magnetic field analysis unit 290 previously calculated the circuit parameters (S402). Then, the time change amount determination unit 240 determines whether or not Δt is equal to or less than a predetermined time change amount Δt ₀ (S404).

When the time change amount determination unit 240 determines that Δt is greater than Δt ₀ in S404, the flux linkage calculation unit 300 calculates the flux linkage [Ψ] linked to each armature circuit (S406). . At this time, the flux linkage calculation unit 300 calculates the magnetic flux density distribution in the motor based on the current value I _m (t _N ) of each armature circuit and the rotor angle θ (t _N ). Based on the magnetic flux density distribution, an interlinkage magnetic flux interlinking with each armature circuit can be calculated.

Then, the current change amount determination unit 230 calculates the difference ΔI _m = | I _m (t _M ) −I _m (t _N ) | from the current value at the time step M for each phase, and the maximum value ΔI _max is calculated as a current change amount ΔI. Then, the current change amount determination unit 230 determines whether ΔI is equal to or less than a predetermined current change amount ΔI ₀ (S408). The current change amount determination unit 230 may determine whether ΔI _m in each phase is equal to or less than ΔI ₀ . ΔI ₀ may be determined for each phase. For example, a value obtained by multiplying I _m (t _N ) by a predetermined constant may be determined as ΔI ₀ of each phase. For example, ΔI ₀ = I _m (t _N ) × 0.3 may be set. In addition, the current change amount determination unit 230 may determine whether or not the current change amount of the phase having the maximum current value in each phase is equal to or less than ΔI ₀ . Further, the current change amount determination unit 230 determines the difference between the total amount of current values for each phase and the total amount of current values at the time step M, that is, | Σ (I _m (t _M ) −I _m (t _N )) | (where m = 1 to n) may be determined as to whether ΔI ₀ or less.

In step S <b> 408, the position change determination unit 250 calculates the difference between the rotor positions θ (t _M ) and θ (tN) at the time step M as the rotor position change Δθ. Then, the position change amount determination unit 250 determines whether or not Δθ is equal to or less than a predetermined position change amount Δθ ₀ . In step S408, the magnetic field analysis control unit 260 determines whether ΔI is ΔI ₀ or less and Δθ is Δθ ₀ or less. Δθ ₀ may be determined in advance according to the number of poles of the rotor.

In S <b> 408, when ΔI is greater than ΔI ₀ or Δθ is greater than Δθ ₀ , the magnetic field analysis control unit 260 causes the magnetic field analysis unit 290 to calculate an inductance. Specifically, the first flux linkage calculation unit 311 increases the current value of the first armature circuit in which the current of I ₁ = I ₁ (t _N ) flows by the minute current change amount ΔI ₁ . In this case, the flux linkage [Ψ ′] of each armature circuit is calculated (S410). The first flux linkage calculation unit 311 has a current value of the first armature circuit of I ₁ + ΔI ₁ , and current values of the other armature circuits of I _m = I _m (t _N ) (m = 2 to 2). n) Calculate the flux linkage [Ψ ′] when the rotor angle is θ (t _N ). For example, the first flux linkage calculation unit 311 sets the magnitude of ΔI _{1 to} a magnitude of several percent of the magnitude of I ₁ . For example, the first flux linkage calculation unit 311, the magnitude of [Delta] I ₁ and 5% of the size of _{I 1 (ΔI 1 = I 1} × 0.05). The first flux linkage calculation unit 311, if _{I 1} + [Delta] I reference current value magnitude predetermined _one (for example, 5% of the magnitude of the rated current _{I r)} is smaller than, _{I 1} + [Delta] I ₁ is the reference current value (e.g., _{I r} × 0.05) may determine the [Delta] I ₁ to be. Thus, the minute current change amount does not necessarily mean a minute current change amount with respect to the absolute value of the current supplied to the motor, and is smaller than the standard magnetic field strength generated in the system. It means a minute change in current that gives a minute change. By providing the reference current value as described above, it is possible to prevent ΔI ₁ from becoming a value that is so small that the calculation accuracy of the numerical calculation is deteriorated. Note that the reference current value is not only defined for the rated current as described above, but may be appropriately set by other methods depending on the system. The first self-inductance calculating unit 321 calculates the self-inductance _{L 11,} the first mutual inductance calculation unit 331 calculates the mutual inductance _L m1 (m = _2~n) of other phases (S412).

Incidentally, the inductance calculation of in S412 includes a flux linkage calculated in S406, the flux linkage calculated in S410, based on the amount of change [Delta] I ₁ of the current value can be performed as follows. First, the flux linkage [Ψ] calculated in S406 is expressed by Expression (8).

  Further, the flux linkage [Ψ ′] calculated in S410 is expressed by Expression (9).

The first self-inductance calculating unit 321 can calculate the self-inductance L ₁₁ by the equation (10) derived from the equations (9) and (8), and the first mutual inductance calculating unit 331 can calculate the mutual inductance L _21. ..L _n1 can be calculated.

Then, the linear region determining unit 340 determines whether L ₁₁ , L ₂₁ ... L _n1 are in the linear region (S414). A method for determining whether or not the linear region determination unit 340 is in the linear region will be described separately. In _S414, if _{L 11,} L 21 ··· L _n1 is that not in the linear region is determined by a method similar to the method of calculating the inductance _{L 11, L} 21 ··· _{L n1} in S410 and S412, The inductance of each armature circuit can be calculated. For example, by calculating the interlinkage magnetic fluxes Ψ ₁ ^{′ to} Ψ n ^′ when the current value of the m-th armature circuit is changed by ΔI _m , the inductances L _1m , L _2m ,. , L _nm can be calculated.

The inductance storage unit 350 stores the inductance [L] calculated by the inductance calculation unit 302 (S440). Then, the induced electromotive force calculation unit 380 calculates the induced electromotive force generated in each armature circuit by the rotor (S450). Specifically, the induced electromotive force calculation unit 380 calculates a change amount Δ [Ψ _mag ] of the linkage flux of each armature circuit by the rotor by the following formula (12) based on formulas (3) and (5). Calculate from

Then, the induced electromotive force calculation unit 380 calculates induced electromotive force [V _i ] = Δ [Ψ _mag ] / Δt based on Δt with respect to the previous time step. Note that, as described above, the inductance calculation unit 302 calculates the influence of nonlinearity in the inductance. Then, the induced electromotive force calculation unit 380 calculates the value obtained by subtracting the magnetic flux expressed by the inductances [L] and [I] from the total interlinkage magnetic flux as represented by the equation (12). And the induced electromotive force [V _i ] is calculated based on this. For this reason, the circuit / control simulator 110 can appropriately perform circuit simulation in which nonlinearity of inductance is incorporated.

Further, the electromagnetic force calculation unit 390 calculates the torque T based on the magnetic flux density distribution calculated in S406 (S452). Then, the output parameter storage unit 270 includes the induced electromotive force [V _i ] calculated by the induced electromotive force calculation unit 380, the torque T calculated by the electromagnetic force calculation unit 390, and the inductance [L] calculated by the inductance calculation unit 302. Store (S454). Then, the parameter output unit 280 outputs the inductance [L], the induced electromotive force [V _i ], and the torque T stored in the output parameter storage unit 270 to the circuit / control simulator 110 (S456).

If Δt is equal to or smaller than Δt ₀ in S404, the magnetic field analysis control unit 260 shifts the process to S456, and the inductance [L] at time step M stored in the output parameter storage unit 270 in S456. The induced electromotive force [V _i ] and the torque T are output to the circuit / control simulator 110. As a result, the electromagnetic field simulator 120 does not need to perform magnetic field analysis such as calculation of magnetic flux, so that the amount of calculation required for the simulation can be reduced.

If the magnetic field analysis control unit 260 determines that ΔI is ΔI ₀ or less and Δθ is Δθ ₀ or less in S408, the process proceeds to S450, and the induced electromotive force calculation unit 380 performs the process in S450. Then, the induced electromotive force [V _i ] is calculated. At this time, the induced electromotive force is based on the inductance [L] and the linkage flux [Ψ] at the time step M stored in the parameter output unit 280 and the linkage flux [Ψ] at the time step N calculated in S406. [V _i ] is calculated. Thus, the electromagnetic field simulator 120 can reduce the calculation amount required for the simulation by appropriately skipping the inductance calculation.

In S414, when it is determined that L ₁₁ , L ₂₁ ... L _n1 are in the linear region, the process proceeds to S450, and the induced electromotive force calculation unit 380 calculates the induced electromotive force [V _i ]. . At this time, the induced electromotive force calculation unit 380 calculates the induced electromotive force [V _i ] using the inductance at time step M stored in S440. As a result, the electromagnetic field simulator 120 can skip the calculation of the remaining inductance only by calculating one inductance of each armature circuit, so that the amount of calculation required for the simulation can be reduced.

Incidentally, the linear region determination unit 340, in _S414, the _{L 11,} L 21 ··· L _n1, and _{L 11, L} 21 ··· _{L n1} at time step M before inductance storage unit 350 stores Each difference may be calculated, and if each difference is equal to or less than a predetermined value, it may be determined that the difference is in the linear region. When the change in inductance from the time step M is small, it is expected that it does not deviate from the linear region. In such a case, the inductance at the time step M is used without calculating the inductance. As a result, the circuit / control simulator 110 can simulate the circuit with sufficient accuracy.

In addition, when comparing the inductance calculated in S412 with the inductance in the previous time step, the linear region determination unit 340 does not compare with the inductance in the time step M but one or more in the previous time step. You may compare with an inductance. Specifically, in S440, the inductance storage unit 350 associates each of the inductances L _m1 (m = 1 to n) at a plurality of time steps with the current values I _m (m = 1 to n) at the time steps. And remember. That is, the inductance storage unit 350 stores the value of _{L m1} in _{I m} space. Then, the linear region determination unit 340 identifies an I _m region where L _m1 is substantially constant, determines whether I _m (t _N ) is included in the I _m region, and is included in the region. On the condition, it may be determined that L _m1 is in the linear region. The inductance storage unit 350 may store the inductance L _m1 (m = 1 to n) in further association with θ, and the linear region determination unit 340 has a region of I _m and θ where L _m1 is substantially constant. If I _m (t _N ) and θ (t _N ) are included in the region, it may be determined that L _m1 is in the linear region. Note that the linear region determination unit 340 can determine the linearity of [Ψ] calculated in S406 with respect to I _m (or I _m and θ) instead of L _m1 . Note that, when the linear region determination unit 340 determines the linearity described above, it may be determined based on [Ψ] calculated in S406 which circuit to determine based on the inductance. That is, the inductance calculating unit 302 may specify a circuit having the largest flux linkage among [Ψ] calculated in S406, and calculate the inductance in S412 for the specified circuit. The linear region determination unit 340 may determine the linearity based on the inductance of the specified circuit.

  Note that the storage function of the output parameter storage unit 270 and the inductance storage unit 350 in S440 and S454 may be included in the memory of the simulator 100. Note that “store” may be synonymous with the memory of the simulator 100 storing values such as inductance.

  The armature circuit is an example of a circuit including the inductance element in the present invention. Further, the rotor angle θ is an example of the position of the movable part in the present invention, and in another example, the position of the movable part may be expressed by three-dimensional coordinates (x, y, z). The position of the movable part may be a relative position of the movable part with respect to a circuit to be simulated. Needless to say, when the simulator 100 simulates a device such as a transformer that has substantially no movable part, the circuit / control simulator 110 does not need to calculate the position of the movable part. In the above description, the operation flow in which the electromagnetic field simulator 120 calculates the resistance [R], which is one of the circuit parameters, is omitted so as not to complicate the description. However, in the electromagnetic field simulator 120, the resistance [R] of the circuit changes according to the frequency or the change frequency of the current or magnetic field (for example, due to the skin effect) or the like. In the case of simulating a device to be operated, the resistance [R] may be calculated according to the position or the frequency of change of current or magnetic field, and the temperature of the circuit is calculated based on the thermal coupling calculation. The resistance [R] corresponding to the calculated temperature may be calculated. Then, the electromagnetic field simulator 120 may output the calculated resistance [R] to the circuit / control simulator 110.

  As described above, the simulator 100 can provide a simulation environment that can reduce the simulation time while accurately simulating the electromagnetic drive system 150.

  FIG. 5 shows an example of the hardware configuration of the simulator 100. The simulator 100 includes a CPU peripheral unit having a CPU 1505, a RAM 1520, a graphic controller 1575, and a display device 1580 connected to each other by a host controller 1582, and a communication interface 1530 connected to the host controller 1582 by an input / output controller 1584. An input / output unit having a hard disk drive 1540 and a CD-ROM drive 1560, and a legacy input / output unit having a ROM 1510, a flexible disk drive 1550, and an input / output chip 1570 connected to the input / output controller 1584.

  The host controller 1582 connects the RAM 1520, the CPU 1505 that accesses the RAM 1520 at a high transfer rate, and the graphic controller 1575. The CPU 1505 operates based on programs stored in the ROM 1510 and the RAM 1520 and controls each unit. The graphic controller 1575 acquires image data generated by the CPU 1505 and the like on a frame buffer provided in the RAM 1520 and displays the image data on the display device 1580. Alternatively, the graphic controller 1575 may include a frame buffer that stores image data generated by the CPU 1505 or the like.

  The input / output controller 1584 connects the host controller 1582 to the hard disk drive 1540, the communication interface 1530, and the CD-ROM drive 1560, which are relatively high-speed input / output devices. The hard disk drive 1540 stores programs and data used by the CPU 1505. The communication interface 1530 is connected to the network communication device 1598 to transmit / receive programs or data. The CD-ROM drive 1560 reads a program or data from the CD-ROM 1595 and provides it to the hard disk drive 1540 and the communication interface 1530 via the RAM 1520.

  The input / output controller 1584 is connected to the ROM 1510, the flexible disk drive 1550, and the relatively low-speed input / output device of the input / output chip 1570. The ROM 1510 stores a boot program that the simulator 100 executes at startup, a program that depends on the hardware of the simulator 100, and the like. The flexible disk drive 1550 reads a program or data from the flexible disk 1590 and provides it to the hard disk drive 1540 and the communication interface 1530 via the RAM 1520. The input / output chip 1570 connects various input / output devices via a flexible disk drive 1550 and, for example, a parallel port, a serial port, a keyboard port, a mouse port, and the like.

  A program executed by the CPU 1505 is stored in a recording medium such as the flexible disk 1590, the CD-ROM 1595, or an IC card and provided by the user. The program stored in the recording medium may be compressed or uncompressed. The program is installed in the hard disk drive 1540 from the recording medium, read into the RAM 1520, and executed by the CPU 1505.

  The program executed by the CPU 1505 causes the simulator 100 to function as the circuit / control simulator 110 and the electromagnetic field simulator 120 described with reference to FIGS. That is, the program causes the circuit / control simulator 110 to function as the current calculation unit 200, the time step determination unit 210, and the position calculation unit 220 described with reference to FIGS. In addition, the program converts the electromagnetic field simulator 120 into the current change amount determination unit 230, the time change amount determination unit 240, the position change amount determination unit 250, the magnetic field analysis control unit 260, which have been described with reference to FIGS. The output parameter storage unit 270, the parameter output unit 280, and the magnetic field analysis unit 290 function. The program also causes the magnetic field analysis unit 290 to perform the interlinkage magnetic flux calculation unit 300, the inductance calculation unit 302, the linear region determination unit 340, the inductance storage unit 350, and the inductance calculation control unit described with reference to FIGS. 360, the induced electromotive force calculation unit 380, and the electromagnetic force calculation unit 390, and the inductance calculation unit 302 is replaced with the first interlinkage magnetic flux calculation unit 311 and the first self-inductance calculation unit described with reference to FIGS. 321, a first mutual inductance calculator 331, a second flux linkage calculator 312, a second self-inductance calculator 322, and a second mutual inductance calculator 332.

  The program shown above may be stored in an external storage medium. As the storage medium, in addition to the flexible disk 1590 and the CD-ROM 1595, an optical recording medium such as a DVD or PD, a magneto-optical recording medium such as an MD, a tape medium, a semiconductor memory such as an IC card, or the like can be used. Further, a storage device such as a hard disk or RAM provided in a server system connected to a dedicated communication network or the Internet may be used as a recording medium, and the program may be provided to the simulator 100 via the network. The program may be provided to the simulator 100 in the form of a library, or may be provided to the simulator 100 in a form that can be linked with other programs statically or dynamically. A program for causing the electromagnetic field simulator 120 to function as the current change amount determination unit 230, the time change amount determination unit 240, the position change amount determination unit 250, the magnetic field analysis control unit 260, the output parameter storage unit 270, and the parameter output unit 280 is The simulator 100 may be provided as a program different from the program that causes the electromagnetic field simulator 120 to function as the magnetic field analysis unit 290.

  As mentioned above, although this invention was demonstrated using embodiment, the technical scope of this invention is not limited to the range as described in the said embodiment. Various modifications or improvements can be added to the above embodiment. It is apparent from the scope of the claims that the embodiments added with such changes or improvements can be included in the technical scope of the present invention.

3 is a diagram illustrating an example of a control target of a simulator 100. FIG. 2 is a diagram illustrating an example of a block configuration of a simulator 100. FIG. 3 is a diagram illustrating an example of a block configuration of a magnetic field analysis unit 290. FIG. It is a figure which shows an example of the operation | movement flow of the simulator. 2 is a diagram illustrating an example of a hardware configuration of a simulator 100. FIG.

Explanation of symbols

100 simulator 110 circuit / control simulator 120 electromagnetic field simulator 150 electromagnetic drive system 160 driver 170 inverter 180 motor 190 equivalent circuit 200 current calculation unit 210 time step determination unit 220 position calculation unit 230 current change amount determination unit 240 time change amount determination unit 250 Position change determination unit 260 Magnetic field analysis control unit 270 Output parameter storage unit 280 Parameter output unit 290 Magnetic field analysis unit 300 Interlinkage magnetic flux calculation unit 302 Inductance calculation unit 311 First interlinkage magnetic flux calculation unit 321 First self-inductance calculation unit 331 1 mutual inductance calculation unit 312 second flux linkage calculation unit 322 second self-inductance calculation unit 332 second mutual inductance calculation unit 340 linear region determination unit 350 inductance storage unit 360 inductor Scan calculation control unit 380 induced electromotive force calculating portion 390 electromagnetic force calculation unit

Claims (20)

A simulator for simulating a circuit including an inductance element,
A current calculation unit for calculating a current value flowing through the inductance element for each time step;
Whether or not a current change amount indicating a difference between a current value at a predetermined time step calculated by the current calculation unit and a current value at a time step before the predetermined time step is larger than a predetermined current change amount A current change amount determination unit for determining whether or not
When the current change amount determination unit determines that the current change amount is larger than a predetermined current change amount, the predetermined amount of time is calculated based on the current value in the predetermined time step calculated by the current calculation unit. A magnetic field analyzer for calculating the inductance of the inductance element in the time step;
An output parameter storage unit for storing inductance in the predetermined time step calculated by the magnetic field analysis unit;
A parameter output unit that outputs the inductance at the predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step; and
When the current change amount determination unit determines that the current change amount is equal to or less than a predetermined current change amount, the magnetic field analysis unit does not calculate the inductance of the inductance element at the predetermined time step, and A simulator comprising: a magnetic field analysis control unit that outputs the inductance in the previous time step stored in the output parameter storage unit from the parameter output unit to the current calculation unit. A time change amount determination unit for determining whether or not a time change amount indicating a difference between the time indicating the predetermined time step and the time indicating the previous time step is larger than a predetermined time change amount; Prepared,
The magnetic field analysis unit determines that the time change amount determination unit determines that the time change amount is larger than a predetermined time change amount, and the current change amount determination unit determines the current change amount. When it is determined that the current change amount is larger, the inductance of the inductance element in the predetermined time step is calculated based on the current value in the predetermined time step,
The output parameter storage unit stores the inductance in the predetermined time step calculated by the magnetic field analysis unit,
The parameter output unit outputs the inductance stored in the output parameter storage unit to the current calculation unit for each time step,
When the time variation determination unit determines that the time variation is equal to or less than a predetermined time variation, the magnetic field analysis control unit causes the magnetic field analysis unit to perform the predetermined time step of the inductance element. 2. The simulator according to claim 1, wherein the inductance in the previous time step stored in the output parameter storage unit is output from the parameter output unit to the current calculation unit without calculating the inductance in step 1. The current calculation unit calculates a current value flowing through each of the plurality of inductance elements magnetically coupled to each other for each time step,
The current change amount determination unit determines whether a maximum value of current change amounts in each of the plurality of inductance elements is larger than a predetermined current change amount,
The magnetic field analysis unit, when the current change amount determination unit determines that the maximum value of the current change amount is larger than a predetermined current change amount, each of the plurality of inductance elements in the predetermined time step. On the basis of the current value of each of the plurality of inductance elements in the predetermined time step,
The output parameter storage unit stores inductances of the plurality of inductance elements in the predetermined time step calculated by the magnetic field analysis unit,
The parameter output unit outputs each inductance of the plurality of inductance elements stored in the output parameter storage unit to the current calculation unit for each time step,
When the current change amount determination unit determines that the maximum value of the current change amount is equal to or less than the current change amount, the magnetic field analysis control unit causes the magnetic field analysis unit to send the plurality of the current change amount determination units in the predetermined time step. The inductance of each of the plurality of inductance elements in the previous time step stored in the output parameter storage unit is output from the parameter output unit to the current calculation unit without calculating the inductance of each inductance element. The simulator according to claim 1. The current calculation unit calculates a current value flowing through each of the plurality of inductance elements magnetically coupled to each other for each time step,
The current change amount determination unit is a total current change amount indicating a difference between a total value of current values passed through each of the plurality of inductance elements in the predetermined time step and a total value of current values in the previous time step. Is greater than a predetermined total current change amount,
The magnetic field analysis unit, when the current change amount determination unit determines that the total current change amount is larger than a predetermined current change amount, each current of the plurality of inductance elements in the predetermined time step. Based on the value, each inductance of the plurality of inductance elements in the predetermined time step is calculated,
The output parameter storage unit stores inductances of the plurality of inductance elements in the predetermined time step calculated by the magnetic field analysis unit,
The parameter output unit outputs the inductances of the plurality of inductance elements in the predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step,
When the current change amount determination unit determines that the total current change amount is equal to or less than a predetermined total current change amount, the magnetic field analysis control unit causes the magnetic field analysis unit to perform the predetermined time step. Without calculating the inductance of each of the plurality of inductance elements, the inductance of each of the plurality of inductance elements in the previous time step stored in the output parameter storage unit is calculated from the parameter output unit to the current calculation unit. The simulator according to claim 1, wherein the simulator is output. A position calculating unit that is magnetically coupled to the inductance element and calculates the position of the movable part driven by the magnetic field generated by the inductance element for each time step;
A position change amount indicating a difference between the position of the movable part at the predetermined time step calculated by the position calculation unit and the position of the movable part at the previous time step is larger than a predetermined position change amount. A position change amount determination unit for determining whether or not
The magnetic field analysis unit determines that the position change amount determination unit determines that the position change amount is greater than a predetermined position change amount, or the current change amount determination unit determines a current change amount in advance. When it is determined that the current change amount is larger than the current change amount, based on the current value calculated by the current calculation unit in the predetermined time step and the position calculated by the position calculation unit in the predetermined time step. Calculate the inductance of the inductance element at the time step,
The output parameter storage unit stores the inductance in the predetermined time step calculated by the magnetic field analysis unit,
The parameter output unit outputs the inductance stored in the output parameter storage unit to the current calculation unit for each time step,
The magnetic field analysis control unit determines that the position change amount determination unit determines that the position change amount is equal to or less than a predetermined position change amount, and the current change amount determination unit determines a current change amount in advance. If it is determined that the current change amount is less than or equal to the current change amount, the magnetic field analysis unit does not calculate the inductance of the inductance element at the predetermined time step, and the output parameter storage unit stores the previous time. The simulator according to claim 1, wherein the inductance in the step is output from the parameter output unit to the current calculation unit. A position calculation unit that calculates the position of the movable unit magnetically coupled to the inductance element and driven by the magnetic field generated by the inductance element for each time step;
The magnetic field analysis unit, when the time change amount determination unit determines that the time change amount is larger than a predetermined time change amount, the current value in the predetermined time step calculated by the current calculation unit and Based on the position in the predetermined time step calculated by the position calculation unit, to calculate the induced electromotive force by the movable part in the inductance element in the predetermined time step,
The output parameter storage unit stores an induced electromotive force generated by the movable unit in the predetermined time step calculated by the magnetic field analysis unit,
When the time variation determination unit determines that the time variation is equal to or less than a predetermined time variation, the magnetic field analysis control unit causes the magnetic field analysis unit to perform the predetermined time step of the inductance element. Without calculating the induced electromotive force due to the movable part in
The simulator according to claim 2, wherein the parameter output unit outputs the induced electromotive force generated by the movable unit stored in the output parameter storage unit to the current calculation unit for each time step. The magnetic field analysis unit is configured to calculate the current value and the position calculation unit in the predetermined time step calculated by the current calculation unit when the time change amount determination unit is larger than a predetermined time change amount. Calculating the electromagnetic force of the movable part at the predetermined time step based on the calculated position at the predetermined time step;
The output parameter storage unit stores the electromagnetic force of the movable unit in the predetermined time step calculated by the magnetic field analysis unit,
When the time variation determination unit determines that the time variation is equal to or less than a predetermined time variation, the magnetic field analysis control unit causes the magnetic field analysis unit to move the movable portion at the predetermined time step. Without calculating the electromagnetic force of the
The simulator according to claim 6, wherein the parameter output unit outputs the electromagnetic force of the movable unit stored in the output parameter storage unit to the position calculation unit for each time step. The magnetic field analysis unit
Based on the current value of the inductance element at the predetermined time step calculated by the current calculation unit, the flux linkage calculation unit that calculates the flux linkage of the inductance element by the current indicated by the current value;
The interlinkage magnetic flux of the inductance element is calculated by the current indicated by the current value changed by a minute current value from the current value calculated by the current calculation unit at the predetermined time step, and the calculated interlinkage magnetic flux and the interlinkage magnetic flux calculation unit The simulator according to claim 1, further comprising: an inductance calculating unit that calculates an inductance of the inductance element based on a value obtained by dividing a difference from the flux linkage calculated by the value by a minute current value. The current calculation unit calculates a current value flowing through each of the plurality of inductance elements magnetically coupled to each other for each time step,
The interlinkage magnetic flux calculation unit calculates the chain of each of the first inductance element and the second inductance element from the current value of each of the first inductance element and the second inductance element in the predetermined time step calculated by the current calculation unit. Calculate the magnetic flux,
The inductance calculator is
The interlinkage magnetic flux of each of the first inductance element and the second inductance element due to the current indicated by the current value changed by a minute current value from the current value of the first inductance element in the predetermined time step calculated by the current calculation unit. A first interlinkage magnetic flux calculating unit for calculating;
Based on the value obtained by dividing the difference between the flux linkage of the first inductance element calculated by the first flux linkage calculator and the flux linkage of the first inductance element calculated by the flux linkage calculator by a minute current value. A first self-inductance calculating unit for calculating a self-inductance of the first inductance element;
Based on a value obtained by dividing the difference between the flux linkage of the second inductance element calculated by the first flux linkage calculator and the flux linkage of the second inductance element calculated by the flux linkage calculator by a minute current value. A first mutual inductance calculation unit for calculating a mutual inductance of the second inductance element between the first inductance element and the first inductance element;
The interlinkage magnetic flux of each of the first inductance element and the second inductance element due to the current indicated by the current value changed by a minute current value from the current value of the second inductance element in the predetermined time step calculated by the current calculation unit. A second flux linkage calculating unit for calculating;
Based on a value obtained by dividing the difference between the flux linkage of the second inductance element calculated by the second flux linkage calculator and the flux linkage of the second inductance element calculated by the flux linkage calculator by a minute current value. A second self-inductance calculating unit for calculating a self-inductance of the second inductance element;
Based on a value obtained by dividing the difference between the flux linkage of the first inductance element calculated by the second flux linkage calculator and the flux linkage of the first inductance element calculated by the flux linkage calculator by a minute current value. The simulator according to claim 8, further comprising: a second mutual inductance calculation unit that calculates a mutual inductance of the first inductance element between the second inductance element and the second inductance element. The inductance calculator is
The self-inductance of the first inductance element in the previous time step calculated by the first self-inductance calculator, the mutual inductance of the second inductance element in the previous time step calculated by the first mutual inductance calculator, and the Inductance storage for storing the self-inductance of the second inductance element in the previous time step calculated by the second self-inductance calculation unit and the mutual inductance of the first inductance element in the previous time step calculated by the second mutual inductance calculation unit. And
The difference between the self-inductance of the first inductance element at the predetermined time step calculated by the first self-inductance calculation unit and the self-inductance of the first inductance element at the previous time step stored in the inductance storage unit is The previous time which is smaller than a predetermined first value and stored in the inductance storage unit and the mutual inductance of the second inductance element in the predetermined time step calculated by the first mutual inductance calculation unit. When the difference from the mutual inductance of the second inductance element in the step is smaller than a predetermined second value, the flux linkage of the first inductance element is in a linear region with respect to the current value of the first inductance element. A linear region determination unit for determining
When the linear region determination unit determines that the linear region is not in the linear region, the second inductance element self-inductance and the mutual inductance of the first inductance element are set in the second self-inductance calculation unit and the second mutual inductance calculation unit, respectively. And when the linear region determining unit determines that the linear region determining unit is in the linear region, the second self-inductance calculating unit and the second mutual inductance calculating unit cause the self-inductance of the second inductance element and the first inductance element to interact with each other. An inductance calculation control unit that outputs the self-inductance of the second inductance element and the mutual inductance of the first inductance element in the previous time step stored in the inductance storage unit without calculating the inductances respectively; Simulator according to claim 9 that. A time change amount determination unit for determining whether or not a time change amount indicating a difference between a time indicating the predetermined time step and a time indicating the previous time step is larger than a predetermined time change amount;
A position calculating unit that is magnetically coupled to the inductance element and calculates the position of the movable part driven by the magnetic field generated by the inductance element for each time step;
A position change amount indicating a difference between the position of the movable part at the predetermined time step calculated by the position calculation unit and the position of the movable part at the previous time step is larger than a predetermined position change amount. A position change amount determination unit for determining whether or not
The magnetic field analysis unit indicates that the position change amount is larger than the predetermined position change amount on the condition that the time change amount determination unit determines that the time change amount is larger than the predetermined time change amount. When the position change amount determination unit determines, or when the current change amount determination unit determines that the current change amount is larger than a predetermined current change amount, the predetermined value calculated by the current calculation unit Based on the current value in the time step and the position in the predetermined time step calculated by the position calculation unit, the inductance of the inductance element in the predetermined time step is calculated,
The output parameter storage unit stores the inductance in the predetermined time step calculated by the magnetic field analysis unit,
The parameter output unit outputs the inductance stored in the output parameter storage unit to the current calculation unit for each time step,
The magnetic field analysis control unit determines that the time change amount determination unit determines that the time change amount is equal to or less than a predetermined time change amount, or the position change amount is equal to or less than a predetermined position change amount. If the position change amount determination unit determines that the current change amount is equal to or less than a predetermined current change amount, the magnetic field analysis unit receives an inductance element. The inductance in the previous time step stored in the output parameter storage unit is output from the parameter output unit to the current calculation unit without calculating the inductance in the predetermined time step. Simulator. A position calculation unit that calculates the position of the movable unit magnetically coupled to the inductance element and driven by the magnetic field generated by the inductance element for each time step;
The magnetic field analysis unit
Based on the current value in the predetermined time step calculated by the current calculation unit and the position in the predetermined time step calculated by the position calculation unit, the induced electromotive force by the movable unit in the inductance element in the predetermined time step And an induced electromotive force calculation unit for calculating
The interlinkage magnetic flux calculation unit is based on a current value flowing through the plurality of inductance elements in the predetermined time step calculated by the current calculation unit and a position of the movable unit in the predetermined time step calculated by the position calculation unit. , The current indicated by the current value and the flux linkage of the inductance element by the movable part,
The inductance calculation unit includes a current value when a minute current value is changed from a current value at the predetermined time step calculated by the current calculation unit and a position of the movable unit at the predetermined time step calculated by the position calculation unit. Based on the current indicated by the current value changed by a minute current value and the flux linkage of the inductance element by the movable part, and the difference between the calculated flux linkage and the flux linkage calculated by the flux linkage calculator Calculate the inductance of the inductance element based on the value divided by the minute current value,
The induced electromotive force calculator is a magnetic flux change value obtained by subtracting the interlinkage magnetic flux in the previous time step from the interlinkage magnetic flux in the predetermined time step calculated by the interlinkage magnetic flux calculator, and the inductance calculator is the predetermined The inductance calculated in the time step, the current change value obtained by subtracting the current value in the previous time step from the current value in the predetermined time step calculated by the current calculation unit, and the time indicated by the predetermined time step The simulator according to claim 8, wherein an induced electromotive force by the movable part is calculated based on a difference from a time indicated by a previous time step. A simulation method for simulating a circuit including an inductance element,
A current calculation stage for calculating a current value flowing through the inductance element for each time step;
Whether the current change amount indicating the difference between the current value at the predetermined time step calculated in the current calculation step and the current value at the time step before the predetermined time step is larger than a predetermined current change amount A current change amount determination stage for determining whether or not,
In the current change amount determination step, when it is determined that the current change amount is larger than a predetermined current change amount, based on the current value in the predetermined time step calculated in the current calculation step, A magnetic field analysis stage for calculating the inductance of the inductance element at a predetermined time step;
An output parameter storage step for storing inductance at the predetermined time step calculated in the magnetic field analysis step;
A parameter output stage for outputting the inductance at the predetermined time step stored in the output parameter storage stage to the current calculation stage for each time step;
In the current change amount determining step, when it is determined that the current change amount is equal to or less than a predetermined current change amount, the magnetic field analysis step does not calculate the inductance of the inductance element at the predetermined time step, and A simulation method comprising: a magnetic field analysis control step of outputting the inductance in the previous time step stored in the output parameter storage step to the current calculation step in the parameter output step. A program for a simulator for simulating a circuit including an inductance element, the simulator
A current calculation unit that calculates a current value flowing through the inductance element for each time step;
Whether or not a current change amount indicating a difference between a current value at a predetermined time step calculated by the current calculation unit and a current value at a time step before the predetermined time step is larger than a predetermined current change amount A current change amount determination unit for determining whether
When the current change amount determination unit determines that the current change amount is larger than a predetermined current change amount, the predetermined amount of time is calculated based on the current value in the predetermined time step calculated by the current calculation unit. A magnetic field analyzer for calculating the inductance of the inductance element in the time step,
An output parameter storage unit for storing inductance in the predetermined time step calculated by the magnetic field analysis unit;
A parameter output unit that outputs the inductance at the predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step;
When the current change amount determination unit determines that the current change amount is equal to or less than a predetermined current change amount, the magnetic field analysis unit does not calculate the inductance of the inductance element at the predetermined time step, and A program that causes the inductance in the previous time step stored in the output parameter storage unit to function as a magnetic field analysis control unit that outputs the inductance from the parameter output unit to the current calculation unit. A simulator for simulating a circuit including an inductance element,
A position calculating unit that is magnetically coupled to the inductance element and calculates the position of the movable part driven by the magnetic field generated by the inductance element for each time step;
A current calculation unit for calculating a current value flowing through the inductance element for each time step;
A position change amount indicating a difference between the position of the movable portion at a predetermined time step calculated by the position calculation portion and the position of the movable portion at a time step before the time step is a predetermined position change amount. A position change amount determination unit for determining whether or not it is greater than,
When the position change amount determination unit determines that the position change amount is larger than a predetermined position change amount, the current value calculated by the current calculation unit and the position calculation unit calculate A magnetic field analyzer that calculates an inductance of the inductance element at the predetermined time step based on the position at the predetermined time step;
An output parameter storage unit for storing inductance in the predetermined time step calculated by the magnetic field analysis unit;
A parameter output unit that outputs the inductance at the predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step; and
When the position change amount determination unit determines that the position change amount is equal to or less than a predetermined position change amount, the magnetic field analysis unit does not calculate the inductance of the inductance element at the predetermined time step, and A simulator comprising: a magnetic field analysis control unit that outputs the inductance in the previous time step stored in the output parameter storage unit from the parameter output unit to the current calculation unit. A simulation method for simulating a circuit including an inductance element,
A position calculating stage that calculates the position of the movable part magnetically coupled to the inductance element and driven by the magnetic field generated by the inductance element for each time step;
A current calculation stage for calculating a current value flowing through the inductance element for each time step;
A position change amount indicating a difference between the position of the movable part in the predetermined time step calculated in the position calculation stage and the position of the movable part in the time step before the time step is a predetermined position change. A position change amount determination stage for determining whether or not the amount is greater than the amount;
In the position change amount determination step, when it is determined that the position change amount is larger than a predetermined position change amount, the current value in the predetermined time step calculated in the current calculation step and the position calculation step A magnetic field analysis step of calculating an inductance of the inductance element in the predetermined time step based on the position in the predetermined time step calculated in
An output parameter storage step for storing inductance at the predetermined time step calculated in the magnetic field analysis step;
A parameter output stage for outputting the inductance at the predetermined time step stored in the output parameter storage stage to the current calculation stage for each time step;
In the position change amount determination step, when it is determined that the position change amount is equal to or less than a predetermined position change amount, the inductance at the predetermined time step of the inductance element is not calculated in the magnetic field analysis step, A simulation method comprising: a magnetic field analysis control step for outputting the inductance in the previous time step stored in the output parameter storage step in the current calculation step from the parameter output step. A program for a simulator for simulating a circuit including an inductance element, the simulator
A position calculation unit that calculates the position of a movable part that is magnetically coupled to the inductance element and is driven by a magnetic field generated by the inductance element, for each time step;
A current calculation unit that calculates a current value flowing through the inductance element for each time step;
A position change amount indicating a difference between the position of the movable portion at a predetermined time step calculated by the position calculation portion and the position of the movable portion at a time step before the time step is a predetermined position change amount. A position change amount determination unit for determining whether or not the value is greater than
When the position change amount determination unit determines that the position change amount is larger than a predetermined position change amount, the current value calculated by the current calculation unit and the position calculation unit calculate A magnetic field analyzer that calculates an inductance of the inductance element at the predetermined time step based on the position at the predetermined time step;
An output parameter storage unit for storing inductance in the predetermined time step calculated by the magnetic field analysis unit;
A parameter output unit that outputs the inductance at the predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step;
When the position change amount determination unit determines that the position change amount is equal to or less than a predetermined position change amount, the magnetic field analysis unit does not calculate the inductance of the inductance element at the predetermined time step, and A program that causes the inductance in the previous time step stored in the output parameter storage unit to function as a magnetic field analysis control unit that outputs the inductance from the parameter output unit to the current calculation unit. A simulator for simulating a circuit including an inductance element,
A current change amount indicating a difference between a current value in a predetermined time step calculated by a current calculation unit that calculates a current value flowing through the inductance element for each time step and a current value in a time step before the predetermined time step is A current change amount determination unit for determining whether or not the current change amount is greater than a predetermined amount;
When the current change amount determination unit determines that the current change amount is larger than a predetermined current change amount, the predetermined amount of time is calculated based on the current value in the predetermined time step calculated by the current calculation unit. A magnetic field analyzer for calculating the inductance of the inductance element in the time step;
An output parameter storage unit for storing inductance in the predetermined time step calculated by the magnetic field analysis unit;
A parameter output unit that outputs the inductance at the predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step; and
When the current change amount determination unit determines that the current change amount is equal to or less than a predetermined current change amount, the magnetic field analysis unit does not calculate the inductance of the inductance element at the predetermined time step, and A simulator comprising: a magnetic field analysis control unit that outputs the inductance in the previous time step stored in the output parameter storage unit from the parameter output unit to the current calculation unit. A simulation method for simulating a circuit including an inductance element,
A current change amount indicating a difference between a current value in a predetermined time step calculated by a current calculation unit that calculates a current value flowing through the inductance element for each time step and a current value in a time step before the predetermined time step is A current change amount determining step for determining whether or not the current change amount is larger than a predetermined current change amount;
In the current change amount determination step, when it is determined that the current change amount is larger than a predetermined current change amount, based on the current value in the predetermined time step calculated by the current calculation unit, A magnetic field analysis stage for calculating the inductance of the inductance element at a predetermined time step;
An output parameter storage step for storing inductance at the predetermined time step calculated in the magnetic field analysis step;
A parameter output step of outputting the inductance at the predetermined time step stored in the output parameter storage step to the current calculation unit for each time step;
In the current change amount determining step, when it is determined that the current change amount is equal to or less than a predetermined current change amount, the inductance of the inductance element in the predetermined time step is not calculated in the magnetic field analysis step, A simulation method comprising: a magnetic field analysis control step for causing the current calculation unit to output the inductance in the previous time step stored in the output parameter storage step in the parameter output step. A program for a simulator for simulating a circuit including an inductance element, the simulator
A current change amount indicating a difference between a current value in a predetermined time step calculated by a current calculation unit that calculates a current value flowing through the inductance element for each time step and a current value in a time step before the predetermined time step is , A current change amount determination unit that determines whether or not the current change amount is greater than a predetermined amount,
When the current change amount determination unit determines that the current change amount is larger than a predetermined current change amount, the predetermined amount of time is calculated based on the current value in the predetermined time step calculated by the current calculation unit. A magnetic field analyzer for calculating the inductance of the inductance element in the time step,
An output parameter storage unit for storing inductance in the predetermined time step calculated by the magnetic field analysis unit;
A parameter output unit that outputs the inductance at the predetermined time step stored in the output parameter storage unit to the current calculation unit for each time step;
When the current change amount determination unit determines that the current change amount is equal to or less than a predetermined current change amount, the magnetic field analysis unit does not calculate the inductance of the inductance element at the predetermined time step, and A program that causes the inductance in the previous time step stored in the output parameter storage unit to function as a magnetic field analysis control unit that outputs the inductance from the parameter output unit to the current calculation unit.

Images