No abstract

SUMMARYDynamic crack microbranching processes in brittle materials are investigated by means of a computational fracture mechanics approach using the finite element method with special interface elements and a topological data structure representation. Experiments indicate presence of a limiting crack speed for dynamic crack in brittle materials as well as increasing fracture resistance with crack speed. These phenomena are numerically investigated by means of a cohesive zone model (CZM) to characterize the fracture process. A critical evaluation of intrinsic versus extrinsic CZMs is briefly presented, which highlights the necessity of adopting an extrinsic approach in the current analysis. A novel topologybased data structure is employed to enable fast and robust manipulation of evolving mesh information when extrinsic cohesive elements are inserted adaptively. Compared to intrinsic CZMs, which include an initial hardening segment in the traction-separation curve, extrinsic CZMs involve additional issues both in implementing the procedure and in interpreting simulation results. These include time discontinuity in stress history, fracture pattern dependence on time step control, and numerical energy balance. These issues are investigated in detail through a 'quasi-steady-state' crack propagation problem in polymethylmethacrylate. The simulation results compare reasonably well with experimental observations both globally and locally, and demonstrate certain advantageous features of the extrinsic CZM with respect to the intrinsic CZM.

SUMMARYThis paper presents a novel compact adjacency-based topological data structure for finite element mesh representation. The proposed data structure is designed to support, under the same framework, both two-and three-dimensional meshes, with any type of elements defined by templates of ordered nodes. When compared to other proposals, our data structure reduces the required storage space while being 'complete', in the sense that it preserves the ability to retrieve all topological adjacency relationships in constant time or in time proportional to the number of retrieved entities. Element and node are the only entities explicitly represented. Other topological entities, which include facet, edge, and vertex, are implicitly represented. In order to simplify accessing topological adjacency relationships, we also define and implicitly represent oriented entities, associated to the use of facets, edges, and vertices by an element. All implicit entities are represented by concrete types, being handled as values, which avoid usual problems encountered in other reduced data structures when performing operations such as entity enumeration and attribute attachment. We also extend the data structure with the use of 'reverse indices', which improves performance for extracting adjacency relationships while maintaining storage space within reasonable limits. The data structure effectiveness is demonstrated by two different applications: for supporting fragmentation simulation and for supporting volume rendering algorithms.

SUMMARYDependence on mesh orientation impacts adversely the quality of computational solutions generated by cohesive zone models. For instance, when considering crack propagation along interfaces between finite elements of 4k structured meshes, both extension of crack length and crack angle are biased according to the mesh configuration. To address mesh orientation dependence in 4k structured meshes and to avoid undesirable crack patterns, we propose the use of nodal perturbation (NP) and edge-swap (ES) topological operation. To this effect, the topological data structure TopS (Int. J. Numer. Meth. Engng 2005; 64: 1529-1556), based on topological entities (node, element, vertex, edge and facet), is utilized so that it is possible to access adjacency information and to manage a consistent data structure in time proportional to the number of retrieved entities. In particular, the data structure allows the ES operation to be done in constant time. Three representative dynamic fracture examples using ES and NP operators are provided: crack propagation in the compact compression specimen, local branching instability, and fragmentation. These examples illustrate the features of the present computational framework in simulating a range of physical phenomena associated with cracking.

No abstract