-
A mixed method for 3D nonlinear elasticity using finite element exterior calculus
Authors:
Bensingh Dhas,
Jamun Kumar,
Debasish Roy,
J N Reddy
Abstract:
This article discusses a mixed FE technique for 3D nonlinear elasticity using a Hu-Washizu (HW) type variational principle. Here, the deformed configuration and sections from its cotangent bundle are taken as additional input arguments. The critical points of the HW functional enforce compatibility of these sections with the configuration, in addition to mechanical equilibrium and constitutive rel…
▽ More
This article discusses a mixed FE technique for 3D nonlinear elasticity using a Hu-Washizu (HW) type variational principle. Here, the deformed configuration and sections from its cotangent bundle are taken as additional input arguments. The critical points of the HW functional enforce compatibility of these sections with the configuration, in addition to mechanical equilibrium and constitutive relations. The present FE approximation distinguishes a vector from a 1-from, a feature not commonly found in FE approximations. This point of view permits us to construct finite elements with vastly superior performance. Discrete approximations for the differential forms appearing in the variational principle are constructed with ideas borrowed from finite element exterior calculus. The discrete equations describing mechanical equilibrium, compatibility and constitutive rule, are obtained by extemizing the discrete functional with respect to appropriate DoF, which are then solved using the Newton's method. This mixed FE technique is then applied to benchmark problems wherein conventional displacement based approximations encounter locking and checker boarding.
△ Less
Submitted 31 August, 2021;
originally announced September 2021.
-
A mixed variational principle in nonlinear elasticity using Cartan's moving frame and implementation with finite element exterior calculus
Authors:
Bensingh Dhas,
Jamun Kumar N,
Debasish Roy,
J N Reddy
Abstract:
This article offers a new perspective for the mechanics of solids using moving Cartan's frame, specifically discussing a mixed variational principle in non-linear elasticity. We treat quantities defined on the co-tangent bundles of reference and deformed configurations as additional unknowns. Such a treatment invites compatibility of the fields with base-space (configurations of the body), so that…
▽ More
This article offers a new perspective for the mechanics of solids using moving Cartan's frame, specifically discussing a mixed variational principle in non-linear elasticity. We treat quantities defined on the co-tangent bundles of reference and deformed configurations as additional unknowns. Such a treatment invites compatibility of the fields with base-space (configurations of the body), so that the configuration can be realised as a subset of the Euclidean space. We first rewrite the metric and connection using differential forms, which are further utilised to write the deformation gradient and Cauchy-Green deformation tensor in terms of frame and co-frame fields. The geometric understanding of stress as a co-vector valued 2-form fits squarely within our overall program. We show that, for a hyperelastic solid, an equation similar to the Doyle-Erciksen formula may be written for the co-vector part of stress. Using these, we write a mixed energy functional in terms of differential forms, whose extremum leads to the compatibility of deformation, constitutive rules and equations of equilibrium. Finite element exterior calculus is then utilised to construct a finite dimensional approximation for the differential forms appearing in the variational principle. These approximations are then used to construct a discrete functional which is then numerically extremised. This discertization leads to a mixed method as it uses independent approximations for differential forms related to stress and deformation gradient. The mixed variational principle is then specialized for 2D case, whose discrete approximation is applied to problems in nonlinear elasticity. An important feature of our FE technique is the lack of additional stabilization. From the numerical study, it is found that the present discretization also does not suffer form locking and related convergence issues.
△ Less
Submitted 11 August, 2021;
originally announced August 2021.
-
FlowMesher: An automatic unstructured mesh generation algorithm with applications from finite element analysis to medical simulations
Authors:
Zhujiang Wang,
Arun R. Srinivasa,
J. N. Reddy,
Adam Dubrowski
Abstract:
In this work, we propose an automatic mesh generation algorithm, FlowMesher, which can be used to generate unstructured meshes for mesh domains in any shape with minimum (or even no) user intervention. The approach can generate high-quality simplex meshes directly from scanned images in OBJ format in 2D and 3D or just from a line drawing in 2-D. Mesh grading can be easily controlled also. The Flow…
▽ More
In this work, we propose an automatic mesh generation algorithm, FlowMesher, which can be used to generate unstructured meshes for mesh domains in any shape with minimum (or even no) user intervention. The approach can generate high-quality simplex meshes directly from scanned images in OBJ format in 2D and 3D or just from a line drawing in 2-D. Mesh grading can be easily controlled also. The FlowMesher is robust and easy to be implemented and is useful for a variety of applications including surgical simulators.
The core idea of the FlowMesher is that a mesh domain is considered as an "airtight container" into which fluid particles are "injected" at one or multiple selected interior points. The particles repel each other and occupy the whole domain somewhat like blowing up a balloon. When the container is full of fluid particles and the flow is stopped, a Delaunay triangulation algorithm is employed to link the fluid particles together to generate an unstructured mesh (which is then optimized using a combination of automated mesh smoothing and element removal in 3D). The performance of the FlowMesher is demonstrated by generating meshes for several 2D and 3D mesh domains including a scanned image of a bone.
△ Less
Submitted 9 March, 2021;
originally announced March 2021.