SEARCH · Search NASA
Results for “Finite deformation theory”
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.
Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation
To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.
Phase-Field Model Incorporating Large Inelastic Strain with Application to the Oxidation of High-Temperature Coating Systems
Structural metallic and ceramic materials at high temperatures often suffer from creep, fatigue and oxidation, which typically involves local large inelastic strain. The large strain is usually approached by Lagrangian-description-based finite deformation theory. Here we develop an efficient, Eulerian-description-based, Incremental Realization of Inelastic Deformation (IRID) algorithm to deal with large inelastic strain with small elastic strain. We apply this algorithm to simulate high temperature oxidation with severe volume expansion using a phase-field approach. The oxidation and deformation of bond coat can lead to spallation and failure of the thermal barrier coatings (TBC) and environmental barrier coatings (EBC), which is a critical concern in design of these coatings. We simulate the growth of thermally grown oxide (TGO) in these coating systems, with phase transformations and microstructure evolution incorporated. The simulation results are compared with relevant experimental reports in the literature.
Exact solution of burst pressure for thick-walled pipes using the flow theory of plasticity
This paper develops an exact solution of burst pressure for defect-free, thick-walled pipes using the flow theory of plasticity in terms of the average shear stress yield criterion (or Zhu-Leis yield criterion). The pipe steel is assumed to obey the power-law strain hardening rule, and large plastic deformation is described by the finite strain theory. On this basis, internal pressure is obtained as a power series function of the effective strains on the inside and outside surfaces of the thick-walled pipe with use of the second Bernoulli numbers. At burst failure, the Zhu-Leis flow solution of burst pressure is determined as a power series solution, and the burst effective strains, the burst effective stresses, and the burst pressure are functions of the diameter ratio (D o / D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). Similarly, the Tresca and von Mises flow solutions of burst pressure are also determined as a power series solution in terms of the Tresca and von Mises yield criteria. A general closed-from exact solution of burst pressure is then proposed for the three yield criteria, and the results showed that the proposed exact solution matches well with the power series solution for each yield criterion. Moreover, the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction that agrees well with the finite element analysis results of burst pressure for thick-walled pipes. Furthermore, two datasets of full-scale burst tests are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin and thick-walled pipes.
Consistent linearization of micromorphic continuum theories
We outline the procedure of consistent linearization and apply it to the micromorphic, microstretch, and micropolar theories of continua. This yields tractable linear theories for nonlinearly elastic microstructured materials undergoing finite deformations. The results may be readily utilized in computational mechanics, stability and bifurcation analyses, and small-deformation problems in the context of these types of continua. Our results generalize those existing in current literature and facilitate their recovery upon incorporating appropriate kinematic and constitutive assumptions.
A large deformation multiphase continuum mechanics model for shock loading of soft porous materials
A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.
Exact solutions of burst pressure for thick-walled cylinders in power-law strain hardening steels
Exact solutions of burst pressure for defect-free, thick-walled cylindrical pressure vessels with capped ends are developed using the flow theory of plasticity in terms of the Tresca, von Mises, and Zhu-Leis yield criteria. The material is assumed to obey a power-law strain hardening rule, and the finite strain theory is adopted to describe large plastic deformation. On this basis, internal pressure is obtained as a complex polynomial function of effective strains on the inside and outside surfaces of the pressure vessel with the use of the Bernoulli numbers. At burst failure, the effective strains and stresses on the inside and outside surfaces are first determined, and then three flow solutions of burst pressure are obtained as a power series function of the diameter ratio (D o /D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). The power series solution is confirmed to agree well with the exact flow solution of burst pressure in a closed-form for each yield criterion, and the results show that the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction of burst pressure. In conclusion, two sets of full-scale burst test data are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin-walled pipes and thick-walled cylindrical pressure vessels.
A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping
A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.
Effective field theories for collective excitations of atomic nuclei
Collective modes emerge as the relevant degrees of freedom that govern low-energy excitations of atomic nuclei. These modes—rotations, pairing rotations, and vibrations—are separated in energy from non-collective excitations, making it possible to describe them in the framework of effective field theory. Rotations and pairing rotations are the remnants of Nambu–Goldstone modes from the emergent breaking of rotational symmetry and phase symmetries in finite deformed and finite superfluid nuclei, respectively. The symmetry breaking severely constrains the structure of low-energy Lagrangians and thereby clarifies what is essential and simplifies the description. The approach via effective field theories exposes the essence of nuclear collective excitations and is defined with a breakdown scale in mind. This permits one to make systematic improvements and to estimate and quantify uncertainties. Effective field theories of collective excitations have been used to compute spectra, transition rates, and other matrix elements of interest. In particular, predictions of the nuclear matrix element for neutrinoless double beta decay then come with quantified uncertainties. This review summarizes these results and also compares the approach via effective field theories to well-known models and ab initio computations.
Wormholes, branes and finite matrices in sine dilaton gravity
We compute the double trumpet in sine dilaton gravity via WdW quantization. The wormhole size is discretized. The wormhole amplitude matches the spectral correlation of a finite-cut matrix integral, where matrices have large but finite dimensions. This strongly suggests an identification of the sine dilaton gravity theory with the q-deformed JT gravity matrix integral. At the very least, it captures all universal content of that matrix model. The disk decomposes into the physical (gauge invariant) solutions of the WdW equation, which are trumpets with discrete sizes. This decomposition modifies the usual no-boundary wavefunction to a normalizable one in sine dilaton gravity.
Instanton Effects in Euclidean vacuum, Real Time Production, and in the Light-front Wave Functions
Nontrivial topological structures of non-Abelian gauge fields were discovered in the 1970s. Instanton solutions, describing vacuum tunneling through topological barriers, have fermionic zero modes which are at the origin of ‘t Hooft effective Lagrangian. In the 1980s, instanton ensembles have been used to explain chiral symmetry breaking. In the 1990s, a large set of numerical simulations were performed deriving Euclidean correlation functions. The special role of scalar diquarks in nucleons and color superconductivity in dense quark matter has been elucidated. In these lectures, we discuss further developments of physics related to gauge topology. We show that the instanton–anti-instanton “streamline” configurations describe “sphaleron transitions” in high-energy collisions, which result in production of hadronic clusters with nontrivial topological/chiral charges. (They are not yet observed, but discussions of dedicated experiments at the LHC and RHIC are ongoing.) Another new direction is instanton effects in hadronic spectroscopy, both in the rest frame and on the light front. We will discuss their role in central and spin-dependent potentials, form factors and antiquark nuclear “sea”. Finally, we summarize the advances in the semiclassical theory of deconfinement, and chiral phase transitions at finite temperature, in QCD, and in some of its “deformed” versions.
Life at the Landau Pole
If a quantum field theory has a Landau pole, the theory is usually called ‘sick’ and dismissed as a candidate for an interacting UV-complete theory. In a recent study on the interacting 4d O(N) model at large N, it was shown that at the Landau pole, observables remain well-defined and finite. In this work, I investigate both relevant and irrelevant deformations of the said model at the Landau pole, finding that physical observables remain unaffected. Apparently, the Landau pole in this theory is benign. As a phenomenological application, I compare the O(N) model to QCD by identifying $Λ_{\overline{MS}}$ with the Landau pole in the O(N) model.
Real-time lattice gauge theory actions: Unitarity, convergence, and path integral contour deformations
Not Available
A spatially-resolved model of neutron-irradiated tungsten coupling stochastic cluster dynamics and finite deformation plasticity
Structural materials used in nuclear reactors face severe degradation in mechanical properties, such as hardening and embrittlement. At the microscopic scale, this occurs due to creation and accumulation of irradiation-induced defects and their interaction with system dislocations. Although techniques exist which can model evolution of irradiation defects, for instance kinetic transport theory-based models, their interaction with mechanical deformation of the bulk material has not been investigated extensively. In this work, we demonstrate a novel spatially-resolved multiscale coupling between microscopic irradiation defect evolution, modeled using Stochastic Cluster Dynamics (SCD) and macroscopic mechanical deformation modeled using a finite-deformation plasticity model. SCD is used to determine the statistically averaged defect cluster spacing, dependent on operating conditions such as irradiation dose and temperature. This acts as an initial condition that governs the critical resolved shear stress of dislocation glide in the macroscopic plasticity model. This framework is used to predict mechanical behavior in post-mortem test of irradiated Tungsten samples, which has found its importance as structural material used in nuclear reactors. The results obtained using the coupled approach are in good agreement with experimental data of uniaxial tension tests. The model is able to capture the effect of temperature and irradiation dose on the material hardening. Two methods are proposed to estimate hardness – using Tabor's Law relating uniaxial yield stress to hardness and from flat-punch simulations. The results are in reasonable agreement with hardness data from micro-indentation experiments of irradiated Tungsten samples. Finally, the model is also able to reveal microstructural details such as spatial variation in defect density and local stress.
An implementation of the phase-field model based on coupled thermomechanical finite element solvers for large-strain twinning, explicit dynamic fracture and the classical Stefan problem
The implementation of a phase-field model in finite elements usually requires significant expertise and involves the development of a user element with additional degrees of freedom. An alternative implementation of the phase-field model within a thermo-mechanical finite element simulation package was presented in (Cho et al 2012 Int. J. Solids Struct. 49 1973–1992), where the phase-field variable is treated as the temperature degree of freedom. However, this approach has only been used for small strain phase-field modelling of martensitic transformations and quasistatic phase-field modelling of fracture. Here, we present a phase-field finite element implementation via the temperature degree of freedom for several additional cases from the literature: (i) the large-strain phase-field description of deformation twinning presented in (Clayton and Knap 2011 Physica D 240 841–858), (ii) phase-field description of brittle fracture with inertial effects based on the theory from (Molnár and Gravouil 2017, Finite Elem. Anal. Des. 130 27–38) and (Miehe et al 2010 Int. J. Numer. Methods Eng. 83 1273–1311) and (iii) the classical Stefan problem of solidification presented in (Mackenzie and Robertson 2002 J. Comput. Phys. 181 526–544). The last problem involves the temperature and phase-field variables as unknowns.
Modeling of surface mechanical behaviors of soft elastic solids: theory and examples
Surfaces of soft solids can have significant surface stress, extensional modulus and bending stiffness. Previous theoretical studies have usually examined cases in which both the surface stress and bending stiffness are constant, assuming small deformation. In this work we consider a general formulation in which the surface can support large deformation and carry both surface stresses and surface bending moments. We show that the large deformation theory can be reduced to the classical linear theory (Shuttleworth equation). We obtain exact solutions for problems of an inflated cylindrical shell and bending of a plate with a finite thickness. Our analysis illustrates the different manners in which surface stiffening and surface bending stabilize these structures. We discuss how the complex surface constitutive behaviors affect the stress field of the bulk. Our calculation provides insights into effects of strain-dependent surface stress and surface bending in the large deformation regime, and can be used as a model to implement surface finite elements to study large deformation of complex structures.
Line operators in Chern-Simons-Matter theories and Bosonization in Three Dimensions II: Perturbative analysis and all-loop resummation
We study mesonic line operators in Chern-Simons theories with bosonic or fermionic matter in the fundamental representation. In this paper, we elaborate on the classification and properties of these operators using all loop resummation of large N perturbation theory. We show that these theories possess two conformal line operators in the fundamental representation. One is a stable renormalization group fixed point, while the other is unstable. They satisfy first-order chiral evolution equations, in which a smooth variation of the path is given by a factorized product of two mesonic line operators. The boundary operators on which the lines can end are classified by their conformal dimension and transverse spin, which we compute explicitly at finite ’t Hooft coupling. We match the operators in the bosonic and fermionic theories. Finally, we extend our findings to the mass deformed theories and discover that the duality still holds true.
Design and Analysis of an Advanced Three-Point Bend Test Approach for Miniature Irradiated Disk Specimens
Candidate tungsten armor materials in a magnetic confinement fusion device must be able to withstand thermal variation that leads to internal stresses caused by the impinging heat load. In addition, the thermomechanical properties of these materials are degraded by irradiation-induced defect accumulation. Fission reactor–based irradiation data are used to predict the fusion neutron damage and property change. This study examines the motivation and design of a custom-designed three-point bend test for neutron-irradiated disk specimens that are 3 mm in diameter to be able to define the flexural strength of advanced tungsten materials, alloys, and composites—and to the extent that embrittlement occurs after neutron irradiation. The theory provided shows a calculation for the flexural deflection and shear deflection due to the small-geometry constraints. A finite element deformation analysis is performed to evaluate the mechanical stress field of disk bend specimens. The stress values above 80% of the maximum stress are concentrated in 2.4 mm of the 3.0-mm length of the centerline across the tungsten disk diameter. A bend test fixture has been designed and fabricated to enable testing of these specimens with precisely engineered tolerance and minimal machine compliance. This fixture will be able to be placed inside a universal testing frame at elevated temperatures for the mechanical property evaluation of future neutron-irradiated disk specimens.