黃 雨，趙紫陽

(1.同濟大學(xué) 巖土及地下工程教育部重點實驗室，上海 200092；2.同濟大學(xué) 土木工程學(xué)院地下建筑與工程系，上海 200092)

人類在探索地外行星過程中發(fā)現(xiàn)諸多外太空微重力環(huán)境下發(fā)生的大型巖質(zhì)滑坡碎屑流。如Apollo 17登月艙觸發(fā)的滑坡巖屑流[1]面積27 km2，在平坦月球表面運動5 km，最大高度落差2 km；火星表面3個大的滑坡寬度范圍近120，最大豎直落差達914 km；火星Gangis Chasma南坡的滑坡堆積物[2]體積100×109 m3，以超過100 km/h速度沿槽狀谷地運動60 km。圖1為NASA伽利略計劃公布的土衛(wèi)八Cassini Regio區(qū)域大型滑坡。坡體由較細土質(zhì)顆粒組成，滑坡距離達600 km[3]。

由此看出，因重力場的變化，微重力下發(fā)生滑坡過程與常規(guī)重力場發(fā)生滑坡地質(zhì)災(zāi)害迥然不同，其超強的動力特性、大規(guī)模高速遠程運動、獨特的失穩(wěn)機理等均遠超現(xiàn)有知識范疇。

圖1 土衛(wèi)八Cassini Regio區(qū)域滑坡[3]

中科院空間科技發(fā)展路線圖提出，2030年實現(xiàn)首次載人登月，2040年建立短期有人值守的月球基地，2050年實現(xiàn)首次載人登陸火星[4]。為此，開展微重力環(huán)境下顆粒物質(zhì)動力學(xué)行為研究，具有明確的國家重大需求背景，亦會促進工程地質(zhì)學(xué)科在空間微重力科學(xué)領(lǐng)域發(fā)展。鑒于高速滑坡等地質(zhì)災(zāi)害給人類生存造成的嚴重后果，需對微(弱)重力環(huán)境下大規(guī)模巖土流動破壞予以關(guān)注、充分考慮。研究微重力下顆粒材料的動力學(xué)行為機理及特點具有極重要現(xiàn)實意義及科學(xué)價值。

顆粒物質(zhì)運動為觸發(fā)滑坡、泥石流等自然災(zāi)害主要原因之一；但早期的顆粒物質(zhì)研究更多關(guān)注顆粒物質(zhì)在正常重力場中的動力學(xué)行為，對顆粒物質(zhì)運動機理研究大多局限于探討顆粒間力的作用及顆粒流動力學(xué)模型，對空間微重力環(huán)境下顆粒物質(zhì)振動、碰撞、流動等動力學(xué)行為機制研究尚少。微重力環(huán)境指重力大小較地球重力加速度980 cm/s2低得多的環(huán)境。顆粒物質(zhì)為具有內(nèi)在聯(lián)系的復(fù)雜體系，整個顆粒介質(zhì)在外力或內(nèi)部應(yīng)力發(fā)生變化時產(chǎn)生流動，表現(xiàn)出類似流體性質(zhì)，形成顆粒流。研究發(fā)現(xiàn)，正常重力與微重力環(huán)境下所得顆粒形態(tài)特征大體相同，形態(tài)呈亞微米球形、橢球形顆粒分布，大小明顯不同，微重力下顆粒是正常重力的兩倍大，顆粒大小在微重力環(huán)境下較正常重力場中變化范圍大得多[5]。微重力環(huán)境中顆粒以懸浮形式存在，顆粒物質(zhì)振動、碰撞、流動等動力學(xué)行為因重力水平極大減小產(chǎn)生明顯變化，使顆粒體系發(fā)生失穩(wěn)的誘因、運動過程及特點與地球重力場中迥然不同。因此，基于微重力獨特的環(huán)境條件，由動力學(xué)角度出發(fā)分析顆粒物質(zhì)振動、碰撞、流動等運動行為是顆粒物質(zhì)研究的重要方向。研究問題有：①顆粒物質(zhì)在振動、碰撞、流動時的基本性質(zhì)、特點；②影響顆粒物質(zhì)振動、碰撞、流動因素；③顆粒物質(zhì)振動、碰撞、流動可能誘發(fā)的后果。此研究有助于了解微重力下顆粒動力學(xué)行為機制，為人類進行空間探索提供有價值的理論依據(jù)。為此，本文對微重力下顆粒物質(zhì)振動、碰撞、流動認識進行綜述，分析總結(jié)顆粒動力學(xué)機制及特點。

