JP5033739B2 - 応力測定方法およびその装置 - Google Patents

Description

本発明は、応力測定方法およびその装置に係り、特に、加工工程における素材および製品（例えば半導体単結晶など）の内部応力を、又発電所，高圧プラント，航空機，車両や各種装置或いはその他の建造物等に組み込んだ部材内の未知の内部応力を、定量的に且つ迅速に測定することを可能とした応力測定方法およびその装置に関する。

金属部材を骨格として使用する各種車両や各種建造物については、一般にその耐用年数は、ほぼ半永久的となっている。これに対し、これら各建造物や車両を構成する各種部材については、部分的な損傷が進行することもあって、多くの場合、予め一定年数の経過とともに補修若しくは部品交換の必要性を定めている。

ここで、各種車両にあっては、その重量を支える車軸をはじめとして、連結用のボルトや動力伝達シャフト等には、静的な荷重或いは動的な荷重が継続して若しくは繰り返して印加される。このため、これらの部材には、多くの場合、上述した各種荷重によって生じる内部応力に充分耐え得る強度を備えた部材が選択され使用されている。そして、かかる場合においても、装置全体の定期点検や損傷し易い箇所の部材の交換の必要性が定められている場合が多い。

一方、上述した各種部材等を対象として、それらに印加される荷重によって実際に生じる内部応力の大きさやその変動を把握することは、当該部材の疲労度等を推測するに際し、又素材の改善や形状構造の研究に際しては重要な情報となる。
このため、従来より、特に超音波を用いてこれら各部材の内部応力を厳密に測定しようとする試みが、例えば下記特許文献１に示すように、比較的多くなされている。
特開平５−２０３５１３

しかしながら、上記特許文献１における特開平５−２０３５１３公報に記載の発明（超音波によるボルト軸力の測定方法）は、音速測定に際しては当該部材の音速を超音波の縦波と横波の両方で測定し特定することが必須の要件となっており、このため多くの手間がかかるばかりでなく、超音波センサとして縦波用超音波センサと横波用超音波センサの両方を常時準備しなければならないという煩わしさを常に伴っていた。
更に、この特許文献１にかかる発明は、縦波と横波の音速比Ｒと応力σとの関係に関して技術的な関連性が未開示であり、産業上の利用性に欠けるという不都合があった。

本発明の目的は、かかる従来例の有する不都合を改善し、特に、縦波又は横波の何れか一方の超音波を用いて被測定物の音速若しくは音の伝搬時間を測定すると共に、これらの情報に基づいて荷重負荷時における当該被測定物の内部に生ずる内部応力をリアルタイムで定量的に測定し得る応力測定方法およびその装置を提供することにある。

本発明にかかる応力測定方法およびその装置は、被測定物の内部応力の測定に際し、当該被測定物内を伝搬する音の伝搬時間を測定してその伝搬時間の変化情報を利用するようにした構成となっている。
そして、この場合、課題解決のための手法として、本発明では、被測定物の素材が備えている固有の常数（応力依存係数α）をまず算定し、これに基づいて被測定物内の内部応力をリアルタイムで具体的に測定するようにした点に特徴を有する。

まず、本発明にかかる応力測定方法（請求項１乃至４）は、それぞれ被測定物内に現実に発生している内部応力を、音の伝搬時間にかかるデータの変化を測定しこれに基づいて算定する場合の、手順にかかる発明であって、各請求項ともそれぞれ以下に示す四つの工程を基本的な構成として備えている。

第１の工程では、被測定物と同一の素材から成るテストピースについて、所定箇所の対向面相互間を伝搬する音の伝搬時間を、荷重を変化させながら印加して当該印加荷重σ_ｎ の大きさに対応した音の伝搬時間ｔ_ｎ と，荷重ゼロの時（無負荷時）の音の伝搬時間ｔ _０ とについて、それぞれ超音波を用いて測定する。

第２の工程では、上記第１の工程で得られる音の伝搬時間ｔ_０ ，ｔ_ｎ にかかる複数情報に基づいて該伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）を算定すると共に、これらの情報及び前述した被測定物の素材が備えている固有の弾性定数等に基づいて当該被測定物特有の応力依存係数αを、予め特定された次式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）により算定する。

第３の工程では、上記第１乃至第２の各工程に相前後して、装置若しくは建造物等に装備され未知の荷重が印加された状態の前記被測定物について前記所定箇所の対向面相互間を伝搬する音の伝搬時間ｔを測定する。
即ち、装備されたままの状態の被測定物について、その現実の音の伝搬時間ｔを測定する。この場合の測定のタイミングは、前述した第１乃至第２の各工程の前でも或いはその後であってもよい。

そして、第４の工程では、上記第１乃至第３の各工程で得られる諸情報に基づいて、被測定物についてその音の伝搬時間ｔの変化分Δｔ（＝ｔ−ｔ_０ ）を算定すると共に、前述した応力依存係数αを組み込んで予め特定された応力σに関する算定式σ＝ｆ（Δｔ，α）を演算し、実際に装備された状態の被測定物についてその現実の内部応力σを算定するようにした。

このため、本発明にかかる応力測定方法によると、第１工程で被測定物における音の伝搬時間ｔ_ｎ ，ｔ _０ 及び第２工程で算定した応力依存係数α等を被測定物の素材にかかる諸情報として導入すると共に、これらと第３工程で測定した被測定物の現実の音の伝搬時間ｔを取り込んで、予め特定された算定式σ＝ｆ（Δｔ，α）を演算し内部応力σを算定するようにした（第４工程）。このため、当該被測定物内に生じている内部応力σを容易に且つ迅速に算定することができる。
更に、測定した音の伝搬時間を直接利用するようにしたので、測定した音の伝搬時間から音速を算出する等の手間（演算）が不要となり内部応力σの算定時間の短縮化を図ることができる。

（削除）

（削除）

又、前述した請求項１記載の応力測定方法にあっては、第１の工程では、所定部材の音の伝搬時間の測定に際しては、所定部材に対する荷重σ_ｎ の印加方向と同方向に超音波を伝搬させるようにする。
この場合、前記第４の工程で設定した未知の内部応力σの算定式は、σ＝ｆ（Δｔ，α）として、σ＝｛Ｅ／（１−α）｝・（Δｔ／ｔ_０ ）が採用されている。
ここで、Δｔは、Δｔ＝ｔ−ｔ_０ であり、記号ｔは未知の荷重が付された状態の被測定物の測定箇所における音伝搬時間を示す。
このようにすると、被測定物内に生じている未知の内部応力σの大きさを、音の伝搬時間情報を使用して具体的に算定することができる。

更に、前述した請求項１に記載の応力測定方法において、第２の工程では、応力依存係数αを算定する算定式として、予め特定されたα＝｛１−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ ）｝、（ここで、ｎは音伝搬時間の測定回数、Ｅは被測定物の素材のヤング率、以下同じ）を採用し、これに基づいて前述した被測定物の素材の応力依存係数αが算定される。
ここで、応力依存係数αを算定する算定式としては、α＝｛１−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ ）｝に代えて、

を採用してもよい（請求項２）。
このようにすると、被測定物の素材固有の応力依存係数αが具体的に算定され、これに基づいて被測定物内に生じている未知の内部応力σの大きさが、音の伝搬時間情報を測定することでリアルタイムで具体的に且つ迅速に算定される。

これに対し、前述した請求項３に記載の応力測定方法にあっては、第１の工程では、被測定物の音の伝搬時間の測定に際しては、被測定物に対する荷重σ_ｎ の印加方向に直交する方向に超音波を伝搬させるようにする。
この場合、第４の工程で設定した内部応力σの算定式σ＝ｆ（Δｔ，α）として、
σ＝−｛Ｅ／（ν＋α）｝・（Δｔ／ｔ_０ ）、（ここで、νは被測定物の素材のポアソン比、以下同じ）が採用されている。ここで、Δｔは、Δｔ＝ｔ−ｔ_０ であり、記号ｔは未知の荷重が付された状態の被測定物の測定箇所における音伝搬時間を示す。

このようにしても、前述した請求項１の場合と同様に、被測定物内に生じている未知の内部応力σの大きさを、音の伝搬時間情報を使用して定量的にしかもリアルタイムで迅速に演算し算定することができる。
更に、荷重σ_ｎ の印加方向に対して直交する方向に超音波を伝搬させるようにしたので、例えば音の伝搬時間の測定に際しては超音波センサ等の装着領域を広範囲から選択することができ、かかる点において作業性能率の向上を期待することができる。

又、前述した請求項３に記載の応力測定方法において、第２の工程では、応力依存係数αを算定する算定式として、
α＝｛−ν−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ ）｝、を採用し、これに基づいて前述した被測定物の素材の応力依存係数αを算定するようにした。
このため、一定条件のもとに応力依存係数αが具体的に数値として特定されることから、前述した請求項１記載の応力測定方法の場合と同様に、これに基づいて被測定物内に生じている未知の内部応力σを、リアルタイムで且つ円滑に算定することが可能となる。
ここで、応力依存係数αを算定する算定式としては、α＝｛−ν−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ ）｝に代えて、

を採用してもよい。

ここで、上述した各請求項にかかる応力測定方法においては、所定部材の音の伝搬時間の測定に際し縦波超音波を使用した場合を前提として説明したが、横波超音波を使用して測定するようにしてもよい。

次に、本発明にかかる応力測定装置（請求項５乃至１０）では、被測定物の内部応力σの測定に際し、当該被測定物と同一の素材から成るテストピースにおける音の伝搬時間を印加荷重σ_ｎ 変化させながら測定しその変化情報ｔ _ｎ ，ｔ _０ （無負荷の場合）を利用するように構成した。
この場合、この請求項５乃至１０に記載の各応力測定装置は、それぞれ音伝搬時間測定手段と、応力依存係数αを算定する応力依存係数算定手段と、音伝搬時間にかかるデータおよび応力依存係数α等に基づいて被測定物の内部応力σを演算し算定する主演算部とを、共通の基本構成として備えている。

この内、音伝搬時間測定手段は、装置若しくは建造物等に装備され未知の荷重が印加された状態の被測定物について所定箇所の対向面相互間を伝搬する音の伝搬時間ｔを、超音波を用いて測定し得る機能（第１の音伝搬時間測定機能）を備えている。又、この音伝搬時間測定手段は、被測定物と同一の素材から成るテストピースに所定箇所の対向面相互間を伝搬する音の伝搬時間を、所定の荷重が変化させながら印加して当該印加荷重σ_ｎ の大きさとその変化に対応した音の伝搬時間ｔ_ｎ ，ｔ _０ （無負荷の場合）を、超音波を用いて測定し得る機能（第２の音伝搬時間測定機能）を備えたものとなっている。

ここで、この音の伝搬時間測定手段は、縦波超音波を使用したものであっても或いは横波超音波を使用したものであってもよい。
又、上記各部材の音の伝搬時間測定に際しては、必要に応じて、荷重σ_ｎ が、予め別に準備された荷重印加手段によって被測定物に対して印加されるようにしてもよい。

