程玲 董景剛

摘要：土工結(jié)構(gòu)數(shù)值模擬中土與結(jié)構(gòu)接觸問(wèn)題一直是學(xué)者研究的熱點(diǎn)。通過(guò)改變普通四邊形單元等效線性黏彈性動(dòng)力本構(gòu)模型中的動(dòng)剪模量和模量指數(shù)，提出一種改進(jìn)薄層接觸單元，其與普通四邊形單元本構(gòu)模型一致，物理意義更為明確，模型計(jì)算參數(shù)更易取得，更適于土與結(jié)構(gòu)接觸面的計(jì)算。對(duì)壩高為200m的面板堆石壩面板與墊層之間設(shè)置改進(jìn)的薄層接觸單元、不加接觸單元、采用無(wú)厚度的Goodman單元及改變薄層接觸單元?jiǎng)恿?shù)等多種工況時(shí)，面板的動(dòng)力響應(yīng)結(jié)果分析比較表明：面板與墊層之間設(shè)置改進(jìn)的薄層接觸單元時(shí)面板的動(dòng)力響應(yīng)介于不設(shè)接觸單元和設(shè)無(wú)厚度的Goodman單元之間，合理選取薄層接觸單元?jiǎng)恿?shù)可以得到理想的結(jié)果。

關(guān)鍵詞：Goodman單元；改進(jìn)薄層單元；動(dòng)力響應(yīng)；接觸面；面板堆石壩

中圖分類號(hào)：TV641.4+3 文獻(xiàn)標(biāo)志碼：A doi：10.3969/j.issn.1000-1379.2018.01.018

水電工程中土與剛性結(jié)構(gòu)之間不存在真實(shí)的過(guò)渡體，由于兩種接觸材料物理力學(xué)性質(zhì)差異巨大，外部荷載作用往往會(huì)使接觸面兩側(cè)結(jié)構(gòu)產(chǎn)生較大的變形差異，并且數(shù)值模擬時(shí)在接觸面處往往會(huì)產(chǎn)生比實(shí)際更大的剪應(yīng)力，導(dǎo)致計(jì)算結(jié)果不盡合理，因此數(shù)值模擬時(shí)在兩種結(jié)構(gòu)接觸面上需設(shè)置特殊的接觸單元進(jìn)行過(guò)渡，以得到更準(zhǔn)確合理的計(jì)算結(jié)果[1]。目前常采用無(wú)厚度Goodman單元和薄層接觸單元等特殊單元來(lái)模擬接觸面之間的相互嵌入、脫離、相對(duì)位移等。常用的接觸單元本構(gòu)關(guān)系都是基于各種假設(shè)建立的，因此接觸單元類型及參數(shù)的選取和應(yīng)用對(duì)于數(shù)值計(jì)算結(jié)果的準(zhǔn)確性極為重要。本文以200m高面板堆石壩為例，分別在面板與墊層之間不設(shè)任何接觸單元、設(shè)置無(wú)厚度的Goodman單元、設(shè)置改進(jìn)薄層接觸單元來(lái)模擬兩者之間的接觸，計(jì)算不同工況下面板壩的動(dòng)力響應(yīng)，并對(duì)計(jì)算結(jié)果及薄層接觸單元參數(shù)敏感性進(jìn)行探討。

1 常用接觸單元簡(jiǎn)介

1.1 無(wú)厚度Goodman單元

該單元為無(wú)厚度一維單元，由4節(jié)點(diǎn)組成。由于其物理意義較為明確且無(wú)厚度，可以較好地反映材料性質(zhì)差異較大的兩種結(jié)構(gòu)接觸面上的張裂和相對(duì)滑動(dòng)，因此得到了廣泛應(yīng)用。然而，只有當(dāng)兩種材料之間的接觸面較為光滑時(shí)，剪切破壞才發(fā)生在這兩種材料的接觸面上。對(duì)于絕大多數(shù)相對(duì)粗糙的兩種材料接觸時(shí)，破壞面一般發(fā)生在土體內(nèi)部，此時(shí)用無(wú)厚度的Goodman單元進(jìn)行數(shù)值模擬存在一定的缺陷。鑒于此，李守德等[2]對(duì)Goodman單元進(jìn)行了修正，使其有限元實(shí)現(xiàn)更為靈活。

1.2 Desai薄層單元

Desai薄層接觸單元[3-5]可以解決無(wú)厚度單元在受壓時(shí)接觸面兩側(cè)相互嵌人及接觸面粗糙時(shí)破壞面發(fā)生在土體內(nèi)部的情況，并且能更好地模擬接觸面的受力變形，所以被廣大學(xué)者所接受。其本構(gòu)關(guān)系為式中；t為單元厚度；σn為法向應(yīng)力；τs為切向剪應(yīng)力；knn為法向勁度系數(shù)；kss為切向勁度系數(shù)；εn為法向應(yīng)變；εs為切向應(yīng)變。

