The least‐squares meshfree method

A new efficient meshfree method is presented in which the first‐order least‐squares method is employed instead of the Galerkin's method. In the meshfree methods based on the Galerkin formulation, the source of many difficulties is in the numerical integration. The current method, in this respect, has different characteristics and is expected to remove some of the integration‐related problems. It is demonstrated through numerical examples that the present formulation is highly robust to integration errors. Therefore, numerical integration can be performed with great ease and effectiveness using very simple algorithms. Copyright © 2001 John Wiley & Sons, Ltd.     

A study on the convergence of least‐squares meshfree method under inaccurate integration

In the authors' previous work, it has been shown through numerical examples that the least‐squares meshfree method (LSMFM) is highly robust to the integration errors while the Galerkin meshfree method is very sensitive to them. A mathematical study on the convergence of the solution of LSMFM under inaccurate integration is presented. New measures are introduced to take into account the integration errors in the error estimates. It is shown that, in LSMFM, solution errors are bounded by approximation errors even when integration is not accurate. Copyright © 2003 John Wiley & Sons, Ltd.     

A posteriori error estimates and an adaptive scheme of least‐squares meshfree method

A posteriori error estimates and an adaptive refinement scheme of first‐order least‐squares meshfree method (LSMFM) are presented. The error indicators are readily computed from the residual. For an elliptic problem, the error indicators are further improved by applying the Aubin–Nitsche method. It is demonstrated, through numerical examples, that the error indicators coherently reflect the actual error. In the proposed refinement scheme, Voronoi cells are used for inserting new nodes at appropriate positions. Numerical examples show that the adaptive first‐order LSMFM, which combines the proposed error indicators and nodal refinement scheme, is effectively applied to the localized problems such as the shock formation in fluid dynamics. Copyright © 2003 John Wiley & Sons, Ltd.     

An improved moving‐least‐squares reconstruction for immersed boundary method

The current work presents an improved immersed boundary method based on the ideas proposed by Vanella and Balaras (M. Vanella, E. Balaras, A moving‐least‐squares reconstruction for embedded‐boundary formulations, J. Comput. Phys. 228 (2009) 6617–6628). In the method, an improved moving‐least‐squares approximation is employed to build the transfer functions between the Lagrangian points and discrete Eulerian grid points. The main advantage of the improved method is that there is no need to obtain the inverse matrix, which effectively eliminates numerical instabilities caused by matrix inversion and reduces the computational cost significantly. Several different flow problems (Taylor‐Green decaying vortices, flows past a stationary circular cylinder and a sphere, and the sedimentation of a free‐falling sphere in viscous fluid) are simulated to validate the accuracy and efficiency of the method proposed in the present paper. The simulation results show good agreement with previous numerical and experimental results, indicating that the improved immersed boundary method is efficient and reliable in dealing with the fluid–solid interaction problems. Copyright © 2015 John Wiley & Sons, Ltd.     

An adaptive domain decomposition coupled finite element–boundary element method for solving problems in elasto‐plasticity

The purpose of this paper is to present an adaptive finite element–boundary element method (FEM–BEM) coupling method that is valid for both two‐ and three‐dimensional elasto‐plastic analyses. The method takes care of the evolution of the elastic and plastic regions. It eliminates the cumbersome of a trial and error process in the identification of the FEM and BEM sub‐domains in the standard FEM–BEM coupling approaches. The method estimates the FEM and BEM sub‐domains and automatically generates/adapts the FEM and BEM meshes/sub‐domains, according to the state of computation. The results for two‐ and three‐dimensional applications in elasto‐plasticity show the practicality and the efficiency of the adaptive FEM–BEM coupling method. Copyright © 2009 John Wiley & Sons, Ltd.     

Least‐squares mixed finite elements for small strain elasto‐viscoplasticity

