林海銘

【摘要】本文通過泥石流溝地表數(shù)據(jù)點(diǎn)在3ds Max建立山體溝壑模型，導(dǎo)出STL網(wǎng)格模型；在HyperMesh中進(jìn)行細(xì)化修改，將生成的有效STL模型導(dǎo)入離散元軟件PFC3D之中；在PFC3D中生成顆粒模型并設(shè)置顆粒及壁面相關(guān)參數(shù)，循環(huán)求解，對泥石流的啟動、堆積和流動三個過程進(jìn)行數(shù)值仿真。

【關(guān)鍵詞】離散元法；泥石流；三維流態(tài)；PFC3D

中圖分類號：K92文獻(xiàn)標(biāo)識碼A文章編號1006-0278（2015）11-137-01

泥石流是由多種因素激發(fā)含有大量泥沙、石塊等的松散土體與水體混合的特殊洪流，是介于崩塌、滑坡等塊體運(yùn)動和高含沙水流運(yùn)動之間的一系列流動過程，具有迸發(fā)突然、來勢兇猛，毀壞性強(qiáng)等特點(diǎn)。我國長期受到泥石流災(zāi)害的困擾，泥石流災(zāi)害分布范圍很廣，發(fā)生頻繁，造成巨大的損失。2010年8月7日甘肅舟曲發(fā)生的特大山洪泥石流造成1748人遇難，直接經(jīng)濟(jì)損失超過14.1多億元；2012年5月10日甘肅岷縣發(fā)生特大冰雹山洪泥石流災(zāi)害共造成49人死亡，35.8萬人受災(zāi)，直接經(jīng)濟(jì)損失超過68億元；2015年8月9日溫州市因臺風(fēng)“蘇迪羅”引發(fā)山洪泥石流造成12人死亡，5人失蹤。

泥石流運(yùn)動屬于非定常的復(fù)雜多相流流動，其研究內(nèi)容屬于多學(xué)科交叉復(fù)合。當(dāng)前對泥石流的研究注重于流動機(jī)理的理論研究，提出了賓漢體模型、膨脹流模型、粘塑流模型等泥石流運(yùn)動模型。但這些模型無法全面描述泥石流的運(yùn)動，并且很難得到解析解，因此數(shù)值計(jì)算仿真顯得很重要。本文基于離散元法對泥石流的三維流動形態(tài)進(jìn)行數(shù)值仿真，對泥石流的理論研究和實(shí)際應(yīng)用都有重要的價(jià)值和意義。

一、離散元法及PFC3D程序

離散單元法（Discrete Element Method， DEM）是一種研究非連續(xù)性顆粒物質(zhì)結(jié)構(gòu)和運(yùn)動規(guī)律的一種數(shù)值方法，由Cundall和Strack于1979年最先提出。其基本思想是把不連續(xù)體或者連續(xù)體離散為具有一定物理意義的獨(dú)立“微元”或“粒子”，相鄰單元之間存在一種或幾種作用力，單元的運(yùn)動受牛頓第二運(yùn)動定律支配，得到各個離散單元的運(yùn)動方程，用時(shí)步迭代的方法求解各單元的運(yùn)動方程，更新單元節(jié)點(diǎn)的速度、位移等物理量，進(jìn)而得到整體模型的形變和運(yùn)動形態(tài)。

PFC3D是由1979年Cundall和Strack推出的三維圓球程序TRUBAL發(fā)展而來，通過離散元方法來模擬球形顆粒介質(zhì)的運(yùn)動及其相互作用，既可直接模擬球形顆粒的運(yùn)動與相互作用問題，也可以通過兩個或多個顆粒與其直接相鄰的顆粒連接形成任意形狀的組合體來模擬塊體結(jié)構(gòu)問題。顆粒單元的直徑可以是一定的，也可按高斯分布規(guī)律分布，單元生成器根據(jù)所描述的單元分布規(guī)律自動進(jìn)行統(tǒng)計(jì)并生成單元。PFC3D既可解決靜態(tài)問題也可解決動態(tài)問題；既可用于參數(shù)預(yù)測，也可用于在原始資料詳細(xì)情況下的實(shí)際模擬；還可以用PFC3D模擬試驗(yàn)代替室內(nèi)試驗(yàn)。PFC3D在模擬顆粒間的相互作用、大變形、斷裂等方面應(yīng)用廣泛。

