蔡曉鴻，蔡勇斌，蔡勇平

（1.吉安市水利局，江西 吉安 343000；2.吉安市水利水電規(guī)劃設(shè)計(jì)院，江西 吉安 343000）

0 引 言

莫爾－庫侖（Mohr－Coulomb）屈服準(zhǔn)則由于其解析表達(dá)式簡(jiǎn)明，且有很好的精確度；加之物理力學(xué)概念直觀明確，參數(shù)易于通過簡(jiǎn)單的試驗(yàn)確定；且未考慮中間主應(yīng)力σ2的影響，從而應(yīng)用于實(shí)際工程偏安全。特別是在巖土力學(xué)與工程仍處于“半理論半經(jīng)驗(yàn)”設(shè)計(jì)水準(zhǔn)的當(dāng)下，受荷載分析統(tǒng)計(jì)、內(nèi)力計(jì)算組合、內(nèi)力控制截面選取的誤差甚或巖土體計(jì)算模型與計(jì)算參數(shù)取用給不準(zhǔn)影響，在真實(shí)的建筑結(jié)構(gòu)計(jì)算分析中不存在精確解，只存在控制解或優(yōu)化解的背景下，莫爾－庫侖屈服準(zhǔn)則為工程界所樂于廣泛采用。

然而，在期刊審稿與項(xiàng)目工程設(shè)計(jì)審查中，常見到不少學(xué)術(shù)、技術(shù)人員在采用莫爾－庫侖屈服準(zhǔn)則的主應(yīng)力表達(dá)式、應(yīng)力不變量表達(dá)式及塑性流動(dòng)剪脹角參數(shù)表達(dá)式，特別是將其應(yīng)用于水工壓力隧洞或隧道工程彈塑性應(yīng)力計(jì)算、流固耦合計(jì)算分析時(shí)，未注意須與所取用應(yīng)力符號(hào)約定相關(guān)聯(lián)適配，常出現(xiàn)誤套誤用現(xiàn)象，致使應(yīng)力計(jì)算成果與后續(xù)研究產(chǎn)生不應(yīng)有的錯(cuò)誤。為避免這一錯(cuò)誤持續(xù)發(fā)生，我們?cè)团c水工壓力隧洞彈塑性應(yīng)力計(jì)算工況相適配的莫爾－庫侖屈服準(zhǔn)則主應(yīng)力表達(dá)式的合理取用做過深入分析討論[1]，并引起水工界的重視。下面，我們將進(jìn)一步探究莫爾－庫侖屈服準(zhǔn)則的應(yīng)力不變量表達(dá)式、塑性流動(dòng)剪脹角參數(shù)表達(dá)式與應(yīng)力符號(hào)約定、大小主應(yīng)力順序間的關(guān)聯(lián)適配性，以使專業(yè)技術(shù)人員對(duì)此有更明晰的認(rèn)知，避免不當(dāng)誤用。

1 與拉應(yīng)力為正、壓應(yīng)力為負(fù)約定相適配的莫爾－庫侖屈服準(zhǔn)則應(yīng)力不變量表達(dá)式

在進(jìn)行水工壓力隧洞等建筑物彈塑性應(yīng)力計(jì)算時(shí)，莫爾－庫侖屈服準(zhǔn)則的主應(yīng)力表達(dá)式，是與其應(yīng)力符號(hào)約定及大小主應(yīng)力順序相關(guān)聯(lián)協(xié)調(diào)的，參考文獻(xiàn)[1]曾詳盡地分析討論過這一問題，本文不予贅述。由于屈服準(zhǔn)則與坐標(biāo)軸方向的旋轉(zhuǎn)無關(guān)，所以屈服準(zhǔn)則除采用主應(yīng)力表示外，常用的另一表示方法是通過應(yīng)力不變量來描述。下面我們來推求與應(yīng)力符號(hào)約定及大小主應(yīng)力順序相適配的莫爾－庫侖屈服準(zhǔn)則應(yīng)力不變量表達(dá)式。

應(yīng)力符號(hào)采用連續(xù)介質(zhì)彈性力學(xué)的約定，即以拉應(yīng)力為正、壓應(yīng)力為負(fù)，相應(yīng)莫爾－庫侖屈服準(zhǔn)則的數(shù)學(xué)表達(dá)式為[2]：

