无网格方法在爆轰数值模拟中的应用

以SPH算法为代表的无网格方法在爆轰波的数值模拟中具有明显的优势,采用SPH算法模拟高能炸药水下爆炸爆轰过程,得到了压力、速度等时历曲线.将数值模拟结果与理论和实验方法得到的结果进行了对比分析,表明SPH算法非常适宜处理高能炸药水下爆炸的极短瞬时具有大变形和高度非均匀的动力学极端情形,且求解结果已达到了较高的精度。     

ICM Method Combined with Meshfree Approximation for Continuum Structure

The independent continuous mapping(ICM) method is integrated into element free Galerkin method and a new implementation of topology optimization for continuum structure is presented.To facilitate the enforcement of the essential boundary condition and derivative of various sensitivities,a singular weight function in element free Galerkin method is introduced.Material point variable is defined to illustrate the condition of material point and its vicinity instead of element or node.The topological variables ...     

An adaptive isogeometric analysis meshfree collocation method for elasticity and frictional contact problems

A collocation method has been recently developed as a powerful alternative to Galerkin's method in the context of isogeometric analysis, characterized by significantly reduced computational cost, but still guaranteeing higher-order convergence rates. In this work, we propose a novel adaptive isogeometric analysis meshfree collocation (IGAM-C) for the two-dimensional (2D) elasticity and frictional contact problems. The concept of the IGAM-C method is based upon the correspondence between the isogeometric collocation and reproducing kernel meshfree approach, which facilitates the robust mesh adaptivity in isogeometric collocation. The proposed method reconciles collocation at the Greville points via the discretization of the strong form of the equilibrium equations. The adaptive refinement in collocation is guided by the gradient-based error estimator. Moreover, the resolution of the nonlinear equations governing the contact problem is derived from a strong form to avoid the disadvantages of numerical integration. Numerical examples are presented to demonstrate the effectiveness, robustness, and straightforward implementation of the present method for adaptive analysis.     

Thin shell analysis from scattered points with maximum‐entropy approximants

We present a method to process embedded smooth manifolds using sets of points alone. This method avoids any global parameterization and hence is applicable to surfaces of any genus. It combines three ingredients: (1) the automatic detection of the local geometric structure of the manifold by statistical learning methods; (2) the local parameterization of the surface using smooth meshfree (here maximum‐entropy) approximants; and (3) patching together the local representations by means of a partition of unity. Mesh‐based methods can deal with surfaces of complex topology, since they rely on the element‐level parameterizations, but cannot handle high‐dimensional manifolds, whereas previous meshfree methods for thin shells consider a global parametric domain, which seriously limits the kinds of surfaces that can be treated. We present the implementation of the method in the context of Kirchhoff–Love shells, but it is applicable to other calculations on manifolds in any dimension. With the smooth approximants, this fourth‐order partial differential equation is treated directly. We show the good performance of the method on the basis of the classical obstacle course. Additional calculations exemplify the flexibility of the proposed approach in treating surfaces of complex topology and geometry. Copyright © 2010 John Wiley & Sons, Ltd.     

对称层合折板结构自由振动分析的移动最小二乘无网格法

提出了一种求解对称层合折板结构自由振动问题的移动最小二乘无网格法。以作者提出的折板无网格模型为基础,将对称层合折板结构视为由不同平面上对称层合板组成的复合结构。先基于一阶剪切变形理论,由移动最小二乘近似推导出各对称层合板的刚度和质量矩阵,再利用板与板间的位移协调条件,将各板的刚度和质量矩阵叠加得到整个结构的刚度和质量矩阵,推导出描述层合折板结构自由振动行为的控制方程。文末算例表明由本文方法得到的解与采用壳单元得到的ANSYS有限元解吻合良好,证明了本文方法的准确性。     

A temperature‐related homogenization technique and its implementation in the meshfree particle method for nanoscale simulations

A new homogenization technique, the temperature‐related Cauchy–Born (TCB) rule, is proposed in this paper with the consideration of the free energy instead of the potential energy. Therefore, temperature effects at the nanoscale can be investigated using continuum approximation with the implementation of the TCB rule. The TCB rule is verified via stress analyses of several crystalline solids. Temperature‐related material instability is also studied. In addition, a new hierarchical multiscale method is developed through implementing the TCB rule into meshfree particle methods. Quasicontinuum meshfree particle simulations are conducted to investigate bending of nanobeams, crack propagation in nanoplates and a three‐dimensional nanoindentation problem. Copyright © 2006 John Wiley & Sons, Ltd.     

