-
Response to David Steigmann's discussion of our paper
Authors:
Thang X. Duong,
Mikhail Itskov,
Roger A. Sauer
Abstract:
We respond to David Steigmann's discussion of our paper "A general theory for anisotropic Kirchhoff-Love shells with in-plane bending of embedded fibers, Math. Mech. Solids, 28(5):1274-1317" (arXiv:2101.03122). His discussion allows us to clarify two misleading statements in our original paper, and confirm that its formulation is fully consistent with the formulation of Steigmann. We also demonstr…
▽ More
We respond to David Steigmann's discussion of our paper "A general theory for anisotropic Kirchhoff-Love shells with in-plane bending of embedded fibers, Math. Mech. Solids, 28(5):1274-1317" (arXiv:2101.03122). His discussion allows us to clarify two misleading statements in our original paper, and confirm that its formulation is fully consistent with the formulation of Steigmann. We also demonstrate that some of our original statements criticized by Steigmann are not wrong.
△ Less
Submitted 29 February, 2024;
originally announced February 2024.
-
A general isogeometric finite element formulation for rotation-free shells with in-plane bending of embedded fibers
Authors:
Thang Xuan Duong,
Mikhail Itskov,
Roger Andrew Sauer
Abstract:
This paper presents a general, nonlinear isogeometric finite element formulation for rotation-free shells with embedded fibers that captures anisotropy in stretching, shearing, twisting and bending -- both in-plane and out-of-plane. These capabilities allow for the simulation of large sheets of heterogeneous and fibrous materials either with or without matrix, such as textiles, composites, and pan…
▽ More
This paper presents a general, nonlinear isogeometric finite element formulation for rotation-free shells with embedded fibers that captures anisotropy in stretching, shearing, twisting and bending -- both in-plane and out-of-plane. These capabilities allow for the simulation of large sheets of heterogeneous and fibrous materials either with or without matrix, such as textiles, composites, and pantographic structures. The work is a computational extension of our earlier theoretical work [1] that extends existing Kirchhoff-Love shell theory to incorporate the in-plane bending resistance of initially straight or curved fibers. The formulation requires only displacement degrees-of-freedom to capture all mentioned modes of deformation. To this end, isogeometric shape functions are used in order to satisfy the required $C^1$-continuity for bending across element boundaries. The proposed formulation can admit a wide range of material models, such as surface hyperelasticity that does not require any explicit thickness integration. To deal with possible material instability due to fiber compression, a stabilization scheme is added. Several benchmark examples are used to demonstrate the robustness and accuracy of the proposed computational formulation.
△ Less
Submitted 5 June, 2023; v1 submitted 1 October, 2021;
originally announced October 2021.
-
A general theory for anisotropic Kirchhoff-Love shells with in-plane bending of embedded fibers
Authors:
Thang Xuan Duong,
Vu Ngoc Khiêm,
Mikhail Itskov,
Roger Andrew Sauer
Abstract:
This work presents a generalized Kirchhoff-Love shell theory that can explicitly capture fiber-induced anisotropy not only in stretching and out-of-plane bending, but also in in-plane bending. This setup is particularly suitable for heterogeneous and fibrous materials such as textiles, biomaterials, composites and pantographic structures. The presented theory is a direct extension of classical Kir…
▽ More
This work presents a generalized Kirchhoff-Love shell theory that can explicitly capture fiber-induced anisotropy not only in stretching and out-of-plane bending, but also in in-plane bending. This setup is particularly suitable for heterogeneous and fibrous materials such as textiles, biomaterials, composites and pantographic structures. The presented theory is a direct extension of classical Kirchhoff-Love shell theory to incorporate the in-plane bending resistance of fibers. It also extends existing second-gradient Kirchhoff-Love shell theory for initially straight fibers to initially curved fibers. To describe the additional kinematics of multiple fiber families, a so-called in-plane curvature tensor -- which is symmetric and of second order -- is proposed. The effective stress tensor and the in-plane and out-of-plane moment tensors are then identified from the mechanical power balance. These tensors are all second order and symmetric in general. Constitutive equations for hyperelastic materials are derived from different expressions of the mechanical power balance. The weak form is also presented as it is required for computational shell formulations based on rotation-free finite element discretizations.
△ Less
Submitted 5 June, 2023; v1 submitted 23 December, 2020;
originally announced January 2021.
-
A chemo-mechano-thermodynamical contact theory for adhesion, friction and (de)bonding reactions
Authors:
Roger A. Sauer,
Thang X. Duong,
Kranthi K. Mandadapu
Abstract:
This work presents a self-contained continuum formulation for coupled chemical, mechanical and thermal contact interactions. The formulation is very general and hence admits arbitrary geometry, deformation and material behavior. All model equations are derived rigorously from the balance laws of mass, momentum, energy and entropy in the framework of irreversible thermodynamics, thus exposing all t…
▽ More
This work presents a self-contained continuum formulation for coupled chemical, mechanical and thermal contact interactions. The formulation is very general and hence admits arbitrary geometry, deformation and material behavior. All model equations are derived rigorously from the balance laws of mass, momentum, energy and entropy in the framework of irreversible thermodynamics, thus exposing all the coupling present in the field equations and constitutive relations. In the process, the conjugated kinematic and kinetic variables for mechanical, thermal and chemical contact are identified, and the analogies between mechanical, thermal and chemical contact are highlighted. Particular focus is placed on the thermodynamics of chemical bonding distinguishing between exothermic and endothermic contact reactions. Distinction is also made between long-range, non-touching surface interactions and short-range, touching contact. For all constitutive relations examples are proposed and discussed comprehensively with particular focus on their coupling. Finally, three analytical test cases are presented that illustrate the thermo-chemo-mechanical contact coupling and are useful for verifying computational models. While the main novelty is the extension of existing contact formulations to chemical contact, the presented formulation also sheds new light on thermo-mechanical contact, since it is consistently derived from basic principles using only few assumptions.
△ Less
Submitted 11 August, 2021; v1 submitted 29 December, 2020;
originally announced December 2020.
-
Isogeometric continuity constraints for multi-patch shells governed by fourth-order deformation and phase field models
Authors:
Karsten Paul,
Christopher Zimmermann,
Thang X. Duong,
Roger A. Sauer
Abstract:
This work presents numerical techniques to enforce continuity constraints on multi-patch surfaces for three distinct problem classes. The first involves structural analysis of thin shells that are described by general Kirchhoff-Love kinematics. Their governing equation is a vector-valued, fourth-order, nonlinear, partial differential equation (PDE) that requires at least $C^1$-continuity within a…
▽ More
This work presents numerical techniques to enforce continuity constraints on multi-patch surfaces for three distinct problem classes. The first involves structural analysis of thin shells that are described by general Kirchhoff-Love kinematics. Their governing equation is a vector-valued, fourth-order, nonlinear, partial differential equation (PDE) that requires at least $C^1$-continuity within a displacement-based finite element formulation. The second class are surface phase separations modeled by a phase field. Their governing equation is the Cahn-Hilliard equation - a scalar, fourth-order, nonlinear PDE - that can be coupled to the thin shell PDE. The third class are brittle fracture processes modeled by a phase field approach. In this work, these are described by a scalar, fourth-order, nonlinear PDE that is similar to the Cahn-Hilliard equation and is also coupled to the thin shell PDE. Using a direct finite element discretization, the two phase field equations also require at least a $C^1$-continuous formulation. Isogeometric surface discretizations - often composed of multiple patches - thus require constraints that enforce the $C^1$-continuity of displacement and phase field. For this, two numerical strategies are presented: For this, two numerical strategies are presented: A Lagrange multiplier formulation and a penalty method. The curvilinear shell model including the geometrical constraints is taken from Duong et al. (2017) and it is extended to model the coupled phase field problems on thin shells of Zimmermann et al. (2019) and Paul et al. (2020) on multi-patches. Their accuracy and convergence are illustrated by several numerical examples considering deforming shells, phase separations on evolving surfaces, and dynamic brittle fracture of thin shells.
△ Less
Submitted 6 October, 2021; v1 submitted 16 January, 2020;
originally announced January 2020.
-
A modified Coulomb's law for the tangential debonding of osseointegrated implants
Authors:
Katharina Immel,
Thang X. Duong,
Vu-Hieu Nguyen,
Guillaume Haiat,
Roger A. Sauer
Abstract:
Cementless implants are widely used in orthopedic and oral surgery. However, debonding-related failure still occurs at the bone-implant interface. It remains difficult to predict such implant failure since the underlying osseointegration phenomena are still poorly understood. Especially in terms of friction and adhesion at the macro-scale, there is a lack of data and reliable models. The aim of th…
▽ More
Cementless implants are widely used in orthopedic and oral surgery. However, debonding-related failure still occurs at the bone-implant interface. It remains difficult to predict such implant failure since the underlying osseointegration phenomena are still poorly understood. Especially in terms of friction and adhesion at the macro-scale, there is a lack of data and reliable models. The aim of this work is to present a new friction formulation that can model the tangential contact behavior between osseointegrated implants and bone tissue, with focus on debonding. The classical Coulomb's law is combined with a state variable friction law to model a displacement-dependent friction coefficient. A smooth state function, based on the sliding distance, is used to model implant debonding. The formulation is implemented in a 3D nonlinear finite element framework, and it is calibrated with experimental data and compared to an analytical model for mode III cleavage of a coin-shaped, titanium implant (Mathieu et al. 2012). Overall, the results show close agreement with the experimental data, especially the peak and the softening part of the torque curve with a relative error of less than 2.25 %. In addition, better estimates of the bone's shear modulus and the adhesion energy are obtained. The proposed model is particularly suitable to account for partial osseointegration, as is also shown.
△ Less
Submitted 21 February, 2020; v1 submitted 12 August, 2019;
originally announced August 2019.
-
A NURBS-based Inverse Analysis for Reconstruction of Nonlinear Deformations of Thin Shell Structures
Authors:
N. Vu-Bac,
T. X. Duong,
T. Lahmer,
X. Zhuang,
R. A. Sauer,
H. S. Parke,
T. Rabczuk
Abstract:
This article presents original work combining a NURBS-based inverse analysis with both kinematic and constitutive nonlinearities to recover the applied loads and deformations of thin shell structures. The inverse formulation is tackled by gradient-based optimization algorithms based on computed and measured displacements at a number of discrete locations. The proposed method allows accurately reco…
▽ More
This article presents original work combining a NURBS-based inverse analysis with both kinematic and constitutive nonlinearities to recover the applied loads and deformations of thin shell structures. The inverse formulation is tackled by gradient-based optimization algorithms based on computed and measured displacements at a number of discrete locations. The proposed method allows accurately recovering the target shape of shell structures such that instabilities due to snapping and buckling are captured. The results obtained show good performance and applicability of the proposed algorithms to computer-aided manufacturing of shell structures.
△ Less
Submitted 5 February, 2019;
originally announced February 2019.
-
A concise frictional contact formulation based on surface potentials and isogeometric discretization
Authors:
Thang X. Duong,
Roger A. Sauer
Abstract:
This work presents a concise theoretical and computational framework for the finite element formulation of frictional contact problems with arbitrarily large deformation and sliding. The aim of this work is to extend the contact theory based on surface potentials (Sauer and De Lorenzis, 2013) to account for friction. Coulomb friction under isothermal conditions is considered here. For a consistent…
▽ More
This work presents a concise theoretical and computational framework for the finite element formulation of frictional contact problems with arbitrarily large deformation and sliding. The aim of this work is to extend the contact theory based on surface potentials (Sauer and De Lorenzis, 2013) to account for friction. Coulomb friction under isothermal conditions is considered here. For a consistent friction formulation, we start with the first and second laws of thermodynamics and derive the governing equations at the contact interface. A so-called interacting gap can then be defined as a kinematic variable unifying both sliding/sticking and normal/tangential contact. A variational principle for the frictional system can then be formulated based on a purely kinematical constraint. The direct elimination approach applied to the tangential part of this constraint leads to the so-called moving friction cone approach of Wriggers and Haraldsson (2003). Compared with existing friction formulations, our approach reduces the theoretical and computational complexity. Several numerical examples are presented to demonstrate the accuracy and robustness of the proposed friction formulation.
△ Less
Submitted 2 April, 2019; v1 submitted 30 August, 2018;
originally announced August 2018.
-
A segmentation-free isogeometric extended mortar contact method
Authors:
Thang Xuan Duong,
Laura De Lorenzis,
Roger A. Sauer
Abstract:
This paper presents a new isogeometric mortar contact formulation based on an extended finite element interpolation to capture physical pressure discontinuities at the contact boundary. The so called two-half-pass algorithm is employed, which leads to an unbiased formulation and, when applied to the mortar setting, has the additional advantage that the mortar coupling term is no longer present in…
▽ More
This paper presents a new isogeometric mortar contact formulation based on an extended finite element interpolation to capture physical pressure discontinuities at the contact boundary. The so called two-half-pass algorithm is employed, which leads to an unbiased formulation and, when applied to the mortar setting, has the additional advantage that the mortar coupling term is no longer present in the contact forces. As a result, the computationally expensive segmentation at overlapping master-slave element boundaries, usually required in mortar methods (although often simplified with loss of accuracy), is not needed from the outset. For the numerical integration of general contact problems, the so-called refined boundary quadrature is employed, which is based on adaptive partitioning of contact elements along the contact boundary. The contact patch test shows that the proposed formulation passes the test without using either segmentation or refined boundary quadrature. Several numerical examples are presented to demonstrate the robustness and accuracy of the proposed formulation.
△ Less
Submitted 12 July, 2018; v1 submitted 4 December, 2017;
originally announced December 2017.
-
A new shell formulation for graphene structures based on existing ab-initio data
Authors:
Reza Ghaffari,
Thang X. Duong,
Roger A. Sauer
Abstract:
An existing hyperelastic membrane model for graphene calibrated from ab-initio data (Kumar and Parks, 2014) is adapted to curvilinear coordinates and extended to a rotation-free shell formulation based on isogeometric finite elements. Therefore, the membrane model is extended by a hyperelastic bending model that reflects the ab-inito data of Kudin et al. (2001). The proposed formulation can be imp…
▽ More
An existing hyperelastic membrane model for graphene calibrated from ab-initio data (Kumar and Parks, 2014) is adapted to curvilinear coordinates and extended to a rotation-free shell formulation based on isogeometric finite elements. Therefore, the membrane model is extended by a hyperelastic bending model that reflects the ab-inito data of Kudin et al. (2001). The proposed formulation can be implemented straight-forwardly into an existing finite element package, since it does not require the description of molecular interactions. It thus circumvents the use of interatomic potentials that tend to be less accurate than ab-initio data. The proposed shell formulation is verified and analyzed by a set of simple test cases. The results are in agreement to analytical solutions and satisfy the FE patch test. The performance of the shell formulation for graphene structures is illustrated by several numerical examples. The considered examples are indentation and peeling of graphene and torsion, bending and axial stretch of carbon nanotubes. Adhesive substrates are modeled by the Lennard-Jones potential and a coarse grained contact model. In principle, the proposed formulation can be extended to other 2D materials.
△ Less
Submitted 1 December, 2017; v1 submitted 28 December, 2016;
originally announced December 2016.
-
A stabilized finite element formulation for liquid shells and its application to lipid bilayers
Authors:
Roger A. Sauer,
Thang X. Duong,
Kranthi K. Mandadapu,
David J. Steigmann
Abstract:
This paper presents a new finite element (FE) formulation for liquid shells that is based on an explicit, 3D surface discretization using $C^1$-continuous finite elements constructed from NURBS interpolation. Both displacement-based and mixed FE formulations are proposed. The latter is needed for area-incompressible material behavior, where penalty-type regularizations can lead to misleading resul…
▽ More
This paper presents a new finite element (FE) formulation for liquid shells that is based on an explicit, 3D surface discretization using $C^1$-continuous finite elements constructed from NURBS interpolation. Both displacement-based and mixed FE formulations are proposed. The latter is needed for area-incompressible material behavior, where penalty-type regularizations can lead to misleading results. In order to obtain quasi-static solutions, several numerical stabilization schemes are proposed based on either stiffness, viscosity or projection. Several numerical examples are considered in order to illustrate the accuracy and the capabilities of the proposed formulation, and to compare the different stabilization schemes. The presented formulation is capable of simulating non-trivial surface shapes associated with tube formation and protein-induced budding of lipid bilayers. In the latter case, the presented formulation yields non-axisymmetric solutions, which have not been observed in previous simulations. It is shown that those non-axisymmetric shapes are preferred over axisymmetric ones.
△ Less
Submitted 17 January, 2016; v1 submitted 15 January, 2016;
originally announced January 2016.
-
A computational formulation for constrained solid and liquid membranes considering isogeometric finite elements
Authors:
Roger A. Sauer,
Thang X. Duong,
Callum J. Corbett
Abstract:
A geometrically exact membrane formulation is presented that is based on curvilinear coordinates and isogeometric finite elements, and is suitable for both solid and liquid membranes. The curvilinear coordinate system is used to describe both the theory and the finite element equations of the membrane. In the latter case this avoids the use of local cartesian coordinates at the element level. Cons…
▽ More
A geometrically exact membrane formulation is presented that is based on curvilinear coordinates and isogeometric finite elements, and is suitable for both solid and liquid membranes. The curvilinear coordinate system is used to describe both the theory and the finite element equations of the membrane. In the latter case this avoids the use of local cartesian coordinates at the element level. Consequently, no transformation of derivatives is required. The formulation considers a split of the in-plane and out-of-plane membrane contributions, which allows the construction of a stable formulation for liquid membranes with constant surface tension. The proposed membrane formulation is general, and accounts for dead and live loading, as well as enclosed volume, area, and contact constraints. The new formulation is illustrated by several challenging examples, considering linear and quadratic Lagrange elements, as well as isogeometric elements based on quadratic NURBS and cubic T-splines. It is seen that the isogeometric elements are much more accurate than standard Lagrange elements. The gain is especially large for the liquid membrane formulation since it depends explicitly on the surface curvature.
△ Less
Submitted 17 October, 2012;
originally announced October 2012.