鋼筋混凝土箱梁瀝青攤鋪溫度場數值分析

鄧祖華，劉其偉

(1.江蘇建研建設工程質量安全鑒定有限公司，江蘇南京 210008;2.東南大學交通學院，江蘇南京 210096)

其中包含n個線性代數方程式，從而解出n個節點溫度值。

將上述方程組寫成矩陣形式，即得到有限元法計算瞬態溫度場的基本方程為:

0 引言

熱傳導問題的數值求解方法主要有差分法和有限元法。其中有限單元法是將一個表示結構或連續體的求解域離散為若干個子域(單元)并通過它們邊界上的結點相互聯結成為組合體，對每一單元假定一個合適的、較簡單的近似解，然后推導求解域的滿足條件，從而得到問題的解，這個解不是準確解，而是近似解。下面具體闡述利用有限單元法分析瀝青攤鋪作用下混凝土箱梁的溫度場。

1 混凝土溫度場計算基本原理

1.1 導熱微分方程

箱梁內部和表面的某一點，在某一瞬間的溫度Ti，可以表示為:

該點的溫度Ti不僅與坐標 x，y，z有關，而且與時間t有關，對于鋼筋混凝土結構可以近似當作各向均為均質、同性的固體，根據Fourier熱傳導理論，可導得下述三維非穩定導熱方程:

式中:λ為混凝土的導熱系數，W/m·℃;c為混凝土的比熱，J/kg·℃;ρ為混凝土的容重，kg/m3;q為單位體積內放出的熱量，J/m3。

當不研究材料的水化熱，即q=0時，有:

僅考慮平面溫度場則:

1.2 初始條件和邊界條件的確定

熱傳導方程建立了物體溫度與時間、空間的關系，但滿足熱傳導方程的解有無限多，還應該給出解決這一具體問題解的附加條件，稱為定解條件。對于非穩態的導熱問題，定解條件有兩個方面，即給出初始時刻溫度分布的初始條件，以及給出箱梁邊界上溫度或換熱情況的邊界條件。

1.2.1 導熱方程的初始條件

為準確對瀝青攤鋪橋梁結構進行分析，取攤鋪到試驗橋梁前橋梁結構各位置的溫度值。

對于平面則:

上式表示的是在初始時刻結構內部溫度場分布為其截面坐標的函數或為常數值，與時間參數無關。

1.2.2 導熱方程的邊界條件

鑒別的時候需要確認魚腥味的程度、純凈度。魚粉應該有較為強烈的魚腥味，但也因原料種類、加工方式及脂肪含量等差別而存在一定的差異，鑒定人員需要鑒別魚腥味中的細微差異。魚腥味的純凈度則是指是否含有魚腥味之外的其他味，非魚腥味應該越少越好。

導熱問題常見的邊界條件可以歸納為以下三類:

1)規定了邊界上的溫度值，稱為第一類邊界條件，它的表達式如下:

式中:TΓ為物體的邊界;Tw為已知的邊界溫度值，℃;f1(x，y，t)為邊界上的溫度函數，隨位置和時間變化。

2)規定了邊界上的熱流密度值，稱為第二類邊界條件。此類邊界條件最簡單的典型例子就是規定邊界上的熱流密度保持特定值。對于非穩態導熱，這類邊界條件可以給出以下關系式:

式中:f2(x，y，t)為邊界上的熱流密度，隨位置和時間變化。

3)規定了邊界上物體與周圍流體間的表面傳熱系數β及周圍流體的溫度Tf，稱為第三類邊界條件。第三類邊界條件可表示為:

瞬態導熱時，式中β及Tf均為時間的函數。

1.3 箱梁溫度場的非線性有限元解法

采用基于微分方程等效積分原理的Galerkin加權余量法，對熱傳導平面方程建立微分方程。

1.3.1 有限單元法熱傳導方程的建立及邊界條件的處理

利用加權余量法及格林公式，可得:

式(11)即為平面溫度場有限單元法計算的基本方程，對于其中的線積分項，可把邊界條件代入，從而滿足邊界條件。

對于第一類邊界條件，上式中的附加線積分項等于0;

對于第三類邊界條件，附加線積分項等于

1.3.2 有限單元法的總剛度矩陣的形成

有限元法的計算最終結果是要求出D中的溫度分布，把域D剖分成E個四邊形單元和n個節點，現在把單元對總體的貢獻綜合起來，最終把各個節點的溫度值求出。

單元總體合成如下:

其中包含n個線性代數方程式，從而解出n個節點溫度值。

將上述方程組寫成矩陣形式，即得到有限元法計算瞬態溫度場的基本方程為:

1.4 溫度場的求解

輸入各種相關數據，劃分單元網格后，依次計算每個單元的熱傳導矩陣、熱容矩陣和荷載列陣，并逐步疊加可形成總的熱傳導陣[K]，熱容矩陣[C]和熱傳導荷載陣{P}。在時間域上進行有限差分就可以初步迭代計算出每一時刻的溫度場。采用向后差分法代入式(13)，化簡可得:

輸入初始溫度場{T}0后，逐步求解式(14)線性方程組，即可求出t時刻各節點的溫度值。

2 箱梁溫度場的理論計算

ANSYS有限元計算軟件熱分析模塊基于以上理論，本文分別采用第一邊界條件及第三邊界條件對箱梁截面進行溫度場分析。

2.1 第一邊界條件求解二維溫度場

采用第一邊界條件，對測試斷面溫度場進行分析。瀝青表層溫度和橋梁結構表面溫度用溫度槍實測獲得。

圖1為根據荷載步列出的主要時刻的橋梁結構溫度分布云圖。

根據理論計算結果做出特定時刻的溫度梯度曲線，并與實測曲線進行比較，見圖2。

從圖2可知，理論計算結果與實測數值吻合情況良好，特別是各個測點所達到的最大溫度值非常相近，最大不超過0.8℃，可以認為理論分析結果能滿足使用精度的要求。

2.2 第三邊界條件求解二維溫度場

采用第三類邊界條件對測試斷面溫度場進行分析。結合測試斷面實測得到的數據和有限元分析結果，對實測值與理論值進行分析比較，以證明理論計算的合理性和正確性，見圖3。

從圖3可以看出，使用第三類邊界條件處理橋梁結構瀝青攤鋪過程的邊界條件與實測結果是吻合的，用第三類邊界條件是可以很好地對模型進行分析，并得到正確的結果。

續圖3 距梁頂不同距離處的溫度實測值與理論值

3 結論

在理論求解定義的初始、邊界條件和材料特性熱工參數準確的前提下，是可以通過理論計算得到符合實際的箱梁溫度場的。

采用第一類邊界條件和第三類邊界條件對試驗橋進行數值分析，兩種方法的計算結果均與試驗結果吻合程度良好。

［1］王勖成.有限單元法［M］.北京:清華大學出版社，2003.

［2］彥啟森，趙慶珠.建筑熱工程［M］.北京:建筑工業出版社，1986.

［3］張朝暉.ANSYS8.0熱分析教程與實例解析［M］.北京:中國鐵道出版社，2005.