A theoretical comparison between inner products in the shift-invert Arnoldi method and the spectral transformation Lanczos method.
A three dimensional finite element method for biological active soft tissue formulation in cylindrical polar coordinates
A hyperelastic constitutive law, for use in anatomically accurate finite element models of living structures, is suggested for the passive and the active mechanical properties of incompressible biological tissues. This law considers the passive and active states as a same hyperelastic continuum medium, and uses an activation function in order to describe the whole contraction phase. The variational and the FE formulations are also presented, and the FE code has been validated and applied to describe...
A three dimensional finite element method for biological active soft tissue Formulation in cylindrical polar coordinates
A hyperelastic constitutive law, for use in anatomically accurate finite element models of living structures, is suggested for the passive and the active mechanical properties of incompressible biological tissues. This law considers the passive and active states as a same hyperelastic continuum medium, and uses an activation function in order to describe the whole contraction phase. The variational and the FE formulations are also presented, and the FE code has been validated and applied to describe...
A three-field augmented Lagrangian formulation of unilateral contact problems with cohesive forces
We investigate unilateral contact problems with cohesive forces, leading to the constrained minimization of a possibly nonconvex functional. We analyze the mathematical structure of the minimization problem. The problem is reformulated in terms of a three-field augmented Lagrangian, and sufficient conditions for the existence of a local saddle-point are derived. Then, we derive and analyze mixed finite element approximations to the stationarity conditions of the three-field augmented Lagrangian....
A variationally consistent generalized variable formulation of the elastoplastic rate problem
The elastoplastic rate problem is formulated as an unconstrained saddle point problem which, in turn, is obtained by the Lagrange multiplier method from a kinematic minimum principle. The finite element discretization and the enforcement of the min-max conditions for the Lagrangean function lead to a set of algebraic governing relations (equilibrium, compatibility and constitutive law). It is shown how important properties of the continuum problem (like, e.g., symmetry, convexity, normality) carry...
Adaptive finite element relaxation schemes for hyperbolic conservation laws
We propose and study semidiscrete and fully discrete finite element schemes based on appropriate relaxation models for systems of Hyperbolic Conservation Laws. These schemes are using piecewise polynomials of arbitrary degree and their consistency error is of high order. The methods are combined with an adaptive strategy that yields fine mesh in shock regions and coarser mesh in the smooth parts of the solution. The computational performance of these methods is demonstrated by considering scalar...
Adaptive Finite Element Relaxation Schemes for Hyperbolic Conservation Laws
We propose and study semidiscrete and fully discrete finite element schemes based on appropriate relaxation models for systems of Hyperbolic Conservation Laws. These schemes are using piecewise polynomials of arbitrary degree and their consistency error is of high order. The methods are combined with an adaptive strategy that yields fine mesh in shock regions and coarser mesh in the smooth parts of the solution. The computational performance of these methods is demonstrated by considering scalar...
Advancing analysis capabilities in ANSYS through solver technology.
ALBERT---Software for scientific computations and applications.
Algebraic domain decomposition solver for linear elasticity
We generalize the overlapping Schwarz domain decomposition method to problems of linear elasticity. The convergence rate independent of the mesh size, coarse-space size, Korn’s constant and essential boundary conditions is proved here. Abstract convergence bounds developed here can be used for an analysis of the method applied to singular perturbations of other elliptic problems.
An a posteriori error analysis for dynamic viscoelastic problems
In this paper, a dynamic viscoelastic problem is numerically studied. The variational problem is written in terms of the velocity field and it leads to a parabolic linear variational equation. A fully discrete scheme is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize time derivatives. An a priori error estimates result is recalled, from which the linear convergence is derived under suitable regularity conditions. Then, an a posteriori...
An a posteriori error analysis for dynamic viscoelastic problems
In this paper, a dynamic viscoelastic problem is numerically studied. The variational problem is written in terms of the velocity field and it leads to a parabolic linear variational equation. A fully discrete scheme is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize time derivatives. An a priori error estimates result is recalled, from which the linear convergence is derived under suitable regularity conditions. Then, an a posteriori error...
An adaptive finite element method for solving a double well problem describing crystalline microstructure
An adaptive finite element method for solving a double well problem describing crystalline microstructure
The minimization of nonconvex functionals naturally arises in materials sciences where deformation gradients in certain alloys exhibit microstructures. For example, minimizing sequences of the nonconvex Ericksen-James energy can be associated with deformations in martensitic materials that are observed in experiments[2,3]. — From the numerical point of view, classical conforming and nonconforming finite element discretizations have been observed to give minimizers with their quality being highly dependent...
An additive Schwarz method for mortar Morley finite element discretizations of 4th order elliptic problem in 2D.
An Algorithm to Compute a Sparse Basis of the Null Space.
An analysis of a contact problem for a cylindrical shell: A primary and dual formulation
In this paper the contact problem for a cylindrical shell and a stiff punch is studied. The existence and uniqueness of a solution is verified. The finite element method is discussed.
An analysis of the boundary layer in the 1D surface Cauchy–Born model
The surface Cauchy–Born (SCB) method is a computational multi-scale method for the simulation of surface-dominated crystalline materials. We present an error analysis of the SCB method, focused on the role of surface relaxation. In a linearized 1D model we show that the error committed by the SCB method is 𝒪(1) in the mesh size; however, we are able to identify an alternative “approximation parameter” – the stiffness of the interaction potential – with respect to which the relative error...
An analysis of the boundary layer in the 1D surface Cauchy–Born model∗
The surface Cauchy–Born (SCB) method is a computational multi-scale method for the simulation of surface-dominated crystalline materials. We present an error analysis of the SCB method, focused on the role of surface relaxation. In a linearized 1D model we show that the error committed by the SCB method is 𝒪(1) in the mesh size; however, we are able to identify an alternative “approximation parameter” – the stiffness of the interaction potential – with respect to which the relative error...
An Analysis of the p-Version of the Finite Element Method for Nearly Incompressible Materials. Uniformly Valid, Optimal Error Estimates