The main aim of this contribution is to provide a mixed finite element for small strain elasto‐viscoplastic material behavior based on the least‐squares method. The L₂‐norm minimization of the residuals of the given first‐order system of differential equations leads to a two‐field functional with displacements and stresses as process variables. For the continuous approximation of the stresses, lowest‐order Raviart–Thomas elements are used, whereas for the displacements, standard conforming elements are employed. It is shown that the non‐linear least‐squares functional provides an a posteriori error estimator, which establishes ellipticity of the proposed variational approach. Further on, details about the implementation of the least‐squares mixed finite elements are given and some numerical examples are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Moving least‐squares approximation with discontinuous derivative basis functions for shell structures with slope discontinuities

Moving least‐squares approximation with discontinuous derivative basis functions (MLSA‐DBF) is introduced for analysis of shell structures with slope discontinuities. To deal with shells with arbitrary slope discontinuities, the Cartesian coordinate is introduced in the construction of MLSA on the shell surface. The possible causes of singularity in the moment matrix of MLSA on the shell surface with slope discontinuities are identified, and the Moore–Penrose pseudoinverse is used to obtain the generalized inverse of the singular moment matrix resulting from linear dependency and insufficient influence nodes in the MLSA. Following the proposed formulations for shear deformable shell structures with slope discontinuities in the Cartesian coordinates, several numerical examples are analyzed to demonstrate the performance, validity, accuracy, and convergence properties of the proposed MLSA‐DBF approach. Copyright © 2007 John Wiley & Sons, Ltd.     

Complex variable moving least‐squares method: a meshless approximation technique

Based on the moving least‐squares (MLS) approximation, we propose a new approximation method—the complex variable moving least‐squares (CVMLS) approximation. With the CVMLS approximation, the trial function of a two‐dimensional problem is formed with a one‐dimensional basis function. The number of unknown coefficients in the trial function of the CVMLS approximation is less than in the trial function of the MLS approximation, and we can thus select fewer nodes in the meshless method that is formed from the CVMLS approximation than are required in the meshless method of the MLS approximation with no loss of precision. The meshless method that is derived from the CVMLS approximation also has a greater computational efficiency. From the CVMLS approximation, we propose a new meshless method for two‐dimensional elasticity problems—the complex variable meshless method (CVMM)—and the formulae of the CVMM for two‐dimensional elasticity problems are obtained. Compared with the conventional meshless method, the CVMM has a greater precision and computational efficiency. For the purposes of demonstration, some selected numerical examples are solved using the CVMM. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of strength difference in elasto‐plasticity for adhesive materials

Experimental evidence of certain adhesive materials reveals elastic strains, plastic strains and hardening. Furthermore, a pronounced strength difference effect between tension, torsion or combined loading is observed. For simulation of these phenomena, a yield function dependent on the first and second basic invariants of the related stress tensor in the framework of elasto‐plasticity is used in this work. A plastic potential with the same mathematical structure is introduced to formulate the evolution equation for the inelastic strains. Furthermore, thermodynamic consistency of the model equations is considered, thus rendering some restrictions on the material parameters. For evolution of the strain like internal variable, two cases are considered, and the consequences on the thermodynamic consistency and the numerical implementation are extensively discussed. The resulting evolution equations are integrated with an implicit Euler scheme. In particular, the reduction of the resulting local problem is performed, and for the finite‐element equilibrium iteration, the algorithmic tangent operator is derived. Two examples are presented. The first example demonstrates the capability of the model equations to simulate the yield strength difference between tension and torsion for the adhesive material Betamate 1496. A second example investigates the deformation evolution of a compact tension specimen with an adhesive zone. Copyright © 2005 John Wiley & Sons, Ltd.     

Floating stress‐point integration meshfree method for large deformation analysis of elastic and elastoplastic materials

