Related papers: Stable numerics for finite-strain elasticity
A method is presented to calculate from first principles the higher-order elastic constants of a solid material. The method relies on finite strain deformations, a density functional theory approach to calculate the Cauchy stress tensor,…
A first-principles method is presented to calculate elastic constants up to the fourth order of crystals with the cubic and hexagonal symmetries. The method relies on the numerical differentiation of the second Piola-Kirchhoff stress tensor…
Modeling microstructural evolution at large strains requires mechanical formulations that remain thermodynamically consistent while capturing significant lattice rotations and transformation-induced stresses. However, most existing…
Elastic constants are central material properties, frequently reported in experimental and theoretical studies. While their computation is straightforward in the absence of thermal fluctuations, finite--temperature methods often suffer from…
Tensor analysis provides a frame-invariant foundation for continuum mechanics, yet numerical implementations rely on matrix representations expressed in user-selected bases. When these bases are non-Cartesian and non-orthonormal, additional…
An efficient and reliable stress computation algorithm is presented, which is based on implicit integration of the local evolution equations of multiplicative finite-strain plasticity/viscoplasticity. The algorithm is illustrated by an…
A new class of semi-implicit numerical schemes for linear advection equation on Cartesian grids is derived that is inspired by so-called $\kappa$-schemes used with fully explicit discretizations for this type of problems. Opposite to fully…
The effect of finite temperature $T$ and finite strain rate $\dot\gamma$ on the statistical physics of plastic deformations in amorphous solids made of $N$ particles is investigated. We recognize three regimes of temperature where the…
We continue our investigation of finite deformation linear viscoelastodynamics by focusing on constructing accurate and reliable numerical schemes. The concrete thermomechanical foundation developed in the previous study paves the way for…
We introduce a new class of mixed finite element methods for 2D and 3D compressible nonlinear elasticity. The independent unknowns of these conformal methods are displacement, displacement gradient, and the first Piola-Kirchhoff stress…
In our recent work on iterative computation in hardware, we showed that arbitrary-precision solvers can perform more favorably than their traditional arithmetic equivalents when the latter's precisions are either under- or over-budgeted for…
In isotropic nonlinear elasticity the corotational stability postulate (CSP) is the requirement that \begin{equation*} \langle\frac{\mathrm{D}^{\circ}}{\mathrm{D} t}[\sigma] , D \rangle > 0 \quad \forall \ D \in \text{Sym}(3)\setminus \{0\}…
Many statistical estimators are defined as the fixed point of a data-dependent operator, with estimators based on minimizing a cost function being an important special case. The limiting performance of such estimators depends on the…
Presented here is a study of well-posedness and asymptotic stability of a "degenerately damped" PDE modeling a vibrating elastic string. The coefficient of the damping may vanish at small amplitudes thus weakening the effect of the…
In order to accelerate implementation of hyperelastic materials for finite element analysis, we developed an automatic numerical algorithm that only requires the strain energy function. This saves the effort on analytical derivation and…
Computational formulations for large strain, polyconvex, nearly incompressible elasticity have been extensively studied, but research on enhancing solution schemes that offer better tradeoffs between accuracy, robustness, and computational…
Numerically simulating deformations in thin elastic sheets is a challenging problem in computational mechanics due to destabilizing compressive stresses that result in wrinkling. Determining the location, structure, and evolution of…
Tangent stabilised large strain isotropic elasticity was recently proposed by Poya et al. [1] wherein by working directly with principal stretches the entire eigenstructure of constitutive and geometric/initial stiffness terms were found in…
We introduce conformal mixed finite element methods for $2$D and $3$D incompressible nonlinear elasticity in terms of displacement, displacement gradient, the first Piola-Kirchhoff stress tensor, and pressure, where finite elements for the…
Fluctuations of the instantaneous local Lagrangian strain $\epsilon_{ij}(\bf{r},t)$, measured with respect to a static ``reference'' lattice, are used to obtain accurate estimates of the elastic constants of model solids from atomistic…