Investigation of prediction methods for the loads and stresses of Apollo type spacecraft parachutes. Volume 2: Stresses
For abstract, see N74-19673.
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.
For abstract, see N74-19673.
A digital computer program has been developed to calculate unsteady loadings caused by motions of lifting surfaces with leading edge or trailing edge controls based on the subsonic kernel function approach. The pressure singularities at hinge line and side edges have been extracted analytically as a preliminary step to solving the integral equation by collocation. The program calculates generalized aerodynamic forces for user supplied deflection modes. Optional intermediate output includes pressure at an array of points, and sectional generalized forces. From one to six controls on the half span can be accommodated.
The vortex-lattice technique for incompressible flow which accounts for separation at sharp edges is modified to account for compressibility. This is accomplished by extending the Prandtl-Glauert transformation to moderate angles of attack. Thus, the aerodynamic characteristics for the compressible case are obtained from the solution of an equivalent incompressible problem. Numerical results are presented for parallelogram and delta wings to assess the effects of compressibility. The results are in good agreement with available experimental data.
Explore the source record for details and available documents.
A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. Using the linearized Euler technique, one decomposes the flow into a mean (or steady) flow plus an unsteady, harmonically varying, small disturbance flow. Linear variable coefficient equations describe the small disturbance behavior of the flow, and are solved using a pseudo-time marching Lax-Wendroff scheme. For the blade motion problem, a harmonically deforming computational rid that conforms to the motion of vibrating blades eliminates large error producing mean flow gradient terms that would otherwise appear in the unsteady flow tangency boundary condition. The paper also presents a new, numerically exact, nonreflecting far-field boundary condition based on an eigenanalysis of the discretized equations. Computed flow solutions demonstrate the computational accuracy and efficiency of the present method. The solution of the linearized Euler equations requires one to two orders of magnitude less computer time than solution of the nonlinear Euler equations using traditional time-accurate time-marching techniques. In addition, the deformable grid significantly improves the accuracy of the solution.
Explore the source record for details and available documents.
A method is proposed for assessing the forces that can delaminate composite laminates and bonded structures. The approach employs a higher-order plate theory to represent portions of the laminate. The ply properties are integrated through the thickness in a manner that is consistent with the displacement assumptions in the plate theory, thus providing for accurate representation of the stacking sequence effects. In addition to stress analysis, the method also provides the strain energy release rate for the growth of existing delaminations. The analysis has been incorporated into a computer code, TTSS (through thickness stretching and shear). Application examples, including a double cantilever beam and a curved laminate, are examined.
Geophysical models of flexural stresses in an elastic lithosphere due to an axisymmetric surface load typically predict a transition with increased distance from the center of the load of radial thrust faults to strike-slip faults to concentric normal faults. These model predictions are in conflict with the absence of annular zones of strike-slip faults around prominent loads such as lunar maria, Martian volcanoes, and the Martian Tharsis rise. We suggest that this paradox arises from difficulties in relating failure criteria for brittle rocks to the stress models. Indications that model stresses are inappropriate for use in fault-type prediction include (1) tensile principal stresses larger than realistic values of rock tensile strength, and/or (2) stress differences significantly larger than those allowed by rock-strength criteria. Predictions of surface faulting that are consistent with observations can be obtained instead by using tensile and shear failure criteria, along with calculated stress differences and trajectories, with model stress states not greatly in excess of the maximum allowed by rock fracture criteria.
A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.
Predicting whether or not vehicle batteries contain sufficient charge to support operations over the remainder of a given flight plan is critical for electric aircraft. This paper describes an approach for identifying upper and lower uncertainty bounds on predictions that aircraft batteries will continue to meet output power and voltage requirements over the remainder of a flight plan. Battery discharge prediction is considered here in terms of the following components; (i) online battery state of charge estimation; (ii) prediction of future battery power demand as a function of an aircraft flight plan; (iii) online estimation of additional parasitic battery loads; and finally, (iv) estimation of flight plan safety. Substantial uncertainty is considered to be an irremovable part of the battery discharge prediction problem. However, high-confidence estimates of flight plan safety or lack of safety are shown to be generated from even highly uncertain prognostic predictions.
The objective of this study is to develop an anlysis to predict buckling loads of ring-stiffened anisotropic cylinders subject to axial compression, torsion, and internal pressure. This structure is modeled as a branched shell. A nonlinear axisymmetric prebuckling equilibrium state is assumed which is amenable to an exact solution within each branch. Axisymmetric geometric imperfections are included. Buckling displacements are represented by a Fourier series in the circumferential coordination and the finite-element method in the axial coordinate. Application of the Trefftz criterion to the second variation of the total potential energy leads to a nonlinear eigenvalue problem for the buckling load and mode. Results are presented for both unstiffened and ring-stiffened cylinders in the form of buckling interaction diagrams. Imperfections can cause an unexpected buckling mode in the ring web which would not occur for the perfect structure, and pressurization diminishes the benefit of adding rings to the unstiffened shell to increase the buckling load. The implementation of the analysis methodology into a structural sizing algorithm is discussed.
Dynamic blade loads, including aerodynamic, gravitational, and inertial effects, are presented for two large horizontal-axis wind turbines: the ERDA-NASA 100 kW Mod-0 and 200 kw Mod-0A wind power systems. Calculated and measured loads are compared for an experimental Mod-0 machine in operation. Predicted blade loads are also given for the higher power Mod-0A wind turbine now being assembled for operation as part of a municipal power plant. Two major structural modifications have been made to the Mod-0 wind turbine for the purpose of reducing blade loads. A stairway within the truss tower was removed to reduce the impulsive aerodynamic loading caused by the tower wake on the downwind rotor blades. Also, the torsional stiffness of the yaw drive mechanism connecting the turbine nacelle to the tower was doubled to reduce rotor-tower interaction loads. Measured reductions in load obtained by means of these two modifications equaled or exceeded predictions.
The problem of predicting aerodynamic loads on the insulating tiles of the space shuttle thermal protection system (TPS) is discussed and seen to require a method for predicting pressure and mass flux in the gaps between tiles. A mathematical model of the tile-gap flow is developed, based upon a slow viscous (Stokes) flow analysis, and is verified against experimental data. The tile-gap pressure field is derived from a solution of the two-dimensional Laplace equation; the mass-flux vector is then calculated from the pressure gradient. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The flow model shows tile-gap mass flux to be very sensitive to the gap width indicating a need for coupling the TPS flow and tile displacement calculation. Analytical and experimental work to improve TPS flow predictions and a possible shuttle TPS hardware modification are recommended.
Robust gear designs consider not only crack initiation, but crack propagation trajectories for a fail-safe design. In actual gear operation, the magnitude as well as the position of the force changes as the gear rotates through the mesh. A study to determine the effect of moving gear tooth load on crack propagation predictions was performed. Two-dimensional analysis of an involute spur gear and three-dimensional analysis of a spiral-bevel pinion gear using the finite element method and boundary element method were studied and compared to experiments. A modified theory for predicting gear crack propagation paths based on the criteria of Erdogan and Sih was investigated. Crack simulation based on calculated stress intensity factors and mixed mode crack angle prediction techniques using a simple static analysis in which the tooth load was located at the highest point of single tooth contact was validated. For three-dimensional analysis, however, the analysis was valid only as long as the crack did not approach the contact region on the tooth.
Described here is the procedure for utilizing a damage constitutive model to predict progressive damage growth in laminated composites. In this model, the effects of the internal damage are represented by strain-like second order tensorial damage variables and enter the analysis through damage dependent ply level and laminate level constitutive equations. The growth of matrix cracks due to fatigue loading is predicted by an experimentally based damage evolutionary relationship. This model is incorporated into a computer code called FLAMSTR. This code is capable of predicting the constitutive response and matrix crack damage accumulation in fatigue loaded laminated composites. The structure and usage of FLAMSTR are presented along with sample input and output files to assist the code user. As an example problem, an analysis of crossply laminates subjected to two stage fatigue loading was conducted and the resulting damage accumulation and stress redistribution were examined to determine the effect of variations in fatigue load amplitude applied during the first stage of the load history. It was found that the model predicts a significant loading history effect on damage evolution.
The effect of rotor blade dynamics on aerodynamic and structural loads is examined for a conventional, main- rotor helicopter using both a comprehensive rotorcraft analysis (CAMRAD) and night test data. The impact of blade dynamics on blade section lift-coefficient time histories is studied by comparing predictions from both a rigid blade analysis and an elastic blade analysis with helicopter flight test data. The elastic blade analysis better predicts high-frequency behavior of section lift. In addition, components of the blade angle of attack, such as elastic blade twist, blade nap rate, blade slope velocity, and inflow, are examined as a function of blade mode. Elastic blade motion affects the blade angle of attack by a few tenths of a degree, and up to the sixth rotor harmonic. A similar study of the influence of blade dynamics on bending and torsion moments was also conducted. The modal analysis of the predicted blade structural loads suggested that five elastic bending deg of freedom (four flap and one lag) and three elastic torsion deg of freedom contributed to calculations of the blade structural loads. However, when structural bending load predictions from several elastic blade analyses were compared with flight test data, an elastic blade model consisting of only three elastic bending modes (first and second flap, and first lag), and two elastic torsion modes was found to be sufficient for maximum correlation.
Buckling loads of thick-walled, orthotropic, simply-supported right circular cylinders are predicted using a new higher-order transverse shear deformation theory. The higher-order theory shows that, by more accurately accounting for transverse shear deformation effects, the predicted buckling load may be reduced by as much as 80 percent compared to predictions based on conventional transverse shear deformation theory. A parametric study of the effect of ply orientation on the buckling load of axially compressed cylinders indicates that laminates containing 0 deg plies are the most sensitive to transverse shear deformation effects. Interaction curves for buckling of cylinders with axial compressive and external pressure loadings indicate that buckling loads due to external pressure loadings are much less sensitive to transverse shear deformation effects than those due to axial compressive loadings.
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.