• Interaction and multiscale mechanics
• /
• v.2 no.2
• /
• pp.129-151
• /
• 2009

This paper presents a meshfree procedure using a convex generalized meshfree (GMF) approximation for the large deformation analysis of particle-reinforced rubber compounds on microscopic level. The convex GMF approximation possesses the weak-Kronecker-delta property that guarantees the continuity of displacement across the material interface in the rubber compounds. The convex approximation also ensures the positive mass in the discrete system and is less sensitive to the meshfree nodal support size and integration order effects. In this study, the convex approximation is generated in the GMF method by choosing the positive and monotonic increasing basis function. In order to impose the periodic boundary condition in the unit cell method for the microscopic analysis, a singular kernel is introduced on the periodic boundary nodes in the construction of GMF approximation. The periodic boundary condition is solved by the transformation method in both explicit and implicit analyses. To simulate the interface de-bonding phenomena in the rubber compound, the cohesive interface element method is employed in corporation with meshfree method in this study. Several numerical examples are presented to demonstrate the effectiveness of the proposed numerical procedure in the large deformation analysis.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.197-210
• /
• 2013

A fundamental trend of processor architecture evolving towards exaflops is fast increasing floating point performance (so-called "free" flops) accompanied by much slowly increasing memory and network bandwidth. In order to fully enjoy the "free" flops, a numerical algorithm of PDEs should request more flops per byte or increase arithmetic intensity. A meshfree/GFEM approximation can be the class of the algorithm. It is shown in a GFEM without extra dof that the kind of approximation takes advantages of the high performance of manycore GPUs by a high accuracy of approximation; the "expensive" method is found to be reversely hardware-efficient on the emerging architecture of manycore.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.137-156
• /
• 2013

In this paper, a meshfree shell adaptive procedure is developed for the applications in the sheet metal forming simulation. The meshfree shell formulation is based on the first-order shear deformable shell theory and utilizes the degenerated continuum and updated Lagrangian approach for the nonlinear analysis. For the sheet metal forming simulation, an h-type adaptivity based on the meshfree background cells is considered and a geometric error indicator is adopted. The enriched nodes in adaptivity are added to the centroids of the adaptive cells and their shape functions are computed using a first-order generalized meshfree (GMF) convex approximation. The GMF convex approximation provides a smooth and non-negative shape function that vanishes at the boundary, thus the enriched nodes have no influence outside the adapted cells and only the shape functions within the adaptive cells need to be re-computed. Based on this concept, a multi-level refinement procedure is developed which does not require the constraint equations to enforce the compatibility. With this approach the adaptive solution maintains the order of meshfree approximation with least computational cost. Two numerical examples are presented to demonstrate the performance of the proposed method in the adaptive shell analysis.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.107-125
• /
• 2013

A strain smoothing meshfree formulation with stabilized conforming nodal integration is presented for modeling the consolidation process in saturated porous media. In the present method, nodal strain smoothing is consistently introduced into the meshfree approximation of strain and pore pressure gradient variables associated with the saturated porous media. Meanwhile, in order to achieve a consistent numerical implementation, a smoothing approximation of the meshfree shape function within a nodal representative domain is also proposed in the stiffness construction. The resulting discrete system of equations is all expressed in smoothed nodal measures that are very efficient for numerical evaluation. Subsequently the space-time fully discrete equations are further established by the generalized trapezoidal rule for time integration. The effectiveness of the proposed meshfree consolidation analysis method is systematically illustrated by several benchmark problems.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2003.06a
• /
• pp.1569-1572
• /
• 2003

Meshfree approximations exhibit significant potential to solve partial differential equations. Meshfree methods have been successfully applied to various problems which the traditional finite element methods have difficulties to handle, including the quasi-static and dynamic fracture. large deformation problems, contact problems, and strain localization problems. A meshfree method based on the reproducing kernel particle approximation(RKPM) is applied to sheet metal forming analysis in this research. Metal forming examples, such as stretch forming and flanging operation, are analyzed to demonstrate the performance of the proposed meshfree method for largely deformed elasto-plastic material.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.83-105
• /
• 2013

In this paper, a multi-scale meshfree-enriched finite element formulation is presented for the analysis of acoustic wave propagation problem. The scale splitting in this formulation is based on the Variational Multi-scale (VMS) method. While the standard finite element polynomials are used to represent the coarse scales, the approximation of fine-scale solution is defined globally using the meshfree enrichments generated from the Generalized Meshfree (GMF) approximation. The resultant fine-scale approximations satisfy the homogenous Dirichlet boundary conditions and behave as the "global residual-free" bubbles for the enrichments in the oscillatory type of Helmholtz solutions. Numerical examples in one dimension and two dimensional cases are analyzed to demonstrate the accuracy of the present formulation and comparison is made to the analytical and two finite element solutions.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.365-372
• /
• 2006

A new meshfree formulation is developed for material discontinuity problems. A local interfacial jump function which is defined as hyperplane function is embedded in the meshless approximation and the approximation accurately models functions with jumps in the displacement and the derivative fields. Diffuse derivative technique copes with difficulty due to complexity of derivative computation of meshfree approximation. Collocation method with diffuse derivative accelerates computing speed for numerical solution. By solving inclusion and composite material problems, the robustness and effectiveness of the method are verified.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.583-590
• /
• 2002

Point collocation method based on the fast derivatives approximation of meshfree shape function is applied to solid mechanics in this study. Enhanced meshfree approximation with approximated derivative of shape function is reviewed, and formulation of linear elastic solid mechanics by point collocation method is presented. It implies that governing equation of solid mechanics with strong form is directly formulated without no numerical integration cells or grid. The regularity of weight function is not required due to a use of approximated derivative, so we propose the exponential type weight function that is discontinuous in first derivative. The convergence and stability of the proposed method is verified by passing the generalized patch test. Also, the efficiency and applicability of the proposed method in solid mechanics is verified by solving types of solid problems. Numerical results show that not only a use of proposed weight function leads lower error and higher convergence rate than that of the conventional weight functions, but also the improved collocation method with derivative approximation enables to compute the derivatives of shape function very fast and accurately enough to replace the classical direct derivative calculation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.495-500
• /
• 2007

A Meshfree is a method used to establish algebraic equations of system for the whole problem domain without the use of a predefined mesh for the domain discretization. A point interpolation method is based on combining radial and polynomial basis functions. Involvement of radial basis functions overcomes possible singularity. Furthermore, the interpolation function passes through all scattered points in an influence domain and thus shape functions are of delta function property. This makes the implementation of essential boundary conditions much easier than the meshfree methods based on the moving least-squares approximation. This study aims to investigate a stress analysis of structural element between a meshfree method and the finite element method. Examples on cantilever type plate and stress concentration problems show that the accuracy and convergence rate of the meshfree methods are high.