發(fā)布時間：2022-02-15 16:23

　　最優(yōu)輸運無網(wǎng)格方法（optimal transportation meshfree method,OTM）相比傳統(tǒng)的SPH方法,可以有效克服現(xiàn)有顯式分析中能量不守恒問題,完全避免SPH方法中存在的拉力不穩(wěn)定性,并可規(guī)避零能模式的出現(xiàn).本文采用OTM方法對鋁制球形彈丸超高速撞擊鋁合金靶板問題開展了數(shù)值模擬研究.為描述鋁在高溫高壓和高應(yīng)變率條件下的動態(tài)力學(xué)響應(yīng),采用Johnson-Cook材料模型和Mie-Grüneisen狀態(tài)方程對超高速撞擊單層靶板進行數(shù)值模擬,得到的彈坑直徑、碎片云長度與寬度以及內(nèi)核碎片云的形態(tài)和分布與實驗吻合較好.此外,對超高速撞擊雙層靶板問題也進行了數(shù)值模擬,模擬得到的碎片云形貌參數(shù)符合實驗結(jié)果,二次碎片云的模擬長度與實驗結(jié)果僅相差2.6%.數(shù)值模擬結(jié)果表明OTM方法非常適合于模擬超高速碰撞問題. 

【文章來源】：北京理工大學(xué)學(xué)報. 2020,40(06)北大核心EICSCD

【文章頁數(shù)】：6 頁

【部分圖文】：

物質(zhì)點與節(jié)點空間離散方案

基于OTM方法所編制程序的主要算法流程如圖2所示.① 幾何模型的初始化. 一般地,將d維的問題域(幾何模型)Ω0進行網(wǎng)格劃分生成初始節(jié)點集{xa,0,a=1,2,…,N},并取每個單元的形心作為一個物質(zhì)點的初始位置,組成物質(zhì)點集{xp,0,p=1,2,…,M},物質(zhì)點的初始體積可設(shè)為該物質(zhì)點所在單元占據(jù)的空間大小.

如圖3所示,鋁制球形彈丸半徑為5 mm,鋁靶板的尺寸為50 mm×50 mm×4 mm,鋁球的撞擊速度設(shè)為與實驗[8] 一致的6.18 km/s. 計算過程中,時間步長由CFL條件給出,每個時間步長的量級約為0.1 μs,室溫設(shè)為300 K.圖4為20 μs時OTM方法計算得到的碎片分布,將其與Hiermaier的實驗[8] 及Hiermaier數(shù)值模擬結(jié)果[8] 進行對比,可以發(fā)現(xiàn)OTM得到的碎片云,整個碎片云的大部分質(zhì)量中位于外泡碎片云的前端,氣泡的大致形貌較為圓潤,在形貌上和分布上都與實驗結(jié)果吻合較好.

【參考文獻】：
期刊論文
[1]超高速碰撞成坑特性分子動力學(xué)模擬[J]. 巨圓圓,張慶明,龍仁榮,武強,龔自正.  北京理工大學(xué)學(xué)報. 2018(12)



本文編號：3626946