水底管道抛石加固过程的DEM-SPH耦合分析

水底管道的抛石加固过程是典型的颗粒-流体耦合问题.采用DEM-SPH耦合方法模拟颗粒-流体系统,其中离散元方法(DEM)用于模拟落石,光滑粒子流体动力学方法(SPH)用于模拟流体.通过三维Voronoi切割算法生成不规则形状的多面体,并基于闵可夫斯基原理构造扩展多面体形态的落石单元.通过SPH的边界排斥力模型计算颗粒与流体间的作用力,从而建立DEM-SPH耦合方法.采用该方法模拟溃坝与楔形块入水的过程,将计算结果与试验结果及其他数值结果进行对比分析,分别验证了 SPH与DEM-SPH耦合方法的合理性.建立锥形结构模拟卸料斗和水底管道,采用DEM-SPH耦合方法模拟落石通过卸料斗入水并与水底管道相互作用的过程,确定了落石和水对管道的作用力,并分析了卸料斗静止与运动时的落石堆积情况.以上研究表明,DEM-SPH方法可有效模拟颗粒材料、水和工程结构的相互作用,可进一步应用于水下抛石过程的结构和参数设计.     

三维自适应FE-SPH耦合算法在多层间隔金属靶侵彻问题中的应用

鉴于有限元算法不能有效地模拟侵彻过程所产生的金属碎片, 本文中基于三维自适应FE-SPH耦合算法的基本理论, 自主开发了模拟多层间隔金属靶侵彻问题的三维FE-SPH耦合计算程序。该程序采用四面体单元对多层间隔金属靶侵彻模型进行初始离散, 计算过程中, 当四面体单元等效塑性应变超过某一设定值时, 单元自动转化为SPH粒子, 并引入有限单元-粒子接触算法和耦合算法, 实现大变形和破碎区域采用SPH方法计算, 克服有限元法单元畸变存在的问题。多层间隔靶侵彻算例分析表明, 三维FE-SPH耦合计算程序采用等效塑性应变作为转化判据计算结果较稳定, 并且能够有效地再现侵彻过程中所产生的碎片, 能够模拟侵彻碎片对后层靶的毁伤效应。     

PTT黏弹性流体的光滑粒子动力学方法模拟

运用光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)方法对基于PHan-Thien-Tanner (PTT)模型的黏弹性流动进行了数值模拟. 首先, 利用SPH方法模拟了基于PTT模型的平板 Poiseuille流, 通过与文献结果的比较, 验证了SPH方法模拟黏弹性流动的准确性和有效性; 随后, 基于PTT模型对黏弹性自由表面流-液滴碰撞问题进行了SPH模拟, 研 究了PTT模型中拉伸参数对碰撞过程的影响. 为了解决张力不稳定问题, 采用简化的 人工应力公式. 数值结果表明, SPH方法可有效而灵活地模拟黏弹自由表面流问题.     

流固耦合破坏分析的多分辨率PD-SPH方法

流固耦合破坏是一类涉及结构变形与破坏以及复杂自由表面现象的强非线性力学问题. 结合近场动力学(peridynamics, PD)与光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)各自的优势并考虑其计算效率问题, 提出一种适用于分析流?固耦合破坏问题的多分辨率PD-SPH混合方法. 分别采用SPH和PD方法以不同的空间和时间分辨率对流体和结构进行离散与求解, 利用具有与流体粒子相同光滑长度的虚粒子处理流?固界面, 以高精度满足界面边界条件. 通过两个经典算例: 液柱静压力下弹性板的变形和溃坝流体冲击弹性闸门的变形问题, 表明提出的多分辨率PD-SPH方法兼具较高的计算精度和计算效率; 对含裂缝的Koyna重力坝水力劈裂问题进行模拟计算, 所得裂缝扩展路径与文献结果吻合, 说明该方法适用于涉及结构破坏的流固耦合问题仿真. 最后尝试采用该方法进行流体冲击作用下含裂纹混凝土板崩塌过程数值仿真, 准确描述混凝土板的断裂破坏和全过程中的流体运动. 多分辨率PD-SPH混合方法或可为流?固耦合作用下的结构损伤破坏仿真提供一种新选择.      

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

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

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

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

