Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSubdivision Shell Elements with Anisotropic Growth
133
0
0.0
(
0
)
Added by
Roman Vetter
Publication date
2012
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
A thin shell finite element approach based on Loops subdivision surfaces is proposed, capable of dealing with large deformations and anisotropic growth. To this end, the Kirchhoff-Love theory of thin shells is derived and extended to allow for arbitrary in-plane growth. The simplicity and computational efficiency of the subdivision thin shell elements is outstanding, which is demonstrated on a few standard loading benchmarks. With this powerful tool at hand, we demonstrate the broad range of possible applications by numerical solution of several growth scenarios, ranging from the uniform growth of a sphere, to boundary instabilities induced by large anisotropic growth. Finally, it is shown that the problem of a slowly and uniformly growing sheet confined in a fixed hollow sphere is equivalent to the inverse process where a sheet of fixed size is slowly crumpled in a shrinking hollow sphere in the frictionless, quasi-static, elastic limit.
rate research
Read More
A method to simulate orthotropic behaviour in thin shell finite elements is proposed. The approach is based on the transformation of shape function derivatives, resulting in a new orthogonal basis aligned to a specified preferred direction for all elements. This transformation is carried out solely in the undeformed state leaving minimal additional impact on the computational effort expended to simulate orthotropic materials compared to isotropic, resulting in a straightforward and highly efficient implementation. This method is implemented for rotation-free triangular shells using the finite element framework built on the Kirchhoff--Love theory employing subdivision surfaces. The accuracy of this approach is demonstrated using the deformation of a pinched hemispherical shell (with a 18{deg} hole) standard benchmark. To showcase the efficiency of this implementation, the wrinkling of orthotropic sheets under shear displacement is analyzed. It is found that orthotropic subdivision shells are able to capture the wrinkling behavior of sheets accurately for coarse meshes without the use of an additional wrinkling model.
A class of peridynamic material models known as constitutive correspondence models provide a bridge between classical continuum mechanics and peridynamics. These models are useful because they allow well-established local constitutive theories to be used within the nonlocal framework of peridynamics. A recent finite deformation correspondence theory (Foster and Xu, 2018) was developed and reported to improve stability properties of the original correspondence model (Silling et al., 2007). This paper presents a stability analysis that indicates the reported advantages of the new theory were overestimated. Homogeneous deformations are analyzed and shown to exibit unstable material behavior at the continuum level. Additionally, the effects of a particle discretization on the stability of the model are reported. Numerical examples demonstrate the large errors induced by the unstable behavior. Stabilization strategies and practical applications of the new finite deformation model are discussed.
In this paper, we present a family of multivariate grid transfer operators appropriate for anisotropic multigrid methods. Our grid transfer operators are derived from a new family of anisotropic interpolatory subdivision schemes. We study the minimality, polynomial reproduction and convergence properties of these interpolatory schemes and link their properties to the convergence and optimality of the corresponding multigrid methods. We compare the performance of our interpolarory grid transfer operators with the ones derived from a family of corresponding approximating subdivision schemes.
This paper addresses some numerical and theoretical aspects of dual Schur domain decomposition methods for linear first-order transient partial differential equations. In this work, we consider the trapezoidal family of schemes for integrating the ordinary differential equations (ODEs) for each subdomain and present four different coupling methods, corresponding to different algebraic constraints, for enforcing kinematic continuity on the interface between the subdomains. Method 1 (d-continuity) is based on the conventional approach using continuity of the primary variable and we show that this method is unstable for a lot of commonly used time integrators including the mid-point rule. To alleviate this difficulty, we propose a new Method 2 (Modified d-continuity) and prove its stability for coupling all time integrators in the trapezoidal family (except the forward Euler). Method 3 (v-continuity) is based on enforcing the continuity of the time derivative of the primary variable. However, this constraint introduces a drift in the primary variable on the interface. We present Method 4 (Baumgarte stabilized) which uses Baumgarte stabilization to limit this drift and we derive bounds for the stabilization parameter to ensure stability. Our stability analysis is based on the ``energy method, and one of the main contributions of this paper is the extension of the energy method (which was previously introduced in the context of numerical methods for ODEs) to assess the stability of numerical formulations for index-2 differential-algebraic equations (DAEs).
Objective: To present the first real-time a posteriori error-driven adaptive finite element approach for real-time simulation and to demonstrate the method on a needle insertion problem. Methods: We use corotational elasticity and a frictional needle/tissue interaction model. The problem is solved using finite elements within SOFA. The refinement strategy relies upon a hexahedron-based finite element method, combined with a posteriori error estimation driven local $h$-refinement, for simulating soft tissue deformation. Results: We control the local and global error level in the mechanical fields (e.g. displacement or stresses) during the simulation. We show the convergence of the algorithm on academic examples, and demonstrate its practical usability on a percutaneous procedure involving needle insertion in a liver. For the latter case, we compare the force displacement curves obtained from the proposed adaptive algorithm with that obtained from a uniform refinement approach. Conclusions: Error control guarantees that a tolerable error level is not exceeded during the simulations. Local mesh refinement accelerates simulations. Significance: Our work provides a first step to discriminate between discretization error and modeling error by providing a robust quantification of discretization error during simulations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
