Finite Element Simulation of a Crack Growth in the Presence of a Hole in the Vicinity of the Crack Trajectory

The aim of this paper was to present a numerical simulation of a crack growth path and associated stress intensity factors (SIFs) for linear elastic material. The influence of the holes’ position and pre-crack locations in the crack growth direction were investigated. For this purpose, ANSYS Mechanical R19.2 was introduced with the use of a new feature known as Separating Morphing and Adaptive Remeshing Technology (SMART) dependent on the Unstructured Mesh Method (UMM), which can reduce the meshing time from up to several days to a few minutes, eliminating long preprocessing sessions. The presence of a hole near a propagating crack causes a deviation in the crack path. If the hole is close enough to the crack path, the crack may stop at the edge of the hole, resulting in crack arrest. The present study was carried out for two geometries, namely a cracked plate with four holes and a plate with a circular hole, and an edge crack with different pre-crack locations. Under linear elastic fracture mechanics (LEFM), the maximum circumferential stress criterion is applied as a direction criterion. Depending on the position of the hole, the results reveal that the crack propagates in the direction of the hole due to the uneven stresses at the crack tip, which are consequences of the hole’s influence. The results of this modeling are validated in terms of crack growth trajectories and SIFs by several crack growth studies reported in the literature that show trustworthy results.

Application of Failure Criteria on Plywood under Bending

In composite materials, the use of failure criteria is necessary to determine the failure forces. Various failure criteria are known, from the simplest ones that compare individual stresses with the corresponding strength, to more complex ones that take into account the sign and direction of the stress, as well as mutual interactions of the acting stresses. This study investigates the application of the maximum stress, Tsai-Hill, Tsai-Wu, Puck, Hoffman and Hashin criteria to beech plywood made from a series of plies of differently oriented beech veneers. Specimens were cut from the manufactured boards at various angles and loaded by bending to failure. The mechanical properties of the beech veneer were also determined. The specimens were modelled using the finite element method with a composite modulus and considering the different failure criteria where the failure forces were calculated and compared with the measured values. It was found that the calculated forces based on all failure criteria were lower than those measured experimentally. The forces determined using the maximum stress criterion showed the best agreement between the calculated and measured forces.

Mixed Mode Crack Propagation in Iliac Bone

Bone is a complex material that can be regarded as an anisotropic elastic composite material. The problem of crack propagation in human bone is analyzed by using a generalization of the maximum tensile stress criterion (MTS). The results concern the critical stress for crack propagation and the direction of the crack path in Iliac bone.

Open hole size effects on tensile properties of 3D braided composites

Owing to the excellent integrated structure, notch-insensitivity, delamination-free characteristics, 3D braided composites have a broad range of engineering applications. In this paper, the notch size effects on two types of 3D braided composites were experimentally examined. Style I incorporated 40% of longitudinal lay-in yarns. Style II was the pure braids. The Point Stress Criterion (PSC) was applied to predict the open-hole strength of 3D braided composites. It is found the 3D braided composites can keep higher proportion residual strength after involving the different circular hole sizes compared to plain woven laminates. The open-hole pure braided specimen shows better performance than that the braids with longitudinal yarns, the lay-in longitudinal yarns improve neither specimens’ un-notched strength, nor the modulus. The predicted open-hole strength were compared with experimental results. The traditional analytical method can predict the open-hole strength of 3D braided composite to some extent. Under uniaxial tensile stress, the failure behaviour of two types of 3D braided specimens are different. For un-notched specimen, clear cracks usually show up on the Style II specimen, while it is not true for Style I coupon. For notched specimen, the crack of both notched specimens will propagate along the notch and finally render the specimen to fail

Modeling of Fatigue Wear of Viscoelastic Coatings

A new model for studying the kinetics of fatigue wear of a viscoelastic coating bonded to a rigid substrate is proposed. The fatigue mechanism is due to the cyclic interaction of the coating with a rough counterbody, which is modeled by a periodic system of smooth indenters. The study includes the solution of the problem of sliding contact of the indenter at a constant velocity along the viscoelastic coating, the calculation of stresses taking into account the mutual effect, and study of the process of damage accumulation in the material. The calculation of the damage function of the surface layer was carried out using the reduced stress criterion. Assuming the possibility of summation of accumulated damage, two processes were considered: delamination of surface layers of the coating and continuous fracture of the surface by the fatigue mechanism. The effect of the sliding velocity and viscoelastic properties of the material on the damage accumulation and the coating wear rate was analyzed. Two types of load, constant and stochastically varying, were used in modeling and analysis. It was found that the rate of fatigue wear of the coating increased and then became constant.

Deformation and Damage Assessments of Two DP1000 Steels Using a Micromechanical Modelling Method