非等温黏弹性复杂流动的改进SPH方法模拟

非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.     

结构爆炸毁伤的浸没多介质有限体积物质点法

结构在爆炸载荷作用下的毁伤现象涉及强非线性激波、固体结构极端变形和破坏破碎、强流固耦合, 给数值计算方法带来了极大的困难与挑战. 针对结构爆炸毁伤问题, 建立了浸没多介质有限体积物质点法(iMMFV-MPM), 采用基于黎曼求解器的多介质有限体积法(MMFVM)模拟爆炸产物和空气的多介质流体, 采用物质点法(MPM)模拟固体结构, 并将提出的基于拉格朗日乘子的连续力浸没边界法(lg-CFIBM)扩展到多介质流体中以处理流固耦合边界条件. 该算法在每个时间步严格满足流固耦合界面处的速度边界条件及动量守恒方程, 不需要重构流固耦合界面, 能够有效地模拟近场爆炸下爆炸产物与结构的相互作用、激波与结构的相互作用和演化以及结构的动态断裂和拓扑变化. 利用iMMFV-MPM对近场爆炸下方形钢筋混凝土靶板的失效模式、外爆载荷下建筑物的毁伤现象以及多腔室内爆炸试验进行了模拟, 模拟结果与相关实验数据吻合良好, 验证了所建立的流固耦合算法的有效性及精度.      

近场水下爆炸瞬态强非线性流固耦合无网格数值模拟研究

近场水下爆炸涉及多相流体的掺杂耦合以及结构的大变形、损伤和断裂等瞬态强非线性现象, 传统的网格算法在模拟近场水下爆炸时面临结构网格畸变、多相界面捕捉精度不足等难题, 鉴于此, 本文建立了完全无网格的近场水下爆炸冲击波和气泡全物理过程瞬态强非线性流固耦合动力学模型. 流体采用基于黎曼求解器的光滑粒子流体动力学(SPH)方法求解, 结构采用重构核粒子法(RKPM)求解, 并基于法向通量边界条件实现流固耦合. 为提高SPH对流场间断的求解精度, 引入黎曼问题思想并结合MUSCL重构算法, 为解决流场粒子体积变化剧烈导致的精度下降问题, 应用了自适应粒子分割与合并方法. 为模拟水下爆炸对结构造成的损伤断裂, 基于退化实体几何表述, 采用Lemaitre损伤算法, 建立了RKPM壳结构断裂损伤模型. 依据所建立的SPH-RKPM流固耦合模型, 对近场水下爆炸冲击波传播、气泡脉动与射流以及结构毁伤进行了模拟, 将得到的冲击波载荷、气泡演化以及结构响应与实验值和其他数值解对比, 验证了当前建立的SPH-RKPM流固耦合模型的有效性和精度, 并给出了水下爆炸载荷特性及其对结构的流固耦合毁伤机制与规律, 旨在为近场水下爆炸载荷预报提供理论和基础性技术支撑, 为毁伤威力评估和舰船防护结构设计提供参考.      

含液颗粒材料液固耦合分析的离散颗粒模型及特征线SPH法

对含液颗粒材料流固耦合分析建议了一个基于离散颗粒模型与特征线SPH法的显式拉格朗日-欧拉无网格方案。在已有的用以模拟固体颗粒集合体的离散颗粒模型[1]基础上,将颗粒间间隙内的流体模型化为连续介质,对其提出并推导了基于特征线的SPH法。数值例题显示了所建议方案在模拟颗粒材料与间隙流相互作用的能力和性能以及间隙流体对颗粒结构承载能力及变形的影响。     

Efficient meshless SPH method for the numerical modeling of thick shell structures undergoing large deformations

