Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP2846764B2 - Thermophysical property measurement method - Google Patents
[go: Go Back, main page]

JP2846764B2 - Thermophysical property measurement method - Google Patents

Thermophysical property measurement method

Info

Publication number
JP2846764B2
JP2846764B2 JP9796492A JP9796492A JP2846764B2 JP 2846764 B2 JP2846764 B2 JP 2846764B2 JP 9796492 A JP9796492 A JP 9796492A JP 9796492 A JP9796492 A JP 9796492A JP 2846764 B2 JP2846764 B2 JP 2846764B2
Authority
JP
Japan
Prior art keywords
thermophysical
thermal expansion
sample
equation
constants
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 - Fee Related
Application number
JP9796492A
Other languages
Japanese (ja)
Other versions
JPH05296953A (en
Inventor
伸吾 住江
弘行 高松
善郎 西元
博文 今中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kobe Steel Ltd
Original Assignee
Kobe Steel Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kobe Steel Ltd filed Critical Kobe Steel Ltd
Priority to JP9796492A priority Critical patent/JP2846764B2/en
Publication of JPH05296953A publication Critical patent/JPH05296953A/en
Application granted granted Critical
Publication of JP2846764B2 publication Critical patent/JP2846764B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【産業上の利用分野】本発明は,試料の体積比熱(ρ
C),熱伝導率(K),熱膨張率(α)の3つの熱物性
定数の内の1つが既知である場合に,他の2つの未知の
熱物性定数を,1つの測定原理を用いて上記1つの既知
の熱物性定数から求める方法に関するものである。
BACKGROUND OF THE INVENTION The present invention relates to a sample having a specific heat of volume (ρ
If one of the three thermophysical constants of C), thermal conductivity (K), and thermal expansion coefficient (α) is known, the other two unknown thermophysical constants are calculated using one measurement principle. And a method for obtaining the same from the one known thermophysical property constant.

【0002】[0002]

【従来の技術】上記したような試料の熱体積比熱,熱伝
導率,熱膨張率は,いずれも熱設計上重要な物理定数で
あり,従来,以下に示すような方法で測定されている。
熱伝導率(K)及び比熱(C)は,例えばレーザフラッ
シュ式熱定数測定法により測定されている。これは,所
定の温度に保持された試料の前面にレーザ光を瞬間的に
照射して加熱し,その瞬間からの試料の裏面の温度変化
を熱伝対や赤外線検出器を用いて測定するものである。
測定装置の概要が,図7に示されている。測定された試
料の裏面の温度は,図8に示すように変化する。上記の
ような温度変化曲線と一次元熱伝導方程式の解である
(1)式により熱拡散率Dが求められ,この熱拡散率D
と予め知られている試料の重量M,密度ρから(2)式
を用いて比熱C及び熱伝導率Kが算出される。 D=1.37(L/π)2 /(t)0.5 …(1) C=Q/(M・ΔT),K=C・D・ρ …(2) ここに,tは測定時間,Lは試料の厚さ,ΔTは温度上
昇,Qは吸収エネルギーである。また,熱膨張率αは,
図9に示すような天秤をてことして使用し,試料を加熱
することによって生じる検出棒の伸縮を差動トランスで
検出する,いわゆる微小定荷重熱膨張方式で測定され
る。
2. Description of the Related Art Thermal volume specific heat, thermal conductivity, and coefficient of thermal expansion of a sample as described above are all important physical constants in thermal design, and are conventionally measured by the following methods.
The thermal conductivity (K) and the specific heat (C) are measured by, for example, a laser flash thermal constant measuring method. This is a method in which a laser beam is instantaneously irradiated and heated on the front surface of a sample held at a predetermined temperature, and the temperature change on the back surface of the sample from that moment is measured using a thermocouple or infrared detector. It is.
The outline of the measuring device is shown in FIG. The measured temperature of the back surface of the sample changes as shown in FIG. The thermal diffusivity D is obtained from the above-mentioned temperature change curve and equation (1) which is a solution of the one-dimensional heat conduction equation.
The specific heat C and the thermal conductivity K are calculated from the weight M and the density ρ of the sample, which are known in advance, using the equation (2). D = 1.37 (L / π) 2 / (t) 0.5 (1) C = Q / (M · ΔT), K = C · D · ρ (2) where t is the measurement time and L Is the thickness of the sample, ΔT is the temperature rise, and Q is the absorbed energy. The coefficient of thermal expansion α is
It is measured by a so-called minute constant load thermal expansion method in which a balance as shown in FIG. 9 is used and the expansion and contraction of the detection rod caused by heating the sample is detected by a differential transformer.

【0003】[0003]

【発明が解決しようとする課題】上記のように,熱伝導
率K,熱膨張率αは,共に熱設計上重要な物理定数であ
るにもかかわらず,前記のような従来技術では別々の測
定装置により,即ち別々の測定原理を用いて測定されて
おり,1つの測定原理の中で同時に得ることはできな
い。そのため,上記熱定数間の一貫性に乏しく,測定精
度についてそれぞれまちまちのデータを使用せざるを得
ない。上記のように,熱伝導率K及び熱膨張率αの測定
原理が別個であるため,それぞれの測定原理に適した形
状に試料を別途加工する必要があり,試料の加工誤差が
測定データのばらつきを促進させると共に,加工コスト
を上昇させる。また,前記したような測定原理では,半
導体ウエハのような薄い材料について測定することが困
難である。また,前記の熱伝導率K及び熱膨張率αの測
定方法では,加工した試料全体の平均値を得るのみであ
って,試料中にK,αの変動があってもそれらを細かく
測定することができない。
As described above, although the thermal conductivity K and the thermal expansion coefficient α are both important physical constants in the thermal design, they are measured separately in the above-mentioned prior art. It is measured by the device, ie using different measurement principles, and cannot be obtained simultaneously in one measurement principle. Therefore, the consistency between the thermal constants is poor, and different data must be used for the measurement accuracy. As described above, since the measurement principles of the thermal conductivity K and the coefficient of thermal expansion α are different, it is necessary to separately process the sample into a shape suitable for each measurement principle, and the processing error of the sample may cause a variation in the measured data. As well as increase processing costs. Further, it is difficult to measure a thin material such as a semiconductor wafer by the above-described measurement principle. In the method of measuring the thermal conductivity K and the coefficient of thermal expansion α described above, only the average value of the entire processed sample is obtained, and even if K and α fluctuate in the sample, they must be measured finely. Can not.

【0004】[0004]

【課題を解決するための手段】上記したような課題を解
決するために,本発明は,試料の体積比熱(ρC),熱
伝導率(K),熱膨張率(α)の3つの熱物性定数の内
の1つが既知である場合に,他の2つの熱物性定数を,
1つの測定原理を用いて上記1つの既知の熱物性定数か
ら求める方法において,上記試料に強度変調した励起光
を断続照射し,その時試料表面に生じた熱膨張振動の振
幅を測定する工程を,異なる変調周波数の励起光を用い
て少なくとも3回繰り返し,得られた3組以上の変調周
波数と熱膨張振幅の関係を前記体積比熱,熱伝導率,熱
膨張率を未知数とする1つの関係式に代入して既知の熱
物性特性と未知の熱物性特性の関係式を求め,上記既知
の熱物性定数から未知の熱物性定数を求めることを特徴
とする熱物性定数測定方法として構成されている。
In order to solve the above-mentioned problems, the present invention provides three thermophysical properties of a sample such as volume specific heat (ρC), thermal conductivity (K), and thermal expansion coefficient (α). If one of the constants is known, the other two thermophysical constants are
In the method of obtaining from the one known thermophysical property constant using one measurement principle, a step of intermittently irradiating the sample with intensity-modulated excitation light and measuring the amplitude of thermal expansion vibration generated on the sample surface at that time includes the steps of: It is repeated at least three times using pump lights of different modulation frequencies, and the relationship between the obtained three or more sets of modulation frequencies and the thermal expansion amplitude is expressed as one relational expression where the volume specific heat, thermal conductivity, and thermal expansion coefficient are unknown. The method is configured as a method for measuring a thermophysical constant characterized by determining a relational expression between a known thermophysical property and an unknown thermophysical property by substituting, and obtaining an unknown thermophysical constant from the known thermophysical constant.

【0005】[0005]

【作用】本発明における測定原理を,図4を用いて説明
する。図4に示すような試料に,励起光を変調角周波数
ωで強度変調して,試料に垂直に照射する。上記励起光
の半径をr,パワーをP0 とする。また,試料の熱伝導
率をK,熱膨張率をα,密度をρ,比熱をCとする。図
5は,上記のような励起光の強度変調状態を示す。上記
のような励起光は,図4に示すように,試料の複素屈折
率で決まる深さdだけ試料に侵入する。侵入状態を斜線
で示す。そこから発生した熱は,熱拡散長μ(=(2K
/ωρC)0.5 )の距離に伝播して体積Aの熱拡散領域
の温度をΔTだけ上昇させる。1回の照射パルス(図5
の網かけ部)のエネルギーはπP0/ω,熱拡散領域の
熱容量はAρCであるから, ΔT=πP0 /(AρC・ω) …(3) また,熱膨張量(図4の2PA,PAは熱膨張振動の振
幅を表す)は, 2PA=α・ΔT・(d+μ) …(4) また,図4より, A=π(r+μ)2 (d+μ) …(5) 上記(3),(5)式を(4)式に代入して整理する
と,(6)式を得る。 2PA=αP0 /(2K+r(8ρCK・ω)0.5 +ρCr2 ω) …(6) (6)式の正当性を評価するため,図1に示す実施例の
装置により予めα,C,ρ,Kの判っている銅系の試料
4種類を測定し,上記(6)式から得られる理論値と比
較した。図6はその結果を示すもので,理論値と実測値
が一致していることを示している。これにより,上記
(6)式の正当性が理解される。上記(6)式における
K,α,ρCを未知数として,これらの関係式を求める
ために,1/(2PA)=y,ω=x2 として(6)式
を変形すると,(7)式のようになる。 y=2K/(αP0 )+(r(8ρCK)0.5 /(αP0 ))x +(ρCr2 /(αP0 ))x2 …(7) 上記(7)式は,上記K,α,ρCの3つの未知数を含
むから,試料に変調周波数ω(x=ω2 )で強度変調し
た励起光を,図5に示す如く断続照射し,その時試料表
面に生じた熱膨張振動の振幅(2PA)(2PA=1/
y)を測定する工程を異なる変調周波数の励起光を用い
て少なくとも3回繰り返し,3組のxとyとの関係を求
めることによって,各係数を決定することができる。そ
の結果,予め判っているP0 やr等も上記決定された係
数のなかに含めると, K/α=X …(8) (ρC・K)0.5 /α=Y …(9) ρC/α=Z …(10) のように書き表される。上記X,Y,Zは,上記のよう
な3回の測定により得られた変調周波数と熱膨張量とを
(7)式に代入することによって演算された係数(定
数)である。上記(8)式及び(10)式より, α=ρC/Z …(11) K=(X/Z)・ρC …(12) が得られる。ここで,例えば上記ρCを予め計測してお
けば,上記(11),(12)式よりαとKを同時に求
めることができる。上記の説明は,ρCを既知として残
りの未知数αとKを(8),(9),(10)式より求
めたが,同様に,例えばαを既知としてρC及びKを求
めたり,Kを既知としてρC及びαを求めることもでき
る。上記のように,この方法によれば,試料の体積比熱
(ρC),熱伝導率(K),熱膨張率(α)の3つの熱
物性定数の内の1つが既知である場合に,他の2つの熱
物性定数を(7)式で表される1つの測定原理を用いて
上記3つの熱物性定数の内の1つから求める方法に係
り,上記のように試料に強度変調した励起光を断続照射
し,その時試料表面に生じた振幅(2PA)を変調周波
数(ω)との組合せで測定する測定工程を異なる変調周
波数の励起光を用いて少なくとも3回繰り返し,得られ
た3組以上の変調周波数と熱膨張振幅の関係を前記体積
比熱,熱伝導率,熱膨張率を未知数とする1つの関係式
((7)式)に代入して既知の熱物性特性と未知の熱物
性特性の関係式((8),(9),(10)式)を求
め,これらの関係式に上記既知の熱物性定数を代入する
ことにより,これらの既知の熱物性定数から未知の熱物
性定数を求めるものである。
The principle of measurement according to the present invention will be described with reference to FIG. Excitation light is intensity-modulated at a modulation angular frequency ω onto a sample as shown in FIG. The radius of the excitation light is r, and the power is P 0 . Further, K is the thermal conductivity of the sample, α is the coefficient of thermal expansion, ρ is the density, and C is the specific heat. FIG. 5 shows the intensity modulation state of the excitation light as described above. The excitation light as described above penetrates into the sample by a depth d determined by the complex refractive index of the sample, as shown in FIG. The intrusion state is indicated by oblique lines. The heat generated therefrom has a heat diffusion length μ (= (2K
/ ΩρC) 0.5 ) to increase the temperature of the heat diffusion region of the volume A by ΔT. One irradiation pulse (Fig. 5
Since the energy of the hatched portion is πP 0 / ω and the heat capacity of the heat diffusion region is AρC, ΔT = πP 0 / (AρC · ω) (3) Also, the amount of thermal expansion (2PA, PA in FIG. 4) Represents the amplitude of the thermal expansion vibration) 2PA = ααΔT ・ (d + μ) (4) From FIG. 4, A = π (r + μ) 2 (d + μ) ... (5) Substituting equation (5) into equation (4) and rearranging yields equation (6). 2PA = αP 0 / (2K + r (8ρCK · ω) 0.5 + ρCr 2 ω) (6) In order to evaluate the validity of the equation (6), α, C, ρ, K Four types of copper-based samples known to have been measured were compared with theoretical values obtained from the above equation (6). FIG. 6 shows the result, and shows that the theoretical value and the measured value match. Thus, the validity of the above equation (6) is understood. If K, α, and ρC in the above equation (6) are set as unknowns, and equation (6) is modified as 1 / (2PA) = y, ω = x 2 to obtain these relational equations, Become like y = 2K / (αP 0 ) + (r (8ρCK) 0.5 / (αP 0 )) x + (ρCr 2 / (αP 0 )) x 2 (7) The above equation (7) expresses the above K, α, Since the three unknowns of ρC are included, the sample is irradiated intermittently with excitation light whose intensity is modulated at the modulation frequency ω (x = ω 2 ) as shown in FIG. 5, and the amplitude of the thermal expansion vibration (2PA ) (2PA = 1 /
Each coefficient can be determined by repeating the step of measuring y) at least three times using pump lights having different modulation frequencies and determining the relationship between three sets of x and y. As a result, if P 0 , r, etc., which are known in advance, are also included in the determined coefficients, K / α = X (8) (ρC · K) 0.5 / α = Y (9) ρC / α = Z (10) X, Y, and Z are coefficients (constants) calculated by substituting the modulation frequency and the amount of thermal expansion obtained by the above three measurements into the equation (7). From the above equations (8) and (10), α = ρC / Z (11) K = (X / Z) · ρC (12) is obtained. Here, for example, if ρC is measured in advance, α and K can be obtained simultaneously from the above equations (11) and (12). In the above description, ρC is known and the remaining unknowns α and K are obtained from equations (8), (9), and (10). Similarly, for example, ρC and K are obtained with α already known, It is also possible to obtain ρC and α as known. As described above, according to this method, when one of the three thermophysical constants, that is, the volume specific heat (ρC), the thermal conductivity (K), and the thermal expansion coefficient (α) of the sample is known, The method of obtaining the two thermophysical constants from one of the above three thermophysical constants using one measurement principle expressed by the equation (7), wherein the excitation light intensity-modulated on the sample as described above Is intermittently illuminated, and the measurement process of measuring the amplitude (2PA) generated on the sample surface at that time in combination with the modulation frequency (ω) is repeated at least three times using the excitation light having different modulation frequencies. The known thermophysical properties and unknown thermophysical properties are obtained by substituting the relationship between the modulation frequency and the thermal expansion amplitude into one relational expression (Equation (7)) where the volume specific heat, thermal conductivity, and thermal expansion coefficient are unknowns. Are obtained (Equations (8), (9), and (10)), and By substituting the thermal property constants, and requests unknown thermal properties constant from these known thermal physical constants.

【0006】[0006]

【実施例】続いて,図1に示した装置を用いて,変調周
波数を変えて3回の測定を行い,これらの測定により得
られたデータを用いて,3つの熱物性定数の内の1つの
既知の熱物性定数から残りの2つの未知の熱物性定数を
求めた実施例につき説明し,本発明の理解に供する。図
1に示した装置の場合,半導体レーザ10は,発振器1
2とドライバ11により,30〜200KHzの変調周
波数で強度変調されたレーザ光を,破線で示すように照
射する。レーザ光は,ダイクロイックミラー5,レンズ
15により,試料16上に直径数ミクロンのスポット光
を照射する(図4)。このレーザ光照射により引き起こ
された試料表面の熱膨張振動は,ヘリウムネオンレーザ
1,偏光ビームスプリッタ2,周波数シフタ18,参照
ミラー3,4分の1波長板4,偏光子6によって構成さ
れるヘテロダイン型の変移計測器で計測される。図1中
のF1 はキャリア周波数,Fb はビート周波数である。
上記ミラー3で反射された光と試料16で反射された光
は,偏光ビームスプリッタ2で合成され,その干渉光は
光検出器7で検出され,光電変換される。こうして電気
信号に置き換えられた変移信号は,バンドパスフィルタ
8,低域通過フィルタ9,乗算器13,ロックインアン
プ14からなる信号処理回路によってノイズを除去され
た後,コンピュータ17で熱膨張振動の振幅(PA)に
変換される。図2は,上記図1に示した測定装置を用い
て,厚さ0.5mmのシリコンウエハを試料として変調周
波数(ω/2π)を6種類に変化させた時の各熱膨張振
動の振幅(PA)の実測値を表したもので,これらの値
は前記(6)式から演算される理論値と合致する。図3
は,上記図2のデータを,前記(6)式を変形するため
に用いたy=1/2PA,x=(ω)0.5 の式によるx
とyの関係に書き改めたものである。こうして得られた
図3の3つの測定値を(7)式に代入し,これら3つの
連立方程式を解くことにより,(8),(10)式を用
いて,係数X及びZを予め判っている値P0 (励起光の
パワー),r(スポット半径)から演算すると, X=7.0×107 〔W/m〕,Z=6.9×1011〔J/m3 〕 が得られた。
EXAMPLE Next, three measurements were performed using the apparatus shown in FIG. 1 while changing the modulation frequency, and one of the three thermophysical constants was determined using the data obtained from these measurements. An example in which the remaining two unknown thermophysical property constants are obtained from one known thermophysical property constant will be described to provide an understanding of the present invention. In the case of the device shown in FIG.
2 and the driver 11 irradiate a laser beam intensity-modulated at a modulation frequency of 30 to 200 KHz as shown by a broken line. The laser beam irradiates a spot light having a diameter of several microns on the sample 16 by the dichroic mirror 5 and the lens 15 (FIG. 4). The thermal expansion vibration of the sample surface caused by this laser beam irradiation is caused by a heterodyne composed of a helium neon laser 1, a polarization beam splitter 2, a frequency shifter 18, a reference mirror 3, a quarter-wave plate 4, and a polarizer 6. It is measured with a type displacement measuring instrument. F 1 in FIG. 1 is the carrier frequency, F b is the beat frequency.
The light reflected by the mirror 3 and the light reflected by the sample 16 are combined by the polarization beam splitter 2, and the interference light is detected by the photodetector 7 and photoelectrically converted. The transition signal thus replaced with the electric signal is subjected to a signal processing circuit including a band-pass filter 8, a low-pass filter 9, a multiplier 13, and a lock-in amplifier 14 to remove noise therefrom. It is converted to amplitude (PA). FIG. 2 shows the amplitudes of the thermal expansion vibrations when the modulation frequency (ω / 2π) is changed to six types using a 0.5 mm thick silicon wafer as a sample using the measuring apparatus shown in FIG. PA), and these values agree with the theoretical values calculated from the above equation (6). FIG.
Is obtained by transforming the data of FIG. 2 into an equation of y = 1 / PA and x = (ω) 0.5 used to transform the equation (6).
And y. By substituting the three measured values in FIG. 3 obtained in this way into equation (7) and solving these three simultaneous equations, the coefficients X and Z are determined in advance using equations (8) and (10). When calculated from the values P 0 (power of the excitation light) and r (spot radius), X = 7.0 × 10 7 [W / m] and Z = 6.9 × 10 11 [J / m 3 ] are obtained. Was done.

