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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Buckling Analysis of Debonded Sandwich Panel Under Compression

A sandwich panel with initial through-the-width debonds is analyzed to study the buckling of its faceskin when subject to an in-plane compressive load. The debonded faceskin is modeled as a beam on a Winkler elastic foundation in which the springs of the elastic foundation represent the sandwich foam. The Rayleigh-Ritz and finite-difference methods are used to predict the critical buckling load for various debond lengths and stiffnesses of the sandwich foam. The accuracy of the methods is assessed with a plane-strain finite-element analysis. Results indicate that the elastic foundation approach underpredicts buckling loads for sandwich panels with isotropic foam cores.

Sleight, David W.↗
NASA NTRS1990

Weight minimization of a support structure

This paper addresses the weight minimization of a circular plate-like structure which resulted in a 26 percent weight reduction. The optimization was performed numerically with the COPES/ADS program using the modified method of feasible directions. Design parameters were the inner thickness and outer thickness of the plate with constraints on maximum yield stress and maximum transverse displacement. Also, constraints were specified for the upper and lower bounds of the fundamental frequency and plate thicknesses. The MSC/NASTRAN finite element program was used for the evaluation of response variables. Original and final designs of the plate were tested using an Instron tension-compression machine to compare finite element results to measured strain data. The difference between finite element strain components and measured strain data was within engineering accuracy.

Kluberdanz, Donald J.↗
DOE OSTI2020

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.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Measurements of elastoresistance under pressure by combining in-situ tunable quasi-uniaxial stress with hydrostatic pressure

Uniaxial stress, as well as hydrostatic pressure are often used to tune material properties in condensed matter physics. Here, we present a setup which allows for the study of the combined effects of quasi-uniaxial stress and hydrostatic pressure. Following earlier designs for measurements under finite stress at ambient pressure (e.g., Chu et al., Science 337, 710 (2012)), the present setup utilizes a piezoelectric actuator to change stress in situ inside the piston-cylinder pressure cell. We show that the actuator can be operated over the full temperature (from 30 K up to 260 K) and pressure range (up to ≈ 2 GPa), resulting in a clear and measurable quasi-uniaxial strain. To demonstrate functionality, measurements of the elastoresistance (i.e., the change of resistance of a sample as a response to quasi-uniaxial strain) under finite hydrostatic pressure on the iron-based compound BaFe 2 As 2 are presented as a proof-of-principle example, and discussed in the framework of electronic nematicity. Overall, this work introduces the combination of in situ tunable quasi-uniaxial stress and large (up to ≈ 2 GPa) hydrostatic pressure as a powerful combination in the study of novel electronic phases. In addition, it also points towards further technical advancements which can be made in the future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Anomalous softening of 3D printed elastomeric foam irradiated under compressive strain

Elastomeric foam is an essential component in many industrial and technological settings, primarily as thermal insulators and as positional/mechanical support cushions. In particular, silicone foam is utilized in harsh environments due to exceptional thermal and chemical stability. Under service conditions within certain applications such material gets exposed to a high dosage of gamma radiation, which can permanently alter the material’s structural and mechanical response properties. Most studies on gamma-exposure under inert or oxidative atmosphere indicate hardening of silicone foam, which is attributed to radiation-induced enhancement in chemical cross-linking. Here we report two contrasting effects depending on whether (non-oxidative) radiation exposure is carried out with the foam under zero or finite compressive strain. While in the former case we observe radiation-hardening consistent with previous studies, in the latter case (50% porous foam under 30% uniaxial compression) we see a monotonic decrease in Young’s modulus with increasing dosage, although solvent swelling experiments on the constituent rubber indicate a net increase in cross-link density independent of the state of strain. We quantitatively model all dose-dependent data using the Ogden Hyperfoam strain-energy function within the framework of Tobolsky two-network scheme and attribute the above anomaly to a combined effect of radiation-induced thickness change (compression set) and inherent nonlinearity in the foam’s stress-strain response.

Coarse-grained models↗
DOE OSTI2023

