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

15-NA17

社会インフラの破壊・非破壊シミュレーションの高度化に資する
大規模数値解析中畑和之（愛媛大学）
コンクリート構造物の老朽化が社会問題となっており，定量的かつ高精度な維持管理技術の構築が望まれている．本課題では，コンクリート部材の劣化・破壊現象の解明か
ら，損傷診断のための非破壊検査手法の構築までを視野に入れた，包括的なモデル化・
シミュレーションの開発を目的とする．ここでは，(A) 複合作用下における破壊現象のモデル化とシミュレーションの実施，(B) 非破壊検査のための波動の順解析・逆解析技
術の高度化，(C) プログラムチューニングおよび数値データ処理の高度化，の研究を 11
名の研究者が分担して遂行している．コンクリートの劣化予測を考える上で，マルチス
ケール・マルチフィジクスを包含した有限要素モデルの導入は不可欠であり，さらに非破壊検査は大規模な波動解析を必要とする．本課題では，計算力学，計算機科学だけでなく，実験力学の研究者が連携した分野横断的体制で研究を行っている．

1. 共同研究に関する情報(1) 共同研究を実施した拠点名

2. 研究の目的と意義

京都大学

国土交通省によれば，建造後 50 年以上経

(2) 共同研究分野

過した社会基盤インフラストラクチャ－（社会インフラ）の割合について，2021 年度には，

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

道路橋で約 28%，河川管理施設で約 40%，海超大規模データ処理系応用分野
超大容量ネットワーク技術分野

岸岸壁で約 29%になるとの集計結果があり，超大規模情報システム関連研究分野

老朽化に伴う事故が懸念されている．近年の(3) 参加研究者の役割分担

財政事情に鑑み，すべて新設の構造物で補う

本研究は，(A) 複合作用下における破壊現ことは困難である．従って，適切な健全度診

象のモデル化とシミュレーションの実施，(B)断による余寿命の評価と，安全を担保した経非破壊検査のための波動の順解析・逆解析技

済的なリノベーションを図ることが現状の術の高度化，を大きな柱とし，研究者が分担最善策であろう．特に，社会インフラの中でして研究を遂行する．また，これらに共通すもコンクリート部材は多用されており，これる計算手法，計算機の課題を改善するために，(C) プログラムチューニングおよび数値デー

らを非破壊的に検査する方法として，X 線透タ処理の高度化を行っている．(C)は主として，

超音波法など多くの手段が提案されている．

構成拠点である京都大学の研究者が中心とな

り，(A)と(B)の研究者と連携して取り組む．(A)しかし，コンクリート部材の内部の状態を定については，FE 解析，破壊力学および材料力いない．

過試験，赤外線，AE，打音法，衝撃弾性波法，

量的に検査する決定的な手法は確立されて学を専門とする以下の研究者，高瀬，猿渡(東コンクリート部材に適した非破壊検査技北大)，車谷(茨城大)，牛島・鳥生(京大)が担術を開発するには，まず，破壊・劣化現象の当する．(B)は波動解析，逆問題および最適化

把握が必要である．コンクリートは，補強用問題を専門とする研究者，京谷・寺田・加藤(東鉄筋と組み合わせて使用される場合がほと

北大)，中畑(愛媛大)，斎藤（群馬大）が担当

んどで，鉄筋とコンクリートの相互作用下にする．(C)，小山田・深沢(京大)が担当する．

おける劣化を考える必要がある．特に，被覆1

学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月するコンクリートは健全であっても内部の並列化コードのチューニング，あるいはプ

鉄筋が腐食した事例も多く，この場合の構造

リ・ポスト処理の最適化など，ソフト・ハー

部材の余耐力を精度よく予測することが肝

ドウエア分野のノウハウも必要である．

要である．また，鉄筋は健全であるが，被覆

これらの技術開発が提示できれば，逆の観するコンクリートが劣化している場合もあ点，すなわち，劣化予想箇所にモニタリング

る．例えば，アルカリ骨材反応，中性化，塩センサを予め設置する，あるいは新設時に維害等が有名である．このような場合，内部ひ持管理が容易になるような構造設計にする

び割れの場合は外側からは見えないため，何等，土木建設工学における全く新しい維持管

らかの方法によってそれを推定することが理技術を提案できる．症状の発生の原因を調必要となる．劣化のメカニズムを把握し，プべ，その原因を取り除く療法，すなわち原因

ロアクティブな非破壊検査技術を確立する療法へと維持管理方針を変えることができことによって効果的な維持管理手法が提案る．従って，現在の維持管理計画に大きく貢できよう．

献するという意味に於いて，本研究の意義は本課題では，破壊現象の解明から長期的な大きい．

劣化予測，さらには損傷度診断のための非破壊検査手法の構築までを視野に入れて，包括

4. 前年度までに得られた研究成果の概要

的な数値解析・シミュレーションを実施する．

H26 年度 JHPCN 課題「社会インフラの破適切なモデルの元で将来起こりうる事象を壊・非破壊シミュレーションの高度化に資す

予測し，プロアクティブな維持管理技術を提る大規模数値解析」に採択して頂いた．以下