【0007】これらの値を(8)及び(10)式に代入
し,試料であるシリコンウエハについての密度ρ=23
30〔kg/m3 〕,比熱C=772.4〔J/kg・
K〕により体積比熱ρC=1.8×106 〔J/m
3 K〕となるので,(11)式より, α=(ρC)/Z=(1.8×106 )/(6.9×1011)=2.6 ×10-6〔/K〕 K=(X/Z)・ρC=((7.0×107 )/(6.9×1011)) ×(1.8×106 )=182〔W/mK〕 が得られる。この値は文献値(α=2.5×10-6〔/
K〕,K=170〔W/mK〕)とよく合致している。
上記したように,この実施例では,試料の体積比熱(ρ
C),熱伝導率(K),熱膨張率(α)の3つの熱物性
定数の内の1つ(体積比熱)が既知であり,他の2つの
熱物性定数(熱伝導率,熱膨張率)が未知である場合に
ついて示したが,熱膨張率(α)が既知であれば,αの
値を(8)式に代入して熱伝導率(K)を求めると共
に,(10)式にαを代入して体積比熱(ρC)を求め
ることも同様の手法により可能である。更に,熱伝導率
(K)が既知であれば,(8)式にKを代入して熱膨張
率(α)を求めると共に,(9)式に既に求めたαとK
とを代入してρCを求めるか,上記(8)式から求めた
αを(10)式に代入してρCを求めることも可能であ
る。要するに,上記実施例における(7)式に,測定に
より求められた強度変調の変調周波数とそれに対応する
熱膨張振動の振幅を代入し,既知の熱物性定数と未知の
熱物性定数の関係式((8),(9),(10)式)を
求めて,これらの式に既知の熱物性定数を代入すること
により,残りの未知の熱物性定数を求めることができる
ものである。なお,上の例では(7)式から3つの未知
数間の関係を得るために3回の測定結果を(7)式に代
入し各係数である定数を求めているが,変調周波数を数
多く変化させてPAとω(即ち,yとx)の関係を求
め,それらより最小二乗法を用いて各係数を算出するよ
うにしてもよい。これにより,更に精度が向上する。
By substituting these values into equations (8) and (10), the density ρ = 23
30 [kg / m 3 ], specific heat C = 772.4 [J / kg ·
K], the specific heat of volume ρC = 1.8 × 10 6 [J / m
Since the 3 K], (11) from the equation, α = (ρC) / Z = (1.8 × 10 6) / (6.9 × 10 11) = 2.6 × 10 -6 [/ K] K = (X / Z) · ρC = ((7.0 × 10 7 ) / (6.9 × 10 11 )) × (1.8 × 10 6 ) = 182 [W / mK] This value is based on the literature value (α = 2.5 × 10 −6 [/
K], K = 170 [W / mK]).
As described above, in this embodiment, the volume specific heat of the sample (ρ
C), thermal conductivity (K), and thermal expansion coefficient (α), one of the three thermophysical constants (specific heat of volume) is known, and the other two thermophysical constants (thermal conductivity, thermal expansion) The coefficient of thermal expansion (α) is unknown, but if the coefficient of thermal expansion (α) is known, the value of α is substituted into equation (8) to determine the thermal conductivity (K), and the equation (10) Can be obtained by the same method by substituting α into the above equation to obtain the volume specific heat (ρC). Further, if the thermal conductivity (K) is known, the thermal expansion coefficient (α) is obtained by substituting K into the equation (8), and α and K already obtained in the equation (9).
Can be substituted to obtain ρC, or α obtained from the above equation (8) can be substituted into the equation (10) to determine ρC. In short, the modulation frequency of the intensity modulation obtained by the measurement and the amplitude of the thermal expansion vibration corresponding thereto are substituted for the equation (7) in the above embodiment, and the relational equation between the known thermophysical property constant and the unknown thermophysical property constant ( (8), (9), and (10) are obtained, and the remaining unknown thermophysical constants can be obtained by substituting the known thermophysical constants into these equations. In the above example, in order to obtain the relationship between the three unknowns from equation (7), the results of the three measurements are substituted into equation (7) to determine the constants as the respective coefficients. Then, the relationship between PA and ω (that is, y and x) may be obtained, and each coefficient may be calculated from them using the least squares method. Thereby, the accuracy is further improved.

【0008】[0008]

【発明の効果】本発明は,以上述べたように構成されて
いるので,試料の体積比熱(ρC),熱伝導率(K),
熱膨張率(α)の3つの内の1つが判れば,他の2つを
1つの測定原理に基づいて得ることができるので,得ら
れた熱物性定数間のばらつきの程度が均質化し,データ
の信頼性が向上する。
The present invention is constructed as described above, so that the sample has a specific volumetric heat (ρC), thermal conductivity (K),
If one of the three coefficients of thermal expansion (α) is known, the other two can be obtained based on one measurement principle, so the degree of variation between the obtained thermophysical constants is homogenized, and the data Reliability is improved.

【図面の簡単な説明】[Brief description of the drawings]

【図1】 本発明を適用することのできる熱膨張振動測
定装置のブロック図。
FIG. 1 is a block diagram of a thermal expansion vibration measuring device to which the present invention can be applied.

【図2】 図1の測定装置で得られた変調周波数と熱膨
張振動量との関係を示すグラフ。
FIG. 2 is a graph showing a relationship between a modulation frequency and a thermal expansion vibration amount obtained by the measurement device of FIG. 1;

【図3】 図2を変形したグラフ。FIG. 3 is a graph obtained by modifying FIG. 2;

【図4】 本発明の理論を説明するための概念図。FIG. 4 is a conceptual diagram for explaining the theory of the present invention.

【図5】 図4に示された励起光の照射状態を示すグラ
フ。
FIG. 5 is a graph showing an irradiation state of the excitation light shown in FIG.

【図6】 (6)式の正当性を評価するための理論値と
実測値との関係を示すグラフ。
FIG. 6 is a graph showing a relationship between a theoretical value and an actually measured value for evaluating the validity of the equation (6).

【図7】 従来の熱伝導率と比熱を測定する装置の概要
図。
FIG. 7 is a schematic diagram of a conventional apparatus for measuring thermal conductivity and specific heat.

【図8】 図7の装置により得られる温度変化を示すグ
ラフ。
FIG. 8 is a graph showing a temperature change obtained by the apparatus of FIG. 7;

【図9】 従来の熱膨張率を測定する装置を示す原理
図。
FIG. 9 is a principle view showing a conventional apparatus for measuring a coefficient of thermal expansion.

【符号の説明】[Explanation of symbols]

10…半導体レーザ(励起光レーザ) 1…ヘリウムネオンレーザ(参照光源) 7…光
検出器
10 semiconductor laser (excitation light laser) 1 helium neon laser (reference light source) 7 photodetector

───────────────────────────────────────────────────── フロントページの続き (72)発明者 今中 博文 神戸市西区高塚台1丁目5−5株式会社 神戸製鋼所 神戸総合技術研究所内 (58)調査した分野(Int.Cl.6,DB名) G01N 25/00 - 25/72 JICSTファイル(JOIS)──────────────────────────────────────────────────続 き Continued on the front page (72) Inventor Hirofumi Imanaka 1-5-5 Takatsukadai, Nishi-ku, Kobe Kobe Steel, Ltd. Kobe Research Institute (58) Field surveyed (Int. Cl. 6 , DB name) G01N 25/00-25/72 JICST file (JOIS)

Claims (1)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 試料の体積比熱(ρC),熱伝導率
(K),熱膨張率(α)の3つの熱物性定数の内の1つ
が既知である場合に,他の2つの熱物性定数を,1つの
測定原理を用いて上記1つの既知の熱物性定数から求め
る方法において, 上記試料に強度変調した励起光を断続照射し,その時試
料表面に生じた熱膨張振動の振幅を測定する工程を,異
なる変調周波数の励起光を用いて少なくとも3回繰り返
し,得られた3組以上の変調周波数と熱膨張振幅の関係
を前記体積比熱,熱伝導率,熱膨張率を未知数とする1
つの関係式に代入して既知の熱物性特性と未知の熱物性
特性の関係式を求め,上記既知の熱物性定数から未知の
熱物性定数を求めることを特徴とする熱物性定数測定方
法。
1. When one of three thermophysical constants of a sample, ie, volume specific heat (ρC), thermal conductivity (K), and coefficient of thermal expansion (α), is known, the other two thermophysical constants are known. Irradiating the sample with an intensity-modulated excitation light intermittently and measuring the amplitude of thermal expansion vibration generated on the surface of the sample at that time. Is repeated at least three times using pump light of different modulation frequencies, and the relationship between the obtained three or more sets of modulation frequencies and the thermal expansion amplitude is defined as the volume specific heat, thermal conductivity, and thermal expansion coefficient as unknowns.
A method for measuring a thermophysical property constant, comprising: obtaining a relational equation between a known thermophysical property property and an unknown thermophysical property property by substituting into two relational equations; and obtaining an unknown thermophysical property constant from the known thermophysical property constant.
JP9796492A 1992-04-17 1992-04-17 Thermophysical property measurement method Expired - Fee Related JP2846764B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP9796492A JP2846764B2 (en) 1992-04-17 1992-04-17 Thermophysical property measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP9796492A JP2846764B2 (en) 1992-04-17 1992-04-17 Thermophysical property measurement method

Publications (2)

Publication Number Publication Date
JPH05296953A JPH05296953A (en) 1993-11-12
JP2846764B2 true JP2846764B2 (en) 1999-01-13

Family

ID=14206360

Family Applications (1)

Application Number Title Priority Date Filing Date
JP9796492A Expired - Fee Related JP2846764B2 (en) 1992-04-17 1992-04-17 Thermophysical property measurement method

Country Status (1)

Country Link
JP (1) JP2846764B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5086863B2 (en) * 2007-08-29 2012-11-28 株式会社神戸製鋼所 Thermophysical property evaluation equipment, measurement method for thermophysical property evaluation

Also Published As

Publication number Publication date
JPH05296953A (en) 1993-11-12

Similar Documents

Publication Publication Date Title
US12092595B2 (en) Steady-state thermo-reflectance method and system to measure thermal conductivity
US4632561A (en) Evaluation of surface and subsurface characteristics of a sample
JPH0347703B2 (en)
WO2015027210A1 (en) Lit method and system for determining material layer parameters of a sample
CN106796899A (en) Optical system and method for determining stress in a substrate
Rani et al. Estimation of defect depth in carbon fibre reinforced polymer using frequency modulated thermal wave imaging: an analytical study
CN117491427A (en) A multilayer wafer bubble defect tomography method based on transmission laser phase locking
Thiel et al. Localization of subsurface defects in uncoated aluminum with structured heating using high-power VCSEL laser arrays
JP2846764B2 (en) Thermophysical property measurement method
JP3029757B2 (en) Sample evaluation method by photothermal displacement measurement
JP2001074637A (en) Dynamic light scattering particle diameter distribution measuring system
Sodnik et al. Photothermal interferometry for nondestructive subsurface defect detection
JP3130708B2 (en) Sample defect evaluation method
CN111307059B (en) Photothermal surface deformation detection and calibration device and method based on wavelength phase-shifting interference
JP4101012B2 (en) Thermal property evaluation method and apparatus for laminated material having thermal resistance
JP2735348B2 (en) Sample evaluation method with a single light source using thermal expansion vibration
EP0772038A2 (en) Method and apparatus for optical measuring by polarization analysis
US3532973A (en) Microwave flaw detector
JP3040140B2 (en) Chromatic aberration measurement method and measurement device
JP3271994B2 (en) Dimension measurement method
JP3662370B2 (en) Optical measuring apparatus and method
RU2114421C1 (en) Method of bispectral pulse-frequency inspection
RU1786411C (en) Method for measuring temperature conductivity coefficient of solid bodies
JPH05226445A (en) Non-contact carrier diffusion length measuring device
Kim et al. Development on the reconstruction of photothermal imaging method for subsurface structure

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees