Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Strain calculation”

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 811 records · Page 45

Behavior of a turbulent boundary layer subjected to sudden transverse strain

Data from two experiments on the development of the components of the Reynolds stress tensor after a sudden application of transverse strain are compared. Computations were based on four different turbulence models: a first-order mixing length model, a second-order two-equation eddy viscosity model, and two second-order Reynolds stress models. The second-order models do not produce dramatic improvements over the simple mixing length model. The Reynolds stress models still need development to represent the physics of shear-strained turbulence well. Finally, it is demonstrated that the assumption of a scalar eddy viscosity, often used in engineering calculations of three-dimensional boundary layers, is quite reasonable.

Higuchi, H.↗

Design of stiffened composite panels with a fracture constraint

An automated procedure for designing minimum-weight composite panels subject to a local damage constraint under tensile loading was developed. A finite element program based on linear elastic fracture mechanics for calculating stress intensity factors (SIF) was incorporated in the design cycle. Panel fracture toughness was obtained by using a strain based criterion. A general purpose mathematical optimization algorithm was used for the weight minimization. Analytical sensitivity derivatives of the SIF employing the adjoint variable technique was used to enhance the computational efficiency of the procedure. Design results for both unstiffened and stiffened plates are presented.

Gurdal, Z.↗

Probabilistic finite elements for fracture mechanics

The probabilistic finite element method (PFEM) is developed for probabilistic fracture mechanics (PFM). A finite element which has the near crack-tip singular strain embedded in the element is used. Probabilistic distributions, such as expectation, covariance and correlation stress intensity factors, are calculated for random load, random material and random crack length. The method is computationally quite efficient and can be expected to determine the probability of fracture or reliability.

Besterfield, Glen↗

Nonlinear analysis and redesign of the mixed-mode bending delamination test

The Mixed Mode Bending (MMB) test uses a lever to simultaneously apply mode I and mode II loading to a split beam specimen. An iterative analysis that accounts for the geometric nonlinearity of the MMB test was developed. The analysis accurately predicted the measured load displacement response and the strain energy release rate, G, of an MMB test specimen made of APC2 (AS4/PEEK). The errors in G when calculated using linear theory were found to be as large as thirty percent in some cases. Because it would be inconvenient to use a nonlinear analysis to analyze MMB data, the MMB apparatus was redesigned to minimize the nonlinearity. The nonlinear analysis was used as a guide in redesigning the MMB apparatus. With the redesigned apparatus, loads were applied through a roller attached to the level and loaded just above the midplane of the test specimen. The redesigned apparatuus has geometric nonlinearity errors of less than three percent, even for materials substantially tougher than APC2. This apparatus was demonstrated by measuring the mixed mode delamination fracture toughness of APC2.

Reeder, J. R.↗

Scale-Similar Models for Large-Eddy Simulations

Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Sarghini, F.↗

Interim Development of new Class B Code Case with Variable Design Lifetime and enhanced design rules to guard against cyclic structural failure mode

This report summarizes updates on the ongoing development of the new American Society of Mechanical Engineers Section III, Division 5, Class B rules to address the potential failure modes of Class B components. An initial design-by-analysis approach and supporting design rules are presented. The elastic-perfectly plastic (EPP) analysis approach has been adopted for primary stress limit and ratcheting check to prevent the component failure against the structural failure mode against primary load and strain accumulation due to cyclic loads. The creep-fatigue damage assessment uses explicitly defined elastic follow-up calculated from the stress concentration region of given component. Damage fraction calculation uses a novel coupled approach to capture the influence of the elastic follow-up and provide adequate conservatism for Class B components. The proposed Class B rules does not use the stress classification approach and uses the combined loads to assess component.

36 MATERIALS SCIENCE↗

Earthquake response of head-mounted equipment in advanced nuclear reactors

