Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite strains”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗
DOE OSTI2024

Pretest Modeling A Spent Nuclear Fuel Seismic Shake Test

The U.S. Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is planning to conduct a series of full-scale shake table tests to simulate hypothetical earthquake conditions and record the response of surrogate spent nuclear fuel (SNF) assemblies in a dry canister storage system mockup. The shake table motions will represent a range of hypothetical earthquake conditions at hypothetical locations in the continental U.S. to generally define the range of mechanical loads that SNF can be expected to experience during extended dry storage periods. This paper describes the pretest predictions made with LS-DYNA models of the mockup storage systems. The test will use two dry storage system configurations, a mockup vertical concrete cask (VCC) and a mockup horizontal storage module (HSM). The test will use a production-quality canister and basket. Within the canister will be four instrumented fuel assemblies with fuel rods containing surrogate mass and 28 instrumented dummy assemblies that are intended to match the mass and outer dimension of a fuel assembly. The finite element models include models of the VCC and HSM on the shake table to calculate the system level dynamic responses and separate single fuel assembly models to calculate stress and strain in fuel assembly components. Both types of models include nonlinear behavior like rod-to-rod contact and the ability for VCC’s to rock and slide. This paper presents the expected response of the VCC, HSM, and fuel assemblies to the shake table testing that is planned to start in April of 2024. The earthquake conditions represent seismic hazards in the 2,000-to-20,000-year return period range. The test data is expected to confirm the expectation that fuel rod cladding will remain intact, fuel assembly structural components like guide tubes will remain intact, and no significant VCC sliding or tipping will occur in the range of conditions to be tested with the shake table.

Klymyshyn, Nicholas A.↗
DOE OSTI2025

Coupled momentum balance and phase-field solver with fenicsx module

Code solves momentum balance and phase-field equations simultaneously. The differential equations are solved on a discretized domain with appropriate boundary and initial conditions using finite element method. Primary purpose of the code is to simulate brittle fracture under dynamic loading. Constitutive equations are that of linear elasticity with degradation of stress due to fracture. Small strain formulation is used.

Zecevic, Milovan↗
DOE OSTI2024

Optimal Geometry for Focused Ion Beam-Milled Samples for Direct-Pull Micro-Tensile Testing Performed In Situ in a Scanning Electron Microscope

A thorough procedure was developed to efficiently manufacture dogbone samples using focused ion beam (FIB) milling for micro-tensile testing. A Bruker PI 89 PicoIndenter, Billerica, MA, USA, was used as a case study, although the analysis and results are applicable to other micro-mechanical testing systems capable of mounting a standard, Ø12.7 mm × Ø3.2 mm pin, scanning electron microscopy (SEM) pin stub (Ted Pella, Redding, CA, USA). Nine dogbones were made from an Fe-45Cu alloy additively manufactured using powder-fed laser-directed energy deposition (DED-LB). Testing showed that fracture was confined to the gauge section for all dogbones and that the fracture mode, ductile vs. brittle, was entirely dependent on the grain orientation relative to the loading direction. The analysis showed that the measured plastic strain to failure can vary from >11% (optimal geometry) to <1% (non-optimal geometry) in micro-tensile testing of high-tensile-strength (>1 GPa) metallic materials. Subsequently, a finite element analysis (FEA) was conducted to identify the improved dogbone geometries. A total of ten thousand dogbone geometries were tested, and their dimensions were defined by a set of four adjustable parameters (corner radius, load surface angle, load surface length, and dogbone head length). The gauge width and gauge length were fixed to 4 µm and 10 µm, respectively. Three-dimensional surface plots of the stress concentration as a function of two parameters were used to identify the optimal ranges of parameter values. The addition of maximum width and length constraints, measuring 25 µm and 30 µm, respectively, allowed us to identify an optimal geometry at load surface angles of 30° and 45°. Their respective dimensions (corner radius, load surface length, and dogbone head length) are, in µm, 12, 6, and 7 and 10, 7, and 7. Testing these two optimal geometries with a range of gauge lengths from 4 to 20 µm showed that smaller gauge lengths only slightly reduced the detrimental stress concentration outside the gauge section. However, smaller gauge lengths will notably improve the FIB surface polishing step as tapering is reduced with smaller dogbone lengths.

Chemistry↗
NASA NTRS1978

The effects of eccentricities on the fracture of off-axis fiber composites

