孫海忠

（1.同濟大學 土木工程學院，上海 200092；2.上海巖土工程勘察設計研究院有限公司，上海 200032）

剪脹性是顆粒材料不同于一般連續(xù)體材料的重要標志.剪脹性是指土體在剪應力作用下體積發(fā)生膨脹的現(xiàn)象，無粘性粒狀土在密實狀態(tài)下易表現(xiàn)出剪脹性.粗粒土為單粒結構，顆粒之間存在咬合作用，在密度較大的條件下受剪切時，位于剪切面附近的顆粒必須以某種方式滑動、轉動或折斷以克服顆粒之間的咬合，尋求剪應變空間，致使體積發(fā)生變化.

土的剪脹性很早被發(fā)現(xiàn)并有了大量的研究工作，不少土的本構模型能夠較好地模擬土的剪脹性，但未考慮實際的土體變形（如顆粒破碎等）是如何發(fā)生的，只是一種宏觀的唯像論方法.一些學者對粗粒土的剪脹性進行了探索，Rowe［1］系統(tǒng)地研究了顆粒材料的變形機理，提出了最小能比原理，推導出了著名的Rowe剪脹模型.Ueng等［2］基于能量守恒原理和Rowe提出的最小能比原理從細觀角度建立了考慮顆粒破碎的剪脹方程，在Rowe剪脹方程的基礎上增加了1項顆粒破碎耗能因子，并通過試驗對耗能因子進行了探討，認為破碎耗能因子與破碎顆粒的表面積成正比.值得一提的是，以目前的試驗條件很難對粗粒土的顆粒破碎表面積進行準確度量，因此，將該剪脹方程應用于實際模型中非常復雜，也不太現(xiàn)實；文獻［3］也在Rowe剪脹方程的基礎上增加了1項顆粒破碎耗能因子，通過粗粒土的大三軸試驗研究提出了試樣破壞時破碎耗能因子與圍壓成線性關系.然而在試驗過程中不同圍壓下的顆粒破碎耗能因子會隨著軸向應變的增加而不斷增大，因此很難對圍壓進行歸一化.如何準確地描述不同圍壓下顆粒破碎耗能因子在試驗過程中的變化趨勢，并將該考慮顆粒破碎的剪脹方程能夠得以合理應用是本文探討的關鍵.本文對Ueng等［2］建立的基于細觀理論的粗粒土剪脹方程進行了改進，對不同圍壓下顆粒破碎耗能因子在試驗過程中的變化趨勢進行了歸一化，并建立了相應的彈塑性本構模型.

1 基于細觀理論的剪脹方程

Rowe［1］根據(jù)緊密粒狀土體的受力平衡條件和變形相容條件考慮變形總是沿著能量耗散最小的路徑進行推導得出著名的Rowe剪脹方程［4］，即

式中：σ1和σ3分別為最大和最小主應力；εv，p和ε1，p分別為塑性體積應變和軸向應變；φμ為顆粒之間的摩擦角.Rowe剪脹方程考慮了土體變形過程中的摩擦和顆粒滑移造成的剪脹耗能.緊密粗粒土在受力變形過程中除了摩擦與顆?；圃斐傻募裘浐哪芡?，顆粒的破碎造成的耗能也不可避免，特別是在高圍壓作用下，土體所受的約束作用增強，顆?；频募裘浶魅酰w粒破碎量變大，若不考慮顆粒破碎耗能效應顯然是不合理的［4］.

在大三軸應力條件下緊密粒狀土體受力及變形如圖1所示，圖中，βi為顆粒之間的接觸面與鉛垂方向的夾角；F1i，F(xiàn)3i分別為顆粒之間接觸力在鉛垂方向和水平方向的分量；Si，Ni分別為顆粒之間接觸力沿顆粒作用面切向和法向分量；ui為顆粒之間的相對位移，yi，xi分別為顆粒之間相對位移在鉛垂方向和水平方向的分量.文獻［2-3］在Rowe剪脹方程的基礎上從細觀角度考慮破碎能耗的修正，得到了以下修正的剪脹方程：

式中：q＝σ1-σ3；p＝（σ1＋2σ3）／3；dεv＝dε1＋2dε3，ε3為圍壓作用產生的應變；WE，p為單位體積粒狀土體的顆粒破碎耗能（有效塑性功）.

在大三軸試驗中由塑性軸向應變、塑性體應變、廣義塑性剪應變三者之間的相互關系得到

式中，εs，p為廣義塑性剪應變.由硬化內變量 M＝得到