值得注意的是，Desai等在研究中[3-5]只指出t=（0.01～0.10）B（B為單元的長(zhǎng)度），并沒(méi)有給出t的具體值，所以其取值有一定的隨意性。

1.3 改進(jìn)的薄層單元

張冬霽等[6]通過(guò)一系列接觸面單剪試驗(yàn)，較為深人地分析了接觸面的破壞機(jī)理，建立了一種接觸面本構(gòu)模型，并將其用于接觸面的有限元計(jì)算中，取得了較滿意的結(jié)果。

殷宗澤等[7]對(duì)Desai薄層單元進(jìn)行改進(jìn)，改進(jìn)后的薄層接觸單元與普通單元的本構(gòu)假設(shè)、應(yīng)變場(chǎng)定義模式相統(tǒng)一，并沒(méi)有引入其他參數(shù)，只是將土體自身的變形分為土體的基本變形和破壞變形兩部分，其概念更為明確，因此具有更高的應(yīng)用價(jià)值。

2 模型及計(jì)算參數(shù)選取

2.1 計(jì)算模型

圖1為面板堆石壩計(jì)算模型。壩頂寬15m，上、下游壩坡坡比均為1：1.4，面板頂部厚度為0.3m、底部厚度按《混凝土面板堆石壩設(shè)計(jì)規(guī)范》（SL228-2013）[8]中的公式計(jì)算后取0.9m。計(jì)算模型網(wǎng)格尺寸取壩高的1/20，在混凝土面板與堆石之間依次設(shè)置5.0m厚的墊層和5.0m厚過(guò)渡層，均上下等寬。計(jì)算壩高為140m，面板在壩高為100m處設(shè)水平施工縫。

2.2 單元本構(gòu)模型及計(jì)算參數(shù)

（1）面板混凝土強(qiáng)度取C30，重度取25.0kN/m3，泊松比取0.165，彈性模量取35GPa，動(dòng)力本構(gòu)模型采用線彈性模型。

（2）壩體墊層料、過(guò)渡料、堆石料采用等效線性黏彈性動(dòng)力本構(gòu)模型[9]，壩體墊層料、過(guò)渡料、堆石料計(jì)算參數(shù)見(jiàn)表1。

（3） Goodman單元的動(dòng)力本構(gòu)模型采用河海大學(xué)的試驗(yàn)結(jié)果[10-13]，計(jì)算參數(shù)見(jiàn)表2。

接觸面的最大動(dòng)剪模量可表示為式中：σn為法向應(yīng)力。

接觸面單元的動(dòng)剪模量K與動(dòng)剪應(yīng)變?chǔ)玫年P(guān)系為式中：τf為破壞剪應(yīng)力，τf=σntanφ，其中φ為接觸面摩擦角。

接觸面單元的阻尼比λ計(jì)算公式為式中：λmax為最大阻尼比，取0.2。

（4）通過(guò)改變普通四邊形單元等效線性黏彈性動(dòng)力本構(gòu)模型中的動(dòng)剪模量K和模量指數(shù)N，提出一種特殊接觸單元。改進(jìn)薄層單元的動(dòng)剪模量K和模量指數(shù)N取值見(jiàn)表3，泊松比μ見(jiàn)表1。

2.3 計(jì)算地震波

動(dòng)力計(jì)算時(shí)地震波選用2008年汶川大地震中茂縣地震波，水平向峰值地震加速度取0.39，豎向峰值地震加速度取0.2g。動(dòng)力響應(yīng)分析時(shí)分別在面板與墊層之間設(shè)置無(wú)厚度Goodman單元、不設(shè)置任何接觸單元、厚度為10cm（計(jì)算模型網(wǎng)格尺寸的1/100）的改進(jìn)薄層單元.具體計(jì)算工況見(jiàn)表3。

3 結(jié)果分析

面板與墊層之間設(shè)置不同接觸單元，薄層單元?jiǎng)蛹裟Ａ縆分別取不設(shè)接觸單元時(shí)的動(dòng)剪模量（即K=1 900kPa/mm2）的0.50倍、0.33倍、0.25倍、0.20倍，模量指數(shù)N改變時(shí)順河向最大加速度放大倍數(shù)和豎向最大加速度放大倍數(shù)見(jiàn)表4。

從表4可以看出，順河向最大加速度放大倍數(shù)小于豎向放大倍數(shù)，設(shè)置接觸面單元工況的最大加速度放大倍數(shù)小于不設(shè)置接觸單元的，說(shuō)明若面板與墊層之間不設(shè)接觸單元，面板與墊層材料剛度的差異就會(huì)導(dǎo)致計(jì)算剛度矩陣分布不合理，從而產(chǎn)生較大的變形和動(dòng)力反應(yīng)誤差。

