Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear elasticity”

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 127 records · Page 7

Computing the effective elasticity of anisotropic porous media from X-ray computed micro-tomography images

Development and optimization of composite materials designed for thermal protection of NASA’s spacecraft requires the understanding of their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the micro-scale on realistic representations of these composites. The digital micro-structures are generated either synthetically or through X-ray micro-tomography reconstructions. One of the main challenges in the prediction of heatshield material structural response is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the micro-scale. These are driven by the fiber layout within the micro-structure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled through the use of a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme. The MPSA, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA), was discussed in a previous presentation at the 15th USCCN and in reference. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Federico Semeraro↗

Investigation of Low-Cycle Bending Fatigue of AISI 9310 Steel Spur Gears

An investigation of the low-cycle bending fatigue of spur gears made from AISI 9310 gear steel was completed. Tests were conducted using the single-tooth bending method to achieve crack initiation and propagation. Tests were conducted on spur gears in a fatigue test machine using a dedicated gear test fixture. Test loads were applied at the highest point of single tooth contact. Gear bending stresses for a given testing load were calculated using a linear-elastic finite element model. Test data were accumulated from 1/4 cycle to several thousand cycles depending on the test stress level. The relationship of stress and cycles for crack initiation was found to be semilogarithmic. The relationship of stress and cycles for crack propagation was found to be linear. For the range of loads investigated, the crack propagation phase is related to the level of load being applied. Very high loads have comparable crack initiation and propagation times whereas lower loads can have a much smaller number of cycles for crack propagation cycles as compared to crack initiation.

Handschuh, Robert F.↗

Investigating Some Technical Issues on Cohesive Zone Modeling of Fracture

This study investigates some technical issues related to the use of cohesive zone models (CZMs) in modeling fracture processes. These issues include: why cohesive laws of different shapes can produce similar fracture predictions; under what conditions CZM predictions have a high degree of agreement with linear elastic fracture mechanics (LEFM) analysis results; when the shape of cohesive laws becomes important in the fracture predictions; and why the opening profile along the cohesive zone length needs to be accurately predicted. Two cohesive models were used in this study to address these technical issues. They are the linear softening cohesive model and the Dugdale perfectly plastic cohesive model. Each cohesive model constitutes five cohesive laws of different maximum tractions. All cohesive laws have the same cohesive work rate (CWR) which is defined by the area under the traction-separation curve. The effects of the maximum traction on the cohesive zone length and the critical remote applied stress are investigated for both models. For a CZM to predict a fracture load similar to that obtained by an LEFM analysis, the cohesive zone length needs to be much smaller than the crack length, which reflects the small scale yielding condition requirement for LEFM analysis to be valid. For large-scale cohesive zone cases, the predicted critical remote applied stresses depend on the shape of cohesive models used and can significantly deviate from LEFM results. Furthermore, this study also reveals the importance of accurately predicting the cohesive zone profile in determining the critical remote applied load.

Wang, John T.↗

Investigation of Low-Cycle Bending Fatigue of AISI 9310 Steel Spur Gears

An investigation of the low-cycle bending fatigue of spur gears made from AISI 9310 gear steel was completed. Tests were conducted using the single-tooth bending method to achieve crack initiation and propagation. Tests were conducted on spur gears in a fatigue test machine using a dedicated gear test fixture. Test loads were applied at the highest point of single tooth contact. Gear bending stresses for a given testing load were calculated using a linear-elastic finite element model. Test data were accumulated from 1/4 cycle to several thousand cycles depending on the test stress level. The relationship of stress and cycles for crack initiation was found to be semi-logarithmic. The relationship of stress and cycles for crack propagation was found to be linear. For the range of loads investigated, the crack propagation phase is related to the level of load being applied. Very high loads have comparable crack initiation and propagation times whereas lower loads can have a much smaller number of cycles for crack propagation cycles as compared to crack initiation.

Handschuh, Robert F.↗

Methods for predicting crack growth in advanced structures

Damage tolerance design and analysis is widely used in fracture critical military aircraft engine components. Linear elastic fracture mechanics techniques have been developed and verified. These are used to predict the crack propagation lives of complex geometries under mission cycling conditions. Research on methods necessary for the prediction of elevated temperature crack growth in advanced structures is reviewed. These include environmentally assisted time-dependent crack growth, non-linear fracture mechanics parameters for thermal mechanical fatigue crack growth, and finite element modeling of crack growth in composites.

