基于光滑粒子流體動力學耦合有限元方法的水射流破煤參數(shù)數(shù)值模擬分析
發(fā)布時間:2020-02-02 00:23
【摘要】:利用非線性沖擊動力學軟件LS-DYNA,采用光滑粒子流體動力學(SPH)耦合有限元單元(FEM)方法,模擬了高壓水射流微觀破煤過程。根據(jù)模擬結(jié)果,紅嶺煤礦二1煤層水力擴孔最佳噴嘴直徑為2.0~2.5 mm,臨界破煤壓力為10 MPa,最佳破煤壓力為20~22 MPa。
【圖文】:
限單元的面M(如粒子i和j),由于這種侵入時非物理侵入,需要做必要的處理。有限元單元與光滑粒子界面模型如圖1,圖1中每個圓代表1個光滑粒子,根據(jù)動量守恒和能量守恒對位于滑移面上的相關(guān)節(jié)點進行速度調(diào)整。圖1有限元單元與光滑粒子界面模型m1v1+m2v2+mivi=m1v'1+m2v'1+miv'im1v1r1+m2v2r2+mivir3=m1v'ir1+m2v'2r2+miv'ir3v'i-v'1=(v'2-v'1)(ri-r1)/(r2-r1{)(1)式中:mi為光滑粒子i的質(zhì)量;m1和m2為與光滑粒子相接觸的有限元的2個節(jié)點M2和M3的質(zhì)量;v1、v2、vi分別為節(jié)點M2、M3及粒子i調(diào)整前沿線段M2、M3法線方向的速度分量;v'1、v'2、v'i為調(diào)整后的法向速度分量;r1、r2、ri為矢徑。根據(jù)式(1)可知M2、M3和i速度分量呈線性分布,通過插值方法對單元節(jié)點和光滑粒子位置進行調(diào)整,使粒子沿滑移面法向方向回到剛好和滑移面相切位置。1.2計算模型建立水射流采用SPH粒子建模,煤體采用FEM建模。選取1/4的水射流和煤體進行模擬(圖2)。水射流長度為50mm,半徑分別設(shè)定1.0~2.5mm,對稱面采用SPH_SYMMETRY_PLANE對稱約束。煤體模型的長寬高均為20mm,底面和外圍側(cè)面采用BOUNDARY_NONREFLECTING透射界面,對稱面采用CONSTRAINED_GLOBAL對稱約束。采用BOUNDARY_SPC_SET約束底面所有方向平動和轉(zhuǎn)動。水射流由11000個SPH粒子組成,煤體結(jié)構(gòu)由8000個正六面體組成,采用CONTACT_ERODING_NOTES_TO_SURFACE耦合接觸。1.3水的狀態(tài)方程和參數(shù)水射流選用Null材料,采用Gruneisen狀?
1和m2為與光滑粒子相接觸的有限元的2個節(jié)點M2和M3的質(zhì)量;v1、v2、vi分別為節(jié)點M2、M3及粒子i調(diào)整前沿線段M2、M3法線方向的速度分量;v'1、v'2、v'i為調(diào)整后的法向速度分量;r1、r2、ri為矢徑。根據(jù)式(1)可知M2、M3和i速度分量呈線性分布,通過插值方法對單元節(jié)點和光滑粒子位置進行調(diào)整,使粒子沿滑移面法向方向回到剛好和滑移面相切位置。1.2計算模型建立水射流采用SPH粒子建模,煤體采用FEM建模。選取1/4的水射流和煤體進行模擬(圖2)。水射流長度為50mm,半徑分別設(shè)定1.0~2.5mm,對稱面采用SPH_SYMMETRY_PLANE對稱約束。煤體模型的長寬高均為20mm,底面和外圍側(cè)面采用BOUNDARY_NONREFLECTING透射界面,對稱面采用CONSTRAINED_GLOBAL對稱約束。采用BOUNDARY_SPC_SET約束底面所有方向平動和轉(zhuǎn)動。水射流由11000個SPH粒子組成,,煤體結(jié)構(gòu)由8000個正六面體組成,采用CONTACT_ERODING_NOTES_TO_SURFACE耦合接觸。1.3水的狀態(tài)方程和參數(shù)水射流選用Null材料,采用Gruneisen狀態(tài)方程,水高速沖擊狀態(tài)下的壓力p的表達式為:圖2水射流破煤模型p=ρoC2μ[1+(1-γ02)μ-a2μ2][1-(S1-1)μ-S2μ2μ+1-S3μ3(μ+1)2]2+(γ0+aμ)E(2)式中:ρ0為水的初始密度;C為聲速;μ為體應(yīng)變,μ=(ρ/ρ0)-1;γ0為Gruneisen系數(shù);a為對γ0的一階體積校正量;S1、S2、S3為vs-vp曲線的斜率系數(shù);E為初始內(nèi)能。水的材料狀態(tài)參數(shù)見表1。表1水的狀態(tài)參數(shù)密度ρ0/(g·cm-3)聲速/(m·s-1)GruneisenC/系數(shù)體積校正
本文編號:2575540
【圖文】:
限單元的面M(如粒子i和j),由于這種侵入時非物理侵入,需要做必要的處理。有限元單元與光滑粒子界面模型如圖1,圖1中每個圓代表1個光滑粒子,根據(jù)動量守恒和能量守恒對位于滑移面上的相關(guān)節(jié)點進行速度調(diào)整。圖1有限元單元與光滑粒子界面模型m1v1+m2v2+mivi=m1v'1+m2v'1+miv'im1v1r1+m2v2r2+mivir3=m1v'ir1+m2v'2r2+miv'ir3v'i-v'1=(v'2-v'1)(ri-r1)/(r2-r1{)(1)式中:mi為光滑粒子i的質(zhì)量;m1和m2為與光滑粒子相接觸的有限元的2個節(jié)點M2和M3的質(zhì)量;v1、v2、vi分別為節(jié)點M2、M3及粒子i調(diào)整前沿線段M2、M3法線方向的速度分量;v'1、v'2、v'i為調(diào)整后的法向速度分量;r1、r2、ri為矢徑。根據(jù)式(1)可知M2、M3和i速度分量呈線性分布,通過插值方法對單元節(jié)點和光滑粒子位置進行調(diào)整,使粒子沿滑移面法向方向回到剛好和滑移面相切位置。1.2計算模型建立水射流采用SPH粒子建模,煤體采用FEM建模。選取1/4的水射流和煤體進行模擬(圖2)。水射流長度為50mm,半徑分別設(shè)定1.0~2.5mm,對稱面采用SPH_SYMMETRY_PLANE對稱約束。煤體模型的長寬高均為20mm,底面和外圍側(cè)面采用BOUNDARY_NONREFLECTING透射界面,對稱面采用CONSTRAINED_GLOBAL對稱約束。采用BOUNDARY_SPC_SET約束底面所有方向平動和轉(zhuǎn)動。水射流由11000個SPH粒子組成,煤體結(jié)構(gòu)由8000個正六面體組成,采用CONTACT_ERODING_NOTES_TO_SURFACE耦合接觸。1.3水的狀態(tài)方程和參數(shù)水射流選用Null材料,采用Gruneisen狀?
1和m2為與光滑粒子相接觸的有限元的2個節(jié)點M2和M3的質(zhì)量;v1、v2、vi分別為節(jié)點M2、M3及粒子i調(diào)整前沿線段M2、M3法線方向的速度分量;v'1、v'2、v'i為調(diào)整后的法向速度分量;r1、r2、ri為矢徑。根據(jù)式(1)可知M2、M3和i速度分量呈線性分布,通過插值方法對單元節(jié)點和光滑粒子位置進行調(diào)整,使粒子沿滑移面法向方向回到剛好和滑移面相切位置。1.2計算模型建立水射流采用SPH粒子建模,煤體采用FEM建模。選取1/4的水射流和煤體進行模擬(圖2)。水射流長度為50mm,半徑分別設(shè)定1.0~2.5mm,對稱面采用SPH_SYMMETRY_PLANE對稱約束。煤體模型的長寬高均為20mm,底面和外圍側(cè)面采用BOUNDARY_NONREFLECTING透射界面,對稱面采用CONSTRAINED_GLOBAL對稱約束。采用BOUNDARY_SPC_SET約束底面所有方向平動和轉(zhuǎn)動。水射流由11000個SPH粒子組成,,煤體結(jié)構(gòu)由8000個正六面體組成,采用CONTACT_ERODING_NOTES_TO_SURFACE耦合接觸。1.3水的狀態(tài)方程和參數(shù)水射流選用Null材料,采用Gruneisen狀態(tài)方程,水高速沖擊狀態(tài)下的壓力p的表達式為:圖2水射流破煤模型p=ρoC2μ[1+(1-γ02)μ-a2μ2][1-(S1-1)μ-S2μ2μ+1-S3μ3(μ+1)2]2+(γ0+aμ)E(2)式中:ρ0為水的初始密度;C為聲速;μ為體應(yīng)變,μ=(ρ/ρ0)-1;γ0為Gruneisen系數(shù);a為對γ0的一階體積校正量;S1、S2、S3為vs-vp曲線的斜率系數(shù);E為初始內(nèi)能。水的材料狀態(tài)參數(shù)見表1。表1水的狀態(tài)參數(shù)密度ρ0/(g·cm-3)聲速/(m·s-1)GruneisenC/系數(shù)體積校正
【相似文獻】
相關(guān)期刊論文 前3條
1 白玉柱,田取珍;沖擊破煤時煤巖中應(yīng)力場變化的實驗分析[J];山西煤炭;2003年02期
2 田取珍,岳龍旺;沖擊破煤參數(shù)的研究[J];太原理工大學學報;2001年02期
3 ;[J];;年期
相關(guān)博士學位論文 前1條
1 梁運培;高壓水射流鉆孔破煤機理研究[D];山東科技大學;2007年
相關(guān)碩士學位論文 前1條
1 李小召;薄煤層開采破煤機理與技術(shù)研究[D];河北工程大學;2013年
本文編號:2575540
本文鏈接:http://sikaile.net/kejilunwen/kuangye/2575540.html
最近更新
教材專著