SEARCH · Search NASA

Results for “von mises”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Elevated temperature axial and torsional fatigue behavior of Haynes 188

The results are reported for high-temperature axial and torsional low-cycle fatigue experiments performed at 760 C in air on thin-walled tubular specimens of Haynes 188, a wrought cobalt-based superalloy. Data are also presented for mean coefficient of thermal expansion, elastic modulus, and shear modulus at various temperatures from room to 1000 C, and monotonic and cyclic stress-strain curves in tension and in shear at 760 C. This data set is used to evaluate several multiaxial fatigue life models (most were originally developed for room temperature multiaxial life prediction) including von Mises equivalent strain range (ASME Boiler and Pressure Code), Manson-Halford, modified multiaxiality factor (proposed in this paper), modified Smith-Watson-Topper, and Fatemi-Socie-Kurath. At von Mises equivalent strain ranges (the torsional strain range divided by the square root of 3, taking the Poisson's ratio to be 0.5), torsionally strained specimens lasted, on average, factors of 2 to 3 times longer than axially strained specimens. The modified multiaxiality factor approach shows promise as a useful method of estimating torsional fatigue life from axial fatigue data at high temperatures. Several difficulties arose with the specimen geometry and extensometry used in these experiments. Cracking at extensometer probe indentations was a problem at smaller strain ranges. Also, as the largest axial and torsional strain range fatigue tests neared completion, a small amount of specimen buckling was observed.

Bonacuse, Peter J.↗

NASA NTRS1991

In-phase and out-of-phase axial-torsional fatigue behavior of Haynes 188 at 760 C

Isothermal, in-phase and out-of-phase axial-torsional fatigue experiments have been conducted at 760 C on uniform gage section, thin-walled tubular specimens of a wrought cobalt-base superalloy, Haynes 188. Test-control and data acquisition were accomplished with a minicomputer. Fatigue lives of the in- and out-of-phase axial-torsional fatigue tests have been estimated with four different multiaxial fatigue life prediction models that were developed primarly for predicting axial-torsional fatigue lives at room temperature. The models investigated were: (1) the von Mises equivalent strain range; (2) the Modified Multiaxiality Factor Approach; (3) the Modified Smith-Watson-Topper Parameter; and (4) the critical shear plane method of Fatemi, Socie, and Kurath. In general, life predictions by the von Mises equivalent strain range model were within a factor of 2 for a majority of the tests and the predictions by the Modified Multiaxiality Factor Approach were within a factor of 2, while predictions of the Modified Smith-Watson-Topper Parameter and of the critical shear plane method of Fatemi, Socie, and Kurath were unconservative and conservative, respectively, by up to factors of 4. In some of the specimens tested under combined axial-torsional loading conditions, fatigue cracks initiated near extensometer indentations. Two design modifications have been proposed to the thin-walled tubular specimen to overcome this problem.

Kalluri, Sreeramesh↗

DOE OSTI2021

Finite Element Analysis (FEA) for Water-Foam Fracturing of Granite Rock

In addition to the foam data that were obtained from literature and that were collected from the current study, simulation data was also generated from finite element analysis (FEA) conducted in this study using COMSOL Multiphysics software. The FEA models were built to simulate the experiments conducted at Oak Ridge National Laboratory (ORNL) on cement and granite samples. In these FEA models, temperature was kept at ambient while the pressure profile resembled the loading conditions during the ORNL experiments, where pressure was either monotonically increased or applied cyclically. The cement material was used as a model material and was used to study Von Mises stress and tensile stress distribution for different bore hole length geometry using a parametric sweep with water as fracturing fluid using solid-fluid interaction module. For the granite material, FEA models were developed for stress analysis of cylindrical samples with water or foam fluids. The solid mechanics module in COMSOL was implemented to solve for Von Mises stress and tensile stress. The fluid-structure interaction module was implemented to solve for water-foam interaction on granite cylinder with addition of fluid-loading on structure, i.e., large deformation in solid mechanics with no impact on fluid deformation. Foam was considered as a pseudo single-phase compressible fluid for which material properties were calculated from water and gas (nitrogen) phases. The density of foam is calculated as a function of the densities of water and nitrogen, while viscosity is a function of temperature. Four types of FEA analyses were modelled: 1. Monotonic injection with water 2. Monotonic injection with foam 3. Cyclic injection with water 4. Cyclic injection with foam All the COMSOL files are converted to a zip file which is save in .mph.

15 GEOTHERMAL ENERGY↗

NASA NTRS1999

Hydrostatic Stress Effects in Metal Plasticity

Since the 1940s, the theory of plasticity has assumed that hydrostatic stress does not affect the yield or postyield behavior of metals. This assumption is based on the early work of Bridgman. Bridgman found that hydrostatic pressure (compressive stress) does not affect yield behavior until a substantial amount of pressure (greater than 100 ksi) is present. The objective of this study was to determine the effect of hydrostatic tension on yield behavior. Two different specimen geometries were examined: an equal-arm bend specimen and a double edge notch specimen. The presence of a notch is sufficient to develop high enough hydrostatic tensile stresses to affect yield. The von Mises yield function, which does not have a hydrostatic component, and the Drucker-Prager yield function, which includes a hydrostatic component, were used in finite element analyses of the two specimen geometries. The analyses were compared to test data from IN 100 specimens. For both geometries, the analyses using the Drucker-Prager yield function more closely simulated the test data. The von Mises yield function lead to 5-10% overprediction of the force-displacement or force-strain response of the test specimens.

Wilson, Christopher D.↗

NASA NTRS1998

An Assessment of the Axial and Radial Dilation of a DPIMS Tantalum Cartridge for Space Shuttle Flight Experiments

Ground-based heat treatment tests are planned on an argon gas-filled tantalum cartridge developed as pan of a Diffusion Processes in Molten Semiconductors (DPIMS) experiment conducted on NASA's Space Shuttle. The possibility that the cartridge may creep during testing and touch the furnace walls is of real concern in this program. The present paper discusses the results of calculations performed to evaluate this possibility. Deformation mechanism maps were constructed using literature data in order to identify the creep mechanism dominant under the appropriate stresses and temperatures corresponding to the test conditions. These results showed that power-law creep was dominant when the grain size of the material exceeded 55 gm but Coble creep was the important mechanism below this value of grain size. Finite element analysis was used to analyze the heat treatment tWs assuming a furnace run away condition (which is a worst case scenario) using the appropriate creep parameters corresponding to grain sizes of 1 and 100 gm. Calculations were also conducted to simulate the effect of an initial 3 tilt of the cartridge assembly, the maximum possible tilt angle. The von Mises stress and su-ain distributions were calculated assuming that the cartridge was fixed at one end as it was heated from ambient temperature to 1823 K in 1.42 h, maintained at 1823 K for 9.5 h and then further heated to an over temperature condition of 2028 K in 0.3 h. The inelastic axial and radial displacements of the cartridge walls were evaluated by resolving the von Mises strain along the corresponding directions. These calculations reveal that the maximum axial and radial displacements are expected to be about 2.9 and 0.25 mm, respectively, for both fine and coarse-grained materials at 2028 K. It was determined that these displacements occur during heat-up to temperature and creep of the cartridge is likely to be relatively insignificant irrespective of grain size. Furthermore, with a 3' tilt of the cartridge, the deflection is increased by only 0.39 gm which is negligible. Since the gap between the furnace heating elements and the cartridge is about 7.5 mm and less than the maximum radial dilation of 0.25 mm at 2028 K, it is concluded that the cartridge is unlikely to touch the furnace walls during the experiments.

Raj, S.V.↗

NASA NTRS2005

Glacial isostatic stress shadowing by the Antarctic ice sheet

Numerous examples of fault slip that offset late Quaternary glacial deposits and bedrock polish support the idea that the glacial loading cycle causes earthquakes in the upper crust. A semianalytical scheme is presented for quantifying glacial and postglacial lithospheric fault reactivation using contemporary rock fracture prediction methods. It extends previous studies by considering differential Mogi-von Mises stresses, in addition to those resulting from a Coulomb analysis. The approach utilizes gravitational viscoelastodynamic theory and explores the relationships between ice mass history and regional seismicity and faulting in a segment of East Antarctica containing the great Antarctic Plate (Balleny Island) earthquake of 25 March 1998 (Mw 8.1). Predictions of the failure stress fields within the seismogenic crust are generated for differing assumptions about background stress orientation, mantle viscosity, lithospheric thickness, and possible late Holocene deglaciation for the D91 Antarctic ice sheet history. Similar stress fracture fields are predicted by Mogi-von Mises and Coulomb theory, thus validating previous rebound Coulomb analysis. A thick lithosphere, of the order of 150-240 km, augments stress shadowing by a late melting (middle-late Holocene) coastal East Antarctic ice complex and could cause present-day earthquakes many hundreds of kilometers seaward of the former Last Glacial Maximum grounding line.

Coulomb stress↗

DOE OSTI2023

Impact of Plate Thickness and Joint Geometry on Residual Stresses in 347H Stainless Steel Welds

Weldments of 347H stainless steel are potentially susceptible to stress relaxation cracking at elevated service temperatures. Mitigation of stress relaxation cracking susceptibility within a multipass weld requires a good understanding of welding practices and manufacturing techniques to control high tensile residual stresses. Here in this work, the dependence of residual stress distribution in 347H stainless steel on base plate thickness, joint geometry design, and preheating condition was systematically investigated by using three-dimensional finite element models. The finite element models were validated through good agreement between neutron diffraction measurements and calculated elastic strains. The single-V-groove welds with and without a preheating step all produced similar peak von Mises residual stresses, above 450 MPa, within both the fusion zone and heat-affected zone (HAZ). In plates thicker than 0.5 in. (12.7 mm), high tensile residual stress could be observed in a relatively large area, from the middle of the plate thickness to underneath the top surface. A double-V groove shifted the high tensile stress area to the middle thickness of the weld. A single-J-groove weld was able to confine the residual stress to a very small region near the middle thickness within the fusion zone and suppressed the von Mises residual stress within the HAZ to below 400 MPa.

14 SOLAR ENERGY↗

DOE OSTI2021

Assessment of asymmetrical rheological behavior of cementitious material for 3D printing application

Among the multiplicity of construction methods based on digital concrete technologies, the most common one remains the extrusion-based process. Two opposed strategies have been developed in order to print a concrete structure using extrusion: the first one emphasizes the ease-of-pumping of material which is massively accelerated to ensure the stability of the structure, when leaving the extruder; the second one uses a firm mortar able to sustain the own weight of the structure during the printing, without requiring any accelerator. Both methods have their pros and cons and don't imply the same rheological requirements. Most of the recent literature assumes that both types of material (fluid and firm) follow the same trend of behavior and obey to the Von Mises plasticity criterion, while it has been demonstrated that firm cement-based materials can exhibit a pressure dependent behavior. Moreover, rheological study is mainly based on shear rheometry, while the extrusion-based printing process requires a description of the rheological behavior in compression to predict the global stability of the structure, and in tension to prevent cracks formation. This paper investigates the fresh behavior of cement paste and mortar under different solicitations (shear, compression, tension) for different water to cement ratios. The analysis of these results allows to describe the transition between a ductile fluid material (symmetrical in tension and compression) that obeys to a Von Mises plasticity criterion, and a brittle and firm material (asymmetrical in tension and compression) that requires a pressure dependent plasticity criterion in order to predict its strength under a given solicitation.

36 MATERIALS SCIENCE↗

NASA NTRS1992

Cyclic axial-torsional deformation behavior of a cobalt-base superalloy

Multiaxial loading, especially at elevated temperature, can cause the inelastic response of a material to differ significantly from that predicted by simple flow rules, i.e., von Mises or Tresca. To quantify some of these differences, the cyclic high-temperature, deformation behavior of a wrought cobalt-based superalloy, Haynes 188, is investigated under combined axial and torsional loads. Haynes 188 is currently used in many aerospace gas turbine and rocket engine applications, e.g., the combustor liner for the T800 turboshaft engine for the RAH-66 Comanche helicopter and the liquid oxygen posts in the main injector of the space shuttle main engine. The deformation behavior of this material is assessed through the examination of hysteresis loops generated from a biaxial fatigue test program. A high-temperature axial, torsional, and combined axial-torsional fatigue data base has been generated on Haynes 188 at 760 C. Cyclic loading tests have been conducted on uniform gauge section tubular specimens in a servohydraulic axial-torsional test rig. Test control and data acquisition were accomplished with a minicomputer. In this paper, the cyclic hardening characteristics and typical hysteresis loops in the axial stress versus axial strain, shear stress versus engineering shear strain, axial strain versus engineering shear strain, and axial stress versus shear stress spaces are presented for cyclic, in-phase and out-of-phase, axial torsional tests. For in-phase tests three different values of the proportionality constant, lambda (ratio of engineering shear strain amplitude to axial strain amplitude), are examined, viz., 0.86, 1.73, and 3.46. In the out-of-phase tests, three different values of the phase angle, phi (between the axial and engineering shear strain waveforms), are studied, viz., 30, 60, and 90 deg with lambda = 1.73. The cyclic hardening behaviors of all the tests conducted on Haynes 188 at 760 C are evaluated using the von Mises equivalent stress-strain and the maximum shear stress-maximum engineering shear strain (Tresca) curves. Comparisons are also made between the hardening behaviors of cyclic axial, torsional, and combined in-phase and out-of-phase axial-torsional fatigue tests. These comparisons are accomplished through simple Ramberg-Osgood type stress-strain functions for cyclic, axial stress-strain and shear stress-engineering shear strain curves.

Bonacuse, Peter J.↗

DOE OSTI2025

Joint Modeling of Wind Speed and Wind Direction Through a Conditional Approach

Atmospheric near surface wind speed and wind direction play an important role in many applications, ranging from air quality modeling, building design, wind turbine placement to climate change research. It is therefore crucial to accurately estimate the joint probability distribution of wind speed and direction. In this work, we develop a conditional approach to model these two variables, where the joint distribution is decomposed into the product of the marginal distribution of wind direction and the conditional distribution of wind speed given wind direction. To accommodate the circular nature of wind direction, a von Mises mixture model is used; the conditional wind speed distribution is modeled as a directional dependent Weibull distribution via a two-stage estimation procedure, consisting of a directional binned Weibull parameter estimation, followed by a harmonic regression to estimate the dependence of the Weibull parameters on wind direction. A Monte Carlo simulation study indicates that our method outperforms two other approaches in estimation efficiency: one that utilizes periodic spline quantile regression and another that generates data from the commonly used Abe-Ley distribution for cylindrical data. We illustrate our method by using the output from a regional climate model to investigate how the joint distribution of wind speed and direction may change under some future climate scenarios. Our method indicates significant changes in the variation of wind speed with respect to some directions.

17 WIND ENERGY↗

NASA NTRS1934

On the theory of laminar boundary layers involving separation

This paper presents a mathematical discussion of the laminar boundary layer, which was developed with a view of facilitating the investigation of those boundary layers in particular for which the phenomenon of separation occurs. The treatment starts with a slight modification of the form of the boundary layer equation first published by Von Mises. Two approximate solutions of this equation are found, one of which is exact at the outer edge of the boundary layer while the other is exact at the wall. The final solution is obtained by joining these two solutions at the inflection points of the velocity profiles. The final solution is given in terms of a series of universal functions for a fairly broad class of potential velocity distributions outside of the boundary layer. Detailed calculations of the boundary layer characteristics are worked out for the case in which the potential velocity is a linear function of the distance from the upstream stagnation point. Finally, the complete separation point characteristics are determined for the boundary layer associated with a potential velocity distribution made up of two linear functions of the distance from the stagnation point. It appears that extensions of the detailed calculations to more complex potential flows can be fairly easily carried out by using the explicit formulae given in the paper. (author)

Von Karman, TH↗

DOE OSTI2023

A new method for detecting abrupt changes in the dependence among multivariate hydrological series based on moving cut total correlation

Knowledge of how to define and estimate the dependence among multivariate hydrological series is essential for detecting abrupt changes in the dependence. Here, in this paper, a new method (BMCTC) is proposed to detect all possible abrupt change points in the dependence among multivariate hydrological series. The total correlation estimated by the matrix-based Renyi's alpha-order entropy functional is firstly introduced to define and measure the dependence strength among multivariate hydrological series. Then, the moving cut total correlation (MCTC) sequence is built by the moving window technique, which is used to measure changes in the dependence strength among multivariate hydrological series. Finally, the Bernaola-Galvan algorithm is used to detect all change points of the MCTC sequence. Simulations are performed to compare the effectiveness of BMCTC with Pearson correlation (BMCPC) and Spearman correlation (BMCSC), Cramer-von Mises (CvM) and copula-based likelihood-ratio (CLR). The results show that all change points are detected by BMCTC regardless of the samples size, but wrong change points or no change points are detected by other methods in most cases. BMCTC is applied to detect change points in the dependence among annual runoff, precipitation and sediment discharge series in the Xiliugou and the Kuyehe River, China. It is found that the dependence among runoff, precipitation and sediment discharge changed abruptly in 1980 and 1996 in the Kuyehe River and in 1999 in the Xiliugou River. These changes are mainly caused by human activities such as construction of water conservancy projects and coal mining.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Accurate Effective Stress Measures: Predicting Creep Life for 3D Stresses Using 2D and 1D Creep Rupture Simulations and Data

Operating structural components experience complex loading conditions resulting in 3D stress states. Current design practice estimates multiaxial creep rupture life by mapping a general state of stress to a uniaxial creep rupture correlation using effective stress measures. The data supporting the development of effective stress measures are nearly always only uniaxial and biaxial, as 3D creep rupture tests are not widely available. This limitation means current effective stress measures must extrapolate from 2D to 3D stress states, potentially introducing extrapolation error. In this work, we use a physics-based, crystal plasticity finite element model to simulate uniaxial, biaxial, and triaxial creep rupture. Here, we use the virtual dataset to assess the accuracy of current and novel effective stress measures in extrapolating from 2D to 3D stresses and also explore how the predictive accuracy of the effective stress measures might change if experimental 3D rupture data was available. We confirm these conclusions, based on simulation data, against multiaxial creep rupture experimental data for several materials, drawn from the literature. The results of the virtual experiments show that calibrating effective stress measures using triaxial test data would significantly improve accuracy and that some effective stress measures are more accurate than others, particularly for highly triaxial stress states. Results obtained using experimental data confirm the numerical findings and suggest that a unified effective stress measure should include an explicit dependence on the first stress invariant, the maximum tensile principal stress, and the von Mises stress.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Anisotropic yield models for lattice unit cell structures exploiting orthotropic symmetry

We report numerical homogenization enables efficient computational analysis and design of multiscale structures made of micro-architected materials including lattice unit cells. However, predicting yield is nontrivial because it requires accurate and efficient predictions of the maximum stress inside the homogenized unit cells. To address this challenge, we develop a macroscale anisotropic yield function for micro-architected materials. The yield function depends on the three-dimensional macroscale stress state and the parameters describing a family of micro-architectures, such as the radii of the struts in a lattice unit cell. To ensure accuracy, we determine the maximum stress using three-dimensional continuum finite-element analysis. To ensure efficiency, we construct surrogate models from the aforementioned high-fidelity results for yield prediction. To reduce simulation costs and surrogate modeling complexity, we leverage orthotropic symmetry commonly found in lattice unit cells. In this paper, we provide a thorough presentation of group representations and the systematic procedure to exploit orthotropic domain symmetry in homogenization and surrogate modeling. We illustrate the use of linear homogenization results to predict yield in specific unit cells without further simulations. We furthermore show that surrogate modeling presents a viable option for anisotropic yield prediction in a continuously parametrized family of micro-architectures. More specifically, despite the dimensionality and degree of nonlinearity in the maximum micro von Mises stress within the unit cell, the surrogate models can predict it with less than 5% error at least 90% of the time. Moreover, the largest under-prediction error, which is more critical than the over-prediction error, is typically less than 10%.

42 ENGINEERING↗

DOE OSTI2022

An extrapolation method to remove spurious stress concentration in micromechanical analyses of composites using pixel-based meshes

It is easier to generate a finite element mesh of complicated geometries with voxel-based meshes compared to conformal meshes. However, the local stresses in voxel-based meshes are higher owing to artificial stress concentration caused due to jagged lines that approximate a curve. Here, we propose a method that accounts for this anomaly which would lead to accurate prediction of maximum von Mises stress in voxel-based meshes. The method is based on the variation of reciprocal of stresses and associated numerical gradient. The prediction of the proposed method is compared against the results of a conformal mesh by considering four two-dimensional example problems. The accuracy of the proposed method is verified using different mesh densities, material properties, and integration schemes. The proposed method was also implemented on a conformal mesh to verify that it filters out just the spurious stress concentration and not the realistic stress concentration.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Nanomechanical investigation of the interplay between pore morphology and crack orientation of amorphous silica

Porous amorphous silica (a-SiO 2 ) is of both fundamental and practical interests, as they exhibit a large specific surface area and tunable porous network. However, the brittle nature of a-SiO 2 and the presence of pre-existing cracks at both micro- and nano-scales lead to complex mechanical behavior. Here, we systematically investigate the effects of pre-existing crack and its orientation on the mechanical properties of a-SiO 2 with varying pore shapes using reactive molecular dynamics simulations. We demonstrate that pore shape will primarily influence the Young’s modulus (E) and critical energy release rate (G IC ). We further investigate the impact of pore shape and crack orientation by local characterization of the structural parameters. By defining the high stress and inter-mediate regions, the overall mechanical properties are found to be greatly influenced by the pore shape which can be reflected through the spatial distribution of von Mises stress. Overall, G IC is found to increase with the increase of ligament length (also known as pore wall thickness). Meanwhile, the effect of the pre-existing crack on the crack propagation process is confirmed by analyzing the density distribution evolution. These results highlight the interplay between pore morphology and crack orientation in controlling the fracture behaviors in brittle porous materials.

42 ENGINEERING↗

DOE OSTI2026

Parametric reduced order models for graded lattice structures

Graded lattice structures, characterized by smoothly varying mechanical properties, hold significant promise for optimizing material distribution in advanced engineering applications. However, accurately modeling these structures poses substantial computational challenges due to the continuous geometric variations within their unit cells. Here, to address these challenges, this paper introduces a novel Efficient Reduced Order Model (EROM) that integrates the Matrix Discrete Empirical Interpolation Method (MDEIM) and Discrete Empirical Interpolation Method (DEIM) with polynomial regression to efficiently manage geometric parametrization in lattice structures. Unlike traditional reduced order models (ROMs) that require extensive precomputed libraries for each geometric configuration, our approach enables continuous geometric variations through a flexible algebraic formulation, significantly reducing computational costs while preserving high accuracy. The method constructs projection matrices for individual unit cells that can be efficiently assembled into global systems, leveraging the repetitive nature of lattice structures. Numerical studies demonstrate that our EROM achieves displacement errors below 1% and von Mises stress prediction errors below 4%, coupled with computational speedups exceeding two orders of magnitude compared to full-order simulations. The proposed method's modularity and scalability make it particularly suitable for design optimization and real-time simulation of functionally graded lattice structures, with applications spanning aerospace to biomedical engineering.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