学際大規模情報基盤共同利用・共同研究拠点 平成 29 年度共同研究 最終報告書 2018 年 5 月

jh170045-NAJ

非均質・異方性材料中を伝搬する弾性波動解析手法の開発と非破壊検査への応用
斎藤隆泰（群馬大学大学院理工学府・環境創生部門）
概要 近年，超音波を用いた非破壊検査法に注目が集まっている．一般に，超音波は固体中で弾性波の性質を示す．そのため，超音波非破壊検査で要求されるシミュレーションは，欠陥や材料界面等による弾性波動散乱問題へと帰着される．しかしながら，超音波の波長は非常に小さいため，その解析は一般的に大規模なものとなる．本研究では，非均質性や異方性といった複雑な性質を示す材料中の波動伝搬を模擬するために必要となる弾性波動解析手法の開発，その高速化，並列化，ならびにそれら材料中の欠陥を特定するための逆問題に取り組んだ．数値解析手法には，境界要素法(BEM)，有限要素法(FEM)，有限積分法(EFIT)のいずれかを，扱う問題の性質毎に選定した。FRP 等の先進材料を対象に，開発した手法を適用することで，提案手法等の妥当性や有効性を検討した。

1. 共同研究に関する情報

FRP 等の異方性材料に対する逆散乱解析を実行(1) 共同研究を実施した拠点名する際，高精度な順解析結果が必要となる．共同

京都大学

研究者の古川は,異方性弾性波動問題に対する境(2) 共同研究分野

界要素法で学位を取得した若手研究者であり,異

超大規模数値計算系応用分野

方性弾性波動問題に詳しい．そのため，異方性を(3) 参加研究者の役割分担

示す FRP 材料中の弾性波動場を,演算子積分時間

平成 29 年度の参加研究者の役割分担は下記

領域境界要素法と呼ばれる最新の時間領域境界要の通りである．申請時の内容と概ね変更はな

素法を用いて解析する役割を担った．またい．

ACA(Adaptive Cross Approximation)を用いた境【斎藤隆泰・群馬大学・研究代表者】

界要素法の高速化や,逆散乱解析手法の構築にも研究代表者として，本研究を統括している．ま携った．

た，時間領域境界要素法の高速化や FRP 中の弾性

【紅露一寛・新潟大学】

波動散乱解析、並びに逆散乱解析を用いた欠陥形

共同研究者の紅露は,長年,境界要素法や有限要状再構成手法の開発，さらには Time-Reversal 法

素法等の数値解析手法の研究に取り組んで来た．の実施にも携わっている．また，演算子積分時間

本研究では,それらの知見を活かし，演算子積分時

領域境界要素法と有限要素法の結合解法の開発に

間領域境界要素法と有限要素法の結合解法の開発も携わった．

に関する研究や，境界要素法自体の高速化に関す【中畑和之・愛媛大学・副研究代表者】る研究にも携った．

研究代表者の斎藤と共に，副研究代表者として【牛島省・京都大学】

本研究を統括している．有限要素法や有限積分法

自由度の大きい大規模弾性波動解析を実施するを用いた大規模弾性波動解析手法の開発や，それには，並列計算は欠かせない．共同研究の拠点大らと境界要素法の結合解法の開発,Time-Reversal学（京都大学）側からの参加研究者として，解析法による逆伝搬解析にも携わった．

コードの並列化等に協力している．また，スーパ【古川陽・東京工業大学】

ーコンピューターのシステムに精通しているため，
1

学際大規模情報基盤共同利用・共同研究拠点 平成 29 年度共同研究 最終報告書 2018 年 5 月

計算実行のアドバイザー的立場の役割も担った．

で多数の繊維から成る非均質材料であり，マクロ【小山田耕二・京都大学】

レベルでは,その影響は音響異方性として波動伝同じく, 共同研究の拠点大学（京都大学）側か搬特性に現れる．すなわち,金属材料中の波動伝搬らの参加研究者として,弾性波動解析結果を可視特性とは全く異なる性質を示す．現状では, FRP に化する方法や，それらデータの効率的なポスト処対する超音波非破壊検査は，金属材料と同様のア理等に携わっている．また,京都大学の計算機を利プローチが適用されているため，検査精度に問題用する上でのアドバイザー的立場の役割も担った． がある．
【一色正晴・愛媛大学】

