Search

Modelling the large deformation of hyperelastic solids under plane stress conditions for arbitrary compressible and nearly incompressible material models is challenging. This is in contrast to the case of full incompressibility where the out-of-plane deformation can be entirely characterised by the in-plane components. A rigorous general procedure for the incorporation of the plane stress condition for the compressible case (including the nearly incompressible case) is provided here, accompanied by a robust and open source finite element code. An isochoric/volumetric decomposition is adopted for nearly incompressible materials yielding a robust single-field finite element formulation. The nonlinear equation for the out-of-plane component of the deformation gradient is solved using a Newton-Raphson procedure nested at the quadrature point level. The model's performance and accuracy are made clear via a series of simulations of benchmark problems. Additional challenging numerical examples of composites reinforced with particles and fibres further demonstrate the capability of this general computational framework.

Cutting-edge smart materials are transforming the domains of soft robotics, actuators, and sensors by harnessing diverse non-mechanical stimuli, such as electric and magnetic fields. Accurately modelling their physical behaviour necessitates an understanding of the complex interactions between the structural deformation and the fields in the surrounding medium. For thin shell structures, this challenge is addressed by developing a shell model that effectively incorporates the three-dimensional field it is embedded in by appropriately accounting for the relevant boundary conditions. This study presents a model for the nonlinear deformation of thin hyperelastic shells, incorporating Kirchhoff-Love assumptions and a rigorous variational approach. The shell theory is derived from 3D nonlinear elasticity by dimension reduction while preserving the boundary conditions at the top and bottom surfaces of the shell. Consequently, unlike classical shell theories, this approach can distinguish between pressure loads applied at the top and bottom surfaces, and delivers a platform to include multi-physics coupling. Numerical examples are presented to illustrate the theory and provide a physical interpretation of the novel mechanical variables of the model., Comment: arXiv admin note: text overlap with arXiv:2308.12300

A geometrically exact dimensionally reduced order model for the nonlinear deformation of thin magnetoelastic shells is presented. The Kirchhoff-Love assumptions for the mechanical fields are generalised to the magnetic variables to derive a consistent two-dimensional theory based on a rigorous variational approach. The general deformation map, as opposed to the mid-surface deformation, is considered as the primary variable resulting in a more accurate description of the nonlinear deformation. The commonly used plane stress assumption is discarded due to the Maxwell stress in the surrounding free-space requiring careful treatment on the upper and lower shell surfaces. The complexity arising from the boundary terms when deriving the Euler-Lagrange governing equations is addressed via a unique application of Green's theorem.The governing equations are solved analytically for the problem of an infinite cylindrical magnetoelastic shell. This clearly demonstrates the model's capabilities and provides a physical interpretation of the new variables in the modified variational approach. This novel formulation for magnetoelastic shells serves as a valuable tool for the accurate design of thin magneto-mechanically coupled devices.

The inflation of hyperelastic thin shells is an important and highly nonlinear problem that arises in multiple engineering applications involving severe kinematic and constitutive nonlinearities in addition to various instabilities. We present an isogeometric approach to compute the inflation of hyperelastic thin shells, following the Kirchhoff-Love hypothesis and associated large deformation. Both the geometry and the deformation field are discretized using Catmull-Clark subdivision bases which provide the C1-continuous finite element framework required for the Kirchhoff-Love shell formulation. To follow the complex nonlinear response of hyperelastic thin shells, the inflation is simulated incrementally, and each incremental step is solved via the Newton-Raphson method enriched with arc-length control. Eigenvalue analysis of the linear system after each incremental step allows for inducing bifurcation to a lower energy mode in case stability of the equilibrium is lost. The proposed method is first validated using benchmarks, and then applied to engineering applications, where we demonstrate the ability to simulate large deformation and associated complex instabilities., Comment: 25 pages, 12 figures

An isogeometric Galerkin approach for analysing the free vibrations of piezoelectric shells is presented. The shell kinematics is specialised to infinitesimal deformations and follow the Kirchhoff-Love hypothesis. Both the geometry and physical fields are discretised using Catmull-Clark subdivision bases. It provides the required C1 continuous discretisation for the Kirchhoff-Love theory. The crystalline structure of piezoelectric materials is described using an anisotropic constitutive relation. Hamilton's variational principle is applied to the dynamic analysis to derive the weak form of the governing equations. The coupled eigenvalue problem is formulated by considering the problem of harmonic vibration in the absence of external load. The formulation for the purely elastic case is verified using a spherical thin shell benchmark. Thereafter, the piezoelectric effect and vibration modes of a transverse isotropic curved plate are analysed and evaluated for the Scordelis-Lo roof problem. Finally, the eigenvalue analysis of a CAD model of a piezoelectric speaker shell structure showcases the ability of the proposed method to handle complex geometries., Comment: 24 pages, 11 figures

A novel implicit integration scheme for the Discrete Element Method (DEM) based on the variational integrator approach is presented. The numerical solver provides a fully dynamical description that, notably, reduces to an energy minimisation scheme in the quasi-static limit. A detailed derivation of the numerical method is presented for the Hookean contact model and tested against an established open source DEM package that uses the velocity-Verlet integration scheme. These tests compare results for a single collision, long-term stability and statistical quantities of ensembles of particles. Numerically, the proposed integration method demonstrates equivalent accuracy to the velocity-Verlet method.

