新型預制U型渠結構的數值模擬分析

孫景路，李欣欣，王正君*，張濱，常俊德

(1.黑龍江省水利科學研究院 黑龍江省季節凍土區工程凍土重點實驗室，哈爾濱 150080； 2.黑龍江大學 水利電力學院，哈爾濱 150080)

新型預制U型渠結構的數值模擬分析

孫景路1，李欣欣2，王正君2*，張濱1，常俊德1

(1.黑龍江省水利科學研究院 黑龍江省季節凍土區工程凍土重點實驗室，哈爾濱 150080； 2.黑龍江大學 水利電力學院，哈爾濱 150080)

本文以水利部“948”計劃項目為背景，在滿足抗凍性能不低于F300，并在給定的混凝土強度等級、受力鋼筋強度等級不變的前提下，借助大型有限元分析軟件ABAQUS對矩型渠Z500、Z1000進行優化設計及數值模擬分析，結合模擬計算結果分析了結構在給定工況下的受力狀態，數值模擬的計算結果表明：在承載力、結構的正常使用及破壞狀態、配筋配置以及經濟等方面優化結構所具有的優勢。同時也證明了本工程優化后的方案具有應用的價值。

預制混凝土矩型渠；數值模擬；渠道防滲；有限元模型

0 引言

渠道防滲是諸多農田灌溉節水措施中經濟合理、技術可行的主要節水措施之一，也是目前我國應用最廣泛的節水工程技術措施[1]，U型渠槽是混凝土渠道防滲中的一種防滲層結構形式，一般采用混凝土現澆或預制而成[2]。混凝土U型槽預制件是90年代開發出來的一種供農業輸水溝道建設需要的特殊產品[3]。預制U型槽新的施工工藝使組件配套安裝更加方便，質量也更容易控制，具有較高的性價比[4]。

U型渠道不僅接近最優輸水斷面，還具有省水、省地、省料和省土等特點，因此在水利工程中被廣泛的采用，這種結構形式也常被應用在中小型閘、壩、涵洞、渡槽和跌水等水工建筑物中[5-6]。目前，中國有許多建于早期的小型農田水利工程受當時的經濟和技術條件所限，大部分渠道均為土渠，長達半個世紀的運行且年久失修，淤塞、倒塌、滲漏嚴重，渠系水利用系數僅為0.3～0.5[7]。在有些抗凍要求較高的寒區，有的U型槽相關指標很難到達使用要求[8]。但有些渠系經過結構改造以后，既提高了灌區水利用系數，也推動了灌區水利事業的發展[9]。因此，應用和推廣U型混凝土預制槽防滲襯砌技術十分適宜且很有必要[10]。

本文已對渠系結構原設計(填方土渠道方案)進行了進一步優化。本文借助大型有限元分析軟件ABAQUS對矩型渠Z500、Z1000進行數值模擬分析與比較，以驗證結構的可靠性和優化方案的可行性。

1 結構的基本情況

1.1 材料

本文擬推薦采用水利部“948”項目《先進渠系建筑物制造關鍵技術引進》新產品——裝配式混凝土組裝式矩形槽替代原規劃設計填方渠道設計方案。

該產品采用德國引進的海斯混凝土成型先進技術設備制造，綜合造價等同國內現有混凝土U型槽，具有相對于填方渠道占地少、土地利用系數高、易于安裝、省時省工，有效提高灌溉水利用系數，增加灌溉面積和耕地面積等優點，是國內傳統混凝土U型槽、梯形槽升級換代產品。

采用C50的混凝土，抗凍性能不低于F300[11-13]；采用HPB235的鋼筋，地基土選用黏土。參照《SL 191-2008水工混凝土結構設計規范》得到相關參數密度為2 000 kg/m3，彈性模量為20MPa，泊松比0.3，摩擦角32°，膨脹角20.5°，粘聚力71 600 Pa。

型號規格斷面尺寸abcdefgrlZ500500×500500460500608080620802000Z10001000×1000100092010008012012011601002000

1.2 工況及建筑等級

本次研究發現，相比對照組，研究組治療有效率對比無顯著差異，P＞0.05，且研究組患者生活質量指標改善情況更為明顯，P＜0.05。分析原因：在MRI影像輔助下實施夾脊穴深刺主要具有以下優勢：首先，可在MR成像輔助下，對病灶位置及腰椎情況進行有效觀測，并可實現對夾脊穴的準確定位，確保治療準確性；其次，通過MRI成像技術，可對治療中針刺深度實現準確把握，在適宜穴位深度實施治療，確保治療有效性[4-5]。

工況1：不考慮凍脹力，渠內滿水；工況2：考慮凍脹力，地基土的凍深為2 m；工況3：考慮冰壓力，冰厚0.3 m；工況4：考慮渠外側填土壓力，渠內滿水。水工建筑物等級為4級。

2 矩型渠的數值模擬分析

ABAQUS 是一套工程模擬的有限元軟件，其解決問題的范圍從相對簡單的線性分析到許多復雜的非線性問題。本文將借助大型有限元分析軟件ABAQUS，對Z500、Z1000進行數值模擬分析，以期得到結構在給定工況下的受力狀態，驗證結構的可靠性。

2.1 構建模型

鋼筋采用的是理想的彈塑性模型。混凝土的模型是根據《混凝土結構設計規范(GB 50010-2010)》得到混凝土單軸受拉、單軸受壓的本構關系，經過計算擬合得到C50混凝土受拉、受壓的應力-應變本構關系模型。地基土則選用Mohr-Coulomb模型。

2.2 建立有限元模型

以Z500矩型渠為例。

(1)幾何模型

結構的部分取2 m長的Z500矩型渠為計算單元。

(2)有限元模型

土體、矩型渠均采用實體單元，網格劃分采用8節點減縮積分單元C3D8R。網格劃分時對渠底圓弧拐角處進行了細化，最終矩型渠結構部分共有4 500個單元；土體部分共有10 000個單元，有限元模型如圖3所示。

(3)荷載及約束條件

矩型渠所受荷載包括自重載荷、渠內水壓力渠兩側土的凍脹力。土體所受荷載包括自重載荷、頂部所受土壓力。因為當結構方案滿足工況1的設計要求時，工況3和工況4自然滿足，故本文的數值分析只考慮工況1和工況2。模型外部約束均為面的法向位移約束。

(4)接觸設置

矩型渠底面與地基土采用面對面接觸，定義接觸時主面為混凝土，從面為土體。根據《建筑地基基礎設計規范(GB50007-2002)》，庫倫摩擦系數為0.3～0.35，本文取0.3。

2.3 計算結果

本文給出結構在工況1和工況2情況下的分析計算結果。以UJ500矩型渠為例，結果匯總見表2。矩形渠在工況1條件下的第一主應力和鋼筋軸向應力模擬如圖5和圖6所示。

表2 混凝土最大應力

注：C50混凝土的抗拉、抗壓強度設計值分別為ft=1.89 MPa、fc=23.1 MPa。

矩形渠在工況2條件下的第一主應力和鋼筋軸向應力模擬如圖7和圖8所示。

圖2 矩型渠的幾何模型Fig.2 Geometric model

圖3 矩型渠有限元模型Fig.3 Integral finite element model

圖4 原方案和優化方案的鋼筋接觸圖Fig.4 Steel contact diagram of the original program and the optimized scheme

圖5 素混凝土、原方案與優化方案的第一主應力Fig.5 The first principal stress of plain concrete，original scheme and optimization scheme

圖6 原方案與優化方案的鋼筋軸向應力Fig.6 Axial stress of reinforcement of original scheme and optimization scheme

圖7 素混凝土、原方案與優化方案的第一主應力Fig.7 The first principal stress of plain concrete，original scheme and optimization scheme

圖8 原方案與優化方案的鋼筋軸向應力Fig.8 Axial stress of reinforcement of original scheme and optimization scheme

3 結論

在給定的三種工況下：

(1)承載能力方面，兩種型號矩型渠的承載能力滿足各自的承載力使用功能要求。

(2)正常使用狀態方面，兩種型號矩型渠在工況1、工況3情況下結構不會出現裂縫。在工況2情況下，對應于兩種型號的矩型渠，混凝土最大壓應力都遠小于相應強度設計值。只有Z1000的混凝土最大拉應力接近于其抗拉強度設計值，且最大拉應力主要分布在矩型渠底部的下表面區域附近，此處矩形渠底會出現裂縫，但沒有向底板厚度方向伸展，能滿足矩型渠正常使用要求。

