Large deformation analysis of rubber based on a reproducing kernel particle method

 A nonlinear formulation of the Reproducing Kernel Particle Method (RKPM) is presented for the large deformation analysis of rubber materials which are considered to be hyperelastic and nearly incompressible. In this approach, the global nodal shape functions derived on␣the basis of RKPM are employed in the Galerkin approximation of the variational equation to formulate the discrete equations of a boundary-value hyperelasticity problem. Existence of a solution in RKPM discretized hyperelasticity problem is discussed. A Lagrange multiplier method and a direct transformation method are presented to impose essential boundary conditions. The characteristics of material and spatial kernel functions are discussed. In the present work, the use of a material kernel function assures reproducing kernel stability under large deformation. Several of numerical examples are presented to study the characteristics of RKPM shape functions and to demonstrate the effectiveness of this method in large deformation analysis. Since the current approach employs global shape functions, the method demonstrates a superior performance to the conventional finite element methods in dealing with large material distortions.     

3D adaptive RKPM method for contact problems with elastic–plastic dynamic large deformation

The adaptive procedure of reproducing kernel particle method (RKPM) for 3D contact problems with elastic–plastic dynamic large deformation is presented. In this study, a modified cell energy error (MCEE) estimate model is constructed to capture the high gradients of stresses behavior in large deformation. Refinement particles with a new proper refinement function are inserted into the high error distribution regions. A domain decomposition method is proposed to determine the support domain size for nodes. A collocation formulation is used in the discretization of the boundary integral of the contact constraint equations formulated by a penalty method. By the use of a particle-to-segment contact algorithm, the contact constraints are imposed directly on the new added contact nodes, consequently the contact forces and their associated stiffness matrices are formulated at the nodal coordinate. For verification of the simulation results, a general benchmark test is applied to justify the accuracy and efficiency of the adaptive RKPM method. Several numerical examples are provided to illustrate the effectiveness and robustness of the suggested approach.     

Quasi-convex reproducing kernel meshfree method

A quasi-convex reproducing kernel approximation is presented for Galerkin meshfree analysis. In the proposed meshfree scheme, the monomial reproducing conditions are relaxed to maximizing the positivity of the meshfree shape functions and the resulting shape functions are referred as the quasi-convex reproducing kernel shape functions. These quasi-convex meshfree shape functions are still established within the framework of the classical reproducing or consistency conditions, namely the shape functions have similar form as that of the conventional reproducing kernel shape functions. Thus this approach can be conveniently implemented in the standard reproducing kernel meshfree formulation without an overmuch increase of computational effort. Meanwhile, the present formulation enables a straightforward construction of arbitrary higher order shape functions. It is shown that the proposed method yields nearly positive shape functions in the interior problem domain, while in the boundary region the negative effect of the shape functions are also reduced compared with the original meshfree shape functions. Subsequently a Galerkin meshfree analysis is carried out by employing the proposed quasi-convex reproducing kernel shape functions. Numerical results reveal that the proposed method has more favorable accuracy than the conventional reproducing kernel meshfree method, especially for structural vibration analysis.     

Synchronized reproducing kernel interpolant via multiple wavelet expansion

In this paper, a new partition of unity – the synchronized reproducing kernel (SRK) interpolant – is derived. It is a class of meshless shape functions that exhibit synchronized convergence phenomenon: the convergence rate of the interpolation error of the higher order derivatives of the shape function can be tuned to be that of the shape function itself. This newly designed synchronized reproducing kernel interpolant is constructed as an series expansion of a scaling function kernel and the associated wavelet functions. These wavelet functions are constructed in a reproducing procedure, simultaneously with the scaling function kernel, by directly enforcing certain orders of vanishing moment conditions. To the authors knowledge, this unique interpolant is the first of its kind to be constructed, and to be used in numerical computations, both in concept and in practice. The new interpolants are in fact a group of special hierarchial meshless bases, and similar counterparts may exist in spline interpolation method, other meshless methods, Galerkin-wavelet method, as well as the finite element method. A detailed account of the subject is presented, and the mathematical principle behind the construction procedure is further elaborated. Another important discovery of this study is that the 1st order wavelet together with the scaling function kernel can be used as a weighting function in Petrov-Galerkin procedures to provide a stable numerical computation in some pathological problems. Benchmark problems in advection-diffusion problems, and Stokes flow problem are solved by using the synchronized reproducing kernel interpolant as the weighting function. Reasonably good results have been obtained. This may open the door for designing well behaved Galerkin procedures for numerical computations in various constrained media.     

板材成形无网格法模拟中的关键技术处理

