李賽，汪優(yōu)，秦志浩,劉建華

(1.中南大學 土木工程學院，湖南 長沙 410075;2.長沙理工大學 交通運輸工程學院，湖南 長沙 410114)

?

基于統(tǒng)計損傷本構(gòu)模型的改進接觸面模型研究

李賽1，汪優(yōu)1，秦志浩1,劉建華2

(1.中南大學 土木工程學院，湖南 長沙 410075;2.長沙理工大學 交通運輸工程學院，湖南 長沙 410114)

考慮樁土界面初始剪切剛度的深度效應，把接觸面剪切應力和剪切應變沿接觸面厚度均勻化，提出基于統(tǒng)計損傷本構(gòu)模型的無厚度單元接觸面本構(gòu)模型，根據(jù)剪切實驗的數(shù)據(jù)對改進的本構(gòu)模型參數(shù)進行擬合，得到了不同法向應力下的模擬曲線。計算結(jié)果表明:本文提出的無厚度單元接觸面本構(gòu)模型是正確的，能夠合理地模擬樁土界面的力學行為，可以描述接觸面剪切剛度的非線性變化,適用性強，有較高的理論和實際應用價值。

樁土界面；無厚度單元；統(tǒng)計損傷本構(gòu)模型；擬合參數(shù)

目前，城市建設(shè)、公路、鐵路以及橋梁等工程日益發(fā)展,實際工程面臨越來越復雜的地質(zhì)條件，樁基礎(chǔ)憑借著處理復雜地基方面的獨特優(yōu)勢而受到廣泛應用。樁土相互作用問題是巖土工程領(lǐng)域的一個研究熱點[1], 樁土界面的荷載傳遞模型對于預測樁的承載變形性狀至關(guān)重要。然而，由于樁體與土體的材料性質(zhì)存在顯著性差異，復雜應力條件下，樁土界面產(chǎn)生不一致的變形，致使樁- 土接觸面上產(chǎn)生錯動、滑移或者開裂等非連續(xù)變形，從而使樁- 土接觸面力學特性變得復雜。因此，建立準確合理的接觸面本構(gòu)模型是研究樁基承載力、沉降及變形等問題的前提。

針對樁土荷載傳遞關(guān)系，不少學者進行了深入的研究。20世紀50年代，Seed等[2]基于荷載傳遞法的雙曲線模型分析了樁土界面荷載傳遞關(guān)系。Desai等[3]首次將損傷力學的基本理論應用于接觸面本構(gòu)關(guān)系，為接觸面的研究提供了新的思路。張忠苗課題組[4]提出了可考慮樁土軟化的樁側(cè)傳遞函數(shù)的統(tǒng)一三折線模型。蒲訶夫[5]基于提出的改進雙曲線模型，成功地在FLAC3D實現(xiàn)了樁土界面模型的二次開發(fā)。彭凱等[6]采用單剪儀對接觸面在泥皮條件下進行剪切試驗，分析了曲線出現(xiàn)的軟化、剪脹和剪縮現(xiàn)象。石熊[7-8]采用大型直剪儀進行靜載及動載下紅黏土與混凝土接觸面的單向直剪試驗，分析了紅黏土與混凝土接觸面的應力、應變規(guī)律及破壞形式。馮大闊等[9]運用80T三維多功能土工試驗機，研究了粗粒土與人造粗糙鋼板形成的接觸面在循環(huán)剪切作用下的動力學特性。楊林德等[10]以接觸面內(nèi)部缺陷分布的隨機性為出發(fā)點，依據(jù)連續(xù)強度理論和統(tǒng)計理論， 建立了土與結(jié)構(gòu)物接觸面的統(tǒng)計損傷本構(gòu)模型。夏紅春[11]提出了基于平面應變條件下Mohr-Coulomb屈服準則的土-結(jié)構(gòu)接觸面微元強度表達式, 基于微元破壞服從雙參數(shù)Weibull分布的假設(shè), 建立了土-結(jié)構(gòu)接觸面統(tǒng)計損傷軟化本構(gòu)模型。