又、応力依存係数算定手段は、前述した音伝搬時間測定手段で測定される複数の各音伝搬時間情報ｔ，ｔ _０ ，ｔ_ｎ および当該音伝搬時間情報ｔ_ｎ に対応した印加荷重情報σ_ｎ を入力すると共にこの内の音伝搬時間情報ｔ_ｎ ，ｔ _０ 又はｔ，ｔ _０ を抽出する時間情報抽出機能と、この抽出された音伝搬時間情報ｔ_ｎ に基づいて当該音の伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）を算定する時間変化分演算機能と、この演算され算出された時間変化分Δｔ_ｎ 及び前述した印加荷重情報σ_ｎ 等に基づいて被測定物を形成する素材特有の応力依存係数αを予め特定された所定の算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）を演算し算定する常数算定機能とを備えている。

この応力依存係数算定手段は、更に、測定された情報の内の未知の荷重が印加された状態の被測定物について測定された音の伝搬時間ｔおよび当該部材にかかる素材の負荷荷重ゼロの場合の音の伝搬時間ｔ_０ に基づいて当該音の伝搬時間の変化分Δｔ（＝ｔ−ｔ_０ ）を算定する機能（他の時間変化分演算機能）を備えている。

更に、主演算部は、前述した音伝搬時間変化情報Δｔおよび前述した応力依存係数算定手段で算定される応力依存係数α等に基づいて、実際に装置若しくは建造物等に装備され未知の荷重が印加された状態の被測定物の内部応力σを、予め特定された算定式
σ＝ｆ（Δｔ，α）により演算しその結果を出力する応力算定機能を備えたものとなっている。

このため、本発明にかかる上記応力測定装置では、被測定物を形成する素材が備えている固有の応力依存係数αを算定することにより、この応力依存係数αおよび前述した現実の被測定物の音伝搬時間の変化分Δｔ等に基づいて当該被測定物内に生じている内部応力σをリアルタイムで迅速に演算することが可能となり、被測定物から測定された音の伝搬時間のデータを直接使用し得るように構成したので（伝搬時間のデータから音速を算定する演算が不要となり）被測定物の内部応力σの算定を、より迅速に成し得ることが可能となっている。

ここで、上述した音伝搬時間測定手段にあっては、前述した所定距離間の音の伝搬時間ｔ _ｎ ，ｔ _０ を、被測定物と同一の素材から成るテストピースについて測定することとした。

又、上述した応力測定装置において、音伝搬時間測定手段にあっては、被測定物に対する音伝搬時間測定用の超音波の伝搬方向を、当該被測定物に印加される荷重σ_ｎ の印加方向と同方向に設定する構成とし、更に、前述した主演算部が、次式
σ＝ｆ（Δｔ，α）＝｛Ｅ／（１−α）｝・（Δｔ／ｔ_０ ）、に基づいて前述した被測定物の内部応力σを具体的に演算する構成とした（請求項５乃至７）。
このようにすると、荷重印加方向と同方向に超音波を伝搬させた場合に得られる音の伝搬時間情報に基づいて、被測定物の未知の内部応力σを、測定条件に合わせて具体的に且つ迅速に演算し算出することが可能となる。

又、前述した応力依存係数算定手段を、具体的には、前述した音伝搬時間測定手段で得られる印加荷重σ_ｎ を含む複数の情報を入力すると共に当該入力された複数の情報から音伝搬時間情報ｔ_ｎ ，ｔ _０ を抽出する情報抽出機能と、この抽出された音伝搬時間情報ｔ_ｎ ，ｔ _０ に基づいて当該印加荷重σ_ｎ に対応した音伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）を算定する伝搬時間変化分算定機能とを有し、この音伝搬時間にかかる情報および前記印加荷重σ_ｎ にかかる情報等に基づいて応力依存係数αを算定する算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）として、次式、
α＝｛１−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ ）｝
を採用すると共に、これに基づいて前記被測定物にかかる素材固有の応力依存係数αを演算し算定する応力依存係数算定機能とを備えて構成した（請求項６）。
このため、荷重の印加方向と同方向に超音波を伝搬させた場合に得られる音の伝搬時間情報に基づいて、当該被測定物の応力依存係数αを有効に且つ具体的に算定することができる。

ここで、前述した応力依存係数算定手段では、応力依存係数αを算定する算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）として、α＝｛１−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ ）｝に代えて、次式

とし、これに基づいて被測定物にかかる素材固有の応力依存係数αを演算し算定するようにしてもよい（請求項７）。
このようにすると、荷重の印加方向と同方向に超音波を伝搬させた場合に得られる音の伝搬時間情報に基づいて当該被測定物の応力依存係数αを高精度に且つ有効に算定することができ、ひいては、被測定物の内部応力σの算定精度を向上させることができる点で都合がよい。

更に、前述した音伝搬時間測定手段にあっては、被測定物に対する音伝搬時間測定用の超音波の伝搬方向を、当該被測定物に印加される荷重σ_ｎ の印加方向に対して直交する方向に設定する構成とすると共に、
前述した主演算部が、次式
σ＝ｆ（Δｔ，α）＝−｛Ｅ／（ν＋α）｝・（Δｔ／ｔ_０ ）
に基づいて当該被測定物の未知の内部応力σを具体的に演算するように構成してもよい（請求項８乃至１０）。

このようにしても、前述した応力測定装置と同等に機能するほか、荷重σ_ｎ の印加方向に対して直交する方向に超音波を伝搬させるようにしたので、例えば音の伝搬時間の測定に際しては超音波センサ等の装着領域を広範囲から選択することができ、かかる点において作業性能率の向上を期待することができる。

ここで、前述した請求項８に記載の応力測定装置において、応力依存係数算定手段を、前述した音伝搬時間測定手段で得られる印加荷重σ_ｎ を含む複数の情報を入力すると共に当該入力された複数の情報から音伝搬時間情報ｔ_ｎ ，ｔ _０ を抽出する情報抽出機能と、この抽出された音伝搬時間情報ｔ_ｎ ，ｔ _０ に基づいて当該印加荷重σ_ｎ に対応した音伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）を算定する伝搬時間変化分算定機能とを有し、
この音伝搬時間にかかる情報および前記印加荷重σ_ｎ にかかる情報等に基づいて応力依存係数αを算定する算定式として、α＝｛−ν−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ ）｝、を採用すると共に、これに基づいて被測定物の素材固有の応力依存係数αを演算し算定する応力依存係数算定機能を備えた構成としてもよい。
このようにしても、荷重印加方向に直交する方向に超音波を伝搬させた場合に得られる音の伝搬時間情報に基づいて当該被測定物の応力依存係数αを、具体的に且つ有効に算定することができる。

ここで、前述した請求項８又は９に記載の応力依存係数算定手段では、応力依存係数αを算定する算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）として、α＝｛−ν−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ ）｝に代えて、次式

を採用し、これに基づいて前記被測定物の素材固有の応力依存係数αを演算し算定するようにしてもよい（請求項１０）。

このようにしても、荷重印加方向に直交する方向に超音波を伝搬させた場合に得られる音の伝搬時間情報に基づいて当該被測定物の応力依存係数αを有効に算定することができ、特に測定した音伝搬時間にかかる情報の全体を対象として平均化したので、被測定物を構成する素材固有の応力依存係数αを音の伝搬時間情報に基づいて迅速に且つ高精度に算定することができ、ひいては、被測定物の内部応力σの算定精度を向上させることができる。

本発明では、上述したように構成され機能するので、これによると、超音波を用いて被測定物の音の伝搬時間を測定すると共に、これらの情報に基づいて当該被測定物の新たに特定された素材固有の応力依存係数αが算定されることから、この算定された応力依存係数αを含むこれらの諸情報に基づいて荷重負荷時における当該被測定物の内部に生ずる内部応力をリアルタイムで定量的に且つ迅速に特定することができ、更に、被測定物の素材が備えている固有の応力依存係数αの算定に際しては、テストピース等の他の部材を不要としたので、その分、容易に且つ迅速に、被測定物の内部に生じている内部応力を特定し得るという従来にない優れた応力測定方法およびその装置を提供することができる。

以下、本発明の最良の実施形態を添付図面に従って説明する。

この実施形態は、図１乃至図２に示すように、被測定物１について、その音の伝搬時間を印加荷重の変化に対応して測定する音伝搬時間測定手段１１と、ここで得られた音の伝搬時間および当該伝搬時間にかかる種々の情報に基づいて応力依存係数αの値を算定する応力依存係数算定手段（応力依存係数算定回路部）１２と、この算定された応力依存係数αの値および前述した音の伝搬時間にかかる種々の情報に基づいて予め特定された所定の関数を演算し被測定物１の内部応力σを算定する主演算部１３とを備えている。

前述した音伝搬時間測定手段１１には、本実施形態では荷重負荷設定手段としての荷重印加装置２１が併設されている。又、図１において、符号１４は情報記憶部を示し、符号１５は必要に応じて所定の情報又は動作指令等を入力する入力部を示す。更に、符号１６は、主演算部１３で得られた被測定物１の内部応力等の情報を表示し又は外部にプリント出力するための出力部を示す。

以下、これを更に詳細に説明する。

〔音伝搬時間測定手段について〕
図１乃至図２において、音伝搬時間測定手段１１は、本実施形態では、製品や素材或いは装置若しくは建造物等に装備され未知の荷重が印加された状態の被測定物１について所定箇所の対向面相互間Ｌを伝搬する音の伝搬時間ｔを超音波を用いて測定する第１の音伝搬時間測定機能を備えている。又、この音伝搬時間測定手段１１は、被測定物１の素材について前記所定箇所の対向面相互間Ｌ_０ を伝搬する音の伝搬時間を無負荷時若しくはこれに近い状態の伝搬時間ｔ_０ と所定の荷重を変化させながら印加して当該印加荷重σ_ｎ の大きさに対応した音の伝搬時間ｔ_ｎ とについてそれぞれ超音波を用いて測定する第２の音伝搬時間測定機能を備えている（図３参照）。

この内、音伝搬時間測定手段１１の第１の音伝搬時間測定機能は、製品や素材或いは装置若しくは建造物等に装備され未知の荷重が印加された状態の被測定物１について、例えば図２に示すように、所定箇所の対向面相互間（一方の面と他方の面の相互間であって、音響的に対向する面のこと、以下同じ）Ｌを伝搬する音の伝搬時間ｔを超音波を用いて測定するという形態をもって実現される。
この図２に示す被測定物１については、本実施形態では２枚の鋼板２Ａ，２Ｂを連結するボルトの場合を想定したもので、当該ボルトである被測定物１が締結された状態では当該被測定物１には未知の荷重は引張荷重となって存在し、その締め付けの度合いによってボルト内に所定の内部応力が生じることとなる。符号３ａ，３ｂはそれぞれ異なったボルト（被測定物１，１）に螺合され前述した２枚の鋼板２Ａ，２Ｂを固定するためのナットを示す。