针对有限元法费时的网格划分/网格重划问题,提出采用无网格法模拟板材成形.采用约束型再生核子法描述板材成形力学方程,仅用一层节点离散板材,同时对形函数进行改造,所得动可容形函数具有插值性,避免了大变形所导致的不稳定形函数及频繁构造形函数,提出了稀疏矩阵数据结构加快接触搜索方法.基于上述措施建立了板材冲压成形无网格数值模拟...     

A meshless collocation method based on the differential reproducing kernel interpolation

A differential reproducing kernel (DRK) interpolation-based collocation method is developed for solving partial differential equations governing a certain physical problem. The novelty of this method is that we construct a set of differential reproducing conditions to determine the shape functions of derivatives of the DRK interpolation function, without directly differentiating the DRK interpolation function. In addition, the shape function of the DRK interpolation function at each sampling node is separated into a primitive function processing Kronecker delta properties and an enrichment function constituting reproducing conditions, so that the nodal interpolation properties are satisfied. A point collocation method based on the present DRK interpolation is developed for the analysis of one-dimensional bar problems, two-dimensional potential problems, and plane problems of elastic solids. It is shown that the present DRK interpolation-based collocation method is indeed a truly meshless approach, with excellent accuracy and fast convergence rate.     

Hermite–Cloud: a novel true meshless method

In this paper, a novel true meshless numerical technique – the Hermite–Cloud method, is developed. This method uses the Hermite interpolation theorem for the construction of the interpolation functions, and the point collocation technique for discretization of the partial differential equations. This technique is based on the classical reproducing kernel particle method except that a fixed reproducing kernel approximation is employed instead. As a true meshless technique, the present method constructs the Hermite-type interpolation functions to directly compute the approximate solutions of both the unknown functions and the first-order derivatives. The necessary auxiliary conditions are also constructed to generate a complete set of partial differential equations with mixed Dirichlet and Neumann boundary conditions. The point collocation technique is then used for discretization of the governing partial differential equations. Numerical results show that the computational accuracy of the Hermite–Cloud method at scattered discrete points in the domain is much refined not only for approximate solutions, but also for the first-order derivative of these solutions.     

Non‐linear version of stabilized conforming nodal integration for Galerkin mesh‐free methods

A stabilized conforming (SC) nodal integration, which meets the integration constraint in the Galerkin mesh‐free approximation, is generalized for non‐linear problems. Using a Lagrangian discretization, the integration constraints for SC nodal integration are imposed in the undeformed configuration. This is accomplished by introducing a Lagrangian strain smoothing to the deformation gradient, and by performing a nodal integration in the undeformed configuration. The proposed method is independent to the path dependency of the materials. An assumed strain method is employed to formulate the discrete equilibrium equations, and the smoothed deformation gradient serves as the stabilization mechanism in the nodally integrated variational equation. Eigenvalue analysis demonstrated that the proposed strain smoothing provides a stabilization to the nodally integrated discrete equations. By employing Lagrangian shape functions, the computation of smoothed gradient matrix for deformation gradient is only necessary in the initial stage, and it can be stored and reused in the subsequent load steps. A significant gain in computational efficiency is achieved, as well as enhanced accuracy, in comparison with the mesh‐free solution using Gauss integration. The performance of the proposed method is shown to be quite robust in dealing with non‐uniform discretization. Copyright © 2002 John Wiley & Sons, Ltd.     

Die shape design optimization of sheet metal stamping process using meshfree method

A die shape design sensitivity analysis (DSA) and optimization for a sheet metal stamping process is proposed based on a Lagrangian formulation. A hyperelasticity‐based elastoplastic material model is used for the constitutive relation that includes a large deformation effect. The contact condition between a workpiece and a rigid die is imposed through the penalty method with a modified Coulomb friction model. The domain of the workpiece is discretized by a meshfree method. A continuum‐based DSA with respect to the rigid die shape parameter is formulated using a design velocity concept. The die shape perturbation has an effect on structural performance through the contact variational form. The effect of the deformation‐dependent pressure load to the design sensitivity is discussed. It is shown that the design sensitivity equation uses the same tangent stiffness matrix as the response analysis. The linear design sensitivity equation is solved at each converged load step without the need of iteration, which is quite efficient in computation. The accuracy of sensitivity information is compared to that of the finite difference method with an excellent agreement. A die shape design optimization problem is solved to obtain the desired shape of the workpiece to minimize spring‐back effect and to show the feasibility of the proposed method. Copyright © 2001 John Wiley & Sons, Ltd.     

