JP7428893B2 - Learning device, abnormal eddy current loss estimation device, learning method, abnormal eddy current loss estimation method, and program - Google Patents

Learning device, abnormal eddy current loss estimation device, learning method, abnormal eddy current loss estimation method, and program Download PDF

Info

Publication number
JP7428893B2
JP7428893B2 JP2020073951A JP2020073951A JP7428893B2 JP 7428893 B2 JP7428893 B2 JP 7428893B2 JP 2020073951 A JP2020073951 A JP 2020073951A JP 2020073951 A JP2020073951 A JP 2020073951A JP 7428893 B2 JP7428893 B2 JP 7428893B2
Authority
JP
Japan
Prior art keywords
eddy current
current loss
time
value
abnormal eddy
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
JP2020073951A
Other languages
Japanese (ja)
Other versions
JP2021170290A (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.)
Nippon Steel Corp
Original Assignee
Nippon Steel Corp
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
Application filed by Nippon Steel Corp filed Critical Nippon Steel Corp
Priority to JP2020073951A priority Critical patent/JP7428893B2/en
Publication of JP2021170290A publication Critical patent/JP2021170290A/en
Application granted granted Critical
Publication of JP7428893B2 publication Critical patent/JP7428893B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)

Description

本発明は、学習装置、異常渦電流損推定装置、学習方法、異常渦電流損推定方法、およびプログラムに関し、特に、軟磁性材料における異常渦電流損を推定するために用いて好適なものである。 The present invention relates to a learning device, an abnormal eddy current loss estimation device, a learning method, an abnormal eddy current loss estimation method, and a program, and is particularly suitable for use in estimating abnormal eddy current loss in soft magnetic materials. .

電気機器の高効率化等のために電気機器の損失を推定することは重要である。例えば、モータで発生する損失には、主に、銅損、鉄損、および機械損がある。モータを高速で運転する場合のモータの全損失に占める鉄損の割合が大きくなる。このため、モータを高速で運転しているときのモータの効率を推定するためには、モータの鉄損を高精度に推定することが必要になる。 It is important to estimate the loss of electrical equipment in order to improve its efficiency. For example, losses occurring in motors mainly include copper loss, iron loss, and mechanical loss. When the motor is operated at high speed, the proportion of iron loss in the total loss of the motor increases. Therefore, in order to estimate the efficiency of the motor when the motor is operating at high speed, it is necessary to estimate the core loss of the motor with high accuracy.

鉄損は、ヒステリシス損、古典渦電流損、および異常渦電流損に分類される。高周波の励磁条件では、鉄損に占める古典渦電流損および異常渦電流損の割合が支配的になる。古典渦電流損はファラデーの電磁誘導の法則により理論的に求められる。これに対し、異常渦電流損は、磁壁の移動に伴い発生する誘導電流による損失であるという仮説は立てられているが、異常渦電流損を定量的に評価するための物理モデルは確立されていない。 Iron loss is classified into hysteresis loss, classical eddy current loss, and abnormal eddy current loss. Under high-frequency excitation conditions, the proportions of classical eddy current loss and abnormal eddy current loss in iron loss become dominant. Classical eddy current loss is theoretically determined by Faraday's law of electromagnetic induction. On the other hand, it has been hypothesized that abnormal eddy current loss is a loss due to induced current that occurs due to the movement of domain walls, but a physical model for quantitatively evaluating abnormal eddy current loss has not been established. do not have.

このような背景の下、特許文献1には、軟磁性材料が励磁されたときに当該軟磁性材料に発生する磁束密度の一周期における大きさの最大値を数値解析の結果から取得し、当該磁束密度の最大値に対応する異常渦電流損係数を導出することが記載されている。また、非特許文献1には、電磁鋼板の磁束密度の時間微分値を用いて、異常渦電流損係数を表すことが記載されている。異常渦電流損係数は、全渦電流損を表現する際に古典渦電流損に乗算される係数である。異常渦電流損係数を用いることにより軟磁性材料の異常渦電流損を導出することができる。 Against this background, Patent Document 1 discloses that when the soft magnetic material is excited, the maximum value of the magnitude of the magnetic flux density generated in the soft magnetic material in one cycle is obtained from the results of numerical analysis, and the It is described that an abnormal eddy current loss coefficient corresponding to the maximum value of magnetic flux density is derived. Further, Non-Patent Document 1 describes that the abnormal eddy current loss coefficient is expressed using the time differential value of the magnetic flux density of the electromagnetic steel sheet. The abnormal eddy current loss coefficient is a coefficient by which the classical eddy current loss is multiplied when expressing the total eddy current loss. By using the abnormal eddy current loss coefficient, the abnormal eddy current loss of the soft magnetic material can be derived.

特許第6206608号公報Patent No. 6206608

光岡隆平、外2名、「プレイモデルと有限要素渦電流解析を用いた電磁鋼板の交流ヒステリシスモデルに関する一検討」、電気学会研究会資料 回転機研究会、SA-12-17/RM-12-17、2012年1月26日、p.95-98Ryuhei Mitsuoka and two others, "Study on AC hysteresis model of electrical steel sheets using play model and finite element eddy current analysis", Institute of Electrical Engineers of Japan Research Group Materials Rotating Machine Research Group, SA-12-17/RM-12- 17, January 26, 2012, p. 95-98 中田高義、高橋則雄著、「電気工学の有限要素法」、第2版、森北出版株式会社、1986年4月Takayoshi Nakata and Norio Takahashi, "Finite Element Method for Electrical Engineering", 2nd edition, Morikita Publishing Co., Ltd., April 1986.

しかしながら、特許文献1に記載の技術では、異常渦電流損係数は、正弦波で励磁した磁束密度の測定の結果から求められる。従って、推定対象の軟磁性材料の磁束密度の時間波形と、異常渦電流損係数を求めた際の磁束密度の時間波形とが異なっていても、それらの磁束密度の最大値が同じであれば、異常渦電流損係数として同じ値が導出される。 However, in the technique described in Patent Document 1, the abnormal eddy current loss coefficient is determined from the result of measuring the magnetic flux density excited by a sine wave. Therefore, even if the time waveform of the magnetic flux density of the soft magnetic material to be estimated is different from the time waveform of the magnetic flux density when calculating the abnormal eddy current loss coefficient, if the maximum value of their magnetic flux densities is the same, then , the same value is derived as the abnormal eddy current loss coefficient.

また、特許文献2に記載の技術では、一周期全体の磁束密度の時間微分値の実効値を用いて異常渦電流損係数を導出しており、磁束密度を時系列のデータとして扱っていない。従って、例えば、磁束密度の時間波形の形が同じであれば、当該時間波形にバイアスが重畳していてもいなくても、異常渦電流損係数として同じ値が導出される。
以上のように、特許文献1、2に記載の技術では、例えば、磁束密度の時間波形に高調波成分が含まれる場合やバイアスが重畳している場合には、磁束密度の時間波形に対応する異常渦電流損を正確に導出することが容易ではない。
また、特許文献1、2に記載の技術では、異常渦電流損の一周期における代表値が導出されるが、一周期内の各時刻における異常渦電流損を導出することができない。
Further, in the technique described in Patent Document 2, the abnormal eddy current loss coefficient is derived using the effective value of the time differential value of the magnetic flux density over one period, and the magnetic flux density is not treated as time-series data. Therefore, for example, if the shape of the time waveform of the magnetic flux density is the same, the same value will be derived as the abnormal eddy current loss coefficient whether or not a bias is superimposed on the time waveform.
As described above, in the technologies described in Patent Documents 1 and 2, for example, when the time waveform of the magnetic flux density contains harmonic components or when bias is superimposed, the time waveform of the magnetic flux density corresponds to the time waveform of the magnetic flux density. It is not easy to accurately derive abnormal eddy current loss.
Further, in the techniques described in Patent Documents 1 and 2, a representative value of the abnormal eddy current loss in one period is derived, but it is not possible to derive the abnormal eddy current loss at each time within one period.

本発明は、以上のような問題点に鑑みてなされたものであり、軟磁性材料における異常渦電流損を正確に導出することができるようにすることを目的とする。 The present invention has been made in view of the above problems, and an object of the present invention is to enable accurate derivation of abnormal eddy current loss in soft magnetic materials.

本発明の学習装置は、軟磁性材料の異常渦電流損を推定するためのニューラルネットワークのパラメータを学習する学習装置であって、前記軟磁性材料を所定の励磁条件で励磁したときの当該軟磁性材料に対する測定値として、当該軟磁性材料の磁化特性から得られる物理量の一周期内の各時刻における測定値と、当該軟磁性材料の鉄損の測定値とを取得する取得手段と、前記軟磁性材料を前記所定の励磁条件で励磁したときの前記軟磁性材料のヒステリシス損および古典渦電流損の一周期における代表値を、マックスウェルの方程式に基づいて導出する第1の電磁場解析手段と、前記取得手段により取得された鉄損の測定値と、前記第1の電磁場解析手段により導出されたヒステリシス損および古典渦電流損の一周期における代表値とに基づいて、前記軟磁性材料における異常渦電流損の一周期における代表値を導出する異常渦電流損導出手段と、前記異常渦電流損導出手段により導出された前記軟磁性材料における異常渦電流損の一周期における代表値と、前記物理量の大きさの一周期内の各時刻における値と、前記物理量の大きさの時間微分値の一周期内の各時刻における値と、に基づいて、教師データを作成する教師データ作成手段と、前記教師データ作成手段により作成された前記教師データを用いて、ニューラルネットワークのパラメータを学習する学習手段と、を有し、前記ニューラルネットワークの入力層には、或る時刻における前記物理量の大きさおよび当該時刻における前記物理量の大きさの時間微分値が少なくとも入力され、前記ニューラルネットワークの出力層から、当該時刻における異常渦電流損が少なくとも出力され、前記学習手段は、前記教師データに含まれる前記軟磁性材料の異常渦電流損の一周期における代表値と、前記ニューラルネットワークの前記出力層から出力される一周期内の各時刻における異常渦電流損の代表値との差の評価値を少なくとも求める目的関数を用いて、前記ニューラルネットワークのパラメータを学習することを特徴とする。
本発明の異常渦電流損推定装置は、前記学習装置でパラメータが学習された前記ニューラルネットワークを用いて軟磁性材料の異常渦電流損を推定する異常渦電流損推定装置であって、前記異常渦電流損の推定対象の軟磁性材料を所定の励磁条件で励磁したときの当該軟磁性材料の前記物理量の大きさおよび当該物理量の大きさの時間微分値を導出する導出手段と、前記導出手段により導出された或る時刻における前記物理量の大きさおよび当該時刻における前記物理量の大きさの時間微分値を前記ニューラルネットワークの前記入力層に与えることにより前記出力層から出力される値に基づいて、当該時刻における前記軟磁性材料の異常渦電流損を推定する推定手段と、を有することを特徴とする。
The learning device of the present invention is a learning device for learning parameters of a neural network for estimating abnormal eddy current loss of a soft magnetic material, and the learning device is a learning device for learning parameters of a neural network for estimating abnormal eddy current loss of a soft magnetic material. an acquisition means for acquiring, as measured values for the material, a measured value at each time within one cycle of a physical quantity obtained from the magnetization property of the soft magnetic material, and a measured value of iron loss of the soft magnetic material; a first electromagnetic field analysis means for deriving representative values in one period of hysteresis loss and classical eddy current loss of the soft magnetic material when the material is excited under the predetermined excitation conditions, based on Maxwell's equation; The abnormal eddy current in the soft magnetic material is determined based on the measured value of iron loss acquired by the acquisition means and the representative values of hysteresis loss and classical eddy current loss in one cycle derived by the first electromagnetic field analysis means. abnormal eddy current loss deriving means for deriving a representative value of loss in one cycle; a representative value of abnormal eddy current loss in one cycle of the soft magnetic material derived by the abnormal eddy current loss deriving means; and a magnitude of the physical quantity. a teacher data creating means for creating teacher data based on a value at each time within one cycle of the magnitude of the physical quantity and a value at each time within one cycle of the time derivative of the magnitude of the physical quantity; a learning means for learning parameters of a neural network using the teacher data created by the creation means, and an input layer of the neural network includes the magnitude of the physical quantity at a certain time and the At least a time differential value of the magnitude of the physical quantity is input, at least an abnormal eddy current loss at the relevant time is output from the output layer of the neural network, and the learning means Using an objective function to at least obtain an evaluation value of the difference between a representative value of abnormal eddy current loss in one cycle and a representative value of abnormal eddy current loss at each time in one cycle output from the output layer of the neural network. The method is characterized in that the parameters of the neural network are learned.
The abnormal eddy current loss estimating device of the present invention is an abnormal eddy current loss estimating device that estimates an abnormal eddy current loss of a soft magnetic material using the neural network whose parameters have been learned by the learning device, A derivation means for deriving the magnitude of the physical quantity of the soft magnetic material whose current loss is to be estimated and a time differential value of the magnitude of the physical quantity when the soft magnetic material whose current loss is to be estimated is excited under predetermined excitation conditions; Based on the value output from the output layer by giving the derived magnitude of the physical quantity at a certain time and the time differential value of the magnitude of the physical quantity at the time to the input layer of the neural network, The method is characterized by comprising an estimating means for estimating the abnormal eddy current loss of the soft magnetic material at a given time.