この場合、測定に際しては、例えば当該ボルト（被測定物１）のナット３ｂを緩めて、その締結力をゼロに近い状態に設定し、かかる状態での音の伝搬時間ｔ_０ を測定する。かかる測定が不可能な場合（ナット３ｂを緩めることができない場合）は、後述する第２の音伝搬時間測定機能によって別途測定される音の伝搬時間にかかる情報に基づいて、当該荷重ゼロの場合の音の伝搬時間ｔ_０ が換算され使用される。
そして、この音伝搬時間測定手段１１の第１の音伝搬時間測定機能が作動して得られる音の伝搬時間にかかる情報ｔ，ｔ_０ に基づいて、後述するように、未知の荷重が印加された状態（実装状態）における被測定物１について、その内部応力σが演算され算定されるようになっている。

又、音伝搬時間測定手段１１の第２の音伝搬時間測定機能は、本発明とは別に、テストピースを用いて測定する場合もあり、参考までに図３に例示する。
即ち、この場合は、例えば図３に示すように、前述した被測定部材１（ボルト）と同一の素材からなる他の部材（テストピース）１Ａについて、所定箇所の対向面相互間の距離Ｌ_ＴＰを伝搬する音の伝搬時間を、無負荷時の伝搬時間ｔ_０ と当該他の部材１Ａに対して所定の荷重を変化しつつ印加された状態の当該印加荷重σ_ｎ の大きさに対応した音の伝搬時間ｔ_ｎ とについて、それぞれ超音波を用いて測定するという形態をもって実現される。
この図３において、他の部材（テストピース）１Ａは、長さがＬ_ＴＰ、両端部の直径がＤ、中央部が一辺の幅Ｗ_ＴＰの正四角形状に形成されたものが使用されている。この場合、長さＬ_ＴＰの内、両端部の断面円形部分がそれぞれ全体の１／５程度を占め、中央部の正四角形状部分が全体の２／５程度を占めている。
ここで、テストピース１Ａの中央部の断面形状を、本実施形態（図３）では正四角形状としたが矩形であってもよい。

この場合、音の伝搬時間の測定に際して特定される対向面相互間の距離の設定箇所としては、この図３では荷重の印加方向と同方向に沿って距離Ｌ_ＴＰを設定し、荷重の印加方向に直交する方向に沿って距離Ｗ_ＴＰを設定し、そのいずれの箇所でも測定し得るようになっている。
そして、後述する応力依存係数算定回路部１２では、この音伝搬時間測定手段１１の第２の音伝搬時間測定機能が作動して得られる音の伝搬時間にかかる情報ｔ_ｎ ，ｔ_０ に基づいて、被測定物１を形成する素材の応力依存係数αを演算され算定されるようになっている。

この音伝搬時間測定手段１１は、具体的には図１乃至図３に示すように、所定周波数の超音波を前述した被測定部材１内に伝搬させてその端面からの反射波を受信する超音波送受波器１１ａ（又は１１ｂ）と、この超音波送受波器１１ａ（又は１１ｂ）で受信される反射超音波の受信のタイミングを音伝搬時間にかかる情報として外部出力する超音波用電気信号送受信回路１１Ａと、この超音波用電気信号送受信回路１１Ａから出力される音伝搬時間にかかる情報を入力して当該超音波の送信時から受信時までの時間を算定する音伝搬時間算定回路部１１Ｂとを備えている。符号１１Ｄは、二つの超音波送受波器１１ａ又は１１ｂと超音波用電気信号送受信回路１１Ａとを連結する切り替えスイッチを示す。

ここで、超音波送受波器１１ａは印加荷重σ_ｎ の印加方向に沿った方向に超音波を送受信するための超音波送受波器を示し、超音波送受波器１１ｂは印加荷重σ_ｎ の印加方向に直交する方向に超音波を送受信するための超音波送受波器を示す。
又、前述した音伝搬時間算定回路部１１Ｂは、更に、後述する荷重負荷設定手段１２から出力され前述した他の部材１Ａに対して印加される荷重情報σ_ｎ を入力すると共に当該荷重の大きさ（ゼロ荷重σ_０ を含む）に対応した音伝搬時間ｔ_０，ｔ_１，ｔ_２，ｔ_３，〜ｔ_Ｋ として集計し整理する伝搬時間整理機能を備えている。

一方、高圧プラントや橋梁等に実装されている状態（所定の荷重が印加された状態）の被測定部材１にあっては、その無負荷時の伝搬時間ｔ_０ の測定が事実上不可能な場合がある。
この場合、当該被測定部材１の長さＬがテストピース１Ａの長さＬ_ＴＰと同一であれば、前述したテストピース１Ａについて測定した無負荷時の伝搬時間ｔ_０ をそのまま流用される。
これに対して、Ｌ≠Ｌ_ＴＰの場合は、前述したテストピース１Ａについて距離Ｌ_ＴＰで測定した無負荷時の伝搬時間ｔ_０ が、比例計算で当該被測定部材１自体の無負荷時の伝搬時間ｔ_０ （対向面相互間Ｌ_ＴＰに対応した無負荷時の伝搬時間ｔ_０ ）に換算して使用されるようになっている。これらの演算処理は、音伝搬時間測定手段１１の音伝搬時間算定回路部１１Ｂが備えている伝搬時間換算機能によって実行されるようになっている。そして、かかる一連の動作は後述する入力部１５からの指令に基づいて実行される。
そして、この種の情報（即ち、被測定部材１１にかかる対向面相互間の距離Ｌおよび被測定部材と同一の素材からなる他の部材１Ａにかかる対向面相互間の距離Ｌ_ＴＰの情報）等は、予め後述する情報記憶部１４に入力され記憶されるようになっている。

尚、音伝搬時間測定手段１１では、被測定物１を形成する素材の応力依存係数αを演算するために必要な印加荷重σ_ｎ の変化と対応する音の伝搬時間に関する情報ｔ_ｎ とを被測定部材と同一の素材からなる他の部材１Ａについて測定した場合を例示したが、図２に例示した被測定物１そのものを（一時的に装備箇所から離脱させて）図３に示す荷重印加装置２１に直接装着して当該被測定物１で直接測定してもよい。
又、新たに所定の建造物等に装備する被測定物１については、装備前に当該被測定物１を前述した図３に示す荷重印加装置２１に直接装着して当該被測定物１で直接測定し、これによって、予め当該素材の応力依存係数αを後述する式（５）又は（１０）に基づいて算定して置いてもよい。

〔荷重負荷設定手段について〕
本実施形態では、前述した音伝搬時間測定手段１１には、荷重負荷設定手段としての荷重印加装置２１が併設されている（図３参照）。この荷重印加装置２１は、前述した被測定物と同一の素材からなる他の部材１Ａに対して（或いは被測定物そのものに対して）所定の引っ張り荷重を印加するためのもので、当該他の部材１Ａにかかる音の伝搬時間の測定に際して印加荷重を変化させることが可能な構造となっている。

即ち、この荷重負荷設定手段（荷重印加装置２１）は、前述した被測定物と同一の素材からなる他の部材（テストピース）１Ａに対して外部から所定の荷重負荷σ_ｎ を連続的に又は段階的に変化させながら印加して、前述した音伝搬時間測定手段１１による音伝搬時間測定時に、当該他の部材（テストピース）１Ａに所定の内部応力を生じさせる機能を備えている。
ここで、この外部からの荷重負荷σ_ｎ に対応して、試料（他の部材１Ａ）内には所定の内部応力が生じる。この場合、試料（他の部材１Ａ）の内部に生じる応力σの大きさは、一般に、印加される荷重負荷σ_ｎ の大きさにほぼ等しいとして扱われる。

ここで、上述した荷重印加装置２１の構成について説明する。
図３において、この荷重印加装置２１は、基台２１Ａと、この基台２１Ａ上に所定間隔を隔てて装備された二本のガイド支持体２１Ｅ，２１Ｆと、この各ガイド支持体２１Ｅ，２１Ｆが備えているガイド支柱２１ｅ，２１ｆに沿って上下動可能に装備された移動体２１Ｇ，２１Ｈと、この各移動体２１Ｇ，２１Ｈをその図３における上端部で一体的に連結する横架支持体２１Ｋと、この横架支持体２１Ｋの中央部に下方に向けて垂下装備された一方の試料保持機構２２と、この一方の試料保持機構２２の下方で当該試料保持機構２２の対向して配設された他方の試料保持機構２３とを備えている。この他方の試料保持機構２３はその下端部が前述した基台２１Ａに固定装備された状態となっている。

この内、一方と他方の各試料保持機構２２，２３は、各中心軸線が同一中心軸線上に沿って配設され、前述した移動体２１Ｇ，２１Ｈが上下動した場合でも、その軸線にぶれが生じない状態で一体的に上下動し得るように装備されている。又、この一方と他方の各試料保持機構２２，２３は、それぞれが前述した他の部材１Ａの端部を同時に又は個別に着脱自在に固着し得るような構造となっている。
ここで、符号１１Ｐ_１ は超音波送受波器１１ａ用の中継端子を示し、符号１１Ｐ_２ は超音波送受波器１１ｂ用の中継端子を示す。

又、上記基台２１Ａには、前述した横架支持体２１Ｋに対してその上下動を付勢する上下動付勢装置（図示せず）が装備されている。
この上下動付勢装置は、本実施形態では油圧を利用した油圧式駆動機構を備え、これによって連続して又は段階的に所定の引っ張り荷重σ_ｎ を各試料保持機構２２，２３を介して前述した他の部材１Ａ等に印加し得るようになっている。
更に、前述した基台２１Ａには、上下動付勢装置（図示せず）の動作を外部操作によって制御する操作制御部２４が装備され、又、この上下動付勢装置によって設定される前述した他の部材１Ａに対する引っ張り荷重の大きさを表示する荷重表示部２５が併設されている。符号２６は荷重情報σ_ｎ を出力するための出力端子を示す。

このため、操作制御部２４を介して上下動付勢装置の動作を制御することにより、前述した横架支持体２１Ｋに対して図３の矢印Ａ方向に向けて所定の荷重を任意に設定制御することか可能となる。これにより、各試料保持機構２２，２３によって保持された他の部材１Ａに対しては、所定の大きさの引っ張り荷重（負荷荷重）を連続的に又は段階的に変化させて設定すると共に、当該設定された負荷荷重σ_ｎ を荷重情報として外部出力し得るようになっている。

