JP2019060740A - 温度測定方法 - Google Patents
温度測定方法 Download PDFInfo
- Publication number
- JP2019060740A JP2019060740A JP2017185945A JP2017185945A JP2019060740A JP 2019060740 A JP2019060740 A JP 2019060740A JP 2017185945 A JP2017185945 A JP 2017185945A JP 2017185945 A JP2017185945 A JP 2017185945A JP 2019060740 A JP2019060740 A JP 2019060740A
- Authority
- JP
- Japan
- Prior art keywords
- sound
- temperature
- probe
- ultrasonic wave
- ultrasonic waves
- 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
Landscapes
- Measuring Temperature Or Quantity Of Heat (AREA)
Abstract
【課題】超音波を用いて物体内の温度分布を測定する温度測定方法を提供する。【解決手段】温度測定方法は、物体10の一側面に超音波の発振用探触子2及び受信用探触子3を一方向に距離Dをあけて配置し、超音波を物体10の一側面から入射させるとともに他側面で反射した超音波を検出して、物体10中を斜めに伝播した超音波の音速を算出するステップと、物体10中の複数の斜角伝播路程ごとに超音波の音速を算出するステップと、複数の斜角伝播路程ごとの超音波の音速に基づいて、物体10の一方向に沿った断面における超音波の音速分布を算出するステップと、物体10の前記断面における超音波の音速分布、物体10の前記断面の所定位置における基準温度及び物体10の温度と超音波の音速との関係式に基づいて、物体10の前記断面における温度分布を算出するステップと、を有する。【選択図】図2
Description
本発明は、超音波を用いて物体内部の温度分布を測定する温度測定方法に関する。
近年、IoTの普及が拡大しており、生産現場において機器の診断や製品の不良を検知するためにビッグデータの活用が期待されている。なかでも、塑性加工の分野では、高強度材料を加工する際の素材温度の上昇や、熱間加工時の金型の温度管理など、加工中の温度計測を望むケースが多い。例えば、加工中の素材の発熱に伴う温度上昇により、成形品の寸法精度が悪くなる、潤滑性が変化して成形不良となる、などの事例が報告されている。また、近年需要の高いハイテンなどのホットプレス成形においても、成形後の製品強度に及ぼす成形温度の影響が著しいため、成形中の温度計測が大変重要となっている。
成形中の金型温度を測定する際、金型に孔をあけて熱電対を挿入して金型内部の温度を測定することが一般的に行われている。しかし、実際に温度を計測したい部分は加工部の直下であり、その位置までの孔を加工することは困難であるうえ、挿入孔をあけることによる強度低下も懸念される。さらに、当該測定方法では、多点測定を行う場合、多数の孔を金型にあける必要があるために、実用上ほぼ不可能である。
一方で、超音波と熱伝導方程式の逆問題解析とを用いて試料内部の温度分布を測定する方法が提案されている(例えば特許文献1,2を参照)。特許文献1,2に記載の温度測定方法は、超音波の音速は温度依存性があるため、音速を求めることで対象材料の温度を得る方法である。
しかしながら、特許文献1,2では、熱伝導方程式の逆問題解析により温度分布を算出しているため、境界条件がない(加熱源がどこか分からない)と温度分布を求めることができない。
本発明は、上記課題に着目してなされたものであり、超音波を用いて物体内部の温度分布を測定することができる温度測定方法を提供することを目的とする。
上記課題を解決する本発明の温度測定方法は、物体の一側面に超音波の発振用探触子及び受信用探触子を一方向に距離をあけて配置し、前記発振用探触子により超音波を前記物体の一側面から入射させるとともに前記物体の他側面で反射した超音波を前記受信用探触子により検出して、前記物体中を斜めに伝播した超音波の音速を算出するステップと、前記発振用探触子及び前記受信用探触子の少なくとも一方について前記物体の一側面における位置を変えて、前記物体中の複数の斜角伝播路程ごとに超音波の音速を算出するステップと、複数の斜角伝播路程ごとの超音波の音速に基づいて、前記物体の一方向に沿った断面における超音波の音速分布を算出するステップと、前記物体の前記断面における超音波の音速分布、前記物体の前記断面の所定位置における基準温度、及び、予め求められた前記物体の温度と超音波の音速との関係式に基づいて、前記物体の前記断面における温度分布を算出するステップと、を有する。
本発明の好ましい実施形態の温度測定方法においては、前記受信用探触子で検出された超音波の波形から、斜角伝播路程ごとに超音波の伝播時間を算出し、斜角伝播路程の路程長さと伝播時間とに基づいて、超音波の音速を算出する。
また、本発明の好ましい実施形態の温度測定方法においては、前記受信用探触子で検出された超音波の波形に対してヒルベルト変換を行う。
また、本発明の好ましい実施形態の温度測定方法においては、前記受信用探触子からの検出信号をプリアンプで増幅する。
また、本発明の好ましい実施形態の温度測定方法においては、複数の斜角伝播路程ごとに、超音波を前記物体の一側面から垂直に入射させた場合の前記物体中を伝播する超音波の基準音速に基づき校正係数を算出して、算出した超音波の音速に当該校正係数を乗ずる。
また、本発明の好ましい実施形態の温度測定方法においては、前記物体の温度と超音波の音速との関係式は、前記物体の温度ごとに超音波を前記物体の一側面から垂直に入射させて前記物体の他側面で反射させた場合の超音波の伝播時間と路程長さとにより音速を算出することで、導出される。なお、この実施態様の温度測定方法においては、前記路程長さを、前記物体の温度に応じて前記物体の熱膨張係数を用いて補正することがさらに好ましい。
本発明の温度測定方法によれば、超音波を用いて物体内部の温度分布を精度よく測定することができる。
以下、本発明の実施形態について添付図面を参照して説明する。本発明の温度測定方法は、超音波を用いて物体内部の温度分布を測定するものである。超音波の音速は、温度依存性を有し、温度が上がるほど音速が下がることが知られており、音速を求めることで測定対象の物体の温度を得ることができる。超音波による温度測定は、応答性がよく、また測定対象の物体に孔などをあける必要がないため、強度低下のデメリットがない。本発明は、超音波による温度測定に、傾斜型CTの断層画像を生成する方法(ラミノグラフィ)を応用することで、物体内部の温度測定を可能としたことに特徴を有する。
測定対象の物体は、特に限定されるものではないが、工学・工業の分野において物体内部の温度分布をモニタリングしたいとの要望が数多くあり、例えば、塑性加工・射出成形におけるワーク(素材)の温度や金型の温度、切削加工・溶接・熱処理などのワーク(素材)の温度などを好適に測定することができる。
図1に、本発明の温度測定方法を実施するための温度測定装置1の概略構成を示す。温度測定装置1は、超音波を発振する発振用探触子2と、超音波を検出する受信用探触子3と、発振用探触子2に電気信号(駆動信号)を送信するパルサ4と、受信用探触子3から電気信号(検出信号)を受信するレシーバ5と、受信用探触子3及びレシーバ5の間に配置されるプリアンプ6と、検出した超音波の波形を測定するオシロスコープ7と、コンピュータなどの制御手段8と、ディスプレイなどの表示手段9と、を備えている。
発振用探触子2は、パルサ4からの電気信号(駆動信号)に基づいて超音波を発振する。発振用探触子2は、物体10の一側面(本実施形態では上面)に配置される。発振用探触子2から発振される超音波は、所定範囲(例えば0.1MHz以上25MHz以下)の周波数を有し、図2に示すように、物体10の一側面から入射される。物体10に入射した超音波は、物体10内を伝播し、一側面とは反対側の他側面(本実施形態では底面)で反射する。
受信用探触子3は、物体10の一側面に発振用探触子2と距離Dをあけて配置され、物体10の一側面から入射した超音波のうち、斜め方向に入射して他側面で反射した超音波を検出する。受信用探触子3は、検出した超音波を電気信号(検出信号)に変換してレシーバ5に送信する。なお、本実施形態では、受信用探触子3及びレシーバ5の間に配置されたプリアンプ6により増幅された電気信号(検出信号)がレシーバ5に送られる。
発振用探触子2及び受信用探触子3は、図3に示すように、物体10の一側面において、図示しない移動手段により、一方向(x方向)に移動可能である。本実施形態では、発振用探触子2及び受信用探触子3は、探触子間の距離Dを一定に維持できる冶具11に取り付けられている。冶具11を一方向に移動することにより、探触子間の距離Dを変化させることなく発振用探触子2及び受信用探触子3を一方向に移動させることができる。なお、冶具11は、探触子間の距離Dを変更することができ、探触子間の距離Dを、例えば10mm、20mm、30mm、40mm、50mm、・・・と10mm間隔で変えて、発振用探触子2及び受信用探触子3を冶具11に取り付けることができる。
受信用探触子3で検出される超音波の入射角θは、探触子間の距離Dに応じて設定される。入射角θは、例えば0度よりも大きく60度以下とすることができ、より広い角度に設定することが好ましい。
パルサ4は、発振用探触子2により超音波を発振させる電気信号(駆動信号)を生成する機能を有する。レシーバ5は、受信用探触子3からの電気信号(検出信号)を増幅する機能と、電気信号(検出信号)をアナログ−デジタル変換する機能とを有する。
オシロスコープ7は、レシーバ5からの電気信号(検出信号)に基づいて、検出した超音波の波形を測定する。オシロスコープ7により得られた波形信号は、コンピュータ8に出力される。
コンピュータ8は、プロセッサ(CPU)より構成される処理部と、RAM、ROM及びハードディスク等により構成される記憶部と、を備える。ROMやハードディスクには、物体10内の温度分布を解析するための各種コンピュータプログラムや、予め求められた物体10の温度と超音波の音速との関係式(物体10の温度に対する超音波の音速データ)などが格納されている。処理部は、詳細は後述するが、オシロスコープ7からの波形信号を解析して、物体10中を伝播した超音波の音速を算出し、これに基づき物体10内の超音波の音速分布を算出する。そして、前記既知の超音波の音速データなどに基づいて、物体10内の温度分布を算出する。そして、物体10内の温度分布をディスプレイ9に表示する。
なお、温度測定装置1は、超音波の発振及び検出が可能な送受信用探触子(図示せず)をさらに備えていてもよい。この探触子は、物体10の一側面より垂直方向に超音波を入射させ、物体10の他側面で反射した超音波を受信することで、物体10を垂直方向に伝播する超音波の音速算出に用いられる。
次に、本発明の温度測定方法の一実施形態について説明する。本発明の温度測定方法は、物体10の一側面に超音波の発振用探触子2及び受信用探触子3を一方向に距離Dをあけて配置し、発振用探触子2により超音波を物体10の一側面から入射させるとともに物体10の他側面で反射した超音波を受信用探触子3により検出して、物体中10を斜めに伝播した超音波の音速を算出するステップS1と、発振用探触子2及び受信用探触子3の少なくとも一方について物体10の一側面における位置を変えて、物体10中の複数の伝播路程ごとに超音波の音速を算出するステップS2と、複数の伝播路程ごとの超音波の音速に基づいて、物体10の一方向に沿った断面における超音波の音速分布を算出するステップS3と、物体10の前記断面における超音波の音速分布、物体10の前記断面の所定位置における基準温度、及び、予め求められた物体10の温度と超音波の音速との関係式に基づいて、物体10の前記断面における温度分布を算出するステップS4と、を有する。
ステップS1では、図2に示すように、発振用探触子2により物体10の一側面から入射した超音波のうち、斜め方向に伝播して物体10の他側面で反射した超音波を受信用探触子3により検出する。これにより、図4(A)のようなパルスエコーを得る。図4(A)のパルスエコーにおいて、最初に現れる波形は発振された超音波を受信したことによるものであり、二つ目の波形が物体10の他側面で反射した超音波(底面エコー)を受信したことによるものである。そのため、物体中10を伝播した超音波の音速を算出する方法としては、1つ目の超音波発振時の波形と2つ目の超音波の底面エコーの波形との時間間隔を発振から受信までの伝播時間とし、物体10の厚みHから斜角伝播路程の路程長さ(≒2H/cosθ)を算出して、路程長さ÷伝播時間より音速を算出することができる。
ここで、1つ目の超音波の底面エコーを拡大した波形を図4(B)に示す。1つ目の発振による波形と2つ目の底面エコーによる波形との時間間隔(伝播時間)を求めるためには、基準となる計測点を決める必要がある。この計測点としては、(i)波形の立ち上がり、(ii)波形のピーク、(iii)波形の横軸との交点が考えられる。ただし、(i)立ち上がりでは、なだらかに振幅が上昇するため測定が難しく、(iii)横軸との交点では、振幅の異なる波が合成された場合に交点がずれるなどといった問題点から、基準となる計測点は(ii)ピークとすることが好ましい。
(ii)波形のピークを基準点として上記時間間隔(伝播時間)を求める際には、受信用探触子3で検出された超音波の波形に対してヒルベルト変換を行うことが好ましい。ヒルベルト変換は、例えば鋼球の打撃試験によるコンクリート内部の欠陥検出の自動化に用いられているものである。ヒルベルト変換を行うことで、波形のピークの位置を抽出し易くなり、上記時間間隔(伝播時間)の測定誤差を低減することができる。上記時間間隔(伝播時間)の測定誤差の低減により、物体10内を伝播する超音波の音速の測定誤差、さらには物体10内の温度分布の測定誤差を低減することができる。
以下の表1は、ヒルベルト変換の有無による音速の測定精度の違いを調べた結果である。室温を23℃の一定の状態で、長さ60mm、幅60mm、厚み30mmのSKD11製の試験片に対して、図5に示すように、試験片の上面より垂直方向に超音波を入射し、試験片の底面で反射した超音波を受信することで、垂直伝播路程での超音波の発振から受信までの時間間隔を算出し、この時間間隔と路程長さ(試験片の厚みの2倍)とから音速を測定した。音速の測定はそれぞれ6回行い、6回の測定の平均値と、6回の測定の標準偏差を求めた。ヒルベルト変換を行わない場合には、音速の平均値は6137.7m/s、標準偏差は2.1m/sであったのに対して、ヒルベルト変換を行うと、音速の平均値は6135.4m/s、標準偏差は0.8m/sとなり、ヒルベルト変換を行うことで、音速測定についての測定誤差を大幅に小さくできることができ、測定精度を向上できることが確認された。
なお、垂直伝播路程で超音波を検出すると、図6のようなパルスエコーを得る。図6のパルスエコーにおいて、最初に現れる波形は発振された超音波を受信したことによるものであり、二つ目の波形が試験片の底面で最初に反射した超音波(第1底面エコー)を受信したことによるものであり、三つ目以降の波形は、第1底面エコーの超音波が試験片の上面でさらに反射し、再び底面で反射した超音波を受信したことによるものである。一般に、探触子と試験片との間には油などの液体を塗布して、微小な空気層を排除して測定を行うが、この油膜厚さが測定における発振から受信までの時間に影響を及ぼす。そのため、垂直伝播路程の超音波の音速を算出する場合には、2つ目の第1底面エコーの波形と三つ目の第2底面エコーの波形との時間間隔Tを発振から受信までの伝播時間として、音速を算出することが好ましい。
ステップS1における波形のピークの位置を抽出し易くする波形処理としては、ヒルベルト変換に限定されるものではなく、例えばフーリエ変換などを行うこともできる。
また、超音波を斜め方向に発振・受信する場合には、図4(A)に示すように、底面エコーの波形の高さが低いために、波形のピークなどを検出して1つ目の発振時の波形との時間間隔を測定することが困難となるおそれがある。そのため、受信用探触子からの電気信号(検出信号)をプリアンプ6にて予め増幅してレシーバ5に入力することで、波形を増大させて、波形のピークなどを検出し易くすることが好ましい。
以下の表2は、プリアンプ6の有無による音速の測定精度の違いを調べた結果である。室温を23℃の一定の状態で、長さ60mm、幅60mm、厚み30mmのSKD11製の試験片に対して、図2に示すように、試験片の上面より斜め方向に超音波を入射し、試験片の底面で反射した超音波を受信することで、超音波の発振から受信までの時間間隔と、路程長さとから音速を測定した。音速の測定は、発振用探触子2及び受信用探触子3の探触子間の距離Dを10mm、20mm、30mm、40mm、50mmとした5つの斜角伝播路程について各2回ずつ行い、その平均値と標準偏差を求めた。プリアンプを用いない場合には、音速の平均値は6125.2m/s、標準偏差は3.2m/sであったのに対して、プリアンプを用いると、音速の平均値は6132.7m/s、標準偏差は1.3m/sとなり、プリアンプを用いることで、斜角伝播路程の音速測定について測定誤差を大幅に小さくできることができ、測定精度を大幅に向上できることが確認された。
さらに、本出願人の検討によると、超音波を斜め方向に発振・受信する場合に、物体10内部の温度が一定であれば、超音波の入射角θに関わらず音速は一定値になるはずであるが、超音波の入射角θによって音速が大きく異なることが分かった。これは、垂直方向に超音波を発射させる探触子を用いて、超音波を物体10に対して斜め方向に入射させると、探触子固有の測定誤差が、超音波の入射角θに応じて音速の測定結果に生じるためである。よって、超音波を斜め方向に発振・受信する場合には、この探触子固有の測定誤差を是正することが好ましい。具体的には、温度が一定の試験片に対して、探触子により垂直方向に超音波を入射して、垂直伝播路程での正確な音速を算出し、これを基準音速VSとする。そして、上記温度が一定の試験片に対して、探触子により斜め方向に超音波を入射して、超音波の入射角θごとに斜角伝播路程での音速VNをそれぞれ算出した後、上記基準音速VSを用いて、超音波の入射角θごとに音速VSの校正係数αN(=VS/VN)を求める。この校正係数αNを、実際の物体10において算出された超音波の音速Vに対して、入射角θに応じて乗じることで、超音波を斜め方向に発振・受信した場合でも、音速の測定精度を向上することができる。
次に、ステップS2では、発振用探触子2及び受信用探触子3の少なくとも一方について、物体10の一側面の一方向における位置を変えて、物体10中の複数の斜角伝播路程について超音波の音速を算出する。例えば、図7に示すように、発振用探触子2を固定し、受信用探触子3を物体10の一側面の一方向に所定のピッチで移動させて、探触子間の距離D(入射角θ)を変えながら超音波を物体10内に伝播させることで、複数の斜角伝播工程A1,A2,A3,A4,A5,・・・について超音波の音速を測定する。なお、受信用探触子3を固定し、発振用探触子2を物体10の一側面の一方向に所定のピッチで移動させてもよい。あるいは、図8に示すように、発振用探触子2及び受信用探触子3の探触子間の距離Dを固定しながら、発振用探触子2及び受信用探触子3を物体10の一側面の一方向に所定のピッチで移動させて、超音波を物体10内に伝播させることで、複数の斜角伝播工程A1,A2,A3,A4,A5,・・・について超音波の音速を測定する。
物体10内の超音波の音速分布をより正確に算出するためには、物体10内のより多くの斜角伝播工程について超音波の音速を測定することが好ましい。例えば、図7においては、受信用探触子3を物体10の一側面の一方向に移動させ終わった後、発振用探触子2を物体10の一側面の一方向にずらして固定し、再度、受信用探触子3を物体10の一側面の一方向に所定のピッチで移動させて、探触子間の距離D(入射角θ)を変えながら超音波を物体10内に伝播させることで、多数の斜角伝播工程について超音波の音速を測定する。また、図8においては、発振用探触子2及び受信用探触子3を物体10の一側面の一方向に移動させ終わった後、探触子間の距離Dを大きくし、再度、発振用探触子2及び受信用探触子3を物体10の一側面の一方向に所定のピッチで移動させて超音波を物体10内に伝播させることで、多数の斜角伝播工程について超音波の音速を測定する。
次に、ステップS3では、複数の斜角伝播路程ごとの超音波の音速に基づいて、物体10の一方向に沿った断面(測定対象断面)における超音波の音速分布を算出する。超音波の音速分布の算出には、傾斜型CTの断層画像を生成する方法(ラミノグラフィ)を応用する。ラミノグラフィに用いられているX線の吸収率、具体的にはX線の通過経路における吸収率の積算値を、超音波の斜角伝播路程における音速の平均値と読み替えることによって、物体10内部の測定対象断面の音速分布を再構成することができる。
音速分布の再構成については、ラミノグラフィの最も基本的な方法として、以下の方法を用いることができる。まず、図9に示すように、物体10内の音速分布を測定する測定対象断面を長さ方向(x方向)m個×厚み方向(y方向)n個の複数のセルに分割する。そして、斜角伝播路程ごとの音速に対して、斜角伝播路程ごとのセルを通過する面積按分値(重み)を乗じて各セルに加算する。具体的には、図9に示された斜角伝播路程A1の超音波の音速V1について説明すると、斜角伝播路程A1について、その全長にわたって、音速の適用範囲を物体10の長さ方向(x方向)にセル幅dだけとり、斜角伝播路程A1が通過するセルごとに、当該セルを含む適用範囲についてその中に含まれる全てのセルの面積比率を算出する。例えば、斜角伝播路程A1が通過するセルC1については、このセルC1を含む適用範囲は、セルC1の他に、隣接するセルC2を含み、この適用範囲におけるセルC1とセルC2との面積比率はu:(1−u)である。この場合には、セルC1の面積按分値(重み)がuとなり、セルC2の面積按分値(重み)が(1−u)となり、セルC1に音速としてV1×uが加算され、セルC2に音速としてV1×(1−u)が加算される。同様にして、斜角伝播路程A1が通過する全てのセルについて、音速V1に面積按分値(重み)を乗じて加算する。また、物体10内のその他の斜角伝播路程についても、音速に各セルに応じた面積按分値(重み)を乗じて各セルに加算する。
そして、全ての斜角伝播路程について各セルに音速を加算した後、各セルに積算された音速の積算値を面積按分値(重み)の積算値で除することで、各セルの音速を求める。具体的には、図9に示されたセルC2を例にして説明すると、セルC2には、セルC2を通過する斜角伝播路程A1に応じた音速V1×(1−u)が加算される。また、セルC2には、セルC2を通過する斜角伝播路程A2の音速V2に面積按分値(重み)vを乗じた音速が加算されるため、音速の積算値は、V1×(1−u)+V2×vとなる。一方で、セルC2の面積按分値(重み)の積算値は、(1−u)+vとなるため、セルC2の音速Vは、以下の式で表される。
V={V1×(1−u)+V2×v}/{(1−u)+v}
これにより、物体10内の測定対象断面の音速分布を得ることができる。なお、図9では、二つの斜角伝播路程A1,A2について説明しているが、これ以上の多数の斜角伝播路程を用いた場合も同様にして各セルの音速を算出できる。
なお、音速分布の算出方法は上記方法に限定されるものではなく、その他のラミノグラフィのアルゴリズム(例えばフィルタ補正逆投影法や逐次近似再構成法など)を応用することによっても、音速分布を再構成することができる。
最後に、ステップS4では、物体10内の測定対象断面の所定位置における基準温度を測定する。そして、この基準温度と、事前に取得した物体10の温度と超音波の音速との関係式(物体10の温度に対する超音波の音速データ)と、ステップS3で算出した物体10内の測定対象断面の音速分布とに基づいて、各セルの音速を温度に換算することで、物体10内の測定対象断面の温度分布を算出する。なお、基準温度とする物体10の位置は、物体10の内部であってもよいし、物体10の表面であってもよい。
物体10の温度と超音波の音速との関係式(物体10の温度に対する超音波の音速データ)については、種々の温度に保持した物体10について、図5と同様に、物体10の一側面(上面)より垂直方向に超音波を入射し、他側面(底面)で反射した超音波を受信することで、垂直伝播路程での超音波の発振から受信までの時間間隔(伝播時間)を算出する。そして、この時間間隔(伝播時間)と路程長さ(物体10の厚みHの2倍)とから音速を測定することで、物体10の各温度に応じた超音波の音速を導出することができる。
なお、物体10の温度と超音波の音速との関係式を導出する際には、物体10の熱膨張を考慮することが好ましい。つまりは、物体10の厚みHは必ずしも一定ではなく、物体10の熱膨張を考慮すると温度に応じて変わる。物体10の厚みHは、超音波の音速を測定する際の路程長さに大きな影響を及ぼす。そのため、例えば室温時の物体10の厚みと、物体10の熱膨張係数と、音速を測定する際の物体10の温度とから、音速を測定する際の物体10の温度に応じた物体10の厚みを算出し、この物体10の熱膨張を考慮した厚みを用いて物体10の各温度における超音波の音速を測定することで、より正確に物体10の温度と超音波の音速との関係式を導出することができる。
以上の通り、本発明の温度測定方法によれば、ラミノグラフィの原理を応用して物体10中を伝播する超音波の音速分布を得て、音速の温度依存性(超音波の物体10内の音速が温度に依存するデータ)を利用して物体10内部の温度分布を測定するから、熱電対などを測定対象の物体10に挿入することなく、また、従来技術のように、熱伝導方程式の逆問題解析により温度分布を算出していないため、境界条件がなくても、物体10内部の温度分布を精度よく測定することができる。
また、測定対象の物体10の熱膨張を考慮することで、温度測定の高精度化を図ることができる。
以上、本発明の実施形態について詳述したが、上述した実施形態は、全ての点で例示であって制限的なものではないため、本発明は上述した実施形態に限定されるものではない。本発明の技術的範囲は、特許請求の範囲によって画定され、また特許請求の範囲の記載と均等の意味及び範囲内での全ての変更を含むものであり、よって、本発明は、本発明の趣旨を逸脱しない限りにおいて種々の変更が可能である。
例えば、上述した実施形態では、物体10内の温度分布の測定の際に、物体10内を斜めに伝播させた斜角伝播路程の超音波の音速だけを計測して物体10内の音速分布を算出し、温度分布に変換している。しかしながら、斜角伝播路程の超音波の音速に加えて、物体10内を垂直方向に超音波を伝播させて垂直伝播路程の超音波の音速も計測することで、物体10内の音速分布を算出し、この垂直伝播路程も含めた音速分布を変換することで温度分布を算出するようにしてもよい。
物体10内の垂直伝播路程の超音波の音速も用いて音速分布を測定する場合、まず、ステップS1及びS2では、物体10の一側面において送受信用探触子を長さ方向に、斜角伝播路程の場合と同様に所定の間隔(例えば10mm間隔)で移動させながら複数箇所で超音波を垂直方向に発振・受信して、複数箇所の垂直伝播路程について音速を測定する。次に、ステップS3において超音波の音速分布を算出するが、例えば図9に示す方法では、垂直伝播路程ごとの音速についても、垂直伝播路程ごとにセルを通過する面積按分値(重み)を乗じて各セルに加算する。図10を用いて説明すると、垂直伝播路程A1〜A8については、その全長にわたり、垂直伝播路程A1〜A8が通過する各セルだけがその音速範囲に含まれる。よって、各セルの面積按分値(重み)は1となり、各セルには、各セルを通過する垂直伝播路程A1〜A8の音速V1〜V8がそのまま往復分、加算される。そして、全ての斜角伝播路程及び垂直伝播路程について各セルに音速を加算した後、各セルに積算された音速の積算値を面積按分値(重み)の積算値で除することで、各セルの音速を求める。最後にステップS4では、物体10内の測定対象断面の所定位置における基準温度と、事前に取得した物体10の温度と超音波の音速との関係式(物体10の温度に対する超音波の音速データ)と、ステップS3で算出した物体10内の測定対象断面の音速分布とに基づいて、各セルの音速を温度に換算することで、物体10内の測定対象断面の温度分布を算出する。
物体10内の温度分布の算出に、斜角伝播路程に加えて垂直伝播路程の超音波の音速を算出することで、より正確な物体10内の音速分布を算出することができるため、温度分布の測定もより高精度とすることができる。
また、上記実施形態では、発振用探触子2及び受信用探触子3を1つずつ用意し、物体10内の複数の斜角伝播路程について超音波の音速を算出する際に、物体10の一側面において、例えば受信用探触子3を発振用探触子2に対して所定の間隔で順次移動させることで、受信用探触子3の位置を変えている(図7)。しかしながら、これに限らず、物体10の一側面に発振用探触子2を1つ設置するとともに複数の受信用探触子3を間隔をあけて設置し、発振用探触子2から物体10内に入射された超音波を複数の受信用探触子3により同時に受信することで、受信用探触子3を移動させなくても、物体10内の複数の斜角伝播路程について超音波の音速を算出することも可能である。また、物体10内の複数の垂直伝播路程における超音波の音速を算出する場合についても、物体10の一側面において、送受信用探触子を所定の間隔で順次移動させなくても(図9)、複数の送受信用探触子を間隔をあけて設置することで、物体10内の複数の垂直伝播路程について超音波の音速を算出することも可能である。
以下に実施例を掲げて本発明をさらに詳しく説明するが、本発明は以下の実施例のみに限定されない。
超音波を用いて内部温度を測定する測定対象物体(試験片)としては、プレス金型で多用されるSKD11を用いた。SKD11製の試験片の寸法は、長さ80mm×幅100mm×厚み30mmであった。超音波の発振や受信には、図1に示す装置を用い、探触子にはオリンパス製の「V1091(5MHz,φ3mm)」を、パルサ・レシーバにはオリンパス製の「5072PR」を、オシロスコープには横河電機製の「DLM2022」を、プリアンプにはオリンパス製の「超音波プリアンプ(型番5676)」を、それぞれを用いた。
試験片には、図11に示すように、一方の端面の幅方向及び厚み方向の中央部にφ10のヒータを挿入して、試験片をヒータにより加熱するとともに、ヒータを挿入した端面と反対側において底面を、ペルチェ素子を用いた冷却板により冷却し、その他の底面を断熱材により断熱した。この状態で試験片を十分な時間保持することで、試験片内部に温度分布を生じさせた。
試験片の一方の側面にφ1のK型シース熱電対を8本挿入し、試験片内部の所定断面における温度分布を計測した。試験片の厚み方向中央部において長さ方向に並べられた6本の熱電対の間隔は10mmとした。図12は、熱電対による試験片内の温度測定結果である。図中、基準温度位置は、超音波による温度分布の算出時に用いた基準温度位置である。
超音波による温度測定位置は、ヒータ位置を対称面として、熱電対により温度計測を行う位置とは反対側にした。試験片内の超音波の音速測定は、発振用探触子及び受信用探触子の探触子間の距離をまず10mmに固定して、発振用探触子及び受信用探触子を試験片の上面の前記温度測定位置を長さ方向に移動させて、複数の斜角伝播路程について音速を測定した。次に、探触子間の距離を、それぞれ20mm、30mm、40mm、50mmに固定して、それぞれ複数の斜角伝播路程について音速を測定した。
試験片内の超音波の音速分布の算出については,ラミノグラフィの最も基本的な方法として、図9に示すように、試験片の音速分布を測定する断面を8×6のセル(各セルの長さd:10mm、厚みh:5mm)に分割する。予め求めた斜角伝播路程ごとのセルを通過する面積案分値(重み)を斜角伝播路程ごとの音速に乗じて各セルに加算し、全ての斜角伝播路程の加算後に重みの積算値で除することで、各セルの音速を求めた。そして、事前に取得した試験片の温度と超音波の音速との関係式から、各セルの音速を温度に換算することで、温度分布を求めた。その測定結果を図13に示す。なお、図13の測定結果は、試験片の熱膨張を考慮しなかった場合の結果である。
図13によれば、熱電対による測定温度と、超音波により測定した温度分布とが、略一致している。よって、試験片内の温度分布を精度よく測定可能であることが確認された。
また、図14は、試験片の熱膨張を考慮した試験片の温度と超音波の音速との関係式を用いて、超音波により試験片内の温度分布を測定した測定結果を示す。熱膨張を考慮した場合には、超音波により測定した温度分布が熱電対による測定温度とさらに良い一致を示しており、試験片左右の高温部及び低温部でも誤差が0.5℃程度である。よって、熱膨張を考慮することで、温度測定の高精度化を図れることが確認された。
2 発振用探触子
3 受信用探触子
6 プリアンプ
10 物体
3 受信用探触子
6 プリアンプ
10 物体
Claims (7)
- 物体の一側面に超音波の発振用探触子及び受信用探触子を一方向に距離をあけて配置し、前記発振用探触子により超音波を前記物体の一側面から入射させるとともに前記物体の他側面で反射した超音波を前記受信用探触子により検出して、前記物体中を斜めに伝播した超音波の音速を算出するステップと、
前記発振用探触子及び前記受信用探触子の少なくとも一方について前記物体の一側面における位置を変えて、前記物体中の複数の斜角伝播路程ごとに超音波の音速を算出するステップと、
複数の斜角伝播路程ごとの超音波の音速に基づいて、前記物体の一方向に沿った断面における超音波の音速分布を算出するステップと、
前記物体の前記断面における超音波の音速分布、前記物体の前記断面の所定位置における基準温度、及び、予め求められた前記物体の温度と超音波の音速との関係式に基づいて、前記物体の前記断面における温度分布を算出するステップと、を有する、温度測定方法。 - 前記受信用探触子で検出された超音波の波形から、斜角伝播路程ごとに超音波の伝播時間を算出し、斜角伝播路程の路程長さと伝播時間とに基づいて、超音波の音速を算出する、請求項1に記載の温度測定方法。
- 前記受信用探触子で検出された超音波の波形に対してヒルベルト変換を行う、請求項2に記載の温度測定方法。
- 前記受信用探触子からの検出信号をプリアンプで増幅する、請求項1〜3のいずれかに記載の温度測定方法。
- 複数の斜角伝播路程ごとに、超音波を前記物体の一側面から垂直に入射させた場合の前記物体中を伝播する超音波の基準音速に基づき校正係数を算出して、算出した超音波の音速に当該校正係数を乗ずる、請求項1〜4のいずれかに記載の温度測定方法。
- 前記物体の温度と超音波の音速との関係式は、前記物体の温度ごとに超音波を前記物体の一側面から垂直に入射させて前記物体の他側面で反射させた場合の超音波の伝播時間と路程長さとにより音速を算出することで、導出される、請求項1〜5のいずれかに記載の温度測定方法。
- 前記路程長さを、前記物体の温度に応じて前記物体の熱膨張係数を用いて補正する、請求項6に記載の温度測定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017185945A JP2019060740A (ja) | 2017-09-27 | 2017-09-27 | 温度測定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017185945A JP2019060740A (ja) | 2017-09-27 | 2017-09-27 | 温度測定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2019060740A true JP2019060740A (ja) | 2019-04-18 |
Family
ID=66176536
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017185945A Pending JP2019060740A (ja) | 2017-09-27 | 2017-09-27 | 温度測定方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2019060740A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020218006A1 (ja) * | 2019-04-24 | 2020-10-29 | 日本電信電話株式会社 | 内部温度測定装置および方法 |
CN113376255A (zh) * | 2021-05-18 | 2021-09-10 | 大庆石油管理局有限公司 | 一种钢管激光电弧复合焊接环焊缝全自动超声波检测方法 |
-
2017
- 2017-09-27 JP JP2017185945A patent/JP2019060740A/ja active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020218006A1 (ja) * | 2019-04-24 | 2020-10-29 | 日本電信電話株式会社 | 内部温度測定装置および方法 |
JP2020180830A (ja) * | 2019-04-24 | 2020-11-05 | 日本電信電話株式会社 | 内部温度測定装置および方法 |
JP7151609B2 (ja) | 2019-04-24 | 2022-10-12 | 日本電信電話株式会社 | 内部温度測定装置および方法 |
CN113376255A (zh) * | 2021-05-18 | 2021-09-10 | 大庆石油管理局有限公司 | 一种钢管激光电弧复合焊接环焊缝全自动超声波检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102297898B (zh) | 金属三阶弹性常数的激光超声测定方法 | |
Zhang et al. | Monte Carlo inversion of ultrasonic array data to map anisotropic weld properties | |
Takahashi et al. | Ultrasonic monitoring of internal temperature distribution in a heated material | |
Cegla et al. | High temperature ultrasonic crack monitoring using SH waves | |
Yamada et al. | Noncontact monitoring of surface temperature distribution by laser ultrasound scanning | |
Hoche et al. | Density, ultrasound velocity, acoustic impedance, reflection and absorption coefficient determination of liquids via multiple reflection method | |
JP2019060740A (ja) | 温度測定方法 | |
CN108918667A (zh) | 一种楔体缺陷检测方法 | |
CN107748205A (zh) | 一种随温度变化的弹性常数测量方法 | |
KR20080018596A (ko) | 초음파 종파와 횡파를 이용한 두께 측정 장치 및 방법 | |
Ihara et al. | New ultrasonic thermometry and its applications to temperature profiling of heated materials | |
BR112014014416B1 (pt) | método para determinar as propriedades elásticas de uma peça esférica | |
JP2010500591A (ja) | 溶接継手の微細構造の検査方法 | |
CN106546353B (zh) | 一种均匀材质构件内部温度场的超声测量方法 | |
US9222918B2 (en) | Sizing of a defect using phased array system | |
Takahashi et al. | Ultrasonic determination of temperature distribution in thick plates during single sided heating | |
JP5575157B2 (ja) | 超音波探傷装置、方法及びプログラム | |
JP3629256B2 (ja) | 超音波tofd法による探傷方法 | |
Barakat et al. | Evaluation of new designed reference blocks for calibration and NDT by optical and ultrasonic techniques | |
Ihara et al. | Non-invasive monitoring of temperature distribution inside materials with ultrasound inversion method | |
Ihara et al. | Ultrasonic thermometry for temperature profiling of heated materials | |
Shamu et al. | Characterization of acoustic beam propagation through high-grade stainless steel pipes for improved pulsed ultrasound velocimetry measurements in complex industrial fluids | |
Hoyle et al. | Ultrasonic algorithms for calculating probe separation distance, combined with full matrix capture with the total focusing method | |
JP2000221076A (ja) | 超音波音速測定方法 | |
Kazys et al. | Online profiling of nonplanar objects by high-resolution air-coupled ultrasonic distance measurements |