案することによって，社会インフラの維持管に，その内容をまとめる．

理の高度化に貢献する手法の開発が本研究(1) 物理・化学作用による破壊現象のモデルの最終目的である．

化とシミュレーションの実施
上述のようにコンクリート部材の劣化は，

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

物理学的要因と化学的要因が複雑に連成し

現在の社会インフラの維持管理は，即時的て生じる．物理的作用によるひび割れ進展をな対応が主であり，対症療法的であると言わ扱うには，コンクリートを非均質材料としてざるを得ない．コンクリートの破壊現象をモ

捉え，メゾスケールのモデル化が必要である．

デル化し，定量的な非破壊検査へと繋げるた象に，X 線 CT 画像を撮影し，断面画像を重

H26 年度は，実際のコンクリート供試体を対間・空間のマルチスケールへ拡張した非常にねることによって 3 次元モデルを作成した．

複雑な問題を解くことになる．関連する物

ここでは，拡張ボクセル有限要素法を応用し，理・化学の専門知識が必要となり，関連する研究者が相互連携して研究を行う必要があ

160x160x160 ボクセル要素で応力解析を行る．また，劣化のメカニズムを熟知し，プロった．非均質な微視構造を有する材料の応力アクティブな非破壊検査技術を確立するた状態を精度良く表現し，現実的な骨材分布でめには，診断・検査を専門とする研究者との破壊シミュレーションを行うことができた．

知識の共有も必要である．さらに，これらのまた，空隙およびクラックに浸透する化学

数理モデリングは必然的に大規模となるこ

物質の流動を精度良く計算するために，多相とは明らかであり，大規模数値解析を効率的場モデルに関する検討を行った．３次元多相

めには，上述のように，物理・化学問題を時に行うためには，計算アルゴリズムの改良や4x4x4 ボクセルを 1 ボクセル要素として，

2場モデルを用いて，コンクリートの亀裂を模学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

擬した多孔質体内の lock-exchange 問題を 3

H25 年度は，ソルバーを大規模計算用にチュ

次元領域で計算し，亀裂内の浸透流についてーニングすることによって，計算速度が格段検討した．間隙に沿った浸透流による熱の移に向上したことが特筆すべき成果である．

動が確認でき，高温流体が(密度小)が上層を，低温流体(密度大)が下層を進行する過程が再5. 今年度の研究成果の詳細(A) 複合作用下における破壊現象のモデル

現された．

(2) 非破壊検査の高度化のためのシミュレー

