Engineering topics
Publications and source records attributed to Witmer, E. A..
The formulation of the problem of finite elastic-viscoplastic deformation of thin beams is described. The finite element analysis of a representative problem is presented and various solution methods are surveyed.
Two computer programs are described for predicting the transient large deflection elastic viscoplastic responses of thin single layer, initially flat unstiffened or integrally stiffened, Kirchhoff-Lov ductile metal panels. The PLATE 1 program pertains to structural responses produced by prescribed externally applied transient loading or prescribed initial velocity distributions. The collision imparted velocity method PLATE 1 program concerns structural responses produced by impact of an idealized nondeformable fragment. Finite elements are used to represent the structure in both programs. Strain hardening and strain rate effects of initially isotropic material are considered.
Finite element modeling alternatives as well as the utility and limitations of the two dimensional structural response computer code CIVM-JET 4B for predicting the transient, large deflection, elastic plastic, structural responses of two dimensional beam and/or ring structures which are subjected to rigid fragment impact were investigated. The applicability of the CIVM-JET 4B analysis and code for the prediction of steel containment ring response to impact by complex deformable fragments from a trihub burst of a T58 turbine rotor was studied. Dimensional analysis considerations were used in a parametric examination of data from engine rotor burst containment experiments and data from sphere beam impact experiments. The use of the CIVM-JET 4B computer code for making parametric structural response studies on both fragment-containment structure and fragment-deflector structure was illustrated. Modifications to the analysis/computation procedure were developed to alleviate restrictions.
The CIVM-JET 4C computer program for the 'finite strain' analysis of 2 d transient structural responses of complete or partial rings and beams subjected to fragment impact stored on tape as a series of individual files. Which subroutines are found in these files are described in detail. All references to the CIVM-JET 4C program are made assuming that the user has a copy of NASA CR-134907 (ASRL TR 154-9) which serves as a user's guide to (1) the CIVM-JET 4B computer code and (2) the CIVM-JET 4C computer code 'with the use of the modified input instructions' attached hereto.
A method of analysis for thin structures that incorporates finite strain, elastic-plastic, strain hardening, time dependent material behavior implemented with respect to a fixed configuration and is consistently valid for finite strains and finite rotations is developed. The theory is formulated systematically in a body fixed system of convected coordinates with materially embedded vectors that deform in common with continuum. Tensors are considered as linear vector functions and use is made of the dyadic representation. The kinematics of a deformable continuum is treated in detail, carefully defining precisely all quantities necessary for the analysis. The finite strain theory developed gives much better predictions and agreement with experiment than does the traditional small strain theory, and at practically no additional cost. This represents a very significant advance in the capability for the reliable prediction of nonlinear transient structural responses, including the reliable prediction of strains large enough to produce ductile metal rupture.
These structural ring deflections lie essentially in one plane and, hence, are called two-dimensional (2-d). The structural rings may be complete or partial; the former may be regarded as representing a fragment containment ring while the latter may be viewed as a 2-d fragment-deflector structure. These two types of rings may be either free or supported in various ways (pinned-fixed, locally clamped, elastic-foundation supported, mounting-bracket supported, etc.). The initial geometry of each ring may be circular or arbitrarily curved; uniform-thickness or variable-thickness rings may be analyzed. Strain-hardening and strain-rate effects of initially-isotropic material are taken into account. An approximate analysis utilizing kinetic energy and momentum conservation relations is used to predict the after-impact velocities of each fragment and of the impact-affected region of the ring; this procedure is termed the collision-imparted velocity method (CIVM) and is used in the CIVM-JET 5 B program. This imparted-velocity information is used in conjunction with a finite-element structural response computation code to predict the transient, large-deflection, elastic-plastic responses of the ring. Similarly, the equations of motion of each fragment are solved in small steps in time. Provisions are made in the CIVM-JET 5B code to analyze structural ring response to impact attack by from 1 to 3 fragments, each with its own size, mass, translational velocity components, and rotational velocity. The effects of friction between each fragment and the impacted ring are included.
The sheet explosive loading technique (SELT) was employed to obtain elastic-plastic, large deflection 3-d transient and/or permanent strain data on simple well defined structural specimens and materials: initially-flat 6061-T651 aluminum panels with all four sides ideally clamped via integral construction. The SELT loading technique was chosen since it is both convenient and provides "forcing function information" of small uncertainty. These data will be useful for evaluating pertinent 3-d structural response prediction methods.
Design considerations, objectives, and approaches used in containing rotor burst debris are discussed. Methods are given for determining the fracture resistance of various materials used in providing lightweight shielding from fragment impact.
Theoretical methods were developed for predicting the large-deflection elastic-plastic transient structural responses of metal containment or deflector (C/D) structures to cope with rotor burst fragment impact attack. For two-dimensional C/D structures both, finite element and finite difference analysis methods were employed to analyze structural response produced by either prescribed transient loads or fragment impact. For the latter category, two time-wise step-by-step analysis procedures were devised to predict the structural responses resulting from a succession of fragment impacts: the collision force method (CFM) which utilizes an approximate prediction of the force applied to the attacked structure during fragment impact, and the collision imparted velocity method (CIVM) in which the impact-induced velocity increment acquired by a region of the impacted structure near the impact point is computed. The merits and limitations of these approaches are discussed. For the analysis of 3-d responses of C/D structures, only the CIVM approach was investigated.
A user-oriented computer program CIVM-JET 4B is described to predict the large-deflection elastic-plastic structural responses of fragment impacted single-layer: (a) partial-ring fragment containment or deflector structure or (b) complete-ring fragment containment structure. These two types of structures may be either free or supported in various ways. Supports accommodated include: (1) point supports such as pinned-fixed, ideally-clamped, or supported by a structural branch simulating mounting-bracket structure and (2) elastic foundation support distributed over selected regions of the structure. The initial geometry of each partial or complete ring may be circular or arbitrarily curved; uniform or variable thicknesses of the structure are accommodated. The structural material is assumed to be initially isotropic; strain hardening and strain rate effects are taken into account.
The sheet explosive loading technique (SELT) was employed to obtain elastic-plastic, large-deflection transient and/or permanent strain data on simple well-defined structural specimens and materials: initially-flat 6061-T651 aluminum beams with both ends ideally clamped via integral construction. The SELT loading technique was chosen since it is both convenient and provides forcing function information of small uncertainty. These data will be useful for evaluating pertinent structural response prediction methods. A second objective was to obtain high-quality transient-strain data for a well-defined structural/material model subjected to impact by a rigid body of known mass, impact velocity, and geometry; large-deflection, elastic-plastic transient response conditions are of primary interest. The beam with both ends clamped and a steel sphere as the impacting body were chosen. The steel sphere was launched vertically by explosive propulsion to achieve various desired impact velocities. The sphere/beam impact tests resulted in producing a wide range of structural responses and permanent deformations, including rupture of the beam from excessive structural response in two cases. The transient and permanent strain data as well as the permanent deflection data obtained are of high quality and should be useful for checking and evaluating methods for predicting the responses of simple 2-d structures to fragment (sphere) impact. Transient strain data very close to the point of impact were not obtained over as long a time as desirable because the gage(s) in that region became detached during the transient response.
An analysis method termed TEJ-JET is described whereby measured transient elastic and inelastic deformations of an engine-rotor fragment-impacted structural ring are analyzed to deduce the transient external forces experienced by that ring as a result of fragment impact and interaction with the ring. Although the theoretical feasibility of the TEJ-JET concept was established, its practical feasibility when utilizing experimental measurements of limited precision and accuracy remains to be established. The experimental equipment and the techniques (high-speed motion photography) employed to measure the transient deformations of fragment-impacted rings are described. Sources of error and data uncertainties are identified. Techniques employed to reduce data reading uncertainties and to correct the data for optical-distortion effects are discussed. These procedures, including spatial smoothing of the deformed ring shape by Fourier series and timewise smoothing by Gram polynomials, are applied illustratively to recent measurements involving the impact of a single T58 turbine rotor blade against an aluminum containment ring. Plausible predictions of the fragment-ring impact/interaction forces are obtained by one branch of this TEJ-JET method; however, a second branch of this method, which provides an independent estimate of these forces, remains to be evaluated.
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.
An approximate analysis, termed the Collision Imparted Velocity Method (CIVM), was employed for predicting the transient structural responses of containment rings or deflector rings which are subjected to impact from turbojet-engine rotor burst fragments. These 2-d structural rings may be initially circular or arbitrarily curved and may have either uniform or variable thickness; elastic, strain hardening, and strain rate material properties are accommodated. This approximate analysis utilizes kinetic energy and momentum conservation relations in order to predict the after-impact velocities of the fragment and the impacted ring segment. This information is then used in conjunction with a finite element structural response computation code to predict the transient, large deflection responses of the ring. Similarly, the equations of motion for each fragment are solved in small steps in time. Also, some comparisons of predictions with experimental data for fragment-impacted free containment rings are presented.
Mathematical optimization of containment/deflection system would save time, effort, and material as well as afford designer greater opportunity to investigate new ideas and variety of materials.
The transient responses of containment and/or deflection rings to impact from an engine rotor-blade fragment are analyzed. Energy and momentum considerations are employed in an approximate analysis to predict the collision-induced velocities which are imparted to the fragment and to the affected ring segment. This collision analysis is combined with the spatial finite-element representation of the ring and a temporal finite-difference solution procedure to predict the resulting large transient elastic-plastic deformations of containment/deflection rings. Some comparisons with experimental data are given.
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.
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.