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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Thermomechanical evolution of oceanic fracture zones

To explain the topographic and gravitational signals of fracture zones (FZs), a realistic model is developed and compared with Seasat altimeter data across three major FZs. The model is based on the thermal, subsidence, and mechanical models of the oceanic lithosphere. The model is used to calculate the FZ topography and the evolution of the geoid signature across a FZ. When lithospheric flexure is included, the model explains the persistence of sharp FZ geoid steps on old seafloor as well as the asymmetry of the geoid steps across major FZs. The model does not explain the disappearance of the thermal part of the FZ geoid signature. The model profiles are used to access the various techniques for estimating the overall geoid offsets across FZs. It is shown that steps may be underestimated by as much as 50 percent. These predictions, if correct, will force reinterpretation of published geoid height versus age data derived from altimeter profiles across FZs.

Sandwell, D. T.↗
NASA NTRS1985

Assessment of global and regional tectonic models for faulting in the ancient terrains of Mars

The structural mapping of tectonic landforms on Mars has provided a great deal of information regarding the planet's tectonic history. The present paper has the objective to test models of regional and global fracturing by evaluating the spatial distribution and orientation of the ancient heavily cratered terrain (HCT) structural landforms. Attention is given to background information concerning earlier studies, landform selection and analysis, rose diagrams, tests for Tharsis generation, and a calculation of thermal stress in the early Martian lithosphere. The contribution of classical Tharsis-centered tectonism to the production of HCT structural landforms was found to be very small. Few structural landforms investigated are either perpendicular to calculated principal stress trajectories or radial to Tharsis.

Schultz, R. A.↗
NASA NTRS1985

Delamination and debonding of materials

The general topics consist of stress analysis, mechanical behavior, and fractography/NDI of composite laminates. Papers are presented on a dynamic hybrid finite-element analysis for interfacial cracks in composites, energy release rate during delamination crack growth in composite laminates, matrix deformation and fracture in graphite-reinforced epoxies, and the role of delamination and damage development on the strength of thick notched laminates. In addition, consideration is given to a new ply model for interlaminar stress analysis, a fracture mechanics approach for designing adhesively bonded joints, the analysis of local delaminations and their influence on composite laminate behavior, and moisture and temperature effects on the mixed-mode delamination fracture of unidirectional graphite/epoxy.

Johnson, W. S.↗
NASA NTRS1988

Application of fracture mechanics to predict the strength of filament-wound composite specimens with surface damage

The strength of a thick surface notched filament wound graphite/epoxy tensile specimen was predicted using fracture mechanics analysis. The prediction model depended upon the depth of the surface notch and upon the type of failure. Specimens failed in two stages: first a sublaminate with thickness equal to the depth of the surface cut, and second, for increasing load, a sublaminate with thickness equal to the initial thickness minus the cut depth. Analytical predictions and experimental data were in good agreement.

Morris, D. H.↗
NASA NTRS1990

Corrosion fatigue crack propagation in metals

This review assesses fracture mechanics data and mechanistic models for corrosion fatigue crack propagation in structural alloys exposed to ambient temperature gases and electrolytes. Extensive stress intensity-crack growth rate data exist for ferrous, aluminum and nickel based alloys in a variety of environments. Interactive variables (viz., stress intensity range, mean stress, alloy composition and microstructure, loading frequency, temperature, gas pressure and electrode potential) strongly affect crack growth kinetics and complicate fatigue control. Mechanistic models to predict crack growth rates were formulated by coupling crack tip mechanics with occluded crack chemistry, and from both the hydrogen embrittlement and anodic dissolution/film rupture perspectives. Research is required to better define: (1) environmental effects near threshold and on crack closure; (2) damage tolerant life prediction codes and the validity of similitude; (3) the behavior of microcrack; (4) probes and improved models of crack tip damage; and (5) the cracking performance of advanced alloys and composites.

Gangloff, Richard P.↗
NASA NTRS1997

Prediction and Verification of Ductile Crack Growth from Simulated Defects in Strength Overmatched Butt Welds

