張 濤，何利軍

（南昌航空大學(xué)土木與建筑學(xué)院，江西 南昌330063）

含分?jǐn)?shù)階導(dǎo)數(shù)元件非線性蠕變模型的二次開發(fā)

張 濤，何利軍

（南昌航空大學(xué)土木與建筑學(xué)院，江西 南昌330063）

通過GDS應(yīng)力三軸儀對湛江黏土在恒定圍壓不同偏應(yīng)力下的三軸蠕變試驗，得到幾組不同偏應(yīng)力下軸向應(yīng)變與時間關(guān)系數(shù)據(jù)，鑒于傳統(tǒng)線性元件組合模型無法精確描述黏土蠕變的非線性問題，提出采用分?jǐn)?shù)階導(dǎo)數(shù)元件分別替換傳統(tǒng)Burger模型中兩個牛頓粘壺線性元件，建立起含兩個分?jǐn)?shù)階導(dǎo)數(shù)元件的非線性蠕變模型，采用C++語言將該模型在FLAC3D中進行二次開發(fā)，通過建立一個三軸數(shù)值算例，將不同偏應(yīng)力下數(shù)值計算結(jié)果與試驗數(shù)據(jù)進行對比，結(jié)果表明提出的非線性蠕變模型相比Burger模型更適合描述湛江黏土的蠕變過程。

分?jǐn)?shù)階導(dǎo)數(shù)；二次開發(fā)；蠕變

巖土的蠕變問題一直是工程基礎(chǔ)處理、邊坡防護和圍巖開挖與支護過程中需要重視的棘手難題，而軟黏土在工程實踐中屬于地質(zhì)不良的土體，其在我國沿海和內(nèi)陸江河沖積流域分布較為廣泛，給工程建設(shè)帶來諸多潛在危害。

軟土的蠕變是巖土流變學(xué)研究的重點內(nèi)容之一，許多研究人員依據(jù)大量的試驗和工程實踐得到了一些可觀的成果：如楊愛武[1]對天津濱海吹填軟土進行三軸不固結(jié)不排水和單軸壓縮蠕變試驗，結(jié)果表明兩類試驗歷時曲線都具有非線性特征；雷華陽[2]對濱海軟土進行不同加荷方式下的蠕變試驗，得到應(yīng)變與應(yīng)力和時間關(guān)系，結(jié)果表明濱海軟土具有明顯的非線性特征且其蠕變變形受其結(jié)構(gòu)制約和影響。這些研究成果使得人們對土體蠕變的認(rèn)識不斷加深，在本構(gòu)模型建立方面也逐漸由單純的線性理論模型試圖描述土體蠕變特性，逐漸發(fā)展為半經(jīng)驗半理論模型和非線性理論模型。

傳統(tǒng)蠕變模型都是通過粘、彈、塑三類線性元件以串并聯(lián)方式組合得到，用以描述對應(yīng)土體的粘彈、粘塑、彈塑等特性，但這些模型具有一定的局限性，只適用于描述線性蠕變過程，而實際軟土蠕變具有很強的非線性特性。本文結(jié)合前人建立非線性蠕變本構(gòu)模型研究方法和在FLAC3D中實現(xiàn)二次開發(fā)經(jīng)驗，提出采用分?jǐn)?shù)階導(dǎo)數(shù)元件（阿貝爾元件）替換傳統(tǒng)Burger模型中的牛頓粘壺，建立起變參數(shù)流變本構(gòu)模型，并在FLAC3D中完成該模型的二次開發(fā)工作，通過建立模擬三軸蠕變數(shù)值算例，比較不同偏應(yīng)力下數(shù)值結(jié)果與試驗數(shù)據(jù)的吻合程度。

