穿甲燃烧弹侵彻陶瓷复合装甲和玻璃复合装甲的FEM-SPH耦合计算模型

为了提高小口径穿甲燃烧弹侵彻陶瓷复合装甲和玻璃复合装甲（透明装甲）的仿真分析精度,本文将传统的FEM（finite element method）-SPH（smooth particle hydrodynamics）耦合计算模型中穿甲燃烧弹弹芯的有限元模型和JC（Johnson-Cook）材料模型分别替换为SPH模型和JH2（Johnson-Holmquist-ceramics）材料模型,提出了新型FEM-SPH耦合计算模型。研究表明,新型FEM-SPH耦合计算模型可以有效模拟弹芯碎裂现象,减少SPH粒子和有限元耦合计算量,进而显著提高仿真模型的计算精度和计算效率,并给出了新型FEM-SPH耦合计算模型的有限元/粒子尺度和建模尺寸的优选结果。     

基于表面力模型的液滴冲击弹性基底SPH模拟

采用光滑粒子动力学SPH方法建立液滴冲击弹性基底的流固耦合数值模型,给出描述粘性流体和弹性固体运动的SPH离散方程和数值处理格式,引入人工耗散项来抑制标准SPH方法的数值震荡。为模拟液滴的表面张力效应,通过精确检测边界粒子,采用拉格朗日插值方法计算表面法向量和曲率,结合界面理论中的连续表面力CSF方法,建立了适用于自由表面液滴的表面力模型,方形液滴变形的模拟结果与拉普拉斯理论解吻合较好。随后,采用SPH流固耦合模型模拟1.0 mm直径水滴以不同速度(0.2 m/s～3.0 m/s)冲击两种薄板型基底,分析了基底弹性变形对液滴铺展、收缩以及回弹行为的影响。     

应用基于GPU的SPH方法模拟二维楔形体入水砰击问题

光滑粒子流体动力学(SPH)法是一种无网格的拉格朗日效值方法,广泛应用于计算流体领域模拟复杂自由表面流问题.SPH方法的主要缺点就是计算量过大,而基于GPU的并行计算方法可使SPH计算得到有效加速.本文应用基于GPU的SPH并行计算方法研究了二维楔形体的入水砰击问题.数值计算结果与文献中对应的解析解比较一致,验证了基于GPU的SPH方法的精度和可靠性.仿真结果同时显示基于GPU的并行计算方法可使SPH计算速度得到显著提高.     

抗差自适应插值粒子滤波及其在组合导航中的应用

针对粒子滤波存在的粒子退化和重要性密度函数难以选取的问题,在吸收抗差自适应滤波、二阶插值滤波和粒子滤波算法优点的基础上,提出了一种新的抗差自适应插值粒子滤波算法。该算法利用二阶插值滤波算法得到重要性密度函数,通过抗差自适应因子实时控制动力学模型误差及观测异常对导航解的影响。将该算法应用于SINS/CNS/SAR组合导航系统进行计算仿真,并与经典的粒子滤波算法进行比较分析。结果表明,提出的滤波算法得到的姿态误差控制在[-0.3′,+0.3′],速度误差控制在[-0.4 m/s,+0.4 m/s],位置误差控制在[-5 m,+5 m],性能明显优于经典的粒子滤波算法。新的滤波算法不但能够有效地抑制粒子退化,而且能够有效地控制动力学模型误差及观测异常的影响,提高了组合导航的滤波精度。     

基于SPH方法的表面张力修正算法及其应用

针对传统SPEI方法中基于CSF模型的表面张力算法,在计算边界、尖角等粒子缺失部位的曲率时存在偏差较大,且粒子秩序较差,对大变形问题表面张力计算精度较低的问题,在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到了表面张力修正方程组.应用本文方法模拟了水溶液中初始...     

高速碰撞数值计算中的SPH自适应粒子分布法