Defects that develop in welds during the fabrication process are frequently manifested as embedded flaws from lack of fusion or lack of penetration. Fracture analyses of welded structures must be able to assess the effect of such defects on the structural integrity of weldments; however, the transferability of R-curves measured in laboratory specimens to defective structural welds has not been fully examined. In the current study, the fracture behavior of an overmatched butt weld containing a simulated buried, lack-of-penetration defect is studied. A specimen designed to simulate pressure vessel butt welds is considered; namely, a center crack panel specimen, of 1.25 inch by 1.25 inch cross section, loaded in tension. The stress-relieved double-V weld has a yield strength 50% higher than that of the plate material, and displays upper shelf fracture behavior at room temperature. Specimens are precracked, loaded monotonically while load-CMOD measurements are made, then stopped and heat tinted to mark the extent of ductile crack growth. These measurements are compared to predictions made using finite element analysis of the specimens using the fracture mechanics code Warp3D, which models void growth using the Gurson-Tvergaard dilitant plasticity formulation within fixed sized computational cells ahead of the crack front. Calibrating data for the finite element analyses, namely cell size and initial material porosities are obtained by matching computational predictions to experimental results from tests of welded compact tension specimens. The R-curves measured in compact tension specimens are compared to those obtained from multi-specimen weld tests, and conclusions as to the transferability of R-curves is discussed.

Nishioka, Owen S.↗
NASA NTRS2000

Effects of Interface Coating and Nitride Enhancing Additive on Properties of Hi-Nicalon SiC Fiber Reinforced Reaction-Bonded Silicon Nitride Composites

Strong and tough Hi-Nicalon SiC fiber reinforced reaction-bonded silicon nitride matrix composites (SiC/ RBSN) have been fabricated by the fiber lay-up approach. Commercially available uncoated and PBN, PBN/Si-rich PBN, and BN/SiC coated SiC Hi-Nicalon fiber tows were used as reinforcement. The composites contained approximately 24 vol % of aligned 14 micron diameter SiC fibers in a porous RBSN matrix. Both one- and two-dimensional composites were characterized. The effects of interface coating composition, and the nitridation enhancing additive, NiO, on the room temperature physical, tensile, and interfacial shear strength properties of SiC/RBSN matrix composites were evaluated. Results indicate that for all three coated fibers, the thickness of the coatings decreased from the outer periphery to the interior of the tows, and that from 10 to 30 percent of the fibers were not covered with the interface coating. In the uncoated regions, chemical reaction between the NiO additive and the SiC fiber occurs causing degradation of tensile properties of the composites. Among the three interface coating combinations investigated, the BN/SiC coated Hi-Nicalon SiC fiber reinforced RBSN matrix composite showed the least amount of uncoated regions and reasonably uniform interface coating thickness. The matrix cracking stress in SiC/RBSN composites was predicted using a fracture mechanics based crack bridging model.

Bhatt, Ramakrishana T.↗
NASA NTRS2005

Toughened Matrix SiC Fiber Reinforced Composites

First matrix cracking stress is a critical parameter for application of Sic fiber reinforced composites in highly stressed, environmentally demanding applications such as turbine blades. High matrix fracture toughness is a key property that contributes to high composite fracture stress. Silicon nitride offers reduced matrix elastic modulus, lower coefficient of thermal expansion, and potentially high fracture toughness compared to Sic matrices. All of these factors can be used to advantage to increase matrix fracture stress. As a first model system we are pursuing toughened silicon nitride matrix composites reinforced with SCS-9 fibers. Fabrication is by tape casting the matrix plies and tape lay-up with fiber plies followed by hot pressing at 1800 C. Progress toward this end will be reported.

Levine, Stanley R.↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A nonlinear elasto-plastic bond model for the discrete element modeling of woody biomass particles

This work presents an experiment-informed, semi-empirical, elasto-plastic bond model for the discrete element modeling of woody biomass particles. The model renders nonlinear plastic deformation of materials when subjected to compression/tension, bending, and twisting, essential for accurately simulating the behavior of biomass in comminution. The model is implemented in an open-source DEM package LIGGGHTS and is assessed in a number of verification tests. The model is applied to simulate the fracturing test of notched loblolly pine blocks. Reduced bond strength is prescribed to the weak planes along the growth rings. It is found that the square-root formulation for describing the elasto-plastic normal bond contact delivers the best agreement with the experimental data over its linear and quadratic counterparts. Furthermore, the model predicts a strong anisotropic behavior of wood blocks with regard to the grain orientations, indicating that the cutting angle should be an important attribute for the optimization of biomass comminution.

