⼟⽊学会全国⼤会 研究討論会（⽇本⼤学津⽥沼，4/Sep/2013)

愛媛⼤学
中畑 和之

nakahata@cee.ehime-u.ac.jp

１．現状認識と問題(センシングと⾮破壊検査の現場から)
１．応⽤⼒学の理論・解析技術は多くの研究実績があるが，実際の現場において必ずしもこれらのアプローチが⼗分活⽤されているとは⾔えない．
２．計算科学や解析データと，思考（解釈）を繋げたメインテナンス技術はまだまだ不⼗分

応⽤⼒学とIT技術の融合によって，新たな⽅法論が創出でき，推定（評価）精度は⾶躍的に伸びる可能性がある2

情報爆発時代における⼟⽊⼯学
データ指向型科学，ビッグデータ，ナレッジマネジメント

推定精度

物理・化学モデリング，最適化・逆解析理論（応⽤⼒学が得意とするところ）

詳細・⼤規模解析

加速？

減速？

現在
モデルベース，1D-CAE，ルール抽出，パターン分類，マイニング，可視化（ITが得意とするところ）

現場への適⽤度3

現場を意識したアプローチ
１．時間軸はともかく，どういう形であれ，結果的に応⽤⼒学をメインテナンスの現場の役に⽴つ学問としたい．
２．⼤量のセンシングデータを計測・記録できるようになった現在，それを安全管理に活かすには，集約化，⾼速処理，情報伝達(共有)が重要なキーワード

現場で使える応⽤⼒学を⽬指して，メインテナンス（⾮破壊検査とセンシング）とITとの融合を模索

4

２．⾮破壊検査へのアプローチ

・HPC技術を探傷の現場へ
・計測データとシミュレーション技術を融合した逆解析

5

２．⾮破壊検査へのアプローチ
⾮破壊検査の種類と近年のニーズ

超⾳波法

電磁波法

波動（振動，⾳波，弾性波，電磁波，X線，⾚外線），磁⼒，電⼒等の物理エネルギーを利⽤．

超⾳波探傷法 (三菱重⼯（株）HP)

電磁波レーダ法(構造総合技術研究所（株）HP)

複雑な3次元形状

配管⾦型

複合材料・異⽅性材料

コンクリート

繊維補強材料(FRP)

配管溶接部

イメージベース波動伝搬シミュレータ
数値シミュレーションをもっと現場で⼿軽に利⽤したい！
・検査対象を正確にモデル化したい．
でも，プリ作業が⾯倒なのは避けたい．
イメージベースモデリング・⼯学的に満⾜できる精度で計算結果がほしい．
有限積分法(FIT)・計算時間が⻑いのはナンセンス
GPUによる並列計算

a.イメージベースモデリング
被検体の写真やCTデータ等のデジタルイメージから数値モデルを作成し，シミュレーションを実⾏する⼿法．
デジカメの写真からClose-up

画像処理
Pixel data
EBSPデータから

表⾯計測データから(光計測)zoom

ボクセル要素積層写真から(X線CT)
zoom

様々なデジタルデータからボクセルデータ(3D)の数値モデルを作成．

b. 有限積分法 FIT*(Finite Integration Technique)・波動伝搬解析法の⼀⼿法．
・構成則と運動⽅程式（ファラデー則とアンペール則）を時間域と空間域で積分し，その離散式を時間ステップごとに逐次計算することで陽的に解を求める⽅法．
⾳波版：A-FIT（空気or液体を伝搬する超⾳波）
弾性波版：E-FIT （固体中を伝搬する超⾳波）
電磁波版：EM-FIT（⾮⾦属を伝搬する電磁波）＜特徴＞
・境界条件の取り扱いが理論的であり，⾮均質材料や異⽅性材料のモデルを計算できる．
・積分セル形状を⼀定とすればステンシル計算となり，陽解法のため並列計算の効率が出やすい．
＜これまでの研究経緯＞
・Validation and verification は，様々な対象で実施済み．
*K. Nakahata, K. Terada, T. Kyoya, M. Tsukino and K. Ishii, Simulation of ultrasonic andelectromagnetic wave propagation for nondestructive testing of concrete using image-based FIT,
Journal of Computational Science and Technology, Vol.6, No.1, pp.28-37, 2012.

計算例

1m
0.06m

