Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Elastic-plastic deformation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Inelastic Strain Analysis of Solder Joint in NASA Fatigue Specimen

The solder fatigue specimen designed by NASA-GSFC/UNISYS is analyzed in order to obtain the inelastic strain history during two different representative temperature cycles specified by UNISYS. In previous reports (dated July 25, 1990, and November 15, 1990), results were presented of the elastic-plastic and creep analysis for delta T = 31 C cycle, respectively. Subsequent results obtained during the current phase, from viscoplastic finite element analysis of the solder fatigue specimen for delta T = 113 C cycle are summarized. Some common information is repeated for self-completeness. Large-deformation continuum formulations in conjunction with a standard linear solid model is utilized for modeling the solder constitutive creep-plasticity behavior. Relevant material properties are obtained from the literature. Strain amplitudes, mean strains, and residual strains (as well as stresses) accumulated due to a representative complete temperature cycle are obtained as a result of this analysis. The partitioning between elastic strains, time-independent inelastic (plastic) strains, and time-dependent inelastic (creep) strains is also explicitly obtained for two representative cycles. Detailed plots are presented for two representative temperature cycles. This information forms an important input for fatigue damage models, when predicting the fatigue life of solder joints under thermal cycling

Dasgupta, Abhijit↗
NASA NTRS1974

A finite element analysis of fatigue crack closure

Experiments have shown that fatigue cracks close at positive loads during constant-amplitude load cycling. The crack-closure phenomenon is caused by residual plastic deformations remaining in the wake of an advancing crack tip. The present paper is concerned with the application of a two-dimensional, nonlinear, finite-element analysis for predicting crack-closure and crack-opening stresses during cyclic loading. A two-dimensional finite-element computer program, which accounts for both elastic-plastic material behavior and changing boundary conditions associated with crack extension and intermittent contact of the crack surfaces under cyclic loading, has been developed. An efficient technique to account for changing boundary conditions was also incorporated into the nonlinear analysis program. This program was subsequently used to study crack extension and crack closure under constant-amplitude and two-level block loading. The calculated crack-closure and crack-opening stresses were qualitatively consistent with experimental observations.

Newman, J. C., Jr.↗
NASA NTRS1988

An elastic-plastic finite-element analysis of the J-resistance curve using a CTOD criterion

The fracture process in compact and bend specimens was simulated using a two-dimensional finite-element analysis of the J-resistance (J-R) curve and a CTOD criterion. The J-R curves were calculated from the numerical results for each specimen type using several different methods. In general, the J-R curves obtained for the bend specimens were found to be higher than those for the compact specimens, especially beyond maximum load. However, below the maximum load, the modified deformation theory of plasticity and the contour-integral J(r) method resulted in very similar J-R curves for both specimen types.

Newman, J. C., Jr.↗
NASA NTRS1990

Three-dimensional elastic-plastic analysis of shallow cracks in single-edge-crack-tension specimens

Three dimensional, elastic-plastic, finite element results are presented for single-edge crack-tension specimens with several shallow crack-length-to-width ratios (0.05 less than or equal to a/W less than or equal to 0.5). Results showed the need to model the initial yield plateau in the stress-strain behavior to accurately model deformation of the A36 steel specimens. The crack-tip-opening-displacement was found to be linearly proportional to the crack-mouth-opening displacement. A new deformation dependent plastic-eta factor equation is presented for calculating the J-integral from test load-displacement records. This equation was shown to be accurate for all crack lengths considered.

Shivakumar, Kunigal N.↗
NASA NTRS1990

Three-dimensional elastic-plastic analysis of shallow cracks in single-edge-crack-tension specimens

Three-dimensional, elastic-plastic, finite-element results are presented for single-edge crack-tension specimens with several shallow crack-length-to-width ratios (0.05 less than or equal to a/W less than or equal to 0.5). Results showed the need to model the initial yield plateau in the stress-strain behavior to accurately model deformation of the A36 steel specimens. The crack-tip-opening-displacement was found to be linearly proportional to the crack-mouth-opening displacement. A new deformation dependent plastic-eta factor equation is presented for calculating the J-integral from test load-displacement records. This equation was shown to be accurate for all crack lengths considered.

Shivakumar, Kunigal N.↗
NASA NTRS1979

An elastic-plastic analysis of fatigue crack closure in modes I and II

In this paper, an efficient elastic-plastic finite element procedure to analyze crack-closure and, its effects on fatigue crack growth under general spectrum loading, is presented. A hybrid-displacement finite element procedure is used to properly treat the stress and strain singularities near the crack-tip; and crack-growth under cyclic loading is simulated by the translation of certain 'core' elements, near the crack-tip, in which proper stress and strain singularities were embedded. Both pure mode I and mode II types of cyclic loading are considered. In the mode I case, four types of cyclic loading, viz., constant amplitude block loading, high-to-low block loading, low-to-high block, and a single overload in an otherwise constant amplitude block loading are considered; whereas in mode II, only a constant amplitude block loading is considered. Detailed results are presented for the crack-closure and opening-stresses, crack-surface deformation profiles, etc., in each case. Certain observations, based on the present numerical results, concerning various factors that cause crack growth acceleration or retardation under general spectrum loading, are presented and discussed.