针对传统光滑粒子法在计算高速碰撞问题时会出现近邻粒子逸出核函数影响域而产生数值破坏这一缺陷,提出了一种根据粒子间距变化自动添加、合并粒子的SPH自适应粒子分布算法。采用该方法对Taylor碰撞和超高速碰撞问题进行了数值模拟,结果表明,该方法可以有效地消除计算中出现的数值破坏,提高计算精度。     

SPH方法在模拟线弹性波传播中的运用

通过对固体中波动问题的模拟建立了一种光滑粒子法的新形式,一种运用SPH的核函数的类似有限体积法的计算方法。通过对统计体积的修正以及对边界粒子的核函数修正,较好地解决了SPH方法中长期以来制约其被广泛应用的主要问题之一边界条件的表述。在此基础上成功地在光滑粒子法中实现了透射边界条件的模拟。同时利用反卷积修正使得较大粒子间距下的计算结果的精度大大提高。这种方法不但保持了SPH的简单性,而且很容易实现应力边界条件。     

完全拉格朗日SPH在冲击问题中的改进和应用

光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)在模拟固体大变形、破碎和裂纹扩展等问题中有天然的优势, 但SPH固有的拉伸不稳定缺陷是SPH在计算固体力学领域进一步应用的一大障碍. 完全拉格朗日SPH (total Lagrangian-SPH, TL-SPH)方法是一种有效的改善拉伸不稳定的措施, 但其仍面临边界区域精度低、界面条件难以施加、损伤裂纹难以模拟等缺陷. 因此, 首先将可达到二阶精度的高阶SPH方法与TL-SPH耦合, 为了节省高阶方法的计算量, 进一步简化粒子选取模式, 提出TL-SFPM (TL-simplified finite particle method)方法; 其次, 将可提高界面精度的DFPM (discontinuous finite particle method)方法与TL-SPH结合, 并提出一种基于黎曼解的界面接触算法, 通过在不同材料粒子间建立黎曼模型求解不同材料间的相互作用, 分别应用于流体?固体接触和固体?固体接触中; 再者, 为了捕捉固体受外载荷后的损伤程度及破坏模式, 提出一种完全拉格朗日框架下的粒子损伤破坏模型; 最后, 通过流?固冲击的带弹性挡板溃坝算例和固?固冲击的子弹撞击靶板算例验证提出的TL-SFPM方法、界面接触算法和损伤破坏模型的合理性和精确性, 进一步扩展TL-SPH方法在计算固体冲击问题中的应用.      

SPH-FEM接触算法在冲击动力学数值计算中的应用

为了充分发挥光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH)易于处理大变形以及有限元(Finite Element Method,FEM)计算精度和效率高的优势,论文基于无网格粒子接触算法,在有限元节点处设置背景粒子,通过接触力的方式计算SPH粒子和有限单元之间的接触问题...     

一种考虑界面不连续的改进的有限粒子法

有限粒子法（finite particle method，FPM）作为SPH（smoothed particle hydrodynamics）方法的重要改进，有效提高了边界区域粒子的近似精度，但是当FPM处理多物理场时，在不连续界面附近的计算精度会大大降低，并且FPM必须满足的矩阵非奇异性也提高了对界面处理的要求。本文中基于DSPH（discontinuous SPH）方法，提出了一种考虑界面不连续的改进FPM—DSFPM（discontinuous special FPM）法，旨在改善FPM在界面不连续处的计算精度，从而进一步提高其计算效率和稳定性。首先，分析了DSFPM的核近似精度。其次，根据不同的工程问题，给出DSFPM处理小变形和大变形问题的算法流程。利用DSFPM、DSPH和FPM等3种方法对弹性铝块小变形碰撞冲击算例进行了模拟，通过对比分析铝块的速度和应力以及计算时间验证了DSFPM算法在非连续界面处计算精度和计算效率的优势。最后，通过结合DSFPM和DFPM（discontinuous FPM）实现了对于大变形问题的模拟。     

Application of particle splitting method for both hydrostatic and hydrodynamic cases in SPH

Smoothed particle hydrodynamics(SPH) method with numerical diffusive terms shows satisfactory stability and accuracy in some violent fluid–solid interaction problems. However, in most simulations, uniform particle distributions are used and the multi-resolution, which can obviously improve the local accuracy and the overall computational efficiency, has seldom been applied. In this paper, a dynamic particle splitting method is applied and it allows for the simulation of both hydrostatic and hydrodynamic problems. The splitting algorithm is that, when a coarse(mother) particle enters the splitting region, it will be split into four daughter particles, which inherit the physical parameters of the mother particle. In the particle splitting process,conservations of mass, momentum and energy are ensured. Based on the error analysis, the splitting technique is designed to allow the optimal accuracy at the interface between the coarse and refined particles and this is particularly important in the simulation of hydrostatic cases. Finally, the scheme is validated by five basic cases, which demonstrate that the present SPH model with a particle splitting technique is of high accuracy and efficiency and is capable for the simulation of a wide range of hydrodynamic problems.Smoothed particle hydrodynamics(SPH)method with numerical diffusive terms shows satisfactory stability and accuracy in some violent fluid–solid interaction problems.However,in most simulations,uniform particle distributions are used and the multi-resolution,which can obviously improve the local accuracy and the overall computational efficiency,has seldom been applied.In this paper,a dynamic particle splitting method is applied and it allows for the simulation of both hydrostatic and hydrodynamic problems.The splitting algorithm is that,when a coarse(mother)particle enters the splitting region,it will be split into four daughter particles,which inherit the physical parameters of the mother particle.In the particle splitting process,conservations of mass,momentum and energy are ensured.Based on the error analysis,the splitting technique is designed to allow the optimal accuracy at the interface between the coarse and refined particles and this is particularly important in the simulation of hydrostatic cases.Finally,the scheme is validated by five basic cases,which demonstrate that the present SPH model with a particle splitting technique is of high accuracy and efficiency and is capable for the simulation of a wide range of hydrodynamic problems.     

一种新型SPH-FEM耦合算法及其在冲击动力学问题中的应用

为了充分发挥光滑粒子流体动力学方法(smoothed particle hydrodynamics,SPH)在处理大变形和有限元(finite element method,FEM)问题时计算精度高的优势,提出了一种新型SPH-FEM耦合算法.该耦合算法在大变形区域使用SPH粒子离散,其余区域使用有限元离散.在耦合界面...     

A robust shell element in meshfree SPH method

With the incorporation of total Lagrangian smoothed particle hydrodynamics(SPH) method equation and moving least square(MLS) function,the traditional SPH method is improved regarding the stability and consistency.Based on Mindlin-Ressiner plate theory,the SPH method simulating dynamic behavior via one layer of particles is applied to plate’s mid-plane,i.e.,a SPH shell model is constructed.Finally,through comparative analyses on the dynamic response of square,stiffened shells and cylindrical shells under various strong impact loads with common finite element software,the feasibility,validity and numerical accuracy of the SPH shell method are verified.Consequently,further researches on SPH shell may well pave the way towards solving problems involving dynamic plastic damage,tearing or even crushing.     

Curvilinear smoothed particle hydrodynamics

We suggest a new set of equations to employ smoothed particle hydrodynamics (SPH) in a curvilinear space, and we refer to it as curvSPH. In classical SPH, the horizontal and vertical resolution of discretization is supposed to be equal for fluid particles. However, curvSPH makes the horizontal and vertical resolutions independent from each other. This is performed by transformation of physical space into an appropriate computational space with a different scale in horizontal and vertical directions. Solving a problem using SPH in a curvilinear space also provides capability to model curved boundaries as straight lines. In classical SPH, special care is needed to reach a uniform mass distribution along curved boundaries; however, producing uniform mass distribution along a line using curvSPH is straight forward. Different simulations, including simulation of a flip bucket are performed to demonstrate the applicability of the proposed method. Good agreement of results with experimental data and classical SPH confirms the capabilities of curvSPH. Copyright © 2016 John Wiley & Sons, Ltd.     

