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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Characterization of the Edge Crack Torsion (ECT) Test for Mode III Fracture Toughness Measurement of Laminated Composites

The edge crack torsion (ECT) test is designed to initiate mode III delamination growth in composite laminates. The test has undergone several design changes during its development. The objective of this paper was to determine the suitability of the current ECT test design a mode III fracture test. To this end, ECT tests were conducted on specimens manufactured from IM7/8552 and S2/8552 tape laminates. Three-dimensional finite element analyses were performed. The analysis results were used to calculate the distribution of mode I, mode II, and mode III strain energy release rate along the delamination front. The results indicated that mode IIIdominated delamination growth would be initiated from the specimen center. However, in specimens of both material types, the measured values of GIIIc exhibited significant dependence on delamination length. Load-displacement response of the specimens exhibited significant deviation from linearity before specimen failure. X-radiographs of a sample of specimens revealed that damage was initiated in the specimens prior to failure. Further inspection of the failure surfaces is required to identify the damage and determine that mode III delamination is initiated in the specimens.

Ratcliffe, James G.↗
NASA NTRS2008

Software Performs Complex Design Analysis

Designers use computational fluid dynamics (CFD) to gain greater understanding of the fluid flow phenomena involved in components being designed. They also use finite element analysis (FEA) as a tool to help gain greater understanding of the structural response of components to loads, stresses and strains, and the prediction of failure modes. Automated CFD and FEA engineering design has centered on shape optimization, which has been hindered by two major problems: 1) inadequate shape parameterization algorithms, and 2) inadequate algorithms for CFD and FEA grid modification. Working with software engineers at Stennis Space Center, a NASA commercial partner, Optimal Solutions Software LLC, was able to utilize its revolutionary, one-of-a-kind arbitrary shape deformation (ASD) capability-a major advancement in solving these two aforementioned problems-to optimize the shapes of complex pipe components that transport highly sensitive fluids. The ASD technology solves the problem of inadequate shape parameterization algorithms by allowing the CFD designers to freely create their own shape parameters, therefore eliminating the restriction of only being able to use the computer-aided design (CAD) parameters. The problem of inadequate algorithms for CFD grid modification is solved by the fact that the new software performs a smooth volumetric deformation. This eliminates the extremely costly process of having to remesh the grid for every shape change desired. The program can perform a design change in a markedly reduced amount of time, a process that would traditionally involve the designer returning to the CAD model to reshape and then remesh the shapes, something that has been known to take hours, days-even weeks or months-depending upon the size of the model.

Source record↗
NASA NTRS2013

Fracture Mechanics Analyses for Interface Crack Problems - A Review

Recent developments in fracture mechanics analyses of the interfacial crack problem are reviewed. The intent of the review is to renew the awareness of the oscillatory singularity at the crack tip of a bimaterial interface and the problems that occur when calculating mode mixity using numerical methods such as the finite element method in conjunction with the virtual crack closure technique. Established approaches to overcome the nonconvergence issue of the individual mode strain energy release rates are reviewed. In the recent literature many attempts to overcome the nonconvergence issue have been developed. Among the many approaches found only a few methods hold the promise of providing practical solutions. These are the resin interlayer method, the method that chooses the crack tip element size greater than the oscillation zone, the crack tip element method that is based on plate theory and the crack surface displacement extrapolation method. Each of the methods is validated on a very limited set of simple interface crack problems. However, their utility for a wide range of interfacial crack problems is yet to be established.

Krueger, Ronald↗
NASA NTRS2019

Through-Thickness Tensile Properties of Triaxial Braided Composites

The through-thickness strength of triaxial braided composites was measured using thin specimens bonded to metallic tabs. Despite providing reasonable engineering estimates of the bulk properties, some test specimen limitations were noted: For strong epoxies, significant shear and tensile stresses occur and induce failure between the first layer and bond-line; the strains over much of the exterior are ~10% lower than interior strains, which are representative of the nominal state; and failure often occurred near or along the adhesive bond line. These interference were partially overcome by adding a large radius to the gage section. Based on the finite element analysis and experimental measurements, the nominal elastic modulus and Poisson's ratio are E(sub 33) = 10.7 GPa and v(sub 31) = 0.08 for the three triaxial braided composites tested. Strengths vary from 14 to 40 MPa depending on the resin. Relative to the in-plane properties, the through thickness modulus, fracture strength, failure strain, and Poisson's ratio are 1/5th, 1/30th, 1/6th, and 1/4th. Interlaminar failure occurred with varying degrees of fiber tear out along with little damage to axial tows, resulting is a smooth undulating fracture surface. Post cure temperature excursions to 177 deg C (350 deg F) had little effect on elastic constants but reduced the strength E862 based composites severely. The strength of PR520 and 5208 based composites increased.

