吳琦，張戈，張碩，肖先煊

(1．華北水利水電大學(xué)，河南 鄭州 450045;2．成都理工大學(xué) 地質(zhì)災(zāi)害防治與環(huán)境保護(hù)國家重點(diǎn)實(shí)驗(yàn)室，四川 成都 610059)

大體積混凝土隨著溫度的變化會(huì)產(chǎn)生微量形變，由于其自身結(jié)構(gòu)體系的復(fù)雜性和所處環(huán)境的多樣性，必須充分考慮溫度形變場對(duì)大體積混凝土產(chǎn)生形變的影響［1］．美國加州大學(xué)的E．L．Wilson［2］最早使用有限元時(shí)間過程分析方法計(jì)算了Dworkshak大壩的溫度形變場． 韓國世宗大學(xué)的Jang-Ho Jay Kim 等［3］應(yīng)用有限元法及差分法計(jì)算了壩體的溫度形變場．在英國和法國，ADINA、ABAQUS、ANSYS 等有限元分析軟件運(yùn)用比較廣泛［4］． 朱伯芳院士和宋敬廷合作編制了我國第一個(gè)混凝土溫度形變有限元分析程序，把有限單元法運(yùn)用于大體積混凝土結(jié)構(gòu)溫度形變場的計(jì)算分析中［5］．

為了進(jìn)一步分析大體積混凝土樁板墻的微變與溫度形變場之間的相關(guān)性，本文基于IBIS －L 地形微變遠(yuǎn)程監(jiān)測系統(tǒng)，監(jiān)測樁板墻隨溫度變化的實(shí)時(shí)形變量，利用有限元理論方法計(jì)算樁板墻產(chǎn)生的溫度形變，并通過ANSYS 有限元軟件進(jìn)行仿真模擬，將溫度形變場引入大體積混凝土的計(jì)算模型中，并結(jié)合實(shí)例總結(jié)兩者之間的對(duì)應(yīng)關(guān)系．

1 IBIS－L 系統(tǒng)變形監(jiān)測分析

1．1 IBIS－L 系統(tǒng)簡介

IBIS－L 地形微變遠(yuǎn)程監(jiān)測系統(tǒng)是集合了步進(jìn)頻率連續(xù)波技術(shù)(SF － CW)、合成孔徑雷達(dá)技術(shù)(SAR)和干涉測量技術(shù)的新型高級(jí)監(jiān)測系統(tǒng)． 該系統(tǒng)通過合成孔徑技術(shù)來獲取監(jiān)測區(qū)域的二維圖像，并將監(jiān)測區(qū)域劃分為許多小單元，每個(gè)小單元都代表著一個(gè)監(jiān)測點(diǎn)的變形，然后應(yīng)用干涉技術(shù)提取其相位變化量［6］．該系統(tǒng)在人造衛(wèi)星和陸地激光掃描的基礎(chǔ)上對(duì)可靠性、精度予以改進(jìn)，可應(yīng)用于對(duì)建筑物和地形的靜態(tài)監(jiān)測． 其最大監(jiān)測距離可達(dá)4 km，測量精度可達(dá)0．1 mm，最小圖像采集時(shí)間為5 min，能夠24 h 不間斷地監(jiān)測整個(gè)目標(biāo)區(qū)域的位移變化量［7］．

1．2 IBIS－L 系統(tǒng)微變監(jiān)測

研究區(qū)樁板墻在每日溫度變化的影響下，整體產(chǎn)生微變位移．應(yīng)用IBIS －L 系統(tǒng)對(duì)樁板墻的微變形量進(jìn)行全方位連續(xù)監(jiān)測，監(jiān)測頻率為1 次/h，周期為20 d．

為了減小監(jiān)測時(shí)的誤差，需要在整個(gè)監(jiān)測區(qū)域范圍內(nèi)選取穩(wěn)定的參照物作為地面穩(wěn)定控制點(diǎn)(Ground Control Points，GCP)，也就是整個(gè)監(jiān)測范圍的基準(zhǔn)點(diǎn)．本次監(jiān)測在儀器位置前方的居民房屋附近選取2 個(gè)地面穩(wěn)定控制點(diǎn)，并在樁板墻上由南向北選取8 個(gè)變形控制點(diǎn)進(jìn)行定點(diǎn)分析． 監(jiān)測時(shí)規(guī)定以IBIS－L 系統(tǒng)的安裝位置為原點(diǎn)，儀器行進(jìn)方向?yàn)閤 軸，與之垂直的雷達(dá)波發(fā)射方向?yàn)閥 軸，GCP點(diǎn)、變形控制點(diǎn)及其位移云圖如圖1 所示． 圖中:正值表示背離儀器方向;負(fù)值表示靠近儀器方向．