本発明の学習方法は、軟磁性材料の異常渦電流損を推定するためのニューラルネットワークのパラメータを学習する学習方法であって、前記軟磁性材料を所定の励磁条件で励磁したときの当該軟磁性材料に対する測定値として、当該軟磁性材料の磁化特性から得られる物理量の一周期内の各時刻における測定値と、当該軟磁性材料の鉄損の測定値とを取得する取得工程と、前記軟磁性材料を前記所定の励磁条件で励磁したときの前記軟磁性材料のヒステリシス損および古典渦電流損の一周期における代表値を、マックスウェルの方程式に基づいて導出する第1の電磁場解析工程と、前記取得工程により取得された鉄損の測定値と、前記第1の電磁場解析工程により導出されたヒステリシス損および古典渦電流損の一周期における代表値とに基づいて、前記軟磁性材料における異常渦電流損の一周期における代表値を導出する異常渦電流損導出工程と、前記異常渦電流損導出工程により導出された前記軟磁性材料における異常渦電流損の一周期における代表値と、前記物理量の大きさの一周期内の各時刻における値と、前記物理量の大きさの時間微分値の一周期内の各時刻における値と、に基づいて、教師データを作成する教師データ作成工程と、前記教師データ作成工程により作成された前記教師データを用いて、ニューラルネットワークのパラメータを学習する学習工程と、を有し、前記ニューラルネットワークの入力層には、或る時刻における前記物理量の大きさおよび当該時刻における前記物理量の大きさの時間微分値が少なくとも入力され、前記ニューラルネットワークの出力層から、当該時刻における異常渦電流損が少なくとも出力され、前記学習工程は、前記教師データに含まれる前記軟磁性材料の異常渦電流損の一周期における代表値と、前記ニューラルネットワークの前記出力層から出力される一周期内の各時刻における異常渦電流損の代表値との差の評価値を少なくとも求める目的関数を用いて、前記ニューラルネットワークのパラメータを学習することを特徴とする。
本発明の異常渦電流損推定方法は、前記学習方法でパラメータが学習された前記ニューラルネットワークを用いて軟磁性材料の異常渦電流損を推定する異常渦電流損推定方法であって、前記異常渦電流損の推定対象の軟磁性材料を所定の励磁条件で励磁したときの当該軟磁性材料の前記物理量の大きさおよび当該物理量の大きさの時間微分値を導出する導出工程と、前記導出工程により導出された或る時刻における前記物理量の大きさおよび当該時刻における前記物理量の大きさの時間微分値を前記ニューラルネットワークの前記入力層に与えることにより前記出力層から出力される値に基づいて、当該時刻における前記軟磁性材料の異常渦電流損を推定する推定工程と、を有することを特徴とする。
The learning method of the present invention is a learning method for learning parameters of a neural network for estimating abnormal eddy current loss of a soft magnetic material, which an acquisition step of acquiring, as measured values for the material, measured values at each time within one cycle of a physical quantity obtained from the magnetization characteristics of the soft magnetic material and a measured value of iron loss of the soft magnetic material; a first electromagnetic field analysis step of deriving representative values in one period of hysteresis loss and classical eddy current loss of the soft magnetic material when the material is excited under the predetermined excitation conditions, based on Maxwell's equation; The abnormal eddy current in the soft magnetic material is calculated based on the measured value of iron loss obtained in the acquisition step and the representative values in one cycle of hysteresis loss and classical eddy current loss derived in the first electromagnetic field analysis step. An abnormal eddy current loss derivation step for deriving a representative value in one cycle of loss, a typical value in one cycle of abnormal eddy current loss in the soft magnetic material derived by the abnormal eddy current loss derivation step, and the magnitude of the physical quantity. a teacher data creation step of creating teacher data based on a value at each time within one cycle of the magnitude of the physical quantity and a value at each time within one cycle of the time derivative of the magnitude of the physical quantity; a learning step of learning parameters of a neural network using the teacher data created in the creation step, and an input layer of the neural network includes the magnitude of the physical quantity at a certain time and the At least a time differential value of the magnitude of the physical quantity is input, at least the abnormal eddy current loss at the relevant time is output from the output layer of the neural network, and the learning step Using an objective function to at least obtain an evaluation value of the difference between a representative value of abnormal eddy current loss in one cycle and a representative value of abnormal eddy current loss at each time in one cycle output from the output layer of the neural network. The method is characterized in that the parameters of the neural network are learned.
The abnormal eddy current loss estimation method of the present invention is an abnormal eddy current loss estimation method for estimating abnormal eddy current loss of a soft magnetic material using the neural network whose parameters have been learned by the learning method, A derivation step of deriving the magnitude of the physical quantity of the soft magnetic material and the time differential value of the magnitude of the physical quantity when the soft magnetic material whose current loss is to be estimated is excited under predetermined excitation conditions; and the derivation step Based on the value output from the output layer by giving the derived magnitude of the physical quantity at a certain time and the time differential value of the magnitude of the physical quantity at the time to the input layer of the neural network, The method is characterized by comprising an estimating step of estimating an abnormal eddy current loss of the soft magnetic material at a given time.

本発明のプログラムの第1の例は、前記学習装置の各手段としてコンピュータを機能させるためのものである。
本発明のプログラムの第2の例は、前記異常渦電流損推定装置の各手段としてコンピュータを機能させるためのものである。
A first example of the program of the present invention is for causing a computer to function as each means of the learning device.
A second example of the program of the present invention is for causing a computer to function as each means of the abnormal eddy current loss estimating device.

本発明によれば、軟磁性材料における異常渦電流損を正確に導出することができる。 According to the present invention, it is possible to accurately derive abnormal eddy current loss in a soft magnetic material.

異常渦電流損推定装置の機能的な構成の一例を示す図である。It is a figure showing an example of the functional composition of an abnormal eddy current loss estimating device. 磁束密度の瞬時値の時間波形と、磁束密度の大きさの時間波形と、磁束密度の大きさの時間微分値の時間波形の一例を示す図である。FIG. 3 is a diagram showing an example of a time waveform of an instantaneous value of magnetic flux density, a time waveform of the magnitude of magnetic flux density, and a time waveform of a time differential value of the magnitude of magnetic flux density. ニューラルネットワークの一例を概念的に示す図である。FIG. 2 is a diagram conceptually showing an example of a neural network. ニューラルネットワークのパラメータの学習方法の一例を説明するフローチャートである。3 is a flowchart illustrating an example of a method for learning parameters of a neural network. 異常渦電流損推定方法の一例を説明するフローチャートである。It is a flow chart explaining an example of an abnormal eddy current loss estimation method. 発明例および比較例における異常渦電流損の推定値と、異常渦電流損の測定値とを比較した結果を示す図である。It is a figure which shows the result of comparing the estimated value of abnormal eddy current loss and the measured value of abnormal eddy current loss in an invention example and a comparative example.

<着想>
本発明の実施形態を説明する前に、本発明者が得た着想について説明する。
異常渦電流損は磁壁の移動に伴って発生する渦電流損である。強磁性体の磁化過程では、まず、磁壁の移動が発生し、且つ、180°磁壁が90°磁壁よりも先に拡大する。このような磁壁の移動の後に磁壁内の磁化方向が変化する過程を経る。そのため、異常渦電流損の発生の程度は磁化状態により変化すると考えられ、また、その大きさは、磁束密度の大きさと、磁束密度の時間変化とによって定まると考えられる。また、磁界強度は磁束密度と対応するものであるため、異常渦電流損の大きさは、磁界強度の大きさと、磁界強度の時間変化によっても定まると考えられる。また、透磁率および微分透磁率も磁化状態により変化するため、異常渦電流損の大きさは、透磁率の大きさと透磁率の時間変化、および、微分透磁率の大きさと微分透磁率の時間変化によっても定まると考えらえる。尚、透磁率は、磁化特性において相互に対応する磁束密度と磁界強度との比(磁束密度を磁界強度で割った値)である。微分透磁率は、磁化特性曲線上の、磁束密度に応じて定まる点における接線の傾きである。
<Idea>
Before describing embodiments of the present invention, an idea obtained by the present inventor will be described.
Abnormal eddy current loss is eddy current loss that occurs due to movement of domain walls. In the magnetization process of a ferromagnetic material, first, movement of the domain wall occurs, and the 180° domain wall expands before the 90° domain wall. After such movement of the domain wall, a process occurs in which the direction of magnetization within the domain wall changes. Therefore, the degree of occurrence of abnormal eddy current loss is considered to vary depending on the magnetization state, and the magnitude is considered to be determined by the magnitude of the magnetic flux density and the temporal change in the magnetic flux density. Furthermore, since the magnetic field strength corresponds to the magnetic flux density, the magnitude of the abnormal eddy current loss is considered to be determined by the magnitude of the magnetic field strength and the temporal change in the magnetic field strength. In addition, since magnetic permeability and differential permeability also change depending on the magnetization state, the magnitude of abnormal eddy current loss is determined by the magnitude of magnetic permeability and time change in permeability, as well as the magnitude of differential permeability and time change in differential permeability. It can be thought that it is also determined by Note that magnetic permeability is the ratio of magnetic flux density and magnetic field strength that correspond to each other in magnetization characteristics (the value obtained by dividing the magnetic flux density by the magnetic field strength). Differential magnetic permeability is the slope of the tangent at a point on the magnetization characteristic curve that is determined according to the magnetic flux density.

磁束密度、磁界強度、透磁率、および微分透磁率は、何れも、磁化特性から得られる物理量であり、磁化特性を表現するパラメータの1つとなる。磁化特性は、磁束密度と磁界強度との関係を表し、磁気ヒステリシス特性とも称される。磁化特性曲線は、磁化特性を表す曲線である。 Magnetic flux density, magnetic field strength, magnetic permeability, and differential magnetic permeability are all physical quantities obtained from magnetization characteristics, and are one of the parameters expressing magnetization characteristics. Magnetization characteristics represent the relationship between magnetic flux density and magnetic field strength, and are also referred to as magnetic hysteresis characteristics. The magnetization characteristic curve is a curve representing magnetization characteristics.

背景技術の欄で説明したように、異常渦電流損は、磁壁の移動に伴い発生する誘導電流による損失であるという仮説は立てられているが、異常渦電流損を定量的に評価するための物理モデルは確立されていない。そこで、本発明者は、前述したように、磁化特性から得られる物理量の大きさおよび時間微分値と、異常渦電流損の大きさとが対応することに着目した。また、本発明者は、異常渦電流損は各時刻において発生することにも着目した。そして、鋭意検討した結果、本発明者は、磁化特性から得られる物理量の大きさと当該物理量の大きさの時間微分値とのそれぞれの或る時刻における値(瞬時値)を説明変数とし、異常渦電流損の当該時刻における値を目的変数とするニューラルネットワークを構築すれば、異常渦電流損を定量的に評価するための物理モデルを用いずに異常渦電流損を正確に推定することができることに想到した。
以下の実施形態は、以上の着想に基づいてなされたものである。
As explained in the background technology section, it is hypothesized that abnormal eddy current loss is a loss due to induced current that occurs due to movement of domain walls, but there is no way to quantitatively evaluate abnormal eddy current loss. A physical model has not been established. Therefore, as described above, the present inventors focused on the fact that the magnitude and time differential value of the physical quantity obtained from the magnetization characteristics correspond to the magnitude of the abnormal eddy current loss. The present inventor also noted that abnormal eddy current loss occurs at each time. As a result of intensive study, the inventors determined that the magnitude of the physical quantity obtained from the magnetization characteristics and the time differential value of the magnitude of the physical quantity at a certain time (instantaneous value) were used as explanatory variables, and the abnormal vortex By constructing a neural network that uses the current loss value at the relevant time as the objective variable, it is possible to accurately estimate abnormal eddy current loss without using a physical model to quantitatively evaluate abnormal eddy current loss. I came up with the idea.
The following embodiments have been made based on the above idea.

以下、図面を参照しながら、本発明の一実施形態を説明する。
<異常渦電流損推定装置の構成>
図1は、異常渦電流損推定装置100の機能的な構成の一例を示す図である。異常渦電流損推定装置100のハードウェアの構成は、例えば、CPU、ROM、RAM、HDD、および各種のインターフェースを有する情報処理装置、または、専用のハードウェアを用いることにより実現される。以下に本実施形態の異常渦電流損推定装置100が有する機能の一例を説明する。尚、本実施形態では、軟磁性材料が電磁鋼板である場合を例示する。電磁鋼板は、例えば、方向性電磁鋼板または無方向性電磁鋼板である。
Hereinafter, one embodiment of the present invention will be described with reference to the drawings.
<Configuration of abnormal eddy current loss estimation device>
FIG. 1 is a diagram showing an example of a functional configuration of an abnormal eddy current loss estimating device 100. The hardware configuration of the abnormal eddy current loss estimating device 100 is realized by using, for example, an information processing device having a CPU, ROM, RAM, HDD, and various interfaces, or dedicated hardware. An example of the functions of the abnormal eddy current loss estimating device 100 of this embodiment will be described below. In addition, in this embodiment, the case where a soft magnetic material is an electromagnetic steel plate is illustrated. The electrical steel sheet is, for example, a grain-oriented electrical steel sheet or a non-oriented electrical steel sheet.

<<取得部101>>
取得部101は、電磁鋼板の磁化特性から得られる物理量の一周期内の各時刻における測定値と、電磁鋼板の鉄損の測定値とを取得する。これらの測定値は、電磁鋼板を同一の励磁条件で励磁したときの測定値である。励磁条件とは、励磁電流の時間波形を定める情報である。励磁電流の時間波形が正弦波である場合、かかる情報は、例えば、実効値と周波数である。また、一周期とは、励磁電流の時間波形の一周期である。また、一周期内の各時刻は、例えば、周期の最初の時刻を基準値(例えば0)とした場合の時刻であり、現実世界における時刻とは必ずしも一致しない。このことは、以降の説明でも同じである。
<<Acquisition unit 101>>
The acquisition unit 101 acquires a measured value at each time within one cycle of a physical quantity obtained from the magnetization characteristics of the electrical steel sheet, and a measured value of iron loss of the electrical steel sheet. These measured values were measured when the electromagnetic steel sheet was excited under the same excitation conditions. The excitation condition is information that determines the time waveform of the excitation current. When the time waveform of the excitation current is a sine wave, such information is, for example, an effective value and a frequency. Moreover, one period is one period of the time waveform of the excitation current. Further, each time within one cycle is, for example, a time when the first time of the cycle is set as a reference value (for example, 0), and does not necessarily match the time in the real world. This also applies to the following explanations.

前述したように、磁化特性から得られる物理量は、磁束密度、磁界強度、透磁率、および微分透磁率である。磁束密度、磁界強度、透磁率、および微分透磁率の少なくとも1つを、磁化特性から得られる物理量として用いることができるが、本実施形態では、磁束密度が、磁化特性から得られる物理量である場合を例示する。 As described above, the physical quantities obtained from the magnetization characteristics are magnetic flux density, magnetic field strength, magnetic permeability, and differential magnetic permeability. At least one of magnetic flux density, magnetic field strength, magnetic permeability, and differential magnetic permeability can be used as a physical quantity obtained from magnetization characteristics, but in this embodiment, when magnetic flux density is a physical quantity obtained from magnetization characteristics exemplify.

電磁鋼板の磁束密度の測定値と鉄損の測定値は、例えば、JIS C 2556:2015「単板試験器による電磁鋼帯の磁気特性の測定方法」や、JIS 2550-1:2011「電磁鋼帯試験方法-第1部:エプスタイン試験器による電磁鋼帯の磁気特性の測定方法」に記載されている公知の手法により得ることができるので、ここでは、その詳細な説明を省略する。 The measured value of magnetic flux density and the measured value of iron loss of electrical steel sheets are, for example, JIS C 2556:2015 "Method for measuring the magnetic properties of electrical steel strip using a veneer tester" and JIS 2550-1:2011 "Method of measuring magnetic properties of electrical steel strip". Since it can be obtained by the known method described in "Strand Test Method - Part 1: Method for Measuring Magnetic Characteristics of Electromagnetic Steel Strip Using Epstein Tester," detailed explanation thereof will be omitted here.

尚、鉄損の測定値は、磁束密度の一周期内の各時刻における測定値を用いて一周期において1つだけ測定されるものである。即ち、鉄損の測定値は、一周期における代表値(平均的な値)である。また、取得部101で取得される各測定値(磁束密度の一周期内の各時刻における測定値と、鉄損の測定値)は、電磁鋼板全体に対する測定値である。従って、或る電磁鋼板を或る励磁条件で励磁した場合に、磁束密度の一周期内の各時刻における測定値と、鉄損の測定値との組が1つ得られる。 Note that the measured value of iron loss is measured only once in one period using the measured value at each time within one period of magnetic flux density. That is, the measured value of iron loss is a representative value (average value) in one cycle. Further, each measurement value (measurement value at each time within one cycle of magnetic flux density and measurement value of iron loss) acquired by the acquisition unit 101 is a measurement value for the entire electrical steel sheet. Therefore, when a certain electromagnetic steel sheet is excited under certain excitation conditions, one set of the measured value at each time within one period of the magnetic flux density and the measured value of iron loss is obtained.

