光滑粒子法中的一种新的核函数

分析了传统的核函数产生压缩失稳现象的原因,提出了消除这种压缩失稳现象的一种新的核函数。采用改进的光滑粒子法,对几种常用的核函数进行了一维应变波的对比计算。结果表明:所提出新的核函数在应力波计算中不但保证了计算精度,还能有效地消除压缩失稳。     

应力波数值计算中的SPH方法

对一维波动方程的SPH(smoothed particle hydrodynamics)格式和有限差分格式进行比较,并采用SPH法模拟了一维应力/应变波, 获得1个可衡量SPH法模拟应力波准确性的重要指标。结果表明,SPH法模拟应力波传播中采用的光滑长度必须不小于粒子间距;采用B-样条核函数和高斯型核函数能够获得良好的应力波图像,而二次型核函数不能,因此二次型核函数不适用于冲击动力学的数值计算。     

弹塑性梁在轴向应力波下的动态屈曲

本文考虑轴向应力波效应，利用分叉理论研究各种支承半无限长弹塑性梁的动态屈曲问题。在轴向阶梯载荷和脉冲载荷冲击下得到了梁的临界屈曲载荷及初始屈曲模态。其结果与实验现象相一致。同时也为研究结构动态屈曲问题提供了有效途径。     

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

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

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

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

基于修正光滑粒子流体动力学算法的低能量耗散数值波浪水池开发

目前, 无网格光滑粒子流体动力学SPH粒子法在波浪与结构物相互作用研究方面得到广泛应用, 但该方法模拟波浪远距离传播时, 常常面临严重的能量耗散问题, 导致波高非物理性降低, 给大范围海域、长时间作用下的波-物耦合作用研究带来一定困难. 对此, 本文采用一种核函数修正算法, 在确保粒子间相互作用对称性的同时, 改进压力梯度离散项的计算精度, 设法解决SPH方法中能量非物理性耗散的难题. 相较于前人减缓能量非物理性衰减的方法, 本文的修正SPH算法避免了自由液面搜索等复杂处理过程, 并能保证动量守恒特性. 数值结果中, 采用振荡液滴、规则波、不规则波等算例, 验证本修正SPH算法的准确性和有效性. 结果表明, 该修正SPH算法能准确模拟振荡液滴形态变化, 且动能保持较好守恒性. 通过数值水池与物理水池两者规则波与不规则波结果的对比分析表明, 基于本文修正SPH算法建立的数值波浪水池具有较好的抗能量衰减效果, 能实现长时间、远距离波浪传播的准确模拟. 此外, 本算法能在低光滑长度系数条件下, 实现精确模拟, 将极大缩减三维SPH模拟的时间, 从而节约计算成本.      

轴向冲击载荷作用下锥壳中弹性应力波传播的计算和实验研究

本文利用有限元分析和模型实验研究了在轴向冲击载荷作用下,锥壳中弹性应力波的传播、计算和实验结果表明,结构中存在着弹性纵波和弹性弯曲波的传播,它们传播的速度各不相同,使壳面承受不同的应力状态;讨论了纵波和弯曲波随壳面的衰减;实验指出,由于边界的影响,即使纵波的反射也会产生新的反射弯曲波沿锥面传播。     

弹性圆柱壳在轴向应力波传播过程中的非对称屈曲

讨论弹性圆柱壳端部受冲击载荷作用，在应力波传播过程中的非对称屈曲问题。通过求解扰动方程得到了动态屈曲的分叉条件、临界载荷和屈曲模态。数值结果表明，当壳壁厚不很薄时，轴对称屈曲临界载荷比非对称临界载荷高；反之，轴对称临界载荷会比非对称临界载荷低。不同的冲击载荷，屈曲模态也将不同。     

光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用

扼要讨论了光滑粒子法的离散思想,充分利用光滑粒子法和有限元方法各自的优点,提出了一种初始时刻用有限元建模,计算过程中大变形单元自动转换为光滑粒子的耦合算法。高速碰撞的系列算例说明,耦合算法不但适宜于计算大变形冲击动力学问题,而且由于集两种方法的优点于一身,可以更高效地模拟一些高速碰撞问题,提高计算效率。     

光滑函数法求解拉压不同弹性模量问题

采用光滑函数技术,对拉压不同弹性模量问题的应力应变关系进行光滑处理,可避免迭代中应力状态的判断,方便计算。同时建立了相应的基于初应力技术的有限元计算模式,仅需对刚度阵三角化一次,避免了考虑剪切刚度带来的不便。文中通过不同算例,对所提算法进行了数值验证,与解析解相比有很好符合。此外,对不同拉压模量的热应力分析进行了初步探讨。     

Effect of rock shapes on brittle fracture using Smoothed Particle Hydrodynamics

Breakage of rocks or particulates plays a major role in various industries, such as mineral and ore processing. Many of the processes used for fracturing materials in these industries have the requirement to produce specified size and/or shape of the products. Numerical modelling can assist in understanding and predicting complex fracture processes, and can be used in designing the equipment and setting the process parameters to ensure desired product quality. In this paper, a mesh-free numerical method, called Smoothed Particle Hydrodynamics (SPH), is extended to predict impact fracture of rocks. SPH is a particle based Lagrangian method which is particularly suited to the analysis of fracture due to its capacity to model large deformation and track the free surfaces generated. A continuum damage model is used to predict the fracture of rocks. Evolution of damage is predicted using the strain history of each particle. Damage inhibits the transmission of tensile stress between particles, and once it reaches unity, the particle is unable to transmit tensile stress, resulting in a macro-crack. Connected macro-cracks lead to complete fragmentation.Firstly, an Unconfined Compressive Strength (UCS) test under uniaxial compression of a rock sample is modelled using SPH and compared against experiments to validate the capability of SPH for prediction of fracture in rocks. The SPH prediction matched the well-known experimentally observed diagonal fracture pattern. SPH is subsequently used to simulate brittle fracture of rocks during impact. Rock specimens of different shapes are examined to determine the effects of shape on both the fracture pattern and the energy dissipation during impact fracture. Rock shape is found to have considerable influence on the fracture process, fragment sizes, energy dissipation, and post-fracture motion of the fragments.     

Smoothed Particle Hydrodynamics Modeling of Transverse Flow in Randomly Aligned Fibrous Porous Media

The Lagrangian smoothed particle hydrodynamics (SPH) method is employed to obtain a meso-/micro-scopic pore-scale insight into the transverse flow across the randomly aligned fibrous porous media in a 2D domain. Fluid is driven by an external body force, and a square domain with periodic boundary conditions imposed at both the streamwise and transverse flow direction is assumed. The porous matrix is established by randomly embedding a certain number of fibers in the square domain. Fibers are represented by position-fixed SPH particles, which exert viscous forces upon, and contribute to the density variations of, the nearby fluid particles. An additional repulsive force, similar in form to the 12-6 Lennard-Jones potential between atoms, is introduced to consider the no-penetrating restraint prescribed by the solid pore structure. This force is initiated from the fixed solid material particle and may act on its neighboring moving fluid particles. Fluid flow is visualized by plotting the local velocity vector field; the meandering fluid flow around the porous microstructures always follow the paths of least resistance. The simulated steady-state flow field is further used to calculate the macroscopic permeability. The dimensionless permeability (normalized by the squared characteristic dimension of the fiber cross section) exhibits an exponential dependence on the porosity within the intermediate porosity range, and the derived dimensionless permeability—porosity relation is found to have only minor dependence on either the relative arrangement condition among fibers or the fiber cross section (shape or area).     

Modeling of droplet collisions using Smoothed Particle Hydrodynamics

We present an approach to model collisions of different droplets using Smoothed Particle Hydrodynamics (SPH). We consider bouncing and coalescence of two droplets. We only discretize the droplets neglecting the gaseous phase and consider a free surface at the boundaries. We use a modified continuum surface force model for the surface tension at a free surface. The transition between bouncing and coalescence is modeled using a critical Weber number and calculating the loss of kinetic energy during the collision to determine the point of coalescence. We demonstrate numerical convergence and analyze the error of the method for the transition of bouncing and coalescence. We show that the proposed approach is applicable to weakly-compressible SPH and incompressible SPH and compare binary collisions of Newtonian droplets with experimental results from the literature. Finally we apply the model to non-Newtonian droplets that show shear-thinning and shear-thickening behavior and discuss the differences to Newtonian droplets.     

Smoothed Particle Hydrodynamics Model for Diffusion through Porous Media

A smoothed particle hydrodynamics (SPH) model is presented for the study of diffusion in spatially periodic porous media. The method of SPH is formulated to solve the convection–diffusion equation for tracer diffusion under steady state and transient conditions. Solutions obtained using SPH are compared with other available solutions and the model is used to calculate diffusion coefficients of spatially periodic porous media for the steady state diffusion problem. Diffusion coefficients are then used to calculate nondimensional diffusivities of the media. The effects of media properties on the values of nondimensional diffusivity are also presented.     

基于光滑粒子流体动力学的滑油泵数值仿真

基于光滑粒子流体动力学方法,构建齿轮泵壳体及内部流体的粒子模型,对内啮合齿轮泵在不同工况下的流量特性进行数值模拟。首先设置了均匀分布的油泵出入口压强,计算得到的流量结果与试验结果吻合,且在中低转速下流量与转速呈线性关系;针对高转速下油泵流量降低的问题,通过适当减小周期模型内部的流体粒子数反映流量降低,获得了与试验值相吻合的连续转速流量的模拟结果;针对影响齿轮泵性能的间隙和空化现象,结合模型特点给出了相应的近似处理方法。通过以上研究,将SPH方法成功地应用于滑油泵问题的分析计算。     

Smoothed Particle Hydrodynamics Model of Non-Aqueous Phase Liquid Flow and Dissolution

A smoothed particle hydrodynamics model was developed to simulate the flow of mixtures of aqueous and non-aqueous phase liquids in porous media and the dissolution of the non-aqueous phase in the aqueous phase. The model was used to study the effects of pore-scale heterogeneity and anisotropy on the steady state dense non-aqueous phase liquid (DNAPL) saturation when gravity driven DNAPL displaces water from initially water saturated porous media. Pore-scale anisotropy was created by using co-oriented non-overlapping elliptically shaped grains to represent the porous media. After a steady state DNAPL saturation was reached, water was injected until a new steady state DNAPL saturation was reached. The amount of trapped DNAPL was found to be greater when DNAPL is displaced in the direction of the major axes of the soil grains than when it is displaced in the direction of the minor axes of the soil grains. The amount of trapped DNAPL was also found to increase with decreasing initial saturation of the continuous DNAPL phase. For the conditions used in our simulations, the saturation of the trapped DNAPL with a smaller initial DNAPL saturation was more than 3 times larger than the amount of trapped DNAPL with a larger initial saturation. These simulations were carried out assuming that the DNAPL did not dissolve in water. Simulations including the effect of dissolution of DNAPL in the aqueous phase were also performed, and effective (macroscopic) mass transfer coefficients were determined. The U.S. Government’s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

The propagation of elastic stress waves in a conical shell under axial impulsive loading

The propagation of elastic stress waves in a conical shell subjected to axial impulsive loading is studied in this paper by means of the finite element calculation and model experiments. It is shown that there are two axisymmetrical elastic stress waves propagating with different velocities, i.e., the longitudinal wave and the bending wave. The attenuation of these waves while propagating along the shell surface is discussed. It is found in experiments that the bending wave is also generated when a longitudinal wave reflects from the fixed end of the shell, and both reflected waves will separate during the propagation due to their different velocities. Southwest Institute of Structural Mechanics