Van Stone, R. H.↗

Fracture mechanics validity limits

Fracture behavior is characteristics of a dramatic loss of strength compared to elastic deformation behavior. Fracture parameters have been developed and exhibit a range within which each is valid for predicting growth. Each is limited by the assumptions made in its development: all are defined within a specific context. For example, the stress intensity parameters, K, and the crack driving force, G, are derived using an assumption of linear elasticity. To use K or G, the zone of plasticity must be small as compared to the physical dimensions of the object being loaded. This insures an elastic response, and in this context, K and G will work well. Rice's J-integral has been used beyond the limits imposed on K and G. J requires an assumption of nonlinear elasticity, which is not characteristic of real material behavior, but is thought to be a reasonable approximation if unloading is kept to a minimum. As well, the constraint cannot change dramatically (typically, the crack extension is limited to ten-percent of the initial remaining ligament length). Rice, et al investigated the properties required of J-type parameters, J(sub x), and showed that the time rate, dJ(sub x)/dt, must not be a function of the crack extension rate, da/dt. Ernst devised the modified-J parameter, J(sub M), that meets this criterion. J(sub M) correlates fracture data to much higher crack growth than does J. Ultimately, a limit of the validity of J(sub M) is anticipated, and this has been estimated to be at a crack extension of about 40-percent of the initial remaining ligament length. None of the various parameters can be expected to describe fracture in an environment of gross plasticity, in which case the process is better described by deformation parameters, e.g., stress and strain. In the current study, various schemes to identify the onset of the plasticity-dominated behavior, i.e., the end of fracture mechanics validity, are presented. Each validity limit parameter is developed in detail, and then data is presented and the various schemes for establishing a limit of the validity are compared. The selected limiting parameter is applied to a set of fracture data showing the improvement of correlation gained.

Lambert, Dennis M.↗

Cell model and elastic moduli of disordered solids - Low temperature limit

The cell theory has been previously employed to compute the equation of state of a disordered condensed system. It is now generalized to include anisotropic stresses. The condition of affine deformation is adopted, transforming an orginally spherical into an ellipsoidal cell. With a Lennard-Jones n-m potential between nonbonded centers, the formal expression for the deformational free energy is derived. It is to be evaluated in the limit of the linear elastic range. Since the bulk modulus in this limit is already known, it is convenient to consider a uniaxial deformation. To begin with, restrictions are made to the low-temperature limit in the absence of entropy contributions. Young's modulus and Poisson's ratio then follow.

Peng, S. T. J.↗

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

Fatigue-Crack-Growth Structural Analysis

Elastic and plastic deformations calculated under variety of loading conditions. Prediction of fatigue-crack-growth lives made with FatigueCrack-Growth Structural Analysis (FASTRAN) computer program. As cyclic loads are applied to initial crack configuration, FASTRAN predicts crack length and other parameters until complete break occurs. Loads are tensile or compressive and of variable or constant amplitude. FASTRAN incorporates linear-elastic fracture mechanics with modifications of load-interaction effects caused by crack closure. FASTRAN considered research tool, because of lengthy calculation times. FASTRAN written in FORTRAN IV for batch execution.

Newman, J. C., Jr.↗

Elastic-Plastic Fracture Mechanics Analysis of Critical Flaw Size in ARES I-X Flange-to-Skin Welds

NASA's Ares 1 Upper Stage Simulator (USS) is being fabricated from welded A516 steel. In order to insure the structural integrity of these welds it is of interest to calculate the critical initial flaw size (CIFS) to establish rational inspection requirements. The CIFS is in turn dependent on the critical final flaw size (CFS), as well as fatigue flaw growth resulting from transportation, handling and service-induced loading. These calculations were made using linear elastic fracture mechanics (LEFM), which are thought to be conservative because they are based on a lower bound, so called elastic, fracture toughness determined from tests that displayed significant plasticity. Nevertheless, there was still concern that the yield magnitude stresses generated in the flange-to-skin weld by the combination of axial stresses due to axial forces, fit-up stresses, and weld residual stresses, could give rise to significant flaw-tip plasticity, which might render the LEFM results to be non-conservative. The objective of the present study was to employ Elastic Plastic Fracture Mechanics (EPFM) to determine CFS values, and then compare these values to CFS values evaluated using LEFM. CFS values were calculated for twelve cases involving surface and embedded flaws, EPFM analyses with and without plastic shakedown of the stresses, LEFM analyses, and various welding residual stress distributions. For the cases examined, the computed CFS values based on elastic analyses were the smallest in all instances where the failures were predicted to be controlled by the fracture toughness. However, in certain cases, the CFS values predicted by the elastic-plastic analyses were smaller than those predicted by the elastic analyses; in these cases the failure criteria were determined by a breakdown in stress intensity factor validity limits for deep flaws (a greater than 0.90t), rather than by the fracture toughness. Plastic relaxation of stresses accompanying shakedown always increases the calculated CFS values compared to the CFS values determined without shakedown. Thus, it is conservative to ignore shakedown effects.