上述研究者根據(jù)接觸面剪切力學試驗研究成果，在接觸面本構(gòu)模型的研究方面取得了一定進展，這些樁土界面荷載傳遞模型能夠較好地分析樁土相互作用，滿足特定的工程應用需要。然而，這些學者不管是基于剪切實驗定性地分析數(shù)據(jù)，還是基于實驗數(shù)據(jù)提出接觸面本構(gòu)模型，都視接觸面厚度為直接的重要影響因素，推導的本構(gòu)模型計算時得考慮接觸面厚度?，F(xiàn)有數(shù)值計算軟件中的無厚度接觸面單元都把界面剪切剛度假定為定值，而實際上接觸面剪切剛度是變化的，很少有人在具有無厚度單元特性的接觸面本構(gòu)模型方面開展進一步研究，提出可以描述接觸面剛度非線性的樁土界面荷載傳遞模型。 本文嘗試基于已有的統(tǒng)計損傷界面本構(gòu)模型，提出可適用于數(shù)值模擬應用的無厚度單元接觸面本構(gòu)模型，用以描述樁土接觸面剛度的非線性變化；為了驗證推導的接觸面本構(gòu)模型的合理性，根據(jù)實驗數(shù)據(jù)對改進的本構(gòu)模型參數(shù)進行擬合，將模擬結(jié)果與實驗結(jié)果進行比較。

1　改進的統(tǒng)計損傷本構(gòu)模型

1.1接觸面統(tǒng)計損傷本構(gòu)模型的介紹

楊林德[10]基于連續(xù)損傷理論和統(tǒng)計強度理論, 建立了接觸面的統(tǒng)計損傷本構(gòu)模型,其表達式為：

(1)

式中：τ為剪應力；G為剪切模量；γ為剪應變，m；F是參數(shù)，可以根據(jù)接觸面剪切實驗數(shù)據(jù)擬合求出。

接觸面統(tǒng)計損傷本構(gòu)模型中的3個參數(shù), 即G、m和F, 可以根據(jù)實驗數(shù)據(jù)確定。首先, 按照材料力學中初始剪切模量的確定方法, 可以求解出式(1)中的參數(shù)G?？紤]到G值受正應力的影響比較大，作者根據(jù)Evgin等( 1996 年) 的研究結(jié)論：正應力的大小對歸一化剪應力-剪切位移(τ/σn-γ)曲線的影響, 除在峰值點較大外, 對曲線其它部分的影響幾乎可以忽略不計。提出根據(jù)歸一化剪應力-剪切位移曲線(τ/σn-γ), 可求得歸一化切向平均剪切模量G’。于是，接觸面在法向應力σn作用下的切向平均剪切模量G，按下式計算

G=G′σn

(2)

則式(1)可以轉(zhuǎn)化為:

Y=A+BX

(3)

由式(2)可知參數(shù)G,再根據(jù)接觸面剪切實驗得到的曲線τ-γ上點的數(shù)據(jù)就可以擬合出式(3)表示的直線，進而得到參數(shù)m、F。

1.2接觸面統(tǒng)計損傷本構(gòu)模型的改進

楊林德推出的接觸面本構(gòu)模型能夠描述接觸面的剪切力學特性，但根據(jù)文獻[10]中的試驗曲線與模擬曲線的對比情況，仍有不小的誤差，因此該模型還需要進一步的改進。

土的粒徑、顆粒級配、密度、接觸面粗糙程度以及施加于接觸面上的法向應力都是影響接觸面厚度t的因素，接觸面間剪切錯動帶的厚度t需要根據(jù)經(jīng)驗或大量的試驗確定[12-14]，考慮到t的值比較小，筆者假設(shè)剪切應力和剪切應變沿接觸面厚度均勻分布，即剪應變γ與剪切位移Δ成線性關(guān)系，剪切應變γ=kΔ，其中，k=1/t,t為接觸面剪切實驗中接觸面的厚度。這樣式(1)就可以直接轉(zhuǎn)化為：

(4)

再根據(jù)式(4)對位移Δ求導，可以得到剪切剛度ks與剪切位移的關(guān)系為：

(5)

