董 勤 喜

(教育科學(xué)研究所，日本 川崎市 214-0033)

工程結(jié)構(gòu)的抗震分析通常采用地震波作為地震荷載輸入，如中國(guó)1975年的海城波、美國(guó)1940年的E1-Centro波及1952年的Taft波等。這些地震波與震源斷層的基巖、覆蓋土層等當(dāng)?shù)氐牡刭|(zhì)情況密切相關(guān)，事實(shí)上不具備通用性。地震是斷層的錯(cuò)動(dòng)所致，因此，可以直接利用斷層模型計(jì)算地震動(dòng)，再利用波動(dòng)理論計(jì)算某一地點(diǎn)的地震動(dòng)。以此為基礎(chǔ)，計(jì)算得到的地震動(dòng)波形將更有應(yīng)用價(jià)值。這是近年來(lái)許多學(xué)者正不斷努力的方向。Nishimura等[1]針對(duì)2011年日本東北部海域發(fā)生的里氏9.0級(jí)大地震出現(xiàn)的同震變形，采用基于開(kāi)源并行計(jì)算環(huán)境FrontISTR[2]的EduS/FrontISTR軟件，直接模擬地殼內(nèi)的大逆沖型震源斷層，成功再現(xiàn)了日本新瀉—神戶間軟弱沉積平原區(qū)域內(nèi)的應(yīng)變集中現(xiàn)象。全球衛(wèi)星觀測(cè)系統(tǒng)和合成孔徑雷達(dá)干涉測(cè)量技術(shù)InSAR的綜合應(yīng)用，使得觀測(cè)同震地震變形的詳細(xì)分布成為可能。這為大規(guī)模地震分析研究提供了良好的實(shí)測(cè)數(shù)據(jù)，使得大規(guī)模地震數(shù)值模擬分析成為可能。嚴(yán)珍珍等[3]基于橫向各向同性PREM地球模型，利用CMT(centroid moment tensor)震源機(jī)制，考慮地表地形及海洋等地球特性，采用Jeroen Tromp等提供的SPECFEM3D_GOLBE譜元法軟件，進(jìn)行了大規(guī)模地震數(shù)值分析，較準(zhǔn)確地重現(xiàn)了長(zhǎng)周期理論頻率值。許才軍等[4]采用彈性半空間矩形位錯(cuò)模型，結(jié)合GPS和InSAR觀測(cè)結(jié)果對(duì)日本東北部海域發(fā)生的里氏9.0級(jí)大地震引起的同震滑動(dòng)分布進(jìn)行了反演分析。

目前，可以進(jìn)行大規(guī)模地震數(shù)值模擬分析的常用軟件[5]見(jiàn)表1所示。

表1 常用大規(guī)模地震波數(shù)值模擬軟件

本文將詳細(xì)介紹基于并行計(jì)算有限元法和斷層模型研制的EduS/FrontISTR軟件及其在地震數(shù)值模擬方面的基本原理，并通過(guò)標(biāo)準(zhǔn)地震模型對(duì)軟件的性能進(jìn)行數(shù)值模擬驗(yàn)證。最后，以2011年日本東北地區(qū)太平洋近海地震為例，應(yīng)用EduS/FrontISTR軟件對(duì)地震發(fā)生后的地震波的傳播過(guò)程等進(jìn)行大規(guī)模地震數(shù)值模擬研究，很好地再現(xiàn)了觀測(cè)結(jié)果。

1 并行計(jì)算有限元地震數(shù)值模擬軟件