Adaptive multiresolution refinement with distance fields

This paper describes a multiresolution approach to field modeling that can be used with any meshfree or mesh‐based method for adaptive solution refinement. The refined solution is represented as a superposition of a coarse (unrefined) solution and a sequence of refinements that provide additional degrees of freedom with higher spatial or functional resolution. Each refinement is treated as a solution to a boundary‐value problem within a specified refinement window. The proposed approach is based on the meshfree method with distance fields (Comput. Mech. 2000; 25 :305–316, Eng. Comput. 2002; 18 (4):295–311) and guarantees C^m continuity of the refined solutions with matching or non‐matching grids. The method does not restrict the shape of the refinement window and does not place any constraints on the type of basis functions, or relative position and resolution of the refinement grids. Combining the proposed approach with hierarchical space decompositions and a posteriori error estimators results in an effective tool for automatic solution refinement. Carefully chosen numerical examples illustrate the power and advantages of the proposed approach. Copyright © 2007 John Wiley & Sons, Ltd.     

Truncation error in mesh‐free particle methods

A truncation error analysis has been developed for the approximation of spatial derivatives in smoothed particle hydrodynamics (SPH) and related first‐order consistent methods such as the first‐order form of the reproducing kernel particle method. Error is shown to depend on both the smoothing length h and the ratio of particle spacing to smoothing length, Δx/h. For uniformly spaced particles in one dimension, analysis shows that as h is reduced while maintaining constant Δx/h, error decays as h² until a limiting discretization error is reached, which is independent of h. If Δx/h is reduced while maintaining constant h (i.e. if the number of neighbours per particle is increased), error decreases at a rate which depends on the kernel function's smoothness. When particles are distributed non‐uniformly, error can grow as h is reduced with constant Δx/h. First‐order consistent methods are shown to remove this divergent behaviour. Numerical experiments confirm the theoretical analysis for one dimension, and indicate that the main results are also true in three dimensions. This investigation highlights the complexity of error behaviour in SPH, and shows that the roles of both h and Δx/h must be considered when choosing particle distributions and smoothing lengths. Copyright © 2005 John Wiley & Sons, Ltd.     

Contact algorithms for the material point method in impact and penetration simulation

The inherent no‐slip contact constraint in the standard material point method (MPM) creates a greater penetration resistance. Therefore, the standard MPM was not able to treat the problems involving impact and penetration very well. To overcome these deficiencies, two contact methods for MPM are presented and implemented in our 3D explicit MPM code, MPM3D. In MPM, the impenetrability condition may not satisfied on the redefined regular grid at the beginning of each time step, even if it has been imposed on the deformed grid at the end of last time step. The impenetrability condition between bodies is only imposed on the deformed grid in the first contact method, while it is imposed both on the deformed grid and redefined regular grid in the second contact method. Furthermore, three methods are proposed for impact and penetration simulation to determine the surface normal vectors that satisfy the collinearity conditions at the contact surface. The contact algorithms are verified by modeling the collision of two elastic rings and sphere rolling problems, and then applied to the simulation of penetration of steel ball and perforation of thick plate with a particle failure model. In the simulation of elastic ring collision, the first contact algorithm introduces significant disturbance into the total energy, but the second contact algorithm can obtain the stable solution by using much larger time step. It seems that both contact algorithms give good results for other problems, such as the sphere rolling and the projectile penetration. Copyright © 2010 John Wiley & Sons, Ltd.     

A weighted least squares particle‐in‐cell method for solid mechanics

A novel meshfree method is proposed that incorporates features of the material point (MPM) and generalized interpolation material point (GIMP) methods and can be used within an existing MPM/GIMP implementation. Weighted least squares kernel functions are centered at stationary grid nodes and used to approximate field values and gradients. Integration is performed over cells of the background grid and material boundaries are approximated with an implicit surface. The proposed method avoids nearest‐neighbor searches while significantly improving accuracy over MPM and GIMP. Implementation is discussed in detail and several example problems are solved, including one manufactured solution which allows measurement of dynamic, non‐linear, large deformation performance. Advantages and disadvantages of the method are discussed. Copyright © 2010 John Wiley & Sons, Ltd.