二、泥石流模型的建立

以甘肅省永靖縣境內(nèi)某一條泥石流溝為模擬研究對象，根據(jù)文獻(xiàn)[3]中提供的泥石流溝地表數(shù)據(jù)點(diǎn)在3dsMAX中建立樣條曲線，再通過樣條曲線生成較為光滑的溝壑曲面，最后導(dǎo)出PFC3D支持的STL三角網(wǎng)格模型（圖1），限于計(jì)算規(guī)模，僅取溝壑段長100m，溝壑寬度為20m。3dsMAX生成的網(wǎng)格較為粗糙，不能滿足計(jì)算要求，再到HyperMesh中細(xì)化網(wǎng)格（圖2），通過節(jié)點(diǎn)耦合排除網(wǎng)格模型中存在的裂隙、孤立邊等幾何缺陷保證STL模型的有效性，并統(tǒng)一三角網(wǎng)格的方向確保能正確計(jì)算顆粒跟網(wǎng)格之間的接觸問題。最后編寫FISH函數(shù)將細(xì)化后的、有效的STL模型讀入PFC3D程序中（圖3）。

編制顆粒生成FISH子程序，生成大小不一的顆粒，并根據(jù)文獻(xiàn)[4]中直剪試驗(yàn)結(jié)果設(shè)置顆粒的摩擦系數(shù)為0.25，密度2.147，法向接觸黏接強(qiáng)度、剪切黏接強(qiáng)度均為6.7Kpa。假設(shè)墻體靜止不動，并賦予其相應(yīng)的接觸剛度和剪切剛度。

三、泥石流流動形態(tài)仿真

顆粒在重力的作用下下落，并受建立在山體上的一個壁面約束，堆積在山體上，不沿溝壑流下，模擬泥石流爆發(fā)前，顆粒的堆積和所建立的攔擋壩對泥石流的攔擋作用。

刪除溝壑模型上加的約束壁面（相當(dāng)于攔擋壩在地震、暴雨沖刷等外力的影響下發(fā)生破壞，失去攔擋作用），顆粒流啟動，開始在溝壑中流動。顆粒顏色代表顆粒速度，從深藍(lán)色到紅色，速度從0 m/s逐一遞增到14 m/s。

由結(jié)果可見，初始3s時(shí)，開始流動的顆粒流頂部的顆粒勢能大，并且受到約束較少，釋放能量快，速度較大；圖a、b中，同一高度，不同尺寸的顆粒速度基本一致，但圖c-e中，同一高度，小顆粒速度明顯高于大顆粒速度，且跑在最前線的是小顆粒，這是由于小顆粒的相對接觸面積小，能量損耗少；顆粒從較為陡峭的坡面流下來后，幾乎在水平面上前行，由于能量的耗散，速度逐漸減慢；最終有少數(shù)顆粒停留在第一個臺階處。

四、結(jié)語

本文將泥石流體簡化為相互之間有粘接效應(yīng)的顆粒，通過PFC3D模擬泥石流的三維流態(tài)。本文先在3ds Max中建立山體溝壑STL模型并導(dǎo)入HyperMesh中細(xì)化，再導(dǎo)入到PFC3D中，然后生成大小不同的顆粒模型，最后模擬泥石流的堆積、啟動和流動等過程。本文采用離散元分析特定情況下的泥石流流態(tài)具有一定的研究價(jià)值，是后續(xù)模擬泥石流防護(hù)措施的基礎(chǔ)。

參考文獻(xiàn)：

[1]袁道先，謝云鶴.巖溶與人類生存、環(huán)境、資源和災(zāi)害[M].桂林：廣西師范大學(xué)出版社，1996.

[2]孫其誠，王光謙.顆粒物質(zhì)力學(xué)導(dǎo)論[M].北京：科學(xué)出版社，200

[3]龔成勇.粘性泥石流運(yùn)動機(jī)理及數(shù)值模擬研究[D].蘭州理工大學(xué)，2009

[4]胡明鑒，汪稔，陳中學(xué)，王志兵.泥石流啟動過程PFC數(shù)值模擬[J].巖土力學(xué)，2010（8）.