筆者基于開(kāi)源并行有限元程序開(kāi)發(fā)平臺(tái)FrontISTR[2]，研究開(kāi)發(fā)了具有自主知識(shí)產(chǎn)權(quán)的EduS/FrontISTR軟件。該軟件是一款通用的高性能并行非線性有限元計(jì)算軟件，適用于個(gè)人PC機(jī)到超級(jí)計(jì)算機(jī)系統(tǒng)。該軟件可以求解常規(guī)的靜態(tài)和動(dòng)態(tài)線性問(wèn)題，也可以求解幾何非線性、材料非線性和邊界非線性(含接觸)等問(wèn)題，還可以有效模擬從微小滑動(dòng)到有限滑動(dòng)問(wèn)題、從無(wú)摩擦接觸到有限摩擦接觸。EduS/FrontISTR軟件內(nèi)的子模塊EduS/FrontISTR/Quake可以進(jìn)行大規(guī)模的地震數(shù)值模擬計(jì)算，包含以下功能模塊：1)斷層模塊和震源模塊；2)邊界能量耗散吸收模塊；3)考慮地表地形和海底形狀的快速網(wǎng)格自動(dòng)生成及地震波后處理模塊；4)通過(guò)大數(shù)據(jù)結(jié)構(gòu)，嵌入非勻質(zhì)地下結(jié)構(gòu)J-SHIS[13]和松原地質(zhì)結(jié)構(gòu)模型[14]等。在EduS/FrontISTR高性能并行計(jì)算支撐平臺(tái)下，運(yùn)行子程序EduS/FrontISTR/Quake，可以快速地實(shí)現(xiàn)地震波的傳播模擬。

1.1 并行有限元法

大規(guī)模地震數(shù)值模擬的計(jì)算區(qū)域很大，需要將計(jì)算區(qū)域分割成若干子區(qū)域。在EduS/FrontISTR軟件中，采用節(jié)點(diǎn)法對(duì)計(jì)算區(qū)域進(jìn)行分割。以圖1(a)所示的計(jì)算區(qū)域?yàn)槔?，將?jì)算區(qū)域分割為4個(gè)子區(qū)域，然后，將各子區(qū)域重新按有限元法進(jìn)行整理，如圖1(b)所示。

計(jì)算區(qū)域分割時(shí)，圖1(a)中各子區(qū)域內(nèi)部的節(jié)點(diǎn)稱為內(nèi)點(diǎn)，圖1(b)中各子區(qū)域重新整理時(shí)各單元引入的其他子區(qū)域的節(jié)點(diǎn)稱為外點(diǎn)，定義一個(gè)子區(qū)域的內(nèi)點(diǎn)是其他子區(qū)域的外點(diǎn)的節(jié)點(diǎn)為邊界點(diǎn)。節(jié)點(diǎn)編號(hào)原則是先內(nèi)點(diǎn)，后外點(diǎn)。以圖1(b)中子域0為例，節(jié)點(diǎn)1到7為內(nèi)點(diǎn)，節(jié)點(diǎn)8到15為外點(diǎn)，節(jié)點(diǎn)1、2、4和6為子域1的邊界點(diǎn)，節(jié)點(diǎn)1到3為子域2的邊界點(diǎn)，節(jié)點(diǎn)1為子域3的邊界點(diǎn)。在軟件中，采用通信表建立邊界和外點(diǎn)的關(guān)系。各子區(qū)域構(gòu)建完成后，采用分布式數(shù)據(jù)存儲(chǔ)，再按常規(guī)的有限元法對(duì)各子區(qū)域獨(dú)立進(jìn)行有限元計(jì)算，各子域間僅利用通信表進(jìn)行數(shù)據(jù)信息交換。

圖1 并行有限元法的計(jì)算區(qū)域分割

1.2 地震斷層模擬

為模擬地震斷層，子程序EduS/FrontISTR/Quake采用分裂節(jié)點(diǎn)法[15](Split-node)處理正斷層和逆斷層處的有摩擦平面剪切裂紋，如圖2所示。通過(guò)分裂節(jié)點(diǎn)法，可以對(duì)斷層的上面和下面指定不同的滑動(dòng)位移，但兩者的相對(duì)位移等于斷層的位錯(cuò)量。

圖2 分裂節(jié)點(diǎn)法

1.3 震源模擬