そして、この場合、この引っ張り荷重（負荷荷重）の設定は、前述した音伝搬時間測定手段１１による音伝搬時間測定を想定しながら行われる。
即ち、音伝搬時間測定手段１１による他の部材（テストピース）１Ａ等に対する音伝搬時間測定は、例えば、負荷荷重がゼロσ_０ の場合の音伝搬時間ｔ_０ 、負荷荷重σ_ｎ （ここで、ｎ＝１，２，３，…，ｋ、以下同じ）を段階的に変化させることによって音伝搬時間の変化を測定値ｔ_ｎ として測定することにより実行される。
ここで、上記荷重印加装置２１はテストピース１Ａに対して引張り荷重を印加する場合について例示したが、圧縮荷重を印加するように構成されたものを使用してもよい。
又、上記荷重印加装置２１にあっては、負荷荷重の発生源として油圧機器を利用した場合を例示したが、同等に機能するものであれば他の荷重発生手段を使用したものであってもよい。

〔応力依存係数算定手段について〕
図１において、符号１２は応力依存係数算定手段としての応力依存係数算定回路部を示す。この応力依存係数算定回路部１２は、前述した音伝搬時間測定手段１１で測定される複数の各音伝搬時間情報ｔ，ｔ_０ ，ｔ_ｎ および当該音伝搬時間情報ｔ_ｎ （ｎ＝１，２，３，〜，ｋ、以下同じ）に対応した印加荷重情報σ_ｎ を入力すると共にこの内の複数の音伝搬時間情報ｔ_０ ，ｔ_ｎ を抽出する時間情報抽出機能と、この抽出された複数の音伝搬時間情報ｔ_０ ，ｔ_ｎ に基づいて当該音の伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）を算定する時間変化分演算機能と、これら音の伝搬時間にかかる諸情報および前述した印加荷重情報等に基づいて被測定物１にかかる素材固有の応力依存係数αを、予め特定された次式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）で演算し算定する応力依存係数算定機能とを備えている。

又、この応力依存係数算定回路部１２は、前述した伝搬時間の変化分Δｔ_ｎ の変化の割合を示す伝搬時間変化比（Δｔ_ｎ ／ｔ_０ ）を演算する時間変化比演算機能を備えている。記号ｔ_０ は印加荷重情報σ_ｎ ＝０の場合の前述した部材１Ａの音伝搬時間を示す。

そして、この応力依存係数算定回路部１２で算定された応力依存係数αおよび使用された各種情報、更には前述した音伝搬時間算定回路部１１Ｂで整理され又は演算された伝搬時間にかかる各種情報は、予め別に装備された情報記憶部１４に送り込まれ記憶されるようになっている。

ここで、素材固有の応力依存係数αを算定する演算式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）に関して、発明者が導き且つ特定した内容を、以下詳細に説明する。

《荷重の印加方向に直交する方向（横方向）に超音波を伝搬させる場合》
この場合の超音波（縦波）の伝搬時間ｔ_ｎ にかかる関数は、Ｗ_ＴＰを図３に示す音の伝搬距離（荷重付加＝０の場合）とすると、次式で表される。
ｔ_ｎ ＝Ｗ_ＴＰ（１−σ_ｎν／Ｅ）／Ｖ_０（１＋α_Ｔ２σ_ｎ／Ｅ）………（１）
この式（１）は、図３では、超音波送受波器１１ｂを作動させて音の伝搬時間を測定した場合に相当する。この式（１）で、分子の部分は横方向の長さを表す。このため、この式（１）の分子の部分は引っ張り応力（負荷）の増加と共に減少する。記号Ｖ_０ は、荷重ゼロのときに当該素材内を伝搬する音の速度を示す。又、σ_ｎ ／Ｅは縦方向のひずみを示す（一般式σ＝Ｅ・εより）。更に、νは当該素材のポアソン比を示し、Ｅはヤング率を示す。縦方向のひずみεにポアソン比νを掛けたものが横方向のひずみとなる。又、α_Ｔ２は被測定物を形成する素材の応力依存係数を示す。

この式（１）で分母は縦波の音速を意味する。この式（１）を展開して微小項を無視し、一部置き換えをし、α_Ｔ２σ_ｎ ／Ｅ≪１、 Ｅ＞σ_ｎ 、Ｗ_ＴＰ／Ｖ_０ ＝ｔ_０ より、次式（２）を得る。
ｔ_ｎ ≒ｔ_０｛１−（σ_ｎ ／Ｅ）（ν＋α_Ｔ２）｝…………………… （２）
これより、
（ｔ_ｎ ／ｔ_０ ）−１＝Δｔ_ｎ ／ｔ_０
＝−（σ_ｎ ／Ｅ）（ν＋α_Ｔ２）…………………… （３）

この式（３）より、応力依存係数α_Ｔ２を算定する演算式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）として、次式を得ることができる。
α＝α_Ｔ２＝−ν−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ）………………… （４）
ここで、ｎは音伝搬時間の測定回数を示す。音伝搬時間の測定は、実際には図３において、荷重σ_ｎ をσ_１ ，σ_２ ，σ_３ ，……σ_ｋ の如く徐々に変化させ、その都度、超音波を伝搬させて音伝搬時間の変化ｔ_１ ，ｔ_２ ，ｔ_３ ，……ｔ_ｋ を測定する。
そして、仮に、精度の高い測定値が得られれば、荷重ゼロ（σ_ｎ ＝０）の場合の音伝搬時間ｔ_０ を基準として、音伝搬時間ｔ_ｎ の何れか一つの音伝搬時間データに基づいて式（４）を演算することにより、どの測定値からでも同一の応力依存係数α_Ｔ２が算定される。

これに対して、測定値と誤差は切り離せないものであることから、多くは下式（５）の手法（若しくはこれと同等の手法）をもって、測定値全体を対象としてその平均値を求め、これにより測定誤差を必要最小限に抑えるようにして、応力依存係数αが算定される。

この応力音速定数α_Ｔ２の具体的な数値計算等については後述する。

《荷重の印加方向と同方向（縦方向）に超音波を伝搬させる場合》
この場合の超音波（縦波）の伝搬時間ｔ_ｎ にかかる関数は、音の伝搬距離をＬ_ＴＰとすると、次式（負荷荷重＝０の場合）で表される。
ｔ_ｎ ＝Ｌ_ＴＰ（１＋σ_ｎ／Ｅ）／｛Ｖ_０（１＋α_Ｔ１σ_ｎ／Ｅ）｝……（６）
この式（６）は、図３では、超音波送受波器１１ａを作動させて音の伝搬時間を測定した場合に相当する。この式（６）の分子の部分は縦方向の長さを表す。このため、応力（負荷）の増加と共に増加する。

又、この式（６）で、分母は縦波の音速を意味する。続いて、式（６）を展開して微小項を無視し、一部置き換えをする。即ち、α_Ｔ１σ_ｎ ／Ｅ≪１、 Ｅ＞σ_ｎ 、Ｌ_ＴＰ／Ｖ_０ ＝ｔ_０ より、次式（７）を得る。
ｔ_ｎ ≒ｔ_０｛（１＋σ_ｎ ／Ｅ）（１−α_Ｔ１）｝…………………… （７）
これより、
（ｔ_ｎ ／ｔ_０ ）−１＝Δｔ_ｎ ／ｔ_０
＝（σ_ｎ ／Ｅ）（１−α_Ｔ１）…………………… （８）
この式（８）より、応力依存係数αを算定する演算式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）として、次式を得ることができる。
α＝α_Ｔ１＝１−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ）…………………… （９）

ここで、この式（９）についても前述した式（４）の場合と同様の理由から、多くは下式（１０）の手法（若しくはこれと同等の手法）をもって、測定値全体を対象としてその平均値を求め、これにより、測定誤差が最小限に抑えられた状態で応力依存係数αが算定される。

………………（１０）
この応力依存係数α_Ｔ１の具体的な数値計算等については後述する。
尚、上述したように、測定方法が異なると同一素材についても応力依存係数αの算定式は異なった形態となる。

〔主演算部（特に内部応力σの算定式）について〕
図１に開示した主演算部１３では、実際に装置若しくは建造物等に装備され未知の荷重が印加された状態の被測定物１の内部応力σを、予め特定された算定式σ＝ｆ（Δｔ，α）によって演算し算定する内部応力算定機能を備えている。
ここで、Δｔは、Δｔ＝ｔ−ｔ_０ を示す。記号ｔは未知の荷重が印加された状態で測定された被測定物１の所定距離間の音の伝搬時間を示す。

この主演算部１３の内部応力算定機能は、前述した応力依存係数算定回路部（応力依存係数算定手段）１２で算定される応力依存係数α_Ｔ２（又はα_Ｔ１）および前述した音伝搬時間測定手段で測定される複数の音伝搬時間情報等に基づいて実行され、未知の荷重が印加された状態の被測定物１の内部応力σが算定されるようになっている。

この内部応力σを算定する演算式σ＝ｆ（Δｔ，α）については、本実施形態では、測定対象たる被測定物１に関してその素材の応力依存係数α_Ｔ２（又はα_Ｔ１）が予め明確に特定されていることを前提として成立したものとなっている。
この場合、図２に例示したように、未知の内部応力σに関する情報を収集するために、被測定物１の現実の音の伝搬時間ｔが超音波を用いて測定される。
ここで、音の伝搬方向の設定を、内部応力σが生じている方向と同方向にするか直交する方向とするかによって、内部応力σを算定する演算式σ＝ｆ（Δｔ，α）が異なったものとなる。かかる演算式の相違は、音の伝搬時間を特定する音の通路長さが、引っ張り方向に沿って伸びる場合、或いはそれと直交して縮む場合との違いによって生じる。

即ち、前述した応力依存係数αの算定式の場合と同様に、未知の荷重の印加方向に直交する方向に超音波を伝搬させて所定距離間の音の伝搬時間ｔを測定する場合と、荷重の印加方向と同方向に超音波を伝搬させて所定距離間の音の伝搬時間ｔを測定する場合とでは、異なった演算式が採用される。以下、これを具体的に説明する。

《未知の荷重の印加方向に直交する方向（横方向）に超音波を伝搬させる場合》
この場合の超音波（縦波）の伝搬時間ｔにかかる関数は、前述した応力依存係数α_Ｔ２を算出する場合（式（１））と同様に、次式で表される。
ｔ＝Ｗ（１−σν／Ｅ）／Ｖ_０（１＋α_Ｔ２σ／Ｅ）………………… （１１）
この式（１１）は、図３で考えると、超音波送受波器１１ｂを作動させて未知の荷重が印加された部材（図示せず）について音の伝搬時間ｔを測定する場合に相当する。
この式（１１）で、分子の部分は横方向の長さを表す。記号Ｗ（図示せず）は被測定物１における負荷応力がゼロの場合の横方向の音伝搬距離を示す（図６における記号Ｗに相当する）。又、α_Ｔ２は当該被測定物１の応力依存係数を、σは被測定物１内に発生している未知の内部応力を示す。その他の構成記号は前述した式（１）の場合と同様である。

