Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computations on changing geometries and structures”

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 19 records

On the Sensitivity of Piezoceramics and Piezopolymers in Structural Integrity Monitoring of Large Trusses

An analytical assessment has been made of the reliability of using integrated microactuators and sensors in the form of piezoceramics and piezopolymers as joint integrity monitors in trussed systems. The concept is first implemented for a simple structure which consists of two truss members with a 45 deg lift angle joined at the apex. A piezoceramic patch (or piezopolymer film) bonded on the surface of one of the members at a location near the joint is used as a collocated actuator/sensor. The overall structural dynamic response under an excitation was modeled by finite element method. Different degrees of nodal constraints at the joints representing various degrees of joint integrity are employed. The resulting dynamic response showed distinct responses for varying joint stiffnesses. Parallel experimental work on a truss model using a multichannel data acquisition system and a digital signal analyzer confirms the results from analysis. We further studied the sensitivity of the micro-sensors to the behavior of joints of large arch truss structure. Results obtained for large trusses with many degrees of freedom indicate optimum locations of sensors for which the dynamic response signatures are distinct and distinguishable for relatively small changes in joint integrity and/or structural geometry. Computations based on finite element modeling show that locating the single actuator/sensor at the joint corresponding to the first loss of static stability appear optimal. Hence, static stability analysis of complex trusses can give us a good indication of the optimum placement of sensors for maximum response. This observation is important if few distributed sensors and actuators are available for placement in constructed facilities made from large trusses with many degrees of freedom. As an extension of this work a dynamic response signature identification technique to monitor in-service degradation of joints is under development for application to the monitoring of the integrity of adhesive joints in composite structures.

Abatan, A. O.↗

Numerical Simulation of Inlet Bleed with Circular Holes on Plate Under Shock-Wave/Boundary-Layer Interactions

A numerical study was performed to investigate the shock-wave/boundary-layer interactions on a flat plate with bleed through one or more circular holes that vent into a plenum. The bleed-hole patterns considered for the study include in-line multiple holes and staggered multiple-row holes that are configured to simulate the patterns used in inlet bleed systems of high performance aircraft. The focus of the study was to examine how the bleed through multiple holes affect bleed rate and the pressure and Mach number distributions. Since the bleed performance was found sensitive to the change in bleed conditions, a computational procedure was developed to give a good turnaround computational time for parametric studies involving changes in bleed hole geometry and the structure of shock-wave/boundary-layer flowfield. The procedure includes the grid-generation methodology and the flow simulation with solutions from the Navier-Stokes equations. The computational techniques permit analysis of complex bleed systems and make possible the investigation of a broader range of design variables associated with inlet bleed operation.

Chyu, Wei J.↗

Optimization of Adaptive Intraply Hybrid Fiber Composites with Reliability Considerations

The reliability with bounded distribution parameters (mean, standard deviation) was maximized and the reliability-based cost was minimized for adaptive intra-ply hybrid fiber composites by using a probabilistic method. The probabilistic method accounts for all naturally occurring uncertainties including those in constituent material properties, fabrication variables, structure geometry, and control-related parameters. Probabilistic sensitivity factors were computed and used in the optimization procedures. For actuated change in the angle of attack of an airfoil-like composite shell structure with an adaptive torque plate, the reliability was maximized to 0.9999 probability, with constraints on the mean and standard deviation of the actuation material volume ratio (percentage of actuation composite material in a ply) and the actuation strain coefficient. The reliability-based cost was minimized for an airfoil-like composite shell structure with an adaptive skin and a mean actuation material volume ratio as the design parameter. At a O.9-mean actuation material volume ratio, the minimum cost was obtained.

Shiao, Michael C.↗

Viscous-flow analysis of a subsonic transport aircraft high-lift system and correlation with flight data

