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

張 濤，何利軍

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

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

張 濤，何利軍

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

圖2 Burger模型Fig.2 Burgermodel

圖3 Sburger模型Fig.3 Sburgermodel

對(duì)于Burger模型其由Maxwell模型和Kelvin模型串聯(lián)而成，分別包含一個(gè)虎克彈簧體和一個(gè)牛頓粘壺，可以用來(lái)描述土體的粘彈性，其一維本構(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ù)同上。

對(duì)于兩元件分?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部分由一個(gè)虎克彈簧體和一個(gè)阿貝爾(Abel)元件串聯(lián)而成，基于前述Abel元件本構(gòu)關(guān)系，很容易給出該部分一維本構(gòu)關(guān)系如下

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

分?jǐn)?shù)階Kelvin部分由一個(gè)虎克彈簧體和一個(gè)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時(shí)在模型兩端施加恒定的單位應(yīng)力，通過(guò)Laplace變換再反演便得到了式（6）的蠕變?nèi)崃勘磉_(dá)式

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

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

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

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

若對(duì)式（11）作個(gè)簡(jiǎn)單變換，令ηK＝ηKtβ，ηM=ηaΓ（1＋r）t1-r，便可進(jìn)一步得到簡(jiǎn)化后的下式

觀察式（12）與Burger三維本構(gòu)關(guān)系式式（3）進(jìn)行對(duì)比發(fā)現(xiàn)，經(jīng)過(guò)替換后的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í)間成指數(shù)形式變化，而B(niǎo)urger模型中這兩類(lèi)粘滯系數(shù)始終為常數(shù)，基于這一特點(diǎn)，采用C++語(yǔ)言將分?jǐn)?shù)階導(dǎo)數(shù)模型在FLAC3D中實(shí)現(xiàn)二次開(kāi)發(fā)成為可能。

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

2.1 二次開(kāi)發(fā)要點(diǎn)

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

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

2.2 二次開(kāi)發(fā)驗(yàn)證

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

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

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

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

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

3 三軸數(shù)值算例

3.1 三軸蠕變?cè)囼?yàn)數(shù)據(jù)

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

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

圖6 三軸蠕變?cè)囼?yàn)曲線Fig.6 Curves of triaxial creep test

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

圖7 經(jīng)“陳氏法”整理后三軸蠕變?cè)囼?yàn)曲線Fig.7 Curves of triaxial creep test under Chen’smethod

3.2 模型建立

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

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

3.3 參數(shù)選取

采用1Stopt數(shù)學(xué)軟件及孫均[18]提出的擬合方法擬合得到兩類(lèi)模型的基本參數(shù)，大致可分為兩部分，第一部分當(dāng)t→o時(shí)，只有Maxwell彈簧存在變形，可擬合出參數(shù)K和GM，第二部分當(dāng)t→∞時(shí)，蠕變達(dá)到穩(wěn)定，將此時(shí)本構(gòu)關(guān)系式減去原本構(gòu)關(guān)系式便得到包含粘性和粘彈性本構(gòu)關(guān)系式，與對(duì)應(yīng)的三軸試驗(yàn)數(shù)據(jù)進(jìn)行回歸擬合。得到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中兩個(gè)模型的相關(guān)擬合系數(shù)R可知，Sburger模型的擬合效果要優(yōu)于Burger模型。其原因在于湛江黏土的蠕變包含有非線性變形部分，而傳統(tǒng)線性模型無(wú)法準(zhǔn)確描述非線性問(wèn)題，所以采用非線性本構(gòu)關(guān)系擬合參數(shù)的效果要比線性本構(gòu)關(guān)系好；其次，Burger模型在30 kPa偏應(yīng)力下，擬合系數(shù)較高，這是在低偏應(yīng)力下，黏土蠕變主要產(chǎn)生線性彈性、粘性和粘彈性變形；而當(dāng)偏應(yīng)力增加，隨著蠕變變形的進(jìn)一步增長(zhǎng)，非線性變形也在不斷上升成為后期變形的主要增長(zhǎng)部分，是導(dǎo)致65～135 kPa偏應(yīng)力下Burger模型擬合效果較低的原因之一。另外，從參數(shù)擬合效果來(lái)看，在數(shù)學(xué)角度上能一定程度說(shuō)明Sburger模型本構(gòu)關(guān)系要比經(jīng)典Burger模型更能描述試驗(yàn)曲線。

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

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

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

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

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

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

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

4 結(jié)論

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

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

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

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

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

[5]何利軍，蔡羽，吳應(yīng)紅.湛江軟黏土的流變本構(gòu)模型研究[J].合肥工業(yè)大學(xué)學(xué)報(bào)：自然科學(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ù)值實(shí)現(xiàn)[J].南昌航空大學(xué)學(xué)報(bào)：自然科學(xué)版，2015（4）：35-39.

[8]李秀.幾類(lèi)粘彈性材料的分?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é)與工程問(wèn)題的分?jǐn)?shù)階導(dǎo)數(shù)建模[M].北京：科學(xué)出版社，2010：124-140.

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

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

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

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

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

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

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

[18]孫鈞著.巖土材料流變及其工程應(yīng)用[M].北京：中國(guó)建筑工業(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é)校級(jí)基金項(xiàng)目（EA201111317）

張濤（1991一），男，碩士研究生，主要研究方向?yàn)槿渥兡Ｐ偷亩伍_(kāi)發(fā)。

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

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