以下の説明では、同一の種類の電磁鋼板を同一の励磁条件で励磁した場合の、磁束密度の一周期内の各時刻における測定値と、鉄損の測定値との組を、必要に応じて測定データと称する。尚、測定データには、当該測定データを得るために用いた励磁条件を示す情報と、当該測定データにおける測定対象の電磁鋼板の種類を示す情報も含まれている。電磁鋼板の種類とは、電磁鋼板の属性に応じた電磁鋼板の分類先を特定する情報である。電磁鋼板の属性は、例えば、電磁鋼板に含まれる成分と当該成分の量である。 In the following explanation, pairs of measured values at each time within one cycle of magnetic flux density and measured values of iron loss when the same type of electrical steel sheet is excited under the same excitation conditions are explained as necessary. This is called measurement data. Note that the measurement data also includes information indicating the excitation conditions used to obtain the measurement data and information indicating the type of electromagnetic steel sheet to be measured in the measurement data. The type of electrical steel sheet is information that specifies the classification of the electrical steel sheet according to the attributes of the electrical steel sheet. The attributes of the electrical steel sheet include, for example, the components contained in the electrical steel sheet and the amounts of the components.

異常渦電流損の推定対象の電磁鋼板と同じ種類の電磁鋼板を、相互に異なる複数の励磁条件で励磁することにより、それぞれの励磁条件に対応する測定データを得る。取得部101は、このようにして得られる複数の測定データを取得する。また、取得部101は、複数種類の電磁鋼板についての測定データを取得してもよい。
測定データの取得の形態として、例えば、外部装置との通信、可搬型記憶媒体に記憶された測定データの読み出し、および異常渦電流損推定装置100のユーザインターフェースに対する測定データの値の入力操作による入力のうち、少なくとも1つを採用することができる。
By exciting an electromagnetic steel sheet of the same type as the electromagnetic steel sheet whose abnormal eddy current loss is to be estimated under a plurality of mutually different excitation conditions, measurement data corresponding to each excitation condition is obtained. The acquisition unit 101 acquires a plurality of measurement data obtained in this way. Further, the acquisition unit 101 may acquire measurement data regarding multiple types of electrical steel sheets.
Forms of acquiring measurement data include, for example, communication with an external device, reading measurement data stored in a portable storage medium, and inputting measurement data values to the user interface of the abnormal eddy current loss estimating device 100. At least one of these can be adopted.

<<第1の電磁場解析部102>>
第1の電磁場解析部102は、電磁鋼板を所定の励磁条件で励磁したときの電磁鋼板のヒステリシス損および古典渦電流損の一周期における代表値を、マックスウェルの方程式に基づいて導出する。第1の電磁場解析部102の処理の具体例を以下に説明する。
<<First electromagnetic field analysis unit 102>>
The first electromagnetic field analysis unit 102 derives representative values of hysteresis loss and classical eddy current loss in one cycle of the electromagnetic steel sheet when the electromagnetic steel sheet is excited under predetermined excitation conditions based on Maxwell's equation. A specific example of the processing of the first electromagnetic field analysis unit 102 will be described below.

第1の電磁場解析部102は、取得部101により取得された測定データに含まれる励磁条件を含む電磁場解析条件に従って、マックスウェルの方程式に基づき、有限要素法を用いて、各微小領域(メッシュ)における磁束密度ベクトルBと渦電流ベクトルJeを一周期内の各時刻において計算する。磁束密度ベクトルBと渦電流ベクトルJeは、予め定められた時間ステップ毎に到来する計算のタイミングにおいて計算される。時間ステップが短ければ、細かい時間ステップで値を導出することができる一方で処理の負荷が増大する。時間ステップは、このようなトレードオフを考慮して予め定められる。 The first electromagnetic field analysis section 102 analyzes each minute region (mesh) using the finite element method based on Maxwell's equations, in accordance with the electromagnetic field analysis conditions including the excitation conditions included in the measurement data acquired by the acquisition section 101. The magnetic flux density vector B and the eddy current vector J e are calculated at each time within one period. The magnetic flux density vector B and the eddy current vector J e are calculated at calculation timings that arrive at every predetermined time step. If the time steps are short, values can be derived at fine time steps, but the processing load increases. The time step is predetermined in consideration of such trade-offs.

尚、各微小領域における磁束密度ベクトルBと渦電流ベクトルJeを計算することができれば、有限要素法以外の方法(差分法等)を用いて電磁場解析を行ってもよい。
磁束密度ベクトルBと渦電流ベクトルJeを計算するための基礎方程式は、一般に、以下の(1)式~(4)式で与えられる。
Note that as long as the magnetic flux density vector B and the eddy current vector J e in each minute region can be calculated, the electromagnetic field analysis may be performed using a method other than the finite element method (such as the finite difference method).
The basic equations for calculating the magnetic flux density vector B and the eddy current vector J e are generally given by the following equations (1) to (4).

Figure 0007428893000001
Figure 0007428893000001

(1)式~(4)式において、μは、透磁率[H/m]であり、Aは、ベクトルポテンシャル[T・m]であり、σは、導電率[S/m]であり、J0は、励磁電流密度[A/m2]であり、φは、スカラーポテンシャル[V]である。
(1)式および(2)式を連立して解いて、ベクトルポテンシャルAとスカラーポテンシャルφを求めた後、(3)式、(4)式から、磁束密度ベクトルBと渦電流ベクトルJeが計算される。
In equations (1) to (4), μ is magnetic permeability [H/m], A is vector potential [T・m], and σ is electrical conductivity [S/m], J 0 is the excitation current density [A/m 2 ], and φ is the scalar potential [V].
After solving equations (1) and (2) simultaneously to obtain vector potential A and scalar potential φ, from equations (3) and (4), magnetic flux density vector B and eddy current vector J e can be calculated. Calculated.

ベクトルポテンシャルAおよびスカラーポテンシャルφを未知変数とした解法として、例えば、ニュートンラプソン法(Newton-Raphson method)がある。
有限要素法により電磁場の解析を実行する手法は、非特許文献2等に詳細に記載されているように、一般的な手法であるので、ここでは、その詳細な説明を省略する。
As a solution method using vector potential A and scalar potential φ as unknown variables, there is, for example, the Newton-Raphson method.
The method of performing electromagnetic field analysis using the finite element method is a general method, as described in detail in Non-Patent Document 2, and therefore a detailed explanation thereof will be omitted here.

第1の電磁場解析部102は、或る微小領域における磁束密度ベクトルBと当該磁束密度ベクトルBに対応する磁界強度ベクトルHの一周期における時間波形から、当該微小領域における磁化特性(磁気ヒステリシス特性)を導出する。そして、第1の電磁場解析部102は、以下の(5)式により、当該導出した磁化特性の面積を、当該微小領域におけるヒステリシス損Whとして導出する。 The first electromagnetic field analysis unit 102 calculates the magnetization characteristic (magnetic hysteresis characteristic) in a certain minute region from the time waveform in one period of the magnetic flux density vector B in a certain minute region and the magnetic field strength vector H corresponding to the magnetic flux density vector B. Derive. Then, the first electromagnetic field analysis unit 102 derives the area of the derived magnetization characteristic as the hysteresis loss W h in the minute region using the following equation (5).

Figure 0007428893000002
Figure 0007428893000002

第1の電磁場解析部102は、以上のような計算を、全ての微小領域について実行する。第1の電磁場解析部102は、全ての微小領域におけるヒステリシス損Whの総和を、電磁鋼板のヒステリシス損の一周期における代表値として導出する。この電磁鋼板のヒステリシス損は、電磁鋼板全体に対するヒステリシス損である。 The first electromagnetic field analysis unit 102 executes the above calculations for all micro regions. The first electromagnetic field analysis unit 102 derives the sum of the hysteresis loss W h in all the micro regions as a representative value of the hysteresis loss of the electrical steel sheet in one cycle. This hysteresis loss of the electrical steel sheet is a hysteresis loss for the entire electrical steel sheet.

また、第1の電磁場解析部102は、或る微小領域における渦電流ベクトルJeと、当該微小領域の大きさと、電磁鋼板の導電率σに基づいて、以下の(6)式により、当該微小領域における古典渦電流損Weを導出する。 In addition, the first electromagnetic field analysis unit 102 calculates the micro-region based on the eddy current vector J e in a certain micro-region, the size of the micro-region, and the electrical conductivity σ of the electromagnetic steel sheet using equation (6) below. Derive the classical eddy current loss W e in the region.

Figure 0007428893000003
Figure 0007428893000003

(6)式において、Tは、渦電流ベクトルJeの周期である。
第1の電磁場解析部102は、以上のような計算を、全ての微小領域について実行する。第1の電磁場解析部102は、全ての微小領域における古典渦電流損Weの総和を、電磁鋼板の古典渦電流損の一周期における代表値として導出する。この電磁鋼板の古典渦電流損は、電磁鋼板全体に対する古典渦電流損である。
In equation (6), T is the period of the eddy current vector J e .
The first electromagnetic field analysis unit 102 performs the above calculations for all micro regions. The first electromagnetic field analysis unit 102 derives the sum total of the classical eddy current loss W e in all minute regions as a representative value of the classical eddy current loss of the electrical steel sheet in one cycle. The classical eddy current loss of this electrical steel sheet is the classical eddy current loss for the entire electrical steel sheet.

第1の電磁場解析部102は、以上のようにして電磁鋼板のヒステリシス損および古典渦電流損の一周期における代表値を導出することを、測定データに含まれる励磁条件のそれぞれについて、電磁鋼板の種類ごとに実行する。
尚、第1の電磁場解析部102が、電磁鋼板のヒステリシス損および古典渦電流損の一周期における代表値を導出する手法は、(5)式および(6)式を用いる手法に限定されず、その他の公知の手法を採用してもよい。例えば、スタインメッツ係数を用いたスタインメッツの式を用いて電磁鋼板のヒステリシス損および古典渦電流損の一周期における代表値を導出してもよい。
以下の説明では、以上のようにして第1の電磁場解析部102により導出される、電磁鋼板のヒステリシス損および古典渦電流損の一周期における代表値を、必要に応じて、ヒステリシス損および古典渦電流損の一周期における代表値と略称する。
The first electromagnetic field analysis unit 102 calculates the representative values of the hysteresis loss and classical eddy current loss of the electrical steel sheet in one cycle for each of the excitation conditions included in the measurement data in the manner described above. Execute by type.
Note that the method by which the first electromagnetic field analysis unit 102 derives the representative values of the hysteresis loss and classical eddy current loss of the electrical steel sheet in one cycle is not limited to the method using equations (5) and (6), Other known methods may also be employed. For example, the Steinmetz equation using the Steinmetz coefficient may be used to derive the representative values of the hysteresis loss and classical eddy current loss in one cycle of the electrical steel sheet.
In the following explanation, representative values of hysteresis loss and classical eddy current loss in one cycle of an electrical steel sheet, which are derived by the first electromagnetic field analysis unit 102 as described above, will be described as needed. It is abbreviated as the representative value of current loss in one cycle.

<<異常渦電流損導出部103>>
異常渦電流損導出部103は、取得部101により取得された測定データに含まれる鉄損の測定値と、第1の電磁場解析部102により導出されたヒステリシス損および古典渦電流損の一周期における代表値とに基づいて、電磁鋼板における異常渦電流損の一周期における代表値を導出する。以下に、異常渦電流損導出部103の処理の具体例を説明する。前述したように取得部101により取得された測定データに含まれる鉄損の測定値は、一周期における平均的な値である。
<<Abnormal eddy current loss deriving unit 103>>
The abnormal eddy current loss deriving unit 103 calculates the iron loss measurement value included in the measurement data acquired by the acquisition unit 101 and the hysteresis loss and classical eddy current loss in one cycle derived by the first electromagnetic field analysis unit 102. Based on the representative values, a representative value in one cycle of abnormal eddy current loss in the electrical steel sheet is derived. A specific example of the processing of the abnormal eddy current loss deriving unit 103 will be described below. As described above, the measured value of iron loss included in the measurement data acquired by the acquisition unit 101 is an average value in one cycle.

異常渦電流損導出部103は、取得部101により取得された測定データの1つを選択する。異常渦電流損導出部103は、選択した測定データに含まれる鉄損の測定値を抽出する。異常渦電流損導出部103は、第1の電磁場解析部102により導出されたヒステリシス損および古典渦電流損の一周期における代表値のうち、選択した測定データに対応するヒステリシス損および古典渦電流損の一周期における代表値を抽出する。測定データに対応するヒステリシス損および古典渦電流損の一周期における代表値とは、当該測定データに含まれる電磁鋼板の種類と同じ種類の電磁鋼板に対して当該測定データに含まれる励磁条件と同じ励磁条件で励磁した場合のヒステリシス損および古典渦電流損の一周期における代表値である。 The abnormal eddy current loss deriving unit 103 selects one of the measurement data acquired by the acquisition unit 101. The abnormal eddy current loss deriving unit 103 extracts the measured value of iron loss included in the selected measurement data. The abnormal eddy current loss deriving unit 103 calculates the hysteresis loss and classical eddy current loss corresponding to the selected measurement data from among the representative values in one cycle of the hysteresis loss and classical eddy current loss derived by the first electromagnetic field analysis unit 102. Extract the representative value in one period. The typical values of hysteresis loss and classical eddy current loss in one cycle corresponding to the measurement data are the same as the excitation conditions included in the measurement data for the same type of electrical steel sheet as the type of electrical steel sheet included in the measurement data. These are representative values for one period of hysteresis loss and classical eddy current loss when excited under excitation conditions.

異常渦電流損導出部103は、抽出したヒステリシス損の一周期における代表値と、抽出した古典渦電流損の一周期における代表値とを加算して、抽出したヒステリシス損および古典渦電流損の一周期における代表値の和を導出する。
異常渦電流損導出部103は、抽出した鉄損の測定値から、抽出したヒステリシス損および古典渦電流損の一周期における代表値の和を減算した値を、電磁鋼板の異常渦電流損の一周期における代表値として導出する。
The abnormal eddy current loss deriving unit 103 adds the representative value of the extracted hysteresis loss in one cycle and the representative value of the extracted classical eddy current loss in one cycle, and calculates the extracted hysteresis loss and classical eddy current loss. Derive the sum of representative values in the period.
The abnormal eddy current loss deriving unit 103 calculates the value obtained by subtracting the sum of the representative values of the extracted hysteresis loss and classical eddy current loss in one cycle from the extracted measured value of iron loss, as the abnormal eddy current loss of the electrical steel sheet. Derived as a representative value in the period.

異常渦電流損導出部103は、以上のようにして電磁鋼板の異常渦電流損の一周期における代表値を導出することを、測定データに含まれる励磁条件のそれぞれについて、電磁鋼板の種類ごとに実行する。電磁鋼板の異常渦電流損の一周期における代表値は、電磁鋼板全体に対する異常渦電流損である。以下の説明では、以上のようにして異常渦電流損導出部103により導出される、電磁鋼板の異常渦電流損の一周期における代表値を、必要に応じて、学習用目的変数データと称する。 The abnormal eddy current loss deriving unit 103 derives the representative value of the abnormal eddy current loss in one cycle of the electrical steel sheet as described above for each type of electrical steel sheet for each excitation condition included in the measurement data. Execute. The typical value of the abnormal eddy current loss of the electrical steel sheet in one cycle is the abnormal eddy current loss for the entire electrical steel sheet. In the following description, the representative value of the abnormal eddy current loss of the electrical steel sheet in one cycle, which is derived by the abnormal eddy current loss deriving unit 103 as described above, will be referred to as learning objective variable data as necessary.