式中：σ為巖土屈服面上的剪應(yīng)力；C為巖土的粘聚力；φ為巖土的摩擦角。

圖1 隧洞圍巖應(yīng)力計(jì)算簡(jiǎn)圖

當(dāng)均勻內(nèi)水壓力p0大于原巖壓力q時(shí)（圖1），圍巖徑向應(yīng)力σr為小主應(yīng)力，圍巖切向應(yīng)力σθ為大主應(yīng)力，其莫爾－庫侖屈服準(zhǔn)則表達(dá)式為（見圖 2（a）、（b））[1]：

當(dāng)均勻內(nèi)水壓力p0小于原巖壓力q時(shí)，則圍巖徑向應(yīng)力σr為大主應(yīng)力，圍巖切向應(yīng)力σθ為小主應(yīng)力，相應(yīng)莫爾－庫侖屈服準(zhǔn)則表達(dá)式為（見圖2（c））：

若以σ1和σ3分別表示最大、最小主應(yīng)力，即有σ1≥σ2≥σ3，則式（2）、（3）可統(tǒng)一于一個(gè)表達(dá)式：

下面我們來推求巖土材料莫爾－庫侖屈服準(zhǔn)則的應(yīng)力不變量表達(dá)式。將式（5）改寫成：

式中：I1為應(yīng)力張量的第一不變量；s1、s2、s3為主應(yīng)力偏量；J2為偏應(yīng)力張量的第二不變量；θ為相似角（Lode角）。

圖2 莫爾-庫倫屈服準(zhǔn)則（拉應(yīng)力為正，壓應(yīng)力為負(fù)）

圖3 莫爾-庫倫屈服準(zhǔn)則（壓應(yīng)力為正，拉應(yīng)力為負(fù)）

將式（8）代入式（7），則得與拉應(yīng)力為正、壓應(yīng)力為負(fù)應(yīng)力符號(hào)約定相適配的巖土材料莫爾－庫侖屈服準(zhǔn)則應(yīng)力不變量表達(dá)式：

2 與壓應(yīng)力為正、拉應(yīng)力為負(fù)相適配的莫爾－庫侖屈服準(zhǔn)則應(yīng)力不變量表達(dá)式

應(yīng)力符號(hào)采用巖土力學(xué)的約定，通常以壓應(yīng)力為正、拉應(yīng)力為負(fù)，莫爾－庫侖屈服準(zhǔn)則的數(shù)學(xué)表達(dá)式為：

此時(shí)，若均勻內(nèi)水壓力p0大于原巖壓力q，圍巖徑向應(yīng)力σr為大主應(yīng)力，圍巖切向應(yīng)力σθ為小主應(yīng)力，其莫爾－庫侖屈服準(zhǔn)則表達(dá)式為（圖 3（a）、（b））：

若均勻內(nèi)水壓力p0小于原巖壓力q，圍巖徑向應(yīng)力σr為小主應(yīng)力，圍巖切向應(yīng)力σθ為大主應(yīng)力，其莫爾－庫侖屈服準(zhǔn)則表達(dá)式為（圖3（c））：

以大、小主應(yīng)力σ1≥σ2≥σ3描述巖土材料莫爾－庫侖屈服準(zhǔn)則，則式（11）、（12）可統(tǒng)一于一個(gè)表達(dá)式：

將式（14）改寫為：

將式（8）代入式（15），整理得與壓應(yīng)力為正、拉應(yīng)力為負(fù)應(yīng)力符號(hào)約定相適配的巖土材料莫爾－庫侖屈服準(zhǔn)則應(yīng)力不變量表達(dá)式：

比較式（9）、式（16）可知，莫爾－庫侖屈服準(zhǔn)則的應(yīng)力不變量表達(dá)式，是與應(yīng)力符號(hào)約定、主應(yīng)力大小順序相關(guān)聯(lián)的。而以往在工程學(xué)術(shù)技術(shù)界，往往疏忽了這一相關(guān)聯(lián)性，存在采用彈性力學(xué)應(yīng)力符號(hào)約定，而誤用與巖土力學(xué)應(yīng)力符號(hào)約定相關(guān)聯(lián)的莫爾－庫侖屈服準(zhǔn)則應(yīng)力不變量表達(dá)式問題，致使后續(xù)理論探求與應(yīng)力計(jì)算分析產(chǎn)生不應(yīng)有的錯(cuò)誤。

