Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

Fatigue Life Methodology for Bonded Composite Skin/Stringer Configurations

A methodology is presented for determining the fatigue life of composite structures based on fatigue characterization data and geometric nonlinear finite element (FE) analyses. To demonstrate the approach, predicted results were compared to fatigue tests performed on specimens which represented a tapered composite flange bonded onto a composite skin. In a first step, tension tests were performed to evaluate the debonding mechanisms between the flange and the skin. In a second step, a 2D FE model was developed to analyze the tests. To predict matrix cracking onset, the relationship between the tension load and the maximum principal stresses transverse to the fiber direction was determined through FE analysis. Transverse tension fatigue life data were used to -enerate an onset fatigue life P-N curve for matrix cracking. The resulting prediction was in good agreement with data from the fatigue tests. In a third step, a fracture mechanics approach based on FE analysis was used to determine the relationship between the tension load and the critical energy release rate. Mixed mode energy release rate fatigue life data were used to create a fatigue life onset G-N curve for delamination. The resulting prediction was in good agreement with data from the fatigue tests. Further, the prediction curve for cumulative life to failure was generated from the previous onset fatigue life curves. The results showed that the methodology offers a significant potential to Predict cumulative fatigue life of composite structures.

Krueger, Ronald↗

Development and Performance Evaluations of HfO2-Si and Rare Earth-Si Based Environmental Barrier Bond Coat Systems for SiC/SiC Ceramic Matrix Composites

Ceramic environmental barrier coatings (EBC) and SiCSiC ceramic matrix composites (CMCs) will play a crucial role in future aircraft propulsion systems because of their ability to significantly increase engine operating temperatures, improve component durability, reduce engine weight and cooling requirements. Advanced EBC systems for SiCSiC CMC turbine and combustor hot section components are currently being developed to meet future turbine engine emission and performance goals. One of the significant material development challenges for the high temperature CMC components is to develop prime-reliant, high strength and high temperature capable environmental barrier coating bond coat systems, since the current silicon bond coat cannot meet the advanced EBC-CMC temperature and stability requirements. In this paper, advanced NASA HfO2-Si based EBC bond coat systems for SiCSiC CMC combustor and turbine airfoil applications are investigated. The coating design approach and stability requirements are specifically emphasized, with the development and implementation focusing on Plasma Sprayed (PS) and Electron Beam-Physic Vapor Deposited (EB-PVD) coating systems and the composition optimizations. High temperature properties of the HfO2-Si based bond coat systems, including the strength, fracture toughness, creep resistance, and oxidation resistance were evaluated in the temperature range of 1200 to 1500 C. Thermal gradient heat flux low cycle fatigue and furnace cyclic oxidation durability tests were also performed at temperatures up to 1500 C. The coating strength improvements, degradation and failure modes of the environmental barrier coating bond coat systems on SiCSiC CMCs tested in simulated stress-environment interactions are briefly discussed and supported by modeling. The performance enhancements of the HfO2-Si bond coat systems with rare earth element dopants and rare earth-silicon based bond coats are also highlighted. The advanced bond coat systems, when integrated with advanced EBC top coats, showed promise to achieve 1500 C temperature capability, helping enable next generation turbine engines with significantly improved engine component temperature capability and long-term durability.

HfO2-Si bond coat↗

The influence of temperature on fatigue-crack growth in a mill-annealed Ti-6Al-4V alloy

To understand the influence of temperature on the rate of fatigue crack growth in high strength metal alloys, constant load amplitude fatigue crack growth experiments were carried out using a 1/4 inch thick (6.35 mm) mill-annealed Ti-6Al-4V alloy plate as a model material. The rates of fatigue crack growth were determined as a function of temperature, ranging from room temperature to about 290 C and as a function of the crack tip, stress intensity factor K, in dehumidified high purity argon environment. The dependence of the rate of fatigue crack growth on K appears to be separable into two regions. The transition correlates with changes in both the microscopic and macroscopic appearances of the fracture surfaces, and suggests a change in the mechanism and the influence of microstructure on fatigue crack growth.

Wei, R. P.↗

Discrete Spring Model for Predicting Delamination Growth in Z-Fiber Reinforced DCB Specimens