そこで，本研究では，材料の非均質性や音響異大規模弾性波動解析結果のデータは膨大なもの

方性を考慮した大規模弾性波動解析手法を開発し，となるため，その計算結果をどのように効率良く現状の超音波非破壊検査で特に問題となってい可視化するかについて検討が必要である．そのたる,FRP 等の異方性や非均質性の性質を有する材料め,様々な弾性波動解析結果のポスト処理に携わをターゲットとした実用的なシミュレーション技った．

術を確立することを目的として研究を行った．

【溝田裕久・愛媛大学大学院生】

FRP は炭素繊維による CFRP，ガラスを母材とし研究副代表者の研究室所属の大学院生として,た GFRP 等，様々な材料が開発されている．しかし有限要素法や有限積分法を用いた大規模波動解析ながら，これら材料の使用頻度は右肩上がりであや Time-Reversal 法を用いた逆解析等の解析コーること，金属のような等方性材料と全く異なる性ド作成補助に携わっている．またプリ・ポスト等質を示すことから考えても，金属材料に対する場の補助としての役割等も担った．

合と異なる新しい超音波非破壊検査法を確立しな【溝上尚弥・愛媛大学大学院生】

ければならない．FRP は，その一部は軍事目的への研究副代表者の研究室所属の大学院生であり，転用等も危惧され，輸出規制も定められている重数値解析のためのプリ・ポスト等の補助としての要な材料である．FRP の利用が進めば進む程，その役割等を担った．

ための維持管理手法の重要性は当然高まる．

【伊藤司・稲垣祐生・群馬大学大学院生】

このように，本研究では，異方性や非均質性と研究代表者の研究室所属の大学院生であり，数いった超音波非破壊検査を複雑にさせる 2 つの要値解析のためのプリ・ポスト等の補助としての役因を考慮できる，弾性波動解析手法の開発に焦点割等を担った．

を当てる．現在，政府の科学技術総合会議は科学技術イノベーションが取り組むべき課題に新しい2. 研究の目的と意義

インフラ点検・診断技術の開発に関するテーマを近年，インフラ構造物等の老朽化に伴い，構造掲げている．本研究は，最新の材料である FRP を材料等の健全性を評価する非破壊検査に注目が集もターゲットにいれた新しい超音波非破壊検査法まっている．特に，波動論を利用した検査法であの開発のための基礎を築くものであり，学術的新る超音波非破壊検査法は，今日まで広く利用され規性のみならず，社会的要請にも応じる時宜を得てきた．超音波は固体中で弾性波として伝搬するた研究である．

ことから，弾性波の特性を把握し，活用すること

3. 当拠点公募型共同研究として実施した意義

が超音波非破壊検査の高度化には重要である．
超音波非破壊検査は,均質な金属材料に対して本研究に参画している研究者は，それぞれ得意は多くの検査実績がある．しかしながら，例えば，な数値解析手法や特徴を有している．互いに非破

近年，航空宇宙や土木分野等で用いられつつある壊検査の研究をしていても，用いる数値解析手法

FRP(Fiber Reinforced Plastic)は,ミクロレベルが異なると，解析対象も無限領域や異方性を対象
2

学際大規模情報基盤共同利用・共同研究拠点 平成 29 年度共同研究 最終報告書 2018 年 5 月

図 2：積層構造を考慮した CFRP 中の超音波伝搬解析モデル．

図 1：円形空洞表面における(a)固体骨格部の変位の絶対値と(b)氷骨格部の変位の絶対値．とする等，異なる場合が多い．また，その可視化方法やデータの処理方法，大規模計算のための並列化もいわば，独自に行ってきた背景がある．さ
らに，研究代表者・副代表者共に，地方国立大学に在籍しているため，大型計算機を使う環境が,現状では本制度以外に存在しない．研究室単位での並列計算機の維持管理は多大な労力も必要となるばかりか,セキュリティの関係上,研究室単位で計算機をシェアし研究を進めることは難しいのが現