Chell, G. Graham↗

Method and models for R-curve instability calculations

This paper presents a simple method for performing elastic R-curve instability calculations. For a single material-structure combination, the calculations can be done on some pocket calculators. On microcomputers and larger, it permits the development of a comprehensive program having libraries of driving force equations for different configurations and R-curve model equations for different materials. The paper also presents several model equations for fitting to experimental R-curve data, both linear elastic and elastoplastic. The models are fit to data from the literature to demonstrate their viability.

Orange, Thomas W.↗

An approach for relating the results of quantitative nondestructive evaluation to intrinsic properties of high-performance materials

One of the most difficult problems the manufacturing community has faced during recent years has been to accurately assess the physical state of anisotropic high-performance materials by nondestructive means. In order to advance the design of ultrasonic nondestructive testing systems, a more fundamental understanding of how ultrasonic waves travel and interact within the anisotropic material is needed. The relationship between the ultrasonic and engineering parameters needs to be explored to understand their mutual dependence. One common denominator is provided by the elastic constants. The preparation of specific graphite/epoxy samples to be used in the experimental investigation of the anisotropic properties (through the measurement of the elastic stiffness constants) is discussed. Accurate measurements of these constants will depend upon knowledge of refraction effects as well as the direction of group velocity propagation. The continuing effort for the development of improved visualization techniques for physical parameters is discussed. Group velocity images are presented and discussed. In order to fully understand the relationship between the ultrasonic and the common engineering parameters, the physical interpretation of the linear elastic coefficients (the quantities that relate applied stresses to resulting strains) are discussed. This discussion builds a more intuitional understanding of how the ultrasonic parameters are related to the traditional engineering parameters.

Miller, James G.↗

Method and models for R-curve instability calculations

This paper presents a simple method for performing elastic R-curve instability calculations. For a single material-structure combination, the calculations can be done on some pocket calculators. On microcomputers and larger, it permits the development of a comprehensive program having libraries of driving force equations for different configurations and R-curve model equations for different materials. The paper also presents several model equations for fitting to experimental R-curve data, both linear elastic and elastoplastic. The models are fit to data from the literature to demonstrate their viability.

Orange, Thomas W.↗

Finite element analysis of large transient elastic-plastic deformations of simple structures, with application to the engine rotor fragment containment/deflection problem

Assumed-displacement versions of the finite-element method are developed to predict large-deformation elastic-plastic transient deformations of structures. Both the conventional and a new improved finite-element variational formulation are derived. These formulations are then developed in detail for straight-beam and curved-beam elements undergoing (1) Bernoulli-Euler-Kirchhoff or (2) Timoshenko deformation behavior, in one plane. For each of these categories, several types of assumed-displacement finite elements are developed, and transient response predictions are compared with available exact solutions for small-deflection, linear-elastic transient responses. The present finite-element predictions for large-deflection elastic-plastic transient responses are evaluated via several beam and ring examples for which experimental measurements of transient strains and large transient deformations and independent finite-difference predictions are available.

Wu, R. W.↗

Controls-structures-interaction dynamics during RCS control of the Orbiter/SRMS/SSF configuration