Semi-Lagrangian reproducing kernel particle method for fragment-impact problems

Fragment-impact problems exhibit excessive material distortion and complex contact conditions that pose considerable challenges in mesh based numerical methods such as the finite element method (FEM). A semi-Lagrangian reproducing kernel particle method (RKPM) is proposed for fragment-impact modeling to alleviate mesh distortion difficulties associated with the Lagrangian FEM and to minimize the convective transport effect in the Eulerian or Arbitrary Lagrangian Eulerian FEM. A stabilized non-conforming nodal integration with boundary correction for the semi-Lagrangian RKPM is also proposed. Under the framework of semi-Lagrangian RKPM, a kernel contact algorithm is introduced to address multi-body contact. Stability analysis shows that temporal stability of the kernel contact algorithm is related to the velocity gradient between two contacting bodies. The performance of the proposed methods is examined by numerical simulation of penetration processes.     

A reproducing kernel method with nodal interpolation property

A general formulation for developing reproducing kernel (RK) interpolation is presented. This is based on the coupling of a primitive function and an enrichment function. The primitive function introduces discrete Kronecker delta properties, while the enrichment function constitutes reproducing conditions. A necessary condition for obtaining a RK interpolation function is an orthogonality condition between the vector of enrichment functions and the vector of shifted monomial functions at the discrete points. A normalized kernel function with relative small support is employed as the primitive function. This approach does not employ a finite element shape function and therefore the interpolation function can be arbitrarily smooth. To maintain the convergence properties of the original RK approximation, a mixed interpolation is introduced. A rigorous error analysis is provided for the proposed method. Optimal order error estimates are shown for the meshfree interpolation in any Sobolev norms. Optimal order convergence is maintained when the proposed method is employed to solve one‐dimensional boundary value problems. Numerical experiments are done demonstrating the theoretical error estimates. The performance of the method is illustrated in several sample problems. Copyright © 2003 John Wiley & Sons, Ltd.     

Green's functions for Kirchhoff plates of irregular shape

A method is proposed for the construction of Green's functions for the Sophie Germain equation in regions of irregular shape with mixed boundary conditions imposed. The method is based on the boundary integral equation approach where a kernel vector function B satisfies the biharmonic equation inside the region. This leads to a regular boundary integral equation where the compensating loads and moments are applied to the boundary. Green's function is consequently expressed in terms of the kernel vector function B, the fundamental solution function of the biharmonic equation, and kernel functions of the inverse regular integral operators. To compute moments and forces, the kernel functions are differentiated under the integral sign. The proposed method appears highly effective in computing both displacements and stress components.     

A Hermite DRK interpolation-based collocation method for the analyses of Bernoulli–Euler beams and Kirchhoff–Love plates

A Hermite differential reproducing kernel (DRK) interpolation-based collocation method is developed for solving fourth-order differential equations where the field variable and its first-order derivatives are regarded as the primary variables. The novelty of this method is that we construct a set of differential reproducing conditions to determine the shape functions of derivatives of the Hermite DRK interpolation, without directly differentiating it. In addition, the shape function of this interpolation at each sampling node is separated into a primitive function possessing Kronecker delta properties and an enrichment function constituting reproducing conditions, so that the nodal interpolation properties are satisfied for the field variable and its first-order derivatives. A weighted least-squares collocation method based on this interpolation is developed for the static analyses of classical beams and plates with fully simple and clamped supports, in which its accuracy and convergence rate are examined, and some guidance for using this method is suggested.     

Finite deformation formulation of a shell element for problems of sheet metal forming

An elastic-plastic thin shell finite element suitable for problems of finite deformation in sheet metal forming is formulated. Hill's yield criterion for sheet materials of normal anisotropy is applied. A nonlinear shell theory in a form of an incremental variational principle and a quasi-conforming element technique are employed in the Lagrangian formulation. The shell element fulfills the inter-element C ¹ continuity condition in a variational sense and has a sufficient rank to allow finite stretching, rotation and bending of the shell element. The accuracy and efficiency of the finite element formulation are illustrated by numerical examples.     

Reproducing kernel collocation method applied to the non-linear dynamics of pipe whip in a plane

