Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You and Your Team.

Learn More →

Trefftz polygonal finite element for linear elasticity: convergence, accuracy, and properties

Trefftz polygonal finite element for linear elasticity: convergence, accuracy, and properties snatarajan@iitm.ac.in Department of Mechanical In this paper, the accuracy and the convergence properties of Trefftz finite element Engineering, Indian Institute method over arbitrary polygons are studied. Within this approach, the unknown of Technology-Madras, Chennai, India displacement field within the polygon is represented by the homogeneous solution Full list of author information to the governing differential equations, also called as the T-complete set. While on the is available at the end of the boundary of the polygon, a conforming displacement field is independently defined to article enforce the continuity of the field variables across the element boundary. An optimal number of T-complete functions are chosen based on the number of nodes of the polygon and the degrees of freedom per node. The stiffness matrix is computed by the hybrid formulation with auxiliary displacement frame. Results from the numerical stud- ies presented for a few benchmark problems in the context of linear elasticity show that the proposed method yields highly accurate results with optimal convergence rates. Keywords: Trefftz finite element, Polytopes, T-complete functions, Boundary integration Background Wachspress [1] introduced the concept of defining basis functions on any wedge form, which also yields interpolants on polytopes of any convex shapes. These shape func - tions are unconventional when compared to the polynomials used in the conventional finite elements. The use of elements with arbitrary number of sides provides flexibility in automatic mesh manipulation. For example, the domain can be discretized without a need to maintain a particular element topology. This is advantageous in adaptive mesh refinement, where a straightforward subdivision of individual elements usually results in hanging nodes. Traditionally, this is eliminated by introducing additional edges/faces to retain conformity. This can be alleviated if the computations are directly on meshes with hanging nodes. However, until recently elements with arbitrary number of sides did not find their applications in the computational mechanics, partly because of the associated difficulties with mesh generation and numerical integration. With the pioneering work of Alwood and Cornes [2], Sukumar and Tabarraei [3], Dasgupta [4], to name a few, now discretization of the domain with finite elements having arbitrary number of sides has gained increased attention [3, 5–10]. This has led to a new area of finite elements called ‘polygonal finite elements .’ There are different techniques to compute the basis functions over arbitrary polygons. Some of them include (a) using length and area measures  [1]; © The Author(s) 2017. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 2 of 17 (b) natural neighbor interpolants [11]; (c) maximum entropy approximant [12]; and (d) harmonic shape functions [8]. As the approximation functions over arbitrary polygonal elements are usually non-polynomial (in particular, rational polynomials) which intro- duces difficulties in the numerical integration, improving numerical integration over polytopes has gained increasing attention [3, 13–15]. It is beyond the scope of this paper to review advances in polygonal finite element methods. Interested readers are referred to the literature  [16, 17] and references therein for detailed discussion. Once the basis functions are constructed, the conventional Galerkin procedure is normally employed to solve the governing equations over the polygonal/polyhedral meshes. Sukumar [18] used Voronoï cells and natural neighbor interpolants to develop a finite difference method on unstructured grids. Rashid and Gullet  [19] proposed a variable element topology finite element method, in which shape functions for convex and non- convex elements are computed in the physical space using constrained minimization procedure. Based on the assumed stress hybrid formulation, Ghosh et al. [20] developed the Voronoï cell finite element method. Tiwary et al. [21] studied the behavior of micro - structures with irregular geometries. Liu et al. [22–24] generalized the concept of strain smoothing technique to arbitrarily shaped polygons. The main idea is to write the strain as the divergence of a spatial average of the compatible strain field. On another front, a fundamental solution less method (Scaled Boundary Method) was introduced by Wolf and Song [25]. It shares the advantages of the FEM and the boundary element method (BEM). Like the FEM, no fundamental solution is required, and like the BEM, the spatial dimension is reduced by one, since only the boundary needs to be discretized, resulting in a decrease in the total degrees of freedom. Ooi et al. [26] employed scaled boundary formulation in polygonal elements to study crack propagation. Apart from the aforementioned formulations, recent studies, among others, include devel- oping polygonal elements based on the virtual nodes  [27] and the virtual element meth- ods [28]. The other possible approach is to employ basis functions that satisfy the differential equation locally  [29, 30]. This method has been studied in detail in  [31, 32] and extended to higher order polygons in [5, 33]. Zienkiewicz [34] presented a concise discussion on dif- ferent approximation procedures to differential equations. It was shown that Trefftz-type approximation is a particular form of weighted residual approximation. This can be used to generate hybrid n fi ite elements. Earlier studies employed boundary-type approximation asso - ciated with Trefftz to develop special type finite elements, for example, elements with holes/ voids [35, 36], for plate analysis [37–39]. Recently, the idea of employing local solutions over arbitrary finite elements has been investigated in [5 , 31–33]. However, its convergence prop- erties and accuracy when applied to linear elasticity need to be investigated. In this paper, hybrid Trefftz arbitrary polygons will be formulated and its convergence properties and accuracy will be numerically studied with a few benchmark problems in the context of linear elasticity. An optimal number of T-complete functions are chosen based on the number of nodes of the polygon and degrees of freedom per node. The sali - ent features of the approach are (a) only the boundary of the element is discretized with 1D finite elements, and (b) explicit form of the shape functions and special numerical integration scheme are not required to compute the stiffness matrix. The paper commences with an overview of the governing equations for elasticity and the corresponding Galerkin form. Section “Overview of hybrid Trefftz finite element Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 3 of 17 method" introduces a hybrid Trefftz-type approximation over arbitrary polytopes. The efficiency, the accuracy, and the convergence properties of the HTFEM (Hybrid Trefftz Finite Element Method) are demonstrated with a few benchmark problems in section “Numerical examples”. The numerical results from the HTFEM are compared with the analytical results and with the polygonal FEM with Laplace/Wachspress interpolants, followed by concluding remarks in the last section. Governing equations and weak form For a 2D static linear elasticity problem defined in the domain  bounded by Ŵ = Ŵ Ŵ , u t Ŵ Ŵ =∅, in the absence of body forces, the governing equation is given by u t ∇· σ = 0 in (1) with the following conditions prescribed on the boundary: u = u in Ŵ (2) σ · n = t on Ŵ where σ is the stress tensor. The discrete equations for this problem are formulated from the Galerkin weak form: T T (∇u) D(∇δu) d� − (δu) t dŴ = 0 (3) � Ŵ where u and δu are the trial and the test functions, respectively, and D is the material con- stitutive matrix. The FEM uses the following trial and test functions: NP NP h h u (x) = N (x)d , δu (x) = N (x)δd , I I I I (4) I=1 I=1 where NP is the total number of nodes in the mesh, N is the shape function matrix, and d is the vector of degrees of freedom associated with node I. Upon substituting Eq. (4) into Eq. (3) and invoking the arbitrariness of δu, we obtain the following discretized algebraic system of equations: Kd = f (5) with K = B DB d� (6) f = N t dŴ, where K is the stiffness matrix and  is the discretized domain formed by the union of elements  . The stiffness matrix is computed over each element and assembled to the global matrix. The size of the stiffness matrix depends on the number of nodes in an element. Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 4 of 17 Generalization to arbitrary polygons The growing interest in the generalization of FE over arbitrary meshes has opened up a new area of finite elements called ‘polygonal finite elements.’ In polygonal finite elements, the number of sides of an element is not restricted to three or four as in the case of 2D. The Voronoï tessellation is a fundamental geometrical construct to generate a polygonal mesh covering a given domain. Polygonal meshes can be generated from Voronoï dia- grams. The Voronoï diagram is a subdivision of the domain into regions V (p ), such that any point in V (p ) is closer to node p than to any other node. Figure 1 shows a Voronoï I I diagram of a point P. The first-order Voronoï V(N) is a subdivision of the Euclidean space R into convex regions, mathematically: T = x ∈ R : d(x, x )< d(x, x )∀J �= I , I I J (7) where d(x , x ), the Euclidean matrix, is the distance between x and x . The quality of I J I J the generated polygonal mesh depends on the randomness in the scattered points. Fig- ure  2 shows a typical Voronoï tessellation of two sets of scattered dataset. The quality of a polygonal mesh determines the accuracy of the solution [3]. To improve the quality of the Voronoï tessellation, the generating point of each Voronoï cell can be used as its center of mass, leading to a special type of Voronoï diagram, called the centroidal Vornoï tessellation (CVT) [40]. Sieger et al. [41] presented an optimizing technique to improve the Voronoï diagrams for use in FE computations. Overview of hybrid Trefftz finite element method The basic idea in the Trefftz FEM is to employ the series of the homogeneous solution to the governing differential equation (see Eq.  1) as the approximation function to model the displacement field within the domain and an independent set of functions to rep - resent the boundary and to satisfy inter-element compatibility (see Fig.  3). The set of h h 5 4 1 S Fig. 1 Voronoï diagram of a point P Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 5 of 17 (a) (b) Fig. 2 Voronoï tessellation, where the red dots are the seed points: a scattered dataset; b polygonal mesh after few iterations functions that are used to represent the displacement field within the domain are also called as T-complete set. The displacement field within an element can be written as u = Nc, x ∈ , (8) where c are the vectors of undetermined coefficients and N are the approximation functions that are selected from the series solution of the homogeneous part of the governing differential equation (see Eq.  1). For linear elastostatics, based on the Mush- elishvili’s complex variable formulation, the N and the corresponding stress fields are given by [29] � � ReZ Jk N = Jk ImZ 2G Jk   (9)  Re(R − S ) Jk Jk T = Re(R + S ) , Jk Jk Jk   ImS Jk k k−1 k k−1 k where J = 1, 2, 3, 4, k = 1, 2, . . .; Z = iκz + iκzz , Z = κz − kzz , Z = iz , 1k 2k 3k k k−1 k−1 k−2 and Z =−z ; R = 2ikz , R = 2kz , R = 0, R = 0 ; and S = ik(k − 1)z 4k 1k 2k 3k 4k 1k k−2 k−1 k−1 z, S = k(k − 1)z z, S = iKz , S = kz . 2k 3k 4k Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 6 of 17 Fig. 3 Hybrid Trefftz polygonal element with a description of inter- and intra-element field However, the intra-element displacement field given by Eq. (8) is non-conforming across the inter-element boundary. The unknown coefficients c are computed from the external boundary conditions and/or from the continuity conditions on the inter-ele- ment boundary. Of the various methods available to enforce these conditions, in this study, we use the hybrid technique. In this technique, the elements are linked through an auxiliary conforming displacement frame which has the same form as in the conven- tional FEM. The displacement field on the element boundary, or otherwise called frame, is given by u˜ = Nq, x ∈ Ŵ, (10) where q are the unknowns of the problem and N are the standard 1D FE shape I I functions. To satisfy the inter-element continuity, a modified variational form is employed [29], given by: T −1 � = σ D σ d� − tδu dŴ + (t − t)u dŴ − tδu dŴ (11) � Ŵ Ŵ Ŵ e1 e2 eI where Ŵ is the inter-element boundary, Ŵ = Ŵ ∩ Ŵ , and Ŵ = Ŵ ∩ Ŵ . The mini - eI e1 u e e2 σ e mization of the modified variational principle given by Eq. (11) leads to the following system of algebraic equations: Kd = f (12) where T −1 K = G H G F = Nt dŴ, and Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 7 of 17 T T H = Q N dŴ; G = Q N dŴ, Ŵ Ŵ where Q = AT, and n 0 n 1 2 A = and T =∇N, 0 n n 2 1 where n and n are the outward normals. The integrals in the above equation can be 1 2 computed by employing standard Gaussian quadrature rules. It is noted that in the com- putation of the stiffness matrix, we have to compute the inverse of the matrix H. The necessary condition for the matrix H to be of full rank is m = N − 1 mim dof (13) where N is the total number of degrees of freedom of the element and m is the dof min minimum number of T-complete functions to be used. Additionally, if does not min guarantee a matrix with full rank, full rank can be achieved by suitably increasing the number of T-complete functions. Numerical examples In this section, we present the convergence and accuracy of the arbitrary polygons with local Trefftz functions using benchmark problems in the context of linear elasticity. The results from the proposed approach are compared with analytical solution where avail- able and with the conventional polygonal finite element method with Laplace interpo - lants. To discuss the results, we employ the following convention: • PFEM Polygonal finite element method with Laplace/Wachspress interpolants (con - ventional approach). The numerical integration within each element is done by sub- dividing the polygon into triangles and employing a sixth-order Dunavant quadra- ture rule. • HT-PFEM Hybrid Trefftz polygonal finite element method. Within each polygon, T-complete functions are employed to compute the stiffness matrix. One-dimen - sional Gaussian quadrature is employed along the boundary of the polygon and the order of the quadrature depends on the number of T-complete functions employed. ® ® The built-in Matlab function voronoin and Matlab functions in PolyTop [42] for build- ing the mesh-connectivity are used to create the polygonal meshes. Interested readers are referred to the corresponding author to obtain a MATLAB code for all the test cases presented in this manuscript. For the purpose of error estimation, we employ the rela- tive error in the displacement norm and in the energy norm, given by: Displacement norm h h (u − u ) · (u − u ) d� ||u − u || =  . (14) L (�) u · u d� � Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 8 of 17 Energy norm h T h (ε − ε ) D(ε − ε ) d� ||u − u || 1 = , (15) H (�) ε Dε d� h h where u, ε are the analytical solution or a reference solution and u , ε are the numerical solution. Cantilever beam bending A two-dimensional cantilever beam was subjected to a parabolic shear load at the free end. The domain is discretized with arbitrary polygonal elements. See Fig.  4 for geom- etry, boundary conditions, and representative polygonal mesh. The geometry of the can - tilever: length L = 10 and height D = 2. The material properties are as follows: Young’s modulus, E = 3e , Poisson’s ratio ν = 0.25, and the parabolic shear force P = 150. The exact solution for displacements is given by (a) (b) Fig. 4 Cantilever beam: a geometry and boundary conditions and b representative polygonal mesh. The domain is discretized with 100 elements Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 9 of 17 Py D u(x, y) = (9L − 3x)x + (2 + ν) y − 6EI 4 P D x (16) 2 2 v(x, y) =− 3νy (L − x) + (4 + 5ν) + (3L − x)x , 6EI 4 3 2 where I = D /12 is the moment of inertia, E = E, ν = ν , and E = E/(1 − ν ), ν = ν/(1 − ν) for plane stress and plane strain, respectively. The numerical convergence of the relative error in the displacement norm and the rel - ative error in the energy norm are shown in Fig. 5. The results from the HT-PFEM and (a) (b) Fig. 5 Bending of thick cantilever beam: convergence results for a the relative error in the displacement norm (L ) and b the relative error in the energy norm. The rate of convergence is also shown, where m is the average slope. In case of the polygonal SBFEM, p denotes the order of the shape functions along each edge of the polygon Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 10 of 17 Polygonal FEM are compared with the available analytical solution. Both the Polygonal 2 1 FEM and the HT-PFEM yield optimal convergence in L and H norm. It is seen that with mesh refinement, both the methods converge to the exact solution. An estimation of the convergence rate is also shown. From Fig. 5, it can be observed that the HT-PFEM yields more accurate results and better convergence rate. Infinite plate with a circular hole In this example, consider an infinite plate with a traction free hole under uni-axial ten - sion (σ = 1) along x-axis. See Fig. 6 for geometry description, boundary description, and representative polygonal mesh. The exact solution of the principal stresses in polar coor - dinates (r , θ) is given by (a) (b) Fig. 6 Infinite plate with a circular hole: a geometry and boundary conditions and b representative polygo - nal mesh L Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 11 of 17 2 4 a 3 3a σ (r, θ) = 1 − (cos 2θ + cos 4θ) + cos 4θ 2 4 r 2 2r 2 4 a 1 3a σ (r, θ) =− (cos 2θ − cos 4θ) − cos 4θ 2 4 r 2 2r (17) 2 4 a 1 3a σ (r, θ) =− (sin 2θ + sin 4θ) + sin 4θ, 2 4 r 2 2r where a is the radius of the hole. Owing to symmetry, only one quarter of the plate is modeled. The material properties are as follows: Young’s modulus E = 10 and Pois- son’s ratio ν = 0.3. In this example, analytical tractions are applied on the boundary. The domain is discretized with polygonal elements, and along each edge of the polygon, the shape function is linear. The convergence rate in terms of the displacement norm is shown in Fig.  7. The relative error in the displacement norm for the PFEM and HT- PFEM is shown in Fig. 7. It can be seen that the HT-PFEM yields more accurate results when compared with the PFEM. The HT-PFEM yields slightly a better convergence rate when compared to the PFEM. Circular beam As a last example, consider a circular cantilevered beam subjected to a prescribed dis- placement u =−0.01 at the free end. The material property, boundary conditions considered for this study, and a representative polygonal mesh are shown in Fig. 8. The material is assumed to linear elastic and in a state of plane stress. The exact solution for the elastic energy is given by U = (ln2 − 0.6) (18) Fig. 7 Infinite plate with a circular hole: convergence results for the relative error in the displacement norm (L ). The rate of convergence is also shown, where m is the average slope Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 12 of 17 (a) (b) Fig. 8 Geometry, boundary conditions, and representative polygonal mesh for a curved cantilever circular beam. The displacements at the left end are constrained and a constant displacement is applied at the other end The convergence in the relative error in the energy norm is shown in Fig.  9. It can be seen that the HT-PFEM yields more accurate results when compared with the PFEM. The HT-PFEM yields a convergence rate of 1.82 and the PFEM yields a convergence rate of 1.0. Both the methods converge to the exact energy with mesh refinement. Next, we consider two problems with complex boundary: (a) a wrench and (b) a two- dimensional crane hook, both subjected to a concentrated force, P = 210 KN. The geometry, loading, and boundary conditions are shown in Figs. 10 and 11 for the wrench and the crane hook, respectively. The material properties are as follows: Young’s modu - lus E = 3e and Poisson’s ratio ν = 0.3. The domain is discretized with arbitrary polygo - nal elements. The appropriate number of T-functions and integration points are chosen based on the number of sides of the polygonal element. As these two problems do not Fig. 9 Circular beam: convergence results for the relative error in the energy norm (L ). The rate of conver- gence is also shown, where m is the average slope Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 13 of 17 r =0.5 d =0.6 d =0.35 r =0.3 l =2 (a) 0.8 0.6 0.4 0.2 -0.2 -0.4 -0.6 -0.8 -1 0 0.5 11.5 2 2.5 (b) Fig. 10 A wrench subjected to point load: a geometry and boundary conditions and b representation polygonal mesh have a closed form solution, we use the results from a FE simulation having 29,016 and 33,104 nodes for wrench and crane hook domain, respectively, as a reference solution. The convergence of the total strain energy with mesh refinement is shown in Fig.  12 for the wrench and crane hook domain. The results from the present method are compared with the results from the conventional PFEM. It is inferred that the present method con- verges faster than the conventional PFEM. Moreover, for the same number of degrees of freedom, the present method is more accurate than the PFEM with triangulation. Concluding remarks In this paper, we studied the convergence and accuracy of hybrid Trefftz polygonal finite elements. The hybrid Trefftz finite elements were constructed by employing the T-complete set of functions and a set of independent auxiliary field on the boundary. From the numerical studies presented, it is seen that hybrid Trefftz finite elements yield Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 14 of 17 R = 20 D = 20 1 R = 35 R = 80 l = 1.86 l = 14.8 R = 50 D = 38 (a) -20 -40 -60 -40 -20 020406080100 (b) Fig. 11 A crane hook subjected to concentrated force: a geometry and boundary conditions and b repre- sentative mesh more accurate results and better convergence rate when compared to the conventional polygonal finite elements with Laplace interpolants. One of the salient features of the hybrid Trefftz approach is that special finite elements with embedded cracks/voids can be constructed. This can then be combined with the extended finite element method to Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 15 of 17 (a) (b) Fig. 12 Convergence of the total strain energy with mesh refinement: a wrench domain and b crane hook model strong and weak discontinuities and singularities within the domain. However, the success of the method relies on the knowledge of T-complete function. This is a topic for future communication. Authors’ contributions H is a Ph.D. student jointly supervised by SN and RKA, Department of Mechanical Engineering, IIT-Madras. SB and EA are external collaborators who helped in clarifying various aspects of the implementation of the method. All authors read and approved the final manuscript. Author details 1 2 Department of Mechanical Engineering, Indian Institute of Technology-Madras, Chennai, India. Faculte des Sciences, de la Technologie et de la Communication, University of Luxembourg, Luxembourg, Luxembourg. Department of Mechanical Engineering, University of Chile, Santiago, Chile. Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 16 of 17 Acknowledgements No funding was received for this project. Competing interests The authors declare that they have no competing interests. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Received: 10 August 2016 Accepted: 3 May 2017 References 1. Wachspress E (1971) A rational basis for function approximation. Springer, New York 2. Alwood R, Cornes G (1969) A polygonal finite element for plate bending problems using the assumed stress approach. Int J Numer Methods Eng 1:135 3. Sukumar N, Tabarraei A (2004) Conforming polygonal finite elements. Int J Numer Methods Eng 61:2045 4. Dasgupta G (2003) Interpolants within convex polygons: Wachspress’ shape functions. J Aerosp Eng (ASCE) 16(1):1 5. Rjasanow S, Weißer S (2012) Higher order BEM-based FEM on polygonal meshes. SIAM J Numer Anal 50:2357 6. Barros FB, de Barcellos CS, Duarte CA (2007) p-Adaptive C k generalized finite element method for arbitrary polygo - nal clouds. Comput Mech 41:175 7. da Veiga LB, Brezzi F, Cangiani A, Manzini G, Marini L, Russo A (2013) Basic principles of virtual element methods. Maths Models Methods Appl Sci 23:199 8. Bishop J (2013) A displacement-based finite element formulation for general polyhedra using harmonic shape func- tions. Int J Numer Methods Eng 97:1. doi:10.1002/nme.4562 9. Biabanaki S, Khoei A (2012) A polygonal finite element method for modeling arbitrary interfaces in large deforma- tion problems. Comput Mech 50:19. doi:10.1007/s00466-011-0668-4 10. Biabanaki S, Khoei A, Wriggers P (2013) Polygonal finite element methods for contact-impact problems on non- conformal meshes. Comput Methods Appl Mech Eng 269:198. doi:10.1016/j.cma.2013.10.025 11. Sukumar N, Moran B, Belytschko T (1998) The natural element method in solid mechanics. Int J Numer Methods Eng 43(5):839 12. Sukumar N (2013) Quadratic maximum-entropy serendipity shape functions for arbitrary planar polygons. Comput Methods Appl Mech Eng 263:27 13. Natarajan S, Bordas S, Mahapatra DR (2009) Numerical integration over arbitrary polygonal domains based on Schwarz-Christoffel conformal mapping. Int J Numer Methods Eng 80:103 14. Mousavi S, Xiao H, Sukumar N (2010) Generalized Gaussian quadrature rules on arbitrary polygons. Int J Numer Methods Eng 82(1):99 15. Talischi C, Paulino GH (2013) http://arxiv.org/pdf/1307.4423v1.pdf. (In review) 16. Fries T, Matthies H (2003) Classification and overview of meshfree methods. Tech. Rep. D38106, Institute of Scientific Computing, Technical University, Braunschweig, Hans-Sommer-Strasse 17. Sukumar N, Malsch E (2006) Recent advances in the construction of polygonal finite element interpolants. Arch Comput Methods Eng 13(1):129 18. Sukumar N (2003) Voronoi cell finite difference method for the diffusion operator on arbitrary unstructured grids. Int J Numer Methods Eng 57:1 19. Rashid M, Gullet P (2000) On a finite element method with variable element topology. Comput Methods Appl Mech Eng 190(11–12):1509 20. Moorthy S, Ghosh S (2000) Adaptivity and convergence in the Voronoi cell finite element model for analyzing heterogeneous materials. Comput Methods Appl Mech Eng 185:37 21. Tiwary A, Hu C, Ghosh S (2007) Numerical conformal mapping method based Voronoi cell finite element model for analyzing microstructures with irregular heterogeneities. Finite Elem Anal Design 43:504 22. Liu G, Nguyen T, Dai K, Lam K (2007) Theoretical aspects of the smoothed finite element method (SFEM). Int J Numer Methods Eng 71(8):902 23. Dai K, Liu G, Nguyen T (2007) An n-sided polygonal smoothed finite element method (nSFEM) for solid mechanics. Finite Elem Anal Design 43:847 24. Nguyen-Thoi T, Liu G, Nguyen-Xuan H (2011) An n-sided polygonal edge-based smoothed finite element method (nES-FEM) for solid mechanics. Int J Numer Methods Eng 27:1446 25. Wolf J, Song C (2001) The scaled boundary finite-element method—a fundamental solution-less boundary—ele - ment method. Comput Methods Appl Mech Eng 190:5551 26. Ooi ET, Song C, Tin-Loi F, Yang Z (2012) Polygon scaled boundary finite elements for crack propagation modelling. Int J Numer Methods Eng 91:319 27. Tang X, Wu S, Zheng C, Zhang J (2009) A novel virtual node method for polygonal elements. Appl Math Mech 30(10):1233 28. da Veiga L, Brezzi F, Cangiani A, Manzini G, Marini L, Russo A (2013) Basic principles of virtual element methods. Math Models Methods Appl Sci 23:199 29. Qin Q (2005) Trefftz finite element method and its applications. Appl Mech Rev 58:316 30. Copeland D, Langer U, Pusch D (2009) In: Bercovier M, Gander M, Komhuber R, Widlund O (eds). Lecture notes in Computational Science and Engineering, vol. 70 Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 17 of 17 31. Hofreither C, Langer U, Pechstein C (2010) Analysis of a non-standard finite element method based on boundary integral operators. Elect Trans Num Anal 37:413 32. Weißer S (2011) Residual error estimate for BEM-based FEM on polygonal meshes. Numerische Math 118:765 33. Weißer S (2012) Finite element methods with local trefftz trial functions. Ph.D. thesis, Universität des Saarlandes, Saarbrücken 34. Zienkiewicz OC (1997) Trefftz type approximation and the generalized finite element method—history and devel- opment. Comput Assis Mech Eng Sci 4:305 35. Piltner R (1985) Special finite elements with holes and internal cracks. Int J Numer Methods Eng 21:1471 36. Qin QH, He XQ (2009) Special elliptic hole elements of Trefftz FEM in stress concentration analysis. J Mech MEMS 1:335 37. Jirousek J, Wróblewski A, Qin Q, He X (1995) A family of quadrilateral hybrid-Trefftz p-elements for thick plate analy- sis. Comput Methods Appl Mech Eng 127:315 38. Qin Q (1995) Hybrid-Trefftz finite element method for Reissner plates on an elastic foundation. Comput Methods Appl Mech Eng 122:379 39. Choo YS, Choi N, Lee BC (2010) A new hybrid-Trefftz triangular and quadrilateral plate elements. Appl Math Model 34:14 40. Du Q, Wang D (2005) Anisotropic centroidal Voronoi tessellations and their applications. SIAM J Sci Comput 26(3):737 41. Sieger D, Alliez P, Botsch M (2010) In: Proceedings of the 19th International Meshing Roundtable, pp. 335–350 42. Talischi C, Paulino GH, Pereira A, Menezes IF (2012) PolyTop: a Matlab implementation of a general topology optimi- zation framework using unstructured polygonal finite element meshes. Struct Multidiscip Optim 45:329 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Asia Pacific Journal on Computational Engineering Springer Journals

Trefftz polygonal finite element for linear elasticity: convergence, accuracy, and properties

Loading next page...
 
/lp/springer-journals/trefftz-polygonal-finite-element-for-linear-elasticity-convergence-ptBpD800Iz
Publisher
Springer Journals
Copyright
Copyright © 2017 by The Author(s)
Subject
Engineering; Theoretical and Applied Mechanics; Computational Science and Engineering; Classical and Continuum Physics; Mathematical Applications in the Physical Sciences
eISSN
2196-1166
DOI
10.1186/s40540-017-0020-3
Publisher site
See Article on Publisher Site

Abstract

snatarajan@iitm.ac.in Department of Mechanical In this paper, the accuracy and the convergence properties of Trefftz finite element Engineering, Indian Institute method over arbitrary polygons are studied. Within this approach, the unknown of Technology-Madras, Chennai, India displacement field within the polygon is represented by the homogeneous solution Full list of author information to the governing differential equations, also called as the T-complete set. While on the is available at the end of the boundary of the polygon, a conforming displacement field is independently defined to article enforce the continuity of the field variables across the element boundary. An optimal number of T-complete functions are chosen based on the number of nodes of the polygon and the degrees of freedom per node. The stiffness matrix is computed by the hybrid formulation with auxiliary displacement frame. Results from the numerical stud- ies presented for a few benchmark problems in the context of linear elasticity show that the proposed method yields highly accurate results with optimal convergence rates. Keywords: Trefftz finite element, Polytopes, T-complete functions, Boundary integration Background Wachspress [1] introduced the concept of defining basis functions on any wedge form, which also yields interpolants on polytopes of any convex shapes. These shape func - tions are unconventional when compared to the polynomials used in the conventional finite elements. The use of elements with arbitrary number of sides provides flexibility in automatic mesh manipulation. For example, the domain can be discretized without a need to maintain a particular element topology. This is advantageous in adaptive mesh refinement, where a straightforward subdivision of individual elements usually results in hanging nodes. Traditionally, this is eliminated by introducing additional edges/faces to retain conformity. This can be alleviated if the computations are directly on meshes with hanging nodes. However, until recently elements with arbitrary number of sides did not find their applications in the computational mechanics, partly because of the associated difficulties with mesh generation and numerical integration. With the pioneering work of Alwood and Cornes [2], Sukumar and Tabarraei [3], Dasgupta [4], to name a few, now discretization of the domain with finite elements having arbitrary number of sides has gained increased attention [3, 5–10]. This has led to a new area of finite elements called ‘polygonal finite elements .’ There are different techniques to compute the basis functions over arbitrary polygons. Some of them include (a) using length and area measures  [1]; © The Author(s) 2017. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 2 of 17 (b) natural neighbor interpolants [11]; (c) maximum entropy approximant [12]; and (d) harmonic shape functions [8]. As the approximation functions over arbitrary polygonal elements are usually non-polynomial (in particular, rational polynomials) which intro- duces difficulties in the numerical integration, improving numerical integration over polytopes has gained increasing attention [3, 13–15]. It is beyond the scope of this paper to review advances in polygonal finite element methods. Interested readers are referred to the literature  [16, 17] and references therein for detailed discussion. Once the basis functions are constructed, the conventional Galerkin procedure is normally employed to solve the governing equations over the polygonal/polyhedral meshes. Sukumar [18] used Voronoï cells and natural neighbor interpolants to develop a finite difference method on unstructured grids. Rashid and Gullet  [19] proposed a variable element topology finite element method, in which shape functions for convex and non- convex elements are computed in the physical space using constrained minimization procedure. Based on the assumed stress hybrid formulation, Ghosh et al. [20] developed the Voronoï cell finite element method. Tiwary et al. [21] studied the behavior of micro - structures with irregular geometries. Liu et al. [22–24] generalized the concept of strain smoothing technique to arbitrarily shaped polygons. The main idea is to write the strain as the divergence of a spatial average of the compatible strain field. On another front, a fundamental solution less method (Scaled Boundary Method) was introduced by Wolf and Song [25]. It shares the advantages of the FEM and the boundary element method (BEM). Like the FEM, no fundamental solution is required, and like the BEM, the spatial dimension is reduced by one, since only the boundary needs to be discretized, resulting in a decrease in the total degrees of freedom. Ooi et al. [26] employed scaled boundary formulation in polygonal elements to study crack propagation. Apart from the aforementioned formulations, recent studies, among others, include devel- oping polygonal elements based on the virtual nodes  [27] and the virtual element meth- ods [28]. The other possible approach is to employ basis functions that satisfy the differential equation locally  [29, 30]. This method has been studied in detail in  [31, 32] and extended to higher order polygons in [5, 33]. Zienkiewicz [34] presented a concise discussion on dif- ferent approximation procedures to differential equations. It was shown that Trefftz-type approximation is a particular form of weighted residual approximation. This can be used to generate hybrid n fi ite elements. Earlier studies employed boundary-type approximation asso - ciated with Trefftz to develop special type finite elements, for example, elements with holes/ voids [35, 36], for plate analysis [37–39]. Recently, the idea of employing local solutions over arbitrary finite elements has been investigated in [5 , 31–33]. However, its convergence prop- erties and accuracy when applied to linear elasticity need to be investigated. In this paper, hybrid Trefftz arbitrary polygons will be formulated and its convergence properties and accuracy will be numerically studied with a few benchmark problems in the context of linear elasticity. An optimal number of T-complete functions are chosen based on the number of nodes of the polygon and degrees of freedom per node. The sali - ent features of the approach are (a) only the boundary of the element is discretized with 1D finite elements, and (b) explicit form of the shape functions and special numerical integration scheme are not required to compute the stiffness matrix. The paper commences with an overview of the governing equations for elasticity and the corresponding Galerkin form. Section “Overview of hybrid Trefftz finite element Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 3 of 17 method" introduces a hybrid Trefftz-type approximation over arbitrary polytopes. The efficiency, the accuracy, and the convergence properties of the HTFEM (Hybrid Trefftz Finite Element Method) are demonstrated with a few benchmark problems in section “Numerical examples”. The numerical results from the HTFEM are compared with the analytical results and with the polygonal FEM with Laplace/Wachspress interpolants, followed by concluding remarks in the last section. Governing equations and weak form For a 2D static linear elasticity problem defined in the domain  bounded by Ŵ = Ŵ Ŵ , u t Ŵ Ŵ =∅, in the absence of body forces, the governing equation is given by u t ∇· σ = 0 in (1) with the following conditions prescribed on the boundary: u = u in Ŵ (2) σ · n = t on Ŵ where σ is the stress tensor. The discrete equations for this problem are formulated from the Galerkin weak form: T T (∇u) D(∇δu) d� − (δu) t dŴ = 0 (3) � Ŵ where u and δu are the trial and the test functions, respectively, and D is the material con- stitutive matrix. The FEM uses the following trial and test functions: NP NP h h u (x) = N (x)d , δu (x) = N (x)δd , I I I I (4) I=1 I=1 where NP is the total number of nodes in the mesh, N is the shape function matrix, and d is the vector of degrees of freedom associated with node I. Upon substituting Eq. (4) into Eq. (3) and invoking the arbitrariness of δu, we obtain the following discretized algebraic system of equations: Kd = f (5) with K = B DB d� (6) f = N t dŴ, where K is the stiffness matrix and  is the discretized domain formed by the union of elements  . The stiffness matrix is computed over each element and assembled to the global matrix. The size of the stiffness matrix depends on the number of nodes in an element. Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 4 of 17 Generalization to arbitrary polygons The growing interest in the generalization of FE over arbitrary meshes has opened up a new area of finite elements called ‘polygonal finite elements.’ In polygonal finite elements, the number of sides of an element is not restricted to three or four as in the case of 2D. The Voronoï tessellation is a fundamental geometrical construct to generate a polygonal mesh covering a given domain. Polygonal meshes can be generated from Voronoï dia- grams. The Voronoï diagram is a subdivision of the domain into regions V (p ), such that any point in V (p ) is closer to node p than to any other node. Figure 1 shows a Voronoï I I diagram of a point P. The first-order Voronoï V(N) is a subdivision of the Euclidean space R into convex regions, mathematically: T = x ∈ R : d(x, x )< d(x, x )∀J �= I , I I J (7) where d(x , x ), the Euclidean matrix, is the distance between x and x . The quality of I J I J the generated polygonal mesh depends on the randomness in the scattered points. Fig- ure  2 shows a typical Voronoï tessellation of two sets of scattered dataset. The quality of a polygonal mesh determines the accuracy of the solution [3]. To improve the quality of the Voronoï tessellation, the generating point of each Voronoï cell can be used as its center of mass, leading to a special type of Voronoï diagram, called the centroidal Vornoï tessellation (CVT) [40]. Sieger et al. [41] presented an optimizing technique to improve the Voronoï diagrams for use in FE computations. Overview of hybrid Trefftz finite element method The basic idea in the Trefftz FEM is to employ the series of the homogeneous solution to the governing differential equation (see Eq.  1) as the approximation function to model the displacement field within the domain and an independent set of functions to rep - resent the boundary and to satisfy inter-element compatibility (see Fig.  3). The set of h h 5 4 1 S Fig. 1 Voronoï diagram of a point P Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 5 of 17 (a) (b) Fig. 2 Voronoï tessellation, where the red dots are the seed points: a scattered dataset; b polygonal mesh after few iterations functions that are used to represent the displacement field within the domain are also called as T-complete set. The displacement field within an element can be written as u = Nc, x ∈ , (8) where c are the vectors of undetermined coefficients and N are the approximation functions that are selected from the series solution of the homogeneous part of the governing differential equation (see Eq.  1). For linear elastostatics, based on the Mush- elishvili’s complex variable formulation, the N and the corresponding stress fields are given by [29] � � ReZ Jk N = Jk ImZ 2G Jk   (9)  Re(R − S ) Jk Jk T = Re(R + S ) , Jk Jk Jk   ImS Jk k k−1 k k−1 k where J = 1, 2, 3, 4, k = 1, 2, . . .; Z = iκz + iκzz , Z = κz − kzz , Z = iz , 1k 2k 3k k k−1 k−1 k−2 and Z =−z ; R = 2ikz , R = 2kz , R = 0, R = 0 ; and S = ik(k − 1)z 4k 1k 2k 3k 4k 1k k−2 k−1 k−1 z, S = k(k − 1)z z, S = iKz , S = kz . 2k 3k 4k Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 6 of 17 Fig. 3 Hybrid Trefftz polygonal element with a description of inter- and intra-element field However, the intra-element displacement field given by Eq. (8) is non-conforming across the inter-element boundary. The unknown coefficients c are computed from the external boundary conditions and/or from the continuity conditions on the inter-ele- ment boundary. Of the various methods available to enforce these conditions, in this study, we use the hybrid technique. In this technique, the elements are linked through an auxiliary conforming displacement frame which has the same form as in the conven- tional FEM. The displacement field on the element boundary, or otherwise called frame, is given by u˜ = Nq, x ∈ Ŵ, (10) where q are the unknowns of the problem and N are the standard 1D FE shape I I functions. To satisfy the inter-element continuity, a modified variational form is employed [29], given by: T −1 � = σ D σ d� − tδu dŴ + (t − t)u dŴ − tδu dŴ (11) � Ŵ Ŵ Ŵ e1 e2 eI where Ŵ is the inter-element boundary, Ŵ = Ŵ ∩ Ŵ , and Ŵ = Ŵ ∩ Ŵ . The mini - eI e1 u e e2 σ e mization of the modified variational principle given by Eq. (11) leads to the following system of algebraic equations: Kd = f (12) where T −1 K = G H G F = Nt dŴ, and Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 7 of 17 T T H = Q N dŴ; G = Q N dŴ, Ŵ Ŵ where Q = AT, and n 0 n 1 2 A = and T =∇N, 0 n n 2 1 where n and n are the outward normals. The integrals in the above equation can be 1 2 computed by employing standard Gaussian quadrature rules. It is noted that in the com- putation of the stiffness matrix, we have to compute the inverse of the matrix H. The necessary condition for the matrix H to be of full rank is m = N − 1 mim dof (13) where N is the total number of degrees of freedom of the element and m is the dof min minimum number of T-complete functions to be used. Additionally, if does not min guarantee a matrix with full rank, full rank can be achieved by suitably increasing the number of T-complete functions. Numerical examples In this section, we present the convergence and accuracy of the arbitrary polygons with local Trefftz functions using benchmark problems in the context of linear elasticity. The results from the proposed approach are compared with analytical solution where avail- able and with the conventional polygonal finite element method with Laplace interpo - lants. To discuss the results, we employ the following convention: • PFEM Polygonal finite element method with Laplace/Wachspress interpolants (con - ventional approach). The numerical integration within each element is done by sub- dividing the polygon into triangles and employing a sixth-order Dunavant quadra- ture rule. • HT-PFEM Hybrid Trefftz polygonal finite element method. Within each polygon, T-complete functions are employed to compute the stiffness matrix. One-dimen - sional Gaussian quadrature is employed along the boundary of the polygon and the order of the quadrature depends on the number of T-complete functions employed. ® ® The built-in Matlab function voronoin and Matlab functions in PolyTop [42] for build- ing the mesh-connectivity are used to create the polygonal meshes. Interested readers are referred to the corresponding author to obtain a MATLAB code for all the test cases presented in this manuscript. For the purpose of error estimation, we employ the rela- tive error in the displacement norm and in the energy norm, given by: Displacement norm h h (u − u ) · (u − u ) d� ||u − u || =  . (14) L (�) u · u d� � Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 8 of 17 Energy norm h T h (ε − ε ) D(ε − ε ) d� ||u − u || 1 = , (15) H (�) ε Dε d� h h where u, ε are the analytical solution or a reference solution and u , ε are the numerical solution. Cantilever beam bending A two-dimensional cantilever beam was subjected to a parabolic shear load at the free end. The domain is discretized with arbitrary polygonal elements. See Fig.  4 for geom- etry, boundary conditions, and representative polygonal mesh. The geometry of the can - tilever: length L = 10 and height D = 2. The material properties are as follows: Young’s modulus, E = 3e , Poisson’s ratio ν = 0.25, and the parabolic shear force P = 150. The exact solution for displacements is given by (a) (b) Fig. 4 Cantilever beam: a geometry and boundary conditions and b representative polygonal mesh. The domain is discretized with 100 elements Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 9 of 17 Py D u(x, y) = (9L − 3x)x + (2 + ν) y − 6EI 4 P D x (16) 2 2 v(x, y) =− 3νy (L − x) + (4 + 5ν) + (3L − x)x , 6EI 4 3 2 where I = D /12 is the moment of inertia, E = E, ν = ν , and E = E/(1 − ν ), ν = ν/(1 − ν) for plane stress and plane strain, respectively. The numerical convergence of the relative error in the displacement norm and the rel - ative error in the energy norm are shown in Fig. 5. The results from the HT-PFEM and (a) (b) Fig. 5 Bending of thick cantilever beam: convergence results for a the relative error in the displacement norm (L ) and b the relative error in the energy norm. The rate of convergence is also shown, where m is the average slope. In case of the polygonal SBFEM, p denotes the order of the shape functions along each edge of the polygon Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 10 of 17 Polygonal FEM are compared with the available analytical solution. Both the Polygonal 2 1 FEM and the HT-PFEM yield optimal convergence in L and H norm. It is seen that with mesh refinement, both the methods converge to the exact solution. An estimation of the convergence rate is also shown. From Fig. 5, it can be observed that the HT-PFEM yields more accurate results and better convergence rate. Infinite plate with a circular hole In this example, consider an infinite plate with a traction free hole under uni-axial ten - sion (σ = 1) along x-axis. See Fig. 6 for geometry description, boundary description, and representative polygonal mesh. The exact solution of the principal stresses in polar coor - dinates (r , θ) is given by (a) (b) Fig. 6 Infinite plate with a circular hole: a geometry and boundary conditions and b representative polygo - nal mesh L Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 11 of 17 2 4 a 3 3a σ (r, θ) = 1 − (cos 2θ + cos 4θ) + cos 4θ 2 4 r 2 2r 2 4 a 1 3a σ (r, θ) =− (cos 2θ − cos 4θ) − cos 4θ 2 4 r 2 2r (17) 2 4 a 1 3a σ (r, θ) =− (sin 2θ + sin 4θ) + sin 4θ, 2 4 r 2 2r where a is the radius of the hole. Owing to symmetry, only one quarter of the plate is modeled. The material properties are as follows: Young’s modulus E = 10 and Pois- son’s ratio ν = 0.3. In this example, analytical tractions are applied on the boundary. The domain is discretized with polygonal elements, and along each edge of the polygon, the shape function is linear. The convergence rate in terms of the displacement norm is shown in Fig.  7. The relative error in the displacement norm for the PFEM and HT- PFEM is shown in Fig. 7. It can be seen that the HT-PFEM yields more accurate results when compared with the PFEM. The HT-PFEM yields slightly a better convergence rate when compared to the PFEM. Circular beam As a last example, consider a circular cantilevered beam subjected to a prescribed dis- placement u =−0.01 at the free end. The material property, boundary conditions considered for this study, and a representative polygonal mesh are shown in Fig. 8. The material is assumed to linear elastic and in a state of plane stress. The exact solution for the elastic energy is given by U = (ln2 − 0.6) (18) Fig. 7 Infinite plate with a circular hole: convergence results for the relative error in the displacement norm (L ). The rate of convergence is also shown, where m is the average slope Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 12 of 17 (a) (b) Fig. 8 Geometry, boundary conditions, and representative polygonal mesh for a curved cantilever circular beam. The displacements at the left end are constrained and a constant displacement is applied at the other end The convergence in the relative error in the energy norm is shown in Fig.  9. It can be seen that the HT-PFEM yields more accurate results when compared with the PFEM. The HT-PFEM yields a convergence rate of 1.82 and the PFEM yields a convergence rate of 1.0. Both the methods converge to the exact energy with mesh refinement. Next, we consider two problems with complex boundary: (a) a wrench and (b) a two- dimensional crane hook, both subjected to a concentrated force, P = 210 KN. The geometry, loading, and boundary conditions are shown in Figs. 10 and 11 for the wrench and the crane hook, respectively. The material properties are as follows: Young’s modu - lus E = 3e and Poisson’s ratio ν = 0.3. The domain is discretized with arbitrary polygo - nal elements. The appropriate number of T-functions and integration points are chosen based on the number of sides of the polygonal element. As these two problems do not Fig. 9 Circular beam: convergence results for the relative error in the energy norm (L ). The rate of conver- gence is also shown, where m is the average slope Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 13 of 17 r =0.5 d =0.6 d =0.35 r =0.3 l =2 (a) 0.8 0.6 0.4 0.2 -0.2 -0.4 -0.6 -0.8 -1 0 0.5 11.5 2 2.5 (b) Fig. 10 A wrench subjected to point load: a geometry and boundary conditions and b representation polygonal mesh have a closed form solution, we use the results from a FE simulation having 29,016 and 33,104 nodes for wrench and crane hook domain, respectively, as a reference solution. The convergence of the total strain energy with mesh refinement is shown in Fig.  12 for the wrench and crane hook domain. The results from the present method are compared with the results from the conventional PFEM. It is inferred that the present method con- verges faster than the conventional PFEM. Moreover, for the same number of degrees of freedom, the present method is more accurate than the PFEM with triangulation. Concluding remarks In this paper, we studied the convergence and accuracy of hybrid Trefftz polygonal finite elements. The hybrid Trefftz finite elements were constructed by employing the T-complete set of functions and a set of independent auxiliary field on the boundary. From the numerical studies presented, it is seen that hybrid Trefftz finite elements yield Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 14 of 17 R = 20 D = 20 1 R = 35 R = 80 l = 1.86 l = 14.8 R = 50 D = 38 (a) -20 -40 -60 -40 -20 020406080100 (b) Fig. 11 A crane hook subjected to concentrated force: a geometry and boundary conditions and b repre- sentative mesh more accurate results and better convergence rate when compared to the conventional polygonal finite elements with Laplace interpolants. One of the salient features of the hybrid Trefftz approach is that special finite elements with embedded cracks/voids can be constructed. This can then be combined with the extended finite element method to Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 15 of 17 (a) (b) Fig. 12 Convergence of the total strain energy with mesh refinement: a wrench domain and b crane hook model strong and weak discontinuities and singularities within the domain. However, the success of the method relies on the knowledge of T-complete function. This is a topic for future communication. Authors’ contributions H is a Ph.D. student jointly supervised by SN and RKA, Department of Mechanical Engineering, IIT-Madras. SB and EA are external collaborators who helped in clarifying various aspects of the implementation of the method. All authors read and approved the final manuscript. Author details 1 2 Department of Mechanical Engineering, Indian Institute of Technology-Madras, Chennai, India. Faculte des Sciences, de la Technologie et de la Communication, University of Luxembourg, Luxembourg, Luxembourg. Department of Mechanical Engineering, University of Chile, Santiago, Chile. Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 16 of 17 Acknowledgements No funding was received for this project. Competing interests The authors declare that they have no competing interests. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Received: 10 August 2016 Accepted: 3 May 2017 References 1. Wachspress E (1971) A rational basis for function approximation. Springer, New York 2. Alwood R, Cornes G (1969) A polygonal finite element for plate bending problems using the assumed stress approach. Int J Numer Methods Eng 1:135 3. Sukumar N, Tabarraei A (2004) Conforming polygonal finite elements. Int J Numer Methods Eng 61:2045 4. Dasgupta G (2003) Interpolants within convex polygons: Wachspress’ shape functions. J Aerosp Eng (ASCE) 16(1):1 5. Rjasanow S, Weißer S (2012) Higher order BEM-based FEM on polygonal meshes. SIAM J Numer Anal 50:2357 6. Barros FB, de Barcellos CS, Duarte CA (2007) p-Adaptive C k generalized finite element method for arbitrary polygo - nal clouds. Comput Mech 41:175 7. da Veiga LB, Brezzi F, Cangiani A, Manzini G, Marini L, Russo A (2013) Basic principles of virtual element methods. Maths Models Methods Appl Sci 23:199 8. Bishop J (2013) A displacement-based finite element formulation for general polyhedra using harmonic shape func- tions. Int J Numer Methods Eng 97:1. doi:10.1002/nme.4562 9. Biabanaki S, Khoei A (2012) A polygonal finite element method for modeling arbitrary interfaces in large deforma- tion problems. Comput Mech 50:19. doi:10.1007/s00466-011-0668-4 10. Biabanaki S, Khoei A, Wriggers P (2013) Polygonal finite element methods for contact-impact problems on non- conformal meshes. Comput Methods Appl Mech Eng 269:198. doi:10.1016/j.cma.2013.10.025 11. Sukumar N, Moran B, Belytschko T (1998) The natural element method in solid mechanics. Int J Numer Methods Eng 43(5):839 12. Sukumar N (2013) Quadratic maximum-entropy serendipity shape functions for arbitrary planar polygons. Comput Methods Appl Mech Eng 263:27 13. Natarajan S, Bordas S, Mahapatra DR (2009) Numerical integration over arbitrary polygonal domains based on Schwarz-Christoffel conformal mapping. Int J Numer Methods Eng 80:103 14. Mousavi S, Xiao H, Sukumar N (2010) Generalized Gaussian quadrature rules on arbitrary polygons. Int J Numer Methods Eng 82(1):99 15. Talischi C, Paulino GH (2013) http://arxiv.org/pdf/1307.4423v1.pdf. (In review) 16. Fries T, Matthies H (2003) Classification and overview of meshfree methods. Tech. Rep. D38106, Institute of Scientific Computing, Technical University, Braunschweig, Hans-Sommer-Strasse 17. Sukumar N, Malsch E (2006) Recent advances in the construction of polygonal finite element interpolants. Arch Comput Methods Eng 13(1):129 18. Sukumar N (2003) Voronoi cell finite difference method for the diffusion operator on arbitrary unstructured grids. Int J Numer Methods Eng 57:1 19. Rashid M, Gullet P (2000) On a finite element method with variable element topology. Comput Methods Appl Mech Eng 190(11–12):1509 20. Moorthy S, Ghosh S (2000) Adaptivity and convergence in the Voronoi cell finite element model for analyzing heterogeneous materials. Comput Methods Appl Mech Eng 185:37 21. Tiwary A, Hu C, Ghosh S (2007) Numerical conformal mapping method based Voronoi cell finite element model for analyzing microstructures with irregular heterogeneities. Finite Elem Anal Design 43:504 22. Liu G, Nguyen T, Dai K, Lam K (2007) Theoretical aspects of the smoothed finite element method (SFEM). Int J Numer Methods Eng 71(8):902 23. Dai K, Liu G, Nguyen T (2007) An n-sided polygonal smoothed finite element method (nSFEM) for solid mechanics. Finite Elem Anal Design 43:847 24. Nguyen-Thoi T, Liu G, Nguyen-Xuan H (2011) An n-sided polygonal edge-based smoothed finite element method (nES-FEM) for solid mechanics. Int J Numer Methods Eng 27:1446 25. Wolf J, Song C (2001) The scaled boundary finite-element method—a fundamental solution-less boundary—ele - ment method. Comput Methods Appl Mech Eng 190:5551 26. Ooi ET, Song C, Tin-Loi F, Yang Z (2012) Polygon scaled boundary finite elements for crack propagation modelling. Int J Numer Methods Eng 91:319 27. Tang X, Wu S, Zheng C, Zhang J (2009) A novel virtual node method for polygonal elements. Appl Math Mech 30(10):1233 28. da Veiga L, Brezzi F, Cangiani A, Manzini G, Marini L, Russo A (2013) Basic principles of virtual element methods. Math Models Methods Appl Sci 23:199 29. Qin Q (2005) Trefftz finite element method and its applications. Appl Mech Rev 58:316 30. Copeland D, Langer U, Pusch D (2009) In: Bercovier M, Gander M, Komhuber R, Widlund O (eds). Lecture notes in Computational Science and Engineering, vol. 70 Hirshikesh et al. Asia Pac. J. Comput. Engin. (2017) 4:3 Page 17 of 17 31. Hofreither C, Langer U, Pechstein C (2010) Analysis of a non-standard finite element method based on boundary integral operators. Elect Trans Num Anal 37:413 32. Weißer S (2011) Residual error estimate for BEM-based FEM on polygonal meshes. Numerische Math 118:765 33. Weißer S (2012) Finite element methods with local trefftz trial functions. Ph.D. thesis, Universität des Saarlandes, Saarbrücken 34. Zienkiewicz OC (1997) Trefftz type approximation and the generalized finite element method—history and devel- opment. Comput Assis Mech Eng Sci 4:305 35. Piltner R (1985) Special finite elements with holes and internal cracks. Int J Numer Methods Eng 21:1471 36. Qin QH, He XQ (2009) Special elliptic hole elements of Trefftz FEM in stress concentration analysis. J Mech MEMS 1:335 37. Jirousek J, Wróblewski A, Qin Q, He X (1995) A family of quadrilateral hybrid-Trefftz p-elements for thick plate analy- sis. Comput Methods Appl Mech Eng 127:315 38. Qin Q (1995) Hybrid-Trefftz finite element method for Reissner plates on an elastic foundation. Comput Methods Appl Mech Eng 122:379 39. Choo YS, Choi N, Lee BC (2010) A new hybrid-Trefftz triangular and quadrilateral plate elements. Appl Math Model 34:14 40. Du Q, Wang D (2005) Anisotropic centroidal Voronoi tessellations and their applications. SIAM J Sci Comput 26(3):737 41. Sieger D, Alliez P, Botsch M (2010) In: Proceedings of the 19th International Meshing Roundtable, pp. 335–350 42. Talischi C, Paulino GH, Pereira A, Menezes IF (2012) PolyTop: a Matlab implementation of a general topology optimi- zation framework using unstructured polygonal finite element meshes. Struct Multidiscip Optim 45:329

Journal

Asia Pacific Journal on Computational EngineeringSpringer Journals

Published: May 25, 2017

References