Finite element analyses were performed to investigate theoretically the effects of in-plane and out-of-plane eccentricities, bending or twisting, and thickness nonuniformity on the axial stress and strain variations across the width of off-axis specimens. The results are compared with measured data and are also used to assess the effects of these eccentricities on the fracture stress of off-axis fiber composites. Guidelines for detecting and minimizing the presence of eccentricities are described.

Chamis, C. C.↗
NASA NTRS1978

The effects of eccentricities on the fracture of off-axis fiber composites

Finite element analyses were performed to investigate theoretically the effects of in-plane and out-of plane eccentricities, bending or twisting, and thickness nonuniformity on the axial stress and strain variations across the width of off-axis specimens. The results are compared with measured data and are also used to access the effects of these eccentricities on the fracture stress of off-axis fiber composites. Guidelines for detecting and minimizing the presence of eccentricities are described.

Chamis, C. C.↗
NASA NTRS1985

On the impact of initially stressed composite laminates

The impact response behavior of initially stressed composite laminates is investigated using the finite element method. An experimentally established contact law is incorporated into the finite element program. The Newmark time integration algorithm is used for solving the time dependent equations of the plate and the impactor. Numerical results, including the contact force history, deflection, and strain in the plate, are presented. Effects of impact velocity, initial stress, and the mass and size of the impactor are discussed.

Sun, C. T.↗
NASA NTRS1996

Vibration Control in Turbomachinery Using Active Magnetic Journal Bearings

The effective use of active magnetic bearings for vibration control in turbomachinery depends on an understanding of the forces available from a magnetic bearing actuator. The purpose of this project was to characterize the forces as functions shaft position. Both numerical and experimental studies were done to determine the characteristics of the forces exerted on a stationary shaft by a magnetic bearing actuator. The numerical studies were based on finite element computations and included both linear and nonlinear magnetization functions. Measurements of the force versus position of a nonrotating shaft were made using two separate measurement rigs, one based on strain gage measurement of forces, the other based on deflections of a calibrated beam. The general trends of the measured principal forces agree with the predictions of the theory while the magnitudes of forces are somewhat smaller than those predicted. Other aspects of theory are not confirmed by the measurements. The measured forces in the normal direction are larger than those predicted by theory when the rotor has a normal eccentricity. Over the ranges of position examined, the data indicate an approximately linear relationship between the normal eccentricity of the shaft and the ratio of normal to principal force. The constant of proportionality seems to be larger at lower currents, but for all cases examined its value is between 0.14 and 0.17. The nonlinear theory predicts the existence of normal forces, but has not predicted such a large constant of proportionality for the ratio. The type of coupling illustrated by these measurements would not tend to cause whirl, because the coupling coefficients have the same sign, unlike the case of a fluid film bearing, where the normal stiffness coefficients often have opposite signs. They might, however, tend to cause other self-excited behavior. This possibility must be considered when designing magnetic bearings for flexible rotor applications, such as gas turbines and other turbomachinery.

Knight, Josiah D.↗
NASA NTRS2020

Lessons Learned from Recent Space Flight Assessments

This paper highlights four space flight assessments, each with an important lesson. The first assessment points to the challenge of material substitution for weight saving. Space Shuttle’s external tank stringers were changed from Al 2024 to Al 2090 without comprehensive analysis. Lack of understanding the low Short-Transverse (S-T) strength associated with Al 2090 led to widespread cracking. The second assessment addresses the failure to understand that time-dependent relationships and loading interactions of multiple disciplines can lead to non-conservative results. An aerodynamic test article of the Orion Multi-Purpose Crew Vehicle (MPCV) Launch Abort System (LAS) was analyzed for worst-case single discipline loads, showing positive structural margins. When disciplines were combined with time dependent loading, the analysis resulted in negative margins. The third assessment highlights the pitfalls of putting blind trust in historical architecture. The MPCV re-entry heat shield adopted the Apollo design approach despite fundamental differences between the two designs (e.g., the underlying substructure design, material differences, scale, and geometry). Manufacturing resulted in cracking along seams of the cured MPCV heat shield and produced lower mechanical properties of the Avcoat witness panels. Subsequently, the design was changed to a block architecture, which introduced unperceived technical, cost, and schedule risk. The fourth assessment shows how the lack of a comprehensive investigation of design modifications and inappropriate Finite Element Analyses (FEA) can lead to trouble. A mass reduction campaign eliminated highly critical internal struts on the ISS Node-1 pressurized module based on inadequate and uncorrelated FEA. High strains were noted and ignored as a modeling aberration. A system level high fidelity FEA discovered large plastic strains in the gussets of all four radial ports. These findings initiated a full-scale test campaign, new test correlated assembly-level FEA, and significant design modifications.