A new class of implicit-explicit (IMEX) methods combined with a p-adaptive mixed finite element formulation is proposed to simulate the diffusion of reacting species. Hierarchical polynomial functions are used to construct an $H(\mathrm{Div})$-conforming base for the flux vectors, and a non-conforming $L^2$ base for the mass concentration of the species. The mixed formulation captures the distinct nonlinearities associated with the constitutive flux equations and the reaction terms. The IMEX method conveniently treats these two sources of nonlinearity implicitly and explicitly, respectively, within a single time-stepping framework. The combination of the p-adaptive mixed formulation and the IMEX method delivers a robust and efficient algorithm. The proposed methods eliminate the coupled effect of mesh size and time step on the algorithmic stability. A residual based a posteriori error estimate that provides an upper bound of the natural error norm is derived. The availability of such estimate which can be obtained with minimal computational effort and the hierarchical construction of the finite element spaces allow for the formulation of an efficient p-adaptive algorithm. A series of numerical examples demonstrate the performance of the approach. It is shown that the method with the p-adaptive strategy accurately solves problems involving travelling waves, and those with discontinuities and singularities. The flexibility of the formulation is also illustrated via selected applications in pattern formation and electrophysiology.

We analyse here the problem of large deformation of dielectric elastomeric membranes under coupled electromechanical loading. Extremely large deformations (enclosed volume changes of 100 times and greater) of a toroidal membrane are studied by the use of a variational formulation that accounts for the total energy due to mechanical and electrical fields. A modified shooting method is adopted to solve the resulting system of coupled and highly nonlinear ordinary differential equations. We demonstrate the occurrence of limit point, wrinkling, and symmetry-breaking buckling instabilities in the solution of this problem. Onset of each of these "reversible" instabilities depends significantly on the ratio of the mechanical load to the electric load, thereby providing a control mechanism for state switching., Comment: 56 pages, 21 figures

Standard low-order finite elements, which perform well for problems involving compressible elastic materials, are known to under-perform when nearly incompressible materials are involved, commonly exhibiting the locking phenomenon. Interior penalty (IP) discontinuous Galerkin methods have been shown to circumvent locking when simplicial elements are used. The same IP methods, however, result in locking on meshes of quadrilaterals. The authors have shown in earlier work that under-integration of specified terms in the IP formulation eliminates the locking problem for rectangular elements. Here it is demonstrated through an extensive numerical investigation that the effect of using under-integration carries over successfully to meshes of more general quadrilateral elements, as would likely be used in practical applications, and results in accurate displacement approximations. Uniform convergence with respect to the compressibility parameter is shown numerically. Additionally, a stress approximation obtained here by postprocessing shows good convergence in the incompressible limit., Comment: 41 pages, 33 figures

An isogeometric approach for solving the Laplace-Beltrami equation on a two-dimensional manifold embedded in three-dimensional space using a Galerkin method based on Catmull-Clark subdivision surfaces is presented and assessed. The scalar-valued Laplace-Beltrami equation requires only C0 continuity and is adopted to elucidate key features and properties of the isogeometric method using Catmull-Clark subdivision surfaces. Catmull-Clark subdivision bases are used to discretise both the geometry and the physical field. A fitting method generates control meshes to approximate any given geometry with Catmull-Clark subdivision surfaces. The performance of the Catmull-Clark subdivision method is compared to the conventional finite element method. Subdivision surfaces without extraordinary vertices show the optimal convergence rate. However, extraordinary vertices introduce error, which decreases the convergence rate. A comparative study shows the effect of the number and valences of the extraordinary vertices on accuracy and convergence. An adaptive quadrature scheme is shown to reduce the error., Comment: 38 pages, 76 figures

Flexoelectricity is characterised by the coupling of the gradient of the deformation and the electrical polarization in a dielectric material. A novel micromorphic approach is presented to accommodate the resulting higher-order gradient contributions arising in this highly-nonlinear and coupled problem within a classical finite element setting. The formulation accounts for all material and geometric nonlinearities, as well as the coupling between the mechanical, electrical and micromorphic fields. The highly-nonlinear system of governing equations are derived using the Dirichlet principle and solved using the finite element method. A series of numerical examples serve to elucidate the theory and to provide insight into this fascinating effect.

A unified classification framework for models of extended plasticity is presented. The models include well known micromorphic and strain gradient plasticity formulations. A unified treatment is possible due to the representation of strain gradient plasticity as an Eringen-type micromorphic continua. The classification is based on the form of the energetic and dissipative model structures and exploits the framework of dissipation-consistent modelling to elucidate the flow relation and yield condition. Models are identified as either serial or parallel. This designation is also applicable to familiar models of classical plasticity. Particular attention is paid to the rate-dependent problem arising from the choice of a smooth dissipation potential. The inability to locally determine the region of admissible stresses for the non-smooth (rate-independent) parallel models of plasticity is made clear.

The potentially significant role of the surface of an elastic body in the overall response of the continuum can be described using the mature theory of surface elasticity. The objective of this contribution is to detail the finite element approximation of the underlying governing equations (both in the volume and on its surface) and their solution using the open-source finite element library deal.II. The fully-nonlinear (geometric and material) setting is considered. The nonlinear problem is solved using a Newton--Raphson procedure wherein the tangent contributions from the volume and surface are computed exactly. The finite element formulation is implemented within the total Lagrangian framework and a Bubnov-Galerkin spatial discretization of the volume and the surface employed. The surface is assumed material. A map between the degrees of freedom on the surface and on the boundary of the volume is used to allocate the contribution from the surface to the global system matrix and residual vector. The deal.II library greatly facilitates the computation of the various surface operators, allowing the numerical implementation to closely match the theory developed in a companion paper. Key features of the theory and the numerical implementation are elucidated using a series of benchmark example problems. The full, documented source code is provided.