<<教師データ作成部104>>
教師データ作成部104は、異常渦電流損導出部103により導出された学習用目的変数データ(異常渦電流損の一周期における代表値)と、電磁鋼板の磁束密度の大きさの一周期内の各時刻における値と、電磁鋼板の磁束密度の大きさの時間微分値の一周期内の各時刻における値と、に基づいて、教師データを作成する。ここで、磁束密度の大きさは、絶対値であるものとする。また、電磁鋼板の磁束密度の大きさの一周期内の各時刻における値と、電磁鋼板の磁束密度の大きさの時間微分値の一周期内の各時刻における値は、取得部101により取得された測定データに含まれる磁束密度の一周期内の各時刻における測定値に基づいて得られる。以下に、教師データ作成部104の処理の具体例を説明する。
<<Teacher data creation unit 104>>
The teacher data creation unit 104 uses the learning objective variable data (representative value in one cycle of the abnormal eddy current loss) derived by the abnormal eddy current loss derivation unit 103 and the magnitude of the magnetic flux density of the electrical steel sheet within one cycle. Teacher data is created based on the value at each time and the value at each time within one period of the time differential value of the magnitude of the magnetic flux density of the electromagnetic steel sheet. Here, the magnitude of the magnetic flux density is assumed to be an absolute value. Further, the value at each time within one cycle of the magnitude of the magnetic flux density of the electrical steel sheet and the value at each time within one cycle of the time differential value of the magnitude of the magnetic flux density of the electrical steel sheet are acquired by the acquisition unit 101. It is obtained based on the measured value at each time within one cycle of the magnetic flux density included in the measured data. A specific example of the processing of the teacher data creation unit 104 will be described below.

教師データ作成部104は、取得部101により取得された測定データの1つを選択する。教師データ作成部104は、選択した測定データに含まれる磁束密度の一周期内の各時刻における測定値に基づいて、磁束密度の大きさの一周期内の各時刻における値を導出する。例えば、磁束密度の一周期内の各時刻における測定値を、磁束密度の大きさの一周期内の各時刻における値としてそのまま用いてもよい。また、例えば、磁束密度の一周期内の各時刻における測定値に対してスムージング処理(平滑化処理)を行ってノイズを除去し、スムージング処理が行われた磁束密度の一周期内の各時刻における値を、磁束密度の大きさの一周期内の各時刻における値としてもよい。 The teacher data creation unit 104 selects one of the measurement data acquired by the acquisition unit 101. The teacher data creation unit 104 derives the value of the magnetic flux density at each time within one cycle based on the measured value at each time within one cycle of the magnetic flux density included in the selected measurement data. For example, the measured value at each time within one cycle of magnetic flux density may be used as is as the value at each time within one cycle of the magnitude of magnetic flux density. In addition, for example, smoothing processing (smoothing processing) is performed on the measured values at each time within one cycle of magnetic flux density to remove noise, and The value may be a value at each time within one cycle of the magnitude of the magnetic flux density.

また、教師データ作成部104は、選択した測定データに含まれる磁束密度の一周期内の各時刻における測定値に基づいて、磁束密度の大きさの時間微分値の一周期内の各時刻における値を導出する。前述したように、例えば、スムージング処理が行われた磁束密度の一周期内の各時刻における値に対して時間微分値を導出してもよい。 Further, the teacher data creation unit 104 generates a value at each time within one cycle of the time differential value of the magnitude of the magnetic flux density, based on the measured value at each time within one cycle of the magnetic flux density included in the selected measurement data. Derive. As described above, for example, a time differential value may be derived for the value at each time within one cycle of the magnetic flux density that has been subjected to the smoothing process.

図2は、磁束密度の瞬時値の時間波形と、磁束密度の大きさの時間波形と、磁束密度の大きさの時間微分値の時間波形の一例を示す図である。図2において、時刻t1から時刻t2までの期間Tが一周期である。
時間波形201は、磁束密度の測定値から得られる磁束密度の瞬時値の時間波形の一例である。時間波形202は、時間波形201から得られる磁束密度の大きさの時間波形の一例である。時間波形202は、時間波形201の各時刻の値の絶対値をとることにより得られる。時間波形203は、時間波形201から得られる磁束密度の大きさの時間微分値の時間波形の一例である。時間波形203は、時間波形202を時間微分することにより得られる。尚、時間微分は、時間に関する一階微分のことを指す。時間微分は、数値微分により実現される。数値微分の手法は、前進差分、後退差分、または中心差分の何れの手法であってもよい。
FIG. 2 is a diagram showing an example of a time waveform of an instantaneous value of magnetic flux density, a time waveform of the magnitude of magnetic flux density, and a time waveform of a time differential value of the magnitude of magnetic flux density. In FIG. 2, the period T from time t1 to time t2 is one cycle.
A time waveform 201 is an example of a time waveform of an instantaneous value of magnetic flux density obtained from a measured value of magnetic flux density. The time waveform 202 is an example of the time waveform of the magnitude of the magnetic flux density obtained from the time waveform 201. The time waveform 202 is obtained by taking the absolute value of each time value of the time waveform 201. The time waveform 203 is an example of a time waveform of the time differential value of the magnitude of the magnetic flux density obtained from the time waveform 201. The time waveform 203 is obtained by time-differentiating the time waveform 202. Note that time differentiation refers to first-order differentiation with respect to time. Time differentiation is realized by numerical differentiation. The numerical differentiation method may be a forward difference method, a backward difference method, or a central difference method.

磁束密度の瞬時値の時間波形201における磁束密度が0の時刻t3、t4、t5においては、磁束密度の大きさの時間波形202は微分不可能となる。そこで、このような時刻t3、t3、t5における、磁束密度の大きさの時間微分値として、例えば、当該時刻t3、t4、t5対して時間的に前後する2つの時刻における、磁束密度の大きさの時間微分値の算術平均値を採用することができる。また、このような時刻t3、t4、t5においては、磁束密度の大きさの時間微分値がないものとしてもよい。
以下の説明では、以上のようにして教師データ作成部104により導出される、磁束密度の大きさの一周期内の各時刻における値と、磁束密度の大きさの時間微分値の一周期内の各時刻における値とを、必要に応じて、学習用説明変数データと称する。
At times t3, t4, and t5 when the magnetic flux density in the instantaneous value of the time waveform 201 of the magnetic flux density is 0, the time waveform 202 of the magnitude of the magnetic flux density becomes non-differentiable. Therefore, as a time differential value of the magnitude of the magnetic flux density at such times t3, t3, and t5, for example, the magnitude of the magnetic flux density at two times temporally preceding and following the times t3, t4, and t5. The arithmetic mean value of the time differential values of can be adopted. Further, at times t3, t4, and t5, there may be no time differential value of the magnitude of the magnetic flux density.
In the following explanation, the value at each time within one cycle of the magnitude of magnetic flux density and the time differential value of the magnitude of magnetic flux density within one cycle, which are derived by the teacher data creation unit 104 as described above, will be described. The values at each time are referred to as learning explanatory variable data as necessary.

また、教師データ作成部104は、選択した測定データに対応する、学習用説明変数データおよび学習用目的変数データが相互に関連付けられたデータを教師データとして作成する。教師データには、選択した測定データに含まれる電磁鋼板の種類を示す情報も含まれる。 Further, the teacher data creation unit 104 creates data in which explanatory variable data for learning and objective variable data for learning are correlated with each other, which corresponds to the selected measurement data, as teacher data. The teacher data also includes information indicating the type of electrical steel sheet included in the selected measurement data.

選択した測定データに対応する学習用説明変数データとは、当該測定データに含まれる磁束密度の一周期内の各時刻における測定値に基づいて教師データ作成部104により導出された、磁束密度の大きさの一周期内の各時刻における値および磁束密度の大きさの時間微分値の一周期内の各時刻における値である。
また、選択した測定データに対応する学習用目的変数データとは、当該測定データに含まれる鉄損の測定値に基づいて異常渦電流損導出部103により導出された、電磁鋼板の異常渦電流損の一周期における代表値である。
The learning explanatory variable data corresponding to the selected measurement data is the magnitude of the magnetic flux density derived by the teacher data creation unit 104 based on the measured values at each time within one cycle of the magnetic flux density included in the measurement data. These are the values at each time within one cycle of the magnetic flux density and the values at each time within one cycle of the time differential value of the magnitude of the magnetic flux density.
Further, the learning objective variable data corresponding to the selected measurement data refers to the abnormal eddy current loss of the electrical steel sheet derived by the abnormal eddy current loss deriving unit 103 based on the measured value of iron loss included in the measurement data. This is a representative value in one cycle.

教師データ作成部104は、測定データを1つずつ選択して以上のようにして電磁鋼板の教師データを作成する。これにより、電磁鋼板の種類別の教師データが作成される。 The teacher data creation unit 104 selects the measurement data one by one and creates the teacher data of the electrical steel sheet as described above. As a result, training data for each type of electrical steel sheet is created.

<<学習部105>>
学習部105は、教師データ作成部104により作成された教師データを用いて、ニューラルネットワークのパラメータを学習する。以下に、学習部105の処理の一例を説明する。
<<Learning Department 105>>
The learning unit 105 uses the teacher data created by the teacher data creation unit 104 to learn the parameters of the neural network. An example of processing by the learning unit 105 will be described below.

学習部105は、教師あり学習を行って、ニューラルネットワークのパラメータを学習する。教師データにおいて学習用目的変数データがラベル(正解データ)となる。
図3は、ニューラルネットワーク300の一例を概念的に示す図である。
ニューラルネットワーク300の入力層301には、或る時刻における磁束密度の大きさ|B|と、当該時刻における磁束密度の大きさの時間微分値d|B|/dtとが少なくとも入力される。また、ニューラルネットネットワークの出力層302から、当該時刻における異常渦電流損の値Waeが少なくとも出力される。ニューラルネットワーク300の中間層(隠れ層)の層数は特に限定されない。
The learning unit 105 performs supervised learning to learn the parameters of the neural network. In the teacher data, the learning objective variable data becomes the label (correct data).
FIG. 3 is a diagram conceptually showing an example of the neural network 300.
The input layer 301 of the neural network 300 receives at least the magnitude of the magnetic flux density |B| at a certain time and the time differential value d|B|/dt of the magnitude of the magnetic flux density at that time. Furthermore, the output layer 302 of the neural network outputs at least the value W ae of the abnormal eddy current loss at the relevant time. The number of intermediate layers (hidden layers) of the neural network 300 is not particularly limited.

学習部105は、教師データに含まれる学習用説明変数データを用いて、或る時刻における磁束密度の大きさ|B|と、当該時刻における磁束密度の大きさの時間微分値d|B|/dtとを入力層301に入力した場合に出力層302から出力される当該時刻における異常渦電流損の値Waeを取得することを、一周期内の各時刻において実行する。このようにして取得される異常渦電流損の値Waeの代表値が、ニューラルネットワーク300から得られる異常渦電流損の一周期における代表値である。代表値としては、例えば、平均値、中央値、実効値を用いることができる。 The learning unit 105 uses learning explanatory variable data included in the teacher data to calculate the magnitude of the magnetic flux density |B| at a certain time and the time differential value d|B|/ of the magnitude of the magnetic flux density at that time. dt is input to the input layer 301, the abnormal eddy current loss value W ae output from the output layer 302 at that time is acquired at each time within one cycle. The representative value of the abnormal eddy current loss value W ae obtained in this way is the representative value in one period of the abnormal eddy current loss obtained from the neural network 300. As the representative value, for example, an average value, a median value, or an effective value can be used.

ニューラルネットワーク300のパラメータを学習する際には、目的関数が用いられる。本実施形態では、目的関数は、教師データに含まれる学習用目的変数データ(異常渦電流損の一周期における代表値)と、ニューラルネットワーク300から得られる異常渦電流損の一周期における代表値との差の評価値を少なくとも導出する関数である。 An objective function is used when learning the parameters of the neural network 300. In this embodiment, the objective function is based on the learning objective variable data (representative value in one cycle of abnormal eddy current loss) included in the teaching data and the representative value in one cycle of abnormal eddy current loss obtained from the neural network 300. is a function that derives at least the evaluation value of the difference between

学習部105は、目的関数の値が最適値になるように、ニューラルネットワーク300のパラメータを学習する。例えば、学習部105は、教師データに含まれる学習用目的変数データ(異常渦電流損の一周期における代表値)と、ニューラルネットワーク300から得られる異常渦電流損の一周期における代表値と、の平均二乗誤差が最小になるときの、ニューラルネットワークのパラメータを、最適なパラメータとして導出する。この場合、前述した評価値は、平均二乗誤差である。 The learning unit 105 learns the parameters of the neural network 300 so that the value of the objective function becomes the optimal value. For example, the learning unit 105 combines the learning objective variable data (representative value in one cycle of abnormal eddy current loss) included in the teacher data and the representative value in one cycle of the abnormal eddy current loss obtained from the neural network 300. The neural network parameters that minimize the mean squared error are derived as the optimal parameters. In this case, the aforementioned evaluation value is the mean square error.

この場合、教師データに含まれる学習用目的変数データ(異常渦電流損の一周期における代表値)と、ニューラルネットワーク300から得られる異常渦電流損の一周期における代表値と、の平均二乗誤差を導出する関数を目的関数とすることができる。目的関数Jの具体例は、以下の(7)式である。 In this case, the mean square error between the learning objective variable data (representative value in one cycle of abnormal eddy current loss) included in the training data and the representative value in one cycle of abnormal eddy current loss obtained from the neural network 300 is calculated. The derived function can be used as an objective function. A specific example of the objective function J is the following equation (7).

Figure 0007428893000004
Figure 0007428893000004

ここで、nは、教師データを特定する変数である。Nは、教師データの数である。tnは、変数nで特定される教師データに含まれる学習用目的変数データ(異常渦電流損の一周期における代表値)である。ynは、変数nで特定される教師データに含まれる学習用説明変数データに基づいてニューラルネットワーク300から得られる異常渦電流損の一周期における代表値である。即ち、ynは、変数nで特定される教師データに含まれる、磁束密度の大きさの一周期内の各時刻における値と、磁束密度の大きさの時間微分値の一周期内の各時刻における値を、ニューラルネットワーク300の入力層301に入力することにより出力層302から出力される一周期内の各時刻における異常渦電流損の値Waeの代表値である。
(7)式においては、目的関数Jの値の最適値は、最小値である。尚、(7)式の右辺に(-1)を掛けたものを目的関数とする場合、目的関数Jの値の最適値は、最大値である。
Here, n is a variable that specifies teacher data. N is the number of training data. t n is learning objective variable data (representative value in one cycle of abnormal eddy current loss) included in the teacher data specified by the variable n. y n is a representative value in one period of the abnormal eddy current loss obtained from the neural network 300 based on the learning explanatory variable data included in the teacher data specified by the variable n. That is, y n is the value at each time within one cycle of the magnitude of magnetic flux density included in the teacher data specified by variable n, and the time differential value of the magnitude of magnetic flux density at each time within one cycle. This is a representative value of the abnormal eddy current loss value W ae at each time within one cycle, which is output from the output layer 302 by inputting the value in the input layer 301 of the neural network 300 .
In equation (7), the optimal value of the objective function J is the minimum value. Note that when the objective function is obtained by multiplying the right side of equation (7) by (-1), the optimal value of the objective function J is the maximum value.

