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

被引:32
作者
杨秀峰
刘谋斌
机构
[1] 中国科学院力学研究所,水动力学与海洋工程重点实验室
关键词
光滑粒子动力学; 稳定性; 核函数; 离散格式;
D O I
暂无
中图分类号
O35 [流体力学];
学科分类号
080103 ; 080704 ;
摘要
光滑粒子动力学方法是一种拉格朗日型无网格粒子方法,在模拟大变形和自由表面流方面具有特殊的优势,已经在工程和科学领域得到了广泛的应用.然而,长期以来,传统光滑粒子动力学方法一直受到应力不稳定性的困扰,从而限制了它的进一步发展和应用.应力不稳定性的根本原因在于应力状态与核函数的不匹配:负压状态下粒子间产生吸引力,吸引力随着粒子间距的减小而增大,导致拉伸不稳定性;正压状态下粒子间产生排斥力,排斥力随着粒子间距的减小而先增大后减小,导致压缩不稳定性.本文通过改进光滑粒子动力学方法的核函数和离散格式,使得无论在正压还是负压状态下粒子间的作用力恒为排斥力,且排斥力随着粒子间距的减小而增大,从而防止粒子聚集等现象,解决应力不稳定问题.分别使用改进前后的光滑粒子动力学方法模拟两个典型的应力不稳定算例,结果表明本文的改进方法能够有效地消除应力不稳定性.
引用
收藏
页码:261 / 268
页数:8
相关论文
共 6 条
  • [1] 耗散粒子动力学处理复杂固体壁面的一种有效方法
    刘谋斌
    常建忠
    [J]. 物理学报, 2010, 59 (11) : 7556 - 7563
  • [2] 光滑粒子流体动力学法拉伸不稳定性改进研究
    贾斌
    马志涛
    张伟
    庞宝君
    [J]. 哈尔滨工业大学学报, 2010, 42 (09) : 1369 - 1373
  • [3] 流体力学中SPH算法张力不稳定性研究
    张嘉钟
    郑俊
    于开平
    魏英杰
    [J]. 工程力学, 2010, 27 (02) : 65 - 72
  • [4] 固体介质中SPH方法的拉伸不稳定性问题研究进展
    傅学金
    强洪夫
    杨月诚
    [J]. 力学进展, 2007, (03) : 375 - 388
  • [5] 光滑粒子法中的一种新的核函数
    张刚明
    王肖钧
    胡秀章
    周钟
    [J]. 爆炸与冲击, 2003, (03) : 219 - 224
  • [6] An improvement for tensile instability in smoothed particle hydrodynamics
    Chen, JK
    Beraun, JE
    Jih, CJ
    [J]. COMPUTATIONAL MECHANICS, 1999, 23 (04) : 279 - 287