Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous finite element”

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.

446 records · Page 25

Lumped Parameter Models of the Central Nervous System for VIIP Research

INTRODUCTION: Current long-duration missions to the International Space Station and future exploration-class missions beyond low-Earth orbit, such as to Mars and asteroids, expose astronauts to increased risk of Visual Impairment and Intracranial Pressure (VIIP) syndrome [1]. It has been hypothesized that the headward shift of cerebral spinal fluid (CSF) and blood in microgravity may cause significant elevation of intracranial pressure (ICP), which in turn induces VIIP syndrome through biomechanical pathways [1, 2]. However, there is insufficient evidence to confirm this hypothesis. In this light, we are developing lumped-parameter models of fluid transport in the central nervous system (CNS) as a means to simulate the influence of microgravity on ICP. The CNS models will also be used in concert with the lumped parameter and finite element models of the eye described in the realted IWS abstracts submitted by Nelson et al., Feola et al. and Ethier et al. METHODS: We have developed a nine compartment CNS model (Figure 1) capable of both time-dependent and steady state fluid transport simulations, based on the works of Stevens et al. [3]. The breakdown of compartments within the model includes: vascular (3), CSF (2), brain (1) and extracranial (3). The boundary pressure in the Central Arteries [A] node is prescribed using an oscillating pressure function PA(t) simulating the carotid pulsatile pressure wave as developed by Linninger et al. [4]. For each time step, pressures are integrated through time using an adaptive-timestep 4th and 5th order Runga-Kutta solver. Once pressures are found, constitutive equations are used to solve for flowrates (Q) between each compartment. In addition to fluid flow between the different compartments, compliance (C) interactions between neighboring compartments are represented. We are also developing a second CNS model based on the works of Linninger et al. [4] which takes a more granular approach to represent the interactions of the intracranial and spinal compartments with the inclusion of arteries, arterioles, capillaries, venules, veins, venous sinus, and ventricles. The flow through the arteries, veins and CSF compartments are governed by continuity, momentum and distensibility balance equations. Furthermore, unlike the Stevens et al. approach, the Monro-Kellie doctrine of constant cranial volume and the bi-phasic nature of the brain parenchyma are implemented. These features appear to be more consistent with the physiologic and anatomical behavior of the CNS, and follow a modeling philosophy similar to the lumped parameter eye model that is intended to be integrated with the CNS model. However, Linninger’s approach has never been implemented to include hydrostatic gradient and microgravity simulation capabilities. Therefore, we aim at implement this modeling approach for spaceflight simulations and assess its overall applicability to VIIP research. OBJECTIVES: We will present verification and validation test results for both models, as well as head-to-head comparison to explore their strengths and limitations with respect to mathematical implementation and physiological significance for VIIP research. In doing so, we hope to provide some guidance to the HRP research community on how to appropriately leverage lumped parameter models for space biomedical research.

Vera, J.↗

Study of Advanced Occupant Models to Quantify Injury Risk for eVTOL Vehicles

Urban transportation is currently evolving from traditional ground-based vehicles (cars, taxis, and buses) to include air-based electric vertical take-off and landing (eVTOL) vehicles which can be utilized for on-demand transportation, cargo transport, and emergency services. These new eVTOL vehicles are designed to be small, lightweight, and able to operate autonomously without user intervention. Safety is a big part of eventual eVTOL adoption, however gaps in the consideration of safety features exist. Anthropomorphic test devices (ATDs) are used in aerospace crashworthiness standards to quantify occupant injury risk and develop improved safety designs for emergency landing situations, but the ATDs currently used in aircraft certification requirements were developed many decades ago. Developments have occurred over the years involving ATD technology, which includes a host of newer and more biofidelic ATDs such as the Test Device for Human Occupant Restraint (THOR). Increased computing power has also allowed for detailed computational human body models (HBMs) to be created, such as the Global Human Body Model Consortium (GHBMC). This study aims to assess the capability of both HBMs and new ATD designs to identify injury mechanisms within eVTOL relevant emergency landing conditions. Finite element (FE) analysis was used to expand upon full-scale and seat level impact testing conducted by researchers at the National Aeronautics and Space Administration (NASA) to look at effects of occupant model configurations on prediction of injury. The GHBMC HBM and THOR ATD models were simulated in the seat level test conditions to characterize differences between these advanced assessment tools and traditional ATDs in the isolated seat loading environment. Results identified key differences in the responses from each of the models utilized and compared their impact response in head, neck, and spinal injury metrics. The THOR model identified potential risks for head injuries due to head impacts on the seat, however it predicted lower spinal loads than the other occupant surrogates. The GHBMC showed distinctly different biomechanical responses compared to the ATD. The GHBMC model is much more deformable than the ATDs and it exhibited higher distribution of forces and increased sensitivity to the duration of acceleration pulses. Both models incorporated into this study identified key mechanisms for injury that should be considered for passenger safety in the development of these novel aircraft. In addition, this study demonstrated the value of FE modeling for running a variety of complex human surrogates to identify potential injury mechanisms for consideration in regulation and development of new aircraft. Continued research in this field to improve validation these models will only lead to safer aircraft and more comprehensive safety measures.