湛江黏土是一種具有一定典型性的強結(jié)構(gòu)性黏土[3]，前期通過GDS應(yīng)力三軸儀對湛江黏土在恒定圍壓不同偏應(yīng)力下的三軸蠕變試驗，得到幾組不同偏應(yīng)力下軸向應(yīng)變與時間關(guān)系數(shù)據(jù)[4]，采用各種模型進行了描述[5]，由于傳統(tǒng)模型都是線性模型，基于分?jǐn)?shù)階軟體元件構(gòu)建非線性模型[6]，在理論上具有更好的描述效果。本文在FLAC3D中基于C++進行二次開發(fā)，相對于基于FISH語言二次開發(fā)而言[7]，其二次開發(fā)效果具有更好的穩(wěn)定性及方便工程技術(shù)人員采用的特點。

1 分?jǐn)?shù)階導(dǎo)數(shù)

1.1 分?jǐn)?shù)階導(dǎo)數(shù)簡介

分?jǐn)?shù)階導(dǎo)數(shù)理論相比于整數(shù)階導(dǎo)數(shù)理論起步較晚。且應(yīng)用范圍并不廣泛，在人們長期的科學(xué)實踐應(yīng)用過程中，研究人員發(fā)現(xiàn)有些問題采用整數(shù)階并不能很好的描述其發(fā)展過程，而分?jǐn)?shù)階能夠給出合理的解釋且參數(shù)簡潔明確。直到上世紀(jì)80年代，分?jǐn)?shù)階導(dǎo)數(shù)才得到廣泛應(yīng)用，主要集中在天氣和氣候研究、醫(yī)學(xué)圖像處理、工程中非線性問題處理等，很好地解釋了自然科學(xué)以及工程領(lǐng)域一些非經(jīng)典現(xiàn)象。由于分?jǐn)?shù)階導(dǎo)數(shù)計算的復(fù)雜性及無法向整數(shù)階那樣得到精確的解析解，現(xiàn)今對分?jǐn)?shù)階導(dǎo)數(shù)還沒有明確的定義，但應(yīng)用廣泛的主要有三種，分別是Riemann-Liouville定義、Grunwald-Letnikov定義和Caputo定義[8-9]，三類定義之間既有聯(lián)系又有區(qū)別。

分?jǐn)?shù)階相比于傳統(tǒng)整數(shù)階理論，優(yōu)點主要在于：① 能夠很好的體現(xiàn)系統(tǒng)函數(shù)發(fā)展過程；② 與非線性模型相比，物理意義更明確，參數(shù)更簡潔；③ 克服了經(jīng)典整數(shù)階理論與試驗數(shù)據(jù)不吻合的嚴(yán)重缺點，使用較少參數(shù)即可獲得很好的效果[10]。

1.2 分?jǐn)?shù)階導(dǎo)數(shù)元件

分?jǐn)?shù)階導(dǎo)數(shù)元件最主要的特點就是非線性[11-12]，相比傳統(tǒng)線性牛頓粘壺元件，更適合描述土體、橡膠、瀝青等具有粘性的復(fù)雜材料體，該元件最早由挪威數(shù)學(xué)家阿貝爾提出，又名Abel（阿貝爾）元件，如圖1。

圖1 阿貝爾元件Fig.1 Abel element

對上述元件，其模型本構(gòu)關(guān)系為

1.3 分?jǐn)?shù)階導(dǎo)數(shù)模型

本文建立的兩元件分?jǐn)?shù)階導(dǎo)數(shù)模型，是將Abel元件分別代替經(jīng)典Burger模型中Maxwell粘壺和Kelvin粘壺而得到（見圖2，圖3）。

圖2 Burger模型Fig.2 Burgermodel

圖3 Sburger模型Fig.3 Sburgermodel

對于Burger模型其由Maxwell模型和Kelvin模型串聯(lián)而成，分別包含一個虎克彈簧體和一個牛頓粘壺，可以用來描述土體的粘彈性，其一維本構(gòu)方程如下

上式及后式中均統(tǒng)一約定，E表示彈性模量，η表示粘滯系數(shù)，G表示剪切模量，下標(biāo)M和K分別代表Mawell部分參數(shù)和Kelvin部分參數(shù)

相應(yīng)三維本構(gòu)方程如

