反力墻結構受力性能有限元分析

張尹， 楊嘉胤

（1.浙江建設職業技術學院，杭州 311231；2.浙江大學建筑設計研究院有限公司，杭州 310028）

0 引言

隨著我國工程建設和結構試驗研究的快速發展，大型復雜工程結構大量涌現，對大型復雜結構進行試驗研究的需求隨之增加，高校和研究機構大型試驗室的建設投入也在增多。而土木工程試驗室中，大型預應力鋼筋混凝土反力墻是可以用于對大型復雜結構進行靜力試驗、反復加載試驗和擬動力試驗的重要試驗設施。

反力墻本身即是一種特種結構，在工作過程中，反力墻臺座和墻體上需要固定被試的結構模型和用于給結構模型加載的加載設備，因此反力墻會承擔試驗過程中結構模型和加載設備的自重和反力。因此反力墻不僅體積龐大，而且受力復雜多變。在極大的靜力或動力荷載作用下，反力墻不僅不能發生強度破壞，其變形控制也極為嚴格，才能不影響結構試驗的精度，這對反力墻的剛度提出了很高的要求。反力墻形式靈活多樣，構造復雜，配筋量大，目前尚未形成系統成熟的分析計算方法，也暫無直接適用設計的規范標準指南。

彭駿[1]針對某高校L型預應力反力墻，采用手算方法驗算反力墻頂點最大撓度和最大裂縫寬度，采用SAP2000軟件利用殼單元考慮剪切變形模擬反力墻進行非線性有限元分析，考慮混凝土的材料非線性，鋼筋考慮為彈性。王德玲等[2]利用ANSYS軟件對反力墻和臺座系統進行了靜力有限元分析，采用帶筋的Solid65單元模擬混凝土和非預應力鋼筋，采用等效載荷法考慮預應力鋼筋的作用，簡化了有限元模型，得出預應力對反力墻和臺座整體剛度和變形的影響。孫柏濤等[3]參考大型反力墻建立了小型鋼架反力墻的ABAQUS有限元模型，采用了C3D8R實體單元，材料選用Q235B鋼材且考慮為彈性。姚琳[4]運用ANSYS和SAP2000軟件對反力墻結構進行有限元分析，在ANSYS中采用整體式模型模擬混凝土和非預應力鋼筋，采用線單元模擬預應力鋼筋，對比不同配筋方案的反力墻的變形、應力、應變分布情況，利用廣廈軟件對反力墻進行結構設計。趙陽等[5]通過預埋件加工和安裝、鋼筋安裝和模板安裝多方面采取措施滿足反力墻對加載孔和表面平整度的高要求。

文中利用Midas軟件建立了反力墻的有限元分析計算模型，采用實體單元模擬混凝土，采用等效荷載法考慮預應力鋼筋的作用?？紤]計算多種復雜工況，通過查看有限元模型在各工況下分析計算所得的應力云圖，分析各工況下反力墻的內力分布規律，對反力墻結構受力特點進行分析，探討了預應力鋼筋對反力墻應力分布的影響和裂縫控制的作用，驗證了建模方法的有效性，可以為反力墻結構布置和配筋設計提供參考。

1 工程概況

文中以某高校土木工程試驗室在建的反力墻試驗系統作為參考案例。除反力墻試驗系統外，該試驗室將設置1000t反力架試驗系統，可與反力墻試驗系統配合使用，完成大型復雜結構的靜力試驗、反復加載試驗和擬動力試驗。反力墻試驗系統中包含兩片獨立反力墻，墻高分別為10m和13m。反力墻墻體底部連接大底盤加載臺座，加載臺座同時是1層地下室，為箱體結構，內設走道、房間、風道等，使空間得到充分利用。

2 有限元模型

2.1 建模與網格劃分

文中參考某高校試驗室在建反力墻試驗系統來確定有限元模型尺寸：反力墻墻高13m，墻體平面尺寸9.6m×4m，臺座高度3m，平面尺寸18.1m×9.8m。

文中采用Midas軟件對反力墻進行有限元分析計算，其建模和劃分網格的方法與ABAQUS軟件及ANSYS軟件不同。在ABAQUS軟件和ANSYS軟件中的建模一般先按照整體尺寸建立模型，然后根據需要的精度進行網格劃分。Midas軟件中建模過程相反，從點到線，從線到面，從面到體，先生成一定尺寸的實體單元，然后對實體單元在三維方向進行復制，組成需要的反力墻整體，文中中所采用實體單元尺寸大部分為100mm，加載點附近采用50mm尺寸的單元，便于精確控制加載點。Midas軟件中反力墻模型如圖1所示。

圖1 反力墻整體建模與單元劃分

2.2 材料與約束

在Midas軟件模型中，采用實體單元模擬混凝土，由于反力墻工作過程中變形很小，材料處于彈性狀態，因此模型中采用混凝土強度等級為C45，僅考慮材料的彈性階段，未輸入塑性階段曲線。反力墻中的配筋包括普通鋼筋和預應力鋼筋。文中將普通鋼筋作為安全余量，在模型中不考慮其對反力墻剛度的貢獻。預應力鋼筋的作用則通過在反力墻頂部施加豎直向下的荷載來考慮，如圖2所示。在臺座底部采用固接，約束全部自由度。

圖2 預應力鋼筋等效荷載

2.3 荷載與工況

在Midas軟件模型中，根據反力墻的工作情況，按照3種工況進行加載，如圖3所示。工況1：在12.000m標高處，施加3個1800kN的集中荷載；工況2：在12.000m標高、9.000m標高、6.000m標高、3.000m標高處，施加2個900kN的集中荷載；工況3：在12.000m標高施加2個1400kN的集中荷載，在9.000m標高施加2個1100kN的集中荷載，在6.000m標高施加2個800kN的集中荷載，在3.000m標高施加2個500kN的集中荷載。加載時，將每個集中力根據加載孔的位置再進行分攤，可以較為精確的模擬反力墻實際的受荷情況，如圖4～圖6所示。

圖3 反力墻加載工況

圖4 工況1加載點

圖5 工況2加載點

圖6 工況3加載點

3 結果分析

文中按照上述建模方法通過Midas軟件建立了反力墻臺座和墻體的有限元模型，得到了反力墻應力分布圖形。在未考慮預應力鋼筋作用的情況下如圖7～圖9所示，在3個工況中，最大拉應力都出現在墻體底部受拉側，此時可以將反力墻視為豎向放置并嵌固在臺座上的懸臂構件，墻體底部受拉側是其關鍵受力部位。當考慮了預應力鋼筋的作用時如圖10～圖12所示，雖然整體應力圖例中出現了拉應力，但是這些拉應力是出現在反力墻頂部預應力鋼筋等效的集中荷載作用點處，這些集中荷載是用于模擬預應力鋼筋的作用，并不是實際存在的加載點，所以無需考慮此處的應力集中，僅需關注在墻體底部的應力分布情況。將墻體底部的單元隔離來看如圖13～圖15，觀察其應力分布情況，可以看出，當考慮了預應力鋼筋的作用時，3種工況下在墻體底部均無拉應力的分布。預應力鋼筋的存在可以有效改善反力墻的應力分布情況，提高反力墻的承載能力，有效的避免關鍵受力部位拉應力的產生，使反力墻滿足一級裂縫控制要求。由于反力墻在工作過程中需要承擔結構模型和加載設備的自重和反力，受力復雜多變，因此反力墻的設計原則是：在極大的靜力或動力荷載作用下，反力墻不僅不能發生強度破壞，還要極為嚴格的控制變形，才能不影響結構試驗的精度。文中的有限元模型計算結果表明，在反力墻配筋設計中，預應力鋼筋的應用對于反力墻變形的控制起到了非常關鍵的作用。

圖7 工況1整體應力（未考慮預應力鋼筋作用）

圖8 工況2整體應力（未考慮預應力鋼筋作用）

圖9 工況3整體應力（未考慮預應力鋼筋作用）

圖10 工況1整體應力（考慮預應力鋼筋作用）

圖11 工況2整體應力（考慮預應力鋼筋作用）

圖12 工況3整體應力（考慮預應力鋼筋作用）

圖13 工況1墻體底部應力（考慮預應力鋼筋作用）

圖14 工況2墻體底部應力（考慮預應力鋼筋作用）

圖15 工況3墻體底部應力（考慮預應力鋼筋作用）

在文中的有限元模型中未考慮普通鋼筋的作用，工作荷載由混凝土和預應力鋼筋共同承擔，在這種情況下，預應力鋼筋發揮了足夠的作用，使墻體底部無拉應力出現，這說明采用文中的建模方法是合理有效的，可將普通鋼筋視為安全余量，使模型更加簡潔，提高了建模和計算效率。文中3個工況中荷載的位置按照反力墻施工圖中實際的加載錨孔定位進行設置，并且考慮了使反力墻關鍵部位處于受力不利狀態的工況，在簡化有限元模型提高計算效率的同時，盡量貼合實際，反映反力墻的受力狀態。

4 結語

文中利用Midas軟件建立了反力墻的有限元分析計算模型，采用實體單元模擬混凝土，采用等效荷載法考慮了預應力鋼筋的作用，通過計算得到不同工況下反力墻的應力分布情況。經過分析得到以下結論：

（1）墻體底部受拉側是反力墻關鍵受力部位。在未考慮預應力鋼筋作用的情況下，在文中所考慮的工作荷載作用下，最大拉應力都出現在墻體底部受拉側。

（2）通過模型對比證明了預應力鋼筋的存在可以有效改善反力墻的應力分布情況，提高反力墻的承載能力，預應力鋼筋的應用對反力墻裂縫控制具有關鍵作用。

（3）文中采用的建模方法簡單有效，計算效率高，可以為反力墻結構布置和配筋設計提供參考。