又、この式（１１）で分母は縦波の音速を意味する。この式（１１）を展開して微小項を無視し、一部置き換えをし、α_Ｔ２σ／Ｅ≪１、Ｅ＞σ、Ｗ／Ｖ_０ ＝ｔ_０ より、次式（１２）を得る。
ｔ≒ｔ_０｛１−（σ／Ｅ）（ν＋α_Ｔ２）｝…………………………… （１２）
これより、
（ｔ／ｔ_０ ）−１＝Δｔ／ｔ_０
＝−（σ／Ｅ）（ν＋α_Ｔ２）………………………… （１３）
この式（１３）より、未知の内部応力σを算定する演算式σ＝ｆ（Δｔ，α）として、次式を得る。
σ＝−｛Ｅ／（ν＋α_Ｔ２）｝・（Δｔ／ｔ_０）……………………… （１４）

これにより、未知の外部荷重が負荷された状態の被測定物１であっても、その時の現実の音伝搬時間ｔを測定することにより（当該素材の固有定数である応力依存係数α_Ｔ２が判っていることを前提）、未知の内部応力σを、当該被測定物１が装備されたままの状態で、リアルタイムで定量的に且つ迅速に算出することが可能となる。

《未知の荷重の印加方向に同方向（縦方向）に超音波を伝搬させる場合》
この場合の超音波（縦波）の伝搬時間ｔにかかる関数は、前述した応力依存係数αを算出する場合（式（６））と同様に、次式で表される。ここで、Ｌ（図示せず）は被測定物１の縦方向の音の伝搬距離を示す。又、α_Ｔ１は当該被測定物１の応力依存係数を、σは被測定物１内に発生している未知の内部応力を示す。その他の構成記号は前述した式（１）の場合と同様である。

ｔ＝Ｌ（１＋σ／Ｅ）／｛Ｖ_０（１＋α_Ｔ１σ／Ｅ）｝……………… （１５）
この式（１５）は、図３では、超音波送受波器１１ａを作動させて音の伝搬時間を測定した場合に相当する。この式（１５）の分子の部分は縦方向の長さを表す。このため、応力（負荷）の増加と共に増加する。

又、この式（１５）で、分母は縦波の音速を意味する。続いて、式（１５）を展開し微小項を無視し、一部置き換えをする。即ち、α_Ｔ１σ／Ｅ≪１、Ｅ＞σ、Ｌ／Ｖ_０ ＝ｔ_０ より、次式（１６）を得る。
ｔ_ｎ ≒ｔ_０｛（１＋σ／Ｅ）（１−α_Ｔ１）｝……………………… （１６）
これより、
（ｔ／ｔ_０ ）−１＝Δｔ／ｔ_０
＝（σ／Ｅ）（１−α_Ｔ１）………………………… （１７）
この式（１７）より、未知の内部応力σ_Ｔ１を算定する演算式σ＝ｆ（Δｔ，α）として、次式を得る。
σ＝｛Ｅ／（１−α_Ｔ１）｝・（Δｔ／ｔ_０ ）……………………… （１８）

これにより、この式（１８）を用いて、前述した式（１４）の場合と同様に、未知の外部荷重が負荷された状態の被測定物１であっても、その時の現実の音伝搬時間ｔを測定することにより、当該素材の固有定数である応力依存係数α_Ｔ１が判っていることを前提として、未知の内部応力σをリアルタイムで定量的に算出することが可能となる。

〔全体的な動作、機能〕
次に、上記実施形態における全体的な機能，動作等を説明する。
本実施形態では、未知の荷重が負荷された状態の被測定物１について、まず、当該被測定物１の素材固有の応力依存係数αを算定し、これに基づいて当該被測定物１の内部応力σを、所定の式σ＝ｆ（Δｔ，α）を演算して特定するようにした。この場合、応力依存係数αの算定に際しては、被測定物１そのものを使用しても、或いは当該被測定物１と同一の素材からなる他の部材１Ａを使用してもよい。本実施形態では後者に依る。

この応力依存係数αの算定は、前述したように当該素材の各種条件下で測定された音の伝搬時間にかかる情報を必要とする。このため、本実施形態では、上述した構成を以下の手順で作動させて、収集された各種データに基づいて応力依存係数αを演算し、これに基づいて内部応力σを算定するようになっている。
ここで、応力依存係数αおよび内部応力σの算定式は、音の伝搬時間の測定条件によって異なることは前述した通りである。このため、最初に図１乃至図４に従って装置の全体的な動作を説明し、その後に、測定条件の相違により生じる細部の動作について説明するものとする。

まず、装置全体を稼働させ、応力依存係数αを測定する。具体的には、応力依存係数αを特定するための手順が、例えば被測定物１と同一の素材からなる他の部材（テストピース）１Ａについて、以下の如く実行される。
この場合、図３に示すように、被測定物１と同一の素材からなる他の部材１Ａに関し、その素材について、まず、所定箇所の対向面相互間を伝搬する音の伝搬時間が測定される。この音の伝搬時間の測定は、超音波送受波器１１ｂ（又は１１ａ）、これを駆動する超音波用電気信号送受信回路１１Ａ、および音伝搬時間算定回路部１１Ｂの連係動作によって実行される。スイッチ１１Ｄは、端子ｂを介して超音波送受波器１１ｂを選択し、又は端子ａを介して超音波送受波器１１ａを選択する。この場合、スイッチ１１Ｄの動作は、入力部１５形の指令によって作動するようになっている。

そして、この超音波用電気信号送受信回路１１Ａで得られる超音波の送受信のタイミングにかかる情報は、音の伝搬時間にかかる情報として、そのまま音伝搬時間算定回路部１１Ｂへ送り込まれ、ここで、音の伝搬時間ｔ_０ ，ｔ_ｎ が算定される（図４、ステップＳ１）。
即ち、上記音伝搬時間算定測定手段１１では、他の部材１Ａに対して外部からの印加荷重σ_ｎ を印加しながら、当該印加荷重がゼロの場合（無負荷時σ_０ ）もしくはこれに近い状態における音の伝搬時間ｔ_０ と、荷重をｎ回変化させながら印加して当該印加荷重σ_ｎ の大きさに対応した音の伝搬時間ｔ_ｎ とについて、それぞれ超音波を用いて測定する（第１の工程）。この印加荷重σ_ｎ はｋ回変化させ、これに対応して音の伝搬時間ｔ_ｎ をｋ回測定する（ｎ＝１，２，３，……，ｋ）。
この音の伝搬時間ｔ_０ ，ｔ_ｎ の測定に際しては、本実施形態では、所定箇所の対向面相互間として、図３における他の部材１Ａの中央部（断面が四角形状）の距離Ｗ_ＴＰ部分が予め設定され、この箇所での多重反射等の測定手法を用いて、音の伝搬時間の測定が実行される。

この第１の工程で音の伝搬時間ｔ_０ ，ｔ_ｎ にかかる複数情報が得られると、次に、これらに基づいて当該伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）、および音伝搬時間変化比Δｔ_ｎ ／ｔ_０ 等が測定値毎にｋ回演算される（図４、ステップＳ２）。続いて、これらの情報および前述した被測定物１を形成する素材固有の弾性定数等に基づいて、当該被測定物固有の応力依存係数αが、予め特定された算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）により算定される（第２の工程）。かかる算定動作は、応力依存係数算定回路部１２で実行される（図４、ステップＳ３）。

ここで、算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）は、音の伝搬時間測定時に、図３で超音波送受波器１１ｂを使用して超音波を印加荷重σ_ｎ の方向に直交する方向に伝搬させた場合には、前述したように、
α＝α_Ｔ２＝−ν−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ ）……………… （４）
又は、

……………（５）
が使用される。

一方、音の伝搬時間測定時に、図３で超音波送受波器１１ａを使用して超音波を印加荷重σ_ｎ の方向に沿った方向に伝搬させた場合には、前述したように、
α＝α_Ｔ１＝１−（Δｔ_ｎ ／ｔ_０）（Ｅ／σ_ｎ）……………………… （９）
又は、

……………（１０）
が使用される。

このため、被測定物の素材固有の応力依存係数αの測定に際しては、印加荷重σ_ｎ の方向に沿った方向又はこれに直交する方向にいずれの方向に超音波を伝搬させても、その音伝搬時間を有効に測定することができ、これに基づいて当該応力依存係数αを定量的に且つ迅速に算出することができ、かかる点においてその汎用化をはかることができる。
又、被測定物の素材固有の応力依存係数αを、被測定物１と同一の素材からなる他の部材１Ａで測定する場合には、装備箇所からの取り外し等が困難な被測定物に対しても、他の部材１Ａによって応力依存係数αを迅速に算定することができ、これにより当該被測定物１の内部応力σの測定をなし得ることから、より一層の汎用化が可能となる。
これに対し、被測定物１のそのものを利用して応力依存係数αを算出する手法では、被測定物１の測定箇所を特定できないものもあって対象範囲が限定されるものの、被測定物と同一の素材からなる他の部材を全く準備する必要がないので、応力依存係数αの特定に際しての作業を迅速になし得るという利点がある。

ここで、上述した各種情報（σ_ｎ ，ｔ_ｎ ，ｔ_０ ，Δｔ_ｎ ，Δｔ_ｎ ／ｔ_０ ）或いは算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）等は、予め情報記憶部１４で記憶され、若しくは必要に応じて新たに記憶されるようになっている。

上述した第１乃至第２の各工程に相前後して、装置若しくは建造物等に装備され未知の荷重が印加された状態の被測定物１について所定箇所の対向面相互間を伝搬する音の伝搬時間ｔが測定される（第３の工程）。
この場合の音の伝搬時間ｔは、未知の荷重が印加された状態の伝搬時間情報であり、かかる情報によって当該被測定物１の内部応力σが演算されることから、再現性のある伝搬時間ｔが得られるまで、実際には数回繰り返し測定される。この音の伝搬時間ｔの測定は、前述した音伝搬時間算定測定手段１１で実行される（図２参照）。
同時に、未知の荷重がゼロ又はゼロに近い状態における音の伝搬時間ｔ_０ も測定される。この場合、この音の伝搬時間ｔ_０（σ_ｎ ≒０の場合）が測定不可能な場合は、前述したように、図３で測定した時間データｔ_０を（比例計算で）距離換算して当該被測定物１における音の伝搬時間ｔ_０（σ_ｎ≒０の場合）として使用される。かかる演算は、前述した入力部１５からの指令に基づいて音伝搬時間算定回路部１１Ｂで実行される（図４、ステップＳ４）。