Crashworthiness↗

Study of Advanced Occupant Models to Quantify Injury Risk for eVTOL Vehicles

Urban transportation is currently evolving from traditional ground-based vehicles (cars, taxis, and buses) to include air-based electric vertical take-off and landing (eVTOL) vehicles which can be utilized for on-demand transportation, cargo transport, and emergency services. These new eVTOL vehicles are designed to be small, lightweight, and able to operate autonomously without user intervention. Safety is a big part of eventual eVTOL adoption, however gaps in the consideration of safety features exist. Anthropomorphic test devices (ATDs) are used in aerospace crashworthiness standards to quantify occupant injury risk and develop improved safety designs for emergency landing situations, but the ATDs currently used in aircraft certification requirements were developed many decades ago. Developments have occurred over the years involving ATD technology, which includes a host of newer and more biofidelic ATDs such as the Test Device for Human Occupant Restraint (THOR). Increased computing power has also allowed for detailed computational human body models (HBMs) to be created, such as the Global Human Body Model Consortium (GHBMC). This study aims to assess the capability of both HBMs and new ATD designs to identify injury mechanisms within eVTOL relevant emergency landing conditions. Finite element (FE) analysis was used to expand upon full-scale and seat level impact testing conducted by researchers at the National Aeronautics and Space Administration (NASA) to look at effects of occupant model configurations on prediction of injury. The GHBMC HBM and THOR ATD models were simulated in the seat level test conditions to characterize differences between these advanced assessment tools and traditional ATDs in the isolated seat loading environment. Results identified key differences in the responses from each of the models utilized and compared their impact response in head, neck, and spinal injury metrics. The THOR model identified potential risks for head injuries due to head impacts on the seat, however it predicted lower spinal loads than the other occupant surrogates. The GHBMC showed distinctly different biomechanical responses compared to the ATD. The GHBMC model is much more deformable than the ATDs and it exhibited higher distribution of forces and increased sensitivity to the duration of acceleration pulses. Both models incorporated into this study identified key mechanisms for injury that should be considered for passenger safety in the development of these novel aircraft. In addition, this study demonstrated the value of FE modeling for running a variety of complex human surrogates to identify potential injury mechanisms for consideration in regulation and development of new aircraft. Continued research in this field to improve validation these models will only lead to safer aircraft and more comprehensive safety measures.

Crashworthiness↗

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗

Efficient Reformulation of the Thermoelastic Higher-order Theory for Fgms