図 3：積層構造を考慮した CFRP 中の超音波伝状である．そのような中で，共通の計算機の利用搬解析結果(上段)シミュレーション結果(下は互いの数値解析手法を共有でき，かつ互いの研段)LDV 計測による可視化結果．

究課題遂行を促進するのに役立っている．実際,以下で述べる FRP 中の超音波伝搬解析も，本課題な研究に関連し，ここでは凍結した固体中の面外波しでは決して実施できなかったと考えられる．

動問題に対する演算子積分時間領域境界要素法にACA の実装を行った．図 1 は，提案手法を用いて4. 前年度までに得られた研究成果の概要

円形空洞に対する入射波の散乱問題を解いた場合

H29 年度新規課題のため,該当しない．の，円形空洞表面における(a)固体骨格部の変位(b)氷骨格部の変位の絶対値を示しており，図1(a)(b)中の実線は，解析解を示してある．いずれ

5. 今年度の研究成果の詳細の要素数𝑁𝑁𝑒𝑒 に対しても，提案手法による解析結果2.で述べた研究目的を達成するために，H29 年
度は，主に以下の内容について実施した．

は，解析解と一致しており，ただしく解析を実行（A-1）境界要素法による大規模弾性波動解析手法できていることが確かめられた．なお，この成果の開発（斎藤・紅露・古川）は，例えば 7 節の研究成果リスト(1)学術論文 f)で既に発表しており，計算効率を向上させた結果境界要素法を用いた弾性波動解析の効率化のた
め に ， 本 研 究 で は ACA(Adaptiveについても言及されている．

CrossApproximation)を適用する方法を試みた．具体的（A-2）有限要素法・有限積分法による大規模弾性には，ACA を演算子積分時間領域境界要素法と呼波動解析手法の開発（斎藤・中畑・溝田）
本課題(A-2)は，課題申込時点では，中畑・溝田

ばれる最新の時間領域境界要素法に適用した．本3

学際大規模情報基盤共同利用・共同研究拠点 平成 29 年度共同研究 最終報告書 2018 年 5 月

図 4：BEM-FEM 結合解法で求めた非均質領域周辺の弾性波動場．の担当であったが，研究代表者の斎藤も参画した．
超音波非破壊検査では，扱う超音波の波長は数 mmのオーダーなので，数 cm 角の試験体に対しても，解くべき問題は大規模となる．ここでは，Imagebased モデリングを施した有限積分法や有限要素法コードを開発した．FRP は，様々な積層角を持って成形されるため，与えられた FRP に対して数値

図 5：逆散乱解析結果(a)CFRP(b)等方性鋼解析を施すには，まず，弾性定数を推定する必要材の場合．

がある．そこで，FRP の弾性定数を推定する方法を

果の一部は，H29 年度安全・安心な社会を築く先進開発し，開発した方法を用いて得られた弾性定数材料・非破壊計測技術シンポジウムにて大学院生を先に開発した解析コードの入力として与えるこが発表し，新進賞を受賞している．とで，FRP 中の超音波の弾性波動解析を行った．図（A-3）境界要素法と有限積分法・有限要素法の結2 は本研究で行った解析モデルの一例であり，プ合による非破壊検査シミュレーション（斎藤・中リプレグの弾性定数を求め，対応する積層角毎に畑・紅露）弾性定数を定めている．それら各層をモデル化し一般的に，境界要素法は無限または半無限領域た場合の，CFRP 上面からの超音波伝搬のシミュレを含む波動解析に威力を発揮し，有限要素法は閉ーション結果を図 3 に示す．図 3（上段）はシミュ領域の解析に威力を発揮する．また，境界要素法

レーション結果を，図 3（下段）は LDV による実験は均質な波動場における散乱問題に強いが，非均結果を示している．解析に用いたパラメータ等の質領域の解析は苦手である．一方，有限要素法は，詳細は省略するが，シミュレーション結果は LDV非均質領域における波動解析を比較的得意としての実験結果と比較的良く一致している．特に，超いる．そこで，本研究では，従来法に比べて数値音波が右下方向へ伝搬する音響異方性の影響を良安定な最新の時間領域境界要素法である，演算子く再現できていることが見て取れる．よって，開積分時間領域境界要素法を有限要素法と結合させ発したコードや FRP の弾性定数推定法の妥当性がる新しい数値解析手法の開発を試みた．具体的に確かめられた．これらの結果は，例えば査読付学