Nakagaki, M.↗
DOE OSTI2025

Non-linear elastic-plastic behavior of the invert glass lithium phosphorous oxynitride (LiPON)

Micro and nano-scale mechanical behavior of network/ionic glasses is dictated by their composition and the number of constraints per structural unit. From this perspective, the glass LiPON presents the opportunity for enhancement of the microscale ductility due to its rather unconstrained orthophosphate structure. Here, we use instrumented nanoindentation with different tip geometries to investigate the mechanical response of LiPON glass. The results reveal that the elastic modulus of LiPON is not constant and depends on pressure. With the method utilizing spherical nanoindentation and continuous stiffness measurement (CSM) we determine the yield point of LiPON. We propose the Drucker-Prager type of yield criterion for LiPON and estimate its yield stress in compression as 2.4 GPa. There exists another stress threshold, however, around 490 MPa, at which the elastic deformation becomes non-linear, and this can be mistaken for the yield stress.

Constitutive model↗
NASA NTRS1973

Elastoplastic transverse properties of a unidirectional fibre reinforced composite

An approximate theoretical study is made to predict the overall transverse elastoplastic properties of a unidirectional fibre reinforced composite. Based on the simplest deformation theory of plasticity, the overall transverse polyaxial elastoplastic constitutive equations for an incompressible composite are obtained. A simple formula is suggested to evaluate the average stress in the fibres. Furthermore, with the prediction of the overall transverse elastic moduli of a compressible composite, the overall transverse uniaxial stress-strain curve for a compressible composite in the entire elastic-plastic range is obtained. A detailed example is given for a Boron-Aluminum composite. Comparison is made with the experimental and theoretical work done by General Dynamics and Adams respectively.

Huang, W.-C.↗
NASA NTRS1983

Application of an Uncoupled Elastic-plastic-creep Constitutive Model to Metals at High Temperature

A uniaxial, uncoupled constitutive model to predict the response of thermal and rate dependent elastic-plastic material behavior is presented. The model is based on an incremental classicial plasticity theory extended to account for thermal, creep, and transient temperature conditions. Revisions to he combined hardening rule of the theory allow for better representation of cyclic phenomenon including the high rate of strain hardening upon cyclic reyield and cyclic saturation. An alternative approach is taken to model the rate dependent inelastic deformation which utilizes hysteresis loops and stress relaxation test data at various temperatures. The model is evaluated and compared to experiments which involve various thermal and mechanical load histories on 5086 aluminum alloy, 304 stainless steel and Hastelloy-X.

Haisler, W. E.↗
NASA NTRS1989

Elevated temperature crack growth

Alloy 718 crack growth experiments were conducted to assess the ability of the selected path-independent (P-I) integrals to describe the elevated temperature crack growth behavior. These tests were performed on single edge notch (SEN) specimens under displacement control with multiple extensometers to monitor the specimen and crack mouth opening displacement (CMOD). The displacements in these tests were sufficiently high to induce bulk cyclic inelastic deformation of the specimen. Under these conditions, the linear elastic fracture mechanics (LEFM) parameter K does not correlate the crack growth data. The experimentally measured displacement gradients at the end of specimen gage length were used as the boundary conditions in elastic-plastic finite element method (FEM) analyses. These analyses were performed with a node release approach using CYANIDE, a GEAE FEM code, which included a gap element which is capable of efficiently simulating crack closure. Excellent correlation was obtained between the experimentally measured and predicted variation of stress and CMOD with crack length and the stress-CMOD loops for Alloy 718 tests conducted at 538 C. This confirmed the accuracy of the FEM crack growth simulation approach. The experimentally measured crack growth rate data correlated well the selected P-I integrals. These investigations have produced significant progress in developing P-I integrals as non-linear fracture mechanics parameters. The results suggest that this methodology has the potential of accurately describing elevated temperature crack growth behavior under the combined influence of thermal cycling and bulk elastic-inelastic deformation states.

Kim, K. S.↗
NASA NTRS1992

Finite-element analyses and fracture simulation in thin-sheet aluminum alloy