Beam theory analysis was applied to predict delamination growth in DCB specimens reinforced in the thickness direction with pultruded pins, known as Z-fibers. The specimen arms were modeled as cantilever beams supported by discrete springs, which were included to represent the pins. A bi-linear, irreversible damage law was used to represent Z-fiber damage, the parameters of which were obtained from previous experiments. Closed-form solutions were developed for specimen compliance and displacements corresponding to Z-fiber row locations. A solution strategy was formulated to predict delamination growth, in which the parent laminate mode I fracture toughness was used as the criterion for delamination growth. The solution procedure was coded into FORTRAN 90, giving a dedicated software tool for performing the delamination prediction. Comparison of analysis results with previous analysis and experiment showed good agreement, yielding an initial verification for the analytical procedure.

Ratcliffe, James G.↗

Irregular Mare Patches (IMPs): 100 Ma or 3 Ga?

IMPs exhibit a perplexing combination of characteristics that are consistent with either an approximately 100 Ma or 3 Ga formation. Dozens of small-area IMPs have crisp morphologies and crater size-frequency distributions (SFDs) that denote relatively recent geologic activity (less than 100 Ma); however, the apparently well-developed regolith on portions of the IMPs are in conflict with such a young age [1]. To test possible formation hypotheses (e.g., [1-5]), which range from ancient volcanism to contemporary outgassing, we examined IMP morphology at the meter-scale with LROC NAC images and derived elevation models. We focused on the largest IMPs (Ina, Sosigenes, Cauchy, Maskelyne, and Nubium), where contacts between deposits are best developed. Most of our observations are consistent with multiple generations of inflation and breakouts (or squeeze-ups) of basaltic lavas that were affected by local slopes. Some of the extrusions coalesced into larger mounds or filled pre-existing craters. We did not observe evidence of large-scale void space (e.g., fissures, fractures, linear depressions, or pits) within or beneath the mounds or rougher deposits (e.g., [5]). But, small-scale voids may be signified by isolated pitted textures. We also did not detect evidence of the cooling fractures or lava plates expected in young lava flows and observed in lunar impact melt deposits. The smooth texture of the mounds is enigmatic. Block-less craters suggest at least 5 m of friable or poorly-cohesive material (such as regolith), yet mound margins exhibit slopes greater than 30 deg requiring significant material strength. Blocks are not common on the mounds, but are sometimes excavated by impacts (usually excavated from beneath the mounds). The uneven deposits are equally enigmatic and texturally varied (blocky, pitted, and crenulated). They are deficient in superposed craters compared to the mounds. If the mounds are indeed of similar age to the rougher units, then their different superposed crater morphologies and SFDs need to be explained by factors other than their ages. Any mounds originally composed of friable surface materials would evolve differently from more coherent deposits (e.g., [6-7]).

Stopar, Julie↗

Rate and time dependent failure of structural adhesives

Studies on two adhesives (Metlbond 1113 and 1113-2) identified as having important applications in the bonding of composite materials are presented. A testing program to ascertain stress-strain, strain-rate, time, yield, and/or failure behavior of these materials in bulk form using uniaxial tensile constant strain-rate, creep, and relaxation tests is described. The stress-strain behavior of each material is shown to be significantly rate dependent. A rate dependent stress whitening (crazing) phenomenon occurs prior to either yield or fracture. A region of linear elasticity, a region of viscoelasticity, and the onset of yielding are identified in the stress-strain behavior. The linear elastic limit and the yield point are shown to be rate dependent and agree well with an empirical equation proposed by Ludwik. A creep to failure phenomenon is shown to exist and is correlated with a delayed yield equation proposed by Crochet. Analytical predictions based on a modified Bingham model are shown to agree well with experimental stress-strain strain-rate data. Analytical predictions based on a modified Ramberg-Osgood equation are also shown for comparison purposes.

Brinson, H. F.↗

Plates and shells containing a surface crack under general loading conditions

