Crack propagation modelling in functionally graded materials using scaled boundary polygons
- Authors: Ooi, Ean Tat , Natarajan, Sundararajan , Song, Chongmin , Tin-Loi, Francis
- Date: 2015
- Type: Text , Journal article
- Relation: International Journal of Fracture Vol. Online first, no. 192 (2015), p. 87-105
- Full Text: false
- Reviewed:
- Description: A recently developed scaled boundary finite element formulation that can model the response of functionally graded materials is further developed to model crack propagation in two-dimensions. This formulation can accurately model the stress singularity at the crack tip in heterogeneous materials. The asymptotic behaviour at the crack tip is analytically represented in the scaled boundary shape functions of a cracked polygon. This enables accurate stress intensity factors to be computed directly from their definitions. Neither local mesh refinement nor asymptotic enrichment functions are required. This novel formulation can be implemented on polygons with an arbitrary number of sides. When modelling crack propagation, the remeshing process is more flexible and leads to only minimal changes to the global mesh structure. Six numerical examples involving crack propagation in functionally graded materials are modelled to demonstrate the salient features of the developed method.
A novel error indicator and an adaptive refinement technique using the scaled boundary finite element method
- Authors: Song, Chongmin , Ooi, Ean Tat , Pramod, Aladurthi , Natarajan, Sundararajan
- Date: 2018
- Type: Text , Journal article
- Relation: Engineering Analysis with Boundary Elements Vol. 94, no. (2018), p. 10-24
- Full Text: false
- Reviewed:
- Description: In this paper, an adaptive refinement strategy based on the scaled boundary finite element method on quadtree meshes for linear elasticity problems is discussed. Within this framework, the elements with hanging nodes are treated as polygonal elements and thus does not require special treatment. The adaptive refinement is supplemented with a novel error indicator. The local error is estimated directly from the solution of the scaled boundary governing equations. The salient feature is that it does not require any stress recovery techniques. The efficacy and the robustness of the proposed approach are demonstrated with a few numerical examples.
A scaled boundary finite element formulation over arbitrary faceted star convex polyhedra
- Authors: Natarajan, Sundararajan , Ooi, Ean Tat , Saputra, Albert , Song, Chongmin
- Date: 2017
- Type: Text , Journal article
- Relation: Engineering Analysis with Boundary Elements Vol. 80, no. (2017), p. 218-229
- Full Text: false
- Reviewed:
- Description: In this paper, a displacement based finite element framework for general three-dimensional convex polyhedra is presented. The method is based on a semi-analytical framework, the scaled boundary finite element method. The method relies on the definition of a scaling center from which the entire boundary is visible. The salient feature of the method is that the discretizations are restricted to the surfaces of the polyhedron, thus reducing the dimensionality of the problem by one. Hence, an explicit form of the shape functions inside the polyhedron is not required. Conforming shape functions defined over arbitrary polygon, such as the Wachpress interpolants are used over each surface of the polyhedron. Analytical integration is employed within the polyhedron. The proposed method passes patch test to machine precision. The convergence and the accuracy properties of the method is discussed by solving few benchmark problems in linear elasticity. © 2017 Elsevier Ltd
Crack propagation modelling in concrete using the scaled boundary finite element method with hybrid polygon-quadtree meshes
- Authors: Ooi, Ean Tat , Natarajan, Sundararajan , Song, Chongmin , Ooi, Ean Hin
- Date: 2017
- Type: Text , Journal article
- Relation: International Journal of Fracture Vol. 203, no. 1-2 (2017), p. 135-157
- Full Text: false
- Reviewed:
- Description: This manuscript presents an extension of the recently-developed hybrid polygon-quadtree-based scaled boundary finite element method to model crack propagation in concrete. This hybrid approach combines the use of quadtree cells with arbitrary sided polygons for domain discretization. The scaled boundary finite element formulation does not distinguish between quadtree cells and arbitrary sided polygons in the mesh. A single formulation is applicable to all types of cells and polygons in the mesh. This eliminates the need to develop transitional elements to bridge the cells belonging to different levels in the quadtree hierarchy. Further to this, the use of arbitrary sided polygons facilitate the accurate discretization of curved boundaries that may result during crack propagation. The fracture process zone that is characteristic in concrete fracture is modelled using zero-thickness interface elements that are coupled to the scaled boundary finite element method using a shadow domain procedure. The scaled boundary finite element method can accurately model the asymptotic stress field in the vicinity of the crack tip with cohesive tractions. This leads to the accurate computation of the stress intensity factors, which is used to determine the condition for crack propagation and the resulting direction. Crack growth can be efficiently resolved using an efficient remeshing algorithm that employs a combination of quadtree decomposition functions and simple Booleans operations. The flexibility of the scaled boundary finite element method to be formulated on arbitrary sided polygons also result in a flexible remeshing algorithm for modelling crack propagation. The developed method is validated using three laboratory experiments of notched concrete beams subjected to different loading conditions.
A scaled boundary finite element formulation with bubble functions for elasto-static analyses of functionally graded materials
- Authors: Ooi, Ean Tat , Song, Chongmin , Natarajan, Sundararajan
- Date: 2017
- Type: Text , Journal article
- Relation: Computational Mechanics Vol. 60, no. 6 (2017), p. 943-967
- Full Text:
- Reviewed:
- Description: This manuscript presents an extension of the recently-developed high order complete scaled boundary shape functions to model elasto-static problems in functionally graded materials. Both isotropic and orthotropic functionally graded materials are modelled. The high order complete properties of the shape functions are realized through the introduction of bubble-like functions derived from the equilibrium condition of a polygon subjected to body loads. The bubble functions preserve the displacement compatibility between the elements in the mesh. The heterogeneity resulting from the material gradient introduces additional terms in the polygon stiffness matrix that are integrated analytically. Few numerical benchmarks were used to validate the developed formulation. The high order completeness property of the bubble functions result in superior accuracy and convergence rates for generic elasto-static and fracture problems involving functionally graded materials. © 2017, Springer-Verlag GmbH Germany.
Crack propagation modelling in functionally graded materials using scaled boundary polygons
- Authors: Ooi, Ean Tat , Natarajan, Sundararajan , Song, Chongmin , Tin-Loi, Francis
- Date: 2015
- Type: Text , Journal article
- Relation: International Journal of Fracture Vol. 192, no. 1 (2015), p. 87-105
- Full Text: false
- Reviewed:
- Description: A recently developed scaled boundary finite element formulation that can model the response of functionally graded materials is further developed to model crack propagation in two-dimensions. This formulation can accurately model the stress singularity at the crack tip in heterogeneous materials. The asymptotic behaviour at the crack tip is analytically represented in the scaled boundary shape functions of a cracked polygon. This enables accurate stress intensity factors to be computed directly from their definitions. Neither local mesh refinement nor asymptotic enrichment functions are required. This novel formulation can be implemented on polygons with an arbitrary number of sides. When modelling crack propagation, the remeshing process is more flexible and leads to only minimal changes to the global mesh structure. Six numerical examples involving crack propagation in functionally graded materials are modelled to demonstrate the salient features of the developed method. © 2015, Springer Science+Business Media Dordrecht.
Dynamic fracture simulations using the scaled boundary finite element method on hybrid polygon-quadtree meshes
- Authors: Ooi, Ean Tat , Natarajan, Sundararajan , Song, Chongmin , Ooi, Ean Hin
- Date: 2016
- Type: Text , Journal article
- Relation: International Journal of Impact Engineering Vol. 90, no. (2016), p. 154-164
- Full Text: false
- Reviewed:
- Description: In this paper, we present an efficient computational procedure to model dynamic fracture within the framework of the scaled boundary finite element method (SBFEM). A quadtree data structure is used to discretise the domain, and 2:1 ratio between the cells is maintained. This limits the number of patterns in the quadtree decomposition and allows for efficient computation of the system matrices. The regions close to the boundary are discretised with arbitrary sided polygons so as to facilitate accurate modelling of the curved boundaries. The stiffness and the mass matrix over all the cells are computed by the SBFEM. Moreover, the semi-analytical nature of the SBFEM enables accurate modelling of the asymptotic stress fields in the vicinity of the crack tip. An efficient remeshing algorithm that combines the quadtree decomposition with simple Boolean operations is proposed to model the crack propagation. The remeshing is restricted only to a small region in the vicinity of the crack tip. The efficiency and the convergence properties of the proposed framework are demonstrated with a few benchmark problems. © 2015 Elsevier Ltd. All rights reserved.
A dual scaled boundary finite element formulation over arbitrary faceted star convex polyhedra
- Authors: Ooi, Ean Tat , Saputra, Albert , Natarajan, Sundararajan , Ooi, Ean Hin , Song, Chongmin
- Date: 2020
- Type: Text , Journal article
- Relation: Computational Mechanics Vol. 66, no. 1 (2020), p. 27-47
- Full Text: false
- Reviewed:
- Description: A novel technique to formulate arbritrary faceted polyhedral elements in three-dimensions is presented. The formulation is applicable for arbitrary faceted polyhedra, provided that a scaling requirement is satisfied and the polyhedron facets are planar. A triangulation process can be applied to non-planar facets to generate an admissible geometry. The formulation adopts two separate scaled boundary coordinate systems with respect to: (i) a scaling centre located within a polyhedron and; (ii) a scaling centre on a polyhedron’s facets. The polyhedron geometry is scaled with respect to both the scaling centres. Polygonal shape functions are derived using the scaled boundary finite element method on the polyhedron facets. The stiffness matrix of a polyhedron is obtained semi-analytically. Numerical integration is required only for the line elements that discretise the polyhedron boundaries. The new formulation passes the patch test. Application of the new formulation in computational solid mechanics is demonstrated using a few numerical benchmarks. © 2020, Springer-Verlag GmbH Germany, part of Springer Nature.
Robust modelling of implicit interfaces by the scaled boundary finite element method
- Authors: Dsouza, Shaima , Pramod, A. L. N. , Ooi, Ean Tat , Song, Chongming , Natarajan, Sundararajan
- Date: 2021
- Type: Text , Journal article
- Relation: Engineering Analysis with Boundary Elements Vol. 124, no. (2021), p. 266-286
- Full Text: false
- Reviewed:
- Description: In this paper, we propose a robust framework based on the scaled boundary finite element method to model implicit interfaces in two-dimensional differential equations in nonhomegeneous media. The salient features of the proposed work are: (a) interfaces can be implicitly defined and need not conform to the background mesh; (b) Dirichlet boundary conditions can be imposed directly along the interface; (c) does not require special numerical integration technique to compute the bilinear and the linear forms and (d) can work with an efficient local mesh refinement using hierarchical background meshes. Numerical examples involving straight interface, circular interface and moving interface problems are solved to validate the proposed technique. Further, the presented technique is compared with conforming finite element method in terms of accuracy and convergence. From the numerical studies, it is seen that the proposed framework yields solutions whose error is O(h2) in L2 norm and O(h) in the H1 semi-norm. Further the condition number increases with the mesh size similar to the FEM. © 2021 Elsevier Ltd