Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 99
PROCEEDINGS OF THE ELEVENTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY Edited by: B.H.V. Topping
Paper 181
A Quadratic Ten Node Tetrahedral Cosserat Point Element for Nonlinear Elasticity M. Jabareen^{1}, E. Hanukah^{2} and M.B. Rubin^{2}
^{1}Faculty of Civil and Environmental Engineering, ^{2}Faculty of Mechanical Engineering,
M. Jabareen, E. Hanukah, M.B. Rubin, "A Quadratic Ten Node Tetrahedral Cosserat Point Element for Nonlinear Elasticity", in B.H.V. Topping, (Editor), "Proceedings of the Eleventh International Conference on Computational Structures Technology", CivilComp Press, Stirlingshire, UK, Paper 181, 2012. doi:10.4203/ccp.99.181
Keywords: Cosserat point element, hyperelasicity, quadratic, tetrahedral element.
Summary
It is well known that generating a mesh for a general threedimensional region is easier using tetrahedral elements than using brick elements. However, homogeneously deformable four node tetrahedral elements exhibit a stiff response to bending relative to eight node brick elements, which admit inhomogeneous deformations. Quadratic tetrahedral elements based on quadratic shape functions with full integration eliminate stiffness to bending modes but are known to exhibit inaccurate response for nearly incompressible materials. Mixed methods can be used for nearly incompressible materials, but they exhibit soft response to bending. Also, it is known that full integration elements exhibit poor response for contact problems with large deformations. In particular, the default quadratic tetrahedral element in the commercial code ABAQUS is an undocumented patented modified element. This modified element is indeed more robust than the full integration element but it is inaccurate in bending.
The objective of this work is to develop a quadratic ten node tetrahedral Cosserat point element (CPE) for nonlinear isotropic hyperelastic materials. In the standard BubnovGalerkin approach the kinematic approximation is assumed to be valid pointwise in the element and the constitutive equation is evaluated at the Gauss points of integration to determine element stiffnesses. In contrast, in the CPE approach [1,2] the kinematic approximation is only used to connect nodal director vectors (i.e. nodal positions) to element director vectors. Specifically, hyperelastic constitutive equations of the CPE are developed by treating the element as a structure with a strain energy function that models the response to all possible modes of deformation. The coefficients in the strain energy function for inhomogeneous deformation modes are determined by matching small deformation solutions for bending problems. Examples show that the resulting CPE exhibits accurate, robust response with a smooth transition from compressible to nearly incompressible material behavior. This new CPE can be used for thin structures and threedimensional bodies. It also outperforms the modified element of ABAQUS both in accuracy and robustness, even for contact problems with large deformations. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