在子程序EduS/FrontISTR/Quake中采用點(diǎn)源法將地震位錯(cuò)表示為一個(gè)集中力，為雙力偶源。斷層參數(shù)的定義如圖3所示。

斷層斷裂時(shí)，數(shù)值模擬計(jì)算采用的矩張量[16]的表達(dá)式為：

式中：M0為地震的標(biāo)量地震矩；Mij為矩張量的分量；其余參數(shù)的定義見(jiàn)圖3。

圖3 斷層參數(shù)的定義

通過(guò)單元的體積積分得到節(jié)點(diǎn)的集中力，即雙力偶源，為

(2)

式中：Pi表示節(jié)點(diǎn)力；N表示位移形狀函數(shù)；xj表示坐標(biāo)方向。

1.4 網(wǎng)格模型生成

利用非均質(zhì)的地下結(jié)構(gòu)模型J-SHIS[13]，得到圖4所示的地層結(jié)構(gòu)形狀，再采用EduS/FrontISTR軟件，可以自動(dòng)生成計(jì)算區(qū)域的有限元網(wǎng)格，最后完成如圖5所示的并行有限元計(jì)算的自動(dòng)建模。

圖4 地層形狀(富士山周邊)

圖5 EduS/FrontISTR軟件生成的網(wǎng)格模型

在EduS/FrontISTR軟件中，可以采用黏性吸收邊界[17]或無(wú)限元處理計(jì)算區(qū)域的邊界。

2 數(shù)值計(jì)算結(jié)果的理論驗(yàn)證

筆者在軟件研制過(guò)程中進(jìn)行了大量數(shù)值模擬驗(yàn)證?，F(xiàn)以地震數(shù)值模擬中典型模型為例，將研制的三維數(shù)值模擬計(jì)算結(jié)果與理論值進(jìn)行對(duì)比驗(yàn)證。

2.1 地殼形變分析的理論驗(yàn)證

數(shù)值計(jì)算區(qū)域取為長(zhǎng)200 km，寬200 km，深100 km，地質(zhì)材料為勻質(zhì)地層。計(jì)算區(qū)域采用六面體單元，網(wǎng)格單元數(shù)約為102.4萬(wàn)個(gè)，網(wǎng)格節(jié)點(diǎn)數(shù)約為106.3萬(wàn)個(gè)。地震斷層的理論模型如圖6所示，相應(yīng)的計(jì)算參數(shù)見(jiàn)表2。使用的程序?yàn)镋duS/FrontISTR/Quake，CPU為Quad Core Intel (R) Xeon X3430 2.40GHz, 采用1個(gè)Core進(jìn)行地殼變動(dòng)數(shù)值模擬的時(shí)間約為3 min，用4個(gè)Core進(jìn)行數(shù)值模擬的時(shí)間約為1.5 min。計(jì)算得到的地面位移與岡田模型[18]的理論位移的對(duì)比如圖7所示?？梢钥闯?，數(shù)值模擬的計(jì)算結(jié)果與理論值吻合得非常好。

圖6 地殼變形的斷層理論模型

斷層長(zhǎng)度L/km斷層寬度W/km傾斜角δ/(°)滑動(dòng)角λ/(°)走向角θ/(°)埋深d/km位錯(cuò)量D0/m201030909027．55

(a)地面東西向分量位移 (b)地面垂向分量位移

2.2 地震波傳播的理論驗(yàn)證

地震波傳播的理論模型如圖8所示，計(jì)算區(qū)域取為長(zhǎng)100 km，寬100 km，深50 km，斷層為水平橫向斷層，地質(zhì)材料常數(shù)為Vs=3 km/s,Vp=5 km/s，ρ=2 500 kg/m3的勻質(zhì)地層，斷層用10×10個(gè)小斷層進(jìn)行數(shù)值模擬。斷層破壞從震中(0,0,-10 km)開(kāi)始，以傳播速度vr=3 km/s同心圓的形式傳播，斷層滑動(dòng)的速度時(shí)間函數(shù)是上升時(shí)間為3 s的三角函數(shù)。地質(zhì)材料的阻尼假設(shè)為無(wú)阻尼和質(zhì)量阻尼，相應(yīng)的計(jì)算參數(shù)見(jiàn)表3。