Damage characterization and micromechanical modelling in dual-phase (DP) steels have recently drawn attention, since any changes in the alloying elements or process route strongly influence the microstructural features, deformation behavior of the phases, and damage to the micro-mechanisms, and subsequently the particular mechanical properties of the material. This approach can be used to stablish microstructure–properties relationships. For instance, the effects of local damage from shear cutting on edge crack sensitivity in the following deformation process can be studied. This work evaluated the deformation and damage behaviors of two DP1000 steels using a microstructure-based approach to estimate the edge cracking resistance. Phase fraction, grain size, phase distribution, and texture were analyzed using electron backscatter diffraction and secondary electron detectors of a scanning electron microscope and employed in 3D representative volume elements. The deformation behavior of the ferrite phase was defined using a crystal plasticity model, which was calibrated through nanoindentation tests. Various loading conditions, including uniaxial tension, equi-biaxial tension, plane strain tension, and shearing, along with the maximum shear stress criterion were applied to investigate the damage initiation and describe the edge cracking sensitivity of the studied steels. The results revealed that a homogenous microstructure leads to homogenous stress–strain partitioning, delayed damage initiation, and high edge cracking resistance.

Energizing ASTM lap joint fracture standards

Several ASTM standards on the fracture of glued and welded joints need attention because they do not consider the Griffith energy criterion of cracking which was proposed a century ago. It is almost as if Griffith never existed because the ASTM definition of failure is the stress criterion postulated by Galileo in 1638 in which stress at failure (i.e. strength = force/area) is defined as the determinant of fracture. Irene Martinez Villegas (Villegas, Rans 2021 Phil. Trans. R. Soc. A 376, 20200296. ( doi:10.1098/rsta.2020.0296 )) shows in this volume that attempts to use ASTM D5868 to standardize welded composite (carbon fibre reinforced polymer, CFRP) lap joints reveal major problems. First, the test is a low angle bend–peel test; not shear. Second, the energy required to break the joint is not emphasized so that joints may have high strength properties but also low toughness; third, the fracture force is not proportional to the lap joint area so the concept of strength independent of sample size is false; fourth, as the CFRP panels are made thicker, the strength rises at constant overlap area so the strength can be any value you want; fifth, the strength of larger joints goes down; this is the size effect noted in many bend-cracking tests, much as Galileo suggested for bent beam fracture in his famous book ‘the larger the machine, the greater its weakness'. The purpose of this paper is to demonstrate that poor ASTM ‘shear strength’ standards should be replaced by a definition of welded lap joint performance based on Griffith's energy conservation argument in which fracture surface energy is the main parameter resisting failure. The foundation of this Griffith-style lap joint analysis for long cracks goes back to 1975 but has been largely ignored until now because it does not fit the Griffith equation for cracked sheets, has no ‘stress intensity factor’, and travels at constant speed, not accelerating like the standard Griffith tension crack. This study of tensile delamination shows that a long lap crack is not driven by stress near the crack but by changes in stored elastic energy in the stretched strips remote from the crack tip, while strain energy release rate is negative. It would be more appropriate to call this lap failure a tensile delamination crack. This article is part of a discussion meeting issue ‘A cracking approach to inventing new tough materials: fracture stranger than friction’.

Comparison of Various Criteria Determining the Direction of Crack Propagation Using the UDMGINI User Procedure Implemented in Abaqus

This paper describes a method of predicting the direction of crack propagation implemented by user subroutines in the Simulia-Abaqus FEA system with the use of the extended finite element method (X-FEM). This method is based on displacements and stresses according to Westergaard’s solution of Griffith’s crack problem. During the calculations, in each crack increment, the algorithm reads the stresses and displacements in the model around the crack tip, calculates the criterion values at the read points, reduces them to a unit distance from the crack tip, fits a polynomial to these points, and finds the minimum of the function closest to the last propagation angle. The algorithm also decides when the crack grows, depending on a chosen criterion. Four criteria have been implemented to predict the direction of failure propagation: the maximum principal stress criterion, the Ottosen–Podgórski criterion, the new criterion described here based on the minimum component values of the displacement vector, and the maximum circumferential tensile stress (MTS). These criteria were verified in two tests: the three-point bending test of the notched beam and the anchor pull-out test. For these tests, the criterion built into Simulia Abaqus does not correctly define the crack path, which causes the crack propagation direction to “rotate” when simulating the fracture. The criteria developed here, in most cases, determine the crack path and the maximum force very well compared to real laboratory tests.

Extension of a basic hypoplastic model for overconsolidated clays

This paper presents a new rate-dependent hypoplastic constitutive model for overconsolidated clays. The model is developed based on a basic hypoplastic model proposed recently for sand. New density and stiffness factors are introduced to account for history dependence. The Matsuoka-Nakai failure surface is incorporated for the limit stress criterion. With six constitutive parameters, the model is capable of predicting the hardening/softening, shear dilation/contraction, and asymptotic state for overconsolidated clays. Comparison between numerical predictions and experimental results shows this model can properly describe the main features of both reconstituted and undisturbed clays with different overconsolidation ratios.