利用ANSYS10.0軟件進行輪軌彈性接觸有限元分析。文章分別對不同載荷條件下,小車輪軌接觸在不同的初始接觸位置處的應力進行分析,得出小車橫向偏移對接觸應力分布的影響。
隨著世界貿易量大幅增長,世界各個港口之間的貨物吞吐量逐年增加。世界各港口,特別是集裝箱港及大型散貨港在最近十幾年發(fā)展迅猛,隨著各個港口碼頭對裝卸效率要求的大幅提高,橋式起重機正趨于大型化、高速化發(fā)展。起重量越來越大。工作速度越來越高,不可避免的小車運行速度也加快。小車運行速度的加快也使得小車車輪發(fā)生和橫向移動,對小車運行的穩(wěn)定性帶來威脅。導向裝置間隙引起的部分偏斜、軌道側面或者水平輪的磨損引起的部分偏斜和軌道水平面上的直線性公差引起的部分偏斜都將造成車輪走偏。
有限元模型的建立與數(shù)值分析
在實際接觸中,由于車輪的橫向移動導致初始位置發(fā)生改變。我們考慮了四種典型的橫截面接觸位置,將有限元方法求得的計算值與赫茲接觸理論值做出比較。分析在不同的接觸位置處應力的分布情況。
在本次模型中采用800t/h的卸船機小車運行軌道進行分析,小車軌道與主軌道相同。
由于車輪的橫向運動,輪廓的每一處都可能發(fā)生接觸。對于車輪和軌道接觸的四個不同橫截面位置處,建立了有限元模型。為了得到滿意的接觸結果,接近接觸區(qū)的輪軌網(wǎng)格對四個模型中任何一個都是適用的。
將有限元計算結果與赫茲理論值進行比較,如下表所示:
表2-1 加載100t時有限元計算結果與赫茲理論值
圖形 | a | b | c | d |
計算值(MPa) | 681 | 719 | 991 | 1001 |
理論值(MPa) | 667 | 736 | 1003 | 1019 |
表2-2相同載荷下有限元計算值的變化率
載荷 | 100t | 變化率 (與a的比值) |
200t | 變化率 (與a的比值) |
300t | 變化率 (與a的比值) |
a | 681 | —— | 1363 | —— | 2043 | —— |
b | 719 | 5.58% | 1439 | 5.57% | 2158 | 5.62% |
c | 991 | 45.52% | 1980 | 45.28% | 2966 | 45.17% |
d | 1001 | 46.98% | 2003 | 46.95% | 3006 | 47.13% |
表2-1是車輪和軌道間的接觸應力和赫茲理論值的比較。從表中可以看到,有限元的計算值和赫茲理論值都很接近。
如果材料假定為線性的,軌道的表面曲率不發(fā)生改變,在輪軌接觸問題中應用赫茲接觸理論可以得到一個合理的準確結果。在小車發(fā)生橫向偏移的過程中,其軌道應力值不斷增大,當偏移到某一位置處,其應力值會有一個大的增量。在重載荷情況下,將會帶來塑性變形并將加劇其蹋面的磨損,影響輪軌的疲勞壽命。降低小車行駛性能與脫軌安全性。為保證小車高速行駛下的安全性,應控制小車在行駛過程中的橫向偏移,盡量減少發(fā)生應力突變。