ニューラルネットワークのパラメータには、重み係数が含まれる。ニューラルネットワークの構造に応じて、その他のパラメータを更に学習してもよい。学習の具体的な手法は、例えば、誤差逆伝搬法や、メタヒューリスティックアルゴリズムを用いた手法等、公知の手法で実現することができる。従って、ここでは、学習の手法の詳細な説明を省略する。尚、目的関数の最適値(最小値、最大値)は、ニューラルネットワークのパラメータを学習するためのアルゴリズムにおいて目的関数の値が最適(最小または最大)とされるときの目的関数の値である。例えば、目的関数Jの導出回数が予め定められた回数になったことや、目的関数Jの変化量が所定値以下になったことを条件として、目的関数の値が最適であるとすることができる。 Neural network parameters include weighting coefficients. Depending on the structure of the neural network, other parameters may also be learned. A specific learning method can be realized by a known method such as an error backpropagation method or a method using a metaheuristic algorithm. Therefore, a detailed explanation of the learning method will be omitted here. Note that the optimal value (minimum value, maximum value) of the objective function is the value of the objective function when the value of the objective function is optimal (minimum or maximum) in an algorithm for learning parameters of a neural network. For example, it is possible to determine that the value of the objective function is optimal on the condition that the number of times the objective function J has been derived has reached a predetermined number, or that the amount of change in the objective function J has become less than or equal to a predetermined value. can.

学習部105は、以上のようにしてニューラルネットワーク300のパラメータを学習することを、電磁鋼板の種類を示す情報として同一の情報が含まれている教師データを用いて実行する。また、電磁鋼板の種類を示す情報として相互に異なる種類を示す情報が含まれる教師データがある場合、学習部105は、それぞれの種類について、ニューラルネットワーク300のパラメータを学習する。このようにして、ニューラルネットワークのパラメータが、電磁鋼板の種類ごとに導出される。 The learning unit 105 learns the parameters of the neural network 300 as described above using teacher data that includes the same information as information indicating the type of electrical steel sheet. Furthermore, if there is training data that includes information indicating mutually different types as information indicating the types of electromagnetic steel sheets, the learning unit 105 learns the parameters of the neural network 300 for each type. In this way, neural network parameters are derived for each type of electrical steel sheet.

以上のように本実施形態では、取得部101、第1の電磁場解析部102、異常渦電流損導出部103、教師データ作成部104、および学習部105を用いることにより、ニューラルネットワークのパラメータを学習する機能が実現される。
<<記憶部106>>
記憶部106は、学習部105により学習されたニューラルネットワーク300のパラメータを示す情報を、電磁鋼板の種類ごとに記憶する。
As described above, in this embodiment, the parameters of the neural network are learned by using the acquisition unit 101, the first electromagnetic field analysis unit 102, the abnormal eddy current loss derivation unit 103, the teacher data creation unit 104, and the learning unit 105. The functions to do this are realized.
<<Storage unit 106>>
The storage unit 106 stores information indicating the parameters of the neural network 300 learned by the learning unit 105 for each type of electrical steel sheet.

<<第2の電磁場解析部111>>
第2の電磁場解析部111は、ニューラルネットワークのパラメータを示す情報が記憶部106に記憶された後に起動する。
第2の電磁場解析部111は、異常渦電流損の推定対象の電磁鋼板を、所定の励磁条件で励磁したときの当該電磁鋼板の磁束密度の大きさおよび磁束密度の大きさの時間微分値を、マックスウェルの方程式に基づいて、複数の微小領域(メッシュ)ごとに導出する。ここでの所定の励磁条件は、電磁鋼板をどのように励磁したときの異常渦電流損を推定するかによって定まるものであり、<<第1の電磁場解析部102>>の項に示した所定の励磁条件と異なる場合がある。例えば、異常渦電流損の推定対象の電磁鋼板の適用先である電気機器(例えばモータ)を駆動する際の励磁条件を所定の励磁条件とすることができる。以下に、第2の電磁場解析部111の処理の具体例を説明する。
<<Second electromagnetic field analysis section 111>>
The second electromagnetic field analysis unit 111 is activated after information indicating the parameters of the neural network is stored in the storage unit 106.
The second electromagnetic field analysis unit 111 calculates the magnitude of the magnetic flux density and the time differential value of the magnitude of the magnetic flux density of the electrical steel sheet whose abnormal eddy current loss is to be estimated when the electrical steel sheet is excited under predetermined excitation conditions. , derived for each of multiple micro regions (meshes) based on Maxwell's equations. The predetermined excitation conditions here are determined by how to estimate the abnormal eddy current loss when the electromagnetic steel sheet is excited, and are based on the predetermined excitation conditions shown in the section of <<First electromagnetic field analysis section 102>>. The excitation conditions may be different from the excitation conditions. For example, the excitation conditions when driving an electrical device (for example, a motor) to which the electromagnetic steel sheet whose abnormal eddy current loss is to be estimated can be set to be the predetermined excitation conditions. A specific example of the processing of the second electromagnetic field analysis section 111 will be described below.

第2の電磁場解析部111は、異常渦電流損の推定対象の電磁鋼板の励磁条件を示す情報を取得する。励磁条件を示す情報の取得の形態として、例えば、外部装置からの送信、可搬型記憶媒体に記憶された測定データの読み出し、および異常渦電流損推定装置100のユーザインターフェースに対する測定データの値の入力操作による入力のうち、少なくとも1つを採用することができる。 The second electromagnetic field analysis unit 111 acquires information indicating the excitation conditions of the electromagnetic steel sheet whose abnormal eddy current loss is to be estimated. Forms of acquiring information indicating excitation conditions include, for example, transmission from an external device, reading of measurement data stored in a portable storage medium, and input of measurement data values to the user interface of the abnormal eddy current loss estimating device 100. At least one of the inputs by operation can be adopted.

第2の電磁場解析部111は、取得した励磁条件を含む電磁場解析条件に従って、マックスウェルの方程式に基づき、有限要素法を用いて、各微小領域における磁束密度ベクトルBと渦電流ベクトルJeを一周期における各時刻において導出する。各微小領域における磁束密度ベクトルBと渦電流ベクトルJeを計算する手法は、第1の電磁場解析部102と同様の手法で実現することができる。
第2の電磁場解析部111は、或る微小領域における磁束密度ベクトルBに基づいて、一周期内の各時刻において、当該微小領域についての、磁束密度の大きさ|B|と、磁束密度の大きさの時間微分値d|B|/dtとを導出する。第2の電磁場解析部111は、このような一周期内の各時刻における、磁束密度の大きさ|B|と、磁束密度の大きさの時間微分値d|B|/dtとの導出を、全ての微小領域のそれぞれについて実行する。磁束密度の一周期内の各時刻における値(瞬時値)から、磁束密度の大きさ|B|と、磁束密度の大きさの時間微分値d|B|/dtとを導出する手法は、<<教師データ作成部104>>の項で説明した手法で実現することができる。従って、ここでは、磁束密度の大きさ|B|と、磁束密度の大きさの時間微分値d|B|/dtとを導出する手法の詳細な説明を省略する。
The second electromagnetic field analysis unit 111 unifies the magnetic flux density vector B and the eddy current vector J e in each minute region using the finite element method based on Maxwell's equations according to the electromagnetic field analysis conditions including the acquired excitation conditions. Derived at each time in the cycle. The method of calculating the magnetic flux density vector B and the eddy current vector J e in each minute region can be realized by the same method as the first electromagnetic field analysis section 102 .
The second electromagnetic field analysis unit 111 calculates the magnitude of the magnetic flux density |B| and the magnitude of the magnetic flux density for the microregion at each time in one cycle based on the magnetic flux density vector B in a certain microregion. The time differential value d|B|/dt is derived. The second electromagnetic field analysis unit 111 derives the magnitude of the magnetic flux density |B| and the time differential value of the magnitude of the magnetic flux density d|B|/dt at each time within one cycle, as follows: Execute for each of all microregions. The method of deriving the magnitude of magnetic flux density |B| and the time differential value of the magnitude of magnetic flux density d|B|/dt from the value (instantaneous value) at each time within one cycle of magnetic flux density is as follows. This can be realized by the method described in the section <Teacher data creation unit 104>>. Therefore, a detailed explanation of the method for deriving the magnitude of the magnetic flux density |B| and the time differential value d|B|/dt of the magnitude of the magnetic flux density will be omitted here.

<<推定部112>>
推定部112は、第2の電磁場解析部111により導出された或る時刻における磁束密度の大きさ|B|および磁束密度の大きさの時間微分値d|B|/dtを、ニューラルネットワーク300の入力層301に入力することにより出力層302から出力される値に基づいて、当該時刻における電磁鋼板の異常渦電流損を推定する。以下に、推定部112の処理の具体例を説明する。
<<Estimation unit 112>>
The estimation unit 112 calculates the magnitude of the magnetic flux density |B| at a certain time derived by the second electromagnetic field analysis unit 111 and the time differential value of the magnitude of the magnetic flux density d|B|/dt using the neural network 300. Based on the value input to the input layer 301 and output from the output layer 302, the abnormal eddy current loss of the electrical steel sheet at the relevant time is estimated. A specific example of the processing of the estimation unit 112 will be described below.

推定部112は、異常渦電流損の推定対象の電磁鋼板の種類を示す情報を取得する。電磁鋼板の種類を示す情報の取得の形態として、例えば、外部装置からの送信、可搬型記憶媒体に記憶された測定データの読み出し、および異常渦電流損推定装置100のユーザインターフェースに対する測定データの値の入力操作による入力のうち、少なくとも1つを採用することができる。以下の説明では、異常渦電流損の推定対象の電磁鋼板を、必要に応じて、推定対象の電磁鋼板と略称する。 The estimation unit 112 acquires information indicating the type of electrical steel sheet whose abnormal eddy current loss is to be estimated. Examples of methods of obtaining information indicating the type of electrical steel sheet include transmission from an external device, reading of measurement data stored in a portable storage medium, and value of measurement data to the user interface of the abnormal eddy current loss estimating device 100. At least one of the input operations can be adopted. In the following description, the electrical steel sheet whose abnormal eddy current loss is to be estimated will be abbreviated as the electrical steel sheet to be estimated, if necessary.

推定部112は、推定対象の電磁鋼板の種類に対応するニューラルネットワーク300のパラメータを示す情報を記憶部106から読み出し、ニューラルネットワーク300のパラメータを、ニューラルネットワークのアルゴリズムに設定する。これにより、パラメータが学習されたニューラルネットワーク300が設定される。 The estimation unit 112 reads information indicating the parameters of the neural network 300 corresponding to the type of electrical steel sheet to be estimated from the storage unit 106, and sets the parameters of the neural network 300 to the algorithm of the neural network. As a result, the neural network 300 whose parameters have been learned is set.

推定部112は、第2の電磁場解析部111により、一周期内の或る時刻において或る微小領域に対して導出された、磁束密度の大きさ|B|および磁束密度の大きさの時間微分値d|B|/dtを、パラメータが学習されたニューラルネットワーク300の入力層301に入力することにより出力層302から出力される異常渦電流損の値Waeを、当該時刻の当該微小領域における異常渦電流損の値Waeとして導出する。 The estimation unit 112 calculates the magnitude of the magnetic flux density |B| and the time differential of the magnitude of the magnetic flux density, which is derived for a certain minute region at a certain time within one cycle by the second electromagnetic field analysis unit 111. By inputting the value d|B|/dt into the input layer 301 of the neural network 300 where the parameters have been learned, the value W ae of the abnormal eddy current loss output from the output layer 302 is calculated in the micro region at the time. It is derived as the value W ae of abnormal eddy current loss.

推定部112は、一周期内の各時刻において、以上のような1つの微小領域における異常渦電流損の値Waeの導出を、全ての微小領域のそれぞれについて実行する。これにより、一周期内の各時刻において、推定対象の電磁鋼板の全ての微小領域における異常渦電流損の値Waeが導出される。 The estimating unit 112 calculates the abnormal eddy current loss value W ae in one microregion as described above for each microregion at each time within one cycle. Thereby, at each time within one period, the value W ae of the abnormal eddy current loss in all the micro regions of the electrical steel sheet to be estimated is derived.

推定部112は、同一の時刻における、全ての微小領域における異常渦電流損の値Waeの総和を、推定対象の電磁鋼板の当該時刻における異常渦電流損として導出する。推定部112は、このような推定対象の電磁鋼板の1つの時刻における異常渦電流損の導出を、一周期内の各時刻のそれぞれについて実行する。これにより、一周期内の各時刻において、推定対象の電磁鋼板の異常渦電流損が導出される。 The estimating unit 112 derives the sum of abnormal eddy current loss values W ae in all minute regions at the same time as the abnormal eddy current loss of the electrical steel sheet to be estimated at the time. The estimating unit 112 calculates the abnormal eddy current loss of the electrical steel sheet to be estimated at one time for each time within one cycle. As a result, the abnormal eddy current loss of the electrical steel sheet to be estimated is derived at each time within one cycle.

そして、推定部112は、推定対象の電磁鋼板の異常渦電流損の一周期における代表値を導出する。推定対象の電磁鋼板の異常渦電流損の一周期における代表値としては、例えば、推定対象の電磁鋼板の異常渦電流損の一周期内の各時刻における値の、平均値、中央値、または実効値を用いることができる。 Then, the estimating unit 112 derives a representative value in one cycle of the abnormal eddy current loss of the electromagnetic steel sheet to be estimated. The representative value of the abnormal eddy current loss of the electrical steel sheet to be estimated in one cycle is, for example, the average value, median value, or effective value of the values at each time within one cycle of the abnormal eddy current loss of the electrical steel sheet to be estimated. value can be used.

以上のように本実施形態では、第2の電磁場解析部111および推定部112を用いることにより、異常渦電流損を推定する機能が実現される。 As described above, in this embodiment, the function of estimating abnormal eddy current loss is realized by using the second electromagnetic field analysis section 111 and the estimation section 112.

<<出力部113>>
出力部113は、推定部112により導出された、推定対象の電磁鋼板の異常渦電流損を示す情報を出力する。本実施形態では、推定部112により、一周期内の各時刻における、推定対象の電磁鋼板の異常渦電流損の値と、その代表値とが導出される。従って、出力部113は、これらの情報を出力する。情報の出力の形態としては、例えば、コンピュータディスプレイへの表示、外部装置への送信、および異常渦電流損推定装置100の内部または外部の記憶媒体への記憶のうち、少なくとも1つを採用することができる。
<<Output section 113>>
The output unit 113 outputs information indicating the abnormal eddy current loss of the electrical steel sheet to be estimated, which has been derived by the estimation unit 112. In this embodiment, the estimating unit 112 derives the abnormal eddy current loss value of the electrical steel sheet to be estimated and its representative value at each time within one period. Therefore, the output unit 113 outputs this information. As the form of outputting the information, for example, at least one of displaying it on a computer display, transmitting it to an external device, and storing it in a storage medium inside or outside the abnormal eddy current loss estimating device 100 is adopted. Can be done.

<動作フローチャート>
次に、図4のフローチャートを参照しながら、異常渦電流損推定装置100によるニューラルネットワーク300のパラメータの学習方法の一例を説明する。図4のフローチャートでは、1つの種類の電磁鋼板に対するニューラルネットワーク300のパラメータを学習する場合を例示する。複数の種類の電磁鋼板に対するニューラルネットワーク300のパラメータを学習する場合には、図4のフローチャートによる処理を、電磁鋼板の種類ごとに実行すればよい。
<Operation flowchart>
Next, an example of a method for learning parameters of the neural network 300 by the abnormal eddy current loss estimating device 100 will be described with reference to the flowchart in FIG. 4. The flowchart in FIG. 4 exemplifies a case where parameters of the neural network 300 for one type of electrical steel sheet are learned. When learning the parameters of the neural network 300 for a plurality of types of electromagnetic steel sheets, the process according to the flowchart of FIG. 4 may be executed for each type of electromagnetic steel sheet.