圖1 緊密粒狀土體的受力和變形［1］Fig.1 Loading and deformation of dense granular aggregates［1］

為了使上述剪脹性方程能夠得到應用，這里關鍵在于dWE，p／dε1的處理.文獻［2］認為 WE，p與破碎顆粒的表面積成正比，但實際上顆粒破碎表面積很難量測.文獻［3］則認為WE，p與圍壓σ3成正比，但WE，p隨軸向應變ε1的增加不斷增大，因此很難直接對圍壓進行歸一化.

對于WE，p的計算采用 Chavez等［5］提出的方法

式中：WR，p為顆粒重新排列所消耗的塑性功，一般假設為某一常數(shù)，該值較小，可以近似為零；WT，p為粗粒土在加載過程中所消耗的總功，即WT，p＝∫dWT，p，因此，dWE，p表示為

文獻［6］開展了圓形、菱形顆粒形狀的粗粒土大三軸排水應力-應變試驗，如圖2所示.

按照式（8）將顆粒破碎所消耗的有效塑性功WE，p分為2個部分進行計算，即軸向應力做功∫σ1dε1，p和圍壓做功∫2σ3dε3，p，相應計算結果如圖3和4，將2個部分塑性功疊加得出WE，p隨軸向應變ε1的變化趨勢，如圖5a和5c所示.經(jīng)過數(shù)值分析，可以得到dWE，p／dε1隨軸向應變ε1的變化關系，如圖5b和5d所示，而ε1與WE，p存在著一一對應關系，因此，dWE，p／dε1可以表示為關 于 有效塑性功WE，p的函數(shù)，即dWE，p／dε1＝f（WE，p）.

通過相關計算，得到了2組土樣剪切時dWE，p／dε1與有效塑性功WE，p的變化關系，如圖6a和6c所示.可以發(fā)現(xiàn)，兩者幾乎呈雙曲線關系，這里建議采用如下關系式進行表述：

式中：a為dWE，p／dε1與WE，p關系曲線初始切線模量的倒數(shù)；b為其關系曲線強度極限值的倒數(shù).用式（9）對2組結果進行最小二乘法非線性擬合，得到如圖6a和6c中虛線所示.相對應的擬合參數(shù)a，b見表1.

圖4 菱形粗粒土塑性功Fig.4 Plastic work associated with angular coarse granular aggregates sample

從表1可見，離差系數(shù)均小于0.05，符合擬合精度要求.參數(shù)a幾乎沒有變化，可以取4種圍壓下的平均值作為模型參數(shù)，圓形和菱形粗粒土的土樣參數(shù)a分別為0.0099和0.0072.而參數(shù)b則隨圍壓的增加而減小，如圖6b和6d所示，兩者幾乎成線性關系，可以表示為

式中：k為比例參數(shù)；c為擬合直線截距.由圖6b和6d的擬合結果得出圓形和菱形粗粒土土樣的k值和c值分別為-0.29，0.46和-0.61，0.92.

將式（9）和式（10）代入式（6），得到最終的剪脹方程

表1 粗粒土回歸系數(shù)Tab.1 Regression coefficient of rounded coarse aggregates

2 本構模型的建立與驗證

2.1 模型的建立

這里通過將前面推導出的剪脹方程代入到相關本構模型來對室內顆粒破碎試驗進行數(shù)值分析，從而間接地驗證該剪脹方程正確與否.

聯(lián)立式（12）和式（13）得到

將式（14）代入式（6）得到塑性勢方程

式中：Mco和Mcres分別為試樣在剪切過程中隨著有效塑性功WE，p的增加而表現(xiàn)出的初始狀態(tài)應力比和殘余狀態(tài)應力比；e為孔隙比；α為有效塑性功參數(shù).

為了應用該本構模型模擬粗粒土在不同加載水平下的應力應變關系和顆粒破碎程度，用FORTRAN語言編制相應的點積分程序并通過顯式積分格式求得各計算步應力增量、塑性剪應變增量、塑性體應變增量，同時進行逐步累加得到應力和應變全量.在對 WE，p、臨界狀態(tài)應力比 Mc（WE，p）數(shù)值計算時其迭代過程為（WE，p）n＋1＝（WE，p）n＋dWE，P，dWE，p＝ qdεs，p＋ pdεv，p， ［Mc（WE，p）］n＋1＝［Mc（WE，p）］n＋d Mc（WE，p），d Mc（WE，p）可通過 對式（16）微分得到

2.2 模型驗證