composite↗
DOE OSTI2022

Modeling and experimental validation of dynamical effects in Bragg coherent x-ray diffractive imaging of finite crystals

Bragg coherent diffractive imaging (BCDI) is a noninvasive microscopy technique that can visualize the shape and internal lattice deviations of crystals with nanoscale spatial resolution and picometer deformation sensitivity. Its strain imaging capability relies on Fourier transform–based iterative phase retrieval algorithms, which are mostly developed under the kinematical approximation. Such approximation prohibits the application of BCDI on larger crystals, which are commonly seen in most emerging functional materials. Understanding the dynamical effect in BCDI, as well as developing a validated method for modeling BCDI at the dynamical diffraction limit, is crucial for applying BCDI to hierarchical systems that contain micron-sized crystals and grains. Thus we report a comparative study on the impact of dynamical diffraction effects by comparing the reconstruction results from two measurements of the same crystal. Forward simulation is implemented to show subtle changes of interference fringes in the diffraction pattern due to the dynamical diffraction, and is compared directly with the experimental data.

36 MATERIALS SCIENCE↗
NASA NTRS1992

An improved plate theory of order (1,2) for thick composite laminates

A new (1,2)-order theory is proposed for the linear elasto-static analysis of laminated composite plates. The basic assumptions are those concerning the distribution through the laminate thickness of the displacements, transverse shear strains and the transverse normal stress, with these quantities regarded as some weighted averages of their exact elasticity theory representations. The displacement expansions are linear for the inplane components and quadratic for the transverse component, whereas the transverse shear strains and transverse normal stress are respectively quadratic and cubic through the thickness. The main distinguishing feature of the theory is that all strain and stress components are expressed in terms of the assumed displacements prior to the application of a variational principle. This is accomplished by an a priori least-square compatibility requirement for the transverse strains and by requiring exact stress boundary conditions at the top and bottom plate surfaces. Equations of equilibrium and associated Poisson boundary conditions are derived from the virtual work principle. It is shown that the theory is particularly suited for finite element discretization as it requires simple C(sup 0)- and C(sup -1)-continuous displacement interpolation fields. Analytic solutions for the problem of cylindrical bending are derived and compared with the exact elasticity solutions and those of our earlier (1,2)-order theory based on the assumed displacements and transverse strains.

Tessler, A.↗
NASA NTRS1991

Improved compression buckling for rectangular composite plates by stiffness tailoring

Buckling resistance is often a controlling criterion in the design of structural elements. Design concepts that lead to increased buckling loads (or strains) can directly lower the structural cost and/or weight by a number of means. This study quantifies the improvements that can be achieved in compression buckling loads of rectangular composite plates by using a simple stiffness tailoring concept. The approach is to position the unidirectional lamina through the thickness and over the planform of the plate so that the buckling load is increased with no loss in in-plane stiffness or increase in weight. Finite element analyses have been used to determine the effects of tailoring on the buckling load of plates with various boundary conditions, aspect ratios, thicknesses, and membrane stiffnesses. Increases in buckling loads (or strains) of nearly 200 percent over the uniform plate buckling loads are shown possible with this tailoring concept.

Biggers, Sherrill B.↗
DOE OSTI2020

In situ measurement and modelling of the growth and length scale of twins in α -uranium

Twinning is an important deformation mechanism in crystalline materials. Twins nucleate, then grow, forming a discrete pattern of deformation bands with a certain thickness. The mechanism for twin growth is not completely understood. We present a simple phase field model which captures the physics of twin growth and reproduces the details observed during in situ electron backscatter diffraction experiments on α-uranium. The interaction between residual dislocations at twin interfaces and mobile dislocations in untwinned regions increases the stress needed for twin growth. This is described by a nonlocal term in the proposed constitutive equations: the nucleation stress of a twin increases proportionally to the twin phase field in a neighborhood. Competition between slip and twinning favors the nucleation of a new twin rather than the growth of a preexisting twin. The phase field model is coupled with a crystal plasticity finite element solver that includes the plastic deformation due to twinning. Simulations are able to reproduce the number of twins and their thicknesses as a function of the strain observed during in situ electron backscatter diffraction experiments. The stress-strain curve is also reproduced. This comparison allows the values of the nucleation stress and the interaction strength between residual and mobile dislocations to be found. It also gives an estimation of the density of residual dislocations that is consistent with the observed twin thickness. Finally, this model can be applied to understand microstructural effects in materials using twinning as a strengthening mechanism.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Global functions in global-local finite-element analysis of localized stresses in prismatic structures

