Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum damage”

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 127 records · Page 7

Mesh objective stochastic simulations of quasibrittle fracture

Continuum finite element (FE) modeling of damage and failure of quasibrittle structures suffers from the spurious mesh sensitivity due to strain localization. Here this issue has been addressed for deterministic analysis through the development of localization limiters. Here this study proposes a mechanism-based model to mitigate the mesh sensitivity in stochastic FE simulations of quasibrittle fracture. The interest is placed on the analysis of large-size structures, where the mesh size is conveniently chosen to be larger than the width of the fracture process zone as well as the correlation length of the random fields of constitutive properties. The present model is formulated within the framework of continuum damage mechanics. Two localization parameters are introduced to describe the evolution of the damage pattern of each finite element. These parameters are used to guide the energy regularization of the constitutive law, as well as to formulate the mesh-dependent probability distributions of constitutive properties. Depending on the prevailing damage pattern, different energy regularization schemes and mesh dependence of the probability distribution functions are used in the constitutive law. The model is applied to simulate the stochastic failure behavior of quasibrittle structures of different geometries featuring different failure processes including damage initiation, localization, and propagation. It is shown that using fixed probability distribution functions of constitutive properties could lead to strong mesh dependence of the prediction of the mean and variance of the structural load capacity. The probability distribution functions of constitutive properties must be linked to the damage pattern, which may evolve during the failure process. Such a mechanism-based modeling of the probability distributions of constitutive properties is essential for mitigating the spurious mesh sensitivity in stochastic FE analysis of quasibrittle fracture.

42 ENGINEERING↗

A study of material damping in large space structures

A constitutive model was developed for predicting damping as a function of damage in continuous fiber reinforced laminated composites. The damage model is a continuum formulation, and uses internal state variables to quantify damage and its subsequent effect on material response. The model is sensitive to the stacking sequence of the laminate. Given appropriate baseline data from unidirectional material, and damping as a function of damage in one crossply laminate, damage can be predicted as a function of damage in other crossply laminates. Agreement between theory and experiment was quite good. A micromechanics model was also developed for examining the influence of damage on damping. This model explicitly includes crack surfaces. The model provides reasonable predictions of bending stiffness as a function of damage. Damping predictions are not in agreement with the experiment. This is thought to be a result of dissipation mechanisms such as friction, which are not presently included in the analysis.

Highsmith, A. L.↗

Creep rupture analysis of a beam resting on high temperature foundation

A simplified uniaxial strain controlled creep damage law is deduced with the use of experimental observation from a more complex strain dependent law. This creep damage law correlates the creep damage, which is interpreted as the density variation in the material, directly with the accumulated creep strain. Based on the deduced uniaxial strain controlled creep damage law, a continuum mechanical creep rupture analysis is carried out for a beam resting on a high temperature elastic (Winkler) foundation. The analysis includes the determination of the nondimensional time for initial rupture, the propagation of the rupture front with the associated thinning of the beam, and the influence of creep damage on the deflection of the beam. Creep damage starts accumulating in the beam as soon as the load is applied, and a creep rupture front develops at and propagates from the point at which the creep damage first reaches its critical value. By introducing a series of fundamental assumptions within the framework of technical Euler-Bernoulli type beam theory, a governing set of integro-differential equations is derived in terms of the nondimensional bending moment and the deflection. These governing equations are subjected to a set of interface conditions at the propagating rupture front. A numerical technique is developed to solve the governing equations together with the interface equations, and the computed results are presented and discussed in detail.

Gu, Randy J.↗

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Unified cohesive zone model (UCZM) for fracturing and fragmenting solids

Here, a Unified Cohesive Zone Model (UCZM), which inherits most of the advantages while overcoming the shortcomings of existing Cohesive Zone Models (CZMs), is proposed. Similar to the traditional extrinsic CZM approach, UCZM dynamically inserts the cohesive elements into the system based on local material states (e.g., stress, strain). However, the transition from continua to discontinua is smoothly achieved, thereby eliminating the “time-discontinuous” issue seen in the extrinsic CZM. Moreover, within the novel UCZM framework, the point of transition from continua to discontinua is controllable through the introduction of crack initialization criteria. As a result, the UCZM allows any material models (e.g., elastic, plastic, damage models) for continuum solids and for discrete fracture behavior to work together. In essence, both an enhanced extrinsic cohesive zone model and an intrinsic cohesive zone model can be represented by the proposed unified model. The proposed UCZM has been verified through different numerical examples. The work demonstrates that the UCZM is a highly effective approach for modeling fracture and fragmentation processes in solids.