3 剪脹角參數(shù)表達(dá)式與應(yīng)力符號(hào)約定的適配性

當(dāng)應(yīng)力符號(hào)采用塑性力學(xué)以拉應(yīng)力為正、壓應(yīng)力為負(fù)約定時(shí)，利用式（5），巖土材料的莫爾-庫侖屈服準(zhǔn)則表達(dá)式可寫成

與之相適應(yīng)的流動(dòng)法則可以用下式表達(dá)[3]：

引入塑性流動(dòng)剪脹角參數(shù)αψ，令，

類似，當(dāng)應(yīng)力符號(hào)采用巖土力學(xué)以壓應(yīng)力為正、拉應(yīng)力為負(fù)約定時(shí)，利用式（14），巖土材料的莫爾-庫侖屈服準(zhǔn)則表達(dá)式可寫成：

相適應(yīng)的流動(dòng)法則可以用下式表達(dá)：

于是據(jù)式（22），可得：

引入塑性流動(dòng)剪脹角參數(shù)βψ，令，

比較式（21）、式（26）可知，與莫爾-庫侖屈服準(zhǔn)則相關(guān)聯(lián)流動(dòng)法則的剪脹角參數(shù)αψ、βψ表達(dá)式，是與應(yīng)力符號(hào)約定及主應(yīng)力大小順序相聯(lián)系的，而以往人們對(duì)此認(rèn)知并不明晰，被許多工程師所疏忽，以致誤套誤用。有必要指出，考慮到遵守非關(guān)聯(lián)流動(dòng)法則的巖土材料模型，在應(yīng)力增量中，塑性應(yīng)變?cè)隽恳彩蔷€性的，而且滿足連續(xù)性條件，從而式（21）、式（26）對(duì)于非關(guān)聯(lián)流動(dòng)法則的巖土材料也是適用的。

4 結(jié)語

通過推求不同應(yīng)力符號(hào)約定與大小順序下的莫爾-庫侖屈服準(zhǔn)則應(yīng)力不變量表達(dá)式及塑性流動(dòng)剪脹角參數(shù)表達(dá)式，揭示了莫爾-庫侖屈服準(zhǔn)則主應(yīng)力表達(dá)式、應(yīng)力不變量表達(dá)式及塑性流動(dòng)剪脹角參數(shù)表達(dá)式與應(yīng)力符號(hào)約定、主應(yīng)力大小順序間的關(guān)聯(lián)適配性，為正確選用莫爾-庫侖屈服準(zhǔn)則表達(dá)式及塑性流動(dòng)剪脹角參數(shù)表達(dá)式提供了指南。顯然，關(guān)于這一問題討論的普適性，可推廣到其他屈服準(zhǔn)則。

有必要指出，將式（1）、（4）、（9）、（20）中 φ 用-φ 代入，則可分別導(dǎo)得式（10）、（13）、（16）、（25），即莫爾-庫侖屈服準(zhǔn)則主應(yīng)力表達(dá)式、應(yīng)力不變量表達(dá)式及塑性流動(dòng)剪脹角參數(shù)表達(dá)式與應(yīng)力符號(hào)約定間的關(guān)聯(lián)式，可轉(zhuǎn)換為與φ的正負(fù)號(hào)選擇相適配，即應(yīng)力符號(hào)以拉應(yīng)力為正，壓應(yīng)力為負(fù)，則相應(yīng)φ取負(fù)號(hào)（即莫爾-庫侖屈服直線與σ負(fù)軸間的夾角）；應(yīng)力符號(hào)以壓應(yīng)力為正，拉應(yīng)力為負(fù)，則相應(yīng)φ取正號(hào)（即莫爾-庫侖屈服直線與σ正軸間的夾角），而這便從理論上揭示了莫爾-庫侖屈服準(zhǔn)則應(yīng)力符號(hào)不同約定間的內(nèi)在聯(lián)系，為采用彈性力學(xué)數(shù)值計(jì)算方法及其應(yīng)力符號(hào)約定求解莫爾-庫侖屈服準(zhǔn)則巖土工程問題提供了一條捷徑，打開了一扇互通之門。