圖8 地震波傳播的斷層理論模型(1/4模型)

斷層長(zhǎng)度L/km斷層寬度W/km傾斜角δ/(°)滑動(dòng)角λ/(°)走向角θ/(°)埋深d/km位錯(cuò)量D0/m10100090102

計(jì)算區(qū)域采用六面體單元，1/4模型的網(wǎng)格單元數(shù)約為51萬(wàn)個(gè)，網(wǎng)格節(jié)點(diǎn)數(shù)約為53萬(wàn)個(gè)。模型的東西向采用無(wú)限元消除邊界波反射的影響。與理論值對(duì)比的位置(觀測(cè)點(diǎn))取為東西向A(-2, 0, 0)，B(-1, 0, 0)，C(0, 0, 0)，D(1, 0, 0)和E(2, 0, 0)等5點(diǎn)(單位為km)。

使用的程序?yàn)镋duS/FrontISTR/Quake，CPU為Quad Core Intel (R) Xeon X3430 2.40GHz,模擬地震波傳播時(shí)間為30 s，時(shí)間步長(zhǎng)為5 ms，采用1個(gè)Core進(jìn)行地震波傳播模擬的計(jì)算時(shí)間約為20 min，用4個(gè)Core的計(jì)算時(shí)間約為10 min。

將5個(gè)觀測(cè)點(diǎn)的理論值與數(shù)值計(jì)算結(jié)果進(jìn)行對(duì)比，以驗(yàn)證數(shù)值模擬結(jié)果的精確性，如圖9所示。數(shù)值計(jì)算的速度波形、峰值和位相同理論解非常一致。數(shù)值模擬的速度波形，在11 s之后出現(xiàn)振蕩現(xiàn)象，這是由于計(jì)算區(qū)域范圍比較狹窄，出現(xiàn)邊界人工波的反射干涉所致。

(a)無(wú)阻尼 (b)質(zhì)量阻尼

對(duì)比無(wú)阻尼和質(zhì)量阻尼的速度波形，可以看出兩者的波形相同，質(zhì)量阻尼波形的位相稍微滯后，峰值是無(wú)阻尼波形的80%左右。綜合以上的結(jié)果，開(kāi)發(fā)的地殼變形和地震波傳播的數(shù)值模擬程序的計(jì)算結(jié)果是精確和可靠的，同時(shí)也驗(yàn)證了阻尼處理方法的正確性。

3 日本東北地區(qū)太平洋近海地震數(shù)值模擬

3.1 再現(xiàn)地殼變形

數(shù)值模擬計(jì)算范圍取為北緯 35°到42°，東經(jīng)138°到145°，深度80 km。再現(xiàn)地殼變形使用4片斷層構(gòu)成的計(jì)算模型[19]，采用勻質(zhì)模型計(jì)算得到的水平位移，與岡田模型及觀測(cè)值的對(duì)比如圖10所示。由圖知，勻質(zhì)模型計(jì)算結(jié)果與岡田模型及觀測(cè)結(jié)果有較好的一致性。進(jìn)一步，基于非均質(zhì)的地下結(jié)構(gòu)模型J-SHIS[13]和松原地震波速度模型[14]相結(jié)合，建立了水平方向1 km分辨率的地殼變形和地震波傳播的大型數(shù)據(jù)庫(kù)。采用非均質(zhì)的地下結(jié)構(gòu)模型計(jì)算得到東西向的應(yīng)變，如圖11所示。非均質(zhì)地下結(jié)構(gòu)模型的數(shù)值模擬結(jié)果很好地揭示了日本東北地區(qū)太平洋近海地震的斷層運(yùn)動(dòng)引起越后平原的東西向伸縮，再現(xiàn)了日本新瀉—神戶地區(qū)間的應(yīng)變集中帶，與觀測(cè)結(jié)果一致。這是國(guó)際上首次使用有限元和非均質(zhì)的地下結(jié)構(gòu)模型再現(xiàn)了應(yīng)變集中帶。