根據(jù)IBIS－L 系統(tǒng)的監(jiān)測結(jié)果，樁板墻20 d 內(nèi)每日發(fā)生的微變與溫度變化趨勢幾乎一致，其相對(duì)趨勢曲線如圖2 所示．

從圖2 中可以看出:①8 個(gè)變形控制點(diǎn)的變化趨勢相近，整體趨勢一致;②樁板墻回彈變形的位移曲線和溫度變化趨勢幾乎一致;③在溫度形變場的影響下，監(jiān)測到的變形回彈位移在4 mm 的范圍內(nèi)波動(dòng)．

圖1 GCP 點(diǎn)、變形控制點(diǎn)及其位移云圖

圖2 樁板墻微變位移與溫度變化趨勢圖

2 大體積混凝土樁板墻溫度形變場的有限元計(jì)算

2．1 溫度形變場熱傳導(dǎo)方程及第三類邊界條件

假設(shè)混凝土為均勻的、各向同性的固體［8］，其熱傳導(dǎo)滿足以下微分方程:

式中:T 為溫度，℃;τ 為時(shí)間，h;a 為導(dǎo)溫系數(shù)，m2/h;θ 為混凝土的絕熱溫升，℃．

第三類邊界條件是大體積混凝土溫度形變場計(jì)算最常用的邊界條件． 混凝土與空氣之間存在熱交換［9］，經(jīng)過混凝土表面的熱流量q 為:

式中:λ 為導(dǎo)熱系數(shù)，kJ/(m·h·℃);n 為邊界外法線方向．

第三類邊界條件假設(shè)經(jīng)過混凝土表面的熱流量q 與混凝土表面溫度T 和氣溫Ta之差成正比，即

式中β 為表面放熱系數(shù)，kJ/(m·h·℃)．

2．2 溫度形變場的有限單元法計(jì)算

有限單元法是比較成熟的分析大體積混凝土溫度場和形變場的數(shù)值計(jì)算方法，需要在空間區(qū)域內(nèi)對(duì)結(jié)構(gòu)體進(jìn)行離散化，同時(shí)也必須在時(shí)間區(qū)域內(nèi)進(jìn)行離散化．當(dāng)結(jié)構(gòu)體的各部分溫度發(fā)生變化時(shí)，結(jié)構(gòu)體將由于熱變形而產(chǎn)生線形變，為α(T－T0)．

大體積混凝土樁板墻由于不均勻溫度場而產(chǎn)生線形變，未產(chǎn)生剪切形變［10］． 這種形變可以看作是結(jié)構(gòu)的初形變?chǔ)?，由ε0可以求得等效節(jié)點(diǎn)溫度荷載Pε0，然后按通常有限元方法求得節(jié)點(diǎn)徑向位移δ．

對(duì)于三維單元e，有

式中:α 為材料的熱膨脹系數(shù)，1/℃;T0為結(jié)構(gòu)的初始溫度場，℃．

等效節(jié)點(diǎn)溫度荷載Pε0的計(jì)算公式為:

式中:B 為單元幾何矩陣;D 為彈性矩陣;ε0為初應(yīng)變矩陣．

根據(jù)虛位移原理，得出增量形式的等效節(jié)點(diǎn)溫度荷載有限元方程，即

式中:K 為混凝土樁板墻的整體剛度矩陣;ΔPε0為溫度增量引起的等效節(jié)點(diǎn)荷載．

研究場地樁板墻高20 m，樁體寬1．5 m，厚2 m，樁間板寬2 m．混凝土強(qiáng)度等級(jí)為C30，混凝土的熱學(xué)參數(shù)見表1．

表1 混凝土的熱學(xué)參數(shù)

由式(6)可求得樁板墻上8 個(gè)變形控制點(diǎn)由溫度形變場引起的節(jié)點(diǎn)徑向位移． 計(jì)算結(jié)果表明，8 個(gè)點(diǎn)的徑向位移基本相同，與實(shí)測位移的變化趨勢基本一致．以變形控制點(diǎn)P4為例，其在20 d 內(nèi)的實(shí)測位移和徑向位移變化趨勢如圖3 所示．由圖3 可見，樁板墻隨著溫度的變化會(huì)產(chǎn)生不同的溫度形變，由溫度引起的徑向位移的計(jì)算值的誤差不超過2%．

圖3 樁板墻P4 點(diǎn)實(shí)測位移與徑向位移變化趨勢圖

3 ANSYS 溫度形變場模擬結(jié)果分析

樁板墻P4點(diǎn)第1 天和第20 天的溫度分布云圖如圖4 所示，不同時(shí)期的溫度形變?nèi)鐖D5 所示．

