WO2015008632A1 - 熱拡散率測定装置 - Google Patents

熱拡散率測定装置 Download PDF

Info

Publication number
WO2015008632A1
WO2015008632A1 PCT/JP2014/067805 JP2014067805W WO2015008632A1 WO 2015008632 A1 WO2015008632 A1 WO 2015008632A1 JP 2014067805 W JP2014067805 W JP 2014067805W WO 2015008632 A1 WO2015008632 A1 WO 2015008632A1
Authority
WO
WIPO (PCT)
Prior art keywords
thermal diffusivity
sample
infrared light
laser beam
phase difference
Prior art date
Application number
PCT/JP2014/067805
Other languages
English (en)
French (fr)
Inventor
羽鳥仁人
大槻哲也
加藤岳雄
Original Assignee
株式会社ベテル
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 株式会社ベテル filed Critical 株式会社ベテル
Priority to US14/409,078 priority Critical patent/US20160033431A1/en
Priority to EP14809582.1A priority patent/EP3012622A4/en
Publication of WO2015008632A1 publication Critical patent/WO2015008632A1/ja

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
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/02—Constructional details
    • G01J5/08—Optical arrangements
    • G01J5/0806—Focusing or collimating elements, e.g. lenses or concave mirrors
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/02—Constructional details
    • G01J5/08—Optical arrangements
    • G01J5/0896—Optical arrangements using a light source, e.g. for illuminating a surface
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited

Definitions

  • the present invention relates to a thermal diffusivity measuring apparatus using a periodic heating radiation temperature measuring method.
  • Patent Document 1 describes an apparatus and a measurement method that use heating of a sample by laser irradiation using an optical fiber with a lens function and temperature measurement by an infrared light detector. In particular, in measurement in the thickness direction, It is described that it is heated by a circular laser beam having a diameter of about 5 mm.
  • Patent Document 2 discloses that an optical fiber with a lens function can be used to reduce the cost of the apparatus and can be applied to various types of samples. In particular, for samples having a large sample thickness and samples having a large thermal diffusivity. In order to fully respond, periodic heating temperature measurement can be provided with optical fiber position adjusting means for controlling the amount of heating per unit area by moving the optical fiber up and down from the sample and heating the surface of the sample. An apparatus for measuring thermophysical properties has been proposed.
  • Patent Document 1 proposes a laser diameter of 5 mm. This is because heat is transmitted two-dimensionally not only in the thickness direction but also in the in-plane direction when the sample is irradiated with a laser. Since the radiation thermometer will detect the two-dimensionally transmitted heat by reducing the irradiation area (laser diameter), the thermal diffusivity in the thickness direction cannot be measured accurately. It is for thinking.
  • Patent Document 1 if an attempt is made to obtain sufficient (detectable) heating performance with a laser diameter of 5 mm, the laser itself becomes expensive, and there is a problem that cannot be adopted due to cost problems.
  • the phase lag is measured while changing the frequency of the heating laser beam, and the thermal diffusivity is obtained from the phase lag.
  • the phase delay decreases as the frequency is increased, and theoretically decreases linearly with respect to the square root of the increasing frequency.
  • the linearly descending part can be used for calculating the thermal diffusivity as reliable data, the data before and after the inflection point are not adopted as unreliable data.
  • the conventional method has a problem that the inflection point is large depending on the laser diameter, and the data is not reliable.
  • the thermal diffusivity measurement device using the periodic heating radiation temperature measurement method consists of laser irradiation (output system) connected to a computer and temperature measurement / data processing system (input system) by a radiation thermometer.
  • a time lag time delay corresponding to the processing time of the network or computer occurs between the input system and the output system. For this reason, in a sample with a short thermal diffusion time, this time delay acts greatly, resulting in a measurement error. Removing this measurement error also becomes a problem.
  • the inventors of the present invention reviewed the lens diameter and laser diameter (irradiation area) of the radiation thermometer.
  • the laser irradiation area is increased to prevent the two-dimensionally transmitted heat from being transmitted to the radiation thermometer. Reduce the diameter).
  • the lens diameter is made as small as possible according to the laser irradiation area so that the temperature measurement range is adjusted to the irradiated portion, but conversely, the measurement noise is reduced by increasing it.
  • the inventors have developed a thermal diffusivity calculation method that takes into account the time delay by measuring how much time delay occurs in materials having different thermal diffusion times.
  • produces with the general apparatus structure of the thermal diffusivity measurement using a period heating radiation temperature measuring method can be reduced. Since noise can be reduced, it no longer depends on the measurement technique of the measurer. That is, in the past, efforts have been made to reduce noise by the technique of the measurer, such as measurement conditions and the number of measurements, and how the measurement noise is evaluated depends on the subjectivity of the measurer. However, according to the present invention, the measurement noise itself can be reduced. Therefore, measurement data can be obtained regardless of the measurement technique and subjectivity of the measurer.
  • FIG. 1 is a conceptual diagram of a thermal diffusion measuring apparatus according to the present invention.
  • FIG. 2 is a scatter diagram showing variations in thermal diffusivity when measured with a detection lens diameter of 50 mm and a laser diameter of 1000 ⁇ m.
  • FIG. 3 is a scatter diagram showing variations in thermal diffusivity when measured with a detection lens diameter of 50 mm and a laser diameter of 150 ⁇ m.
  • FIG. 4 is a table showing a guide for the ratio of the literature values of the thermal diffusivity when the detection lens diameter and the laser diameter are changed.
  • FIG. 5 is a table showing the phase delay when the laser diameter is changed.
  • FIG. 6 shows the literature value ratio between the thermal diffusivity and the thermal diffusivity for each type of metal. Time delay is not taken into consideration.
  • FIG. 2 is a scatter diagram showing variations in thermal diffusivity when measured with a detection lens diameter of 50 mm and a laser diameter of 1000 ⁇ m.
  • FIG. 3 is a scatter diagram showing variations in thermal diffusivity when measured
  • FIG. 7 shows the literature value ratio of thermal diffusivity to thermal diffusivity for each type of metal.
  • the time delay is 14 ⁇ s.
  • FIG. 8 shows an example of the thickness corresponding k calculation position storage means. The vertical axis represents the sample thickness, and the horizontal axis represents the square root of the applied laser frequency.
  • FIG. 9 is a block diagram of a thermal diffusivity measuring apparatus according to the present invention.
  • FIG. 1 is a conceptual diagram of a thermal diffusivity measuring apparatus according to the present invention.
  • the heating laser 2 modulates the light intensity at the frequency f by the modulation signal emitted from the function generator 1.
  • the heating laser 2 is composed of a semiconductor laser having a wavelength of 808 nm and an output of 5 W, for example.
  • the heating laser beam emitted from the heating laser 2 passes through the condenser lens 4 and is irradiated on the upper surface of the sample 5.
  • the thermal diffusivity measurement error in the thickness direction is reduced by adjusting the laser diameter through the condenser lens 4.
  • the infrared light condensing optical system 6 is installed on the side opposite to the heating laser beam irradiation side with the sample 5 interposed therebetween.
  • An infrared fiber 8 is connected to the radiation thermometer 9.
  • the infrared condensing optical system 6 is composed of two CaF2 plano-convex lenses combined with the convex sides back to back, condensing the infrared rays radiated from the sample surface, and applying infrared rays to the mid-infrared optical fiber. Have a role to input.
  • the radiation thermometer 9 detects the intensity of infrared light emitted from a region having a diameter of 250 ⁇ m on the lower surface of the sample 5.
  • the infrared light converging optical system 6 includes two lenses, such as CaF2, Si, Ge and ZnSe, which are finite correction systems and materials coated with gold, using a material having good transmittance for light having a wavelength of 1 ⁇ m to 10 ⁇ m.
  • a parabolic mirror or the like can be used.
  • the radiation thermometer a radiation thermometer using InSb or the like as a detection element can be used. According to the present invention, a sufficient effect can be obtained even with a single element.
  • the infrared light detected by the radiation thermometer 9 is converted into an electric signal, and the phase difference ⁇ of the temperature change that vibrates at the frequency f is measured by the lock-in amplifier 10.
  • the value of ⁇ is a phase difference with the function generator 1 branching a modulation signal for modulating the heating laser 2 and inputting it to the lock-in amplifier 10, and using the modulation signal as a reference.
  • the PC 11 records the phase difference ⁇ , for example, every f0. 5 (every 0.5 Hz).
  • time delay correction In the present invention, the time delay generated in the lock-in amplifier 10 is corrected. A correction method will be described. After measuring the phase difference with the lock-in amplifier, the PC 11 records the measured phase, the modulation frequency, and the phase difference between the two.
  • the PC 11 converts the phase difference into time. Further, the PC 11 has a setting unit for arbitrarily setting a time delay according to the characteristics of the lock-in amplifier. After correcting the time delay set by the setting unit, the phase difference is recalculated and stored as a new phase difference ⁇ .
  • the method of obtaining the slope k is, for example, as follows.
  • the PC 11 obtains the slope by the least square method by taking the frequency and phase difference at the 10 front and rear positions for each frequency with the square root of the frequency as the horizontal axis and the phase difference ⁇ as the vertical axis. That is, the slope for each frequency can be calculated.
  • This slope is compared with the slope corresponding to the before and after frequencies measured every f0.5. As a result of the comparison, the smallest part is adopted as the slope k where the difference falls within 5%.
  • the above calculation of k is performed because when the square root of the frequency is plotted on the horizontal axis and the phase difference ⁇ is plotted on the vertical axis, the linear region changes as the thickness of the sample changes. According to the experiment, when the sample is thinned, the range in which the linear region appears is changed to a higher frequency (for example, see the experimental result of 50 ⁇ m in FIG. 8). If the thickness is equal to or greater than a certain value, the difference in slope by the least square method is small regardless of the frequency (for example, see the experimental result of 1000 ⁇ m in FIG. 8).
  • the frequency and phase difference to be selected for calculating the slope k can be specified in consideration of the above situation. That is, the frequency to be selected by the thickness-corresponding k calculation position storage means as shown in FIG. 8 is specified in advance, and an inclination corresponding to a suitable frequency is adopted as k according to the thickness input by the thickness input means. I can do it. In addition, it is not limited to the Example of FIG. It is also possible to memorize whether the horizontal axis is a phase difference and whether it is appropriate for calculating k in the same manner.
  • the lens 4 When measuring the thermal diffusivity in the thickness direction, the lens 4 is used, and the surface of the sample 5 is heated by circular irradiation.
  • the infrared light collecting optical system 6 is arranged so as to measure the temperature of the center of the heated circle with the sample having a thickness d interposed therebetween.
  • the phase difference ⁇ taking into account the time delay is measured, and the PC 11 plots and displays the phase difference ⁇ with respect to f0.5.
  • the inclination k calculation means calculates the inclination k and calculates the thermal diffusivity.
  • the thermal diffusivity ⁇ in the thickness direction of the sample 1 is expressed by the following equation.
  • d is the thickness of the sample.
  • the apparatus according to the present invention can measure the thermal diffusivity in the thickness direction of the sample.
  • the horizontal axis represents the thermal diffusion time ( ⁇ seconds)
  • the vertical axis represents the ratio of the thermal diffusivity to the literature value.
  • the vertical axis indicates that there is no measurement error as it is closer to 1.
  • the detection system lens diameter was 50 mm and the laser diameter was 1000 ⁇ m
  • the result was as shown in FIG.
  • the average value of the absolute values of errors from the literature values was 20.8%.
  • the samples used in FIG. 2 are pure materials of various metals. Including Ag, Cu, Mo, Ta, Ti, SUS.
  • the sample had an arbitrary thickness of 50-1000 ⁇ m.
  • the vertical axis and the horizontal axis are the same as those in FIG. 2, but the measurement is performed with a detection system lens diameter of 50 mm and a laser diameter of 150 ⁇ m.
  • the average of the absolute values of errors from the literature values was 8%.
  • FIG. 4 shows the document value ratio error of the thermal diffusivity when the combination of the laser diameter and the detection lens diameter is changed.
  • the document value ratio error of the thermal diffusivity is an average value of absolute values of errors from the document value.
  • the result of the literature value ratio error of thermal diffusivity was within 15%.
  • FIG. 5 is a plot for calculating the slope k.
  • the slope is linear, but when the error is large, the slope is curvilinear, so surrounding data cannot be adopted. Therefore, data that is not curvilinear is desirable, but as shown in FIG. 5, an inflection point is generated at 1000 ⁇ m, a short straight line is generated before and after that, and it is not clear in which part the slope k should be measured. It cannot be adopted as data.
  • 150 ⁇ m and 300 ⁇ m have a linear distribution, and it is easy to calculate the slope k.
  • FIG. 6 as in FIGS. 2 and 3, the horizontal axis represents the thermal diffusion time ( ⁇ sec), and the vertical axis represents the ratio of the thermal diffusivity to the literature value.
  • the measurement sample and thickness are the same as in FIG. The vertical axis indicates that there is no measurement error as it is closer to 1.
  • FIG. 6 shows a case where the diameter of the detection system lens is 50 mm and the laser diameter is 150 ⁇ m while changing the metal to be measured. In FIG. 6, time delay is not considered. Therefore, it can be seen that the thinner the sample (the shorter the thermal diffusion time), the greater the effect of time delay.
  • FIG. 8 shows an example of the thickness corresponding k calculation position storage means shown in the flowchart of FIG.
  • the vertical axis represents the thickness of the sample
  • the horizontal axis represents the square root of the frequency of the applied laser.
  • the thickness-corresponding k calculation position storage means may take a phase difference on the horizontal axis. This is because it is possible to determine whether or not the phase difference is appropriate for calculating k in the same manner.
  • FIG. 9 is a flowchart of the calculation means.
  • the PC 11 records the irradiation frequency and the phase difference for each frequency while changing the frequency. After the recording is completed, the time delay correction considering the time delay parameter is performed.
  • the correction method is as described above.
  • the PC 11 calculates the inclination k based on the thickness input from the sample thickness input unit. That is, the phase difference corresponding to the input thickness is read from the frequency stored in the thickness correspondence storage means, and the slopes are calculated and calculated as k using the least square method at the 10 points before and after the thickness correspondence frequency.
  • the present invention can be used to measure thermophysical values of materials widely used in advanced industries.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

