Unified one-dimensional finite element for the analysis of hyperelastic soft materials and structures.

Based on the Carrera unified formulation (CUF) and first-invariant hyperelasticity, this work proposes a displacement-based high order one-dimensional (1 D) finite element model for the geometrical and physical nonlinear analysis of isotropic, slightly compressible soft material structures. Different strain energy functions are considered and they are decomposed in a volumetric and an isochoric part, the former acting as penalization of incompressibility. Given the material Jacobian tensor, the nonlinear governing equations are derived in a unified, total Lagrangian form by expanding the three-dimensional displacement field with arbitrary cross-section polynomials and using the virtual work principle. The exact analytical expressions of the elemental internal force vector and tangent matrix of the unified beam model are also provided. Several problems are addressed, including uniaxial tension, bending of a slender structure, compression of a three-dimensional block, and a thick pinched cylinder. The proposed model is compared with analytical solutions and literature results whenever possible. It is demonstrated that, although 1 D, the present CUF-based finite element can address simple to complex nonlinear hyperelastic phenomena, depending on the theory approximation order.

Quasi-static fracture analysis by coupled three-dimensional peridynamics and high order one-dimensional finite elements based on local elasticity.

This work investigates quasi-static crack propagation in specimens made of brittle materials by combining local and non-local elasticity models. The portion of the domain where the failure initiates and then propagates is modeled via three-dimensional bond-based peridynamics (PD). On the other hand, the remaining regions of the structure are analyzed with high order one-dimensional finite elements based on the Carrera unified formulation (CUF). The coupling between the two zones is realized by using Lagrange multipliers. Static solutions of different fracture problems are provided by a sequential linear analysis. The proposed approach is demonstrated to combine the advantages of the CUF-based classical continuum mechanics models and PD by providing, in an efficient manner, both the failure load and the shape of the crack pattern, even for three-dimensional problems.

Multi-level hp-finite cell method for embedded interface problems with application in biomechanics.

This work presents a numerical discretization technique for solving 3-dimensional material interface problems involving complex geometry without conforming mesh generation. The finite cell method (FCM), which is a high-order fictitious domain approach, is used for the numerical approximation of the solution without a boundary-conforming mesh. Weak discontinuities at material interfaces are resolved by using separate FCM meshes for each material sub-domain and weakly enforcing the interface conditions between the different meshes. Additionally, a recently developed hierarchical hp-refinement scheme is used to locally refine the FCM meshes to resolve singularities and local solution features at the interfaces. Thereby, higher convergence rates are achievable for nonsmooth problems. A series of numerical experiments with 2- and 3-dimensional benchmark problems is presented, showing that the proposed hp-refinement scheme in conjunction with the weak enforcement of the interface conditions leads to a significant improvement of the convergence rates, even in the presence of singularities. Finally, the proposed technique is applied to simulate a vertebra-implant model. The application showcases the method's potential as an accurate simulation tool for biomechanical problems involving complex geometry, and it demonstrates its flexibility in dealing with different types of geometric description.

High-order spectral/hp element discretisation for reaction-diffusion problems on surfaces: Application to cardiac electrophysiology.

We present a numerical discretisation of an embedded two-dimensional manifold using high-order continuous Galerkin spectral/hp elements, which provide exponential convergence of the solution with increasing polynomial order, while retaining geometric flexibility in the representation of the domain. Our work is motivated by applications in cardiac electrophysiology where sharp gradients in the solution benefit from the high-order discretisation, while the computational cost of anatomically-realistic models can be significantly reduced through the surface representation and use of high-order methods. We describe and validate our discretisation and provide a demonstration of its application to modelling electrochemical propagation across a human left atrium.