1 振 動

已有諸多對在微重力環(huán)境下以振動作為驅(qū)動源的顆粒運動行為規(guī)律進行的研究。如美國宇航局在太空探針火箭上進行的顆粒物質(zhì)力學(xué)實驗，用空間微重力環(huán)境中振動干、濕顆粒模擬地震條件，獲得動力學(xué)過程圖像[6]。Pierrer等[7-8]進行的火箭實驗中發(fā)現(xiàn)在粒子數(shù)密度足夠大區(qū)域或容器中顆粒按晶格緊密排列，容器中只有兩個顆粒時會出現(xiàn)周期性“超聲壓縮”最終達到共振狀態(tài)。中科學(xué)院物理研究所在八號育種衛(wèi)星上首次開展的對微重力下顆粒物質(zhì)在小幅振動驅(qū)動下運動行為研究[9]發(fā)現(xiàn)氣體分布狀態(tài)的顆粒速度分布滿足非高斯指數(shù)分布規(guī)律，非通常的分子氣體麥克斯韋分布。Garrabos等[10-12]在空中客機A300-0g拋物線飛行計劃實驗中研究氣體分布狀態(tài)的顆粒受振動時對容器壁的沖量分布，顆粒數(shù)N=36,48的實驗結(jié)果見圖2，其中PDF指幾率分布函數(shù)，A為振幅，I為撞擊強度，f為容器振動頻率。由2圖看出，在不同實驗條件下撞擊強度I標度化后的分布函數(shù)p(I)滿足指數(shù)分布[12]。

圖2 撞擊強度I的幾率分布函數(shù)[12]

顆粒體系為多尺度耗散體系，因密度或體積變化存在相轉(zhuǎn)變過程[6]。Opsomer等[13]利用Mini-Texus5試驗參數(shù)在三維模擬中發(fā)現(xiàn)微重力下振動顆粒體系存在4種不同動態(tài)機制：Ⅰ氣態(tài)、Ⅱ局部聚集、Ⅲ完全聚集、Ⅳ彈跳聚集。Ⅰ～Ⅳ動態(tài)機制相變圖見圖3，其中，r為顆粒半徑，η為顆粒密度。由圖3看出，r，η不同取值顆粒體系分別對應(yīng)4種不同狀態(tài)，振動會使r，η變化，導(dǎo)致顆粒體系發(fā)生相變。

Saadatmand等[14]將一球形顆粒置于微重力狀態(tài)下充滿液體的容器中，平行于z軸方向以一定振幅、頻率施加振動，由于顆粒與周圍液體密度不同，顆粒被迫產(chǎn)生類似鄰近容器壁振動。實驗結(jié)果表明，液體黏稠度為1 Pa·s時，顆粒與容器壁間壓力小于液體內(nèi)部其它區(qū)域壓力，與最近容器壁間流體速度最高，且在其與最近容器壁間空隙中振蕩作往返運動。可見，顆粒體系的穩(wěn)定性與激發(fā)振動大小密切相關(guān)。受到振動時體系內(nèi)部顆粒分布及密度會發(fā)生變化，從而改變顆粒間摩擦力，力鏈網(wǎng)絡(luò)中力的分布亦會變化。力鏈上任何局部或微小位置變動均可能引起顆粒體系較大變化，造成“螻蟻之穴，潰堤千里”結(jié)果[15]。即空間環(huán)境中較小振動均可能觸發(fā)大型滑坡。

