Non-destructive evaluation of 3D plastic strain and residual stress distributions induced by earthquakes using non-linear inverse analysis

2022 Fiscal Year Research-status Report

• Project/Area Number: 22K03831
• Research Institution: Kogakuin University
• Principal Investigator: 小川 雅 工学院大学, 工学部, 准教授 (90635236)
• Project Period (FY): 2022-04-01 – 2025-03-31
• Keywords: 塑性ひずみ / 変位 / 非破壊 / 逆問題 / 非線形

Outline of Annual Research Achievements

地震発生後、発電プラントの迅速な再稼働の可否の判断が求められる。本研究では、デジタルカメラを用いた画像相関法により非破壊計測した変位とX線回折法による表面ひずみから、それらの原因である３次元の塑性ひずみを逆問題解析により推定できるようにすることを目指す。ここで、地震荷重により生じる塑性変形は比較的大きくなることがあるため、その場合には発生する塑性ひずみと変位との関係が非線形になる。
２０２２年度には、まず１回の外力によって塑性変形が発生する場合、その結果として発生する変位から逆問題解析により原因である塑性ひずみを正しく推定する方法について検証した。ここで、塑性ひずみと変位との関係が線形になる場合と非線形になる場合の２つについて、単純な３点曲げ解析により得られた塑性ひずみ分布を正解とし、その結果として生じる変位から逆問題解析を試みた。その結果、線形問題の場合には従来の一般的な逆解析手法により正しく推定できることを確認した。また、非線形の場合にも、提案する繰返し計算手法により正解の値へと収束した。
次に、２回の外力の発生により塑性変形が生じる場合について検証する際、ソフトウェアの問題から、せん断成分の塑性ひずみを推定することができないため、本数値解析では、せん断成分の塑性ひずみが生じないように材料特性を設定した。ただし、この結果により生じた変位から塑性ひずみを逆解析する際、数値解析プログラムの問題で塑性ひずみと変位との関係行列を正しく得ることができていないことが明らかとなった。２０２３年度は、まずこのプログラムの問題を修正する必要がある。

Current Status of Research Progress

2: Research has progressed on the whole more than it was originally planned.

Reason

２０２２年度には、原因である塑性ひずみとその結果として生じる変位との関係が非線形となる場合にも、本提案手法により、原理的に部材表面の２次元変位から正しく３次元の塑性ひずみ分布を推定できることを数値解析により示すことができた。これにより、１回の外力によって塑性変形が発生する場合の塑性ひずみの推定方法を確立することができた。研究を進めていく中で、利用している汎用の有限要素解析ソフトではせん断成分の塑性ひずみを考慮することができないことがわかったが、せん断成分を無視できるように、材料特性の設定を行うことができた。まだ、複数回の塑性変形が生じる場合について、本手法の有効性を検証することができていないが、検証を行うための数値解析プログラムを概ね設定することができている。引き続き、数値解析により、手法の原理的な有効性の評価を実施する。

Strategy for Future Research Activity

２０２２年度には、概ね本提案手法の原理的な検証を行うことができたものの、複数回の塑性ひずみが生じる場合について、検証することが課題となっている。
２０２３年度には、引き続き数値解析プログラムを修正し、２回の塑性ひずみが生じる場合に、非破壊に取得可能な表面変位と表面弾性ひずみの値から、逆問題解析により、正しく塑性ひずみ分布を評価できるかどうかを検証する。
次に、実際の構造物に対する適用に向けて、まずは数値解析により計測誤差の影響について考察する。変位や弾性ひずみにおける現実的な計測誤差を想定した場合の本手法の推定精度を評価し、比較的推定精度が低い場合には、本手法の推定精度を向上させる方法について検討する必要がある。
実際の構造物に対する本提案手法の適用を試みる際には、１回の塑性変形が生じた場合と複数回の塑性変形が発生する場合について、本手法の推定精度を評価する。ここで、正解の塑性ひずみ分布を知るために、外力の入力履歴を数値解析上で再現し、その際に得られる塑性ひずみ分布を正解として、本手法による残留応力の推定値との比較を行う。変位の計測は、デジタルカメラによる画像相関法により非破壊に計測する場合も考えられるが、計測誤差の大きさによっては、レーザー変位計を用いる場合も考えられる。なお、構想物にはひずみゲージを貼り付け、意図した変形を再現できているかを確認する予定である。

Causes of Carryover

当初、２０２２年度は比較的性能のよいデスクトップパソコンを計上し、また有限要素解析ソフトのライセンスやメッシュ分割ソフトのソフトウェアを計上した。ただ、予算申請時より全体の予算が縮小された関係で、当該年度は有限要素解析ソフトのグレードを下げ、本学で利用可能なソフトを利用することとした。また、当該年度は手法の原理検証を行うため、それほど計算時間がかからない。原理検証の後、モデルを発展させる際には高性能な計算機が必要であるものの、現時点では既存の計算機を利用することとした。
また、次年度には実験計測において多くの費用を要することになる。デジタル画像相関法による変位計測を実施するだけでなく、表面残留応力の計測も必要となる。残留応力の測定には本研究室にある可搬型のＸ線回折装置を利用するが、経年劣化のため、X線の管球を交換しなければならない可能性があることもわかった。さらには、国際会議や各種学会も予定されているため、次年度に予算を盛り込む必要があると考え、次年度使用額を生じさせることとなった。