SEARCH · Search NASA

Results for “finite strain”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1974

On numerically accurate finite element solutions in the fully plastic range

A general criterion for testing a mesh with topologically similar repeat units is given, and it is shown that only a few conventional element types and arrangements are suitable for computations in the fully plastic range. Further, a new variational principle, which can easily be incorporated into an existing finite element program, is presented. This makes it possible to perform accurate computations even for element designs that would not normally be suitable. An appendix extends the discussion to elastic-plastic computation at finite strain.

Nagtegaal, J. C.↗

NASA NTRS1977

The relationship between crustal tectonics and internal evolution in the moon and Mercury

The relationship between crustal tectonics and thermal evolution is discussed in terms of the moon and Mercury. Finite strain theory and depth and temperature-dependent thermal expansion are used to evaluate previous conclusions about early lunar history. Factors bringing about core differentiation in the first 0.6 b.y. of Mercurian evolution are described. The influence of concentrating radioactive heat sources located in Mercury's crust on the predicted contraction is outlined. The predicted planetary volume change is explored with regard to quantitative limits on the extent of Mercurian core solidification. Lunar and Mercurian thermal stresses involved in thermal evolution are reviewed, noting the history of surface volcanism. It is concluded that surface faulting and volcanism are closely associated with the thermal evolution of the whole planetary volume. As the planet cools or is heated, several types of tectonic and volcanic effects may be produced by thermal stress occurring in the lithosphere.

Solomon, S. C.↗

NASA NTRS1979

Energy absorption mechanisms during crack propagation in metal matrix composites

The stress distributions around individual fibers in a unidirectional boron/aluminum composite material subjected to axial and transverse loadings are being studied utilizing a generalized plane strain finite element analysis. This micromechanics analysis was modified to permit the analysis of longitudinal sections, and also to incorporate crack initiation and propagation. The analysis fully models the elastoplastic response of the aluminum matrix, as well as temperature dependent material properties and thermal stress effects. The micromechanics analysis modifications are described, and numerical results are given for both longitudinal and transverse models loaded into the inelastic range, to first failure. Included are initially cracked fiber models.

Murphy, D. P.↗

NASA NTRS1981

Analysis of crack propagation as an energy absorption mechanism in metal matrix composites

The crack initiation and crack propagation capability was extended to the previously developed generalized plane strain, finite element micromechanics analysis. Also, an axisymmetric analysis was developed, which contains all of the general features of the plane analysis, including elastoplastic material behavior, temperature-dependent material properties, and crack propagation. These analyses were used to generate various example problems demonstrating the inelastic response of, and crack initiation and propagation in, a boron/aluminum composite.

Adams, D. F.↗

NASA NTRS1983

Characteristics of thermally-induced transverse cracks in graphite epoxy composite laminates

The characteristics of thermally induced transverse cracks in T300/5208 graphite-epoxy cross-ply and quasi-isotropic laminates were investigated both experimentally and analytically. The formation of transverse cracks and the subsequent crack spacing present during cool down to -250 F (116K) and thermal cycling between 250 and -250 F (116 and 394K) was investigated. The state of stress in the vicinity of a transverse crack and the influence of transverse cracking on the laminate coefficient of thermal expansion (CTE) was predicted using a generalized plane strain finite element analysis and a modified shear lag analysis. A majority of the cross-ply laminates experienced transverse cracking during the initial cool down to -250 F whereas the quasi-isotropic laminates remained uncracked. The in situ transverse strength of the 90 degree layers was more than 1.9 times greater than the transverse strength of the unidirectional 90 degree material for all laminates investigated.

Adams, D. S.↗

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.↗