"熱伝導もしくは熱応力の数値解析システムであって、解析領域を少なくとも2つ以上の領域に分け、少なくとも一つの領域を、有限要素法もしくは境界要素法にて解析し、他の少なくとも一つの領域を、等価回路近似に基づいた手法にて解析することを特徴とする数値解析システム。請求項1において、前記等価回路近似に基づいた手法にて解析する方法に対して、熱伝導解析の場合は、前記解析領域の熱抵抗Rもしくは熱抵抗Rの逆数であるアドミッタンスY=R1を予め求めておき、さらに、前記解析領域の境界における熱量Qから、熱等価回路方程式ΔT=RQに基づいて、前記解析領域の境界における温度変化ΔTを導出することを特徴とする数値解析システム。請求項2において、前記解析領域の熱抵抗Rもしくは熱抵抗の逆数アドミッタンスY=R1を予め求める方法に対して、前記解析領域の境界において、熱量が流出入する箇所もしくは温度変化を観測したい箇所が2つ以上ある場合は、前記解析領域の熱抵抗行列[R]もしくは熱抵抗行列[R]の逆行列であるアドミッタンス行列[Y]=[R]1を予め求めておき、さらに、前記解析領域の境界における熱量列[Q]から、熱等価回路方程式[ΔT]=[R][Q]に基づいて、前記解析領域の境界における温度変化列[ΔT]を導出することを特徴とする数値解析システム。請求項1において、前記等価回路近似に基づいた手法にて解析する方法に対して、熱応力解析の場合は、請求項2または請求項3に記載した熱等価回路方程式にて、前記解析領域の境界における温度変化ΔTもしくは温度変化列[ΔT]を導出した後、その温度変化を用いて、更に熱膨張あるいは熱収縮により発生する応力を導出することを特徴とする数値解析システム。請求項1において、解析領域を少なくとも2つ以上の領域に分割する方法に対して、形状を変更しない領域と形状を変更し得る領域とに分割し、形状を変更しない領域に対しては、等価回路近似に基づいた手法の解析対象とし、形状を変更し得る領域に対しては、有限要素法もしくは境界要素法の解析対象とすることを特徴とする数値解析システム。請求項1において、解析領域を少なくとも2つ以上の領域に分割する方法に対して、熱源が一つ以上存在する領域と熱源が存在しない領域とに分割し、熱源が存在しない領域に対しては、等価回路近似に基づいた手法の解析対象とし、熱源が一つ以上存在する領域に対しては、有限要素法もしくは境界要素法の解析対象とすることを特徴とする数値解析システム。請求項2または請求項3において、前記解析領域で求めた熱抵抗値Rまたは熱抵抗行列[R]、或いは、熱抵抗Rの逆数であるアドミッタンスYまたは熱抵抗行列[R]の逆行列であるアドミッタンス行列[Y]に対して、それらの値を一度導出した後は、それをデータベースDB化して活用することを特徴とする数値解析システム。請求項1において、前記等価回路近似に基づいた手法にて解析する方法に対して、等価回路近似に基づいた手法にて解析実行の対象とする解析領域に対して、更に細かく、複数に細分した解析領域に分割し、各細分した解析領域に対して、請求項2または請求項3に記載した解析領域で求めた熱抵抗値Rまたは熱抵抗行列[R]、或いは、熱抵抗Rの逆数であるアドミッタンスYまたは熱抵抗行列[R]の逆行列であるアドミッタンス行列[Y]を求めて、それらの値や行列を合成することで、等価回路近似に基づいた手法にて解析実行する解析領域の熱抵抗値Rまたは熱抵抗行列[R]、或いは、熱抵抗Rの逆数であるアドミッタンスYまたは熱抵抗行列[R]の逆行列であるアドミッタンス行列[Y]を求めることを特徴とする数値解析システム。請求項8において、等価回路近似に基づいた手法にて解析実行の対象とする解析領域を更に細分した解析領域に分割する方法に対して、細分した解析領域に分割する指標乃至は単位として、部品乃至は周辺をも含めた部品を一つの細分した解析領域の単位として分割する、または、形状を変更しない領域と形状を変更し得る領域との細分した解析領域に分割することを特徴とする数値解析システム。請求項8において、等価回路近似に基づいた手法にて解析実行の対象とする解析領域を更に細分領域に分割する方法に対して、細分領域に分割する指標として、等価回路近似に基づいた手法にて解析実行の対象とする解析領域に対して、主な熱の流入部から主な熱の流出部までの主な熱の流れの方向に、平行に近い細分領域分割線の数に比べ、主な熱の流れの方向に、垂直に近い細分領域分割線の数が等しい乃至はそれ以上となることを特徴とする数値解析システム。請求項1において、解析領域を少なくとも2つ以上の領域に分け、少なくとも一つの領域を有限要素法もしくは境界要素法にて解析し、他の少なくとも一つの領域を等価回路近似に基づいた手法にて解析する方法に対して、一方の領域で得られた解析結果を、他の領域と接する境界にて、他の領域における次の解析の境界値として受け渡し、物理量の保存及び整合性が合うように連成解析することを特徴とする数値解析システム。請求項1において、解析領域を少なくとも2つ以上の領域に分け、少なくとも一つの領域を有限要素法もしくは境界要素法にて解析し、他の少なくとも一つの領域を等価回路近似に基づいた手法にて解析する方法に対して、まず、有限要素法FEMもしくは境界要素法BEMで解析する領域に対して、等価回路近似に基づいた手法にて解析する対象の領域との境界における温度変化をΔT=0として初期設定して有限要素法もしくは境界要素法にて対象領域の内部及び境界の熱流束q分布、温度変化ΔT分布を求め、それに基づき、等価回路近似に基づいた手法の対象領域との境界における熱流束分布から、前記境界の熱量QBを計算し、次に、前記境界の熱量QBを、等価回路近似に基づいた手法の対象領域の熱等価回路方程式[ΔT]=[R][Q]に基づいて、解析領域の境界における温度変化ΔTBを導出し、次に、境界における温度変化の初期設定としてΔT=0として有限要素法もしくは境界要素法で解析して得られているΔT分布に対して、解析領域の境界における温度変化ΔTB分の値を加算して、温度変化ΔT分布を導出することを特徴とする数値解析システム。請求項1において、解析領域を少なくとも2つ以上の領域に分け、少なくとも一つの領域を有限要素法もしくは境界要素法にて解析し、他の少なくとも一つの領域を等価回路近似に基づいた手法にて解析する方法に対して、有限要素法もしくは境界要素法で解析する領域と等価回路近似に基づいた手法で解析する領域との境界における熱量Qを、全体系に存在する全ての熱源から予測した熱量Qを設定し、それにより、まず、等価回路近似に基づいた手法の対象領域の熱等価回路方程式ΔT=RQに基づいて、解析領域の境界における温度変化ΔTを導出し、次に、境界の温度変化ΔTを境界初期条件として、有限要素法もしくは境界要素法にて対象領域の内部及び境界の熱流束q分布、温度変化ΔT分布を導出することを特徴とする数値解析システム。請求項11において、一方の領域で得られた解析結果を他の領域と接する境界にて、他の領域における次の解析の境界値として受け渡し、物理量の保存及び整合性が合うように連成解析する方法に対して、有限要素法もしくは境界要素法の解析、等価回路近似に基づいた解析の連成解析を、少なくとも2回以上繰り返し反復計算することを特徴とする数値解析システム。請求項14において、有限要素法もしくは境界要素法の解析、等価回路近似に基づいた解析の連成解析を、少なくとも2回以上繰り返し反復計算する方法に対して、各領域の各手法で求めた、境界における最新の温度変化ΔTの残差を求め、その残差に緩和係数ω≦1を掛けた値を、前回の解析あるいは前々回の解析結果の値に加算した値を次の解析のための境界値として受け渡すことを特徴とする数値解析システム。請求項15において、各領域の各手法で求めた、境界における最新の温度変化ΔTの残差を求め、その残差に緩和係数ω≦1を掛けた値を、前回の解析あるいは前々回の解析結果の値に加算した値を次の解析のための境界値として受け渡す方法に対して、前記緩和係数ωを前記連成解析の反復計算の途中で、値を変動させることを特徴とする数値解析システム。請求項16において、前記緩和係数ωを前記連成解析の反復計算の途中で、値を変動させる方法に対して、反復計算回数が少ない時は小さい値ω≦05に設定し、反復計算回数が多くなった際に、緩和係数ωとして大きな値<05<ω≦10になるように前記緩和係数ωを反復計算の途中で、値を変動させることを特徴とする数値解析システム。請求項1において、解析計算する手段として、PCクラスターまたはマルチコアPC乃至はマルチスレッドPCを用いる場合は、PC毎またはコア乃至はスレッド毎に各分割領域、各部分領域毎の計算を割り当て実施し、全体計算を高速化することを特徴とする数値解析システム。請求項1において、熱伝導もしくは熱応力解析に対して、前記解析対象の物理量以外でも、熱解析と同様に、スカラーポテンシャルにより場の記述が可能な物理量の解析を特徴とする数値解析システム。請求項1において、解析領域を少なくとも2つ以上の領域に分ける方法、乃至は、請求項8に記載の等価回路近似に基づいた解析の対象とする解析領域に対して、更に細かく、複数に細分した解析領域に分割する方法に対して、前記領域分割、前記細分領域分割に関して、解析を実行するユーザが入力、乃至は、設定するユーザインターフェース機能を備えることを特徴とする数値解析システム。コンポーネント毎に領域を切り出して、コンポーネント毎に熱評価する為、他の領域の影響が考慮されず、全体として熱評価する際の精度が不足する。また、パワエレシステム全体系をFEMで熱解析するためには、非常に大規模な計算コスト時間やメモリを必要とし、現実的ではない。上記の課題を解決するためには、パワエレシステム全体系などの大規模な対象に対して、現実的な計算コストの範囲内で、高精度な熱伝導解析システムを提供することが必要である。本発明では、前記の課題を解決するため、主に以下の構成に基づいた数値解析システムである。1解析領域を少なくとも2つ以上の領域に分け、2少なくとも一つの領域を、有限要素法もしくは境界要素法で解析し、3他の少なくとも一つの等価回路近似に基づいた解析にて解析する数値解析システム。本発明の熱伝導もしくは熱応力の数値解析システムは、特に、等価回路近似に基づいた手法にて解析する方法に対して、熱伝導解析の場合は、前記解析領域の熱抵抗Rもしくは熱抵抗Rの逆数であるアドミッタンスY=R1を予め求めておき、さらに、前記解析領域の境界における熱量Qから、熱等価回路方程式ΔT=RQに基づいて、前記解析領域の境界における温度変化ΔTを導出することが好ましい。また、解析領域の熱抵抗Rもしくは熱抵抗の逆数アドミッタンスY=R1を予め求める方法に対して、前記解析領域の境界において、熱量が流出入する箇所もしくは温度変化を観測したい箇所が2つ以上ある場合は、前記解析領域の熱抵抗行列[R]もしくは熱抵抗行列[R]の逆行列であるアドミッタンス行列[Y]=[R]1を予め求めておき、さらに、前記解析領域の境界における熱量列[Q]から、熱等価回路方程式[ΔT]=[R][Q]に基づいて、前記解析領域の境界における温度変化列[ΔT]を導出するが好ましい。また、前記等価回路近似に基づいた手法にて解析する方法に対して、熱応力解析の場合は、上述した熱等価回路方程式にて、前記解析領域の境界における温度変化ΔTもしくは温度変化列[ΔT]を導出した後、その温度変化を用いて、更に熱膨張あるいは熱収縮により発生する応力を導出することが好ましい。また、解析領域を少なくとも2つ以上の領域に分割する方法に対して、形状を変更しない領域と形状を変更し得る領域とに分割し、形状を変更しない領域に対しては、等価回路近似に基づいた手法の解析対象とし、形状を変更し得る領域に対しては、有限要素法もしくは境界要素法の解析対象とすることが好ましい。また、解析領域を少なくとも2つ以上の領域に分割する方法に対して、熱源が一つ以上存在する領域と熱源が存在しない領域とに分割し、熱源が存在しない領域に対しては、等価回路近似に基づいた手法の解析対象とし、熱源が一つ以上存在する領域に対しては、有限要素法もしくは境界要素法の解析対象とすることが好ましい。また、前記解析領域で求めた熱抵抗値Rまたは熱抵抗行列[R]、或いは、熱抵抗Rの逆数であるアドミッタンスYまたは熱抵抗行列[R]の逆行列であるアドミッタンス行列[Y]に対して、それらの値を一度導出した後は、それをデータベースDB化して活用することが好ましい。また、前記等価回路近似に基づいた手法にて解析する方法に対して、等価回路近似に基づいた手法にて解析実行の対象とする解析領域に対して、更に細かく、複数に細分した解析領域細分領域に分割し、各細分した解析領域に対して、前述の解析領域で求めた熱抵抗値Rまたは熱抵抗行列[R]、或いは、熱抵抗Rの逆数であるアドミッタンスYまたは熱抵抗行列[R]の逆行列であるアドミッタンス行列[Y]を求めて、それらの値や行列を合成することで、等価回路近似に基づいた手法にて解析実行する解析領域の熱抵抗値Rまたは熱抵抗行列[R]、或いは、熱抵抗Rの逆数であるアドミッタンスYまたは熱抵抗行列[R]の逆行列であるアドミッタンス行列[Y]を求めることが好ましい。また、等価回路近似に基づいた手法にて解析実行の対象とする解析領域を更に細分した解析領域に分割する方法に対して、細分した解析領域に分割する指標乃至は単位として、部品乃至は周辺をも含めた部品を一つの細分した解析領域の単位として分割する、または、形状を変更しない領域と形状を変更し得る領域との細分した解析領域に分割することが好ましい。また、等価回路近似に基づいた手法にて解析実行の対象とする解析領域を更に細分領域に分割する方法に対して、細分領域に分割する指標として、等価回路近似に基づいた手法にて解析実行の対象とする解析領域に対して、主な熱の流入部から主な熱の流出部までの主な熱の流れの方向に、平行に近い細分領域分割線の数に比べ、主な熱の流れの方向に、垂直に近い細分領域分割線の数が等しい乃至はそれ以上となることが好ましい。また、解析領域を少なくとも2つ以上の領域に分け、少なくとも一つの領域を有限要素法もしくは境界要素法にて解析し、他の少なくとも一つの領域を等価回路近似に基づいた手法にて解析する方法に対して、一方の領域で得られた解析結果を、他の領域と接する境界にて、他の領域における次の解析の境界値として受け渡し、物理量の保存及び整合性が合うように連成解析することが好ましい。また、解析領域を少なくとも2つ以上の領域に分け、少なくとも一つの領域を有限要素法もしくは境界要素法にて解析し、他の少なくとも一つの領域を等価回路近似に基づいた手法にて解析する方法に対して、まず、有限要素法FEMもしくは境界要素法BEMで解析する領域に対して、等価回路近似に基づいた手法にて解析する対象の領域との境界における温度変化をΔT=0として初期設定して有限要素法もしくは境界要素法にて対象領域の内部及び境界の熱流束q分布、温度変化ΔT分布を求め、それに基づき、等価回路近似に基づいた手法の対象領域との境界における熱流束分布から、前記境界の熱量QBを計算し、次に、前記境界の熱量QBを、等価回路近似に基づいた手法の対象領域の熱等価回路方程式[ΔT]=[R][Q]に基づいて、解析領域の境界における温度変化ΔTBを導出し、次に、境界における温度変化の初期設定としてΔT=0として有限要素法もしくは境界要素法で解析して得られているΔT分布に対して、解析領域の境界における温度変化ΔTB分の値を加算して、温度変化ΔT分布を導出することが好ましい。また、解析領域を少なくとも2つ以上の領域に分け、少なくとも一つの領域を有限要素法もしくは境界要素法にて解析し、他の少なくとも一つの領域を等価回路近似に基づいた手法にて解析する方法に対して、有限要素法もしくは境界要素法で解析する領域と等価回路近似に基づいた手法で解析する領域との境界における熱量Qを、全体系に存在する全ての熱源から予測した熱量Qを設定し、それにより、まず、等価回路近似に基づいた手法の対象領域の熱等価回路方程式ΔT=RQに基づいて、解析領域の境界における温度変化ΔTを導出し、次に、境界の温度変化ΔTを境界初期条件として、有限要素法もしくは境界要素法にて対象領域の内部及び境界の熱流束q分布、温度変化ΔT分布を導出することが好ましい。また、一方の領域で得られた解析結果を他の領域と接する境界にて、他の領域における次の解析の境界値として受け渡し、物理量の保存及び整合性が合うように連成解析する方法に対して、有限要素法もしくは境界要素法の解析、等価回路近似に基づいた解析の連成解析を、少なくとも2回以上繰り返し反復計算することが好ましい。また、有限要素法もしくは境界要素法の解析、等価回路近似に基づいた解析の連成解析を、少なくとも2回以上繰り返し反復計算する方法に対して、各領域の各手法で求めた、境界における最新の温度変化ΔTの残差を求め、その残差に緩和係数ω≦1を掛けた値を、前回の解析あるいは前々回の解析結果の値に加算した値を次の解析のための境界値として受け渡すことが好ましい。また、各領域の各手法で求めた、境界における最新の温度変化ΔTの残差を求め、その残差に緩和係数ω≦1を掛けた値を、前回の解析あるいは前々回の解析結果の値に加算した値を次の解析のための境界値として受け渡す方法に対して、前記緩和係数ωを前記連成解析の反復計算の途中で、値を変動させることが好ましい。また、緩和係数ωを前記連成解析の反復計算の途中で、値を変動させる方法に対して、反復計算回数が少ない時は小さい値ω≦05に設定し、反復計算回数が多くなった際に、緩和係数ωとして大きな値<05<ω≦10になるように前記緩和係数ωを反復計算の途中で、値を変動させることが好ましい。また、解析計算する手段として、PCクラスターまたはマルチコアPC乃至はマルチスレッドPCを用いる場合は、PC毎またはコア乃至はスレッド毎に各分割領域、各部分領域毎の計算を割り当て実施し、全体計算を高速化することが好ましい。また、熱伝導もしくは熱応力解析に対して、前記解析対象の物理量以外でも、熱解析と同様に、スカラーポテンシャルにより場の記述が可能な物理量の解析を特徴とする。また、解析領域を少なくとも2つ以上の領域に分ける方法、乃至は、前述した等価回路近似に基づいた解析の対象とする解析領域に対して、更に細かく、複数に細分した解析領域に分割する方法に対して、前記領域分割、前記細分領域分割に関して、解析を実行するユーザが入力、乃至は、設定するユーザインターフェース機能を備えることが好ましい。本発明は、パワエレシステム全体系などの大規模な対象に対して、現実的な計算コストの範囲内で、高精度な熱伝導解析することが可能となる。以下、図面を参照し、本発明の実施例について説明する。まず、第一の実施例について説明する。図1は本発明の数値解析プログラムの処理フローを表し、図2は、本発明の特徴を示す、数値解析システム全体の処理フローを表す。また、図3に本発明の数値解析システムのハードウェア構成を、図4に本発明の概要を示す。まず、図3を用いて、本実施例のハードウェア構成を説明する。本実施例の数値解析システム1は、ハードウェア構成としては、計算機2、入力部3、出力部4、表示部5、記録部6、データベース7から構成される。ここで、出力部4と表示部5は併用されても良い。本発明に基づいた処理により本発明の数値解析を実行するプログラム10は、計算機2に付随する、ハードディスク等の記録部6の中に保存されており、計算機2、乃至は、計算機2にネットワーク9で接続されている他の計算機8の演算能力を用いて、本発明に基づいた本発明の数値解析を実行するプログラム10による数値解析が実行される。ここで、計算機2或いは計算機8は、PC、乃至は、PCクラスター、乃至は、マルチコアPC、乃至は、マルチスレッドPC、乃至は、スーパーコンピュータなどの数値演算可能な計算機全般を示す。数値解析するユーザが入力部3を用いて、数値解析対象のモデルを作成し、必要な解析条件を入力して、数値解析され、その結果が出力部4により出力され、表示部5により表示される。次に、本実施例の数値解析システム1、及び、本実施例に基づいた処理により本発明の数値解析を実行するプログラム10の処理フローを各々、図1、図2を用いて説明する。まず、図2を用いて、数値解析システム全体の大きな処理の流れを説明する。ユーザ入力18として、ユーザが解析対象領域全領域の指定をし、全領域を複数の領域に大分割する為の諸条件例えば、領域分割ラインの指定などを設定入力し、また、大分割した各領域内の解析手法の選択及び諸条件を設定し、等価回路近似に基づいた解析対象領域を、更に細分割化する場合の諸設定、及び、各分割領域の解析手法間の選択及びその諸設定を入力し、本発明の数値解析を実行するプログラム10を実行し、その結果を表示するというのが大きな処理の流れである。次に、本発明の数値解析を実行するプログラム10の詳細処理について図1を用いて説明する。ユーザ入力18により指定された諸条件に基づいて、全体の領域を領域分割11し、分割された各領域に対してユーザが指定した解析手法を割り当てる12。このとき、少なくとも一つの領域は、等価回路近似に基づく解析法が割り当てられ、他の領域は、FEMもしくはBEMの解析手法が選択される。ここで、図4に示すような熱源とヒートシンクが一つずつ存在し、銅配線をFR4などの絶縁材で挟んだ体系を適用例として示す。この例では、全体系の領域が図4に示すように2領域に分割し、熱源が存在する方の領域をFEMもしくはBEM対象領域、ヒートシンクが存在する方の領域を価回路近似に基づく解析対象領域とする。本実施例では、FEMもしくはBEM解析対象の分割領域内の処理13を先に実行し、その後、等価回路近似に基づく解析対象の分割領域内の処理15を実施することを特徴とする。具体的には、FEMもしくはBEM対象領域領域Aの解析13では、対象の分割領域に対する計算モデル、メッシュ作成などの準備13-1を実施するが、このとき、特に、大分割した各領域間の境界条件として、温度変化ΔTB=0と仮に初期設定する。13-1の準備に基づき、FEMもしくはBEMの解析13-2を実施し、得られた結果13-3に対して、大分割した各領域間の境界上の熱流束qB分布から、境界面での熱量QBを熱流束qBの面積分などの方法で算出し、もう一方の領域領域Bでの解析の境界条件として受け渡す14。ここで、等価回路近似に基づく解析対象の分割領域領域B内の処理15の手法について、図5を用いて説明する。等価回路近似に基づく解析対象として選定した領域領域Bの境界面に、図4、図5に示すように、銅配線断面が表出している場合、銅配線断面がその周囲の絶縁材より熱伝導率が3桁程度大きいことから、この銅配線断面を通してのみ、この領域から他の領域への熱伝導は行われると仮定し、これらの銅配線断面をポート19と称する。他の領域と熱伝導可能な端子のような役割であり、このポートに対して、その領域内の熱回路方程式[ΔT]=[R][Q]を求める。ここで、[]は数列または、行列を表し、その大きさは、ポートの数をNpとすると、数列の場合、Np、或いは、Np-1の数列、行列の場合、Np×Np、或いはNp-1×Np-1の行列を表す。つまり、該当する領域の熱抵抗行列[R]、あるいは、その逆行列である熱アドミッタンス行列[Y]=[R]1を求めることで、熱回路方程式[ΔT]=[R][Q]から、各ポートの熱量Qを入力すれば、各ポートのΔTが算出されるという解析手法である。この方法によると、ブロック内部の複雑さに関係なく、ポートの熱流を忠実に表現できる一方、一度、該当する領域の熱抵抗行列[R]あるいは、その逆行列である熱アドミッタンス行列[Y]=[R]1を求めてしまえば、熱回路方程式[ΔT]=[R][Q]による単なる行列演算で温度変化ΔTが算出される為、非常に高速に温度変化ΔTが評価可能となる。上記に説明した解析方法に基づいて、次に、等価回路近似に基づく解析対象の分割領域領域B内の処理15を実施する。このときの準備15-1としては、該当する領域の計算モデル作成や、メッシュ生成などの準備であり、次に、解析15-2としては、該当する領域の熱抵抗行列[R]あるいは、その逆行列である熱アドミッタンス行列[Y]=[R]1を一度だけ求め、その熱抵抗行列[R]、及び、領域Aの解析結果として既に得られている境界面のポートの熱量QBを用いて、熱回路方程式[ΔT]=[R][Q]に基づき、境界面ポートのΔTBを算出する。ここで、特記すべき点は、解析15-2での、該当する領域の熱抵抗行列[R]あるいは、その逆行列である熱アドミッタンス行列[Y]=[R]1を求めるのは、当該領域の形状が変わらない限り、一度だけでよいということである。従って、一度得られた熱抵抗行列[R]あるいは、その逆行列である熱アドミッタンス行列[Y]=[R]1をデータベース7に格納し、同じ形状の領域で等価回路近似に基づく解析を行う際に、データベース7から、該当する熱抵抗行列[R]あるいは、その逆行列である熱アドミッタンス行列[Y]=[R]1を呼び出し、何度も利用することで、解析の高速化が期待できる。以上により得られた結果15-3である境界面のポートの温度変化ΔTBから、物理量の保存及び整合性が合うように、前記境界における温度変化の初期設定としてΔTB=0として有限要素法もしくは境界要素法で解析して得られていたΔT分布に対して、当該解析領域の境界における温度変化ΔTB分の値を加算して、温度変化分布を加工調整処理16を実施して、温度変化ΔT′=ΔTΔTBである全体系の連成解析結果17が得られる。このとき、実際に本実施例を用いて得られた結果を図6に示す。全体系をFEMで計算した結果と誤差10以内で一致していることが確認できた。次に、本発明の第二の実施例について、図7、図8、図9を用いて説明する。ここでは、図7、図8に示すように、第一実施例と上流解析と下流解析の関係が逆転している。つまり、等価回路近似に基づく解析対象の分割領域内の処理15を先に実行し、その後、FEMもしくはBEM解析対象の分割領域内の処理13を実施することを特徴とする。具体的には、有限要素法もしくは境界要素法で解析する領域領域Aと、等価回路近似に基づいた手法の対象領域領域Bとの境界における熱量Qを、全体系に存在する全ての熱源から予測した熱量Qを設定し、それにより、まず、前記等価回路近似に基づく解析対象の分割領域内の処理15を対象領域領域Bに対して、熱等価回路方程式ΔT=RQに基づいて、当該解析領域の境界における温度変化ΔTBを導出し、もう一方の領域領域Aでの解析の境界条件として受け渡す20。次に、前記境界の温度変化ΔTBを境界初期条件として、有限要素法もしくは境界要素法にてFEMもしくはBEM解析対象の分割領域内の処理13を実施し、対象領域領域Aの内部及び境界の熱流束q分布、温度変化ΔT分布を導出する。図7の処理フローに示すように、本実施例によると、FEMもしくはBEM解析対象の分割領域内の処理13が、境界の初期条件ΔTBが既に等価回路近似に基づく解析手法の結果である為、第一実施例の処理フローで必要であった温度変化分布の加工調整処理16が不要となる。また、このとき、実際に本実施例を用いて得られた結果を図9に示す。全体系をFEMで計算した結果と誤差10以内で一致していることが確認できた。次に、本発明の第三の実施例について、図10、図11、図12を用いて説明する。ここでは、図4、図8に示した、等価回路近似に基づく解析対象の領域領域Bに対して、図10に示すように、更に細分割することを特徴とする。例えば、体系が複雑になり、それに伴い熱流も複雑な場合など、等価回路近似に基づく解析対象の領域領域Bを細分割せずに、内部の熱流を熱抵抗行列[R]やアドミッタンス行列[Y]で表現することが困難な場合がある。そのような場合は、図10に示すように、等価回路近似に基づく解析対象の領域領域Bに対して、更に、細分割すると解析精度が向上する場合がある。細分割した場合の等価回路近似に基づく解析手法の概要を図11に示す。図11に示すように、細分割した領域i毎に、これまでと同様に熱抵抗行列[R]iやアドミッタンス行列[Y]iを求め、それを足し合わせて合成行列Σ[R]i、または、Σ[Y]iを生成することで、等価回路近似に基づく解析対象の領域領域Bの全体系の熱抵抗行列[R]やアドミッタンス行列[Y]を求める方法である。この合成行列を用いて、これまでと同様、熱回路方程式[ΔT]=[R][Q]に基づき、各ポートのΔTBを算出する。但し、ここでは、ポートは、FEMもしくはBEM対象領域領域Aに対する境界のみではなく、等価回路近似に基づく解析対象の領域領域Bの細分領域間の境界にも存在する。また、細分領域分割に伴い、準備15-1としては、該当する領域の計算モデル作成や、メッシュ生成の準備の他に、細分領域分割に関する設定、細分割などの処理が必要になる。ここで、細分領域に分割する指標として、等価回路近似に基づいた手法にて解析実行の対象とする当該解析領域に対して、主な熱の流入部から主な熱の流出部までの主な熱の流れの方向に、平行に近い細分領域分割線の数に比べ、前記主な熱の流れの方向に垂直に近い細分領域分割線の数が等しい乃至はそれ以上となるように分割すると、精度が更に向上することが期待できる。以上、本実施例に依れば、等価回路近似に基づく解析対象の領域領域B内のポートが増える分、内部の温度変化の評価ポイントが増え、更に、精度が向上することが期待できる。また、このとき、実際に本実施例を用いて得られた結果を図12に示す。全体系をFEMで計算した結果と誤差10以内で一致していることが確認できた。但し、ここでは、FEM解析から等価回路近似に基づいた解析の流れとなっているが、実施例2のように、解析順番が逆になっても良い。次に、本発明の第四の実施例について、図13を用いて説明する。ここでは、細分領域に分割する指標乃至は細分割領域の単位として、部品乃至はその周辺をも含めた部品を一つの細分領域の単位として分割する、乃至は、形状を変更しない領域と形状を変更し得る領域への細分領域に分割することを特徴とする。つまり、図13に示すような、部品単位を規準として、細分割した場合、部品毎に熱抵抗行列[R]やアドミッタンス行列[Y]を一度導出してしまえば、部品毎に熱抵抗行列[R]やアドミッタンス行列[Y]がデータベース7に格納される為、同じ部品が実装されている場合は、その領域の解析は非常に高速になることが期待できる。次に、本発明の第五の実施例について、図14を用いて説明する。ここでは、第一実施例及び第二実施例で示したような、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15を連成解析するフローに対して、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15の連成解析を、FEMもしくはBEM解析対象の分割領域内の処理13と等価回路近似に基づく解析対象の分割領域内の処理15との連成解析を2回以上繰り返し反復計算処理21することを特徴とする。本実施例の処理フローを図14に示す。FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15の連成サイクルを2回以上繰り返し21、次の式のような判定式に基づき、反復計算の継続か中止かを判定する。-ΔTFEM/ΔTY<ε...判定22または-ΔTFEM/ΔTY<ε...判定23ここで、ΔTY、ΔTFEMは、各々、等価回路近似に基づく解析対象の分割領域内の処理15の解析結果の温度変化、FEMもしくはBEM解析対象の分割領域内の処理13の解析結果の温度変化を示し、上付きnは繰り返し回数nを表す。つまり、2つの手法による結果の違いや誤差を評価している式となる。さらに、εは反復計算継続化中止かを判定する基準値であり、システムが予め準備した値を採用しても良く、また、ユーザ自身に指定してもらっても良い。例えば、ε=01、或いは、ε=10などの1未満の値が良い。但し、ここで、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15における、解析準備処理13-1,15-1のうち、解析モデル作成やメッシュ作成に関しては、n=1回目の反復回数時のみでよい。つまり、一度だけ処理をすれば、その後の処理は不要となる。したがって、本実施例によれば、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15の結果に乖離がある場合等、繰り返し計算できる為、計算結果の信憑性が向上し、また、指定した尤度内に収束したと判定した時点で計算を中止できる為、余計な計算をしなくて良い。次に、本発明の第六の実施例について、図15、図16を用いて説明する。ここでは、第五の実施例に対して、特に、等価回路近似に基づく解析対象の分割領域内の処理15の結果である、領域境界における解析結果ΔTYに対して、FEMで得られた結果ΔTFEMとの残差に緩和係数Cωをかけた値を等価回路近似に基づく解析結果をFEMもしくはBEM解析の境界条件として設定する処理20において、具体的には、次の式に基づき、FEMもしくはBEM解析の境界条件の設定処理24することを特徴とする。ΔTFEM=ΔTFEMCω-ΔTFEM...判定24但し、ここで緩和係数はCω<1とする。つまり、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15で、結果を徐々に近づけるといった方法である。本実施例の効果を実際に計算により確認した結果を図16に示す。FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15の結果に乖離がある場合等、緩和係数はCω=1では解ΔTが発散振動してしまうが、緩和係数はCω=05<1では、解を徐々に近づけることにより、安定に収束している。つまり、本実施例に依れば、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15の結果に乖離がある場合等に、解が安定収束することが期待できる。次に、本発明の第七の実施例について、図17、図18を用いて説明する。ここでは、第六の実施例に対して、等価回路近似に基づく解析対象の分割領域内の処理15の結果である、領域境界における解析結果ΔTYに対して、FEMで得られた結果ΔTFEMとの残差に緩和係数Cω<1をかけた値を反復計算の継続か中止の判定処理22において、途中で緩和係数Cωの値を変えることを特徴とする。図181に示すように、緩和係数はCωの値が小さいほど、安定確実に解は収束するが、反復回数を多く必要とするため、非常に計算時間がかかってしまう。一方、図182に示すように、緩和係数はCωの値がCω<1の範囲内で大きくすれば、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15で、結果が早く近づくが、解が振動してしまう。そこで、本実施例では、等価回路近似に基づく解析対象の分割領域内の処理15の結果が早く収束することに着目して、次式の緩和係数Cω値の変更判定処理25に基づき、等価回路近似に基づく解析結果の残差を求めて、その値が、ある決められた値α以下になった場合に、緩和係数Cω値の変更設定処理26し直す。-ΔTY/ΔTY<α...判定25但し、上記判定後に、設定する緩和係数Cωの値はCω=1でよい。つまり、Cω=1と設定し直しても良い。本実施例の処理フローを図17に示す。また、本実施例の効果を実際に計算により確認した結果を図183に示す。最初、緩和係数はCω=05で、安定に解を収束させ、上記の反復計算の継続か中止の判定処理23で以下になった際に、緩和係数はCω=10に設定変更し、収束を加速できているのがわかる。従って、本実施例によると、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15の結果に乖離がある場合等でも、解が早く、安定収束することが期待できる。次に、本発明の第八の実施例について説明する。これまでは、熱伝導解析を中心に、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15の連成解析について記載してきたが、ここでは、熱応力も同様に、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15の連成解析により算出することを特徴とする。熱応力の一般式を下記する。σ=Eεr=EαΔTここで、Eはヤング率、εrは熱ひずみ、αは線膨張係数である。本実施例では、温度変化ΔTまではこれまでの実施例と同様に求め、境界面に対して、垂直な熱応力として、上記式により求める。熱応力解析は、一般的に3次元解析が必要となり、計算コストが多大となるため、本実施例に依れば、高速高精度な熱応力解析が可能となる。次に、本発明の第九の実施例について説明する。これまでは、熱伝導解析を中心に、FEMもしくはBEM解析対象の分割領域内の処理13と、等価回路近似に基づく解析対象の分割領域内の処理15の連成解析について記載してきたが、ここでは、熱伝導解析と同様に、スカラーポテンシャルにより場の記述が可能な物理量の解析を対称とすることを特徴とする。例えば、電場などがその一例である。スカラーポテンシャルによって、場の記述が可能であれば、本発明のような領域を分割して、それを組み合わせるといった線形性が保たれる為、基本的には本発明を適用可能と考える。次に、本発明の第十の実施例について、図3を用いて説明する。ここでは、図3に示す計算機2に対して、PCクラスターまたはマルチコアPC乃至はマルチスレッドPCを用いる場合は、PC毎またはコア乃至はスレッド毎に各分割領域、各細分割領域の計算を割り当て実施し、全体計算を高速化することを特徴とする。互いに独立な計算ほど本実施例に効果が大きい為、例えば、等価回路近似に基づく解析対象の領域領域Bに対して、細分割した領域毎の熱抵抗行列[R]やアドミッタンス行列[Y]を求める計算を、本実施例のようなマルチ演算処理化することで大幅に計算速度が向上する。"