流冰碰撞下閘門的有限元分析

摘 要：在黃河開河時期，經常發生冰塞、冰壩等災害，對水工建筑物造成破壞。文章運用ANSYS/LS-DYNA軟件對流冰撞擊閘門的過程進行了仿真模擬，模擬了流冰撞擊閘門時冰和閘門的受力過程，同時繪制出不同速度下流冰的壓強曲線。文章從有限元仿真模擬研究上，對流冰碰撞水工建筑物問題進行了分析，為以后更深入的研究工作提供了參考。

關鍵詞：ANSYS/LS-DYNA；仿真模擬；流冰；閘門

黃河內蒙古河段地處黃河流域最北端，由于緯度較高，因此造成氣溫上暖下寒，并且冬季嚴寒而漫長，氣溫在0℃以下的時間可持續4～5個月，是穩定封河河段。在開河階段，極其容易發生冰壩、冰塞等險情，并且大塊的流冰對水工建筑物會造成不同程度的破壞，因此，開展流冰對建筑物的破壞研究意義重大。目前利用有限元軟件ANSYS/LS-DYNA對撞擊問題進行模擬分析越來越普遍，文章利用此軟件對流冰撞擊閘門問題進行分析，希望可以為冰凌的防治工作提供一定參考。

1 流冰碰撞閘門的運動方程

在文章中，流冰碰撞閘門屬于低速度的非線性問題，在碰撞的瞬間流冰將發生變形，整個碰撞的過程中伴隨著能量交換和吸收，在有限元方法中，其運動方程可以表示為

2 數值模型的建立

文章采用顯示六面體實體單元：SOLID164，建立三維實體模型，并采用沙漏粘性阻尼，以節省存儲量，增快計算速度，克服零能模式。

文章建立的有限元模型分為平板剛閘門和流冰兩部分，其中平板剛閘門部分由五根縱梁，兩根主梁，兩根次梁和面板組成，面板寬為2640mm，長為3000mm，厚為20mm，主梁長、寬、高分別為2640mm、80mm、320mm。冰的尺寸為2m×2m×0.2m，水深5m，根據冰撞擊閘門速度的不同，分別建立了流冰速度為0.5m/s、0.7m/s、1.0m/s、1.5m/s四種碰撞類型。

碰撞模型如圖1所示。

閘門跟流冰的材料參數如表1、表2所示。

3 流冰與閘門的碰撞分析

3.1 碰撞能量分析

流冰撞擊閘門是一個瞬間的過程，其能量交換也是在瞬間完成的，在動能轉化為內能的過程中，整體的總能量是守恒的。具體的能量轉換為：開始只有流冰具有動能，碰撞后，轉化為冰塊變形破碎的能量、剩余的動能、各物體之間因摩擦產生的內能、閘門的動能和內能和沙漏能。

碰撞過程中能量變化如圖2所示，曲線A表示動能，曲線B表示內能，曲線C表示總能量。由圖可以看出，當流冰與閘門碰撞瞬間，流冰動能突然減小，在0.002s時降到了最低，隨后在一個范圍內有限的波動。在整個撞擊過程中，開始內能為零，動能最大，后來內能大于動能，說明整個過程動能大部分轉化成了內能。

3.2 碰撞時流冰應力云圖

流冰在碰撞瞬間的第一主應力云圖和第三主應力云圖如圖3、4所示。

從圖中可以看出，在碰撞瞬間冰最大拉應力為1.300×10-5Mpa，最大壓應力為2.858×106pa，在碰撞后，冰的拉應力和壓應力都在減小。

3.3 碰撞時閘門的應力云圖

在碰撞瞬間閘門的第一主應力云圖和第三主應力云圖5、圖6所示。

由圖可以看出，在碰撞瞬間，閘門的最大拉應力為2.700×10-4Mpa，最大壓應力為3.457×10-5Mpa。

分析流冰和閘門接觸部分的單元，找到應力最大的單元，如圖7所示，單元所受到的應力在碰撞的瞬間急劇增加，0.07s左右達到最大壓應力值，約1.49×107Pa，小于鋼材的屈服應力2.9×107Pa，所以閘門在碰撞中沒有發生破壞。

4 不同流冰速度對冰壓強的影響

在流冰尺寸為2m×2m×0.2m，流冰速度分別為0.5m/s、0.7m/s、1.0m/s、1.5m/s時繪制出冰體壓強趨勢如圖8所示，從圖中可以看出，隨著流冰速度的增加，冰體壓強呈增大趨勢。

5 結束語

作者還做了當流冰的長寬都為2m×2m時，冰的厚度分別為0.1m、0.3m時候閘門的分析，分析得出當流冰速度分別為0.5m/s、0.7m/s、1.0m/s、1.5m/s時候，閘門所受的最大壓強都小于鋼材的屈服應力，閘門不會受到破壞。但是當流冰的長寬繼續同步增大時，閘門所受的壓強也隨著增大，因此應提前對大體積流冰進行爆破，以防止對水工建筑物造成破壞。

參考文獻

[1]時黨勇，李裕春，張勝民.基于ANSYS/LS-DYNA8.1進行顯示動力分析[M].北京：清華大學出版社，2005.

[2]孟聞遠，張蕊，王俊峰.黃河防凌減災方案新探索[J].華北水利水電學院學報， 2012，04.

[3]可素娟，王敏，饒素秋，等.黃河冰凌研究[M].鄭州：黃河水利出版社，2002，12.

[4]熊安平.基于LS-DYNA仿真模擬對船橋撞擊的研究[D].華東交通大學，2012.

作者簡介：楊再云，男，漢族，云南麗江，工程碩士，華北水利水電大學，研究方向：水工結構。