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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Shock-wave properties and high-pressure equations of state of geophysically important materials

Shock wave (Hugoniot), shock temperature, and release data are presented for several geophysically important, refractory materials. A sensitive multichannel optical pyrometer was developed to measure shock temperatures (2500 to 5600 K at pressures from 48 to 117 GPa) in anorthite (CaAl2Si2O8) glass. Shock temperatures of 3750 to 6000 K at pressures from 140 to 182 GPa were measured in calcium oxide (CaO). Temperature data were used to constrain the energetics of the B1-B2 phase transition at 70 GPa in CaO, and to construct a finite strain equation of state for CaO consistent with previous Hugoniot data. The CaO equation of state was used with equation of state parameters of other oxides to construct a theoretical mixed oxide Hugoniot of anorthite, which is in agreement with new Hugoniot data above about 50 GPa, determined using experimental techniques developed. The mixed oxide model, however, overestimates the shock temperatures, and does not accurately predict measured release paths.

Boslough, M. B.↗
NASA NTRS1985

Adiabatic thermal models for planetary bodies

In a number of recent experiments it was found that the logarithmic derivative with respect to volume of the adiabatic temperature increase with pressure P to be an approximately constant quantity n. It was found that n decreases slightly with temperature, to be virtually unaffected by increasing pressure and to take values between 4 and 8 for a wide variety of materials. It is shown that these findings can be substantiated from thermodynamic arguments, finite strain theory, atomic potential theory and experimental data on the thermal expansion coefficient and the bulk modulus B. It will be shown that n is independent of pressure if it is exactly equal to dB/dP + 1. For these materials d log gamma/d log v = -1, where gamma is the thermodynamic Gruneisenparameter. It will increase with P during an isothermal transformation if n dB/dP + 1 and decrease of n dB/dP + 1. For most materials n is close to dB/dP and the changes will be slight if pressures do not become too extreme. During an adiabatic transformation n is virtually constant. Adiabatic thermal models for planetary bodies were calculated and are presented.

Spohn, T.↗
NASA NTRS1986

Reinforcement of composite laminate free edges with U-shaped caps

Generalized plane strain finite element analysis is used to predict reduction of interlaminar normal stresses when a U-shaped cap is bonded to the edge of a laminate. Three-dimensional composite material failure criteria are used in a progressive laminate failure analysis to predict failure loads of laminates with different edge cap designs. In an experimental program, symmetric 11-layer graphite-epoxy laminates with a one-layer cap of Kevlar-epoxy cloth are shown to be 130 to 140 percent stronger than uncapped laminates under static tensile and tension-tension fatigue loading. In addition, the coefficient of variation of the static tensile failure load decreases from 24 to 8 percent when edge caps are added. The predicted failure load calculated with the finite element results is 10 percent lower than the actual failure load. For both capped and uncapped laminates, actual failure loads are much lower than those predicted using classical lamination theory stresses and a two-dimensional failure criterion. Possible applications of the free edge reinforcement concept are described, and future research is suggested.

Howard, W. E.↗
NASA NTRS1986

Computational solid mechanics (finite elements and boundary elements) - Present status and future directions

Computational finite-element and boundary-element methods are reviewed, and their application to the mechanics of solids is discussed. Stability conditions for general FEMs are considered in addition to the use of least-order, stable, invariant, or hybrid/mixed isoparametric elements as alternatives to the displacement-based isoparametric elements. The use of symbolic manipulation, adaptive mesh refinement, transient dynamic response, and boundary-element methods for linear elaslticity and finite-strain problems of inelastic materials are also discussed.

Atluri, S. N.↗
NASA NTRS1987

A nonlinear viscoelastic constitutive equation - Yield predictions in multiaxial deformations

Yield stress predictions of a nonlinear viscoelastic constitutive equation for amorphous polymer solids have been obtained and are compared with the phenomenological von Mises yield criterion. Linear viscoelasticity theory has been extended to include finite strains and a material timescale that depends on the instantaneous temperature, volume, and pressure. Results are presented for yield and the correct temperature and strain-rate dependence in a variety of multiaxial deformations. The present nonlinear viscoelastic constitutive equation can be formulated in terms of either a Cauchy or second Piola-Kirchhoff stress tensor, and in terms of either atmospheric or hydrostatic pressure.

Shay, R. M., Jr.↗
NASA NTRS1988

Origin of petrofabrics and magnetic anisotropy in ordinary chondrites

Three-dimensional finite strain and magnetic susceptibility anisotropy have been determined for 15 ordinary chondrites. The axes of strain and magnetic ellipsoids roughly correlate in both magnitude and orientation. The shapes of these ellipsoids are generally oblate spheroids that define a dominant foliation and a weak lineation. These characteristics suggest deformation involving uniaxial compaction. The degree of uniaxial deformation correlates with intensity of shock, as indicated by optical, TEM and chemical criteria. These data, plus the lack of a relationship between foliation and metamorphic history, indicate that dynamic processes, i.e., impacts, produced planar deformation fabrics in chondrites.