James P Smith↗
DOE OSTI2021

Initial Developments in Modeling Graphite Behavior

Graphite is used in High Temperature Reactor designs as a moderator, reflector and core material responsible for protecting the fuel and maintaining structural stability. From a safety perspective, it is imperative to be able to predict how graphite will behave under reactor operating conditions which may compromise the core’s structural integrity or fuel’s safety performance. The following report summarizes recent modeling work undertaken at Idaho National Laboratory which focuses on graphite degradation behavior. One phenomenon which is explored in this report and can cause significant degradation to graphite properties is oxidation. Oxidation behavior in graphite is highly temperature dependent. At low temperatures, oxygen can fully penetrate a component and cause a homogeneous damage profile. At high temperatures, only the graphite near the surface of a component is affected. Understanding and modeling this temperature dependence is essential to predicting graphite behavior. In this report, three graphite models are discussed. The first model is used to investigate strength loss after low temperature oxidation. It does this by generating an approximate graphite microstructure then determining the required load to cause crack propagation. It has been shown that the model reproduces the quasi brittle stress versus strain relationship observed in graphite. Strength loss results from the model are shown to match experimental values. The second model investigates strength loss under high temperature oxidation conditions. This model is more applicable to full scale graphite components where the oxidation damage is often inhomogeneous throughout a component. The third model computes stresses in graphite under reactor conditions. Specifically, the model incorporates the effect of creep, irradiation dose, thermal strains, and oxidation in order to predict internal stresses in a component. All of these models are implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open source, parallel finite element framework.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2009

Further Development of Ko Displacement Theory for Deformed Shape Predictions of Nonuniform Aerospace Structures

The Ko displacement theory previously formulated for deformed shape predictions of nonuniform beam structures is further developed mathematically. The further-developed displacement equations are expressed explicitly in terms of geometrical parameters of the beam and bending strains at equally spaced strain-sensing stations along the multiplexed fiber-optic sensor line installed on the bottom surface of the beam. The bending strain data can then be input into the displacement equations for calculations of local slopes, deflections, and cross-sectional twist angles for generating the overall deformed shapes of the nonuniform beam. The further-developed displacement theory can also be applied to the deformed shape predictions of nonuniform two-point supported beams, nonuniform panels, nonuniform aircraft wings and fuselages, and so forth. The high degree of accuracy of the further-developed displacement theory for nonuniform beams is validated by finite-element analysis of various nonuniform beam structures. Such structures include tapered tubular beams, depth-tapered unswept and swept wing boxes, width-tapered wing boxes, and double-tapered wing boxes, all under combined bending and torsional loads. The Ko displacement theory, combined with the fiber-optic strain-sensing system, provide a powerful tool for in-flight deformed shape monitoring of unmanned aerospace vehicles by ground-based pilots to maintain safe flights.

Ko, William L.↗
NASA NTRS1985

Finite Element Prediction of Loss Factors for Structures with Frequency-dependent Damping Treatments

A finite element procedure is described for calculating the loss factors for elastic structures to which frequency-dependent viscoelastic damping treatments were applied. The frequency dependence of the viscoelastic damping material is treated by approximating its shear modulus with a second-order polynomial so that the stiffnesses associated with the constant, linear, and quadratic terms can be combined, respectively, with the stiffness, damping, and mass matrices assembled for the rest of the structure. A single complex eigenvalue analysis is then performed in which the eigenvalues are purely imaginary. The loss factor is computed by the modal strain energy (MSE) approach. In the the MSE approach, the loss factor of a composite structure vibrating in one of its natural modes may be visualized as a weighted average of the loss factors of the component parts, with the relative stored energies as weighting constants. The finite element procedure, which can treat very general geometries, is illustrated for the case of a vibrating constrained-layer damped plate.

Everstine, G. C.↗
NASA NTRS2000

Vascular mechanics of the coronary artery