上式中σ1=σ1-σ3，σ11=σ1+2σ3，其它參數(shù)同上。

對于兩元件分?jǐn)?shù)階導(dǎo)數(shù)模型本構(gòu)方程的推導(dǎo)，分兩部分，第一部分導(dǎo)出分?jǐn)?shù)階Maxwell部分本構(gòu)關(guān)系；第二部分導(dǎo)出分?jǐn)?shù)階Kelvin部分本構(gòu)關(guān)系；上述兩部分串聯(lián)便得到了兩元件分?jǐn)?shù)階導(dǎo)數(shù)Burger模型，所以這兩部分本構(gòu)關(guān)系的和即為兩元件分?jǐn)?shù)階導(dǎo)數(shù)Burger模型本構(gòu)關(guān)系，具體推導(dǎo)如下：

分?jǐn)?shù)階Maxwell部分由一個虎克彈簧體和一個阿貝爾(Abel)元件串聯(lián)而成，基于前述Abel元件本構(gòu)關(guān)系，很容易給出該部分一維本構(gòu)關(guān)系如下

式中：ηa為阿貝爾粘壺基礎(chǔ)粘滯系數(shù)，r為分?jǐn)?shù)階階數(shù)。

分?jǐn)?shù)階Kelvin部分由一個虎克彈簧體和一個Abel元件串聯(lián)而成，其本構(gòu)關(guān)系推導(dǎo)如下[13]

式中：Q＝EK／ηK，D為微分算子，r為分?jǐn)?shù)階階數(shù)，若r＝1即為整數(shù)階。

將Heaviside單位階躍函數(shù)θ（t）代入式（5）并令θ（t）＝σ（t），表示在t＞0時在模型兩端施加恒定的單位應(yīng)力，通過Laplace變換再反演便得到了式（6）的蠕變?nèi)崃勘磉_式

上式中Y1＝1／EK，由展開級數(shù)再逐項反演可得，因此蠕變?nèi)崃靠蛇M一步寫為

由于上式含有無窮級數(shù)和Gamma函數(shù)，形式復(fù)雜，不便于直接應(yīng)用，需對上式進一步簡化。在上式無窮級數(shù)內(nèi)添加n=-1項，同時在無窮級數(shù)外減去此項，并進一步令r=1-β，0≤β＜1，實際表明軟粘土蠕變過程與時間密切相關(guān)，為了保留蠕變時間的指數(shù)關(guān)系，需對無窮級數(shù)內(nèi)的Gamma函數(shù)進行簡化，令β=0，則上式變?yōu)?/p>

式（9）即為分?jǐn)?shù)階Kelvin模型的蠕變?nèi)崃?，若β?即成為線性粘彈性模型，由式（9）和式（4）相加便得到兩元件分?jǐn)?shù)階Burger模型一維本構(gòu)方程，如下式

相應(yīng)三維本構(gòu)方程形式如下

若對式（11）作個簡單變換，令ηK＝ηKtβ，ηM=ηaΓ（1＋r）t1-r，便可進一步得到簡化后的下式

觀察式（12）與Burger三維本構(gòu)關(guān)系式式（3）進行對比發(fā)現(xiàn)，經(jīng)過替換后的Sburger模型本構(gòu)關(guān)系同Burger模型三維本構(gòu)關(guān)系完全一樣，兩者的區(qū)別之一在于分?jǐn)?shù)階導(dǎo)數(shù)模型中粘滯系數(shù)ηK＝ηKtβ和ηM=ηaΓ（1＋r）t1-r是隨時間成指數(shù)形式變化，而Burger模型中這兩類粘滯系數(shù)始終為常數(shù)，基于這一特點，采用C++語言將分?jǐn)?shù)階導(dǎo)數(shù)模型在FLAC3D中實現(xiàn)二次開發(fā)成為可能。

2 分?jǐn)?shù)階導(dǎo)數(shù)模型二次開發(fā)

2.1 二次開發(fā)要點