Various through and part-through crack problems in plates and shells are considered. The line-spring model of Rice and Levy is generalized to the skew-symmetric case to solve surface crack problems involving mixed-mode, coplanar crack growth. Compliance functions are introduced which are valid for crack depth to thickness ratios at least up to .95. This includes expressions for tension and bending as well as expressions for in-plane shear, out-of-plane shear, and twisting. Transverse shear deformation is taken into account in the plate and shell theories and this effect is shown to be important in comparing stress intensity factors obtained from the plate theory with three-dimensional solutions. Stress intensity factors for cylinders obtained by the line-spring model also compare well with three-dimensional solution. By using the line-spring approach, stress intensity factors can be obtained for the through crack and for part-through crack of any crack front shape, without recalculation integrals that take up the bulk of the computer time. Therefore, parameter studies involving crack length, crack depth, shell type, and shell curvature are made in some detail. The results will be useful in brittle fracture and in fatigue crack propagation studies. All problems considered are of the mixed boundary value type and are reducted to strongly singular integral equations which make use of the finite-part integrals of Hadamard. The equations are solved numerically in a manner that is very efficient.

Joseph, Paul F.↗

Characterization of Low-Velocity Impact Damage in Toughened Non-Crimp Fabric Laminates

The damage response of resin-infused composites was investigated using quasi-static indentation (QSI) and instrumented drop-weight impact testing. The test specimens were composed of six layers of non-crimp fabrics where each layer contained +45°, 90°, −45°, and 0° plies. The goals of the testing were to identify damage responses that may differ from traditional thermoset materials and provide damage characterization data to evaluate current damage prediction analyses and aid improvements to current damage models. The damage was characterized using X-ray computed tomography (CT) and ultrasonic testing (UT). Incremental QSI loading was used on a single specimen to investigate damage growth in nine steps at progressively higher loads. Realtime acoustic emission data were used to determine the unloading points. Unlike typical thermoset composites, the force-displacement was smooth and showed no obvious indications of damage. The primary damage consisted of large lines of fiber fracture on the plies near the impact surface and delaminations through the thickness. Nine specimens were impacted at nominal impact energies of 25 ft-lbs, 30 ft-lbs, and 35 ft-lbs. The force histories at each energy followed nearly identical paths. Again, no obvious signs of damage were observed in any of the force histories. The general damage response in the dynamic tests was similar to the QSI-loaded specimens. However, repeat specimens impacted with a given impact energy showed variations in the damage. Some characteristics of the damage response were different than those typically observed in thermoset composites and will need to be evaluated using current analysis methods.

Composite Materials↗

Characterization of Low-Velocity Impact Damage in Toughened Non-Crimp Fabric Laminates

The damage response of resin-infused composites was investigated using quasi-static indentation (QSI) and instrumented drop-weight impact testing. The test specimens were composed of six layers of non-crimp fabrics where each layer contained +45°, 90°, −45°, and 0° plies. The goals of the testing were to identify damage responses that may differ from traditional thermoset materials and provide damage characterization data to evaluate current damage prediction analyses and aid improvements to current damage models. The damage was characterized using X-ray computed tomography (CT) and ultrasonic testing (UT). Incremental QSI loading was used on a single specimen to investigate damage growth in nine steps at progressively higher loads. Realtime acoustic emission data were used to determine the unloading points. Unlike typical thermoset composites, the force-displacement was smooth and showed no obvious indications of damage. The primary damage consisted of large lines of fiber fracture on the plies near the impact surface and delaminations through the thickness. Nine specimens were impacted at nominal impact energies of 25 ft-lbs, 30 ft-lbs, and 35 ft-lbs. The force histories at each energy followed nearly identical paths. Again, no obvious signs of damage were observed in any of the force histories. The general damage response in the dynamic tests was similar to the QSI-loaded specimens. However, repeat specimens impacted with a given impact energy showed variations in the damage. Some characteristics of the damage response were different than those typically observed in thermoset composites and will need to be evaluated using current analysis methods.

Composite Materials↗

NASCRAC - A computer code for fracture mechanics analysis of crack growth

NASCRAC - a computer code for fracture mechanics analysis of crack growth - is described in this paper. The need for such a code is increasing as requirements grow for high reliability and low weight in aerospace components. The code is comprehensive and versatile, as well as user friendly. The major purpose of the code is calculation of fatigue, corrosion fatigue, or stress corrosion crack growth, and a variety of crack growth relations can be selected by the user. Additionally, crack retardation models are included. A very wide variety of stress intensity factor solutions are contained in the code, and extensive use is made of influence functions. This allows complex stress gradients in three-dimensional crack problems to be treated easily and economically. In cases where previous stress intensity factor solutions are not adequate, new influence functions can be calculated by the code. Additional features include incorporation of J-integral solutions from the literature and a capability for estimating elastic-plastic stress redistribution from the results of a corresponding elastic analysis. An example problem is presented which shows typical outputs from the code.