令位移Δ=0，則初始剪切剛度為：

ksi=Gk

(6)

為了考慮樁土之間的滑移，采用Alonso[15]提出的方法來確定樁土界面初始剛度, 即等于極限剪應力和極限相對位移的比值。大量實驗表明，混凝土樁-土界面剪切試驗得到的極限相對剪切位移受法向應力σn影響不大。根據(jù)剪切實驗擬合出的極限剪應力τult-σn直線可以求出接觸面的粘聚力c和內(nèi)摩擦角φ，則極限剪應力為：

τult=c+σntanφ

(7)

由于樁土界面初始剛度等于極限剪應力和極限相對位移的比值，故可求出樁土界面初始剪切剛度ksi。

式(4)可以轉(zhuǎn)化為

Y=C+DX

(8)

由式(6)可知

因此只需求出初始剪切剛度ksi，結(jié)合接觸面實驗τ-Δ曲線上的點就可求出式(8)中的C和D。參數(shù)m和F可以用C和D表示，即：

(9)

將式(6)和式(9)代入式(4)中得：

τ=ksiΔexp[-ΔDexpC]

(10)

同理式(5)變?yōu)椋?/p>

ks=ksiexp[-ΔDexpC]-ksiDΔDexpCexp[-ΔDexpC]

(11)

由式(10)和(11)可知接觸面剪切實驗中接觸面的厚度t對剪切應力和剪切剛度無直接影響，也就是說本文開發(fā)的接觸面本構(gòu)模型計算時可忽略接觸面厚度。產(chǎn)生這種結(jié)果是因為公式推導過程中的兩個因素：其一，剪應變γ與剪切位移Δ線性關(guān)系的假設(shè)把接觸面剪切應力和剪切應變沿厚度均化；其二，通過式(4)對位移Δ求導，把初始剪切剛度ksi引入到公式(10)，但初始剪切剛度ksi的解法并沒有直接考慮接觸面厚度。與文獻[10]中的統(tǒng)計損傷界面本構(gòu)模型相比，改進的本構(gòu)模型是一個與接觸面厚度無關(guān)的單元，本文定義為“無厚度單元”，且無厚度單元接觸面本構(gòu)模型能描述接觸面剪切剛度的非線性，可以考慮應用于接觸面數(shù)值模擬。 因此，本文提出的接觸面本構(gòu)模型適用性強，有較高的理論和實際應用價值。

2　不同法向應力下的參數(shù)擬合

為了驗證接觸面本構(gòu)模型的合理性和有效性，本文選取文獻[16]中樁-土接觸面剪切性質(zhì)的室內(nèi)單剪試驗數(shù)據(jù)進行曲線擬合，含水量17%的黃色黏土與混凝土接觸面試樣單剪試驗數(shù)據(jù)如圖1所示。筆者將取法向應力σn分別等于22，55，88和107kPa時實驗曲線上的點，對本文的無厚度單元接觸面本構(gòu)模型參數(shù)進行擬合。

圖1　含水量17%接觸面試驗曲線Fig.1 Experimental curve on the interface of moisture content 17%

2.1法向應力σn=22 kPa情況下參數(shù)擬合的過程

由式(10)和(11)可知ksi對曲線的擬合影響很大，選擇合理的求解方法尤為重要。為了考慮樁土之間的滑移，可采用Alonso[15]提出的方法來確定樁土界面初始剛度, 即極限剪應力和極限相對位移的比值。如圖2所示, 取樁土界面的τ-Δ曲線的初始切線與τ=τult線的交點對應的樁土相對位移為表觀極限相對位移Δau, 也就是說Δau=τult/ksi。當取破壞比Rf= 1. 0時,Δau=τult/ksi，為了將Δau與Δu建立聯(lián)系, 引入一個參數(shù)χ=Δu/Δau, 通過調(diào)整參數(shù)χ使曲線模型在樁土相對位移等于Δu時的剪應力τ盡可能接近極限剪應力τult。這樣就建立起樁土界面的初始剪切剛度ksi與極限剪應力之間的關(guān)系,即ksi=χτult/Δu，其中極限剪應力越大,ksi越大。