得出兩類模型差異后，將Sburger模型在FLAC3D中進行二次開發(fā)，主要工作是通過C++語言建立循環(huán)結(jié)構(gòu)[14]，實現(xiàn)參數(shù)ηK和ηM隨蠕變時間呈指數(shù)關(guān)系變化。借鑒FLAC3D中內(nèi)置的Burger模型源程序(包含. cpp源文件和.h頭文件)編寫規(guī)則和各函數(shù)功能，以該源程序為藍(lán)本通過if..else條件選擇結(jié)構(gòu)和for循環(huán)語句實現(xiàn)兩類粘滯系數(shù)隨時間的變化關(guān)系。并將更改后的分?jǐn)?shù)階模型程序文件以項目形式添加入visual studio 2008中生成Dll(動態(tài)鏈接庫)文件[15]，并命名為userSburger.dll。然后將該鏈接復(fù)制到FLAC3D內(nèi)置模型程序文件目錄下，通過在FLAC3D窗口輸入相關(guān)命令，完成所開發(fā)模型同F(xiàn)LAC3D的接口、載入，初步實現(xiàn)該模型二次開發(fā)。

模型開發(fā)過程中對程序的更改主要集中在基本參數(shù)的定義、常量與過程量的定義、基類的描述、模型注冊ID號和Run（）函數(shù)中本構(gòu)模型計算關(guān)系的修改等。

2.2 二次開發(fā)驗證

在開發(fā)模型同F(xiàn)LAC3D進行接口、載入后，需要對模型開發(fā)的正確性及用于數(shù)值計算的可靠性進行驗證。驗證方法分為兩部分，一是判斷所開發(fā)模型程序編寫是否正確，即該模型能否用于數(shù)值計算；其次判斷所開發(fā)模型是否實現(xiàn)兩類粘滯系數(shù)隨時間成指數(shù)形式變化，這是開發(fā)成功的關(guān)鍵。首先通過建立一個簡單算例，在Sburger模型中對兩個分?jǐn)?shù)階階數(shù)γ和β取值0退化成Burger模型形式，得ηK=ηK和ηM=ηa；保持其它參數(shù)不變，計算結(jié)果如圖4。

圖4 Sburger模型退化后計算結(jié)果Fig.4 Calculation of deteriorated Sburgermodel

從圖4可以看出，對Sburger模型的分?jǐn)?shù)階階數(shù)取值為0后，保持其他參數(shù)相同，其數(shù)值計算結(jié)果與內(nèi)置Burger模型完全吻合。說明兩類分?jǐn)?shù)階模型按照上述方式通過取特殊值可退化為經(jīng)典Burger模型，也進一步表明開發(fā)的兩類分?jǐn)?shù)階模型是可以在FLAC3D中作為本構(gòu)模型來選取并參與數(shù)值計算。

圖5 Sburger模型中不同階數(shù)數(shù)值計算曲線Fig.5 Numerical calculation curves of different orders in Sburgermodel

從圖5可得：對Sburger模型階數(shù)取不同值，數(shù)值計算曲線存在明顯差異，其基本規(guī)律為當(dāng)β取小，γ取大時，試樣在等速蠕變階段的變形速率較大；而β取大，γ取小時，試樣在等速蠕變階段的變形速率較小；出現(xiàn)這一規(guī)律的原因在于β取值越小使得ηK隨時間增加的變化率下降，即單位時間內(nèi)的粘滯系數(shù)增長率降低，導(dǎo)致在恒定偏應(yīng)力作用下，其對應(yīng)蠕變變形量增加；反之β取值越大，單位時間內(nèi)粘滯系數(shù)增長率較高，使得對應(yīng)蠕變變形量降低。同理對于另一階數(shù)γ，其取相同值所得ηM的結(jié)果正好與β相反。綜上，兩個分?jǐn)?shù)階階數(shù)的不同取值，將影響兩類粘滯系數(shù)隨時間變化增長速率的快慢，從而在恒定偏應(yīng)力下，使得等速蠕變階段蠕變變形存在顯著差異。