Sneyd, Deana S.↗
NASA NTRS1988

Deep crustal deformation by sheath folding in the Adirondack Mountains, USA

As described by McLelland and Isachsen, the southern half of the Adirondacks are underlain by major isoclinal (F sub 1) and open-upright (F sub 2) folds whose axes are parallel, trend approximately E-W, and plunge gently about the horizontal. These large structures are themselves folded by open upright folds trending NNE (F sub 3). It is pointed out that elongation lineations in these rocks are parallel to X of the finite strain ellipsoid developed during progressive rotational strain. The parallelism between F sub 1 and F sub 2 fold axes and elongation lineations led to the hypothesis that progressive rotational strain, with a west-directed tectonic transport, rotated earlier F sub 1-folds into parallelism with the evolving elongation lineation. Rotation is accomplished by ductile, passive flow of F sub 1-axes into extremely arcuate, E-W hinges. In order to test these hypotheses a number of large folds were mapped in the eastern Adirondacks. Other evidence supporting the existence of sheath folds in the Adirondacks is the presence, on a map scale, of synforms whose limbs pass through the vertical and into antiforms. This type of outcrop pattern is best explained by intersecting a horizontal plane with the double curvature of sheath folds. It is proposed that sheath folding is a common response of hot, ductile rocks to rotational strain at deep crustal levels. The recognition of sheath folds in the Adirondacks reconciles the E-W orientation of fold axes with an E-W elongation lineation.

Mclelland, J. M.↗
NASA NTRS1990

Static compression of Ca(OH)2 at room temperature - Observations of amorphization and equation of state measurements to 10.7 GPa

X-ray diffraction measurements are reported for Ca(OH)2 portlandite as it is compressed to 37.6 GPa in the diamond cell at room temperature. Between 10.7 and 15.4 GPa crystalline Ca(OH)2 transforms to a glass, and on decompression the glass recrystallizes between 3.6 and 5.1 GPa. Below pressures of 10.7 GPa the elastic compression of crystalline Ca(OH)2 was measured. A finite strain analysis of these data shows that the isothermal bulk modulus and its pressure derivative are 37.8 + or - 1.8 GPa and 5.2 + or - 0.7 at zero pressure. The change in the unit cell dimensions indicates that the linear incompressibilities of Ca(OH)2 differ by a factor of three.

Meade, Charles↗
NASA NTRS1993

Computation of strain energy release rates for skin-stiffener debonds modeled with plate elements

An efficient method for predicting the strength of debonded composite skin-stiffener configurations is presented. This method, which is based on fracture mechanics, models the skin and the stiffener with two-dimensional (2D) plate elements instead of three-dimensional (3D) solid elements. The skin and stiffener flange nodes are tied together by two modeling techniques. In one technique, the corresponding flange and skin nodes are required to have identical translational and rotational degrees-of-freedom. In the other technique, the corresponding flange and skin nodes are only required to have identical translational degrees-of-freedom. Strain energy release rate formulas are proposed for both modeling techniques. These formulas are used for skin-stiffener debond cases with and without cylindrical bending deformations. The cylindrical bending results are compared with plane-strain finite element results. Excellent agreement between the two sets of results is obtained when the second technique is used. Thus, from these limited studies, a preferable modeling technique for skin-stiffener debond analysis using plate elements is established.

Wang, J. T.↗
NASA NTRS1993

Inelastic deformation of metal matrix composites

A theoretical model capable of predicting the thermomechanical response of continuously reinforced metal matrix composite laminates subjected to multiaxial loading was developed. A micromechanical model is used in conjunction with nonlinear lamination theory to determine inelastic laminae response. Matrix viscoplasticity, residual stresses, and damage to the fiber/matrix interfacial zone are explicitly included in the model. The representative cell of the micromechanical model is considered to be in a state of generalized plane strain, enabling a quasi two-dimensional analysis to be performed. Constant strain finite elements are formulated with elastic-viscoplastic constitutive equations. Interfacial debonding is incorporated into the model through interface elements based on the interfacial debonding theory originally presented by Needleman, and modified by Tvergaard. Nonlinear interfacial constitutive equations relate interfacial tractions to displacement discontinuities at the interface. Theoretical predictions are compared with the results of an experimental program conducted on silicon carbide/titanium (SiC/Ti) unidirectional, (O4), and angle-ply, (+34)(sub s), tubular specimens. Multiaxial loading included increments of axial tension, compression, torque, and internal pressure. Loadings were chosen in an effort to distinguish inelastic deformation due to damage from matrix plasticity and separate time-dependent effects from time-independent effects. Results show that fiber/matrix debonding is nonuniform throughout the composite and is a major factor in the effective response. Also, significant creep behavior occurs at relatively low applied stress levels at room temperature.