The objective of this paper is to present an extension of the Lagrangian Smoothed Particle Hydrodynamics (SPH) method to solve three-dimensional shell-like structures undergoing large deformations. The present method is an enhancement of the classical stabilized SPH commonly used for 3D continua, by introducing a Reissner–Mindlin shell formulation, allowing the modeling of moderately thin structure using only one layer of particles in the shell mid-surface. The proposed Shell-based SPH method is efficient and very fast compared to the classical continuum SPH method. The Total Lagrangian Formulation valid for large deformations is adopted using a strong formulation of the differential equilibrium equations based on the principle of collocation. The resulting non-linear dynamic problem is solved incrementally using the explicit time integration scheme, suited to highly dynamic applications. To validate the reliability and accuracy of the proposed Shell-based SPH method in solving shell-like structure problems, several numerical applications including geometrically non-linear behavior are performed and the results are compared with analytical solutions when available and also with numerical reference solutions available in the literature or obtained using the Finite Element method by means of ABAQUS^© commercial software.     

叶轮机械迴转面流场的有限元近似解方法的研究和数值计算

本文提出一种新的概念:即在应用有限元处理不同类型的流体力学、传热学问题时,应该根据各自问题的物理本质,采用不同型式的基函数(shape function),并在文[1]的基础上,应用有限元近似解方法(FEAM)对直角弯管流动进行了数值计算和精度分析,还对带长短叶片的前向离心通风机的迴转面流场进行了数值计算和分析,在迴转面流场计算中采用了一种新的简化方法,解决了在非矩形单元上求取近似解的困难,计算结果表明在稀疏网格下能保持较高精度,从而明显地提高了计算速度。     

Fluctuation Splitting Schemes for the Compressible and Incompressible Euler and Navier-Stokes Equations

Introduced in the late eighties by Roe, fluctuation splitting (or residual distribution) schemes have recently emerged as a viable alternative to Finite Volume and Finite Element methods for PDE based, fluid dynamics simulations using unstructured meshes. Their application to the numerical approximation of the compressible and incompressible Euler and Navier-Stokes equations is described, emphasizing low Mach number and incompressible applications. The advantages provided by time-preconditioning techniques are discussed and details of the implementation are given.     

A modification of Taylor-Galerkin Finite Element Method and its application

Two basic hypothesises of Taylor-Galerkin Finite Element Method are studied in this paper.One of them which is unreasonable is redefined.The only hypothesis becomes the standpoint of Generalized Finite Element.We use this idea to analysis stream function-vorticity equations with Modified Taylor-Galerkin Finite Element Method,and give the two-step solving method,which makes the solving process more reasonable than ever before.Several computational examples reveal that the results of this new method are satisfied.     

Development and validation of off-road tire-gravelly soil interaction using advanced computational techniques

This paper presents a novel modelling technique to compute the interaction between an 8x4 off-road truck and gravelly soil (sand with gravel soil). The off-road truck tire size 315/80R22.5 is modelled using the Finite Element Analysis (FEA) technique and validated using manufacturer-provided data in static and dynamic responses. The gravelly soil is modelled using Smoothed-Particle Hydrodynamics (SPH) technique and calibrated against physical measurements using pressure-sinkage and direct shear-strength tests. The tire-gravelly soil interaction is captured using the node symmetric node to segment with edge treatment algorithm deployed for interaction between FEA and SPH elements. The model setup consists of four tires presenting the four axles of the truck, the first tire is a free-rolling steering tire, the second and third tires are driven tires and the fourth tire is a free-rolling push tire. The truck tires-gravelly soil interaction is computed and validated against physical measurements performed in Göteborg, Sweden. The effect of gravelly soil compaction and truck loading on the tire performance is discussed and investigated.     

Examining a modified algorithm of smoothed particle hydrodynamics for a high velocity perforation of an aluminum beam

The wide spreading of utilizing of smoothed particle hydrodynamics (SPH) for numerical studies of the complex and high rate deformations of continuums, led the current study to gain a more reliable simulation by employing a modified compressible smoothed particle hydrodynamics (MCSPH) algorithm which could be a more accurate and stable technique in high tension regions, in despite of incompressible standard SPH. The main feature of the modified compressible SPH algorithm relies on a three steps solution procedure to calculate the pressure gradient, the deviatoric stress tensor, and the body forces separately. This algorithm is free of any artificial viscosity in its formulations, as well as welcoming to compressible effects which permits the pressure shock waves in high rate plastic deformation. To examine the accuracy of the algorithm, a benchmark problem of colliding rubber cylinders was simulated first and then a high velocity perforation process of an aluminum beam struck by a rigid projectile was simulated in various projectile speeds, and the failure response of the beam in each case was accompanied by crack propagation process. The prominent capability of the utilized MCSPH can be more illustrated when it was used in simulation of thickness crack propagation a tiny crack paths and defragmentation which can be encountered as a not easy numerical case study. The adequate assurance has been more fortified when the results were compared to those reported from a Finite Element method study.     

