鄧新征

(湖南省高速百通建設(shè)投資有限公司，湖南 長沙 410000)

0 前言

ANSYS 中內(nèi)置了基于D－P 準則的彈塑性計算，可以通過降低強度參數(shù)的方法來確定塑性區(qū)開展的范圍，并通過塑性區(qū)是否貫通來確定其安全系數(shù)，但D－P 準則屈服面的選取方式有很多，有必要對不同屈服面下確定的安全系數(shù)合理性進行研究。

關(guān)于穩(wěn)定性的計算方式，許多學(xué)者進行了研究，劉志勇等［1］基于極限平衡方法，通過Geo-studio 軟件對不同比例脫硫石膏摻入量的灰渣子壩穩(wěn)定性的數(shù)值分析。劉懷忠等［2］提出了利用剪應(yīng)變增量的方式搜索滑面，并用矢量和的方式計算邊坡的安全系數(shù)。在強度折減方面，王培濤等［3］將強度折減法引入到邊坡穩(wěn)定性研究中，并基于PFC2D 對巖質(zhì)邊坡穩(wěn)定性分析。陳星［4］基于FLCA3D 運行平臺，通過fish 語言二次開發(fā)了強度折減法，并對某邊坡進行了三維穩(wěn)定性分析。施建勇［5］等認為在對材料強度進行折減的同時，也應(yīng)該對彈性模量E 進行折減，從而提出了修正的強度折減法。遲世春等［6］通過研究發(fā)現(xiàn)對于強度折減有限元法分析土坡穩(wěn)定性，彈塑性模型更適合。

上述研究內(nèi)容對基于D－P 準則的彈塑性分析較少，對D－P 準則屈服面選擇的研究更是少見報道?；诖?，本文基于ANSYS 的二次開發(fā)語言APDL，分別采用D－P 準則的內(nèi)接圓和外接圓作為屈服面，對某土坡的穩(wěn)定性進行了研究，并將研究成果與經(jīng)典的瑞典條分法計算結(jié)果進行比較，分析了土坡穩(wěn)定計算中屈服面的選擇，以求為類似的工程研究提供參考。

1 基本理論

1.1 D－P強度準則基本理論

D－P 準則是 Drucker 和 Prager 于 1952年提出的，其基本思想是用內(nèi)切于M－C 準則的六棱錐的圓錐作為屈服面，從而避免了M－C 準則屈服面的棱角，使得屈服面變得光滑，基本函數(shù)形式為［7］:

式中，I1，J2表示應(yīng)力第1 不變量和應(yīng)力偏張量第2不變量，如式(2);α，k 為D－P 準則與材料有關(guān)的常數(shù)。對于平面應(yīng)變，α，k 與 c，φ 之間的關(guān)系可以表示為式(3)和式(4):

式中，c 為粘結(jié)力;φ 為內(nèi)摩擦角。當以外接圓為屈服面時取“－”，當以內(nèi)接圓為屈服面時取“+”。

對于實際工程，常常采用中間主應(yīng)力系數(shù)n 來反映中主應(yīng)力 σ2與大、小主應(yīng)力 σ1、σ3的關(guān)系。如式(5):

根據(jù)式(5)，σ2=nσ1+(1－n)σ3，將 σ2代入I1，J2中，將其轉(zhuǎn)化為 σ1+ σ3，σ1－σ3及 n 的關(guān)系，可得:

將式(5)和式(6)代入式(1)，可得D－P 準則的數(shù)學(xué)表達式為:

1.2 瑞典條分法

瑞典條分法由W.Fellenious 等人于1927年提出，認為滑動面為圓弧面，將邊坡劃分為若干土條，不計土條間的相互作用力，計算每個土體的滑動力矩與抗滑力矩，并將抗滑力矩與滑動力矩的比值定義為安全系數(shù)，對于整個滑面的安全系數(shù)則按照式(1)進行計算，其基本計算簡圖如圖1。

圖1 瑞典條分法計算簡圖

整個土坡相應(yīng)于滑動面(AD)的穩(wěn)定性系數(shù)為:

式中:αi為土條i 圓弧半徑與豎直線的夾角;li為土條i 滑動面ef 的弧長;ci、φi為滑動面上土的粘聚力及內(nèi)摩擦角。

在實際情況，對于比較均質(zhì)的土質(zhì)邊坡，可以近似地認為其滑裂面為圓弧面，因此可以較好地解決這類問題。

2 工程分析

某土坡坡高為10 m，坡腳為45°(如圖2所示)，室內(nèi)常規(guī)試驗表明，該土坡的力學(xué)參數(shù)如表1。結(jié)合工程斷面較長的特征，在有限元計算中采用平面應(yīng)變進行計算。

圖2 土坡模型(單位:m)

表1 邊坡力學(xué)參數(shù)表

2.1 外接圓為屈服面時的安全系數(shù)

ANSYS 的內(nèi)置D－P 準則計算方式為外接圓為屈服面，可直接進行強度折減計算。

結(jié)合下式:

本文分別對c，tanφ 進行強度折減并計算，折減系數(shù)分別為1.1，1.2 ，1.3(對應(yīng)的強度參數(shù)如表2)。并以塑性區(qū)是貫通或者計算不收斂對應(yīng)的折減系數(shù)作為安全系數(shù)。圖3列出了不同折減系數(shù)下以外接圓作為屈服面時塑性開展圖。

