Normalized SPH without boundary deficiency and its application to transient solid mechanics problems

Smoothed particle hydrodynamics (SPH) method is a powerful tool for modeling solid mechanics problems, especially for large deformation problems. However, it suffers from boundary deficiency and difficulty of boundary condition treatment. In this work, a normalized SPH method is proposed to overcome these problems. The method is based on a newly developed normalized particle approximation. To derive this particle approximation, a normalized kernel approximation which is accurate for derivatives of linear functions everywhere in a problem domain is constructed, and all integral terms of the normalized kernel approximation including boundary terms are discretized by particle summations. The normalized particle approximation is free of matrix inversion, consequently attractive in computational stability and simplicity compared with other corrective particle approximations. Its approximation accuracy is demonstrated by calculating derivatives of test functions. Based on this particle approximation, the formulation of the normalized SPH method for transient solid mechanics problems is derived. Moreover, a direct method of treating traction boundary conditions is presented by making use of the boundary term of the normalized particle approximation. The accuracy and capability of the normalized SPH method are validated by the calculation of elastic wave propagation in solids and compared with commonly used SPH method.

     

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.     

Adaptive smoothed particle hydrodynamics for high strain hydrodynamics with material strength

This paper presents the implementation of an adaptive smoothed particle hydrodynamics (ASPH) method for high strain Lagrangian hydrodynamics with material strength. In ASPH, the isotropic kernel in the standard SPH is replaced with an anisotropic kernel whose axes evolve automatically to follow the mean particle spacing as it varies in time, space, and direction around each particle. Except for the features inherited from the standard SPH, ASPH can capture dimension-dependent features such as anisotropic deformations with a more generalized elliptical or ellipsoidal influence domain. Two numerical examples, the impact of a plate against a rigid surface and the penetration of a cylinder through a plate, are investigated using both SPH and ASPH. The comparative studies show that ASPH has better accuracy than the standard SPH when being used for high strain hydrodynamic problems with inherent anisotropic deformations. PACS 46.15.-x, 83.10.Rs, 83.50.-v     

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.     

The Kirchhoff–Helmholtz integral for anisotropic elastic media

The Kirchhoff–Helmholtz integral is a powerful tool to model the scattered wavefield from a smooth interface in acoustic or isotropic elastic media due to a given incident wavefield and observation points sufficiently far away from the interface. This integral makes use of the Kirchhoff approximation of the unknown scattered wavefield and its normal derivative at the interface in terms of the corresponding quantities of the known incident field. An attractive property of the Kirchhoff–Helmholtz integral is that its asymptotic evaluation recovers the zero-order ray theory approximation of the reflected wavefield at all observation points where that theory is valid. Here, we extend the Kirchhoff–Helmholtz modeling integral to general anisotropic elastic media. It uses the natural extension of the Kirchhoff approximation of the scattered wavefield and its normal derivative for those media. The anisotropic Kirchhoff–Helmholtz integral also asymptotically provides the zero-order ray theory approximation of the reflected response from the interface. In connection with the asymptotic evaluation of the Kirchhoff–Helmholtz integral, we also derive an extension to anisotropic media of a useful decomposition formula of the geometrical spreading of a primary reflection ray.     

Directional dependent Green’s function and Kelvin images

The inverse of a linear differential operator is an integral operator with a kernel which is commonly known as the Green’s function of the differential operator. Therefore, the knowledge of the Green’s function of a linear problem leads directly to an integral representation of its solution. Any Green’s function is split into a singular part that carries the localized singularity of the Dirac measure and a regular part that is controlled by the Dirichlet boundary condition. In some relatively simple cases, this regular part can be interpreted as the contribution of imaginary sources which lie in the complement of the fundamental domain. If a problem is associated with the Laplace operator, such as the biharmonic operator or the Papkovitch potentials, which both govern Linear Elastostatics, the construction of such Green’s functions are of extremely large importance. All these are well-behaving procedures as long as we live in the highly symmetric geometry represented by the spherical system. But, if we live in a directional-dependent environment, such as the one imposed by the ellipsoidal geometry, the above procedures become extremely complicated, if not impossible. In the present work, the Green’s functions and their Kelvin image systems are obtained for the interior and the exterior regions of an ellipsoid. It is amazing, although not unjustified, that besides the point image source, that is needed for the isotropic spherical case, in the case of ellipsoidal domains, the necessary image system involves a full two-dimensional distribution of imaginary sources to account for the anisotropic character of the ellipsoidal domains.     

On linking n-dimensional anisotropic and isotropic Green’s functions for infinite space,half-space,bimaterial, and multilayer for conduction

We establish exact mathematical links between the n-dimensional anisotropic and isotropic Green’s functions for diffusion phenomena for an infinite space, a half-space, a bimaterial and a multilayered space. The purpose of this work is not to attempt to present a solution procedure, but to focus on the general conditions and situations in which the anisotropic physical problems can be directly linked with the Green’s functions of a similar configuration with isotropic constituents. We show that, for Green’s functions of an infinite and a half-space and for all two-dimensional configurations, the exact correspondences between the anisotropic and isotropic ones can always be established without any regard to the constituent conductivities or any other information. And thus knowing the isotropic Green’s functions will readily provide explicit expressions for anisotropic Green’s functions upon back transformation. For three- and higher-dimensional bimaterials and layered spaces, the correspondence can also be found but the constituent conductivities need to satisfy further algebraic constraints. When these constraints are fully satisfied, then the anisotropic Green’s functions can also be obtained from those of the isotropic ones, or at least in principle.     

Nonlinear surface waves in media with complex rheology

Weak nonlinear waves in a generalized viscoelastic medium with internal oscillators are considered. The rheological relations contain higher time derivatives of the stresses and strains as well as their tensor products. The method of expansion in a small parameter with the introduction of slow time and a running space coordinate is employed. The first approximation gives wave velocities and relations between the parameters equivalent to the results of an acoustic analysis at elastic wave fronts [1]. The second approximation leads to an evolution equation for the displacement velocity. For this a Fourier-Laplace double integral transformation is used. Reversion to the inverse transforms of the unknown functions leads to an integrodifferential evolution equation, which contains a Hubert transform and is a generalization of the Benjamin-Ono equation of deep water theory.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 95–103, September–October, 1990.     

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.     

A kernel gradient free (KGF) SPH method

The finite particle method (FPM) is a modified SPH method with high order accuracy while retaining the advantages of SPH in modeling problems with free surfaces, moving interfaces, and large deformations. In both SPH and FPM, kernel gradient is necessary in kernel and particle approximation of a field function and its derivatives. In this paper, a new FPM is presented, which only involves kernel function itself in kernel and particle approximation. The kernel gradient is not necessary in the whole computation, and this approach is thus referred to as a kernel gradient free (KGF) SPH method. This is helpful when a kernel function is not differentiable or the resultant kernel gradient is not sufficiently smooth, and thus it is more general in selecting a kernel function. Moreover, different from the original FPM with an asymmetric corrective matrix, in the new FPM, the resultant corrective matrix is symmetric, and this is advantageous in particle approximations. A series of numerical examples have been conducted to show the efficiencies of KGF‐SPH including one‐dimensional mathematical tests of polynomial functions with equal or variable smoothing length and two‐dimensional incompressible fluid flow of shear cavity. It is found that KGF‐SPH is comparable with FPM in accuracy and is flexible as SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

Partial Uniqueness and Obstruction to Uniqueness in Inverse Problems for Anisotropic Elastic Media

We consider the inverse problem of identifying the density and elastic moduli for three-dimensional anisotropic elastic bodies, given displacement and traction measurements made at their surface. These surface measurements are modelled by the dynamic Dirichlet-to-Neumann map on a finite time interval. For linear or nonlinear anisotropic hyperelastic bodies we show that the displacement-to-traction surface measurements do not change when the density and elasticity tensor in the interior are transformed tensorially by a change of coordinates fixing the surface of the body to first order. Our main tool, a new approach in inverse problems for elastic media, is the representation of the equations of motion in a covariant form (following Marsden and Hughes, 1983) that preserves the underlying physics.In the case of classical linear elastodynamics we then investigate how the type of anisotropy changes under coordinate transformations. That is, we analyze the orbits of general linear, anisotropic elasticity tensors under the action by pull-back of diffeomorphisms that fix the surface of the elastic body to first order, and derive a pointwise characterization of parts of the orbits under this action. For example, we show that the orbit of isotropic elastic media, at any point in the body, consists of some transversely isotropic and some orthotropic elastic media. We then derive the first uniqueness result in the inverse problem for anisotropic media using surface displacement-traction data: uniqueness of three elastic moduli for tensors in the orbit of isotropic elasticity tensors. ^†Partially supported by an MSRI Postdoctoral Fellowship. Research at MSRI is supported in part by NSF grant DMS-9850361. This work was conducted while the first author was a Gibbs Instructor at Yale University. ^‡Partially supported by an MSRI Postdoctoral Fellowship, and by NSF grant DMS-9801664 (9996350).     

Superposing scheme for the three-dimensional Green’s functions of an anisotropic half-space

This paper presents a method of superposition for the half-space Green’s functions of a generally anisotropic material subjected to an interior point loading. The mathematical concept is based on the addition of a complementary term to the Green’s function in an anisotropic infinite domain. With the two-dimensional Fourier transformation, the complementary term is derived by solving the generalized Stroh eigenrelation and satisfying the boundary conditions on the free surface with the use of Green’s functions in the full-space case. The inverse Fourier transform leads to the contour integrals, which can be evaluated with the application of Cauchy residue theorem. Application of the present results is made to obtain analytical expression for the orthotropic materials which were not reported previously. The closed-form solutions for the transversely isotropic and isotropic materials derived directly from the solutions as being a special case are also given in this paper.     

A finite volume scheme for solving anisotropic diffusion on ALE‐AMR grids

In the context of High Energy Density Physics and more precisely in the field of laser plasma interaction, Lagrangian schemes are commonly used. The lack of robustness due to strong grid deformations requires the regularization of the mesh through the use of Arbitrary Lagrangian Eulerian methods. Theses methods usually add some diffusion and a loss of precision is observed. We propose to use Adaptive Mesh Refinement (AMR) techniques to reduce this loss of accuracy. This work focuses on the resolution of the anisotropic diffusion operator on Arbitrary Lagrangian Eulerian‐AMR grids. In this paper, we describe a second‐order accurate cell‐centered finite volume method for solving anisotropic diffusion on AMR type grids. The scheme described here is based on local flux approximation which can be derived through the use of a finite difference approximation, leading to the CCLADNS scheme. We present here the 2D and 3D extension of the CCLADNS scheme to AMR meshes. Copyright © 2017 John Wiley & Sons, Ltd.     

固体介质中SPH方法的拉伸不稳定性问题研究进展

光滑粒子流体动力学法(smoothed particle hydrodynamics, SPH)是一种基于核估计的无网格Lagrange数值方法.它用粒子方程离散流体动力学的连续方程, 既可以处理有限元难于处理的大变形和严重扭曲问题, 又可以处理有限差分法不易处理的自由边界和材料界面的问题, 在固体力学中的冲击、爆炸和裂纹模拟中具有广阔的发展前景.但是, 该算法的拉伸不稳定性(tensile instability)问题是它在固体力学领域中应用的最大障碍.对SPH稳定性分析表明, 算法不稳定性的条件仅与应力状态和核函数的2阶导数有关.目前, 应力点法(stress points)、Lagrange核函数法、人工应力法(artificialstress)、修正光滑粒子法(corrective smoothed particle method, CSPM)和守恒光滑法(conservativesmoothing)以及其他一些方法成功地改善了SPH的拉伸不稳定性, 但是每一种方法都不能彻底解决SPH的拉伸不稳定性问题.本文介绍了SPH法的方程和Von Neumann稳定性分析的思想, 以及国内外在这几个方面的研究成果及其最新进展, 同时指出目前研究中存在的问题和研究的方向.      

各向异性位势平面问题的规则化边界元法

基于转化域方程为边界积分方程的极限定理及一个新颖的基本解分解技术, 建立间接变量规则化边界积分方程, 它有效地避免了奇异积分的直接计算. 与已有方法比,该方法不将问题变换为各向同性的问题去处理, 因而无需反演运算, 也有别于Galerkin方法, 无需计算重积分. 可计算任意边界位势梯度, 而不仅限于法向通量. 针对椭圆边界的边值问题, 提交一种精确单元来描述边界几何. 数值算例表明, 所提算法稳定且效率高, 所得数值结果与精确解吻合较好.      

层状横观各向同性饱和土的非轴对称动力响应

通过方位角的Fourier变换,将圆柱坐标系下横观各向同性饱和土的Biot非轴对称波动方 程转化为一组一阶常微分方程组. 然后基于径向Hankel变换,建立问题的状态方程;求解状态方程后,得到传递矩阵. 进而利用传递矩阵,结合饱和层状地基的边界条件、排水条件及层间接触和连续条件,求解 了任意震源力作用下层状横观各向同性饱和地基频域动力响应问题. 时域解可通过频率的Fourier积分得到.     

Transient thermoelastic waves in an anisotropic hollow cylinder due to localized heating

This paper presents a theoretical study of transient ultrasonic guided waves generated by concentrated heating of the outer surface of an infinite anisotropic hollow circular cylinder. Generalized thermoelastic theory proposed by Lord and Shulman is adopted to model the dynamic thermoelastic behavior of the cylinder. The concentrated heat source model used is to represent heating due to a pulsed laser beam, which is focused on the outer surface of the cylinder. A semi-analytical finite element (SAFE) method is employed to evaluate guided wave modes in the cylinder. Using integral transform techniques, the modal wave forms are obtained in frequency and wave number domains. Time histories of the propagating modes are then calculated by applying inverse Fourier transformation in the time domain. Numerical results showing the dispersion curves for the group velocities of the propagating modes and transient radial displacements are presented. For this purpose it is assumed that the cylinder is made of transversely isotropic silicon nitride (Si₃N₄). Attention is focused on the propagation characteristics of longitudinal and flexural modes separately.     

Anisotropic yield criterion for polycrystalline metals using texture and crystal symmetries

An anisotropic yield criterion for polycrystalline metals which uses texture data and takes advantage of crystal symmetries is presented. A linear transformation is developed to map an anisotropic yield surface for a polycrystal to an appropriate isotropic yield surface. The transformation developed reflects the symmetry of the material being modeled. First, the transformation is determined. Then, information regarding the orientation distribution (texture) of the crystals in a polycrystalline aggregate is used to determine, via averaging, the transformation for the polycrystal. The transformation, along with appropriate isotropic yield surface, provides a phenomenological approach to modeling yield, yet accounts for microstructural texture. The approach reduces to the Hill (1950) anisotropic plasticity theory under certain conditions. The yield surfaces and R-values for various face-centered-cubic ( fcc) polycrystalline textures are computed by this method. Results compare favorably with those given by other theories, and with experiment. The method proves to have the computational efficiency of phenomenological approaches to modeling yield, while effectively incorporating the physics of more complex crystallographic approaches.     

多圆形刚性夹杂的反平面问题

构造任意分布且相互影响的多个圆形刚性夹杂模型的复应力函数，采用复变函数方法，达到满足各个夹杂的边界条件，利用坐标变换和围线积分将求解方程组化为线性代数方程组，推导出了圆形刚性夹杂任意分布的界面应力解析表达式，算例对多夹杂与单夹杂两种模型的界面应力最大值进行了对比，同时还给出了界面应力最大值随夹杂间距的变化规律，求出了刚性夹杂的合理间距。本文发展的分析方法为研究夹杂材料的细观机理探索了一条有效的分析途径。