Harris, D. O.↗

Evidence for Hesperian Impact-Induced Hydrothermalism on Mars

Several hydrated silicate deposits on Mars are observed within craters and are interpreted as excavated Noachian material. Toro crater (71.8 deg E, 17.0 deg N), located on the northern edge of the Syrtis Major Volcanic Plains, shows spectral and morphologic evidence of impact-induced hydrothermal activity. Spectroscopic observations were used to identify extensive hydrated silicate deposits, including prehnite, chlorites, smectites, and opaline material, a suite of phases that frequently results from hydrothermal alteration in terrestrial craters and also expected on Mars from geochemical modeling of hydrothermal environments. When combined with altimetry and high-resolution imaging data, these deposits appear associated predominantly with the central uplift and with portions of the northern part of the crater floor. Detailed geologic mapping of these deposits reveals geomorphic features that are consistent with hydrothermal activity that followed the impact event, including vent-like and conical mound structures, and a complex network of tectonic structures caused by fluid interactions such as fractures and joints. The crater age has been calculated from the cumulative crater size-frequency distributions and is found to be Early Hesperian. The evidence presented here provides support for impact-induced hydrothermal activity in Toro crater, that extends phyllosilicate formation processes beyond the Noachian era.

Marzo, Giuseppe A.↗

The Komplast Experiment: Space Environmental Effects after 12 Years in LEO (and Counting)

The Komplast materials experiment was designed by the Khrunichev State Research and Production Space Center, together with other Russian scientific institutes, and has been carried out by Mission Control Moscow since 18. Komplast panels fitted with material samples and sensors were located on the International Space Station (ISS) Functional Cargo Block (FGB) module exterior surface. Within the framework of this experiment, the purpose was to study the effect of the low earth orbit (LEO) environment on exposed samples of various materials. The panels were sent into orbit with the FGB when it launched on November 20, 1998. Panels #2 and #10 were retrieved during Russian extravehicular activity in February 2011 and sealed within cases to temporarily protect the samples from exposure to air until they could be studied on the ground. Panel #2 contained an experiment to detect micrometeoroid and orbital debris (MMOD) impacts, radiation and UV sensors, several pieces of electrical cable, and samples made from elastomeric and fluoroplastic materials. Panel #10 contained a temperature sensor, and both carbon composite and adhesive-bonded samples. A figure shows the location of panels #2 and #10 on the FGB module aft endcone. The panels were subsequently returned to Earth by Space Shuttle Discovery on the STS-133/ULF-5 mission after 12 years of LEO exposure and opened in an argon chamber at the Institute of Nuclear Physics at Moscow State University in July 2011. Based on the results of analyzing the readings from sensors located on Komplast panels and in studying material samples from the panels, the comprehensive effect of spaceflight factors on the FGB (at the locations of Panels #2 and #10) was evaluated. Total solar exposure was determined to be 960 +/- 200 kJ/square cm or 21,000 equivalent solar hours. Because of location of these two panels and the ISS flight attitude, atomic oxygen (AO) fluence was relatively low for such a long duration exposure, approximately 1.5x10 (exp 21) O atoms/sq. cm. Most of the AO fluence occurred during the early phase of ISS assembly. Temperatures ranged from a maximum of 107 C to a minimum of -80 C. The MMOD environment was determined, as is shown. Interestingly, the distribution of craters and low-velocity impact particles observed in the 5-50 micron size range was approx 2-3 orders of magnitude higher than that predicted by the ORDEM2000 model of natural and artificial origin microparticles in the ISS orbit. Contamination observations were also made, through the evaluation of optical properties on thermal control coatings and by the examination of visible deposits located near several samples on the panels, primarily from elastomeric samples. These elastomer samples were extensively investigated for their physical property changes and sealing capability. In addition, 16 samples of adhesive-bonded joints using three types of epoxy adhesive were evaluated studying fracture toughness, failure surface, and adhesive volume properties. Both rubber and adhesive-bonded specimens also underwent additional (post-flight) ground-based exposures to simulate a total 30 years in LEO prior to their evaluation. Overall, results indicate that space environmental effects will not adversely impact the service life of the FGB through 2028. Our investigation is complete and a summary of the results obtained from this uniquely long-duration exposure experiment will be presented.

