Acessibilidade / Reportar erro

Imposition of Dirichlet Boundary Conditions in Element Free Galerkin Method through an Object-Oriented Implementation

Abstract

One of the main drawbacks of Element Free Galerkin (EFG) method is its dependence on moving least square shape functions which don’t satisfy the Kronecker Delta property, so in this method it’s not possible to apply Dirichlet boundary conditions directly. The aim of the present paper is to discuss different aspects of three widely used methods of applying Dirichlet boundary conditions in EFG method, called Lagrange multipliers, penalty method, and coupling with finite element method. Numerical simulations are presented to compare the results of these methods form the perspective of accuracy, convergence and computational expense. These methods have been implemented in an object oriented programing environment, called INSANE, and the results are presented and compared with the analytical solutions.

Keywords:
Element Free Galerkin method; EFG-FE Coupling; Lagrange Multipliers; Penalty Method; Object-Oriented Programming

1 INTRODUCTION

Conventional numerical methods like Finite Element Method (FEM) are very powerful in solving different engineering problems, and they are efficiently being widely used through commercial software packages. But due to their reliance on mesh, these methods show limitations or difficulties when dealing with problems having large deformations or arbitrary path discontinuities.

Meshfree methods were introduced to eliminate part of those difficulties such as distorted elements, need for re-meshing and the other limitations which arise because of element connectivity. A large variety of meshfree methods have been developed up to now for solving different engineering problems, such as smooth particle hydrodynamics (SPH) by (Lucy, 1997Lucy L. B (1997) A numerical approach to the testing of fission hypothesis. The Astronomical Journal 82: 1013-1024.), diffuse elements, (Nayroles et al., 1992Nayroles B, Touzot G, Villon P. (1992) Generalizing the finite element method: diffuse approximation and diffuse elements. Computational Mechanics 10: 307-318.), Element-Free Galerkin method (EFG), (Belytschko et al., 1994Belytschko T, Lu Y. Y, Gu L. (1994) Element-free Galerkin methods.. Int J. Numerical Methods in Engineering 37: 229-256.), reproducing kernel particle method (RKPM), (Liu et al., 1995Liu W. K, lun S, Li S, (1995) Reproducing kernel particle methods for structural dynamics. Int. J. Numerical Methods in Engineering 38: 1655-1679.), and hp-clouds methods, (Duarte and Oden, 1995Duarte C. A, Oden J. T. (1995) Hp clouds - a meshless method to solve boundary-value problems. Technical Report, Texas Institute for Computational and Applied Mathematics, University of Texas at Austin, 95-105.).

Among introduced meshfree methods, element-free Galerkin is one of the most common methods which is based on the Moving Least Square (MLS) approximation, (Lancaster et al., 1981Lancaster P, Salkauskas K, (1981) Surface generated by moving least square methods. J. Computational Mathematics 37: 141-58.). Since being introduced by (Belytschko et al., 1994Belytschko T, Lu Y. Y, Gu L. (1994) Element-free Galerkin methods.. Int J. Numerical Methods in Engineering 37: 229-256.), this method has been used for solving different engineering problems like elasticity, static and dynamic fracture, (Fleming et al., 1997Fleming M, Chu YA, Moran B, (1997) Enriched element-free Galerkin methods for crack-tip fields. Int. J. Numerical Methods in Engineering 40: 1483-504.), (Belytschko et al., 1995aBelytschko T, Lu YY, Gu L. (1995a) Element-free Galerkin methods for static and dynamic fracture, Int. J. of Solids and Structures 32: 2547-70, 1196Belytschko T., Tabbara M. (1996) Dynamic fracture using element-free Galerkin methods. Int. J. Numerical Methods in Engineering 39: 923-38., 1997), (Guiamatsia et al., 2009Guiamatsia I, Falzon B.G, Davies G.A.O, (2009) Element-Free Galerkin modeling of composite damage. Compos. Sci. Technol. 69: 2640-648.), manufacturing processes, (Yonghui et al., 2008Yonghui L., Jun C., Song Y. (2008) Numerical simulation of three-dimensional bulk forming processes by the element-free Galerkin method. Int. J. Advanced Manufacturing Technology 36(5): 442-450.), and heat transfer, (Singh et al., 2002Singh I.V, Sandeep K, Prakash R. (2002) The element free Galerkin method in three dimensional steady state heat conduction. Int. J. Computational Engineering Science 3(3): 291-303.).

The EFG method is a very promising method for the treatment of partial differential equations. Because of the absence of element connectivity, node particles can be added without the need for re-meshing, so adaptive refinement of the discretization can be performed easily. This makes the EFG method more flexible than finite element method, and prominent for solving problems involving large deformation, fracture and fragmentation. But EFG method still has some drawbacks compared with FE method. The most difficulties of this method are: 1) as a meshfree method, this method is also computationally expensive because of the complex meshfree interpolations and the complex implementation algorithm, 2) since this method is based on the moving least squares shape functions, it is difficult to impose Dirichlet boundary conditions. Considerable researches have been devoted to overcome these difficulties, such as Lagrange multipliers method, (Belytschko et al., 1994Belytschko T, Lu Y. Y, Gu L. (1994) Element-free Galerkin methods.. Int J. Numerical Methods in Engineering 37: 229-256.), and the penalty method (Zhu and Atluri, 1998Zhu T., Atluri S.N. (1998) A modified collocation method and a penalty formulation for enforcing the essential boundary conditions in the element free Galerkin method. Computational Mechanics 21(3): 211-222.), based on a modification of weak form. Other methods were proposed that can be interpreted as a modification of the interpolation shape functions, see for instance, (Belytschko et al., 1995bBelytschko T, Organ D, Krongauz Y. (1995b) A coupled finite element-element-free Galerkin method. Computational Mechanics 17: 186-195.), (Huerta et al., 2000Huerta A, Fernandez-Mendez S. (2000) Enrichment and coupling of the finite element and meshless method. Int. J. Numerical Methods in Engineering 48: 1615 1636.). The Lagrange multipliers method, firstly introduced by (Babuška, 1973Babuška I. (1973) The finite element method with Lagrange multipliers. J. Numerical Mathematics 20: 179-192.), has been thoroughly studied by several researchers in (Pitkäranta, 1979Pitkäranta J. (1979) Boundary subspaces for the finite element method of Lagrange multipliers. J. Numerical Mathematics 33: 273-289.), (Barbosa et al., 1992Barbosa H.J.C, Hughes T.J.R. (1992) Boundary Lagrange multipliers in finite element methods: error analysis in natural norms. J. Numerical Mathematics 62: 1-15.), (Stenberg, 1995Stenberg R. (1995) On some techniques for approximating boundary conditions in the finite element methods. J. Computer in Applied Mathematics 63: 139-148.). This method is used to approximate the Dirichlet boundary conditions by using some unknown parameters called Lagrange multipliers. It was proven that this method had optimal convergence rates for the finite element spaces. The penalty method requires only the choice of the penalty parameter. The large values of this parameter must be used in order to impose the Dirichlet boundary condition in a proper manner. In practice, that might lead to ill-conditioned systems of equations, reducing the applicability of this method. Further information on penalty method can be found in (Zhu and Atluri, 1998Zhu T., Atluri S.N. (1998) A modified collocation method and a penalty formulation for enforcing the essential boundary conditions in the element free Galerkin method. Computational Mechanics 21(3): 211-222.).

Coupling EFG with finite element in order to apply Dirichlet boundary conditions is one the methods that can be interpreted as a modification of the interpolation shape functions. In order to do so, several coupling techniques are proposed in different studies. (Belytschko et al., 1995bBelytschko T, Organ D, Krongauz Y. (1995b) A coupled finite element-element-free Galerkin method. Computational Mechanics 17: 186-195.), (Krongauz and Belytschko, 1996Krongauz Y, Belytschko T. (1996) Enforcement of essential boundary conditions in meshless approximations using finite elements. Computer Methods in Applied Mechanics and Engineering 131: 133-145.) developed an algorithm for coupling EFG and FEM by a mixed interpolation in the transition domain, where FE nodes are substituted by particles and connected via ramp functions to the EFG nodes. The drawback of this method is that the derivatives are discontinuous along the interface of EFG and FE domains. (Huerta et al., 2000Huerta A, Fernandez-Mendez S. (2000) Enrichment and coupling of the finite element and meshless method. Int. J. Numerical Methods in Engineering 48: 1615 1636.) developed a hierarchical approximation based on finite element and mesh-free methods, which is able to remove the discontinuities in the derivatives of shape functions across the inner boundaries that was the major drawback of the method proposed by Belytschko[14]. Another approach for coupling EFG and FE is by using Lagrange multipliers. In this approach the substitution of FE nodes by particles is not necessary, (Hegen, 1996Hegen D. (1996) Element free Galerkin methods in combination with finite element approaches. Computer Methods in Applied Mechanics and Engineering 135: 143-166.), (Rabczuk and Belytschko, 2006Rabczuk T, Belytschko T. (2006) Application of mesh free methods to static fracture of reinforced concrete structures. Int. J. Fracture 137(1-4): 19-49.). Another coupling approach using Lagrange multipliers was proposed by Belytschko and Xiao, called ‘bridging domain coupling method, (Belytschko and Xiao, 2004). In their approach, FE nodes and meshfree particles could overlap. (Gu and Zhang, 2008Gu Y.T, Zhang L.C. (2008) Coupling of the meshfree and finite element methods for determination of the crack tip fields. Engineering Fracture Mechanics 75: 986-1004.) developed a new concurrent simulation technique to couple the EFG method with the finite element method for the analysis of crack tip fields, using the ‘bridging domain coupling method’ concept and Lagrange multipliers. (Liu et al., 2008Liu L.C, Dong X.H, Li C.X. (2008) An Adaptive EFG-FE Coupling Method for the Numerical Simulation Of Extrusion Processes. Acta Metallurgica Sinica (English Letters), 21: 380-388.) presented an adaptive EFG-FE coupling method oriented for bulk metal forming. In their method the simulation is performed using finite elements initially. During the process, the FE region with severe deformation occurrence is automatically converted to an EFG region.

The application of FEM object-oriented programming has been receiving great attention over the last two decades. As a result, a bunch of FEM-based codes used object-oriented concept as their implementation strategy, such as: a C++ based finite element code developed by (Patzak et al., 2001Patzaak B., Bittnar Z., (2001) Design of object oriented finite element code. Advances in Engineering Software 32: 759-767.), a meshless Object-oriented programming (OOP) code using non-uniform rational b-spline, (Zhang et al., 2006Zhang X., Subbarayan G, (2006) jNURBS: An object-oriented, symbolic framework for integrated, meshless analysis and optimal design. Advances in Engineering Software 37: 287-311), Freefem++, (Hecht, 2012Hecht F., (2012) New development in freefem++. Journal of Numerical Mathematics 20(3-4): 251-265.), a fem code to solve two- and three-dimensional problems, a general framework for weak-form meshfree methods by (Hsieh et al., 2014Hsieh Y-M, Pan M-S (2014) ESFM: An essential software framework for meshfree methods. Advances in Engineering Software 76: 133-147.), and an educational FE software for truss structures presented by (Zou et al., 2014Zou W., Li X., Guo G. (2014) EFESTS, Educational finite element software for truss structure. Part II: Linear static analysis. International Journal of Mechanical Engineering Education 42(4): 307-319.).

In the present work, the three introduced methods of applying Dirichlet boundary conditions are investigated through developing MATLAB codes and also implementation of an object-oriented programming environment, called INSANE (INteractive Structural ANalysis Environment- http://www.insane.dees.ufmg.br) platform, (Alves et al., 2013Alves P.D, Barros F.B, Pitangueira R.L.S (2013) An object oriented approach to the generalized finite element method. Advances in Engineering Software, 59: 1-18.), (Fonseca and Pitangueira, 2007Fonseca F.T, Pitangueira R.L.S (2007) An object oriented class organization for dynamic geometrically non-linear. CMNE (Congress on Numerical Methods in Engineering)/CILAMCE 2007, Porto, 13-15 June, Portugal). The main contribution of the present work is regarding to implementation which, to the best knowledge of the author, is the first object oriented program for meshfree methods and also EFG/FE coupling approach. The implemented EFG program is highly efficient with a wide range of applications since it' s built on the basis of a well-developed finite element program, not only in the sense of using integration cells but, indeed the ability of using several features of FEM code which are already implemented, or will be further implemented in INSANE computational platform. The most significant example for this is the capability of using existing physically nonlinear models of FEM code, such as smeared crack and damage models in order to solve continuum damage problems in meshfree framework. This paper is arranged as follows: first a brief review of EFG formulation will be presented, and the methods of applying Dirichlet boundary conditions which are studied in this work will be explained. Then the INSANE platform will be briefly introduced with emphasis on its meshfree modules. In Sect. 5 the numerical examples will be presented and the results will be compared with the analytical solutions.

2 ELEMENT FREE GALERKIN METHOD

2.1 Moving Least Square Approximation

The EFG method utilizes the moving least square approximants, which are constructed in terms of nodes only:

(1a)

where p(x) is a complete polynomial basis of arbitrary order and ai (x) are coefficients which are obtained at any point x by minimizing the following:

(1b)

in which nsn is the number of nodes in the support domain of point x, and w(x-xI) is considered as a cubic spline weight function defined by:

(2)

Minimizing Eq. (1b) leads to:

(3)

where:

(4)

and,

(5)

(6)

By substituting Eq. (3) into Eq. (1a), the MLS approximants are defined as:

(7)

where ∅I (x) are the EFG shape functions and are defined as:

(8)

2.2 EFG with Lagrange Multipliers

For any trial function, uS, and test function, vU, the Galerkin weak form of the EFG method with Lagrange multipliers is:

(9)

with the following boundary conditions:

u - g = 0 on Γg

-Ku.n = h on Γh

where U and S are the spaces of test and trial functions, respectively, Γg and Γh are the Dirichlet and Neumann boundaries, λ is Lagrange multiplier, and f is the excitation vector.

Applying the moving least square approximations (Eq. (7)), the resulting system of equations for EFG method with Lagrange multipliers is:

(10)

where:

(11)

As can be seen in Eq. (10), this method introduces additional unknown variables, λk , which results in a larger system of equations compared with FEM, which means increasing the computational cost.

2.3 EFG with Penalty Method

The Galerkin weak form for EFG with penalty method is:

(12)

in which αis the penalty parameter used to apply Dirichlet boundary condition. Using MLS approximation leads to the following system of equations:

(13)

where:

(14)

This method produces a system of equations with the same dimensions of FEM, but choosing the proper value for the penalty parameter is critical for this method. Penalty parameter should be as large as possible, but not too large to produce numerical problems due to increasing the condition number of the system matrix.

2.4 EFG Coupled with FEM

In order to couple EFG with FEM using ramp function, (Belytschko et al., 1995bBelytschko T, Organ D, Krongauz Y. (1995b) A coupled finite element-element-free Galerkin method. Computational Mechanics 17: 186-195.), domain of the problem is divided into three subdomains, as shown in Figure 1.

Figure 1
Domain discretization for EFG-FE coupled method.

In ΩEFG , the solution value at each point is approximated using the MLS approximants as in Eq. (7). In ΩFE , the standard finite element interpolation functions are used to approximate the solution, and in ΩIE which the interface elements are located, the following expression is used:

(15)

where uFE (x) and uEFG (x) are approximations for u in Ωi given by the FE and EFG approximations, respectively, and R(x) is the ramp function defined by (Belytschko et al., 1995bBelytschko T, Organ D, Krongauz Y. (1995b) A coupled finite element-element-free Galerkin method. Computational Mechanics 17: 186-195.):

(16)

In other words, the ramp function is equal to the sum of the FE shape functions of interface element nodes that are on the EFG boundary. By substituting the FE and EFG solution approximations in Eq. (15) it obtains:

(17)

which the ÑI are the shape functions of nodes of interface elements and are defined by [14]:

(18)

And the derivatives of the shape functions are:

(19)

The coupled shape functions and their derivatives for a typical one dimensional and two dimensional cases are shown in Figure 2 and Figure 3 respectively. It can be clearly seen that the derivatives of shape functions are discontinues on the borders of interface element.

Figure 2
Coupled Shape functions in 1-D: (a) Shape function, (b) Shape function derivative.

Figure 3
Coupled shape functions in 2-D: (a) Shape function, (b) Shape function derivative with respect to x, (c) Shape function derivative with respect to y.

3 INTERACTIVE STRUCTURAL ANALYSIS ENVIRONMENT (INSANE)

The INSANE environment is an open source software implemented in Java, an OOP language. It is composed of three main applications: pre-processor, processor and postprocessor. The pre- and post-processor are interactive graphical applications that provide tools to build representations of the structural problem as well as its results visualization. The processor is the numerical core of the system, responsible for obtaining the results of the analysis. There are different computational methods available for the processor, such as FEM, (Fonseca and Pitangueira, 2007Fonseca F.T, Pitangueira R.L.S (2007) An object oriented class organization for dynamic geometrically non-linear. CMNE (Congress on Numerical Methods in Engineering)/CILAMCE 2007, Porto, 13-15 June, Portugal), Generalized/Extended FEM, (Alves et al., 2013Alves P.D, Barros F.B, Pitangueira R.L.S (2013) An object oriented approach to the generalized finite element method. Advances in Engineering Software, 59: 1-18.), (Malekan and Barros, 2016Malekan M, Barros FB (2016) Well-conditioning global-local analysis using stable generalized/extended finite element method for linear elastic fracture mechanics. Computational Mechanics 58(5): 819-831.), (Malekan et al., 2016aMalekan M, Barros F.B, Pitangueira R.L.S, Alves P.D (2016a) An object-oriented class organization for global-local generalized finite element method. Latin American Journal of Solids and Structures 13(13):2529-2551., 2016bMalekan M, Barros F.B, Silva R.P, (2016b) Development and implementation of a well-conditioning approach toward generalized/extended finite element method into an object-oriented platform. Revista Interdisciplinar de Pesquisa em Engenharia 2(14):1-20., 2017Malekan M, Barros F.B, Pitangueira R.L.S, Alves P.D, Penna S.S (2017) A computational framework for a two-scale generalized/extended finite element method: generic imposition of boundary conditions. Engineering Computations, 34(3), pp 1-32.), and Boundary element method, (Anacleto et al., 2013Anacleto F.E.S, Ribeiro T.S.A, Ribeiro G.O, Pitangueira R.L.S, Penna S.S. (2013) An object-oriented tridimensional self-regular boundary element method implementation, Engineering Analysis with Boundary Elements 37(10):1276-1284), (Peixoto et al., 2016Peixoto R.G, Anacleto F.E.S, Ribeiro G.O, Pitangueira R.L.S (2016) A solution strategy for non-linear implicit BEM formulation using a unified constitutive modelling framework, Engineering Analysis with Boundary Elements, 64:294-310.).

3.1 Main Structure

The numerical core of INSANE is composed of two interfaces: Assembler, and Persistence and two abstract classes: Model and Solution. Figure 4 shows the Unified Modeling Language (UML) form of the numerical core.

Figure 4
Organization of numerical core of INSANE.

The Assembler interface builds the matrix of system of equations given by the discretization of a boundary or initial value problem, being implemented with following general representation:

(20)

Where X is the solution vector; A, B and C are matrices and D is a vector, where may or may not be dependent on the state variable or its derivatives.

Solution class has the necessary methods for solving the matrix system of Eq. (20), and the Model abstract class contains the data of the discrete model and provides Assembler with all necessary information to assemble Eq. (20).

Both Model and Solution communicate with the Persistence interface, which treats the input data and persists the output data to the other applications, any time it observes a modification of the discrete model state. This observation strategy occurs according to the Observer-Observable design pattern, (Horstmann et al., 2008Horstmann CS, Cornell G. (2008) Core java - fundamentals, vol. 1, 8thed. Sun Microsystems Press.). In this mechanism when an object of type observer (which implements the interface java.util.observer) is created, it will be added to a list of observers and will have a list of observables (which extends the class java.util.observable). In case of any modification in the state of any observable object, the change propagation mechanism is triggered, and the observer objects are notified to update themselves. This process guarantees the consistency between the observer and the observed components. In INSANE, the observer component is Persistence and the observed components are Solution and Model abstract classes.

3.2 Assembler

The MfreeAssembler class has the necessary methods for assembling the matrix system of the Eq. (20), with three ways to apply essential boundary conditions: coupling with finite element method, using the method presented in (Fernández-Méndez et al., 2004Fernández-Méndez S, Huerta A. (2004) Imposing essential boundary conditions in mesh-free methods. Computer Methods in Applied Mechanics and Engineering 193: 1257-1275.), Lagrange multipliers, and penalty method. In static analysis using penalty method or EFG-FE coupling technique, the matrix system of equations is:

(21)

where the matrix C represents the stiffness matrix of the model, X is the vector of nodal parameters, and D is the vector of forces. The indices u and p indicate unknown and prescribed values, respectively.

In static analysis using Lagrange multipliers the system of equations to be solved is:

(22)

This class has methods to build different parts of the matrix system, according to Eq. (21) and Eq. (22). The organization of the classes which implement the Assembler interface is illustrated in Figure 5. In this figure the UML diagram of the MfreeAssembler is showing the main methods of this class.

Figure 5
Organization of Assembler interface.

3.3 Model

The objective of the Model abstract class is to represent the discrete model of the problem. It is implemented by the object of class FemModel which itself consists of objects of different classes specifying the attributes of the model. Figure 6 shows the UML diagram of the MfreeModel class and its related classes. MfreeModel uses objects of Element class as EFG integration cells.

Figure 6
UML diagram of MfreeModel class.

3.4 Solution

The abstract class Solution, manages the solution process of solving matrix system of Eq. (20). In case of linear static problem a class named SteadyState implements the class Solution.

Each object of type SteadyState has attributes of two other objects of type Assembler and LinearEquationSystems .The class LinearEquationSystems has different methods for standard techniques of solving linear equation systems. This class receives the matrix system of equations (Eq. (20)) built by class Assembler, and computes the solution vector. For solving a linear static meshfree problem this class returns the vector of nodal values.

3.5 Analysis Flowchart

Figure 7 shows a detailed flowchart for meshfree analysis described in this section. The process starts with creating a meshfree input file, then defining analysis parameters such as size and shape of the influence domain. Two types of solution procedures are available, either linear elastic or nonlinear analysis. Following selection of the solver type, the problem will be solved and user can work on the results using the available post-processor part. Similar flowchart can be drawn for EFG/FE coupling approach, except the coupling domain have to be defined before solving the model.

Figure 7
Flowchart for Meshfree analysis.

4 NUMERICAL RESULTS

In order to investigate the difference between the methods of applying Dirichlet boundary conditions in EFG method, two well-known cases are studied: Timoshenko beam and Laplace equation. A computational code is developed in order to solve these problems with EFG meshfree method, using the three previously mentioned methods to apply Dirichlet boundary conditions. The case of Timoshenko beam is also implemented in INSANE platform.

4.1 Cantilever Beam

For a Timeshenko beam as an elastostatic problem, the differential equation is as follows:

(23)

In which b is the body force, σ =Dε, and ε =u, Γu and Γt are the essential and traction boundaries, respectively. The resulting system of linear equations is:

(24)

in which:

(25)

(26)

The model of the cantilever beam is illustrated in Figure 8. The following parameters are considered for the model: L=48, D=12, E=30e6, υ=0.3, and P=-1000.

Figure 8
Timoshenko beam.

Regular meshes were used, with 4 × 4 quadrature in each EFG cell and also for the coupled FE elements. A linear basis was used for EFG and standard Q4 elements were used in ΩFE .

Two cases are investigated for coupled EFG-FE problem. In the first case, the left half of the beam is modeled using FE elements, and the right half by EFG nodes. In the last strip of the FE domain the interface elements are located. The results for the deformation of the beam, in the y direction are illustrated in Figure 9 and are compared with Lagrange multipliers and analytical solution. The exact solution for the vertical displacement of Timoshenko beam is:

(27)

where I is the moment of inertia for a beam with rectangular cross-section and unit thickness.

Figure 9
Deflection of Timoshenko beam in y direction.

In the second case only a strip of the interface elements is located in the position of Γu and the remaining parts of the beam are modeled by EFG nodes. This strip of interface elements makes it possible to impose the essential boundary condition accurately.

The L2 norm of error for the displacement in the y direction is obtained from the Eq. (28) and compared for the two cases, and with the FE results in Figure 10. In this figure h is the distance between the nodes.

(28)

Figure 10
Variations of error versus mesh size, for 36, 64, and 240 nodes.

As can be seen in this figure, the errors of the coupled method are smaller than that of finite element. Moreover, in case 2, which only a strip of interface elements are inserted at the essential boundary, the coupled method has the best accuracy.

The Central Processing Unit (CPU) times of the solution for these two cases are compared with each other and with FE in Figure 11. As can be expected, case 2 involves the largest CPU time (because of the larger number of meshfree particles), but in cases of coarser nodal distributions, it has resulted in better accuracy, using almost the same CPU time.

Figure 11
Variation of CPU time versus error, for 36, 64, and 240 nodes.

The results of solving Timoshenko beam problem with INSANE software are illustrated in Figure 12 for vertical displacements of points on the beam’s mid-plane. Regular meshes were used, with 5 × 5 quadrature points in each EFG cell and also for the coupled FE elements. A linear basis was used for EFG and standard eight-node quadrilateral (Q8) elements were used in ΩFE in order to use EFG-FE coupling method. Only a strip of interface elements are placed at Γu and also at x=L, where the load is applied. The results show a very promising accuracy even using a coarse nodal distribution.

Figure 12
Deflection of Timoshenko beam in y direction using INSANE implementation.

4.2 Laplace's Equation with Mixed Boundary Conditions

The second numerical example is the Laplace's equation with mixed boundary conditions on a rectangle. The differential equation for this problem is:

(29)

The analytical solution for this problem is:

(30)

In order to investigate the EFG-FE coupling results in this problem only a strip of interface elements are located at the essential boundary locations, as shown in Figure 13.

Figure 13
Domain discretization for Laplace problem.

The obtained solutions over the whole domain are demonstrated in figures 14-16 for EFG-FE coupled, Lagrange multipliers and penalty method (α =1e4). The absolute values of error are shown for each method in two cases of 200 and 800 nodes.

Figure 14
Solution of Laplace equation using EFG-FE coupling method, a) solution, b) absolute error of approximation using 200 nodes, c) absolute error of approximation using 800 nodes.

Figure 15
Solution of Laplace equation using Lagrange multipliers, a) solution, b) absolute error of approximation using 200 nodes, c) absolute error of approximation using 800 nodes.

Figure 16
Solution of Laplace equation using penalty method, a) solution, b) absolute error of approximation using 200 nodes, c) absolute error of approximation using 800 nodes.

In Figure 14, the values of error are equal to zero at Dirichlet boundaries. Comparing these three figures, penalty method shows the largest error at Dirichlet boundaries. On the other hand, by increasing the number of nodes, the error of approximation at Dirichlet boundaries using Lagrange multipliers shows more significant decrease than the penalty method.

The values of solution at Neumann boundary of domain, x=5, are compared with the analytical solution in Figure 17 for three methods of Lagrange multipliers, penalty, and EFG-FE coupling.

Figure 17
Values of solution for Laplace equation at x=5.

4.3 Effect of Penalty Parameter on Convergence

For the case of Timoshenko beam the logarithmic values of energy norm against nodal spacing, h, are plotted in Figure 18. This plot is obtained using the EFG method for different penalty parameters, α = 1e2 = 1e8 The nodal spacing was chosen from h = 0.25 to h = 0.1.

Figure 18
Energy norm against nodal spacing for Timoshenko beam using different penalty parameters, α = 1e2 = 1e8.

It is clearly seen that the error in energy norm for penalty parameter α = 1e4 is less than the other cases, so, this value has been used for obtaining EFG results and been compared with other approaches like Lagrange multipliers and EFG-FE coupling.

Figure 19 compares the rate of convergence considering different values of nodal spacing, h, for the solution of Timoshenko beam using EFG method based on two different approaches in order to apply Dirichlet boundary conditions: Lagrange multipliers method and penalty method. In this figure, the convergence rate of EFG approximation using Lagrange multipliers and penalty method are about 0.94269 and 1.1762, respectively, which shows that EFG solution with penalty method converges faster than EFG with Lagrange multipliers.

Figure 19
Rate of convergence in energy norm for the 2-D problem using Lagrange and penalty (α = 1e4) methods.

L2-norm of error against nodal spacing for the 2D Laplace problem is illustrated in Figure 20. This figure contains L2-norm for different approaches, like FE, EFG using Lagrange multipliers, EFG using penalty method (for α = 1e2 and α = 1e4) and EFG-FE coupling approach. Although the rate of convergence for FE method, EFG using Lagrange multipliers method and EFG-FE coupled method are the same (equal to 1.99), but the error in L2-norm for EFG using Lagrange multipliers is less than the other two methods. In addition, it can be seen that EFG using penalty method with α = 1e4 shows the best results considering both convergence rate (equal to 4.449) and the accuracy.

Figure 20
Rate of convergence in energy norm for the 2D EFG problem using Lagrange and penalty (α = 1e4) methods.

5 CONCLUSION

Three methods of applying Dirichlet boundary conditions in element-free Galerkin method were investigated: Lagrange multipliers and penalty method as methods based on modification of weak form, and coupling EFG with finite element in order to directly enforce the essential boundary conditions. Numerical examples were studied to investigate different aspects of these methods such as accuracy of the results, convergence, and computational cost. Implementation of aforementioned methods of applying Dirichlet boundary conditions in EFG meshfree method in INSANE platform, an Object-oriented environment, were explained, showing the capability of this code in solving linear problems using EFG method. The coupling between EFG and FE approach works perfectly and it can be used to analyze complex problems with different boundary conditions. This capability of whole implementation in a single code helps the user to migrate easily from FE to EFG analysis, or using both methodologies in a single problem. Also, object-oriented based implementation allow us to use the advantages of other methods, such as the boundary element, hp-cloud, or generalized finite element methods that are currently available in INSANE code.

Future studies will focus on solving nonlinear problems such as damage propagation using meshfree part of the INSANE computational platform. In addition, the implemented meshfree program could also communicate with the GFEM module of the INSANE code, and using the implemented enrichment functions and methods of treating crack singularities in those classes, become an extended EFG (XEFG) method and be able to explicitly model the problems with discrete crack propagation in structures; which the last two aspects are also the subjects of ongoing works of the authors.

Acknowledgments

This research is funded by Brazilian research agencies CNPq (National Council for Scientific and Technological Developments - Grant 308785/2014-2), FAPEMIG (Foundation for Research Support of the State of Minas Gerais - Grant PPM-00669-15) and CAPES (Coordination for the Improvement of Higher Education Personnel).

References

  • Alves P.D, Barros F.B, Pitangueira R.L.S (2013) An object oriented approach to the generalized finite element method. Advances in Engineering Software, 59: 1-18.
  • Anacleto F.E.S, Ribeiro T.S.A, Ribeiro G.O, Pitangueira R.L.S, Penna S.S. (2013) An object-oriented tridimensional self-regular boundary element method implementation, Engineering Analysis with Boundary Elements 37(10):1276-1284
  • Babuška I. (1973) The finite element method with Lagrange multipliers. J. Numerical Mathematics 20: 179-192.
  • Barbosa H.J.C, Hughes T.J.R. (1992) Boundary Lagrange multipliers in finite element methods: error analysis in natural norms. J. Numerical Mathematics 62: 1-15.
  • Belytschko T, Lu Y. Y, Gu L. (1994) Element-free Galerkin methods.. Int J. Numerical Methods in Engineering 37: 229-256.
  • Belytschko T, Lu YY, Gu L. (1995a) Element-free Galerkin methods for static and dynamic fracture, Int. J. of Solids and Structures 32: 2547-70
  • Belytschko T, Organ D, Krongauz Y. (1995b) A coupled finite element-element-free Galerkin method. Computational Mechanics 17: 186-195.
  • Belytschko T, Xiao SP. (2004) A bridging domain method for coupling continua with molecular dynamics. Computer Methods in Applied Mechanics and Engineering 193: 1645-1669.
  • Belytschko T., Tabbara M. (1996) Dynamic fracture using element-free Galerkin methods. Int. J. Numerical Methods in Engineering 39: 923-38.
  • Duarte C. A, Oden J. T. (1995) Hp clouds - a meshless method to solve boundary-value problems. Technical Report, Texas Institute for Computational and Applied Mathematics, University of Texas at Austin, 95-105.
  • Fernández-Méndez S, Huerta A. (2004) Imposing essential boundary conditions in mesh-free methods. Computer Methods in Applied Mechanics and Engineering 193: 1257-1275.
  • Fleming M, Chu YA, Moran B, (1997) Enriched element-free Galerkin methods for crack-tip fields. Int. J. Numerical Methods in Engineering 40: 1483-504.
  • Fonseca F.T, Pitangueira R.L.S (2007) An object oriented class organization for dynamic geometrically non-linear. CMNE (Congress on Numerical Methods in Engineering)/CILAMCE 2007, Porto, 13-15 June, Portugal
  • Gu Y.T, Zhang L.C. (2008) Coupling of the meshfree and finite element methods for determination of the crack tip fields. Engineering Fracture Mechanics 75: 986-1004.
  • Guiamatsia I, Falzon B.G, Davies G.A.O, (2009) Element-Free Galerkin modeling of composite damage. Compos. Sci. Technol. 69: 2640-648.
  • Hecht F., (2012) New development in freefem++. Journal of Numerical Mathematics 20(3-4): 251-265.
  • Hegen D. (1996) Element free Galerkin methods in combination with finite element approaches. Computer Methods in Applied Mechanics and Engineering 135: 143-166.
  • Horstmann CS, Cornell G. (2008) Core java - fundamentals, vol. 1, 8thed. Sun Microsystems Press.
  • Hsieh Y-M, Pan M-S (2014) ESFM: An essential software framework for meshfree methods. Advances in Engineering Software 76: 133-147.
  • Huerta A, Fernandez-Mendez S. (2000) Enrichment and coupling of the finite element and meshless method. Int. J. Numerical Methods in Engineering 48: 1615 1636.
  • Krongauz Y, Belytschko T. (1996) Enforcement of essential boundary conditions in meshless approximations using finite elements. Computer Methods in Applied Mechanics and Engineering 131: 133-145.
  • Lancaster P, Salkauskas K, (1981) Surface generated by moving least square methods. J. Computational Mathematics 37: 141-58.
  • Liu L.C, Dong X.H, Li C.X. (2008) An Adaptive EFG-FE Coupling Method for the Numerical Simulation Of Extrusion Processes. Acta Metallurgica Sinica (English Letters), 21: 380-388.
  • Liu W. K, lun S, Li S, (1995) Reproducing kernel particle methods for structural dynamics. Int. J. Numerical Methods in Engineering 38: 1655-1679.
  • Lucy L. B (1997) A numerical approach to the testing of fission hypothesis. The Astronomical Journal 82: 1013-1024.
  • Malekan M, Barros F.B, Pitangueira R.L.S, Alves P.D (2016a) An object-oriented class organization for global-local generalized finite element method. Latin American Journal of Solids and Structures 13(13):2529-2551.
  • Malekan M, Barros F.B, Pitangueira R.L.S, Alves P.D, Penna S.S (2017) A computational framework for a two-scale generalized/extended finite element method: generic imposition of boundary conditions. Engineering Computations, 34(3), pp 1-32.
  • Malekan M, Barros F.B, Silva R.P, (2016b) Development and implementation of a well-conditioning approach toward generalized/extended finite element method into an object-oriented platform. Revista Interdisciplinar de Pesquisa em Engenharia 2(14):1-20.
  • Malekan M, Barros FB (2016) Well-conditioning global-local analysis using stable generalized/extended finite element method for linear elastic fracture mechanics. Computational Mechanics 58(5): 819-831.
  • Nayroles B, Touzot G, Villon P. (1992) Generalizing the finite element method: diffuse approximation and diffuse elements. Computational Mechanics 10: 307-318.
  • Patzaak B., Bittnar Z., (2001) Design of object oriented finite element code. Advances in Engineering Software 32: 759-767.
  • Peixoto R.G, Anacleto F.E.S, Ribeiro G.O, Pitangueira R.L.S (2016) A solution strategy for non-linear implicit BEM formulation using a unified constitutive modelling framework, Engineering Analysis with Boundary Elements, 64:294-310.
  • Pitkäranta J. (1979) Boundary subspaces for the finite element method of Lagrange multipliers. J. Numerical Mathematics 33: 273-289.
  • Rabczuk T, Belytschko T. (2006) Application of mesh free methods to static fracture of reinforced concrete structures. Int. J. Fracture 137(1-4): 19-49.
  • Singh I.V, Sandeep K, Prakash R. (2002) The element free Galerkin method in three dimensional steady state heat conduction. Int. J. Computational Engineering Science 3(3): 291-303.
  • Stenberg R. (1995) On some techniques for approximating boundary conditions in the finite element methods. J. Computer in Applied Mathematics 63: 139-148.
  • Yonghui L., Jun C., Song Y. (2008) Numerical simulation of three-dimensional bulk forming processes by the element-free Galerkin method. Int. J. Advanced Manufacturing Technology 36(5): 442-450.
  • Zhang X., Subbarayan G, (2006) jNURBS: An object-oriented, symbolic framework for integrated, meshless analysis and optimal design. Advances in Engineering Software 37: 287-311
  • Zhu T., Atluri S.N. (1998) A modified collocation method and a penalty formulation for enforcing the essential boundary conditions in the element free Galerkin method. Computational Mechanics 21(3): 211-222.
  • Zou W., Li X., Guo G. (2014) EFESTS, Educational finite element software for truss structure. Part I: Preprocess. International Journal of Mechanical Engineering Education 42(4): 298-306.
  • Zou W., Li X., Guo G. (2014) EFESTS, Educational finite element software for truss structure. Part II: Linear static analysis. International Journal of Mechanical Engineering Education 42(4): 307-319.

Publication Dates

  • Publication in this collection
    June 2017

History

  • Received
    13 Dec 2016
  • Reviewed
    16 Feb 2017
  • Accepted
    10 Apr 2017
Individual owner www.lajss.org - São Paulo - SP - Brazil
E-mail: lajsssecretary@gmsie.usp.br