この第３の工程で実行される被測定物１についての音の伝搬時間ｔ，ｔ_０の測定は、図３に示すように、前述した第１の工程で測定される音の伝搬時間ｔ_ｎ，ｔ_０の測定とは異なった時間帯で測定される。この場合、測定のタイミングとしては、前述した第１乃至第２の各工程の前でも後でもよい。
そして、この第３の工程で得られる音の伝搬時間ｔおよびこれに対応する荷重ゼロ（若しくはゼロに近い荷重）の場合の音の伝搬時間にかかる情報ｔ_０については、前述した第２の工程の場合と同様にその変化分Δｔ（＝ｔ−ｔ_０）、および音伝搬時間変化比Δｔ／ｔ_０ 等が音伝搬時間算定回路部１１Ｂで算定され（図４、ステップＳ５）、その算定結果の情報Δｔ，Δｔ／ｔ_０ は、この第３の工程で得られる音の伝搬時間ｔ，ｔ_０ と共に、情報記憶部１３で記憶されると共に主演算部１３へ直接送り込まれるようになっている。

そして、上記第１乃至第３の各工程で得られる諸情報に基づいて、予め特定された内部応力σに関する次式 σ＝ｆ（Δｔ，α）、（ここで、Δｔは Δｔ＝ｔ−ｔ_０ ）を演算し、これによって実際に装備され未知の荷重が印加された状態の被測定物１について、その現実の内部応力σが特定される（第４の工程）。

この内部応力σの演算は、前述したように主演算部１３で実行される。この主演算部１３では、前述したように応力依存係数算定回路部１２で算定された応力依存係数αおよび被測定物１について測定された音の伝搬時間ｔ，ｔ_０ 等が入力され、これらに基づいて前述した所定の演算式が演算され、未知の荷重が印加された被測定物１の内部応力が算定される（図４、ステップＳ６）。
ここで、被測定物１の内部応力にかかる算定式σ＝ｆ（Δｔ，α）は、例えば音の伝搬時間測定時に、図３で超音波送受波器１１ｂを使用して超音波を荷重の印加方向に直交する方向に伝搬させた場合と同等に設定した場合には、前述したように、
σ＝−｛Ｅ／（ν＋α_Ｔ２）｝・（Δｔ／ｔ_０）……………………… （１４）
が使用される。

一方、音の伝搬時間測定時に、図３で超音波送受波器１１ａを使用して超音波を荷重の印加方向に沿った方向に伝搬させた場合と同等に設定した場合には、前述したように、
σ＝｛Ｅ／（１−α_Ｔ１）｝・（Δｔ／ｔ_０）………………………… （１８）
が使用される。
本実施形態では、図２に示すように超音波を荷重の印加方向に沿った方向に伝搬させて未知の荷重が印加された被測定物１の音の伝搬時間を測定するようにしたので、内部応力にかかわる算定式σ＝ｆ（Δｔ，α）としては式（１８）が使用される。

このように、本実施形態は上述したように構成され機能するので、これによると、被測定物１を形成する素材固有の応力依存係数αを算定することにより、この応力依存係数αおよび前述した現実の被測定物の音伝搬時間の変化分Δｔ等に基づいて当該被測定物内に生じている内部応力σをリアルタイムで迅速に演算することが可能となり、更に、被測定物から測定された音の伝搬時間のデータを直接利用し得るように構成したので（伝搬時間のデータから音速を算定する演算が不要となり）被測定物の内部応力σを、より迅速に算定することができ、応力依存係数αを具体的に算定し得るので、未知の荷重が印加された状態の被測定物の内部応力σを迅速に且つ定量的に算定し得るという優れた応力測定方法及びその装置を得ることができる。

上述した実施形態において、未知の荷重が印加された状態の被測定物の内部応力σの測定に際しては、事前に、被測定物の素材固有の応力依存係数αを的確に且つ迅速に特定することが重要となっている。
ここでは、応力依存係数αを特定する場合の具体例を示す。
最初に、図３に開示した手法を用いて音の伝搬時間を測定し、これに基づいて当該部材の応力依存係数αを、式（４），（５）および式（９），（１０）に基づいて算定してみた。以下、これを説明する。

〔応力依存係数αの算定式の内容〕
（１）．印加荷重σ_ｎ の印加方向に直交する方向（横方向）に超音波を伝搬させて測定した場合
この場合の応力依存係数αは、前述したように、
α＝α_Ｔ２＝−ν−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ）………………… （４）
この式（４）をＹ＝ａＸ±ｂの形態に直してＸーＹ座標と同等の座標上に表示し、応力依存係数αの性質（定数か否か）をまず確認する。
式（４）より、Ｅ（Δｔ_ｎ ／ｔ_０ ）＝−（ν＋α）σ_ｎ …………………（４ａ）
具体的には、まず、ＸーＹ直交座標を想定して上式（４ａ）の左辺Ｅ（Δｔ_ｎ ／ｔ_０ ）をＹ座標上に、左辺のσ_ｎ をＸ座標上にとると、図５（Ａ）に示す線図となる。
この図５（Ａ）から明らかのように、「−（ν＋α）」は方向係数である。又、νはポアッソン比、αは２次および３次の弾性定数からなる定数（応力依存係数）である。従って、この方向係数−（ν＋α）は材料固有の値をとる。このため、今後の課題としては多くの実用材料についてα値の実測が望まれる。

（２）．印加荷重σ_ｎの印加方向に沿った方向（縦方向）に超音波を伝搬させて測定した場合
この場合の応力依存係数αは、前述したように、
α＝α_Ｔ１＝１−（Δｔ_ｎ ／ｔ_０ ）（Ｅ／σ_ｎ）…………………………（９）
この式（９）をＹ＝ａＸ±ｂの形態に直してＸーＹ座標と同等の座標上に表示して、応力依存係数αの性質（定数か否か）をまず確認する。
式（９）より、Ｅ（Δｔ_ｎ ／ｔ_０ ）＝（１−α）σ_ｎ ……………………（９ａ）
具体的には、上述した式（４ａ）の場合と同様に、まずＸーＹ直交座標を想定して上式（９ａ）の左辺Ｅ（Δｔ_ｎ ／ｔ_０ ）をＹ座標上に、左辺のσ_ｎ をＸ座標上にとると、図５（Ｂ）に示す線図となる。
この図５（Ｂ）から明らかのように、「（１−α）」は方向係数である。又、αは２次および３次の弾性定数からなる定数（応力依存係数）である。従って、この方向係数（１−α）は材料固有の値をとる。このため、今後の課題としては多くの実用材料につき、この場合のα値の実測が望まれる。

〔応力依存係数αの算定式における「Δｔ_ｎ ／ｔ_０ 」と「σ_ｎ」の関連性〕
（１）．測定条件
被測定物と同一の素材からなる部材１Ａの形状および寸法を図６（Ａ）に示す。長さＬは１８０〔ｍｍ〕、中央部の四角形状の断面は一辺が３０〔ｍｍ〕×１５〔ｍｍ〕角に精密加工されている。又、この部材の素材としてＳ２０Ｃの鋼材を使用した。この鋼材（Ｓ２０Ｃ）固有のヤング率Ｅ等の弾性定数を図６（Ｂ）の図表に示す。

（２）．図表化
「Δｔ_ｎ ／ｔ_０ 」と「σ_ｎ」との対応関係を図表化する。図７にこれを示す。
この図７で、表示「横方向」の欄は、印加荷重σ_ｎ の印加方向に直交する方向に縦波超音波を伝搬させて得られる時間データに基づいて算定した「Δｔ_ｎ ／ｔ_０ 」の値を示す。又、表示「縦方向」の欄は、印加荷重σ_ｎ の印加方向に沿って縦波超音波を伝搬させて得られる時間データに基づいて算定した「Δｔ_ｎ ／ｔ_０ 」の値を示す。この場合、印加荷重（負荷応力）σ_ｎ については、これを５．４４〔ＭＰａ〕ずつ変化させて縦波超音波の伝搬時間ｔ_ｎ を測定し、Δｔ_ｎ ＝ｔ_ｎ −ｔ_０ ，Δｔ_ｎ ／ｔ_０ を算定した。

（３）．グラフ化
「Δｔ_ｎ ／ｔ_０ 」と「σ_ｎ」との関係をグラフ化する。図８にこれを示す。
この内、図８（Ａ）は、印加荷重σ_ｎ の印加方向に直交する方向に縦波超音波を伝搬させて音の伝搬時間を測定した場合で、測定データとしては図７の「横方向」の欄に記載のデータがこれに対応する。この図８（Ａ）では、一つ一つのドットがσ_ｎとΔｔ_ｎ ／ｔ_０ との対のデータである。これらのドットを繋いだのがジグザグの線である。そして、このジグザグの線からなるσ_ｎとΔｔ_ｎ ／ｔ_０ とのデータ群を最小二乗法で処理した結果の直線Ｆ_２がジグザグの線からなるデータ群を図示の如く左上から右下に貫いている。そして、この図８（Ａ）の実測値に基づいた線図Ｆ_２が図５（Ａ）の線図に対応したものとなっている。

又、図８（Ｂ）は、印加荷重σ_ｎの印加方向に沿って縦波超音波を伝搬させて音の伝搬時間を測定した場合で、測定データとしては図７の「縦方向」の欄に記載のデータがこれに対応する。これらのドットを繋いだのが比較的滑らかなジグザグの線である。そして、このジグザグの線からなるσ_ｎとΔｔ_ｎ ／ｔ_０ とのデータ群を最小二乗法で処理した結果の直線Ｆ_１がジグザグの線からなるデータ群を図示の如く左下から右上に貫いている。そして、この図８（Ｂ）の実測値に基づいた線図Ｆ_１が図５（Ｂ）の線図に対応したものとなっている。この図８（Ｂ）の実測値では方向係数が正である。
この図８（Ａ）（Ｂ）で得られる傾斜直線の直線Ｆ_２ ，Ｆ_１ が、当該測定値の平均値として得られる実験式（一次関数）を示す。

（４）．Ｓ２０Ｃの鋼材にかかる応力依存係数αの算定例
４−１） ．応力依存係数α_Ｔ２（図８（Ａ）参照）の場合
図８（Ａ）で明らかのように、測定値の全体的な傾向（最小二乗法で補正された直線）より実験式（一次関数）として次式（４ｂ）を得る。
ｙ＝−２．７７７９ｘ−２８５．９ …………………………………（４ｂ）
この式（４ｂ）と前述した式（４ａ）とを対応させる。
方向係数＝ｄｙ／ｄｘ
＝−２．７７７９
＝（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ……………………………（４ｃ）
ここで、図６（Ｂ）の図表より、Ｓ２０Ｃの鋼材のヤング率とポアッソン比はそれぞれ下記値をとる。
Ｅ＝２０６７８３ＭＰａ、 ν＝０．２８８５ ……………………（４ｄ）
この式（４ｃ）（４ｄ）と式（４）より、次式α_Ｔ２を得る。
α＝α_Ｔ２＝−ν−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ
＝−０．２８８５−２０６７８３×（−２．７７７９）×１０^−６
＝０．２８６ …………………………………………………（４ｅ）

即ち、印加荷重σ_ｎ の印加方向に直交する方向に超音波を伝搬させて音の伝搬時間を測定した場合、Ｓ２０Ｃの鋼材の応力依存係数αは「０．２８６」となる。
ここで、前述した式（５）についても測定値に基づいて同様の演算を試みた。その結果、応力依存係数αの値として上式（４ｅ）の場合とほぼ同等の値を得ることができた。

このため、内部応力σの算定に際し、この算定された応力依存係数αとＳ２０Ｃの鋼材の弾性係数とを、前述した内部応力σの式（１４）に代入すると、内部応力σの算定式は、
σ＝σ_Ｔ２
＝−｛Ｅ／（ν＋α_Ｔ２）｝・（Δｔ／ｔ_０ ）
＝−｛２０６７８３／（０．２８８５＋０．２８６）｝・（Δｔ／ｔ_０）×１０^−６
＝−３５９９３５．５９（Δｔ／ｔ_０ ）×１０^−６ ………………（４ｆ）
となる。

この式（４ｆ）で明らかのように、被測定物の応力依存係数αが事前に特定されていると内部応力σの演算式が極めて単純化される。このため、当該未知の荷重が印加された被測定物についてその時の音伝搬時間ｔと、当該被測定物について負荷荷重ゼロの時の音伝搬時間ｔ_０ とを測定することにより、当該被測定物についてその内部応力を瞬時にしかも定量的に算出することが可能となる。

４−２） ．応力依存係数α_Ｔ１（図８（Ｂ）参照）の場合
図８（Ｂ）で明らかのように、測定値の全体的な傾向（最小二乗法で補正された直線）より実験式（一次関数）として次式（９ｂ）を得る。
ｙ＝１１．５２５ｘ−２４１．４５ ………………………………（９ｂ）
この式（９ｂ）と前述した式（９ａ）とを対応させる。
方向係数＝ｄｙ／ｄｘ
＝１１．５２５
＝（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ …………………………（９ｃ）
ここで、図６（Ｂ）の図表より、Ｓ２０Ｃの鋼材のヤング率として、
Ｅ＝２０６７８３ＭＰａ ……………………………………………（９ｄ）
式（９ｃ）（９ｄ）および式（９）より次式α_Ｔ１を得る。
α＝α_Ｔ１＝１−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ
＝１−２０６７８３×１１．５２５×１０^−６
＝−１．３８３２……………………………………………（９ｅ）
即ち、印加荷重σ_ｎ の印加方向に沿った方向に超音波を伝搬させて音の伝搬時間を測定した場合、Ｓ２０Ｃの鋼材の応力依存係数αは「−１．３８３２」となる。
ここで、前述した式（１０）についても測定値に基づいて同様の演算を試みた。その結果、応力依存係数αの値として上式（９ｅ）の場合とほぼ同等の値を得ることができた。

このため、この算定された応力依存係数αとＳ２０Ｃの鋼材の弾性係数とを、前述した内部応力σの式（１８）に代入すると、
σ＝σ_Ｔ１＝｛Ｅ／（１−α_Ｔ１）｝・（Δｔ／ｔ_０ ）
＝｛２０６７８３／（１＋１．３８３２）｝・（Δｔ／ｔ_０）×１０^−６
＝８６７６６．９５１（Δｔ／ｔ_０ ）×１０^−６………（９ｆ）
となる。
この式（９ｆ）で明らかのように、被測定物の応力依存係数αが事前に特定されていると内部応力σの演算式が極めて単純化される。これがため、未知の荷重が印加された被測定物についてその時の音伝搬時間ｔと、当該被測定物について負荷荷重ゼロの時の音伝搬時間ｔ_０ とを測定することにより、被測定物についてその内部応力を瞬時にしかも定量的に算出することが可能となる。

尚、上記実施形態および実施例は、音伝搬時間の測定に際して縦波超音波を使用した場合について例示したが、横波超音波を使用してもよい。横波超音波を使用した場合でも、同等の手順により、これと同等の応力依存係数αおよび内部応力σにかかる算定式を得ることができる。
又、上記実施形態は、負荷荷重σ_ｎ として引張荷重の場合について例示したが、圧縮荷重であってもよい。この圧縮荷重の場合、内部応力σはベクトル量であるためその符号が変わる。

更に、上記発明およびその実施形態等は、音の伝搬時間を測定して応力依存係数αを算定し未知内部応力σを算出する手法について説明したが、音速情報によっても、これと同等の応力依存係数αおよび内部応力σにかかる各算定式を得ることができる。
この場合、応力依存係数αの算定式は、
α＝Ｅ（ΔＶ_ｎ ／Ｖ_０ ）／σ_ｎ ……………………（イ）
となり、
内部応力σの算定式は、
σ＝Ｅ（ΔＶ／Ｖ_０ ）／α ……………………………（ロ）
となる。これにより、固有常数αが判明すると、内部応力σが極く容易に算定することができる。

ここで、ΔＶ_ｎ ＝Ｖ_ｎ −Ｖ_０ 、ΔＶ＝Ｖ−Ｖ_０ を示す。又、記号Ｖ，Ｖ_ｎ は被測定物の音速で、Ｖは未知の荷重が負荷された状態下で測定された音速を、Ｖ_ｎは応力依存係数αを算定するために印加される荷重負荷σ_ｎ に対応して変化する音速を、それぞれ示す。記号Ｅは被測定物のヤング率を示す。
又、この式（イ）（ロ）については、被測定物に印加する（又は印加されている）荷重σ_ｎ の向きに対して同方向又は直交する方向のいずれの方向に超音波を伝搬させても、適用あるものとする。

本発明は、上述したように、未知の荷重が負荷されている箇所の部材の内部応力を迅速に且つ定量的に測定するためのものであり、過度の荷重が印加された場合に事前にこれを検知して当該部材の破壊を有効に回避することを可能としたものである。このため、本発明は、音速又は音の伝搬時間が測定可能な部材であれば、橋梁の連結部材や各種ボルトをはじめ、電車や自動車の車軸、その他の回転体の支軸等、更には飛行体の部品連結部等、あらゆる箇所について、その内部応力を有効に測定することができ、需要範囲は広い。

本発明の最良の一実施形態を示すブロック図である。 図１における音伝搬時間測定手段による被測定物そのものに対する音の伝搬時間の測定状態を示す説明図である。 図１における音伝搬時間測定手段による被測定物と同一の素材について音の伝搬時間を測定する場合（引張荷重を変化させつつ印加して）を示す説明図である。 図１に示す実施形態の装置全体の動作を示すフローチャートである。 応力依存係数αの性質（印加荷重と音伝搬時間変化比との関連性）を示す図で、図５（Ａ）は音伝搬時間の測定時に超音波を横方向に伝搬させた場合を示し、図５（Ｂ）は音伝搬時間の測定時に超音波を縦方向に伝搬させた場合を示す。 図３における音の伝搬時間の測定に使用される被測定物と同一の素材からなる部材（テストピース）の情報を示す図で、図６（Ａ）は当該部材（テストピース）の形状および概略寸法を示す参考図、図６（Ａ）は当該部材（テストピース）の弾性係数を示す図表である。 図６に開示した部材（Ｓ２０Ｃの鋼材で形成した）を図２に開示した手法で測定し算出された音伝搬時間変化比（Δｔ_ｎ／ｔ_０）とこれに対応する印加荷重（引張荷重σ_ｎ）との関係を示す図表である。 図７の測定情報をグラフ化しその傾斜角情報（測定値の平均値）から応力依存係数αを算出することを意図したもので、図８（Ａ）は音伝搬時間の測定時に超音波を横方向に伝搬させた場合を示す図表、図８（Ｂ）は音伝搬時間の測定時に超音波を縦方向に伝搬させた場合を示す図表である。

符号の説明

１ 被測定物
１Ａ 被測定物と同一の素材からなる他の部材（テストピース）
１１ 音伝搬時間測定手段
１１ａ，１１ｂ 超音波送受波器
１１Ａ 超音波用信号送受信回路部
１１Ｂ 音伝搬時間算定回路部
１２ 応力依存係数算定回路部（応力依存係数算定手段）
１３ 主演算部
１４ 情報記憶部
１５ 入力部
１６ 出力部
２１ 荷重印加装置

Claims (10)

1. 被測定物と同一の素材から成るテストピースについて、その所定箇所の対向面相互間に荷重を変化させながら印加して、当該印加荷重σ_ｎ の大きさに応じて変化する音の伝搬時間ｔ_ｎ ，及び荷重ゼロの時（無負荷時）の音の伝搬時間ｔ _０ を、超音波を用いて且つ当該超音波を前記印加荷重σ_ｎ の印加方向と同方向に伝搬させて測定する第１の工程と、
   この第１の工程で得られる音の伝搬時間ｔ_ｎ とｔ _０ とに基づいて、当該伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）を算定すると共に、これらの情報および前記音伝搬時間情報ｔ_ｎ に対応した印加荷重情報σ_ｎ 並びに前記被測定物を形成する素材のヤング率Ｅを含む固有の弾性定数等を含んで予め特定された応力依存係数αに関する下式、
   α＝ｆ（Δｔ_ｎ ，σ_ｎ ）＝｛１−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝
   により当該被測定物の素材固有の応力依存係数αを算定する第２の工程と、
   この第１乃至第２の各工程に相前後して、装置若しくは建造物等に装備され未知の荷重が印加された状態の前記被測定物について前記所定箇所の対向面相互間を伝搬する音の伝搬時間ｔを測定する第３の工程と、
   この第１乃至第３の各工程で得られる諸情報に基づいて、前記第３工程における被測定物の実際の内部応力σを、前記応力依存係数αを含んで予め特定された内部応力σに関する次式、
   σ＝ｆ（Δｔ，α）＝Ｅ（Δｔ／ｔ_０ ）／（１−α）
   （ここで、Δｔ＝ｔ−ｔ_０ ）
   、を演算し、実際に装備された状態下にある前記被測定物の現実の内部応力σを特定する第４の工程とを備え、
   前記第１の工程では、前記テストピースにおける前記対向面相互間の距離Ｌ _ｒｐ は前記第３の工程における前記被測定物の前記対向面相互間の距離Ｌと同一距離であることを前提として、当該対向面相互間の距離Ｌ _ｒｐ が設定されており、
   前記第２の工程では、更に、前記応力依存係数αの算定時に成される伝搬時間の変化分Δｔ_ｎ にかかる伝搬時間変化比「Δｔ_ｎ ／ｔ_０ 」とこれに対応する応力依存係数「α」との対応関係を順次演算処理しそのデータを記憶する応力依存係数特定工程を同時に実行する構成としたことを特徴とする応力測定方法。
2. 前記請求項１に記載の応力測定方法において、
   前記第２の工程における前記応力依存係数αの算定式、α＝ｆ（Δｔ_ｎ ，σ_ｎ ）を、前記α＝｛１−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝に代えて、下式、
   

   

   としたことを特徴とする応力測定方法。
3. 被測定物と同一の素材から成るテストピースについて、その所定箇所の対向面相互間に荷重を変化させながら印加して、当該印加荷重σ_ｎ の大きさに応じて変化する音の伝搬時間ｔ_ｎ ，及び荷重ゼロの時（無負荷時）の音の伝搬時間ｔ _０ を、超音波を用いて且つ当該超音波を前記印加荷重σ_ｎ の印加方向に直交する方向に伝搬させて測定する第１の工程と、
   この第１の工程で得られる音の伝搬時間ｔ_ｎ とｔ _０ とに基づいて、当該伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）を算定すると共に、これらの情報および前記音伝搬時間情報ｔ_ｎ に対応した印加荷重情報σ_ｎ 並びに前記被測定物を形成する素材のヤング率Ｅとポアソン比νとを含む固有の弾性定数等を含んで予め特定された応力依存係数αに関する下式、
   α＝ｆ（Δｔ_ｎ ，σ_ｎ ）＝｛−ν−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝
   を用いて当該被測定物の素材固有の応力依存係数αを算定する第２の工程と、
   この第１乃至第２の各工程に相前後して、装置若しくは建造物等に装備され未知の荷重が印加された状態の前記被測定物について前記所定箇所の対向面相互間を伝搬する音の伝搬時間ｔを測定する第３の工程と、
   この第１乃至第３の各工程で得られる諸情報に基づいて、前記第３の工程における被測定物の実際の内部応力σを、前記応力依存係数αを含んで予め特定された内部応力σに関する次式、
   σ＝ｆ（Δｔ，α）＝−Ｅ（Δｔ／ｔ_０ ）／（ν＋α）、
   （ここで、Δｔ＝ｔ−ｔ_０ ）
   を演算し、実際に装備された状態下にある前記被測定物の現実の内部応力σを特定する第４の工程とを備え、
   前記第１の工程では、前記テストピースにおける前記対向面相互間の距離Ｌ _ｒｐ は前記第３の工程における前記被測定物の前記対向面相互間の距離Ｌと同一距離であることを前提として、当該対向面相互間の距離Ｌ _ｒｐ が設定されており、
   前記第２の工程では、更に、前記応力依存係数αの算定時に成される伝搬時間の変化分Δｔ_ｎ にかかる伝搬時間変化比「Δｔ_ｎ ／ｔ_０ 」とこれに対応する応力依存係数「α」との対応関係を順次演算処理しそのデータを記憶する応力依存係数特定工程を同時に実行する構成としたことを特徴とする応力測定方法。
4. 前記請求項３に記載の応力測定方法において、
   前記第２の工程における前記応力依存係数αの算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）を、前記α＝｛−ν−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝に代えて、下式、
   

   

   としたことを特徴とする応力測定方法。
5. 装置若しくは建造物等に装備され未知の荷重が印加された状態の被測定物について所定箇所の対向面相互間を伝搬する音の伝搬時間ｔを超音波を用いて測定すると共に、前記被測定物と同一の素材から成るテストピースについて前記被測定物の前記所定箇所の対向面相互間に対応した距離について，荷重を変化させながら印加して当該印加荷重σ_ｎ の大きさに応じて変化する音の伝搬時間ｔ_ｎ 及び荷重ゼロの時（無負荷時）の音の伝搬時間ｔ _０ を、超音波を用いて且つ当該超音波を前記印加荷重σ_ｎ の印加方向と同方向に伝搬させて測定する音伝搬時間測定手段を有し、
   この音伝搬時間測定手段で得られる音の伝搬時間ｔ_ｎ とｔ _０ とにより、当該伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）を算定すると共に、これらの情報および前記音伝搬時間情報ｔ_ｎ に対応した印加荷重情報σ_ｎ 並びに前記被測定物を形成する素材のヤング率Ｅを含む固有の弾性定数等を含んで予め特定された応力依存係数αに関する下式、
   α＝ｆ（Δｔ_ｎ ，σ_ｎ ）＝｛１−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝
   によって、当該被測定物の素材固有の応力依存係数αを算定する応力依存係数算定手段を設け、
   この応力依存係数算定手段で算定される応力依存係数αおよび前記音伝搬時間測定手段で測定される複数の音伝搬時間情報等に基づいて、実際に装置若しくは建造物等に装備され未知の荷重が印加された状態の前記被測定物の未知の内部応力σを、前記応力依存係数αを含んで予め特定された次式、
   σ＝ｆ（Δｔ，α）＝Ｅ（Δｔ／ｔ_０ ）／（１−α）、
   （ここで、Δｔ＝ｔ−ｔ_０ ）
   を演算することにより算定する主演算部を備え、
   前記音伝搬時間測定手段で測定の対象とする前記テストピースにおける前記対向面相互間の距離Ｌ _ｒｐ は前記被測定物の前記対向面相互間の距離Ｌと同一距離であることを前提として、当該対向面相互間の距離Ｌ _ｒｐ が設定されており、
   前記応力依存係数算定手段が、前記応力依存係数αの算定時に成される伝搬時間の変化分Δｔ_ｎ にかかる伝搬時間変化比「Δｔ_ｎ ／ｔ_０ 」とこれに対応する応力依存係数「α」との対応関係を順次演算処理しそのデータを記憶する応力依存係数特定機能を備えていることを特徴とする応力測定装置。
6. 前記請求項５に記載の応力測定装置において、
   前記応力依存係数算定手段は、時間情報抽出機能と、伝搬時間変化分算定機能と、応力依存係数算定機能とを有すると共に、
   前記時間情報抽出機能が、前記音伝搬時間測定手段で得られる印加荷重σ_ｎ を含む前記複数の情報から音伝搬時間情報ｔ_ｎ ，ｔ_０ を抽出するように機能し、
   前記伝搬時間変化分算定機能が、前記時間情報抽出機能で抽出された音伝搬時間情報ｔ_ｎ ，ｔ_０ に基づいて前記印加荷重σ_ｎ に対応した音伝搬時間の変化分Δｔ_ｎ を算定するように機能し、
   前記応力依存係数算定機能が、前記音伝搬時間にかかる情報および前記印加荷重σ_ｎ にかかる情報等に基づいて応力依存係数αにかかる前記算定式
   α＝ｆ（Δｔ_ｎ ，σ_ｎ ）＝｛１−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝
   を演算し前記被測定物を形成する素材固有の応力依存係数αを算定するように機能する構成としたことを特徴とする応力測定装置。
7. 前記請求項５又は６に記載の応力測定装置において、
   前記応力依存係数算定手段では、応力依存係数αを算定する算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）として、前記α＝｛１−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝に代えて、次式、
   

   

   としたことを特徴とする応力測定装置。
8. 装置若しくは建造物等に装備され未知の荷重が印加された状態の被測定物について所定箇所の対向面相互間を伝搬する音の伝搬時間ｔを超音波を用いて測定すると共に、前記被測定物と同一の素材から成るテストピースについて前記被測定物の前記所定箇所の対向面相互間に対応した距離について，荷重を変化させながら印加して当該印加荷重σ_ｎ の大きさに応じて変化する音の伝搬時間ｔ_ｎ 及び荷重ゼロの時（無負荷時）の音の伝搬時間ｔ _０ を、超音波を用いて且つ当該超音波を前記印加荷重σ_ｎ の印加方向と直交する方向に伝搬させて測定する音伝搬時間測定手段を有し、
   この音伝搬時間測定手段で得られる音の伝搬時間ｔ_ｎ とこれに対応した距離のテストピース（被測定物と同一の素材）を用いて測定され換算された無負荷時の音の伝搬時間ｔ_０ とにより当該伝搬時間の変化分Δｔ_ｎ （＝ｔ_ｎ −ｔ_０ ）を算定すると共に、これらの情報および前記音伝搬時間情報ｔ_ｎ に対応した印加荷重情報σ_ｎ 並びに前記被測定物を形成する素材のヤング率Ｅとポアソン比νとを含む固有の弾性定数等を含んで予め特定された応力依存係数αに関する下式、
   α＝ｆ（Δｔ_ｎ ，σ_ｎ ）＝｛−ν−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝
   を演算することにより当該被測定物の素材固有の応力依存係数αを算定する応力依存係数算定手段を設け、
   この応力依存係数算定手段で算定される応力依存係数αおよび前記音伝搬時間測定手段で測定される複数の音伝搬時間情報等に基づいて、実際に装置若しくは建造物等に装備され未知の荷重が印加された状態の前記被測定物の未知の内部応力σを、前記応力依存係数αを含んで予め特定された次式、
   σ＝ｆ（Δｔ，α）＝−Ｅ（Δｔ／ｔ_０ ）／（ν＋α）、（ここで、Δｔ＝ｔ−ｔ_０ ）
   を演算することにより算定する主演算部を備え、
   前記音伝搬時間測定手段で測定の対象とする前記テストピースにおける前記対向面相互間の距離Ｌ _ｒｐ は前記被測定物の前記対向面相互間の距離Ｌと同一距離であることを前提として、当該対向面相互間の距離Ｌ _ｒｐ が設定されており、
   前記応力依存係数算定手段が、前記応力依存係数αの算定時に成される伝搬時間の変化分Δｔ_ｎ にかかる伝搬時間変化比「Δｔ_ｎ ／ｔ_０ 」とこれに対応する応力依存係数「α」との対応関係を順次演算処理しそのデータを記憶する応力依存係数特定機能を備えていることを特徴とする応力測定装置。
9. 前記請求項８に記載の応力測定装置において、
   前記応力依存係数算定手段が、時間情報抽出機能と、伝搬時間変化分算定機能と、応力依存係数算定機能とを有すると共に、
   前記時間情報抽出機能が、前記音伝搬時間測定手段で得られる印加荷重σ_ｎ を含む複数の情報から音伝搬時間情報ｔ_ｎ ，ｔ_０ を抽出するように機能し、
   前記伝搬時間変化分算定機能が、前記時間情報抽出機能が稼働して抽出された音伝搬時間情報ｔ_ｎ ，ｔ_０ に基づいて前記印加荷重σ_ｎ に対応した音伝搬時間の変化分Δｔ_ｎ を算定するように機能し、
   前記応力依存係数算定機能が、前記音伝搬時間にかかる情報および前記印加荷重σ_ｎ にかかる情報に基づいて応力依存係数αにかかる前記算定式
   α＝ｆ（Δｔ_ｎ ，σ_ｎ ）＝｛−ν−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝
   を演算し前記被測定物を形成する素材固有の応力依存係数αを算定するように機能する構成としたことを特徴とする応力測定装置。
10. 前記請求項８又は９に記載の応力測定装置において、
    前記応力依存係数算定手段では、応力依存係数αを算定する算定式α＝ｆ（Δｔ_ｎ ，σ_ｎ ）として、前記α＝｛−ν−Ｅ（Δｔ_ｎ ／ｔ_０ ）／σ_ｎ ｝に代えて、次式、
    

    

    としたことを特徴とする応力測定装置。

Images