JP4097176B2 - 熱解析方法および装置 - Google Patents
熱解析方法および装置 Download PDFInfo
- Publication number
- JP4097176B2 JP4097176B2 JP2001235556A JP2001235556A JP4097176B2 JP 4097176 B2 JP4097176 B2 JP 4097176B2 JP 2001235556 A JP2001235556 A JP 2001235556A JP 2001235556 A JP2001235556 A JP 2001235556A JP 4097176 B2 JP4097176 B2 JP 4097176B2
- Authority
- JP
- Japan
- Prior art keywords
- heat transfer
- thermal analysis
- transfer coefficient
- deviation
- temperature
- 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.)
- Expired - Lifetime
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/18—Investigating or analyzing materials by the use of thermal means by investigating thermal conductivity
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Complex Calculations (AREA)
Description
【発明の属する技術分野】
本発明は、例えば、内燃機関のピストンの温度分布解析などに利用することができる熱解析方法および装置に関し、さらに詳細には応答曲面法を用いた熱解析方法および装置に関する。
【0002】
【従来の技術】
熱解析を行う際には、例えば、特開平11−118740号公報に示されているように、自然対流熱伝達が生じる熱伝達境界面における熱伝達率を計算する実験式を、諸条件に応じて自動的に選択し、熱伝達率を決定するために仮想した熱伝達境界面の温度値を、有限要素法(FEM)あるいは有限差分法などの手段によって得られた熱伝達境界面の温度との整合性を取りながら繰り返し計算によって自動的に求めている。
【0003】
【発明が解決しようとする課題】
しかし、上記したように、自然対流熱伝達が生じる熱伝達境界面における熱伝達率を計算するための実験式を、諸条件に応じて自動的に選択するにしても、対象物の形状モデルにあった適切な実験式を選択するためと熱伝達率を決定するためには膨大な繰り返し計算が必要であり、膨大な手数と時間を必要とするという問題があった。
【0004】
また、実験式は解析される物体形状毎に必要であって、実験式の選択が自動化されても、対象物の細部が変更された場合には、変更の度毎に実験式を作成し、さらに熱伝達率を決定するために繰り返し計算を実行しなければならず、膨大な手数と時間を必要とするという問題があった。
【0005】
本発明は上記の問題点に鑑みなされたもので、対象物体の形状にとらわれずに、汎用的に温度分布を得ることができる熱解析方法および装置を提供することを目的とする。
【0006】
【課題を解決するための手段】
本発明の請求項1および3にかかる熱解析方法および装置は、熱解析の対象物形状を有限個の有限要素あるいはセルに分割する対象領域分割を行い、この対象領域分割結果を用いて前記対象領域の熱解析を行う熱解析方法および装置において、
領域分割された対象形状表面の複数位置における熱伝達率を用いて対象物中の温度分布を求め、求めた温度分布に基づいて対象形状中における特定点の温度を求め、求めた前記特定点の温度と該特定点における予め与えた目標温度との偏差を求め、前記偏差が所定範囲内か否かを判別し、前記偏差が所定範囲外と判別されたときは熱伝達率を変更し、変更された熱伝達率で対象物の温度分布を求めるために熱伝達率を更新して、前記偏差が前記所定範囲内になるまで繰り返すことを特徴とする。
【0007】
本発明の請求項1および3にかかる熱解析方法および装置によれば、領域分割された対象形状表面の複数位置における熱伝達率を用いて対象物中の温度分布が求められ、求められた温度分布に基づいて対象形状中における特定点の温度が求められ、特定点の求められた温度と対象物が所望強度を得るように前記特定点における予め与えた目標温度との偏差が求められ、前記偏差が所定範囲内か否かが判別されて、前記偏差が所定範囲外と判別されたとき熱伝達率が変更され、変更された熱伝達率に前記温度分布を求めるための熱伝達率が更新されて、前記偏差が前記所定範囲内になるまで繰り返される。
【0008】
上記のように本発明の請求項1および3にかかる熱解析方法および装置によれば、自動的に繰り返して熱伝達率が演算によって求められるために、実験式の選択を行う必要もなく、温度分布と熱伝達率を効率よく得ることができる。
【0009】
さらに、熱伝達率と温度分布とを連続して得ることができるため、最終的な温度分布と熱伝達率とがほぼ同時に求められて作業効率が向上する。
【0010】
また、本発明の請求項1および3にかかる熱解析方法および装置において、第1の行程および手段は熱伝達率を予測変数とし、熱伝達率を所定の範囲内で変化させて得た複数の近似式を連立方程式とし、該連立方程式を解いて応答曲面を求め、該応答曲面を対象物の温度分布とするようにしてもよい。
【0011】
また、このようにしたときは、実験式によるわけではなく、近似式の連立方程式を熱伝達率に関して解くことによって得られた熱伝達率を用いて温度分布が得られ、偏差の収束が時間的に早く得られて、最終的な熱伝達率や温度分布が求まるまでの繰り返し回数は少なくて済む。
【0012】
【発明の実施の形態】
以下、本発明にかかる熱解析方法および装置を実施の形態によって説明する。
【0013】
図1は、本発明の実施の一形態にかかる熱解析方法が適用される熱解析装置の構成を示すブロック図である。
【0014】
図1に示した熱解析装置10を内燃機関のピストンに対して適用した場合を例に説明する。
【0015】
熱解析装置10は、入力装置1と、図示しないROMと、作業領域を有するテンポラリメモリ3と、ROMおよびテンポラリメモリ3と協働する制御装置2と、出力装置4とを備えている。
【0016】
入力装置1からピストンの3次元データ、熱伝達率の初期値、予め定めた特定点の温度を入力する。
【0017】
制御装置2は、3次元データからピストンの3次元モデルを作成し、FEM解析のために3次元モデルの領域分割を行う形状算出部21と、領域分割されたピストン表面の複数位置における熱伝達率を用いてピストンの温度分布を求める温度分布算出部22と、求めた温度分布からピストンの対象形状中における特定点の温度を求める特定点温度算出部23と、求めた特定点の温度と該特定点に対して予め与えた特定点の温度とから偏差を求め、求めた偏差が所定範囲内か否かを判別する偏差算出・判別部24と、偏差が所定範囲外と判別されたとき熱伝達率を変更し、変更された熱伝達率を温度分布算出部22にて用いる熱伝達率として、偏差が所定範囲内に入るまで繰り返させる熱伝達率更新部25とを機能的に備えている。
【0018】
さらに詳細には、温度分布算出部22は、応答曲面を決定するための解析を含み、FEM解析による解析による温度分布解析を行う。
【0019】
先ず、熱解析装置10の作用の説明に先立って、熱解析装置10が適用される解析システムについての概略を説明する。
【0020】
対象物であるピストンのモデルが作成され、続いてFEM解析のためのメッシュ、またはセルに分割される。次いでピストンの諸元、計算条件が入力され、温度分布解析が行われ、温度分布解析結果に基づいて爆発荷重と慣性力と熱応力等に対する強度解析が行われ、続いて疲労計算、安全率評価などが行われる。
【0021】
熱解析装置10は、上記の温度分布解析において用いられる。
【0022】
ここで、温度分布解析を行う理由は次の如くである。
【0023】
最高300℃にも達する高温にさらされるピストンに用いられるアルミニウム合金材料は温度による強度低下が大きい。このため、ピストン各部の温度分布を求めておくことにより、爆発力と慣性力とを組み合わせることによる温度の影響を考慮した疲労限界線図によってピストンの強度評価、応力評価が行えるためである。
【0024】
上記のように構成された熱解析装置10の作用について図2のフローチャートに基づいて説明する。
【0025】
プログラムが開始されると、入力装置1から入力された3次元データに基づいて形状算出部21により、対象物であるピストンの3次元モデルが作成され、続いてFEM解析のためのメッシュ、またはセルに分割される(ステップS1)。続いて、入力装置1を介して熱境界条件が設定される(ステップS2)。
【0026】
熱境界条件として、図3に模式的に示すように、ピストン51の天井52からの入熱、ピストンリング53からの出熱、ピストン51の内側からの出熱、ピストン51の外側からの出熱、ピストンピン55の内側からの出熱およびコネクチングロッド(コンロッド)56端からの出熱に対応する熱伝達率が初期条件として、入力装置1において設定される。図3において符号54はピストンピン55が挿入されるピストンピンボスを示している。
【0027】
実際の内燃機関においては、ガスの燃焼によって熱が発生し、ピストン51、ピストンピン55、コンロッド56、図示しないシリンダなどの部品を介して放熱し、最終的には冷却水、オイルなどを介して放熱される。
【0028】
この現象を完全にシミュレートすることができないために、上記のように熱の出入りをパラメータとして表し、熱解析を行うのである。本発明の実施の形態では、熱の出入りのパラメータとして熱伝達率を用いる。
【0029】
具体的には、パラメータとしての各部熱伝達率が図4に示すように、ピストン51の天井52の熱伝達率P1、ピストン51の側面の熱伝達率P2、トップリング53−1の熱伝達率P3、セカンドリング53−2の熱伝達率P4、オイルリング53−3の熱伝達率P5、ピストン51の内面の熱伝達率P6、ピストンピン55の熱伝達率P7の7つとした場合を例に説明する。
【0030】
ピストン天井52の熱伝達率P1、ピストン51の側面の熱伝達率P2、トップリング53−1の熱伝達率P3、セカンドリング53−2の熱伝達率P4、オイルリング53−3の熱伝達率P5、ピストン51の内面の熱伝達率P6、ピストンピン55の熱伝達率P7の7つが初期条件として設定される。
【0031】
さらに、ピストン51の断面において図5に黒丸で示す予め定めた特定点における目標温度が設定される。この目標温度は所望強度のピストン51を得るためにそれぞれの強度着目点である特定点における目標温度であり、予めステップS2において、熱伝達率P1〜P7と共に初期条件として入力装置1を介して設定される。
【0032】
初期条件の設定に続いて、初期条件として設定された熱伝達率P1〜P7が読み込まれ(ステップS3)、対象形状、すなわちピストン51中における予め定めた特定点の目標温度が読み込まれる(ステップS4)。
【0033】
次いで、熱伝達率P1〜P7が用いられて温度分布算出部22において応答曲面法に基づいて応答曲面が求められる(ステップS5)。
【0034】
ステップS5における応答曲面モデルの作成方法について説明する。
【0035】
いま、各パラメータである7つの熱伝達率P1〜P7を予測変数とし、応答関数として2次多項式を採用した場合、応答曲線法に基づき応答曲面式は下記の式(1)にて示される。
【0036】
【数1】
【0037】
ここで、式(1)におけるa0、bi、cj、dkはそれぞれ未知の値であり、熱伝達率P1〜P7を、所定の範囲、例えば(0.1〜1.0)の範囲内で何回か変動させて、変動の結果に基づいて7つ以上の方程式を生成し、そのうちの7つの方程式を連立方程式として、熱伝達率P1〜P7について解いて、応答曲面を求める。
【0038】
ステップS5をさらに説明すれば、例えば、図7に示す如く熱伝達率P1〜P7(各パラメータ)の各々に対して水準の決定が行われる(ステップS51)。各パラメータを2水準とした場合を例示すれば、図8に示すように、例えば各パラメータの水準を2に揃え、例えば下記の表1に示すようにL8直交表に割り付けられる。
【0039】
【表1】
【0040】
直交表に割り付けられた水準に基づき熱伝達率P1〜P7の寄与が同等になるようにN回の解析が設定される(ステップS52)。
【0041】
ステップS52に続いて設定されたN回の解析が繰り返し行われて(ステップS53)、解析結果に基づいて最小自乗法などを用いて応答曲面が作成される(ステップS55)。
【0042】
ステップS55において作成された応答曲面から特定点の温度が目標温度となるように応答曲面モデルより熱伝達率P1〜P7が決定される(ステップS6)。その熱伝達率P1〜P7を用いて熱伝達解析を行うことで対象領域であるピストン51の温度分布が得られる(ステップS7)。
【0043】
ステップS7において算出された対象領域の温度分布は模式的に図5において破線で示す如くである。具体的には濃淡図形で示される。図5はピストンピン55と直交する方向でのピストン51の断面図であり、二点鎖線はピストンピンボス54の位置を示している。図6もピストン51の断面図であって、図5と直交する方向で切断したときのピストン51の中心線に対して左の一方側のみを示し、二点鎖線はピストンピンボス54の中心線の位置を示している。
【0044】
ステップS7において得られた温度分布から、特定点温度算出部23によって、前記特定点における温度が算出され(ステップS8)、特定点における算出温度と対応する特定点の目標温度との偏差△Tが偏差算出・判別部24により求められる(ステップS9)。
【0045】
ステップS9において求められた偏差△Tが予め定めた所定範囲内であるか否かが、偏差算出・判別部24によってチェックされる(ステップS10)。
【0046】
ステップS10において偏差△Tが予め定めた範囲内でないと判別されたときは、熱伝達率更新部25により、これまでに得られた解析結果とステップS8の解析結果より応答曲面モデルが更新され(ステップS11)、ステップS6から再び実行される。
【0047】
この実行により偏差△Tが予め定めた所定範囲内に入ったと判別されると、偏差△Tが予め定めた所定範囲内に入ったときの熱伝達率P1〜P7が、ピストン51の熱伝達率として決定され、その熱伝達率P1〜P7に基づくピストン51の温度分布が得られる。
【0048】
上記のステップS6からステップS11の繰り返しによって、熱伝達率P1〜P7(絶対値)の変動の様子の一例は図8に示す如くであり、熱伝達率P1〜P7が目標値(偏差△Tが収束したときの値)へ収束していく様子の一例は図9に示す如くであり、偏差△Tが収束していく様子の一例は図10に示す如くであった。この例では5回の繰り返しによって偏差△Tは所定範囲内に入り、従来の手法による場合に比較して5分の1の時間で所定の設定目標温度範囲に収束した。
【0049】
図10において、オイルリング中央部、セカンドリング中央部、トップリング中央部、ピストン天井肉抜き部(a)、ピストン天井上中央部(b)、ピストン下中央部(c)、リブ中央部(d)、ピストンピンボス面取り内側上部(e)、ピストンピンボス面取り内側油溝上部(f)およびピストンピンボス面取り内側下部(g)における偏差△Tが収束していく様子を示している。
【0050】
ここで、オイルリング中央部、セカンドリング中央部、トップリング中央部を除いて、ピストン天井肉抜き部(a)、ピストン天井上中央部(b)、ピストン天井下中央部(c)、リブ中央部(d)、ピストンピンボス面取り内側上部(e)、ピストンピンボス面取り内側油溝上部(f)およびピストンピンボス面取り内側下部(g)の各部は、ピストンピン55と直交する方向でのピストン51の断面図である図11(a)および図11(a)の底部からみたピストン51の模式図である図11(b)において、図10における符号と同一の符号(a)〜(g)にて示してある。
【0051】
図11(a)および(b)において、符号59はピストンの内周側と外周側とを連通させる油溝であり、符号60はピストンピンボス54とピストン外周壁との間をつないでいるリブを示している。また、aはピストン天井裏面側において、ピストン軽量化のために肉厚を薄くされた部分である。
【0052】
以上説明したように、熱解析装置10によるときは、強度上着目点での目標温度に対して該着目点での温度変化の履歴から最適熱伝達率が得られ、最適熱伝達率による対象物における温度分布が迅速に得られる。
【0053】
ちなみに、強度上着目点の実測温度が欲しいときにも、その着目点の温度測定が可能とは限らず、着目点の周辺の温度を測定して測定結果に基づいて温度分布を推定するために、実測データに合うように熱伝達率を変えながら計算を繰り返す場合に比較して、熱解析装置10によるときは時間的に1/5から1/10程度の時間で済む。
【0054】
【発明の効果】
以上説明したように本発明にかかる熱解析方法および装置によれば、自動的に繰り返して熱伝達率が演算によって求められるために、実験式の選択を行う必要がなく、熱伝達率を効率よく得ることができるという効果が得られる。
【0055】
また、偏差の収束に至るまでの時間が少なくて済み、作業効率が向上するという効果も得られる。
【図面の簡単な説明】
【図1】本発明の実施の一形態にかかる熱解析装置の構成を示すブロック図である。
【図2】本発明の実施の一形態にかかる熱解析装置の説明に供するフローチャートである。
【図3】本発明の実施の一形態にかかる熱解析装置が適用されるピストンの模式断面図である。
【図4】本発明の実施の一形態にかかる熱解析装置における熱伝達率の説明に供する一部断面図である。
【図5】本発明の実施の一形態にかかる熱解析装置によるときの温度分布の説明図である。
【図6】本発明の実施の一形態にかかる熱解析装置によるときの温度分布の説明図である。
【図7】本発明の実施の一形態にかかる熱解析装置における温度分布の計算の説明に供するフローチャートである。
【図8】本発明の実施の一形態にかかる熱解析装置における繰り返し数対熱伝達率の変動の様子を示す説明図である。
【図9】本発明の実施の一形態にかかる熱解析装置における繰り返し数対熱伝達率の目標値への収束の様子を示す説明図である。
【図10】本発明の実施の一形態にかかる熱解析装置における繰り返し数対偏差の収束の様子を示す説明図である。
【図11】図10における測定個所を示す概略図である。
【符号の説明】
1…入力装置 2…制御装置
3…テンポラリメモリ 4…出力装置
10…熱解析装置 21…形状算出部
22…温度分布算出部 23…特定点温度算出部
24…偏差算出・判定部 25…熱伝達率更新部
Claims (4)
- 熱解析の対象物形状を有限個の有限要素あるいはセルに分割する対象領域分割を行い、この対象領域分割結果を用いて前記対象領域の熱解析を行う熱解析実行行程を有する熱解析方法において、
前記熱解析実行行程は、
領域分割された対象形状表面の複数位置における熱伝達率を用いて対象物中の温度分布を求める第1の行程と、
前記第1の行程によって求めた温度分布に基づいて対象形状中における特定点の温度を求める第2の行程と、
前記第2の行程によって求めた前記特定点の温度と対象物が所望強度を得るように前記特定点における予め与えた目標温度との偏差を求める第3の行程と、
前記偏差が所定範囲内か否かを判別する第4の行程と、
前記偏差が所定範囲外と判別されたとき熱伝達率を変更する第5の行程と、
を備え、
前記第5の行程によって変更された熱伝達率に前記第1の行程にて用いる熱伝達率を更新して、前記偏差が前記所定範囲内になるまで前記第1の行程から繰り返すことを特徴とする熱解析方法。 - 請求項1記載の熱解析方法において、第1の行程は熱伝達率を予測変数とし、熱伝達率を所定の範囲内で変化させて得た複数の近似式を連立方程式として、該連立方程式を解いて応答曲面を求め、該応答曲面を対象物中の温度分布とすることを特徴とする熱解析方法。
- 熱解析の対象物形状を有限個の有限要素あるいはセルに分割する対象領域分割を行い、この対象領域分割結果を用いて前記対象領域の熱解析を行う熱解析装置において、
領域分割された対象形状表面の複数位置における熱伝達率を用いて対象物中の温度分布を求める第1の手段と、
前記第1の手段によって求めた温度分布に基づいて対象形状中における特定点の温度を求める第2の手段と、
前記第2の手段によって求めた前記特定点の温度と対象物が所望強度を得るように前記特定点における予め与えた目標温度との偏差を求める第3の手段と、
前記偏差が所定範囲内か否かを判別する第4の手段と、
前記偏差が所定範囲外と判別されたとき熱伝達率を変更する第5の手段と、
を備え、
第5の手段によって変更された熱伝達率に前記第1の手段にて用いる熱伝達率を更新して、前記偏差が前記所定範囲内になるまで前記第1の手段から繰り返すことを特徴とする熱解析装置。 - 請求項3記載の熱解析装置において、第1の手段は熱伝達率を予測変数とし、熱伝達率を所定の範囲内で変化させて得た複数の近似式を連立方程式として、該連立方程式を解いて応答曲面を求め、該応答曲面を対象物中の温度分布とすることを特徴とする熱解析装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001235556A JP4097176B2 (ja) | 2001-08-02 | 2001-08-02 | 熱解析方法および装置 |
FR0209429A FR2828283B1 (fr) | 2001-08-02 | 2002-07-25 | Procede et appareil pour analyse thermique |
US10/202,884 US6885964B2 (en) | 2001-08-02 | 2002-07-26 | Method and apparatus for thermal analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001235556A JP4097176B2 (ja) | 2001-08-02 | 2001-08-02 | 熱解析方法および装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003042984A JP2003042984A (ja) | 2003-02-13 |
JP4097176B2 true JP4097176B2 (ja) | 2008-06-11 |
Family
ID=19066981
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001235556A Expired - Lifetime JP4097176B2 (ja) | 2001-08-02 | 2001-08-02 | 熱解析方法および装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US6885964B2 (ja) |
JP (1) | JP4097176B2 (ja) |
FR (1) | FR2828283B1 (ja) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004188382A (ja) * | 2002-12-13 | 2004-07-08 | Fuji Heavy Ind Ltd | 塗装乾燥予測方法、塗装乾燥予測システムおよび記録媒体 |
KR100796926B1 (ko) * | 2005-09-27 | 2008-01-22 | 동부제강주식회사 | 칼라 도장 강판 생산라인에서 강판온도 예측 방법 |
US8065122B2 (en) * | 2007-07-16 | 2011-11-22 | Durr Systems, Inc. | Method of designing or evaluating a bake oven |
US8180598B2 (en) * | 2009-02-10 | 2012-05-15 | Owens Corning Intellectual Capital, Llc | System and method for performing thermal analysis on a building through universal meshing |
US8214182B2 (en) * | 2009-05-12 | 2012-07-03 | GM Global Technology Operations LLC | Methods of predicting residual stresses and distortion in quenched aluminum castings |
JP5409497B2 (ja) * | 2010-04-16 | 2014-02-05 | 住友重機械工業株式会社 | 伝熱解析装置及び伝熱解析方法 |
CN102592024B (zh) * | 2012-01-06 | 2013-04-03 | 北京航空航天大学 | 一种确定径向热传导稳态温度极大值的热网络建模方法 |
FR2993360B1 (fr) * | 2012-07-13 | 2015-02-20 | Commissariat Energie Atomique | Procede de mesure de la conductivite thermique d’un materiau mince anisotrope |
DE102015007455B4 (de) | 2015-06-05 | 2018-10-31 | Audi Ag | Verfahren zum Betreiben einer Brennkraftmaschine mit einer Verringerung oder Deaktivierung einer Kolbenkühlung in Abhängigkeit von einer modellierten Kolbentemperatur sowie entsprechende Brennkraftmaschine |
JP6296045B2 (ja) * | 2015-12-08 | 2018-03-20 | トヨタ自動車株式会社 | 内燃機関の制御装置 |
CN106446364B (zh) * | 2016-09-08 | 2019-06-25 | 东南大学 | 一种温度场-热路直接耦合的电机热分析方法 |
WO2018225243A1 (ja) * | 2017-06-09 | 2018-12-13 | 富士通株式会社 | 熱解析装置、熱解析方法及び熱解析プログラム |
US11783099B2 (en) * | 2018-08-01 | 2023-10-10 | General Electric Company | Autonomous surrogate model creation platform |
CN109344460B (zh) * | 2018-09-14 | 2023-04-14 | 上海理工大学 | 一种传热结构仿生分层生长方法 |
CN109238724B (zh) * | 2018-10-24 | 2023-07-07 | 浙江大学城市学院 | 一种可调节燃烧室样件传热特性研究系统 |
CN114117852B (zh) * | 2021-11-18 | 2022-10-21 | 华北电力大学 | 一种基于有限差分工作域划分的区域热负荷滚动预测方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5202843A (en) * | 1990-04-25 | 1993-04-13 | Oki Electric Industry Co., Ltd. | CAE system for automatic heat analysis of electronic equipment |
JP3170375B2 (ja) * | 1993-01-25 | 2001-05-28 | 日新製鋼株式会社 | 熱延鋼板の温度予測方法 |
JP3042601B2 (ja) * | 1996-10-31 | 2000-05-15 | ファイラックインターナショナル株式会社 | 流体燃料の改質用セラミックス触媒を用いた内燃機関及 びこれを使用した交通輸送用又は発電用手段 |
US5931140A (en) * | 1997-05-22 | 1999-08-03 | General Motors Corporation | Internal combustion engine thermal state model |
JP3274095B2 (ja) * | 1997-07-18 | 2002-04-15 | 富士通株式会社 | 加熱炉内の被加熱物の熱解析装置及びそれを用いたリフロー炉の制御装置並びにそのプログラムを記録したコンピュータ読み取り可能な記録媒体 |
JPH11118740A (ja) * | 1997-10-17 | 1999-04-30 | Canon Inc | 熱解析方法及び熱解析装置 |
-
2001
- 2001-08-02 JP JP2001235556A patent/JP4097176B2/ja not_active Expired - Lifetime
-
2002
- 2002-07-25 FR FR0209429A patent/FR2828283B1/fr not_active Expired - Fee Related
- 2002-07-26 US US10/202,884 patent/US6885964B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2003042984A (ja) | 2003-02-13 |
FR2828283B1 (fr) | 2005-06-03 |
US20030060999A1 (en) | 2003-03-27 |
US6885964B2 (en) | 2005-04-26 |
FR2828283A1 (fr) | 2003-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4097176B2 (ja) | 熱解析方法および装置 | |
US8666706B2 (en) | Material property distribution determination for fatigue life calculation using dendrite arm spacing and porosity-based models | |
US8306790B2 (en) | Apparatus, method, and program for acoustic characteristics optimization model analysis | |
US20160034614A1 (en) | Materials property predictor for cast aluminum alloys | |
EP2486377B1 (en) | Rotor behaviour determination | |
GB2395821A (en) | Control process calibration combining initital map with correction function | |
KR101544457B1 (ko) | 최적 설계 파라미터 탐색을 위한 최적화 방법 | |
US20020177985A1 (en) | Computer system and method for radial cooled bucket optimization | |
JP6516308B1 (ja) | 応力解析装置 | |
JP7074015B2 (ja) | 探索装置 | |
JP2005306174A (ja) | タイヤ性能予測方法 | |
JP2004210106A (ja) | タイヤ経時変化予測方法、タイヤ特性予測方法、タイヤ設計方法、タイヤ製造方法、空気入りタイヤおよびプログラム | |
CN109871642B (zh) | 一种用于蒸汽发生器的内部结构的温度场建模方法 | |
CN110750936B (zh) | 一种确定销孔型线的方法及装置 | |
JP6065472B2 (ja) | タイヤモデル作成方法、タイヤ断面形状決定方法、タイヤモデル作成装置、及びプログラム | |
CN115720658A (zh) | 用于物理零件的3d打印的热量感知工具路径重新排序 | |
CN116050202A (zh) | 适应多种平均应力表现的Chaboche叶片疲劳寿命预测模型及方法 | |
KR102397107B1 (ko) | 반복적 전역최적화에 의해 도출된 인공오차 공분산 최적해에 기반한 분산점 칼만필터 연산 구성의 구조물 손상여부 모니터링 장치 | |
CN112782974B (zh) | 一种控制参数标定方法 | |
JP4998765B2 (ja) | エンジン性能の予測解析システム、予測解析方法及び予測解析プログラム | |
JP2006247746A (ja) | 溶接解析方法 | |
CN117790300B (zh) | 一种精细线路的动态蚀刻补偿方法 | |
US6871105B2 (en) | Methods and apparatus for reducing hyperplanes in a control space | |
CN115906338B (zh) | 一种烧结炉炉门端盖的优化设计评定方法 | |
JP2008031916A (ja) | 内燃機関の自着火予測装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041203 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20051003 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071211 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080206 |
|
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: 20080304 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080307 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110321 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110321 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120321 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130321 Year of fee payment: 5 |