化とシミュレーションの実施(高瀬，猿渡，車ションと逆解析技術の開発

谷，牛島，鳥生）

超音波アレイ探触子を用いて，きずの位置

上述のようにコンクリート部材の劣化は，や形状を 3 次元的に再構成する手法について

物理学的要因と化学的要因が複雑に連成し

シミュ レー シ ョ ン を 行 っ た ． こ の 方 法 はDelay and Sum 方式による映像化方法であり，て生じる．本年度は，鉄筋とコンクリートの

動弾性有限積分法(EFIT)を用いて，超音波の的なサイズで破壊シミュレーションを行う．

受信波を模擬した．アレイ素子の数(100 個)

また，固体中の温度変化に伴う内圧変化によだけ順解析の回数を要する．また，ボクセル

数は 1000×1000×750(合計 7 億 5000 万)の規模

るコンクリートの破壊現象を想定して，圧縮で，3500 ステップの陽的更新を行った．計算について検討を行った．

並列(Flat MPI)で約 15 分であった．きずを模材の高精度等価フレームモデルの提案

複合材料としての破壊モデルを構築し，現実性流体-固体の連成場に対する数値解析手法

時間は，京都大学スパコンシステム A の 4096

[1] 腐食劣化した鉄筋・鉄骨コンクリート部
鋼材の腐食を考慮した RC 梁の繰り返し載した５つの人工欠陥が再構成できた．
また，粒子フィルタを用いてコンクリート

荷３点曲げ解析を行った．計算は商用コード

内部の構造同定を行った．粒子フィルタでは，(LS-DYNA)を利用した．図 1 に解析モデルを

状態量を表す多数の粒子を数値計算し，計測示す．RC 梁はコンクリート及び主鉄筋を値との尤度に基づいて粒子の生成・消滅を行

SOLID 要素でモデル化し，せん断補強筋は梁う．本研究では，コンクリート中の骨材の体要素でモデル化した．有限要素メッシュ数は
約 220 万要素となった．鉄筋は応力ひずみ曲積含有率(骨材率)を超音波エコーから推定することを試みた．Intel Xeon Phi コプロセッサを利用した MIC 並列計算を EFIT に導入し，線を多点近似によって表現した等方弾塑性体とした．降伏応力ρy は 393MPa，引張強さ

超音波エコー計算を高速化した．粒子フィル

ρt は 582MPa である．コンクリートは圧縮強タの更新に伴い，平均値が真値に近づくこと度を元に各パラメータを設定できる連続局を示した．

(3) プログラム高度化および数値データの効面キャップモデルを適用した．コンクリート

グループメンバによる対面会議を 1 回開催

繰返載荷はモデル中央において強制変位
によって行った．鋼材の腐食は主鉄筋の中央スによって，(1)有限要素解析の大規模化および高速化，(2)Intel Xeon Phi コプロセッサをとし，腐食による断面欠損割合に応じて断面用いた波動解析の高速化が可能となった．

重曲線，図 3 にひび割れ進展状況を示す．

の圧縮強度ρc は 36MPa とした．

率的な処理部から軸方向両側 175mm の区間を腐食区間した．また，拠点側からの技術的なアドバイ

積を減じることでモデル化を行った．腐食率

は 0%，14%，24.5%とした．図 2 に変位－荷(1)の有限要素法は陰的に更新する方法であり，これまで計算効率が芳しくなかった．

3

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

図 1: モデル概略図

図 4: 強制加振モデル

図 2: 変位－荷重曲線
図 5: 横軸：空隙（模擬損傷）長さ，縦軸：共振周
波数[2] 鉄筋の弾塑性挙動とそれに起因するコンクリートのひび割れ進展のシミュレーショ
図 3: ひび割れ進展状況

ン

計算時間は 512 並列で 32 時間であった．

鉄筋とコンクリートの非線形力学挙動を

腐食領域の高精度なモデル化や損傷の詳細モデル化した数値解析手法を用いて，せん断な表現にはさらなるモデルの微細化が必要

補強筋の異なる RC はりの破壊挙動を 3 次元であり，今後並列数の大規模化により解決をで詳細にシミュレートした．前年度も同じ内図る．

容を報告しているが，本年度は実験と詳細な

また，模擬損傷を導入した試験体モデルに比較を行ったので，それについて述べる．

対する強制加振試験のシミュレーションを

破壊シミュレーションには，
鉄筋に von-Mises
塑性モデル，コンクリートに修正 von-Mises

領域拡大に伴う振動応答の変化を精度よく損傷モデルを用いた非線形有限要素法に基予測可能であることを確認した．

づく陰的な解析手法を採用している．鉄筋の行った(図 4)．図 5 に示すように，梁の損傷

ヤング率は 200GPa，初期降伏応力は 380MPaとし，コンクリートのヤング率は 29GPa，破壊発生ひずみは，2～3 MPa で引張破壊が生4

じるように定めている．
圧縮引張強度比 k は，

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

損傷モデルにおいて，理想的な一軸引張と一を比較すると，目視で確認されたひび割れの

軸圧縮に対する強度比である．理想的な一軸

本数や分布を精度よく再現できている．ひび

圧縮の要素試験を行うことができれば，円柱割れの発生パターンに着目すると，最初に曲供試体の圧縮強度よりも大きな値になる可げひび割れが発生した後，せん断補強筋が少能性があるため，ここでは，圧縮引張強度比
k を 10 と 20 で設定することとする．破壊シなくなるにつれて，せん断破壊を示す典型的

ミュレーションに用いる有限要素モデルをートされている．

な斜めひび割れに発展する様子がシミュレ

図 6 に示す．ここでは，せん断補強筋 10 本の RC モデル（スパン 1200mm）を作成した．

Experiment S10: 120 kN

各モデルの要素数は，四面体要素を用いて，約 300 万要素である．非線形解析における増分幅がほぼ一致するように，各モデルの強制変位とステップ数を設定し，各モデルで材料0.00010.10

Equivalent strain distribution at load step: 160/500 (120 kN)

パラメータを変えずに破壊シミュレーショ
ンを行い，数値解析と実験結果を比較する．

S10yInternal damage distribution at load step: 160/500 (120 kN)

図 7: 実験と数値解析のひび割れ分布の比較

x
z

Es

200 GPaEc

29 GPaMaterial parameter of steel
s y0
ns
Q

0.3

380 MPa

150 MPab
15

0.2

0.0001

0.1 N/mm

10, 20Material parameter of concrete
k0
Gf
nc

k

図 6: RC はりの数値モデル

図 8: せん断補強筋 0,6,10 本の RC はりの荷重－変位関係の比較

図 7 に，実験と数値解析におけるひび割れ
分布の比較を示す．実験結果は，目視で確認

また，図 8 は，実験と数値解析における荷されたひび割れの分布である．数値解析の結重－変位関係の比較を示す．ここでは，せん

果は，モデル表面での等価ひずみ分布とモデ

断補強筋が 0 本と 6 本の場合の数値解析も行

ル内部の損傷変数の分布を示している．本研

い，実験結果と比較している．実験結果と数

究では，コンクリートのひび割れを損傷によ値解析結果ともに，せん断補強筋の数が減る

るひずみの局所化でモデル化しているため，につれて，延性的な曲げ破壊から脆性的なせひずみの分布と損傷変数の分布をひび割れん断破壊を示している．材料パラメータを変と見なすことができる．解析結果と実験結果

5

えずに，せん断補強筋のモデル化の違いのみ

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

10 倍に設定したため，特に t = 0.012 [s] (図で，数値解析結果は実験結果を定量的に再現

10 (a)) から分かるように，固体部分におけるできている．

[3] 固体中の温度変化に伴う内圧変化のシ

等温線の間隔は空隙内部に比べて大きくなっミュレーション

ており，熱伝導率の差異を反映した妥当な結

２次元の圧縮性流体・固体熱連成計算手法

果が得られたと考えられる．

を用いて，温度変化に伴う固体中の空隙における内圧変化の評価シミュレーションを行った．本計算手法では流体の圧縮性が考慮されるため，固体の温度変化によって内部の空隙で生じる圧力変化を適切に計算することが可能である．また，本計算手法では固体と流体間に境界条件を設定する必要がなく，単純な直交構造格子を用いて固体内部の熱伝導と，任意形状の空隙内部における流体の対流を統一的に計算することが可能である．

計算領域を図 9 に示す．図 9 に示されるように，対象とする固体の空隙率は 7.33%であ

図 9: 計算領域 (長さの単位：[mm])る．各空隙については，計算領域を 2.5 [mm]幅の正方形格子で分割し，各格子内に直径 0.5〜1.0[mm]の空隙をランダムに配置した (合計 2,048 個)．空隙内の流体は理想気体とし，プラントル数は 0.717，比熱比は 1.4 とした．
固体の熱伝導率は流体の 10 倍とし，その他の物性値については初期状態の流体と同じ(a) t = 0.012 [s]とした．計算領域の左側面および右側面の温度 Th，Tc はそれぞれ 500 [K]，300 [K]で一定(b) t = 0.180 [s]図 10: 各時刻における等温線(領域：0.0≦x1≦1.0 [mm]，79.0≦x2≦80.0 [mm]，等温線の間隔：10 [K])とし，上部および底部は断熱条件とした．また，計算領域内の初期温度 Tc は 300 [K]，空隙内の初期圧力 P0 は 1.01325×105 [Pa] とし温度変化に伴う空隙内部の圧力変化を確た．計算格子数は各方向に 3,200×6,400，並列

認するために，計算領域を図 11 のように間利用し，並列プロセス数は各方向に 32×64 の領域内の無次元平均温度 T*，空隙内の無次元計算手法として Flat MPI による領域分割法を

隔 2.5 [mm]で x1 方向に分割し，分割された各

合計 2,048 並列とした．

平均圧力 P*を算出した．なお，T*，P*は，

図 10 に 0.0≦x1≦1.0 [mm]，79.0≦x2≦80.0，[mm]の領域 (図 9 で拡大された領域) におけ，と定義される．ここで，る各時刻の等温線を示す．なお，図 10 の点線は分割された各領域の平均温度 (単位：[K])，は流体と固体の境界面を表す．図 10 に示されは各領域に含まれる空隙内圧力の平均値(単位：[Pa])である．また，T*，P*の定義点

るように，計算開始と共に左側面で固体が加は各領域の中心点とした (図 11 の黒丸)．

熱され，領域内の温度が上昇していく状況が確認できる．また，固体の熱伝導率を流体の6

図 12 に t = 0.18 [s]における T*および P*の水

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

平方向分布を示す．なお，図 12 において x1*

微視構造に内在するマイクロクラックを予測

は水平方向座標 x1 を計算領域の水平方向長する手法（逆解析）についても示す．

さ 80 [mm] で無次元化した無次元座標であ[1] Biot モデルを用いた波動の散逸を考慮しる．図 12 に示されるように，加熱面 (x1* = 0)

た大規模波動解析

付近では温度，空隙の内圧が共に上昇してい固体中の間隙に流体が飽和した材料の波動

る．各空隙はそれぞれが完全に閉じており，伝搬問題は，Biot によりモデル化されており，外部への流体の移動は生じない．したがって，

飽和多孔質弾性波動問題と呼ばれる．飽和多各空隙内部における流体の状態変化は定積

孔質弾性波動問題を波動問題に有効な手法の変化であり，温度上昇に伴って空隙内部の圧

1 つである時間領域境界要素法を用いて解析力が上昇するという計算結果は妥当である

する場合，時間領域境界要素解析に必要な基と考えられる．

本解を閉じた形で得ることが難しいという問以上，本研究では圧縮性流体と固体の熱連題点が生じる．そこで，本研究では，時間領成計算手法を用いて，固体を加熱した際に生域境界要素法の定式化で必須な境界積分方程じる空隙の内圧上昇を計算できることを確式に現れる畳み込み積分に，演算子積分法と認した．

呼ばれる畳み込み積分の近似解法を適用する
ことで，この問題を解決した．演算子積分法
を用いることで，時間領域基本解を直接使用することなく，ラプラス変換域基本解を用いることで，数値解を求めることに成功した．

(a) T*の分布

図 11: 一定間隔で x1 方向に分割された計算領域 (単位：[mm]，各領域の x1 方向：2.5 [mm] )(B) 非破壊検査のための波動の順解析・逆解
析技術の高度化（京谷，寺田，加藤，中畑，齋藤）
ここでは，コンクリート中の弾性波の伝搬および散乱を模擬するために，Biot モデルを(b) P*の分布用いた大規模境界要素解析（順解析）を行っ

図 12 : x1 方向に分割された各領域内における T*およた．また，コンクリート中のきずを再構成すび P*の分布 (t = 0.18 [s])るためのアレイ映像化（逆解析）についても検討を行っている．また，コンクリートのある位置における異方性材料剛性から，材料の7

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

ノードに 1MPI プロセスを割当て，合計 64 ノ

一方，時間領域境界要素法の欠点の 1 つでード使用し，1 ノード当たり 32OpenMP スレ

ある計算負荷の低減についても検討を行ってッドで並列化を実行している．図 13(a)より入いる．まず，MPI-OpenMP のハイブリッド並

射縦波が先頭列の球形空洞により散乱され，列化を開発した手法に適用することで，3 次

図 13(b)より 2 列目の空洞にも散乱されている

元飽和多孔質弾性体中に存在する球形空洞に様子を見て取れる．なお，この場合の境界要よる入射波の散乱問題を解析した．解析結果の一部を図 13 に示す．空洞は全部で 28 個で素解析に要した計算時間(内点計算を除く)はあり，1 つの空洞当たり 512 三角形要素で離散化した．また，全時間 step 数は 64 とした．

速化を目指し、高速多重極法を適用すること

64=3670016（空洞数×要素数×1 要素当たりの[2] 全波形サンプリング処理によるコンクリ

およそ 1.7 時間だった．現在は，さらなる高そ の た め ， 全 自 由 度 は 28 × 512 × 4 ×

自由度×全時間 step 数）である．境界要素法を行っている．

ートの内部鉄筋・空隙の再構成を使用しているため，空洞の表面のみを離散超音波アレイ探触子を用いて，きずの位置化していることに注意されたい．並列化は 1や形状を 3 次元的に再構成する手法について検討を行った．再構成原理は，全波形サンプ
リング処理方式(FSAP 方式)をベースとして
おり，平面アレイ状に配列した素子で得られる散乱波を基に，散乱体(きずや鉄筋)を再構成する手法である．ここでは，動弾性有限積
分法(EFIT)を用いて，シミュレーションで生成した波動から，コンクリート中の鉄筋の形状を再構成する．FSAP 方式では，アレイの素子数分，超音波の送受信をする必要があり，EFIT のシミュレーションを繰り返す．

図 14 に示すように，
500×500×350[mm]のコ

ンクリートの上部にマトリクスアレイ探触子(16×16=256 素子，素子ピッチ 14mm)の設置を想定する．この数値モデルのボクセル数は1000×1000×700(合計 7 億)である．モデル内には直径 20 mm の鉄筋が深さ方向に平行に 5つ存在する．数値解析では，探触子が接触するボクセルの垂直応力グリッドを，200kHzの入射波形で励起することで超音波を発生させている．EFIT は，3600 ステップの陽的更新を行う．

図 13: 飽和多孔質弾性体中の球形空洞による入射縦波の散乱(a)20step，(b)62step における空洞周辺の変位場8学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

Array elements (256ch)
223

500
200
150

x3

x2
x1

500

0.0

223

50

50

250

x2

x1 Cross section N-2

R

-1.0

(a)

Cross section N-1

(mm)

14mm

16 elements

250

Cross section N-1

x3

20

1.0

(b)
Cross section N-2

16 elements
図 14: 3 次元コンクリートモデル(上)とマトリクスア
レイ探触子のモデル(下)(c)

図 15: 5 つの鉄筋の再構成シミュレーション，（a）等

ある１素子から超音波を送信し，各素子で散乱波を受信する．送信素子を変えながら，

値面，(b)と(c)断面表示

み合わ せの 散 乱 波 形 を 得 る ． こ の 波 形 を[3] 有限要素法を用いた逆均質化法による内これを繰り返して，計 256 素子×256 素子の組
FSAP 方式に入力し，開口合成処理を行うこ

部ひび割れ分布の同定とで，鉄筋の形状を再現した結果を図 15 に劣化の主な要因であるひび割れが発現す示す．図 15 の断面図からも，5 つの鉄筋は識る初期段階にはマイクロクラックという状

別できることがわかるが，コンクリート境界

態があり，将来的な構造物の劣化を引き起こ面付近の鉄筋は再構成されていない．これは，す要因となるため，早期発見が重要である．

マトリクスアレイ探触子から発振した超音

ここでは，コンクリート材料のある位置にお波ビームが，鉄筋に鋭角で入射するために，けるマクロ的な異方性材料剛性（観測剛性

D*）が既知である仮定とした場合に，その材

探触子への反射波が小さくなるためである．
京都大学スパコンシステム A の 2048 並列(Flat MPI)では約 25 分を要した．従って 256

料微視構造に内在するマイクロクラックのション，計 107 時間程度を要した．

を開発し，数値シミュレーションを実施した．

劣化状況を予測するトポロジー最適化手法

素子すべての波形を得るのに 256 シミュレー

解析には京都大学スパコンシステム A を用9いた．材料微細構造は周期境界を有し，材料

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

はコンクリートを想定してヤング率 25 [GPa]，

もので，f =1 のときはマイクロクラックの存

ョンでは，一方向の単純なマイクロクラックのθ=30°, 45°の同定結果について観察す

ポアソン比 0.167 とした． このシミュレーシ
を 対 象 に ひ び 割 れ 角 度 θ =0, 30, 45, 120,

在しない健全な状態を意味する．まず，図 17135°の 5 ケースについて同定を行った．同る．同定によって得られた均質化材料剛性

定精度の検証が目的であるため，ターゲットとなる観測剛性 D*については，図 16 に示す

非常に近い値となっており，高精度で同定さ5 つのマイクロクラック分布を準備し，それれたことが確認できる．また，得られたトポ

ぞれの観測剛性 D*を均質化法によって求め，

ロジーを見ても，いずれもマイクロクラック

これを所与として扱うものとした．

当該同定問題は，一意に解が求まらない逆問

CH は，いずれの成分も観測剛性 D*のそれとの方向は再現できているといえる．ちなみに題を扱っているため，例え数値的な同定精度が高くても様々なクッラク分布が解として得られることは分かっている．

図 16：所与のクラック分布と観測剛性

図 17 は，その同定結果である．また，f は
目的関数値である材料体積量を正規化した

図 17：同定されたクラック分布と均質化剛性10学際大規模情報基盤共同利用・共同研究拠点 平成 27 年度共同研究 最終報告書 2016 年 5 月

次に，θ=120°, 135°について同定を行っ[2]鉄筋コンクリートの大規模有限要素モ

た結果を考察する．左右対称となるθ=120°

デルを作成し，鉄筋の塑性変形とコンクリーとθ=30°，あるいはθ=135°とθ=45°の観トのひび割れ進展挙動をそれぞれ解析することで，鉄筋コンクリートの破壊挙動を 3 次

測剛性および均質化剛性の成分をそれぞれ比較すると，(1, 3)成分と(2, 3)成分の符号が異元で詳細に再現することを予定している．

なっているだけで，その成分の大きさはほぼ[3] 今年度は，圧縮性流体と固体の熱連成

同じとなっていることが分かる．これは，力計算手法を用いて，固体を加熱した際に生じ学的にも正しいもので均質化解析の妥当性

る空隙の内圧上昇を計算できることを確認と同定精度の高さを裏付けるものであると言える．その結果，θ=120°, 135°のマイクした．今後は，鉄筋の有無など，内部の構造

ロクラック分布も所与の角度を再現できていて検討していく予定である．

が温度や空隙の圧力上昇に与える影響についると言える．一方，θ=0°の同定結果をみ(B) 非破壊検査のための波動の順解析・逆解

ると，得られたマイクロクラック分布の方向

析技術の高度化[1]MPI-OpenMP を用いたハードウェアには再現できているといえるが，材料剛性の同定精度については他の 4 ケースの結果と比べよる高速化のみならず，高速畳み込み演算をてやや慎重に観察する必要がある．具体的には，均質化剛性 CH の(1, 1)成分は，観測剛性

用いた時間方向の高速化，積分核に依存しな

D*のそれに対し精度よく同定できているが，い高速多重極法を用いた空間方向の高速化

その他の成分はそもそも観測剛性 D*の成分

実現についても検討していく．

が小さ過ぎるため，十分零に近い数値を示せ

レーションへの拡張を行いたい．また，コン

たものの，所与の値のオーダー（桁数）に対クリート用のリニアアレイ探触子を試作ししてはやや隔たりがあるように思える．これており，このシミュレーション結果と実験結は，観測剛性 D*の(2, 2)成分が極めて小さい，[2]骨材による波動減衰を考慮したシミュ

果を比較する予定である．

すなわち，ひび割れに対して垂直である鉛直[3]については，これまで一方向のマイクロ方向に対して力学的にほとんど剛性がない

クラックという単純なモデルを対象としてような極限的な状態を数値材料試験で同定いたがマイクロクラックが途中で分岐するしようとしているので，この結果に対する見ような複雑なモデルに対して，どの程度の精解としては，同定精度が悪化したというより度が得られるか，またその精度を担保するたも数値的に限界の状態を表している．

めの手法の開発と大規模計算にチャレンジ
していく．

(C) プログラムチューニングおよび数値デー

6. 今年度の進捗状況と今後の展望(A) 合作用下における破壊現象のモデル化

タ処理の高度化とシミュレーションの実施

拠点側からの技術的なアドバイスによっ[1] 今後，載荷試験により損傷を導入したて，(A)の有限要素解析の大規模化および高実機試験結果のシミュレーションを行い，試速化に取り組んだ．今後は，チューニングし験結果と比較することで，RC 構造物の破壊たコードを用いて大規模な計算にチャレン

現象及びその評価手法の精度を確認する．まジしていく予定である．他の問題における大た，腐食鉄筋の膨張など，コンクリートと鉄規模計算の効率化についても継続的に議論筋の相互作用について現実的なモデル化をを行っていく予定である．

取り入れていきたい．

11

以上より，本年度は(A)，(B)，(C)について，

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

個々の研究者が独立した成果を上げている．

破壊力学的損傷を考慮したコンクリート中の個々の研究は予定通り遂行されているが，い物質移動解析手法の開発とその性能評価，土ずれは(A)で予測した破壊モデルを用いて(B)pp.I_161-I_170, 2015．

木学会論文集 A2（応用力学），Vol.71, No.2,

で非破壊検査シミュレーションを行う等，(A)と(B)を有機的に連携していくことを予定

8. M. Kurumatani, K. Terada, J. Kato, T. Kyoya, K.

している．そのためにはより大きな計算スケ

fracture mechanics for concrete, Engineering

Kashiyama, An isotropic damage model based onールとなると予想されるが，現在の計算機資Fracture Mechanics, Vol.155, pp.49-66, 2016.

源で効率的な計算を行えるよう(C)との連携も図りたい．

(2) 国際会議プロシーディングス

7.

9.

研究成果リストD.

Toriu,S.

Ushijima,

Applicationofcomputational method for thermal interactions(1) 学術論文

1. K. Nakahata, S. Tokumasu, A. Sakai, Y. Iwata, K.

between compressible fluids and solids to natural

Ohira, Y. Ogura, High-speed ultrasonic imaging

convection

using scattering amplitude with a flexible arrayaroundcircularcylinders,

WCCM&APCOM 2016, Seoul, Korea, 2016.(掲

transducer, NDT & E International, Vol.82,

載予定)

pp.13-25, 2016.

2. 車谷麻緒，根本優輝，相馬悠人，寺田賢二郎，(3) 国際会議発表

コンクリートの破壊力学を考慮した鉄筋コン

10. T. Saitoh, M. Shimoda, A. Furukawa and S.

性能評価，日本計算工学会論文集，Vol.2016,

for elastic wave scattering in 2-D anisotropic

3. 中畑和之, 堀口貴志, GPU 計算を実装した全computational methods (ICCM2015), Pullman

Hirose, Fast multipole boundary element method

クリートの 3 次元破壊シミュレーションとそのsolids,pp.20160004, 2016．

6th

International

conferenceonHotel, Auckland, New Zealand, 2015.

波形サンプリング処理方式による内部欠陥の三次元超音波映像 化, 日本機械学会論文 集,

11. A. Furukawa, T. Saitoh and S. Hirose: A

4. A. Ibrahim, A. Zabri, K. Nakahata, Identification

fluid-saturated porous

Vol.81, No.832, 15-00471, 2015.

boundary element method for wave scattering inof elastic parameters of an inclusion by a particle

rocks,

BEM/MRM38,

EM/MRM38, Balmer Lawn Hotel, New Forest,

filter using ultrasonic waves, International Journal

England, 2015.

of Mechanical and Materials Engineering, Vol.10,

12. M. Kurumatani, S. Abe, Y. Nemoto, Y. Soma, N.

2015.

concrete using an isotropic damage model, 13th

No.1, Article 23, DOI:10.1186/s40712-015-0050-y,

Henmi, A method for fracture simulation ofUS

5. 丸山泰蔵, 斎藤隆泰, 廣瀬壮一, 時間方向に

National

Mechanics, 2015.

高速化された 3 次元音響波動問題に対するCongressonComputational

CQ-BEM, 計 算 数 理 工 学 論 文 集 ,Vol.15,

13. M. Kurumatani, S. Abe, Y. Nemoto, Y. Soma, N.

6. 車谷麻緒, 川瀬晴香, 拡張ボクセル有限要素

reinforced concrete using a damage model based

12-151204, 2015.

Henmi, A method for fracture simulation of法の開発とその性能評価, 日本計算工学会論文

onfracture

7. 小林賢司，車谷麻緒，岡崎慎一郎，廣瀬壮一，Computational Mechanics, 2015.

International

集, Vol.2015, pp.20150011, 2015.

12mechanics

Workshops

foronconcrete,

Advances

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

23. 中畑和之, 平岡拓也, 金井淳, 廣瀬壮一, メ(4) 国内会議発表

14. 中野葵，車谷麻緒，鉄筋コンクリートの破ゾスケール構造に着目したコンクリート中の壊シミュレーション結果の可視化方法に関する研究，第 43 回土木学会関東支部技術研究発

24. 古川陽, 斎藤隆泰, 廣瀬壮一, 周波数領域

表会，CD-ROM，2016．

15.相馬悠人，阿部俊逸，車谷麻緒，幾何学的非境界要素法を用いた異方性飽和多孔質弾性体

線形性を考慮した損傷モデルと塑性モデルに

学術講演会, CD-ROM, 2015.

よる鉄筋コンクリートの破壊シミュレーショ

25. 丸山泰蔵, 斎藤隆泰, 廣瀬壮一, 3 次元音響

CD-ROM，2016．

高速化について, 土木学会第 70 回年次学術講

超音波伝搬特性に関する考察, 土木学会第 70回年次学術講演会, CD-ROM, 2015.

の 3 次元波動散乱解析, 土木学会第 70 回年次

ン，第 43 回土木学会関東支部技術研究発表会，
16.星裕斗，相馬悠人，車谷麻緒，鉄筋の幾何形

波動問題に対する CQ-FMBEM の時間方向の演会, CD-ROM, 2015.

状を考慮した RC はりの 3 次元破壊シミュレ

26. 相馬悠人, 根本優輝, 車谷麻緒, 損傷モデ

支部技術研究発表会，CD-ROM，2016．

壊シミュレーション, 土木学会第 70 回年次学

数値解析による RC はりの破壊進展挙動の比27. 邊見哲一, 車谷麻緒, 実験・計測・解析に

表会，CD-ROM，2016．

礎的研究, 土木学会第 70 回年次学術講演会

内部ひび割れを考慮した鉄筋コンクリート中

28. 車谷麻緒, 小林賢司, 岡崎慎一郎, 損傷モの塩化物イオン濃度の定量的予測手法に関する研究，第 43 回土木学会関東支部技術研究発

デルを用いたコンクリート中の物質移動解析

表会，CD-ROM，2016．

19.邊見哲一，車谷麻緒，画像解析を用いた RC

年次学術講演会, CD-ROM, 2015.

29. 加藤準治, 河西亮輔, 京谷孝史, 寺田賢二

はりに発生する内部ひび割れ進展挙動の計測，
第 43 回土木学会関東支部技術研究発表会，郎, 逆均質化法を用いた微細亀裂分布の同定,

CD-ROM，2016．

20. A. Zabri, K. Nakahata, Study on identification

概要集, CD-ROM, 2015.

30. 根本優輝, 相馬悠人, 車谷麻緒, 鉄筋コン

particle filter, 土木学会第 70 回年次学術講演

第 20 回計算工学講演会, CD-ROM, 2015.

ーションに関する研究，第 43 回土木学会関東

ルと塑性モデルを用いた RC はりの 3 次元破17.小圷祐輔，邊見哲一，車谷麻緒，画像計測と術講演会, CD-ROM, 2015.

較考察，第 43 回土木学会関東支部技術研究発

よる RC はりの破壊力学挙動評価に関する基18.安藏尚，小林賢司，車谷麻緒，岡崎慎一郎，

CD-ROM, 2015.

手法の開発とその性能評価, 土木学会第 70 回

of wave velocity and density of defect using

土木学会第 18 回応用力学シンポジウム講演

クリートの非線形有限要素解析とその V&V,

会, CD-ROM, 2015.

31. 古川陽, 斎藤隆泰, 廣瀬壮一, 異方性飽和

レイアンテナを用いた電磁波レーダ法による要素法, 第 20 回計算工学講演会, CD-ROM,

CD-ROM, 2015.

32.小林賢司, 車谷麻緒, 岡崎慎一郎, 損傷モデ

21. 岸岡大樹, 中畑和之, 清 良平, 無指向型ア

多孔質弾性体の 3 次元波動問題に対する境界

鉄筋再構成, 土木学会第 70 回年次学術講演会,
22. 泉英輝, 堀口貴志, 中畑和之, アレイ映像

ルを用いた損傷後の鉄筋コンクリート中の物
質 移 動 解 析 , 第 20 回 計 算 工 学 講 演 会 ,化のための全波形サンプリング処理方式の性
能検証,

CD-ROM, 2015.

土木学会第 70 回年次学術講演会,

CD-ROM, 2015.

2015.

13

33. 阿部俊逸, 車谷麻緒, 寺田賢二郎, 幾何学

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

的非線形性を考慮した準脆性材料の損傷解析,
第 20 回計算工学講演会, CD-ROM, 2015.

34. 相馬悠人, 根本優輝, 車谷麻緒, 損傷モデ
ルと弾塑性モデルを用いた鉄筋コンクリートの破壊シミュレーション, 第 20 回計算工学講
演会，CD-ROM, 2015.

35. 邊見哲一, 根本忍, 車谷麻緒, 画像解析を用いたコンクリートの破壊シミュレーションの V&V, 第 20 回計算工学講演会, CD-ROM,
2015.

36. 車谷麻緒, 川瀬晴香, 非均質構造解析のための拡張ボクセル有限要素法, 第 20 回計算工学講演会, CD-ROM, 2015.

37. 番場良平, 新宅勇一, 村松眞由, 高瀬慎介,
森口周二, 寺田賢二郎, 堤成一郎, 地震荷重を受ける鋼製橋脚のき裂進展評価, 第 20 回計算工学講演会, CD-ROM, 2015.

38. 新宅勇一, 村松眞由, 高瀬慎介, 森口周二,
寺田賢二郎, 堤成一郎, 車谷麻緒, 有限被覆法による自動き裂進展解析, 第 20 回計算工学
講演会, CD-ROM, 2015.

39. 新宅勇一, 村松眞由, 高瀬慎介, 森口周二,
寺田賢二郎, 堤成一郎, 埋込型結合力モデルを導入した損傷構成則, 第 20 回計算工学講演会, CD-ROM, 2015.

(5) その他（特許，プレス発表，著書等）
なし

14