"測定対象物の熱歪を測定する歪測定装置であって、前記測定対象物の第1ステレオ画像を生成する第1ステレオカメラと、前記第1ステレオカメラと離間した位置に配設され、前記測定対象物の第2ステレオ画像を生成する第2ステレオカメラと、前記第1ステレオ画像及び前記第2ステレオ画像から前記測定対象物の三次元形状を求めて、前記実歪を求める実歪算出手段と、前記測定対象物の温度分布を検出するサーモビュアと、前記サーモビュアによって検出された温度分布から前記熱自由歪を求める熱自由歪算出手段と、前記実歪算出手段によって求められた前記実歪から、前記熱自由歪算出手段によって求められた前記熱自由歪を減じた差を、拘束歪として求める拘束歪算出手段と、を備えることを特徴とする歪測定装置。請求項1に記載の歪測定装置において、前記サーモビュアは、前記第1ステレオカメラ又は前記第2ステレオカメラと近接して配置されていることを特徴とする歪測定装置。請求項1又は請求項2に記載の歪測定装置において、前記測定対象物の線膨張係数と前記測定対象物の温度とを対応付けて記憶する膨張係数記憶手段を更に備え、前記熱自由歪算出手段は、前記サーモビュアによって検出された温度分布に含まれる各温度に対応する線膨張係数を、前記膨張係数記憶手段から読み出して、前記熱自由歪を求めることを特徴とする歪測定装置。請求項1から請求項3のいずれか1つに記載の歪測定装置において、前記サーモビュアから前記測定対象物に向かう視線方向と直交する面と前記測定対象物の表面とのなす角である傾斜角に対応付けて、前記サーモビュアによって検出された温度分布に含まれる各温度を補正する補正係数を記憶する補正係数記憶手段と、前記実歪算出手段によって求められた前記測定対象物の三次元形状から前記傾斜角を求め、求められた傾斜角に対応する補正係数を前記補正係数記憶手段から読み出して、読み出された補正係数を用いて、前記サーモビュアによって検出された温度分布を補正する温度補正手段と、を備え、前記熱自由歪算出手段は、前記温度補正手段によって補正された温度分布から前記熱自由歪を求めることを特徴とする歪測定装置。請求項1から請求項4のいずれか1つに記載の歪測定装置を用いて、前記測定対象物の温度と線膨張係数との関係を求める線膨張係数測定方法であって、前記測定対象物と同一の材料からなり、少なくとも2つの位置にマークが描かれたテストピースを、予め設定された温度に均一に加熱して、前記測定対象物と略同一の位置に配置する均一加熱工程と、前記テストピースに描かれた2つのマーク間の距離を前記第1ステレオカメラ及び前記第2ステレオカメラの少なくとも一方を用いて求める測定工程と、前記測定工程において求められた2つのマーク間の距離から線膨張係数を求める係数算出工程と、を繰り返し実行することによって、前記測定対象物の温度と線膨張係数との関係を求めることを特徴とする線膨張係数測定方法。請求項1から請求項4のいずれか1つに記載の歪測定装置を用いて、前記サーモビュアから前記測定対象物に向かう視線方向と直交する面と前記測定対象物の表面とのなす角である傾斜角と前記サーモビュアの補正係数との関係を求めるサーモビュアの補正係数測定方法であって、前記測定対象物と同一の材料からなり、円柱形のテストピースを、予め設定された温度に均一に加熱して、前記測定対象物と略同一の位置に配置する均一加熱工程と、前記サーモビュアによって前記テストピースの表面の温度分布を検出する検出工程と、前記検出工程において検出された前記テストピースの表面の温度分布から補正係数を求める係数算出工程と、を実行することによって、前記測定対象物の傾斜角と補正係数との関係を求めることを特徴とするサーモビュアの補正係数測定方法。しかしながら、上記特許文献1に記載の三次元形状測定装置では、測定対象物の歪みを測定することはできるものの、熱歪によって発生する応力を求めることはできない。すなわち、例え、エンジンの排気系部品が高温になったときであっても、エンジンの排気系部品が自由膨張する場合には、上記応力は作用せず、上記応力は、エンジンとの接続部等における外部からの拘束を受けること、及び、エンジンの排気系部品の温度分布等によって作用するものであるからである。すなわち、上記三次元形状測定装置では、測定対象物の寿命に影響する上記応力熱歪によって発生する応力に対応する歪この歪を、本願においては、拘束歪というを測定することはできない。本発明は、上記課題に鑑みてなされたものであって、拘束歪を容易に測定することの可能な歪測定装置、線膨張係数測定方法、及び、サーモビュアの補正係数測定方法を提供することを目的としている。上記課題を解決するために、本発明に係る歪測定装置は、以下のように構成されている。すなわち、本発明に係る歪測定装置は、測定対象物の熱歪を測定する歪測定装置であって、前記測定対象物の第1ステレオ画像を生成する第1ステレオカメラと、前記第1ステレオカメラと離間した位置に配設され、前記測定対象物の第2ステレオ画像を生成する第2ステレオカメラと、前記第1ステレオ画像及び前記第2ステレオ画像から前記測定対象物の三次元形状を求めて、前記実歪を求める実歪算出手段と、前記測定対象物の温度分布を検出するサーモビュアと、前記サーモビュアによって検出された温度分布から前記熱自由歪を求める熱自由歪算出手段と、前記実歪算出手段によって求められた前記実歪から、前記熱自由歪算出手段によって求められた前記熱自由歪を減じた差を、拘束歪として求める拘束歪算出手段と、を備えることを特徴としている。かかる構成を備える歪測定装置によれば、第1ステレオカメラによって前記測定対象物の第1ステレオ画像が生成される。また、前記第1ステレオカメラと離間した位置に配設された第2ステレオカメラによって、前記測定対象物の第2ステレオ画像が生成される。そして、前記第1ステレオ画像及び前記第2ステレオ画像から前記測定対象物の三次元形状が求められて、前記実歪が求められる。また、サーモビュアによって、前記測定対象物の温度分布が検出され、検出された温度分布から前記熱自由歪が求められる。更に、求められた前記実歪から、求められた前記熱自由歪を減じた差が、拘束歪として求められるため、拘束歪を容易に測定することができる。すなわち、本発明に係る歪測定装置は、測定対象物の実際の歪である実歪から、測定対象物の熱による自由膨張の歪である熱自由歪を減じた差を拘束歪として定義して、測定対象物の熱疲労に伴う寿命を評価する指標として、この拘束歪を測定するものである。そして、第1ステレオカメラ及び第2ステレオカメラによって、それぞれ、生成された前記第1ステレオ画像及び前記第2ステレオ画像から前記実歪が求められ、サーモビュアによって検出された前記測定対象物の温度分布から前記熱自由歪が求められるため、前記実歪から前記熱自由歪を減じた差である拘束歪を容易に測定することができるのである。また、本発明に係る歪測定装置は、前記サーモビュアが、前記第1ステレオカメラ又は前記第2ステレオカメラと近接して配置されていることが好ましい。かかる構成を備える歪測定装置によれば、前記サーモビュアが、前記第1ステレオカメラ又は前記第2ステレオカメラと近接して配置されているため、拘束歪を更に容易に測定することができる。すなわち、前記サーモビュアと近接配置されているカメラを基準として実歪を測定することによって、前記実歪を求める格子点と、前記熱自由歪を求める格子点との位置合わせが容易となるため、拘束歪を更に容易に測定することができるのである。また、本発明に係る歪測定装置は、前記測定対象物の線膨張係数と前記測定対象物の温度とを対応付けて記憶する膨張係数記憶手段を更に備え、前記熱自由歪算出手段が、前記サーモビュアによって検出された温度分布に含まれる各温度に対応する線膨張係数を、前記膨張係数記憶手段から読み出して、前記熱自由歪を求めることが好ましい。かかる構成を備える歪測定装置によれば、膨張係数記憶手段に、前記測定対象物の線膨張係数と前記測定対象物の温度とが対応付けて記憶されており、前記サーモビュアによって検出された温度分布に含まれる各温度に対応する線膨張係数図9参照が、前記膨張係数記憶手段から読み出されて、前記熱自由歪が求められるため、前記熱自由歪を容易に求めることができる。また、本発明に係る歪測定装置は、前記サーモビュアから前記測定対象物に向かう視線方向と直交する面と前記測定対象物の表面とのなす角である傾斜角に対応付けて、前記サーモビュアによって検出された温度分布に含まれる各温度を補正する補正係数を記憶する補正係数記憶手段と、前記実歪算出手段によって求められた前記測定対象物の三次元形状から前記傾斜角を求め、求められた傾斜角に対応する補正係数を前記補正係数記憶手段から読み出して、読み出された補正係数を用いて、前記サーモビュアによって検出された温度分布を補正する温度補正手段とを備え、前記熱自由歪算出手段が、前記温度補正手段によって補正された温度分布から前記熱自由歪を求めることが好ましい。かかる構成を備える歪測定装置によれば、前記サーモビュアから前記測定対象物に向かう視線方向と直交する面と前記測定対象物の表面とのなす角である傾斜角に対応付けて、前記サーモビュアによって検出された温度分布に含まれる各温度を補正する補正係数が、補正係数記憶手段に記憶されており、前記測定対象物の三次元形状から前記傾斜角が求められ、求められた傾斜角に対応する補正係数図12参照が前記補正係数記憶手段から読み出されて、読み出された補正係数を用いて、前記サーモビュアによって検出された温度分布が補正される。そして、補正された温度分布から前記熱自由歪が求められるため、前記熱自由歪を正確に且つ容易に求めることができる。また、本発明に係る線膨張係数測定方法は、以下のように構成されている。すなわち、本発明に係る線膨張係数測定方法は、上記いずれかの歪測定装置を用いて、前記測定対象物の温度と線膨張係数との関係を求める線膨張係数測定方法であって、前記測定対象物と同一の材料からなり、少なくとも2つの位置にマークが描かれたテストピースを、予め設定された温度に均一に加熱して、前記測定対象物と略同一の位置に配置する均一加熱工程と、前記テストピースに描かれた2つのマーク間の距離を前記第1ステレオカメラ及び前記第2ステレオカメラの少なくとも一方を用いて求める測定工程と、前記測定工程において求められた2つのマーク間の距離から線膨張係数を求める係数算出工程と、を繰り返し実行することによって、前記測定対象物の温度と線膨張係数との関係を求めることを特徴としている。かかる構成を備える線膨張係数測定方法によれば、前記測定対象物と同一の材料からなり、少なくとも2つの位置にマークが描かれたテストピースが、予め設定された温度に均一に加熱されて、前記測定対象物と略同一の位置に配置される。そして、前記テストピースに描かれた2つのマーク間の距離が求められ、求められた2つのマーク間の距離から線膨張係数が求められるため、上記いずれかの歪測定装置を用いて、前記測定対象物の温度と線膨張係数との関係を容易に求めることができる図7~図9参照。更に、本発明に係るサーモビュアの補正係数測定方法は、以下のように構成されている。すなわち、本発明に係るサーモビュアの補正係数測定方法は、上記いずれかの歪測定装置を用いて、前記サーモビュアから前記測定対象物に向かう視線方向と直交する面と前記測定対象物の表面とのなす角である傾斜角と前記サーモビュアの補正係数との関係を求めるサーモビュアの補正係数測定方法であって、前記測定対象物と同一の材料からなり、円柱形のテストピースを、予め設定された温度に均一に加熱して、前記測定対象物と略同一の位置に配置する均一加熱工程と、前記サーモビュアによって前記テストピースの表面の温度分布を検出する検出工程と、前記検出工程において検出された前記テストピースの表面の温度分布から補正係数を求める係数算出工程と、を実行することによって、前記測定対象物の傾斜角と補正係数との関係を求めることを特徴としている。かかる構成を備えるサーモビュアの補正係数測定方法によれば、前記測定対象物と同一の材料からなる円柱形のテストピースが、予め設定された温度に均一に加熱されて、前記測定対象物と略同一の位置に配置される。そして、前記サーモビュアによって前記テストピースの表面の温度分布が検出され、検出された前記テストピースの表面の温度分布から補正係数が求められるため、上記いずれかの歪測定装置を用いて、前記測定対象物の傾斜角と補正係数との関係を容易に求めることができる図10~図12参照。本発明に係る歪測定装置によれば、第1ステレオカメラによって前記測定対象物の第1ステレオ画像が生成される。また、前記第1ステレオカメラと離間した位置に配設された第2ステレオカメラによって、前記測定対象物の第2ステレオ画像が生成される。そして、前記第1ステレオ画像及び前記第2ステレオ画像から前記測定対象物の三次元形状が求められて、前記実歪が求められる。また、サーモビュアによって、前記測定対象物の温度分布が検出され、検出された温度分布から前記熱自由歪が求められる。更に、求められた前記実歪から、求められた前記熱自由歪を減じた差が、拘束歪として求められるため、拘束歪を容易に測定することができる。また、本発明に係る線膨張係数測定方法によれば、前記測定対象物と同一の材料からなり、少なくとも2つの位置にマークが描かれたテストピースが、予め設定された温度に均一に加熱して、前記測定対象物と略同一の位置に配置される。そして、前記テストピースに描かれた2つのマーク間の距離が求められ、求められた2つのマーク間の距離から線膨張係数が求められるため、上記歪測定装置を用いて、前記測定対象物の温度と線膨張係数との関係を容易に求めることができる。更に、本発明に係るサーモビュアの補正係数測定方法によれば、前記測定対象物と同一の材料からなる円柱形のテストピースが、予め設定された温度に均一に加熱されて、前記測定対象物と略同一の位置に配置される。そして、前記サーモビュアによって前記テストピースの表面の温度が測定され、測定された前記テストピースの表面の温度から補正係数が求められるため、上記歪測定装置を用いて、前記測定対象物の傾斜角と補正係数との関係を求容易にめることができる。以下、本発明の実施形態を、図面を参照して説明する。-歪測定装置のハード構成-まず、図1、図2を参照して、本発明に係る歪測定装置の構成について説明する。図1は、本発明に係る歪測定装置100の一例を示す構成図である。図2は、図1に示す歪測定装置100の平面図及び側面図である。図2aは、図1に示す歪測定装置100の平面図であり、図2bは、図1に示す歪測定装置100の側面図である。歪測定装置100は、測定対象物5の実際の歪である実歪εcから、測定対象物5の熱による自由膨張の歪である熱自由歪εtを減じた差を拘束歪εrとして測定する装置であって、図1に示すように、コンピュータ1、第1ステレオカメラ2、第2ステレオカメラ3、及び、サーモビュア4を備えている。また、第1ステレオカメラ2、第2ステレオカメラ3、及び、サーモビュア4の視界の略中心位置には、測定対象物5が配置されている。なお、ここでは、測定対象物5は、長方形状の板状体であって、その左右の両端は壁Wによって熱膨張が拘束されている。また、測定対象物5は、図略の熱源によって、その温度が第1温度T1例えば、20°Cから第2温度T2例えば、900°Cまで加熱される。本実施形態においては、便宜上、測定対象物5のz軸方向の変形は考えない測定対象としないものとする。測定対象物5の第1ステレオカメラ2第2ステレオカメラ3、及び、サーモビュア4側の正面51には、図1に示すように格子線が描かれている。図3は、図1に示す測定対象物5の正面51に描かれた格子の一例を示す図である。図3に示すように、測定対象物5の正面51に描かれた格子線は、x軸方向にN個例えば、50個、y軸方向にM個例えば、10個の領域に正面51を均等に分割している。ここでは、各格子をi,jで表す。添え字iは、x軸方向の格子の順番を表し、左端が1であり、右端がNである。同様に、添え字jは、y軸方向の格子の順番を表し、上端が1であり、下端がMである。また、以下の説明においては、各格子点i,jに対応する正面51の温度を、温度Ti,jと表記する。更に、各格子点i,jに対応する正面51の実歪εc、熱自由歪εt、及び、拘束歪εrを、それぞれ、実歪εci,j、熱自由歪εti,j、及び、拘束歪εri,jと表記する。ここで、添え字iは、1からNまでのいずれかの整数であり、添え字jは、1からMまでのいずれかの整数である。コンピュータ1は、いわゆるパーソナルコンピュータ等からなり、歪測定装置100全体の動作を制御するものであって、CPUCentralProcessingUnit、ROMReadOnlyMemory、RAMRandomAccessMemory等を備えている。ROMは、種々の制御プログラム等を記憶する。CPUは、ROMに記憶された種々の制御プログラムを読み出して各種処理を実行する。RAMは、CPUでの演算結果等を一時的に記憶するメモリである。また、コンピュータ1は、操作入力部、表示部、HDDHarddiskdrive等を備えている。操作入力部は、外部からの操作を受け付けるものであって、キーボード、マウス等からなる。表示部は、LCDLiquidCrystalDisplay等からなり、CPUによる演算結果等を外部から視認可能に表示するものである。HDDは、種々のデータを記憶するものである。第1ステレオカメラ2は、図2bに示すように、サーモビュア4の上に載置され、測定対象物5の第1ステレオ画像PS1を生成するカメラである。第1ステレオカメラ2によって生成された第1ステレオ画像PS1は、コンピュータ1図4に示す画像取得部11へ出力される。なお、第1ステレオ画像PS1は、後述する実歪算出部12によって実歪εci,jを求めるときの基準とされる画像である。第2ステレオカメラ3は、第1ステレオカメラ2と離間した位置に配設され、測定対象物5の第2ステレオ画像PS2を生成するカメラである。第2ステレオカメラ3によって生成された第2ステレオ画像PS2は、コンピュータ1図4に示す画像取得部11へ出力される。ここで、第2ステレオカメラ3は、図2aに示すように、測定対象物5の中心線CLについて第1ステレオカメラ2と線対称となる位置に配置されている。サーモビュア4は、第1ステレオカメラ2の下に配設され、測定対象物5の温度分布を検出する装置である。サーモビュア4によって検出された温度分布を示す情報は、コンピュータ1図4に示す温度取得部13へ出力される。上述のように、実歪εci,jを求めるときの基準とされる第1ステレオ画像PS1を生成する第1ステレオカメラ2が、サーモビュア4の上側に載置されているため、拘束歪εri,jを更に容易に測定することができる。すなわち、サーモビュア4と近接配置されている第1ステレオカメラ2を基準として実歪εci,jを測定することによって、実歪εci,jを求める格子点i,jと、熱自由歪εti,jを求める格子点i,jとの位置合わせが容易となるため、拘束歪εri,jを更に容易に測定することができるのである。本実施形態では、第1ステレオカメラ2が、サーモビュア4の上側に載置されている場合について説明するが、第1ステレオカメラ2又は第2ステレオカメラ3がサーモビュア4と近接して配置されている形態であればよい。例えば、第2ステレオカメラ3がサーモビュア4の真横に配置されている形態でもよい。この場合には、第2ステレオカメラ3によって生成される第2ステレオ画像PS2が、後述する実歪算出部12によって実歪εci,jを求めるときの基準となる画像として用いられる。-コンピュータの機能構成-次に、図4を参照してコンピュータ1の構成について説明する。図4は、図1にコンピュータの機能構成の一例を示す機能構成図である。コンピュータ1は、ROM等に記憶された制御プログラムを読み出して実行することによって、画像取得部11、実歪算出部12、温度取得部13、補正係数記憶部14、温度補正部15、膨張係数記憶部16、熱自由歪算出部17、及び、拘束歪算出部18等の機能部として機能する。ここで、実歪算出部12、補正係数記憶部14、温度補正部15、膨張係数記憶部16、熱自由歪算出部17、及び、拘束歪算出部18は、本発明に係る歪測定装置100の一部に相当する。画像取得部11は、第1ステレオカメラ2によって生成された第1ステレオ画像PS1、及び、第2ステレオカメラ3によって生成された第2ステレオ画像PS2を取得する機能部である。実歪算出部12は、画像取得部11によって取得された第1ステレオ画像PS1及び第2ステレオ画像PS2に基づいて、測定対象物5の正面51における実歪εci,jを求める機能部である。ここで、実歪算出部12は、特許請求の範囲に記載の実歪算出手段に相当する。具体的には、実歪算出部12は、画像取得部11によって取得された第1ステレオ画像PS1及び第2ステレオ画像PS2に基づいて、第1温度T1例えば、20°C及び第2温度T2例えば、900°Cの、それぞれの場合に、測定対象物5の正面51における各格子点i,jの中心位置のx-y座標x1ij、y1ij及びx2ij、y2ijを求める。ここで、添え字iは、1からNのいずれかの整数であり、添え字jは、1からMのいずれかの整数である。このように、2つのステレオ画像に基づいて各格子点i,jの中心位置のx-y座標を求める方法は、公知であるImageCorrelationforDeformationandShapeMeasurementsBasicConceptsTheoryandApplications、Chapter4:TwoDimensionalandThreeDimensionalComputerVision、P65P80、著者:SuttonMA、OrteuJ、2009年4月出版、出版社:Springer、ISBN:9780387787466、及び、特開2009-270915号公報、特開2001-241928号公報参照から、ここでは、その説明を省略する。次に、実歪算出部12は、第1温度T1及び第2温度T2の、それぞれの場合に、測定対象物5の正面51における各格子点i,jの中心位置の座標x1ij、y1ij及びx2ij、y2ijから、各格子点i,jにおける実歪εci,jを、幅方向x軸方向の実歪εcxi,j、及び、高さ方向y軸方向の実歪εcyi,jとして求める。幅方向の実歪εcxi,jは、例えば、次の1~4式によって求められる。εcxi,j=ΔLxi,j/L1xi,j1ΔLxi,j=L2xi,j-L1xi,j2L1xi,j=x1ij-x1i-1j2y1ij-y1i-1j2123L2xi,j=x2ij-x2i-1j2y2ij-y2i-1j2124また、高さ方向の実歪εcyi,jは、例えば、次の5~8式によって求められる。εcyi,j=ΔLyi,j/L1yi,j5ΔLyi,j=L2yi,j-L1yi,j6L1yi,j=x1ij-x1ij-12y1ij-y1ij-12127L2yi,j=x2ij-x2ij-12y2ij-y2ij-12128すなわち、実歪εci,jは、第1温度T1から第2温度T2に変化した場合の、各格子点i,jの中心位置から隣接する格子点の中心位置との距離の変化L1→L2に基づいて求められる。より具体的には、幅方向の実歪εcxi,jは、第1温度T1から第2温度T2に変化した場合の、各格子点i,jの中心位置から左側に隣接する格子点の中心位置までの距離の変化L1x→L2xに基づいて求められる。また、高さ方向の実歪εcyi,jは、第1温度T1から第2温度T2に変化した場合の、各格子点i,jの中心位置から上側に隣接する格子点の中心位置までの距離の変化L1y→L2yに基づいて求められる。なお、以下の説明においては、幅方向の実歪εcxi,j及び高さ方向の実歪εcyi,jを、便宜上、実歪εci,jと総称する。温度取得部13は、サーモビュア4によって検出された温度分布情報を取得する機能部である。また、温度取得部13は、取得した温度分布情報に基づいて、測定対象物5の正面51における各格子点i,jに対応する位置の温度Tijを求める。補正係数記憶部14は、サーモビュア4から測定対象物5に向かう視線方向と直交する面と測定対象物5の表面とのなす角である傾斜角θに対応付けて、サーモビュア4によって検出された温度分布に含まれる各温度を補正する補正係数βθ図9参照を記憶する機能部である。ここで、補正係数記憶部14は、特許請求の範囲に記載の補正係数記憶手段に相当する。具体的には、補正係数記憶部14は、ROM等にマップ又はLUTルックアップテーブルとして、傾斜角θに対応付けて補正係数βθを記憶する。また、本実施形態においては、補正係数βθは、図10~図12を用いて後述する本発明に係るサーモビュア4の補正係数測定方法によって求められ、補正係数記憶部14に記憶されているものとする。温度補正部15は、実歪算出部12によって求められた測定対象物5の三次元形状から傾斜角θを求め、求められた傾斜角θに対応する補正係数βθを補正係数記憶部14から読み出して、読み出された補正係数βθを用いて、サーモビュア4によって検出された温度分布を補正する機能部である。ここで、温度補正部15は、特許請求の範囲に記載の温度補正手段に相当する。具体的には、温度補正部15は、次の9式によって、サーモビュア4によって検出された各格子点i,jの温度TSi,jを補正して、補正後の温度TAi,jを求める。TAi,j=TSi,j/βi,j9ここで、βi,jは、各格子点i,jの傾斜角に対応する補正係数である。ただし、上述のように、本実施形態では、便宜上、測定対象物5のz軸方向の変形は考えない測定対象としないため、傾斜角θは0として補正が行われる。上述のように、サーモビュア4から測定対象物5に向かう視線方向と直交する面と測定対象物5の表面とのなす角である傾斜角θに対応付けて、サーモビュア4によって検出された温度分布Tに含まれる各温度TSi,jを補正する補正係数βθが、補正係数記憶部14に記憶されており、温度補正部15によって、測定対象物5の三次元形状から傾斜角θが求められ、求められた傾斜角θに対応する補正係数βθ図12参照が補正係数記憶部14から読み出されて、読み出された補正係数βθを用いて、サーモビュア4によって検出された温度TSi,jが補正される。そして、補正された温度TAi,jから熱自由歪εti,jが求められるため、熱自由歪εti,jを正確に且つ容易に求めることができる。膨張係数記憶部16は、測定対象物5の線膨張係数αと測定対象物5の温度Tとを対応付けて記憶する機能部である。ここで、膨張係数記憶部16は、特許請求の範囲に記載の膨張係数記憶手段に相当する。具体的には、膨張係数記憶部16は、ROM等にマップ又はLUTルックアップテーブルとして、温度Tに対応付けて線膨張係数αTを記憶する。また、本実施形態においては、線膨張係数αTは、図7~図9を用いて後述する本発明に係る線膨張係数測定方法によって求められ、膨張係数記憶部16に記憶されているものとする。熱自由歪算出部17は、サーモビュア4によって検出された測定対象物5の温度分布から熱自由歪εtを求める機能部である。熱自由歪算出部17は、特許請求の範囲に記載の熱自由歪算出手段に相当する。具体的には、熱自由歪算出部17は、サーモビュア4によって検出された温度分布に含まれる各格子点i,jの各温度Ti,jに対応する線膨張係数αTi,jを、膨張係数記憶部16から読み出して、次の10、11式によって熱自由歪εti,jを、それぞれ、幅方向の熱自由歪εtxi,j、及び、高さ方向の熱自由歪εtyi,jとして求める。εtxi,j=αTi,j-α0×L1xi,j10εtyi,j=αTi,j-α0×L1yi,j11ここで、α0は、測定対象物5の初期温度である第1温度T1における線膨張係数αの値である。また、10式のL1xi,j、及び、11式のL1yi,jは、それぞれ、上記3式、及び、上記7式で与えられる第1温度T1における測定対象物5の幅方向の格子間距離及び高さ方向の格子間距離である。なお、以下の説明においては、幅方向の熱自由歪εtxi,j及び高さ方向のεtyi,jを、便宜上、熱自由歪εti,jと総称する。上述のように、膨張係数記憶部16に、測定対象物5の線膨張係数αTと測定対象物5の温度Tとが対応付けて記憶されており、サーモビュア4によって検出された温度分布に含まれる各温度Ti,jに対応する線膨張係数αTi,j図9参照が、膨張係数記憶部16から読み出されて、熱自由歪εti,jが求められるため、熱自由歪εti,jを容易に求めることができる。拘束歪算出部18は、実歪算出部12によって求められた実歪εci,jから、熱自由歪算出部17によって求められた熱自由歪εti,jを減じた差を、拘束歪εri,jとして求める機能部である。ここで、拘束歪算出部18は、特許請求の範囲に記載の拘束歪算出手段に相当する。すなわち、拘束歪算出部18は、次の12式によって、各格子点i,jの拘束歪εri,jを求めるεri,j=εci,j-εti,j12ここで、図5を参照して、拘束歪εrについて説明する。図5は、本発明に係る拘束歪の考え方を示す概念図である。円板状の測定対象物5Aは、左右両端が壁WAに移動を拘束された状態で、第1温度T1例えば、20°Cから第2温度T2例えば、900°Cまで加熱されて熱膨張する。この場合に、温度上昇に伴って、測定対象物5Aは熱膨張するが、左右両端が壁WAに移動を拘束されているため、右側の図で示すように、第2温度T2まで加熱された状態では、縦方向に長い長円又は、楕円状の形状となる。本願発明に係る歪測定装置100は、図5の左端の図である初期状態と、図5の右側の図である加熱後の状態との間に、図5の中央の図で示す、測定対象物5Aが第2温度T2で自由膨張した状態を想定している。そして、拘束歪εrを、右端の図で示す状態における歪みである実歪εcと、中央の図で示す状態である熱自由歪εtとの差と定義して、この拘束歪εrを、測定対象物5Aに作用する応力に関連する情報として使用するものである。換言すれば、本願発明は、測定対象物5Aの加熱、冷却の繰り返しによる測定対象物5Aの寿命に影響する上記応力を、拘束歪εrを用いて推定することをその目的としている。-歪測定装置コンピュータ1の動作-次に、図6を参照して、本発明に係る歪測定装置100主に、コンピュータ1の動作を説明する。図6は、図4に示す歪測定装置100の動作の一例を示すフローチャートである。まず、画像取得部11によって、第1ステレオカメラ2及び第2ステレオカメラ3から、それぞれ、第1ステレオ画像PS1、及び、第2ステレオ画像PS2が取得されるステップS101。そして、実歪算出部12によって、ステップS101で取得された第1ステレオ画像PS1、及び、第2ステレオ画像PS2に基づいて、測定対象物5の正面51における実歪εci,jが求められるステップS103。次いで、温度取得部13によって、サーモビュア4で検出された各格子点i,jに対応する位置の温度TSijが取得されるステップS105。次に、温度補正部15によって、各格子点i,jの温度TSi,jが補正されて、補正後の温度TAi,jが求められるステップS107。そして、熱自由歪算出部17によって、ステップS107で補正された温度TAi,jに基づいて、熱自由歪εti,jが求められるステップS109。次に、拘束歪算出部18によって、ステップS103で求められた実歪εci,jからステップS109で求められた熱自由歪εti,jを減じた差として拘束歪εri,jが求められてステップS111、処理が終了される。上述のように、第1ステレオカメラ2によって測定対象物5の第1ステレオ画像PS1が生成される。また、第1ステレオカメラ2と離間した位置に配設された第2ステレオカメラ3によって、測定対象物5の第2ステレオ画像PS2が生成される。そして、第1ステレオ画像PS1及び第2ステレオ画像PS2から測定対象物5の三次元形状が求められて、実歪εci,jが求められる。また、サーモビュア4によって、測定対象物5の温度分布Tijが検出され、検出された温度分布Tijから熱自由歪εti,jが求められる。更に、求められた実歪εci,jから、求められた熱自由歪εti,jを減じた差が、拘束歪εri,jとして求められるため、拘束歪εri,jを容易に測定することができる。すなわち、本発明に係る歪測定装置100は、測定対象物5の実際の歪である実歪εci,jから、測定対象物5の熱による自由膨張の歪である熱自由歪εti,jを減じた差を拘束歪εri,jとして定義して、測定対象物5の熱疲労に伴う寿命を評価する指標として、この拘束歪εri,jを測定するものである。そして、第1ステレオカメラ2及び第1ステレオカメラ3によって、それぞれ、生成された第1ステレオ画像PS1及び第2ステレオ画像PS2から実歪εci,jが求められ、サーモビュア4によって検出された測定対象物5の温度分布Tijから熱自由歪εti,jが求められるため、実歪εci,jから熱自由歪εti,jを減じた差である拘束歪εri,jを容易に測定することができるのである。-線膨張係数測定方法-次に、図7~図9を参照して、本発明に係る線膨張係数測定方法について説明する。図7は、本発明に係る線膨張係数測定方法の一例を示す構成図である。図7に示すように、本発明に係る線膨張係数測定方法は、上述の歪測定装置100を用いて温度Tと線膨張係数αTとの関係を求めるものである。図1に示す測定対象物5に換えて、図7に示すように、テストピース6を、図1に示す測定対象物5と略同一の位置に配置される。テストピース6は、測定対象物5と同一の材料からなり、その正面61に、4つのマークM1~M4が描かれている。マークM1~M4は、正方形の4つの各頂点に対応する位置に描かれており、初期温度T0例えば、20°Cにおいて、上記正方形の一辺の長さはL0である。そして、テストピース6は、予め設定された温度TP例えば、第1温度T1=20°Cから第2温度T2=900°Cまでの間の10°C間隔の温度に均一に加熱されて、図1に示す測定対象物5と略同一の位置に配置される。この工程が、特許請求の範囲に記載の均一加熱工程に相当する。次に、第1ステレオカメラ2によって、テストピース6の正面61に描かれた2つのマーク間ここでは、マークM3とマークM4との間の距離LPが求められる。この工程が、特許請求の範囲に記載の測定工程に相当する。なお、距離LPは、温度TPが設定される度に例えば、第1温度T1=20°Cから第2温度T2=900°Cまでの間の10°C間隔の温度に設定される度に、求められる。そして、求められた2つのマーク間の距離LPから、次の13式によって線膨張係数αが求められる。この工程が、特許請求の範囲に記載の係数算出工程に相当する。α=LP-L0/L013次に、図8を参照して本発明に係る線膨張係数測定方法の手順を説明する。図8は、本発明に係る線膨張係数測定方法の処理手順の一例を示すフローチャートである。まず、テストピース6を加熱する温度TPが設定されるステップS201。そして、テストピース6がステップS201で設定された温度TPに均一に加熱されるステップS203。次いで、ステップS203で温度TPに均一に加熱されたテストピース6が図1に示す測定対象物5と略同一の位置に配置されるステップS205。ここで、ステップS203及びステップS205が、特許請求の範囲に記載の均一加熱工程に相当する。次に、テストピース6の2つのマーク間ここでは、マークM3とマークM4との間の距離LPが求められるステップS207。このステップS207が、特許請求の範囲に記載の測定工程に相当する。そして、上記13式によって、線膨張係数αが求められるステップS209。このステップS209が、特許請求の範囲に記載の係数算出工程に相当する。次に、測定を終了するか否かの判定が行われるステップS211。この判定は、例えば、20°Cから900°Cまでの間の10°C間隔の温度について、全て、線膨張係数αの測定が完了した場合に、終了すると判定される。ステップS211でYESの場合には、処理が終了される。ステップS211でNOの場合には、処理がステップS201に戻され、ステップS201以降の処理が繰り返し実行される。図9は、線膨張係数αの測定結果の一例を示すグラフG1である。横軸は、テストピース6又は、測定対象物5の温度TPであって、縦軸は、線膨張係数αである。このようにして、測定対象物5と同一の材料からなり、少なくとも2つの位置にマークが描かれたテストピース6が、予め設定された温度TPに均一に加熱されて、測定対象物3と略同一の位置に配置される。そして、テストピース6の2つのマークM3、M4間の距離LPが求められ、求められた2つのマークM3、M4間の距離LPから線膨張係数αが求められるため、上記の歪測定装置100を用いて、容易に測定対象物5の温度TPと線膨張係数αとの関係を求めることができる。本実施形態では、第1ステレオカメラ2を用いて距離LPが求められる場合について説明するが、第1ステレオカメラ2及び第2ステレオカメラ3の少なくとも一方を用いて距離LPを求める形態であればよい。-サーモビュアの補正係数測定方法-次に、図10~図12を参照して、本発明に係るサーモビュア4の補正係数βθの測定方法について説明する。図10は、本発明に係るサーモビュア4の補正係数βθの測定方法の一例を示す構成図である。図10に示すように、本発明に係るサーモビュア4の補正係数βθの測定方法は、上述の歪測定装置100を用いて、サーモビュア4から測定対象物5に向かう視線方向と直交する面と測定対象物5の表面とのなす角である傾斜角θと補正係数βθとの関係を求めるものである。図1に示す測定対象物5に換えて、図10に示すように、テストピース7を、図1に示す測定対象物5と略同一の位置に配置される。テストピース7は、測定対象物5と同一の材料からなり、円柱形に形成されている。そして、テストピース7は、予め設定された温度TQ例えば、600°Cに均一に加熱されて、図1に示す測定対象物5と略同一の位置に配置される。この工程が、特許請求の範囲に記載の均一加熱工程に相当する。次に、サーモビュア4によって、テストピース6の表面の温度TSθの分布が検出される。この工程が、特許請求の範囲に記載の検出工程に相当する。ここで、傾斜角θは、サーモビュア4から測定対象物5に向かう視線方向と直交する面とテストピース6の表面とのなす角である。そして、検出された温度TSθの分布から、次の14式によって補正係数βθが求められる。この工程が、特許請求の範囲に記載の係数算出工程に相当する。βθ=TSθ/TQ14次に、図11を参照して本発明に係るサーモビュア4の補正係数βθの測定方法の手順を説明する。図11は、本発明に係るサーモビュア4の補正係数βθの測定方法における処理手順の一例を示すフローチャートである。まず、テストピース7が温度TQに均一に加熱されるステップS301。次いで、ステップS301で温度TQに均一に加熱されたテストピース7が図1に示す測定対象物5と略同一の位置に配置されるステップS303。ここで、ステップS301及びステップS303が、特許請求の範囲に記載の均一加熱工程に相当する。次に、サーモビュア4によって、テストピース7の表面における温度TSθの分布が検出されるステップS305。このステップS305が、特許請求の範囲に記載の検出工程に相当する。そして、上記14式によって、サーモビュア4の補正係数βθが求められてステップS307、処理が終了される。このステップS307が、特許請求の範囲に記載の係数算出工程に相当する。図12は、サーモビュア4の補正係数βθの測定結果の一例を示すグラフである。図12aは、傾斜角θと検出温度TSθとの関係を示すグラフG2であり、図12bは、傾斜角θと補正係数βθとの関係を示すグラフG3である。図12aのグラフG2に示すように、傾斜角θの絶対値が大きい程、検出温度TSθは低下する。また、図12bのグラフG3に示すように、傾斜角θの絶対値が大きい程、補正係数βθは小さくなる。このようにして、測定対象物5と同一の材料からなる円柱形のテストピース7が、予め設定された温度TQに均一に加熱されて、測定対象物5と略同一の位置に配置される。そして、サーモビュア4によってテストピース7の表面における温度TSθの分布が検出され、検出されたテストピース7の表面における温度TSθの分布から補正係数βθが求められるため、歪測定装置100を用いて、容易に測定対象物5の傾斜角θと補正係数βθとの関係を求めることができる。本実施形態では、テストピース7が円柱形に形成されている場合について説明するが、テストピース7がその他の形状に形成されている形態でもよい。例えば、テストピース7が、平板状に形成されており、y軸方向図10の紙面と直交する方向と平行なテストピース7の中心軸を中心として、テストピース7を回転させて、傾斜角θを変化させながら、傾斜角θ毎のテストピース7の表面における温度TSθを求める形態でもよい。この場合には、傾斜角θと補正係数βθとの関係を更に正確に求めることができる。-他の実施形態-本実施形態では、歪測定装置100がコンピュータ1において実歪算出部12、補正係数記憶部14、温度補正部15、膨張係数記憶部16、熱自由歪算出部17、及び、拘束歪算出部18の機能部として構成されている場合について説明したが、実歪算出部12、補正係数記憶部14、温度補正部15、膨張係数記憶部16、熱自由歪算出部17、及び、拘束歪算出部18の機能部のうち、少なくとも1つの機能部が、電子回路等のハードウェアで構成されている形態でもよい。"