Implementation and experimental validation of nonlocal damage in a large-strain elasto-viscoplastic FFT-based framework for predicting ductile fracture in 3D polycrystalline materials

Ductile materials, such as metal alloys, can undergo substantial deformation before failure. Additionally, these materials are usually of polycrystalline composition and exhibit strongly anisotropic behavior at small length scales. Previously developed fast Fourier transform (FFT)-based models can model ductile fracture of isotropic materials or the elastic–plastic behavior of anisotropic polycrystalline materials; however, there remains a need to couple both capabilities. This work extends a large-strain FFT-based crystal plasticity model to simulate ductile fracture of polycrystalline materials. In this work, a triaxiality-based continuum damage mechanics (CDM) formulation is incorporated into a large-strain elasto-viscoplastic FFT (LS-EVPFFT) framework. The CDM formulation is augmented with an integral-based nonlocal regularization approach that correctly handles gas-phase material necessary to model unconstrained surfaces. To validate the damage-enabled LS-EVPFFT framework, mesoscale copper tensile coupons were machined using microwire electrical discharge machining and experimentally characterized using electron backscatter diffraction. In-situ optical digital image correlation was performed during uniaxial testing to provide a side-by-side comparison of the experimental and computational strain fields and stress–strain responses. The damage-enabled LS-EVPFFT framework can simulate the complete macroscopic stress–strain response of ductile polycrystals to failure. The model reproduces necking behavior that qualitatively agrees with experimental observations. By leveraging the relatively low computational cost of the damage-enabled LS-EVPFFT framework, the framework presented here allows the ductile fracture response of 3D polycrystalline materials to be tractably predicted.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Gravitational spreading of high terrain in Ishtar Terra, Venus

A plane-strain viscoelastic finite element algorithm derived by Melosh and Rafesky (1980, 1983) is used to model the gravitational spreading of the high topography of Ishtar Terra. Brittle failure and relaxation of topographic relief are predicted as functions of time for ranges of crustal thickness and thermal gradient and for observed ranges in topographic elevation and slope. For a crustal thickness greater than 10 km, a thermal gradient of 15 K/km or more, and average values of relief (3 km) and slope (3 deg), the topography relaxes to 25 percent of its original height in less than 10 m.y. Although initial failure, in the form of normal faulting on the highlands and margin and occasional shallow thrusting in the lowlands, is predicted to occur much earlier, values of horizontal surface strain that are likely to be observable (about 1 percent) do not accumulate until significant relaxation of relief begins.

Smrekar, Suzanne E.↗
DOE OSTI2023

Computational design of metamaterials with self contact

Inverse homogenization in combination with contact modeling, topology optimization and shape optimization is used to design metamaterials with optimized macroscopic response. The homogenization assumes length scale separation which allows the non-linear macroscopic behavior to be obtained by analyzing a single unit cell in a lattice structure. Self contact in the unit cell, which is modeled using a third medium contact method, is leveraged to obtain a complex homogenized response. The inverse homogenization problem is initially formulated as a topology optimization problem, where the macroscopic stress–strain behavior is tuned to our liking. However, it is well known that boundary phenomena are difficult to model in topology optimization and that interface modeling is crucial to accurately analyze contact. For that reason, the boundary representation of the topology optimized design is extracted and used as initial design in a subsequent shape optimization. The behaviors of our designs are verified by performing rigorous post-processing analyzes using conforming meshes and conventional contact formulations.

42 ENGINEERING↗
DOE OSTI2022

Tunable phononic bandgap materials designed via topology optimization

Topology optimization is used to design phononic bandgap materials that are tunable by mechanical deformation. A periodic media is considered, which due to the assumption of length scale separation, allows the dispersion relations to be obtained by analyzing a single unit cell subjected to Floquet–Bloch boundary conditions. A finite macroscopic deformation is applied to the unit cell to affect its geometry and hence dispersion. We tune the dispersion–deformation relation to our liking by solving a topology optimization problem using nonlinear programming. The adjoint method is employed to compute the sensitivities, and the non-differentiability of degenerate eigenvalues is avoided using symmetric polynomials. Several tunable phononic crystal designs are presented. Also, a verification analysis is performed, wherein the optimized design is interpreted and analyzed using a conforming finite element mesh.