由此可知，微重力下顆粒物質(zhì)在振動驅(qū)動下的基本性質(zhì)為：① 顆粒物質(zhì)數(shù)量較少時，顆粒之間發(fā)生碰撞概率較小，顆粒體系在足夠大區(qū)域中可達共振狀態(tài)；② 顆粒速度分布為非高斯指數(shù)分布；③振動顆粒體系隨顆粒半徑、密度不同存在4種相變。

2 碰 撞

空間微重力環(huán)境下以懸浮形式存在的顆粒由于非彈性碰撞致能量損耗，顆粒會產(chǎn)生聚集或粘結(jié)，而氣體分子系統(tǒng)則不會[6]。微重力下微米大小顆粒碰撞速度足夠小，碰撞后顆粒相互粘結(jié)形成不規(guī)則顆粒集合體，粘結(jié)后兩顆粒會繞粘結(jié)體質(zhì)心旋轉(zhuǎn)，集合體隨顆粒的逐漸變大，會更快從周圍氣體中分離[16]，見圖4。由圖4可知，兩灰塵顆粒以9 mm/s速度發(fā)生碰撞后粘結(jié)在一起并繞其中心旋轉(zhuǎn)。

圖4 兩灰塵顆粒以9 mm/s速度碰撞的粘結(jié)次序圖[16]

圖5 單個灰塵顆粒以17 mm/s速度與2個灰塵顆粒碰撞的粘結(jié)次序圖[16]

單個與2個灰塵顆粒發(fā)生碰撞后粘結(jié)一起的次序見圖5。Weidling等[16]在微重力實驗中用新裝置研究大小不同灰塵顆粒集合體間的自由碰撞，推斷顆粒碰撞后由粘結(jié)至彈跳的分離變化并非突然，存在速度相差兩數(shù)量級的過渡帶，在過渡帶中顆粒碰撞后可能發(fā)生粘結(jié)或彈跳。發(fā)生粘結(jié)概率會隨碰撞速度減小而增大。顆粒越小，碰撞時發(fā)生粘結(jié)概率越大。

單顆粒與集合體間發(fā)生碰撞后較單顆粒間更易發(fā)生粘結(jié)。顆粒碰撞后粘結(jié)概率隨碰撞速度的減小而增大，隨顆粒尺寸的減小而增大[17]?；覊m顆粒集合體碰撞后粘結(jié)會導(dǎo)致更大結(jié)構(gòu)的形成[18]。Blum[19]在落塔中進行多顆粒體系以超低速碰撞試驗發(fā)現(xiàn)，碰撞速度大于一定值時碰撞顆粒碰撞后會粘結(jié)在一起，入射顆粒貫入深度與碰撞速度有關(guān)，深度大于其半徑時由目標顆粒彈跳進入射軌道。顆粒碰撞時粘結(jié)、彈跳可同時發(fā)生。由此可見，空間微重力環(huán)境下顆粒間碰撞后是否粘結(jié)取決于碰撞時顆粒速度大小，顆粒碰撞時粘結(jié)概率隨碰撞速度的減小而增大，但速度大于某臨界值時顆粒碰撞后又黏結(jié)在一起。在速度相差兩個數(shù)量級過渡帶中顆粒碰撞后可能粘結(jié)亦可能分離。

3 流 動