は，局所的に非均質性を有する固体（非均質部分術雑誌に投稿され，既に採択されている(7 節研究はコンクリートを想定）に対する弾性波動解析手成果リスト(1)学術論文 d),g))．また，本研究の成法を開発した．なお，非均質なコンクリート部分
4

学際大規模情報基盤共同利用・共同研究拠点 平成 29 年度共同研究 最終報告書 2018 年 5 月にはボクセル有限要素法を採用し，コンクリートのモデル化には X 線 CT スキャン画像を用いた．また，これら非均質領域における散乱波の計算には，演算子積分時間領域境界要素法を用いた．図 4 は，本研究で行った解析結果の一例を示している（例えば７節研究成果リスト(1)学術論文 c)）．図 4 中の赤点線は，有限要素法と境界要素法の結合境界を示している．図 4 より，水平方向に伝搬する入射平面波は，骨材等に散乱されながら伝搬していることがわかる．また，非均質領域により散乱された超音波が，無限遠方に伝搬している様子を見て取ることができる．すなわち，開領域を境界要素法で解析しているため，無限遠方を容易に扱う事が可能となった．
（B-1）非均質・異方性材料中の欠陥に対する逆散図 6：時間反転法による欠陥の同定（上段）

乱解析（斎藤・古川）

解析モデル（下段）解析結果．

本研究では，まず(A-1)で開発している演算子積分時間領域境界要素法を 2 次元面内異方性弾性波のような強異方性材料に対しても，欠陥形状を再動問題へ拡張した．解析対象として，現在，原子構成できていることがわかる．なお，図 5(a)に示

力機器等に用いられているオーステナイト系鋼材

す CFRP の場合の結果は，真円の空洞がやや横長におよび CFRP を選定した．これらはいずれもミクロ再構成されている．これは，本解析で取り扱ったな支点から見れば非均質な性質を示しているが，CFRP は，一方向繊維強化 CFRP であり，水平方向超音波の波長によっては均一な異方性と扱って問の弾性波速度が極めて速いためであり，散乱波の題ないであろう．そこで，いわば均質化された弾受信方向によって再構成に用いるデータ数が異な性定数を用いて，これら材料中の異方性弾性波動る事が原因の一つとして考えられる．なお，この解析を実施した．現状ではこれら材料に対する定成果は，既に査読付き論文に投稿している（7 節研量的な超音波非破壊検査法は確立されていない．

究成果リスト(1)学術論文 a））．

まず，オーステナイト系鋼材，CFRP 材料を想定し，（B-2）Time-Reversal 法に基づく逆伝搬解析（斎入射平面波に対する欠陥の散乱問題を解析した．

藤・中畑・溝田）次に，これら異方性材料中の欠陥に対する逆散乱

課題申請時点では，Time-Reversal 法は中畑・溝解析の定式化を行った．定式化は，一般の異方性田の分担であったが，H29 年度開始前に，本研究の材料に適用できるよう工夫した．逆散乱解析にお前段階として行った斎藤の Time-Reversal 法に関ける欠陥形状再構成には，欠陥からの散乱波形がする研究が第 44 回土木学会関東支部で講演賞を必要となる．ここでは，先の解析で求めた平面波受賞（発表は大学院生）したため，本研究には研を想定した場合の数値シミュレーションによる散究代表者の斎藤も参画している．図 6 に Time-乱模擬波形を用いて欠陥形状を再構成することをreversal 法を用いて等方弾性体内部の 4 つの空洞行った．図 5(a)に CFRP 中の欠陥に対する逆散乱を再構成した結果を示す．解析では図 6（上段）で解析結果を示す．比較のため等方性材料に対する示すように，64 素子から成るリニアアレイ探触子結果も図 5(b)に示してある．CFRP は繊維方向に強を用いた場合の超音波非破壊検査を想定したシミ

い異方性の性質を示すことで知られているが，そュレーションを実施している．図 6（下段）は，図5

学際大規模情報基盤共同利用・共同研究拠点 平成 29 年度共同研究 最終報告書 2018 年 5 月