表2 不同折減系數(shù)下對應(yīng)的強度參數(shù)

從圖3的比變化過程中可以看出，該土坡的塑性區(qū)從無到有發(fā)展，當強度降低1.3 倍時，在坡腳處出現(xiàn)了小范圍的塑性開展區(qū)，這與實際工程中塑性區(qū)應(yīng)出現(xiàn)的是相吻合的:邊坡在自重荷載作用下，在坡腳處應(yīng)力最大，塑性區(qū)應(yīng)從此處開展。故認為此時強度折減倍數(shù)就是安全系數(shù)，即該邊坡的安全以外接圓為屈服面時對應(yīng)到安全系數(shù)為1.3。

圖3 外接圓作為屈服面時塑性區(qū)

2.2 內(nèi)接圓為屈服面時的安全系數(shù)

ANSYS 中默認的D－P 準則以外接圓作為屈服面，欲將內(nèi)接圓作為屈服面，可通過兩種途徑實現(xiàn):①通過ANSYS 內(nèi)置語言APDL 進行二次開發(fā);②基于式(3)和式(4)對強度參數(shù)進行換算，本節(jié)思路為在給定的強度參數(shù)下，按外接圓計算出相應(yīng)α 和k，式(3)和式(4)按內(nèi)接圓計算出相應(yīng)的強度參數(shù)。本文采用第2 種方法進行計算。

由于強度折減應(yīng)從材料本身的力學(xué)參數(shù)向進行折減，故采用外接圓作為屈服面時強度應(yīng)按以下方式進行折減:參照上節(jié)折減方式，先對其強度參數(shù)進行折減，折減系數(shù)分別為 1.1，1.2，1.25 ，1.3，1.35，1.4，再按照上述換算方式計算出內(nèi)接圓為屈服面時的強度參數(shù)。同樣以塑性區(qū)是貫通或者計算不收斂對應(yīng)的折減系數(shù)作為安全系數(shù)。其具不同折減系數(shù)下的塑性開展區(qū)如圖4。

計算結(jié)果表明，采用內(nèi)接圓作為屈服面的強度折減法的計算結(jié)果與外接圓計算成果相似，塑性區(qū)也是從無到有發(fā)展，當強度折減系數(shù)為1.2 時，在坡腳出現(xiàn)塑性貫通區(qū)，認為此邊坡的安全系數(shù)為1.2。

2.3 不同屈服面下安全系數(shù)分析

采用瑞典條分法按式(9)對該邊坡的穩(wěn)定性進行計算，通過自編的程序進行迭代計算，通過計算表明該邊坡的安全系數(shù)為1.29。這與外接圓屈服面計算的安全系數(shù)1.3 比較接近。因此，采用該ANSYS 內(nèi)置外接圓屈服面的計算方式是合理的。

圖4 內(nèi)接圓作為屈服面時塑性區(qū)

3 結(jié)論

分別將D－P 準則外接圓和內(nèi)接圓作為屈服面，采用強度折減法對某邊坡的穩(wěn)定性進行了研究，并將計算結(jié)果與瑞典條分法的計算成果相對比。主要得到以下結(jié)論:

1)基于D－P 準則的強度折減法在邊坡穩(wěn)定性計算中是合理的，計算結(jié)果能夠代表工程的實際情況。

2)采用外接圓作為屈服面計算的安全系數(shù)較內(nèi)接圓屈服面大。且采用外接圓計算的安全系數(shù)與瑞典條分法計算的成果相近，說明在工程計算中采用ANSYS 內(nèi)置外接圓為屈服面的彈塑性計算方式是合理的。

［1］劉志勇，陳 群，董 渝.脫硫石膏的摻入對灰渣子壩穩(wěn)定性的影響［J］.水電站設(shè)計，2015，31(1):16 －19.

［2］劉懷忠，謝紅強，何江達，等.基于FLAC3D 的邊坡滑面搜索法及動力穩(wěn)定性分析［J］.四川大學(xué)學(xué)報(工程科學(xué)版)，2014，16(S1):111 －115.

［3］王培濤，楊天鴻，于慶磊，等.節(jié)理邊坡巖體參數(shù)獲取與 PFC2D應(yīng)用研究［J］.采礦與安全工程學(xué)報，2013，30(4):560 －565.

［4］陳 星.某開挖邊坡的三維強度折減和點安全度評價［J］.人民黃河，2010，32(5):138 －139.

［5］施建勇，曹秋榮，周璐翡.修正有限元強度折減法與失穩(wěn)判據(jù)在邊坡穩(wěn)定分析中的應(yīng)用［J］.巖土力學(xué)，2013，32(S2):237 －241.

［6］遲世春，關(guān)立軍.應(yīng)用強度折減有限元法分析土坡穩(wěn)定的適應(yīng)性［J］.哈爾濱工業(yè)大學(xué)學(xué)報，2005，37(9):1298 －1302.

［7］張學(xué)言，閆澍旺.巖土塑性力學(xué)基礎(chǔ)(第2 版)［M］.天津:天津大學(xué)出版社，2004.