【課題】 周期加熱放射測温法熱物性測定装置において、厚み方向の測定誤差を解消する。 【解決手段】 試料にレーザービームを周波数fで周期的に照射するための加熱レーザービーム照射手段と、試料のある一点から放射される赤外光を集光するための赤外光集光手段とが、試料を挟みそれぞれ対向する位置に配置され、試料の周期的な温度変化から熱拡散率の測定を行う 周期加熱放射測温法熱物性測定装置において、 前記赤外光集光手段により集光された赤外光を放射温度計まで導く赤外用ファイバを備え、 前記放射温度計の温度変化の周期と前記加熱レーザービームの周期との位相差θを測定し、当該位相差θと前記周波数fから熱拡散率を演算する制御手段を備え、 前記加熱レーザービーム照射手段は、最大直径300μmの円形に相当する面積で照射する手段を備え、 前記赤外光集光手段は、最大50ｍｍの集光レンズを備えることに特徴を有する。

Description

熱拡散率測定装置
　本発明は、周期加熱放射測温法を用いた熱拡散率測定装置に関するものである。
従来、材料の熱拡散率を測定するために試料の一部に周期的に加熱を加え加熱した位置からある距離を離れた部分の温度測定を行い、その温度変化の周期より熱拡散率を得る方法および装置が知られている。
　特許文献１には、レンズ機能付き光ファイバを用いたレーザ照射による試料の加熱と赤外光検知器による温度測定を用いる装置および測定法が記載されており、特に、厚み方向の測定においては、直径5mm程度の円形のレーザービームによって加熱されることが記載されている。
　また、特許文献２には、レンズ機能付光ファイバを用いて装置の低価格を図り、各種の試料に適用することができ、特に試料厚さの大小の試料および熱拡散率大小の試料に対して充分に対応するために、前記光ファイバを試料から上下動させ、試料の面を加熱する単位面積当りの加熱量を制御するようにした光ファイバ位置調整手段を備えることが出来る周期加熱測温法熱物性測定装置が提案されている。
特開2011-185852号公報 特開2009-139163号公報
　一般的に、放射温度計を用いる厚み方向の熱拡散率測定装置においては、レーザー径が1000μm以上の照射面積で加熱されることが多い。例えば、特許文献１では、レーザー径5mmで提案されている。何故ならば、試料にレーザーを照射すれば、厚み方向だけでなく、面内方向も含めて二次元的に熱が伝わる。照射面積（レーザー径）を小さくすることで、放射温度計には二次元的に伝達された熱が検出されることになると考えるため、厚み方向の熱拡散率を正確に測定することが出来ないと考えるためである。
しかし、特許文献１のように、５mmのレーザ径で十分な（検出できる程度の）加熱性能を得ようとするとレーザー自体が高価になるため、コストの問題から採用できない課題がある。
別の課題もある。周期加熱放射測温法を用いた熱拡散率測定方法においては、加熱レーザー光の周波数を変化させながら、位相遅れを測定し、位相遅れから熱拡散率を求める方法をとる。一般的に、周波数を上げれば位相遅れは小さくなるし、理論的には、増加する周波数の平方根に対して、直線的に下がっていく。直線的に下がっている部分は信頼できるデータとして熱拡散率の算出に使用出来るものの、変曲点の前後のデータは信頼出来ないデータとして採用しない。従来の方法では、レーザー径によっては変曲点が大きく出てしまい、データが信頼出来ないという問題があった。
さらに別の課題もある。すなわち、周期加熱放射測温法を用いた熱拡散率測定装置では、コンピュータに接続されたレーザー照射（出力系）、及び放射温度計による測温・データ処理系（入力系）からなるが、この入力系と出力系にネットワークやコンピュータの処理時間分のタイムラグ（時間遅れ）が生じる。そのため、熱拡散時間の短い試料では、この時間遅れが大きく作用して、測定誤差となる。この測定誤差を除去することも課題となる。
上記課題に鑑みて、本発明の発明者らは、放射温度計のレンズ径、レーザー径（照射面積）を見直した。
すなわち、一般的な厚み方向の測定ではレーザー照射面積を大きくして放射温度計に対して、二次元的に伝達された熱が伝わらないようにするところ、逆にレーザー照射面積を小さくした（レーザー径を小さくする）。
また、一般的な厚み方向の測定では、できるだけ、レンズ径をレーザー照射面積に合わせて小さくして測温範囲を照射部分に合わせるが、逆に大きくして測定ノイズを低減させた。
さらに、発明者らは、熱拡散時間の異なる材料で時間遅れがどの程度生じるかを測定して、時間遅れを考慮した熱拡散率算出方法を開発した。
本発明によれば、周期加熱放射測温法を用いた熱拡散率測定の一般的な装置構成で発生するノイズを低減させることができる。ノイズを低減させることができるため、測定者の測定技術に依存しなくなった。すなわち、従来、測定条件や測定回数等、測定者の技術によってノイズを減らす努力が行われていたし、測定ノイズをどのように評価するかについて、測定者の主観に依存していた。しかし、本発明によれば、測定ノイズそれ自体を低減させることが出来る。そのため、測定者の測定技術や主観によらずに測定データを求めることが出来る。
図１は本発明の熱拡散測定装置の概念図である。 図２は、検出レンズ径が50mm、レーザー径1000μmで測定した場合の熱拡散率のばらつきを示す散布図である。 図３は、検出レンズ径が50mm、レーザー径150μmで測定した場合の熱拡散率のばらつきを示す散布図である。 図４は、検出レンズ径とレーザー径を変化させた場合に、熱拡散率の文献値比の目安を示す表である。 図５は、レーザー径を変化させた場合の位相遅れを示した表である。 図６は、金属の種類別に、熱拡散率と熱拡散率の文献値比を示したものである。なお、時間遅れは考慮していない。 図７は、金属の種類別に、熱拡散率と熱拡散率の文献値比を示したものである。時間遅れを１４μ秒考慮済みのものである。 図８は、厚み対応k算出位置記憶手段の例を示したものである。縦軸に試料の厚み、横軸には印加したレーザのー周波数の平方根をとっている。 図９は本発明に係る熱拡散率測定装置のブロック図である。
本発明の実施例について図面を参照して説明する。図1は本発明の熱拡散率測定装置の概念図である。加熱レーザー２は、ファンクションジェネレータ１から発した変調用信号により周波数fで光強度の変調が行われる。加熱レーザー２は、例えば波長808nm、出力5Wの半導体レーザーで構成される。加熱レーザー２より発した加熱レーザービームは、集光レンズ４を通り、試料５の上面に照射される。本発明では、集光レンズ４を通過してレーザー径を調整することによって、厚み方向の熱拡散率測定誤差を低減するものである。
赤外光集光光学系６は、試料５を挟んで加熱レーザービームの照射側とは反対側に設置される。赤外ファイバ８が放射温度計９に接続されている。例えば、赤外集光光学系６は２枚のCaF2からなる平凸レンズが凸側を背中合わせにして組み合わされており、試料面から放射された赤外線を集光するとともに中赤外線用光ファイバに赤外線を入力する役目を持っている。
放射温度計９は、試料５の下面の直径250μmの領域から発せられる赤外光の強度が検知される。ここで赤外光集光光学系６には、CaF2、Si、Ge、ZnSeなど波長1μmから10μmの光に良好な透過性を有する材料を用いた有限補正系のレンズや金コートされた2つの放物面ミラーなどを用いることができる。放射温度計は、InSbなどを検出素子とする放射温度計を用いることが出来る。なお、本発明によれば単素子のものでも十分に効果が得られる。
放射温度計９によって検知された赤外光は電気信号へ変換され、ロックインアンプ１０により周波数fで振動する温度変化の位相差θが測定される。ここで前記θの値は、ファンクションジェネレータ１が加熱レーザー2を変調するための変調信号を分岐してロックインアンプ１０に入力し、前記変調信号を基準とした位相差である。PC１１は、位相差θを、例えば、f0 . 5ごとに（0.5Hzごとに）記録する。
（時間遅れの補正）
本発明では、ロックインアンプ１０で生じる時間遅れを補正する。補正する方法を説明する。ロックインアンプで位相差を測定した後、ＰＣ１１は計測した位相、変調周波数、及び両者の位相差を記録する。
PC１１は位相差を時間に換算する。また、PC１１は、ロックインアンプの特性に応じて、時間遅れを任意で設定するための設定部を有する。前記設定部で設定された時間遅れを補正した上で、位相差を再計算し、新たな位相差θとして記憶する。
（厚みに応じた傾きkの算出）
傾きkの求め方は、例えば、次のようなものである。PC１１は、周波数の平方根を横軸、位相差θを縦軸として、周波数ごとに前後１０箇所の周波数及び位相差をとって、最小二乗法による傾きを求める。すなわち、周波数ごとの傾きが計算できる。この傾きを、f0 . 5ごとに測定した前後の周波数に対応した傾きと比較する。比較した結果、差異が５％以内に収まる部分で最小の部分を傾きkとして採用する。
上記のようなkの計算は、周波数の平方根を横軸、位相差θを縦軸としてプロットした場合、試料の厚みが変化すると直線領域が変化するために行われる。実験によれば、試料が薄くなると、直線領域がでる範囲が周波数が高いほうに変化する（例えば、図８の５０μmの実験結果を参照）。厚みが一定以上であれば、周波数に関わらず、前記最小二乗法による傾きに差異は小さい（例えば、図８の１０００μmの実験結果を参照）。
本発明による装置では、上記事情を考慮して、傾きkを演算するための選択すべき周波数や位相差を特定することもできる。すなわち、図８のような厚み対応k算出位置記憶手段によって選択すべき周波数を予め特定しておき、厚み入力手段によって入力された厚みに応じて、好適な周波数に対応した傾きをkとして採用することが出来る。なお、図８の実施例に限定されない。横軸を位相差として、同じようにkを算出するのに適切かどうかを記憶しておくことも出来る。
厚さ方向の熱拡散率を測定する場合レンズ４を使用し、試料５ の表面を円形状の照射により加熱を行う。赤外光集光光学系６は、厚さdの試料を挟んで 、加熱された円の中心の温度を測定するように配置される。その後、周波数fを変更するステップごとに、時間遅れを考慮した位相差θを測定し、PC１１は、位相差θをf0 . 5 に対してプロットして表示する。傾きk算出手段は、傾きkを算出し、熱拡散率を算出する。このとき前記試料1の厚さ方向の熱拡散率αは次の式で表される。ここで、dは試料の厚みである。
Figure JPOXMLDOC01-appb-M000001
　　上述の通り、本発明に係る装置では、試料の厚み方向の熱拡散率を測定することが出来る。
図２は、横軸に熱拡散時間（μ秒）、縦軸に熱拡散率の文献値との比をとったものである。縦軸は１に近いほど測定誤差がないことを示す。検出系レンズの直径が50mm、レーザー径1000μmで測定した場合、図２のようになった。誤差の水準を定量化するため、文献値からの誤差の絶対値の平均値をとったところ、２０．８％となった。
なお、図２で使用した試料は各種の金属の純物質である。Ag,Cu,Mo,Ta,Ti,SUSを含む。試料は50-1000μmの任意の厚みとした。
図３は、縦軸と横軸は図２と同じであるが、検出系レンズの直径が50mm、レーザー径150μmで測定したものである。文献値からの誤差の絶対値の平均値をとったところ、８％となった。
図４は、レーザー径と検出レンズ径の組み合わせを変化させた場合の、熱拡散率の文献値比誤差がどの程度かを示したものである。熱拡散率の文献値比誤差とは、文献値からの誤差の絶対値の平均値である。図４に示された通り、レーザー径300μmまでで、かつ検出系レンズの直径が50mmまでの場合、熱拡散率の文献値比誤差が１５％以内の結果を得られた。
図５は、前記傾きkを算出するためのプロット図である。誤差が少ない場合には、傾きは直線的となるが、誤差が多い場合には、傾きは曲線的となるため、周辺のデータは採用できない。そのため、曲線的でないデータが望ましいが、図５が示す通り、1000μmでは変曲点を生じてその前後で短い直線が生じており、傾きkをどの部分で計測すればよいか判然としないため、データとして採用できない。一方、150μm、300μmでは、直線的な分布をしており、傾きkを算出するのが容易である。
図６は、図２及び図３と同様横軸に熱拡散時間（μ秒）、縦軸に熱拡散率の文献値との比をとったものである。測定試料や厚みは図２と同じである。縦軸は１に近いほど測定誤差がないことを示す。測定する金属を変えながら、検出系レンズの直径が50mm、レーザー径150μmで測定した場合、図６のようになった。図６では、時間遅れを考慮していない。そのため、試料が薄いほど（熱拡散時間が短いほど）、時間遅れの影響が大きく出ていることが見て取れる。
図７は、縦軸と横軸は図６と同じであるが、時間遅れを14μ秒として設定して計算したものである。熱拡散時間による測定誤差が少なくなっていることが分かる。
図８は、図９のフローチャートに示す厚み対応k算出位置記憶手段の例を示したものである。縦軸に試料の厚み、横軸には印加したレーザーの周波数の平方根をとっている。なお、厚み対応k算出位置記憶手段は、横軸に位相差をとっても良い。位相差によっても、同じようにkを算出するのに適切か否かを判定できるからである。
図９は、演算手段のフローチャートである。PC１１は周波数を変更しながら、周波数ごとに照射周波数と位相差を記録していく。記録が終了してから、時間遅れパラメータを考慮した時間遅れ補正を行う。補正の方法は前述のとおりである。
また、PC１１は試料厚み入力部から入力された厚みに基づいて、傾きkを演算する。すなわち、厚み対応記憶手段に記憶された周波数から、入力された厚みに対応した位相差を読み出し、前記厚み対応周波数の前後１０点を最小二乗法を用いて傾きを求めてkとして演算する。
本発明は、先端産業で広く用いられている材料の熱物性値を計測するために利用可能である。
１　ファンクションジェネレータ
２　加熱レーザー
３　光ファイバ
４　集光レンズ
５　試料
６　赤外光集光光学系
７　集光レンズ
８　赤外ファイバ
９　放射温度計
１０　ロックインアンプ
１１　PC