The seismic response of safety-related equipment mounted on the head of an advanced reactor, including pumps, control rod drive mechanisms, and reactor monitoring devices, will affect the design and layout of many advanced reactors. High earthquake-induced accelerations in such equipment may challenge their seismic qualification and trigger the need for additional support framing on the reactor head. Base isolation is a design solution that can drastically reduce seismic demands on equipment. This article describes a set of earthquake-simulator experiments conducted on a scale-model of a base-isolated reactor vessel including four representations of head-mounted equipment, with frequencies spanning from 4.5 to 27 Hz. Dynamic responses of the head-mounted equipment, including displacements, accelerations, and strains, were measured in the experiments for three support conditions: conventional, and seismically isolated using single concave Friction Pendulum (SFP) bearings and triple Friction Pendulum (TFP) bearings. Seismic isolation was effective at reducing equipment responses (accelerations, displacements, and strains) with respect to those in the conventionally supported vessel across a range of seismic inputs. Companion numerical studies highlight the accuracy to be expected in the calculation of different response quantities for lightly damped equipment. The importance of characterizing damping in head-mounted, safety-related equipment through physical experiments to support design and risk assessment is made clear through the numerical simulations.

Engineering↗

Prediction of room-temperature multiferroicity in strained MoCr 2 S 6 monolayer

The contrasting d-orbital occupation required for ferroelectricity vs ferromagnetism makes it difficult for their coexistence in two-dimensional materials, especially at high temperature. To resolve this intrinsic contradiction, we propose a layered MoCr 2 S 6 multiferroics by alloying magnetic Cr element into the ferroelectric 1T phase of the MoS 2 matrix. First-principles calculations disclose that a spontaneous symmetry breaking, depending on the Mo atom displacement, leads to a robust ferroelectricity, which coexists with a ferromagnetic order originated from two neighboring Cr atoms. The effect can be further enhanced by tensile strain to bring about a room-temperature multiferroicity. Our results shed new light on the fundamental understanding of multiferroics and display promising applications in spintronics and multistate data storage.

2D materials↗

Chemical control of polarization in thin strained films of a multiaxial ferroelectric: Phase diagrams and polarization rotation

The emergent behaviors in thin films of a multiaxial ferroelectric (FE) due to electrochemical coupling between the rotating polarization and surface ions are explored within the framework of the 2–4 Landau-Ginzburg-Devonshire (LGD) thermodynamic potential combined with the Stephenson-Highland (SH) approach. The combined LGD-SH approach allows us to describe the electrochemical switching and rotation of a polarization vector in a multiaxial ferroelectric film covered by surface ions with a charge density defined by the oxygen pressure. In this work, we calculate phase diagrams, analyze the dependence of polarization components on the applied voltage, and discuss the peculiarities of quasistatic ferroelectric, dielectric, and piezoelectric hysteresis loops in thin strained multiaxial ferroelectric films. The nonlinear surface screening by oxygen ions makes the diagrams very different from the known diagrams of, e.g., strained BaTiO 3 films. Quite unexpectedly, we predict the appearance of ferroelectric reentrant phases. The obtained results point to the possibility to control the appearance and features of ferroelectric, dielectric, and piezoelectric hysteresis in multiaxial FE films covered with surface ions by varying their concentration via the partial oxygen pressure. The LGD-SH description of a multiaxial FE film can be further implemented within the Bayesian optimization framework, paving the way toward predictive materials optimization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ferroelectric HfO 2 and the importance of strain

Ferroelectric oxides based on HfO 2 show tremendous promise for the next generation of memory and logic devices. The ferroelectric polymorph is one of several that can be derived from the high symmetry cubic fluorite structure of HfO 2 . A single grain of HfO 2 may consist of a coherent mixture of multiple orientational and translational variants of different polymorphs. Here, we use symmetry-adapted strain-order parameters to elucidate the relationship between the different HfO 2 polymorphs and their symmetrically equivalent variants. We use first-principles electronic structure methods to identify minimum energy pathways and map them in subspaces of the symmetry-adapted strain order parameters. We next investigate the atomic structure of domain boundaries that separate coexisting variants of ferroelectric HfO 2 . Further, we rely on Gibbsian excess quantities and a precise specification of mechanical boundary conditions to describe the thermodynamic properties of domain boundaries. Our first-principles calculations show that the O and Hf shuffle arrangement within a domain boundary is closely related to the intermediate shuffle patterns of the homogeneous pathways between ferroelectric variants. Furthermore, the preferred structure within a boundary is very sensitive to local strain constraints imposed by the adjacent ferroelectric variants, leading to highly anisotropic domain boundary energies.