High-lift system aerodynamics has been gaining attention in recent years. In an effort to improve aircraft performance, comprehensive studies of multi-element airfoil systems are being undertaken in wind-tunnel and flight experiments. Recent developments in Computational Fluid Dynamics (CFD) offer a relatively inexpensive alternative for studying complex viscous flows by numerically solving the Navier-Stokes (N-S) equations. Current limitations in computer resources restrict practical high-lift N-S computations to two dimensions, but CFD predictions can yield tremendous insight into flow structure, interactions between airfoil elements, and effects of changes in airfoil geometry or free-stream conditions. These codes are very accurate when compared to strictly 2D data provided by wind-tunnel testing, as will be shown here. Yet, additional challenges must be faced in the analysis of a production aircraft wing section, such as that of the NASA Langley Transport Systems Research Vehicle (TSRV). A primary issue is the sweep theory used to correlate 2D predictions with 3D flight results, accounting for sweep, taper, and finite wing effects. Other computational issues addressed here include the effects of surface roughness of the geometry, cove shape modeling, grid topology, and transition specification. The sensitivity of the flow to changing free-stream conditions is investigated. In addition, the effects of Gurney flaps on the aerodynamic characteristics of the airfoil system are predicted.

Potter, R. C.↗

Numerical analysis of stiffened shells of revolution. Volume 4: Engineer's program manual for STARS-2S shell theory automated for rotational structures - 2 (statics) digital computer program

The engineering programming information for the digital computer program for analyzing shell structures is presented. The program is designed to permit small changes such as altering the geometry or a table size to fit the specific requirements. Each major subroutine is discussed and the following subjects are included: (1) subroutine description, (2) pertinent engineering symbols and the FORTRAN coded counterparts, (3) subroutine flow chart, and (4) subroutine FORTRAN listing.

Svalbonas, V.↗

Ion propulsion

An ion engine is a plasma thruster which produces thrust by extracting ions from the plasma and accelerating them to high velocity with an electrostatic field. The ions are then neutralized and leave the engine as high velocity neutral particles. The advantages of ion engines are high specific impulse and efficiency and their ability to operate over a wide range of input powers. In comparison with other electric thrusters, the ion engine has higher efficiency and specific impulse than thermal electric devices such as the arcjet, microwave, radiofrequency and laser heated thrusters and can operate at much lower current levels than the MPD thruster. However, the thrust level for an ion engine may be lower than a thermal electric thruster of the same operating power, consistent with its higher specific impulse, and therefore ion engines are best suited for missions which can tolerate longer duration propulsive phases. The critical issue for the ion engine is lifetime, since the prospective missions may require operation for several thousands of hours. The critical components of the ion engine, with respect to engine lifetime, are the screen and accelerating grid structures. Typically, these are large metal screens that must support a large voltage difference and maintain a small gap between them. Metallic whisker growth, distortion and vibration can lead to arcing, and over a long period of time ion sputtering will erode the grid structures and change their geometry. In order to study the effects of long time operation of the grid structure, we are developing computer codes based on the Particle-In-Cell (PIC) technique and Laser Induced Fluorescence (LIF) diagnostic techniques to study the physical processes which control the performance and lifetime of the grid structures.

Meserole, J. S.↗

Ion propulsion

An ion engine is a plasma thruster which produces thrust by extracting ions from the plasma and accelerating them to high velocity with an electrostatic field. The ions are then neutralized and leave the engine as high velocity neutral particles. The advantages of ion engines are high specific impulse and efficiency and their ability to operate over a wide range of input powers. In comparison with other electric thrusters, the ion engine has higher efficiency and specific impulse than thermal electric devices such as the arcjet, microwave, radio frequency, and laser heated thrusters and can operate at much lower current levels than the MPD thruster. However, the thrust level for an ion engine may be lower than a thermal electric thruster of the same operating power, consistent with its higher specific impulse, and therefore ion engines are best suited for missions which can tolerate longer duration propulsive phases. The critical issue for the ion engine is lifetime, since the prospective missions may require operation for several thousands of hours. The critical components of the ion engine, with respect to engine lifetime, are the screen and accelerating grid structures. Typically, these are large metal screens that must support a large voltage difference and maintain a small gap between them. Metallic whisker growth, distortion, and vibration can lead to arcing, and over a long period of time ion sputtering will erode the grid structures and change their geometry. In order to study the effects of long time operation of the grid structure, we are developing computer codes based on the particle-in-cell (PIC) technique and laser induced fluorescence (LIF) diagnostic techniques to study the physical processes which control the performance and lifetime of the grid structures.

Meserole, J. S.↗

Mechanical properties of triaxially braided composites: Experimental and analytical results

