Elastic-plastic torsion problem for strain-hardening materials
Successive approximation procedure for solving elastic-plastic problem for torsion of solid prismatic bar made of strain hardening material
SEARCH · Search NASA
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.
Successive approximation procedure for solving elastic-plastic problem for torsion of solid prismatic bar made of strain hardening material
Phase diagrams supported by density functional theory methods can be crucial for designing high-entropy alloys that are subset of multi-principal-element alloys. We present phase and property analysis of quinary (MoW) x Zr y (TaTi) 1-x-y refractory high-entropy alloys from combined Calculation of Phase Diagram (CALPHAD) and density-functional theory results, supplemented by molecular dynamics simulations. Both CALPHAD and density-functional theory analysis of phase stability indicates a Mo-W-rich region of this quinary has a stable single-phase body-centered-cubic structure. We report first quinary composition from Mo-W-Ta-Ti-Zr family of alloy with pseudo-elastic behavior, i.e., hysteresis in stress-strain. Our analysis shows that only Mo-W-rich compositions of Mo-W-Ta-Ti-Zr, i.e., Mo+W, show reproducible hysteresis in stress-strain responsible for pseudo-elastic behavior. The (MoW) 85 Zr 7.5 (TaTi) 7.5 was down-selected based on temperature-dependent phase diagram analysis and molecular dynamics simulations predicted elastic behavior that reveals twinning-assisted pseudoelastic behavior. While mostly unexplored in body-centered-cubic crystals, twinning is a fundamental deformation mechanism that competes against dislocation slip in crystalline solids. Furthermore, this alloy shows identical cyclic deformation characteristics during uniaxial < 100 > loading, i.e., the pseudoelasticity is isotropic in loading direction. Additionally, a temperature increase from 77 to 1,500 K enhances the elastic strain recovery in load-unload cycles, offering possibly control to tune the pseudoelastic behavior.
Application of isotropic pressure or uniaxial strain alters the elastic properties of materials; sufficiently large strains can drive structural transformations. Linear elasticity describes stability against infinitesimal strains, while nonlinear elasticity describes the response to finite deformations. Here, it was previously shown that uniaxial strain along [100] drives refractory metals and alloys towards mechanical instabilities. These include an extensional instability, and a symmetry-breaking orthorhombic distortion caused by a Jahn-Teller-Peierls instability that splays the cubic lattice vectors. Here we analyze these transitions in depth. Eigenvalues and eigenvectors of the Wallace tensor identify and classify linear instabilities in the presence of strain. We show that both instabilities are discontinuous, leading to discrete jumps in the lattice parameters. We provide physical intuition for the instabilities by analyzing the changes in first-principles energy, stress, bond lengths, and angles upon application of strain. Electronic band structure calculations show differential occupation of bonding and antibonding orbitals, driven by the changing bond lengths and leading to the structural transformations. Strain thresholds for these instabilities depend on the valence electron count.
Numerical implementation of an isotropic creep-damage model for concrete in multiphysics object-oriented simulation environment (MOOSE) finite element framework is presented in this paper. The constitutive model considers the combined effect of instantaneous and delayed strains on damage propagation. The implementation considers creep using generalized Maxwell or Kelvin-Voigt models. Using strain splitting assumptions, the total mechanical strains are split into elastic and creep components. Damage is considered to evolve as a function of the elastic and creep strains. This work considers damage as a function of fracture energy using the characteristic length of each finite element. This approach preserves the energy release rate of each element and avoids vanishing energy dissipation as the mesh is refined. Further, a creep-damage parameter is used to quantify the effect of creep strain on damage. The model is tested against published results on notched three-point bending specimens involving non-linear creep and predicts that about a third of the creep strain contributes towards damage for the experiments simulated. Results show that the proposed framework has predictive capabilities, and the model can be extended for more complex systems.
Contact problem for elastic foundations subjected to bending strains
The solution of small strain elastic-plastic stress analysis problems by the p-version of the finite element method is discussed. The formulation is based on the deformation theory of plasticity and the displacement method. Practical realization of controlling discretization errors for elastic-plastic problems is the main focus. Numerical examples which include comparisons between the deformation and incremental theories of plasticity under tight control of discretization errors are presented.
The strain formulation in elasticity and the compatibility condition in structural mechanics have neither been understood nor have they been utilized. This shortcoming prevented the formulation of a direct method to calculate stress. We have researched and understood the compatibility condition for linear problems in elasticity and in finite element analysis. This has lead to the completion of the method of force with stress (or stress resultant) as the primary unknown. The method in elasticity is referred to as the completed Beltrami-Michell formulation (CBMF), and it is the integrated force method (IFM) in structures. The dual integrated force method (IFMD) with displacement as the primary unknown has been formulated. IFM and IFMD produce identical responses. The variational derivation of the CBMF yielded the new boundary compatibility conditions. The CBMF can be used to solve stress, displacement, and mixed boundary value problems. The IFM in structures produced high-fidelity response even with a modest finite element model. The IFM has influenced structural design considerably. A fully utilized design method for strength and stiffness limitation has been developed. The singularity condition in optimization has been identified. The CBMF and IFM tensorial approaches are robust formulations because of simultaneous emphasis on the equilibrium equation and the compatibility condition.
Diamond is not only the hardest material in nature, but is also an extreme electronic material with an ultrawide bandgap, exceptional carrier mobilities, and thermal conductivity. Straining diamond can push such extreme figures of merit for device applications. We microfabricated single-crystalline diamond bridge structures with ~1 micrometer length by ~100 nanometer width and achieved sample-wide uniform elastic strains under uniaxial tensile loading along the [100], [101], and [111] directions at room temperature. We also demonstrated deep elastic straining of diamond microbridge arrays. The ultralarge, highly controllable elastic strains can fundamentally change the bulk band structures of diamond, including a substantial calculated bandgap reduction as much as ~2 electron volts. Our demonstration highlights the immense application potential of deep elastic strain engineering for photonics, electronics, and quantum information technologies.
An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is ideally suited to isotropically hardening Prandtl-Reuss materials. Further, the formulation is given in a manner which allows any conventional finite element program, for 'small strain' elastic-plastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension. The paper closes with a unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures. Further, a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain, and the inadequacies of some of these are commented upon.
Strain redistribution corrections were developed for a simplified inelastic analysis procedure to economically calculate material cyclic response at the critical location of a structure for life prediction proposes. The method was based on the assumption that the plastic region in the structure is local and the total strain history required for input can be defined from elastic finite-element analyses. Cyclic stress-strain behavior was represented by a bilinear kinematic hardening model. The simplified procedure predicts stress-strain response with reasonable accuracy for thermally cycled problems but needs improvement for mechanically load-cycled problems. Neuber-type corrections were derived and incorporated in the simplified procedure to account for local total strain redistribution under cyclic mechanical loading. The corrected simplified method was used on a mechanically load-cycled benchmark notched-plate problem. The predicted material response agrees well with the nonlinear finite-element solutions for the problem. The simplified analysis computer program was 0.3% of the central processor unit time required for a nonlinear finite-element analysis.
Explore the source record for details and available documents.
Extended x-ray absorption fine structure (EXAFS) measurements of single-crystal Ge/GeSn radial heterostructure nanowires are used to examine the effects of composition on both short-range order (SRO) and longer-range disorder in GeSn alloys. GeSn has prompted significant interest because it can achieve a direct band gap for sufficient Sn concentrations beyond the equilibrium solid solubility limit in an all-group IV system. Short-range order in this material is particularly interesting as it has been predicted to affect the band gap independent of average composition or strain effects. By independently controlling the Sn composition and GeSn thickness during chemical vapor deposition of misfitting GeSn shells around ultrathin, elastically compliant, Ge core nanowires, the elastic misfit strain in the GeSn is minimized for Sn compositions over the studied range ≈Ge 0.96 Sn 0.04 to Ge 0.88 Sn 0.12 . The degree of SRO was found to decrease with increasing Sn composition. Additionally, damping of the EXAFS signal was observed as the Sn content increased, particularly for increasingly distant atomic shells about the absorbing atom, even for scattering paths not involving Sn atoms. This result is quantified as an increase in the mean-squared relative displacement parameters of the shells. These measurements reveal the accommodation of local strain due to the presence of the highly size-mismatched Sn atoms in the Ge diamond cubic lattice (≈14%), which may have effects on the band structure of the material in addition to the influence of short-range atomic order. Comparison among the nanowire samples allows for calculation of the topological rigidity parameter, a ∗∗ , for the first-neighbor bond lengths. Furthermore, these exhibit chemically distinct values for Ge-Ge, Ge-Sn, and Sn-Sn, and they are consistent with the value a ∗∗ = 0.75 ± 0.07 confirming the general applicability of the model to alloys with both large amounts of natural misfit strain and the potential for short-range order.
An elastic large displacement thick-shell mixed finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at the element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses is included in the mixed variational functional. The nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the nonlinear elastic equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to computationally simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.
An elastic large displacement thick-shell hybrid finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at he element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses in included in the mixed variational functional. Nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to numerically simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.
The development of two new shell finite elements for applications to large deflection problems is considered. The elements in question are doubly curved and of triangular and quadrilateral planform. They are restricted to small strains of elastic materials, and can accommodate large rotations. The elements described, which are based on relatively simple linear elements, make use of a new displacement function approach specifically designed for strongly nonlinear problems. The displacement function development for nonlinear applications is based on certain beam element formulations, and the strain-displacement equations are of a shallow shell type. Additional terms were included in these equations in an attempt to avoid the large errors characteristic of shallow shell elements in certain types of problems. An incremental nonlinear solution procedure specifically adopted to the element formulation was developed. The solution procedure is of combined incremental and total Lagrangian type, and uses a new updating scheme. A computer program was written to evaluate the developed formulations. This program can accommodate small element groups in arbitrary arrangements. Two simple programs were successfully solved. The results indicate that this new type of element has definite promise and should be a fruitful area for further research.
Geodetic work done over the past century is examined to investigate block motion in areas of intense tectonic deformation, with special attention to geodetic releveling data obtained for Japan. Problems in interpreting strainmeter and tiltmeter to study the behavior of crustal blocks and block boundaries are discussed. Block dimensions of 5 to 50 km seem to occur frequently in regions of intense tectonic activity. Block boundaries are often considerably weaker than contiguous crustal blocks, resulting in a concentration of strains at boundaries and decrease of strains within blocks. Block to block variations in tilt phase and magnitude are observed. There is some evidence that block boundaries may exhibit strain-dependent elastic properties or respond viscoelastically, possibly accounting for the slow transmission of tectonic deformation reported in Japan. If nonlinear behavior is characteristic of regions fragmented by crustal blocks, it will generally not be possible to apply a site correction factor based on the observed distortion of tidal or seismic strains in interpreting secular strains.
The self-consistent (SC) theory is the most commonly used mean-field homogenization method to estimate the mechanical response behavior of polycrystals based on the knowledge of the properties and orientation distribution of constituent single-crystal grains. The original elastic SC method can be extended to thermo-elasticity by adding a stress-free strain to an elastic constitutive relation that expresses stress as a linear function of strain. With the addition of this independent term, the problem remains linear. Although the thermo-elastic self-consistent (TESC) model has important theoretical implications for the development of self-consistent homogenization of non-linear polycrystals, in this paper, we focus on TESC applications to actual thermo-elastic problems involving non-cubic (i.e. thermally anisotropic) materials. To achieve this aim, here we provide a thorough description of the TESC theory, which is followed by illustrative examples involving cooling of polycrystalline non-cubic metals. The TESC model allows studying the effect of crystallographic texture and single-crystal elastic and thermal anisotropy on the effective thermo-elastic response of the aggregate and on the internal stresses that develop at the local level.
Nanoindentation is a critical technique to probe mechanical properties at the micrometer and sub-micrometer scales, accompanied by challenges from indentation size effect, pile-up/sink-in effect, and strain rate sensitivity. In this work, different nanoindentation techniques have been employed to explore Ni-based concentrated solid solution alloys (CSAs) with the addition of 3d transition metal elements including Co, Cr, Mn, and Fe, including unique single-crystal Ni, NiCo, NiFe, Ni 80 Cr 20 , and NiCoFeCr samples with (100) surfaces. A procedure of nanoindentation tests and data analysis/correction have been developed, and a data set of hardness, elastic modulus, strain rate sensitivity, and activation volume for Ni-based CSAs are provided, including the less explored binary alloys such as Ni 80 Cr 20 and Ni 80 Mn 20 . The results show that the type of alloying elements is more critical than the number of elements in strengthening: Co does not provide strengthening in NiCo, while Cr, Mn, and Fe are effective strengthening elements. Cr is the most effective among all the 3d transition metal elements. Furthermore, atomic-level lattice distortion is responsible for the strengthening and the role of stacking fault energy is insignificant in Ni-based CSAs at room temperature. In summary, nanoindentation shows increasing promise as a reliable and fast tool to provide comprehensive mechanical information for new alloy design and development.