Shaevich, S. K.↗

Skin-Stiffener Debond Prediction Based on Computational Fracture Analysis

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used with limited success primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities however, requires the successful demonstration of the methodology on structural level. For this purpose a panel was selected that is reinforced with stringers. Shear loading causes the panel to buckle and the resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. For finite element analysis, the panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot and the panel in the vicinity of the embedded defect were modeled with a local 3D solid model. Across the width of the stringer foot the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with the mixed-mode failure criterion of the graphite/epoxy material. For small applied loads the failure index is well below one across the entire width. With increasing load the failure index approaches one first near the edge of the stringer foot from which delamination is expected to grow. With increasing delamination lengths the buckling pattern of the panel changes and the failure index increases which suggests that rapid delamination growth from the initial defect is to be expected.

Krueger, Ronald↗

Episodic plate tectonics on Venus

Studies of impact craters on Venus from the Magellan images have placed important constraints on surface volcanism. Some 840 impact craters have been identified with diameters ranging from 2 to 280 km. Correlations of this impact flux with craters on the Moon, Earth, and Mars indicate a mean surface age of 0.5 +/- 0.3 Ga. Another important observation is that 52 percent of the craters are slightly fractured and only 4.5 percent are embayed by lava flows. These observations led researchers to hypothesize that a pervasive resurfacing event occurred about 500 m.y. ago and that relatively little surface volcanism has occurred since. Other researchers have pointed out that a global resurfacing event that ceased about 500 MYBP is consistent with the results given by a recent study. These authors carried out a series of numerical calculations of mantle convection in Venus yielding thermal evolution results. Their model considered crustal recycling and gave rapid planetary cooling. They, in fact, suggested that prior to 500 MYBP plate tectonics was active in Venus and since 500 MYBP the lithosphere has stabilized and only hot-spot volcanism has reached the surface. We propose an alternative hypothesis for the inferred cessation of surface volcanism on Venus. We hypothesize that plate tectonics on Venus is episodic. Periods of rapid plate tectonics result in high rates of subduction that cool the interior resulting in more sluggish mantle convection.

Turcotte, Donald↗

Design, analysis, and fabrication of a pressure box test fixture for tension damage tolerance testing of curved fuselage panels

A pressure box test fixture was designed and fabricated to evaluate the effects of internal pressure, biaxial tension loads, curvature, and damage on the fracture response of composite fuselage structure. Previous work in composite fuselage tension damage tolerance, performed during NASA contract NAS1-17740, evaluated the above effects on unstiffened panels only. This work extends the tension damage tolerance testing to curved stiffened fuselage crown structure that contains longitudinal stringers and circumferential frame elements. The pressure box fixture was designed to apply internal pressure up to 20 psi, and axial tension loads up to 5000 lb/in, either separately or simultaneously. A NASTRAN finite element model of the pressure box fixture and composite stiffened panel was used to help design the test fixture, and was compared to a finite element model of a full composite stiffened fuselage shell. This was done to ensure that the test panel was loaded in a similar way to a panel in the full fuselage shell, and that the fixture and its attachment plates did not adversely affect the panel.

Smith, P. J.↗

Isotopic Constraints on the Genesis of Carbonates in Martian Meteorite ALH 84001