During the assembly flights of the Space Station Freedom (SSF), the Orbiter will either dock with the SSF and retract to the final berthed position, or will grapple the SSF using the Shuttle Remote Manipulator System (SRMS) and maneuver the SRMS coupled vehicles to their final berthed position. The SRMS method is expected to take approximately one to one and a half hours to complete and require periodic attitude corrections by either the Orbiter or the SSF reaction control system (RCS) or continuous control by a control moment gyro (CMG) system with RCS desaturation as required. Free drift of the attached vehicles is not currently thought to be acceptable because the desired system attitude will quickly deteriorate due to unbalanced gravity gradient and aerodynamic torques resulting in power generation problems, thermodynamic control problems, and communications problems. This paper deals with the simulation and control of the SRMS during trunnion/latch interaction dynamics and during RCS maneuvers. The SRMS servo drive joints have highly non-linear elastic characteristics which tend to degrade sensitive control strategies. In addition the system natural frequencies are extremely low and depend on the drive joint deflections and SRMS geometric position. The lowest mean period of oscillation for the Orbiter/SRMS/SSF(MB6) system in brakes hold mode positioned near the final berthed position is approximately 120 seconds. A detailed finite element model of the SRMS has been developed and used in a newly developed SRMS systems dynamics simulation to investigate the non-linear transient response dynamics of the Orbiter/SRMS/SSF systems. The present SRMS control strategy of brakes only recommended by the Charles Draper Labs is contrasted with a robust controller developed by the authors. The robust controller uses an optimal inear quadratic regulator (LQR) to optimally place the closed-loop poles of a multivariable continuous-time system within the common region of an open sector with the sector angle plus or minus 45 degrees from the negative real axis, and the left-hand side of a parallel to the imaginary axis in the complex s-plane. This guarantees that the critical damping ratio for the desired control modes is equal to or in excess of 0.707. The matrix sign function is used for solving the Riccati equations which appear in the controller design procedure. Fast and stable algorithms have recently been developed for the computation of the matrix sign function. Simulation results are given which demonstrate the potential CSI involvement for the current SRMS control system and the proposed control system.

Schliesing, J. A.↗

Physical interpretation and application of principles of ultrasonic nondestructive evaluation of high-performance materials

An ultrasonic measurement system employed in the experimental interrogation of the anisotropic properties (through the measurement of the elastic stiffness constants) of the uniaxial graphite-epoxy composites is presented. The continuing effort for the development of improved visualization techniques for physical parameters is discussed. The background is set for the understanding and visualization of the relationship between the phase and energy/group velocity for propagation in high-performance anisotropic materials by investigating the general requirements imposed by the classical wave equation. The consequences are considered when the physical parameters of the anisotropic material are inserted into the classical wave equation by a linear elastic model. The relationship is described between the phase velocity and the energy/group velocity three dimensional surfaces through graphical techniques.

Miller, James G.↗

Sandwich beam compressive test method

Tests presented show that the sandwich beam in four-point bending could be used to obtain reliable compressive elastic constants for graphite/polyimide laminates although some difficulties were encountered measuring composite compressive strengths. Data were obtained for ultimate stress, ultimate strain, Young's modulus, and Poisson's ratio values in 117K, room temperature, and 589K test environments for the HTS1/PMR-15 material system. A total of 36 compressive tests on various laminates were performed. Also, 24 tensile tests on the same laminate orientations were performed to obtain input to the analytical portion of this program. A portion of the beam test section was analyzed using a linear elastic finite element computer program to predict the stress state in the graphite/polyimide composite. The influence of the honeycomb core of this stress state was analytically determined.

Shuart, M. J.↗

A geometrically nonlinear analysis of interlaminar stresses in unsymmetrically laminated plates subjected to uniform thermal loading

An analytical study of interlaminar stresses in unsymmetrically laminated plates is presented. The study examines the linear elastic, large deflection response of square laminated composite plates subjected to uniform thermal loading. Both cross-ply and angle-ply, symmetric and unsymmetric, laminates are examined to evaluate the effects of mismatch between adjacent layers in elastic constants and coefficients of thermal expansion. A geometrically nonlinear kinematic description is used to predict the large out-of-plane (transverse) deflections. The nonlinear, three-dimensional boundary value problems are formulated from elasticity theory and approximate solutions are determined using the finite element method. A global/local analysis procedure is used to obtain improved free edge stress predictions. For the laminates and loading conditions considered, the results indicate that the out-of-plane deflections of the unsymmetric laminates reduce interlaminar shear stresses. These deflections also reduce interlaminar normal stresses in some laminates and increase these stresses for others. The results indicate that as the out-of-plane deflections become large, the differences in interlaminar stress predictions between linear and nonlinear theory can become quite large.

Norwood, D. S.↗