JP2005140509A - Equivalent thermal conductivity calculation device and equivalent thermal conductivity calculation method for composite material - Google Patents
Equivalent thermal conductivity calculation device and equivalent thermal conductivity calculation method for composite material Download PDFInfo
- Publication number
- JP2005140509A JP2005140509A JP2003374010A JP2003374010A JP2005140509A JP 2005140509 A JP2005140509 A JP 2005140509A JP 2003374010 A JP2003374010 A JP 2003374010A JP 2003374010 A JP2003374010 A JP 2003374010A JP 2005140509 A JP2005140509 A JP 2005140509A
- Authority
- JP
- Japan
- Prior art keywords
- thermal conductivity
- equivalent thermal
- input
- boundary
- node
- 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.)
- Pending
Links
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
Description
本発明は、伝導境界から被伝導境界への熱伝導を単一部材として構成した場合に等価な熱伝導率とみなせる等価熱伝導率を出力する出力手段を備えている複合材の等価熱伝導率算出装置解析値により等価な熱伝導率を算出する装置および方法に関する。 The present invention provides an equivalent thermal conductivity of a composite material having output means for outputting an equivalent thermal conductivity that can be regarded as an equivalent thermal conductivity when the heat conduction from the conduction boundary to the conduction boundary is configured as a single member. The present invention relates to an apparatus and a method for calculating an equivalent thermal conductivity based on a calculated apparatus analysis value.
従来の技術では、等価な熱伝導率を算出する場合、複合材の体積(ただし、形状が特定の形状方向で同一の平面的形状の場合は面積)の割合により等価な熱伝導率算出を算出していた。この様な方法および装置としては、特許文献1の様な等価熱伝導率を算出するものがある。
しかし、上記方法では計算精度の点で問題がある場合がある。 However, the above method may have a problem in terms of calculation accuracy.
そこで、外形や物性値が違う部材間の境界が曲面や曲線である場合の様な複雑な形状でも、精度良く算出を実行することが可能な等価熱伝導率算出装置および等価熱伝導率算出方法が求められていた。 Therefore, an equivalent thermal conductivity calculation device and an equivalent thermal conductivity calculation method capable of accurately calculating even a complicated shape such as a curved surface or a curved boundary between members having different external shapes and physical property values Was demanded.
上記課題の解決を目的とした本発明の第一の発明は、条件を入力する入力手段と、伝導境界から被伝導境界への熱伝導を単一部材として構成した場合に等価な熱伝導率とみなせる等価熱伝導率を出力する出力手段を備えている複合材の等価熱伝導率算出装置において、
寸法、形状、物性値の特定により条件決定可能なモデルをそなえており、
(a)寸法、形状、物性値の入力を促す物性入力手段と、
(b)入力された形状をモデルに適合した上で要素分割する要素分割手段と、
(c)境界条件の入力を促す境界条件入力手段と、
(d)節点の温度を算出する節点温度決定手段と、
(e)節点の温度に基づき等価熱伝導率を算出する熱伝導計算手段と、
を具備することを特徴とする複合材の等価熱伝導率算出装置である。
The first invention of the present invention aimed at solving the above-mentioned problems is an input means for inputting conditions, and an equivalent thermal conductivity when the heat conduction from the conduction boundary to the conduction boundary is configured as a single member. In an equivalent thermal conductivity calculation device for a composite material provided with output means for outputting an equivalent thermal conductivity that can be considered,
We have a model that can determine conditions by specifying dimensions, shapes and physical properties.
(A) physical property input means for prompting input of dimensions, shapes, and physical property values;
(B) element dividing means for dividing the input shape into elements after conforming to the model;
(C) boundary condition input means for prompting input of boundary conditions;
(D) a node temperature determining means for calculating the temperature of the node;
(E) a heat conduction calculating means for calculating an equivalent thermal conductivity based on the temperature of the node;
An equivalent thermal conductivity calculation device for a composite material, comprising:
また、本発明の第二の発明は、伝導境界から被伝導境界への熱伝導を単一部材として構成した場合に等価な熱伝導率とみなせる等価熱伝導率を求める等価熱伝導率算出方法において、
(f)寸法、形状、物性値の入力を受け付ける物性入力処理と、
(g)入力された形状をモデルに適合する様に要素分割する要素分割処理と、
(h)境界条件の入力を受け付ける境界条件入力処理と、
(i)節点の温度を演算処理する節点温度決定処理と、
(j)節点の温度に基づき等価熱伝導率を演算する熱伝導計算処理と、
を具備することを特徴とする複合材の等価熱伝導率算出方法である。
Further, the second invention of the present invention is an equivalent thermal conductivity calculation method for obtaining an equivalent thermal conductivity that can be regarded as an equivalent thermal conductivity when the heat conduction from the conduction boundary to the conduction boundary is configured as a single member. ,
(F) physical property input processing for receiving input of dimensions, shapes, and physical property values;
(G) element division processing for dividing the input shape so as to conform to the model;
(H) boundary condition input processing for receiving input of boundary conditions;
(I) a nodal temperature determination process for computing the nodal temperature;
(J) a heat conduction calculation process for calculating an equivalent thermal conductivity based on the temperature of the node;
It is the equivalent thermal conductivity calculation method of the composite material characterized by comprising.
本発明によれば、外形や物性値が違う部材間の境界が曲面や曲線である場合の様な複雑
な形状でも、精度良く算出を実行することが可能な等価熱伝導率算出装置および等価熱伝導率算出方法を提供することが可能になった。
According to the present invention, an equivalent thermal conductivity calculation device and an equivalent heat that can be accurately calculated even in a complicated shape such as a case where a boundary between members having different external shapes and physical property values is a curved surface or a curved line. It has become possible to provide a conductivity calculation method.
以下に図面を用いて本発明の1実施形態を説明する。 Hereinafter, an embodiment of the present invention will be described with reference to the drawings.
図1において1は入力装置、2は処理装置(またはCPU)、3は表示装置(またはディスプレイ)、4は出力装置(またはプリンタ)、5は記憶装置である。 In FIG. 1, 1 is an input device, 2 is a processing device (or CPU), 3 is a display device (or display), 4 is an output device (or printer), and 5 is a storage device.
入力装置1は、寸法、形状、熱伝導率などの物性値を入力するものである。処理装置2はソフトウェアの働きに基づいて各装置を制御するとともに、下記の処理を行う。
(1)寸法、形状、物性値の情報の入力を促し、入力値を受け付ける処理
(2)入力された形状を要素分割する処理
(3)境界条件の設定の入力を促し、入力値を受け付ける処理
(4)有限要素法、差分法などによる節点温度を計算する処理
(5)等価熱伝導率の算出する処理
表示装置3は計算結果の表示を行う。出力装置4は計算結果によって得られたデータを記憶装置5に記憶した後、等価熱伝導率値を出力するものである。記憶装置5は上記の処理を行うプログラムや算出結果を記憶するものである。
The
(1) Process for prompting input of dimension, shape, and physical property value information and receiving input value (2) Process for dividing input shape into elements (3) Process for prompting input of setting of boundary condition and receiving input value (4) Process for calculating the node temperature by the finite element method, the difference method, etc. (5) Process for calculating the equivalent thermal conductivity The
図2のフローチャートに基づいて一実施例を説明する。 An embodiment will be described based on the flowchart of FIG.
ここでは、図3のようなPGA(Pin Grid Array)のピンを例に説明する。
等価熱伝導率を算出するため、ピン数と基板サイズから図4のような等価な計算領域を設定する。計算領域はピンの占める割合より求める。また今回はピン数9としたが、違うピン数の場合でもよい。
Here, a description will be given taking as an example a pin of a PGA (Pin Grid Array) as shown in FIG.
In order to calculate the equivalent thermal conductivity, an equivalent calculation region as shown in FIG. 4 is set from the number of pins and the substrate size. The calculation area is obtained from the ratio of pins. Although the number of pins is 9 this time, a different number of pins may be used.
計算モデルは熱伝導率の測定で使用される絶対法、比較法の2通りの方法を参考に、解析モデルを作成し等価熱伝導率を算出する。 The calculation model is based on the absolute method and the comparison method used in the measurement of thermal conductivity, and an analytical model is created to calculate the equivalent thermal conductivity.
まず絶対法を参考にした等価熱伝導率算出方法について説明する。 First, an equivalent thermal conductivity calculation method referring to the absolute method will be described.
図4のようなモデル形状、ピン、空気および樹脂の熱伝導率を入力する(ステップ1)。 The model shape as shown in FIG. 4, pins, air, and the thermal conductivity of the resin are input (step 1).
次に、計算領域をメッシュ分割する(ステップ2)。このメッシュ分割はピンが円筒形状であることを利用して、ピンの径に応じて計算精度が上がる様に自動的に分割する。形状により平面方向の分割だけで厚み方向が近似計算で十分で厚み方向の分割が必要でない場合もあるし、逆に厚み方向の分割だけで平面方向が近似計算で十分で平面方向の分割が必要でない場合もあるが、一般的には両方の分割が必要である。また、分割はピンが円筒形状であるため、境界部分で細かく分割するのが通常である。 Next, the calculation area is divided into meshes (step 2). This mesh division utilizes the fact that the pins are cylindrical, and automatically divides them so that the calculation accuracy increases according to the diameter of the pins. Depending on the shape, the thickness direction can be approximated only by dividing in the plane direction, and division in the thickness direction may not be necessary. Conversely, the division in the thickness direction is sufficient only by dividing in the thickness direction, and division in the plane direction is necessary. In some cases, both divisions are generally required. Further, since the pins are cylindrical, the division is usually finely divided at the boundary portion.
図5のような境界条件を設定する。ここでqは熱流量、t2は温度である(ステップ3)。これらの入力方法としては、ステップ1と同様に画面に数値の入力を促す枠を表示し、その枠にカーソルが移動したときに入力する様にするのが一般的である。なお、l2は厚みであり、Sは面積である。
A boundary condition as shown in FIG. 5 is set. Here, q is the heat flow rate, and t2 is the temperature (step 3). As these input methods, it is common to display a frame for prompting numerical input on the screen in the same manner as in
次に、有限要素法、差分法等により熱解析、若しくは熱流体解析を実施し、節点温度を求める(ステップ4)。このとき節点の温度は適宜間隔にて求める。この間隔はモデル形
状により適宜精度が一定以上になる様な間隔である。
Next, thermal analysis or thermal fluid analysis is performed by a finite element method, a difference method, or the like to obtain a node temperature (step 4). At this time, the temperature of the node is obtained at appropriate intervals. This interval is such that the accuracy becomes a certain level or more as appropriate depending on the model shape.
次式より等価熱伝導率λを算出する(ステップ5)。
λ=−(q/S){l2/(t1−t2)}
ここでt1は計算より求めた節点における平均温度である(図6)。
The equivalent thermal conductivity λ is calculated from the following equation (step 5).
λ = − (q / S) {l2 / (t1-t2)}
Here, t1 is the average temperature at the node obtained by calculation (FIG. 6).
次に、比較法を参考にした計算方法について述べる。なお、特段説明がないものは実施例1と同様である。図8のようなモデル作成をする。また標準試料の熱伝導率をλ1とする。
図9のように計算領域をメッシュ分割する(ステップ22)。
図8のように境界条件を設定する(ステップ32)。
有限要素法,差分法等により熱解析、若しくは熱流体解析を実施し、節点温度を求める(ステップ42)。
Next, a calculation method with reference to the comparison method will be described. Note that those not specifically described are the same as those in the first embodiment. A model is created as shown in FIG. The thermal conductivity of the standard sample is λ1.
The calculation area is divided into meshes as shown in FIG. 9 (step 22).
A boundary condition is set as shown in FIG. 8 (step 32).
Thermal analysis or thermal fluid analysis is performed by a finite element method, a difference method, or the like to obtain a node temperature (step 42).
そして次式より等価熱伝導率λを算出する(ステップ52)。なおl1は標準試料の厚みであり、t3は下面温度である。
λ=λ1{(t1−t)/(t−t3)}(l2/l1)
ここで温度tは、ステップ42で求めた標準試料との境界における節点温度の平均値である。
Then, the equivalent thermal conductivity λ is calculated from the following equation (step 52). Here, l1 is the thickness of the standard sample, and t3 is the lower surface temperature.
[lambda] = [lambda] 1 {(t1-t) / (t-t3)} (l2 / l1)
Here, the temperature t is an average value of the nodal temperatures at the boundary with the standard sample obtained in step 42.
本発明は、はんだバンプ、ピン、ビアなどがある様な、複合材の等価熱伝導率算出装置および等価熱伝導率算出方法に利用可能である。 The present invention can be used for an equivalent thermal conductivity calculation device and an equivalent thermal conductivity calculation method for composite materials such as solder bumps, pins, and vias.
1 入力装置
2 CPU
3 ディスプレイ
4 プリンタ
5 記憶装置
1
3
Claims (2)
寸法、形状、物性値の特定により条件決定可能なモデルをそなえており、
寸法、形状、物性値の入力を促す物性入力手段と、
入力された形状をモデルに適合した上で要素分割する要素分割手段と、
境界条件の入力を促す境界条件入力手段と、
節点の温度を算出する節点温度決定手段と、
節点の温度に基づき等価熱伝導率を算出する熱伝導計算手段と、
を具備することを特徴とする複合材の等価熱伝導率算出装置。 A composite material having input means for inputting conditions and output means for outputting equivalent thermal conductivity that can be regarded as equivalent thermal conductivity when heat conduction from the conduction boundary to the conduction boundary is configured as a single member. In the equivalent thermal conductivity calculation device,
We have a model that can determine conditions by specifying dimensions, shapes and physical properties.
Physical property input means for prompting input of dimensions, shapes, physical property values;
An element dividing means for dividing an input shape in conformity with a model,
Boundary condition input means for prompting input of boundary conditions;
A node temperature determining means for calculating the temperature of the node;
Heat conduction calculation means for calculating the equivalent thermal conductivity based on the temperature of the node;
A device for calculating an equivalent thermal conductivity of a composite material.
寸法、形状、物性値の入力を受け付ける物性入力処理と、
入力された形状をモデルに適合する様に要素分割する要素分割処理と、
境界条件の入力を受け付ける境界条件入力処理と、
節点の温度を演算処理する節点温度決定処理と、
節点の温度に基づき等価熱伝導率を演算する熱伝導計算処理と、
を具備することを特徴とする複合材の等価熱伝導率算出方法。 In the equivalent thermal conductivity calculation method to obtain the equivalent thermal conductivity that can be regarded as equivalent thermal conductivity when the heat conduction from the conduction boundary to the conduction boundary is configured as a single member,
Physical property input processing that accepts input of dimensions, shapes, physical property values,
Element division processing to divide the input shape to fit the model,
Boundary condition input processing for receiving boundary condition input;
A node temperature determination process for computing the temperature of the node;
Heat conduction calculation processing for calculating the equivalent thermal conductivity based on the temperature of the node;
A method for calculating an equivalent thermal conductivity of a composite material, comprising:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003374010A JP2005140509A (en) | 2003-11-04 | 2003-11-04 | Equivalent thermal conductivity calculation device and equivalent thermal conductivity calculation method for composite material |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003374010A JP2005140509A (en) | 2003-11-04 | 2003-11-04 | Equivalent thermal conductivity calculation device and equivalent thermal conductivity calculation method for composite material |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005140509A true JP2005140509A (en) | 2005-06-02 |
Family
ID=34685867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003374010A Pending JP2005140509A (en) | 2003-11-04 | 2003-11-04 | Equivalent thermal conductivity calculation device and equivalent thermal conductivity calculation method for composite material |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2005140509A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009282020A (en) * | 2008-04-23 | 2009-12-03 | Central Res Inst Of Electric Power Ind | Method, apparatus and program for evaluating heat shielding performance of coating layer |
JP2010008403A (en) * | 2008-05-27 | 2010-01-14 | Central Res Inst Of Electric Power Ind | Method, apparatus and program for evaluating heat shield performance of coating layer |
JP2011164032A (en) * | 2010-02-12 | 2011-08-25 | Eiji Nemoto | Method for measuring principal-axis thermal constant of two-dimensional anisotropic heat conductor using multi-point temperature measurement by pulse/period method, and measurement apparatus thereof |
CN110362935A (en) * | 2019-07-18 | 2019-10-22 | 济南大学 | A kind of modeling of plastic-steel door and window frame heat transfer and calculation method |
-
2003
- 2003-11-04 JP JP2003374010A patent/JP2005140509A/en active Pending
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009282020A (en) * | 2008-04-23 | 2009-12-03 | Central Res Inst Of Electric Power Ind | Method, apparatus and program for evaluating heat shielding performance of coating layer |
JP2010008403A (en) * | 2008-05-27 | 2010-01-14 | Central Res Inst Of Electric Power Ind | Method, apparatus and program for evaluating heat shield performance of coating layer |
JP2011164032A (en) * | 2010-02-12 | 2011-08-25 | Eiji Nemoto | Method for measuring principal-axis thermal constant of two-dimensional anisotropic heat conductor using multi-point temperature measurement by pulse/period method, and measurement apparatus thereof |
CN110362935A (en) * | 2019-07-18 | 2019-10-22 | 济南大学 | A kind of modeling of plastic-steel door and window frame heat transfer and calculation method |
CN110362935B (en) * | 2019-07-18 | 2022-09-09 | 济南大学 | Heat conduction modeling and calculating method for plastic-steel door and window frame |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101935148B1 (en) | Method for estimating a variation in preload applied to linear guideway | |
CN109002632B (en) | Liquid cooling system simulation method and device | |
JPWO2016056642A1 (en) | Blood flow analysis device for blood flow simulation, method and computer software program | |
JP2006343311A (en) | Three-dimensional measurement data inspection method using parametric tolerance | |
JP2005140509A (en) | Equivalent thermal conductivity calculation device and equivalent thermal conductivity calculation method for composite material | |
JP2007316032A (en) | Analytical data producing device, analytical data producing method, analytical data producing program, and analyzer | |
Gong et al. | Determining local residual stresses from high resolution wafer geometry measurements | |
JP2013054657A (en) | Information processing apparatus and method thereof | |
JP2006350504A (en) | Analytic result display method, analytic result display program and analytic result display device | |
US11526637B2 (en) | Information processing apparatus, information processing method, and storage medium for creating a thermal network model in a short time | |
JP2014203195A (en) | Design support device | |
Harari et al. | A note on a recent study of stabilized finite element computations for heat conduction | |
Frey et al. | An adaptive numerical scheme for solving incompressible 2‐phase and free‐surface flows | |
JP2009020625A (en) | Analysis model creation method, program, recording medium and analysis model creation device | |
JP3023969B2 (en) | Method for analyzing temperature of cooling / heating cycle structure and design apparatus for mold apparatus system | |
JP5235573B2 (en) | Strength analysis method, strength analysis apparatus, and strength analysis program | |
JP2008052380A (en) | Element split system of finite element method analysis model, element split apparatus, element split method, and its program | |
Ely et al. | A comparative study of metamodeling techniques for predictive process control of welding applications | |
Calderon et al. | Approximating nonequilibrium processes using a collection of surrogate diffusion models | |
KR102100113B1 (en) | Pressure sensor and calibration method pressure distribution measuring device | |
JP7247845B2 (en) | Heat transfer analyzer | |
JP2006215711A (en) | Method and device for creating mesh from three-dimensional cad model | |
JPH1185812A (en) | Method and system for analyzing curved state of structure using finite element method, and record medium recorded with curved state analysis program of the structure | |
JP6826919B2 (en) | Data distribution device and data distribution ratio determination method | |
JP2007108020A (en) | Inflow calculation program, apparatus and method therefor |