6（上段）の点線内の領域について解析結果を可視並列化を施した大規模弾性波動解析手法を開発化した結果を示しており，リニアアレイ探触子のし，FRP に適用することで一定の成果を得ている．

各素子からの超音波発信時間を制御することで，しかしながら，GPU を活かした超大規模計算や最斜め 30°に超音波を送信している．図 6（下段）

新の XeonPhi 等を活かした高速計算には至っていより，4 つの空洞の位置を概ね特定できているこない．引き続き開発したコードの改良を行いながとが見て取れる．

ら，これらハードウェアを利用した超大規模弾性（C）波動伝搬解析の並列化及びポスト処理の効率波動解析への応用課題が残されていると考えられ化（牛島・小山田・一色）る．

上 記 (A-1)-(B-2) に 至 る 数 値 解 析 コ ー ド の（A-3）境界要素法と有限積分法・有限要素法の結OpenMP, MPI 並列化や，その計算結果に対するポ合による非破壊検査シミュレーション（斎藤・中スト処理等で連携し，解析自体や可視化技術に関畑・紅露）
時間領域境界要素法と有限要素法の結合解法のする効率化を行った．

開発には成功しているため，一定の成果を残せて6. 今年度の進捗状況と今後の展望いると言える．国内学会発表も行っているため，今年度の進捗状況は概ね良好であったと考えら

次回のステップとして解析結果をまとめ，海外のれる．5 節で述べた(A-1)-(B-2)はいずれも課題申査読付論文に投稿することを検討している．ただ込時点での予定研究内容であり，いずれの項目もし，解析は，現状では 2 次元弾性波動問題に留ま中間報告の時点で既にいくつかの結果を導くことっている．そのため，3 次元問題への拡張も残されができており，予定通り研究を進めることができた課題である．特に，3 次元問題を扱う場合は，問た．今後は，それぞれの研究内容に対して，下記題の自由度が大きいため，計算負荷の大きい境界の点に留意して研究を進める必要がある．

要素法の解析部分には，高速多重極法や ACA 等を（A-1）境界要素法による大規模弾性波動解析手法取り入れて計算を高速化する検討も必要であるとの開発（斎藤・紅露・古川）考えられる．一方で，有限積分法と時間領域境界H29 年度は境界要素法を高速化するための ACA

要素法の結合解法の開発については，引き続き検コードの開発を行い，いくつかのパラメータ解析討を行っていく予定である．

や検証を行った上で査読付論文にも投稿している． （B-1）非均質・異方性材料中の欠陥に対する逆散乱解析（斎藤・古川）

しかしながら，解析対象とした波動場は，2 次元であること，弾性波動問題への適用や異方性弾性波

2 次元面外波動問題，面内波動問題を扱った場動問題への適用等，今後適用すべき問題は数多い．

合の異方性材料中の空洞欠陥に対する逆散乱解析

そのため，今後は解析対象範囲を広げ，それぞれは，既に定式化を行い，一定の成果を示したこと定式化を行い，解析を行っていく必要がある．特で，査読付論文にもアクセプトされた．しかしなに，異方性弾性波動問題の高速化は，今後の大きがら，一般的に CFRP は層間剥離が問題となるため，な課題である．高速多重極法や ACA を適用するだ欠陥の形状はき裂に近いものが対象となる．逆散けでなく，例えばプログラム高度化共同研究課題

乱解析は，積分方程式を基礎とした定式化である

プログラム等への公募申請についても検討し，解が，欠陥の対象がき裂であるか，空洞であるかに析コードを高度化することも検討している．将来

よって，定式化が異なる．そのため，今後は，異方的には MPI や OpenMP の並列化のみならず，GPU を性材料中のき裂を対象とした欠陥形状再構成手法用いた高速化についても検討が必要であろう．

を開発することが必要となると考えられる．また，（A-2）有限要素法・有限積分法による大規模弾性3 次元問題への拡張も今後，取組むべき課題であ波動解析手法の開発（斎藤・中畑・溝田）る．そのためには，3 次元異方性弾性波動問題に対6

学際大規模情報基盤共同利用・共同研究拠点 平成 29 年度共同研究 最終報告書 2018 年 5 月