Functionally graded materials (FGMs) are characterized by spatially variable microstructures which are introduced to satisfy given performance requirements. The microstructural gradation gives rise to continuously or discretely changing material properties which complicate FGM analysis. Various techniques have been developed during the past several decades for analyzing traditional composites and many of these have been adapted for the analysis of FGMs. Most of the available techniques use the so-called uncoupled approach in order to analyze graded structures. These techniques ignore the effect of microstructural gradation by employing specific spatial material property variations that are either assumed or obtained by local homogenization. The higher-order theory for functionally graded materials (HOTFGM) is a coupled approach developed by Aboudi et al. (1999) which takes the effect of microstructural gradation into consideration and does not ignore the local-global interaction of the spatially variable inclusion phase(s). Despite its demonstrated utility, however, the original formulation of the higher-order theory is computationally intensive. Herein, an efficient reformulation of the original higher-order theory for two-dimensional elastic problems is developed and validated. The use of the local-global conductivity and local-global stiffness matrix approach is made in order to reduce the number of equations involved. In this approach, surface-averaged quantities are the primary variables which replace volume-averaged quantities employed in the original formulation. The reformulation decreases the size of the global conductivity and stiffness matrices by approximately sixty percent. Various thermal, mechanical, and combined thermomechanical problems are analyzed in order to validate the accuracy of the reformulated theory through comparison with analytical and finite-element solutions. The presented results illustrate the efficiency of the reformulation and its advantages in analyzing functionally graded materials.

Bansal, Yogesh↗

Silicon Nitride Plates for Turbine Blade Application: FEA and NDE Assessment

Engine manufacturers are continually attempting to improve the performance and the overall efficiency of internal combustion engines. The thermal efficiency is typically improved by raising the operating temperature of essential engine components in the combustion area. This reduces the heat loss to a cooling system and allows a greater portion of the heat to be used for propulsion. Further improvements can be achieved by diverting part of the air from the compressor, which would have been used in the combustor for combustion purposes, into the turbine components. Such a process is called active cooling. Increasing the operating temperature, decreasing the cooling air, or both can improve the efficiency of the engine. Furthermore, lightweight, strong, tough hightemperature materials are required to complement efficiency improvement for nextgeneration gas turbine engines that can operate with minimum cooling. Because of their low-density, high-temperature strength, and thermal conductivity, ceramics are being investigated as potential materials for replacing ordinary metals that are currently used for engine hot section components. Ceramic structures can withstand higher operating temperatures and other harsh environmental factors. In addition, their low densities relative to metals helps condense component mass (ref. 1). The objectives of this program at the NASA Glenn Research Center are to develop manufacturing technology, a thermal barrier coating/environmental barrier coating (TBC/EBC), and an analytical modeling capability to predict thermomechanical stresses, and to do minimal burner rig tests of silicon nitride (Si3N4) and SiC/SiC turbine nozzle vanes under simulated engine conditions. Furthermore, and in support of the latter objectives, an optimization exercise using finite element analysis and nondestructive evaluation (NDE) was carried out to characterize and evaluate silicon nitride plates with cooling channels.

Abdul-Aziz, Ali↗

Superalloy Disk With Dual-Grain Structure Spin Tested

Advanced nickel-base disk alloys for future gas turbine engines will require greater temperature capability than current alloys, but they must also continue to deliver safe, reliable operation. An advanced, nickel-base disk alloy, designated Alloy 10, was selected for evaluation in NASA s Ultra Safe Propulsion Project. Early studies on small test specimens showed that heat treatments that produced a fine grain microstructure promoted high strength and long fatigue life in the bore of a disk, whereas heat treatments that produced a coarse grain microstructure promoted optimal creep and crack growth resistance in the rim of a disk. On the basis of these results, the optimal combination of performance and safety might be achieved by utilizing a heat-treatment technology that could produce a fine grain bore and coarse grain rim in a nickel-base disk. Alloy 10 disks that were given a dual microstructure heat treatment (DMHT) were obtained from NASA s Ultra-Efficient Engine Technology (UEET) Program for preliminary evaluation. Data on small test specimens machined from a DMHT disk were encouraging. However, the benefit of the dual grain structure on the performance and reliability of the entire disk still needed to be demonstrated. For this reason, a high temperature spin test of a DMHT disk was run at 20 000 rpm and 1500 F at the Balancing Company of Dayton, Ohio, under the direction of NASA Glenn Research Center personnel. The results of that test showed that the DMHT disk exhibited significantly lower crack growth than a disk with a fine grain microstructure. In addition, the results of these tests could be accurately predicted using a two-dimensional, axisymmetric finite element analysis of the DMHT disk. Although the first spin test demonstrated a significant performance advantage associated with the DMHT technology, a second spin test on the DMHT disk was run to determine burst margin. The disk burst in the web at a very high speed, over 39 000 rpm, in line with the predicted location and speed. Furthermore, significant growth of the disk was observed before failure, in line with predictions, clearly demonstrating the reliability and safety of the DMHT technology. Although successful spin testing in Ultra Safe's Nickel Disk Program represents a significant milestone for DMHT technology, realistic engine operation will require repeated loading of a DMHT nickel disk. For this reason, a cyclic spin test study of DMHT nickel disk technology has been proposed to start in fiscal year 2003 under NASA's Aviation Safety Program. The goal of this program will be to determine the fatigue failure mechanism in DMHT nickel disks, thereby demonstrating the reliability and safety of DMHT technology under repetitive loading conditions encountered in realistic engine operation.

