JP2013050377A - 寿命・余寿命評価機能付き非破壊検査装置 - Google Patents
寿命・余寿命評価機能付き非破壊検査装置 Download PDFInfo
- Publication number
- JP2013050377A JP2013050377A JP2011188390A JP2011188390A JP2013050377A JP 2013050377 A JP2013050377 A JP 2013050377A JP 2011188390 A JP2011188390 A JP 2011188390A JP 2011188390 A JP2011188390 A JP 2011188390A JP 2013050377 A JP2013050377 A JP 2013050377A
- Authority
- JP
- Japan
- Prior art keywords
- crack
- shape
- life
- remaining life
- defect
- 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.)
- Withdrawn
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
Abstract
【課題】数値シミュレーションに特化した専門の知識や計算機資源を必要とせず、簡易な操作、ノートPCレベルの計算速度と記憶容量だけでき裂進展挙動をシミュレーションし、破壊に到るまでの余寿命の算出を非破壊検査の現場で行うことを可能とする。
【解決手段】非破壊検査手段と寿命・余寿命評価手段とを備えた非破壊検査装置であって、寿命・余寿命評価手段は、非破壊検査手段により推定されたき裂面形状に基づいて、拡張型有限要素法により直交格子型の固定メッシュにて板厚方向に進展するき裂の応力解析を行い、構造部材が破断するまでの寿命・余寿命を評価する。
【選択図】図11
【解決手段】非破壊検査手段と寿命・余寿命評価手段とを備えた非破壊検査装置であって、寿命・余寿命評価手段は、非破壊検査手段により推定されたき裂面形状に基づいて、拡張型有限要素法により直交格子型の固定メッシュにて板厚方向に進展するき裂の応力解析を行い、構造部材が破断するまでの寿命・余寿命を評価する。
【選択図】図11
Description
本発明は、非破壊検査に関するものであって、特に、寿命・余寿命評価機能を備えた非破壊検査装置に関し、さらに、数値構造解析や信頼性工学の技術分野にも関するものである。
漏洩磁束密度計測等の非破壊検査手法によりき裂の検出を行おうとする試みはあるものの幾何形状を捉えるまでであり、部材の寿命・余寿命評価を直接与えるものではない(本発明者等による特許文献1及び非特許文献2参照)。
一方、寿命・余寿命評価に必要な構造解析を行うためには複雑な計算モデル(メッシュ)の作成が必要であり、モデル化・解析の専門家、あるいはソフトウェア等の人的、計算機的資源が必要とされる(特許文献2、3参照)。このため非破壊検査の現場にて寿命・余寿命評価や保守計画の立案を行うことが難しかった。
一方、寿命・余寿命評価に必要な構造解析を行うためには複雑な計算モデル(メッシュ)の作成が必要であり、モデル化・解析の専門家、あるいはソフトウェア等の人的、計算機的資源が必要とされる(特許文献2、3参照)。このため非破壊検査の現場にて寿命・余寿命評価や保守計画の立案を行うことが難しかった。
中住、鈴木、「磁気センサの分解能を考慮した逆解析に基づく欠陥形状の評価」、保全学、日本保全学会、vol.9、No.4(2011年)
非破壊検査の現場から得られたき裂形状のデータから直接的に、構造解析を経て寿命・余寿命評価まで実行する手段が存在しないため、現場から解析部署にデータを移した後に応力解析及び寿命・余寿命評価を行う手続きが現状となっている。そのことが保守計画全体の効率低下を招いている。
上記課題を解決するために、本発明では、拡張型有限要素法(XFEM)を用い、直交型メッシュにて有限要素解析を行う。これによりき裂形状を考慮したメッシュを作成する必要がなくなり、部材寸法・形状、き裂形状、荷重負荷条件等を入力することで、初心者でも容易に応力拡大係数解析と寿命・余寿命評価が可能となる。
本発明は、繰り返し荷重の負荷を受ける構造部材の欠陥を検出する非破壊検査手段と、欠陥形状入力手段と、寿命・余寿命評価手段を備えた寿命・余寿命評価機能付き非破壊検査装置であって、前記欠陥形状入力手段は、前記非破壊検査手段により検出された厚板状の構造部材の表面に発生した表面き裂の形状の入力、及び、構造部材の形状等の寸法・材料諸元、負荷条件の入力を行うものであり、前記寿命・余寿命評価手段は、入力された前記表面き裂形状を初期き裂縁形状とし、固定メッシュによる拡張型有限要素法を用いて、3次元応力解析を行い、その結果得られる周辺部の応力分布よりき裂縁上での応力拡大係数を求め、求めた応力拡大係数が予め定めた上限値を超えてなければ、Paris則により1回の負荷によるき裂縁上でのき裂進展量を求め、求めたき裂進展量によりき裂縁形状を更新し、以下、更新された裂縁形状に基づいてき裂縁上での応力拡大係数を求め、求めた応力拡大係数が予め定めた上限値を超えるまで繰り返し、当該上限値を超えるまでに要した繰り返し回数に基づいて寿命・余寿命を評価することを特徴とする。
また、本発明は、上記寿命・余寿命評価機能付き非破壊検査装置において、さらに、前記非破壊検査手段は、FGセンサにより測定した漏洩磁束密度分布と欠陥上に存在する磁荷との幾何学的対応関係を示す応答関数を用いて逆解析を行って磁荷の分布を復元させることにより磁性体の構造部材の欠陥形状を検出するものであって、前記応答関数として、FGセンサがセンサ長2Lを有することを考慮して修正した修正応答関数を用いることを特徴とする。
また、本発明は、非破壊検査手段と、欠陥形状入力手段と、寿命・余寿命評価手段を備えた寿命・余寿命評価機能付き非破壊検査装置を用いて寿命・余寿命を評価する方法であって、繰り返し荷重の負荷を受ける厚板状の構造部材の表面き裂形状を前記非破壊検査手段により検出するステップと、前記欠陥形状入力手段により、検出された前記表面き裂形状に基づいて初期き裂縁形状の入力、及び、前記構造部材の形状等の寸法・材料諸元、負荷条件の入力を行うステップと、
前記寿命・余寿命評価手段により固定メッシュを採用した拡張型有限要素法を用いて、
(a)前記初期き裂縁形状に基づき、3次元応力解析を行い、その結果得られる周辺部の応力分布よりき裂縁上での応力拡大係数を求めるステップと、
(b)求めた応力拡大係数が予め定めた上限値を超えてなければ、Paris則により1回の負荷によるき裂縁上でのき裂進展量を求めるステップと、
(c)求めたき裂進展量によりき裂縁形状を更新するステップと、
(d)以下、更新されたき裂縁形状に基づいて、上記(a)〜(c)を、応力拡大係数が予め定めた上限値を超えるまで繰り返し、当該上限値を超えるまでに要した繰り返し回数に基づいて寿命・余寿命を評価するステップと、
からなることを特徴とする。
また、本発明は、上記寿命・余寿命評価方法において、さらに、前記構造部材の欠陥形状を前記非破壊検査手段により検出するステップは、FGセンサにより測定した漏洩磁束密度分布と欠陥上に存在する磁荷との幾何学的対応関係を示す応答関数を用いて逆解析を行って磁荷の分布を復元させることにより磁性体の構造物の欠陥形状を評価し、前記応答関数として、FGセンサがセンサ長2Lを有することを考慮して修正した修正応答関数を用いることを特徴とする。
また、本発明は、非破壊検査装置の非破壊検査手段で厚板状の構造部材から検出された表面き裂形状に基づき欠陥形状入力手段により入力された初期き裂縁形状、構造部材の形状等の寸法・材料諸元、及び繰り返し荷重の負荷条件を読み取るステップと、
固定メッシュを採用した拡張型有限要素法を用いて、
(a)前記初期き裂縁形状に基づき、3次元応力解析を行い、その結果得られる周辺部の応力分布よりき裂縁上での応力拡大係数を求めるステップと、
(b)求めた応力拡大係数が予め定めた上限値を超えてなければ、Paris則により1回の負荷によるき裂縁上でのき裂進展量を求めるステップと、
(c)求めたき裂進展量によりき裂縁形状を更新するステップと、
(d)以下、更新されたき裂縁形状に基づいて、上記(a)〜(c)を、応力拡大係数が予め定めた上限値を超えるまで繰り返し、当該上限値を超えるまでに要した繰り返し回数に基づいて寿命・余寿命を評価するステップと、
を非破壊検査装置に行わせることを特徴とする寿命・余寿命評価プログラム。
本発明は、繰り返し荷重の負荷を受ける構造部材の欠陥を検出する非破壊検査手段と、欠陥形状入力手段と、寿命・余寿命評価手段を備えた寿命・余寿命評価機能付き非破壊検査装置であって、前記欠陥形状入力手段は、前記非破壊検査手段により検出された厚板状の構造部材の表面に発生した表面き裂の形状の入力、及び、構造部材の形状等の寸法・材料諸元、負荷条件の入力を行うものであり、前記寿命・余寿命評価手段は、入力された前記表面き裂形状を初期き裂縁形状とし、固定メッシュによる拡張型有限要素法を用いて、3次元応力解析を行い、その結果得られる周辺部の応力分布よりき裂縁上での応力拡大係数を求め、求めた応力拡大係数が予め定めた上限値を超えてなければ、Paris則により1回の負荷によるき裂縁上でのき裂進展量を求め、求めたき裂進展量によりき裂縁形状を更新し、以下、更新された裂縁形状に基づいてき裂縁上での応力拡大係数を求め、求めた応力拡大係数が予め定めた上限値を超えるまで繰り返し、当該上限値を超えるまでに要した繰り返し回数に基づいて寿命・余寿命を評価することを特徴とする。
また、本発明は、上記寿命・余寿命評価機能付き非破壊検査装置において、さらに、前記非破壊検査手段は、FGセンサにより測定した漏洩磁束密度分布と欠陥上に存在する磁荷との幾何学的対応関係を示す応答関数を用いて逆解析を行って磁荷の分布を復元させることにより磁性体の構造部材の欠陥形状を検出するものであって、前記応答関数として、FGセンサがセンサ長2Lを有することを考慮して修正した修正応答関数を用いることを特徴とする。
また、本発明は、非破壊検査手段と、欠陥形状入力手段と、寿命・余寿命評価手段を備えた寿命・余寿命評価機能付き非破壊検査装置を用いて寿命・余寿命を評価する方法であって、繰り返し荷重の負荷を受ける厚板状の構造部材の表面き裂形状を前記非破壊検査手段により検出するステップと、前記欠陥形状入力手段により、検出された前記表面き裂形状に基づいて初期き裂縁形状の入力、及び、前記構造部材の形状等の寸法・材料諸元、負荷条件の入力を行うステップと、
前記寿命・余寿命評価手段により固定メッシュを採用した拡張型有限要素法を用いて、
(a)前記初期き裂縁形状に基づき、3次元応力解析を行い、その結果得られる周辺部の応力分布よりき裂縁上での応力拡大係数を求めるステップと、
(b)求めた応力拡大係数が予め定めた上限値を超えてなければ、Paris則により1回の負荷によるき裂縁上でのき裂進展量を求めるステップと、
(c)求めたき裂進展量によりき裂縁形状を更新するステップと、
(d)以下、更新されたき裂縁形状に基づいて、上記(a)〜(c)を、応力拡大係数が予め定めた上限値を超えるまで繰り返し、当該上限値を超えるまでに要した繰り返し回数に基づいて寿命・余寿命を評価するステップと、
からなることを特徴とする。
また、本発明は、上記寿命・余寿命評価方法において、さらに、前記構造部材の欠陥形状を前記非破壊検査手段により検出するステップは、FGセンサにより測定した漏洩磁束密度分布と欠陥上に存在する磁荷との幾何学的対応関係を示す応答関数を用いて逆解析を行って磁荷の分布を復元させることにより磁性体の構造物の欠陥形状を評価し、前記応答関数として、FGセンサがセンサ長2Lを有することを考慮して修正した修正応答関数を用いることを特徴とする。
また、本発明は、非破壊検査装置の非破壊検査手段で厚板状の構造部材から検出された表面き裂形状に基づき欠陥形状入力手段により入力された初期き裂縁形状、構造部材の形状等の寸法・材料諸元、及び繰り返し荷重の負荷条件を読み取るステップと、
固定メッシュを採用した拡張型有限要素法を用いて、
(a)前記初期き裂縁形状に基づき、3次元応力解析を行い、その結果得られる周辺部の応力分布よりき裂縁上での応力拡大係数を求めるステップと、
(b)求めた応力拡大係数が予め定めた上限値を超えてなければ、Paris則により1回の負荷によるき裂縁上でのき裂進展量を求めるステップと、
(c)求めたき裂進展量によりき裂縁形状を更新するステップと、
(d)以下、更新されたき裂縁形状に基づいて、上記(a)〜(c)を、応力拡大係数が予め定めた上限値を超えるまで繰り返し、当該上限値を超えるまでに要した繰り返し回数に基づいて寿命・余寿命を評価するステップと、
を非破壊検査装置に行わせることを特徴とする寿命・余寿命評価プログラム。
現場において非破壊検査だけでなく、短時間で応力解析と寿命・余寿命評価まで行うことで、迅速な保守計画立案と実行が可能となる。
本発明は非破壊検査の現場で寿命・余寿命評価を行うものである。非破壊検査方法は特に限定しないが、例えば漏洩磁束密度法等が挙げられる。これによりき裂形状を推定することができる。
図1は、その概念図であり、繰り返し引張り荷重を受ける厚板状の構造部材に発生した表面き裂の形状概観を示す。
図1は、その概念図であり、繰り返し引張り荷重を受ける厚板状の構造部材に発生した表面き裂の形状概観を示す。
本発明の寿命・余寿命評価機能付き非破壊検査装置では、非破壊検査手段によって求まったき裂形状を、寿命・余寿命評価手段に入力する。本発明で対象とするき裂形状は、厚板状構造部材に発生する表面き裂(表面欠陥)である。通常これは半楕円形のき裂縁(き裂先端を連ねた曲線)を持つことが多いが、本発明では一般性を高めるため必ずしも半楕円形に限定せず、き裂縁を有限個の線分の結合で近似する。すなわちき裂縁は、有限個の線分と頂点とで構成される。例えば、非破壊検査手段によるき裂形状画像を見ながら、き裂縁上の点をマウス等の入力装置で順次選択していく。
図2は、その入力例を示すものであり、き裂面の存在する断面上において線分の結合によるき裂縁の近似モデル化を示す。このモデルを用いてき裂面及びき裂縁周辺の応力状態を3次元的に解析する。
図2は、その入力例を示すものであり、き裂面の存在する断面上において線分の結合によるき裂縁の近似モデル化を示す。このモデルを用いてき裂面及びき裂縁周辺の応力状態を3次元的に解析する。
従来、寿命・余寿命予測を行うためには、き裂形状を考慮した数値計算モデル(メッシュ)を作成し、これを有限要素法による構造解析を行うことが必要である。しかし一般に3次元構造物では、き裂形状を考慮したメッシュを作成するためには専門的な知識を持った人間による修正が不可欠であり完全な自動化が難しい。そのため構造解析を非破壊検査の現場で行うことは非常に困難であった。
そこで、本発明では拡張型有限要素法(XFEM)を使用することでこの問題を解決する。この拡張型有限要素法ではき裂の幾何学的情報を計算機内部で保有するためき裂形状を考慮したメッシュを作成する必要がなく、自動的に作成可能な直交格子型のメッシュを使用することができる。またき裂伝播解析ではメッシュの再作成が不要である。図3に有限要素法と本発明の拡張型有限要素法における解析メッシュの形状(2次元モデルの場合)を示す。
そこで、本発明では拡張型有限要素法(XFEM)を使用することでこの問題を解決する。この拡張型有限要素法ではき裂の幾何学的情報を計算機内部で保有するためき裂形状を考慮したメッシュを作成する必要がなく、自動的に作成可能な直交格子型のメッシュを使用することができる。またき裂伝播解析ではメッシュの再作成が不要である。図3に有限要素法と本発明の拡張型有限要素法における解析メッシュの形状(2次元モデルの場合)を示す。
拡張型有限要素法(XFEM)では要素内部の変位を複雑な関数で近似できる。これをエンリッチ(Enrich)と呼ぶ。本発明ではこの関数にヘビサイドのステップ関数と特異関数を用いる。
以下、本発明で用いるXFEMの定式化を説明する。エンリッチされる節点とき裂の位置関係を図4に示す。ある節点について、結合する形状関数がき裂によって完全に切断される場合をJ属性、形状関数が完全には切断されず内部にき裂端を含む場合をC属性と表記する。またき裂端を原点とし、き裂方向及びき裂に垂直な方向をそれぞれ
とする局所座標系
を定義し、全体座標系
と区別する。使用する有限要素を4節点アイソパラメトリック四角形要素とする。要素e内の変位場
は次式で表現される。
以下、本発明で用いるXFEMの定式化を説明する。エンリッチされる節点とき裂の位置関係を図4に示す。ある節点について、結合する形状関数がき裂によって完全に切断される場合をJ属性、形状関数が完全には切断されず内部にき裂端を含む場合をC属性と表記する。またき裂端を原点とし、き裂方向及びき裂に垂直な方向をそれぞれ
とする局所座標系
を定義し、全体座標系
と区別する。使用する有限要素を4節点アイソパラメトリック四角形要素とする。要素e内の変位場
は次式で表現される。
ここで、r及びθは、局所座標系
と図4に示す位置関係にある極座標系の座標値である。式(1)右辺の
は通常の有限要素法で現れる節点自由度ベクトル、
は、それぞれヘビサイド関数、特異関数がエンリッチされたことに対応してそれらの係数として付加されるベクトルである。
と図4に示す位置関係にある極座標系の座標値である。式(1)右辺の
は通常の有限要素法で現れる節点自由度ベクトル、
は、それぞれヘビサイド関数、特異関数がエンリッチされたことに対応してそれらの係数として付加されるベクトルである。
本発明では構造部材表面に発生した半楕円型表面き裂を解析する。解析を行うためにはき裂縁上の各点で、局所座標系を構築することが必要となる。その定義を図5に示す。また、図2に示したき裂モデルに対して、本発明装置内で作成された局所座標系ベクトル(主法線方向及び接線方向)をそれぞれ図6(主法線方向)、図7(接線方向)に示す。
本発明では構造解析に拡張型有限要素法を用いるため直交格子型のメッシュを扱う。図8に本発明の実施例で用いるメッシュを示す。き裂の形状と独立に作成するため完全に自動的に作成することができる。
本発明では構造解析に拡張型有限要素法を用いるため直交格子型のメッシュを扱う。図8に本発明の実施例で用いるメッシュを示す。き裂の形状と独立に作成するため完全に自動的に作成することができる。
部材に作用する荷重(引き張り荷重,曲げモーメント)の値を与え、応力解析を行う。すなわち、以下の方程式を解く。
ここで、
は上述の方法で定義した変位に基づいて計算した剛性行列
は節点自由度ベクトル、そして
は入力した荷重情報を元に節点荷重に換算した荷重ベクトルである。
ここで、
は上述の方法で定義した変位に基づいて計算した剛性行列
は節点自由度ベクトル、そして
は入力した荷重情報を元に節点荷重に換算した荷重ベクトルである。
応力拡大係数はJ積分を3次元に拡張した領域積分法で求める。J積分の被積分関数に、ある重み関数を掛け、その後にGaussの発散定理を適用し、そして境界がき裂端に収縮した極限を考えるものとすると、次式のように変形される。
ただし、σij,ui,xiはそれぞれ応力,変位,座標であり、またW,δ1jはそれぞれ仕事量,クロネッカーのデルタ関数である。またqは重み関数と呼ばれ、J積分の領域Vに分布する0から1までの値をとる関数である。J積分領域の形状を図9に示す。
上式で計算されたJ値を用いて応力拡大係数は次式のKIとして求まる。ただし、Eは構造材料のヤング率、LはJ積分領域のき裂縁に沿った長さである.
上述した方法により、初期のき裂形状に対してき裂縁上の各離散点での応力拡大係数の値を逐次算出し、き裂縁全体での分布を求める。この方法で求めた応力拡大係数は、通常、き裂縁に沿って多少変動する分布となる。そのため曲線による近似(カーブフィッティング)を行い、応力拡大係数分布の平滑化を図る。
本発明では、平滑化に用いる曲線として放物線を使用する。図2に示した初期き裂形状に対して、構造解析の結果求まる応力拡大係数とそのカーブフィッティングを図10に示す。この処理により、き裂縁上の任意の点で、応力拡大係数を算出することができる。
ただし、σij,ui,xiはそれぞれ応力,変位,座標であり、またW,δ1jはそれぞれ仕事量,クロネッカーのデルタ関数である。またqは重み関数と呼ばれ、J積分の領域Vに分布する0から1までの値をとる関数である。J積分領域の形状を図9に示す。
上式で計算されたJ値を用いて応力拡大係数は次式のKIとして求まる。ただし、Eは構造材料のヤング率、LはJ積分領域のき裂縁に沿った長さである.
上述した方法により、初期のき裂形状に対してき裂縁上の各離散点での応力拡大係数の値を逐次算出し、き裂縁全体での分布を求める。この方法で求めた応力拡大係数は、通常、き裂縁に沿って多少変動する分布となる。そのため曲線による近似(カーブフィッティング)を行い、応力拡大係数分布の平滑化を図る。
本発明では、平滑化に用いる曲線として放物線を使用する。図2に示した初期き裂形状に対して、構造解析の結果求まる応力拡大係数とそのカーブフィッティングを図10に示す。この処理により、き裂縁上の任意の点で、応力拡大係数を算出することができる。
得られた応力拡大係数に基づき、き裂の進展量を算出する。き裂の進展量の算出には、Paris則を用いる。すわなち、繰り返し荷重一回当たりのき裂進展量(き裂伝播速度)da/dNは次式で与えられる。
da/dN=C(ΔK)m
ここでaはき裂長さ、Nは繰り返し荷重の反復回数、ΔKは応力拡大係数範囲、C,mは材料パラメータである。応力拡大係数範囲ΔKは次式で求まる。
ΔK=Kmax−Kmin
Kmax,Kminはそれぞれき裂の開口時,閉口時における応力拡大係数である。
da/dN=C(ΔK)m
ここでaはき裂長さ、Nは繰り返し荷重の反復回数、ΔKは応力拡大係数範囲、C,mは材料パラメータである。応力拡大係数範囲ΔKは次式で求まる。
ΔK=Kmax−Kmin
Kmax,Kminはそれぞれき裂の開口時,閉口時における応力拡大係数である。
上述の方法により求めた裂進展量を元にき裂を進展させる。すなわち、き裂縁上の各離散点をその局所座標系の主法線方向(図6に示した方向)に進展量だけ移動させる。そして更新されたき裂形状に対して再度構造解析を行う。拡張型有限要素法を用いるためメッシュの更新は不要である。応力拡大係数が予め定まる限界値Kfcを超えた時点で構造部材が破断すると考え、K>Kfcを破断条件とする。
き裂進展解析のフローチャートを図11に示す。また、き裂進展解析の結果得られた、き裂縁の進展挙動を図12に示す。
き裂進展解析のフローチャートを図11に示す。また、き裂進展解析の結果得られた、き裂縁の進展挙動を図12に示す。
本発明では、き裂縁を線分結合で表現しており、それにより任意形状のき裂面への応用が可能である。図13に、2個の半円が部分的に結合した形状のき裂を示す。本発明では、このような複雑な形状のき裂であっても、応力拡大係数の算出が可能である。図14にその結果を示す。き裂の結合部分で高い応力拡大係数が生じていることが確認できる。
初期き裂縁形状の入力は、非破壊検査手段によるき裂形状画像を見ながら、き裂縁上の点をマウス等の入力装置で順次選択していく例を示したが、他の入力手段を用いることもできる。
構造部材の欠陥を検出する非破壊検査手段は、欠陥形状を検出できるものであればどのようなものでもよい。特に、FGセンサを用いた漏洩磁束密度分布法による場合には、本発明者等が先に提案(特許文献1参照)したFGセンサ長を考慮した修正応答関数を用いると、逆解析による欠陥形状がより正確に検出できる。
初期き裂縁形状の入力は、非破壊検査手段によるき裂形状画像を見ながら、き裂縁上の点をマウス等の入力装置で順次選択していく例を示したが、他の入力手段を用いることもできる。
構造部材の欠陥を検出する非破壊検査手段は、欠陥形状を検出できるものであればどのようなものでもよい。特に、FGセンサを用いた漏洩磁束密度分布法による場合には、本発明者等が先に提案(特許文献1参照)したFGセンサ長を考慮した修正応答関数を用いると、逆解析による欠陥形状がより正確に検出できる。
構造部材が受ける繰り返し荷重としては、例えば、始業と終業により毎日1回ずつ繰り返されるものや、自動車等が通過する毎に繰り返されるもの等、繰り返し負荷される荷重であれば適用できる。また、繰り返し時間間隔がわかっていれば、繰り返し時間間隔に繰り返し回数をかけることで寿命・余寿命を時間長さで評価することもできる。
Claims (5)
- 繰り返し荷重の負荷を受ける構造部材の欠陥を検出する非破壊検査手段と、欠陥形状入力手段と、寿命・余寿命評価手段を備えた寿命・余寿命評価機能付き非破壊検査装置であって、
前記欠陥形状入力手段は、前記非破壊検査手段により検出された厚板状の構造部材の表面に発生した表面き裂の形状の入力、及び、構造部材の形状等の寸法・材料諸元、負荷条件の入力を行うものであり、
前記寿命・余寿命評価手段は、入力された前記表面き裂形状を初期き裂縁形状とし、固定メッシュによる拡張型有限要素法を用いて、3次元応力解析を行い、その結果得られる周辺部の応力分布よりき裂縁上での応力拡大係数を求め、求めた応力拡大係数が予め定めた上限値を超えてなければ、Paris則により1回の負荷によるき裂縁上でのき裂進展量を求め、求めたき裂進展量によりき裂縁形状を更新し、以下、更新された裂縁形状に基づいてき裂縁上での応力拡大係数を求め、求めた応力拡大係数が予め定めた上限値を超えるまで繰り返し、当該上限値を超えるまでに要した繰り返し回数に基づいて寿命・余寿命を評価することを特徴とする寿命・余寿命評価機能付き非破壊検査装置。 - 前記非破壊検査手段は、FGセンサにより測定した漏洩磁束密度分布と欠陥上に存在する磁荷との幾何学的対応関係を示す応答関数を用いて逆解析を行って磁荷の分布を復元させることにより磁性体の構造部材の欠陥形状を検出するものであって、前記応答関数として、FGセンサがセンサ長2Lを有することを考慮して修正した修正応答関数を用いることを特徴とする請求項1記載の寿命・余寿命評価機能付き非破壊検査装置。
- 非破壊検査手段と、欠陥形状入力手段と、寿命・余寿命評価手段を備えた寿命・余寿命評価機能付き非破壊検査装置を用いて寿命・余寿命を評価する方法であって、
繰り返し荷重の負荷を受ける厚板状の構造部材の表面き裂形状を前記非破壊検査手段により検出するステップと、
前記欠陥形状入力手段により、検出された前記表面き裂形状に基づいて初期き裂縁形状の入力、及び、前記構造部材の形状等の寸法・材料諸元、負荷条件の入力を行うステップと、
前記寿命・余寿命評価手段により固定メッシュを採用した拡張型有限要素法を用いて、
(a)前記初期き裂縁形状に基づき、3次元応力解析を行い、その結果得られる周辺部の応力分布よりき裂縁上での応力拡大係数を求めるステップと、
(b)求めた応力拡大係数が予め定めた上限値を超えてなければ、Paris則により1回の負荷によるき裂縁上でのき裂進展量を求めるステップと、
(c)求めたき裂進展量によりき裂縁形状を更新するステップと、
(d)以下、更新されたき裂縁形状に基づいて、上記(a)〜(c)を、応力拡大係数が予め定めた上限値を超えるまで繰り返し、当該上限値を超えるまでに要した繰り返し回数に基づいて寿命・余寿命を評価するステップと、
からなることを特徴とする寿命・余寿命評価方法。 - 前記構造部材の欠陥形状を前記非破壊検査手段により検出するステップは、FGセンサにより測定した漏洩磁束密度分布と欠陥上に存在する磁荷との幾何学的対応関係を示す応答関数を用いて逆解析を行って磁荷の分布を復元させることにより磁性体の構造物の欠陥形状を評価し、前記応答関数として、FGセンサがセンサ長2Lを有することを考慮して修正した修正応答関数を用いることを特徴とする請求項3記載の寿命・余寿命評価方法。
- 非破壊検査装置の非破壊検査手段で厚板状の構造部材から検出された表面き裂形状に基づき欠陥形状入力手段により入力された初期き裂縁形状、構造部材の形状等の寸法・材料諸元、及び繰り返し荷重の負荷条件を読み取るステップと、
固定メッシュを採用した拡張型有限要素法を用いて、
(a)前記初期き裂縁形状に基づき、3次元応力解析を行い、その結果得られる周辺部の応力分布よりき裂縁上での応力拡大係数を求めるステップと、
(b)求めた応力拡大係数が予め定めた上限値を超えてなければ、Paris則により1回の負荷によるき裂縁上でのき裂進展量を求めるステップと、
(c)求めたき裂進展量によりき裂縁形状を更新するステップと、
(d)以下、更新されたき裂縁形状に基づいて、上記(a)〜(c)を、応力拡大係数が予め定めた上限値を超えるまで繰り返し、当該上限値を超えるまでに要した繰り返し回数に基づいて寿命・余寿命を評価するステップと、
を非破壊検査装置に行わせることを特徴とする寿命・余寿命評価プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011188390A JP2013050377A (ja) | 2011-08-31 | 2011-08-31 | 寿命・余寿命評価機能付き非破壊検査装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011188390A JP2013050377A (ja) | 2011-08-31 | 2011-08-31 | 寿命・余寿命評価機能付き非破壊検査装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2013050377A true JP2013050377A (ja) | 2013-03-14 |
Family
ID=48012524
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011188390A Withdrawn JP2013050377A (ja) | 2011-08-31 | 2011-08-31 | 寿命・余寿命評価機能付き非破壊検査装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2013050377A (ja) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107748198A (zh) * | 2017-09-19 | 2018-03-02 | 国家电网公司 | 基于感应视磁阻抗法的变电站接地网腐蚀诊断系统及方法 |
WO2018084011A1 (ja) * | 2016-11-04 | 2018-05-11 | 横河電機株式会社 | 減肉検出装置、減肉検出システム、減肉検出方法及びプログラム |
JP2018081070A (ja) * | 2016-11-04 | 2018-05-24 | 横河電機株式会社 | 減肉検出装置、減肉検出システム、減肉検出方法及びプログラム。 |
CN109238892A (zh) * | 2018-11-20 | 2019-01-18 | 东北大学 | 一种转子系统钢环式弹支结构强度设计及在线监测方法 |
CN109308556A (zh) * | 2017-07-28 | 2019-02-05 | 中国石油化工股份有限公司 | 一种基于内检测数据的管道评价方法及装置 |
CN109783850A (zh) * | 2018-12-08 | 2019-05-21 | 西安电子科技大学 | 高加速应力筛选试验的剩余寿命评估与可靠性分析方法 |
CN110147643A (zh) * | 2019-06-12 | 2019-08-20 | 中国神华能源股份有限公司 | 车钩钩体剩余寿命确定方法和装置 |
CN110930055A (zh) * | 2019-12-03 | 2020-03-27 | 上海核工程研究设计院有限公司 | 一种含缺陷管道破损安全期的评估系统及方法 |
CN112765840A (zh) * | 2020-12-29 | 2021-05-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种快速评估表面划痕对结构裂纹扩展寿命影响的方法 |
-
2011
- 2011-08-31 JP JP2011188390A patent/JP2013050377A/ja not_active Withdrawn
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018084011A1 (ja) * | 2016-11-04 | 2018-05-11 | 横河電機株式会社 | 減肉検出装置、減肉検出システム、減肉検出方法及びプログラム |
JP2018081070A (ja) * | 2016-11-04 | 2018-05-24 | 横河電機株式会社 | 減肉検出装置、減肉検出システム、減肉検出方法及びプログラム。 |
CN109308556B (zh) * | 2017-07-28 | 2023-03-24 | 中国石油化工股份有限公司 | 一种基于内检测数据的管道评价方法及装置 |
CN109308556A (zh) * | 2017-07-28 | 2019-02-05 | 中国石油化工股份有限公司 | 一种基于内检测数据的管道评价方法及装置 |
CN107748198A (zh) * | 2017-09-19 | 2018-03-02 | 国家电网公司 | 基于感应视磁阻抗法的变电站接地网腐蚀诊断系统及方法 |
CN109238892A (zh) * | 2018-11-20 | 2019-01-18 | 东北大学 | 一种转子系统钢环式弹支结构强度设计及在线监测方法 |
CN109783850A (zh) * | 2018-12-08 | 2019-05-21 | 西安电子科技大学 | 高加速应力筛选试验的剩余寿命评估与可靠性分析方法 |
CN109783850B (zh) * | 2018-12-08 | 2023-04-07 | 西安电子科技大学 | 高加速应力筛选试验的剩余寿命评估与可靠性分析方法 |
CN110147643A (zh) * | 2019-06-12 | 2019-08-20 | 中国神华能源股份有限公司 | 车钩钩体剩余寿命确定方法和装置 |
CN110147643B (zh) * | 2019-06-12 | 2023-11-14 | 中国神华能源股份有限公司 | 车钩钩体剩余寿命确定方法和装置 |
CN110930055A (zh) * | 2019-12-03 | 2020-03-27 | 上海核工程研究设计院有限公司 | 一种含缺陷管道破损安全期的评估系统及方法 |
CN112765840B (zh) * | 2020-12-29 | 2022-10-11 | 中国航空工业集团公司西安飞机设计研究所 | 一种快速评估表面划痕对结构裂纹扩展寿命影响的方法 |
CN112765840A (zh) * | 2020-12-29 | 2021-05-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种快速评估表面划痕对结构裂纹扩展寿命影响的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2013050377A (ja) | 寿命・余寿命評価機能付き非破壊検査装置 | |
Zenzen et al. | A modified transmissibility indicator and Artificial Neural Network for damage identification and quantification in laminated composite structures | |
Nanda et al. | Vibration based structural damage detection technique using particle swarm optimization with incremental swarm size | |
Michopoulos et al. | Modeling and simulation of multiphysics systems | |
Zhang et al. | Multivariable wavelet finite element-based vibration model for quantitative crack identification by using particle swarm optimization | |
Holl et al. | 3D multiscale crack propagation using the XFEM applied to a gas turbine blade | |
Pereira et al. | Three-dimensional crack growth with hp-generalized finite element and face offsetting methods | |
Nanda et al. | Modal parameter based inverse approach for structural joint damage assessment using unified particle swarm optimization | |
RU2554676C2 (ru) | Способ проектирования подводного оборудования, подверженного вызванному водородом растрескиванию под напряжением | |
Khatir et al. | Crack identification using eXtended IsoGeometric analysis and particle swarm optimization | |
Qinami et al. | Variational eigenerosion for rate‐dependent plasticity in concrete modeling at small strain | |
Katam et al. | A review on structural health monitoring: past to present | |
Kourehli | Damage diagnosis of structures using modal data and static response | |
Dincal et al. | Nondestructive damage detection in Euler–Bernoulli beams using nodal curvatures—Part I: Theory and numerical verification | |
Nanda et al. | Damage assessment from curvature mode shape using unified particle swarm optimization | |
JP6982847B2 (ja) | 破壊現象の解析装置、破壊現象の解析方法、及びプログラム | |
JP2019144765A (ja) | 柱状ハニカム構造体の最大引張応力推定方法及び最大引張応力推定装置 | |
Lee et al. | Application of XFEM to model stationary crack and crack propagation for pressure containing subsea equipment | |
Chadaram et al. | Identification of cracks length by XFEM and machine learning algorithm | |
Okereke et al. | Computational mechanics and the finite element method | |
Krishnamurthy et al. | Damage characterization using the extended finite element method for structural health management | |
Tian et al. | Fracture and fatigue analyses of cracked structures using the iterative method | |
Nanda et al. | Element damage assessment in semi rigid connected structures using modal domain data | |
Meng et al. | Prediction of optimal inspection time for structural fatigue life | |
JP2004171190A (ja) | ひび割れ表示方法、プログラム、記録媒体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A300 | Withdrawal of application because of no request for examination |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20141104 |