A windowed collocation method, based on a moving least squares reproducing kernel particle approximation of functions, is explored for spatial discretization of the strongly non-linear system of partial differential equations governing large, planar whipping motion of a cantilever pipe subjected to a follower force pulse (the blow-down force) normal to the deflected centreline at its tip. This problem was discussed by Reid et al. [An elastic–plastic hardening–softening cantilever beam subjected to a force pulse at its tip: a model for pipe whip. Proc R Soc London A1998;454:997–1029] where a space–time finite difference discretization was employed to solve the governing partial differential equation of motion. It was shown that, despite the deflected shape predictions being accurate, numerical solutions of these equations might exhibit problematic (possibly spurious) steep localized gradients. The resolution of this problem in the context of structural mechanics is novel and is the subject of this paper. In particular, it is demonstrated that it is possible to reduce significantly such spurious and localized numerical instabilities through a windowed collocation approach with a suitable choice of the window size. The collocation procedure presently adopted is based on the moving least squares reproducing kernel particle method. Material and structural non-linearity in the beam (pipe) model is incorporated via an elastic–plastic-hardening–softening moment–curvature relationship. The projected ordinary differential equations are then integrated in time through a fifth order, explicit Runge–Kutta method with adaptive step sizes.     

Reproducing Kernel Particle Method in Plasticity of Pressure-Sensitive Material with Reference to Powder Forming Process

In this paper, an application of the reproducing kernel particle method (RKPM) is presented in plasticity behavior of pressure-sensitive material. The RKPM technique is implemented in large deformation analysis of powder compaction process. The RKPM shape function and its derivatives are constructed by imposing the consistency conditions. The essential boundary conditions are enforced by the use of the penalty approach. The support of the RKPM shape function covers the same set of particles during powder compaction, hence no instability is encountered in the large deformation computation. A double-surface plasticity model is developed in numerical simulation of pressure-sensitive material. The plasticity model includes a failure surface and an elliptical cap, which closes the open space between the failure surface and hydrostatic axis. The moving cap expands in the stress space according to a specified hardening rule. The cap model is presented within the framework of large deformation RKPM analysis in order to predict the non-uniform relative density distribution during powder die pressing. Numerical computations are performed to demonstrate the applicability of the algorithm in modeling of powder forming processes and the results are compared to those obtained from finite element simulation to demonstrate the accuracy of the proposed model.     

On finite element modelling of surface tension Variational formulation and applications – Part I: Quasistatic problems

This paper describes variational formulation and finite element discretization of surface tension. The finite element formulation is cast in the Lagrangian framework, which describes explicitly the interface evolution. In this context surface tension formulation emerges naturally through the weak form of the Laplace–Young equation.The constitutive equations describing the behaviour of Newtonian fluids are approximated over a finite time step, leaving the governing equations for the free surface flow function of geometry change rather than velocities. These nonlinear equations are then solved by using Newton-Raphson procedure.Numerical examples have been executed and verified against the solution of the ordinary differential equation resulting from a parameterization of the Laplace-Young equation for equilibrium shapes of drops and liquid bridges under the influence of gravity and for various contact angle boundary conditions.     

Solving 2D transient rolling contact problems using the BEM and mathematical programming techniques

This work presents a new approach to the transient rolling contact of two‐dimensional elastic bodies. A solution will be obtained by minimizing a general B‐differentiable function representing the equilibrium equations and the contact conditions at each time step. Inertial effects are not taken into account and the boundary element method is used to compute the elastic influence coefficients of the surface points involved in contact (equilibrium equations). The contact conditions are represented with the help of variational inequalities and projection functions. Finally, the minimization problem is solved using the Generalized Newton's Method with line search. The results are compared with some example problems and the influence of discretization and integration time step on the results is discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

Potential fields caused by point sources in regions with apertures and foreign inclusions

A semi-analytic approach is developed for obtaining potential fields generated by point sources in regions embedded with foreign inclusions and containing holes of various shape. This study promotes an extension of the range of effective implementation of the Green's function version of the boundary integral equation method in mechanics of contemporary materials. Equivalents of Green's functions are obtained to boundary-contact value problems posed for Laplace equation on piecewise homogeneous multiply connected regions. Dirichlet, Neumann and Robin boundary conditions can be imposed on the outer boundary of the region and on contours of the apertures, while the ideal contact conditions are assumed on the interfacial contours. Source points could be located either outside or inside the inclusion. A Green's function modification of the method of functional equations is applied. A number of illustrative examples included to show the potential of the approach.     

基于重构核的最小二乘配点法求解封闭声腔响应

基于重构核思想,应用无网格配点法构造近似函数,并利用最小二乘方法的原理解决边界问题,离散控制微分方程,建立求解的代数方程.并将此方法应用于封闭声腔响应的求解,即对亥姆霍兹方程进行离散,建立其最小二乘无网格配点格式.该方法的系数矩阵是对称正定的,因而保证了解的稳定性.通过数值算例分别验证了配点均匀分布与随机分布时此方法的...