3 三軸數(shù)值算例

3.1 三軸蠕變試驗數(shù)據(jù)

兩類模型的基本參數(shù)均依據(jù)三軸蠕變試驗數(shù)據(jù)通過各模型三維本構(gòu)關(guān)系擬合得到，表1為湛江軟粘土試樣基本物理力學(xué)性質(zhì)[16]，圖6為該試樣在50 kPa恒定圍壓不同偏應(yīng)力下，柱體軸向應(yīng)變量隨時間變化關(guān)系。

表1 湛江軟粘土試樣基本物理力學(xué)性質(zhì)Tab.1 Basic physical andmechanical properties of Zhanjiang clay samp le

圖6 三軸蠕變試驗曲線Fig.6 Curves of triaxial creep test

圖6是對試樣進行分級加載條件下得到，而只有分別加載才能得到更加真實的土體蠕變曲線，且考慮到巖土材料流變的非線性特性。需將分級加載下的蠕變數(shù)據(jù)轉(zhuǎn)化為分別加載下的蠕變數(shù)據(jù)，這里采用陳宗基先生提出并由其學(xué)生發(fā)展的“陳氏法”[17]進行處理，整理后的蠕變數(shù)據(jù)如圖7。

圖7 經(jīng)“陳氏法”整理后三軸蠕變試驗曲線Fig.7 Curves of triaxial creep test under Chen’smethod

3.2 模型建立

為了對比分?jǐn)?shù)階導(dǎo)數(shù)模型與Burger模型描述湛江軟粘土蠕變過程的可靠性，在FLAC3D中建立一個圓柱體模型(如圖8)算例模擬三軸蠕變過程。圓柱高為2個單位，長和寬均為0．5個單位，半徑為0．25個單位，先對模型施加50 kPa的圍壓和重力加速度，然后模型各方向位移清零，再在上頂面施加豎直方向壓力，蠕變時步取為1e-4，分別采用Burger模型和本文開發(fā)的分?jǐn)?shù)階導(dǎo)數(shù)Sburger模型進行計算。模型參數(shù)取值依據(jù)相應(yīng)本構(gòu)關(guān)系基于三軸蠕變試驗數(shù)據(jù)擬合得到。

圖8 柱體模型示意圖Fig.8 Schematic diagram of cylindermodel

3.3 參數(shù)選取

采用1Stopt數(shù)學(xué)軟件及孫均[18]提出的擬合方法擬合得到兩類模型的基本參數(shù)，大致可分為兩部分，第一部分當(dāng)t→o時，只有Maxwell彈簧存在變形，可擬合出參數(shù)K和GM，第二部分當(dāng)t→∞時，蠕變達到穩(wěn)定，將此時本構(gòu)關(guān)系式減去原本構(gòu)關(guān)系式便得到包含粘性和粘彈性本構(gòu)關(guān)系式，與對應(yīng)的三軸試驗數(shù)據(jù)進行回歸擬合。得到Burger模型在不同偏應(yīng)力下參數(shù)如表2，同理可得分?jǐn)?shù)階Sburger模型參數(shù)如表3。并由擬合相關(guān)系數(shù)R可以看出，各參數(shù)擬合情況較好。

表2 恒定圍壓不同偏應(yīng)力下burger模型相關(guān)參數(shù)Tab.2 Parameters of Burgermodel under different partial stress of constant confining pressure

表3 恒定圍壓不同偏應(yīng)力下Sburger模型相關(guān)參數(shù)Tab.3 Parameters of Sburgermodel under different partial stress of constant confining pressure