圖10 數(shù)值計(jì)算的地面水平位移與岡田模型及觀測(cè)值的對(duì)比

圖11 數(shù)值計(jì)算得到的越后平原東西向應(yīng)變分布[1]

3.2 再現(xiàn)地震波傳播過(guò)程

數(shù)值模擬計(jì)算范圍和非均質(zhì)的地下結(jié)構(gòu)模型同3.1節(jié)相同，地震斷層位錯(cuò)分布采用日本氣象廳提供的原始數(shù)據(jù)[20]，該斷層的位錯(cuò)分布是利用震源附近觀測(cè)點(diǎn)的時(shí)序擬合曲線得到，最大位錯(cuò)為37 m，斷層的范圍為475 km×175 km，由走行方向19個(gè)，傾斜方向7個(gè)，共133個(gè)小斷層構(gòu)成，如圖12所示。

圖12 斷層的位錯(cuò)分布

計(jì)算區(qū)域采用六面體和四面體單元，深度方向根據(jù)非均質(zhì)地下結(jié)構(gòu)地震波速結(jié)構(gòu)，通過(guò)四面體過(guò)渡單元(遷移層)，網(wǎng)格尺寸逐漸增大，以便節(jié)省內(nèi)存和節(jié)約計(jì)算時(shí)間。整個(gè)模擬計(jì)算區(qū)域的網(wǎng)格節(jié)點(diǎn)數(shù)約為5 000萬(wàn)個(gè)，約1.5億自由度。計(jì)算區(qū)域的網(wǎng)格如圖13所示(由于網(wǎng)格太多，未能詳細(xì)顯示)。網(wǎng)格分割過(guò)程中，每個(gè)小斷層再進(jìn)一步分割為400個(gè)更小的斷層，總共使用5萬(wàn)3 200個(gè)斷層來(lái)模擬斷層滑動(dòng)，斷層滑動(dòng)的速度時(shí)間函數(shù)采用8 s上升的三角函數(shù)來(lái)模擬。

圖13 數(shù)值模擬計(jì)算區(qū)域網(wǎng)格示意圖

模擬地震波傳播時(shí)間為300 s，時(shí)間步長(zhǎng)為25 ms，時(shí)間步為12 000，采用24個(gè)Core進(jìn)行地震波傳播模擬的計(jì)算時(shí)間約為3 h。圖14所示的計(jì)算結(jié)果清晰地顯示了地震發(fā)生后，地震波向日本全土傳播的過(guò)程。

圖14 地震波傳播過(guò)程

由圖14可以看出，數(shù)值模擬在150 s時(shí)刻地震波的峰值剛好到達(dá)福島核發(fā)電廠附近。實(shí)際上，地震發(fā)生后，第2個(gè)波峰到達(dá)日本福島核發(fā)電廠的時(shí)間約為150 s，正是該波峰引起海底地殼的變動(dòng)，產(chǎn)生的海嘯將核發(fā)電廠摧毀。數(shù)值模擬與實(shí)際發(fā)生海嘯的時(shí)間吻合的很好。這種利用地震斷層位錯(cuò)模型計(jì)算得到的地震波更接近實(shí)際情況，有利于工程結(jié)構(gòu)的災(zāi)害評(píng)估。

4 結(jié)論