0.5m
0.1m
コンクリートコア中の超⾳波伝搬(E-FIT，Flat MPI, 約0.5億ボクセル)

コンクリート床板中の電磁波伝搬(EM-FIT, Flat MPI, 約3.5億ボクセル)

c.

GPUによる並列計算
SMX:ストリーミングプロセッサエクストリーム

SP:ストリーミング
プロセッサ
SMX1つあたり192個シェアードメモリ(超⾼速)

ビデオメモリ（グローバルメモリ）低速

レジスタ・メモリ(超⾼速)

NVIDIA製GPUの構造(Kepler)

GPC:グラフィックプロセッシング
クラスタ

GeForce GTX680 (GPU)
コア数
プロセッサ性能

メモリ

1536(192x2x4)

Core i7 965 (CPU)4クロック (MHz)

1006

3200

ピーク演算性能

4121 Glops

51.2Gflops

標準メモリ設定

2GB GDDR5

各⾃PCによる

192.2625.6

メモリバンド幅 (GB/sec)

MPI-OpenMPハイブリッド↓
MPI-CUDA ハイブリッド（HPC技術の利⽤！）
GPUメモリで計算
start

ボクセル
データ
⼊⼒

(FITの体積積分１セルが1スレッドに対応）

速度場の計算(電場)

応⼒場の計算
(磁場)
陽的に更新

袖部通信・
時刻歴応答の出⼒

endIT企業・⾮破壊検査会社と共同で「シミュレータ SWAN21」を開発
GPUによる超⾼速計算(愛媛⼤開発のGPU化したFITのコード)

PCの空きスロットに差し込むだけで使⽤開始．ノートPCで駆動可能．

⾼機能GUIの設計・開発(ジャパンプローブ(株)，eCompute(株)と共同開発)

⾼速シミュレータを利⽤した新しい逆解析の提案？
被検体内部が⾮均質，異⽅性材料

路程の推定が困難（速度分布が不均⼀）
時間反転法*?
界⾯形状が複雑（凹凸，ギザギザも）逆解析（⽋陥の再構成）
が困難
SWAN-21• K. Nakahata and K. Kimoto, Reconstruction of flaws in heterogeneous media using image-based FIT• and time reversal approach, Journal of Solid Mechanics and Materials Engineering, Vol.6, No.6,• pp.771-781, 2012.

時間反転法の原理
sourceアレイ探触⼦計測されたエコーを時間反転して“⼊射”線形波動の時間に対する可逆性を仮定すれば，Sourceに集束する（定式化済み）

実験による検証

電⼦スキャン装置
・⽇⽴エンジニアリングアンドサービス製

アレイ探触⼦・ジャパンプローブ製
・中⼼周波数： 3MHz(400KHz)
・N=64 (24)・1素⼦毎のパルス波の送受信が可能

・P=0.8mm(5mm)

くさび・ジャパンプローブ製
・ポリスチレン

検証実験１
アルミニウム中の⽋陥（スリット）の映像化
送信：素⼦No.1〜24の合成波
アレイ探触⼦（72素⼦）

受信：No.25〜72素⼦（計48個）

48個の受信波を時間反転し，「SWAN21」に⼊⼒する．

くさびスリット(4mm)

SWAN21上に作成した数値モデル(スリットの位置は未知)

Re-emission oftime-reversed waves

1

2

3

エコーはもともとスリットから来たものであるから．スリットに集束するはずである．
可視化プロセスで，集束点(最も振幅が⼤きかったとき)が，実際の空洞の位置をあらわす．

Mises応⼒で可視化

３．センシングへのアプローチ

・ITプラットフォームを⽤いたソフトウエア設計
・計測データと可視化技術を融合した構造物の特性評価

26

３．センシングへのアプローチ
振動の多点計測および可視化システム（3軸加速度センサ，MEMS）
データ転送ノード

センサノード（⾃作）
センサノード
加速度計測
多点でセンサ配置
データ転送ノード
A/D変換，ワイヤレス通信
基地局PC
加速度→変位（数値積分）
3次元可視化，
データ処理

基地局PCへ転送
ワイヤレス通信ナショナルインスツルメンツ（NI）社の開発プラットフォームを利⽤して⼀元管理（データ取得→転送→可視化）

歩道橋で実験松⼭市某所の歩道橋

実験条件
加振：⼈の跳躍による衝撃
サンプリングレート：1.0kHz
バンドパスフィルタ(2 to 400Hz)
データ計測時間：50 sec.

16点で得られた変形データをCAD上にマッピング．
衝撃による橋梁床版振動の3次元リアルタイム可視化

しゃがむ1.48 [s]，跳躍開始2.08 [s]，着地後は⾃由振動

3次元可視化から分かること
z7番ノードの周波数解析

1357911131546810121416y2Frequency[Hz]Power spectrumx

ウエーブレット変換

3.9Hz
3.9Hz

Frequency[Hz]24.7HzTime[sec]

3.9Hz：1次の曲げモード
24.7Hz : 曲げモードなのか，ねじりモードか？，何次モードなのか？

[dB]

狭帯域のバンドパスフィルタを作⽤して，再び3次元可視化(24.07Hz – 25.43Hz）

Power spectrum

卓越周波数の把握

24.7Hzは何？

3D visualization

曲げ
1次モード
Frequency[Hz]これは3次の曲げモードであることが分かる（2次では無い）

松⼭市某所の歩道橋

実験条件
加振：⾼覧を2Hzで揺する
サンプリングレート：1.0kHz
バンドパスフィルタ(2 to 400Hz)
データ計測時間：50 sec.

⾼覧の⽚側を2Hzで揺する狭帯域のバンドパスフィルタを作⽤して，再び3次元可視化(15.07Hz – 16.17Hz）卓越周波数の把握

3D visualization

Power spectrum

15.7Hzは何？

⼤量のデータを集約し，⼈に伝達するためには結果の⾒える化(可視化)は⾮常に有効．
また，損傷が⽣じた場合には，可視化結果の異常から，損傷箇所を視覚的にローカライズできる可能性がある．

Frequency[Hz]曲げ1次モード

ねじり
1次モード

これは2次のねじりモードであることが分かる

IT基盤の恩恵の享受
学際⼤規模情報基盤共同利⽤・共同研究拠点
jh130013-NA08
構造物の劣化のモデル化とメンテナンス技術の向上に資する⼤規模数値解析
研究拠点：京都⼤学
研究体制：京都⼤（⽜島，⼩⼭⽥，岩下），中央⼤（樫⼭），東北⼤（京⾕，寺⽥，加藤，⾼瀬），岐⾩⼤（永井），九州⼤（浅井），茨城⼤（⾞⾕），産総研（松本），⽔資源機構（藤岡），愛媛⼤（中畑），
敬称略京都⼤学スパコン 32ノード(1024コア 通常期），128ノード(4096コア 集中期）の利⽤
グランドチャレンジ的な問題について⼤規模計算を実施．研究者の計算技術の向上と情報交換．

情報通信研究機構

サイエンスクラウド

研究者が個別の研究環境を⽤意することなく，低コストで科学研究環境（データ処理計算機やデータストレージ）を利⽤できる．ネットワーク管理，システム（サーバ・ストレージ・可視化環境）管理などから解放されて，研究に専念できる．

35

４．将来予測と研究展望
検査装置はコンパクト化，解析はクラウドで今後，中⼩規模のインフラは，スケジュールに基づいて，オンサイト検査が相変わらず実施されるはず．現場における⾼速・⾼精度計算ツールの利⽤は拡⼤．
ただ，検査装置⾃体は必要最低限の機能（センサ部やデータ転送部）を残して，簡素化される．あとはクラウドコンピュータ上でデータ処理・解析を⾏い，検査現場にはアウトプットが表⽰される時代が来ることを想定している．
検査員の技量によらない信頼性のある評価システムの開発
検査員の技量（個⼈知）よりも，インテグレートされたシステムとしての評価⼿法の信頼性（集合知）がキーとなる．逆解析の要素技術⾃体はもちろん重要．さらに，⼤量データの集約，⾼速処理，さらにはアウトプットをクラウドから現場に転送し検査員が状況判断するための情報伝達技術（可視化など）がキー．
技術者教育だからといって個⼈技量は要らないわけでなく，クラウドが出⼒した検査結果に対して，思考（解釈）できるか，それらを活⽤できるか，ある種の想像⼒やデザイン能⼒も必要だと考える．そのためには，⼟⽊⼯学の根幹をなす学問（構造⼒学，⽔理学，⼟質⼒学，etc）の習得だけでなく，IT技術を活⽤するための次世代の教育カリキュラムが不可避だと考える．