With the traditional fracture mechanics approaches, an initial crack and self-similar progression of cracks are assumed. In this treatise, theoretical and numerical tools are developed to mathematically describe non-self-similar progression of cracks without specifying an initial crack. A cohesive-decohesive zone model, similar to the cohesive zone model known in fracture mechanics as Dugdale-Barenblatt model, is adopted to represent the degradation of the material ahead of the crack tip. This model unifies strength-based crack initiation and fracture based crack progression. The cohesive-decohesive zone model is implemented with an interfacial surface material that consists of an upper and lower surface connected by a continuous distribution of normal and tangential nonlinear elastic springs that act to resist either Mode I opening, Mode II sliding, Mode III sliding, or mixed mode. The initiation of fracture is determined by the interfacial strength and the progression of fracture is determined by the critical energy release rate. The material between two adjacent laminae of a laminated composite structure or the material between the adherend and the adhesive is idealized with an interfacial surface material to predict interfacial fracture. The interfacial surface material is positioned within the bulk material to predict discrete cohesive cracks. The proper work-conjugacy relations between the stress and deformation measures are identified for the interfacial surface theory. In the principle of virtual work, the interfacial cohesive-decohesive tractions are conjugate to the displacement jumps across the upper and lower surfaces. A finite deformation kinematics theory is developed for the description of the upper and lower surface such that the deformation measures are invariant with respect to superposed rigid body translation and rotation. Various mechanical softening constitutive laws thermodynamically consistent with damage mechanics are postulated that relate the interfacial tractions to the displacement jump. An exponential function is used for the constitutive law such that it satisfies a multi-axial stress criterion for the onset of delamination, and satisfies a mixed mode fracture criterion for the progression of delamination. A damage parameter is included to prevent the restoration of the previous cohesive state between the interfacial surfaces. In addition, interfacial constitutive laws are developed to describe the contact-friction behavior. Interface elements applicable to two dimensional and three dimensional analyses are formulated for the analyses of contact, friction, and delamination problems. The consistent form of the interface element internal force vector and the tangent stiffness matrix are considered in the formulation. We investigate computational issues related to interfacial interpenetration, mesh sensitivity, the number of integrations points and the integration scheme, mathematical form of the softening constitutive law, and the convergence characteristics of the nonlinear solution procedure when cohesive-decohesive constitutive laws are used. To demonstrate the predictive capability of the interface finite element formulation, steadystate crack growth is simulated for quasi-static loading of various fracture test configurations loaded under Mode I, Mode II, Mode III, and mixed-mode loading. The finite element results are in agreement with the analytical results available in the literature and those developed in this work. A progressive failure methodology is developed and demonstrated to simulate the initiation and material degradation of a laminated panel due to intralaminar and interlaminar failures. Initiation of intralaminar failure can be by a matrix-cracking mode, a fiber-matrix shear mode, and a fiber failure mode. Subsequent material degradation is modeled using damage parameters for each mode to selectively reduce lamina material properties. The interlaminar failure mechanism such as delamination is simulated by positioning interface elements between adjacent sublaminates. The methodology is validated with respect to experimental data available in the literature on the response and failure of quasi-isotropic panels with centrally located circular cutouts. Very good agreement between the progressive failure analysis and the experiments is achieved if the failure analyses includes the interaction of intralaminar and interlaminar failures in the postbuckling response of the panels. In addition, ideas concerning the implementation of a fatigue model incorporated with a cohesive zone model are discussed. / Ph. D.

Using experimental data obtained from standard fracture test configurations, theoretical and numerical tools are developed to mathematically describe non-self-similar progression of cracks without specifying an initial crack. A cohesive-decohesive zone model, similar to the cohesive zone model known in the fracture mechanics literature as the Dugdale-Barenblatt model, is adopted to represent the degradation of the material ahead of the crack tip. This model unifies strength-based crack initiation and fracture-mechanics-based crack progression. The cohesive-decohesive zone model is implemented with an interfacial surface material that consists of an upper and a lower surface that are connected by a continuous distribution of normal and tangential nonlinear elastic springs that act to resist either Mode I opening, Mode II sliding, Mode III sliding, or a mixed mode. The initiation of fracture is determined by the interfacial strength and the progression of the crack is determined by the critical energy release rate. The adhesive is idealized with an interfacial surface material to predict interfacial fracture. The interfacial surface material is positioned within the bulk material to predict discrete cohesive cracks. The interfacial surface material is implemented through an interface element, which is incorporated in ABAQUS using the user defined element (UEL) option. A procedure is established to formulate a rate dependent model based on experiments carried out on compact tension test specimens. The rate dependent model is incorporated into the interface element approach to capture the unstable crack growth observed in experiments under quasi-static loading conditions. The compact tension test gives the variation of the fracture toughness with the rate of loading, this information is processed and a relationship between the fracture toughness and the rate of the opening displacement is established. The cohesive-decohesive zone model is implemented through a material model to be used in an explicit code (LS-DYNA). Dynamic simulations of the standard test configurations for Mode I (Double Cantilever Beam) and Mode II (End Load Split) are carried out using the explicit code. Verification of these coupon tests leads to the crash analysis of realistic structures like the square composite tube. Analyses of bonded and unbonded square tubes are presented. These tubes shows a very uncharacteristic failure mode: the composite material disintegrates on impact, and this has been captured in the analysis. Disadvantages of the interface element approach are well documented in the literature. An alternative method, known as the Extended Finite Element Method (XFEM), is implemented here through an eight-noded quadrilateral plane strain element. The method, based on the partition-of-unity, is used to study simple test configuration like the three-point bend problem and a double cantilever beam. Functionally graded materials are also simulated and the results are compared to the experimental results available in the literature. / Ph. D.