由表2和表3中兩個模型的相關(guān)擬合系數(shù)R可知，Sburger模型的擬合效果要優(yōu)于Burger模型。其原因在于湛江黏土的蠕變包含有非線性變形部分，而傳統(tǒng)線性模型無法準(zhǔn)確描述非線性問題，所以采用非線性本構(gòu)關(guān)系擬合參數(shù)的效果要比線性本構(gòu)關(guān)系好；其次，Burger模型在30 kPa偏應(yīng)力下，擬合系數(shù)較高，這是在低偏應(yīng)力下，黏土蠕變主要產(chǎn)生線性彈性、粘性和粘彈性變形；而當(dāng)偏應(yīng)力增加，隨著蠕變變形的進一步增長，非線性變形也在不斷上升成為后期變形的主要增長部分，是導(dǎo)致65～135 kPa偏應(yīng)力下Burger模型擬合效果較低的原因之一。另外，從參數(shù)擬合效果來看，在數(shù)學(xué)角度上能一定程度說明Sburger模型本構(gòu)關(guān)系要比經(jīng)典Burger模型更能描述試驗曲線。

將上述兩模型采用對應(yīng)參數(shù)分別在不同偏應(yīng)力下進行三軸數(shù)值計算，并將兩者數(shù)值計算結(jié)果同三軸蠕變試驗數(shù)據(jù)進行對比，如圖9～圖12。

由圖9～圖12不難看出，開發(fā)兩元件分?jǐn)?shù)階導(dǎo)數(shù)模型比傳統(tǒng)經(jīng)典Burger模型，更適合描述湛江軟粘土的初始蠕變和穩(wěn)定蠕變階段。從而驗證了本文提出采用分?jǐn)?shù)階元件替換經(jīng)典整數(shù)階模型中牛頓粘壺元件，建立變參數(shù)流變模型，描述軟土非線性蠕變思路的正確性；也進一步說明，本文對分?jǐn)?shù)階導(dǎo)數(shù)模型在

FLAC3D中實現(xiàn)二次開發(fā)的正確性和可行性。

圖9 偏應(yīng)力30 kPa下數(shù)值計算曲線Fig.9 Numerical calculation curve under the partial stress 30 kPa

圖10 偏應(yīng)力65 kPa下數(shù)值計算曲線Fig.10 Numerical calculation curve under the partial stress 65 kPa

圖11 偏應(yīng)力100 kPa下數(shù)值計算曲線Fig.11 Numerical calculation curve under the partial stress 100 kPa

圖12 偏應(yīng)力135 kPa下數(shù)值計算曲線Fig.12 Numerical calculation curve under the partial stress 135 kPa

4 結(jié)論

從本文提出分?jǐn)?shù)階元件替換傳統(tǒng)線性粘壺元件，建立變參數(shù)蠕變模型，并在FLAC3D中完成二次開發(fā)和驗證的結(jié)果來看，分?jǐn)?shù)階元件相對于傳統(tǒng)元件，對蠕變曲線的刻畫能力更強，預(yù)示著對工程的數(shù)值模擬結(jié)果更加準(zhǔn)確。由于土體蠕變變形具有廣泛的非線性特性，分?jǐn)?shù)階元件相對于整數(shù)階元件，將更具備普遍意義。伴隨著工程建設(shè)的需要和蠕變模型自身的發(fā)展，本文研究成果將在未來的巖土工程設(shè)計計算中被越來越多地接納和采用。

[1]楊愛武，張兆杰，孔令偉.不同應(yīng)力狀態(tài)下軟黏土蠕變特性試驗研究[J].巖土力學(xué)，2014，35（S2）：53-60.

[2]雷華陽，賈亞芳，李肖.濱海軟土非線性蠕變特性的試驗研究[J].巖石力學(xué)與工程學(xué)報，2013，31（1）：2806-2816.

[3]拓勇飛，孔令偉，郭愛國，等.湛江地區(qū)結(jié)構(gòu)性軟土的賦存規(guī)律及其工程特性[J].巖土力學(xué)，2004（12）：1879-1884.

[4]蔡羽，孔令偉，郭愛國，等.剪應(yīng)變率對湛江強結(jié)構(gòu)性黏土力學(xué)性狀的影響[J].巖土力學(xué)，2006（8）：1235-1240.

[5]何利軍，蔡羽，吳應(yīng)紅.湛江軟黏土的流變本構(gòu)模型研究[J].合肥工業(yè)大學(xué)學(xué)報：自然科學(xué)版，2011（4）：565-570.