まず、ステップS401において、取得部101は、測定データを1つ取得する。測定データには、同一の種類の電磁鋼板を同一の励磁条件で励磁したときの、磁束密度の一周期における測定値と鉄損の測定値との組が含まれる。また、測定データには、当該測定データを得るために用いた励磁条件を示す情報と、当該測定データにおける測定対象の電磁鋼板の種類を示す情報も含まれる。 First, in step S401, the acquisition unit 101 acquires one piece of measurement data. The measurement data includes a set of a measured value of magnetic flux density in one cycle and a measured value of iron loss when the same type of electromagnetic steel sheet is excited under the same excitation conditions. The measurement data also includes information indicating the excitation conditions used to obtain the measurement data and information indicating the type of electrical steel sheet to be measured in the measurement data.

次に、ステップS402において、第1の電磁場解析部102は、ステップS401で取得された測定データに含まれる励磁条件を含む電磁場解析条件に従って、マックスウェルの方程式に基づき、有限要素法を用いて、各微小領域における磁束密度ベクトルBと渦電流ベクトルJeを一周期における各時刻において計算する。第1の電磁場解析部102は、一周期内の各時刻における各微小領域の磁束密度ベクトルBと渦電流ベクトルJeに基づいて、電磁鋼板のヒステリシス損および古典渦電流損の一周期における代表値を導出する。 Next, in step S402, the first electromagnetic field analysis unit 102 uses the finite element method based on Maxwell's equations according to the electromagnetic field analysis conditions including the excitation conditions included in the measurement data acquired in step S401. The magnetic flux density vector B and eddy current vector J e in each minute region are calculated at each time in one cycle. The first electromagnetic field analysis unit 102 calculates representative values of hysteresis loss and classical eddy current loss in one cycle of the electrical steel sheet based on the magnetic flux density vector B and eddy current vector J e of each minute region at each time in one cycle. Derive.

次に、ステップS403において、異常渦電流損導出部103は、ステップS401で取得された測定データに含まれる鉄損の測定値から、ステップS402で導出されたヒステリシス損および古典渦電流損の一周期における代表値の和を減算した値を、電磁鋼板の異常渦電流損の一周期における代表値(学習用目的変数データ)として導出する。 Next, in step S403, the abnormal eddy current loss deriving unit 103 calculates one period of the hysteresis loss and classical eddy current loss derived in step S402 from the measured value of iron loss included in the measurement data acquired in step S401. The value obtained by subtracting the sum of the representative values in is derived as the representative value (learning objective variable data) in one cycle of the abnormal eddy current loss of the electrical steel sheet.

次に、ステップS404において、教師データ作成部104は、ステップS401で取得された測定データに含まれる磁束密度の一周期内の各時刻における測定値に基づいて、磁束密度の大きさの一周期内の各時刻における値と、磁束密度の大きさの時間微分値の一周期内の各時刻における値とを、学習用説明変数データとして導出する。教師データ作成部104は、このようにして導出した学習用説明変数データと、ステップS403で導出された学習用目的変数データと、ステップS401で取得された測定データに含まれる電磁鋼板の種類を示す情報とを含む教師データを1つ作成する。 Next, in step S404, the teacher data creation unit 104 calculates the magnitude of the magnetic flux density within one cycle based on the measured values at each time within one cycle of the magnetic flux density included in the measurement data acquired in step S401. The value at each time and the value at each time within one cycle of the time differential value of the magnitude of the magnetic flux density are derived as explanatory variable data for learning. The teacher data creation unit 104 indicates the type of electromagnetic steel sheet included in the learning explanatory variable data derived in this way, the learning objective variable data derived in step S403, and the measurement data acquired in step S401. Create one piece of teacher data including information.

次に、ステップS405において、異常渦電流損推定装置100は、ニューラルネットワーク300のパラメータの学習に必要な数の教師データが作成されたか否かを判定する。この数は、例えば、異常渦電流損推定装置100のユーザインターフェースの操作に基づいてオペレータにより予め設定される。様々な励磁条件における教師データを用いることによりニューラルネットワークの推定精度は向上する。従って、教師データの数は多い方が好ましい。一方で、教師データを作成すればするほど作業の負担が増大する。このような観点から、教師データの数を定めることができる。 Next, in step S405, the abnormal eddy current loss estimating device 100 determines whether the necessary number of training data for learning the parameters of the neural network 300 has been created. This number is preset by an operator based on the operation of the user interface of the abnormal eddy current loss estimating device 100, for example. The estimation accuracy of the neural network is improved by using training data under various excitation conditions. Therefore, it is preferable to have a large number of training data. On the other hand, the more training data is created, the more the work burden increases. From this perspective, the number of training data can be determined.

例えば、測定データを用いずに任意の教師データを用意する。当該教師データを用いてニューラルネットワーク300のパラメータを学習することを、教師データの数を変えて実行する。学習後のニューラルネットワークから出力される値が略一定となる学習データの数の最小値を調査する。この学習データの数の最小値を、ニューラルネットワーク300のパラメータの学習に必要な数とすることができる。 For example, arbitrary teacher data may be prepared without using measurement data. Learning the parameters of the neural network 300 using the teacher data is performed by changing the number of teacher data. Investigate the minimum number of learning data at which the value output from the neural network after learning is approximately constant. The minimum value of the number of learning data can be set to the number necessary for learning the parameters of the neural network 300.

ステップS405の判定の結果、ニューラルネットワーク300のパラメータの学習に必要な数の教師データが作成されていない場合、処理は、ステップS401に戻る。ステップS401において、取得部101は、ステップS401で既に取得している測定データとは異なる測定データを1つ取得する。ステップS401で既に取得している測定データとは異なる測定データは、ステップS401で既に取得している測定データに含まれる励磁条件とは異なる励磁条件を示す情報が含まれる測定データである。そして、新たに取得した測定データに基づいて、ステップS402~S404の処理が実行される。 As a result of the determination in step S405, if the number of teacher data necessary for learning the parameters of the neural network 300 has not been created, the process returns to step S401. In step S401, the acquisition unit 101 acquires one measurement data different from the measurement data already acquired in step S401. The measurement data different from the measurement data already acquired in step S401 is measurement data that includes information indicating an excitation condition different from the excitation condition included in the measurement data already acquired in step S401. Then, the processes of steps S402 to S404 are executed based on the newly acquired measurement data.

以上のようにしてステップS405において、ニューラルネットワーク300のパラメータの学習に必要な数の教師データが作成されたと判定されるまで、ステップS401~S405の処理が繰り返し実行される。そして、ステップS405において、ニューラルネットワーク300のパラメータの学習に必要な数の教師データが作成されたと判定されると、処理は、ステップS406に進む。 As described above, the processes of steps S401 to S405 are repeatedly executed until it is determined in step S405 that the necessary number of teacher data for learning the parameters of the neural network 300 has been created. If it is determined in step S405 that the necessary number of teacher data for learning the parameters of the neural network 300 has been created, the process proceeds to step S406.

処理がステップS406に進むと、学習部105は、ステップS404で作成された教師データを用いて、ニューラルネットワーク300のパラメータを学習する。
最後に、ステップS407において、記憶部106は、ステップS406で学習されたニューラルネットワーク300のパラメータを示す情報を記憶する。このとき、記憶部106は、ニューラルネットワーク300のパラメータを示す情報と、当該ニューラルネットワーク300による推定対象の電磁鋼板の種類を示す情報とを相互に関連付けて記憶する。ステップS407の処理が終了すると、図4のフローチャートによる処理が終了する。
When the process proceeds to step S406, the learning unit 105 learns the parameters of the neural network 300 using the teacher data created in step S404.
Finally, in step S407, the storage unit 106 stores information indicating the parameters of the neural network 300 learned in step S406. At this time, the storage unit 106 stores information indicating the parameters of the neural network 300 and information indicating the type of electrical steel sheet to be estimated by the neural network 300 in association with each other. When the process in step S407 ends, the process according to the flowchart in FIG. 4 ends.

次に、図5のフローチャートを参照しながら、異常渦電流損推定装置100による異常渦電流損推定方法の一例を説明する。図5のフローチャートは、図4のフローチャートによる処理が終了した後に開始する。
まず、ステップS501において、第2の電磁場解析部111は、推定対象の電磁鋼板の励磁条件を示す情報を取得し、取得した励磁条件を含む電磁場解析条件に従って、マックスウェルの方程式に基づき、有限要素法を用いて、各微小領域における磁束密度ベクトルBと渦電流ベクトルJeを一周期における各時刻において計算する。第2の電磁場解析部111は、各微小領域における磁束密度ベクトルBに基づいて、一周期内の各時刻において、磁束密度の大きさ|B|と、磁束密度の大きさの時間微分値d|B|/dtとを導出することを、各微小領域のそれぞれについて実行する。
Next, an example of a method for estimating abnormal eddy current loss by the abnormal eddy current loss estimating device 100 will be described with reference to the flowchart in FIG. 5 . The flowchart in FIG. 5 starts after the process according to the flowchart in FIG. 4 is completed.
First, in step S501, the second electromagnetic field analysis unit 111 acquires information indicating the excitation conditions of the electromagnetic steel sheet to be estimated, and in accordance with the electromagnetic field analysis conditions including the acquired excitation conditions, calculates the finite element based on Maxwell's equation. Using the method, the magnetic flux density vector B and eddy current vector J e in each minute region are calculated at each time in one cycle. The second electromagnetic field analysis unit 111 calculates the magnitude of the magnetic flux density |B| and the time differential value d| of the magnitude of the magnetic flux density at each time within one cycle based on the magnetic flux density vector B in each minute region. B|/dt is derived for each minute region.

次に、ステップS502において、推定部112は、ステップS501で導出された、磁束密度の大きさ|B|および磁束密度の大きさの時間微分値d|B|/dtを、図4のフローチャートによりパラメータが学習されたニューラルネットワーク300の入力層301に入力することにより出力層302から出力された異常渦電流損の値Waeに基づいて、推定対象の電磁鋼板の一周期内の各時刻における異常渦電流損と、推定対象の電磁鋼板の異常渦電流損の一周期における代表値とを導出する。 Next, in step S502, the estimation unit 112 calculates the magnitude of the magnetic flux density |B| and the time differential value of the magnitude of the magnetic flux density d|B|/dt derived in step S501 according to the flowchart of FIG. Based on the abnormal eddy current loss value W ae output from the output layer 302 by inputting it to the input layer 301 of the neural network 300 where the parameters have been learned, the abnormality at each time within one cycle of the electrical steel sheet to be estimated is determined. The eddy current loss and the representative value in one cycle of the abnormal eddy current loss of the electrical steel sheet to be estimated are derived.

ステップS502においてニューラルネットワーク300の入力層301に入力される、磁束密度の大きさ|B|および磁束密度の大きさの時間微分値d|B|/dtは、或る1つ時刻および或る1つの微小領域におけるものである。従って、ニューラルネットワーク300の出力層302から出力される異常渦電流損の値Waeは、当該時刻および当該微小領域における値である。従って、ニューラルネットワーク300によって、1つの時刻ごと、および、1つの微小領域ごとに、異常渦電流損の値Waeが導出される。 The magnitude of the magnetic flux density |B| and the time differential value of the magnitude of the magnetic flux density d|B|/dt, which are input to the input layer 301 of the neural network 300 in step S502, are determined at a certain time and a certain time. This is in two microscopic areas. Therefore, the abnormal eddy current loss value W ae output from the output layer 302 of the neural network 300 is the value at the time and in the minute area. Therefore, the value W ae of the abnormal eddy current loss is derived by the neural network 300 for each time and for each minute region.

推定部112は、同一の時刻における異常渦電流損の値Waeの総和を、当該時刻における異常渦電流損の値として導出する。推定部112は、このような導出を、一周期内の各時刻において実行することにより、一周期内の各時刻において、推定対象の電磁鋼板の異常渦電流損を導出する。推定部112は、推定対象の電磁鋼板の異常渦電流損の一周期内の各時刻における値の代表値を、推定対象の電磁鋼板の異常渦電流損の一周期における代表値として導出する。 The estimation unit 112 derives the sum of the abnormal eddy current loss values W ae at the same time as the abnormal eddy current loss value at the time. The estimation unit 112 derives the abnormal eddy current loss of the electrical steel sheet to be estimated at each time within one cycle by performing such derivation at each time within one cycle. The estimation unit 112 derives a representative value of the values at each time within one cycle of the abnormal eddy current loss of the electrical steel sheet to be estimated as a representative value of the abnormal eddy current loss of the electrical steel sheet to be estimated in one cycle.

最後に、ステップS503において、出力部113は、ステップS502で導出された、推定対象の電磁鋼板の異常渦電流損を示す情報を出力する。ステップS503の処理が終了すると、図5のフローチャートによる処理が終了する。 Finally, in step S503, the output unit 113 outputs information indicating the abnormal eddy current loss of the electrical steel sheet to be estimated, which was derived in step S502. When the process of step S503 ends, the process according to the flowchart of FIG. 5 ends.

<実施例>
次に、実施例を説明する。
発明例として、単板の電磁鋼板と同じ種類の電磁鋼板に対する測定データを用いて、本実施形態で説明した手法によりニューラルネットワーク300のパラメータを学習した。当該単板の電磁鋼板について、パラメータが学習されたニューラルネットワーク300を用いて、本実施形態で説明した手法により異常渦電流損の一周期における代表値を推定した。
<Example>
Next, an example will be described.
As an example of the invention, the parameters of the neural network 300 were learned by the method described in the present embodiment using measurement data for an electrical steel sheet of the same type as a single electromagnetic steel sheet. Regarding the single electromagnetic steel sheet, a representative value of the abnormal eddy current loss in one cycle was estimated by the method described in this embodiment using the neural network 300 whose parameters were learned.

比較例として、発明例の単板の電磁鋼板と同一の単板の電磁鋼板について、特許文献1に記載の手法を用いて、異常渦電流損を推定した。
比較例においては、具体的に以下のようにして異常渦電流損を導出した。
当該単板の電磁鋼板について、特許文献1に記載の手法により、異常渦電流損係数を導出し、異常渦電流損係数を用いて導電率を補正する。そして、当該単板の電磁鋼板について、有限要素法を用いて、各微小領域(メッシュ)における磁束密度ベクトルBと渦電流ベクトルJeを計算する。そして、(6)式のσに補正後の導電率を与えると共に(6)式の右辺全体に異常渦電流損係数を乗算して全渦電流損を導出することと、(6)式のσに補正前の導電率を与えて古典渦電流損を導出することとを実行する。そして、全渦電流損から古典渦電流損を減算した値を異常渦電流損とする。この異常渦電流損は、一周期における代表値である。
As a comparative example, abnormal eddy current loss was estimated using the method described in Patent Document 1 for a single electromagnetic steel sheet that is the same as the single electromagnetic steel sheet of the invention example.
In the comparative example, the abnormal eddy current loss was specifically derived as follows.
Regarding the single electromagnetic steel sheet, an abnormal eddy current loss coefficient is derived by the method described in Patent Document 1, and the conductivity is corrected using the abnormal eddy current loss coefficient. Then, for the single electromagnetic steel sheet, the magnetic flux density vector B and eddy current vector J e in each minute region (mesh) are calculated using the finite element method. Then, the corrected conductivity is given to σ in equation (6), and the entire right side of equation (6) is multiplied by the abnormal eddy current loss coefficient to derive the total eddy current loss. The conductivity before correction is given to derive the classical eddy current loss. Then, the value obtained by subtracting the classical eddy current loss from the total eddy current loss is defined as the abnormal eddy current loss. This abnormal eddy current loss is a typical value in one cycle.

