RBF-based meshless modeling of strain localization and fracture

This work attempts to contribute further knowledge and understanding in the discipline of computational science in general and numerical modeling of discontinuity problems in particular. Of particular interest is numerical simulation of dynamic strain localization and fracture problems. The distinguishing feature in this study is the employment of neural-networks-(RBF)-based meshfree methods, which differentiates the present approach from many other computational approaches for numerical simulation of strain localization and fracture mechanics. As a result, new meshfree methods based on RBF networks, namely moving RBF-based meshless methods, have been devised and developed for solving PDEs. Unlike the conventional RBF methods, the present moving RBF is locally supported and yields sparse, banded resultant matrices, and better condition numbers. The shape functions of the new method satisfy the Kroneckerdelta property, which facilitates the imposition of the essential boundary conditions. In addition, the method is applicable to arbitrary domain and scattered nodes. To capture the characteristics of discontinuous problems, the method is further improved by special techniques including coordinate mapping and local partition of unity enrichment. Results of simulation of strain localization and fracture, presented in the latter chapters of the thesis, indicate that the proposed meshless methods have been successfully applied to model such problems

Error estimate and adaptive refinement in mixed discrete least squares meshless method

The node moving and multistage node enrichment adaptive refinement procedures are extended in mixed discrete least squares meshless (MDLSM) method for efficient analysis of elasticity problems. In the formulation of MDLSM method, mixed formulation is accepted to avoid second-order differentiation of shape functions and to obtain displacements and stresses simultaneously. In the refinement procedures, a robust error estimator based on the value of the least square residuals functional of the governing differential equations and its boundaries at nodal points is used which is inherently available from the MDLSM formulation and can efficiently identify the zones with higher numerical errors. The results are compared with the refinement procedures in the irreducible formulation of discrete least squares meshless (DLSM) method and show the accuracy and efficiency of the proposed procedures. Also, the comparison of the error norms and convergence rate show the fidelity of the proposed adaptive refinement procedures in the MDLSM method

Computational fracture modelling by an adaptive cracking particle method

Conventional element-based methods for crack modelling suffer from remeshing and mesh distortion, while the cracking particle method is meshless and requires only nodal data to discretise the problem domain so these issues are addressed. This method uses a set of crack segments to model crack paths, and crack discontinuities are obtained using the visibility criterion. It has simple implementation and is suitable for complex crack problems, but suffers from spurious cracking results and requires a large number of particles to maintain good accuracy. In this thesis, a modified cracking particle method has been developed for modelling fracture problems in 2D and 3D. To improve crack description quality, the orientations of crack segments are modified to record angular changes of crack paths, e.g. in 2D, bilinear segments replacing straight segments in the original method and in 3D, nonplanar triangular facets instead of planar circular segments, so continuous crack paths are obtained. An adaptivity approach is introduced to optimise the particle distribution, which is refined to capture high stress gradients around the crack tip and is coarsened when the crack propagates away to improve the efficiency. Based on the modified method, a multi-cracked particle method is proposed for problems with branched cracks or multiple cracks, where crack discontinuities at crack intersections are modelled by multi-split particles rather than complex enrichment functions. Different crack propagation criteria are discussed and a configurational-force-driven cracking particle method has been developed, where the crack propagating angle is directly given by the configuration force, and no decomposition of displacement and stress fields for mixed-mode fracture is required. The modified method has been applied to thermo-elastic crack problems, where the adaptivity approach is employed to capture the temperature gradients around the crack tip without using enrichment functions. Several numerical examples are used to validate the proposed methodology

Adaptive meshless point collocation methods: investigation and application to geometrically non-linear solid mechanics

Conventional mesh-based methods for solid mechanics problems suffer from issues resulting from the use of a mesh, therefore, various meshless methods that can be grouped into those based on weak or strong forms of the underlying problem have been proposed to address these problems by using only points for discretisation. Compared to weak form meshless methods, strong form meshless methods have some attractive features because of the absence of any background mesh and avoidance of the need for numerical integration, making the implementation straightforward. The objective of this thesis is to develop a novel numerical method based on strong form point collocation methods for solving problems with geometric non-linearity including membrane problems. To address some issues in existing strong form meshless methods, the local maximum entropy point collocation method is developed, where the basis functions possess some advantages such as the weak Kronecker-Delta property on boundaries. r- and h-adaptive strategies are investigated in the proposed method and are further combined into a novel rh-adaptive approach, achieving the prescribed accuracy with the optimised locations and limited number of points. The proposed meshless method with h-adaptivity is then extended to solve geometrically non-linear problems described in a Total Lagrangian formulation, where h-adaptivity is again employed after the initial calculation to improve the accuracy of the solution effciently. This geometrically non-linear method is finally developed to analyse membrane problems, in which the out-of-plane deformation for membranes complicates the governing PDEs and the use of hyperelastic materials makes the computational modelling of membrane problems challenging. The Newton-Raphson arc-length method is adopted here to capture the snap-through behaviour in hyperelastic membrane problems. Several numerical examples are presented for each proposed algorithm to validate the proposed methodology and suggestions are made for future work leading on from the findings of this thesis

POD for real-time simulation of hyperelastic soft biological tissue using the point collocation method of finite spheres

The point collocation method of finite spheres (PCMFS) is used to model the hyperelastic response of soft biological tissue in real time within the framework of virtual surgery simulation. The proper orthogonal decomposition (POD) model order reduction (MOR) technique was used to achieve reduced-order model of the problem, minimizing computational cost. The PCMFS is a physics-based meshfree numerical technique for real-time simulation of surgical procedures where the approximation functions are applied directly on the strong form of the boundary value problem without the need for integration, increasing computational efficiency. Since computational speed has a significant role in simulation of surgical procedures, the proposed technique was able to model realistic nonlinear behavior of organs in real time. Numerical results are shown to demonstrate the effectiveness of the new methodology through a comparison between full and reduced analyses for several nonlinear problems. It is shown that the proposed technique was able to achieve good agreement with the full model; moreover, the computational and data storage costs were significantly reduced

Isogeometric analysis: an overview and computer implementation aspects

Isogeometric analysis (IGA) represents a recently developed technology in computational mechanics that offers the possibility of integrating methods for analysis and Computer Aided Design (CAD) into a single, unified process. The implications to practical engineering design scenarios are profound, since the time taken from design to analysis is greatly reduced, leading to dramatic gains in efficiency. The tight coupling of CAD and analysis within IGA requires knowledge from both fields and it is one of the goals of the present paper to outline much of the commonly used notation. In this manuscript, through a clear and simple Matlab implementation, we present an introduction to IGA applied to the Finite Element (FE) method and related computer implementation aspects. Furthermore, implemen- tation of the extended IGA which incorporates enrichment functions through the partition of unity method (PUM) is also presented, where several examples for both two-dimensional and three-dimensional fracture are illustrated. The open source Matlab code which accompanies the present paper can be applied to one, two and three-dimensional problems for linear elasticity, linear elastic fracture mechanics, structural mechanics (beams/plates/shells including large displacements and rotations) and Poisson problems with or without enrichment. The Bezier extraction concept that allows FE analysis to be performed efficiently on T-spline geometries is also incorporated. The article includes a summary of recent trends and developments within the field of IGA

Coupling of adaptive refinement with variational multiscale element free Galerkin method for high gradient problems

In this thesis, a new adaptive refinement coupled with variational multiscale element free Galerkin method (EFGM) is developed for solving high gradient problems. The aim of this thesis is to propose a new framework of moving least squares (MLS) approximation with coupling method based on the variational multiscale concept. Additional new nodes will be inserted automatically at high gradient regions by adaptive algorithm based on refinement criteria. An enrichment function is embedded in the MLS approximation for the fine scale part of the problem. Besides, this new technique will be parallelized by using OpenMP which is based on shared memory architecture. The proposed new approach is first applied in two-dimensional large localized gradient problem, transient heat conduction problem as well as Burgers' equation in order to analyze the accuracy of the proposed method and validated with an available analytic solutions. The obtained numerical results show a very good agreement with the analytic solutions and is able to obtain more accurate results than the standard EFGM. It is found that the average relative error of this new method is reduced in the range of 15% to 70%. Besides, this new method is also extended to solve two-dimensional sine-Gordon solitons. The results obtained show good agreement with the published results. Moreover, the parallelization of adaptive variational multiscale EFGM can improve the computational efficiency by reducing the execution time without loss of accuracy. Therefore, the capability and robustness of this new method has the potential to investigate more complicated problems in order to produce higher precision solutions with shorter computational time

Semi-Infinite Structure Analysis with Bimodular Materials with Infinite Element

The modulus of elasticity of some materials changes under tensile and compressive states is simulated by constructing a typical material nonlinearity in a numerical analysis in this paper. The meshless Finite Block Method (FBM) has been developed to deal with 3D semi-infinite structures in the bimodular materials in this paper. The Lagrange polynomial interpolation is utilized to construct the meshless shape function with the mapping technique to transform the irregular finite domain or semi-infinite physical solids into a normalized domain. A shear modulus strategy is developed to present the nonlinear characteristics of bimodular material. In order to verify the efficiency and accuracy of FBM, the numerical results are compared with both analytical and numerical solutions provided by Finite Element Method (FEM) in four examples