An algorithm for two-dimensional mesh generation based on the pinwheel tiling

We propose a new two-dimensional meshing algorithm called PINW able to generate meshes that accurately approximate the distance between any two domain points by paths composed only of cell edges. This technique is based on an extension of pinwheel tilings proposed by Radin and Conway. We prove that the algorithm produces triangles of bounded aspect ratio. This kind of mesh would be useful in cohesive interface finite element modeling when the crack propagation pathis an outcome of a simulation process.Comment: Short version appears in Proceedings of 2004 International Meshing Roundtable at http://www.imr.sandia.go

Time continuity in cohesive finite element modeling

We introduce the notion of time continuity for the analysis of cohesive zone interface finite element models. We focus on “initially rigid ” models in which an interface is inactive until the traction across it reaches a critical level. We argue that methods in this class are time discontinuous, unless special provision is made for the opposite. Time discontinuity leads to pitfalls in numerical implementations: oscillatory behavior, non-convergence in time and dependence on nonphysical regularization parameters. These problems arise at least partly from the attempt to extend uniaxial traction-displacement relationships to multiaxial loading. We also argue that any formulation of a time-continuous functional traction-displacement cohesive model entails encoding the value of the traction components at incipient softening into the model. We exhibit an example of such a model. Most of our numerical experiments concern explicit dynamics

Eulerian framework for inelasticity based on the Jaumann rate and a hyperelastic constitutive relation–part II: finite strain elastoplasticity

An integrable Eulerian rate formulation of finite deformation elasticity is developed, which relates the Jaumann or other objective corotational rate of the Kirchhoff stress with material spin to the same rate of the left Cauchy-Green deformation measure through a deformation dependent constitutive tensor. The proposed constitutive relationship can be written in terms of the rate of deformation tensor in the form of a hypoelastic material model. Integrability conditions, under which the proposed formulation yields (a) a Cauchy elastic and (b) a Green elastic material model are derived for the isotropic case. These determine the deformation dependent instantaneous elasticity tensor of the material. In particular, when the Cauchy integrability criterion is applied to the stressstrain relationship of a hyperelastic material model, an Eulerian rate formulation of hyperelasticity is obtained. This formulation proves crucial for the Eulerian finite strain elastoplastic model developed in part II of this work. The proposed model is formulated and integrated in the fixed background and extends the notion of an integrable hypoelastic model to arbitrary corotational objective rates and coordinates. Integrability was previously shown for the grade-zero hypoelastic model with use of the logarithmic (D) rate, the spin of which is formulated in principal coordinates. Uniform deformation examples of rectilinear shear, closed path four-step loading, and cyclic elliptical loading are presented. Contrary to classical grade-zero hypoelasticity, no shear oscillation, elastic dissipation, or ratcheting under cyclic load is observed when the simple Zaremba-Jaumann rate of stress is employed

An algorithm for two-dimensional mesh generation based on the pinwheel tiling ∗

We propose a new two-dimensional meshing algorithm called PINW able to generate meshes that accurately approximate the distance between any two domain points by paths composed only of cell edges. This technique is based on an extension of pinwheel tilings proposed by Radin and Conway. We prove that the algorithm produces triangles of bounded aspect ratio. This kind of mesh would be useful in cohesive interface finite element modeling when the crack propagation path is an outcome of a simulation process.