圖2　樁土界面剪應力與相對位移的關(guān)系曲線Fig.2 Pile soil interface shear stress and relative displacement curve

式中：y為模擬應力值；yi為實驗應力值；n為實驗曲線上點的個數(shù)。

圖3　誤差系數(shù)變化圖Fig.3 Error coefficient curve

Table1Fittingparametervaluesunderthenormalstress22kPa

法向應力/kPaχ△u/mmDC相關(guān)系數(shù)R22212.52.50.36483.12360.9991

圖4　法向應力σn =22 kPa情況下的模擬效果圖Fig.4 Simulation results under the normal stress 22 kPa

由圖3可知χ= 12.5時誤差系數(shù)最小，為6.134 5，擬合的參數(shù)見表1，將擬合出的參數(shù)代入公式(10)得到模擬曲線見圖4。

2.2σn=55 kPa、σn=88 kPa和σn=107 kPa情況下參數(shù)擬合的結(jié)果

基于上文法向應力為22 kPa時的模擬過程，分別擬合出法向應力σn=55，σn=88和σn=107 kPa時的參數(shù)值，具體參數(shù)見表2 。根據(jù)ksi=χτult/Δu求出不同法向應力下的接觸面初始剪切剛度見表3。不同法向應力情況下的模擬曲線跟實驗曲線的對比結(jié)果見圖5，圖6和圖7。

表2　擬合獲取的參數(shù)值

表3不同法向應力下的初始剪切剛度值

Table3Initialshearstiffnessvaluesunderdifferentnormalstresses

法向應力/kPa225588107初始剪切剛度/(MPa·m-1)135.9141.03167.4234.28

由表3可知隨著法向應力的增加，接觸面初始剪切剛度剛度不斷增大，與剪切實驗曲線上的變化規(guī)律相同。根據(jù)圖4，圖5，圖6和圖7中不同法向應力下模擬曲線與實驗曲線的對比情況可以看出，模擬曲線和實驗曲線能夠比較準確地吻合，擬合參數(shù)的相關(guān)系數(shù)均在0.99以上。因此，本文提出的改進統(tǒng)計損傷本構(gòu)模型是可行的，能夠合理地模擬出樁土界面的力學行為。同時本文改進的接觸面本構(gòu)模型參數(shù)少，擬合過程簡便且具有無厚度接觸面單元的特性，在數(shù)值模擬分析樁土界面力學行為方面占有一定的優(yōu)勢，有較高的工程應用價值。

圖5　法向應力σn =55 kPa情況下的模擬效果圖Fig.5 Simulation results under the normal stress 55 kPa

圖6　法向應力σn =88 kPa情況下的模擬效果圖Fig.6 Simulation results under the normal stress 88 kPa

圖7　法向應力σn =107 kPa情況下的模擬效果圖Fig.7 Simulation results under the normal stress 107 kPa

3　結(jié)論

1)本文提出了基于統(tǒng)計損傷本構(gòu)模型的無厚度單元接觸面本構(gòu)模型,能描述接觸面剛度的非線性變化，適用性強，具有理論和實際應用價值。

2)根據(jù)剪切實驗的數(shù)據(jù)對改進的接觸面本構(gòu)模型參數(shù)進行擬合，得到了不同法向應力下的模擬曲線，將模擬曲線與實驗曲線進行對比，結(jié)果表明：本文提出的改進統(tǒng)計損傷本構(gòu)模型是正確的，能夠合理地模擬出樁土界面的力學行為。在擬合參數(shù)的過程中，考慮到剪切實驗曲線中的極限相對位移Δu基本不變，本文統(tǒng)一取Δu為2.5mm。但是，這并不代表擬合時不同法向應力下的Δu必須為定值，在擬合不同土性的剪切實驗曲線時，若Δu變化較大，就應該為不同的值。

[1] 劉蓓． 紅黏土與混凝土結(jié)構(gòu)接觸面力學特性的試驗研究[D]． 長沙: 中南大學，2013．

LIUBei．Thetextstudyonmechanicalpropertiesofinterfacebetweenredclayandconcretestructure[D]．Changsha:CentralSouthUniversity， 2013．