42 ENGINEERING↗
NASA NTRS1983

Microfissuring of Inconel 718

A tentative mathematical computer model of the microfissuring process during electron beam welding of Inconel 718 has been constructed. Predictions of the model are compatible with microfissuring tests on eight 0.25-in. thick test plates. The model takes into account weld power and speed, weld loss (efficiency), parameters and material characteristics. Besides the usual material characteristics (thermal and strength properties), a temperature and grain size dependent critical fracture strain is required by the model. The model is based upon fundamental physical theory (i.e., it is not a mere data interpolation system), and can be extended to other metals by suitable parameter changes.

Nunes, A. C., Jr.↗
DOE OSTI2020

Relationship between ductility and the porosity of additively manufactured AlSi10Mg

Additive manufacturing via selective laser melting can result in variable levels of internal porosity both between build plates and within components from the same build. In this work, sample porosity levels were compared to tensile properties for 176 samples spanning eight different build plates. Sample porosity was measured both by Archimedes density, which provided an estimation of overall porosity, and by observation of voids in the fracture surface, which provided an estimation of the porosity at the failure plane. The porosity observed at the fracture surface consistently demonstrated higher porosity than that suggested by Archimedes density. The porosity values obtained from both methods were compared against the mechanical results. Sample porosity appears to have some correlation to the ultimate tensile strength, yield strength, and modulus, but the strongest relationship is observed between porosity and ductility. Three different models were used to relate the fracture surface porosity to the ductility. The first method was a simple linear regression analysis, while the other two models have been used to relate porosity to ductility in cast alloys. It is shown that all three models fit the data well over the observed porosity ranges, suggesting that the models taken from casting theory can extend to additively manufactured metals. Finally, it is proposed that the non-destructive Archimedes method could be used to estimate an approximate sample ductility through the use of correlations realized here. Such a relationship could prove useful for design and for a deeper understanding of the impact of pores on tensile behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2019

EGS Collab Testbed 1 Earth Model Input Files

The EGS Collab is conducting experiments in hydraulic fracturing at a depth of 1.5 km in the Sanford Underground Research Facility (SURF) on the 4850 Level. A total of eight ~60m-long subhorizontal boreholes were drilled at that depth on the western rib of the West Access Drift. Six of these holes are used for geophysical monitoring, one is used for hydraulic fracturing and the remaining hole was designed as a production borehole. In addition to these eight boreholes, 4 5-m Jack leg boreholes were drilled for housing geophones. This submission package includes various data type that were assembled to create Earth Models of the testbed. Note: The coordinate system used is local Homestake Mine Coordinate (HMC) system from an old gold mine that was in operation for over 100 years.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Modeling heat transport processes in enhanced geothermal systems: Validation study from EGS Collab Experiment 1

Heat recovery from enhanced geothermal systems (EGS) is a complex process involving heat transport in both fracture networks and rock formations. A comprehensive understanding of and the ability to model the underlying heat transport mechanisms is important for the success of EGS but remains challenging in practice due to the generally insufficient characterization of EGS reservoirs. In the present study, we analyze an extensively monitored intermediate-scale EGS field experiment performed in a well-characterized testbed. The high-resolution, high-quality measurements from the field experiment enable the development of a high-fidelity model incorporating a well-constrained fracture network. Based on the field experiment, we investigate the complex heat transport processes in an EGS-relevant environment and validate the capability of a numerical approach in simulating these inherently coupled heat transport processes. A series of numerical simulations were performed to study the effects of different heat transport mechanisms, including thermal convection with fracture flow, thermal conduction in rock formations, and the Joule-Thomson effect. The agreement of thermal responses between field measurements and simulation results indicates that our numerical approach can appropriately model the heat transport processes pertaining to heat recovery from EGS reservoirs.

Wu, Hui↗
DOE OSTI2023