42 ENGINEERING↗

Coupled diffusion-deformation-damage model for polymers used in hydrogen infrastructure

The soft materials used in the infrastructure of hydrogen storage and distribution systems are vulnerable because exposure to high-pressure hydrogen can lead to mechanical damage and property degradation. Some materials, such as polymers, can undergo mechanical failure due to a phenomenon wherein hydrogen gas diffuses through the polymer chains and occupies preexisting cavities or voids inside the polymer material. If the hydrogen gas pressure is reduced more quickly than hydrogen can diffuse back out of the polymer which is necessary for some applications, the trapped hydrogen gas instead escapes by rupturing the material, causing surface blistering or permanent damage. In this study, a continuum mechanics-based fully coupled diffusion-deformation model with damage is developed to predict the stress distribution and damage propagation while the polymer undergoes rapid decompression failure. The hyperelastic material model, along with the maximum principal strain failure theory, was chosen for this study as it represents the nonlinear material response with brittle failure observed in uniaxial tensile tests perfectly. EPDM polymer was chosen for this study because of its commercial availability and common use in hydrogen storage and distribution system. It has superior mechanical properties, high and low temperature resistance, and certain compounds work well in hydrogen gas. This work is useful for design engineers to alter the parameters while manufacturing polymer composites to increase their performance in a high-pressure hydrogen environment.

Kulkarni, Shank S.↗

A cumulative damage model for continuous fiber composite laminates with matrix cracking and interply delaminations

The present cumulative damage model for the prediction of stiffness loss in graphite/epoxy laminates applies a thermomechanical constitutive theory for elastic composites with distributed damage. The model proceeds from a continuum mechanics and thermodynamics approach in which the distributed damage is characterized by a set of second-order tensor-valued internal state variables. A set of damage-dependent laminated plate equations is obtained; this is developed by modifying classical Kirchhoff plate theory.

Allen, David H.↗

Fracture Dynamics in Silicon Anode Solid-State Batteries

Solid-state batteries (SSBs) with silicon anodes could enable improved safety and energy density compared to lithium-ion batteries. However, degradation arising from the massive volumetric changes of silicon anodes during cycling is not well understood in solid-state systems. Here, we use operando X-ray computed microtomography to reveal micro- to macro-scale chemo-mechanical degradation processes of silicon anodes in SSBs. Mud-type channel cracks driven by biaxial tensile stress form across the electrode during delithiation. We also find detrimental cracks at the silicon/ solid electrolyte interface that form due to local reaction competition between neighboring domains of different sizes. Continuum phase-field damage modeling quantifies stress-driven channel cracking and shows that the lithiated silicon stress state is critical for determining the extent of interfacial fracture. This work reveals the mechanisms that govern SSBs compared to conventional lithium-ion batteries and provides guidelines for engineering chemomechanically resilient electrodes for high-energy batteries.

25 ENERGY STORAGE↗

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗

Sound speed measurements in shock-compressed cemented Tungsten carbides at pressures up to 100 GPa

To gain insights into thermodynamic states attained during shock compression of cemented tungsten carbide with 3.7 wt.% cobalt binder, we present results of longitudinal sound (release wave) speed measurements and their analysis at peak stresses up to 100 GPa (volumetric compression ratio ~ 15%). The sound speeds are determined using front-surface impact and release-wave overtake plate impact experimental configurations using laser interferometry. The measured sound speed data along with estimates for bulk sound speeds obtained using the fourth-order Birch-Murnaghan EoS and thermodynamics are used to determine the longitudinal moduli and shear moduli of shocked tungsten carbide at the various peak compression states attained in the experiments. Here, the longitudinal sound speeds were found to increase linearly with volume compression ratio from 6.97 ± 0.010 km/s at ambient conditions to 8.26 ± 0.156 km/s at a volume compression ratio of ~ 15%. The corresponding longitudinal elastic moduli also increase nearly linearly with the volume compression ratio but remain consistently lower than their theoretical predictions based on continuum models with no damage. Also, the sensitivity of shear moduli to pressure, as predicted by the Steinberg-Guinan model, is reduced substantially and the shear moduli of cemented WC with 3.7 wt.% Co remains nearly constant at ~ 310 GPa at the various peak compression stress states investigated in the present study.