[2]SeedHB,ReeseLC.Theactionofsoftclayalongfrictionpiles[J].Transaction,ASCE, 1955(122):731-754.

[3]DesaicsZ．Thinlayerelementforinterface[J]．InternationalJournalforNumericalandAnalyticalMethodsinGeomechanics， 1984，8( 1) : 19-43.

[4] 張忠苗. 樁基工程[M]. 北京: 中國建筑工業(yè)出版社, 2007.

ZHANGZhongmiao.Pilefoundationengineering[M].Beijing:ChinaBuildingIndustryPress, 2007.

[5] 蒲訶夫, 鄭俊杰, 章榮軍. 樁土界面荷載傳遞模型的改進及 其數(shù)值實現(xiàn)[J]. 華中科技大學學報(城市科學版), 2010, 27(1):84-88.PUHefu,ZHENGJunjie,ZHANGRongjun.Improvementontheloadtransfermodelofpile-soilinterfaceanditsnumericalsimulation[J].JournalofHuazhongUniversityofScienceandTechnology(UrbanScienceEdition), 2010, 27(1): 84-88.

[6] 彭凱, 朱俊高, 伍小玉, 等. 不同泥皮粗粒土與結(jié)構(gòu)接觸面力學特性實驗[J]. 重慶大學學報(自然科學版), 2011, 34(1): 110-115.

PENGKai,ZHUJungao,WUXiaoyu,etal.Experimentsoninfluenceofslurrykindsonmechanicalbehaviorofinterfacebetweengravelandconcrete[J].JournalofChongqingUniversity(NaturalScienceEdition), 2011, 34(1): 110-115.

[7] 石熊，張家生，劉蓓，等． 紅黏土與混凝土接觸面剪切特性試驗研究[J]． 中南大學學報( 自然科學版) ，2015，46( 5) : 1826-1831．

SHIXiong，ZHANGJiasheng，LIUBei，etal．Experimentalresearchonshearingpropertiesofinterfacebetweenredclayandconcrete[J]．JournalofCentralSouthUniversity(ScienceandTechnology) ，2015，46 ( 5 ): 1826-1831．

[8] 石熊, 張家生, 鄧國棟, 等. 循環(huán)荷載作用下紅黏土與混凝土接觸面剪切特性試驗研究[J]. 鐵道科學與工程學報, 2014,11(3): 88-93.

SHIXiong,ZHANGJiasheng,DENGGuodong,etal.Experimentalresearchonshearingpropertyofredclay-concreteinterfaceundercyclicloading[J].JournalofRailwayScienceandEngineering, 2014, 11(3): 88-93.

[9] 馮大闊，張建民． 粗粒土與結(jié)構(gòu)接觸面靜動力學特性的大型單剪試驗研究[J]． 巖土工程學報，2012，34 ( 7) : 1201- 1208．

FENGDakuo，ZHANGJianmin．Monotonicandcyclicbehaviorsofcoarse-grainedsoilstructureinterfaceusinglarge-scalesimplesheardevice[J].ChineseJournalofGeotechnicalEngineering，2012，34( 7) : 1201-1208．

[10] 楊林德，劉齊建． 土-結(jié)構(gòu)物接觸面統(tǒng)計損傷本構(gòu)模型[J]． 地下空間與工程學報，2006，2 ( 1) : 79-83．

YANGLinde，LIUQijian．Researchonstatisticaldamagemodelforsoilstructureinterface[J]．ChineseJournalofUndergroundSpaceandEngineering，2006，2(1) : 79-83．

[11] 夏紅春, 周國慶, 商翔宇. 基于Weibull分布的土-結(jié)構(gòu)接觸面統(tǒng)計損傷軟化本構(gòu)模型[J]. 中國礦業(yè)大學學報, 2007, 36(6): 734-737.

XIAHongchun,ZHOUGuoqing,SHANGXiangyu.StatisticaldamagesofteningconstitutivemodelofsoilstructureinterfacebasedonWeibullRandomdistribution[J].JournalofChinaUniversityofMining&Technology, 2007, 36(6): 734-737.