Claims (3)

  1. 試料にレーザービームを周波数fで周期的に照射するための加熱レーザービーム照射手段と、試料のある一点から放射される赤外光を集光するための赤外光集光手段とが、試料を挟みそれぞれ対向する位置に配置され、試料の周期的な温度変化から熱拡散率の測定を行う 周期加熱放射測温法熱物性測定装置において、
    前記赤外光集光手段により集光された赤外光を放射温度計まで導く赤外用ファイバを備え、
    前記放射温度計の温度変化の周期と前記加熱レーザービームの周期との位相差θを測定し、当該位相差θと前記周波数fから熱拡散率を演算する制御手段を備え、
    前記加熱レーザービーム照射手段は、最大直径300μmの円形に相当する照射面積で照射する手段を備え、
    前記赤外光集光手段は、最大直径50ｍｍの集光レンズを備えることに特徴を有する熱拡散率測定装置。
  2. 時間遅れパラメータを入力する入力手段を備え、
    前記熱拡散率を演算する制御手段が前記位相差θを計算する際に、前記入力手段により入力された時間遅れパラメータを考慮して熱拡散率を演算することに特徴を有する請求項１に記載の熱拡散率測定装置。
  3. 試料厚みを入力する入力手段を備え、
    前記熱拡散率を演算する制御手段は、前記入力された厚みに基づいて熱拡散率を演算することに特徴を有する請求項１又は２に記載の熱拡散率測定装置。