動(dòng)剪模量K對(duì)加速度放大倍數(shù)的影響要大于模量指數(shù)N的；隨著改進(jìn)的薄層單元?jiǎng)恿Ρ緲?gòu)模型中參數(shù)K和N取值的減小，相應(yīng)的加速度放大倍數(shù)也減小。說(shuō)明隨著動(dòng)剪模量K和模量指數(shù)N取值的減小使面板與墊層之間的材料差異得到了較好的過(guò)渡，面板上一定的剪力轉(zhuǎn)移到了墊層上，從而使得面板與墊層之間實(shí)現(xiàn)更好的變形協(xié)調(diào)。

圖2為面板與墊層之間不設(shè)接觸單元、設(shè)置無(wú)厚度Goodman單元、設(shè)置改進(jìn)的薄層單元（N=0.50，動(dòng)剪模量取0.50K）時(shí)面板順河向和豎向加速度放大倍數(shù)沿壩高的分布情況。

由圖2可以看出，加速度放大倍數(shù)隨著壩高的增高呈現(xiàn)先緩慢減小后迅速增大的趨勢(shì)，在約0.8壩高處突然大幅增大（鞭梢效應(yīng)），所以在實(shí)際工程設(shè)計(jì)施工時(shí)相對(duì)壩高超過(guò)0.8時(shí)要采取一定的抗震加固措施。

4 結(jié)語(yǔ)

在進(jìn)行土與結(jié)構(gòu)數(shù)值模擬時(shí)采用接觸單元是十分必要的。本文通過(guò)改變普通四邊形單元計(jì)算參數(shù)得到一種特殊的薄層接觸單元，其與普通四邊形單元本構(gòu)模型一致，物理意義更為明確，模型計(jì)算參數(shù)更易取得，更適于土與結(jié)構(gòu)接觸面的計(jì)算。對(duì)壩高為200m的面板堆石壩動(dòng)力響應(yīng)分析表明：面板與墊層之間不采用接觸單元時(shí)面板動(dòng)力響應(yīng)最大，采用無(wú)厚度的Goodman單元時(shí)面板動(dòng)力響應(yīng)最小，改變薄層單元的動(dòng)力參數(shù)對(duì)動(dòng)力響應(yīng)有一定的影響，當(dāng)模量指數(shù)N取0.4、動(dòng)剪模量取墊層料動(dòng)剪模量的0.2倍時(shí)，面板動(dòng)力響應(yīng)與采用Goodman單元很接近；合理選取改進(jìn)薄層接觸單元的計(jì)算參數(shù)，可以很好地模擬面板與墊層之間的過(guò)渡。

參考文獻(xiàn)：

[1]朱泓，殷宗澤.土與結(jié)構(gòu)材料接觸面性能研究綜述[J].河海科技進(jìn)展，1994（4）：1-8.

[2]李守德，俞洪良.Goodman接觸面單元的修正與探討[J].巖石力學(xué)與工程學(xué)報(bào)，2004，25（4）：2628-2631.

[3]DESAI C S，MA Y.Modelling of Joints and Interfaces Usingthe Disturbed-State Concept[J].International Journal for Nu-merical and Analytical Methods in Geomechanics，1992，16（9）：623-653.

[4]DESAI C S，RIGBY D B.Cyclic Interface and Joint ShearDevice Including Pore Pressure Effects[J].Journal ofGeotechnical and Geoenvironmental Engineering，1997，123（6）：568-579.

[5]DESAI C S，PRADHAN S K，COHEN D.Cyclic Testingand Constitutive Modeling of Saturated Sand-Concrete Inter-faces Using the Disturbed State Concept[J].InternationalJournal of Geomechanics，2005，5（4）：286-294.

[6]張冬霽，盧廷浩一種土與結(jié)構(gòu)接觸面模型的建立及應(yīng)用[J].巖土工程學(xué)報(bào)，1998，20（6）：62-66.

[7]殷宗澤，朱泓，許國(guó)華.土與結(jié)構(gòu)材料接觸面的變形及其數(shù)學(xué)模擬[J].巖土工程學(xué)報(bào)，1994，16（2）：14-22.

[8]中華人民共和國(guó)水利部.混凝土面板堆石壩設(shè)計(jì)規(guī)范：SL 228-2013[S].北京：中國(guó)水利水電出版社，2013：16-17.

[9]韓國(guó)城，孔憲京，李俊杰.面板堆石壩動(dòng)力破壞性態(tài)及抗震措施試驗(yàn)研究[J].水利學(xué)報(bào)，1990，21（5）：61-67.

[10]丁勇春，王建華，陳錦劍.接觸面單元計(jì)算方法探討[J].巖土力學(xué)，2006（S1）：211-213.

[11]胡黎明，濮家騮.土與結(jié)構(gòu)物接觸面物理力學(xué)特性試驗(yàn)研究[J].巖土工程學(xué)報(bào)，2001，23（4）：431-435.

[12]胡黎明.土與結(jié)構(gòu)物接觸面數(shù)值模擬[J].工程力學(xué)，2001（增刊）：467-471.

[13]欒茂田，武亞軍.土與結(jié)構(gòu)間接觸面的非線性彈性-理想塑性模型及其應(yīng)用[J].巖石力學(xué)，2004，25（4）：507-513.