Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7198334B2 - Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program, and recording medium recording the program - Google Patents
[go: Go Back, main page]

JP7198334B2 - Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program, and recording medium recording the program - Google Patents

Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program, and recording medium recording the program Download PDF

Info

Publication number
JP7198334B2
JP7198334B2 JP2021198366A JP2021198366A JP7198334B2 JP 7198334 B2 JP7198334 B2 JP 7198334B2 JP 2021198366 A JP2021198366 A JP 2021198366A JP 2021198366 A JP2021198366 A JP 2021198366A JP 7198334 B2 JP7198334 B2 JP 7198334B2
Authority
JP
Japan
Prior art keywords
power supply
supply voltage
voltage waveform
coil current
waveform
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2021198366A
Other languages
Japanese (ja)
Other versions
JP2022033920A (en
Inventor
健治 宮田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from JP2017252916A external-priority patent/JP6991061B2/en
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2021198366A priority Critical patent/JP7198334B2/en
Publication of JP2022033920A publication Critical patent/JP2022033920A/en
Application granted granted Critical
Publication of JP7198334B2 publication Critical patent/JP7198334B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Electric Motors In General (AREA)

Description

本発明は、電源電圧波形算出方法、回路連成磁界解析方法、ならびにこれらの方法を実行するためのプログラム及びプログラムを記録した記録媒体に関する。 The present invention relates to a power supply voltage waveform calculation method, a circuit coupled magnetic field analysis method, a program for executing these methods, and a recording medium recording the program.

従来、モータ等の電磁機器を駆動する磁場発生源として、磁石及びコイル電流がある。電磁機器に所定の動作をさせるためには、コイル電流波形を設計通りに流すことが望まれる。通常、電磁機器とそれを駆動する電源回路とを組み合わせた回路連成磁界解析においては、電源電圧波形を入力して、磁場を記述する磁気ベクトルポテンシャルに関する量ならびにコイル電流を未知変数とした解析が実行される。 Conventionally, there are magnets and coil currents as magnetic field generating sources for driving electromagnetic devices such as motors. In order to cause an electromagnetic device to perform a predetermined operation, it is desirable to apply a coil current waveform as designed. Normally, in a circuit coupled magnetic field analysis that combines an electromagnetic device and a power supply circuit that drives it, the power supply voltage waveform is input, and the amount related to the magnetic vector potential that describes the magnetic field and the coil current are the unknown variables. executed.

この場合、所望のコイル電流波形にするために、電源電圧波形の振幅や位相を調整することになる。しかし、これでは、数回の試行錯誤が必要になる。また、1回の計算に関しても、コイル電流波形が未知のため、定常のコイル電流波形を得るのに、多くの計算時間を要する。通常、定常解に到達するのに、100~数100ステップ程度の過渡解析が必要になるため、三次元解析では数日間の計算時間が必要になることもある。 In this case, the amplitude and phase of the power supply voltage waveform are adjusted to obtain a desired coil current waveform. However, this requires several rounds of trial and error. Moreover, since the coil current waveform is unknown even for one calculation, it takes a long calculation time to obtain a steady coil current waveform. Usually, 100 to several 100 steps of transient analysis are required to arrive at a steady-state solution, so three-dimensional analysis may require several days of calculation time.

これを補助するために、所望のコイル電流波形を入力した磁界解析からコイル誘起電圧波形を平均磁気ベクトルポテンシャルの後退差分により求め、巻線抵抗とコイル電流により、コイル線間電圧波形を算出する方法がある(非特許文献1参照)。 In order to assist this, the coil induced voltage waveform is obtained from the backward difference of the average magnetic vector potential from the magnetic field analysis in which the desired coil current waveform is input, and the coil line voltage waveform is calculated from the winding resistance and the coil current. There is (see Non-Patent Document 1).

高橋則雄著 「磁気工学の有限要素法」 朝倉書店 (2013年) pp. 251-253Norio Takahashi, Finite Element Method for Magnetic Engineering, Asakura Shoten (2013) pp. 251-253

しかしながら、非特許文献1に記載された従来法では、所望のコイル電流波形を入力としたコイル線間電圧波形の算出にとどまっており、電源電圧波形の算出には至っていない。
また、定常場への収束を早める方法として、TP-EEC(Time Periodic-Explicit Error Correction)法または簡易TP-EEC法といった手法が知られているが、この方法を用いても、目的の定常場を求めるのに、ある程度の時間ステップによる過渡解析を必要とするため、定常場に達するまである程度時間がかかるという問題があった。
However, in the conventional method described in Non-Patent Document 1, the calculation of the coil line-to-line voltage waveform is limited to the input of the desired coil current waveform, and the power supply voltage waveform is not calculated.
In addition, as a method of hastening the convergence to a stationary field, the TP-EEC (Time Periodic-Explicit Error Correction) method or the simplified TP-EEC method are known. Since a transient analysis with a certain amount of time steps is required to obtain

本発明の目的は、まず第一番目の問題(所望のコイル電流波形になるように電源電圧波形入力による回路連成磁界解析を数回試行錯誤的に実施する必要があるという問題)を解決するための電源電圧波形算出方法を提供することにある。また二番目の問題(電源電圧波形を最終決定した後のコイル電流波形の定常解を得るのにある程度時間を要する)を解決するために、回路連成磁界解析方法、ならびに上記両方法に関するプログラム及びプログラムを記録した記録媒体を提供することにある。 The object of the present invention is to solve the first problem (the problem that it is necessary to perform circuit coupled magnetic field analysis by inputting a power supply voltage waveform several times by trial and error so as to obtain a desired coil current waveform). The object of the present invention is to provide a method for calculating a power supply voltage waveform. In order to solve the second problem (it takes a certain amount of time to obtain a stationary solution of the coil current waveform after the final determination of the power supply voltage waveform), a circuit coupled magnetic field analysis method, a program and a program related to both the above methods, To provide a recording medium recording a program.

上記課題を解決するために、例えば特許請求の範囲に記載の構成を採用する。
本発明は上記課題を解決する手段を複数含んでいるが、その一例を挙げるならば、本発明の電源電圧波形算出方法は、以下のステップ(a)~(d)から構成される。
(a)微分方程式を数値的に解くための離散化データと、解析プロセスを制御するためのコイル電流波形データを含むコントロールデータからなる入力データを計算機に読み込ませるステップ、
(b)入力データを用いて、計算機によりコイル電流波形データを強制電流とした過渡磁界解析を行うステップ、
(c)電源回路における回路方程式を用いて、過渡磁界解析を行うステップで得られたコイル誘起電圧波形と、入力データとして用いたコイル電流波形データから電源電圧波形を計算機により算出するステップ。
In order to solve the above problems, for example, the configurations described in the claims are adopted.
The present invention includes a plurality of means for solving the above problems. To give one example, the power supply voltage waveform calculation method of the present invention comprises the following steps (a) to (d).
(a) causing a computer to read input data consisting of discretized data for numerically solving differential equations and control data including coil current waveform data for controlling the analysis process;
(b) using the input data to perform a transient magnetic field analysis with a computer using the coil current waveform data as a forced current;
(c) A step of calculating a power supply voltage waveform from the coil induced voltage waveform obtained in the transient magnetic field analysis step and the coil current waveform data used as input data using a circuit equation in the power supply circuit.

また、本発明の回路連成磁界解析方法は、上記ステップ(a)~(c)に加えて,以下の(d)~(f)を含む。
(d)電源電圧波形から基本正弦波電圧成分を計算機により算出するステップ、
(e)電源電圧基本正弦波電圧成分を電源電圧波形として計算機により入力するステップ、
(f)コントロールデータとして入力されたコイル電流波形の初期値を未知数であるコイル電流の初期値として、計算機により回路連成磁界解析を行うステップ。
In addition to the above steps (a) to (c), the circuit coupling magnetic field analysis method of the present invention includes the following (d) to (f).
(d) computing a fundamental sinusoidal voltage component from the power supply voltage waveform;
(e) inputting the power supply voltage fundamental sinusoidal voltage component as a power supply voltage waveform into a calculator;
(f) A step of performing a circuit coupling magnetic field analysis with a computer using the initial value of the coil current waveform input as control data as the initial value of the coil current, which is an unknown quantity.

本発明の電源電圧波形算出方法によれば、所望のコイル電流波形を実現するための電源電圧波形を直接算出することができる。また、前述のステップ(a)~(d)はモータ設計プロセスで必ず実施する工程なので、電源電圧波形が決定した段階では、実質、前述のステップ(e)~(g)のみを実行すれば、最終的なコイル電流波形を求めることができる。したがって、本発明の回路連成磁界解析方法によれば、回路連成磁界解析による解析時間を短縮できるという効果がある。その結果、モータ等の電気機器設計開発のための設計開発を短期間化・低コスト化することが可能となる。
上記した以外の課題、構成及び効果は、以下の実施形態例の説明により明らかにされる。
According to the power supply voltage waveform calculation method of the present invention, it is possible to directly calculate the power supply voltage waveform for realizing a desired coil current waveform. In addition, since steps (a) to (d) described above are always performed in the motor design process, when the power supply voltage waveform is determined, substantially only steps (e) to (g) described above are performed. A final coil current waveform can be obtained. Therefore, according to the circuit coupled magnetic field analysis method of the present invention, there is an effect that the analysis time for the circuit coupled magnetic field analysis can be shortened. As a result, it is possible to reduce the time and cost of design development for designing and developing electrical equipment such as motors.
Problems, configurations, and effects other than those described above will be clarified by the following description of exemplary embodiments.

本発明の第1の実施の形態例を説明するためのフローチャートである。4 is a flow chart for explaining the first embodiment of the present invention; 本発明の第2の実施の形態例を説明するためのフローチャートである。9 is a flow chart for explaining a second embodiment of the present invention; 本発明の第3の実施の形態例を説明するためのフローチャートである。9 is a flow chart for explaining a third embodiment of the present invention; 本発明を実装するための装置の概略構成図である。1 is a schematic block diagram of an apparatus for implementing the present invention; FIG. 本発明の第4の実施の形態例を説明するためのフローチャートである。FIG. 14 is a flow chart for explaining a fourth embodiment of the present invention; FIG. 本発明の第5の実施の形態例を説明するためのフローチャートである。FIG. 11 is a flow chart for explaining a fifth embodiment of the present invention; FIG. 本発明の第1~第5の実施の形態例の解析対象を示す同期電動機のコイルと電源の関係を示す図である。FIG. 4 is a diagram showing the relationship between the coils of the synchronous motor and the power supply, which are the objects of analysis in the first to fifth embodiments of the present invention; 従来法において、電源波形からコイル電流波形を安定化させる手順を示す波形図である。FIG. 4 is a waveform diagram showing a procedure for stabilizing a coil current waveform from a power supply waveform in a conventional method; 本発明の方法を用いた回路連成磁界解析から少ない補正回数で安定化したコイル電流波形を得る手順を示す波形図である。FIG. 4 is a waveform diagram showing a procedure for obtaining a coil current waveform stabilized with a small number of corrections from circuit coupled magnetic field analysis using the method of the present invention; 第5の実施形態例における補正を1回行った回路連成磁界解析から得られたコイル電流波形と補正を行わないコイル電流波形の違いを示す波形図である。FIG. 11 is a waveform diagram showing the difference between a coil current waveform obtained from circuit coupled magnetic field analysis with one correction and a coil current waveform without correction in the fifth embodiment.

最初に、図1~3、図5~6に示す本発明の実施の形態例1~5について説明する前に、本発明の課題を解決するための数学的な根拠を従来法と比較して、図7~図9を参照して説明する。 First, before describing Embodiments 1 to 5 of the present invention shown in FIGS. 1 to 3 and FIGS. , with reference to FIGS.

<従来の回路連成磁界解析手法の説明>
図7は、同機モータを含む磁界解析領域100とモータコイル101ならびにモータコイル101に接続された電源回路200を示す。
<Explanation of conventional circuit coupling magnetic field analysis method>
FIG. 7 shows a magnetic field analysis area 100 including the machine motor, a motor coil 101, and a power supply circuit 200 connected to the motor coil 101. FIG.

図7に示すように、同期モータを含む磁界解析領域100のモータコイル101(U、V、W)は、それぞれ抵抗RとインダクタンスLを介して、三相交流電源200(UV、VW、WU)に接続されている。
モータコイル101における設計上適正なコイル電流波形を実現するための電源電圧波形を高速に求めることが第一目的である。
As shown in FIG. 7, the motor coils 101 (U, V, W) in the magnetic field analysis area 100 including the synchronous motor are connected to the three-phase AC power supply 200 (UV, VW, WU) through the resistance R and the inductance L, respectively. It is connected to the.
The first purpose is to quickly obtain a power supply voltage waveform for realizing a coil current waveform that is appropriate for design in the motor coil 101 .

モータコイル101に流すコイル電流波形は、モータの運転特性を最適化する観点で、その波形が定まるのであるが、その際に、どのような電源電圧波形を供給すればよいかは事前にはわからない。そこで、従来法では、ユーザは、三相交流の電源電圧波形を仮設定して、定常のコイル電流波形を高速に求めるために三相交流TP-EEC法または三相交流簡易TP-EEC法という手法を採用していた。 The waveform of the coil current flowing through the motor coil 101 is determined from the viewpoint of optimizing the operating characteristics of the motor. . Therefore, in the conventional method, the user temporarily sets the power supply voltage waveform of the three-phase AC and obtains the steady coil current waveform at a high speed. used the method.

図8は、この三相交流TP-EEC法を用いてコイル電流波形を求める従来の方法を説明するための図である。図8Aに示すように、この従来の手法は、三相交流電源の電圧波形がベースになっている。図8Bは、図8Aに示した三相交流電源の電圧波形を電源200に加え、かつ三相交流簡易TP-EEC法による補正をしないときのコイル電流波形を示している。
これに対して、図8Cは、三相交流電源の電圧波形を元にした回路連成の過渡磁界解析において、三相交流簡易TP-EEC法により、電気角60°相当の時間間隔で3回補正した場合のコイル電流波形を示している。
FIG. 8 is a diagram for explaining a conventional method for obtaining a coil current waveform using this three-phase AC TP-EEC method. As shown in FIG. 8A, this conventional technique is based on the voltage waveform of a three-phase AC power supply. FIG. 8B shows the coil current waveform when the voltage waveform of the three-phase AC power supply shown in FIG. 8A is applied to the power supply 200 and no correction is made by the three-phase AC simple TP-EEC method.
On the other hand, in FIG. 8C, in the transient magnetic field analysis of the circuit coupling based on the voltage waveform of the three-phase AC power supply, three times at a time interval equivalent to an electrical angle of 60° by the three-phase AC simple TP-EEC method It shows the coil current waveform after correction.

図8Bでは、三相交流簡易TP-EEC法による補正をしていないため、モータコイル101のコイル電流波形は、25ms後でも定常状態に達していない。しかし、図8Cに示すように、三相交流簡易TP-EEC法により3回の補正を行った場合には、5ms後にほぼ定常状態のコイル電流波形に達していることが分かる。図8Cのコイル電流波形が不連続に変化している箇所が、5msまでに3箇所認められるが、この不連続な変化点の3箇所が三相交流簡易TP-EEC法による3回の補正に対応している。 In FIG. 8B, the coil current waveform of the motor coil 101 has not reached a steady state even after 25 ms because correction by the three-phase AC simple TP-EEC method is not performed. However, as shown in FIG. 8C, when correction is performed three times by the three-phase AC simplified TP-EEC method, the coil current waveform reaches a substantially steady state after 5 ms. There are three places where the coil current waveform in FIG. 8C changes discontinuously within 5 ms. Yes.

<本発明の回路連成磁界解析手法の説明>
次に、図9を参照して、本発明の回路連成磁界解析手法を用いて、定常的なコイル電流波形を高速に求める手法について説明する。図9Aに示すように、本発明では、設計仕様で決めた三相のコイル電流基本正弦波波形から出発する点が従来の方法とは異なっている。後述するように電源電圧波形を基準にとると、コイル電流波形はこの基本正弦波波形からやや歪んだ波形になる。本発明の手法は、元のコイル電流波形から、図9Bに示すような電源電圧波形を直接求める。図9Bから分かるように、このとき、電源電圧波形はきれいな正弦波になっておらず、モータを構成する磁性体の非線形磁気特性により高調波成分を含んでいる。
<Description of circuit coupling magnetic field analysis method of the present invention>
Next, with reference to FIG. 9, a technique for obtaining a stationary coil current waveform at high speed using the circuit coupled magnetic field analysis technique of the present invention will be described. As shown in FIG. 9A, the present invention differs from the conventional method in that it starts from a three-phase coil current basic sinusoidal waveform determined by design specifications. As will be described later, when the power supply voltage waveform is used as a reference, the coil current waveform becomes a slightly distorted waveform from this basic sinusoidal waveform. The technique of the present invention directly determines the source voltage waveform as shown in FIG. 9B from the original coil current waveform. As can be seen from FIG. 9B, at this time, the power supply voltage waveform is not a clean sine wave, and contains harmonic components due to the nonlinear magnetic properties of the magnetic material that constitutes the motor.

次に図9Bの電源電圧波形から高調波成分を取り除き、基本正弦波電圧成分を抽出する。その結果を図9Cに示す。
最後に、図9Cに示す電源電圧波形の基本正弦波電圧成分を使って回路連成磁界解析を実施する。図9Bと図9Cに示す電源電圧波形は互いに近い波形をしているので、この回路連成磁界解析で求まるコイル電流波形は、図9Aに示した元のコイル電流波形に近い波形になる。このため、この回路連成磁界解析においては、モータコイル101のコイル電流の初期値として、図9Aに示した元のコイル電流波形の初期値を用いることで、過渡解析の初期の段階で定常波形に近い波形を形成できる。より高精度な定常場を得たい場合は、さらにここで三相交流ETF(Error-correction Time-interval Flexible)法による補正を1回実施することで図9Dに示されるように初期の数ステップ目から(図9Dでは2ステップ目から)定常のコイル電流波形が得られる。
Next, the harmonic components are removed from the power supply voltage waveform of FIG. 9B to extract the fundamental sinusoidal voltage component. The results are shown in FIG. 9C.
Finally, circuit coupled magnetic field analysis is performed using the fundamental sinusoidal voltage component of the power supply voltage waveform shown in FIG. 9C. Since the power supply voltage waveforms shown in FIGS. 9B and 9C are close to each other, the coil current waveform obtained by this circuit coupled magnetic field analysis is close to the original coil current waveform shown in FIG. 9A. Therefore, in this circuit coupled magnetic field analysis, by using the initial value of the original coil current waveform shown in FIG. 9A as the initial value of the coil current of the motor coil 101, the steady waveform can form a waveform close to If you want to obtain a steady -state field with higher accuracy , you can further correct the initial field as shown in FIG . A steady coil current waveform is obtained from several steps (from the second step in FIG. 9D).

図8に示す従来の方法と比べて、図9に示す本発明による回路連成磁界解析では、瞬時にコイル電流が定常状態になっていることが分かる。つまり、図8の従来法と比較して、本発明の手法を用いると、定常状態になるまでの過渡解析の時間ステップ数が大幅に削減されている。 Compared to the conventional method shown in FIG. 8, the circuit coupled magnetic field analysis according to the present invention shown in FIG. 9 shows that the coil current is in a steady state instantaneously. That is, compared with the conventional method of FIG. 8, the method of the present invention significantly reduces the number of time steps for transient analysis until the steady state is reached.

<回路方程式に基づいた回路連成磁界解析の説明>
以上、本発明の回路連成磁界解析法の原理の概略を説明したが、本発明では、最初にコイル電流波形の初期値を入力して通常の磁界解析を実施する点に特徴がある。回路系は回路理論により複数の基本ループで構成される。例えば、図7の例では、以下の3つの基本ループが形成されている。
(a)モータコイルU→抵抗R→インダクタンスL→電源WU→インダクタンスL→抵抗R→モータコイルWの基本ループ。
(b)モータコイルU→抵抗R→インダクタンスL→電源UV→インダクタンスL→抵抗R→モータコイルVの基本ループ。
(c)モータコイルV→抵抗R→インダクタンスL→電源VW→インダクタンスL→抵抗R→モータコイルWの基本ループ。
<Explanation of circuit coupling magnetic field analysis based on circuit equation>
The outline of the principle of the circuit coupling magnetic field analysis method of the present invention has been described above, but the present invention is characterized in that the initial value of the coil current waveform is first inputted to perform a normal magnetic field analysis. A circuit system consists of a plurality of basic loops according to circuit theory. For example, in the example of FIG. 7, the following three basic loops are formed.
(a) Basic loop of motor coil U→resistor R→inductance L→power supply WU→inductance L→resistor R→motor coil W;
(b) A basic loop of motor coil U→resistor R→inductance L→power source UV→inductance L→resistor R→motor coil V;
(c) Basic loop of motor coil V→resistor R→inductance L→power source VW→inductance L→resistor R→motor coil W;

このため、回路を流れる電流は複数の基本ループ電流Ik ioop(k=1~n 但し、ここではn=3)から形成され、電源電圧も複数の基本ループ電圧V loopで表現することができる。基本ループ電圧V loopは基本ループ電流I loopを用いて次式で求められる。 Therefore, the current flowing through the circuit is formed from a plurality of basic loop currents I k ioop (k = 1 to n, where n = 3 here), and the power supply voltage can also be expressed by a plurality of basic loop voltages V k loop . can. The basic loop voltage V k loop is obtained by the following equation using the basic loop current I k loop .

Figure 0007198334000001
Figure 0007198334000001
Figure 0007198334000002
Figure 0007198334000002
Figure 0007198334000003
Figure 0007198334000003

ここで、式(1)の右辺の第1項は基本ループkにおける複数のコイル誘起電圧項、第2項は基本ループkにおける複数の抵抗による電圧降下量である。また、第3項は基本ループkにおける複数のインダクタンスによる電圧降下量である。図7に示すような、三相交流同期モータの回路連成磁界解析の場合、kは1~3である。 Here, the first term on the right side of equation (1) is a plurality of coil induced voltage terms in the basic loop k, and the second term is a voltage drop due to a plurality of resistances in the basic loop k. The third term is the amount of voltage drop due to multiple inductances in the basic loop k. In the case of circuit coupled magnetic field analysis of a three-phase AC synchronous motor as shown in FIG. 7, k is 1-3.

また、式(1)におけるaは有限要素領域の辺jに割り当てられた磁気ベクトルポテンシャルに関する未知変数である。式(2)におけるR、Lは、回路上におけるk番目の抵抗及びインダクタンスであり、Tikはi番目の回路要素を流れる電流Iiとk番目の基本ループ電流Ik loopとを連結する行列で、行列成分は「+1」か「-1」か「0」の値を持っている。 Also, a j in equation (1) is an unknown variable related to the magnetic vector potential assigned to side j of the finite element region. R k and L k in equation (2) are the k-th resistance and inductance on the circuit, and T ik is the current I i flowing through the i-th circuit element and the k-th basic loop current I k loop . matrix, where the matrix elements have values of '+1', '-1', or '0'.

また、式(3)におけるNjは有限要素領域における辺jに関する辺要素ベクトル基底関数で、npはp番目のコイルの電流の向きを表す単位ベクトルである。i番目の回路要素を流れる電流Iiは、式(4)で表すことができる。 Also, N j in Equation (3) is a side element vector basis function for side j in the finite element domain, and n p is a unit vector representing the current direction of the p-th coil. A current I i flowing through the i-th circuit element can be expressed by Equation (4).

Figure 0007198334000004
Figure 0007198334000004

基本ループ電圧Vk loopは基本ループkに電源が1個しかなければ、電源電圧そのものであり、複数の電源が存在していれば、向きを考慮した総和量になる。
図7に示す回路構成であれば、電源200(WU、VW、UV)の中の1個の電源がそれぞれの基本ループの電源に対応しているので、電源電圧は基本ループ電圧から直接算出される。
The basic loop voltage V k loop is the power supply voltage itself if there is only one power supply in the basic loop k, and if there are a plurality of power supplies, it becomes the total amount considering the direction.
With the circuit configuration shown in FIG. 7, one power supply in the power supply 200 (WU, VW, UV) corresponds to the power supply of each basic loop, so the power supply voltage is calculated directly from the basic loop voltage. be.

但し、三相交流電源系のY結線において、中性点の連結線がない場合は、基本ループ上に2個の電源が存在することになる。しかし、この場合には、三相交流電源の総和電圧がゼロとみなせる場合、この条件式を用いて、基本ループ電圧から各電源の電圧を算出することができる。 However, in the Y-connection of the three-phase AC power supply system, if there is no neutral point connection line, two power supplies exist on the basic loop. However, in this case, if the total voltage of the three-phase AC power supply can be regarded as zero, the voltage of each power supply can be calculated from the basic loop voltage using this conditional expression.

なお、式(1)は、回路上の基本ループ電流を用いた式であるが、本発明の連成磁界解析は、この式に限定されるものではなく、回路上の各点の電位を用いた回路方程式としてもよい。本発明のポイントは、回路方程式を解いてコイル電流波形から電源電圧波形を算出することにあり、用いる回路方程式の種類を限定するものではない。 Note that the equation (1) is an equation using the basic loop current on the circuit, but the coupled magnetic field analysis of the present invention is not limited to this equation, and the potential at each point on the circuit is used. It is also possible to use a circuit equation with The point of the present invention is to solve the circuit equation and calculate the power supply voltage waveform from the coil current waveform, and the type of circuit equation to be used is not limited.

さらに、上記計算式で得られる電源電圧波形から基本波形成分の振幅及び初期位相を抽出し、これを回路連成磁界解析の入力電圧波形の主成分として、通常の磁界解析を行う。このとき、前述したコイル電流波形の初期値を磁界解析の初期値に用いる。 Further, the amplitude and initial phase of the fundamental waveform component are extracted from the power supply voltage waveform obtained by the above formula, and normal magnetic field analysis is performed using these as the main components of the input voltage waveform for circuit coupled magnetic field analysis. At this time, the initial value of the coil current waveform described above is used as the initial value of the magnetic field analysis.

図9で説明したように、この場合、得られるコイル電流波形は、元のコイル電流波形に近い波形になる。したがって、定常解に近いコイル電流波形の初期値から回路連成磁界解析を始めるので、定常解への収束が高速になり、定常解に到達できるまでの解析時間を大幅に短縮することが可能になる。 As described with reference to FIG. 9, in this case, the obtained coil current waveform becomes a waveform close to the original coil current waveform. Therefore, since the circuit coupled magnetic field analysis starts from the initial value of the coil current waveform that is close to the steady-state solution, the convergence to the steady-state solution becomes faster, and the analysis time required to reach the steady-state solution can be greatly shortened. Become.

さらに、少なくとも2タイムステップの過渡解析を実施後、1回の三相交流ETF法による補正を実施してもよい。ここでタイムステップとは、過渡解析における時間軸分割における1刻みを意味する。 Furthermore, after performing transient analysis for at least two time steps, one correction by the three-phase AC ETF method may be performed. Here, the time step means one increment in the time axis division in the transient analysis.

また、本発明に係る電源電圧波形算出法に関するプログラム、回路連成磁界解析に関するプログラムは、上記の一連のプロセスをコーディングすることを特徴としている。 Further, the program relating to the power supply voltage waveform calculation method and the program relating to circuit coupled magnetic field analysis according to the present invention are characterized by coding the series of processes described above.

以下、図1~6を参照して本発明の実施の形態例について詳細に説明する。
<第1の実施の形態例の説明>
図1は、本発明の第1の実施の形態例の解析のプロセスを説明するためのフローチャートである。
図1に示すように、最初に、微分方程式を数値的に解くための離散化データ11(メッシュデータ)及び解析プロセスをコントロールするためのコントロールデータ12(コイル電流波形データを含む)から構成される入力データファイル10を計算機に読み込ませる(ステップS13、S14)。
Embodiments of the present invention will be described in detail below with reference to FIGS. 1 to 6. FIG.
<Description of First Embodiment>
FIG. 1 is a flow chart for explaining the analysis process of the first embodiment of the present invention.
As shown in FIG. 1, first, it consists of discretized data 11 (mesh data) for numerically solving differential equations and control data 12 (including coil current waveform data) for controlling the analysis process. The input data file 10 is read into the computer (steps S13 and S14).

ここで離散化データ11(メッシュデータ)は、構造物(ここでは、図7の磁界解析領域にある同期モータ100をいう)を複数の有限個の要素(メッシュ)に分解したデータである。一般に解析空間の方程式は、連続では解くことができないので、メッシュに分割して解く有限要素法(FEM:Finite Element Method)という手法が用いられる。つまり、有限要素法とは、解析空間を分割して数値解析を行う手法であり、メッシュが細かいほど解析の精度は上がる。コントロールデータ12は、解析条件に関する入力データである。このコントロールデータの中には、コイル電流波形の初期値データが含まれている。 Here, the discretized data 11 (mesh data) is data obtained by decomposing a structure (here, the synchronous motor 100 in the magnetic field analysis area of FIG. 7) into a plurality of finite elements (mesh). In general, equations in an analytical space cannot be solved continuously, so a method called the finite element method (FEM) is used to solve the equations by dividing them into meshes. In other words, the finite element method is a method of performing numerical analysis by dividing the analysis space, and the finer the mesh, the higher the precision of the analysis. Control data 12 is input data relating to analysis conditions. This control data includes initial value data of the coil current waveform.

これら入力データは、解析実行モジュール20に供給される。解析実行モジュール20は、これらの入力データを用いて、上記コイル電流波形を強制電流とした過渡磁界解析を実行する(ステップS21)。 These input data are supplied to the analysis execution module 20 . The analysis executing module 20 uses these input data to execute a transient magnetic field analysis with the coil current waveform as a forced current (step S21).

次に、タイムステップ毎に得られた未知変数ajの値を、式(1)に代入することにより、すべての基本ループ電圧波形を求めるプロセスを実行する(ステップS22)。そして、得られた基本ループ電圧波形から、電源電圧波形を算出する(ステップS23)。その後、得られた電源電圧波形を記憶装置3(図4参照)に記憶し(ステップS31)、この電源電圧波形を表示装置2(図4参照)に表示する(ステップS32)。 Next, by substituting the value of the unknown variable a j obtained at each time step into the equation (1), a process of obtaining all basic loop voltage waveforms is executed (step S22). Then, the power supply voltage waveform is calculated from the obtained basic loop voltage waveform (step S23). After that, the obtained power supply voltage waveform is stored in the storage device 3 (see FIG. 4) (step S31), and this power supply voltage waveform is displayed on the display device 2 (see FIG. 4) (step S32).

上述の第1の実施形態例によれば、コントロールデータとして入力されたコイル電流波形の初期値から直接電源電圧波形を算出して、電源電圧波形を表示装置2に表示することができるので、ユーザは、電源電圧波形を直接把握することができる。 According to the first embodiment described above, the power supply voltage waveform can be directly calculated from the initial value of the coil current waveform input as control data, and the power supply voltage waveform can be displayed on the display device 2. can directly grasp the power supply voltage waveform.

<第2の実施形態例の説明>
次に、図2を参照して、本発明の第2の実施形態例について説明する。図2において、電源電圧波形を算出するプロセスであるステップS23までの処理は第1の実施形態例と同じなので説明は省略する。
<Description of Second Embodiment>
A second embodiment of the present invention will now be described with reference to FIG. In FIG. 2, the processing up to step S23, which is the process of calculating the power supply voltage waveform, is the same as in the first embodiment, so the description is omitted.

第2の実施形態例では、ステップS23で求めた電源電圧波形から基本正弦波電圧成分を算出している(ステップS24)。基本正弦波電圧成分は、ステップS23で求めた電源電圧波形から高調波を除くことで容易に算出することができる。ステップS23で得られた電源電圧波形及びステップS24で求めた基本正弦波電圧成分は、記憶装置3に記憶されるとともに(ステップ31)、表示装置2に表示される(ステップS32)。 In the second embodiment, the basic sinusoidal voltage component is calculated from the power supply voltage waveform obtained in step S23 (step S24). The fundamental sinusoidal voltage component can be easily calculated by removing harmonics from the power supply voltage waveform obtained in step S23. The power supply voltage waveform obtained in step S23 and the fundamental sinusoidal voltage component obtained in step S24 are stored in the storage device 3 (step S31) and displayed on the display device 2 (step S32).

本発明の第2の実施形態例によれば,初期に入力したコイル電流波形から算出した電源電圧波形から、基本正弦波電圧成分を算出し、電源電圧波形とともに基本正弦波電圧成分を表示装置2に表示することができるので、ユーザは電源電圧波形と基本正弦波電圧成分をともに直接把握することが可能となる。 According to the second embodiment of the present invention, the basic sine wave voltage component is calculated from the power supply voltage waveform calculated from the initially input coil current waveform, and the basic sine wave voltage component is displayed on the display device 2 together with the power supply voltage waveform. Therefore, the user can directly grasp both the power supply voltage waveform and the fundamental sine wave voltage component.

<第3の実施形態例の説明>
次に、図3の解析プロセスを参照して、本発明の第3の実施形態例について説明する。この第3の実施形態例でも、電源電圧波形を算出するステップS23までは第1または第2の実施の形態例と同じあるから、その説明は省略する。
<Description of Third Embodiment>
A third exemplary embodiment of the present invention will now be described with reference to the analysis process of FIG. This third embodiment is the same as the first or second embodiment up to step S23 for calculating the power supply voltage waveform, so the description thereof will be omitted.

第3の実施形態例では、ステップS23で求めた電源電圧波形からパルス幅変調(PWM:Pulse Width Modulation)電圧波形を算出する(ステップS25)。そして、ステップS23とステップS25で得られた電源電圧波形及びPWM電圧波形を記憶装置3に記憶し(ステップS31)、この電源電圧波形及びPWM電圧波形を表示装置2に表示する(ステップS32)。 In the third embodiment, a pulse width modulation (PWM) voltage waveform is calculated from the power supply voltage waveform obtained in step S23 (step S25). Then, the power supply voltage waveform and PWM voltage waveform obtained in steps S23 and S25 are stored in the storage device 3 (step S31), and the power supply voltage waveform and PWM voltage waveform are displayed on the display device 2 (step S32).

第3の実施形態例は、電源電圧波形の基本正弦波電圧波形の代わりに、PWM電圧波形を用いて以降の解析を行うものであり、例えば、数kHz~20kHzの搬送波を使ってPWM電圧波形を作ることができる。PWM電圧波形を使うとオンとオフの繰り返しでスイッチング制御を行うことができるので、多少のノイズは発生するものの、周波数の決まった基本正弦波電圧駆動に比べてモータの回転数制御がやりやすくなるというメリットがある。 In the third embodiment, instead of the basic sine wave voltage waveform of the power supply voltage waveform, the PWM voltage waveform is used for the subsequent analysis. can be made. By using a PWM voltage waveform, switching control can be performed by repeatedly turning on and off, so although some noise is generated, it is easier to control the motor rotation speed than with a basic sine wave voltage drive with a fixed frequency. There is an advantage.

本発明の第1~第3の実施の形態例を実現する解析システムの一例を図4に示す。
図4に示すように、本発明の第1~第3の実施の形態例の解析システムは、計算機1、表示装置2、記憶装置3、及び入力装置4から構成される。
FIG. 4 shows an example of an analysis system that implements the first to third embodiments of the present invention.
As shown in FIG. 4, the analysis system of the first to third embodiments of the present invention comprises a computer 1, a display device 2, a storage device 3, and an input device 4. FIG.

図4では、記憶装置3は、それを明示するために計算機1の外に出しているが、計算機1の内部に設置してもよい。計算機1には、上述した第1~第3の実施の形態例の解析プロセスのうち、少なくともいずれか1つのアルゴリズムに基づくプログラムが格納されている。 In FIG. 4, the storage device 3 is shown outside the computer 1 for clarity, but it may be installed inside the computer 1 as well. The computer 1 stores a program based on at least one of the analysis processes of the first to third embodiments described above.

入力装置4は、例えばキーボードやマウスであり、この入力装置4から、前述した磁界解析領域のメッシュデータや、その他の操作に必要なコントロールデータを含む入力データが入力される。入力データが入力されると、計算機1は、格納されているプログラムに従い、入力データの読み取りやメッシュデータに関連した演算処理を実行する。そして、処理結果を表示装置2に表示し、あるいは、処理結果をデータファイルとして記憶装置3に記憶する。なお、得られた処理結果の全てではなく、その一部を表示したり記憶したりしてもよい。 The input device 4 is, for example, a keyboard or a mouse, and input data including mesh data of the magnetic field analysis region and control data required for other operations are input from the input device 4 . When the input data is input, the computer 1 reads the input data and performs arithmetic processing related to the mesh data according to the stored program. Then, the processing result is displayed on the display device 2, or the processing result is stored in the storage device 3 as a data file. It should be noted that not all of the obtained processing results, but some of them may be displayed or stored.

<第4の実施形態例の説明>
次に、図5に示す解析プロセスを参照して、本発明の第4の実施形態例について説明する。図5に示すように、電源電圧波形から基本正弦波電圧成分を算出するステップS24までのプロセスは第2の実施形態例と同じなので、その説明は省略する。
<Description of Fourth Embodiment>
A fourth example embodiment of the present invention will now be described with reference to the analysis process shown in FIG. As shown in FIG. 5, the process up to step S24 for calculating the basic sine wave voltage component from the power supply voltage waveform is the same as in the second embodiment, so the description thereof will be omitted.

ステップS24で求めた電源電圧波形の基本正弦波電圧成分は、電源電圧波形として計算機1(図4)に入力される(ステップS41)。計算機1は、コントロールデータとして入力したコイル電流波形の初期値を、これから求めるコイル電流(未知変数)の初期値として、回路連成磁界解析を実行する(ステップS42)。
そして、ステップS42で得られたコイル電流波形及びトルク波形等のモータ特性データを記憶装置3に記憶する(ステップS31)。このコイル電流波形及びトルク波形等のモータ特性データは表示装置2に表示される(ステップS32)。
The basic sinusoidal voltage component of the power supply voltage waveform obtained in step S24 is input to the calculator 1 (FIG. 4) as a power supply voltage waveform (step S41). Calculator 1 executes the circuit coupling magnetic field analysis using the initial value of the coil current waveform input as control data as the initial value of the coil current (unknown variable) to be obtained (step S42).
Then, the motor characteristic data such as the coil current waveform and the torque waveform obtained in step S42 are stored in the storage device 3 (step S31). The motor characteristic data such as the coil current waveform and torque waveform are displayed on the display device 2 (step S32).

ここで、ステップS24で得られる電源電圧波形の基本正弦波電圧成分は、一般的に電源電圧波形に近い波形をしているので、回路連成磁界解析で得られるコイル電流波形もステップS14で用いられるコントロールデータ12に入力したコイル電流の初期値に近い波形をしている。 Here, since the fundamental sinusoidal voltage component of the power supply voltage waveform obtained in step S24 generally has a waveform close to the power supply voltage waveform, the coil current waveform obtained by circuit coupled magnetic field analysis is also used in step S14. It has a waveform close to the initial value of the coil current input to the control data 12 to be supplied.

このため、本発明の第4の実施形態例によれば,ステップS42で、未知変数であるコイル電流を、近似的なコイル電流の初期値から出発して回路連成磁界解析を実行できるので、従来法に比べて格段と高速に準定常解が得られる。 For this reason, according to the fourth embodiment of the present invention, in step S42, the coil current, which is an unknown variable, can be started from an approximate initial value of the coil current to perform the circuit coupled magnetic field analysis. A quasi-stationary solution can be obtained much faster than the conventional method.

<第5の実施形態例の説明>
次に、図6を参照して本発明の第5の実施形態例を説明する。
第5の実施形態例は、電源電圧波形の基本正弦波電圧成分を電源電圧波形として入力するプロセス41までは、第4の実施形態例と同じなので、その説明は省略する。
<Description of Fifth Embodiment>
Next, a fifth embodiment of the present invention will be described with reference to FIG.
The fifth embodiment is the same as the fourth embodiment up to the process 41 for inputting the basic sinusoidal voltage component of the power supply voltage waveform as the power supply voltage waveform, so the explanation thereof will be omitted.

上述したように、ステップS42では、コントロールデータに含まれるコイル電流波形の初期値を、求める未知変数であるコイル電流の初期値にして回路連成磁界解析を実行した。このステップS42において、タイムステップ1及びタイムステップ2では通常の回路連成磁界解析が実施される。この回路連成磁界解析におけるタイムステップ2が終了した段階で、時間微分項が解析に関与する三相成分(代表的には三相コイル鎖交磁束)に後述する三相交流ETF法による補正を施す(ステップS43)。 As described above, in step S42, the initial value of the coil current waveform included in the control data is set to the initial value of the coil current, which is the unknown variable to be obtained, and the circuit coupling magnetic field analysis is executed. In step S42, normal circuit coupled magnetic field analysis is performed in time step 1 and time step 2. FIG. At the end of time step 2 in this circuit coupled magnetic field analysis, the three-phase component (typically three-phase coil interlinkage magnetic flux) whose time differential term is involved in the analysis is corrected by the three-phase AC ETF method described later. (step S43).

図10は、第5の実施形態例の解析結果を示す図である。図10Aは、図6のステップS43の補正しない場合の例であり、図10Bは、ステップS43で2タイムステップ目の解析終了後に、三相交流ETF法による補正を1回実施した場合の例である。
図10Aの左図に示すように、元のコイル電流波形の初期値を用いた場合、もともと定常解に近い初期値から出発したことになるので、補正がなくともほぼ定常解に近いコイル電流波形が早期に得られる。
FIG. 10 is a diagram showing analysis results of the fifth embodiment. FIG. 10A is an example in which correction is not performed in step S43 of FIG. 6, and FIG. 10B is an example in which correction is performed once by the three-phase AC ETF method after the analysis of the second time step is completed in step S43. be.
As shown in the left diagram of FIG. 10A, when the initial value of the original coil current waveform is used, the starting value is originally close to the stationary solution. is obtained early.

図10の右側の図は、図10の左側の波形図のピーク付近を拡大して示した図である。この図10Aの右側の図を見ると分かるように、ステップS43の補正なしの場合、三相コイル電流の振幅誤差は±1.6%である。一方、図10Bの右側の図に示すように、三相交流ETF法による1回の補正を行った場合には、コイル電流の振幅誤差が±0.04%になり、さらに高精度の定常解を得ることができることが分かる。 The diagram on the right side of FIG. 10 is an enlarged diagram showing the vicinity of the peak of the waveform diagram on the left side of FIG. As can be seen from the diagram on the right side of FIG. 10A, the amplitude error of the three-phase coil current is ±1.6% without the correction in step S43. On the other hand, as shown in the diagram on the right side of FIG. 10B, when one correction is performed by the three-phase AC ETF method, the amplitude error of the coil current becomes ±0.04%, and a steady solution with higher accuracy is obtained. I know it can be done.

<三相交流ETF法の説明>
この三相交流ETF法は、既に特開2014-142701号公報「時間周期非線形場の解析方法及び時間周期非線形場の解析プログラム」に開示されている方法であるが、具体的には次式で示す補正を実施している。
<Explanation of three-phase AC ETF method>
This three-phase AC ETF method is a method already disclosed in Japanese Unexamined Patent Application Publication No. 2014-142701 "Method for analyzing time-period nonlinear field and analysis program for time-period nonlinear field". The correction shown is implemented.

Figure 0007198334000005
Figure 0007198334000005

Figure 0007198334000006
Figure 0007198334000006

Figure 0007198334000007
Figure 0007198334000007

Figure 0007198334000008
Figure 0007198334000008

Figure 0007198334000009
Figure 0007198334000009

ここに、U、V、Wという記号は補正対象である三相成分を意味する量であり、Unew,Vnew,Wnewは補正後のU、V、Wの値である。dU、dV、dWは、それぞれタイムステップ1からタイムステップ2にかけてのU、V、Wの変動量である。また、dθは時間ステップ幅を電気角に換算した量である。なお、αは任意の値でよいが、通常は1/2の値が用いられる。 Here, the symbols U, V, and W are quantities that mean the three-phase components to be corrected, and U new , V new , and W new are the values of U, V, and W after correction. dU, dV, and dW are the amounts of variation of U, V, and W from time step 1 to time step 2, respectively. Also, dθ is the amount obtained by converting the time step width into an electrical angle. Although α may be any value, a value of 1/2 is normally used.

以上、本発明の第1~第5の実施形態例を順に説明してきた。ここで説明した例はあくまで実施の形態例であり、本発明はこれらに限定されるものではない。本発明は、特許請求の範囲に記載した内容を逸脱しない限りにおいて、他の応用例及び変形例を含むことは勿論である。 So far, the first to fifth embodiments of the present invention have been described in order. The examples described here are merely embodiments, and the present invention is not limited to these. The present invention naturally includes other applications and modifications without departing from the scope of the claims.

1…計算機、2…表示装置、3…記憶媒体、4…入力装置、10…入力データファイル、11…離散化データ、12…コントロールデータ、20 …解析プロセス、100…同期モータを含む磁界解析領域、101…モータコイル、200…電源 1 computer, 2 display device, 3 storage medium, 4 input device, 10 input data file, 11 discretized data, 12 control data, 20 analysis process, 100 magnetic field analysis area including synchronous motor , 101... Motor coil, 200... Power supply

Claims (11)

微分方程式を数値的に解くための離散化データと、解析プロセスを制御するためのコイル電流波形データを含むコントロールデータからなる入力データを計算機に読み込ませるステップと、
前記入力データを用いて、前記計算機により、前記コイル電流波形データを強制電流とした過渡磁界解析を行うステップと、
電源回路における回路方程式を用いて、前記過渡磁界解析を行うステップで得られたコイル誘起電圧波形と、入力データとして用いた前記コイル電流波形データから電源電圧波形を前記計算機により算出するステップと、
を含む電源電圧波形算出方法。
causing a computer to read input data consisting of discretized data for numerically solving a differential equation and control data including coil current waveform data for controlling the analysis process;
using the input data to perform a transient magnetic field analysis with the computer using the coil current waveform data as a forced current;
A step of calculating a power supply voltage waveform from the coil induced voltage waveform obtained in the step of performing the transient magnetic field analysis and the coil current waveform data used as input data by the calculator, using a circuit equation in the power supply circuit;
Power supply voltage waveform calculation method including.
前記コイル電流波形データから電源電圧波形を前記計算機により算出するステップは、回路上の各点の電位を用いた回路方程式を用いてコイル電流波形から電源電圧波形を算出するステップを含む
請求項1に記載の電源電圧波形算出方法。
2. The step of calculating the power supply voltage waveform from the coil current waveform data by the calculator includes the step of calculating the power supply voltage waveform from the coil current waveform using a circuit equation using the potential of each point on the circuit. The described power supply voltage waveform calculation method.
さらに、前記計算機により算出された解析結果である前記電源電圧波形を記憶装置に記憶するステップと、
前記計算機により、前記解析結果を表示装置に表示するステップと、
を含む請求項1または2に記載の電源電圧波形算出方法。
Furthermore, a step of storing the power supply voltage waveform, which is the analysis result calculated by the computer, in a storage device;
a step of displaying the analysis result on a display device by the computer;
The power supply voltage waveform calculation method according to claim 1 or 2, comprising:
さらに、前記計算機により、前記電源電圧波形から電源電圧基本正弦波電圧成分またはPWM電圧波形を算出するステップを含む
請求項1または2に記載の電源電圧波形算出方法。
3. The power supply voltage waveform calculation method according to claim 1, further comprising a step of calculating a power supply voltage basic sine wave voltage component or a PWM voltage waveform from the power supply voltage waveform by the calculator.
微分方程式を数値的に解くための離散化データと、解析プロセスを制御するためのコイル電流波形データを含むコントロールデータからなる入力データを計算機に読み込ませるステップと、
前記入力データを用いて、前記計算機により、前記コイル電流波形データを強制電流とした過渡磁界解析を行うステ
ップと、
電源回路における回路方程式を用いて、前記過渡磁界解析を行うステップで得られたコイル誘起電圧波形と、入力データとして用いた前記コイル電流波形データから電源電圧波形を前記計算機により算出するステップと、
前記電源電圧波形から電源電圧基本正弦波電圧成分を前記計算機により算出するステップと、
前記電源電圧基本正弦波電圧成分を電源電圧波形として前記計算機により入力するステップと、
前記コントロールデータとして入力されたコイル電流波形の初期値を未知変数であるコイル電流の初期値として、前記計算機により回路連成磁界解析を行うステップと、
を含む回路連成磁界解析方法。
causing a computer to read input data consisting of discretized data for numerically solving a differential equation and control data including coil current waveform data for controlling the analysis process;
using the input data to perform a transient magnetic field analysis with the computer using the coil current waveform data as a forced current;
A step of calculating a power supply voltage waveform from the coil induced voltage waveform obtained in the step of performing the transient magnetic field analysis and the coil current waveform data used as input data by the calculator, using a circuit equation in the power supply circuit;
a step of calculating a power supply voltage fundamental sinusoidal voltage component from the power supply voltage waveform by the calculator;
inputting the power supply voltage fundamental sinusoidal voltage component as a power supply voltage waveform into the calculator;
a step of performing a circuit coupled magnetic field analysis with the computer using the initial value of the coil current waveform input as the control data as the initial value of the coil current, which is an unknown variable;
A circuit coupled magnetic field analysis method including
前記コイル電流波形データから電源電圧波形を前記計算機により算出するステップは、回路上の各点の電位を用いた回路方程式を用いてコイル電流波形から電源電圧波形を算出するステップを含む
請求項5に記載の回路連成磁界解析方法。
6. The step of calculating the power supply voltage waveform from the coil current waveform data by the calculator includes the step of calculating the power supply voltage waveform from the coil current waveform using a circuit equation using the potential of each point on the circuit. The circuit coupled magnetic field analysis method described.
さらに、前記計算機により、少なくとも2タイムステップの過渡解析を実施した後に、三相交流ETF法による補正を1回行うステップを、
含む請求項5または6に記載の回路連成磁界解析方法。
Furthermore, after performing a transient analysis for at least two time steps with the computer, a step of performing correction once by the three-phase AC ETF method,
7. The circuit coupled magnetic field analysis method according to claim 5 or 6.
微分方程式を数値的に解くための離散化データと、解析プロセスを制御するためのコントロールデータからなる入力データを計算機に読み込ませる手順と、
前記入力データを用いて、コイル電流波形を強制電流とした過渡磁界解析を前記計算機により行う手順と、
電源回路における回路方程式を用いて、前記過渡磁界解析を行う手順で得られたコイル誘起電圧波形と、入力データとして用いたコイル電流波形データから電源電圧波形を前記計算機により算出する手順と、
を前記計算機に実行させるためのプログラム。
a procedure for reading input data into a computer, which consists of discretized data for numerically solving a differential equation and control data for controlling the analysis process;
A procedure for performing a transient magnetic field analysis with the computer using the input data with the coil current waveform as a forced current;
Using circuit equations in a power supply circuit, a coil induced voltage waveform obtained in the procedure for performing the transient magnetic field analysis, and a procedure for calculating the power supply voltage waveform by the computer from the coil current waveform data used as input data;
is executed by the computer.
前記コイル電流波形データから電源電圧波形を前記計算機により算出する手順は、回路上の各点の電位を用いた回路方程式を用いてコイル電流波形から電源電圧波形を算出する手順を含む
請求項8に記載のプログラム。
9. The step of calculating the power supply voltage waveform from the coil current waveform data by the computer includes the step of calculating the power supply voltage waveform from the coil current waveform using a circuit equation using the potential of each point on the circuit. program as described.
さらに、前記電源電圧波形から電源電圧基本正弦波電圧成分を前記計算機により算出する手順と、
前記電源電圧基本正弦波電圧成分を、前記計算機により電源電圧波形として入力する手順と、
前記コントロールデータとして入力されたコイル電流波形の初期値を未知変数であるコイル電流の初期値として、前記計算機により回路連成磁界解析を行う手順と、
を前記計算機に実行させるための請求項8または9に記載のプログラム。
Further, a step of calculating a power supply voltage fundamental sine wave voltage component from the power supply voltage waveform by the calculator;
a step of inputting the power supply voltage fundamental sinusoidal voltage component as a power supply voltage waveform using the calculator;
A step of performing circuit coupled magnetic field analysis by the computer using the initial value of the coil current waveform input as the control data as the initial value of the coil current, which is an unknown variable;
10. The program according to claim 8 or 9, for causing said computer to execute:
請求項8~10のいずれか1項に記載のプログラムを記録した記録媒体。


A recording medium recording the program according to any one of claims 8 to 10.


JP2021198366A 2017-12-28 2021-12-07 Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program, and recording medium recording the program Active JP7198334B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021198366A JP7198334B2 (en) 2017-12-28 2021-12-07 Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program, and recording medium recording the program

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017252916A JP6991061B2 (en) 2017-12-28 2017-12-28 Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program and recording medium on which the program is recorded
JP2021198366A JP7198334B2 (en) 2017-12-28 2021-12-07 Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program, and recording medium recording the program

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2017252916A Division JP6991061B2 (en) 2017-12-28 2017-12-28 Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program and recording medium on which the program is recorded

Publications (2)

Publication Number Publication Date
JP2022033920A JP2022033920A (en) 2022-03-02
JP7198334B2 true JP7198334B2 (en) 2022-12-28

Family

ID=87885122

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021198366A Active JP7198334B2 (en) 2017-12-28 2021-12-07 Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program, and recording medium recording the program

Country Status (1)

Country Link
JP (1) JP7198334B2 (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014142701A (en) 2013-01-22 2014-08-07 Hitachi Ltd Analysis method of time period nonlinear field and analysis program of time period nonlinear field
JP2016003971A (en) 2014-06-18 2016-01-12 株式会社日立製作所 Magnetic field / circuit coupled field analysis method and circuit coupled analyzer or synchronous rotating machine designed using the same

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014142701A (en) 2013-01-22 2014-08-07 Hitachi Ltd Analysis method of time period nonlinear field and analysis program of time period nonlinear field
JP2016003971A (en) 2014-06-18 2016-01-12 株式会社日立製作所 Magnetic field / circuit coupled field analysis method and circuit coupled analyzer or synchronous rotating machine designed using the same

Also Published As

Publication number Publication date
JP2022033920A (en) 2022-03-02

Similar Documents

Publication Publication Date Title
JP6987304B2 (en) Power converters, machine learners, and how to generate trained models
EP3806311A1 (en) Model predictive pulse pattern control based on small-signal pulse pattern optimization
Putek et al. Shape and topology optimization of a permanent-magnet machine under uncertainties
Stolze Advanced finite-set model predictive control for power electronics and electrical drives
JP6991061B2 (en) Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program and recording medium on which the program is recorded
JP5523265B2 (en) High-speed steady-field analysis method, high-speed steady-field analysis apparatus, steady-state high-speed analysis program, and computer-readable recording medium recording this program
JP7198334B2 (en) Power supply voltage waveform calculation method, circuit coupled magnetic field analysis method, program, and recording medium recording the program
Sobczyk et al. A new approach to steady state analysis of nonlinear electrical circuits
CN103226655B (en) Be suitable for the improvement EMTP algorithm of PWM converter averaging model
Sobczyk et al. Discrete differential operators for periodic and two-periodic time functions
JP7170601B2 (en) Motor simulator and its program
CN117856651B (en) Finite set sliding mode predictive control method and system for single-phase LC filter inverter
JP6650333B2 (en) Frequency characteristic analyzer, method, and program
JP5061175B2 (en) Steady-field fast solution, steady-field fast analysis program, and recording medium
JP5948258B2 (en) Time-period nonlinear field analysis method and time-period nonlinear field analysis program
JP4500155B2 (en) Simulation method
Šmídl et al. Frequency spectrum shaping using finite control set MPC with LQ lookahead
CN107863772B (en) A High-Order Trajectory Sensitivity Calculation Method Based on Generalized Galerkin
Agobert et al. Optimization on frequency constraints with FFT using automatic differentiation on hybrid ODE applications
Nowak et al. Field‐circuit model of the dynamics of electromechanical device supplied by electronic power converters
Iqbal Finite-Time Stability-Based Chaos Suppression of Brushless DC Motor
JP7159068B2 (en) Inductance analysis method, inductance analysis program, and inductance analysis system
Qu et al. Extraction of dynamic low-order nonlinear inductor models based on steady state solutions
Pimentel et al. High-order $ A $-stable and $ L $-stable state-space discrete modeling of continuous systems
Lei Efficient Power System Transient Simulation for Stability Studies Based on Frequency Response Optimized Approximation

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20211207

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20221209

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20221213

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221216

R150 Certificate of patent or registration of utility model

Ref document number: 7198334

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150