36 MATERIALS SCIENCE↗

Emergent multiferroism with magnetodielectric coupling in EuTiO 3 created by a negative pressure control of strong spin-phonon coupling

Negative pressure has emerged as a powerful tool to tailor the physical properties of functional materials. However, a negative pressure control of spin-phonon coupling for engineering magnetism and multiferroicity has not been explored to date. Here, using uniform three-dimensional strain-induced negative pressure in nanocomposite films of (EuTiO 3 ) 0.5 :(MgO) 0.5 , we demonstrate an emergent multiferroicity with magnetodielectric coupling in EuTiO 3 , matching exactly with density functional theory calculations. Density functional theory calculations are further used to explore the underlying physics of antiferromagnetic-paraelectric to ferromagnetic-ferroelectric phase transitions, the spin-phonon coupling, and its correlation with negative pressures. The observation of magnetodielectric coupling in the EuTiO 3 reveals that an enhanced spin-phonon coupling originates from a negative pressure induced by uniform three-dimensional strain. Our work provides a route to creating multiferroicity and magnetoelectric coupling in single-phase oxides using a negative pressure approach.

36 MATERIALS SCIENCE↗

The direct effects of strain on burning rates of composite solid propellants

A mathematical model is developed to predict burn rate augmentation due to strain in a composite solid propellant. The model assumes the effect is due to the ability of the flame to penetrate the small fissures and voids that form when a propellant is strained. The number and size of these fissures is obtained by applying a flaw propagation analysis to randomly distributed flaws that form when the binder-oxidizer particle bonds break under stress. A flame height is calculated with Summerfield's burn rate equation and is used to compute the burn rate augmentation based upon the additional burn area created when the flame penetrates the fissures. Comparisons are made with data obtained from published sources. The existence of threshold pressure and strains, above which augmentation occurs, is verified although the model predicts a lower threshold pressure and higher threshold strain than expected. Further results and applications of the model are discussed.

Langhenry, M. T.↗

Interfacial Reactivity-Triggered Oscillatory Lattice Strains of Nanoalloys

Understanding the structure evolution of nanoalloys under reaction conditions is vital to the design of active and durable catalysts. Herein, we report an operando measurement of the dynamic lattice strains of dual-noble-metal alloyed with an earth-abundant metal as a model electrocatalyst in a working proton-exchange membrane fuel cell using synchrotron high-energy X-ray diffraction coupled with pair distribution function analysis. Additionally, the results reveal an interfacial reaction-triggered oscillatory lattice strain in the alloy nanoparticles upon surface dealloying. Analysis of the lattice strains with an apparent oscillatory irregularity in terms of frequency and amplitude using time-frequency domain transformation and theoretical calculation reveals its origin from a metal atom vacancy diffusion pathway to facilitate realloying upon dealloying. This process, coupled with surface metal partial oxidation, constitutes a key factor for the nanoalloy's durability under the electrocatalytic oxygen reduction reaction condition, which serves as a new guiding principle for engineering durable or self-healable electrocatalysts for sustainable fuel cell energy conversion.