This paper describes our research into the vascular mechanics of the coronary artery and plaque. The three sections describe the determination of arterial mechanical properties using intravascular ultrasound (IVUS), a constitutive relation for the arterial wall, and finite element method (FEM) models of the arterial wall and atheroma. METHODS: Inflation testing of porcine left anterior descending coronary arteries was conducted. The changes in the vessel geometry were monitored using IVUS, and intracoronary pressure was recorded using a pressure transducer. The creep and quasistatic stress/strain responses were determined. A Standard Linear Solid (SLS) was modified to reproduce the non-linear elastic behavior of the arterial wall. This Standard Non-linear Solid (SNS) was implemented into an axisymetric thick-walled cylinder numerical model. Finite element analysis models were created for five age groups and four levels of stenosis using the Pathobiological Determinants of Atherosclerosis Youth (PDAY) database. RESULTS: The arteries exhibited non-linear elastic behavior. The total tissue creep strain was epsilon creep = 0.082 +/- 0.018 mm/mm. The numerical model could reproduce both the non-linearity of the porcine data and time dependent behavior of the arterial wall found in the literature with a correlation coefficient of 0.985. Increasing age had a strong positive correlation with the shoulder stress level, (r = 0.95). The 30% stenosis had the highest shoulder stress due to the combination of a fully formed lipid pool and a thin cap. CONCLUSIONS: Studying the solid mechanics of the arterial wall and the atheroma provide important insights into the mechanisms involved in plaque rupture.

NASA Program Biomedical Research and Countermeasur↗
DOE OSTI2020

On the mechanical behavior of a Nb 3 Sn superconducting coil during a quench: Two-dimensional finite element analysis of a quench heater protected magnet

New high-field accelerator magnets based on $Nb_3Sn$ superconductors are pushing the boundaries of magnet design and quench protection towards new limits. While their large stored energies and current densities result in a very challenging scenario for magnet protection, their great electromagnetic forces create also new requirements in terms of magnet design and stress management techniques. Furthermore, the strain sensitivity of $Nb_3Sn$ cables turns the electro-mechanical limits of the conductor into a parameter of the highest importance, where conductor degradation becomes a critical aspect in magnet operation. The coupling of all the above-mentioned considerations during quench is a case of special interest that adds further complexity to the design of $Nb_3Sn$ magnets. The objective of this paper is to provide a complete two-dimensional investigation of the coil and magnet structure mechanics during a quench event. Here, the analysis is performed using a combination of finite element codes that provide the necessary input for the mechanical study. The core of the modelling strategy relies on a thermal-electric model, whose results are directly used as loads in the mechanical simulation. In doing so, the stress evolution during and after a quench is obtained. We focus for this time in the analysis of a quench heater protected magnet, where dynamic effects are of less importance in contrast to other protection systems like CLIQ.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Fundamental Mechanisms Controlling Dislocation-Obstacle Interactions in Metals and Alloys

This project combined modeling and experimentation across micro-meso-macro length scales to investigate the role of long-range internal stresses associated with obstacles and dislocation populations, in the deformation response of cubic metals. In addition to unique insights into the physical behavior of metals, advances were also made in modeling techniques and tools at the meso and atomistic level, and new characterization methods were developed and disseminated to the scientific community. Twenty eight journal papers, two dissertations and three theses have served as vehicles to share the gained scientific knowledge, and two opensource software packages have been developed, maintained and distributed. Additional publications are in progress.

36 MATERIALS SCIENCE↗
NASA NTRS1984

A comparison of several methods of calculating fracture toughness from composite laminate test data

Critical fracture toughness was determined by two different techniques for graphite/epoxy laminates of three stacking sequences and several thicknesses. Critical fracture toughness was determined by a finite element stress analysis of a center-cracked tension specimen which yielded the stress intensity factor as a function of specimen thickness and applied load. As an alternative approach the critical strain energy release rate was determined from the compliance calibration technique. Test results from both procedures are compared and discussed.

Harris, C. E.↗
NASA NTRS1985

Axisymmetric solid elements by a rational hybrid stress method

Four-node axisymmetric solid elements are derived by a new version of hybrid method for which the assumed stresses are expressed in complete polynomials in natural coordinates. The stress equilibrium conditions are introduced through the use of additional displacements as Lagrange multipliers. A rational procedure is to choose the displacement terms such that the resulting strains are also of complete polynomials of the same order. Example problems all indicate that elements obtained by this procedure lead to better results in displacements and stresses than that by other finite elements.

Tian, Z.↗