[6]何利軍，孔令偉，吳文軍，張先偉，蔡羽.采用分?jǐn)?shù)階導(dǎo)數(shù)描述軟黏土蠕變的模型[J].巖土力學(xué)，2011（S2）：239-243.

[7]何利軍.分?jǐn)?shù)階導(dǎo)數(shù)蠕變模型及在FLAC3D中的數(shù)值實現(xiàn)[J].南昌航空大學(xué)學(xué)報：自然科學(xué)版，2015（4）：35-39.

[8]李秀.幾類粘彈性材料的分?jǐn)?shù)階積分模型及其數(shù)值解法[D].寧夏：寧夏大學(xué)，2013：8-24.

[9]王玉嬌.分?jǐn)?shù)階導(dǎo)數(shù)及其應(yīng)用[D].太原：太原理工大學(xué)，2014：10-44.

[10]陳文，孫洪廣，李西成.力學(xué)與工程問題的分?jǐn)?shù)階導(dǎo)數(shù)建模[M].北京：科學(xué)出版社，2010：124-140.

[11]郭佳奇，喬春生，徐沖，等.基于分?jǐn)?shù)階微積分的Kelvin-Voigt流變模型[J].中國鐵道科學(xué)，2009，30（4）：1-6.

[12]何志磊，朱珍德，朱明禮，等.基于分?jǐn)?shù)階導(dǎo)數(shù)的非定常蠕變本構(gòu)模型研究[J].巖土力學(xué)，2016，37（3）：737-744.

[13]王志方，張國忠，劉剛.膠凝原油的黏彈性流變模型[J].高校化學(xué)工程學(xué)報，2008，22（2）：351-355.

[14]古萬榮著.Visual C++完全自學(xué)手冊[M].北京：機械工業(yè)出版社，2009：29-111.

[15]陳育民，徐鼎平.FLAC/FLAC3D基礎(chǔ)與工程實例[M].2版.北京：中國水利水電出版社，2013：50-165.

[16]孔令偉，何利軍，張先偉.湛江黏土的蠕變模型與變參數(shù)塑性元件[J].巖土力學(xué)，2012，33（8）：2241-2246.

[17]陳宗基，康文法.巖石的封閉應(yīng)力、蠕變和擴容及本構(gòu)方程[J].巖石力學(xué)與工程學(xué)報，1991，10（4）：299-212.

[18]孫鈞著.巖土材料流變及其工程應(yīng)用[M].北京：中國建筑工業(yè)出版社，1999：146-158.

Secondary Development of Non-linear Creep Modelw ith Fractional Order Derivative Elements

Zhang Tao，He Lijun
（School of Civil Engineering and Architecture，Nanchang Hangkong University，Nanchang 330063，China）

Triaxial creep test of Zhanjiang clay was conducted at constant confining pressure and under different deviatoric stress through GDS,which obtained several groups of data concerning the relationship between axial strain and time under different deviatoric stress.In view of the fact that traditional linear elements of nonlinear creepmodel does not accurately describe the clay problem,a fractional order derivative elementwas proposed to replace the two linear elements of the Newton’s pot in the Burgermodel respectively.Then a non-linear creepmodel with two fractional derivatives was established.By using C++language,the secondary development of themodel in the FLAC3D was completed.Finally,comparison wasmade between numerical results under different stress and experimental data by establishing a triaxial numerical example.The results show that compared with the Burgermodel,the non-linear creepmodel ismore favorable for the description of creep process for Zhanjiang clay.

fractional derivative；secondary development；creep

1005-0523（2017）05-0021-08

T U433

A

2017－05－09

南昌航空大學(xué)校級基金項目（EA201111317）

張濤（1991一），男，碩士研究生，主要研究方向為蠕變模型的二次開發(fā)。

何利軍（1977一），男，講師，博士，研究方向為結(jié)構(gòu)性土體本構(gòu)模型研究。

（責(zé)任編輯 王建華）