Lissenden, C. J.↗
NASA NTRS1992

Modeling of Thermal Barrier Coatings

The project examined the effectiveness of studying the creep behavior of thermal barrier coating system through the use of a general purpose, large strain finite element program, NIKE2D. Constitutive models implemented in this code were applied to simulate thermal-elastic and creep behavior. Four separate ceramic-bond coat interface geometries were examined in combination with a variety of constitutive models and material properties. The reason for focusing attention on the ceramic-bond coat interface is that prior studies have shown that cracking occurs in the ceramic near interface features which act as stress concentration points. The model conditions examined include: (1) two bond coat coefficient of thermal expansion curves; (2) the creep coefficient and creep exponent of the bond coat for steady state creep; (3) the interface geometry; and (4) the material model employed to represent the bond coat, ceramic, and superalloy base.

Ferguson, B. L.↗
NASA NTRS1999

Testing and Analysis of Composite Skin/Stringer Debonding Under Multi-Axial Loading

Damage mechanisms in composite bonded skin/stringer constructions under uniaxial and biaxial (in-plane/out- of-plane) loading conditions were examined. Specimens consisted of a tapered composite flange bonded onto a composite skin. Tests were performed under monotonic loading conditions in tension, three-point bending, and combined tension/bending . For combined tension/bending testing, a unique servohydraulic load frame was used that was capable of applying both in-plane tension and out-of-plane bending loads simultaneously. Specimen edges were examined on the microscope to document the damage occurrence and to identify typical damage patterns. The observations showed that, for all three load cases, failure initiated in the flange, near the flange tip, causing the flange to almost fully debond from the skin. A two-dimensional plane-strain finite element model was developed to analyze the different test cases using a geometrically nonlinear solution. For all three loading conditions, principal stresses exceeded the transverse strength of the material in the flange area. Additionally, delaminations of various lengths were simulated in two locations where delaminations were observed. The analyses showed that unstable delamination propagation is likely to occur in one location at the loads corresponding to matrix ply crack initiation for all three load cases.

Krueger, Ronald↗
NASA NTRS1998

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗
NASA NTRS2003

Tensor Fields for Use in Fractional-Order Viscoelasticity

To be able to construct viscoelastic material models from fractional0order differentegral equations that are applicable for 3D finite-strain analysis requires definitions for fractional derivatives and integrals for symmetric tensor fields, like stress and strain. We define these fields in the body manifold. We then map them ito spatial fields expressed in terms of an Eulerian or Lagrangian reference frame where most analysts prefer to solve boundary problems.

Freed, Alan D.↗
NASA NTRS2004

Debonding in Composite Skin/Stringer Configurations Under Multi-Axial Loading

The objective of this work was to investigate the damage mechanisms in composite bonded skin/stringer constructions under uniaxial and biaxial (in-plane/out-of-plane) loading conditions as typically experienced by aircraft crown fuselage panels. The specimens for all tests were identical and consisted of a tapered composite flange, representing a stringer or frame, bonded onto a composite skin. Tests were performed under monotonic loading conditions in tension, three-point bending, and combined tension/bending to evaluate the debonding mechanisms between the skin and the bonded stringer. For combined tension/bending testing, a unique servohydraulic load frame was used that was capable of applying both loads simultaneously. Microscopic investigations of the specimen edges were used to document the damage occurrence and to identify typical damage patterns. The observations showed that, for all three load cases, failure initiated in the flange near the flange tip causing the flange to almost fully debond from the skin. A two-dimensional plain-strain finite element model was developed to analyze the different test cases using a geometrically nonlinear solution. For all three loading conditions, principal stresses exceeded the transverse strength of the material in the flange area. Additionally, delaminations of various lengths were simulated in the locations where delaminations were experimentally observed. The analyses showed that unstable delamination propagation is likely to occur at the loads corresponding to matrix ply crack initiation for all three loadings.

Cvitkovich, Michael K.↗
NASA NTRS1983

Two-dimensional Finite Element Modeling for Modeling Tectonic Stress and Strain

Techniques of finite element analysis in two dimensional plane strain were applied to problems of geophysics and tectonics. More specifically, the flexibility of the finite element method was employed to address problems involving geological complexity and fault interactions. The modeling of effective anisotropy in material elastic properties proved useful in describing the deformation of faulted crustal blocks. The applications of this modeling work to problems of actual tectonics in southern California was explored. Preliminary models show encouraging agreement with measured tectonic strain in this region, and modeling work was done to gain an understanding of the stress state in a locked fault region with future seismic potential.

Lyzenga, G. A.↗