A two-dimensional, elastic-plastic finite-element analysis was used with a critical crack-tip-opening angle (CTOA) fracture criterion to model stable crack growth in thin-sheet 2024-T3 aluminum alloy under monotonic loading after precracking at different cyclic stress levels. Tests were conducted on three types of specimens: middle-crack, three-hole-crack and blunt-notch tensile specimens. An experiment technique was developed to measure CTOA during crack growth initiation and stable tearing using a high-resolution video camera and recorder. Crack front shapes were also measured during initiation and stable tearing using a fatigue marker-load technique. Three-dimensional elastic-plastic finite-element analyses of these crack shapes for stationary cracks were conducted to study the crack-front opening displacements. Predicted load against crack extension on middle-crack tension specimens agreed well with test results even for large-scale plastic deformations. The analyses were able to predict the effects of specimen size and precracking stress history on stable tearing. Predicted load against load-line displacements agreed well with test results up to maximum load bu the analyses tended to overpredict displacements as crack grew beyond the maximum load under displacement-controlled conditions. During the initiation phase, the measured CTOA values were high but decreased and remained nearly constant after a small amount of stable tearing. The constant value of CTOA agree well with the calculated value from the finite-element analysis. The larger CTOA values measured at the sheet surface during the initiation phase may be associated with the crack tunneling observed in the tests. Three-dimensional analyses for nonstraight crack fronts predicted much higher displacements near the free surface than in the interior.

Newman, J. C., Jr.↗
NASA NTRS1993

A numerical model for predicting crack path and modes of damage in unidirectional metal matrix composites

A finite element-based numerical technique has been developed to simulate damage growth in unidirectional composites. This technique incorporates elastic-plastic analysis, micromechanics analysis, failure criteria, and a node splitting and node force relaxation algorithm to create crack surfaces. Any combination of fiber and matrix properties can be used. One of the salient features of this technique is that damage growth can be simulated without pre-specifying a crack path. In addition, multiple damage mechanisms in the forms of matrix cracking, fiber breakage, fiber-matrix debonding and plastic deformation are capable of occurring simultaneously. The prevailing failure mechanism and the damage (crack) growth direction are dictated by the instantaneous near-tip stress and strain fields. Once the failure mechanism and crack direction are determined, the crack is advanced via the node splitting and node force relaxation algorithm. Simulations of the damage growth process in center-slit boron/aluminum and silicon carbide/titanium unidirectional specimens were performed. The simulation results agreed quite well with the experimental observations.

Bakuckas, J. G.↗
NASA NTRS1993

A fracture criterion for widespread cracking in thin-sheet aluminum alloys

An elastic-plastic finite-element analysis was used with a critical crack-tip-opening angle (CTOA) fracture criterion to model stable crack growth in thin-sheet 2024-T3 aluminum alloy panels with single and multiple-site damage (MSD) cracks. Comparisons were made between critical angles determined from the analyses and those measured with photographic methods. Calculated load against crack extension and load against crack-tip displacement on single crack specimens agreed well with test data even for large-scale plastic deformations. The analyses were also able to predict the stable tearing behavior of large lead cracks in the presence of stably tearing MSD cracks. Small MSD cracks significantly reduced the residual strength for large lead cracks.

Newman, J. C., Jr.↗
NASA NTRS1993

Elastic-plastic finite element analyses of an unidirectional, 9 vol percent tungsten fiber reinforced copper matrix composite

Micromechanical modeling via elastic-plastic finite element analyses were performed to investigate the effects that the residual stresses and the degree of matrix work hardening (i.e., cold-worked, annealed) have upon the behavior of a 9 vol percent, unidirectional W/Cu composite, undergoing tensile loading. The inclusion of the residual stress-containing state as well as the simulated matrix material conditions proved to be significant since the Cu matrix material exhibited plastic deformation, which affected the subsequent tensile response of the composite system. The stresses generated during cooldown to room temperature from the manufacturing temperature were more of a factor on the annealed-matrix composite, since they induced the softened matrix to plastically flow. This event limited the total load-carrying capacity of this matrix-dominated, ductile-ductile type material system. Plastic deformation of the hardened-matrix composite during the thermal cooldown stage was not considerable, therefore, the composite was able to sustain a higher stress before showing any appreciable matrix plasticity. The predicted room temperature, stress-strain response, and deformation stages under both material conditions represented upper and lower bounds characteristic of the composite's tensile behavior. The initial deformation stage for the hardened material condition showed negligible matrix plastic deformation while for the annealed state, its initial deformation stage showed extensive matrix plasticity. Both material conditions exhibited a final deformation stage where the fiber and matrix were straining plastically. The predicted stress-strain results were compared to the experimental, room temperature, tensile stress-strain curve generated from this particular composite system. The analyses indicated that the actual thermal-mechanical state of the composite's Cu matrix, represented by the experimental data, followed the annealed material condition.

Sanfeliz, Jose G.↗
NASA NTRS1993

Current capabilities for simulating the extreme distortion of thin structures subjected to severe impacts

The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Key, Samuel W.↗