期刊
  出版年
  关键词
结果中检索 Open Search
Please wait a minute...
选择: 显示/隐藏图片
改进的物理粘性SPH方法及其在溃坝问题中的应用
韩亚伟, 强洪夫
计算物理    2012, 29 (5): 693-699.  
摘要307)      PDF (1285KB)(1281)   
在低雷诺数物理粘性SPH方法基础上引入再生核粒子法进行密度重构,既避免了用人工粘性所导致的数值耗散问题,又提高了低雷诺数物理粘性SPH方法的数值稳定性;以溃坝问题为例,对比分析低雷诺数物理粘性SPH方法和本文方法的仿真结果表明,本文方法可有效消除数值不稳定,压强和速度分布更加光滑,粒子秩序更好,可应用于雷诺数较高或粘性不可忽略的流动问题.
相关文章 | 多维度评价
基于SPH-FEM转换算法的7.62mm步枪弹冲击30CrMnSiA钢板的数值计算
张志春, 强洪夫, 傅学金, 夏薇
计算物理    2012, 29 (1): 73-81.  
摘要559)      PDF (7972KB)(840)   
从转换判据和时间步长控制两个方面对光滑粒子流体动力学-有限元(SPH-FEM)转换算法进行改进,采用改进的SPH-FEM转换算法对7.62 mm步枪弹冲击特殊热处理的30CrMnSiA钢板进行全尺寸三维数值计算,将子弹在冲击过程中发生严重畸变的有限单元转换为SPH粒子继续参与计算.子弹采用弹塑性模型,靶板采用修正的Johnson-Cook强度模型和Gruneisen状态方程.针对不同的子弹初速,计算靶板盘形凹陷和冲塞两种破坏模式,计算结果与实验结果吻合较好.表明改进的SPH-FEM转换算法能够合理高效地利用有限单元和SPH粒子,为低强度子弹冲击高强度靶板的计算提供一种有效途径.
相关文章 | 多维度评价
修正表面张力算法的SPH方法及其实现
强洪夫, 陈福振, 高巍然
计算物理    2011, 28 (3): 375-384.  
摘要333)      PDF (745KB)(1423)   
在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到表面张力修正方程组;通过半圆形算例测试方法和Morris方法在边界定位、法向计算和曲率计算等影响表面张力关键因素的求解精度,研究曲率计算中应采用的光滑长度值.模拟初始方形液滴在表面张力作用下的自然变化过程,并与Morris方法及VOF有限体积法进行对比,表明方法精度较高,稳定性较好.最后,模拟水溶液中两个油滴的互溶过程.
相关文章 | 多维度评价
基于粘结界面模型的三维裂纹扩展研究
张志春, 强洪夫, 周伟
计算物理    2010, 27 (4): 586-592.  
摘要260)      PDF (342KB)(1154)   
根据界面断裂的特点,将粘结界面模型应用于有限元,采用非线性显式动力学算法,编写粘结-体积单元计算程序.其中体积和粘结单元分别采用四面体和三棱柱单元;粘结单元本构关系选择双线性内聚力-位移曲线,避免了界面完全断裂后发生穿透;以界面相对位移的形式定义线性损伤准则;以Benzeggagh-Kenane模式的扩展准则确定界面失效位移;以Turon模式的初始损伤准则确定界面启裂位移.该程序能够实现混合加载模式下材料界面断裂问题的三维数值模拟.用所编写的程序(CVFEM)分别对Ⅰ,Ⅱ,Ⅲ型裂纹扩展问题进行数值模拟,并且与Abaqus6.7计算结果进行对比.
相关文章 | 多维度评价
完全变光滑长度SPH法及其实现
强洪夫, 高巍然
计算物理    2008, 25 (5): 569-575.  
摘要399)      PDF (506KB)(1500)   
提出完全变光滑长度SPH法及其算法实现.方程组基于对称形式核函数近似,SPH密度演化方程与变光滑长度方程隐式关联;在Springel提出的全守恒SPH方程组基础上,通过将分散核近似形式改进为对称核近似形式得到SPH动量方程和能量方程.采用迭代求解密度演化方程和变光滑长度方程,显式求解SPH动量方程和能量方程,增加的计算量相对很少.给出三个1D激波管算例和2D Sedov算例验证方法的有效性.数值结果表明,算法保持动量和能量的守恒律,解决了传统SPH法中由于变光滑长度影响带来的计算误差,且在模拟2DSedov问题时能得到比Springel方法更准确的压强峰值位置和中心压强值.特别适合于模拟爆炸与冲击、大变形大扭曲等密度梯度和光滑长度梯度剧烈变化的问题.
相关文章 | 多维度评价