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.
Report presents study of mathematical modeling of stresses and strains, reaching beyond limits of elasticity, in bars and plates. Study oriented toward development of capability to predict stresses and resulting elastic and plastic strains in butt welds.
We present an estimate of the optical performance of the Thirty Meter Telescope (TMT) after execution of the full telescope alignment plan. The TMT alignment is performed by the Global Metrology System (GMS) and the Alignment and Phasing System (APS). The GMS first measures the locations of the telescope optics and instruments as a function of elevation angle. These initial measurements will be used to adjust the optics positions and build initial elevation look-up tables. Then the telescope is aligned using starlight as the input for the APS at multiple elevation angles. APS measurements are used to refine the telescope alignment to build elevation and temperature dependent look-up tables. Due to the number of degrees of freedom in the telescope (over 10,000), the ability of the primary mirror to correct aberrations on other optics, the tight optical performance requirements and the multiple instrument locations, it is challenging to develop, test and validate these alignment procedures. In this paper, we consider several GMS and APS operational scenarios. We apply the alignment procedures to the model-generated TMT, which consists of various quasi-static errors such as polishing errors, passive supports errors, thermal and gravity deformations and installation position errors. Using an integrated optical model and Monte-Carlo framework, we evaluate the TMT’s aligned states using optical performance metrics at multiple instrument and field of view locations. The optical performance metrics include the Normalized Point Source Sensitivity (PSSN), RMS wavefront error before and after Adaptive Optics (AO) correction, pupil position change, and plate scale distortion.
Poisson's ratio in materials is governed by the following aspects of the microstructure: the presence of rotational degrees of freedom, non-affine deformation kinematics, or anisotropic structure. Several structural models are examined. The non-affine kinematics are seen to be essential for the production of negative Poisson's ratios for isotropic materials containing central force linkages of positive stiffness. Non-central forces combined with pre-load can also give rise to a negative Poisson's ratio in isotropic materials. A chiral microstructure with non-central force interaction or non-affine deformation can also exhibit a negative Poisson's ratio. Toughness and damage resistance in these materials may be affected by the Poisson's ratio itself, as well as by generalized continuum aspects associated with the microstructure.
This report describes the computational steps involved in executing a finite-element-based perturbation model for computing the rotor dynamic coefficients of a shrouded pump impeller or a simple seal. These arise from the fluid/rotor interaction in the clearance gap. In addition to the sample cases, the computational procedure also applies to a separate category of problems referred to as the 'seal-like' category. The problem, in this case, concerns a shrouded impeller, with the exception that the secondary, or leakage, passage is totally isolated from the primary-flow passage. The difference between this and the pump problem is that the former is analytically of the simple 'seal-like' configuration, with two (inlet and exit) flow-permeable stations, while the latter constitutes a double-entry / double-discharge flow problem. In all cases, the problem is that of a rotor clearance gap. The problem here is that of a rotor excitation in the form of a cylindrical whirl around the housing centerline for a smooth annular seal. In its centered operation mode, the rotor is assumed to give rise to an axisymmetric flow field in the clearance gap. As a result, problems involving longitudinal or helical grooves, in the rotor or housing surfaces, go beyond the code capabilities. Discarding, for the moment, the pre- and post-processing phases, the bulk of the computational procedure consists of two main steps. The first is aimed at producing the axisymmetric 'zeroth-order' flow solution in the given flow domain. Detailed description of this problem, including the flow-governing equations, turbulence closure, boundary conditions, and the finite-element formulation, was covered by Baskharone and Hensel. The second main step is where the perturbation model is implemented, with the input being the centered-rotor 'zeroth-order' flow solution and a prescribed whirl frequency ratio (whirl frequency divided by the impeller speed). The computational domain, in the latter case, is treated as three dimensional, with the number of computational planes in the circumferential direction being specified a priori. The reader is reminded that the deformations in the finite elements are all infinitesimally small because the rotor eccentricity itself is a virtual displacement. This explains why we have generically termed the perturbation model the 'virtually' deformable finite-element category. The primary outcome of implementing the perturbation model is the tangential and radial components, F(sub theta)(sup *) and F(sub r)(sup *) of the fluid-exerted force on the rotor surface due to the whirling motion. Repetitive execution of the perturbation model subprogram over a sufficient range of whirl frequency ratios, and subsequent interpolation of these fluid forces, using the least-square method, finally enable the user to compute the impeller rotor dynamic coefficients of the fluid/rotor interaction. These are the direct and cross-coupled stiffness, damping, and inertia effects of the fluid/rotor interaction.
A computer model has been constructed to simulate the compliance and load sharing in a spur gear mesh. The model adds the effect of rim deflections to previously developed state-of-the-art gear tooth deflection models. The effects of deflections on mesh compliance and load sharing are examined. The model can treat gear meshes composed to two external gears or an external gear driving an internal gear. The model includes deflection contributions from the bending and shear in the teeth, the Hertzian contact deformations, and primary and secondary rotations of the gear rims. The model shows that rimmed gears increase mesh compliance and, in some cases, improve load sharing.
A computer model has been constructed to simulate the compliance and load sharing in a spur gear mesh. The model adds the effect of rim deflections to previously developed state-of-the-art gear tooth deflection models. The effects of deflections on mesh compliance and load sharing are examined. The model can treat gear meshes composed of two external gears or an external gear driving an internal gear. The model includes deflection contributions from the bending and shear in the teeth, the Hertzian contact deformations, and primary and secondary rotations of the gear rims. The model shows that rimmed gears increase mesh compliance and, in some cases, improve load sharing.
Lightweight aerospace structures, such as low areal density composite space reflectors, are highly flexible and may undergo large deflection under applied loading, especially during the launch phase. Accordingly, geometrically nonlinear analysis that takes into account the effect of finite rotation may be needed to determine the deformed shape for a clearance check and the stress and strain state to ensure structural integrity. In this study, deformation of the space reflector is determined under static conditions using a geometrically nonlinear solid shell finite element model. For the solid shell element formulation, the kinematics of deformation is described by six variables that are purely vector components. Because rotational angles are not used, this approach is free of the limitations of small angle increments. This also allows easy connections between substructures and large load increments with respect to the conventional shell formulation using rotational parameters. Geometrically nonlinear analyses were carried out for three cases of static point loads applied at selected points. A chart shows results for a case when the load is applied at the center point of the reflector dish. The computed results capture the nonlinear behavior of the composite reflector as the applied load increases. Also, they are in good agreement with the data obtained by experiments.
Thermoplastic Composites (TPCs) are being increasingly considered for aerospace applications given their faster manufacturing cycles and lower cost. The NASA High-Rate Composite Aircraft Manufacturing (HiCAM) project aims to evaluate and mature thermoplastic manufacturing technologies to achieve a 4X – 6X increase in the production rate of Next Generation Single Aisle (NGSA) commercial aircraft. Within the purview of the NASA HiCAM project, this work investigates the response of thermoplastic composites under High Energy Dynamic Impact (HEDI) conditions. HEDI tests were conducted on panels fabricated with a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) material system. A blunt metal projectile was chosen as the impactor which impacted the T700/LMPAEK panels over a range of impact velocities imparted by a gas gun test setup. As expected, lower velocity impacts caused the projectile to rebound, whereas higher velocity impacts resulted in the projectile penetrating the test panels. For the test cases in which the projectile rebounded, the predominant damage modes recorded with Ultrasonic (UT) scans included interlaminar delaminations which exhibited a “rotating fan” like structure when viewed in the through-thickness direction. Additionally, the dynamic deflection of the center point of the back face of the test panels and the velocity of the projectile was also measured during the test. Finite element models were developed to model a test case in which the projectile rebounded. These models were developed using LS-DYNA® wherein the primary objective was to evaluate the ability of the material model MAT299 to capture the dynamic response and damage modes in the T700/LMPAEK test panels. MAT299 is a Deformation Gradient Decomposition (DGD) based Continuum Damage Mechanics (CDM) material model. The interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA, wherein the mixed-mode fracture is captured via the Benzegaggh-Kenane (B-K) mode-dependent fracture energy interpolation law. The above-mentioned modeling methods have mostly been applied to thermoset composites and adapting them to TPCs involves addressing challenges associated with appropriately representing material behavior. Therefore, the proposed paper would discuss the systematic approach taken to adapt thermoset composites modeling practices for applications to TPCs, while appropriately addressing material behavior. The proposed paper would show that the predicted back face center point deflection correlated reasonably well with experimental data and the peak deflection was predicted to be within 5% of the experimental measurements. The projectile rebound velocity, while predicted to be higher than the experimental measurement, was within reasonable bounds. Additionally, it will be shown that the predicted delamination shapes were in reasonable agreement with experimental data.
The mountain belts of Ishtar Terra are unique on Venus. Models for their formation include mantle upwelling, mantle downwelling, and horizontal convergence. The present forms of these models are too simple to predict surface strain, topography, or gravity. More detailed models will require specific constraints as imposed by geologic relations. In order to develop specific constraints for geodynamic models, we examine the geology of Ishtar Terra as viewed in Magellan SAR imagery in an attempt to interpret regional surface strain patterns. In this paper, we present geologic and structural relations that leads us to postulate that Ishtar deformed belts result from shear forces within the mantle acting on the lithosphere, and not by horizontal forces from colliding plates. We propose that the surface strains result from differential strain and displacement of domains within the upper mantle, and that further analysis of Ishtar deformation may allow us to identify individual domains within the mantle, and to constrain displacement trajectories between domains.
The reusable launch vehicle (RLV) thrust cell liner, or thrust chamber, is a critical component of the Space Shuttle Main Engine (SSME). It is designed to operate in some of the most severe conditions seen in engineering practice. This requirement, in conjunction with experimentally observed 'dog-house' failure modes characterized by bulging and thinning of the cooling channel wall, provides the motivation to study the factors that influence RLV thrust cell liner performance. Factors or parameters believed to be directly related to the observed characteristic deformation modes leading to failure under in-service loading conditions are identified, and subsequently investigated using the cylindrical version of the higher-order theory for functionally graded materials in conjunction with the Robinson's unified viscoplasticity theory and the power-law creep model for modeling the response of the liner s constituents. Configurations are analyzed in which specific modifications in cooling channel wall thickness or constituent materials are made to determine the influence of these parameters on the deformations resulting in the observed failure modes in the outer walls of the cooling channel. The application of thermal barrier coatings and functional grading are also investigated within this context. Comparison of the higher-order theory results based on the Robinson and power-law creep model predictions has demonstrated that, using the available material parameters, the power-law creep model predicts more precisely the experimentally observed deformation leading to the 'dog-house' failure mode for multiple short cycles, while also providing much improved computational efficiency. However, for a single long cycle, both models predict virtually identical deformations. Increasing the power-law creep model coefficients produces appreciable deformations after just one long cycle that would normally be obtained after multiple cycles, thereby enhancing the efficiency of the analysis. This provides a basis for the development of an accelerated modeling procedure to further characterize dog-house deformation modes in RLV thrust cell liners. Additionally, the results presented herein have demonstrated that the mechanism responsible for deformation leading to 'dog-house' failure modes is driven by pressure, creep/relaxation and geometric effects.
The rate-of-slip vector on the San Andreas fault in central California estimated from geodetic and Holocene geological data is inconsistent with the prediction of rigid plate models such as RM2. This well-known 'San Andreas discrepancy' is diagnostic of plate deformation distributed both east of the fault in the Basin and Range and west of the fault along the California continental margin. Constraints on the integrated deformation rates across these two regions consistent with: (1) the kinematical boundary conditions imposed by the rigid plate model; (2) neotectonic and paleoseismic estimates of deformation rates; (3) ground-based geodetic measurements; and (4) rates of change observed by very long baseline interferometry along seven baselines to western U.S. sites are constructed. The space-geodetic data on Basin and Range extension taken over a 4-year interval are compatible with geological observations averaged over the Holocene; the best estimate of its integrated deformation rate, provided by the joint inversion of both data types, is 9.7 + or - 2.1 mm/yr, N 56 deg W + or - 10 deg, too small and in the wrong direction to account entirely for the San Andreas discrepancy. The integral of this deformation, estimated by subtracting the Basin and Range contribution from the discrepancy vector, requires significant right-lateral shear parallel to the San Andreas (13 + or - 5 mm/yr) and some compression perpendicular to it (9 + or - 3 mm/yr).
This paper will discuss the objectives, approach, and anticipated results of a study of earthquakes in Costa Rica. GPS measurements will be taken and field surveys will be made. Assessments of seismic strain accumulation and post-seismic deformation will be made in an effort to understand the effect these processes have on regional tectonic models.
Spin forming is an advanced manufacturing process widely used in the aerospace and defense sectors to produce lightweight, high-strength cylindrical components with tight dimensional tolerances. This study explores the applicability of the path-dependent Mechanical Threshold Stress (MTS) constitutive model by simulating the evolution of geometry, machining forces, and plastic deformation during the spin forming of a 10-mm thick 6061-O aluminum cylinder. While numerical modeling of spin forming has advanced substantially over the past decade, systematic verification and experimental validation of material models remain limited, particularly in predicting through-thickness process evolution. The MTS model, incorporating a Voce hardening rule, is employed for its ability to represent cyclic loading, rapidly varying temperature fields, and strain rates characteristic of spin forming. Numerical convergence analysis indicates discretization uncertainties between 0.3% and 9.2% for key quantities of interest. Experimental validation demonstrates that the MTS model, when implemented with a verified mesh, accurately reproduces both elastic and plastic behavior of 6061-O aluminum, predicting peak roller loads within 11–18% of measurements, geometric tolerances within 3%, and plastic strain distributions within 10% of experimental values. Collectively, these results establish a validated computational framework for predictive spin-forming simulations with quantified confidence, providing a foundation for extension to other alloys, geometries, and forming conditions.
A direct method for determining the effective elastic stiffnesses and deformation behavior of composite box-beam (BB) structures is developed analytically, validated, and demonstrated. The BB walls are modeled as orthotropic-ply laminated plates, so that the elastic properties vary both through the thickness and around the BB contour; deformation is described in terms of extension, bending, twisting, shearing, and torsion-related out-of-plane warping. Numerical results for several BB configurations are presented in tables and graphs and compared with experimental data and FEM computations (Stemple and Lee, 1989): good general agreement is obtained for cross-ply, antisymmetric, and symmetric layups, except for symmetric layups with ply orientation angles theta greater than 30 deg. Both out-of-plane warping and transverse shear coupling are found to have significant effects on BB elastic response.
Simple models of the Martian crust are summarized that predict extensional deformation style on the basis of depth, material friction and strength, and hydraulic conditions appropriate to the planet. These models indicate that tension fractures may be common features on Mars, given adequate differential stress conditions. Examples of tension fractures on Mars inferred from morphological criteria are examined based on the probable geologic conditions in which they formed and on model constraints. It is proposed that the grabens and collapse features of Valles Marineris are controlled by tension fractures in intact basement rocks that lie below impact ejecta.
Increasingly, carbon composite structures are being used in aerospace applications. Their highstrength, high-stiffness, and low-weight properties make them good candidates for replacing many aerospace structures currently made of aluminum or steel. Recently, many of the aircraft engine manufacturers have developed new commercial jet engines that will use composite fan cases. Instead of using traditional composite layup techniques, these new fan cases will use a triaxially braided pattern, which improves case performance. The impact characteristics of composite materials for jet engine fan case applications have been an important research topic because Federal regulations require that an engine case be able to contain a blade and blade fragments during an engine blade-out event. Once the impact characteristics of these triaxial braided materials become known, computer models can be developed to simulate a jet engine blade-out event, thus reducing cost and time in the development of these composite jet engine cases. The two main problems that have arisen in this area of research are that the properties for these materials have not been fully determined and computationally efficient computer models, which incorporate much of the microscale deformation and failure mechanisms, are not available. The research reported herein addresses some of the deficiencies present in previous research regarding these triaxial braided composite materials. The current research develops new techniques to accurately quantify the material properties of the triaxial braided composite materials. New test methods are developed for the polymer resin composite constituent and representative composite coupons. These methods expand previous research by using novel specimen designs along with using a noncontact measuring system that is also capable of identifying and quantifying many of the microscale failure mechanisms present in the materials. Finally, using the data gathered, a new hybrid micromacromechanical computer model is created to simulate the behavior of these composite material systems under static and ballistic impact loading using the test data acquired. The model also quantifies the way in which the fiber/matrix interface affects material response under static and impact loading. The results show that the test methods are capable of accurately quantifying the polymer resin under a variety of strain rates and temperature for three loading conditions. The resin strength and stiffness data show a clear rate and temperature dependence. The data also show the hydrostatic stress effects and hysteresis, all of which can be used by researchers developing composite constitutive models for the resins. The results for the composite data reveal noticeable differences in strength, failure strain, and stiffness in the different material systems presented. The investigations into the microscale failure mechanisms provide information about the nature of the different material system behaviors. Finally, the developed computer model predicts composite static strength and stiffness to within 10 percent of the gathered test data and also agrees with composite impact data, where available.
The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.