顆粒流研究開始于20世紀50年代。Bagnold[20]通過對蠟球在甘油-酒精-水溶液中進行同軸剪切流變儀實驗研究，揭示出顆粒的摩擦應(yīng)力本構(gòu)關(guān)系。Ogawa[21]引入“顆粒溫度”概念，定義平動顆粒溫度為T=/3，其中無規(guī)速度c=(v-u)，u=，尖括號表示整體平均。Cundall等開發(fā)出第1個研究顆粒介質(zhì)力學(xué)行為的二維離散元模擬程序，獲得與動光彈實驗極吻合結(jié)果，成為離散元模擬顆粒運動開端。Spencer[22]發(fā)展了以Mohr-Coulumb條件為基礎(chǔ)的“雙剪切”理論并用于顆粒流動研究。實驗證明該模型較有效。有研究亦觀察到顆粒流自發(fā)無序混合、指進、非均勻尺寸顆粒偏聚等[23-26]。

Cundall等首創(chuàng)的離散元(DEM)成為重要顆粒流動數(shù)值模擬方法，即將介質(zhì)離散成獨立元或粒子，單個顆粒為1個元，相鄰間存在幾種相互作用。自然界中多數(shù)顆粒形狀不規(guī)則，大部分呈橢圓形，用球形近似不會產(chǎn)生明顯誤差，故圓球顆粒間相互作用成為離散顆粒動力學(xué)基礎(chǔ)[27]。孫其誠等[28]據(jù)對顆粒接觸作用方式不同分為3種模型：① 基于接觸力學(xué)的離散模型；② 軟球模型；③硬球模型。其中①為基礎(chǔ)，②、③為基于接觸力學(xué)離散模型的簡化。工程中將以軟球模型為主。

Hofmeister等[29]用玄武巖球在真空條件下準二維沙漏裝置中模擬顆粒物質(zhì)的流動行為，實驗測得不同重力水平下顆粒流動體積流量、休止角、摩擦角、速度、顆粒流厚度變化見圖6。據(jù)理論，顆粒流體積流量為重力的平方根[30]，但實驗結(jié)果表明指數(shù)為0.60，較理論預(yù)測值偏大。顆粒物質(zhì)流動速度隨深度的增加而減??；重力水平越高顆粒流動速度減小越快；顆粒流厚度隨重力水平增大而增大。由圖6可知，摩擦角隨重力的變化幅度較休止角大。

Huang等[31]通過利用微重力落塔設(shè)備分別采集微重力狀態(tài)及1 g重力狀態(tài)下顆粒材料沿不同傾角斜坡流動的相關(guān)數(shù)據(jù)發(fā)現(xiàn)，模型箱傾角對巖土顆粒材料流動特性有較大影響。小傾角時僅部分砂土流動，且持時較短；大傾角時全部砂土流動，持時較長。給定初始速度時巖土顆粒材料在微重力環(huán)境中流動持時更長、范圍更大。進入微重力環(huán)境后原未獲得初速度的巖土顆粒在整個實驗中處于靜止狀態(tài)。已獲得初速度的巖土顆粒則在減弱的重力場中繼續(xù)流動，且持時更長，流動范圍更廣。微重力環(huán)境下運動中的巖土顆粒能量衰減較慢，流動至模型箱尾部時能順模型箱壁面改變方向繼續(xù)流動，最終在殘余重力及顆粒間相互作用下流動停止。30°傾角下常規(guī)重力環(huán)境與微重力環(huán)境巖土顆粒流動構(gòu)型對比見圖7。

圖6 休止角、摩擦角與重力水平函數(shù)關(guān)系[29]

綜上可知，顆粒物質(zhì)流動時隨重力水平增大顆粒體積流量呈指數(shù)增加，速度隨深度增加下降越快，表面流厚度呈增大趨勢；摩擦角、休止角與重力水平大小相關(guān)。微重力環(huán)境下顆粒沿斜坡流動隨斜坡傾角的變化而不同，運動的巖土顆粒微重力下流動時間更長。由于重力極大減小，顆粒物質(zhì)流動特點被放大。

4 結(jié) 論

顆粒物質(zhì)空間動力學(xué)行為特點、機理、動力學(xué)模型、觸發(fā)因素及可能產(chǎn)生的后果等直接關(guān)系到未來對空間的探索;微重力條件提供了直接觀察顆粒運動規(guī)律的環(huán)境，使對顆粒體系的基本動力學(xué)機制研究成為可能。通過綜述，結(jié)論如下：

