JP6645623B2 - 積層鉄心の弾性マトリックス決定方法および振動解析方法 - Google Patents
積層鉄心の弾性マトリックス決定方法および振動解析方法 Download PDFInfo
- Publication number
- JP6645623B2 JP6645623B2 JP2019524096A JP2019524096A JP6645623B2 JP 6645623 B2 JP6645623 B2 JP 6645623B2 JP 2019524096 A JP2019524096 A JP 2019524096A JP 2019524096 A JP2019524096 A JP 2019524096A JP 6645623 B2 JP6645623 B2 JP 6645623B2
- Authority
- JP
- Japan
- Prior art keywords
- elastic modulus
- core
- laminated core
- natural frequency
- laminated
- 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
Links
- 238000004458 analytical method Methods 0.000 title claims description 49
- 239000011159 matrix material Substances 0.000 title claims description 45
- 238000000034 method Methods 0.000 title claims description 31
- 238000012360 testing method Methods 0.000 claims description 49
- 229910000831 Steel Inorganic materials 0.000 claims description 37
- 239000010959 steel Substances 0.000 claims description 37
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 36
- 238000010030 laminating Methods 0.000 claims description 32
- 238000012916 structural analysis Methods 0.000 claims description 21
- 238000003475 lamination Methods 0.000 claims description 10
- 238000005316 response function Methods 0.000 claims description 6
- 238000005452 bending Methods 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 239000011162 core material Substances 0.000 description 73
- 229910001224 Grain-oriented electrical steel Inorganic materials 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 229910000576 Laminated steel Inorganic materials 0.000 description 7
- 238000003860 storage Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 239000011521 glass Substances 0.000 description 3
- 238000004804 winding Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 2
- 229910000976 Electrical steel Inorganic materials 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 239000000696 magnetic material Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/20—Metals
- G01N33/204—Structure thereof, e.g. crystal structure
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/12—Analysing solids by measuring frequency or resonance of acoustic waves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01F—MAGNETS; INDUCTANCES; TRANSFORMERS; SELECTION OF MATERIALS FOR THEIR MAGNETIC PROPERTIES
- H01F27/00—Details of transformers or inductances, in general
- H01F27/24—Magnetic cores
- H01F27/245—Magnetic cores made from sheets, e.g. grain-oriented
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01F—MAGNETS; INDUCTANCES; TRANSFORMERS; SELECTION OF MATERIALS FOR THEIR MAGNETIC PROPERTIES
- H01F27/00—Details of transformers or inductances, in general
- H01F27/33—Arrangements for noise damping
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/014—Resonance or resonant frequency
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0231—Composite or layered materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0234—Metals, e.g. steel
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02827—Elastic parameters, strength or force
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/26—Scanned objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01F—MAGNETS; INDUCTANCES; TRANSFORMERS; SELECTION OF MATERIALS FOR THEIR MAGNETIC PROPERTIES
- H01F1/00—Magnets or magnetic bodies characterised by the magnetic materials therefor; Selection of materials for their magnetic properties
- H01F1/01—Magnets or magnetic bodies characterised by the magnetic materials therefor; Selection of materials for their magnetic properties of inorganic materials
- H01F1/03—Magnets or magnetic bodies characterised by the magnetic materials therefor; Selection of materials for their magnetic properties of inorganic materials characterised by their coercivity
- H01F1/12—Magnets or magnetic bodies characterised by the magnetic materials therefor; Selection of materials for their magnetic properties of inorganic materials characterised by their coercivity of soft-magnetic materials
- H01F1/14—Magnets or magnetic bodies characterised by the magnetic materials therefor; Selection of materials for their magnetic properties of inorganic materials characterised by their coercivity of soft-magnetic materials metals or alloys
- H01F1/16—Magnets or magnetic bodies characterised by the magnetic materials therefor; Selection of materials for their magnetic properties of inorganic materials characterised by their coercivity of soft-magnetic materials metals or alloys in the form of sheets
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01F—MAGNETS; INDUCTANCES; TRANSFORMERS; SELECTION OF MATERIALS FOR THEIR MAGNETIC PROPERTIES
- H01F41/00—Apparatus or processes specially adapted for manufacturing or assembling magnets, inductances or transformers; Apparatus or processes specially adapted for manufacturing materials characterised by their magnetic properties
- H01F41/02—Apparatus or processes specially adapted for manufacturing or assembling magnets, inductances or transformers; Apparatus or processes specially adapted for manufacturing materials characterised by their magnetic properties for manufacturing cores, coils, or magnets
- H01F41/0206—Manufacturing of magnetic cores by mechanical means
- H01F41/0233—Manufacturing of magnetic circuits made from sheets
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Crystallography & Structural Chemistry (AREA)
- Theoretical Computer Science (AREA)
- Acoustics & Sound (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Manufacturing Cores, Coils, And Magnets (AREA)
Description
本発明は、変圧器鉄心などの積層鉄心の振動解析を実施する際に適用する鉄心の弾性変形における応力と歪の関係を示す構成式中の弾性マトリックス決定方法および振動解析方法に関する。
配電用変圧器などの変圧器は電磁鋼板を積層した積層鉄心にコイルを巻装することにより構成されている。変圧器として重要とされる性能には鉄損（無負荷損）特性、励磁電流特性、騒音特性などがある。
配電用変圧器は、様々な場所に設置されているが、特に市街地に設置される変圧器には騒音が小さいことが強く求められる。このように昨今では、変圧器が設置される周辺環境への配慮などから、特に騒音特性がますます重要となっている。
変圧器の鉄心材料として、方向性電磁鋼板が多く使用されている。方向性電磁鋼板には磁歪と称される励磁に伴う材料伸縮があり、この励磁磁歪振動が変圧器騒音の主な原因と言われている。このため変圧器騒音性能は使用する電磁鋼板の磁歪性能に強く依存するとされ、低騒音変圧器を製造するに際しては低磁歪特性を有する電磁鋼板が鉄心材料として使用される。
配電用変圧器は、様々な場所に設置されているが、特に市街地に設置される変圧器には騒音が小さいことが強く求められる。このように昨今では、変圧器が設置される周辺環境への配慮などから、特に騒音特性がますます重要となっている。
変圧器の鉄心材料として、方向性電磁鋼板が多く使用されている。方向性電磁鋼板には磁歪と称される励磁に伴う材料伸縮があり、この励磁磁歪振動が変圧器騒音の主な原因と言われている。このため変圧器騒音性能は使用する電磁鋼板の磁歪性能に強く依存するとされ、低騒音変圧器を製造するに際しては低磁歪特性を有する電磁鋼板が鉄心材料として使用される。
しかし、磁歪性能の優れた電磁鋼板を実際に使用して鉄心を製造したにもかかわらず、十分な変圧器低騒音特性が得られない場合がしばしばみられる。このようなことが起こる原因を調査してみると、変圧器鉄心の固有振動数と電磁鋼板の磁歪振動との共鳴現象であると考えられるケースが多くみられる。このため、変圧器鉄心の固有振動をはじめとする機械振動特性を計算予測することは変圧器を設計・製造する上で極めて重要である。
そこで、磁歪が生ずる磁性体を含む電磁部品を有限要素解析における複数の有限要素の組み合わせで表現した数値解析モデルに基づいて電磁部品に与えられる磁束密度に応じた有限要素の各節点又は各有限要素の歪みと等価な節点力を算出する解析装置および解析方法が提案されている（例えば、特許文献１参照）。
そこで、磁歪が生ずる磁性体を含む電磁部品を有限要素解析における複数の有限要素の組み合わせで表現した数値解析モデルに基づいて電磁部品に与えられる磁束密度に応じた有限要素の各節点又は各有限要素の歪みと等価な節点力を算出する解析装置および解析方法が提案されている（例えば、特許文献１参照）。
しかしながら、上記特許文献１に記載された先行技術では、力の釣合いの式、応力と歪みとの関係を示した構成式、および変位と歪みの関係式によって構成される構造解析の支配方程式を用いて、準静的構造解析を行うようにしている。
このうち、応力テンソル｛σ｝と歪みテンソル｛ε｝との関係を示した構成式は、
｛σ｝＝｛Ｄ｝｛ε｝ （｛ ｝はテンソルを示す）
で表されている。
ここで、｛Ｄ｝は歪みと応力の関係を表したテンソルである。成分表示すると、（１）式のようになる。
このうち、応力テンソル｛σ｝と歪みテンソル｛ε｝との関係を示した構成式は、
｛σ｝＝｛Ｄ｝｛ε｝ （｛ ｝はテンソルを示す）
で表されている。
ここで、｛Ｄ｝は歪みと応力の関係を表したテンソルである。成分表示すると、（１）式のようになる。
式中、｛Ｄ｝は８１成分、｛σ｝と｛ε｝はそれぞれ９成分ある。ここで、物理量としてのテンソルは対称テンソルとなるので、｛σ｝と｛ε｝の独立成分はそれぞれ６成分となる。したがって、構成式を
［σ］＝［Ｄ］［ε］ （［ ］は行列であることを示す。［Ｄ］を弾性マトリックスと称する。）
と行列表示して、成分で表示すると、（２）式のようになる。
［σ］＝［Ｄ］［ε］ （［ ］は行列であることを示す。［Ｄ］を弾性マトリックスと称する。）
と行列表示して、成分で表示すると、（２）式のようになる。
垂直応力σｉと垂直歪みεｉとの関係が弾性マトリックス［Ｄ］を使用して表され、剪断歪みτｉｊとせん断歪みγｉｊとの関係がせん断弾性係数Ｇを用いて表されている。
しかしながら、上記先行技術では、弾性マトリックス［Ｄ］やせん断弾性係数Ｇの設定については何ら記載がなく、弾性マトリックス［Ｄ］やせん断弾性係数Ｇをどのように決定するかについては記載がない。但し、一般的には、構造解析の対象となる電磁部品を構成する電磁鋼板等の部材自体の弾性係数をそのまま適用する場合が多い。
しかしながら、上記先行技術では、弾性マトリックス［Ｄ］やせん断弾性係数Ｇの設定については何ら記載がなく、弾性マトリックス［Ｄ］やせん断弾性係数Ｇをどのように決定するかについては記載がない。但し、一般的には、構造解析の対象となる電磁部品を構成する電磁鋼板等の部材自体の弾性係数をそのまま適用する場合が多い。
この場合には、構造解析の対象となる電磁部品の機械振動計算結果と実際に機械振動を測定した結果とを比較した場合に、計算値と実測値との間に大きな乖離があることが知られている。
計算値と実測値との乖離を抑制するためには、構造解析の対象となる電磁部品を試作して弾性係数を測定する必要があるが、それでは本末転倒であり、構造解析プログラムによる構造解析を電磁部品の設計に反映することは困難である。
そこで、本発明は、上述した従来技術の課題に着目してなされたものであり、積層鉄心の弾性係数を最適に決定することができる積層鉄心の弾性マトリックス決定方法および振動解析方法を提供することを目的としている。
計算値と実測値との乖離を抑制するためには、構造解析の対象となる電磁部品を試作して弾性係数を測定する必要があるが、それでは本末転倒であり、構造解析プログラムによる構造解析を電磁部品の設計に反映することは困難である。
そこで、本発明は、上述した従来技術の課題に着目してなされたものであり、積層鉄心の弾性係数を最適に決定することができる積層鉄心の弾性マトリックス決定方法および振動解析方法を提供することを目的としている。
上記の課題を解決するために、本発明の趣旨は以下のとおりである。
［１］応力と歪みとの関係を行列表示で表した構成式を使用して鋼板を積層した積層鉄心の振動解析を行うに際し、前記構成式中の弾性マトリックスに含まれる、前記積層鉄心の積層方向を含む面の横弾性係数、及び、前記積層鉄心の積層方向の縦弾性係数を決定する弾性マトリックス決定方法であって、振動解析の対象となる前記積層鉄心と同じ鋼板を積層してなる試験用鉄心を製作する第１ステップと、前記第１ステップで製作した前記試験用鉄心を支持して積層方向に加振することにより曲げ振動の固有振動数を実測する第２ステップと、前記試験用鉄心の構造解析モデルについて、積層方向の縦弾性係数として前記試験用鉄心に使用した鋼板１枚の縦弾性係数を近似値として与え、積層方向を含む面の横弾性係数をパラメータとして値を変えながら複数回構造解析を行って、前記積層方向を含む面の横弾性係数と固有振動数との関係を推定する第３ステップと、前記第３ステップで推定した前記積層方向を含む面の横弾性係数と固有振動数との関係に、前記第２ステップで実測した前記試験用鉄心の固有振動数の実測値を与えて前記試験用鉄心の前記積層方向を含む面の横弾性係数の近似値を決定する第４ステップと、前記試験用鉄心の構造解析モデルについて、前記積層方向を含む面の横弾性係数として前記第４ステップで決定した前記積層方向を含む面の横弾性係数の近似値を与え、前記積層方向の縦弾性係数をパラメータとして値を変えながら複数回構造解析を行って、前記積層方向の縦弾性係数と固有振動数との関係を推定する第５ステップと、前記第５ステップで推定した前記積層方向の縦弾性係数と固有振動数との関係に、前記第２ステップで実測した前記試験用鉄心の固有振動数の実測値を与えて前記試験用鉄心の前記積層方向の縦弾性係数の近似値を決定する第６ステップと、前記試験用鉄心の構造解析モデルについて、前記積層方向の縦弾性係数として前記第６ステップで決定した前記積層方向の縦弾性係数の近似値を与えて、前記第３ステップから前記第６ステップの処理を実行することを繰り返すことにより、前記積層鉄心の積層方向を含む面の横弾性係数、及び、前記積層鉄心の積層方向の縦弾性係数を決定する第７ステップと、を含むことを特徴とする弾性マトリックス決定方法。
［１］応力と歪みとの関係を行列表示で表した構成式を使用して鋼板を積層した積層鉄心の振動解析を行うに際し、前記構成式中の弾性マトリックスに含まれる、前記積層鉄心の積層方向を含む面の横弾性係数、及び、前記積層鉄心の積層方向の縦弾性係数を決定する弾性マトリックス決定方法であって、振動解析の対象となる前記積層鉄心と同じ鋼板を積層してなる試験用鉄心を製作する第１ステップと、前記第１ステップで製作した前記試験用鉄心を支持して積層方向に加振することにより曲げ振動の固有振動数を実測する第２ステップと、前記試験用鉄心の構造解析モデルについて、積層方向の縦弾性係数として前記試験用鉄心に使用した鋼板１枚の縦弾性係数を近似値として与え、積層方向を含む面の横弾性係数をパラメータとして値を変えながら複数回構造解析を行って、前記積層方向を含む面の横弾性係数と固有振動数との関係を推定する第３ステップと、前記第３ステップで推定した前記積層方向を含む面の横弾性係数と固有振動数との関係に、前記第２ステップで実測した前記試験用鉄心の固有振動数の実測値を与えて前記試験用鉄心の前記積層方向を含む面の横弾性係数の近似値を決定する第４ステップと、前記試験用鉄心の構造解析モデルについて、前記積層方向を含む面の横弾性係数として前記第４ステップで決定した前記積層方向を含む面の横弾性係数の近似値を与え、前記積層方向の縦弾性係数をパラメータとして値を変えながら複数回構造解析を行って、前記積層方向の縦弾性係数と固有振動数との関係を推定する第５ステップと、前記第５ステップで推定した前記積層方向の縦弾性係数と固有振動数との関係に、前記第２ステップで実測した前記試験用鉄心の固有振動数の実測値を与えて前記試験用鉄心の前記積層方向の縦弾性係数の近似値を決定する第６ステップと、前記試験用鉄心の構造解析モデルについて、前記積層方向の縦弾性係数として前記第６ステップで決定した前記積層方向の縦弾性係数の近似値を与えて、前記第３ステップから前記第６ステップの処理を実行することを繰り返すことにより、前記積層鉄心の積層方向を含む面の横弾性係数、及び、前記積層鉄心の積層方向の縦弾性係数を決定する第７ステップと、を含むことを特徴とする弾性マトリックス決定方法。
［２］前記積層鉄心は、電磁鋼板を複数枚積層して構成されていることを特徴とする［１］に記載の弾性マトリックス決定方法。
［３］前記積層鉄心は、変圧器用積層鉄心であることを特徴とする［２］に記載の弾性マトリックス決定方法。
［３］前記積層鉄心は、変圧器用積層鉄心であることを特徴とする［２］に記載の弾性マトリックス決定方法。
［４］応答関数を用いて応力と歪みとの関係をテンソル表示で表した構成式の弾性マトリックスに、［１］に記載の弾性マトリックス決定方法で決定した前記積層鉄心の積層方向を含む面の横弾性係数、及び、前記積層鉄心の積層方向の縦弾性係数を組み込んで鋼板を積層した積層鉄心の振動解析を行うことを特徴とする積層鉄心の振動解析方法。
［５］前記積層鉄心は、電磁鋼板を複数枚積層した変圧器用鉄心であることを特徴とする［４］に記載の積層鉄心の振動解析方法。
［５］前記積層鉄心は、電磁鋼板を複数枚積層した変圧器用鉄心であることを特徴とする［４］に記載の積層鉄心の振動解析方法。
本発明に係る積層鉄心の弾性マトリックス決定方法によれば、振動解析に使用する応力と歪みとの関係を表す構成式の弾性係数を積層鋼板間の滑りを考慮して適正に設定して振動特性の実測値と計算値との乖離を抑制することができる。
また、本発明に係る積層鉄心の振動解析方法は、適正に設定した弾性係数を構成式に組み込んで振動解析を行うことにより、振動解析精度を向上させることができる。
また、本発明に係る積層鉄心の振動解析方法は、適正に設定した弾性係数を構成式に組み込んで振動解析を行うことにより、振動解析精度を向上させることができる。
以下、本発明に係る積層鉄心の弾性マトリックス決定方法および積層鉄心の振動解析方法の実施形態を図面に基づいて説明する。なお、各図面は模式的なものであって、現実のものとは異なる場合がある。また、以下の実施形態は、本発明の技術的思想を具体化するための装置や方法を例示するものであり、構成を下記のものに特定するものでない。すなわち、本発明の技術的思想は、特許請求の範囲に記載された技術的範囲内において、種々の変更を加えることができる。
振動解析装置１０は、図１に示すように、ＣＰＵ１１を備えた演算処理装置１２で構成されている。ＣＰＵ１１には、内部バス１３を介してＲＡＭ，ＲＯＭ等の内部記憶装置１４、外部記憶装置１５、キーボード、マウス等の入力装置１６およびディスプレイに画像データを出力する出力装置１７が接続されている。
外部記憶装置１５は、ハードディスクドライブやソリッドステートドライブ等の読み出しが可能なディスクドライブと、記録媒体からのデータを読み出すＣＤ、ＤＶＤ、ＢＤ等のドライブ装置を含んで構成されている。この外部記憶装置１５に市販の振動解析プログラムを格納した記録媒体１８をセットし、読み出した振動解析プログラムをディスクドライブにインストールする。なお、振動解析プログラムのインストールは、記録媒体１８を使用する場合に限らず、ネットワークを介して振動解析プログラムをダウンロードするようにしてもよい。
外部記憶装置１５は、ハードディスクドライブやソリッドステートドライブ等の読み出しが可能なディスクドライブと、記録媒体からのデータを読み出すＣＤ、ＤＶＤ、ＢＤ等のドライブ装置を含んで構成されている。この外部記憶装置１５に市販の振動解析プログラムを格納した記録媒体１８をセットし、読み出した振動解析プログラムをディスクドライブにインストールする。なお、振動解析プログラムのインストールは、記録媒体１８を使用する場合に限らず、ネットワークを介して振動解析プログラムをダウンロードするようにしてもよい。
ＣＰＵ１１は、インストールされた振動解析プログラムにしたがって入力された解析用入力データを用いて有限要素法を用いて振動解析を行い、振動モードに応じた固有振動数等の解析結果を出力装置１７からディスプレイに出力して表示する。この解析結果は、ディスプレイに表示する場合に限らず、プリンタで印刷したり、ネットワーク経由で送信したりすることができる。
そして、本実施形態で解析対象とする積層鉄心２１は、例えば配電用変圧器として使用する三相三脚変圧器用の積層鉄心であって、図２に示すように、上ヨーク２２ａおよび下ヨーク２２ｂ間に三本の脚部２２ｃを連結した板厚０．３ｍｍの方向性電磁鋼板２２を例えば３３３枚積層してガラステープを巻き付けて固定されている。
そして、本実施形態で解析対象とする積層鉄心２１は、例えば配電用変圧器として使用する三相三脚変圧器用の積層鉄心であって、図２に示すように、上ヨーク２２ａおよび下ヨーク２２ｂ間に三本の脚部２２ｃを連結した板厚０．３ｍｍの方向性電磁鋼板２２を例えば３３３枚積層してガラステープを巻き付けて固定されている。
一例として、上ヨーク２２ａおよび下ヨーク２２ｂの寸法は、幅１００ｍｍ×長さ５００ｍｍに設定されている。また、三本の脚部２２ｃの寸法は、幅１００ｍｍ×長さ３００ｍｍに設定され、上ヨーク２２ａおよび下ヨーク２２ｂ間に１００ｍｍ間隔で連結されている。
このような三相三脚変圧器の積層鉄心２１の振動解析を行う場合には、構造解析の支配方程式となる応力と歪みとの関係を示した構成式が使用される。
この構成式は、積層物を等価均質体に置き換え、積層の影響を弾性マトリックスの値で表現すると下記（３）式のようになる。
このような三相三脚変圧器の積層鉄心２１の振動解析を行う場合には、構造解析の支配方程式となる応力と歪みとの関係を示した構成式が使用される。
この構成式は、積層物を等価均質体に置き換え、積層の影響を弾性マトリックスの値で表現すると下記（３）式のようになる。
［σ］＝［Ｃ］［ε］ ・・・（３）
ここで、［σ］は応力マトリックス、［Ｃ］は応答関数としての弾性マトリックス（ステフィネスマトリックス）、［ε］は歪みマトリックスである。
ここで、鋼板の積層方向をＺ方向とし、このＺ方向と直行する２次元平面の一方をＸ方向、他方をＹ方向とすると、応力マトリックス［σ］は、図３に示すように、垂直成分がＸ方向の垂直応力σｘ、Ｙ方向の垂直応力σｙおよびＺ方向の垂直応力σｚで表され、せん断成分が、ＺＸ平面のせん断応力τｚｘ、ＹＺ平面のせん断応力τｙｚおよびＸＹ平面のせん断応力τｘｙで表される。
ここで、［σ］は応力マトリックス、［Ｃ］は応答関数としての弾性マトリックス（ステフィネスマトリックス）、［ε］は歪みマトリックスである。
ここで、鋼板の積層方向をＺ方向とし、このＺ方向と直行する２次元平面の一方をＸ方向、他方をＹ方向とすると、応力マトリックス［σ］は、図３に示すように、垂直成分がＸ方向の垂直応力σｘ、Ｙ方向の垂直応力σｙおよびＺ方向の垂直応力σｚで表され、せん断成分が、ＺＸ平面のせん断応力τｚｘ、ＹＺ平面のせん断応力τｙｚおよびＸＹ平面のせん断応力τｘｙで表される。
同様に、歪みマトリックス［ε］は、垂直成分がＸ方向の垂直歪みεｘ、Ｙ方向の垂直歪みεｙおよびＺ方向の垂直歪みεｚで表される、せん断成分がＺＸ平面のせん断歪みはγｚｘ、ＹＺ平面のせん断歪みγｙｚおよびＸＹ平面のせん断歪みγｘｙで表される。
また、弾性マトリックス［Ｃ］は、３６個の弾性係数Ｃｉｊ（ｉ＝１〜６，ｊ＝１〜６）で表される。
これらをマトリックス表示すると下記（４）式となる。
また、弾性マトリックス［Ｃ］は、３６個の弾性係数Ｃｉｊ（ｉ＝１〜６，ｊ＝１〜６）で表される。
これらをマトリックス表示すると下記（４）式となる。
積層鉄心は、方向性電磁鋼板を積層して製造するので、積層鉄心の機械的対称性を有する他、積層する鋼板の長手方向とその直角方向にも１８０度対称性を有するので、異方性分類としては直交異方性を有することになる。このため、直交異方性を有する物体については基本的に、下記（５）式のようにＣ１１、Ｃ１２、Ｃ１３、Ｃ２２、Ｃ２３、Ｃ３３、Ｃ４４、Ｃ５５およびＣ６６の計９個の弾性係数で表すことができる。
このうち、弾性係数Ｃ１１、Ｃ１２、Ｃ１３、Ｃ２２、Ｃ２３、Ｃ３３については縦弾性係数Ｅｘ、ＥｙおよびＥｚとポアソン比νｘｙ、νｙｘ、νｙｚ、νｚｙ、νｚｘ、νｘｚとによって（６）〜（１２）式で算出することができる。
ここで、ＥｘはＸ方向縦弾性係数（ヤング率）、ＥｙはＹ方向縦弾性係数（ヤング率）、νｘｙはＸＹ平面のポアソン比、νｙｘはＹＸ平面のポアソン比、νｙｚはＹＺ平面のポアソン比、νｚｙはＺＹ平面のポアソン比、νｚｘはＺＸ平面のポアソン比、νｘｚはＸＺ平面のポアソン比である。また、弾性係数Ｃ４４はＹＺ平面の横弾性係数Ｇｙｚであり、弾性係数Ｃ５５はＺＸ平面の横弾性係数Ｇｚｘであり、弾性係数Ｃ６６はＸＹ平面の横弾性係数Ｇｘｙである。 そして、縦弾性係数ＥｘおよびＥｙは１枚の鋼板の縦弾性係数Ｅｘ０、Ｅｙ０と等しく設定することができるが、縦弾性係数Ｅｚは１枚の鋼板の縦弾性係数Ｅｚ０と略等しく設定できない。その理由は、積層鋼板間には、わずかながら隙間があるからである。
ここで、ＥｘはＸ方向縦弾性係数（ヤング率）、ＥｙはＹ方向縦弾性係数（ヤング率）、νｘｙはＸＹ平面のポアソン比、νｙｘはＹＸ平面のポアソン比、νｙｚはＹＺ平面のポアソン比、νｚｙはＺＹ平面のポアソン比、νｚｘはＺＸ平面のポアソン比、νｘｚはＸＺ平面のポアソン比である。また、弾性係数Ｃ４４はＹＺ平面の横弾性係数Ｇｙｚであり、弾性係数Ｃ５５はＺＸ平面の横弾性係数Ｇｚｘであり、弾性係数Ｃ６６はＸＹ平面の横弾性係数Ｇｘｙである。 そして、縦弾性係数ＥｘおよびＥｙは１枚の鋼板の縦弾性係数Ｅｘ０、Ｅｙ０と等しく設定することができるが、縦弾性係数Ｅｚは１枚の鋼板の縦弾性係数Ｅｚ０と略等しく設定できない。その理由は、積層鋼板間には、わずかながら隙間があるからである。
また、ＸＹ平面の横弾性係数Ｇｘｙも１枚の鋼板の横弾性係数Ｇｘｙ０と等しく設定できるが、ＺＸ平面の横弾性係数ＧｚｘおよびＹＺ平面の横弾性係数Ｇｙｚは１枚の鋼板の横弾性係数Ｇｚｘ０およびＧｙｚをそのまま設定することはできない。
その理由は、積層鋼板は、積層された各鋼板の境界面で積層方向と直交するＸ方向およびＹ方向に滑りが生じることから、鋼板間の滑りの影響を横弾性係数ＧｚｘおよびＧｙｚに反映する必要があるからである。
したがって、弾性マトリックスにおいて、縦弾性係数ＥｚおよびＺＸ平面の横弾性係数ＧｚｘおよびＹＺ平面の横弾性係数Ｇｙｚを決定することが、積層鉄心の応力と歪みの関係式を使用する振動解析において重要事項となる。
その理由は、積層鋼板は、積層された各鋼板の境界面で積層方向と直交するＸ方向およびＹ方向に滑りが生じることから、鋼板間の滑りの影響を横弾性係数ＧｚｘおよびＧｙｚに反映する必要があるからである。
したがって、弾性マトリックスにおいて、縦弾性係数ＥｚおよびＺＸ平面の横弾性係数ＧｚｘおよびＹＺ平面の横弾性係数Ｇｙｚを決定することが、積層鉄心の応力と歪みの関係式を使用する振動解析において重要事項となる。
ところが、ＺＸ平面の横弾性係数ＧｚｘおよびＹＺ平面の横弾性係数Ｇｙｚを鋼板間の滑りの影響を反映した値とするには、実際に三相三脚変圧器用の積層鉄心を製作して、正確な横弾性係数ＧｚｘおよびＧｙｚを測定しなければならない。しかしながら、製作した三相三脚変圧器用の積層鉄心の実測された固有振動数が、目標とする固有振動数と合わない場合には、鉄心を再度製作する必要がある。この過程を繰り返して固有振動数が目標とする固有振動数に一致する鉄心を製作することは困難である。
そこで、本実施形態では、図４に示す手順で弾性マトリックス決定方法を実施することにより、積層された電磁鋼板間の積層方向の縦弾性率および滑りを考慮した横弾性係数を決定する。
そこで、本実施形態では、図４に示す手順で弾性マトリックス決定方法を実施することにより、積層された電磁鋼板間の積層方向の縦弾性率および滑りを考慮した横弾性係数を決定する。
弾性マトリックスを決定するには、鉄心そのものを製作するのではなく、まず鉄心と同じ電磁鋼板を用いた試験用鉄心を製作する（ステップＳ１）。例えば厚さ０．３ｍｍ、幅３０ｍｍ、長さ２８０ｍｍの方向性電磁鋼板を、例えば２０枚積層してガラステープを巻き付けて固定した試験用直方体鉄心を製作する。
ここで、試験用鉄心の形状は直方体に限られるものではない。直方体だけでなく、角が丸みを帯びていたり、一部にノッチや切り欠けがあったりしてもよく、必ずしも直方体となっていなくてもよい。また、楕円状の電磁鋼板を積層した鉄心でもよい。
ここで、試験用鉄心の形状は直方体に限られるものではない。直方体だけでなく、角が丸みを帯びていたり、一部にノッチや切り欠けがあったりしてもよく、必ずしも直方体となっていなくてもよい。また、楕円状の電磁鋼板を積層した鉄心でもよい。
以降、試験用鉄心を直方体鉄心として、本発明の詳細を説明する。
次いで、作成した試験用直方体鉄心と加振器との間に加速度と力を測定するためのインピーダンスヘッドを取り付け、試験用直方体鉄心を中央支持して固有振動数測定を実施する（ステップＳ２）。加振方向を鋼板積層方向として、中央支持（両端自由）で曲げ振動の１次モードおよび高次モードの固有振動数を求める。中央支持の場合は、高次モードは３次モードになる。
この時の支持位置は、中央支持（両端自由）としてもよいし、両端支持（中央自由）としてもよい。さらに、試験用鉄心の一方の端を固定して、もう一方の端を自由としてもよい。
次いで、作成した試験用直方体鉄心と加振器との間に加速度と力を測定するためのインピーダンスヘッドを取り付け、試験用直方体鉄心を中央支持して固有振動数測定を実施する（ステップＳ２）。加振方向を鋼板積層方向として、中央支持（両端自由）で曲げ振動の１次モードおよび高次モードの固有振動数を求める。中央支持の場合は、高次モードは３次モードになる。
この時の支持位置は、中央支持（両端自由）としてもよいし、両端支持（中央自由）としてもよい。さらに、試験用鉄心の一方の端を固定して、もう一方の端を自由としてもよい。
次に、振動解析装置１０にインストールされた構造解析ソフトを用いて試験用鉄心の固有振動解析を実施する。まず縦弾性係数Ｅｚを、試験用鉄心に使用する1枚の鋼板のそれと等しいと仮定して、縦弾性係数Ｅｚの１次近似値をＥｚ１とする（ステップＳ３）。
そして、鉄心の横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇとして横弾性係数Ｇをパラメータとして値を変えながら複数回構造解析を実施した結果に基づいて横弾性係数Ｇと固有振動数との関係を推定し、図６（ａ）に示す特性線図を作成する。
この特性線図は、横軸に横弾性係数Ｇをとり、縦軸に固有振動数をとったもので、両端自由の曲げ振動モードの１次モードおよび３次モードの試験用鉄心の横弾性係数Ｇと固有振動数との関係を表している。
そして、鉄心の横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇとして横弾性係数Ｇをパラメータとして値を変えながら複数回構造解析を実施した結果に基づいて横弾性係数Ｇと固有振動数との関係を推定し、図６（ａ）に示す特性線図を作成する。
この特性線図は、横軸に横弾性係数Ｇをとり、縦軸に固有振動数をとったもので、両端自由の曲げ振動モードの１次モードおよび３次モードの試験用鉄心の横弾性係数Ｇと固有振動数との関係を表している。
次いで、図６（ａ）の特性線図に第１モードおよび第３モードの固有振動数実測値を与えて横弾性係数Ｇをそれぞれ算出することにより、横弾性係数Ｇの1次近似値を決定することができる（ステップＳ４）。
例えば、第１モードおよび第３モードのそれぞれの横弾性係数Ｇ１ａおよびＧ１ｂの平均値を横弾性係数Ｇの1次近似値Ｇ１（Ｇ１＝（Ｇ１ａ＋Ｇ１ｂ）／２）とすることができる。
次に、ステップＳ４で求められた横弾性係数Ｇの１次近似値Ｇ１を使用して、縦弾性係数Ｅｚの２次近似値を求める。これは横弾性係数Ｇの近似値を求める場合と同様であり、横弾性係数Ｇを１次近似値Ｇ１として、縦弾性係数Ｅｚをパラメータとして値を変えながら複数回構造解析を実施した結果に基づいて縦弾性係数Ｅｚと固有振動数との関係を推定し、図６（ｂ）に示す特性線図を作成する。図６（ｂ）の特性線図に第１モードおよび第３モードの固有振動数実測値を与えて縦弾性係数Ｅｚの２次近似値を決定することができる（ステップＳ５）。
例えば、第１モードおよび第３モードのそれぞれの横弾性係数Ｇ１ａおよびＧ１ｂの平均値を横弾性係数Ｇの1次近似値Ｇ１（Ｇ１＝（Ｇ１ａ＋Ｇ１ｂ）／２）とすることができる。
次に、ステップＳ４で求められた横弾性係数Ｇの１次近似値Ｇ１を使用して、縦弾性係数Ｅｚの２次近似値を求める。これは横弾性係数Ｇの近似値を求める場合と同様であり、横弾性係数Ｇを１次近似値Ｇ１として、縦弾性係数Ｅｚをパラメータとして値を変えながら複数回構造解析を実施した結果に基づいて縦弾性係数Ｅｚと固有振動数との関係を推定し、図６（ｂ）に示す特性線図を作成する。図６（ｂ）の特性線図に第１モードおよび第３モードの固有振動数実測値を与えて縦弾性係数Ｅｚの２次近似値を決定することができる（ステップＳ５）。
例えば、第１モードおよび第３モードのそれぞれの縦弾性係数Ｅｚ２ａおよびＥｚ２ｂの平均値を縦弾性係数Ｅｚの２次近似値Ｅｚ２（Ｅｚ２＝（Ｅｚ２ａ＋Ｅｚ２ｂ）／２）とすることができる。
そして、ステップＳ５で決定した縦弾性係数Ｅｚのｎ次近似値をステップＳ４の構造解析の縦弾性係数Ｅｚとして横弾性係数Ｇの（ｎ＋１）次近似値を決定し、決定した横弾性係数Ｇの（ｎ＋１）次近似値を使用してステップＳ５で縦弾性係数Ｅｚの（ｎ＋２）次近似値を決定する操作を繰り返すことで、固有振動数測定値は不変であるので、第１モードおよび第３モードでの固有振動数の実測値に一致する又は略一致する縦弾性係数Ｅｚおよび横弾性係数Ｇを決定することができる（ステップＳ６）。
そして、ステップＳ５で決定した縦弾性係数Ｅｚのｎ次近似値をステップＳ４の構造解析の縦弾性係数Ｅｚとして横弾性係数Ｇの（ｎ＋１）次近似値を決定し、決定した横弾性係数Ｇの（ｎ＋１）次近似値を使用してステップＳ５で縦弾性係数Ｅｚの（ｎ＋２）次近似値を決定する操作を繰り返すことで、固有振動数測定値は不変であるので、第１モードおよび第３モードでの固有振動数の実測値に一致する又は略一致する縦弾性係数Ｅｚおよび横弾性係数Ｇを決定することができる（ステップＳ６）。
次に、振動解析装置１０にインストールされた構造解析ソフトを用いて図２に示す形状の三相三脚変圧器鉄心の固有振動解析を実施する（ステップＳ７）。三相三脚変圧器鉄心の機械的弾性係数のうち縦弾性係数Ｅｚ、横弾性係数ＧｙｚおよびＧｚｘを上記弾性マトリックス決定方法で決定した値とし、その他の機械的弾性係数の値は方向性電磁鋼板の機械的弾性係数として前述した（５）式に組み込んで計算を行う。
このようにして、図５に示す中央磁脚に対して左右磁脚が逆方向に捩じれるモードＡと、３つの磁脚が座屈するモードＢと、上ヨークが後方に湾曲し、下ヨークが前方に座屈するモードＣの３つの固有振動モードが存在することが分かり、各モードでの固有振動数を算出することができる。
このようにして、図５に示す中央磁脚に対して左右磁脚が逆方向に捩じれるモードＡと、３つの磁脚が座屈するモードＢと、上ヨークが後方に湾曲し、下ヨークが前方に座屈するモードＣの３つの固有振動モードが存在することが分かり、各モードでの固有振動数を算出することができる。
上記の説明では、試験用直方体鉄心の固有振動数測定を中央支持で実施したので、固有振動は第１モードおよび第３モードとしたが、支持の方法によって振動モードは異なる。たとえば、両端支持であれば、第１モードおよび第２モードとなる。したがって、固有振動の２次以上の高次モードは、支持の方法によって適宜選択してよい。
また、上記の説明では、縦弾性係数Ｅｚの１次近似値を仮定して、横弾性係数Ｇと固有振動数の関係から横弾性係数Ｇの１次近似値を求め、それをもとにして縦弾性係数Ｅｚの２次近似値を求めるようにしたが、この順序を逆にして、横弾性係数Ｇの１次近似値を固定して縦弾性係数Ｅｚの１次近似値を求め、それをもとに横弾性係数Ｇの２次近似値を求めるようにしてもよいことはいうまでもない。
また、上記の説明では、縦弾性係数Ｅｚの１次近似値を仮定して、横弾性係数Ｇと固有振動数の関係から横弾性係数Ｇの１次近似値を求め、それをもとにして縦弾性係数Ｅｚの２次近似値を求めるようにしたが、この順序を逆にして、横弾性係数Ｇの１次近似値を固定して縦弾性係数Ｅｚの１次近似値を求め、それをもとに横弾性係数Ｇの２次近似値を求めるようにしてもよいことはいうまでもない。
先ず、解析対象となる三相三脚変圧器鉄心に使用する厚さ０．３ｍｍの方向性電磁鋼板と同じ方向性電磁鋼板を、幅３０ｍｍ、長さ２８０ｍｍの長方形に成形し、この方向性電磁鋼板を２０枚積層してガラステープを巻き付けて固定した試験用直方体鉄心を製作した。
作成した試験用直方体鉄心と加振器との間にインピーダンスヘッドを取り付け、試験用直方体鉄心を中央支持して固有振動数測定を実施した。加振方向は鋼板積層方向として中央支持（両端自由）の曲げ振動の１次モードおよび３次モードの固有振動数実測値Ｆ１およびＦ３を求めた。
作成した試験用直方体鉄心と加振器との間にインピーダンスヘッドを取り付け、試験用直方体鉄心を中央支持して固有振動数測定を実施した。加振方向は鋼板積層方向として中央支持（両端自由）の曲げ振動の１次モードおよび３次モードの固有振動数実測値Ｆ１およびＦ３を求めた。
一方、試験用直方体鉄心の縦弾性係数Ｅｚを、前記試験用直方体鉄心に使用する１枚の鋼板の縦弾性係数（Ｅｚ１＝２２０ＧＰａ）に設定し、前記試験用直方体の横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇとして、横弾性係数Ｇをパラメータとして値を変えながら複数回構造解析を実施して、図６（ａ）および（ｂ）に示す特性曲線を作成した。
作成した特性曲線を使用して１次モードの固有振動数実測値Ｆ１および３次モードの固有振動数実測値Ｆ３から横弾性係数Ｇの１次近似値Ｇ１を決定した。この場合、Ｇ１ａ＝０．２５、Ｇ１ｂ＝０．４５であった。この結果、方向性電磁鋼板間の滑りを考慮した前記試験用直方体の横弾性係数Ｇの１次近似値としてＧ１＝０．３５［ＧＰａ］を決定することができた。
作成した特性曲線を使用して１次モードの固有振動数実測値Ｆ１および３次モードの固有振動数実測値Ｆ３から横弾性係数Ｇの１次近似値Ｇ１を決定した。この場合、Ｇ１ａ＝０．２５、Ｇ１ｂ＝０．４５であった。この結果、方向性電磁鋼板間の滑りを考慮した前記試験用直方体の横弾性係数Ｇの１次近似値としてＧ１＝０．３５［ＧＰａ］を決定することができた。
次に、試験用直方体鉄心の横弾性係数Ｇの１次近似値をＧ１＝０．３５［ＧＰａ］として、前記試験用直方体の縦弾性係数Ｅｚをパラメータとして値を変えながら複数回構造解析を実施して、１次モードおよび３次モードを示す特性曲線（図６（ｂ））を作成した。Ｅｚ２ａ＝８、Ｅｚ２ｂ＝１４であった。そして、方向性電磁鋼板間の滑りを考慮した前記試験用直方体の縦弾性係数Ｅｚの２次近似値として、Ｅｚ２＝１１［ＧＰａ］を得ることができた。
次に、試験用直方体鉄心の縦弾性係数Ｅｚの２次近似値をＥｚ２＝１１［ＧＰａ］として、横弾性係数Ｇの近似値を決定し、ここで求められた新たな横弾性係数Ｇの近似値を用いて、縦弾性係数Ｅｚの近似値を求めることを繰り返し、その結果、前記試験用直方体の固有振動数算出値が固有振動数実測値に一致する縦弾性係数および横弾性係数として、Ｅｚ＝１０［ＧＰａ］およびＧ＝０．２９［ＧＰａ］をそれぞれ得た。
次に、試験用直方体鉄心の縦弾性係数Ｅｚの２次近似値をＥｚ２＝１１［ＧＰａ］として、横弾性係数Ｇの近似値を決定し、ここで求められた新たな横弾性係数Ｇの近似値を用いて、縦弾性係数Ｅｚの近似値を求めることを繰り返し、その結果、前記試験用直方体の固有振動数算出値が固有振動数実測値に一致する縦弾性係数および横弾性係数として、Ｅｚ＝１０［ＧＰａ］およびＧ＝０．２９［ＧＰａ］をそれぞれ得た。
ところで、変圧器鉄心の固有振動解析を実施すると、図５に示すような振動モードが出現する。
ここで、従来例のように、変圧器鉄心の機械的弾性係数を、方向性電磁鋼板の機械的弾性係数と等しいとして図２に示す形状の三相三脚変圧器鉄心の固有振動解析を行った結果、固有振動数はモードＡが１００６Ｈｚ、モードＢが２００４Ｈｚとなり、モードＣの振動モードは出現しなかった。
ここで、従来例のように、変圧器鉄心の機械的弾性係数を、方向性電磁鋼板の機械的弾性係数と等しいとして図２に示す形状の三相三脚変圧器鉄心の固有振動解析を行った結果、固有振動数はモードＡが１００６Ｈｚ、モードＢが２００４Ｈｚとなり、モードＣの振動モードは出現しなかった。
次に、図２に示す形状の三相三脚変圧器の積層鉄心２１を板厚０．３ｍｍの方向性電磁鋼板を用いて実際に製造後、インパクトハンマーによる加振を行って鉄心の固有振動数測定を実施した。その結果、図５に示した３つの振動モードＡ、ＢおよびＣが観察され、計測された固有振動数は、モードＡが１４１Ｈｚ、モードＢが２５５Ｈｚ、モードＣが２９３Ｈｚであった。
このように、従来の固有振動解析で得られる固有振動数は、実測された固有振動数とは大きく異なる結果となる。これは、従来の固有振動解析が、積層鉄心の積層方向を含む二面の横弾性係数を積層された前記鋼板間の滑りを考慮していないことに起因すると考えられる。
このように、従来の固有振動解析で得られる固有振動数は、実測された固有振動数とは大きく異なる結果となる。これは、従来の固有振動解析が、積層鉄心の積層方向を含む二面の横弾性係数を積層された前記鋼板間の滑りを考慮していないことに起因すると考えられる。
そこで、三相三脚変圧器の横弾性係数を上記弾性マトリックス決定方法で決定した値であるＧ＝０．２９［ＧＰａ］とし、縦弾性係数をＥｚ＝１０［ＧＰａ］として、振動解析装置１０を使用して振動解析を実施した。その結果、固有振動数は、図５に示すモードＡが１６６Ｈｚ、モードＢが２０７Ｈｚ、モードＣが３５５Ｈｚとなった。これらの値は、従来の解析方法による結果と比較すると、はるかによく一致している。
以上の結果から、本実施形態による弾性マトリックス決定方法によって、弾性マトリックス［Ｃ］における横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを決定し、決定した横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを前述した（３）式の構成方程式に組み込んで、固有振動解析を実施することにより、振動特性を従来方法に比較して高精度で計算することができることが確認できた。
以上の結果から、本実施形態による弾性マトリックス決定方法によって、弾性マトリックス［Ｃ］における横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを決定し、決定した横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを前述した（３）式の構成方程式に組み込んで、固有振動解析を実施することにより、振動特性を従来方法に比較して高精度で計算することができることが確認できた。
このように、本実施形態によると、試験用鉄心を用いて横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇをパラメータとして複数回固有振動解析を実施することにより、横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを積層された方向性電磁鋼板間の滑りを考慮した値に設定することが可能となる。このため、決定した横弾性係数を弾性マトリックスに組み込んで振動解析を実施することにより、積層鉄心の振動特性を高精度で計算することができる。
したがって、種々の電磁鋼板を使用して三相三脚変圧器鉄心を設計した場合に、同一の電磁鋼板を使用した試験用鉄心を製作し、この試験用鉄心を使用して、縦弾性係数Ｅｚおよび横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを決定することにより、積層鋼板間の滑りを考慮した適正な弾性マトリックスすなわち応答関数を決定することができる。
したがって、種々の電磁鋼板を使用して三相三脚変圧器鉄心を設計した場合に、同一の電磁鋼板を使用した試験用鉄心を製作し、この試験用鉄心を使用して、縦弾性係数Ｅｚおよび横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを決定することにより、積層鋼板間の滑りを考慮した適正な弾性マトリックスすなわち応答関数を決定することができる。
また、決定した横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを弾性マトリックスに組み込んで固有振動解析を行うことにより、三相三脚変圧器鉄心の各振動モードでの固有振動数を高精度で算出することができる。
なお、上記実施形態では、試験用鉄心について、縦弾性係数Ｅｚと横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇの近似値を逐次交互に決めることを繰り返して、最終的なＥｚとＧの値を決定する場合について説明したが、これに限定されるものではなく、試験用鉄心について縦弾性係数Ｅｚおよび横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇの２つをパラメータとして固有振動数計算値が固有振動数実測値により一致するように固有振動解析を複数回実施して固有振動数に対する縦弾性係数Ｅｚおよび横弾性係数Ｇの三次元マップを作成し、この三次元マップから固有振動数実測値に相当する縦弾性係数Ｅｚおよび横弾性係数Ｇを決定する、等の方法でもよい。
なお、上記実施形態では、試験用鉄心について、縦弾性係数Ｅｚと横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇの近似値を逐次交互に決めることを繰り返して、最終的なＥｚとＧの値を決定する場合について説明したが、これに限定されるものではなく、試験用鉄心について縦弾性係数Ｅｚおよび横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇの２つをパラメータとして固有振動数計算値が固有振動数実測値により一致するように固有振動解析を複数回実施して固有振動数に対する縦弾性係数Ｅｚおよび横弾性係数Ｇの三次元マップを作成し、この三次元マップから固有振動数実測値に相当する縦弾性係数Ｅｚおよび横弾性係数Ｇを決定する、等の方法でもよい。
この結果、試験用鉄心を使用して、より固有振動数実測値に近い固有振動数計算値を得ることができる縦弾性係数Ｅｚおよび横弾性係数Ｇを決定することができる。
したがって、決定した縦弾性係数Ｅｚおよび横弾性係数Ｇを応力−歪み構成式に組み込んで三相三脚変圧器鉄心等の積層鉄心の固有振動解析を行うことにより、より高精度の固有振動解析を行うことができる。
また、上記実施形態では、試験用鉄心として、直方体鉄心を使用した場合について説明したが、これに限定されるものではなく、試験用三相三脚変圧器鉄心を製作し、この試験用三相三脚変圧器鉄心について固有振動数が測定値と一致する横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇ、又は縦弾性係数Ｅｚおよび横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを決定するようにしてもよい。
したがって、決定した縦弾性係数Ｅｚおよび横弾性係数Ｇを応力−歪み構成式に組み込んで三相三脚変圧器鉄心等の積層鉄心の固有振動解析を行うことにより、より高精度の固有振動解析を行うことができる。
また、上記実施形態では、試験用鉄心として、直方体鉄心を使用した場合について説明したが、これに限定されるものではなく、試験用三相三脚変圧器鉄心を製作し、この試験用三相三脚変圧器鉄心について固有振動数が測定値と一致する横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇ、又は縦弾性係数Ｅｚおよび横弾性係数Ｇｙｚ＝Ｇｚｘ＝Ｇを決定するようにしてもよい。
また、上記実施形態では、応答関数としてステフィネスマトリックスである弾性マトリックス［Ｃ］を使用する応力と歪みの関係式を使用して振動解析を行う場合について説明したが、これに限定されるものではなく、応答関数としてコンプライアンスマトリックス［Ｓ］を使用する歪みと応力の関係式｛ε｝＝［Ｓ］｛σ｝を使用して振動解析を行うようにしてもよい。この場合には、Ｓ４４＝１／Ｇｙｚ、Ｓ５５＝１／Ｇｚｘ、Ｓ６６＝１／Ｇｚｙとなる。
さらに、上記実施形態では、三相三脚変圧器についての振動解析について説明したが、これに限定されるものではなく三相五脚変圧器や他の変圧器における積層鉄心の振動解析にも本発明を適用することができる。
さらに、上記実施形態では、三相三脚変圧器についての振動解析について説明したが、これに限定されるものではなく三相五脚変圧器や他の変圧器における積層鉄心の振動解析にも本発明を適用することができる。
１０ 振動解析装置
１１ ＣＰＵ
１２ 演算処理装置
１３ 内部バス
１４ 内部記憶装置
１５ 外部記憶装置
１６ 入力装置
１７ 出力装置
１８ 記録媒体
２１ 積層鉄心
２２ 電磁鋼板
２２ａ 上ヨーク
２２ｂ 下ヨーク
２２ｃ 脚部
１１ ＣＰＵ
１２ 演算処理装置
１３ 内部バス
１４ 内部記憶装置
１５ 外部記憶装置
１６ 入力装置
１７ 出力装置
１８ 記録媒体
２１ 積層鉄心
２２ 電磁鋼板
２２ａ 上ヨーク
２２ｂ 下ヨーク
２２ｃ 脚部
Claims (5)
- 応力と歪みとの関係を行列表示で表した構成式を使用して鋼板を積層した積層鉄心の振動解析を行うに際し、前記構成式中の弾性マトリックスに含まれる、前記積層鉄心の積層方向を含む面の横弾性係数、及び、前記積層鉄心の積層方向の縦弾性係数を決定する弾性マトリックス決定方法であって、
振動解析の対象となる前記積層鉄心と同じ鋼板を積層してなる試験用鉄心を製作する第１ステップと、
前記第１ステップで製作した前記試験用鉄心を支持して積層方向に加振することにより曲げ振動の固有振動数を実測する第２ステップと、
前記試験用鉄心の構造解析モデルについて、積層方向の縦弾性係数として前記試験用鉄心に使用した鋼板１枚の縦弾性係数を近似値として与え、積層方向を含む面の横弾性係数をパラメータとして値を変えながら複数回構造解析を行って、前記積層方向を含む面の横弾性係数と固有振動数との関係を推定する第３ステップと、
前記第３ステップで推定した前記積層方向を含む面の横弾性係数と固有振動数との関係に、前記第２ステップで実測した前記試験用鉄心の固有振動数の実測値を与えて前記試験用鉄心の前記積層方向を含む面の横弾性係数の近似値を決定する第４ステップと、
前記試験用鉄心の構造解析モデルについて、前記積層方向を含む面の横弾性係数として前記第４ステップで決定した前記積層方向を含む面の横弾性係数の近似値を与え、前記積層方向の縦弾性係数をパラメータとして値を変えながら複数回構造解析を行って、前記積層方向の縦弾性係数と固有振動数との関係を推定する第５ステップと、
前記第５ステップで推定した前記積層方向の縦弾性係数と固有振動数との関係に、前記第２ステップで実測した前記試験用鉄心の固有振動数の実測値を与えて前記試験用鉄心の前記積層方向の縦弾性係数の近似値を決定する第６ステップと、
前記試験用鉄心の構造解析モデルについて、前記積層方向の縦弾性係数として前記第６ステップで決定した前記積層方向の縦弾性係数の近似値を与えて、前記第３ステップから前記第６ステップの処理を実行することを繰り返すことにより、前記積層鉄心の積層方向を含む面の横弾性係数、及び、前記積層鉄心の積層方向の縦弾性係数を決定する第７ステップと、
を含むことを特徴とする弾性マトリックス決定方法。 - 前記積層鉄心は、電磁鋼板を複数枚積層して構成されていることを特徴とする請求項１に記載の弾性マトリックス決定方法。
- 前記積層鉄心は、変圧器用積層鉄心であることを特徴とする請求項２に記載の弾性マトリックス決定方法。
- 応答関数を用いて応力と歪みとの関係をテンソル表示で表した構成式の弾性マトリックスに、請求項１に記載の弾性マトリックス決定方法で決定した前記積層鉄心の積層方向を含む面の横弾性係数、及び、前記積層鉄心の積層方向の縦弾性係数を組み込んで鋼板を積層した積層鉄心の振動解析を行うことを特徴とする積層鉄心の振動解析方法。
- 前記積層鉄心は、電磁鋼板を複数枚積層した変圧器用鉄心であることを特徴とする請求項４に記載の積層鉄心の振動解析方法。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018009662 | 2018-01-24 | ||
| JP2018009662 | 2018-01-24 | ||
| PCT/JP2019/001704 WO2019146554A1 (ja) | 2018-01-24 | 2019-01-21 | 積層鉄心の弾性マトリックス決定方法および振動解析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2019146554A1 JPWO2019146554A1 (ja) | 2020-02-06 |
| JP6645623B2 true JP6645623B2 (ja) | 2020-02-14 |
Family
ID=67394921
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019524096A Active JP6645623B2 (ja) | 2018-01-24 | 2019-01-21 | 積層鉄心の弾性マトリックス決定方法および振動解析方法 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US11513101B2 (ja) |
| EP (1) | EP3745286A4 (ja) |
| JP (1) | JP6645623B2 (ja) |
| CA (1) | CA3089090C (ja) |
| MX (1) | MX2020007739A (ja) |
| WO (1) | WO2019146554A1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023027444A1 (ko) * | 2021-08-27 | 2023-03-02 | 주식회사 엘지화학 | 다층 소재의 물성 예측 시스템 및 방법 |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7147831B2 (ja) * | 2020-02-21 | 2022-10-05 | Ｊｆｅスチール株式会社 | 変圧器の積層鉄心の弾性マトリックス決定方法および振動解析方法 |
| MX2024001685A (es) * | 2021-08-05 | 2024-03-04 | Jfe Steel Corp | Metodo de determinacion de matriz de elasticidad de nucleo de hierro apilado de transformador y metodo de analisis de vibracion. |
| CN114818148B (zh) * | 2022-06-24 | 2022-09-13 | 中国飞机强度研究所 | 飞机关键构件在振动试验中的等效设计方法 |
Family Cites Families (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3943755A (en) * | 1974-10-22 | 1976-03-16 | Tokyo Shibaura Electric Co., Ltd. | Method and apparatus for measuring the magnitude of a clamping load applied to a laminated iron core of an electric machine |
| US4499770A (en) * | 1982-07-22 | 1985-02-19 | The United States Of America As Represented By The Secretary Of Commerce | Systems for monitoring changes in elastic stiffness in composite materials |
| CN85102855B (zh) * | 1985-04-01 | 1987-09-23 | 曹培生 | 一种非线性电磁振动装置 |
| GB9201880D0 (en) * | 1992-01-29 | 1992-03-18 | Rasmussen O B | Laminated films |
| US7027048B2 (en) * | 2002-05-31 | 2006-04-11 | Ugs Corp. | Computerized deformation analyzer |
| JP5309431B2 (ja) * | 2006-08-04 | 2013-10-09 | 新日鐵住金株式会社 | 鋼板剪断面の鋼板間抵抗が高い電磁鋼の積層鋼板およびそのカシメ方法 |
| JP5003609B2 (ja) * | 2008-06-20 | 2012-08-15 | 株式会社Ｉｈｉ | 異方性材の等価板厚決定方法 |
| JP5953686B2 (ja) * | 2011-09-22 | 2016-07-20 | Ｊｆｅスチール株式会社 | 変圧器の騒音予測方法 |
| JP6009305B2 (ja) | 2012-09-28 | 2016-10-19 | 株式会社Ｊｓｏｌ | 解析装置、解析方法及びコンピュータプログラム |
| WO2015021182A1 (en) * | 2013-08-06 | 2015-02-12 | Bp Corporation North America Inc. | Image-based direct numerical simulation of petrophysical properties under simulated stress and strain conditions |
| JP2015032295A (ja) * | 2013-08-07 | 2015-02-16 | ヤマハ株式会社 | 減衰振動解析方法 |
| CN106461614B (zh) * | 2014-05-28 | 2019-04-26 | 汉阳大学校产学协力团 | 利用超声波的热化评价及强度估算装置及方法 |
| WO2016027874A1 (ja) * | 2014-08-21 | 2016-02-25 | 公立大学法人大阪市立大学 | 応力可視化装置および力学物性値可視化装置 |
| CN104359549B (zh) * | 2014-11-27 | 2015-07-08 | 国家电网公司 | 一种变压器铁芯振动噪声分析方法 |
| US11965858B2 (en) * | 2018-11-06 | 2024-04-23 | Jfe Steel Corporation | Laminated iron core elastic matrix determination method and vibration analysis method |
-
2019
- 2019-01-21 US US16/960,931 patent/US11513101B2/en active Active
- 2019-01-21 CA CA3089090A patent/CA3089090C/en active Active
- 2019-01-21 JP JP2019524096A patent/JP6645623B2/ja active Active
- 2019-01-21 MX MX2020007739A patent/MX2020007739A/es unknown
- 2019-01-21 WO PCT/JP2019/001704 patent/WO2019146554A1/ja unknown
- 2019-01-21 EP EP19743839.3A patent/EP3745286A4/en active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023027444A1 (ko) * | 2021-08-27 | 2023-03-02 | 주식회사 엘지화학 | 다층 소재의 물성 예측 시스템 및 방법 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPWO2019146554A1 (ja) | 2020-02-06 |
| MX2020007739A (es) | 2020-09-25 |
| EP3745286A1 (en) | 2020-12-02 |
| US20200340952A1 (en) | 2020-10-29 |
| CA3089090C (en) | 2022-05-24 |
| CA3089090A1 (en) | 2019-08-01 |
| US11513101B2 (en) | 2022-11-29 |
| WO2019146554A1 (ja) | 2019-08-01 |
| EP3745286A4 (en) | 2021-03-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6645623B2 (ja) | 積層鉄心の弾性マトリックス決定方法および振動解析方法 | |
| Kitagawa et al. | Analysis of structural deformation and vibration of a transformer core by using magnetic property of magnetostriction | |
| JP5556882B2 (ja) | 磁性体のシミュレーション方法及びプログラム | |
| JP6519810B2 (ja) | 変圧器内部異常および劣化の診断方法と診断装置 | |
| JP6009305B2 (ja) | 解析装置、解析方法及びコンピュータプログラム | |
| JP7147831B2 (ja) | 変圧器の積層鉄心の弾性マトリックス決定方法および振動解析方法 | |
| JP6729837B1 (ja) | 積層鉄心の弾性マトリックス決定方法および振動解析方法 | |
| Case | Numerical analysis of the vibration and acoustic characteristics of large power transformers | |
| JP7218839B1 (ja) | 変圧器の積層鉄心の弾性マトリックス決定方法及び振動解析方法 | |
| JP7207623B1 (ja) | 変圧器の騒音性能診断方法及び低騒音化方法 | |
| Masti et al. | On the influence of core laminations upon power transformer noise | |
| Monturet et al. | Optimal dimensioning of a piezoelectric bimorph actuator | |
| JP2023144273A (ja) | 積層鉄心の弾性マトリックス決定装置、積層鉄心の弾性マトリックス決定方法およびコンピュータプログラム | |
| Rens et al. | Numerical and experimental evaluation of magnetostriction and magnetic forces on transformer stacks and joints for the assessment of core vibrations | |
| Haettel et al. | Prediction of transformer core noise | |
| Kaltenbacher et al. | Anisotropic model for the numerical computation of magnetostriction in grain-oriented electrical steel sheets | |
| Dobrucki et al. | Generalized Oberst beam method for measuring viscoelastic parameters of layered composite components | |
| Lahyaoui et al. | Macroscopic model of magnetostriction based on energy minimization | |
| CN115526079A (zh) | 分步等效的变压器振动噪声仿真方法 | |
| Gao et al. | Homogenization Technique of Magnetostriction Calculation in Laminated Iron Core | |
| Villalva et al. | Study of Correlation between Vibration Tests and Finite Element Simulations by Means of Impulse Excitation Technique Applied to Components of an Automatized Transmission System | |
| Gao | Improved Methods for Measuring Characteristic Functions of Underwater Transducers | |
| Butler et al. | Transducer models | |
| MURGIDA | Vibroacoustic methodology for transformers air-borne noise computation based on FEM multiphysics approach |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190508 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20190508 |
|
| A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20190807 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190813 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20191004 |
|
| 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: 20191210 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20191223 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6645623 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |