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.

     

BOUNDARY ELEMENT METHOD FOR BUCKLING EIGENVALUE PROBLEM AND ITS CONVERGENCE ANALYSIS

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂｕｃｋｌｉｎｇｅｉｇｅｎｖａｌｕｅｐｒｏｂｌｅｍｈａｓｉｍｐｏｒｔａｎｔｓｉｇｎｉｆｉｃａｎｃｅｉｎｔｈｅｓｔａｂｉｌｉｔｙａｎａｌｙｓｉｓｏｆｅｎｇｉｎｅｅｒｉｎｇｓｔｒｕｃｔｕｒｅ.Ｈｅｎｃｅｔｈｅｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｆｏｒｔｈｅｓｅｐｒｏｂｌｅｍｓｉｓｅｘｔｒｅｍｅｌｙｍｅａｎｉｎｇｆｕｌｉｎｃｏｍｐｕｔａｔｉｏｎａｌｍｅｃｈａｎｉｃｓ.ＴｈｅｐｒｅｓｅｎｔｃｏｍｐｕｔａｔｉｏｎａｌｍｅｔｈｏｄｓｆｏｃｕｓｏｎＦＥＭ ,ｄｉｆｆｅｒｅｎｃｅｍ…     

Boundary element method for buckling eigenvalue problem and its convergence analysis

The conditions for determining solution of buckling eigenvalue problem are discussed. The corresponding system of integral equations with constraint conditions and boundary variational equations with Lagrange multiplier are established. The theorems on the existence and uniqueness of the solution for these problems are given. The corresponding boundary element method is constructed and the error estimation for the approximation solution is obtained. Finally the numerical example is given. Foundation item: the National Natural Science Foundation Pre-research Project (T4107015) Biography: Ding Rui (1969-)     

Time‐accurate stabilized finite‐element model for weakly nonlinear and weakly dispersive water waves

Introduction of a time‐accurate stabilized finite‐element approximation for the numerical investigation of weakly nonlinear and weakly dispersive water waves is presented in this paper. To make the time approximation match the order of accuracy of the spatial representation of the linear triangular elements by the Galerkin finite‐element method, the fourth‐order time integration of implicit multistage Padé method is used for the development of the numerical scheme. The streamline‐upwind Petrov–Galerkin (SUPG) method with crosswind diffusion is employed to stabilize the scheme and suppress the spurious oscillations, usually common in the numerical computation of convection‐dominated flow problems. The performance of numerical stabilization and accuracy is addressed. Treatments of various boundary conditions, including the open boundary conditions, the perfect reflecting boundary conditions along boundaries with irregular geometry, are also described. Numerical results showing the comparisons with analytical solutions, experimental measurements, and other published numerical results are presented and discussed. Copyright © 2007 John Wiley & Sons, Ltd.     

The improved element-free Galerkin method for three-dimensional wave equation

The paper presents the improved element-freeGalerkin(IEFG) method for three-dimensional wave propagation.The improved moving least-squares(IMLS) approximation is employed to construct the shape function,whichuses an orthogonal function system with a weight function asthe basis function.Compared with the conventional movingleast-squares(MLS) approximation,the algebraic equationsystem in the IMLS approximation is not ill-conditioned,andcan be solved directly without deriving the inverse matrix.Because there are fewer coefficients in the IMLS than in theMLS approximation,fewer nodes are selected in the IEFGmethod than in the element-free Galerkin method.Thus,theIEFG method has a higher computing speed.In the IEFGmethod,the Galerkin weak form is employed to obtain a discretized system equation,and the penalty method is appliedto impose the essential boundary condition.The traditionaldifference method for two-point boundary value problems isselected for the time discretization.As the wave equationsand the boundary-initial conditions depend on time,the scaling parameter,number of nodes and the time step length areconsidered for the convergence study.     

弹性力学的一种边界无单元法

首先对移动最小二乘副近法进行了研究，针对其容易形成病态方程的缺点，提出了以带权的正交函数作为基函数的方法－改进的移动最小二乘副近法，改进的移动最小二乘逼近法比原方法计算量小，精度高，且不会形成病态方程组，然后，将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合，提出了弹性力学的一种边界无单元法，这种边界无单元法法是边界积分方程的无网格方法，与原有的边界积分方程的无网格方法相比，该方法直接采用节点变量的真实解为基本未知量，是边界积分方程无网格方法的直接解法，更容易引入界条件，且具有更高的精度，最后给出了弹性力学的边界无单元法的数值算例，并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。     

Effective Flux Boundary Conditions for Upscaling Porous Media Equations

We introduce a new algorithm for setting pressure boundary conditions in subgrid simulations of porous media flow. The algorithm approximates the flux in the boundary cell as the flux through a homogeneous inclusion in a homogeneous background, where the permeability of the inclusion is given by the cell permeability and the permeability of the background is given by the ambient effective permeability. With this approximation, the flux in the boundary cell scales with the cell permeability when that permeability is small, and saturates at a constant value when that permeability is large. The flux conditions provide Neumann boundary conditions for the subgrid pressure. We call these boundary conditions effective flux boundary conditions (EFBCs). We give solutions for the flux through ellipsoidal inclusions in two and three dimensions, assuming symmetric tensor permeabilities whose principal axes align with the axes of the ellipse. We then discuss the considerations involved in applying these equations to scale up problems in geological porous media. The key complications are heterogeneity, fluctuations at all length scales, and boundary conditions at finite scales.     

Dynamic stability of rotating cylindrical shells subjected to periodic axial loads

In this paper, the dynamic stability of rotating cylindrical shells under static and periodic axial forces is investigated using a combination of the Ritz method and Bolotin’s first approximation. The kernel particle estimate is employed in hybridized form with harmonic functions, to approximate the 2-D transverse displacement field. A system of Mathieu–Hill equations is obtained through the application of the Ritz energy minimization procedure. The principal instability regions are then obtained via Bolotin’s first approximation. In this formulation, both the hoop tension and Coriolis effects due to the rotation are accounted for. Various boundary conditions are considered, and the present results represent the first instance in which, the effects of boundary conditions for this class of problems, have been reported in open literature. Effects of rotational speeds on the instability regions for different modes are also examined in detail.     

Appropriate boundary conditions for the solution of the equations of unsteady one-dimensional gas flow by the Lax-Wendroff method

This paper presents a new characteristic approximation to the boundary conditions, required in the solution of gas flow problems by the Law-Wendroff method. The accuracy of this and other currently used methods is assessed by a comparison with the exact solutions of two test problems     

Wave propagation and diffraction in bodies with noncircular cylindrical boundaries

Conclusion The foregoing results make it possible to formulate approximate solutions of linear problems (in four approximations) for the steady-state diffraction and propagation of waves of different physical natures in bodies bounded by noncircular cylindrical surfaces. The problem is reduced in each approximation to corresponding problems in a circular cylindrical coordinate system with identical homogeneous equations in all approximations and different right-hand sides of the boundary conditions in each approximation. The results are obtained in a general form for all linear problems of continuum mechanics.Institute of Mechanics, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Prikladnaya Mekhanika, Vol. 9, No. 9, pp. 3–11, September, 1973.     

Deformation problem for an elastically fixed plate modeling a coating partially delaminated from the substrate (Plane Strain)

An asymptotic solution of the problemindicated in the title is obtained at distances large compared with the plate width and some promising methods for its use, in particular, for calculating the coefficients in the boundary conditions of the plate elastic fixation which models a coating partially delaminated from the substrate, are outlined. The possibility of considering the delamination in the approximation of the plate weak bending (the plate approximation) and the possibility of neglecting the tangential stress action along the contact boundary are implemented. The substrate is considered as a half-infinite elastic solid. This solution was obtained by using the Fourier transform and the solution of the resulting equation by the Wiener-Hopf method. The obtained asymptotic solution can be used to study problems related to coating delamination, especially on soft thick substrata.     

Equation of small transverse vibrations of an elastic pipeline filled with a transported fluid

The integro-differential equation of small transverse vibrations of a rectilinear elastic pipeline filled with a transported fluid is obtained. The pipeline vibrations are described in the linear setting in the beam approximation. The mutual dynamic influence of motions of the pipeline and the filling fluid is taken into account. A complete trigonometric series method is presented for solving problems with various boundary and initial conditions for the pipeline deflection.     

不可压缩N-S方程的稳定分步算法及耦合离散

在有限增量微积分(finite increment calculus, FIC)的理论框架下,通过引入一个附加变量,发展了压力稳定型分步算法,有效改善了经典 分步算法的压力稳定性,同时还避免了标准FIC方法中存在的空间高阶导数的计算. 为保证 数值方法同时具有较快的计算速度和较好的健壮性,发展了有限元与无网格的耦合空间离散 方法. 该方案可在网格发生扭曲的区域采用无网格法空间离散以保证求解的精度和稳定性, 而在网格质量较好的区域以及本质边界上保留使用有限元法空间离散以提高计算效率和便于 施加本质边界条件. 方腔流考题的数值模拟结果突出地显示了所发展的压力稳定型分步算 法比经典分步算法具有更好的压力稳定性,能够有效消除速度-压力插值空间违反LBB条件而 导致的压力场的虚假数值振荡. 平面Poisseuille流动和一个典型型腔充填过程的数值模拟 结果, 表明了发展的耦合离散方案相对于单一的有限元法和单一的无网格法在综合考虑计 算效率和算法健壮性方面的突出优点.     

Nonlinear equations of elastic deformation of plates

A method for constructing nonlinear equations of elastic deformation of plates with boundary conditions for stresses and displacements at the face surfaces in an arbitrary coordinate system is proposed. The initial three–dimensional problem of the nonlinear theory of elasticity is reduced to a one–parameter sequence of two–dimensional problems by approximating the unknown functions by truncated series in Legendre polynomials. The same unknowns are approximated by different truncated series. In each approximation, a linearized system of equations whose differential order does not depend on the boundary conditions at the face surfaces which can be formulated in terms of stresses or displacements is obtained.     

Non-element method of solving 2D boundary problems defined on polygonal domains modeled by Navier equation

The paper presents a non-element method of solving boundary problems defined on polygonal domains modeled by corner points. To solve these problems a parametric integral equation system (PIES) is used. The system is characterized by a separation of the approximation of boundary geometry from the approximation of boundary functions. This feature makes it possible to effectively investigate the convergence of the obtained solutions with no need of performing the approximation of boundary geometry. The testing examples included confirm high accuracy of the solutions.     

Improved non-singular local boundary integral equation method

When the source nodes are on the global boundary in the implementation of local boundary integral equation method (LBIEM),singularities in the local boundary integrals need to be treated specially. In the current paper,local integral equations are adopted for the nodes inside the domain trod moving least square approximation (MLSA) for the nodes on the global boundary,thus singularities will not occur in the new al- gorithm.At the same time,approximation errors of boundary integrals are reduced significantly.As applications and numerical tests,Laplace equation and Helmholtz equa- tion problems are considered and excellent numerical results are obtained.Furthermore, when solving the Hehnholtz problems,the modified basis functions with wave solutions are adapted to replace the usually-used monomial basis functions.Numerical results show that this treatment is simple and effective and its application is promising in solutions for the wave propagation problem with high wave number.     

Analyzing the 2D fracture problems via the enriched boundary element-free method

In this paper, the enriched boundary element-free method for two-dimensional fracture problems is presented. An improved moving least-squares (IMLS) approximation, in which the orthogonal function system with a weight function is used as the basis function, is used to obtain the shape functions. The IMLS approximation has greater computational efficiency and precision than the existing moving least-squares (MLS) approximation, and does not lead to an ill-conditioned system of equations. Combining the boundary integral equation (BIE) method and the IMLS approximation, a boundary element-free method (BEFM), for two-dimensional fracture problems is obtained. For two-dimensional fracture problems, the enriched basis function is used at the tip of the crack, and then the enriched BEFM is presented. In comparison with other existing meshless boundary integral equation methods, the BEFM is a direct numerical method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be implemented easily, which leads to a greater computational precision. When the enriched BEFM is used, the singularity of the stresses at the tip of the crack can be shown better than that in the BEFM. For the purposes of demonstration, some selected numerical examples are solved using the enriched BEFM.     

弹性静力问题的无网格弱-强形式结合法

基于改进的移动最小二乘(MLS)二阶导数近似,建立了一种求解弹性静力问题的无网格弱-强形式结合法(MLS-MWS)。该方法采用节点离散求解域,通过MLS构造形函数,将求解域划分为边界域和内部域,并分别使用控制方程的局部弱形式和强形式来建立离散系统方程。对强形式中涉及的近似函数二阶导数计算,提出了一种将其转化为求两次一阶导数的方法,与传统方法相比,该方法计算简单、精度高。MLS-MWS法结合了弱、强形式无网格法的优点,Neumann边界条件容易满足,并且只需在边界区域进行积分。文中应用该方法分析了两个弹性力学平面问题,分析结果表明本文方法具有良好的精度和收敛性。     

A stream function–vorticity formulation‐based immersed boundary method and its applications

A new stream function–vorticity formulation‐based immersed boundary method is presented in this paper. Different from the conventional immersed boundary method, the main feature of the present model is to accurately satisfy both governing equations and boundary conditions through velocity correction and vorticity correction procedures. The velocity correction process is performed implicitly based on the requirement that velocity at the immersed boundary interpolated from the corrected velocity field accurately satisfies the nonslip boundary condition. The vorticity correction is made through the stream function formulation rather than the vorticity transport equation. It is evaluated from the firstorder derivatives of velocity correction. Two simple and efficient ways are presented for approximation of velocity‐correction derivatives. One is based on finite difference approximation, while the other is based on derivative expressions of Dirac delta function and velocity correction. It was found that both ways can work very well. The main advantage of the proposed method lies in its simple concept, easy implementation, and robustness in stability. Numerical experiments for both stationary and moving boundary problems were conducted to validate the capability and efficiency of the present method. Good agreements with available data in the literature were achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

On the no-slip boundary conditions for squeeze flow of viscous fluids

A typical class of boundary conditions for squeeze flow problems in lubrication approximation is the one in which the squeezing rate and the width between the squeezing plates are constant. This hypothesis is justified by claiming that the plates moves so slowly that they can be considered static. In this short note we prove that this assumption leads to a contradiction and hence cannot be used.