A regular boundary element method for fluid flow

The Boundary Element Method is now well established as a valid numerical technique for the solution of field problems, equal to the Finite Element Method in generality and surpassing it in computational efficiency in some cases.¹ In this paper is presented a 'Regular Boundary Element Method' as applied to inviscid laminar fluid flow problems. It involves the formation of a system of regular integral equations obtained by moving the singularity outside the domain of the given problem. It is also shown that non-conforming elements may be used whereby freedoms are not defined at the geometric nodes under the boundary element discretization. A linear element is developed here; higher order variants could easily be defined. Satisfactory numerical results have been obtained using the proposed regular method with both conventional (continuous across the boundary) and non-conforming boundary elements for two-dimensional inviscid laminar fluid flow problems having regular and singular solutions.     

Immersed finite element method for fluid-structure interactions

In this paper, we present a detailed derivation of the numerical method, Immersed Finite Element Method (IFEM), for the solution of fluid-structure interaction problems. This method is developed based on the Immersed Boundary (IB) method that was initiated by Peskin, with additional capabilities in handling nonuniform and independent meshes and applying arbitrary boundary conditions on both fluid and solid domains. A higher order interpolation function is adopted from one of the mesh-free methods, the Reproducing Kernel Particle Method (RKPM), which relieves the uniformity constraint of fluid meshes. Two 2-D example problems are presented to illustrate the capabilities of the algorithm. The accuracy in the numerical analysis demonstrates that the IFEM algorithm is a reliable and robust numerical approach to solve fluid and deformable solid interactions.     

Stability of a spinning liquid-filled spacecraft

Summary The stability of a spinning liquid-filled spacecraft has been investigated in the present paper. Using Galerkin's method, the attitude dynamic equations have been given. The Liapunov direct method was employed to obtain a sufficient condition for stability. Three kinds of characteristic modals were investigated: free motion of inviscid fluid, slosh motion and non-slosh motion. All characteristic problems can be solved numerically by the Finite Element Method or the Boundary Element Method. It has been demonstrated that the viscosity of the fluid has a dissipative effect at large Reynolds number, while the slosh motion plays a destabilizing role. The non-slosh model of fluid does not affect the stability criterion. Accepted for publication 19 October 1996     

An improved KGF‐SPH with a novel discrete scheme of Laplacian operator for viscous incompressible fluid flows

The kernel gradient free (KGF) smoothed particle hydrodynamics (SPH) method is a modified finite particle method (FPM) which has higher order accuracy than the conventional SPH method. In KGF‐SPH, no kernel gradient is required in the whole computation, and this leads to good flexibility in the selection of smoothing functions and it is also associated with a symmetric corrective matrix. When modeling viscous incompressible flows with SPH, FPM or KGF‐SPH, it is usual to approximate the Laplacian term with nested approximation on velocity, and this may introduce numerical errors from the nested approximation, and also cause difficulties in dealing with boundary conditions. In this paper, an improved KGF‐SPH method is presented for modeling viscous, incompressible fluid flows with a novel discrete scheme of Laplacian operator. The improved KGF‐SPH method avoids nested approximation of first order derivatives, and keeps the good feature of ‘kernel gradient free’. The two‐dimensional incompressible fluid flow of shear cavity, both in Euler frame and Lagrangian frame, are simulated by SPH, FPM, the original KGF‐SPH and improved KGF‐SPH. The numerical results show that the improved KGF‐SPH with the novel discrete scheme of Laplacian operator are more accurate than SPH, and more stable than FPM and the original KGF‐SPH. Copyright © 2015 John Wiley & Sons, Ltd.