Analysis of thermally induced structural vibrations by finite element techniques
Analysis of thermally induced structural vibrations by finite element techniques
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.
Analysis of thermally induced structural vibrations by finite element techniques
A finite element computer program (TAP 2) for steady-state and transient thermal analyses of convectively cooled structures is presented. The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature-dependent thermal parameters is performed using the Newton-Raphson iteration method. Transient analyses are performed using an implicit Crank-Nicolson time integration scheme with consistent or lumped capacitance matrices as an option. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. User instructions and sample problems are presented in appendixes.
The finite element ablation and thermal response (FEAtR, hence forth called FEAR) design and analysis program simulates the one, two, or three-dimensional ablation, internal heat conduction, thermal decomposition, and pyrolysis gas flow of thermal protection system materials. As part of a code validation study, two-dimensional axisymmetric results from FEAR are compared to thermal response data obtained from an arc-jet stagnation test in this paper. The results from FEAR are also compared to the two-dimensional axisymmetric computations from the two-dimensional implicit thermal response and ablation program under the same arcjet conditions. The ablating material being used in this arcjet test is phenolic impregnated carbon ablator with an LI-2200 insulator as backup material. The test is performed at the NASA, Ames Research Center Interaction Heating Facility. Spatially distributed computational fluid dynamics solutions for the flow field around the test article are used for the surface boundary conditions.
A new finite element is described for performing detailed thermal stress analysis of thin orthotropic shells of revolution. The element provides for temperature loadings which may vary over the surface of the shell as well as through the thickness. In a number of sample calculations, results from the present method are compared with analytical solutions as well as with independent numerical analyses. Such calculations are carried out for two cylinders, a conical frustum, a truncated hemisphere, and an annular plate. Generally, the agreement between the present solution and the other solutions is excellent.
A hierarchical finite element approach for thermal-structural analysis is presented. The approach employs a common nodal discritization and seeks improvements in the accuracy of the analyses by using hierarchical interpolation functions with nodeless variables. The effectiveness of the integrated approach is assessed for three applications with two-dimensional elements by comparison with conventional finite element thermal-structural solutions. Improvements in the accuracy of temperatures and thermal-stresses are demonstrated. The applications demonstrate the practical importance of having flexibility in refining each analysis independently while maintaining a common discretization, and show that the hierarchical approach offers potential for the development of a general method for integrated thermal-structural analysis.
A hierarchical flux-based finite element method is developed for both a one and two dimensional thermal structural analyses. Derivation of the finite element equations is presented. The resulting finite element matrices associated with the flux based formulation are evaluated in a closed form. The hierarchical finite elements include additional degrees of freedom in the approximation of the element variable distributions by the use of nodeless variables. The nodeless variables offer increased solution accuracy without the need for defining actual nodes and rediscretizing the finite element model. Thermal and structural responses are obtained from a conventional linear finite element method and exact solutions. Results show that the hierarchical flux-based method can provide improved thermal and structural solution accuracy with fewer elements when compared to results for the conventional linear element method.
The near-term objective is to develop an understanding of the states of stresses and strains in a Zirconia-yttria thermal barrier coating (TBC) experiencing a given temperature drop. Results so obtained are expected to facilitate experimental work. In order to gain realistic insights into the distribution of stresses and strains in a complex TBC, the finite element approach was selected to model a cylindrical TBC specimen. Experimental evidence reported in the literature indicated the presence of rough interface between the ceramic coat and bond coat. Oxidation of the bond coat at ceramic-bond interface was observed, as was a small amount of cracking in the ceramics near the ceramic-bond interface. To account for these complex features, a plane-strain finite element computer program known as TBCOC was developed, taking advantage of a generic computer code known as MARC. This generic code was made available through the use of a supercomputer (Cray 1). The TBCOC model contains 1316 nodal points and 2140 finite elements. It is capable of a uniform isothermal loading. Results of a sample computer run are presented.
Finite-element, transient heat transfer analyses were performed for the first-stage blades of the space shuttle main engine (SSME) high-pressure fuel turbopump. The analyses were based on test engine data provided by Rocketdyne. Heat transfer coefficients were predicted by performing a boundary-layer analysis at steady-state conditions with the STAN5 boundary-layer code. Two different peak-temperature overshoots were evaluated for the startup transient. Cutoff transient conditions were also analyzed. A reduced gas temperature profile based on actual thermocouple data was also considered. Transient heat transfer analyses were conducted with the MARC finite-element computer code.
Explore the source record for details and available documents.
The purpose of this research was to perform sensitivity studies and develop a process to perform thermal and structural analysis and sizing of the latest Metallic Thermal Protection System (TPS) developed at NASA LaRC (Langley Research Center). Metallic TPS is a key technology for reducing the cost of reusable launch vehicles (RLV), offering the combination of increased durability and competitive weights when compared to other systems. Accurate sizing of metallic TPS requires combined thermal and structural analysis. Initial sensitivity studies were conducted using transient one-dimensional finite element thermal analysis to determine the influence of various TPS and analysis parameters on TPS weight. The thermal analysis model was then used in combination with static deflection and failure mode analysis of the sandwich panel outer surface of the TPS to obtain minimum weight TPS configurations at three vehicle stations on the windward centerline of a representative RLV. The coupled nature of the analysis requires an iterative analysis process, which will be described herein. Findings from the sensitivity analysis are reported, along with TPS designs at the three RLV vehicle stations considered.
The application of adaptive finite element methods to the solution of transient heat conduction problems in two dimensions is investigated. The computational domain is represented by an unstructured assembly of linear triangular elements and the mesh adaptation is achieved by local regeneration of the grid, using an error estimation procedure coupled to an automatic triangular mesh generator. Two alternative solution procedures are considered. In the first procedure, the solution is advanced by explicit timestepping, with domain decomposition being used to improve the computational efficiency of the method. In the second procedure, an algorithm for constructing continuous lines which pass only once through each node of the mesh is employed. The lines are used as the basis of a fully implicit method, in which the equation system is solved by line relaxation using a block tridiagonal equation solver. The numerical performance of the two procedures is compared for the analysis of a problem involving a moving heat source applied to a convectively cooled cylindrical leading edge.
A structural performance and resizing (SPAR) finite-element thermal analysis computer program was used in the heat-transfer analysis of the space shuttle orbiter subjected to reentry aerodynamic heating. Three wing cross sections and one midfuselage cross section were selected for the thermal analysis. The predicted thermal protection system temperatures were found to agree well with flight-measured temperatures. The calculated aluminum structural temperatures also agreed reasonably well with the flight data from reentry to touchdown. The effects of internal radiation and of internal convection were found to be significant. The SPAR finite-element solutions agreed reasonably well with those obtained from the conventional finite-difference method.
A structural performance and resizing finite element thermal analysis computer program was used in the reentry heat transfer analysis of the space shuttle. Two typical wing cross sections and a midfuselage cross section were selected for the analysis. The surface heat inputs to the thermal models were obtained from aerodynamic heating analyses, which assumed a purely turbulent boundary layer, a purely laminar boundary layer, separated flow, and transition from laminar to turbulent flow. The effect of internal radiation was found to be quite significant. With the effect of the internal radiation considered, the wing lower skin temperature became about 39 C (70 F) lower. The results were compared with fight data for space transportation system, trajectory 1. The calculated and measured temperatures compared well for the wing if laminar flow was assumed for the lower surface and bay one upper surface and if separated flow was assumed for the upper surfaces of bays other than bay one. For the fuselage, good agreement between the calculated and measured data was obtained if laminar flow was assumed for the bottom surface. The structural temperatures were found to reach their peak values shortly before touchdown. In addition, the finite element solutions were compared with those obtained from the conventional finite difference solutions.
Papers entitled, An Adaptive Finite Element Procedure for Compressible Flows and Strong Viscous-Inviscid Interactions, and An Adaptive Remeshing Method for Finite Element Thermal Analysis, were presented at the June 27 to 29, 1988, meeting of the AIAA Thermophysics, Plasma Dynamics and Lasers Conference, San Antonio, Texas. The papers describe research work supported under NASA/Langley Research Grant NsG-1321, and are submitted in fulfillment of the progress report requirement on the grant for the period ending February 29, 1988.
An integrated thermal-structural finite element approach for efficient coupling of thermal and structural analysis is presented. New thermal finite elements which yield exact nodal and element temperatures for one dimensional linear steady state heat transfer problems are developed. A nodeless variable formulation is used to establish improved thermal finite elements for one dimensional nonlinear transient and two dimensional linear transient heat transfer problems. The thermal finite elements provide detailed temperature distributions without using additional element nodes and permit a common discretization with lower order congruent structural finite elements. The accuracy of the integrated approach is evaluated by comparisons with analytical solutions and conventional finite element thermal structural analyses for a number of academic and more realistic problems. Results indicate that the approach provides a significant improvement in the accuracy and efficiency of thermal stress analysis for structures with complex temperature distributions.
Preflight predictions of the structural temperature distributions during entry are compared with data from the initial Shuttle flight. Finite element thermal analysis programming was used to model the heat flow on Shuttle structures and actual gas properties of air were employed in the analyses of aerodynamic heating. Laminar, separated, and turbulent heat fluxes were calculated for varying locations on the craft using velocity-attitude and angle-of-attack projections taken from the nominal STS-1 trajectory. Temperature time histories of the first flight are compared with laminar and turbulent flow assumptions and an unpredicted rapid cooling 1800 sec into entry is credited to inaccurate assumptions of structural heat dissipative properties or flow conditions in that time phase of the flight; additional discrepancies in descriptions of heating of the upper fuselage are attributed to a lack of knowledge of the complex flow patterns existing over that area of the Shuttle body.
NASA is studying the next generation of space suit for safe operations in space, from the International Space Station to the lunar surface as part of the Exploration Extravehicular Mobility Unit project. The backplate structure of the portable life support system (PLSS) serves as the attachment point for assembly of the PLSS and contains embedded cooling channels. The thermal loop plumbing assembly and backplate are combined into a single component such that the fluid lines run through the midplane of the backplate. The manufacturing of the Ti-6Al-4V backplate consists of machining channels into a plate, welding cover plates on each channel using an electron beam, and additional machining for excess material removal and final detail. Careful control of the electron beam weld process is necessary to completely seal cover plates, produce consistent weld tracks, and produce channels with uniform cross sections free of debris and defects. The backplate contains multiple weld channels with unique geometries. The weld process parameters include choosing the beam current and focus, accelerating voltage, travel speed, pause times between weld tracks, and sequence of weld tracks. A transient finite element thermal analysis was developed to provide support for choosing and controlling the weld process parameters to maintain a consistent melt pool along each of the weld tracks. A moving distributed heat source is modeled to represent the heat input generated by the electron beam. The heat source model is a critical input to the analysis for accurate temperature predictions during the weld process. A modified conical heat source was developed to account for the change in working distance between the backplate and cover plates. The geometry of the material being welded must be accounted for by the heat source model to define the correct energy input. The weld process control is designed according to the predicted thermal history and melt pool dimensions for each of the weld channels. The analysis model, modified heat source, temperature and melt pool predictions, application of the thermal data to the weld parameters, and comparison with micrograph cross-sections from weld trials will be discussed in this presentation.
Future space exploration missions will require the development of more advanced in-space radiators. These radiators should be highly efficient and lightweight, deployable heat rejection systems. Typical radiators for in-space heat mitigation commonly comprise a substantial portion of the total vehicle mass. A small mass savings of even 5-10% can greatly improve vehicle performance. The objective of this paper is to present the development of detailed tools for the analysis and design of in-space radiators using evolutionary computation techniques. The optimality criterion is defined as a two-dimensional radiator with a shape demonstrating the smallest mass for the greatest overall heat transfer, thus the end result is a set of highly functional radiator designs. This cross-disciplinary work combines topology optimization and thermal analysis design by means of a genetic algorithm The proposed design tool consists of the following steps; design parameterization based on the exterior boundary of the radiator, objective function definition (mass minimization and heat loss maximization), objective function evaluation via finite element analysis (thermal radiation analysis) and optimization based on evolutionary algorithms. The radiator design problem is defined as follows: the input force is a driving temperature and the output reaction is heat loss. Appropriate modeling of the space environment is added to capture its effect on the radiator. The design parameters chosen for this radiator shape optimization problem fall into two classes, variable height along the width of the radiator and a spline curve defining the -material boundary of the radiator. The implementation of multiple design parameter schemes allows the user to have more confidence in the radiator optimization tool upon demonstration of convergence between the two design parameter schemes. This tool easily allows the user to manipulate the driving temperature regions thus permitting detailed design of in-space radiators for unique situations. Preliminary results indicate an optimized shape following that of the temperature distribution regions in the "cooler" portions of the radiator. The results closely follow the expected radiator shape.