為了驗證模型的適用性，擬針對文獻［6］開展的圓形（包括亞圓形）和菱形（包括亞菱形）2種顆粒形狀的粗粒土大三軸顆粒破碎試驗結果（2種粗粒土的臨界狀態(tài)線在e-ln p平面中均為直線型），用該模型對粗粒土在不同圍壓下的強度和變形特性進行數(shù)值模擬.首先分別將這2種粗粒土的試驗結果繪制在e-ln p平面中，得出它們的臨界狀態(tài)線，如圖7所示，圖中ec為臨界狀態(tài)孔隙比.

圖7 粗粒土的臨界狀態(tài)線Fig.7 Critical state line of coarse aggregates

根據(jù)文獻［8］提出的參數(shù)確定方法得出其他模型參數(shù)，模型參數(shù)如表2所示，表中，Γ，λ為臨界狀態(tài)線ec＝Γ-λln p中擬合參數(shù).模擬結果與試驗結果的對比如圖8所示.

表2 模型參數(shù)Tab.2 Model parameters

圖8 粗粒土試驗數(shù)據(jù)與模擬結果Fig.8 Comparison between data and simulation for coarse granular aggregates sample

由圖8可知，本文模型不僅可以比較準確地預測剪應力q和軸向應變ε1的關系，還可以反映剪縮和剪脹特性對平均應力的依存性，即隨著圍壓σ3的增大，剪脹性減小，剪縮性變大.

另外可以發(fā)現(xiàn)，圓形粗粒土的剪脹性較小，而菱形粗粒土的剪脹性較為明顯，表明顆粒形狀對于粗粒土的剪脹性具有一定的影響.值得一提的是，本文模型能夠在較大應變范圍內較好地描述粗粒土隨著應力水平的增長由體縮向體脹發(fā)展的體變特性，證明了本文所推導出的考慮顆粒破碎的剪脹方程能夠較好地描述粗粒土的剪脹特性.

3 結語

基于細觀理論本文在Ueng和Chen提出的考慮顆粒破碎粗粒土剪脹方程的基礎上對不同圍壓下顆粒破碎耗能因子在試驗過程中的變化趨勢進行了歸一化，建立了僅與顆粒破碎所消耗的有效塑性功直接相關的剪脹方程.為了驗證該剪脹方程的正確性，引用筆者曾提出的修正硬化準則，建立了考慮顆粒破碎的粗粒土彈塑性本構模型，對2種不同顆粒形狀粗粒土室內大三軸試驗的數(shù)值模擬表明，該模型能夠在較大的應變范圍內很好地描述粗粒土的應力-應變及體變特性，從而間接地證明了本文所推導出的考慮顆粒破碎剪脹方程的正確性.

［1］Rowe P W. The stress-dilatancy relation for the static equilibrium of an assembly of particles in contact ［J］／／Proceeding Royal Society，1962，269（1339）：500.

［2］Ueng T S，Chen T J.Energy aspects of particle breakage in drained shear of sands［J］.Geotechnique，2000，50（1）：65.

［3］Indrartna B，Salim W.Modeling of particle breakage of coarse aggregates incorporating strength and dilatancy ［J］.Geotechnical Engineering，2002，4：243.

［4］遲世春，賈宇峰.土顆粒破損耗能對羅維剪脹模型的修正［J］.巖土工程學報，2005，27（11）：1266.CHI Shichun，JIA Yufeng.Rowe’s stress-dilatancy model modified for energy dissipation of particle breakage［J］.Chinese Journal of Geotechnical Engineering，2005，27（11）：1266.

［5］Chavez C，Alonso E E.A constitutive model for crushed granular aggregates which includes suction effects［J］.Soils and Foundations，2003，43（4）：215.

［6］Varadarajan A，Sharma K G，Venkatachalam K.Testing and modeling two rockfill materials［J］.Journal of Geotechnical and Geoenvironmental Engineering，ASCE，2003，129（3）：206.

［7］Pietruszczak S，Stolle D.Modeling of sand behavior under earthquake excitation［J］.International Journal for Numerical and Analytical Method in geomechanics，1987，11：221.

［8］孫海忠，黃茂松.考慮顆粒破碎的粗粒土臨界狀態(tài)彈塑性本構模型［J］.巖土工程學報，2010，32（8）：1284.SUN Haizhong，HUANG Maosong.A critical state elastoplastic model for coarse granular material incorporating particle crushing［J］.Chinese Journal of Geotechnical Engineering，2010，32（8）：1284.

［9］Li X S，Dafalias Y F.Dilatancy for cohesionless soils［J］.Geotechnique，2000，50（4）：449.