静磁場解析
まず、静磁場解析から始めます。以下の多くの機能は交流定常解析や過渡解析でも共通して使用できます。
ELMCUR（要素電流ソース）を用いた静磁場解析
EMSolutionで使われるELMCUR（要素電流ソース）は、要素それぞれに電流密度あるいは要素面を通過する電流を与え、電流ソースとするものです。
要素の面や方向がそろっている場合に使用することができます。その例として、電気学会三次元静磁界検証モデル（電気学会技術報告（II部）第286号、三次元静磁界数値計算技術、昭和63年12月）の解析を行います。
本モデルでは、100mm×100mm×200mm（高さ）の鉄芯に四角形コイルによる磁場を加え、その磁場分布を求めるものです。
コイルは高さ100mm,内幅150mm,外幅200mmで、R25（内側）,R50（外側）の角取りがあります。鉄芯は比透磁率1000とし、コイルには3000ATの電流が流れるとします。対称性より1/8領域（x,y,z≧0）をモデル化します。
まず、Fig.1のように二次元メッシュを定義します。FEMAP入力の場合、物性番号としてFEMAPのProperty IDを用います。Fig.1の様に各領域にproperty IDを定義して下さい。
このメッシュデータをファイルpre_geom2D.neuとしてEMSolutionに入力します。ｚ方向への拡張は、ファイル2D_to_3Dを入力することにより行います。上の二つのファイルと実行制御ファイルinputをを同じディレクトリに置き、EMSolutionを実行します。
EMSolution実行画面で、実行メニューを選択しますとファイルダイアログが現れますので、input.1ファイルを選択して実行を開始してください。
EMSolutionは、pre_geom2D.neuと2D_to_3Dをもとに、三次元メッシュ（Fig.2）を作製します。そのメッシュデータは、post_geomとして出力されます。もちろん、三次元メッシュを直接入力して計算を始めることもできます。ELMCURでソース電流を与えるためには、要素の向きをそろえておく必要があります。今の場合、電流は面3から入り面5から出る様に六面体要素を定義します。
これは、FEMAPにおいて、各Surfaceを定義する時に、第一番目のEdgeを電流が流入する面に来るようにすれば可能です。
 |  |
| Fig.1 モデル二次元メッシュ | Fig.2 三次元メッシュ |
出力ファイルoutput（output.1）が出力されます。outputには、入力メッシュのサマリー、計算の過程、計算結果等が出力されます。また、ファイルcheckには入力のエコー、ICCG法の詳細収束過程が出力されます。実行の異常が見られるときは、まずこの両ファイルを確認下さい。
また、stderrファイルにエラーメッセージが出力されることがあります。実行ディレクトリには各種のファイルが出力されますが、リスタートを予定している場合にはファイルはそのまま残しておいて下さい。実行後、ポスト出力を行います。FEMAPを実行できる方は、まず、ファイル名post_geom, magnetic, currentに.neuを付け、FEMAPで読みます。
post_geomにはProperty, Materialの定義は入っていませんので、必要ならFEMAP上で定義して下さい。Fig.3に磁束密度分布出力図、Fig.4にコイル内の電流分布出力図を示します。磁束密度はT単位、電流密度はA/m2単位で出力れます。電流密度はコイル内でほぼ一定になっています。
 |  |
| Fig.3 磁束密度分布（Ｔ） | Fig.4 コイル内電流分布（A/m2） |
次に、リスタートにより電磁力を出力します。このときポスト処理プロセスのみを行えば済みます。再実行時には、ファイルoutput等出力ファイルが上書きされますので、必要な場合はファイル名の変更等を行い保存してください。まず、節点力分布を求めるため、ファイルinputの一部を変更し（input.2）、同じディレクトリでEMSolutionを再実行します。
ファイルoutputの電磁力出力部分をList.1に示します。モデル部分1/8領域に働く鉄芯部とコイル部に働く各方向のトータル電磁力です。電磁力分布はファイルforceに出力されます。Fig.5に鉄芯部に加わる節点電磁力を示します。各節点に働く電磁力（N）で示しています。透磁率一定の部材には表面しか電磁力が現れません。
同様に、コイル部に対してローレンツ力を求めます（input.3）。ファイルoutputの、各物性毎の電磁力（List.1）をみると、上の節点力求めたコイルに働く電磁力とここで求めたものがほぼ一致します。この場合、電磁力分布は、ファイルforce_J_Bに出力されます。Fig.6にコイル部のローレンツ力分布（N/m3）を示します。
 |  |
| Fig.5 鉄芯節点力分布（N） | Fig.6 コイルローレンツ力分布（N/m3） |
List.1
| 節点力により求めた電磁力 *** Force to nodes （N） and torque（Nm） ********************************************** |
| No. | Fx（N） | Fy（N） | Fz（N） | Mx（Nm） | My（Nm） | Mz（Nm） |
| MAT No. | 1 | 1.5726e+000 | 1.5806e+000 | 1.0975e+000 | -8.0615e-002 | 8.0065e-002 | 2.9609e-004 |
| MAT No. | 3 | -2.8444e-001 | -2.8693e-001 | -2.1822e+000 | -1.2893e-001 | 1.2901e-001 | -1.1002e-004 |
| Total | | 1.4487e+000 | 1.4488e+000 | -1.0733e+000 | -2.1717e-001 | 2.1717e-001 | -5.6200e-006 |
| ローレンツ力より求めたコイル電磁力 *** Force J*B （N） （Volume total） ********************************************** |
| No. | Fx | Fy | Fz |
| MAT No. | 3 | -2.6849e-001 | -2.6853e-001 | -2.2211e+000 |
| Total（N） | | -2.6849e-001 | -2.6853e-001 | -2.2211e+000 |
使用データ：pre_geom2D.neu,2D_to_3D,input.1,output.1, input.2,input.3
静磁場解析