An important consideration in the global local finite-element method (GLFEM) is the availability of global functions for the given problem. The role and mathematical requirements of these global functions in a GLFEM analysis of localized stress states in prismatic structures are discussed. A method is described for determining these global functions. Underlying this method are theorems due to Toupin and Knowles on strain energy decay rates, which are related to a quantitative expression of Saint-Venant's principle. It is mentioned that a mathematically complete set of global functions can be generated, so that any arbitrary interface condition between the finite element and global subregions can be represented. Convergence to the true behavior can be achieved with increasing global functions and finite-element degrees of freedom. Specific attention is devoted to mathematically two-dimensional and three-dimensional prismatic structures. Comments are offered on the GLFEM analysis of NASA flat panel with a discontinuous stiffener. Methods for determining global functions for other effects are also indicated, such as steady-state dynamics and bodies under initial stress.

Dong, Stanley B.↗
DOE OSTI2020

Constitutive modeling of δ-phase zircaloy hydride based on strain rate dependent nanoindentation and nano-scale impact dataset

In this work, a constitutive model for zircaloy-4 embedded with δ-phase hydride rim/blister structures (δ-phase zircaloy hydride) is presented. Nanoindentation and nano-scale impact experiments were performed on the δ-hydride blisters as a function of temperature in order to obtain strain rate dependent location specific constitutive response. A power law viscoplastic constitutive model was chosen to represent the strain rate dependent behavior. The strength coefficient and the strain hardening exponent in the power law relation were found to decrease with increase in temperature indicating softening behavior at higher temperatures. The strain rate exponent remained constant with increase in temperature indicating that the rate dependent ductility does not change significantly at higher temperatures. A finite element method based model was created to analyze the failure behavior and the ductile to brittle transition in zircaloy hydride. The model predicted that fracture strength of δ-phase zircaloy hydride decreases linearly with increase in temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reconfigurable quantum phononic circuits via piezo-acoustomechanical interactions

Abstract We show that piezoelectric strain actuation of acoustomechanical interactions can produce large phase velocity changes in an existing quantum phononic platform: aluminum nitride on suspended silicon. Using finite element analysis, we demonstrate a piezo-acoustomechanical phase shifter waveguide capable of producing ± π phase shifts for GHz frequency phonons in 10s of μm with 10s of volts applied. Then, using the phase shifter as a building block, we demonstrate several phononic integrated circuit elements useful for quantum information processing. In particular, we show how to construct programmable multi-mode interferometers for linear phononic processing and a dynamically reconfigurable phononic memory that can switch between an ultra-long-lifetime state and a state strongly coupled to its bus waveguide. From the master equation for the full open quantum system of the reconfigurable phononic memory, we show that it is possible to perform read and write operations with over 90% quantum state transfer fidelity for an exponentially decaying pulse.

Taylor, Jeffrey C. (ORCID:0000000215982436)↗
DOE OSTI2020

Thermal / structural analysis of the HB 650 thermal shield

Fermilab’s PIP-II project’s superconducting linear accelerator will drive the next generation of particle accelerators through a revolution in beam intensity. Key to beam intensification are the high efficiencies of niobium-tin superconducting radiofrequency (SRF) cavities operating at cryogenic temperatures near 5 Kelvin. A multifaceted approach is employed to achieve and maintain the extreme temperature. Vacuum provides the first barrier to thermal convection. Physical thermal intercept zones further isolate the exterior shell vacuum vessel at 300 Kelvin from the supercooled beamline. The first thermal intercept lies just inside the exterior vacuum vessel forming the 40 Kelvin zone. This 40K thermal shield is the focus of this investigation. In preparation for operation, the thermal shield is cooled from 300K to 40K in a slow process over 2-1/2 days. The factors limiting the cooling rate are a high thermal gradient which produces thermal strain and high mechanical stress. D ecreasing the time required to reach operational temperatures is the desired research goal. Physical testing data and finite element method (FEM) computer models from the Single Spoke Resonator – 1 (SSR1) and Linac Coherent Light Source (LCLS-II) cryomodules will be evaluated to determine needed design changes. The thermal shield design for the High Beta 650 MHz (HB650) will be modified and analyzed using FEM computer models to determine maximum cooling rate with a target of 50 thermal cycles of life. Maximum allowable stress will be determined through low cycle fatigue calculations.

43 PARTICLE ACCELERATORS↗
NASA NTRS1974