(1) 微重力環(huán)境中顆粒發(fā)生振動、碰撞、流動條件、觸發(fā)因素、機制、特點及可能產(chǎn)生的后果需綜合考慮巖土顆粒性質(zhì)(大小、形狀、材料)、荷載特征、外界環(huán)境等多種因素。

(2) 顆粒振動、碰撞、流動三種動力學(xué)行為間相互作用，可能存在相互關(guān)系，且直接影響顆粒體系穩(wěn)定性，故需研究顆粒動力學(xué)行為與穩(wěn)定性關(guān)系。

(3) 微重力條件下顆粒物質(zhì)動力學(xué)模型相對滯后，已有適用于重力場顆粒運動的多種模型是否適用于微重力下顆粒運動尚不清楚。因此需探尋適用于空間微重力環(huán)境中多顆粒相互作用及非均勻顆粒動力學(xué)特征的新理論、新模型，指導(dǎo)登月、火星探索。

[1] Masursky H, Colton G W, El-Baz F. Apollo over the moon: a view from orbit[M]. Washington, D.C: National Aeronautics and Space Administration, 1978,SP-362.

[2] Malin space science systems. mars orbiter sees landslide [EB/OL]. http://spaceflightnow.com/news/n0111/05marslandslide/,2001.11.05.

[3] NASA. Giant landslide on iapetus [EB/OL]. http://photojournal.jpl.nasa.gov/catalog/PIA06171,2005.01.07.

[4] 人民網(wǎng). 中科院發(fā)布未來五十年中國科技發(fā)展路線圖[EB/OL]. http://scitech.people.com.cn/GB/9451859.html, 2009.06.11.

[5] 徐宇. 煙塵顆?；咎匦缘臏y試與分析[J]. 化學(xué)推進劑與高分子材料, 2002 (1): 18-20.

XU Yu. Test and analysis of basic characteristics of soot particles [J]. Chemical propellants and polymeric materials, 2002 (1): 18-20.

[6] 厚美瑛. 空間環(huán)境顆粒物質(zhì)運動行為的研究[J]. 物理, 2008, 37(10): 729-732.

HOU Mei-ying. Research on the motor behavior of particulate matter in the space [J]. Physical, 2008, 37(10): 729-732.

[7] Evesque P, Adjémian F. Stress fluctuations and macroscopic stick-slip in granular materials[J]. The European Physical Journal E, 2002, 9(3): 253-259.

[8] Evesqu P. Granular materials : fundamentals and applications[M]. Ed by S. Antony, Royal Society of Chemistry, 2004:29-62.

[9] Hou M, Liu R, Zhai G, et al. Velocity distribution of vibration-driven granular gas in Knudsen regime in microgravity[J]. Microgravity Science and Technology, 2008, 20(2): 73-80.

[10] Garrabos Y, Evesque P, Palencia F, et al. Coherent behavior of balls in a vibrated box[EB/OL]. http://arxiv.org/pdf/cond-mat/0611613v1.pdf, 2006.11.23.

[11] Leconte M, Garrabos Y, Palencia F, et al. Inelastic ball-plane impact: an accurate way to measure the normal restitution coefficient[J]. Applied Physics Letters, 2006, 89(24): 1-3.

[12] Falcon E, Aumatre S, Evesque P, et al. Collision statistics in a dilute granular gas fluidized by vibrations in low gravity[J]. Europhysics Letters, 2006, 74(5): 830-836.

[13] Opsomer E, Ludewig F, Vandewalle N. Phase transitions in vibrated granular systems in microgravity[J]. Physical Review E, 2011, 84(5): 1-5.

[14] Saadatmand M, Kawaji M, Hu H H. Vibration-Induced attraction of a particle towards a wall in microgravity-the mechanism of attraction force[J]. Microgravity Science and Technology, 2012, 24(1): 53-64.