圖4 樁板墻P4 點(diǎn)的溫度分布云圖

圖5 樁板墻P4 點(diǎn)的溫度形變分布圖

從圖5 中可以看出，樁板墻P4點(diǎn)在20 d 的監(jiān)測周期內(nèi)，當(dāng)日溫差變化較大時(shí)，其形變隨之增大;反之則亦然．如第1 天溫差為11 ℃，第20 天溫差為9 ℃，二者隨之產(chǎn)生的溫度形變相差0．2 mm． ANSYS 有限元軟件數(shù)值仿真模擬的溫度形變在第1天、第5 天、第15 天、第20 天時(shí)分別為3．61、2．63、2．96、3．41 mm，這與IBIS －L 系統(tǒng)監(jiān)測的實(shí)際數(shù)據(jù)基本相同．

4 結(jié) 語

1)IBIS－L 地形微變遠(yuǎn)程監(jiān)測系統(tǒng)可以全方位精確地監(jiān)測研究區(qū)樁板墻每天隨溫度變化而產(chǎn)生的微變位移，監(jiān)測數(shù)據(jù)形成的微變位移曲線與實(shí)時(shí)溫度變化曲線的變化趨勢一致．

2)ANSYS 有限單元軟件能夠較好地仿真大體積混凝土樁板墻受溫度形變場影響而產(chǎn)生溫度形變的過程，其計(jì)算和模擬出來的溫度形變與IBIS －L 系統(tǒng)實(shí)際監(jiān)測的微變位移基本相同． 但是由于使用ANSYS 有限單元軟件數(shù)值模擬樁板墻溫度形變時(shí)對(duì)彈性模量、熱力學(xué)公式等參數(shù)的選擇過于理想化，導(dǎo)致數(shù)值分析的結(jié)果與現(xiàn)場實(shí)測數(shù)據(jù)之間存在一定的誤差．

3)通過對(duì)現(xiàn)場20 d 的連續(xù)監(jiān)測與分析，8 個(gè)變形控制點(diǎn)對(duì)應(yīng)的大體積混凝土樁板墻的溫度形變有如下規(guī)律:連續(xù)高溫或連續(xù)低溫時(shí)，其產(chǎn)生的變形變化量較小，曲線走勢較平穩(wěn);連續(xù)高低溫交錯(cuò)，即相對(duì)溫差較大時(shí)，產(chǎn)生的變形變化量較大，曲線走勢會(huì)出現(xiàn)突增現(xiàn)象．

［1］朱伯芳． 大體積混凝土溫度應(yīng)力與溫度控制［M］． 北京:中國電力出版社，1999:1 －3．

［2］Malkawl Abdallah I Husein，Mutasher Saad A，Qiu Tony J．Thermal-structural modeling and temperature control of roller compacted the concrete gravity dam［J］． Journal of Performance of Constructed Facilities，2003，17(4):177－187．

［3］Jang-Ho Jay Kim，Sang-Eun Jeon，Jin-Keun Kim． Development of new device for measuring thermal stresses［J］．Cement and Concrete Researsh，2002，32(10):1645－1651．

［4］張朝暉．ANSYS12．0 熱分析工程應(yīng)用［M］． 北京:中國鐵道出版社，2010:1 －101．

［5］胡少偉，魯文妍．基于XFEM 的混凝土三點(diǎn)彎曲梁開裂數(shù)值模擬研究［J］． 華北水利水電大學(xué)學(xué)報(bào)(自然科學(xué)版)，2014，35(4):48 －51．

［6］邱志偉，汪學(xué)琴．地基雷達(dá)系統(tǒng)IBIS－L 在大壩變形監(jiān)測中的應(yīng)用［J］．長江科學(xué)院院報(bào)，2014，31(10):104－107．

［7］王峰，許強(qiáng)，肖先煊，等．基于地形微變遠(yuǎn)程監(jiān)測系統(tǒng)的滑坡變形監(jiān)測［J］．長春工程學(xué)院學(xué)報(bào)(自然科學(xué)版)，2014，15(1):66 －71．

［8］康明睿，薛偉辰．配FRP 箍筋混凝土梁的抗剪性能研究進(jìn)展［J］．華北水利水電大學(xué)學(xué)報(bào)(自然科學(xué)版)，2015，36(2):15 －20．

［9］劉亞基，蘇志敏．大體積混凝土澆筑塊溫度應(yīng)力場直接耦合法仿真計(jì)算［J］． 科學(xué)技術(shù)與工程，2010，10(16):57 －60．

［10］黃君寶． 大體積混凝土溫度場應(yīng)力場仿真方法研究［D］．天津:天津大學(xué)，2009:21 －22．