SEARCH · Engineering Papers
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.
Program valuable in analyzing and sizing filamentary composite panels. Panel Analysis and Sizing code (PASCO) for buckling and vibration analysis and sizing of prismatic structure having arbitrary cross section primarily intended for analyzing and sizing stiffened panels made of laminated orthotropic materials. PASCO written in FORTRAN IV.
Phenolic Impregnated Carbon Ablator was the heatshield material for the Stardust probe and is also a candidate heatshield material for the Orion Crew Module. As part of the heatshield qualification for Orion, physical and thermal properties were measured for newly manufactured material, included emissivity, heat capacity, thermal conductivity, elemental composition, and thermal decomposition rates. Based on these properties, an ablation and thermal-response model was developed for temperatures up to 3500 K and pressures up to 100 kPa. The model includes orthotropic and pressure-dependent thermal conductivity. In this work, model validation is accomplished by comparison of predictions with data from many arcjet tests conducted over a range of stagnation heat flux and pressure from 107 Watts per square centimeter at 2.3 kPa to 1100 Watts per square centimeter at 84 kPa. Over the entire range of test conditions, model predictions compare well with measured recession, maximum surface temperatures, and in depth temperatures.
Stress-intensity factors are calculated for a cracked infinite sheet adhesively bonded to a stringer, and debonding of the adhesive layer is predicted. The stringer is modeled as a semiinfinite sheet. Adhesive nonlinearity is also included. Both the sheet and stringer are treated as homogeneous, orthotropic materials, a set of integral equations is formulated and solved to obtain the adhesive shear stresses and crack-tip stress-intensity factors. Adhesive debonding is predicted using a rupture criterion based on the combined adhesive stresses. A through-the-thickness crack is located in the infinite sheet perpendicular to the edge of the stringer. When the crack is not under the stringer, the debond extends along the edge of the stringer. When the crack tip is beneath the stringer, the debond extends to the crack tip, then along the edge of the stringer. Stress levels required for debond initiation decrease as the crack tip is moved beneath the stringer. With a nonlinear adhesive, the debond initiates at higher applied stress levels than in linear adhesive cases. Compared with the linear adhesive solution, modeling a nonlinear adhesive causes the stress-intensity factor to decrease when debonding is included.
Understanding the high velocity impact response of polymer matrix composites with complex architectures is critical to many aerospace applications, including engine fan blade containment systems where the structure must be able to completely contain fan blades in the event of a blade-out. Despite the benefits offered by these materials, the complex nature of textile composites presents a significant challenge for the prediction of deformation and damage under both quasi-static and impact loading conditions. The relatively large mesoscale repeating unit cell (in comparison to the size of structural components) causes the material to behave like a structure rather than a homogeneous material. Impact experiments conducted at NASA Glenn Research Center have shown the damage patterns to be a function of the underlying material architecture. Traditional computational techniques that involve modeling these materials using smeared homogeneous, orthotropic material properties at the macroscale result in simulated damage patterns that are a function of the structural geometry, but not the material architecture. In order to preserve heterogeneity at the highest length scale in a robust yet computationally efficient manner, and capture the architecturally dependent damage patterns, a previously-developed subcell modeling approach where the braided composite unit cell is approximated as a series of four adjacent laminated composites is utilized. This work discusses the implementation of the subcell methodology into the commercial transient dynamic finite element code LS-DYNA (Livermore Software Technology Corp.). Verification and validation studies are also presented, including simulation of the tensile response of straight-sided and notched quasi-static coupons composed of a T700/PR520 triaxially braided [0deg/60deg/-60deg] composite. Based on the results of the verification and validation studies, advantages and limitations of the methodology as well as plans for future work are discussed.
The results of an experimental and numerical study of the effects of initial imperfections on the buckling response and failure of unstiffened thin-walled compression-loaded graphite-epoxy cylindrical shells are presented. The shells considered in the study have six different orthotropic or quasi-isotropic shell-wall laminates and two different shell-radius-to-thickness ratios. The numerical results include the effects of geometric shell-wall mid-surface imperfections, shell-wall thickness variations, local shell-wall ply-gaps associated with the fabrication process, shell-end geometric imperfections, nonuniform end loads, and the effects of elastic boundary conditions. Selected cylinder parameter uncertainties were also considered. Results that illustrate the effects of imperfections and uncertainties on the nonlinear response characteristics, buckling loads and failure the shells are presented. In addition, a common failure analysis is used to predict material failures in the shells.
An optimization study of composite stiffened cylinders is discussed. The mathematical model for the buckling has been coupled successfully with the optimization program AESOP. The buckling analysis is based on the use of the smeared theory for the buckling of stiffened orthotropic cylindrical shells. The loading, radius, and length of the cylinder are assumed to be known parameters. An optimum solution gives the value of cross-sectional dimensions and laminate orientations. The different types of buckling modes are identified. Mathematical models are developed to show the relationships of the parameters.
A detailed evaluation of the suitability of the Iosipescu specimen tested in the modified Wyoming fixture is presented. An experimental investigation using conventional strain gage instrumentation and moire interferometry is performed. A finite element analysis of the Iosipescu shear test for unidirectional and cross-ply composites is used to assess the uniformity of the shear stress field in the vicinity of the notch, and demonstrate the effect of the nonuniform stress field upon the strain gage measurements used for the determination of composite shear moduli. From the test results for graphite-epoxy laminates, it is shown that the proximity of the load introduction point to the test section greatly influences the individual gage readings for certain fiber orientations but the effect upon shear modulus measurement is relatively unimportant. A numerical study of the load contact effect shows the sensitivity of some fiber configurations to the specimen/fixture contact mechanism and may account for the variations in the measured shear moduli. A comparison of the strain gage readings from one surface of a specimen with corresponding data from moire interferometry on the opposite face documented an extreme sensitivity of some fiber orientations to eccentric loading which induced twisting and yielded spurious shear stress-strain curves. In the numerical analysis, it is shown that the Iosipescu specimens for different fiber orientations have to be modeled differently in order to closely approximate the true loading conditions. Correction factors are needed to allow for the nonuniformity of the strain field and the use of the average shear stress in the shear modulus evaluation. The correction factors, which are determined for the region occupied by the strain gage rosette, are found to be dependent upon the material orthotropic ratio and the finite element models. Based upon the experimental and numerical results, recommendations for improving the reliability and accuracy of the shear modulus values are made, and the implications for shear strength measurement discussed. Further application of the Iosipescu shear test to woven fabric composites is presented. The limitations of the traditional strain gage instrumentation on the satin weave and high tow plain weave fabrics is discussed. Test results of a epoxy based aluminum particulate composite is also presented. A modification of the Iosipescu specimen is proposed and investigated experimentally and numerically. It is shown that the proposed new specimen design provides a more uniform shear stress field in the test section and greatly reduces the normal and shear stress concentrations in the vicinity of the notches. While the fabrication and the material cost of the proposed specimen is tremendously reduced, it is shown the accuracy of the shear modulus measurement is not sacrificed.
The dynamic response of a composite laminate plate to transverse impact loading is simulated numerically. An FEM scheme based on an eight-point brick element and three-dimensional linear elasticity theory is employed, assuming that the laminate layers are homogeneous and orthotropic. The derivation of the model is given in detail; the numerical implementation is briefly outlined; and results for the impact of 0.5-inch-diameter Al spheres at 500, 1000, and 1500 inch/sec on clamped or simply supported 3 x 3-inch plates made of T300/934 graphite/epoxy prepreg with ply orientation (0/-45/45/90)2S are presented in extensive graphs. The model is shown to provide accurate information on the histories of impactor-plate contact force; impactor displacement and velocity; and displacement, velocity, strain, and stress throughout the plate.
An analytical approach to determine an optimum laminate for a variety of thermal and mechanical loading combinations is presented. The analysis is performed for a linear elastic material under static mechanical and uniform thermal loadings. The problem is restricted to a unit width and length laminate with angle orientations resulting in an orthotropic, symmetric, and balanced configuration. An objective function defining total strain energy, is formulated and an optimum laminate design determined subject to constraints on stiffness, average coefficient of thermal expansion, and strength. The objective function is formulated in terms of the orientation angles, number of plies, and material properties. The method presented has, in varying degrees, shown that the design of a laminate can be accomplished using strain energy minimization as the primary criteria. The results of various combinations of applied constraints in the optimized design process are presented and discussed.
A numerical study of the postbuckling behavior of open-section stiffened composite compression panels is presented with emphasis on the effects of an anisotropic attached flange on results, the strain distribution near collapse, and the change of buckle pattern during postbuckling response. Results are obtained for a blade stiffened panel with orthotropic or anisotropic attached flanges from a new verison of STAGS, a general, branched shell, nonlinear computer program. Comparisons between these results indicate the effects of anisotropic flanges on the results. Elastic strain distributions are obtained from STAGS and indications are given about the collapse loads and modes of a particular panel. Change of buckle pattern is studied using STAGS and a special purpose computer program for long plates. Comparisons are made with experimental results.
A methodology is presented for the tension fatigue analysis and life prediction of composite laminates subjected to tension fatigue loading. The methodology incorporates both the generic fracture mechanics characterization of delamination and the assessment of the infuence of damage on laminate fatigue life. Tension fatigue tests were conducted on quasi-isotropic and orthotropic glass epoxy, graphite epoxy, and glass/graphite epoxy hybrid laminates, demonstrating good agreement between measured and predicted lives.
This paper deals with a dynamic stability analysis of unidirectional fiber-reinforced composite viscoelastic plates subjected to compressive edge loads. The integrodifferential equations governing the stability problem are obtained by using, in conjunction with a Boltzmann hereditary constitutive law for a three-dimensional viscoelastic medium, a higher-order shear deformation theory of orthotropic plates. Such a theory incorporates transverse shear deformation, transverse normal stress, and rotatory inertia effects. The solution of the stability problem as considered within this paper concerns the determination of the critical in-plane edge loads yielding the asymptotic instability. Numerical applications, based on material properties derived within the framework of Aboudi's micromechanical model, are presented and pertinent conclusions concerning the nature of the loss of stability and the influence of various parameters are outlined.
Room temperature compressive buckling strengths of eight cylinders, four boron-epoxy and four boron-epoxy reinforced-titanium, with diameter to thickness ratios ranging between 40 and 67 are determined experimentally and compared with analytical predictions. Numerical buckling strengths are presented for Donnell's, Flugge's and Sanders' shell theories for anisotropic and orthotropic material cases. Comparison of analytical predictions with experimental results indicates good agreement and the recommended correlation factor suggested in the literature is applicable for design. For the cylinders tested, the correlation between experiment and theory ranged from 0.73 to 0.97.
The vibration of an initially compressed plate having a circular geometry and orthotropy is examined in an affine space. The classical linear plate theory and Hamilton's principle are employed. The equations of motion of the plate are particularly simple in the chosen affine space, permitting a free-vibration study of the entire spectrum of composite materials with polar orthotropy. Approximate but very accurate standing-wave-type mode shapes are used in solving the essentially double eigenvalue problem to determine the effects of midplane forces on the vibration frequencies of the plate. The results indicate that the affine-space frequency increases with increasing stiffness ratio but decreases with increasing midplane compression. It is also discovered that, contrary to the trends observed for rectangular geometry and orthotropy by Oyibo (1981), Brunelle (1982), and Brunelle and Oyibo (1983), the affine-space frequency increases with increasing generalized Poisson's ratio.
Heat conduction throughout the blade and the distribution of thermal stresses caused by the temperature distribution were determined for a laminated wood wind turbine blade in both the horizontal and vertical positions. Results show that blade cracking is not due to thermal stresses induced by insulation. A method and practical example of thermal stress analysis for an engineering body of orthotropic materials is presented.
Predicted engineering stiffness parameters of filament-wound composite-material cylinders verified with respect to experimental data, by use of equations developed straightforwardly from applicable formulation of Hooke's law. Equations derived in engineering study of filament-wound rocket-motor cases, also applicable to other cylindrical pressure vessels made of orthotropic materials.
The viscoplasticity theory based on overstress, a simplified version of Sutcu (1985) theory, is presented for the case of orthotropy. The theory permits asymptotic solutions for checking numerical experiments and for identification of material properties. It is applicable to materials with rate(time) dependence, hysteresis, and a permanent set, such as cubic crystals, directionally solidified alloys, and metal matrix composites. Though it is a unified type theory, with effects of creep and plasticity not separated, it allows to identify viscous and plastic contributions to the stress in the asymptotic state.
The results of an experimental and analytical study of the effects of initial imperfections on the buckling and postbuckling response of unstiffened thin-walled compression-loaded graphite-epoxy cylindrical shells are presented. The shells considered in the study have four different shell-wall laminates and two different shell-radius-to-thickness ratios. The shell-wall laminates include two different orthotropic laminates and two different quasi-isotropic laminates. The shell-radius-to-thickness ratios include shell-radius-to-thickness ratios equal to 100 and 200. The results identify the effects of traditional and nontraditional initial imperfections on the nonlinear response characteristics and buckling loads of the shells. The traditional imperfections include the geometric shell-wall mid-surface imperfections that are commonly discussed in the literature on thin shell buckling. The nontraditional imperfections include shell-wall thickness variations, local shell-wall ply-gaps associated with the fabrication process, shell-end geometric imperfections, nonuniform applied end loads, and variations in the boundary conditions including the effects of elastic boundary conditions. A high-fidelity nonlinear shell analysis procedure that accurately accounts for the effects of these traditional and nontraditional imperfections on the nonlinear response characteristics and buckling loads of the shells is described. The analysis procedure includes a nonlinear static analysis that predicts the stable response characteristics of the shells, and a nonlinear transient analysis that predicts the unstable response characteristics. The results of a local shell-wall stress analysis used to estimate failure stresses are also described.