Discrete element model for powder grain interactions under high compressive stress

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. Finally, as validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Crack Development in Cross-Ply Laminates Under Uniaxial Tension

This study addresses matrix-dominated failures in carbon fiber/polymer matrix composite laminates in a cross-ply lay-up. The events of interest are interlaminar fracture in the form of transverse cracks in the 90' plies and longitudinal splitting in the 0 deg plies and interlaminar fracture in the form of 0/90 delamination. These events were observed using various nondestructive evaluation (NDE) techniques during static tensile tests. Acoustic emission (AE) x radiography, and edge view microscopy were the principal ones utilized in a real-time environment. A comparison of the NDE results with an analytical model based on the classical linear fracture mechanics concept of strain energy release rate as a criterion for crack growth was performed. The virtual crack closure theory was incorporated with a finite element model to generate strain energy release rate curves for the analytical case. Celion carbon fiber/polyimide matrix (G30-500/PMR-15) was the material tested with cross-ply lay-ups of (0(2)/90(6))s and (0(4)/90(4))s. The test specimens contained thermally induced cracks caused by the high-temperature processing. The analytical model was updated to compensate for the initial damage and to study further accumulation by taking into account the crack interactions. By correlating the experimental and analytical data, the critical energy release rates were found for the observable events of interest.

Gyekenyesi, Andrew L.↗
DOE OSTI2021

Simulation of hardened cement degradation and estimation of uncertainty in predicted failure times with peridynamics

Modeling the degradation of cement-based infrastructure due to aqueous environmental conditions continues to be a challenge. In order to develop a capability to predict concrete infrastructure failure due to chemical degradation, here we created a chemomechanical model of the effects of long-term water exposure on cement paste. The model couples the mechanical static equilibrium balance with reactive–diffusive transport and incorporates fracture and failure via peridynamics (a meshless simulation method). The model includes fundamental aspects of degradation of ordinary Portland cement (OPC) paste, including the observed softening, reduced toughness, and shrinkage of the cement paste, and increased reactivity and transport with water induced degradation. This version of the model focuses on the first stage of cement paste decalcification, the dissolution of portlandite. Given unknowns in the cement paste degradation process and the cost of uncertainty quantification (UQ), we adopt a minimally complex model in two dimensions (2D) in order to perform sensitivity analysis and UQ. We calibrate the model to existing experimental data using simulations of common tests such as flexure, compression and diffusion. Then we calculate the global sensitivity and uncertainty of predicted failure times based on variation of eleven unique and fundamental material properties. We observed particularly strong sensitivities to the diffusion coefficient, the reaction rate, and the shrinkage with degradation. Also, the predicted time of first fracture is highly correlated with the time to total failure in compression, which implies fracture can indicate impending degradation induced failure; however, the distributions of the two events overlap so the lead time may be minimal. Extension of the model to include the multiple reactions that describe complete degradation, viscous relaxation, post-peak load mechanisms, and to three dimensions to explore the interactions of complex fracture patterns evoked by more realistic geometry is straightforward and ongoing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Residual Stresses Affect Cell Fragment Movement

Predictive modeling tools such as the finite element method can be of tremendous help in assessing the reliability and long-term performance of photovoltaic modules. In order to obtain accurate results, the proper modeling of materials and manufacturing processes are of utmost importance. Module fabrication introduces thermo-mechanical stresses inside the module laminate, which need to be accounted for as residual stresses in finite element simulations. We found that cell fragment movement and crack opening displacements of fractured silicon cells within modules are affected by those residual stresses. Cell cracking remains a challenging topic in assessing the reliability and durability of damaged modules. Hence, accurately quantifying the separation and movement between cell fragments creates the foundation for establishing reliable lifetime and performance assessments of fractured silicon modules. Here, we present a modeling approach that uses upper and lower bounds to accurately account for the residual stresses introduced by the module lamination process. We designed a four-point flexure coupon test of a laminated, fractured silicon strip to validate our numerical results and found good agreement between our modeling methodology and the experimental data. Finally, we discuss the implications of the residual stresses on the normal crack opening and metallization wear-out of fractured silicon cells.

14 SOLAR ENERGY↗