多相流动的光滑粒子流体动力学方法研究综述

光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)具有粒子方法的无网格和全拉格朗日特征, 适用于具有界面大变形、不连续性和多物理场的多相流的高精度模拟. SPH方法模拟多相流已有大量报道, 具体的实现方式也大不相同. 本文首先阐述了采用SPH方法模拟流体的基本控制方程, 以及求解过程中需要考虑的流体压力求解、表面张力、固体边界等问题. 整理和总结了基于SPH方法进行多相流模拟的主要实现方式: (1)双流体模型的拉格朗日求解器: 两相离散为两组独立SPH粒子, 并用显式相间作用耦合两相; (2)多相SPH方法: SPH方法对多相流模拟的自然延伸, 相间作用由SPH参数隐式描述; (3) SPH与其他离散方法的耦合: 差异较大的两相各自采用不同离散方法, 发挥不同拉格朗日方法的优点; (4) SPH和基于网格方法的耦合: 网格方法处理简单的单相流动主体, 获得精度和效率间的平衡. 另外, 还在模拟参数物理化等方面论述了与SPH方法模拟多相流相关的一些改进和修正方法, 并在最后讨论和建议了提高多相流SPH模拟效率和精度的措施.      

一种SPH应力修正算法及自由表面流中的应用

提出了一种SPH应力修正算法,即模型中的拉应力和压应力分别采用不同的插值核函数和状态方程来处理,改善应力稳定性问题。介绍了一种改进的Quintic核函数,用于改善模型中压应力的稳定性。通过增加钟型核函数的光滑长度,改善模型中拉应力的稳定性。采用该应力修正算法模拟了无重力条件下方形液滴的震荡变形过程,对比分析了不同算法的模拟结果。此外,为进一步验证算法的适用性,模拟了溃坝算例。研究表明,改进的Quintic型核函数明显改善了粒子聚集现象,该SPH应力修正方法可以使液滴具有更均匀的粒子分布以及更光滑的自由表面,有效改善了SPH方法中的压应力不稳定作用以及自由表面流的模拟精度。     

An assessment of particle methods for approximating anisotropic dispersion

We derive a smoothed particle hydrodynamics (SPH) approximation for anisotropic dispersion that only depends upon the first derivative of the kernel function and study its numerical properties. In addition, we compare the performance of the newly derived SPH approximation versus an implementation of the particle strength exchange (PSE) method and a standard finite volume method for simulating multiple scenarios defined by different combinations of physical and numerical parameters. We show that, for regularly spaced particles, given an adequate selection of numerical parameters such as kernel function and smoothing length, the new SPH approximation is comparable with the PSE method in terms of convergence and accuracy and similar to the finite volume method. On other hand, the performance of both particle methods (SPH and PSE) decreases as the degree of disorder of the particle increases. However, we demonstrate that in these situations the accuracy and convergence properties of both particle methods can be improved by an adequate choice of some numerical parameters such as kernel core size and kernel function. Copyright © 2012 John Wiley & Sons, Ltd.     

Application of B‐splines and curved geometries to boundaries in SPH

Smoothed particle hydrodynamics (SPH) has been increasing in popularity rapidly in recent years and is being used for an ever wider range of applications. Central to almost all of these application is the inclusion of accurate wall boundaries. We present here a discussion of boundaries in SPH, in particular, focusing on reflected ghost‐particle boundaries. We show how one can include curved shapes as geometrical objects and more generally as parametric non‐uniform rational B‐spline (NURBS) curves. By properly considering the reflection operation, we derive a correction factor that demonstrably improves the accuracy of the SPH solution and present examples to confirm this. NURBS are standard for representing both 2D curves and 3D surfaces. We detail how they can be practically included in an SPH implementation, including how to calculate various required quantities and reflect particles in the NURBS object. Copyright © 2014 John Wiley & Sons, Ltd.