高巍然 作品数:14 被引量:135 H指数:8 供职机构: 第二炮兵装备研究院 更多>> 发文基金: 国家重点基础研究发展计划 教育部“新世纪优秀人才支持计划” 教育部跨世纪优秀人才培养计划 更多>> 相关领域: 理学 化学工程 兵器科学与技术 金属学及工艺 更多>>
基于罚函数SPH新方法的水模拟充型过程的数值分析 被引量:14 2011年 同传统的网格法相比,光滑粒子流体动力学(SPH)方法非常适合于求解大变形和自由表面流动问题。阐述了SPH理论及其应用,并用罚函数处理流体与壁面的相互作用,以解决传统SPH本质边界条件不易施加的问题。对水模拟的充型过程实验进行数值分析,并和文献实验结果以及传统SPH进行对比,最终表明仿真结果与实验非常吻合,比传统SPH方法捕捉到了更多的流动细节和特征;同时由于使用了XSPH算法,使得粒子更加整齐有序。 强洪夫 韩亚伟 王坤鹏 高巍然关键词:罚函数法 本质边界条件 光滑粒子流体动力学-有限元法接触算法研究 被引量:3 2011年 耦合光滑粒子流体动力学(SPH)方法和有限元法(FEM),是研究冲击动力学问题的一种有效途径。为解决SPH粒子和有限单元间的接触问题,提出了SPH-FEM接触算法。该算法是在有限元节点处设置背景粒子,采用SPH接触算法的思想,计算施加在SPH粒子和有限元节点上的接触力,并且以外力的形式分别加入到SPH动量方程和有限元动力学方程中。利用SPH-FEM接触算法,对两杆撞击以及圆柱形钢弹正冲击钢板发生的冲塞破坏过程进行了三维数值模拟,靶板采用含损伤的Johnson-Cook模型和Grüneisen状态方程,模拟结果与实验结果吻合较好。 张志春 强洪夫 高巍然关键词:有限元法 修正表面张力算法的SPH方法及其实现 被引量:14 2011年 在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到表面张力修正方程组;通过半圆形算例测试方法和Morris方法在边界定位、法向计算和曲率计算等影响表面张力关键因素的求解精度,研究曲率计算中应采用的光滑长度值.模拟初始方形液滴在表面张力作用下的自然变化过程,并与Morris方法及VOF有限体积法进行对比,表明方法精度较高,稳定性较好.最后,模拟水溶液中两个油滴的互溶过程. 强洪夫 陈福振 高巍然关键词:光滑粒子流体动力学 曲率 聚能射流侵彻有限尺寸混凝土靶板的SPH数值模拟 聚能射流侵彻混凝土靶板不仅包括聚能炸药爆炸挤压药型罩形成射流的过程还包括混凝土在高速冲击载荷下的破坏损伤过程,这些过程都涉及到大变形问题.SPH 方法作为一种无网格纯Lagrangian 粒子方法非常适合求解这类问题.本... 高巍然 张义忠 陈福振关键词:光滑粒子流体动力学 聚能射流 侵彻 数值模拟 文献传递 SPH-FEM接触算法在冲击动力学数值计算中的应用 被引量:27 2011年 为了充分发挥光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH)易于处理大变形以及有限元(Finite Element Method,FEM)计算精度和效率高的优势,论文基于无网格粒子接触算法,在有限元节点处设置背景粒子,通过接触力的方式计算SPH粒子和有限单元之间的接触问题.使用SPH-FEM接触算法分别对球头钢弹斜冲击钢板和平头钢弹正冲击钢板的情况进行了三维数值计算,SPH采用完全变光滑长度算法,FEM采用矩阵向量积的EBE(Element-By-Element)算法.将SPH-FEM接触算法的计算结果与LS-DYNA的计算结果以及实验结果进行了对比验证. 张志春 强洪夫 高巍然关键词:有限元方法 基于SPH方法的表面张力修正算法及其应用 被引量:9 2011年 针对传统SPH方法中基于CSF模型的表面张力算法,在计算边界、尖角等粒子缺失部位的曲率时存在偏差较大,且粒子秩序较差,对大变形问题表面张力计算精度较低的问题,在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到了表面张力修正方程组。应用本文方法模拟了水溶液中初始方形油液滴在表面张力作用下的自然变化过程,与VOF有限体积法进行对比,结果非常吻合;最后,基于本文方法模拟了二维水溶液中两个油滴的碰撞破碎过程。 强洪夫 陈福振 高巍然关键词:光滑粒子流体动力学 曲率 基于SPH方法的低韦伯数下三维液滴碰撞聚合与反弹数值模拟研究 被引量:4 2012年 光滑粒子流体动力学方法(SPH方法)作为纯拉格朗日粒子方法,可以有效避免网格法在模拟大变形过程中带来的网格扭曲等缺陷,适合模拟含大变形的液滴碰撞聚合与反弹过程。该文基于Ott和Schnetter提出的修正SPH方法,利用有限差分与SPH一阶导数相结合的方法处理粘性项中的二阶导数问题,进行Couette流算例验证,数值解与理论解基本吻合;利用范德华模型中的分子内聚力项处理表面张力问题,进行初始方形液滴在表面张力作用下自然变化过程的算例验证,变化过程符合宏观实际;运用阐述的SPH方法模拟三维低韦伯数条件下相同尺寸的两个正十四烷液滴正碰与偏心碰撞动力学问题,同时对相同条件下的油与水等难溶液滴间正碰与偏心碰撞问题进行了模拟研究,仿真结果均与文献实验结果进行对比,两者非常吻合,表明该方法适合处理诸如液滴碰撞等大变形的问题。 强洪夫 陈福振 高巍然关键词:粘性 完全变光滑长度SPH法及其实现 被引量:22 2008年 提出完全变光滑长度SPH法及其算法实现.方程组基于对称形式核函数近似,SPH密度演化方程与变光滑长度方程隐式关联;在Springel提出的全守恒SPH方程组基础上,通过将分散核近似形式改进为对称核近似形式得到SPH动量方程和能量方程.采用迭代求解密度演化方程和变光滑长度方程,显式求解SPH动量方程和能量方程,增加的计算量相对很少.给出三个1D激波管算例和2D Sedov算例验证方法的有效性.数值结果表明,算法保持动量和能量的守恒律,解决了传统SPH法中由于变光滑长度影响带来的计算误差,且在模拟2DSedov问题时能得到比Springel方法更准确的压强峰值位置和中心压强值.特别适合于模拟爆炸与冲击、大变形大扭曲等密度梯度和光滑长度梯度剧烈变化的问题. 强洪夫 高巍然关键词:光滑粒子流体动力学 无网格法 激波 修正变光滑长度SPH方法及其应用 被引量:4 2007年 提出了修正的变光滑长度SPH方法及其算法实现。与传统变光滑长度SPH法相比,新方程组基于对称形式核函数近似,并从根本上考虑了变光滑长度影响。方程组中,密度演化方程与变光滑长度方程隐式关联;而SPH动量方程和能量方程在Springel提出的全守恒SPH方程组基础上,改进其分散核近似形式为对称核近似形式得到。由于新方程组密度演化方程和变光滑长度方程隐式耦合,采用迭代求解密度演化方程和变光滑长度方程,显式求解SPH动量方程和能量方程的办法,迭代过程增加的计算量相对很少。给出了2个一维激波管算例和二维Sedov算例来验证方法的有效性。结果表明,新方法很好地解决了传统SPH法中变光滑长度影响,特别是在模拟二维Sedov问题时新方法能得到比Springel方法更准确的压强峰值位置,中心压强也更接近理论值。特别适合模拟爆炸与冲击、大变形大扭曲等密度和光滑长度变化剧烈的问题。 强洪夫 高巍然关键词:光滑粒子流体动力学 无网格法 激波管 基于修正SPH方法的爆轰波绕射传播的数值模拟 被引量:2 2009年 为研究内置金属球体的球壳装药在一点起爆情况下的非理想爆轰波的传播行为,基于修正的光滑粒子流体动力学(SPH)方法进行绕射爆轰波传播的数值模拟,并应用一种工程上简便的近似计算方法对该球壳装药模型的爆轰波传播时间进行计算。结果表明,数值模拟的特征爆轰时间与理论计算结果非常吻合,验证了修正SPH方法在爆炸数值模拟领域的可行性。 周建辉 孙新利 高巍然 韩亚伟 聂虎关键词:爆轰波 数值模拟