筆者基于開(kāi)源并行有限元計(jì)算環(huán)境FrontISTR，研究開(kāi)發(fā)了高性能的地殼變形和地震波傳播并行計(jì)算軟件EduS/FrontISTR，同時(shí)構(gòu)建了非均質(zhì)的地下結(jié)構(gòu)數(shù)據(jù)庫(kù)。該軟件可以高效地對(duì)大規(guī)模地震數(shù)值計(jì)算進(jìn)行建模和并行計(jì)算，數(shù)值模擬結(jié)果很好地再現(xiàn)了地震后的地表觀測(cè)結(jié)果或理論結(jié)果，能夠很好地再現(xiàn)地震發(fā)生后地震波的持續(xù)傳播過(guò)程?；贓duS/FrontISTR軟件計(jì)算得到的地震波將更有實(shí)際應(yīng)用價(jià)值，將有助于地震后工程結(jié)構(gòu)的災(zāi)害評(píng)估。

日本京都大學(xué)西村卓也教授(原日本國(guó)土地理院主任研究員)為本文繪制了圖10，加工了圖11，日本計(jì)算力學(xué)研究中心吉見(jiàn)顯一朗為本文繪制了地震波傳播圖片等，在此一并致謝。

[1]NISHIMURA T, SUITO H, KOBAYASHI T, et al. Excess strain in the Echigo Plain sedimentary basin, NE Japan: evidence from coseismic deformation of the 2011 Tohoku-Oki earthquake[J].Geophysical Journal Inertnation, 2016, 205(3): 1613.

[2]日本計(jì)算工學(xué)會(huì)GreenCAE研究會(huì)共催ユーザ會(huì)資料[EB/OL].http://www.multi.k.u-tokyo.ac.jp/FrontISTR/.

[3]嚴(yán)珍珍，張懷，范湘濤，等. 日本Mw9.0 大地震激發(fā)全球地震波傳播特性的數(shù)值模擬研究與對(duì)比分析[J]. 中國(guó)科學(xué)：地球科學(xué)，2014，44(2)：259.