Wu, Zhi-Peng [State University of New York (SUNY),↗

Predicting Failure Using Deep Learning SAND Report

Accurate prediction of ductile failure is critical to Sandia’s NW mission, but the models are computationally heavy. The costs of including high-fidelity physics and mechanics that are germane to the failure mechanisms are often too burdensome for analysts either because of the person-hours it requires to input them or because of the additional computational time, or both. In an effort to deliver analysts a tool for representing these phenomena with minimal impact to their existing workflow, our project sought to develop modern data-driven methods that would add microstructural information to business-as-usual calculations and expedite failure predictions. The goal is a tool that receives as input a structural model with stress and strain fields, as well as a machine-learned model, and output predictions of structural response in time, including failure. As such, our project spent substantial time performing high-fidelity, three-dimensional experiments to elucidate materials mechanisms of void nucleation and evolution. We developed crystal-plasticity finite-element models from the experimental observations to enrich the findings with fields not readily measured. We developed engineering length-scale simulations of replicated test specimens to understand how the engineering fields evolve in the presence of fine-scale defects. Finally, we developed deep learning convolutional neural networks, and graph-based neural networks to encode the findings of the experiments and simulations and make forward predictions in time for structural performance. This project demonstrated the power of data-driven methods for model development, which have the potential to vastly increase both the accuracy and speed of failure predictions. These benefits and the methods necessary to develop them are highlighted in this report. However, many challenges remain to implementing these in real applications, and these are discussed along with potential methods for overcoming them.

97 MATHEMATICS AND COMPUTING↗

Thermophysical characterization of composite materials under transient heating conditions

Thermophysical property measurements were made under transient heating conditions on several materials being considered for use in SCOUT rocket motors. The materials included were ATJ graphite, MX 2600 silica phenolic, FM 5272 cellulose phenolic, and two carbon-carbon composites: CARBITEX 700 and RPP-4. The ATJ was included as a reference or base line material to check performance of the transient tests as it was not expected to be sensitive to heating rate. Measurements included in the program were thermal conductivity, strength, compressive stress-strain (carbon-carbon only), thermal expansion and the effective thermal expansion under partially restrained conditions. Development of this latter measurement was a major part of the program. It consisted of partially restraining the expansion of a specimen as it was heated, measuring the load and strain which occurred (together with a simultaneous modulus determination by superimposing a small cyclic load) and using these quantities to calculate what the effective thermal expansion would have to be to produce the observed stress and deformation. For materials which are sensitive to heating rate, such as reinforced phenolics, it was believed that this would provide a more realistic determination of the thermal expansion as it more nearly simulates the conditions experienced in end use.

Roetling, J.↗

Stress-intensity factors for internal surface cracks in cylindrical pressure vessels

The paper presents stress-intensity factors for a wide range of semi-elliptical surface cracks on the inside of pressurized cylinders. The ratio of crack depth to crack length ranged from 0.2 to 1; the ratio of crack depth to wall thickness ranged from 0.2 to 0.8; and the ratio of wall thickness to vessel radius was 0.1 to to 0.25. The stress-intensity factors were calculated by a three-dimensional finite-element method using singularity elements along the crack front and linear-strain elements elsewhere. An equation for the stress-intensity factors was obtained which applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was also made between the results and other analyses of internal surface cracks in cylinders. The results from a boundary-integral equation method were in agreement and those from another finite-element method were in fair agreement (+ or - 8 percent) with the results.

Newman, J. C., Jr.↗

On composites with periodic structure

The overall moduli of a composite with an isotropic elastic matrix containing periodically distributed (anisotropic) inclusions or voids, can be expressed in terms of several infinite series which only depend on the geometry of the inclusions or voids, and hence can be computed once and for all for given geometries. For solids with periodic structures these infinite series play exactly the same role as does Eshelby's tensor for a single inclusion or void in an unbounded elastic medium. For spherical and circular-cylindrical geometries, the required infinite series are calculated and the results are tabulated. These are then used to estimate the overall elastic moduli when either the overall strains or the overall stresses are prescribed, obtaining the same results. These results are compared with other estimates and with experimental data. It is found that the model of composites with periodic structure yields estimates in excellent agreement with the experimental observations.

Nemat-Nasser, S.↗

Convergence of strain energy release rate components for edge-delaminated composite laminates

Strain energy release rates for edge delaminated composite laminates were obtained using quasi 3 dimensional finite element analysis. The problem of edge delamination at the -35/90 interfaces of an 8-ply composite laminate subjected to uniform axial strain was studied. The individual components of the strain energy release rates did not show convergence as the delamination tip elements were made smaller. In contrast, the total strain energy release rate converged and remained unchanged as the delamination tip elements were made smaller and agreed with that calculated using a classical laminated plate theory. The studies of the near field solutions for a delamination at an interface between two dissimilar isotropic or orthotropic plates showed that the imaginary part of the singularity is the cause of the nonconvergent behavior of the individual components. To evaluate the accuracy of the results, an 8-ply laminate with the delamination modeled in a thin resin layer, that exists between the -35 and 90 plies, was analyzed. Because the delamination exists in a homogeneous isotropic material, the oscillatory component of the singularity vanishes.

Raju, I. S.↗