[15] 孫其誠, 王光謙, 胡凱衡. 顆粒物質(zhì)力學(xué)幾個關(guān)鍵問題的思考[J]. 自然科學(xué)進展, 2008, 18(10): 1104-1110.

SUN Qi-cheng, WANG Guang-qian, HU Kai-heng. Several key issues on granular materials mechanics[J]. Progress in Natural Science, 2008, 18(10): 1104-1110.

[16] Weidling R, Güttler C, Blum J. Free collisions in a microgravity many-particle experiment. I. dust aggregate sticking at low velocities[J]. Icarus, 2012, 218(1): 688-700.

[17] Beitz E, Güttler C, Weidling R, et al. Free collisions in a microgravity many-particle experiment-II: the collision dynamics of dust-coated chondrules[J]. Icarus, 2012, 218(1): 701-706.

[18] Kothe S, Blum F, Weidling R, et al. Free collisions in a microgravity many-particle experiment.III:the collision behavior of sub-millimeter-sized dust aggregates[J]. Icarus, 2013, 225(1): 75-85.

[19] Blum J. Astrophysical microgravity experiments with dust particles[J]. Microgravity Science and Technology, 2010, 22(4): 517-527.

[20] Bagnold R A. Experiments on a gravity-free dispersion of large solid particles in a Newtonian fluid under shear[J]. Proc R Soc Lond A, 1954, 225(1160): 49-63.

[21] Ogawa S. Multitempreture theory of granular materials[C].// Cowin S C, Satake M. Proceedings of continuum mechanics and statistical approaches to the mechanics of granular media of US-Japan[A].Tokyo: Gukunjuntsu bunken Fukyukai, 1978: 208-217.

[22] Spencer A J M. Plane deformations of ideal granular materials[J]. Journal of the Mechanics and Physics of Solids, 1983, 31(2): 105-122.

[23] Shinbrot T, Alexander A, Muzzio F J. Spontaneous chaotic granular mixing[J]. Nature, 1999, 397(6721): 675-678.

[24] Pouliquen O, Delour J, Savage S B. Fingering in granular flows[J]. Nature, 1997, 386(24): 816-817.

[25] Maske H A,Havlin S,King P R,et al. Spontaneous stratification in granular mixtures[J]. Nature,1997,386(27): 379-381.

[26] Baxter J, Tüzün U, Heyes D, et al. Stratification in poured granular heaps[J]. Nature, 1998, 391(6663): 136-136.

[27] 吳清松, 胡茂彬. 顆粒流的動力學(xué)模型和實驗研究進展[J]. 力學(xué)進展, 2002, 32(2): 250-258.

WU Qing-song,HU Mao-bin. Progress on kinetic model and experimental research of granular flow[J]. Advances in mechanics, 2002, 32(2): 250-258.

[28] 孫其誠, 王光謙. 顆粒流動力學(xué)及其離散模型評述[J]. 力學(xué)進展. 2008, 38(1): 87-100.

SUN Qi-cheng, WANG Guang-qian. Commentary on granular flow dynamics and discrete model[J]. Advances in mechanics, 2008, 38(1): 87-100.

[29] Hofmeister P G, Blum J, Hei?elmann D. The flow of granular matter under reduced-gravity conditions[C].// Nakagawa M, Luding S. Proceedings of the 6th International Conference on Micromechanics of Granular Media[A]. New York: AIP Publishing, 2009:71-74.

[30] Hagen G. Druckund bewegung des trockenen sandes[J]. Bericht über Diezur Bekanntmachung Geeigneten Verhandlungen der K?niglich Preu?ischen Akademie der Wissenschaften zu Berlin, 1852, 35(2): 35-42.

[31] Huang Y, Mao W. First results derived from a drop-tower testing system for granular flow in a microgravity environment[J]. Landslides, 2013, 10(4):493-501.