Gayda, John↗

Risk of Bone Fracture due to Spaceflight-induced Changes to Bone

The Human Research Program (HRP) is taking a biomechanical approach to assessing subsequent fracture risk in active astronauts as a consequence of spaceflight exposure. Triennial testing of active and retired astronauts continues as Space Medicine monitors for a premature diagnosis of primary osteoporosis, which is associated with age-related bone loss and skeletal fragility. This updated 2023 Fracture Evidence Report has expanded the description of skeletal changes to capture the full effects of spaceflight on bone: • Recent analysis of fractures from health records of the full astronaut cohort, suggests an increased incidence rate of hip and spine fractures in astronauts following longer spaceflight duration flights compared to incident rates found in shorter duration flights. • Routine preflight-to-postflight surveillance by DXA (dual-energy x-ray absorptiometry) does not provide the full detection of loss and recovery in the long-duration (LD) astronaut nor the full recovery of hip trabecular bone. • The inclusion of hip quantitative computed tomography (QCT) in flight studies delineates effects of spaceflight, of countermeasures, and of post-flight recovery on cortical and trabecular bone parameters—some of which are verified predictors of hip fracture in the aged. • QCT detects and compares the distinct countermeasure effects of the pharmaceutical alendronate and of resistive exercise (on the Advanced Resistive Exercise Device, ARED) in specific cortical and trabecular bone sub-regions during spaceflight. • A published comparison of QCT-determined loss rates of hip trabecular bone in LD astronauts compared to terrestrial cohorts suggest that accelerated loss rates in trabecular volumetric BMD (vBMD) during spaceflight might be analogous to skeletal effects of accelerated loss rates in females due to menopause, potentially leading to disruptions in trabecular microarchitecture. • A dataset of finite element (FE) estimates of hip bone strength in aging terrestrial cohorts (spanning astronaut age-range) provides comparative context for changes in the FE of hip bone strength in LD astronauts, including comparison (force unit of newton) to percentiles (50th, 75th, and 100th) of sex-matched aged humans with hip fractures. • Risk of fracture, due to the mechanical overloading of bones, is being updated with IMPACT, the next-generation tool suite for probabilistic risk assessment (PRA) for exploration missions being created by the Exploration and Medical Capability Element (ExMC) at NASA. IMPACT is currently in development and will not be included in this 2024 update. The risk for fracture necessitates understanding the relationship between applied loads to bone and the biomechanical competence of bone. The Risk for Early Onset Osteoporosis focuses on the weakened condition of bone (including development of new technologies, measurements of novel skeletal attributes, translation of multiple measures to an index of bone fragility, and interpretations of data used to reflect a weakened bone), while the Risk for Fracture assesses factors that influence the probability an astronaut would encounter applied loads exceeding the biomechanical competence of bones, resulting in fracture. This Evidence Report combines the research gaps and tasks associated with both risks.

Fracture↗

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗

Probabilistic Risk Assessment for Astronaut Post Flight Bone Fracture