A new meshfree formulation of stress‐point integration, called the floating stress‐point integration meshfree method, is proposed for the large deformation analysis of elastic and elastoplastic materials. This method is a Galerkin meshfree method with an updated Lagrangian procedure and a quasi‐implicit time‐advancing scheme without any background cell for domain integration. Its new formulation is based on incremental equilibrium equations derived from the incremental virtual work equation, which is not generally used in meshfree formulations. Hence, this technique allows the temporal continuity of the mechanical equilibrium to be naturally achieved. The details of the new formulation and several examples of the large deformation analysis of elastic and elastoplastic materials are presented to show the validity and accuracy of the proposed method in comparison with those of the finite element method. Copyright © 2012 John Wiley & Sons, Ltd.     

Optimization of thermo‐elasto‐plastic processes using Eulerian sensitivity analysis

A computational scheme for the analysis and optimization of quasi‐static thermo‐mechanical processes is presented in this paper. In order to obtain desirable mechanical transformations in a workpiece using a thermal treatment process, the optimal control parameters need to be determined. The problem is addressed by posing the process as a decoupled thermo‐mechanical finite element problem and performing an optimization using gradient methods. The forward problem is solved using the Eulerian formulation since it is computationally more efficient compared to an equivalent Lagrangian formulation. The design sensitivities required for the optimization are developed analytically using direct differentiation. This systematic design approach is applied to optimize a laser forming process. The objective is to maximize the angular distortion of a specimen subject to the constraint that the phase transition temperature is not exceeded at any point in the model. Copyright © 2003 John Wiley & Sons, Ltd.     

Mixed plate bending elements based on least‐squares formulation

A finite element formulation for the bending of thin and thick plates based on least‐squares variational principles is presented. Finite element models for both the classical plate theory and the first‐order shear deformation plate theory (also known as the Kirchhoff and Mindlin plate theories, respectively) are considered. High‐order nodal expansions are used to construct the discrete finite element model based on the least‐squares formulation. Exponentially fast decay of the least‐squares functional, which is constructed using the L₂ norms of the equations residuals, is verified for increasing order of the nodal expansions. Numerical examples for the bending of circular, rectangular and skew plates with various boundary conditions and plate thickness are presented to demonstrate the predictive capability and robustness of the new plate bending elements. Plate bending elements based on this formulation are shown to be insensitive to both shear‐locking and geometric distortions. Copyright © 2004 John Wiley & Sons, Ltd.     

A modified least‐squares mixed finite element with improved momentum balance

The main goal of this contribution is to provide an improved mixed finite element for quasi‐incompressible linear elasticity. Based on a classical least‐squares formulation, a modified weak form with displacements and stresses as process variables is derived. This weak form is the basis for a finite element with an advanced fulfillment of the momentum balance and therefore with a better performance. For the continuous approximation of stresses and displacements on the triangular and tetrahedral elements, lowest‐order Raviart–Thomas and linear standard Lagrange interpolations can be used. It is shown that coercivity and continuity of the resulting asymmetric bilinear form could be established with respect to appropriate norms. Further on, details about the implementation of the least‐squares mixed finite elements are given and some numerical examples are presented in order to demonstrate the performance of the proposed formulation. Copyright © 2009 John Wiley & Sons, Ltd.     

Space–time meshfree collocation method: Methodology and application to initial‐boundary value problems

A novel space–time meshfree collocation method (STMCM) for solving systems of non‐linear ordinary and partial differential equations by a consistent discretization in both space and time is proposed as an alternative to established mesh‐based methods. The STMCM belongs to the class of truly meshfree methods, i.e. the methods that do not have any underlying mesh, but work on a set of nodes only without any a priori node‐to‐node connectivity. Instead, the neighbouring information is established on‐the‐fly. The STMCM is constructed using the Interpolating Moving Least‐squares technique, which allows a simplified implementation of boundary conditions due to fulfillment of the Kronecker delta property by the kernel functions, which is not the case for the major part of other meshfree methods. The method is validated by several examples ranging from interpolation problems to the solution of PDEs, whereas the STMCM solutions are compared with either analytical or reference ones. Copyright © 2009 John Wiley & Sons, Ltd.     

A computational approach for multi‐scale analysis of heterogeneous elasto‐plastic media subjected to short duration loads

The asymptotic expansion homogenization (AEH) approach has found wide acceptance for the study of heterogeneous structures due to its ability to account for multi‐scale features. The emphasis of the present study is to develop consistent AEH numerical formulations to address elasto‐plastic material response of structures subjected to short‐duration transient loading. A second‐order accurate velocity‐based explicit time integration method, in conjunction with the AEH approach, is currently developed that accounts for large deformation non‐linear material response. The approach is verified under degenerate homogeneous conditions using existing experimental data in the literature and its ability to account for heterogeneous conditions is demonstrated for a number of test problems. Copyright © 2004 John Wiley & Sons, Ltd.     

The least‐squares finite element method in elasticity. Part II: Bending of thin plates

A least‐squares finite element method (LSFEM) for bending problems of thin plates is developed. This LSFEM is based on the first‐order deflection‐slope‐moment‐shear force formulation. Four compatibility conditions are added into the first‐order system; thus, the method can accommodate all kinds of equal‐order interpolations. Numerical experiments on various examples show that the method achieves an optimal rate of convergence for all eight variables. Copyright © 2002 John Wiley & Sons, Ltd.     

A continuum sensitivity method for finite thermo‐inelastic deformations with applications to the design of hot forming processes

A computational framework is presented to evaluate the shape as well as non‐shape (parameter) sensitivity of finite thermo‐inelastic deformations using the continuum sensitivity method (CSM). Weak sensitivity equations are developed for the large thermo‐mechanical deformation of hyperelastic thermo‐viscoplastic materials that are consistent with the kinematic, constitutive, contact and thermal analyses used in the solution of the direct deformation problem. The sensitivities are defined in a rigorous sense and the sensitivity analysis is performed in an infinite‐dimensional continuum framework. The effects of perturbation in the preform, die surface, or other process parameters are carefully considered in the CSM development for the computation of the die temperature sensitivity fields. The direct deformation and sensitivity deformation problems are solved using the finite element method. The results of the continuum sensitivity analysis are validated extensively by a comparison with those obtained by finite difference approximations (i.e. using the solution of a deformation problem with perturbed design variables). The effectiveness of the method is demonstrated with a number of applications in the design optimization of metal forming processes. Copyright © 2002 John Wiley & Sons, Ltd.     

A least‐squares coupling method between a finite element code and a discrete element code

This paper presents a coupling method between a discrete element code CeaMka3D and a finite element code Sem. The coupling is based on a least‐squares method, which adds terms of forces to finite element code and imposes the velocity at coupling particles. For each coupling face, a small linear system with a constant matrix is solved. This method remains conservative in energy and shows good results in applications. Copyright © 2014 John Wiley & Sons, Ltd.     

Augmented Lagrangian and Galerkin least‐squares methods for membrane contact

In this paper, we propose a stabilized finite element method for the numerical solution of contact between a small deformation elastic membrane and a rigid obstacle. We limit ourselves to friction‐free contact, but the formulation is readily extendable to more complex situations.     

A simple integration algorithm for a non‐associated anisotropic plasticity model for sheet metal forming

In this paper, an anisotropic material model based on a non‐associated flow rule and nonlinear mixed isotropic‐kinematic hardening is developed. The quadratic Hill48 yield criterion is considered in the non‐associated model for both yield function and plastic potential to account for anisotropic behavior. The developed model is integrated based on fully implicit backward Euler's method. The resulting problem is reduced to only two simple scalar equations. The consistent local tangent modulus is obtained by exact linearization of the algorithm. All numerical development was implemented into user‐defined material subroutine for the commercial finite element code ABAQUS/Standard. The performance of the present algorithm is demonstrated by numerical examples. Copyright © 2015 John Wiley & Sons, Ltd.