する基本解の遠方場近似を求める等，解決すべきc)森川光・藤縄和宏・市川諒・斎藤隆泰：2 次元弾

いくつかの問題が残されているので，これら課題

性波動問題に対する演算子積分時間領域境界要素を 1 つずつ解決していく予定である．

法・イメージベース有限要素法結合解法の開発,（B-2）Time-Reversal 法に基づく逆伝搬解析（斎第 22 回計算工学講演会論文集, Vol.22,CD-ROM 収藤・中畑・溝田）録,(2017)

2 次元異方性弾性波動問題を対象に成果を出し，d) 溝上尚弥・中畑和之・黄木景二・堤三佳・森亜

いくつかの口頭発表を行っている．しかしながら，也華・斎藤隆泰：レーザースキャンによる音響異時間反転法は，受信波形を解析対象の波動場に再方性を有する CFRP 中の超音波の可視化と弾性ス入射させて，如何にして定量的に欠陥位置を示す

ティフネスの推定,土木学会論文集 A2(応用力学),かに課題が残されている．すなわち，時間反転さ

vol.73,pp.I-343-I354,(2017)

れた波動は時間領域であるため，再入射された複e) T. Maruyama, T. Saitoh and S. Hirose：Numerical数の波動が収束する箇所を時間と共に追跡する必

study on sub-harmonic generation due to interior and要があるが，現状では，その追跡方法は，必ずしsurface breaking cracks with contact boundaryも定量化されていない．今後は，例えば時間反転conditions using time-domain boundary element法にトポロジー感度の概念を取り入れて，欠陥位method, International Journal of Solids and Structures,置を定量的に推定する方法等を検討している．まvol.126-127, pp.74-89, (2017)

た，2 次元異方性弾性波動問題へと拡張し，査読付f) 古川陽・斎藤隆泰・廣瀬壮一：周波数領域境界論文に投稿できるよう，解析結果をまとめていく要素法と ACA による凍結した多孔質体内部の介在必要もあるだろう．本研究では，超音波非破壊評

物による波動散乱解析，計算数理工学論文集,価法を想定し，時間反転法の利用を検討したが，vol.17,pp.83-88,(2017)

時間反転法を用いて AE(アコースティックエミッg) T. Saitoh and A. Mori: Development of new elasticション)法等，弾性波計測の類似分野における時間

constant estimation method using the laser ultrasonic

反転法の応用を行うことも考えている．

visualization（C）波動伝搬解析の並列化及びポスト処理の効率

Environmental Engineering, vol.1, pp. 669-674, (2017),

化（牛島・小山田・一色）

DOI: 10.1201/9781315116259-118

testing,

Civil,Architectureand(A-1)-(B-2)に渡る関連項目について，引き続き， (2) 国際会議プロシーディングス
並列化およびポスト処理の効率化について検討すa)Numerical simulation of ultrasonic wave propagation

る予定である．特に，Xeon Phi や GPU 等の利用とin fiber reinforced plastic using image-based modeling,

超大規模問題への応用は今後，検討すべき事項でN. Mizokami, K. Nakahata, K. Ogi, H. Yamawaki, M.

ある．

Shiwa, AIP Conference Proceedings, Vol.1806, (2017)
pp.150005, DOI: 10.1063/1.4974729

7.

研究成果リストb) K. Nakahata, K. Ogi, T. Namita, K. Ohira, M.

(1) 学術論文

Maruyama, T. Shiina: Photoacoustic microscopica)稲垣祐生・斎藤隆泰・古川陽・廣瀬壮一：一方向

imaging of surface and subsurface damages in CFRP,

炭素繊維強化 CFRP 中の欠陥に対する逆散乱解析，

Review of Progress in Quantitative Nondestructive

計算数理工学論文集,vol.17,pp.13-18,(2017)

Testing, AIP Conference Proceedings, now printing.

b)An ultrasonic testing method for homogeneousanisotropic materials, H. Mizota, Y. Nagashima, K.

(3) 国際会議発表

Nakahata, Insight, Vol.59, No.7,pp.351-357,(2017),

a) T. Saitoh, Y. Inagaki, A. Furukawa and S. Hirose: 2-DOI:10.1784/insi.2017.59.7.351.

