概要
過渡解析において、発熱の時間平均を出力できるようになりました。ここでは、交流解析と同じ場合を過渡解析で解き、発熱を比較します。
解説
ここでは線形解析の場合を示しますが、非線形過渡解析で周期性ののある場合の時間平均を取る際も有効です。
まず、交流定常解析を行います。縦横0.1mで1cm厚の銅板に50Hzの磁場が板に垂直に加わるとします(input.AC)。これを実行しますと、outputに次のように一周期平均発熱量が出力されます。この量は解析領域(全体の1/8)に対する発熱量です。
また、発熱分布はheatファイルに出力され、Fig.1の様になります。
***********************************************************************************
*** Total Joule heatin loss ( 1 cycle average ) in regions ***
***********************************************************************************
MAT No. Q (W)
1 1.3973e+04
Total 1.3973e+04
次に過渡解析で同じ場合を計算してみます。一周期を40分割し、二周期分計算します。半周期毎の平均発熱を出力します(input.transient)。
outputには以下のように出力されます。
二周期目でほぼ定常に達し、交流定常解析ともほぼ一致します。
また、最後の半周期での平均発熱分布をFig.2に示します。
交流定常解析と良く一致していることが解ります
Fig.1 一周期平均発熱密度分布
$(W/m^3)$
Fig.2 過渡解析による
半周期平均発熱密度分布$(W/m^3)$
*************************************************************************************
* Average from Step No. 1 to Step No. 20. Period 2.0000000000e-02 sec *
*************************************************************************************
*** Total Joule heatin loss in regions **********************
MAT No. Q (W)
1 1.4769e+04
Total 1.4769e+04
*************************************************************************************
* Average from Step No. 21 to Step No. 40. Period 2.0000000000e-02 sec *
*************************************************************************************
*** Total Joule heatin loss in regions **********************
MAT No. Q (W)
1 1.3602e+04
Total 1.3602e+04
*************************************************************************************
* Average from Step No. 41 to Step No. 60. Period 2.0000000000e-02 sec *
*************************************************************************************
*** Total Joule heatin loss in regions **********************
MAT No. Q (W)
1 1.3519e+04
Total 1.3519e+04
*************************************************************************************
* Average from Step No. 61 to Step No. 80. Period 2.0000000000e-02 sec *
*************************************************************************************
*** Total Joule heatin loss in regions **********************
MAT No. Q (W)
1 1.3530e+04
Total 1.3530e+04
************************************************************************************