Oxygen isotopic analyses in approximately 20 micrometer spots in a chemically diverse suite of carbonates from ALH 84001 show highly variable delta(exp 18)O values from +5.4 to +25.3%. The isotopic data are correlated with the major element composition of the carbonate. The earliest forming (Ca-rich) carbonates have the lowest delta(exp 18)O values and the late-forming Mg-rich carbonates have the highest delta(exp 18)O values. Two models that can explain the isotopic variation were investigated. The carbonates could have formed in a water-rich environment at relatively low, but highly variable temperatures. In this open-system case the lower limit to the temperature variation is approximately 125 C, with fluctuations of over 250 C possible within the constraints of the model, depending on fluid composition. Alternatively the data can be explained by a closed-system model in which carbonates precipitated from a limited amount of a CO2-rich fluid. This scenario can reproduce the range of isotopic values observed, even at relatively high temperatures (greater than 500 C). Thus, the oxygen isotopic compositions do not provide unequivocal evidence for formation of the carbonates at low temperature. Neither of these scenarios is consistent with a biological origin of the carbonates and their associated features. Olivine from ALH 84001 occurs as clusters within orthopyroxene adjacent to fractures containing disrupted carbonate globules and feldspathic shock glass. The inclusions are irregular in shape and range in size from approximately 40 micrometers to submicrometer. The olivine exhibits a limited range of chemical composition from Fo(sub 65) to Fo(sub 66). We measured delta(exp 18)O values of the olivine to be +5.1 +/- 1.4%, indistinguishable within uncertainty from the host orthopyroxene. The data suggest that the olivine formed at high temperature (greater than 800 C), and is probably unrelated to carbonate formation. Instead the olivine probably formed by metamorphic reactions involving either dehydration of hydrous silicates or reduction of opx + spinel. If the reactions took place after carbonate formation, they could have caused devolatilization of siderite, producing the magnetite in the ALH carbonate globules. This scenario is also inconsistent with a biological origin for the features in ALH 84001.

Leshin, Laurie A.↗

Three-dimensional finite-element analysis of chevron-notched fracture specimens

Stress-intensity factors and load-line displacements were calculated for chevron-notched bar and rod fracture specimens using a three-dimensional finite-element analysis. Both specimens were subjected to simulated wedge loading (either uniform applied displacement or uniform applied load). The chevron-notch sides and crack front were assumed to be straight. Crack-length-to-specimen width ratios (a/w) ranged from 0.4 to 0.7. The width-to-thickness ratio (w/B) was 1.45 or 2. The bar specimens had a height-to-width ratio of 0.435 or 0.5. Finite-element models were composed of singularity elements around the crack front and 8-noded isoparametric elements elsewhere. The models had about 11,000 degrees of freedom. Stress-intensity factors were calculated by using a nodal-force method for distribution along the crack front and by using a compliance method for average values. The stress intensity factors and load-line displacements are presented and compared with experimental solutions from the literature. The stress intensity factors and load-line displacements were about 2.5 and 5 percent lower than the reported experimental values, respectively.

Raju, I. S.↗

GENOA-PFA: Progressive Fracture in Composites Simulated Computationally

GENOA-PFA is a commercial version of the Composite Durability Structural Analysis (CODSTRAN) computer program that simulates the progression of damage ultimately leading to fracture in polymer-matrix-composite (PMC) material structures under various loading and environmental conditions. GENOA-PFA offers several capabilities not available in other programs developed for this purpose, making it preferable for use in analyzing the durability and damage tolerance of complex PMC structures in which the fiber reinforcements occur in two- and three-dimensional weaves and braids. GENOA-PFA implements a progressive-fracture methodology based on the idea that a structure fails when flaws that may initially be small (even microscopic) grow and/or coalesce to a critical dimension where the structure no longer has an adequate safety margin to avoid catastrophic global fracture. Damage is considered to progress through five stages: (1) initiation, (2) growth, (3) accumulation (coalescence of propagating flaws), (4) stable propagation (up to the critical dimension), and (5) unstable or very rapid propagation (beyond the critical dimension) to catastrophic failure. The computational simulation of progressive failure involves formal procedures for identifying the five different stages of damage and for relating the amount of damage at each stage to the overall behavior of the deteriorating structure. In GENOA-PFA, mathematical modeling of the composite physical behavior involves an integration of simulations at multiple, hierarchical scales ranging from the macroscopic (lamina, laminate, and structure) to the microscopic (fiber, matrix, and fiber/matrix interface), as shown in the figure. The code includes algorithms to simulate the progression of damage from various source defects, including (1) through-the-thickness cracks and (2) voids with edge, pocket, internal, or mixed-mode delaminations.

Murthy, Pappu L. N.↗