Introduction: Space flight potentially reduces the loading that bone can resist before fracture. This reduction in bone integrity may result from a combination of factors, the most common reported as reduction in astronaut BMD. Although evaluating the condition of bones continues to be a critical aspect of understanding space flight fracture risk, defining the loading regime, whether on earth, in microgravity, or in reduced gravity on a planetary surface, remains a significant component of estimating the fracture risks to astronauts. This presentation summarizes the concepts, development, and application of NASA's Bone Fracture Risk Module (BFxRM) to understanding pre-, post, and in mission astronaut bone fracture risk. The overview includes an assessment of contributing factors utilized in the BFxRM and illustrates how new information, such as biomechanics of space suit design or better understanding of post flight activities may influence astronaut fracture risk. Opportunities for the bone mineral research community to contribute to future model development are also discussed. Methods: To investigate the conditions in which spaceflight induced changes to bone plays a critical role in post-flight fracture probability, we implement a modified version of the NASA Bone Fracture Risk Model (BFxRM). Modifications included incorporation of variations in physiological characteristics, post-flight recovery rate, and variations in lateral fall conditions within the probabilistic simulation parameter space. The modeled fracture probability estimates for different loading scenarios at preflight and at 0 and 365 days post-flight time periods are compared. Results: For simple lateral side falls, mean post-flight fracture probability is elevated over mean preflight fracture probability due to spaceflight induced BMD loss and is not fully recovered at 365 days post-flight. In the case of more energetic falls, such as from elevated heights or with the addition of lateral movement, the contribution of space flight quality changes is much less clear, indicating more granular assessments, such as Finite Element modeling, may be needed to further assess the risks in these scenarios.

Risk Analysis↗

Local Ultrasonic Resonance Spectroscopy of Lithium Metal Batteries for Aerospace Applications

As next-generation aircraft and vehicles continue to develop, so do their associated energy demands. Lithium metal batteries are a leading candidate to fulfill this energy requirement, but these batteries are prone to internal dendrite defects that can lead to catastrophic thermal runaway events. Current battery management systems are capable of mitigating such risks, but are unable to detect such defects until thermal runaway has already begun. Various nondestructive evaluation (NDE) techniques, particularly ultrasonic NDE, can directly monitor internal battery parameters giving them the potential to detect critical defects prior to catastrophic failure. However, most of the current battery NDE research has focused on improved battery state-of-charge (SOC) and stateof- health (SOH) monitoring with little emphasis on critical defect detection. Thus, a measurement technique sensitive to subtle battery defects is needed. In addition, the complex mechanics of ultrasound in porous, thin, multilayered batteries prompt the use of physics-based simulation to guide inspections. In this work, an ultrasonic NDE technique has been developed utilizing frequency domain analysis of local battery resonances to detect the presence of battery defects. This technique is a practical extension of local ultrasonic resonance spectroscopy (LURS) – which previously required non-contact laser ultrasonics – to measurements with piezoelectric contact and immersion scan transducers. To extend the technique to work with piezoelectric transducers, ultrasonic battery measurements were compared to a sans-battery calibration measurement. Then, a linear systems deconvolution was used to eliminate the transfer functions of extraneous factors such as the transducer and electronics, leaving only the frequency-dependent battery reflection coefficient. The LURS technique was first validated on stainless steel and aluminum plates, producing reflection coefficients in line with analytical and numerical finite element modeling (FEM) results. Functioning Li-metal pouch cells were then seeded with lithium chip defects prior to LURS measurements. The presence of these defects is shown to cause a measurable shift in the battery’s through-thickness local resonances. 2D, frequency-domain poroelastic models of ultrasonic propagation in a single-cell lithium metal pouch battery were created and corroborated these findings. Thus, this work has both extended and proven the feasibility of the LURS technique in the detection of local battery defects.

Ultrasound↗

A new look at the simultaneous analysis and design of structures

