光滑粒子动力学SPH方法应力不稳定性的一种改进方案

杨秀峰; 刘谋斌

doi:10.7498/aps.61.224701

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名

邮箱

手机号码

标题

留言内容

验证码

摘要

光滑粒子动力学方法是一种拉格朗日型无网格粒子方法, 在模拟大变形和自由表面流方面具有特殊的优势, 已经在工程和科学领域得到了广泛的应用. 然而, 长期以来, 传统光滑粒子动力学方法一直受到应力不稳定性的困扰, 从而限制了它的进一步发展和应用. 应力不稳定性的根本原因在于应力状态与核函数的不匹配:负压状态下粒子间产生吸引力, 吸引力随着粒子间距的减小而增大, 导致拉伸不稳定性;正压状态下粒子间产生排斥力, 排斥力随着粒子间距的减小而先增大后减小, 导致压缩不稳定性. 本文通过改进光滑粒子动力学方法的核函数和离散格式, 使得无论在正压还是负压状态下粒子间的作用力恒为排斥力, 且排斥力随着粒子间距的减小而增大, 从而防止粒子聚集等现象, 解决应力不稳定问题. 分别使用改进前后的光滑粒子动力学方法模拟两个典型的应力不稳定算例, 结果表明本文的改进方法能够有效地消除应力不稳定性.

关键词:

作者及机构信息

1.
中国科学院力学研究所, 水动力学与海洋工程重点实验室, 北京 100190

参考文献

[1]	Lucy L B 1977 Astron. J. 82 1013
[2]	Gingold R A, Monaghan J J 1977 Mon. Not. R. Astron. Soc. 181 375
[3]	Liu G R, Liu M B 2003 Smoothed particle hydrodynamics: a meshfree particle method (Singapore: World Scientific) p36
[4]	Liu M B, Liu G R, Zong Z 2008 Int. J. Comput. Methods 5 135
[5]	Monaghan J J 2005 Rep. Prog. Phys. 68 1703
[6]	Liu M B, Chang J Z 2010 Acta Phys. Sin. 59 7556 (in Chinese) [刘谋斌, 常建忠 2010 物理学报 59 7556]
[7]	Swegle J W, Hicks D L, Attaway S W 1995 J. Comput. Phys. 116 123
[8]	Dyka C T, Ingel R P 1995 Comput. Struct. 57 573
[9]	Zhang J Z, Zheng J, Yu K P, Wei Y J 2010 Eng. Mech. 27 65 (in Chinese) [张嘉钟, 郑俊, 于开平, 魏英杰 2010 工程力学 27 65]
[10]	Zheng J 2010 Ph. D. Dissertation (Harbin: Harbin Institute of Technoloty) (in Chinese) [郑俊 2010 博士学位论文 (哈尔滨: 哈尔滨工业大学)]
[11]	Morris J P 1995 Pulb. Astron. Soc. Aust. 13 97
[12]	Dyka C T, Randles P W, Ingel R P 1997 Int. J. Num. Methods Eng. 40 2325
[13]	Randles P W, Libersky L D 2000 Int. J. Num. Methods Eng. 48 1445
[14]	Morris J P 1996 Ph.D. Dissertation (Melbourne: Monash University)
[15]	Johnson G R, Stryk R A, Beissel S R 1996 Comput. Methods Appl. Mech. Eng. 139 347
[16]	Johnson G R, Beissel S R 1996 Int. J. Num. Methods Eng. 39 2725
[17]	Hicks D L, Liebrock L M 1999 Comput. Math. Appl. 38 1
[18]	Belytschko T, Guo Y, Liu W K, Xiao S P 2000 Int. J. Num. Methods Eng. 48 1359
[19]	Zhang G M, Wang X J, Hu X Z, Zhou Z 2003 Explosion and Shock Waves 23 219 (in Chinese) [张刚明, 王肖钧, 胡秀章, 周钟 2003 爆炸与冲击 23 219]
[20]	Hicks D L, Swegle J W, Attaway S W 1997 Appl. Math. Comput. 85 209
[21]	Wen Y, Hicks D L, Swegle J W 1994 Sandia Report SAND94-1932
[22]	Chen J K, Beraun J E, Jih C J 1999 Comput. Mech. 23 279
[23]	Monaghan J J 2000 J. Comput. Phys. 159 290
[24]	Gray J P, Monaghan J J, Swift R P 2001 Comput. Methods Appl. Mech. Eng. 190 6641
[25]	Jia B, Ma Z T, Zhang W, Pang B J 2010 J. of Harbin Institute of Technology 42 1369 (in Chinese) [贾斌, 马志涛, 张伟, 庞宝君 2010 哈尔滨工业大学学报 42 1369]
[26]	Liu W K, Jun S, Zhang Y F 1995 Int. J. Num. Methods in Fluids 20 1081
[27]	Dilts G A 1999 Int. J. Num. Methods Eng. 44 1115
[28]	Fu X J, Qiang H F, Yang Y C 2007 Advances in Mech. 37 375 (in Chinese) [傅学金, 强洪夫, 杨月诚 2007 力学进展 37 375]
[29]	Morris J P, Fox P J, Zhu Y 1997 J. Comput. Phys. 136 214
[30]	Fulk D A, Quinn D W 1996 J. Comput. Phys. 126 165

施引文献

[1]	Lucy L B 1977 Astron. J. 82 1013
[2]	Gingold R A, Monaghan J J 1977 Mon. Not. R. Astron. Soc. 181 375
[3]	Liu G R, Liu M B 2003 Smoothed particle hydrodynamics: a meshfree particle method (Singapore: World Scientific) p36
[4]	Liu M B, Liu G R, Zong Z 2008 Int. J. Comput. Methods 5 135
[5]	Monaghan J J 2005 Rep. Prog. Phys. 68 1703
[6]	Liu M B, Chang J Z 2010 Acta Phys. Sin. 59 7556 (in Chinese) [刘谋斌, 常建忠 2010 物理学报 59 7556]
[7]	Swegle J W, Hicks D L, Attaway S W 1995 J. Comput. Phys. 116 123
[8]	Dyka C T, Ingel R P 1995 Comput. Struct. 57 573
[9]	Zhang J Z, Zheng J, Yu K P, Wei Y J 2010 Eng. Mech. 27 65 (in Chinese) [张嘉钟, 郑俊, 于开平, 魏英杰 2010 工程力学 27 65]
[10]	Zheng J 2010 Ph. D. Dissertation (Harbin: Harbin Institute of Technoloty) (in Chinese) [郑俊 2010 博士学位论文 (哈尔滨: 哈尔滨工业大学)]
[11]	Morris J P 1995 Pulb. Astron. Soc. Aust. 13 97
[12]	Dyka C T, Randles P W, Ingel R P 1997 Int. J. Num. Methods Eng. 40 2325
[13]	Randles P W, Libersky L D 2000 Int. J. Num. Methods Eng. 48 1445
[14]	Morris J P 1996 Ph.D. Dissertation (Melbourne: Monash University)
[15]	Johnson G R, Stryk R A, Beissel S R 1996 Comput. Methods Appl. Mech. Eng. 139 347
[16]	Johnson G R, Beissel S R 1996 Int. J. Num. Methods Eng. 39 2725
[17]	Hicks D L, Liebrock L M 1999 Comput. Math. Appl. 38 1
[18]	Belytschko T, Guo Y, Liu W K, Xiao S P 2000 Int. J. Num. Methods Eng. 48 1359
[19]	Zhang G M, Wang X J, Hu X Z, Zhou Z 2003 Explosion and Shock Waves 23 219 (in Chinese) [张刚明, 王肖钧, 胡秀章, 周钟 2003 爆炸与冲击 23 219]
[20]	Hicks D L, Swegle J W, Attaway S W 1997 Appl. Math. Comput. 85 209
[21]	Wen Y, Hicks D L, Swegle J W 1994 Sandia Report SAND94-1932
[22]	Chen J K, Beraun J E, Jih C J 1999 Comput. Mech. 23 279
[23]	Monaghan J J 2000 J. Comput. Phys. 159 290
[24]	Gray J P, Monaghan J J, Swift R P 2001 Comput. Methods Appl. Mech. Eng. 190 6641
[25]	Jia B, Ma Z T, Zhang W, Pang B J 2010 J. of Harbin Institute of Technology 42 1369 (in Chinese) [贾斌, 马志涛, 张伟, 庞宝君 2010 哈尔滨工业大学学报 42 1369]
[26]	Liu W K, Jun S, Zhang Y F 1995 Int. J. Num. Methods in Fluids 20 1081
[27]	Dilts G A 1999 Int. J. Num. Methods Eng. 44 1115
[28]	Fu X J, Qiang H F, Yang Y C 2007 Advances in Mech. 37 375 (in Chinese) [傅学金, 强洪夫, 杨月诚 2007 力学进展 37 375]
[29]	Morris J P, Fox P J, Zhu Y 1997 J. Comput. Phys. 136 214
[30]	Fulk D A, Quinn D W 1996 J. Comput. Phys. 126 165

引用本文:

Citation:

计量

文章访问数: 5597
PDF下载量: 1793
被引次数: 0

光滑粒子动力学SPH方法应力不稳定性的一种改进方案

1. 中国科学院力学研究所, 水动力学与海洋工程重点实验室, 北京 100190

收稿日期: 2012-03-16

修回日期: 2012-06-12

刊出日期: 2012-11-05

摘要: 光滑粒子动力学方法是一种拉格朗日型无网格粒子方法, 在模拟大变形和自由表面流方面具有特殊的优势, 已经在工程和科学领域得到了广泛的应用. 然而, 长期以来, 传统光滑粒子动力学方法一直受到应力不稳定性的困扰, 从而限制了它的进一步发展和应用. 应力不稳定性的根本原因在于应力状态与核函数的不匹配:负压状态下粒子间产生吸引力, 吸引力随着粒子间距的减小而增大, 导致拉伸不稳定性;正压状态下粒子间产生排斥力, 排斥力随着粒子间距的减小而先增大后减小, 导致压缩不稳定性. 本文通过改进光滑粒子动力学方法的核函数和离散格式, 使得无论在正压还是负压状态下粒子间的作用力恒为排斥力, 且排斥力随着粒子间距的减小而增大, 从而防止粒子聚集等现象, 解决应力不稳定问题. 分别使用改进前后的光滑粒子动力学方法模拟两个典型的应力不稳定算例, 结果表明本文的改进方法能够有效地消除应力不稳定性.

关键词:

全文HTML

参考文献 (30)

搜索

留言板