Browsing by Author "Reddy, Batmanathan Daya"

Now showing 1 - 2 of 2
  • Thumbnail Image
    Item
    Open Access
    A virtual element method for transversely isotropic elasticity
    (2018) Van Huyssteen, Daniel; Reddy, Batmanathan Daya
    This work studies the approximation of plane problems concerning transversely isotropic elasticity, using a low-order virtual element method (VEM). The VEM is an alternative finite element method characterised by complete freedom in determining element geometries that are otherwise polygonal in two dimensions, or polyhedral in three. Transversely isotropic materials are characterised by an axis of symmetry perpendicular to a plane of isotropy, and have applications ranging from fibre reinforcement to biological materials. The governing equations of the transversely isotropic elasticity problem are derived and a virtual element formulation of the problem is presented along with a sample implementation of the method. This work focuses on the treatment of near-incompressibility and near-inextensibility. These are explored both for homogeneous problems, in which the plane of isotropy is fixed; and non-homogeneous problems, in which the fibre directions defining the plane of isotropy vary with position. In the latter case various options are explored for approximating the non-homogeneous terms at an element level. The VEM approximations are shown through a range of numerical examples to be robust and locking-free, for a selection of element geometries, and fibre directions corresponding to mild and strong inhomogeneity.
  • Thumbnail Image
    Item
    Open Access
    Convergent finite element approximations for problems of near-incompressible and near-inextensible transversely isotropic linear elasticity
    (2019) Rasolofoson, Faraniaina; Reddy, Batmanathan Daya
    This work comprises a detailed theoretical and computational study of the boundary value problem for transversely isotropic linear elastic bodies. The main objective is the development and implementation of low-order finite element methods that are uniformly convergent in the incompressible and inextensible limits. The first step in the investigation is a study of the constitutive relation for transversely isotropic elasticity, and establishment of conditions on the five material parameters under which the relation is pointwise stable. This forms the basis for a study of well-posedness of the weak displacement-based formulation. Conforming finite element approximations are studied. The error estimate indicates the possibility of extensional locking; on the other hand, anisotropy, measured as the ratio of Young’s moduli in the fibre and transverse directions, plays a role in minimizing or even eliminating volumetric locking behaviour. Extensional locking is circumvented with the use of selective under-integration, in the context of low-order quadrilateral elements. Its equivalence with mixed and perturbed Lagrangian methods are shown. A series of numerical results illustrates the various features of the formulations considered. In a second approach, interior penalty or discontinuous Galerkin (DG) formulations of the problem are considered. Low-order approximations on triangles are adopted, with the use of three interior penalty discontinuous Galerkin methods, viz. nonsymmetric, symmetric and incomplete. It is known that these methods are uniformly convergent in the incompressible limit for the case of isotropy. This property carries over to the transversely isotropic case for moderate anisotropy. An error estimate suggests the possibility of extensional locking, and under-integration of the extensional edge terms is proposed as a remedy. This modification is shown to lead to an error estimate that is consistent with locking-free behaviour. Numerical tests confirm the uniformly convergent behaviour, at an optimal rate, of the under-integrated scheme.