[4] 許才軍，何平，溫?fù)P茂，等. 日本2011 Tohoku-Oki Mw 9.0級(jí)地震的同震形變及其滑動(dòng)分布反演：GPS和InSAR約束[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，2012，37(12)：1387.

[5]張林波，鄭偉英，盧本卓，等. 并行自適應(yīng)有限元軟件平臺(tái)PHG及其應(yīng)用[J]. 中國(guó)科學(xué)：信息科學(xué)，2016，46(10)：1442.

[6]董勤喜，柴山恭.基于開(kāi)源并行計(jì)算平臺(tái)FrontISTR的非均質(zhì)地下結(jié)構(gòu)的地殼變形和地震波傳播的數(shù)值模擬[C]//日本文部科學(xué)省次世代革新數(shù)值模擬軟件的研究開(kāi)發(fā)最終成果報(bào)告會(huì).東京都:東京大學(xué)生產(chǎn)技術(shù)研究所，2013.

[7]KOMATITSCH D, TSUBOI S, CHEN J, et al. A 14.6 billion degree of freedom, 5teraflops, 2.5 terabyte earthquake simulation on the Earth simulator[C]//Proceedings of the 2003 ACM/IEEE International Conference on Supercomputing (SC03). New York: ACM, 2003:4-11.

[8]CUI Y, POYRAZ E, OLSEN K B, et al. Physics-based seismic hazard analysis on petascale heterogenous supercomputers[C]// Proceeding of 2013 ACM/IEEE International Conference on High Performance Computing, Networking, Storage and Analysis (SC13). New York: ACM, 2013:1-12.

[9]HEINECKE A, BREUER A, RETTENBERGER S, et al. Petascale high order dynamic rupture earthquake simulations on heterogenous supercomputers[C]//Proceeding of the ACM/IEEE International Conference on High Performance Computing, Networking, Storage and Analysis (SC14).Piscataway: IEEE Press, 2014:3-14.

[10]ICHIMURA T, FUJITA K, QUINAY P E B, et al.Implicit nonlinear wave simulation with 1.08T DOF and 0.270T unstructured finite elements to enhance comprehensive earthquake simulation[C]// Proceeding of the ACM/IEEE International Conference on High Performance Computing, Networking, Storage and Analysis (SC15). New York: ACM, 2015:1-12.

[11] MAZZIERI I, STUPAZZINI M, GUIDOTTI R, et al. SPEED: spectral elements in elastodynamics with discontinuous Galerkin: a non-conforming approach for 3D multi-scale problems[J].Int J for Numerical Meth in Eng, 2013, 95(12): 991.

[12]BAO H, BIELAK J, GHATTAS O, et al. Large-scale simulation of elastic wave propagation in heterogeneous media on parallel computers[J].Computer Methods in Applied Mechanics & Engineering, 1998, 152:85.

[13]J-SHIS 地震ハザードステーション[EB/OL].http://www.j-shis.bosai.go.jp/.

[14]MATSUBARA M, OBARA K,KASAHARA K. Three-dimensional P-and S-wave velocity structures beneath the Japan islands obtained by high-density seismic stations by seismic tomography[J].Tectonophysics, 2008, 454(1): 86.

[15]MELOSH H J, RAEFSKY A. A simple and efficient method or introducing faults into finite element computations[J]. Bulletin of the Seismological Society of America, 1981, 71： 1391.

[16] AKI K, RICHARDS P G. Quantitative seismology theory and methods[M].New York : W H Freeman and Company, 1980:61.

[17] LYSMER J, KUHLEMEYER R L.Finite dynamic model for infinite media[J]. ASCE Journal of Engineering Mechanics Division, 1969, 95： 859.

[18]OKADA Y. Internal deformation due to shear and tensile faults in a half-space[J]. Bulletin of the Seismological Society of America, 1992, 82： 1018.

[19]NISHIMURA T, MUNEKANE H, YARAI H. The 2011 off the Pacific coast of Tohoku earthquake and its aftershocks observed by GEONET[J]. Earth Planets Space, 2011, 63(7):631.

[20]平成23年(2011年)東北地方太平洋沖地震調(diào)查報(bào)告(日文)，第I編:日本氣象廳技術(shù)報(bào)告:第133號(hào)[R].東京都:日本氣象廳,2012.

特約專家介紹

董勤喜(1962—)，男，漢族，山東蘭陵人，中共黨員，科學(xué)技術(shù)博士。1997年博士畢業(yè)于奧地利因斯布魯克大學(xué)土木與建筑工程學(xué)院。1998年8月至2000年7月日本學(xué)術(shù)振興會(huì)(JSPS)博士后。主要從事隧道開(kāi)挖數(shù)值模擬，材料參數(shù)反演，道路及機(jī)場(chǎng)跑道設(shè)計(jì)修復(fù)軟件的研發(fā)，流固耦合并行計(jì)算的研發(fā)，高性能并行計(jì)算及大規(guī)模地震數(shù)值模擬的研究，日本下一代超級(jí)計(jì)算機(jī)京[K-computer；E級(jí)]上先進(jìn)制造及熱塑性復(fù)合材料數(shù)值模擬軟件等的研究。已在《應(yīng)用數(shù)學(xué)和力學(xué)》《日本土木學(xué)會(huì)論文集》《日本道路工程學(xué)會(huì)論文集》《J. Beijing Inst. of Techn.》《Appl. Math. Mech.》《Nonlinear Analysis-Theory, Methods & Applications》《Phys. Rev. E》《Phys. Lett. A》《J. Trauma-Injury Infection and Critical Care》《Int. J. Geomech., ASCE》《Int. J. Numer. Anal. Meth. Geomech》《Finite Elements in Analysis and Design》《Geophys. J. Int.》及《Int. J. Pavement》等重要學(xué)術(shù)刊物上發(fā)表論文60余篇。主持和參與日本鐵道運(yùn)輸基礎(chǔ)研究基金、JST自然科學(xué)基金、日本文部科學(xué)省火山研究基金、日本文部科學(xué)省面向E級(jí)(百億億次級(jí)；旗艦2020計(jì)劃)超級(jí)計(jì)算機(jī)上復(fù)合材料數(shù)值模擬研究基金資助的研究課題。