また、当該単板の電磁鋼板の鉄損の測定値から、当該単板の電磁鋼板のヒステリシス損および古典渦電流損の和を減算した値を、異常渦電流損の測定値として導出した。この異常渦電流損は、一周期における代表値である。単板の電磁鋼板のヒステリシス損、単板の古典渦電流損は、それぞれ、(5)式、(6)式を用いて導出した。 In addition, the value obtained by subtracting the sum of the hysteresis loss and classical eddy current loss of the single electromagnetic steel sheet from the measured value of the iron loss of the single electromagnetic steel sheet was derived as the measured value of abnormal eddy current loss. This abnormal eddy current loss is a typical value in one cycle. The hysteresis loss of a single electromagnetic steel sheet and the classical eddy current loss of a single sheet were derived using equations (5) and (6), respectively.

以上のような異常渦電流損の推定を、励磁条件を異ならせて行った。図6は、発明例および比較例における異常渦電流損の推定値と、異常渦電流損の測定値とを比較した結果を示す。
図6に示す値は、異常渦電流損の推定値の測定値に対する誤差率(=推定値÷測定値-1)を百分率で表したものである。
波形1、2、3は、磁束密度の時間波形を示す。波形1は、基本周波数が50Hzであり、最大値が1.0Tであり、高調波が相対的に少ない時間波形であることを示す。波形2は、基本周波数が50Hzであり、最大値が1.5Tであり、高調波が相対的に多い時間波形であることを示す。波形3は、基本周波数が1kHzであり、最大値が1.5Tであり、高調波が相対的に多い時間波形であることを示す。
The abnormal eddy current loss as described above was estimated under different excitation conditions. FIG. 6 shows the results of comparing estimated values of abnormal eddy current loss and measured values of abnormal eddy current loss in the invention example and comparative example.
The values shown in FIG. 6 are the error rate (=estimated value÷measured value−1) of the estimated value of abnormal eddy current loss with respect to the measured value, expressed as a percentage.
Waveforms 1, 2, and 3 indicate time waveforms of magnetic flux density. Waveform 1 has a fundamental frequency of 50 Hz and a maximum value of 1.0 T, indicating that it is a temporal waveform with relatively few harmonics. Waveform 2 has a fundamental frequency of 50 Hz and a maximum value of 1.5 T, indicating that it is a time waveform with relatively many harmonics. Waveform 3 has a fundamental frequency of 1 kHz and a maximum value of 1.5 T, indicating that it is a time waveform with relatively many harmonics.

図6に示すように、発明例では、比較例に比べて、異常渦電流損を精度よく推定することができることが分かる。特に、発明例では、磁束密度の時間波形の基本周波数が高い場合や、磁束密度の時間波形に高調波成分が重畳される場合でも、異常渦電流損を精度よく推定することができる。
また、比較例では、異常渦電流損として、一周期における代表値しか導出することができない。これに対し、発明例では、異常渦電流損として、一周期内の各時刻における値を導出することができる。
As shown in FIG. 6, it can be seen that in the invention example, abnormal eddy current loss can be estimated with higher accuracy than in the comparative example. In particular, in the example of the invention, abnormal eddy current loss can be accurately estimated even when the fundamental frequency of the time waveform of magnetic flux density is high or when harmonic components are superimposed on the time waveform of magnetic flux density.
Further, in the comparative example, only a representative value in one cycle can be derived as the abnormal eddy current loss. In contrast, in the example of the invention, a value at each time within one period can be derived as the abnormal eddy current loss.

<まとめ>
以上のように本実施形態では、異常渦電流損推定装置100は、或る時刻における磁束密度の大きさ|B|と、当該時刻における磁束密度の大きさの時間微分値d|B|/dtとが入力される入力層301と、当該時刻における異常渦電流損の値Waeが出力される出力層302とを有するニューラルネットワーク300のパラメータを学習する。異常渦電流損推定装置100は、学習後のニューラルネットワーク300を用いて、各時刻における異常渦電流損の値Waeを推定する。従って、磁束密度の時間波形に高調波が重畳される場合や、磁束密度の時間波形が高周波数である場合であっても、異常渦電流損を精度よく導出することができる。また、一周期における各時刻における異常渦電流損を導出することができる。よって、電磁鋼板における異常渦電流損を正確に導出することができる。
<Summary>
As described above, in this embodiment, the abnormal eddy current loss estimating device 100 calculates the magnitude of magnetic flux density at a certain time |B| and the time differential value d|B|/dt of the magnitude of magnetic flux density at that time. The parameters of the neural network 300 are learned, which has an input layer 301 to which the values are input, and an output layer 302 to which the abnormal eddy current loss value W ae at the time is output. The abnormal eddy current loss estimating device 100 uses the learned neural network 300 to estimate the abnormal eddy current loss value W ae at each time. Therefore, even when harmonics are superimposed on the time waveform of magnetic flux density or when the time waveform of magnetic flux density has a high frequency, abnormal eddy current loss can be derived with high accuracy. Furthermore, the abnormal eddy current loss at each time in one cycle can be derived. Therefore, the abnormal eddy current loss in the electrical steel sheet can be accurately derived.

<変形例>
本実施形態では、磁束密度が磁化特性から得られる物理量である場合を例示した。しかしながら、前述したように、磁束密度、磁界強度、透磁率、および微分透磁率の少なくとも1つを、磁化特性から得られる物理量として用いることができる。ただし、磁束密度と磁界強度は、相互に対応する物理量であり、一方を用い、他方を用いなくても、異常渦電流損の推定精度に大きな影響を与えないと考えられる。また、磁束密度と磁界強度との双方を用いると、ニューラルネットワークのパラメータの学習、およびニューラルネットワークを用いた異常渦電流損の推定の際の処理負荷が高くなる虞がある。従って、磁化特性から得られる物理量として磁束密度を採用する場合には磁界強度を採用せず、逆に、磁界強度を採用する場合には磁束密度を採用しないようにするのが好ましい。
<Modified example>
In this embodiment, the case where the magnetic flux density is a physical quantity obtained from magnetization characteristics is exemplified. However, as described above, at least one of magnetic flux density, magnetic field strength, magnetic permeability, and differential magnetic permeability can be used as the physical quantity obtained from the magnetization characteristics. However, the magnetic flux density and the magnetic field strength are physical quantities that correspond to each other, and it is thought that using one and not using the other will not have a large effect on the estimation accuracy of abnormal eddy current loss. Furthermore, if both magnetic flux density and magnetic field strength are used, there is a risk that the processing load will increase when learning the parameters of the neural network and estimating abnormal eddy current loss using the neural network. Therefore, when magnetic flux density is used as a physical quantity obtained from magnetization characteristics, it is preferable not to use magnetic field strength, and conversely, when magnetic field strength is used, it is preferable not to use magnetic flux density.

本実施形態では、ニューラルネットワークのパラメータを学習する機能(取得部101、第1の電磁場解析部102、異常渦電流損導出部103、教師データ作成部104、および学習部105)を実現する学習装置と、ニューラルネットワークのパラメータを記憶する機能(記憶部106)を実現する記憶する記憶装置が、異常渦電流損を推定する機能(第2の電磁場解析部111および推定部112)を実現する異常渦電流損推定装置の内部に含まれる場合を例示した。しかしながら、これらの機能を複数の装置で実現してもよい。 In this embodiment, a learning device that realizes the function of learning neural network parameters (acquisition unit 101, first electromagnetic field analysis unit 102, abnormal eddy current loss derivation unit 103, teacher data creation unit 104, and learning unit 105) and a storage device that realizes the function of storing the neural network parameters (storage unit 106), and a storage device that realizes the function of storing the parameters of the neural network (the storage unit 106), and the abnormal eddy current loss that realizes the function of estimating the abnormal eddy current loss (the second electromagnetic field analysis unit 111 and the estimation unit 112). The case where it is included inside the current loss estimating device is illustrated. However, these functions may be realized by multiple devices.

例えば、ニューラルネットワークのパラメータを学習する機能(取得部101、第1の電磁場解析部102、異常渦電流損導出部103、教師データ作成部104、および学習部105)を実現する学習装置と、ニューラルネットワークのパラメータを記憶する機能(記憶部106)を実現する記憶装置と、異常渦電流損を推定する機能(第2の電磁場解析部111および推定部112)を実現する異常渦電流損推定装置とを有する異常渦電流損推定システムを構築してもよい。また、学習装置の内部に記憶装置が含まれていてもよい。 For example, a learning device that realizes a function of learning the parameters of a neural network (an acquisition unit 101, a first electromagnetic field analysis unit 102, an abnormal eddy current loss derivation unit 103, a teacher data creation unit 104, and a learning unit 105) and a neural A storage device that realizes a function of storing network parameters (storage unit 106), and an abnormal eddy current loss estimation device that realizes a function of estimating abnormal eddy current loss (second electromagnetic field analysis unit 111 and estimation unit 112). An abnormal eddy current loss estimation system may be constructed. Further, a storage device may be included inside the learning device.

本実施形態では、異常渦電流損の推定対象の電磁鋼板を、所定の励磁条件で励磁したときの当該電磁鋼板の磁束密度の大きさおよび磁束密度の大きさの時間微分値を導出する導出部として第2の電磁場解析部111を適用した場合について説明したが、導出部は第2の電磁場解析部111に限られない。即ち、異常渦電流損の推定対象の電磁鋼板を、所定の励磁条件で励磁したときの当該電磁鋼板の磁束密度の大きさおよび磁束密度の大きさの時間微分値を導出できればよく、例えば実験を行い、実験結果からこれらの値を導出することとしてもよい。例えば、電磁鋼板に励磁コイルと検出コイルとを巻き回し、励磁コイルに励磁電流を流したときに検出コイルに発生する誘導起電力に基づいて磁束密度を測定してもよい。また、第2の電磁場解析部111を適用した場合として、電磁鋼板の磁束密度の大きさおよび磁束密度の大きさの時間微分値を、マックスウェルの方程式に基づいて、複数の微小領域(メッシュ)ごとに導出することを示したが、これに限られない。例えば、所定の励磁条件で励磁される電磁鋼板の等価回路を用いてメッシュ分割無しに磁束密度の大きさおよび磁束密度の大きさの時間微分値を導出することとしてもよい。尚、以上のようにするに際し、磁束密度を用いなくてもよいことは、<変形例>の項の冒頭で説明した通りである。実験により磁界強度を検出する際には、例えば、Hコイルを用いてもよい。 In this embodiment, a derivation unit that derives the magnitude of the magnetic flux density and the time differential value of the magnitude of the magnetic flux density of the electrical steel sheet whose abnormal eddy current loss is to be estimated is excited under predetermined excitation conditions. Although the case where the second electromagnetic field analysis section 111 is applied has been described, the derivation section is not limited to the second electromagnetic field analysis section 111. In other words, it is only necessary to derive the magnitude of the magnetic flux density of the electrical steel sheet whose abnormal eddy current loss is to be estimated and the time differential value of the magnitude of the magnetic flux density when the electrical steel sheet is excited under predetermined excitation conditions. These values may be derived from the experimental results. For example, an excitation coil and a detection coil may be wound around an electromagnetic steel plate, and the magnetic flux density may be measured based on the induced electromotive force generated in the detection coil when an excitation current is passed through the excitation coil. In addition, when the second electromagnetic field analysis unit 111 is applied, the magnitude of the magnetic flux density of the electromagnetic steel sheet and the time differential value of the magnitude of the magnetic flux density are divided into multiple micro regions (mesh) based on Maxwell's equation. Although it has been shown that it can be derived for each, it is not limited to this. For example, the magnitude of the magnetic flux density and the time differential value of the magnitude of the magnetic flux density may be derived without mesh division using an equivalent circuit of an electromagnetic steel sheet that is excited under predetermined excitation conditions. As explained at the beginning of the <Modifications> section, it is not necessary to use the magnetic flux density when doing the above. When detecting the magnetic field strength through experiments, for example, an H coil may be used.

尚、以上説明した本発明の実施形態は、コンピュータがプログラムを実行することによって実現することができる。また、前記プログラムを記録したコンピュータ読み取り可能な記録媒体及び前記プログラム等のコンピュータプログラムプロダクトも本発明の実施形態として適用することができる。記録媒体としては、例えば、フレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、CD-ROM、磁気テープ、不揮発性のメモリカード、ROM等を用いることができる。
また、以上説明した本発明の実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
Note that the embodiments of the present invention described above can be realized by a computer executing a program. Furthermore, a computer-readable recording medium on which the program is recorded and a computer program product such as the program can also be applied as an embodiment of the present invention. As the recording medium, for example, a flexible disk, hard disk, optical disk, magneto-optical disk, CD-ROM, magnetic tape, nonvolatile memory card, ROM, etc. can be used.
Furthermore, the embodiments of the present invention described above are merely examples of implementation of the present invention, and the technical scope of the present invention should not be construed as limited by these. It is something. That is, the present invention can be implemented in various forms without departing from its technical idea or its main features.

100:異常渦電流損推定装置、101:取得部、102:第1の電磁場解析部、103:異常渦電流損導出部、104:教師データ作成部、105:学習部、106:記憶部、111:第2の電磁場解析部、112:推定部、113:出力部 100: Abnormal eddy current loss estimation device, 101: Acquisition unit, 102: First electromagnetic field analysis unit, 103: Abnormal eddy current loss derivation unit, 104: Teacher data creation unit, 105: Learning unit, 106: Storage unit, 111 : second electromagnetic field analysis section, 112: estimation section, 113: output section

Claims (9)