42 ENGINEERING↗
DOE OSTI2021

Thermomechanical Microstructural Predictions of Fracture Nucleation of Zircaloy-4 Alloys With δ and ε Hydride Distributions

Abstract A crystalline dislocation-density formulation that was incorporated with a nonlinear finite-element (FE) method was utilized to understand and to predict the thermomechanical behavior of an hexagonal closest packed (h.c.p.) zircaloy system with hydrides with either face-centered cubic (f.c.c.) or body-centered cubic (b.c.c.) hydrides. This formulation was then used with a recently developed fracture methodology that is adapted for finite inelastic strains and multiphase crystalline systems to understand how different microstructurally based fracture modes nucleate and propagate. The interrelated microstructural characteristics of the different crystalline hydride and matrix phases with the necessary orientation relationships (ORs) have been represented, such that a detailed physical understanding of fracture nucleation and propagation can be predicted for the simultaneous thermomechanical failure modes of hydride populations and the matrix. The effects of volume fraction, morphology, crystalline structure, and orientation and distribution of the hydrides on simultaneous and multiple fracture modes were investigated for radial, circumferential, and mixed distributions. Another key aspect was accounting for temperatures changes due to the effects of thermal conduction and dissipated plastic work and their collective effects on fracture. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and growth. The predictions provide validated insights into why circumferential hydrides are more fracture-resistant than radial hydrides for different volume fractions and thermomechanical loading conditions.

Engineering↗
NASA NTRS1975

Application of NASTRAN for stress analysis of left ventricle of the heart

Knowing the stress and strain distributions in the left ventricular wall of the heart is a prerequisite for the determination of the muscle elasticity and contractility in the process of assessing the functional status of the heart. NASTRAN was applied for the calculation of these stresses and strains and to help in verifying the results obtained by the computer program FEAMPS which was specifically designed for the plane-strain finite-element analysis of the left ventricular cross sections. Adopted for the analysis are the true shape and dimensions of the cross sections reconstructed from multiplanar X-ray views of a left ventricle which was surgically isolated from a dog's heart but metabolically supported to sustain its beating. A preprocessor was prepared to accommodate both FEAMPS and NASTRAN, and it has also facilitated the application of both the triangular element and isoparameteric quadrilateral element versions of NASTRAN. The stresses in several crucial regions of the left ventricular wall calculated by these two independently developed computer programs are found to be in good agreement. Such confirmation of the results is essential in the development of a method which assesses the heart performance.

Pao, Y. C.↗
NASA NTRS1987

Impact-induced fracture in a quasi-isotropic laminate

A systematic study of impact-induced fracture in a quasi-isotropic laminated composite is carried out. The main focus of the study is to understand damage initiation when a laminate is subjected to the impact of a foreign object. The total incipient damage of a laminate subjected to impact at higher than threshold velocity is also presented. The incipient damage is restricted to small growth from the initiation of the damage. The experimentally collected data are interpreted using a two-dimensional plane-strain finite-element analysis. The qualitative comparison suggests that the skew cracks in the proximal layers are primarily due to transverse shear stress. Skew cracks in the middle layers are also due to transverse shear stress. Vertical cracks in the distal layer are due to flexural stress transverse to the fiber orientation of the layer.

Joshi, S. P.↗
NASA NTRS1985

Finite element analysis of thermal barrier coatings

