SEARCH · Search NASA
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.
Stress-strain distributions in thick walled cylinder of strain-hardening material, elastic- plastically strained by internal pressure
The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.
Effect of curvature and elastic-plastic material behavior of cylindrical thin shell under critical load
Extension of the classical theory of fracture to viscoelastic and elastic-plastic materials in which the plasticity effects are confined to a narrow band encompassing the crack front. It is suggested that the Griffith-Irwin criterion of fracture, which requires that the energy release rate computed for a given boundary value problem equal the critical threshold, ought to be replaced by a differential equation governing the slow growth of a crack prior to the onset of rapid propagation. A new term which enters the equation of motion in the dissipative media is proportional to the energy lost within the end sections of the crack, and thus it reflects the extent of inelastic behavior of a solid. A concept of 'apparent' surface energy is introduced to account for the geometry-dependent and the rate-dependent phenomena which influence toughness of an inelastic solid. Three hypotheses regarding the condition for fracture in the subcritical range of load are compared.
The Neuber stress-concentration relation for notches in an elastic-plastic material subjected to shear loading was generalized for a crack in a finite plate subjected to tensile loading. An equation was derived which related the linear elastic stress-intensity factor, the applied stress, and two material parameters. The equation was then used as a two-parameter fracture criterion for surface- and through-cracked specimens. Fracture data from the literature on surface- and through-cracked sheet and plate specimens of steel, titanium alloy, titanium weldment, and aluminum alloy tested at room and cryogenic temperature were analyzed according to the proposed equation. For surface cracks, wide ranges of crack-depth to crack-length ratio and crack-depth to specimen-thickness ratio were considered. For through cracks, wide ranges of crack length and specimen width were also considered.-
An extension is proposed of the classical theory of fracture to viscoelastic and elastic-plastic materials in which the plasticity effects are confined to a narrow band encompassing the crack front. It is suggested that the Griffith-Irwin criterion of fracture, which requires that the energy release rate computed for a given boundary value problem equals the critical threshold, ought to be replaced by a differential equation governing the slow growth of a crack prior to the onset of rapid propagation. A new term which enters the equation of motion in the dissipative media is proportional to the energy lost within the end sections of the crack, and thus reflects the extent of inelastic behavior of a solid. A concept of apparent surface energy is introduced to account for the geometry dependent and the rate dependent phenomena which influence toughness of an inelastic solid. Three hypotheses regarding the condition for fracture in the subcritical range of load are compared. These are: (1) constant fracture energy (Cherepanov), (2) constant opening displacement at instability (Morozov) and (3) final stretch criterion (Wnuk).
The elastic-plastic deformation behavior under cyclic stress of a number of different engineering materials was experimentally investigated with the aid of high-precision methods of measuring, some of which had been newly developed. Experiments made with a variety of steels, the titanium alloy Ti-A16-V4, a cobalt (tungsten) alloy, the high-temperature material Nimonic 90 and Dural (A1-Cu) are reported. The theory given in an attempt to explain these experiments is aimed at finding general formulas for the cyclic stress-strain behavior materials.
Successive approximation procedure for solving elastic-plastic problem for torsion of solid prismatic bar made of strain hardening material
Nonlinear boundary value problems of an infinite elastic-plastic plate with a circular hole subjected to pure tension and pure shear at infinity are solved by a method involving Fourier series and finite difference. On the basis of these solutions, the validity of Neuber's relationship between the stress and strain concentration factors for the plane stress problems is examined and a generalized Stowell formula for the stress concentration factor is proposed for problems in which the applied loading may be pure shear as well as pure tension and, furthermore, other stress states. By the same method of solution, the stress distributions around a rigid circular cylindrical inclusion embedded in an infinite rigid-plastic matrix subjected to uniform transverse pure shear and tension are obtained.
Based upon the Principle of Virtual Work and D'Alembert's Principle, the assumed-displacement version of the spatial finite-element method is developed to predict the large deflection transient responses of structures including elastic-plastic, strain hardening, and strain-rate material behavior. The formulations are developed in detail for curved beamlike structures undergoing planar (1) Bernouilli-Euler-type or (2) Timoshenko-type deformation behavior. The resulting equations of motion are solved timewise by a finite-difference numerical procedure. The present predictions are evaluated via several beam and ring examples for which experimental measurements and independent finite-difference predictions in both space and time are available; very good agreement is noted. The consequences of employing several types of timewise finite-difference operators are examined. Also, some comparisons between finite-element predictions and finite-difference predictions are shown to illustrate 'typical comparisons' of efficiency for a given prediction accuracy.
Elastic-plastic shear cracks propagation in two phase composite materials using dislocation model, determining critical fracture strength
The following details calibration of a material model for Al7075-T6511. This aluminum alloy is commonly used across a host of engineering applications. Owing to its widespread prevalence, there is great benefit in improving simulation predictions for this alloy. In the present effort, a calibration is performed of its elastic-plastic response accounting for both rate and temperature dependence. The calibration is informed by a series of tests that include specimens of different geometries tested at different rates and temperatures. All specimens are derived from the same barstock, 3.5 inches in diameter. The fitted model itself uses an anisotropic, Hill yield surface coupled with a Johnson-Cook hardening model. Failure predictions are had by means of a modified Wilkins failure criterion. Following calibration of the material model, a validation exercise is performed against platepuncture experiments. These experiments include multiple probe shapes, probe diameters, and plate thicknesses. The puncture experiments are replicated in simulation with mesh studies performed to assess uncertainty. Key quantities of interest, notably the absorbed energy up to failure, are compared between simulation and experiment providing a means to assess the suitability of the calibration in puncture simulations.
The formulation and check out porblems for a computer code DYNAPLAS, which analyzes the large deflection elastic-plastic dynamic response of stiffened shells of revolution, are presented. The formulation for special discretization is by the finite element method with finite differences being used for the evaluation of the pseudo forces due to material and geometric nonlinearities. Time integration is by the Houbolt method. The stiffeners may be due to concentrated or distributed eccentric rings and spring supports at arbitrary angles around the circumference of the elements. Check out porblems include the comparison of solutions from DYNAPLAS with experimental and other computer solutions for rings, conical and cylindrical shells and a curved panel. A hypothetical submarine including stiffeners and missile tube is studied under a combination of hydrostatic and dynamically applied asymmetrical pressure loadings.
A user-oriented FORTRAN 4 computer program, called JET 3, is presented. The JET 3 program, which employs the spatial finite-element and timewise finite-difference method, can be used to predict the large two-dimensional elastic-plastic transient Kirchhoff-type deformations of a complete or partial structural ring, with various support conditions and restraints, subjected to a variety of initial velocity distributions and externally-applied transient forcing functions. The geometric shapes of the structural ring can be circular or arbitrarily curved and with variable thickness. Strain-hardening and strain-rate effects of the material are taken into account.
A study was made to estimate the magnitudes of elastic stresses and elastic-plastic stresses and strains in sheet or foil laminar composites. Using a model tungsten/80Ni + 20Cr laminar composite and assuming cooling or heating through a temperature range of 80-2000 deg F (26.5-1093.5 deg C), calculated elastic stresses exceeded published or estimated strengths of the constituents. Elastic-plastic stress-strain solutions resulted in lower estimated stress levels but with the concomitant occurrence of sufficiently large strain ranges to suggest possible thermal fatigue problems. Limited experimental studies using tungsten/80Ni + 20Cr foil and sheet laminar composites, slowly cycled between 80 and 2000 deg F (26.5-1093.5 deg C) or rapidly cycled between 80 and 1600 deg F (26.5-871 deg C) produced varying degrees of observable structural damage in from 1 to 11 cycles depending upon temperature transition rate and laminae thickness; these particular results might not occur with other combinations of materials.
An approximate collision analysis, termed the collision-force method, was developed for studying impact-interaction of an engine rotor blade fragment with an initially circular containment ring. This collision analysis utilizes basic mass, material property, geometry, and pre-impact velocity information for the fragment, together with any one of three postulated patterns of blade deformation behavior: (1) the elastic straight blade model, (2) the elastic-plastic straight shortening blade model, and (3) the elastic-plastic curling blade model. The collision-induced forces are used to predict the resulting motions of both the blade fragment and the containment ring. Containment ring transient responses are predicted by a finite element computer code which accommodates the large deformation, elastic-plastic planar deformation behavior of simple structures such as beams and/or rings. The effects of varying the values of certain parameters in each blade-behavior model were studied. Comparisons of predictions with experimental data indicate that of the three postulated blade-behavior models, the elastic-plastic curling blade model appears to be the most plausible and satisfactory for predicting the impact-induced motions of a ductile engine rotor blade and a containment ring against which the blade impacts.
The uniaxial deformation of an elastic-plastic medium containing a doubly periodic square array of circular cylindrical voids is studied under plane-strain conditions. Both the effects of geometrical nonlinearities resulting from large deformation and physical nonlinearities arising from plastic material behavior are included in formulating the problem. A variational principle is used as the basis for implementing a finite-element solution. Results are obtained for the change in void shape and size under increasing overall strain, the overall tensile behavior of the material with voids, and the development of the plastic zone about a void.
Numerical procedures for accurate determination of elastic stress intensity factors for the general two-dimensional crack problem are reviewed. The elastic perfectly plastic state of crack tip deformation is studied by a finite element procedure. Elastic-plastic fields in the immediate vicinity of a crack tip are determined numerically by finite element procedures based on asymptotic studies of crack tip singularities in plastic materials. The small-scale yielding problem is modeled, and expressions for crack tip opening displacement, shear singularity amplitude, and plastic zone extent are derived. A finite element solution to the large-scale yielding of a circumferentially cracked round tension bar is obtained. The three-dimensional aspects of flawed structures and numerical methods of treating them are studied. Ductile fracture mechanisms, in particular crack tip fracture on the microscale, are discussed.