非線形定常熱伝導解析

材料特性が温度の関数であるシステムにおける温度分布を計算します。

(1)
( K C +H )T+R ( T+ T abs ) 4 =f MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaadaqadaqaaiaahUeadaWgaaWcbaGaam 4qaaqabaGccqGHRaWkcaWHibaacaGLOaGaayzkaaGaaCivaiabgUca RiaahkfadaqadaqaaiaahsfacaWHRaGaaCivamaaBaaaleaacaWHHb GaaCOyaiaahohaaeqaaaGccaGLOaGaayzkaaWaaWbaaSqabeaacaaI 0aaaaOGaeyypa0JaaCOzaaaa@434B@
ここで、
K C MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHlbWaaSbaaSqaaiaadoeaaeqaaa aa@3383@
熱伝導マトリックス
H
自由伝達による境界伝達マトリックス
T
未知の節点温度
R MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWGsbaaaa@3292@
放射交換マトリックス
T abs MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHubWaaSbaaSqaaiaahggacaWHIb GaaC4Caaqabaaaaa@3595@
次で定義した、絶対温度のスケール: PARAM, TABS
f
温度荷重ベクトル

方程式の系は節点温度 T を見つけるために解かれます。

温度荷重ベクトルは次のように表されます:(2)
f= f B + f H + f Q + f R MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbGaeyypa0JaaCOzamaaBaaale aacaWGcbaabeaakiabgUcaRiaahAgadaWgaaWcbaGaamisaaqabaGc cqGHRaWkcaWHMbWaaSbaaSqaaiaadgfaaeqaaOGaey4kaSIaaCOzam aaBaaaleaacaWGsbaabeaaaaa@3E21@
ここで、
f B MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbWaaSbaaSqaaiaadkeaaeqaaa aa@339D@
QBDY1で指定される境界での熱流束によるパワー
f H MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbWaaSbaaSqaaiaadkeaaeqaaa aa@339D@
CONVで指定される伝達による境界伝達ベクトル(CONVGバルク/サブケースの組み合わせで、自由対流の自動定義を有効にすることができます。)
f Q MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbWaaSbaaSqaaiaadkeaaeqaaa aa@339D@
QVOLで指定される内部熱生成によるパワーのベクトル
f R MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbWaaSbaaSqaaiaadkeaaeqaaa aa@339D@
RADBCで指定される放射による境界放射ベクトル

式 1の左側のマトリックスは、温度境界条件が指定されない限り特異となります。

方程式のシステム(上記)はニュートン法を使用して解かれます。解のコントロールは、NLPARMバルクデータエントリでパラメータを定義することで与えられます。TEMPERATURE(INITIAL)を使用することにより、初期温度分布の設定が可能です。非線形熱伝導解析による温度結果は後続の構造解析で使用できます。

非線形定常熱伝導解析の設定

次の手順は、非線形定常熱伝導解析を設定するためのガイドラインです:
  1. サブケース情報セクションのソリューションシーケンスの識別子(ANALYSIS)を使用し、ANALYSIS=NLHEATを用いて非線形定常熱伝導解析を選択します。
  2. 初期温度分布はTEMPERATUREサブケース情報エントリ(type=INITIAL)を使用して定義できます。初期温度分布の予測が優れていれば、ソルバーの収束が改善されます。
  3. 温度依存の熱材料特性は、MATT4バルクデータエントリを使って定義することが可能です。
  4. サブケースに対して非線形解法が必要であることを示すには、NLPARMサブケース情報エントリが必要です。このサブケースエントリは、収束トレランスおよびその他の非線形パラメーターを指定するNLPARMバルクデータエントリを指します。
  5. 荷重条件および境界条件は、入力デックのバルクデータセクションで定義します。これらをサブケース情報エントリセクション内で参照する必要があり、これにはサブケースでSPCおよびLOADエントリを使用します。各サブケースは、1つの荷重ベクトルを定義します。

放射の入力

  1. 放射(Radiation)は1つの境界条件と見なされます(伝達(Convection)と同様)。
  2. 放射交換は、サーフェス要素と黒体空間節点との間で為されます。
  3. 環境空間節点はSPCバルクデータエントリによって定義され、RADBCバルクデータエントリのNODAMBフィールドを介して認識されます。
  4. RADBCはサーフェス要素に適用される際にのみ使用されます(CHBDYEバルクデータ)。
  5. 環境節点は、サーフェス要素と空間との間の放射交換について自身の温度を擁する黒体として扱われます。
  6. 放射ビュー係数( RADBCバルクデータエントリのFAMB)は、サーフェス要素と環境節点との間の幾何学的関係を定義します。
  7. 輻射率(EMIS1フィールド)と吸収係数(ABSORPフィールド)は、RADMバルクデータエントリ上で定義されることが可能です。
  8. RADMエントリは、サーフェス要素(CHBDYEエントリ)によって直接参照されます。
  9. 放射の計算には、2つのパラメータが必要です:
    • PARAM, TABS: 壁温度スケールファクターの絶対値
    • PARAM, SIGMA: シュテファンボルツマン定数

例:ANALYSISNLPARMの使用方法

SUBCASE  5
ANALYSIS=NLHEAT
SPC=10
LOAD=20
NLPARM=30
…
BEGIN BULK
…
NLPARM, 30
…
ENDDATA
注:
  1. 非線形熱伝導解析に基づく最適化は、現時点ではサポートされていません。
  2. CNTNLSUBは現時点では、非線形熱伝導解析用としてサポートされていません。
  3. シェル要素は、熱伝導解析においては膜として扱われます。複合材プロパティ(PCOMP)は、均質化によってNLHEATでサポートされています。シェル要素の板厚を通した温度分布は計算されません。節点温度のみが決定されます(節点ごとに1つの自由度)
  4. 複合プレートまたはシェル要素の流束は、要素全体の均質化された伝導係数で計算されます。