The near-term objective is to develop an understanding of the states of stresses and strains in a Zirconia-yttria thermal barrier coating (TBC) experiencing a given temperature drop. Results so obtained are expected to facilitate experimental work. In order to gain realistic insights into the distribution of stresses and strains in a complex TBC, the finite element approach was selected to model a cylindrical TBC specimen. Experimental evidence reported in the literature indicated the presence of rough interface between the ceramic coat and bond coat. Oxidation of the bond coat at ceramic-bond interface was observed, as was a small amount of cracking in the ceramics near the ceramic-bond interface. To account for these complex features, a plane-strain finite element computer program known as TBCOC was developed, taking advantage of a generic computer code known as MARC. This generic code was made available through the use of a supercomputer (Cray 1). The TBCOC model contains 1316 nodal points and 2140 finite elements. It is capable of a uniform isothermal loading. Results of a sample computer run are presented.

Chang, G. C.↗
NASA NTRS1984

Silicon ribbon stress-strain activities

The finite element method is used to investigate stress/strain in silicon ribbon. Failure considerations such as residual stress, buckling material non-linearity and creep are discussed. Temperature profiles are presented.

Wada, B. K.↗
NASA NTRS1987

Calculation of strain-energy release rates with higher order and singular finite elements

A general finite element procedure for obtaining strain-energy release rates for crack growth in isotropic materials is presented. The procedure is applicable to two-dimensional finite element analyses and uses the virtual crack-closure method. The procedure was applied to nonsingular 4-noded (linear), 8-noded (parabolic), and 12-noded (cubic) elements and to quarter-point and cubic singularity elements. Simple formulas for strain-energy release rates were obtained with this procedure for both nonsingular and singularity elements. The formulas were evaluated by applying them to two mode I and two mixed mode problems. Comparisons with results from the literature for these problems showed that the formulas give accurate strain-energy release rates.

Raju, I. S.↗
NASA NTRS1986

Simple formulas for strain-energy release rates with higher order and singular finite elements

A general finite element procedure for obtaining strain-energy release rates for crack growth in isotropic materials is presented. The procedure is applicable to two-dimensional finite element analyses and uses the virtual crack-closure method. The procedure was applied to non-singular 4-noded (linear), 8-noded (parabolic), and 12-noded (cubic) elements and to quarter-point and cubic singularity elements. Simple formulas for strain-energy release rates were obtained with this procedure for both non-singular and singularity elements. The formulas were evaluated by applying them to two mode I and two mixed mode problems. Comparisons with results from the literature for these problems showed that the formulas give accurate strain-energy release rates.

Raju, I. S.↗
NASA NTRS1980

Computational strategy for the solution of large strain nonlinear problems using the Wilkins explicit finite-difference approach

The STEALTH code system, which solves large strain, nonlinear continuum mechanics problems, was rigorously structured in both overall design and programming standards. The design is based on the theoretical elements of analysis while the programming standards attempt to establish a parallelism between physical theory, programming structure, and documentation. These features have made it easy to maintain, modify, and transport the codes. It has also guaranteed users a high level of quality control and quality assurance.

Hofmann, R.↗
NASA NTRS1982

Creep crack-growth: A new path-independent integral (T sub c), and computational studies

The development of valid creep fracture criteria is considered. Two path-independent integral parameters which show some degree of promise are the C* and (Delta T)sub c integrals. The mathematical aspects of these parameters are reviewed by deriving generalized vector forms of the parameters using conservation laws which are valid for arbitrary, three dimensional, cracked bodies with crack surface tractions (or applied displacements), body forces, inertial effects, and large deformations. Two principal conclusions are that (Delta T)sub c has an energy rate interpretation whereas C* does not. The development and application of fracture criteria often involves the solution of boundary/initial value problems associated with deformation and stresses. The finite element method is used for this purpose. An efficient, small displacement, infinitesimal strain, displacement based finite element model is specialized to two dimensional plane stress and plane strain and to power law creep constitutive relations. A mesh shifting/remeshing procedure is used for simulating crack growth. The model is implemented with the quartz-point node technique and also with specially developed, conforming, crack-tip singularity elements which provide for the r to the n-(1+n) power strain singularity associated with the HRR crack-tip field. Comparisons are made with a variety of analytical solutions and alternate numerical solutions for a number of problems.

Stonesifer, R. B.↗