王文龍,張杰勝 (中鐵四局集團第一工程有限公司,安徽 合肥 230041)
瀝青路面攤鋪質(zhì)量的好壞與攤鋪機熨平板平直性存在直接關(guān)系。熨平板在攤鋪過程中要接觸150℃的瀝青混合料必然受熱膨脹使其發(fā)生變形,而攤鋪單幅的無拱路面要求熨平板必須是平直的,因此在攤鋪前應(yīng)對熨平板進行進行熱變形分析,通過設(shè)置合理的預(yù)拱度使其在攤鋪受熱變形時達到平直的要求。
攤鋪機熨平板的熱變形的計算與求解,實際上是一個熱彈性力學(xué)問題。分析熨平板的熱變形,需要根據(jù)熨平板的工作狀況、施工方法以及材料的傳熱學(xué)特性,按照熱彈性力學(xué)原理進行計算。為了確定熨平板的熱變形,需進行兩方面的計算:一是按照熱傳導(dǎo)理論,根據(jù)熨平板的熱學(xué)性質(zhì)、內(nèi)部熱源、初始條件和邊界條件,計算彈性體內(nèi)各點在各瞬時的溫度,而前后兩個溫度場之差就是彈性體的變溫;二是按照“熱彈性力學(xué)”,根據(jù)彈性體的變溫來求出熨平板各點的溫度應(yīng)變。
一般而言,熱分析可分為穩(wěn)態(tài)熱分析和瞬態(tài)熱分析兩種。其中,溫度場不隨時間而發(fā)生變化的傳熱過程稱為穩(wěn)態(tài)傳熱,反之則為瞬態(tài)傳熱。考慮到攤鋪機熨平板工作時,首先要對其進行充分的預(yù)熱,因此,可以認為熨平板的分析類型為穩(wěn)態(tài)熱分析。根據(jù)能量守恒原理,穩(wěn)態(tài)熱分析的熱平衡方程可用矩陣形式:
式中:[K]——傳導(dǎo)矩陣,包含導(dǎo)熱系數(shù),對流系數(shù)及輻射率和形狀系數(shù)
{T}——節(jié)點的溫度向量
{Q}——節(jié)點熱流率向量,包含熱生成。在穩(wěn)態(tài)熱分析時,任一節(jié)點的溫度都不隨時間變化
ANSYS在熱分析的過程中,如果材料的熱性能、邊界條件隨溫度變化,或是含有非線性單元和考慮熱輻射,則該分析為非線性分析3。在此次分析中選用了非線性單元solid98,因此分析類型屬于非線性穩(wěn)態(tài)輻射,則該分析為非線性分析3。在此次分析中選用了非線性單元solid98,因此分析類型屬于非線性穩(wěn)態(tài)熱分析。其熱平衡方程為:
式中:[c(T)]——比熱容矩陣,考慮系統(tǒng)內(nèi)能的增加
{T}——溫度對時間的導(dǎo)數(shù)
當熨平板結(jié)構(gòu)的溫度場求得后,就可以根據(jù)彈性體內(nèi)的已知變溫來決定體內(nèi)的溫度應(yīng)變。其計算公式為:
式中:a——彈性體的線膨脹系數(shù)
T——變溫,即溫度的改變量
ANSYS作為有限元分析軟件在熱分析問題方面具有強大的功能,而且界面友好,便于應(yīng)用。ANSYS熱分析基于能量守恒原理的熱平衡方程,用有限元法計算物體內(nèi)部各節(jié)點的溫度,并導(dǎo)出其他熱物理參數(shù)。
ANSYS進行熱分析的基本原理是將所處理的對象首先劃分成有限個單元(每個單元包含若干個節(jié)點),然后根據(jù)能量守恒原理求解一定邊界條件和初始條件下每一節(jié)點處的熱平衡方程,由此計算出各節(jié)點溫度值,繼而進一步求解出其他相關(guān)量。
2.2.1 幾何模型的建立
在對結(jié)構(gòu)進行熱分析時,為了提高有限元計算的精度,立體模型與實際零件結(jié)構(gòu)的形狀、尺寸應(yīng)一致,但立體模型過于復(fù)雜,必然增加有限元單元網(wǎng)格劃分的難度和計算工作量。因此對熨平板的形狀作了必要的簡化,例如,直徑小于50mm的孔在立體模型中一律省略,對熨平板某些部位的特殊結(jié)構(gòu),如尖角、圓角等在建模時可以忽略。對于用螺栓固定的零部件結(jié)構(gòu),假設(shè)螺栓已充分擰緊,把它們結(jié)合成一個實體零件。本文采用三維設(shè)計主流軟件Pro/E來建立熨平板實體模型,將建立好的模型通過Pro與ANSYS的專用接口導(dǎo)入ANSYS。
2.2.2 在ANSYS中進行前處理
首先需要定義材料屬性,包括彈性模量、泊松比、比熱容、導(dǎo)熱系數(shù)、線膨脹系數(shù)和密度。熨平板共有兩組不同材料屬性,一是底板材料通常為中碳鋼,二是箱體材料低碳鋼,如表1所示。
溫度與位移分布統(tǒng)計表 表2
攤鋪機熨平板材料屬性 表1
本文采用ANSYS中的直接耦合法進行熨平板熱分析,因此采用耦合單元solid98進行分網(wǎng)。鑒于熨平板結(jié)構(gòu)的復(fù)雜性,為減少網(wǎng)格劃分的困難,先對所有實體零件進行粘結(jié)操作,再使用智能網(wǎng)格劃分技術(shù)分網(wǎng)。
2.2.3 施加載荷,求解計算
考慮到模型左側(cè)面為對稱平面,因此在此接觸面上施加對稱邊界條件,根據(jù)熨平板在工作狀態(tài)下的實際安裝方式,模型與調(diào)拱裝置和大臂連接處均采用固定約束。模型的底面在工作時與瀝青混合料相接觸,因此在ANSYS中直接對其施加熱載荷,溫度為130℃(此值為瀝青溫度的平均值)。熨平板工作時與空氣發(fā)生自然對流,屬于熱分析中的第三類邊界條件,對流邊界條件可以作為面荷載施加于實體的表面來計算固體和流體間的熱交換。對流系數(shù)按照經(jīng)驗取6J,空氣溫度取20℃,忽略熱輻射,然后進入求解器求解計算。
2.2.4 進行通用后處理(POST1)并顯示分析結(jié)果
求解完成后,在通用后處理器(POsr1)中可以觀察溫度分布云圖、位移分布云圖與應(yīng)力分布云圖。溫度與位移分布如圖1和圖2所示。
①由圖1可見,溫度逐漸升高,表明熨平板底板處溫度較高(達到130℃),沿底板向上溫度迅速下降,噴火口向上溫度幾乎接近環(huán)境溫度。因此,安裝在噴火口附近的隔熱氈的厚度不必過大,以免浪費。
②由圖2可見,膨脹量逐漸加大,表明底板處與熱瀝青相接觸,膨脹量較大,上板處溫度接近環(huán)境溫度,膨脹量很小,因此熨平板下端明顯向上翹曲,顯然,熨平板熱變形后的橫向形狀為上底小、下底大的梯形,在底板溫度是130℃時,最大變形量是1.336mm,如表2所示。
在施工中采取預(yù)留變形量的方法,以此來補償熨平板熱變形對施工質(zhì)量的影響。
本文提出的方法在東海博世試車場進行運用,為檢測運用效果,現(xiàn)場檢測采取預(yù)留變形量熨平板在施工過程中對瀝青面層的影響,檢測方法是在熨平板加熱前,檢測熨平板的平整度,在攤鋪機勻速行駛5m后檢測橫斷面平整度,檢測結(jié)果如表3所示。
從表中數(shù)據(jù)顯示,檢測斷面的寬度是6m,其中平整度最差的是斷面1#,其最大值與最小值差值是1.62mm,精度完全滿足國內(nèi)最高攤鋪精度2mm/4m的平整度精度要求,取得的良好的效果。
[1]石偉.有限元分析基礎(chǔ)與應(yīng)用教程[M].北京:機械工業(yè)出版社,2010.
[2]米國華.淺談如何提高路面施工的平整度[J].中國高新技術(shù)企業(yè),2010(6).
[3]田國富,閆維佳,張大衛(wèi).應(yīng)用于攤鋪機熨平板的綜合設(shè)計方法[J].沈陽工業(yè)大學(xué)學(xué)報,2005(5).