This paper investigates the unnotched tensile properties of two-dimensional triaxial braid reinforced composites from both an experimental and analytical viewpoint. The materials are graphite fibers in an epoxy matrix. Three different reinforcing fiber architectures were considered. Specimens were cut from resin transfer molded (RTM) composite panels made from each braid. There were considerable differences in the observed elastic constants from different size strain gage and extensometer readings. Larger strain gages gave more consistent results and correlated better with the extensometer readings. Experimental strains correlated reasonably well with analytical predictions in the longitudinal, zero degree, fiber direction but not in the transverse direction. Tensile strength results were not always predictable even in reinforcing directions. Minor changes in braid geometry led to disproportionate strength variations. The unit cell structure of the triaxial braid was discussed with the assistence of computer analysis of the microgeometry. Photomicrographs of the braid geometry were used to improve upon the computer graphics representations of unit cells. These unit cells were used to predict the elastic moduli with various degrees of sophistication. The simple and the complex analyses were generally in agreement but none adequately matched the experimental results for all the braids.

Masters, John E.↗

Mechanical properties of triaxially braided composites - Experimental and analytical results

This paper investigates the unnotched tensile properties of two-dimensional (2D) triaxial braid-reinforced composites both experimentally and analytically. The materials are graphite fibers in an epoxy matrix. Three different reinforcing fiber architectures were considered. There were considerable differences in the observed elastic constants from different size strain gage and extensometer readings. Larger strain gages gave more consistent results and correlated better with the extensometer readings. Experimental strains correlated reasonably well with analytical predictions in the longitudinal, 0 deg, fiber direction but not in the transverse direction. Tensile strength results were not always predictable even in reinforcing directions. Minor changes in braid geometry led to disproportionate strength variations. The unit cell structure of the triaxial braid was discussed with the assistance of computer analysis of the microgeometry. Photomicrographs of braid geometry were used to improve upon the computer graphics representations of unit cells. These unit cells were used to predict the elastic moduli with various degrees of sophistication. The simple and the complex analyses were generally in agreement, but none adequately matched the experimental results for all the braids.

Masters, John E.↗

Mechanical properties of triaxially braided composites: Experimental and analytical results

The unnotched tensile properties of 2-D triaxial braid reinforced composites from both an experimental and an analytical viewpoint are studied. The materials are graphite fibers in an epoxy matrix. Three different reinforcing fiber architectures were considered. Specimens were cut from resin transfer molded (RTM) composite panels made from each braid. There were considerable differences in the observed elastic constants from different size strain gage and extensometer reading. Larger strain gages gave more consistent results and correlated better with the extensometer reading. Experimental strains correlated reasonably well with analytical predictions in the longitudinal, 0 degrees, fiber direction but not in the transverse direction. Tensile strength results were not always predictable even in reinforcing directions. Minor changes in braid geometry led to disproportionate strength variations. The unit cell structure of the triaxial braid was discussed with the assistance of computer analysis of the microgeometry. Photomicrographs of braid geometry were used to improve upon the computer graphics representations of unit cells. These unit cells were used to predict the elastic moduli with various degrees of sophistication. The simple and the complex analyses were generally in agreement but none adequately matched the experimental results for all the braids.

Masters, John E.↗

Producing Science-Ready Radar Datasets for the Retrieval of Forest Structure Parameters from Backscatter: Correcting for Terrain Topography and Changes in Vegetation Reflectivity

Radar backscatter data contain both geometric and radiometric distortions due to underlying topography and the radar viewing geometry. Our objective is to develop a radiometric correction algorithm specific to the UAVSAR system configuration that would improve retrieval of forest structure parameters. UAVSAR is an airborne Lband radar capable of repeat?pass interferometry producing images with a spatial resolution of 5m. It is characterized by an electronically steerable antenna to compensate for aircraft attitude. Thus, the computation of viewing angles (i.e. look, incidence and projection) must include aircraft attitude angles (i.e. yaw, pitch and roll) in addition to the antenna steering angle. In this presentation, we address two components of radiometric correction: area projection and vegetation reflectivity. The first correction is applied by normalization of the radar backscatter by the local ground area illuminated by the radar beam. The second is a correction due to changes in vegetation reflectivity with viewing geometry.

polarimetry↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Technique for Calculating Solution Derivatives With Respect to Geometry Parameters in a CFD Code