PCT/JP2014/067805 2013-07-17 2014-07-03 熱拡散率測定装置 WO2015008632A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US14/409,078 US20160033431A1 (en) 2013-07-17 2014-07-03 Thermal diffusivity measuring device
EP14809582.1A EP3012622A4 (en) 2013-07-17 2014-07-03 Thermal diffusivity measurement device

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2013-148743 2013-07-17
JP2013148743A JP2015021791A (ja) 2013-07-17 2013-07-17 熱拡散率測定装置

Publications (1)

Publication Number Publication Date
WO2015008632A1 true WO2015008632A1 (ja) 2015-01-22

Family

ID=52346101

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/067805 WO2015008632A1 (ja) 2013-07-17 2014-07-03 熱拡散率測定装置

Country Status (4)

Country Link
US (1) US20160033431A1 (ja)
EP (1) EP3012622A4 (ja)
JP (1) JP2015021791A (ja)
WO (1) WO2015008632A1 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3372340A1 (de) 2017-03-10 2018-09-12 Scheppach Fabrikation von Holzbearbeitungsmaschinen GmbH Verfahren zur steuerung eines elektrowerkzeugs und zur durchführung des verfahrens eingerichtetes elektrowerkzeug
WO2019204140A1 (en) * 2018-04-17 2019-10-24 The Government Of The United States Of America, As Represented By The Secretary Of The Navy Rapid multiplexed infrared 3d nano-tomography
CN109557129B (zh) * 2018-10-29 2021-05-11 同济大学 一种薄膜热扩散系数的测量方法
WO2020162121A1 (ja) * 2019-02-06 2020-08-13 パナソニックＩｐマネジメント株式会社 厚み計測方法及び厚み計測装置、並びに欠陥検出方法及び欠陥検出装置
DE102021127596A1 (de) 2021-10-22 2023-04-27 Linseis Messgeräte Gesellschaft mit beschränkter Haftung Temperaturleitfähigkeitsmessgerät
WO2023180906A1 (en) * 2022-03-20 2023-09-28 Nova Ltd. Combined photothermal and ocd for semi-opaque structures

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03156351A (ja) * 1989-08-30 1991-07-04 Mitsui Toatsu Chem Inc 交流加熱による熱拡散率測定方法、熱伝導率測定方法および熱拡散率測定装置
JPH06221913A (ja) * 1992-10-06 1994-08-12 Hamamatsu Photonics Kk 散乱吸収体の光学情報計測装置及び方法
JP2005091119A (ja) * 2003-09-17 2005-04-07 Ricoh Co Ltd 微小領域熱物性解析方法及び装置
JP2009139163A (ja) 2007-12-05 2009-06-25 Beteru:Kk 周期加熱放射測温法熱物性測定装置およびそれを用いた熱物性測定方法
JP2010169453A (ja) * 2009-01-20 2010-08-05 Yamanashi Gijutsu Kobo:Kk 異物検査装置及び検査方法
JP2011185852A (ja) 2010-03-10 2011-09-22 National Institute Of Advanced Industrial Science & Technology 熱拡散率測定装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1422518A1 (en) * 2002-11-20 2004-05-26 Communaute Europeenne Method and system for measuring the thermal diffusivity

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03156351A (ja) * 1989-08-30 1991-07-04 Mitsui Toatsu Chem Inc 交流加熱による熱拡散率測定方法、熱伝導率測定方法および熱拡散率測定装置
JPH06221913A (ja) * 1992-10-06 1994-08-12 Hamamatsu Photonics Kk 散乱吸収体の光学情報計測装置及び方法
JP2005091119A (ja) * 2003-09-17 2005-04-07 Ricoh Co Ltd 微小領域熱物性解析方法及び装置
JP2009139163A (ja) 2007-12-05 2009-06-25 Beteru:Kk 周期加熱放射測温法熱物性測定装置およびそれを用いた熱物性測定方法
JP2010169453A (ja) * 2009-01-20 2010-08-05 Yamanashi Gijutsu Kobo:Kk 異物検査装置及び検査方法
JP2011185852A (ja) 2010-03-10 2011-09-22 National Institute Of Advanced Industrial Science & Technology 熱拡散率測定装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
See also references of EP3012622A4
SHUN NAKAYAMA ET AL.: "Spot Shuki Kanetsuho ni yoru Banjo Shiryo no Laser -kei Oyobi Kando Bunpu o Koryo shita Netsu Kakusanritsu Sokuteiho", THE 33RD JAPAN SYMPOSIUM ON THERMOPHYSICAL PROPERTIES KOEN RONBUNSHU, 3 October 2012 (2012-10-03), pages 164 - 166, XP008175675 *