Small-scale yielding analysis of mixed mode plane-strain crack problems

The small-scale yielding analysis of an elastic-plastic body with a line crack under plane-strain conditions subject to combinations of Mode I and II loadings is examined. The analysis of the near-tip field follows the works of Hutchinson and Rice and Rosengren. Dominant singularity solutions governing the asymptotic behavior of the stresses and strains at the crack tip are obtained for the complete range of loadings between Mode I and II. The results of an accurate finite element technique, which imbeds the dominant singularity solutions, directly relates the near-tip behavior to the elastic stress intensity factors K(I) and K(II). Implications of this study to mixed mode fracture mechanics is also discussed, particularly with respect to the direction of crack initiation and the relation of fracture toughness under mixed modes to that in Mode I. Details of the mixed mode plastic zone sizes and shapes are also given.

Shih, C. F.↗
NASA NTRS1980

Three dimensional inelastic finite element analysis of laminated composites

Formulations of the inelastic response of laminated composites to thermal and mechanical loading are used as the basis for development of the computer NALCOM (Nonlinear Analysis of Laminated Composites) computer program which uses a fully three dimensional isoparametric finite element with 24 nodes and 72 degrees of freedom. An incremental solution is performed with nonlinearities introduced as pseudoloads computed for initial strains. Equilibrium iteration may be performed at every step. Elastic and elastic-plastic response of boron/epoxy and graphite/epoxy graphite/epoxy and problems of curing 0/90 sub s Gr/Ep laminates with and without circular holes are analyzed. Mechanical loading of + or - 45sub s Gr/Ep laminates is modeled and symmetry conditions which exist in angle-ply laminates are discussed. Results are compared to experiments and other analytical models when possible. All models are seen to agree reasonably well with experimetnal results for off-axis tensile coupons. The laminate analyses show the three dimensional effects which are present near holes and free corners.

Griffin, O. H., Jr.↗
NASA NTRS1984

Biaxial constitutive equation development for single crystals

Current gas turbine engines utilize large single crystal superalloy components in the hot section. Structural analysis of these components requires a valid stress strain temperature constitutive equation. The goal of the program described is to create one or more models and verify these models. A constitutive equation based on an assumed slip behavior of a single slip system was formulated, programmed, and debugged. Specifically, the basic theory for a model based on aggravating slip behavior on individual slip systems was formulated and programmed and some simulations were run using assumed values of constants. In addition, a formulation allowing strain controlled simulations was completed. An approach to structural analysis of the specimen was developed. This approach uses long tube consistancy conditions and finite elements specially formulated to take advantage of the symmetry of 100 oriented specimens.

Jordan, E. H.↗
NASA NTRS1985

Use of statical indentation laws in the impact analysis of laminated composite plates

The low-velocity impact response of graphite/epoxy laminates was investigated theoretically and experimentally. A nine-node isoparametric plate finite element in conjunction with an empirical contact law was used for the theoretical investigation. Theoretical results are in good agreement with strain-gage experimental data. The results of the investigation indicate that the present theoretical procedure describes the impact response of laminate for low-impact velocities.

Tan, T. M.↗
NASA NTRS1986

A refined shear deformation theory for the analysis of laminated plates

A refined, third-order plate theory that accounts for the transverse shear strains is presented, the Navier solutions are derived for certain simply supported cross-ply and antisymmetric angle-ply laminates, and finite-element models are developed for general laminates. The new theory does not require the shear correction factors of the first-order theory (i.e., the Reissner-Mindlin plate theory) because the transverse shear stresses are represented parabolically in the present theory. A mixed finite-element model that uses independent approximations of the generalized displacements and generalized moments, and a displacement model that uses only the generalized displacements as degrees of freedom are developed. The displacement model requires C sup 1-continuity of the transverse deflection across the inter-element boundaries, whereas the mixed model requires a C sup 0-element. Also, the mixed model does not require continuous approximations (between elements) of the bending moments. Numerical results are presented to show the accuracy of the present theory in predicting the transverse stresses. Numerical results are also presented for the nonlinear bending of plates, and the results compare well with the experimental results available in the literature.

Reddy, J. N.↗
NASA NTRS1988

Finite Element Composite Analysis Program (FECAP) for a microcomputer

A special purpose finite element composite analysis program for analyzing composite material behavior with a microcomputer is described. The formulation assumes a state of generalized plane strain in a material consisting of two or more orthotropic phases. Loading can be mechanical and/or thermal. The theoretical background, computer implementation, and program users guide are described in detail. A sample program is solved showing the required user input and computer generated output.

Bowles, David E.↗