The minimum weight optimization of structural systems, subject to strength and displacement constraints as well as size side constraints, was investigated by the Simultaneous ANalysis and Design (SAND) approach. As an optimizer, the code NPSOL was used which is based on a sequential quadratic programming (SQP) algorithm. The structures were modeled by the finite element method. The finite element related input to NPSOL was automatically generated from the input decks of such standard FEM/optimization codes as NASTRAN or ASTROS, with the stiffness matrices, at present, extracted from the FEM code ANALYZE. In order to avoid ill-conditioned matrices that can be encountered when the global stiffness equations are used as additional nonlinear equality constraints in the SAND approach (with the displacements as additional variables), the matrix displacement method was applied. In this approach, the element stiffness equations are used as constraints instead of the global stiffness equations, in conjunction with the nodal force equilibrium equations. This approach adds the element forces as variables to the system. Since, for complex structures and the associated large and very sparce matrices, the execution times of the optimization code became excessive due to the large number of required constraint gradient evaluations, the Kreisselmeier-Steinhauser function approach was used to decrease the computational effort by reducing the nonlinear equality constraint system to essentially a single combined constraint equation. As the linear equality and inequality constraints require much less computational effort to evaluate, they were kept in their previous form to limit the complexity of the KS function evaluation. To date, the standard three-bar, ten-bar, and 72-bar trusses have been tested. For the standard SAND approach, correct results were obtained for all three trusses although convergence became slower for the 72-bar truss. When the matrix displacement method was used, correct results were still obtained, but the execution times became excessive due to the large number of constraint gradient evaluations required. Using the KS function, the computational effort dropped, but the optimization seemed to become less robust. The investigation of this phenomenon is continuing. As an alternate approach, the code MINOS for the optimization of sparse matrices can be applied to the problem in lieu of the Kreisselmeier-Steinhauser function. This investigation is underway.

Striz, Alfred G.↗

Continuation of tailored composite structures of ordered staple thermoplastic material

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Santare, Michael H.↗

Carbon Nanofiber-Based, High-Frequency, High-Q, Miniaturized Mechanical Resonators

High Q resonators are a critical component of stable, low-noise communication systems, radar, and precise timing applications such as atomic clocks. In electronic resonators based on Si integrated circuits, resistive losses increase as a result of the continued reduction in device dimensions, which decreases their Q values. On the other hand, due to the mechanical construct of bulk acoustic wave (BAW) and surface acoustic wave (SAW) resonators, such loss mechanisms are absent, enabling higher Q-values for both BAW and SAW resonators compared to their electronic counterparts. The other advantages of mechanical resonators are their inherently higher radiation tolerance, a factor that makes them attractive for NASA s extreme environment planetary missions, for example to the Jovian environments where the radiation doses are at hostile levels. Despite these advantages, both BAW and SAW resonators suffer from low resonant frequencies and they are also physically large, which precludes their integration into miniaturized electronic systems. Because there is a need to move the resonant frequency of oscillators to the order of gigahertz, new technologies and materials are being investigated that will make performance at those frequencies attainable. By moving to nanoscale structures, in this case vertically oriented, cantilevered carbon nanotubes (CNTs), that have larger aspect ratios (length/thickness) and extremely high elastic moduli, it is possible to overcome the two disadvantages of both bulk acoustic wave (BAW) and surface acoustic wave (SAW) resonators. Nano-electro-mechanical systems (NEMS) that utilize high aspect ratio nanomaterials exhibiting high elastic moduli (e.g., carbon-based nanomaterials) benefit from high Qs, operate at high frequency, and have small force constants that translate to high responsivity that results in improved sensitivity, lower power consumption, and im - proved tunablity. NEMS resonators have recently been demonstrated using topdown, lithographically fabricated ap - proaches to form cantilever or bridgetype structures. Top-down approaches, however, rely on complicated and expensive e-beam lithography, and often require a release mechanism. Reso - nance effects in structures synthesized using bottom-up approaches have also recently been reported based on carbon nanotubes, but such approaches have relied on a planar two-dimensional (2D) geometry. In this innovation, vertically aligned tubes synthesized using a bottom- up approach have been considered, where the vertical orientation of the tubes has the potential to increase integration density even further. The simulation of a vertically oriented, cantilevered carbon nanotube was performed using COMSOL Multi - physics, a finite element simulation package. All simulations were performed in a 2D geometry that provided consistent results and minimized computational complexity. The simulations assumed a vertically oriented, cantilevered nanotube of uniform density (1.5 g/cu cm). An elastic modulus was assumed to be 600 GPa, relative permittivity of the nanotube was assumed to be 5.0, and Poisson s ratio was assumed to be 0.2. It should be noted that the relative permittivity and Poisson s ratio for the nanotubes of interest are not known accurately. However, as in previous simulations, the relative permittivity and Poisson s ratios were treated as weak variables in the simulation, and no significant changes were recognized when these variables were varied.

Kaul, Anupama B.↗