(3)鋼筋配置方面，三種工況下，原方案及優化方案結構中鋼筋的應力都遠小于其強度設計值，表明鋼筋的強度沒能得到充分發揮且原方案的配筋遠遠不滿足現行設計規范的“適筋破壞”構造要求。

(4)結構的破壞狀態，優化后結構的延性得到明顯提高。由于不滿足設計規范的“適筋破壞”構造要求，原方案的結構具有明顯的脆性特征，一旦超載將發生不可修復的破壞。

(5)經濟效益方面，從制造成本方面考慮，素混凝土及原方案結構的制造成本要遠遠低于優化后的結構。

[1] 周維博，李立新，何武權，等.我國渠道防滲技術研究與進展[J].水利水電科技進展，2004，24(5)：60-63.

[2] 施全祖.現澆 U 型渠槽在渠道防滲工程中的推廣應用[J].甘肅農業，2004(6)：85-86.

[3] 李光正.砼U型槽預制件在山區引水溝渠建設中的應用[J].云南農業，1997(12)：15

[4] 劉道維，孫景路，張濱.渠道防滲新型 U(矩)型槽施工技術探討[J].森林工程，2015，31(5)：130-133.

[5] 王中海.混凝土 U 型渠在中小型灌區中的應用[J].現代農業科技，2011(9)：261-262.

[6] 邵文華，劉鳳林.U型槽防凍害的有效措施[J].吉林水利，2004，12(4)：35-36

[7] 劉木蘭.淺談預制混凝土U型槽在重點縣工程建設中的應用[J].科技論壇，2010(8)：9-10

[8] 孫景路，張濱，張守杰，等.渠道防滲新型結構體系的開發與應用[J].森林工程，2015，31(4)：158-160.

[9] 郭志龍.預制混凝土 U 型渠在釜山灌區渠系改造中的應用[J].山西水利，2013，29(9)：29-30.

[10] 王海霞.“U”型混凝土預制槽在渠道防滲襯砌工程中的應用[J].水利水電建設，2013(5)：57-58

[11] 中華人民共和國水利部.GB50288—99灌溉與排水渠系 工程設計規范[S].北京：水利水電出版社，1999.

[12] 中華人民共和國水利部.SL482—201l灌溉與排水渠系建筑物設計規范[S].北京：中國水利水電出版社，2011.

[13] 中華人民共和國水利部，SL/T50363—2006節水灌溉技術規范[S].北京：中國水利水電出版社，2006.

NumericalSimulationAnalysisofNewPrecastUChannel

Sun Jinglu1，Li Xinxin2，Wang Zhengjun2*，Zhang Bin1，Chuang Junde1

(1.Heilongjiang Water Conservancy Institute， Heilongjiang Seasonal Frozen Soil Region Engineering Frozen Soil Key Laboratory，Harbin 150080； 2.Institute of Hydraulic and Electric Engineering，Heilongjiang University，Harbin 150080

Based on the “948” project of the Ministry of Water Resources，when the frost resistance performance is not less than F300，and under the condition of unchanged strength grade of concrete and reinforced steel bar，the optimal design and numerical simulation of U channel Z500，Z1000 is carried out by finite element analysis software ABAQUS.The stress state of the structure under given conditions is analyzed combined with the simulation results.The numerical simulation results show that the optimization structure has advantages in the bearing capacity，the normal use of the structure and the destruction of the state，reinforced steel bar configuration and economy.At the same time，it proves that the optimization scheme of this project has the value of application.

Precast concrete U channel；numerical simulation；canal seepage control；finite element model

S 277； TV 91

：A

：1001-005X(2017)05-0113-04

2017-05-05

水利部“948”計劃項目(201316)

孫景路，碩士，高級工程師。研究方向：工程凍土和寒區水工技術研究。E-mail：sunjinglu74@163.com

王正君，博士，教授。研究方向：建筑材料及結構檢測。E-mail：wzjsir@163.com

孫景路，李欣欣，王正君，等.新型預制U型渠結構的數值模擬分析[J].森林工程，2017，33(5)：113-116.