[12] 殷宗澤，朱泓，許國華． 土與結(jié)構(gòu)材料接觸面的變形及其數(shù)學模擬[J].巖土工程學報，1994，16( 3) : 14-22．

YINZongze，ZHUHong，XUGuohua．Numericalsimulationofthedeformationintheinterfacebetweensoilandstructuralmaterial[J].ChineseJournalofGeotechnicalEngineering，1994， 16( 3) : 14-22.

[13] 盧廷浩, 鮑伏波. 接觸面薄層單元耦合本構(gòu)模型[J]. 水利學報, 2000, 52(2): 71-75.

LUTinghao,BAOFubo.Acoupledconstitutivemodelforinterfacethin-layerelement[J].JournalofHydraulicEngineering, 2000, 52(2): 71-75.

[14] 胡黎明． 土與結(jié)構(gòu)物接觸面力學特性研究和工程應用[D].北京: 清華大學，2000．

HULiming.Mechanicalcharacteristicsofsoil-structureinterfaceandengineeringapplications[D].Beijing:TsinghuaUniversity，2000．

[15]AlonsoE,JosaA,LedesmaA.Negativeskinfrictiononpiles:Asimplifiedanalysisandpredictionprocedure[J].Geo-technique, 1984, 34 ( 3 ): 341-357.

[16] 楊麗君, 王偉, 盧廷浩. 樁-土接觸面剪切性質(zhì)室內(nèi)單剪試驗研究[J]. 公路, 2008, 8(8): 209-212.

YANGLijun,WANGWei,LUTinghao.Astudyonlaboratorysimpleshearingtestofpile-soilinterface[J].Highway, 2008, 8(8): 209-212.

[17] 陳仁朋, 周萬歡, 曹衛(wèi)平,等.改進的樁土界面荷載傳遞雙曲線模型及其在單樁負摩阻力時間效應研究中的應用[J]. 巖土工程學報, 2007, 29 ( 6):824-830.

CHENRenPeng,ZHOUWanhuan,CAOWeiping,etal.Improvedhyperbolicmodelofload-transferforPile-soilinterfaceanditsapplicationinstudyofnegativefrictionofsinglepilesconsideringtimeeffect[J].ChineseJournalofGeotechnicalEngineering, 2007, 29 ( 6):824-830.

Research on improved contact surface constitutive modelbased on statistical damage constitutive model

LISAI1,WANGYou1,QINZHIHAO1，LIUJianhua2

(1. School of Civil Engineering and Architecture, Central South University, Changsha 410075, China;2.SchoolofTrafficandTransPortationengineering,ChangshaUniversityofScienceandTechnology,Changsha410076,China)

Consideringtheinitialshearstiffnessdeptheffectofthepile-soilinterfaceandkeepingtheuniformityalongthecontactsurfacethicknessdirectionoftheshearstressandshearstrain,thepaperputforwardthezero-thicknesszonecontactsurfaceconstitutivemodelbasedonthestatisticaldamageconstitutivemodel.Accordingtotheshearexperimentdata,weobtainedthesimulationcurvesfittingparametersoftheimprovedconstitutivemodelunderdifferentnormalstresses.Calculationresultsindicatedthatthezero-thicknesszonecontactsurfaceconstitutivemodeldevelopedwascorrect,whichcouldsimulatereasonablythemechanicalbehaviorofthepile-soilinterfaceanddescribethenonlinearityofthepile-soilcontactsurfaceshearstiffness.Italsohasstrongapplicability,relativelyhightheoreticalandpracticalapplicationvalue.

pile-soilinterface;zero-thicknesszone;statisticaldamageconstitutivemodel;fittingparameters

2015-12-07

國家自然科學基金資助項目(51308552, 51408060)；湖南省高校創(chuàng)新平臺開放基金資助項目(14K105)

汪優(yōu)(1978-)，女，湖南湘潭人，副教授,博士，從事樁基礎(chǔ)及路基工程研究；E-mail：ywang1920@csu.edu.cn

TU473.1

A

1672-7029(2016)07-1247-06