Also Published As

Publication number Publication date
EP3012622A1 (en) 2016-04-27
US20160033431A1 (en) 2016-02-04
EP3012622A4 (en) 2017-03-22
JP2015021791A (ja) 2015-02-02

Similar Documents

Publication Publication Date Title
WO2015008632A1 (ja) 熱拡散率測定装置
JP5558629B2 (ja) レーザ加工装置
JP5414058B2 (ja) 熱拡散率測定装置
US5713665A (en) Method and apparatus for thermal diffusivity measurement
JP6382912B2 (ja) 試料を光熱分析するための方法及び装置
WO2000058700A1 (fr) Systeme de mesure de temperature
JP5288533B2 (ja) 周期加熱放射測温法熱物性測定装置およびそれを用いた熱物性測定方法
JP6724543B2 (ja) レーザ光吸収率測定方法、レーザ光吸収率測定装置及びレーザ加工方法
JP5662742B2 (ja) 粒径計測装置及び粒径計測方法
US20090274191A1 (en) Single beam optical apparatus and method
EP0429081A2 (en) Silicon wafer temperature measurement by optical transmission monitoring
KR101557878B1 (ko) 유동 유체의 스루풋을 결정하기 위한 방법 및 장치
CN109421275A (zh) 用于制造三维物体的设备
CN109416332B (zh) 热扩散率测定装置、热扩散率测定方法以及程序
Bouchut et al. Fused silica thermal conductivity dispersion at high temperature
JP5015581B2 (ja) 熱物性測定方法
JP4761547B2 (ja) レーザフラッシュ法を用いた熱定数の測定方法
JP7331732B2 (ja) アバランシェフォトダイオードの評価方法
JP2013092502A (ja) 温度測定装置及び放射率測定装置
RU2591227C2 (ru) Способ и система для бесконтактного определения радиуса пучка излучения
Tong et al. Direct and quantitative photothermal absorption spectroscopy of individual particulates
Pordzik et al. Frontal pyrometric snapshot measurements of the keyhole wall temperature in laser welding of pure aluminium
Schneeberg et al. Optical properties of a new IR converter for laser beam analysis
CN106909083B (zh) 基于拉曼散射的测量区域温度实时控制系统和方法
JP6421500B2 (ja) ガス分析装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 2014809582

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 14409078

Country of ref document: US

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14809582

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE