JP4768927B2 - 打音検査物の非破壊検査方法、及び品質管理方法 - Google Patents
打音検査物の非破壊検査方法、及び品質管理方法 Download PDFInfo
- Publication number
- JP4768927B2 JP4768927B2 JP2001141731A JP2001141731A JP4768927B2 JP 4768927 B2 JP4768927 B2 JP 4768927B2 JP 2001141731 A JP2001141731 A JP 2001141731A JP 2001141731 A JP2001141731 A JP 2001141731A JP 4768927 B2 JP4768927 B2 JP 4768927B2
- Authority
- JP
- Japan
- Prior art keywords
- dynamic
- vibration
- elastic modulus
- test object
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、打音検査物（例えば、製鉄、製鋼において使用される転炉、取鍋等の工業窯炉および樋など、またアルミニューム、ダクタイル鋳鉄など溶融金属を扱う産業での取鍋等の工業窯炉および樋などの耐火物、特に不定形耐火物、或いは、土木、建築構造物の建設に用いられるシールドセグメント、ヒューム管、コンクリート杭、消波ブロックなどのプレキャストコンクリート製品の如きコンクリート部材）の健全度を検査、管理する打音検査物の非破壊検査方法、その非破壊検査装置及び品質管理方法に関する。
【０００２】
【従来の技術】
打音検査物の非破壊検査に関し、例えば、コンクリート表面を打撃してコンクリート自体を伝搬する打撃音を採取することでコンクリート構造物内部の亀裂、空隙等の反映を健全度として検査する方法が特開平５−３２２８６１号公報に記述されている。
【０００３】
これによると、打撃音の受信方法を工夫することで内部空隙のある状況を受信音の周波数スペクトル解析によって周波数の振幅の大小で健全度を把握できるものとされている。空隙の有無を判定できる限界例は２０ｃｍの幅の空間を実施例で示している。また、この方法は健全度の判定方法として振幅の大小が例示的に示されているのみであり、健全度がどの程度かという指数化については言及されていない。更に、コンクリート構造物はベースコンクリート上に打設されたとの記述があるように、地上に一体固定された構造体を念頭に置いた非破壊検査方法であった。以上のように、構造体が固定されて、空隙の大きさもセンチメートルオーダーで壁状であり、健全性は指数化されていない。
【０００４】
一方、耐火物、特に不定形耐火物は、耐火煉瓦に比べ使用現場での施工作業性が良好であることから、多くの場合、現場での流し込みなどで現場一体型施工が実施されてきたが、近年、品質向上の観点から、現場一体型施工ではなく、分割された成形体を事前に製造し、それを組立て、施工する方法が取り入れられて来ている。組立部品としての不定形耐火物成形体の重量は１０００ｋｇ程度となる場合がある。
【０００５】
現場一体型施工を実施していた従来から、不定形耐火物は成形時にラミネーション、空隙、素材の偏析などによる欠陥が発生しやすく、その解決策の一つとして組立施工が導入され、前記問題点は減少の傾向にある中で、組立部品の製造技術の向上、出荷管理による現場使用時の補修方法への反映などの観点から、組立部品としての不定形耐火物成形体個々の健全度を指数化して把握する方法が望まれていた。
【０００６】
不定形耐火物の欠陥を把握する方法は特開平７−２１８４７９号公報に不定形耐火物施工体の非破壊試験方法と装置が示されている。これによると、施工体に超音波を送信し、施工体内部にある欠陥による反射波を受信することが記述されている。
【０００７】
しかしながら、この方法では不定形耐火物成形体としての健全度を把握するには計測に時間が掛かり健全度を指数化する方法としては問題があった。前記特開平５−３２２８６１号公報に記述されている技術は、前述の如く、健全度を指数化する方法が不明確で、把握できる健全性も空隙の大きさがセンチメートルオーダーとみられ、不定形耐火物成形体に適用するには不適当と判断されること、また、固定されたことを前提としたコンクリート構造体での健全度把握方法であり、固定されていないと言う自由度をもつ不定形耐火物成形体へ適用する技術としては問題があった。
【０００８】
【発明が解決しようとする課題】
そこで本発明の目的は、打音検査物（例えば、耐火物、特に不定形耐火物、コンクリート部材、その他の非固定物）の非破壊検査において、センチメートルオーダー以下の空隙、亀裂まで短時間に検出可能な検査方法、その装置を提供し、更に不定形耐火物成形体の如き非固定打音検査物の健全度を指数化することであり、これにより組立部品等としての製造・出荷の品質を管理する手段を提供することにある。
【０００９】
【課題を解決するための手段】
本発明者らは、これら従来の問題を解決するため、不定形耐火物やコンクリート部材が固定されていないことに着想を得て工夫された耐火物、特に不定形耐火物成形体やプレキャストコンクリート製品の如きコンクリート部材、その他の非固定打音検査物の非破壊検査方法とその装置を完成した。前記目的を達成するために、本発明は次のように構成する。
【００１０】
１．打音検査物表面を打撃した際の振幅が最大となる振動数を測定し、当該振動数と打音検査物の理想固有振動数と比較し、打音検査物の健全度を指数化する非破壊検査方法であって、
前記打音検査物の振動数の測定が、前記打音検査物を有限要素に分割し、有限要素法による固有値解析を実施することにより打音検査物の振動形式を算出し、当該振動形式が得られるように、打音検査物を支持する場所、打撃する場所及び受信する場所を特定して行われる構成であり、
前記理想固有振動数が、前記有限要素法による固有値解析において、打音検査物の理想状態での動弾性率、動ポアソン比及び密度を用いて固有値解析し振動数を算出した当該振動数を理想固有振動数とする構成であり、
前記理想状態での動弾性率及び動ポアソン比を求める手段が、仮の動ポアソン比から仮の動弾性率を求め、前記打音検査物の形状から動剪断弾性率を求め、当該仮の動弾性率と当該動剪断弾性率から、新たな仮の動ポアソン比を求め、当該新たな仮の動ポアソン比から新たな仮の動弾性率を求めて、これを第１回目の動弾性率とし、これを複数回繰り返して、前の回に求めた動弾性率との差が認められなかったときの動弾性率及び動ポアソン比を理想状態での動弾性率及び動ポアソン比とする構成であり、
打音検査物の振動数の受信において、受信振動のうち、打撃による打音検査物の振動開始から１〜１００ミリ秒の受信振動を解析の対象から除くこと、
を特徴とする打音検査物の非破壊検査方法。
【００１６】
３．打音検査物の品質を管理する品質管理方法において、上記１又は２に記載の方法により得られた打音検査物の前記健全度を表す指数を用いて、所定値を外れる場合に異常と判定することを特徴とする打音検査物の品質管理方法。
【００１７】
４．打音検査物が、耐火物又はコンクリート部材の如き非固定物であることを特徴とする上記１又は２に記載の打音検査物の非破壊検査方法。
【００１９】
５．打音検査物が、耐火物又はコンクリート部材の如き非固定物であることを特徴とする上記４に記載の打音検査物の品質管理方法。
【００２０】
本発明の参考発明は、下記の通りである。
６．打音検査物表面を打撃した際の振幅が最大となる振動数を測定し、当該振動数、打音検査物の密度及び理想形状常数から打音検査物の動弾性率を算出し、当該動弾性率と打音検査物の理想状態での動弾性率と比較し、打音検査物の健全度を指数化すること。
７．前記打音検査物表面を打撃する際に、打音検査物を支持する場所、打撃する場所、及び打音検査物の振動を受信する場所を特定すること。
８．打音検査物の健全度を指数化する非破壊検査装置において、打音検査物を支持する支持体、打音検査物表面を打撃する打撃体、打音検査物の振動を受信する１又は２以上の受信体及び前記受信体が受信した受信振動を解析する解析装置を含むこと。
【００２１】
９．前記打撃に使用する打撃体の振動数が２０Ｈｚ〜２０ＫＨｚであること。
【００２２】
１０．前記打撃による打音検査物の振動が２０Ｈｚ〜２０ＫＨｚであること。
【００２３】
１１．前記理想形状常数が打音検査物の理想固有振動数、動弾性率及び密度から算出されること。
【００２４】
１２．前記打音検査物の理想状態での動弾性率、動ポアソン比及び密度を、前記打音検査物と同一材質で製造された小形体の撓み振動数、捻り振動数、寸法及び重量の測定値に基づいて算出し、当該動弾性率、動ポアソン比及び密度を前記理想状態での動弾性率、動ポアソン比及び密度とすること。
【００２５】
１３．前記小形体の体積が０．０３ｍ3以下であること。
【００２９】
１４．前記複数の受信体は、打音検査物の振動が同位相、または１８０度位相の異なると特定される支持体上の該打音検査物に設置されていること。
【００３０】
本発明において、振動形式とは以下のように定義する。
【００３１】
均一な物質で与えられた形状において衝撃などの極短時間の外力により、振動する、所謂固有振動（非減衰自由振動）において、振動形式とは、その多数ある固有振動を形状の変化（所謂、振動の節、腹の位置などの変位（移動）の大小）する様式に着目し分類したものである。その振動形式を振動モード等と称することがある。
【００３２】
尚、或る形状の振動形式（モード）を求める一般的な方法は、例えば、その形状を有限要素に分割し、有限要素法による固有値解析を実施することにより得られる。この際、均一な物質として各要素に対して、密度、動弾性率、動ポアソン比の物性値を一定値として入れるが、振動形式はこの物性値に依らず、形状により定まる。
【００３３】
【発明の実施形態】
次に、本発明の実施態様について図を参照して詳細に説明する。併せて、本発明の参考発明についても説明する。
【００３４】
先ず、本発明の打音検査物が耐火物である場合について、主に説明する。
図１は本発明に係る不定形耐火物成形体３の非破壊検査装置１の構成を示す一実施例である。図１によると、不定形耐火物成形体３は支持体２により地面から隔てられ、打撃体４による打撃により不定形耐火物成形体３を振動させる。この打撃による振動は特定された振動形式で固有振動を大きく振幅させるよう不定形耐火物成形体３の支持体２の場所を特定し、また打撃する場所を特定する。更に、受信体５により受信する場所も固有振動を大きく振幅させた場所におく。このようにして得られた受信信号を解析装置６で解析することで不定形耐火物成形体３の健全度を指数化できる。
【００３５】
不定形耐火物成形体３の製造方法の例としては水、ピッチ等の流動材と不定形耐火原料を混合し、流し込み・振動などの方法で型へ入れ成形し、養生及び乾燥して、その後の型から取り出したものである。この成形体の重量は１０００ｋｇ程度となる場合がある。
【００３６】
図２は本発明に係る不定形耐火物成形体３に関し別の非破壊検査装置１の構成を示す一実施例で、受信体５が複数ある場合の例である。図２では受信体５が２個の場合（５ａ，５ｂ）を示している。複数の受信体５ａ，５ｂを支持体２上の不定形耐火物成形体３に設置するに際して、特定された振動形式で固有振動が同位相、または１８０度位相の異なると特定される場所へ設置する。このようにして得られた受信信号を解析装置６で解析することで不定形耐火物成形体３の健全度を指数化する際に、測定される振動数の確認ができる。
【００３７】
次に、実施例として溶鋼鍋の敷部を不定形耐火物成形体で分割した場合の例を示す。図３は溶鋼鍋の敷部を組み立てる不定形耐火物成形体の一つである。（ａ）は平面図、（ｂ）は側面図である。この容積は０．３３ｍ3であり、嵩密度３．０ｔ／ｍ3の材料を使用すると、この不定形耐火物成形体の重量は約１トンになる。この不定形耐火物成形体は水と不定形耐火原料を混合し、流し込み・振動の方法で型へ入れ成形し、室温で３日養生、１５０℃で５日間乾燥して、その後の型から取り出したものである。
【００３８】
耐火物を有限要素に分割し、有限要素法による固有値解析を実施することにより耐火物の振動形式を算出し、当該振動形式が得られるように、耐火物を支持する場所、打撃する場所及び受信する場所を特定する方法を説明する。
【００３９】
図４は図３の形状を有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した場合を示す第一の振動形式の例図である。振動形式（モード）の出力には形状のみが重要な入力データであり、密度、動弾性率、動ポアソン比は均一体と仮定した各々につき或る値（一定値）を入力することで求められる。固有値解析の計算アルゴリズムは逆べき乗法で計算した。形状は解析モデル形状とするため三次元ソリッド要素により要素分割を行った。実施例では要素数２４２４、節点数３１７６で実施した。要素・節点の相対的変位（移動）は色分けでき、変位（移動）量を見極める判断ができる。
【００４０】
図４（ａ）は平面図、（ｂ）は立面図、（ｃ）は側面図であり、特に図４（ｂ）により不定形耐火物成形体が中央を対象点に左右の両端部が捻れている状況が分かる。従って、中央部分底部は変位（移動）が小さい。また、左右の端部（角）の変位（移動）が大きいことが分かる。このように変位（移動）の大小から、有限要素法により算出された振動形式が得られるように、耐火物を支持する場所、打撃する場所、受信体を置く場所を特定すればよい。この結果から図５に示す如く、図４（ａ）〜（ｃ）で特定された振動形式で大きい振幅が得られるように、耐火物を支持する支持体２を置く場所（△印）、打撃体４により打撃する場所（□印）、また、受信体５ａ、５ｂを置く場所（○印）の例を示した。また、この図は受信体を２個置く場合を示している。この場合は、受信体を置く２つの場所の振動が１８０度位相の異なると特定される場所へ受信体を設置する場合を示している。図４の（ａ）、（ｂ）、（ｃ）から両受信体を置く場所の変位（移動）を観察すると、一方の受信体の場所が上昇方向になると、もう一方の受信体の場所が下降方向になることから、受信体の受信信号の位相が１８０度異なると特定される。
【００４１】
図５は例であり、これにより上記の特定する考え方を制限するものではなく、図４（ａ）、（ｂ）、（ｃ）から得られた変位（移動）の小さい、大きいと判断される場所から実施しやすい場所・方法を選べる。具体的に選ぶ場所を図４（ａ）、（ｂ）、（ｃ）に示す。耐火物を支持する支持体２を置く場所は図４でＡ、Ｂの線分で示される範囲で変位（移動）が小さく適しており、この範囲で選べる。また、打撃体４により打撃する場所、受信体５を置く場所は図４でＣ、Ｄ、Ｅ、Ｆで示される範囲で変位（移動）が大きく、各々打撃する場所、受信体５を置く場所として適している。
【００４２】
以上のように有限要素による固有値解析により得られた捻れを生じる振動形式の例を説明したが、耐火物の振動が点対称に捻れを生じる振動形式があることは、有限要素法による解析を行わなくとも幾何学的な対象性から判断できる場合もある。
【００４３】
しかし、次の例のように、より高次の振動形式で振動を大きく振幅させるように支持する場所、打撃する場所、また、その振動を大きい振幅で受信するように受信体を置く場所が特定される方法として有限要素法による固有値解析が好ましい。
【００４４】
図６はこの形状を有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した場合を示す第二の振動形式の例図である。振動形式を求める計算方法は図４、図５の場合と同様である。（ａ）は平面図、（ｂ）は立面図、（ｃ）は側面図である。有限要素法による解析によれば、耐火物内周側の端部であるＲの付いた角（底部）Ａ、Ｂを一つの固定端とし、外周部の端部からそれぞれ１／４の場所Ｃ、Ｄに変位（移動）の最小のところがあることがわかる。これは側面図（ｃ）のＲ付き角（底部）Ａ、Ｂから外周部に水平に線を引くと外周部の曲線をほぼ２分割することからも分かる。また、この点を結ぶ線上（Ａ−Ｄ、Ｂ−Ｃ）の変位（移動）も小さい。一方、中央部分と外周部の端部（角）の変位（移動）が大きいことが分かる。このように変位（移動）の大小から、算出された振動形式が得られるように、耐火物を支持する場所、打撃する場所、また、受信体を置く場所が特定される。この結果から図７に示す如く、図６（ａ）〜（ｃ）で特定された振動形式が得られるように、耐火物を支持する支持体２を置く場所（△印）、打撃体４により打撃する場所（□印）、また、受信体５ａ，５ｂを置く場所（○印）を示した。また、この図は受信体を２個置く場合を示している。この場合は、受信体を置く２つの場所の振動が同位相となると特定される場所へ受信体を設置する場合を示している。図６の（ａ）、（ｂ）、（ｃ）から両受信体を置く場所の変位（移動）を観察すると、一方の受信体の場所が上昇方向になると、もう一方の受信体の場所も上昇方向になることから、受信体の受信信号の位相が同位相になると特定される。図７は例であり、これにより上記の特定する考え方を制限するものではなく、図６（ａ）、（ｂ）、（ｃ）から得られた変位（移動）の小さい、大きいと判断される場所から実施しやすい場所・方法を選べる。具体的に選ぶ場所を図６（ａ）、（ｂ）、（ｃ）に示す。耐火物を支持する支持体２を置く場所は図６でＡ−Ｄ、Ｂ−Ｄの線分で示される範囲で変位（移動）が小さく適しており、この範囲で選べる。また、打撃体４により打撃する場所、受信体５を置く場所は図６でＥ、Ｆ、Ｇ、Ｈで示される範囲で変位（移動）が大きく、各々打撃する場所、受信体５を置く場所として適している。
【００４５】
次に、耐火物の理想状態での動弾性率、動ポアソン比、密度を算出する方法について説明する。
【００４６】
前記成形体は０．３３ｍ3と大きく、嵩密度３．０ｔ／ｍ3の材料を使用すると、この不定形耐火物成形体の重量は約１トンになる。このように大形状では成形上のラミネーション、空隙、素材の偏析などによる欠陥が生じる可能性が回避出来ないが、成形容積を小さくすることで、例えば１／１０倍の容積以下（０．０３ｍ3）に小さくすることで、成形上の前記欠陥を大きく回避できることに着目した。大きさ以外は組立部品とする不定形耐火物成形体の成形方法と同一とすることが好ましく、大きさの違いから流し込み成形の振動時間が異なることや、また、養生、乾燥の時間が異なることがあることなど、成形体の品質を変えない範囲で選定することを制限するものではない。本発明はこの方法により成形されたものを理想状態での成形体、理想状態成形体とした。
【００４７】
この理想状態成形体を打撃し、打撃による成形体の振動数を測定し、理想状態での動弾性率、動ポアソン比を算出した。理想状態成形体の一例の形状として、縦１５ｃｍ×横１５ｃｍ×長さ５３ｃｍの０．０１２ｍ3矩形体を作製した。この形状を複数個作製することで、理想状態でのバラツキを明確にすることができる。ここでは、３個（不定形耐火物Ａ）、あるいは１０個（不定形耐火物Ｂ）の成形体を作製し、重量を測定して寸法から密度を算出し、また撓み振動数、捻り振動数を測定し、それにより動弾性率(撓み)、動剪断弾性率、動ポアソン比を算出した。その後、平均値の動弾性率（撓み）、動ポアソン比、密度を求めた。この動弾性率（撓み）、動ポアソン比、密度が不定形耐火物成形体の理想状態での動弾性率、動ポアソン比、密度である。尚、ＪＩＳ Ａ １１２７−１９９３には矩形体など簡単な形状での撓み振動数、捻り振動数の測定方法及び動弾性率(撓み)、動剪断弾性率の計算方法が示されており、今回の成形体を打撃し、打撃による成形体の振動を測定、弾性率を計算する場合においてそれに倣っており、以下に要点を説明する。
【００４８】
表１に不定形耐火物Ａに関し縦１５ｃｍ×横１５ｃｍ×長さ５３ｃｍの矩形体の撓み振動数、捻り振動数を測定した結果と計測値の密度を示す。この結果から動弾性率（撓み）Ｅd（Ｎ／ｍｍ2）、動剪断弾性率Ｇd（Ｎ／ｍｍ2）と動ポアソン比νd［−］を求める方法を説明する。
【００４９】
動弾性率（撓み）Ｅdは、長さＬ、幅さｂ、高さｔ、（本例では長さＬ＝５３０ｍｍ、幅さｂ＝１５０ｍｍ、高さｔ＝１５０ｍｍ）、質量ｍ（ｋｇ）、密度［ｋｇ／ｍ3］で、（密度は計測値、振動数ｆ（１／ｓｅｃ）は測定値）を用いて次式で求められる。
【００５０】
尚、Ｔは修正係数である。本実施例では矩形体であるので、t／（３．４６４Ｌ）と動ポアソン比νdの関数で与えられる。動ポアソン比νdを１／６とした場合の修正係数Tは表２で示されている。例えば、ｔ＝１５０［ｍｍ］、Ｌ＝５３０［ｍｍ］、νd＝１／６の場合、Ｔ＝１．５である。
【００５１】
更に、そのＴを用いて、動ポアソン比がνdである場合の修正係数Ｔ'は次式で求められる。
T'＝T(1+(0.26νd +3.22 (νd)2)(t/(3.464L))/(1+0.1328(t/(3.464L))
【００５２】
動弾性率（撓み）Ｅdを求める実際の手順は、先ず、動ポアソン比νdを１／６と仮定して、仮の動弾性率（撓み）Ｅdを求め、後に述べる手順で求める動剪断弾性率Ｇdから得られる動ポアソン比νdを用いて再度計算し直す。これが第１回目の値となり、これを繰り返すことで精度のある動弾性率（撓み）Ｅdが求められる。
【００５３】
本例の場合は動ポアソン比を求める試行を６回繰り返した。１回目の値と２回目は有意な差が認められ、２回目から６回目までは差が認められなかった。その結果を表３に示す。
【００５４】
動剪断弾性率Ｇd（Ｎ／ｍｍ2）は次式で与えられる。
【００５５】
ここで、Ｒは形状係数と呼ばれ、ｔ／ｂの関数（ｔ≦ｂ）であり、正方形断面（本実施例は縦１５０ｍｍ×横１５０ｍｍ）の角柱の場合は、Ｒ＝１．１８３である。そこで、動弾性率（歪み）Ｅd、動剪断弾性率Ｇdより、動ポアソン比νdは
νd＝｛Ｅd／（２Ｇd）｝−１で求められる。
【００５６】
上記と同様に、不定形耐火物Ｂについても密度、撓み振動数、捻り振動数を測定し、動弾性率（歪み）Ｅd、動剪断弾性率Ｇd、動ポアソン比νdを算出し、更にその平均値を求め、理想状態での動弾性率、動ポアソン比、密度を算出した。
【００５７】
前記の理想状態での動弾性率、動ポアソン比、密度を表４に示す。
【００５８】
尚、前記の成形体を打撃し、打撃による成形体の振動を測定する方法によらず、例えば超音波の透過時間から動弾性率を求める方法もあり、このような他の方法によって決めることを制限しない。但し、この場合は動ポアソン比を経験値などで仮定することになる。また、前記の成形体を打撃し、打撃による成形体の振動を測定し、動弾性率（撓み）を算出する実施例を示したが、縦振動を測定して動弾性率（縦）を算出し動弾性率とすることを制限するものではない。本発明の視点は小形成形体が理想状態との観点であり、その形状からの動弾性率を測定する方法を制限するものではない。
【００５９】
【表１】
【００６０】
【表２】
【００６１】
【表３】
【００６２】
次に、耐火物の理想固有振動数を算出する方法を説明する。
有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した。形状の入力データのほか、表４の理想状態での密度、動弾性率、動ポアソン比を各々要素につき入力することで理想固有振動数が求められる。固有値解析の計算アルゴリズムは逆べき乗法で計算した。形状は図４と同様に三次元ソリッド要素により分割を行い、実施例では要素数２４２４、節点数３１７６で実施した。不定形耐火物Aでは、実際のデータを密度は要素No.１から要素No.２４２４の全要素を２２７０（ｋｇ／ｍ3）、動弾性率は同じく全要素を２９．０Ｅ９（意味は２９．０×１０9）（Ｐａ）、動ポアソン比は同じく全要素を０．２５とした。不定形耐火物Bでは、実際のデータを密度は要素No.１から要素No.２４２４の全要素を３３４０（ｋｇ／ｍ3）、動弾性率は同じく全要素を１４９Ｅ９（意味は１４９×１０9）（Ｐａ）、動ポアソン比は同じく全要素を０．１４とした。
【００６３】
このようにして有限要素法（ＦＥＭ）による固有値解析を実施して得られた理想固有振動数を表５に示す。このように、理想状態での諸物性（密度、動弾性率、動ポアソン比）と図４、図６など振動形式を特定することで理想固有振動数が決まる。
【００６４】
この表５により分かるように図４、図６の理想固有振動数は不定形耐火物の種類Ａ、Ｂにより大きく異なる。一方で、耐火物の種類Ａ、Ｂによる理想固有振動数の比（Ａ／Ｂ）は図４、図６のそれぞれの振動形式でほぼ一定（０．５３、０．５２）であり、また図４、図６の振動形式による振動数の比（（１）／（２））はＡ、Ｂでほぼ一定（０．８８、０．８９）の関係が得られる。
【００６５】
表５の同一の振動形式でＡ、Ｂの理想状態での物性値に基づく理想固有振動数の変化（例えば図６の振動形式で４９４Ｈｚ、９５０Ｈｚ）は、見方を変えると、これは耐火物の振動数を測定した結果に対する解釈を与えている。例えば、耐火物Ｂについて、振動数を測定した結果、理想固有振動数（９５０Ｈｚ）より小さい振動数（８３０Ｈｚ）になる場合、密度、動弾性率、動ポアソン比が反映されて、理想状態でのＢの諸物性からＡのほうへずれたためと解釈される。この考えに基づき、該耐火物の振動数を測定し、該耐火物の理想固有振動数と比較し健全度として、測定振動数を指数化する。
【００６６】
Ａ、Ｂ間の固有振動数の比がほぼ一定（０．５３、０．５２）であるので図４、図６のいずれの振動形式で評価しても、同等の健全度を表すことが分かる。
【００６７】
測定された振動数と理想固有振動数との比で健全度を指数化する方法としては、上述のようであるが、測定された振動数から算出された動弾性率と理想状態での動弾性率との比が健全度の指数として好ましい。
【００６８】
そこで、測定された振動数から動弾性率の算出方法について説明する。
理想状態での動弾性率Ｅ0［Pa］と密度ρ0［ｋｇ／ｍ3］、理想固有振動数ｆ0［１／sec］が次式の関係にある。
（Ｅ0）=(Ｋ)×（ρ0）×（ｆ0）2
【００６９】
ここでＫを理想形状常数［ｍ2］と称する。この理想形状常数は形状、動ポアソン比に関係する常数である。
【００７０】
この関係式により、理想形状常数Ｋを算出し、耐火物の実測された密度ρm、測定された振動数ｆmを入れて、弾性率Ｅmが求められる。
【００７１】
以上の説明をまとめると
【００７２】
不定形材料Ｂの図６における理想形状常数は次式より、
(理想形状常数Ｋ[m2])＝（動弾性率Ｅ0：１４９×１０9）／｛（密度ρ0：３３４０）（理想固有振動数ｆ0：９５０）2｝
Ｋ：４９．４[m2]である。
【００７３】
弾性率Ｅmの一例を計算すると、密度ρm＝３２５０［ｋｇ／ｍ3］、測定された振動数ｆm＝８３２［１／sec］を代入して、
弾性率Ｅmが１１１[ＧＰａ]となる。
【００７４】
健全度の指数として測定された振動数と理想固有振動数との比、動弾性率と理想状態での動弾性率との比ではなく、測定され振動数や算出された動弾性率そのものを使うことを制限するものではない。いずれにしても、理想固有振動数、理想状態での弾性率が想定されることが好ましい。
【００７５】
【表４】
【００７６】
【表５】
【００７７】
次に、耐火物の表面を打撃し、該打撃による耐火物の振動数を測定する範囲、打撃体の振動数について説明する。
【００７８】
不定形耐火物成形体の振動形式と理想固有振動数は有限要素法（ＦＥＭ）により固有値解析を実施して求めることは既に説明した。そこでは理想状態での物性値（密度、動弾性率、動ポアソン比）を適宜変更することで、予測される振動が算出され、表４に示すような代表的な理想状態での物性からは、図４、図６の場合で大型形状の理想固有振動数は表５に示すようにせいぜい１ＫＨｚである。一般的に、大型形状即ち、長尺なものになると周波数が小さくなる。一方、耐火物で製造する代わりに、金属で代替し、例としてアルミニューム金属製で図３の１／４サイズで約１４Ｋｇの場合で実測すると、図４、図６の振動形式の理想固有振動数はせいぜい３ＫＨｚである。
【００７９】
また、振動には倍振動、３倍振動などもあり、打撃による成形体の振動が２０Ｈｚから２０ＫＨｚにあると限定した。
【００８０】
また、この振動を大きく振幅させるためには、打撃に使用する打撃体の加振振動数が２０Ｈｚから２０ＫＨｚにあることが必要であり、ゴムは２０Ｈｚ〜５００Ｈｚ、プラスチックは５００Ｈｚ〜３Ｋｚ、３ＫＨｚ以上はスチールがその振動数の振幅を大きくさせるために適当な材質の例である。
【００８１】
従って、その振動を測定するための受信体は上記２０Ｈｚから２０ＫＨｚの大きく振幅された振動数が受信できる装置で有ればよく、特開平５−３２２８６１号公報に記述されている工夫された受信体であることに制限するものではなく、また通常の振動ピックアップ、加速度計であってもよい。
【００８２】
次にどの振動形式を測定の対象とするかの選択方法について説明する。
理想状態でのデータから得られる有限要素法による固有値解析で得た振動形式の例を図４、図６で説明したが、有限要素法による固有値解析から他にも多くの振動形式がある。振動形式と理想固有振動数は一対一の関係があるので、簡単のために、理想固有振動数で示すと、表４の不定形耐火物Ａについての理想固有振動数は表６に示すようになる。表６に示す以外にも高次の振動数があるが、省略する。表６で示す振動数のうちどれを選択するかとの判断が必要であるが、振動解析の精度を重視するために、比較的小さい振動数の振動を選択することが望ましい。
【００８３】
また、比較的小さい振動形式の中からは図４、図６に示すように、図４は支持する場所が２点、図６では、支持する場所が４点であり、成形体の重量が１トンにも及ぶ場合があり、安全上の配慮からも選択される。
【００８４】
【表６】
【００８５】
このような観点から、以下の実施例では図６で特定される振動を大きく振幅させ、その振動を解析した例を説明する。
【００８６】
図６に、表４の種類Ａにより製造した不定形耐火物成形体を、図７で特定されるように支持及び打撃し、大きく振幅させ、振動を受信した場合の解析結果を示す。 解析装置は市販の高速フーリエ変換装置（ＦＦＴ）を内蔵した解析装置である。
【００８７】
図６（ａ）は、打撃体がゴムの場合、図６（ｂ）は、打撃体がプラスチックの場合、図６（ｃ）は、打撃体が鋼球の場合で、横軸が周波数、縦軸がその振幅を示す。理想固有振動数４９４Ｈｚ付近の測定値は５００Ｈｚで、これは解析装置の分解能による誤差が原因で５００Ｈｚと測定された。理想固有振動数４９４Ｈｚを大きく振幅させる打撃方法としてゴム、プラスチックが適当な材質であることが分かる。何れの場合も理想固有振動数付近に明確な振動数が解析されるが、その他にも１ＫＨｚ、１．４ＫＨｚなどに共通のピークがある。
【００８８】
これら高次の振動数は、不定形耐火物成形体の健全度を指数化する本非破壊検査方法にとって意味がなく、除外することで解析精度が向上したり、測定振動が特定された振動を測定しているか判断する際にノイズになるので除外する。
【００８９】
この高次の振動数を除外する方法としては、打撃による成形体の振動開始から１〜１００ミリ秒の受信振動を解析の対象から除き、その後の受信振動を解析することを発明した。また、請求項４に記載の発明において、成形体の理想固有振動数を含みながら最も小さい測定周波数を選択する機能を有することが好ましい。
【００９０】
前記請求項４に記載の発明においては、ＦＦＴ解析の以前に高次の振動数を測定しないことにより、ＦＦＴ解析によりスペクトル解析の精度が向上し、最も小さい測定周波数を選択する機能に１００Ｈｚ、２００Ｈｚ、５００Ｈｚ、１ＫＨｚ、２ＫＨｚ、５ＫＨｚ、１０ＫＨｚ、２０ＫＨｚから測定周波数を選択する機能を有することで解析装置が容易に使用できる。また、このような目的で特定の周波数以上のフィルターを掛けＦＦＴ解析の以前に高次の振動を測定しない方法もあり、これを制限するものではない。
【００９１】
高次の振動数を除外する方法の実施例として、Ａ材質の不定形耐火物成形体を図７で示すように支持及び打撃し、大きく振幅させ、受信した場合での実施例としての解析結果を図９に示す。この解析に当たり、成形体の振動開始から５０ミリ秒の受信振動を解析の対象から除外した。この実施例は、また２個の受信体を配置した例であり（前記８の発明）、また前記１４の好ましい実施態様として、成形体の振動が同位相になる場所へ受信体を配置した場合の実施例でもある。
【００９２】
図９（ａ）で５００Ｈｚの振動数が大きく振幅していることが確認できた。また、図８の１ＫＨｚ、１．４ＫＨｚの振動数が図９では小さく押さえられていることが分かる。この結果、成形体の振動が同位相と特定される場所へ設置した受信体の受信信号を解析した結果、振動が、図９（ｂ）で示すように同位相を主体にした状態が確認できるまでに向上した。
【００９３】
このように、打撃による成形体の振動開始から５０ミリ秒の受信振動を解析の対象から除き、その後の受信振動を解析することで、受信した振動が測定対象の振動であることを確認することが出来た。実施例では、打撃による成形体の振動開始から５０ミリ秒の受信振動を解析の対象から除いたが、この時間は、理想固有振動数により適宜変動させ、適正な値を選定することができる。だたし、受信振動を解析の対象から除く期間が長期間になると、解析対象の振動が小さくなり、誤差が生じ易いので、上限を１００ミリ秒とした。また、１ミリ秒未満では、外乱の受信振動を解析の対象から除くことが不充分である。従って、本発明において、打撃による成形体の振動開始から１〜１００ミリ秒の受信振動を解析の対象から除き、その後の受信振動を解析する機能を有することとした。
【００９４】
次に、不定形耐火物の成形体表面を打撃し、該打撃による不定形耐火物の成形体の振動を解析することで該不定形耐火物成形体の健全度を指数化することを実施例で説明する。
【００９５】
表４の種類Ｂの不定形耐火物成形体を、図７で特定されるように支持及び打撃し、大きく振幅させ、振動を受信した。種類Ｂの不定形耐火物成形体は５個製造した。理想状態での動弾性率、密度は表４Ｂに示す１４９ＧＰａ、３３４０ｋｇ／ｍ3、理想固有振動数は表５に示すとおり９５０Ｈｚである。この場合の不定形材料Ｂの理想形状常数は(理想形状常数Ｋ[m2])＝（動弾性率Ｅ0：１４９×１０9[Pa]）／｛（密度ρ0：３３４０[ｋｇ／m3]）（理想固有振動数ｆ0：９５０[1／sec]）2｝より、
Ｋ＝４９．４[m2]であり、これを用いて、実測した振動数から既に述べた次式で弾性率を算出した。Ｅm＝Ｋ・ρm・ｆm 2
尚、密度は何れも３２５０[ｋｇ／m3]であった。
【００９６】
表７に測定された振動数と上述の式で算出した動弾性率を示す。
【００９７】
【表７】
【００９８】
表７より、振動数、動弾性率が何れも理想固有振動数、理想状態での動弾性率より小さく、異なった値であった。この値を指数化した結果を表８に示す。理想固有振動数あるいは理想状態での動弾性率を１とした比の値で示されている。
【００９９】
【表８】
【０１００】
表８に示すように耐火物の表面を打撃し、該打撃による耐火物の振動を最大となるように耐火物の振動数を測定する本発明によると、容易に該不定形耐火物成形体の健全度を指数化することができる。
【０１０１】
表８では理想固有振動数、理想状態での動弾性率を１とした比で示したが、製造の技術管理上はいくつか成形体を測定した結果で、最大のものを１とするような比で示すことや、測定された振動数そのものや、算出した弾性率そのものの値を指数にすることを制限するものではない。尚、弾性率に比べ振動数の場合は健全度が圧縮されて表現される傾向であり、理想状態での動弾性率との比を採用することが好ましい。
【０１０２】
本発明の目的である耐火物の非破壊検査により健全度を指数化することは表８に示された。次に、これにより耐火物の品質を管理する方法の一例を示す。
【０１０３】
図１０は、表８の結果を母数５個としてワイブルプロットしたものである。図１０中の（ａ）、（ｂ）は表８中の（ａ）、（ｂ）と対応している。動弾性率に比べ振動数の場合は健全度が圧縮されて表現される傾向であり、理想固有振動数との比よりも理想弾性率との比を採用することが好ましいことが分かる。
【０１０４】
右側上方にあるほど理想に近く、左側下方にあるほど健全度が悪化していることを示す。理想から離れた品質がこのように順位付けられ、製造技術との関連性を容易に把握でき、優れた製造技術の管理が実施できるばかりでなく、使用上の品質を把握して使用中の補修箇所、時期を予測する方法へ活用できるなどの面で有益な作用を奏することができる。
【０１０５】
この他にも組立部品として品質を管理する統計上の方法はあり、それを制限するものではない。
【０１０６】
以上の実施態様においては、不定形耐火物成形体の非破壊検査方法、その非破壊検査装置及びそれによる品質管理方法を実施例に説明しているが、本発明は不定形耐火物成形体に限定されるものではなく、定形耐火物、ノズル煉瓦、炉底用カーボン煉瓦、炭化珪素、窒化珪素などのセラミック焼結体、各種セラミックスなどの固定されていない耐火物、成形体へも適用可能なものである。
【０１０７】
次に、本発明の打音検査物がコンクリート部材である場合について説明する。
尚、図は前記耐火物の説明で用いたものを使用する。
【０１０８】
図１は本発明に係るコンクリート部材３の非破壊検査装置１の構成を示す一実施例である。図１によると、コンクリート部材３は支持体２により地面から隔てられ、打撃体４による打撃によりコンクリート部材３を振動させる。この打撃による振動は特定された振動形式で固有振動を大きく振幅させるようコンクリート部材３の支持体２の場所を特定し、また打撃する場所を特定する。更に、受信体５により受信する場所も固有振動を大きく振幅させた場所におく。このようにして得られた受信信号を解析装置６で解析することでコンクリート部材３の健全度を指数化できる。
【０１０９】
図２は本発明に係るコンクリート部材３に関し別の非破壊検査装置１の構成を示す一実施例で、受信体５が複数ある場合の例である。図２では受信体５が２個の場合（５ａ，５ｂ）を示している。複数の受信体５ａ，５ｂを支持体２上のコンクリート部材３に設置するに際して、特定された振動形式で固有振動が同位相、または１８０度位相の異なると特定される場所へ設置する。このようにして得られた受信信号を解析装置６で解析することでコンクリート部材３の健全度を指数化する際に、測定される振動数の確認ができる。
【０１１０】
次に、実施例として或る形状のコンクリート部材を示す。（ａ）は平面図、（ｂ）は側面図である。この容積は０．３３ｍ3であり、嵩密度２．３ｔ／ｍ3の材料を使用すると、このコンクリート部材の重量は約０．７６トンになる。
【０１１１】
コンクリート部材を有限要素に分割し、有限要素法による固有値解析を実施することによりコンクリート部材の振動形式を算出し、当該振動形式が得られるように、コンクリート部材を支持する場所、打撃する場所及び受信する場所を特定する方法を説明する。
【０１１２】
図４は図３の形状を有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した場合を示す第一の振動形式の例図である。振動形式（モード）の出力には形状のみが重要な入力データであり、密度、動弾性率、動ポアソン比は均一体と仮定した各々につき或る値（一定値）を入力することで求められる。固有値解析の計算アルゴリズムは逆べき乗法で計算した。形状は解析モデル形状とするため三次元ソリッド要素により要素分割を行った。実施例では要素数２４２４、節点数３１７６で実施した。要素・節点の相対的変位（移動）は色分けでき、変位（移動）量を見極める判断ができる。
【０１１３】
図４（ａ）は平面図、（ｂ）は立面図、（ｃ）は側面図であり、特に図４（ｂ）によりコンクリート部材が中央を対象点に左右の両端部が捻れている状況が分かる。従って、中央部分底部は変位（移動）が小さい。また、左右の端部（角）の変位（移動）が大きいことが分かる。このように変位（移動）の大小から、有限要素法により算出された振動形式が得られるように、コンクリート部材を支持する場所、打撃する場所、受信体を置く場所を特定すればよい。
【０１１４】
この結果から図５に示す如く、図４（ａ）〜（ｃ）で特定された振動形式で大きい振幅が得られるように、コンクリート部材を支持する支持体２を置く場所（△印）、打撃体４により打撃する場所（□印）、また、受信体５ａ、５ｂを置く場所（○印）の例を示した。
【０１１５】
また、この図は受信体を２個置く場合を示している。この場合は、受信体を置く２つの場所の振動が１８０度位相の異なると特定される場所へ受信体を設置する場合を示している。図４の（ａ）、（ｂ）、（ｃ）から両受信体を置く場所の変位（移動）を観察すると、一方の受信体の場所が上昇方向になると、もう一方の受信体の場所が下降方向になることから、受信体の受信信号の位相が１８０度異なると特定される。
【０１１６】
図５は例であり、これにより上記の特定する考え方を制限するものではなく、図４（ａ）、（ｂ）、（ｃ）から得られた変位（移動）の小さい、大きいと判断される場所から実施しやすい場所・方法を選べる。具体的に選ぶ場所を図４（ａ）、（ｂ）、（ｃ）に示す。コンクリート部材を支持する支持体２を置く場所は図４でＡ、Ｂの線分で示される範囲で変位（移動）が小さく適しており、この範囲で選べる。また、打撃体４により打撃する場所、受信体５を置く場所は図４でＣ、Ｄ、Ｅ、Ｆで示される範囲で変位（移動）が大きく、各々打撃する場所、受信体５を置く場所として適している。
【０１１７】
以上のように有限要素による固有値解析により得られた捻れを生じる振動形式の例を説明したが、コンクリート部材の振動が点対称に捻れを生じる振動形式があることは、有限要素法による解析を行わなくとも幾何学的な対象性から判断できる場合もある。
【０１１８】
しかし、次の例のように、より高次の振動形式で振動を大きく振幅させるように支持する場所、打撃する場所、また、その振動を大きい振幅で受信するように受信体を置く場所が特定される方法として有限要素法による固有値解析が好ましい。
【０１１９】
図６はこの形状を有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した場合を示す第二の振動形式の例図である。振動形式を求める計算方法は図４、図５の場合と同様である。（ａ）は平面図、（ｂ）は立面図、（ｃ）は側面図である。有限要素法による解析によれば、コンクリート部材内周側の端部であるＲの付いた角（底部）Ａ、Ｂを一つの固定端とし、外周部の端部からそれぞれ１／４の場所Ｃ、Ｄに変位（移動）の最小のところがあることがわかる。これは側面図（ｃ）のＲ付き角（底部）Ａ、Ｂから外周部に水平に線を引くと外周部の曲線をほぼ２分割することからも分かる。また、この点を結ぶ線上（Ａ−Ｄ、Ｂ−Ｃ）の変位（移動）も小さい。一方、中央部分と外周部の端部（角）の変位（移動）が大きいことが分かる。
【０１２０】
このように変位（移動）の大小から、算出された振動形式が得られるように、コンクリート部材を支持する場所、打撃する場所、また、受信体を置く場所が特定される。この結果から図７に示す如く、図６（ａ）〜（ｃ）で特定された振動形式が得られるように、コンクリート部材を支持する支持体２を置く場所（△印）、打撃体４により打撃する場所（□印）、また、受信体５ａ，５ｂを置く場所（○印）を示した。また、この図は受信体を２個置く場合を示している。この場合は、受信体を置く２つの場所の振動が同位相となると特定される場所へ受信体を設置する場合を示している。
【０１２１】
図６の（ａ）、（ｂ）、（ｃ）から両受信体を置く場所の変位（移動）を観察すると、一方の受信体の場所が上昇方向になると、もう一方の受信体の場所も上昇方向になることから、受信体の受信信号の位相が同位相になると特定される。図７は例であり、これにより上記の特定する考え方を制限するものではなく、図６（ａ）、（ｂ）、（ｃ）から得られた変位（移動）の小さい、大きいと判断される場所から実施しやすい場所・方法を選べる。具体的に選ぶ場所を図６（ａ）、（ｂ）、（ｃ）に示す。
【０１２２】
コンクリート部材を支持する支持体２を置く場所は図６でＡ−Ｄ、Ｂ−Ｄの線分で示される範囲で変位（移動）が小さく適しており、この範囲で選べる。また、打撃体４により打撃する場所、受信体５を置く場所は図６でＥ、Ｆ、Ｇ、Ｈで示される範囲で変位（移動）が大きく、各々打撃する場所、受信体５を置く場所として適している。
【０１２３】
次に、コンクリート部材の理想固有振動数を算出する方法を説明する。
有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した。形状の入力データのほか、表９の理想状態での密度、動弾性率、動ポアソン比を各々要素につき入力することで理想固有振動数が求められる。固有値解析の計算アルゴリズムは逆べき乗法で計算した。形状は図４と同様に三次元ソリッド要素により分割を行い、実施例では要素数２４２４、節点数３１７６で実施した。コンクリート部材では、実際のデータを密度は要素Ｎｏ．１から要素Ｎｏ．２４２４の全要素を２２７０（ｋｇ／ｍ3）、動弾性率は同じく全要素を２９．０Ｅ９（意味は２９．０×１０9）（Ｐａ）、動ポアソン比は同じく全要素を０．２５とした。
【０１２４】
このようにして有限要素法（ＦＥＭ）による固有値解析を実施して得られた理想固有振動数を表１０に示す。このように、理想状態での諸物性（密度、動弾性率、動ポアソン比）と図４、図６など振動形式を特定することで理想固有振動数が決まる。
【０１２５】
そこで、測定された振動数から動弾性率の算出方法について説明する。
理想状態での動弾性率Ｅ0［Ｐａ］と密度ρ0［ｋｇ／ｍ3］、理想固有振動数ｆ0［１／ｓｅｃ］が次式の関係にある。
（Ｅ0）＝（Ｋ）×（ρ0）×（ｆ0）2
【０１２６】
ここでＫを理想形状常数［ｍ2］と称する。この理想形状常数は形状、動ポアソン比に関係する常数である。
【０１２７】
この関係式により、理想形状常数Ｋを算出し、コンクリート部材の実測された密度ρm、測定された振動数ｆmを入れて、弾性率Ｅmが求められる。
【０１２８】
以上の説明をまとめると
【０１２９】
コンクリート部材の図６における理想形状常数は次式より、
（理想形状常数Ｋ［ｍ2］）＝（動弾性率Ｅ0：２９×１０9）／｛（密度ρ0：２２７０）（理想固有振動数ｆ0：４９４）2｝
Ｋ：５２．３５［ｍ2］である。
【０１３０】
弾性率Ｅmの一例を計算すると、密度ρm＝２０００［ｋｇ／ｍ3］、測定された振動数ｆm＝３００［１／ｓｅｃ］を代入して、
弾性率Ｅmが９．４２［ＧＰａ］となる。
【０１３１】
健全度の指数として測定された振動数と理想固有振動数との比、動弾性率と理想状態での動弾性率との比ではなく、測定され振動数や算出された動弾性率そのものを使うことを制限するものではない。いずれにしても、理想固有振動数、理想状態での弾性率が想定されることが好ましい。
【０１３２】
次に、コンクリート部材の表面を打撃し、該打撃によるコンクリート部材の振動数を測定する範囲、打撃体の振動数について説明する。
【０１３３】
コンクリート部材の振動形式と理想固有振動数は有限要素法（ＦＥＭ）により固有値解析を実施して求めることは既に説明した。そこでは理想状態での物性値（密度、動弾性率、動ポアソン比）を適宜変更することで、予測される振動が算出され、表９に示すような代表的な理想状態での物性からは、図４、図６の場合で大型形状の理想固有振動数は表１０に示すようにせいぜい１ＫＨｚである。一般的に、大型形状即ち、長尺なものになると周波数が小さくなる。一方、コンクリートで製造する代わりに、金属で代替し、例としてアルミニューム金属製で図３の１／４サイズで約１４Ｋｇの場合で実測すると、図４、図６の振動形式の理想固有振動数はせいぜい３ＫＨｚである。
【０１３４】
また、振動には倍振動、３倍振動などもあり、打撃による成形体の振動が２０Ｈｚから２０ＫＨｚにあると限定した。
【０１３５】
また、この振動を大きく振幅させるためには、打撃に使用する打撃体の加振振動数が２０Ｈｚから２０ＫＨｚにあることが必要であり、ゴムは２０Ｈｚ〜５００Ｈｚ、プラスチックは５００Ｈｚ〜３Ｋｚ、３ＫＨｚ以上はスチールがその振動数の振幅を大きくさせるために適当な材質の例である。
【０１３６】
従って、その振動を測定するための受信体は上記２０Ｈｚから２０ＫＨｚの大きく振幅された振動数が受信できる装置で有ればよく、特開平５−３２２８６１号公報に記述されている工夫された受信体であることに制限するものではなく、また通常の振動ピックアップ、加速度計であってもよい。
【０１３７】
次にどの振動形式を測定の対象とするかの選択方法について説明する。
理想状態でのデータから得られる有限要素法による固有値解析で得た振動形式の例を図４、図６で説明したが、有限要素法による固有値解析から他にも多くの振動形式がある。振動形式と理想固有振動数は一対一の関係があるので、簡単のために、理想固有振動数で示すと、表９のコンクリート部材についての理想固有振動数は表１１に示すようになる。表１１に示す以外にも高次の振動数があるが、省略する。表１１で示す振動数のうちどれを選択するかとの判断が必要であるが、振動解析の精度を重視するために、比較的小さい振動数の振動を選択することが望ましい。
【０１３８】
また、比較的小さい振動形式の中からは図４、図６に示すように、図４は支持する場所が２点、図６では、支持する場所が４点であり、コンクリート部材の重量が１トンにも及ぶ場合があり、安全上の配慮からも選択される。
【０１３９】
次に、コンクリート部材の成形体表面を打撃し、該打撃によるコンクリート部材の成形体の振動を解析することで該コンクリート部材の健全度を指数化することを実施例で説明する。
【０１４０】
表９のコンクリート部材を、図６で特定されるように支持及び打撃し、大きく振幅させ、振動を受信した。コンクリート部材は５個製造した。理想状態での動弾性率、密度は表９に示す２９ＧＰａ、２２７０ｋｇ／ｍ3、理想固有振動数は表１０に示すとおり４９４Ｈｚである。この場合のコンクリート部材の理想形状常数は（理想形状常数Ｋ［ｍ2］）＝（動弾性率Ｅ0：２９×１０9［Ｐａ］）／｛（密度ρ0：２２７０［ｋｇ／ｍ3］）（理想固有振動数ｆ0：４９４［１／ｓｅｃ］）2｝より、Ｋ＝５２．３５［ｍ2］であり、これを用いて、実測した振動数から既に述べた次式で弾性率を算出した。Ｅm＝Ｋ・ρm・ｆm 2
尚、密度は何れも２２００［ｋｇ／ｍ3］であった。
【０１４１】
表１２に測定された振動数と上述の式で算出した動弾性率を示す。
表１２より、振動数、動弾性率が何れも理想固有振動数、理想状態での動弾性率より小さく、異なった値であった。この値を指数化した結果を表１３に示す。理想固有振動数あるいは理想状態での動弾性率を１とした比の値で示されている。
【０１４２】
表１３に示すようにコンクリート部材の表面を打撃し、該打撃によるコンクリート部材の振動を最大となるようにコンクリート部材の振動数を測定する本発明によると、容易に該コンクリート部材の健全度を指数化することができる。
【０１４３】
表１３では理想固有振動数、理想状態での動弾性率を１とした比で示したが、製造の技術管理上はいくつか成形体を測定した結果で、最大のものを１とするような比で示すことや、測定された振動数そのものや、算出した弾性率そのものの値を指数にすることを制限するものではない。尚、弾性率に比べ振動数の場合は健全度が圧縮されて表現される傾向であり、理想状態での動弾性率との比を採用することが好ましい。
【０１４４】
本発明の目的であるコンクリート部材の非破壊検査により健全度を指数化することは表１３に示された。
【０１４５】
【表９】
【０１４６】
【表１０】
【０１４７】
【表１１】
【０１４８】
【表１２】
【０１４９】
【表１３】
【０１５０】
以上、耐火物及びコンクリート部材に関して本発明の説明を行ったが、本発明は、非固定の打音検査物であればこれらに限らず、いずれのものにも適用できる。
【０１５１】
【発明の効果】
本発明によれば、耐火物、コンクリート部材、その他の非固定打音検査物の非破壊検査により、耐火物、コンクリート部材、その他の非固定打音検査物の健全度を容易に指数化することでき、これにより組立部品等としての製造・出荷の品質を管理することで、製造技術との関連性を容易に把握でき、優れた製造技術の管理が実施できるばかりでなく、使用上の品質を把握して使用中の補修箇所、時期を予測する方法へ活用できるなどの面で有益な効果を奏する。
【図面の簡単な説明】
【図１】本発明に係る不定形耐火物成形体（又はコンクリート部材）の非破壊検査装置の一実施例を示す概略図
【図２】受信体が複数ある場合の本発明に係る不定形耐火物成形体（又はコンクリート部材）の非破壊検査装置の一実施例を示す概略図
【図３】（ａ）不定形耐火物成形体（又はコンクリート部材）の実施例の平面図
（ｂ）不定形耐火物成形体（又はコンクリート部材）の実施例の側面図
【図４】固有値解析（第一の例）
（ａ）有限要素法による固有値解析で第一の例（支持点２ヶ）の平面図
（ｂ）有限要素法による固有値解析で第一の例（支持点２ヶ）の立面図
（ｃ）有限要素法による固有値解析で第一の例（支持点２ヶ）の側面図
【図５】固有値解析（第一の例）において、特定された支持する場所、打撃する場所、受信する場所を示す図
【図６】固有値解析（第二の例）
（ａ）有限要素法による固有値解析で第二の例（支持点４ヶ）の平面図
（ｂ）有限要素法による固有値解析で第二の例（支持点４ヶ）の立面図
（ｃ）有限要素法による固有値解析で第二の例（支持点４ヶ）の側面図
【図７】固有値解析（第二の例）において、特定された支持する場所、打撃する場所、受信する場所を示す図
【図８】Ａ材質の不定形耐火物成形体を図７に示すように支持及び打撃し、振動を受信した場合の解析結果を示す図
（ａ）打撃体がゴムの場合
（ｂ）打撃体がプラスチックの場合
（ｃ）打撃体が鋼球の場合
【図９】Ａ材質の不定形耐火物成形体を図７に示すように支持及び打撃し、振動を受信した場合の解析結果を示す図
（ａ）受信から５０ミリ秒の受信振動を解析対象から除いた解析結果
（ｂ）２個の受信体が受信した振動の位相関係を示す図
【図１０】打音による成形体の健全度を表す指数で品質を管理する一例の図
【符号の説明】
１ 非破壊検査装置
２ 支持体
３ 不定形耐火物成形体（又はコンクリート部材）
４ 打撃体
５ 受信体
５ａ 受信体
５ｂ 受信体
６ 解析装置
【発明の属する技術分野】
本発明は、打音検査物（例えば、製鉄、製鋼において使用される転炉、取鍋等の工業窯炉および樋など、またアルミニューム、ダクタイル鋳鉄など溶融金属を扱う産業での取鍋等の工業窯炉および樋などの耐火物、特に不定形耐火物、或いは、土木、建築構造物の建設に用いられるシールドセグメント、ヒューム管、コンクリート杭、消波ブロックなどのプレキャストコンクリート製品の如きコンクリート部材）の健全度を検査、管理する打音検査物の非破壊検査方法、その非破壊検査装置及び品質管理方法に関する。
【０００２】
【従来の技術】
打音検査物の非破壊検査に関し、例えば、コンクリート表面を打撃してコンクリート自体を伝搬する打撃音を採取することでコンクリート構造物内部の亀裂、空隙等の反映を健全度として検査する方法が特開平５−３２２８６１号公報に記述されている。
【０００３】
これによると、打撃音の受信方法を工夫することで内部空隙のある状況を受信音の周波数スペクトル解析によって周波数の振幅の大小で健全度を把握できるものとされている。空隙の有無を判定できる限界例は２０ｃｍの幅の空間を実施例で示している。また、この方法は健全度の判定方法として振幅の大小が例示的に示されているのみであり、健全度がどの程度かという指数化については言及されていない。更に、コンクリート構造物はベースコンクリート上に打設されたとの記述があるように、地上に一体固定された構造体を念頭に置いた非破壊検査方法であった。以上のように、構造体が固定されて、空隙の大きさもセンチメートルオーダーで壁状であり、健全性は指数化されていない。
【０００４】
一方、耐火物、特に不定形耐火物は、耐火煉瓦に比べ使用現場での施工作業性が良好であることから、多くの場合、現場での流し込みなどで現場一体型施工が実施されてきたが、近年、品質向上の観点から、現場一体型施工ではなく、分割された成形体を事前に製造し、それを組立て、施工する方法が取り入れられて来ている。組立部品としての不定形耐火物成形体の重量は１０００ｋｇ程度となる場合がある。
【０００５】
現場一体型施工を実施していた従来から、不定形耐火物は成形時にラミネーション、空隙、素材の偏析などによる欠陥が発生しやすく、その解決策の一つとして組立施工が導入され、前記問題点は減少の傾向にある中で、組立部品の製造技術の向上、出荷管理による現場使用時の補修方法への反映などの観点から、組立部品としての不定形耐火物成形体個々の健全度を指数化して把握する方法が望まれていた。
【０００６】
不定形耐火物の欠陥を把握する方法は特開平７−２１８４７９号公報に不定形耐火物施工体の非破壊試験方法と装置が示されている。これによると、施工体に超音波を送信し、施工体内部にある欠陥による反射波を受信することが記述されている。
【０００７】
しかしながら、この方法では不定形耐火物成形体としての健全度を把握するには計測に時間が掛かり健全度を指数化する方法としては問題があった。前記特開平５−３２２８６１号公報に記述されている技術は、前述の如く、健全度を指数化する方法が不明確で、把握できる健全性も空隙の大きさがセンチメートルオーダーとみられ、不定形耐火物成形体に適用するには不適当と判断されること、また、固定されたことを前提としたコンクリート構造体での健全度把握方法であり、固定されていないと言う自由度をもつ不定形耐火物成形体へ適用する技術としては問題があった。
【０００８】
【発明が解決しようとする課題】
そこで本発明の目的は、打音検査物（例えば、耐火物、特に不定形耐火物、コンクリート部材、その他の非固定物）の非破壊検査において、センチメートルオーダー以下の空隙、亀裂まで短時間に検出可能な検査方法、その装置を提供し、更に不定形耐火物成形体の如き非固定打音検査物の健全度を指数化することであり、これにより組立部品等としての製造・出荷の品質を管理する手段を提供することにある。
【０００９】
【課題を解決するための手段】
本発明者らは、これら従来の問題を解決するため、不定形耐火物やコンクリート部材が固定されていないことに着想を得て工夫された耐火物、特に不定形耐火物成形体やプレキャストコンクリート製品の如きコンクリート部材、その他の非固定打音検査物の非破壊検査方法とその装置を完成した。前記目的を達成するために、本発明は次のように構成する。
【００１０】
１．打音検査物表面を打撃した際の振幅が最大となる振動数を測定し、当該振動数と打音検査物の理想固有振動数と比較し、打音検査物の健全度を指数化する非破壊検査方法であって、
前記打音検査物の振動数の測定が、前記打音検査物を有限要素に分割し、有限要素法による固有値解析を実施することにより打音検査物の振動形式を算出し、当該振動形式が得られるように、打音検査物を支持する場所、打撃する場所及び受信する場所を特定して行われる構成であり、
前記理想固有振動数が、前記有限要素法による固有値解析において、打音検査物の理想状態での動弾性率、動ポアソン比及び密度を用いて固有値解析し振動数を算出した当該振動数を理想固有振動数とする構成であり、
前記理想状態での動弾性率及び動ポアソン比を求める手段が、仮の動ポアソン比から仮の動弾性率を求め、前記打音検査物の形状から動剪断弾性率を求め、当該仮の動弾性率と当該動剪断弾性率から、新たな仮の動ポアソン比を求め、当該新たな仮の動ポアソン比から新たな仮の動弾性率を求めて、これを第１回目の動弾性率とし、これを複数回繰り返して、前の回に求めた動弾性率との差が認められなかったときの動弾性率及び動ポアソン比を理想状態での動弾性率及び動ポアソン比とする構成であり、
打音検査物の振動数の受信において、受信振動のうち、打撃による打音検査物の振動開始から１〜１００ミリ秒の受信振動を解析の対象から除くこと、
を特徴とする打音検査物の非破壊検査方法。
【００１６】
３．打音検査物の品質を管理する品質管理方法において、上記１又は２に記載の方法により得られた打音検査物の前記健全度を表す指数を用いて、所定値を外れる場合に異常と判定することを特徴とする打音検査物の品質管理方法。
【００１７】
４．打音検査物が、耐火物又はコンクリート部材の如き非固定物であることを特徴とする上記１又は２に記載の打音検査物の非破壊検査方法。
【００１９】
５．打音検査物が、耐火物又はコンクリート部材の如き非固定物であることを特徴とする上記４に記載の打音検査物の品質管理方法。
【００２０】
本発明の参考発明は、下記の通りである。
６．打音検査物表面を打撃した際の振幅が最大となる振動数を測定し、当該振動数、打音検査物の密度及び理想形状常数から打音検査物の動弾性率を算出し、当該動弾性率と打音検査物の理想状態での動弾性率と比較し、打音検査物の健全度を指数化すること。
７．前記打音検査物表面を打撃する際に、打音検査物を支持する場所、打撃する場所、及び打音検査物の振動を受信する場所を特定すること。
８．打音検査物の健全度を指数化する非破壊検査装置において、打音検査物を支持する支持体、打音検査物表面を打撃する打撃体、打音検査物の振動を受信する１又は２以上の受信体及び前記受信体が受信した受信振動を解析する解析装置を含むこと。
【００２１】
９．前記打撃に使用する打撃体の振動数が２０Ｈｚ〜２０ＫＨｚであること。
【００２２】
１０．前記打撃による打音検査物の振動が２０Ｈｚ〜２０ＫＨｚであること。
【００２３】
１１．前記理想形状常数が打音検査物の理想固有振動数、動弾性率及び密度から算出されること。
【００２４】
１２．前記打音検査物の理想状態での動弾性率、動ポアソン比及び密度を、前記打音検査物と同一材質で製造された小形体の撓み振動数、捻り振動数、寸法及び重量の測定値に基づいて算出し、当該動弾性率、動ポアソン比及び密度を前記理想状態での動弾性率、動ポアソン比及び密度とすること。
【００２５】
１３．前記小形体の体積が０．０３ｍ3以下であること。
【００２９】
１４．前記複数の受信体は、打音検査物の振動が同位相、または１８０度位相の異なると特定される支持体上の該打音検査物に設置されていること。
【００３０】
本発明において、振動形式とは以下のように定義する。
【００３１】
均一な物質で与えられた形状において衝撃などの極短時間の外力により、振動する、所謂固有振動（非減衰自由振動）において、振動形式とは、その多数ある固有振動を形状の変化（所謂、振動の節、腹の位置などの変位（移動）の大小）する様式に着目し分類したものである。その振動形式を振動モード等と称することがある。
【００３２】
尚、或る形状の振動形式（モード）を求める一般的な方法は、例えば、その形状を有限要素に分割し、有限要素法による固有値解析を実施することにより得られる。この際、均一な物質として各要素に対して、密度、動弾性率、動ポアソン比の物性値を一定値として入れるが、振動形式はこの物性値に依らず、形状により定まる。
【００３３】
【発明の実施形態】
次に、本発明の実施態様について図を参照して詳細に説明する。併せて、本発明の参考発明についても説明する。
【００３４】
先ず、本発明の打音検査物が耐火物である場合について、主に説明する。
図１は本発明に係る不定形耐火物成形体３の非破壊検査装置１の構成を示す一実施例である。図１によると、不定形耐火物成形体３は支持体２により地面から隔てられ、打撃体４による打撃により不定形耐火物成形体３を振動させる。この打撃による振動は特定された振動形式で固有振動を大きく振幅させるよう不定形耐火物成形体３の支持体２の場所を特定し、また打撃する場所を特定する。更に、受信体５により受信する場所も固有振動を大きく振幅させた場所におく。このようにして得られた受信信号を解析装置６で解析することで不定形耐火物成形体３の健全度を指数化できる。
【００３５】
不定形耐火物成形体３の製造方法の例としては水、ピッチ等の流動材と不定形耐火原料を混合し、流し込み・振動などの方法で型へ入れ成形し、養生及び乾燥して、その後の型から取り出したものである。この成形体の重量は１０００ｋｇ程度となる場合がある。
【００３６】
図２は本発明に係る不定形耐火物成形体３に関し別の非破壊検査装置１の構成を示す一実施例で、受信体５が複数ある場合の例である。図２では受信体５が２個の場合（５ａ，５ｂ）を示している。複数の受信体５ａ，５ｂを支持体２上の不定形耐火物成形体３に設置するに際して、特定された振動形式で固有振動が同位相、または１８０度位相の異なると特定される場所へ設置する。このようにして得られた受信信号を解析装置６で解析することで不定形耐火物成形体３の健全度を指数化する際に、測定される振動数の確認ができる。
【００３７】
次に、実施例として溶鋼鍋の敷部を不定形耐火物成形体で分割した場合の例を示す。図３は溶鋼鍋の敷部を組み立てる不定形耐火物成形体の一つである。（ａ）は平面図、（ｂ）は側面図である。この容積は０．３３ｍ3であり、嵩密度３．０ｔ／ｍ3の材料を使用すると、この不定形耐火物成形体の重量は約１トンになる。この不定形耐火物成形体は水と不定形耐火原料を混合し、流し込み・振動の方法で型へ入れ成形し、室温で３日養生、１５０℃で５日間乾燥して、その後の型から取り出したものである。
【００３８】
耐火物を有限要素に分割し、有限要素法による固有値解析を実施することにより耐火物の振動形式を算出し、当該振動形式が得られるように、耐火物を支持する場所、打撃する場所及び受信する場所を特定する方法を説明する。
【００３９】
図４は図３の形状を有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した場合を示す第一の振動形式の例図である。振動形式（モード）の出力には形状のみが重要な入力データであり、密度、動弾性率、動ポアソン比は均一体と仮定した各々につき或る値（一定値）を入力することで求められる。固有値解析の計算アルゴリズムは逆べき乗法で計算した。形状は解析モデル形状とするため三次元ソリッド要素により要素分割を行った。実施例では要素数２４２４、節点数３１７６で実施した。要素・節点の相対的変位（移動）は色分けでき、変位（移動）量を見極める判断ができる。
【００４０】
図４（ａ）は平面図、（ｂ）は立面図、（ｃ）は側面図であり、特に図４（ｂ）により不定形耐火物成形体が中央を対象点に左右の両端部が捻れている状況が分かる。従って、中央部分底部は変位（移動）が小さい。また、左右の端部（角）の変位（移動）が大きいことが分かる。このように変位（移動）の大小から、有限要素法により算出された振動形式が得られるように、耐火物を支持する場所、打撃する場所、受信体を置く場所を特定すればよい。この結果から図５に示す如く、図４（ａ）〜（ｃ）で特定された振動形式で大きい振幅が得られるように、耐火物を支持する支持体２を置く場所（△印）、打撃体４により打撃する場所（□印）、また、受信体５ａ、５ｂを置く場所（○印）の例を示した。また、この図は受信体を２個置く場合を示している。この場合は、受信体を置く２つの場所の振動が１８０度位相の異なると特定される場所へ受信体を設置する場合を示している。図４の（ａ）、（ｂ）、（ｃ）から両受信体を置く場所の変位（移動）を観察すると、一方の受信体の場所が上昇方向になると、もう一方の受信体の場所が下降方向になることから、受信体の受信信号の位相が１８０度異なると特定される。
【００４１】
図５は例であり、これにより上記の特定する考え方を制限するものではなく、図４（ａ）、（ｂ）、（ｃ）から得られた変位（移動）の小さい、大きいと判断される場所から実施しやすい場所・方法を選べる。具体的に選ぶ場所を図４（ａ）、（ｂ）、（ｃ）に示す。耐火物を支持する支持体２を置く場所は図４でＡ、Ｂの線分で示される範囲で変位（移動）が小さく適しており、この範囲で選べる。また、打撃体４により打撃する場所、受信体５を置く場所は図４でＣ、Ｄ、Ｅ、Ｆで示される範囲で変位（移動）が大きく、各々打撃する場所、受信体５を置く場所として適している。
【００４２】
以上のように有限要素による固有値解析により得られた捻れを生じる振動形式の例を説明したが、耐火物の振動が点対称に捻れを生じる振動形式があることは、有限要素法による解析を行わなくとも幾何学的な対象性から判断できる場合もある。
【００４３】
しかし、次の例のように、より高次の振動形式で振動を大きく振幅させるように支持する場所、打撃する場所、また、その振動を大きい振幅で受信するように受信体を置く場所が特定される方法として有限要素法による固有値解析が好ましい。
【００４４】
図６はこの形状を有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した場合を示す第二の振動形式の例図である。振動形式を求める計算方法は図４、図５の場合と同様である。（ａ）は平面図、（ｂ）は立面図、（ｃ）は側面図である。有限要素法による解析によれば、耐火物内周側の端部であるＲの付いた角（底部）Ａ、Ｂを一つの固定端とし、外周部の端部からそれぞれ１／４の場所Ｃ、Ｄに変位（移動）の最小のところがあることがわかる。これは側面図（ｃ）のＲ付き角（底部）Ａ、Ｂから外周部に水平に線を引くと外周部の曲線をほぼ２分割することからも分かる。また、この点を結ぶ線上（Ａ−Ｄ、Ｂ−Ｃ）の変位（移動）も小さい。一方、中央部分と外周部の端部（角）の変位（移動）が大きいことが分かる。このように変位（移動）の大小から、算出された振動形式が得られるように、耐火物を支持する場所、打撃する場所、また、受信体を置く場所が特定される。この結果から図７に示す如く、図６（ａ）〜（ｃ）で特定された振動形式が得られるように、耐火物を支持する支持体２を置く場所（△印）、打撃体４により打撃する場所（□印）、また、受信体５ａ，５ｂを置く場所（○印）を示した。また、この図は受信体を２個置く場合を示している。この場合は、受信体を置く２つの場所の振動が同位相となると特定される場所へ受信体を設置する場合を示している。図６の（ａ）、（ｂ）、（ｃ）から両受信体を置く場所の変位（移動）を観察すると、一方の受信体の場所が上昇方向になると、もう一方の受信体の場所も上昇方向になることから、受信体の受信信号の位相が同位相になると特定される。図７は例であり、これにより上記の特定する考え方を制限するものではなく、図６（ａ）、（ｂ）、（ｃ）から得られた変位（移動）の小さい、大きいと判断される場所から実施しやすい場所・方法を選べる。具体的に選ぶ場所を図６（ａ）、（ｂ）、（ｃ）に示す。耐火物を支持する支持体２を置く場所は図６でＡ−Ｄ、Ｂ−Ｄの線分で示される範囲で変位（移動）が小さく適しており、この範囲で選べる。また、打撃体４により打撃する場所、受信体５を置く場所は図６でＥ、Ｆ、Ｇ、Ｈで示される範囲で変位（移動）が大きく、各々打撃する場所、受信体５を置く場所として適している。
【００４５】
次に、耐火物の理想状態での動弾性率、動ポアソン比、密度を算出する方法について説明する。
【００４６】
前記成形体は０．３３ｍ3と大きく、嵩密度３．０ｔ／ｍ3の材料を使用すると、この不定形耐火物成形体の重量は約１トンになる。このように大形状では成形上のラミネーション、空隙、素材の偏析などによる欠陥が生じる可能性が回避出来ないが、成形容積を小さくすることで、例えば１／１０倍の容積以下（０．０３ｍ3）に小さくすることで、成形上の前記欠陥を大きく回避できることに着目した。大きさ以外は組立部品とする不定形耐火物成形体の成形方法と同一とすることが好ましく、大きさの違いから流し込み成形の振動時間が異なることや、また、養生、乾燥の時間が異なることがあることなど、成形体の品質を変えない範囲で選定することを制限するものではない。本発明はこの方法により成形されたものを理想状態での成形体、理想状態成形体とした。
【００４７】
この理想状態成形体を打撃し、打撃による成形体の振動数を測定し、理想状態での動弾性率、動ポアソン比を算出した。理想状態成形体の一例の形状として、縦１５ｃｍ×横１５ｃｍ×長さ５３ｃｍの０．０１２ｍ3矩形体を作製した。この形状を複数個作製することで、理想状態でのバラツキを明確にすることができる。ここでは、３個（不定形耐火物Ａ）、あるいは１０個（不定形耐火物Ｂ）の成形体を作製し、重量を測定して寸法から密度を算出し、また撓み振動数、捻り振動数を測定し、それにより動弾性率(撓み)、動剪断弾性率、動ポアソン比を算出した。その後、平均値の動弾性率（撓み）、動ポアソン比、密度を求めた。この動弾性率（撓み）、動ポアソン比、密度が不定形耐火物成形体の理想状態での動弾性率、動ポアソン比、密度である。尚、ＪＩＳ Ａ １１２７−１９９３には矩形体など簡単な形状での撓み振動数、捻り振動数の測定方法及び動弾性率(撓み)、動剪断弾性率の計算方法が示されており、今回の成形体を打撃し、打撃による成形体の振動を測定、弾性率を計算する場合においてそれに倣っており、以下に要点を説明する。
【００４８】
表１に不定形耐火物Ａに関し縦１５ｃｍ×横１５ｃｍ×長さ５３ｃｍの矩形体の撓み振動数、捻り振動数を測定した結果と計測値の密度を示す。この結果から動弾性率（撓み）Ｅd（Ｎ／ｍｍ2）、動剪断弾性率Ｇd（Ｎ／ｍｍ2）と動ポアソン比νd［−］を求める方法を説明する。
【００４９】
動弾性率（撓み）Ｅdは、長さＬ、幅さｂ、高さｔ、（本例では長さＬ＝５３０ｍｍ、幅さｂ＝１５０ｍｍ、高さｔ＝１５０ｍｍ）、質量ｍ（ｋｇ）、密度［ｋｇ／ｍ3］で、（密度は計測値、振動数ｆ（１／ｓｅｃ）は測定値）を用いて次式で求められる。
尚、Ｔは修正係数である。本実施例では矩形体であるので、t／（３．４６４Ｌ）と動ポアソン比νdの関数で与えられる。動ポアソン比νdを１／６とした場合の修正係数Tは表２で示されている。例えば、ｔ＝１５０［ｍｍ］、Ｌ＝５３０［ｍｍ］、νd＝１／６の場合、Ｔ＝１．５である。
【００５１】
更に、そのＴを用いて、動ポアソン比がνdである場合の修正係数Ｔ'は次式で求められる。
T'＝T(1+(0.26νd +3.22 (νd)2)(t/(3.464L))/(1+0.1328(t/(3.464L))
【００５２】
動弾性率（撓み）Ｅdを求める実際の手順は、先ず、動ポアソン比νdを１／６と仮定して、仮の動弾性率（撓み）Ｅdを求め、後に述べる手順で求める動剪断弾性率Ｇdから得られる動ポアソン比νdを用いて再度計算し直す。これが第１回目の値となり、これを繰り返すことで精度のある動弾性率（撓み）Ｅdが求められる。
【００５３】
本例の場合は動ポアソン比を求める試行を６回繰り返した。１回目の値と２回目は有意な差が認められ、２回目から６回目までは差が認められなかった。その結果を表３に示す。
【００５４】
動剪断弾性率Ｇd（Ｎ／ｍｍ2）は次式で与えられる。
ここで、Ｒは形状係数と呼ばれ、ｔ／ｂの関数（ｔ≦ｂ）であり、正方形断面（本実施例は縦１５０ｍｍ×横１５０ｍｍ）の角柱の場合は、Ｒ＝１．１８３である。そこで、動弾性率（歪み）Ｅd、動剪断弾性率Ｇdより、動ポアソン比νdは
νd＝｛Ｅd／（２Ｇd）｝−１で求められる。
【００５６】
上記と同様に、不定形耐火物Ｂについても密度、撓み振動数、捻り振動数を測定し、動弾性率（歪み）Ｅd、動剪断弾性率Ｇd、動ポアソン比νdを算出し、更にその平均値を求め、理想状態での動弾性率、動ポアソン比、密度を算出した。
【００５７】
前記の理想状態での動弾性率、動ポアソン比、密度を表４に示す。
【００５８】
尚、前記の成形体を打撃し、打撃による成形体の振動を測定する方法によらず、例えば超音波の透過時間から動弾性率を求める方法もあり、このような他の方法によって決めることを制限しない。但し、この場合は動ポアソン比を経験値などで仮定することになる。また、前記の成形体を打撃し、打撃による成形体の振動を測定し、動弾性率（撓み）を算出する実施例を示したが、縦振動を測定して動弾性率（縦）を算出し動弾性率とすることを制限するものではない。本発明の視点は小形成形体が理想状態との観点であり、その形状からの動弾性率を測定する方法を制限するものではない。
【００５９】
【表１】
【表２】
【表３】
次に、耐火物の理想固有振動数を算出する方法を説明する。
有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した。形状の入力データのほか、表４の理想状態での密度、動弾性率、動ポアソン比を各々要素につき入力することで理想固有振動数が求められる。固有値解析の計算アルゴリズムは逆べき乗法で計算した。形状は図４と同様に三次元ソリッド要素により分割を行い、実施例では要素数２４２４、節点数３１７６で実施した。不定形耐火物Aでは、実際のデータを密度は要素No.１から要素No.２４２４の全要素を２２７０（ｋｇ／ｍ3）、動弾性率は同じく全要素を２９．０Ｅ９（意味は２９．０×１０9）（Ｐａ）、動ポアソン比は同じく全要素を０．２５とした。不定形耐火物Bでは、実際のデータを密度は要素No.１から要素No.２４２４の全要素を３３４０（ｋｇ／ｍ3）、動弾性率は同じく全要素を１４９Ｅ９（意味は１４９×１０9）（Ｐａ）、動ポアソン比は同じく全要素を０．１４とした。
【００６３】
このようにして有限要素法（ＦＥＭ）による固有値解析を実施して得られた理想固有振動数を表５に示す。このように、理想状態での諸物性（密度、動弾性率、動ポアソン比）と図４、図６など振動形式を特定することで理想固有振動数が決まる。
【００６４】
この表５により分かるように図４、図６の理想固有振動数は不定形耐火物の種類Ａ、Ｂにより大きく異なる。一方で、耐火物の種類Ａ、Ｂによる理想固有振動数の比（Ａ／Ｂ）は図４、図６のそれぞれの振動形式でほぼ一定（０．５３、０．５２）であり、また図４、図６の振動形式による振動数の比（（１）／（２））はＡ、Ｂでほぼ一定（０．８８、０．８９）の関係が得られる。
【００６５】
表５の同一の振動形式でＡ、Ｂの理想状態での物性値に基づく理想固有振動数の変化（例えば図６の振動形式で４９４Ｈｚ、９５０Ｈｚ）は、見方を変えると、これは耐火物の振動数を測定した結果に対する解釈を与えている。例えば、耐火物Ｂについて、振動数を測定した結果、理想固有振動数（９５０Ｈｚ）より小さい振動数（８３０Ｈｚ）になる場合、密度、動弾性率、動ポアソン比が反映されて、理想状態でのＢの諸物性からＡのほうへずれたためと解釈される。この考えに基づき、該耐火物の振動数を測定し、該耐火物の理想固有振動数と比較し健全度として、測定振動数を指数化する。
【００６６】
Ａ、Ｂ間の固有振動数の比がほぼ一定（０．５３、０．５２）であるので図４、図６のいずれの振動形式で評価しても、同等の健全度を表すことが分かる。
【００６７】
測定された振動数と理想固有振動数との比で健全度を指数化する方法としては、上述のようであるが、測定された振動数から算出された動弾性率と理想状態での動弾性率との比が健全度の指数として好ましい。
【００６８】
そこで、測定された振動数から動弾性率の算出方法について説明する。
理想状態での動弾性率Ｅ0［Pa］と密度ρ0［ｋｇ／ｍ3］、理想固有振動数ｆ0［１／sec］が次式の関係にある。
（Ｅ0）=(Ｋ)×（ρ0）×（ｆ0）2
【００６９】
ここでＫを理想形状常数［ｍ2］と称する。この理想形状常数は形状、動ポアソン比に関係する常数である。
【００７０】
この関係式により、理想形状常数Ｋを算出し、耐火物の実測された密度ρm、測定された振動数ｆmを入れて、弾性率Ｅmが求められる。
【００７１】
以上の説明をまとめると
不定形材料Ｂの図６における理想形状常数は次式より、
(理想形状常数Ｋ[m2])＝（動弾性率Ｅ0：１４９×１０9）／｛（密度ρ0：３３４０）（理想固有振動数ｆ0：９５０）2｝
Ｋ：４９．４[m2]である。
【００７３】
弾性率Ｅmの一例を計算すると、密度ρm＝３２５０［ｋｇ／ｍ3］、測定された振動数ｆm＝８３２［１／sec］を代入して、
【００７４】
健全度の指数として測定された振動数と理想固有振動数との比、動弾性率と理想状態での動弾性率との比ではなく、測定され振動数や算出された動弾性率そのものを使うことを制限するものではない。いずれにしても、理想固有振動数、理想状態での弾性率が想定されることが好ましい。
【００７５】
【表４】
【表５】
次に、耐火物の表面を打撃し、該打撃による耐火物の振動数を測定する範囲、打撃体の振動数について説明する。
【００７８】
不定形耐火物成形体の振動形式と理想固有振動数は有限要素法（ＦＥＭ）により固有値解析を実施して求めることは既に説明した。そこでは理想状態での物性値（密度、動弾性率、動ポアソン比）を適宜変更することで、予測される振動が算出され、表４に示すような代表的な理想状態での物性からは、図４、図６の場合で大型形状の理想固有振動数は表５に示すようにせいぜい１ＫＨｚである。一般的に、大型形状即ち、長尺なものになると周波数が小さくなる。一方、耐火物で製造する代わりに、金属で代替し、例としてアルミニューム金属製で図３の１／４サイズで約１４Ｋｇの場合で実測すると、図４、図６の振動形式の理想固有振動数はせいぜい３ＫＨｚである。
【００７９】
また、振動には倍振動、３倍振動などもあり、打撃による成形体の振動が２０Ｈｚから２０ＫＨｚにあると限定した。
【００８０】
また、この振動を大きく振幅させるためには、打撃に使用する打撃体の加振振動数が２０Ｈｚから２０ＫＨｚにあることが必要であり、ゴムは２０Ｈｚ〜５００Ｈｚ、プラスチックは５００Ｈｚ〜３Ｋｚ、３ＫＨｚ以上はスチールがその振動数の振幅を大きくさせるために適当な材質の例である。
【００８１】
従って、その振動を測定するための受信体は上記２０Ｈｚから２０ＫＨｚの大きく振幅された振動数が受信できる装置で有ればよく、特開平５−３２２８６１号公報に記述されている工夫された受信体であることに制限するものではなく、また通常の振動ピックアップ、加速度計であってもよい。
【００８２】
次にどの振動形式を測定の対象とするかの選択方法について説明する。
理想状態でのデータから得られる有限要素法による固有値解析で得た振動形式の例を図４、図６で説明したが、有限要素法による固有値解析から他にも多くの振動形式がある。振動形式と理想固有振動数は一対一の関係があるので、簡単のために、理想固有振動数で示すと、表４の不定形耐火物Ａについての理想固有振動数は表６に示すようになる。表６に示す以外にも高次の振動数があるが、省略する。表６で示す振動数のうちどれを選択するかとの判断が必要であるが、振動解析の精度を重視するために、比較的小さい振動数の振動を選択することが望ましい。
【００８３】
また、比較的小さい振動形式の中からは図４、図６に示すように、図４は支持する場所が２点、図６では、支持する場所が４点であり、成形体の重量が１トンにも及ぶ場合があり、安全上の配慮からも選択される。
【００８４】
【表６】
このような観点から、以下の実施例では図６で特定される振動を大きく振幅させ、その振動を解析した例を説明する。
【００８６】
図６に、表４の種類Ａにより製造した不定形耐火物成形体を、図７で特定されるように支持及び打撃し、大きく振幅させ、振動を受信した場合の解析結果を示す。 解析装置は市販の高速フーリエ変換装置（ＦＦＴ）を内蔵した解析装置である。
【００８７】
図６（ａ）は、打撃体がゴムの場合、図６（ｂ）は、打撃体がプラスチックの場合、図６（ｃ）は、打撃体が鋼球の場合で、横軸が周波数、縦軸がその振幅を示す。理想固有振動数４９４Ｈｚ付近の測定値は５００Ｈｚで、これは解析装置の分解能による誤差が原因で５００Ｈｚと測定された。理想固有振動数４９４Ｈｚを大きく振幅させる打撃方法としてゴム、プラスチックが適当な材質であることが分かる。何れの場合も理想固有振動数付近に明確な振動数が解析されるが、その他にも１ＫＨｚ、１．４ＫＨｚなどに共通のピークがある。
【００８８】
これら高次の振動数は、不定形耐火物成形体の健全度を指数化する本非破壊検査方法にとって意味がなく、除外することで解析精度が向上したり、測定振動が特定された振動を測定しているか判断する際にノイズになるので除外する。
【００８９】
この高次の振動数を除外する方法としては、打撃による成形体の振動開始から１〜１００ミリ秒の受信振動を解析の対象から除き、その後の受信振動を解析することを発明した。また、請求項４に記載の発明において、成形体の理想固有振動数を含みながら最も小さい測定周波数を選択する機能を有することが好ましい。
【００９０】
前記請求項４に記載の発明においては、ＦＦＴ解析の以前に高次の振動数を測定しないことにより、ＦＦＴ解析によりスペクトル解析の精度が向上し、最も小さい測定周波数を選択する機能に１００Ｈｚ、２００Ｈｚ、５００Ｈｚ、１ＫＨｚ、２ＫＨｚ、５ＫＨｚ、１０ＫＨｚ、２０ＫＨｚから測定周波数を選択する機能を有することで解析装置が容易に使用できる。また、このような目的で特定の周波数以上のフィルターを掛けＦＦＴ解析の以前に高次の振動を測定しない方法もあり、これを制限するものではない。
【００９１】
高次の振動数を除外する方法の実施例として、Ａ材質の不定形耐火物成形体を図７で示すように支持及び打撃し、大きく振幅させ、受信した場合での実施例としての解析結果を図９に示す。この解析に当たり、成形体の振動開始から５０ミリ秒の受信振動を解析の対象から除外した。この実施例は、また２個の受信体を配置した例であり（前記８の発明）、また前記１４の好ましい実施態様として、成形体の振動が同位相になる場所へ受信体を配置した場合の実施例でもある。
【００９２】
図９（ａ）で５００Ｈｚの振動数が大きく振幅していることが確認できた。また、図８の１ＫＨｚ、１．４ＫＨｚの振動数が図９では小さく押さえられていることが分かる。この結果、成形体の振動が同位相と特定される場所へ設置した受信体の受信信号を解析した結果、振動が、図９（ｂ）で示すように同位相を主体にした状態が確認できるまでに向上した。
【００９３】
このように、打撃による成形体の振動開始から５０ミリ秒の受信振動を解析の対象から除き、その後の受信振動を解析することで、受信した振動が測定対象の振動であることを確認することが出来た。実施例では、打撃による成形体の振動開始から５０ミリ秒の受信振動を解析の対象から除いたが、この時間は、理想固有振動数により適宜変動させ、適正な値を選定することができる。だたし、受信振動を解析の対象から除く期間が長期間になると、解析対象の振動が小さくなり、誤差が生じ易いので、上限を１００ミリ秒とした。また、１ミリ秒未満では、外乱の受信振動を解析の対象から除くことが不充分である。従って、本発明において、打撃による成形体の振動開始から１〜１００ミリ秒の受信振動を解析の対象から除き、その後の受信振動を解析する機能を有することとした。
【００９４】
次に、不定形耐火物の成形体表面を打撃し、該打撃による不定形耐火物の成形体の振動を解析することで該不定形耐火物成形体の健全度を指数化することを実施例で説明する。
【００９５】
表４の種類Ｂの不定形耐火物成形体を、図７で特定されるように支持及び打撃し、大きく振幅させ、振動を受信した。種類Ｂの不定形耐火物成形体は５個製造した。理想状態での動弾性率、密度は表４Ｂに示す１４９ＧＰａ、３３４０ｋｇ／ｍ3、理想固有振動数は表５に示すとおり９５０Ｈｚである。この場合の不定形材料Ｂの理想形状常数は(理想形状常数Ｋ[m2])＝（動弾性率Ｅ0：１４９×１０9[Pa]）／｛（密度ρ0：３３４０[ｋｇ／m3]）（理想固有振動数ｆ0：９５０[1／sec]）2｝より、
Ｋ＝４９．４[m2]であり、これを用いて、実測した振動数から既に述べた次式で弾性率を算出した。Ｅm＝Ｋ・ρm・ｆm 2
尚、密度は何れも３２５０[ｋｇ／m3]であった。
【００９６】
表７に測定された振動数と上述の式で算出した動弾性率を示す。
【００９７】
【表７】
表７より、振動数、動弾性率が何れも理想固有振動数、理想状態での動弾性率より小さく、異なった値であった。この値を指数化した結果を表８に示す。理想固有振動数あるいは理想状態での動弾性率を１とした比の値で示されている。
【００９９】
【表８】
表８に示すように耐火物の表面を打撃し、該打撃による耐火物の振動を最大となるように耐火物の振動数を測定する本発明によると、容易に該不定形耐火物成形体の健全度を指数化することができる。
【０１０１】
表８では理想固有振動数、理想状態での動弾性率を１とした比で示したが、製造の技術管理上はいくつか成形体を測定した結果で、最大のものを１とするような比で示すことや、測定された振動数そのものや、算出した弾性率そのものの値を指数にすることを制限するものではない。尚、弾性率に比べ振動数の場合は健全度が圧縮されて表現される傾向であり、理想状態での動弾性率との比を採用することが好ましい。
【０１０２】
本発明の目的である耐火物の非破壊検査により健全度を指数化することは表８に示された。次に、これにより耐火物の品質を管理する方法の一例を示す。
【０１０３】
図１０は、表８の結果を母数５個としてワイブルプロットしたものである。図１０中の（ａ）、（ｂ）は表８中の（ａ）、（ｂ）と対応している。動弾性率に比べ振動数の場合は健全度が圧縮されて表現される傾向であり、理想固有振動数との比よりも理想弾性率との比を採用することが好ましいことが分かる。
【０１０４】
右側上方にあるほど理想に近く、左側下方にあるほど健全度が悪化していることを示す。理想から離れた品質がこのように順位付けられ、製造技術との関連性を容易に把握でき、優れた製造技術の管理が実施できるばかりでなく、使用上の品質を把握して使用中の補修箇所、時期を予測する方法へ活用できるなどの面で有益な作用を奏することができる。
【０１０５】
この他にも組立部品として品質を管理する統計上の方法はあり、それを制限するものではない。
【０１０６】
以上の実施態様においては、不定形耐火物成形体の非破壊検査方法、その非破壊検査装置及びそれによる品質管理方法を実施例に説明しているが、本発明は不定形耐火物成形体に限定されるものではなく、定形耐火物、ノズル煉瓦、炉底用カーボン煉瓦、炭化珪素、窒化珪素などのセラミック焼結体、各種セラミックスなどの固定されていない耐火物、成形体へも適用可能なものである。
【０１０７】
次に、本発明の打音検査物がコンクリート部材である場合について説明する。
尚、図は前記耐火物の説明で用いたものを使用する。
【０１０８】
図１は本発明に係るコンクリート部材３の非破壊検査装置１の構成を示す一実施例である。図１によると、コンクリート部材３は支持体２により地面から隔てられ、打撃体４による打撃によりコンクリート部材３を振動させる。この打撃による振動は特定された振動形式で固有振動を大きく振幅させるようコンクリート部材３の支持体２の場所を特定し、また打撃する場所を特定する。更に、受信体５により受信する場所も固有振動を大きく振幅させた場所におく。このようにして得られた受信信号を解析装置６で解析することでコンクリート部材３の健全度を指数化できる。
【０１０９】
図２は本発明に係るコンクリート部材３に関し別の非破壊検査装置１の構成を示す一実施例で、受信体５が複数ある場合の例である。図２では受信体５が２個の場合（５ａ，５ｂ）を示している。複数の受信体５ａ，５ｂを支持体２上のコンクリート部材３に設置するに際して、特定された振動形式で固有振動が同位相、または１８０度位相の異なると特定される場所へ設置する。このようにして得られた受信信号を解析装置６で解析することでコンクリート部材３の健全度を指数化する際に、測定される振動数の確認ができる。
【０１１０】
次に、実施例として或る形状のコンクリート部材を示す。（ａ）は平面図、（ｂ）は側面図である。この容積は０．３３ｍ3であり、嵩密度２．３ｔ／ｍ3の材料を使用すると、このコンクリート部材の重量は約０．７６トンになる。
【０１１１】
コンクリート部材を有限要素に分割し、有限要素法による固有値解析を実施することによりコンクリート部材の振動形式を算出し、当該振動形式が得られるように、コンクリート部材を支持する場所、打撃する場所及び受信する場所を特定する方法を説明する。
【０１１２】
図４は図３の形状を有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した場合を示す第一の振動形式の例図である。振動形式（モード）の出力には形状のみが重要な入力データであり、密度、動弾性率、動ポアソン比は均一体と仮定した各々につき或る値（一定値）を入力することで求められる。固有値解析の計算アルゴリズムは逆べき乗法で計算した。形状は解析モデル形状とするため三次元ソリッド要素により要素分割を行った。実施例では要素数２４２４、節点数３１７６で実施した。要素・節点の相対的変位（移動）は色分けでき、変位（移動）量を見極める判断ができる。
【０１１３】
図４（ａ）は平面図、（ｂ）は立面図、（ｃ）は側面図であり、特に図４（ｂ）によりコンクリート部材が中央を対象点に左右の両端部が捻れている状況が分かる。従って、中央部分底部は変位（移動）が小さい。また、左右の端部（角）の変位（移動）が大きいことが分かる。このように変位（移動）の大小から、有限要素法により算出された振動形式が得られるように、コンクリート部材を支持する場所、打撃する場所、受信体を置く場所を特定すればよい。
【０１１４】
この結果から図５に示す如く、図４（ａ）〜（ｃ）で特定された振動形式で大きい振幅が得られるように、コンクリート部材を支持する支持体２を置く場所（△印）、打撃体４により打撃する場所（□印）、また、受信体５ａ、５ｂを置く場所（○印）の例を示した。
【０１１５】
また、この図は受信体を２個置く場合を示している。この場合は、受信体を置く２つの場所の振動が１８０度位相の異なると特定される場所へ受信体を設置する場合を示している。図４の（ａ）、（ｂ）、（ｃ）から両受信体を置く場所の変位（移動）を観察すると、一方の受信体の場所が上昇方向になると、もう一方の受信体の場所が下降方向になることから、受信体の受信信号の位相が１８０度異なると特定される。
【０１１６】
図５は例であり、これにより上記の特定する考え方を制限するものではなく、図４（ａ）、（ｂ）、（ｃ）から得られた変位（移動）の小さい、大きいと判断される場所から実施しやすい場所・方法を選べる。具体的に選ぶ場所を図４（ａ）、（ｂ）、（ｃ）に示す。コンクリート部材を支持する支持体２を置く場所は図４でＡ、Ｂの線分で示される範囲で変位（移動）が小さく適しており、この範囲で選べる。また、打撃体４により打撃する場所、受信体５を置く場所は図４でＣ、Ｄ、Ｅ、Ｆで示される範囲で変位（移動）が大きく、各々打撃する場所、受信体５を置く場所として適している。
【０１１７】
以上のように有限要素による固有値解析により得られた捻れを生じる振動形式の例を説明したが、コンクリート部材の振動が点対称に捻れを生じる振動形式があることは、有限要素法による解析を行わなくとも幾何学的な対象性から判断できる場合もある。
【０１１８】
しかし、次の例のように、より高次の振動形式で振動を大きく振幅させるように支持する場所、打撃する場所、また、その振動を大きい振幅で受信するように受信体を置く場所が特定される方法として有限要素法による固有値解析が好ましい。
【０１１９】
図６はこの形状を有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した場合を示す第二の振動形式の例図である。振動形式を求める計算方法は図４、図５の場合と同様である。（ａ）は平面図、（ｂ）は立面図、（ｃ）は側面図である。有限要素法による解析によれば、コンクリート部材内周側の端部であるＲの付いた角（底部）Ａ、Ｂを一つの固定端とし、外周部の端部からそれぞれ１／４の場所Ｃ、Ｄに変位（移動）の最小のところがあることがわかる。これは側面図（ｃ）のＲ付き角（底部）Ａ、Ｂから外周部に水平に線を引くと外周部の曲線をほぼ２分割することからも分かる。また、この点を結ぶ線上（Ａ−Ｄ、Ｂ−Ｃ）の変位（移動）も小さい。一方、中央部分と外周部の端部（角）の変位（移動）が大きいことが分かる。
【０１２０】
このように変位（移動）の大小から、算出された振動形式が得られるように、コンクリート部材を支持する場所、打撃する場所、また、受信体を置く場所が特定される。この結果から図７に示す如く、図６（ａ）〜（ｃ）で特定された振動形式が得られるように、コンクリート部材を支持する支持体２を置く場所（△印）、打撃体４により打撃する場所（□印）、また、受信体５ａ，５ｂを置く場所（○印）を示した。また、この図は受信体を２個置く場合を示している。この場合は、受信体を置く２つの場所の振動が同位相となると特定される場所へ受信体を設置する場合を示している。
【０１２１】
図６の（ａ）、（ｂ）、（ｃ）から両受信体を置く場所の変位（移動）を観察すると、一方の受信体の場所が上昇方向になると、もう一方の受信体の場所も上昇方向になることから、受信体の受信信号の位相が同位相になると特定される。図７は例であり、これにより上記の特定する考え方を制限するものではなく、図６（ａ）、（ｂ）、（ｃ）から得られた変位（移動）の小さい、大きいと判断される場所から実施しやすい場所・方法を選べる。具体的に選ぶ場所を図６（ａ）、（ｂ）、（ｃ）に示す。
【０１２２】
コンクリート部材を支持する支持体２を置く場所は図６でＡ−Ｄ、Ｂ−Ｄの線分で示される範囲で変位（移動）が小さく適しており、この範囲で選べる。また、打撃体４により打撃する場所、受信体５を置く場所は図６でＥ、Ｆ、Ｇ、Ｈで示される範囲で変位（移動）が大きく、各々打撃する場所、受信体５を置く場所として適している。
【０１２３】
次に、コンクリート部材の理想固有振動数を算出する方法を説明する。
有限要素法（ＦＥＭ）の汎用プログラムを用いて（例えば、計算プログラム商品名：ＭＳＣ−ＭＡＲＣ）固有値解析のコマンドを実施した。形状の入力データのほか、表９の理想状態での密度、動弾性率、動ポアソン比を各々要素につき入力することで理想固有振動数が求められる。固有値解析の計算アルゴリズムは逆べき乗法で計算した。形状は図４と同様に三次元ソリッド要素により分割を行い、実施例では要素数２４２４、節点数３１７６で実施した。コンクリート部材では、実際のデータを密度は要素Ｎｏ．１から要素Ｎｏ．２４２４の全要素を２２７０（ｋｇ／ｍ3）、動弾性率は同じく全要素を２９．０Ｅ９（意味は２９．０×１０9）（Ｐａ）、動ポアソン比は同じく全要素を０．２５とした。
【０１２４】
このようにして有限要素法（ＦＥＭ）による固有値解析を実施して得られた理想固有振動数を表１０に示す。このように、理想状態での諸物性（密度、動弾性率、動ポアソン比）と図４、図６など振動形式を特定することで理想固有振動数が決まる。
【０１２５】
そこで、測定された振動数から動弾性率の算出方法について説明する。
理想状態での動弾性率Ｅ0［Ｐａ］と密度ρ0［ｋｇ／ｍ3］、理想固有振動数ｆ0［１／ｓｅｃ］が次式の関係にある。
（Ｅ0）＝（Ｋ）×（ρ0）×（ｆ0）2
【０１２６】
ここでＫを理想形状常数［ｍ2］と称する。この理想形状常数は形状、動ポアソン比に関係する常数である。
【０１２７】
この関係式により、理想形状常数Ｋを算出し、コンクリート部材の実測された密度ρm、測定された振動数ｆmを入れて、弾性率Ｅmが求められる。
【０１２８】
以上の説明をまとめると
コンクリート部材の図６における理想形状常数は次式より、
（理想形状常数Ｋ［ｍ2］）＝（動弾性率Ｅ0：２９×１０9）／｛（密度ρ0：２２７０）（理想固有振動数ｆ0：４９４）2｝
Ｋ：５２．３５［ｍ2］である。
【０１３０】
弾性率Ｅmの一例を計算すると、密度ρm＝２０００［ｋｇ／ｍ3］、測定された振動数ｆm＝３００［１／ｓｅｃ］を代入して、
【０１３１】
健全度の指数として測定された振動数と理想固有振動数との比、動弾性率と理想状態での動弾性率との比ではなく、測定され振動数や算出された動弾性率そのものを使うことを制限するものではない。いずれにしても、理想固有振動数、理想状態での弾性率が想定されることが好ましい。
【０１３２】
次に、コンクリート部材の表面を打撃し、該打撃によるコンクリート部材の振動数を測定する範囲、打撃体の振動数について説明する。
【０１３３】
コンクリート部材の振動形式と理想固有振動数は有限要素法（ＦＥＭ）により固有値解析を実施して求めることは既に説明した。そこでは理想状態での物性値（密度、動弾性率、動ポアソン比）を適宜変更することで、予測される振動が算出され、表９に示すような代表的な理想状態での物性からは、図４、図６の場合で大型形状の理想固有振動数は表１０に示すようにせいぜい１ＫＨｚである。一般的に、大型形状即ち、長尺なものになると周波数が小さくなる。一方、コンクリートで製造する代わりに、金属で代替し、例としてアルミニューム金属製で図３の１／４サイズで約１４Ｋｇの場合で実測すると、図４、図６の振動形式の理想固有振動数はせいぜい３ＫＨｚである。
【０１３４】
また、振動には倍振動、３倍振動などもあり、打撃による成形体の振動が２０Ｈｚから２０ＫＨｚにあると限定した。
【０１３５】
また、この振動を大きく振幅させるためには、打撃に使用する打撃体の加振振動数が２０Ｈｚから２０ＫＨｚにあることが必要であり、ゴムは２０Ｈｚ〜５００Ｈｚ、プラスチックは５００Ｈｚ〜３Ｋｚ、３ＫＨｚ以上はスチールがその振動数の振幅を大きくさせるために適当な材質の例である。
【０１３６】
従って、その振動を測定するための受信体は上記２０Ｈｚから２０ＫＨｚの大きく振幅された振動数が受信できる装置で有ればよく、特開平５−３２２８６１号公報に記述されている工夫された受信体であることに制限するものではなく、また通常の振動ピックアップ、加速度計であってもよい。
【０１３７】
次にどの振動形式を測定の対象とするかの選択方法について説明する。
理想状態でのデータから得られる有限要素法による固有値解析で得た振動形式の例を図４、図６で説明したが、有限要素法による固有値解析から他にも多くの振動形式がある。振動形式と理想固有振動数は一対一の関係があるので、簡単のために、理想固有振動数で示すと、表９のコンクリート部材についての理想固有振動数は表１１に示すようになる。表１１に示す以外にも高次の振動数があるが、省略する。表１１で示す振動数のうちどれを選択するかとの判断が必要であるが、振動解析の精度を重視するために、比較的小さい振動数の振動を選択することが望ましい。
【０１３８】
また、比較的小さい振動形式の中からは図４、図６に示すように、図４は支持する場所が２点、図６では、支持する場所が４点であり、コンクリート部材の重量が１トンにも及ぶ場合があり、安全上の配慮からも選択される。
【０１３９】
次に、コンクリート部材の成形体表面を打撃し、該打撃によるコンクリート部材の成形体の振動を解析することで該コンクリート部材の健全度を指数化することを実施例で説明する。
【０１４０】
表９のコンクリート部材を、図６で特定されるように支持及び打撃し、大きく振幅させ、振動を受信した。コンクリート部材は５個製造した。理想状態での動弾性率、密度は表９に示す２９ＧＰａ、２２７０ｋｇ／ｍ3、理想固有振動数は表１０に示すとおり４９４Ｈｚである。この場合のコンクリート部材の理想形状常数は（理想形状常数Ｋ［ｍ2］）＝（動弾性率Ｅ0：２９×１０9［Ｐａ］）／｛（密度ρ0：２２７０［ｋｇ／ｍ3］）（理想固有振動数ｆ0：４９４［１／ｓｅｃ］）2｝より、Ｋ＝５２．３５［ｍ2］であり、これを用いて、実測した振動数から既に述べた次式で弾性率を算出した。Ｅm＝Ｋ・ρm・ｆm 2
尚、密度は何れも２２００［ｋｇ／ｍ3］であった。
【０１４１】
表１２に測定された振動数と上述の式で算出した動弾性率を示す。
表１２より、振動数、動弾性率が何れも理想固有振動数、理想状態での動弾性率より小さく、異なった値であった。この値を指数化した結果を表１３に示す。理想固有振動数あるいは理想状態での動弾性率を１とした比の値で示されている。
【０１４２】
表１３に示すようにコンクリート部材の表面を打撃し、該打撃によるコンクリート部材の振動を最大となるようにコンクリート部材の振動数を測定する本発明によると、容易に該コンクリート部材の健全度を指数化することができる。
【０１４３】
表１３では理想固有振動数、理想状態での動弾性率を１とした比で示したが、製造の技術管理上はいくつか成形体を測定した結果で、最大のものを１とするような比で示すことや、測定された振動数そのものや、算出した弾性率そのものの値を指数にすることを制限するものではない。尚、弾性率に比べ振動数の場合は健全度が圧縮されて表現される傾向であり、理想状態での動弾性率との比を採用することが好ましい。
【０１４４】
本発明の目的であるコンクリート部材の非破壊検査により健全度を指数化することは表１３に示された。
【０１４５】
【表９】
【表１０】
【表１１】
【表１２】
【表１３】
以上、耐火物及びコンクリート部材に関して本発明の説明を行ったが、本発明は、非固定の打音検査物であればこれらに限らず、いずれのものにも適用できる。
【０１５１】
【発明の効果】
本発明によれば、耐火物、コンクリート部材、その他の非固定打音検査物の非破壊検査により、耐火物、コンクリート部材、その他の非固定打音検査物の健全度を容易に指数化することでき、これにより組立部品等としての製造・出荷の品質を管理することで、製造技術との関連性を容易に把握でき、優れた製造技術の管理が実施できるばかりでなく、使用上の品質を把握して使用中の補修箇所、時期を予測する方法へ活用できるなどの面で有益な効果を奏する。
【図面の簡単な説明】
【図１】本発明に係る不定形耐火物成形体（又はコンクリート部材）の非破壊検査装置の一実施例を示す概略図
【図２】受信体が複数ある場合の本発明に係る不定形耐火物成形体（又はコンクリート部材）の非破壊検査装置の一実施例を示す概略図
【図３】（ａ）不定形耐火物成形体（又はコンクリート部材）の実施例の平面図
（ｂ）不定形耐火物成形体（又はコンクリート部材）の実施例の側面図
【図４】固有値解析（第一の例）
（ａ）有限要素法による固有値解析で第一の例（支持点２ヶ）の平面図
（ｂ）有限要素法による固有値解析で第一の例（支持点２ヶ）の立面図
（ｃ）有限要素法による固有値解析で第一の例（支持点２ヶ）の側面図
【図５】固有値解析（第一の例）において、特定された支持する場所、打撃する場所、受信する場所を示す図
【図６】固有値解析（第二の例）
（ａ）有限要素法による固有値解析で第二の例（支持点４ヶ）の平面図
（ｂ）有限要素法による固有値解析で第二の例（支持点４ヶ）の立面図
（ｃ）有限要素法による固有値解析で第二の例（支持点４ヶ）の側面図
【図７】固有値解析（第二の例）において、特定された支持する場所、打撃する場所、受信する場所を示す図
【図８】Ａ材質の不定形耐火物成形体を図７に示すように支持及び打撃し、振動を受信した場合の解析結果を示す図
（ａ）打撃体がゴムの場合
（ｂ）打撃体がプラスチックの場合
（ｃ）打撃体が鋼球の場合
【図９】Ａ材質の不定形耐火物成形体を図７に示すように支持及び打撃し、振動を受信した場合の解析結果を示す図
（ａ）受信から５０ミリ秒の受信振動を解析対象から除いた解析結果
（ｂ）２個の受信体が受信した振動の位相関係を示す図
【図１０】打音による成形体の健全度を表す指数で品質を管理する一例の図
【符号の説明】
１ 非破壊検査装置
２ 支持体
３ 不定形耐火物成形体（又はコンクリート部材）
４ 打撃体
５ 受信体
５ａ 受信体
５ｂ 受信体
６ 解析装置
Claims (5)
- 打音検査物表面を打撃した際の振幅が最大となる振動数を測定し、当該振動数と打音検査物の理想固有振動数と比較し、打音検査物の健全度を指数化する非破壊検査方法であって、
前記打音検査物の振動数の測定が、前記打音検査物を有限要素に分割し、有限要素法による固有値解析を実施することにより打音検査物の振動形式を算出し、当該振動形式が得られるように、打音検査物を支持する場所、打撃する場所及び受信する場所を特定して行われる構成であり、
前記理想固有振動数が、前記有限要素法による固有値解析において、打音検査物の理想状態での動弾性率、動ポアソン比及び密度を用いて固有値解析し振動数を算出した当該振動数を理想固有振動数とする構成であり、
前記理想状態での動弾性率及び動ポアソン比を求める手段が、仮の動ポアソン比から仮の動弾性率を求め、前記打音検査物の形状から動剪断弾性率を求め、当該仮の動弾性率と当該動剪断弾性率から、新たな仮の動ポアソン比を求め、当該新たな仮の動ポアソン比から新たな仮の動弾性率を求めて、これを第１回目の動弾性率とし、これを複数回繰り返して、前の回に求めた動弾性率との差が認められなかったときの動弾性率及び動ポアソン比を理想状態での動弾性率及び動ポアソン比とする構成であり、
打音検査物の振動数の受信において、受信振動のうち、打撃による打音検査物の振動開始から１〜１００ミリ秒の受信振動を解析の対象から除くこと、
を特徴とする打音検査物の非破壊検査方法。 - 前記打音検査物の振動数の受信において、受信振動のうち、前記打音検査物の理想固有振動数を含みながら、最も小さい測定周波数を選択し、該選択の際に、１００Ｈｚ、２００Ｈｚ、５００Ｈｚ、１ＫＨｚ、２ＫＨｚ、５ＫＨｚ、１０ＫＨｚ、２０ＫＨｚの何れかから測定周波数を選択することを特徴とする請求項１に記載の打音検査物の非破壊検査方法。
- 打音検査物の品質を管理する品質管理方法において、請求項１又は２に記載の方法により得られた打音検査物の前記健全度を表す指数を用いて、所定値を外れる場合に異常と判定することを特徴とする打音検査物の品質管理方法。
- 打音検査物が、耐火物又はコンクリート部材の如き非固定物であることを特徴とする請求項１又は２に記載の打音検査物の非破壊検査方法。
- 打音検査物が、耐火物又はコンクリート部材の如き非固定物であることを特徴とする請求項３に記載の段検査物の品質管理方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001141731A JP4768927B2 (ja) | 2001-05-11 | 2001-05-11 | 打音検査物の非破壊検査方法、及び品質管理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001141731A JP4768927B2 (ja) | 2001-05-11 | 2001-05-11 | 打音検査物の非破壊検査方法、及び品質管理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002340869A JP2002340869A (ja) | 2002-11-27 |
| JP4768927B2 true JP4768927B2 (ja) | 2011-09-07 |
Family
ID=18988135
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001141731A Expired - Fee Related JP4768927B2 (ja) | 2001-05-11 | 2001-05-11 | 打音検査物の非破壊検査方法、及び品質管理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4768927B2 (ja) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4740748B2 (ja) * | 2006-01-16 | 2011-08-03 | 新日鉄エンジニアリング株式会社 | 耐火物煉瓦間の空隙検査方法 |
| JP2008020425A (ja) * | 2006-07-10 | 2008-01-31 | Marui:Kk | コンクリート構造物の非破壊診断方法と装置 |
| JP2010071748A (ja) * | 2008-09-17 | 2010-04-02 | Nippon Telegr & Teleph Corp <Ntt> | コンクリートポールの損傷検知方法 |
| JP2012017998A (ja) * | 2010-07-06 | 2012-01-26 | Shin Etsu Chem Co Ltd | 多結晶シリコン棒および多結晶シリコン棒の検査方法ならびに多結晶シリコン棒の製造方法 |
| JP5238762B2 (ja) | 2010-07-06 | 2013-07-17 | 信越化学工業株式会社 | 多結晶シリコン棒および多結晶シリコン棒の製造方法 |
| CN103862006B (zh) * | 2014-03-07 | 2016-06-08 | 攀钢集团攀枝花钢钒有限公司 | 一种判断连铸板坯的皮下裂纹缺陷的方法 |
| CN109208956A (zh) * | 2018-10-10 | 2019-01-15 | 中国水利水电科学研究院 | 一种抽水蓄能电站主厂房立柱动力设计方法 |
| JP6855646B1 (ja) * | 2019-12-10 | 2021-04-07 | 黒崎播磨株式会社 | スライディングノズルプレート用耐火物及びその製造方法 |
| US20240199095A1 (en) * | 2021-05-07 | 2024-06-20 | Nec Corporation | Abnormality detection device, abnormality detection method, and storage medium |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0356851A (ja) * | 1989-07-25 | 1991-03-12 | Toyota Motor Corp | 割れ検出方法 |
-
2001
- 2001-05-11 JP JP2001141731A patent/JP4768927B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002340869A (ja) | 2002-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10247718B2 (en) | Non-destructive apparatus, system and method for determining pull-out capacity of anchor bolts | |
| Ma et al. | Real-time quantification of fresh and hardened mechanical property for 3D printing material by intellectualization with piezoelectric transducers | |
| JP4768927B2 (ja) | 打音検査物の非破壊検査方法、及び品質管理方法 | |
| JP4980245B2 (ja) | プロセス容器のライニングの非破壊試験 | |
| JP4195176B2 (ja) | 高強度の鉄筋コンクリート部材の若材齢時におけるひび割れの判定方法、高強度の鉄筋コンクリート部材の若材齢時におけるひび割れ判定装置、高強度の鉄筋コンクリート部材の養生方法と養生期間を決定する装置および高強度の鉄筋コンクリートの打設方法 | |
| KR101291504B1 (ko) | 콘크리트 구속 수축 시험 장치 및 방법 | |
| KR101513067B1 (ko) | 비파괴 동탄성 계수 측정장치 및 이를 이용한 동탄성 계수 측정방법 | |
| CN109208657B (zh) | 一种低应变质量检测测具及其测试方法 | |
| CN111561970A (zh) | 一种钢混叠合柱的混凝土浇筑监测方法 | |
| CN115901457A (zh) | 筒压法测定抹灰砂浆抗压强度测强曲线的建立方法 | |
| Erdem et al. | Impact effect on different sized reinforced concrete specimens | |
| CN208334278U (zh) | 一种基于钻芯取样冲击回波的混凝土结构弹性模量检测装置 | |
| CN112729082A (zh) | 基于整体变形监测的实体构件外约束度评价方法 | |
| US8485033B2 (en) | Testing technique for determination of the dynamic poisson's ratio of a material with longitudinal and cross-sectional resonant frequencies of a circular solid rod specimen | |
| JP2016205832A (ja) | 推定方法 | |
| Kaur et al. | Feasibility of using resin-jacketed piezo sensors for monitoring of concrete strength using electro-mechanical impedance technique | |
| CN108931576B (zh) | 一种基于钻芯取样冲击回波的混凝土结构弹性模量检测方法 | |
| Wang et al. | Experimental study of effects of water-cement ratio on the acoustic emission rate “a” values in concrete | |
| CN217358821U (zh) | 筒柱式混凝土结构应力测量装置 | |
| JP6910891B2 (ja) | 乾式二重床の実建物における重量床衝撃音レベル低減量の推定方法 | |
| CN217466923U (zh) | 一种混凝土约束收缩环形新型抗裂试验装置 | |
| CN111829911B (zh) | 浮筑地面填充找平层的性能测定装置 | |
| KR101031165B1 (ko) | 초음파속도에 기초한 슬립폼 인양여부 평가방법 | |
| Staquet et al. | Ultrasonic monitoring of setting and autogenous shrinkage development of high performance concrete | |
| Buratti et al. | Experimental study on the flexural behaviour of fibre reinforced concretes strengthened with steel and macro-synthetic fibres |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070910 |
|
| RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20091020 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100427 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100623 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100810 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101012 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101214 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101216 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110614 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110617 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4768927 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140624 Year of fee payment: 3 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140624 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |