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

TLDR: A novel numerical method based on strong form point collocation methods for solving problems with geometric non-linearity including membrane problems is developed, where the basis functions possess some advantages such as the weak Kronecker-Delta property on boundaries.

Abstract: 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.