軟磁性材料の異常渦電流損を推定するためのニューラルネットワークのパラメータを学習する学習装置であって、
前記軟磁性材料を所定の励磁条件で励磁したときの当該軟磁性材料に対する測定値として、当該軟磁性材料の磁化特性から得られる物理量の一周期内の各時刻における測定値と、当該軟磁性材料の鉄損の測定値とを取得する取得手段と、
前記軟磁性材料を前記所定の励磁条件で励磁したときの前記軟磁性材料のヒステリシス損および古典渦電流損の一周期における代表値を、マックスウェルの方程式に基づいて導出する第1の電磁場解析手段と、
前記取得手段により取得された鉄損の測定値と、前記第1の電磁場解析手段により導出されたヒステリシス損および古典渦電流損の一周期における代表値とに基づいて、前記軟磁性材料における異常渦電流損の一周期における代表値を導出する異常渦電流損導出手段と、
前記異常渦電流損導出手段により導出された前記軟磁性材料における異常渦電流損の一周期における代表値と、前記物理量の大きさの一周期内の各時刻における値と、前記物理量の大きさの時間微分値の一周期内の各時刻における値と、に基づいて、教師データを作成する教師データ作成手段と、
前記教師データ作成手段により作成された前記教師データを用いて、ニューラルネットワークのパラメータを学習する学習手段と、を有し、
前記ニューラルネットワークの入力層には、或る時刻における前記物理量の大きさおよび当該時刻における前記物理量の大きさの時間微分値が少なくとも入力され、
前記ニューラルネットワークの出力層から、当該時刻における異常渦電流損が少なくとも出力され、
前記学習手段は、前記教師データに含まれる前記軟磁性材料の異常渦電流損の一周期における代表値と、前記ニューラルネットワークの前記出力層から出力される一周期内の各時刻における異常渦電流損の代表値との差の評価値を少なくとも求める目的関数を用いて、前記ニューラルネットワークのパラメータを学習することを特徴とする学習装置。
A learning device for learning neural network parameters for estimating abnormal eddy current loss in soft magnetic materials,
Measured values for the soft magnetic material when the soft magnetic material is excited under predetermined excitation conditions include measured values at each time within one period of a physical quantity obtained from the magnetization characteristics of the soft magnetic material, and the soft magnetic material. an acquisition means for acquiring a measured value of iron loss;
A first electromagnetic field analysis means for deriving representative values in one cycle of hysteresis loss and classical eddy current loss of the soft magnetic material when the soft magnetic material is excited under the predetermined excitation conditions based on Maxwell's equations. and,
Based on the measured value of iron loss acquired by the acquisition means and the representative values of hysteresis loss and classical eddy current loss in one cycle derived by the first electromagnetic field analysis means, abnormal eddy in the soft magnetic material is determined. Abnormal eddy current loss derivation means for deriving a representative value of current loss in one cycle;
A typical value of the abnormal eddy current loss in the soft magnetic material derived by the abnormal eddy current loss deriving means in one cycle, a value of the magnitude of the physical quantity at each time within one cycle, and a value of the magnitude of the physical quantity at each time within one cycle. a value at each time within one cycle of the time differential value, and a teacher data creation means for creating teacher data based on the time differential value at each time within one cycle;
a learning means for learning parameters of a neural network using the teacher data created by the teacher data creation means,
At least the magnitude of the physical quantity at a certain time and the time differential value of the magnitude of the physical quantity at the time are input to the input layer of the neural network,
At least the abnormal eddy current loss at the time is output from the output layer of the neural network,
The learning means is configured to calculate a representative value of the abnormal eddy current loss of the soft magnetic material in one period included in the training data and an abnormal eddy current loss at each time in one period output from the output layer of the neural network. A learning device, characterized in that the parameters of the neural network are learned using an objective function that at least obtains an evaluation value of a difference between a representative value and a representative value.
前記物理量は、磁束密度、磁界強度、透磁率、および微分透磁率の少なくとも1つを含むことを特徴とする請求項1に記載の学習装置。 The learning device according to claim 1, wherein the physical quantity includes at least one of magnetic flux density, magnetic field strength, magnetic permeability, and differential magnetic permeability. 請求項1または2に記載の学習装置でパラメータが学習された前記ニューラルネットワークを用いて軟磁性材料の異常渦電流損を推定する異常渦電流損推定装置であって、
前記異常渦電流損の推定対象の軟磁性材料を所定の励磁条件で励磁したときの当該軟磁性材料の前記物理量の大きさおよび当該物理量の大きさの時間微分値を導出する導出手段と、
前記導出手段により導出された或る時刻における前記物理量の大きさおよび当該時刻における前記物理量の大きさの時間微分値を前記ニューラルネットワークの前記入力層に与えることにより前記出力層から出力される値に基づいて、当該時刻における前記軟磁性材料の異常渦電流損を推定する推定手段と、を有することを特徴とする異常渦電流損推定装置。
An abnormal eddy current loss estimating device that estimates abnormal eddy current loss of a soft magnetic material using the neural network whose parameters have been learned by the learning device according to claim 1 or 2,
Derivation means for deriving the magnitude of the physical quantity of the soft magnetic material and the time differential value of the magnitude of the physical quantity when the soft magnetic material whose abnormal eddy current loss is to be estimated is excited under predetermined excitation conditions;
The magnitude of the physical quantity at a certain time derived by the derivation means and the time differential value of the magnitude of the physical quantity at the time are given to the input layer of the neural network to obtain the value output from the output layer. an estimation means for estimating the abnormal eddy current loss of the soft magnetic material at the relevant time based on the above information.
前記導出手段は、前記異常渦電流損の推定対象の軟磁性材料を所定の励磁条件で励磁したときの当該軟磁性材料の前記物理量の大きさおよび当該物理量の大きさの時間微分値を、マックスウェルの方程式に基づいて、複数の微小領域ごとに導出する第2の電磁場解析手段であることを特徴とする請求項3に記載の異常渦電流損推定装置。 The derivation means maximizes the magnitude of the physical quantity of the soft magnetic material and the time differential value of the magnitude of the physical quantity when the soft magnetic material whose abnormal eddy current loss is to be estimated is excited under predetermined excitation conditions. The abnormal eddy current loss estimating device according to claim 3, characterized in that the second electromagnetic field analysis means derives the electromagnetic field for each of a plurality of micro regions based on Well's equation. 前記推定手段は、前記導出手段により導出された或る時刻における前記物理量の大きさおよび当該時刻における前記物理量の大きさの時間微分値を前記ニューラルネットワークの前記入力層に与えることにより前記出力層から出力される当該時刻における異常渦電流損を取得することを、一周期内の各時刻において実行し、当該取得した各時刻における異常渦電流損に基づいて、前記軟磁性材料の異常渦電流損の一周期における代表値を推定することを特徴とする請求項3または4に記載の異常渦電流損推定装置。 The estimating means provides the magnitude of the physical quantity at a certain time derived by the deriving means and the time differential value of the magnitude of the physical quantity at the time to the input layer of the neural network, thereby estimating the magnitude of the physical quantity from the output layer. Obtaining the output abnormal eddy current loss at the relevant time is executed at each time within one cycle, and based on the abnormal eddy current loss at each acquired time, the abnormal eddy current loss of the soft magnetic material is calculated. The abnormal eddy current loss estimating device according to claim 3 or 4, wherein the abnormal eddy current loss estimating device estimates a representative value in one cycle. 軟磁性材料の異常渦電流損を推定するためのニューラルネットワークのパラメータを学習する学習方法であって、
前記軟磁性材料を所定の励磁条件で励磁したときの当該軟磁性材料に対する測定値として、当該軟磁性材料の磁化特性から得られる物理量の一周期内の各時刻における測定値と、当該軟磁性材料の鉄損の測定値とを取得する取得工程と、
前記軟磁性材料を前記所定の励磁条件で励磁したときの前記軟磁性材料のヒステリシス損および古典渦電流損の一周期における代表値を、マックスウェルの方程式に基づいて導出する第1の電磁場解析工程と、
前記取得工程により取得された鉄損の測定値と、前記第1の電磁場解析工程により導出されたヒステリシス損および古典渦電流損の一周期における代表値とに基づいて、前記軟磁性材料における異常渦電流損の一周期における代表値を導出する異常渦電流損導出工程と、
前記異常渦電流損導出工程により導出された前記軟磁性材料における異常渦電流損の一周期における代表値と、前記物理量の大きさの一周期内の各時刻における値と、前記物理量の大きさの時間微分値の一周期内の各時刻における値と、に基づいて、教師データを作成する教師データ作成工程と、
前記教師データ作成工程により作成された前記教師データを用いて、ニューラルネットワークのパラメータを学習する学習工程と、を有し、
前記ニューラルネットワークの入力層には、或る時刻における前記物理量の大きさおよび当該時刻における前記物理量の大きさの時間微分値が少なくとも入力され、
前記ニューラルネットワークの出力層から、当該時刻における異常渦電流損が少なくとも出力され、
前記学習工程は、前記教師データに含まれる前記軟磁性材料の異常渦電流損の一周期における代表値と、前記ニューラルネットワークの前記出力層から出力される一周期内の各時刻における異常渦電流損の代表値との差の評価値を少なくとも求める目的関数を用いて、前記ニューラルネットワークのパラメータを学習することを特徴とする学習方法。
A learning method for learning neural network parameters for estimating abnormal eddy current loss in soft magnetic materials, the method comprising:
Measured values for the soft magnetic material when the soft magnetic material is excited under predetermined excitation conditions include measured values at each time within one period of a physical quantity obtained from the magnetization characteristics of the soft magnetic material, and the soft magnetic material. an acquisition step of acquiring a measured value of iron loss;
A first electromagnetic field analysis step of deriving representative values of hysteresis loss and classical eddy current loss in one cycle of the soft magnetic material when the soft magnetic material is excited under the predetermined excitation conditions based on Maxwell's equations. and,
Based on the measured value of iron loss acquired in the acquisition step and the representative values in one cycle of hysteresis loss and classical eddy current loss derived in the first electromagnetic field analysis step, abnormal eddy in the soft magnetic material is determined. An abnormal eddy current loss derivation step of deriving a representative value of current loss in one cycle;
A representative value in one cycle of the abnormal eddy current loss in the soft magnetic material derived by the abnormal eddy current loss deriving step, a value at each time within one cycle of the magnitude of the physical quantity, and a value of the magnitude of the physical quantity at each time within one cycle. A teacher data creation step of creating teacher data based on the value at each time within one cycle of the time differential value;
a learning step of learning parameters of a neural network using the teacher data created in the teacher data creation step;
At least the magnitude of the physical quantity at a certain time and the time differential value of the magnitude of the physical quantity at the time are input to the input layer of the neural network,
At least the abnormal eddy current loss at the time is output from the output layer of the neural network,
The learning step includes determining a representative value of the abnormal eddy current loss of the soft magnetic material in one cycle included in the teacher data and an abnormal eddy current loss at each time in one cycle output from the output layer of the neural network. A learning method characterized in that the parameters of the neural network are learned using an objective function that at least obtains an evaluation value of the difference between the 1 and the representative value.
請求項6に記載の学習方法でパラメータが学習された前記ニューラルネットワークを用いて軟磁性材料の異常渦電流損を推定する異常渦電流損推定方法であって、
前記異常渦電流損の推定対象の軟磁性材料を所定の励磁条件で励磁したときの当該軟磁性材料の前記物理量の大きさおよび当該物理量の大きさの時間微分値を導出する導出工程と、
前記導出工程により導出された或る時刻における前記物理量の大きさおよび当該時刻における前記物理量の大きさの時間微分値を前記ニューラルネットワークの前記入力層に与えることにより前記出力層から出力される値に基づいて、当該時刻における前記軟磁性材料の異常渦電流損を推定する推定工程と、を有することを特徴とする異常渦電流損推定方法。
An abnormal eddy current loss estimation method for estimating abnormal eddy current loss of a soft magnetic material using the neural network whose parameters have been learned by the learning method according to claim 6,
a derivation step of deriving the magnitude of the physical quantity of the soft magnetic material and the time differential value of the magnitude of the physical quantity when the soft magnetic material whose abnormal eddy current loss is to be estimated is excited under predetermined excitation conditions;
The magnitude of the physical quantity at a certain time derived in the derivation step and the time differential value of the magnitude of the physical quantity at the time are given to the input layer of the neural network to obtain the value output from the output layer. an estimation step of estimating the abnormal eddy current loss of the soft magnetic material at the relevant time based on the above.
請求項1または2に記載の学習装置の各手段としてコンピュータを機能させるためのプログラム。 A program for causing a computer to function as each means of the learning device according to claim 1 or 2. 請求項3~5の何れか1項に記載の異常渦電流損推定装置の各手段としてコンピュータを機能させるためのプログラム。 A program for causing a computer to function as each means of the abnormal eddy current loss estimating device according to any one of claims 3 to 5.
JP2020073951A 2020-04-17 2020-04-17 Learning device, abnormal eddy current loss estimation device, learning method, abnormal eddy current loss estimation method, and program Active JP7428893B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020073951A JP7428893B2 (en) 2020-04-17 2020-04-17 Learning device, abnormal eddy current loss estimation device, learning method, abnormal eddy current loss estimation method, and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020073951A JP7428893B2 (en) 2020-04-17 2020-04-17 Learning device, abnormal eddy current loss estimation device, learning method, abnormal eddy current loss estimation method, and program

Publications (2)

Publication Number Publication Date
JP2021170290A JP2021170290A (en) 2021-10-28
JP7428893B2 true JP7428893B2 (en) 2024-02-07

Family

ID=78119569

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020073951A Active JP7428893B2 (en) 2020-04-17 2020-04-17 Learning device, abnormal eddy current loss estimation device, learning method, abnormal eddy current loss estimation method, and program

Country Status (1)

Country Link
JP (1) JP7428893B2 (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002122571A (en) 2000-10-12 2002-04-26 Kenzo Miya Method and apparatus for inspection of defect
JP2008309573A (en) 2007-06-13 2008-12-25 Mitsubishi Heavy Ind Ltd Eddy current flaw detector and eddy current flaw detection method
JP2018109592A (en) 2017-01-06 2018-07-12 新日鐵住金株式会社 Device, method, and program for electromagnetic field analysis

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002122571A (en) 2000-10-12 2002-04-26 Kenzo Miya Method and apparatus for inspection of defect
JP2008309573A (en) 2007-06-13 2008-12-25 Mitsubishi Heavy Ind Ltd Eddy current flaw detector and eddy current flaw detection method
JP2018109592A (en) 2017-01-06 2018-07-12 新日鐵住金株式会社 Device, method, and program for electromagnetic field analysis

Also Published As

Publication number Publication date
JP2021170290A (en) 2021-10-28

Similar Documents

Publication Publication Date Title
JP6206608B1 (en) Electromagnetic field analysis apparatus, electromagnetic field analysis method, and program
Kowal et al. Comparison of iron loss models for electrical machines with different frequency domain and time domain methods for excess loss prediction
JP6834670B2 (en) Magnetic property measurement system and magnetic property measurement method
Gabi et al. FEM simulations of incremental permeability signals of a multi-layer steel with consideration of the hysteretic behavior of each layer
Knypiński et al. Application of a PSO algorithm for identification of the parameters of Jiles-Atherton hysteresis model
EP2420853B1 (en) Analysis apparatus and analysis method
JP7428893B2 (en) Learning device, abnormal eddy current loss estimation device, learning method, abnormal eddy current loss estimation method, and program
Kuczmann Dynamic Preisach hysteresis model
CN114580235A (en) Design method of vibration scaling model of converter transformer
JP6578845B2 (en) Electromagnetic field analysis apparatus, electromagnetic field analysis method, and program
He et al. An improved model for anomalous loss utilizing loss separation and comparison with ANN model in electrical steel sheet
Lorenz Electrical machine iron loss predictions-A unique engineering approach utilizing transient finite element methods-Part 1: Theory and calculation method
JP6610346B2 (en) Magnetic property analyzing apparatus, magnetic property analyzing method, and computer program
JP7124441B2 (en) Electromagnetic field analysis system, electromagnetic field analysis method, and program
Chothani et al. Convolution Neural Network and XGBoost-Based Fault Identification in Power Transformer
JP6984426B2 (en) Electromagnetic field analyzer, electromagnetic field analysis method, and program
Li et al. A posteriori error estimation and adaptive mesh refinement controlling in finite element analysis of 3-D steady state eddy current fields
Zurek et al. Extended T (x) hysteresis model with magnetic anisotropy
Gürbüz et al. A new methodology for incorporating the cutting deterioration of electrical sheets into electromagnetic finite-element simulation
Zhang et al. An improved modeling of vector magnetic properties of electrical steel sheet for FEM application and its experimental verification
Longhitano et al. Temperature-dependent hysteresis model for soft magnetic materials
JP2019125097A (en) Electromagnetic field analysis device, electromagnetic field analysis method, and program
Reljić et al. The estimation of iron losses in a non-oriented electrical steel sheet based on the artificial neural network and the genetic algorithm approaches
JP6417152B2 (en) Device power estimation method and system
Sofiane et al. Improving the Swelling Phenomenon in the Dynamic Jiles-Atherton Hysteresis Model Using Magnetic Viscosity.

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20221205

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231130

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20231226

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240108

R151 Written notification of patent or utility model registration

Ref document number: 7428893

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151