36 MATERIALS SCIENCE↗

Physical, resource supply, and biological controls on nutrient processing along the river continuum

Nutrient impairment has led to damages to US surface and groundwater systems in excess of 100 billion dollars per year. Therefore, there is a strong need to develop methods to predict the transport, uptake, and export of nutrients along fluvial networks. We present results that are based on a data-driven mechanistic understanding of three factors that largely control nutrient uptake and export: 1) interactions between transport-related processes (mass transfer to metabolically active zones), 2) resource supply dynamics (nutrient concentration, stoichiometric constraints, etc.), and 3) biological controls (microbial community structure and function). Our results were generated from column experiments conducted along the Jemez River-Rio Grande continuum, which spans four orders of magnitude in mean annual discharge, more than 2000 m in altitude, and more than 500 km of stream longitude. Two resource supply injections were performed on each of the columns, i.e., a nitrate only addition, followed by a stoichiometrically ‘balanced’ 106Carbon:16Nitrogen:1Phosphorus addition. We quantified NO3-N uptake kinetics while constraining three variables: stream order, sediment type and type of injection (N vs stoichiometrically ‘balanced’ C:N:P). Following the laboratory nutrient uptake experiments, the columns were destructively sampled and the contents were homogenized to collect subsamples for DNA sequencing. Amplicon analysis was carried out as described by the Earth Microbiome Protocol for 16s and ITS sequencing.

54 ENVIRONMENTAL SCIENCES↗

A Continuum Dislocation Dynamics Crystal Plasticity Approach to Irradiated Body-Centered Cubic α-Iron

Radiation-induced embrittlement of reactor pressure vessel (RPV) steels can potentially limit the operating life of nuclear power plants. Over extended exposure to radiation doses, these body-centered cubic (BCC) irons demonstrate irradiation damage. Here, we present a continuum dislocation density (CDD) crystal plasticity model to capture the interaction among dislocations and self-interstitial atom (SIA) loops in α-iron. We demonstrate the importance of modeling cross slip using a combined stochastic Monte Carlo approach and the role of slip system strength anisotropy in capturing stochastic cross slip interactions. Through these captured interactions, the CDD crystal plasticity model can capture both the stress response and the physical evolution of dislocations on different slip system planes. Single-crystal verification experiments are used to calibrate the CDD crystal plasticity model, and a set of simplified polycrystalline simulations demonstrates the model’s ability to capture the stress response from tensile experiments on α-iron.

36 MATERIALS SCIENCE↗

Capturing the fractocohesive length scale in elastomers through a statistical mechanics-based gradient enhanced damage model

Here, this study aims to examine modeling flaw sensitivity in elastomers. The direct incorporation of polymer chain statistical mechanics considerations into a continuum stretch-based gradient-enhanced damage formulation, in turn, allows a representation of diffuse chain damage and fracture events that align with known micromechanical mechanisms. Through a series of numerical experiments, we simulate crack propagation and extract the fracture energy as an output of the model, while keeping track of the micromechanical signatures of diffuse chain damage that accommodate fracture propagation and eventually influence flaw sensitivity. Finally, by combining the fracture toughness and the work to rupture, we identify a fractocohesive length of the material, corresponding to the full width of the damage process zone. As the damage-to-fracture cascade in the proposed GED model is influenced by the introduction of a length scale associated with network imperfection and long-range load transfer, the emerging relationship of the two length scales is discussed, providing a potential link between microscopic damage mechanisms and the observed macroscopic fracture response.

36 MATERIALS SCIENCE↗

Effect of surface defects on damage and failure of dynamically driven cylindrical shells