A solution has been developed to the challenges of computation of derivatives with respect to geometry, which is not straightforward because these are not typically direct inputs to the computational fluid dynamics (CFD) solver. To overcome these issues, a procedure has been devised that can be used without having access to the mesh generator, while still being applicable to all types of meshes. The basic approach is inspired by the mesh motion algorithms used to deform the interior mesh nodes in a smooth manner when the surface nodes, for example, are in a fluid structure interaction problem. The general idea is to model the mesh edges and nodes as constituting a spring-mass system. Changes to boundary node locations are propagated to interior nodes by allowing them to assume their new equilibrium positions, for instance, one where the forces on each node are in balance. The main advantage of the technique is that it is independent of the volumetric mesh generator, and can be applied to structured, unstructured, single- and multi-block meshes. It essentially reduces the problem down to defining the surface mesh node derivatives with respect to the geometry parameters of interest. For analytical geometries, this is quite straightforward. In the more general case, one would need to be able to interrogate the underlying parametric CAD (computer aided design) model and to evaluate the derivatives either analytically, or by a finite difference technique. Because the technique is based on a partial differential equation (PDE), it is applicable not only to forward mode problems (where derivatives of all the output quantities are computed with respect to a single input), but it could also be extended to the adjoint problem, either by using an analytical adjoint of the PDE or a discrete analog.

Mathur, Sanjay↗

User's manual to the ACTION computer code

The form and interpretation of input and output data are defined. The Analysis of Crash Transients in Inelastic and Geometrically Nonlinear structures program (ACTION) performs nonlinear transient response analysis of structures subjected to time varying loads, allowing for nonlinear, time independent material properties and large geometry changes.

Kamat, M. P.↗

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗

NDE and SHM Simulation for CFRP Composites

Ultrasound-based nondestructive evaluation (NDE) is a common technique for damage detection in composite materials. There is a need for advanced NDE that goes beyond damage detection to damage quantification and characterization in order to enable data driven prognostics. The damage types that exist in carbon fiber-reinforced polymer (CFRP) composites include microcracking and delaminations, and can be initiated and grown via impact forces (due to ground vehicles, tool drops, bird strikes, etc), fatigue, and extreme environmental changes. X-ray microfocus computed tomography data, among other methods, have shown that these damage types often result in voids/discontinuities of a complex volumetric shape. The specific damage geometry and location within ply layers affect damage growth. Realistic threedimensional NDE and structural health monitoring (SHM) simulations can aid in the development and optimization of damage quantification and characterization techniques. This paper is an overview of ongoing work towards realistic NDE and SHM simulation tools for composites, and also discusses NASA's need for such simulation tools in aeronautics and spaceflight. The paper describes the development and implementation of a custom ultrasound simulation tool that is used to model ultrasonic wave interaction with realistic 3-dimensional damage in CFRP composites. The custom code uses elastodynamic finite integration technique and is parallelized to run efficiently on computing cluster or multicore machines.

Leckey, Cara A. C.↗

Lightning Return-Stroke Current Waveforms Aloft, From Measured Field Change, Current, and Channel Geometry

Direct current measurements are available near the attachment point from both natural cloud-to-ground lightning and rocket-triggered lightning, but little is known about the rise time and peak amplitude of return-stroke currents aloft. We present, as functions of height, current amplitudes, rise times, and effective propagation velocities that have been estimated with a novel remote-sensing technique from data on 24 subsequent return strokes in six different lightning flashes that were triggering at the NASA Kennedy Space Center, FL, during 1987. The unique feature of this data set is the stereo pairs of still photographs, from which three-dimensional channel geometries were determined previously. This has permitted us to calculate the fine structure of the electric-field-change (E) waveforms produced by these strokes, using the current waveforms measured at the channel base together with physically reasonable assumptions about the current distributions aloft. The computed waveforms have been compared with observed E waveforms from the same strokes, and our assumptions have been adjusted to maximize agreement. In spite of the non-uniqueness of solutions derived by this technique, several conclusions seem inescapable: 1) The effective propagation speed of the current up the channel is usually significantly (but not unreasonably) faster than the two-dimensional velocity measured by a streak camera for 14 of these strokes. 2) Given the deduced propagation speed, the peak amplitude of the current waveform often must decrease dramatically with height to prevent the electric field from being over-predicted. 3) The rise time of the current wave front must always increase rapidly with height in order to keep the fine structure of the calculated field consistent with the observations.

Willett, J. C.↗