D inverse scattering analysis for a defect in authentic
7

学際大規模情報基盤共同利用・共同研究拠点 平成 29 年度共同研究 最終報告書 2018 年 5 月

stainless steels, The 8th International conference on術シンポジウム,愛媛大学,pp.77-78,(2018)

computational methods (ICCM2017), ID2747, (2017)j) 前原佑・斎藤隆泰・今井済・佐藤明良:有限要素法を用いた L 字型 CFRP 中の三次元超音波伝搬(4) 国内会議発表

シミュレーション,第 25 回超音波による非破壊評a)溝上尚弥・中畑和之・斎藤隆泰：レーザー超音

価 シ ン ポ ジ ウ ム , 東 京 都立 産 業 技 術 セ ン タ ー ,波法による超音波伝搬の可視化と弾性定数推定のpp.133-134, (2018)

試 み , 平 成 29 年 度 土 木 学 会 全 国 大 会 , 九 州 大k)中畑和之：FRP を伝搬する超音波のモデル化と学,CD-ROM 収録,番号 I-376,pp.751-752,(2017)

表層きずの光音響イメージング,H29 年度安全・安b)稲垣祐生・斎藤隆泰:オーステナイト系鋼材中の心な社会を築く先進材料・非破壊計測技術シンポ

欠陥に対する 2 次元順解析および逆散乱解析,平ジウム,pp,i-ⅷ,(2018)

成 29 年度土木学会全国大会,九州大学,CD-ROM 収l) 森川光・斎藤隆泰：トポロジー感度を用いた 2録,番号 I-119,pp.237-238,(2017)

次元等方性材料中の欠陥形状再構成に関する研究,c)大芦健太・斎藤隆泰:時間反転法を用いた FRP 中

第 45 回土木学会関東支部技術発表会 講演概要の欠陥形状再構成に関する研究,平成 29 年度土木集,CD-ROM 収録, I-34, ,山梨大学

学会全 国大会,九 州大学,CD-ROM 収 録,番 号 I126,pp.251-252,(2017)(5) その他（特許，プレス発表，著書等）

d)中畑和之：革新的な非破壊検査技術開発への取a)中畑和之・小澤耀生・大平克己・小川健三：低周り組み,日本保全学会 第 14 回学術講演会要旨集， 波アレイ探触子を用いたコンクリート深部の埋設特別企画, pp.5-12, (2017)

物の高速映像化, 検査技術, 第 22 巻，第 10 号,e) 斎藤隆泰・今井済・佐藤明良：有限要素法を用

pp.6-10, (2017)

いた L 字型 CFRP 中の超音波伝搬シミュレーショb)中畑 和之・黄木 景二・斎藤 隆泰:炭素繊維強ン,平成 29 年度非破壊検査協会秋季講演大会,福化樹脂中を伝搬する超音波の可視化計測とシミュ

岡国際会議場,pp.9-10,(2017)

レーション, 超音波 TECHNO, 第 29 巻, 第 4 号,f) 溝上尚弥・中畑和之・黄木景二・堤三佳・森亜

pp.69-74, (2017)

也華・斎藤隆泰：レーザー超音波法を用いた CFRPc)斎藤隆泰：逆散乱解析法を用いたオーステナイ中の波動場の可視化と数理モデルの構築,第 64 回

ト系鋼材中の欠陥形状再構成,超音波 TECHNO,第理論応用力学講演会講演論文集,機械振興会

30 巻, No.2, pp.85-89, (2018)

館,USB 収録,(2017)g) 大芦健太・斎藤隆泰：時間反転法を用いたオーステナイト系鋼材中の欠陥形状再構成に関する研究,第 64 回理論応用力学講演会講演論文集,機械振興会館,USB 収録,(2017)h) 稲垣祐生・斎藤隆泰：逆散乱解析によるオーステナイト系鋼材中の欠陥形状再構成,第 64 回理論
応用力学講演会講演論文集,機械振興会館,USB 収録,(2017)i) 大芦健太・斎藤隆泰・中畑和之：レーザー超音波可視化試験を用いた擬似等方積層板に対する弾性定数の推定と FDTD 法による検証, 平成 29 年度
安全・安心な社会を築く先進材料・非破壊計測技8