Here, in this work, expanding cylinder experiments are performed using a gas gun on the aluminum 6061 alloy to understand the dynamic deformation and failure behavior of this material in the presence of surface defects in the form of scratches. The experiments are then modeled with the modified Tepla (for tensile plasticity) model, which is equipped with the damage yield function of Gurson [“Continuum theory of ductile rupture by void nucleation and growth: Part I—Yield criteria and flow rules for porous ductile media,” (1977)] and Tvergaard and Needleman [Acta Metall. 32(1), 157–169 (1984)], over-stress viscosity [F. Addessio and J. Johnson, J. Appl. Phys. 74(3), 1640–1648 (1993)], and the underlying viscoplastic model of Preston et al. [J. Appl. Phys. 93(1), 211–220 (2003)] (Preston–Tonks–Wallace model) for the dense material. The Tepla model parameters were calibrated using the Bayesian approach against experimentally measured free-surface velocity data obtained from plate impact experiments on the aluminum 6061 alloy. The evolution of radial expansion and shapes of the deformed cylinder from the simulations on a pristine cylinder wall and a cylinder wall with a longitudinal surface defect shows excellent agreement with the corresponding experimental results. Finally, the role of surface defect size and cylinder wall geometry on the dynamic strength of the cylinder wall was also investigated through a series of simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mechanics of composite materials: Recent advances; Proceedings of the Symposium, Virginia Polytechnic Institute and State University, Blacksburg, VA, August 16-19, 1982

The present conference on the mechanics of composites discusses microstructure's influence on particulate and short fiber composites' thermoelastic and transport properties, the elastoplastic deformation of composites, constitutive equations for viscoplastic composites, the plasticity and fatigue of metal matrix composites, laminate damping mechanisms, the micromechanical modeling of Kevlar/epoxy composites' time-dependent failure, the variational characterization of waves in composites, and computational methods for eigenvalue problems in composite design. Also discussed are the elastic response of laminates, elastic coupling nonlinear effects in unsymmetrical laminates, elasticity solutions for laminate problems having stress singularities, the mechanics of bimodular composite structures, the optimization of laminated plates and shells, NDE for laminates, the role of matrix cracking in the continuum constitutive behavior of a damaged composite ply, and the energy release rates of various microcracks in short fiber composites.

Hashin, Z.↗

Crack layer theory

A damage parameter is introduced in addition to conventional parameters of continuum mechanics and consider a crack surrounded by an array of microdefects within the continuum mechanics framework. A system consisting of the main crack and surrounding damage is called crack layer (CL). Crack layer propagation is an irreversible process. The general framework of the thermodynamics of irreversible processes are employed to identify the driving forces (causes) and to derive the constitutive equation of CL propagation, that is, the relationship between the rates of the crack growth and damage dissemination from one side and the conjugated thermodynamic forces from another. The proposed law of CL propagation is in good agreement with the experimental data on fatigue CL propagation in various materials. The theory also elaborates material toughness characterization.

Chudnovsky, A.↗

Crack layer theory

A damage parameter is introduced in addition to conventional parameters of continuum mechanics and consider a crack surrounded by an array of microdefects within the continuum mechanics framework. A system consisting of the main crack and surrounding damage is called crack layer (CL). Crack layer propagation is an irreversible process. The general framework of the thermodynamics of irreversible processes are employed to identify the driving forces (causes) and to derive the constitutive equation of CL propagation, that is, the relationship between the rates of the crack growth and damage dissemination from one side and the conjugated thermodynamic forces from another. The proposed law of CL propagation is in good agreement with the experimental data on fatigue CL propagation in various materials. The theory also elaborates material toughness characterization.

Chudnovsky, A.↗

Inelastic Deformation of Metal Matrix Composites: Plasticity and Damage Mechanisms - Part 1

The deformation mechanisms of a Ti 15-3/SCS6 (SiC fiber) metal matrix composite (MMC) were investigated using a combination of mechanical measurements and microstructural analysis. The objectives were to evaluate the contributions of plasticity and damage to the overall inelastic response, and to confirm the mechanisms by rigorous microstructural evaluations. The results of room temperature experiments performed on 0 degree and 90 degree systems primarily are reported in this report. Results of experiments performed on other laminate systems and at high temperatures will be provided in a forthcoming report. Inelastic deformation of the 0 degree MMC (fibers parallel to load direction) was dominated by the plasticity of the matrix. In contrast, inelastic deformations of the 90 degree composite (fibers perpendicular to loading direction) occurred by both damage and plasticity. The predictions of a continuum elastic plastic model were compared with experimental data. The model was adequate for predicting the 0 degree response; however, it was inadequate for predicting the 90 degree response largely because it neglected damage. The importance of validating constitutive models using a combination of mechanical measurements and microstructural analysis is pointed out. The deformation mechanisms, and the likely sequence of events associated with the inelastic deformation of MMCs, are indicated in this paper.

Majumdar, B. S.↗