Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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 811 records · Page 45

Methods Developed by the Tools for Engine Diagnostics Task to Monitor and Predict Rotor Damage in Real Time

Tools for Engine Diagnostics is a major task in the Propulsion System Health Management area of the Single Aircraft Accident Prevention project under NASA s Aviation Safety Program. The major goal of the Aviation Safety Program is to reduce fatal aircraft accidents by 80 percent within 10 years and by 90 percent within 25 years. The goal of the Propulsion System Health Management area is to eliminate propulsion system malfunctions as a primary or contributing factor to the cause of aircraft accidents. The purpose of Tools for Engine Diagnostics, a 2-yr-old task, is to establish and improve tools for engine diagnostics and prognostics that measure the deformation and damage of rotating engine components at the ground level and that perform intermittent or continuous monitoring on the engine wing. In this work, nondestructive-evaluation- (NDE-) based technology is combined with model-dependent disk spin experimental simulation systems, like finite element modeling (FEM) and modal norms, to monitor and predict rotor damage in real time. Fracture mechanics time-dependent fatigue crack growth and damage-mechanics-based life estimation are being developed, and their potential use investigated. In addition, wireless eddy current and advanced acoustics are being developed for on-wing and just-in-time NDE engine inspection to provide deeper access and higher sensitivity to extend on-wing capabilities and improve inspection readiness. In the long run, these methods could establish a base for prognostic sensing while an engine is running, without any overt actions, like inspections. This damage-detection strategy includes experimentally acquired vibration-, eddy-current- and capacitance-based displacement measurements and analytically computed FEM-, modal norms-, and conventional rotordynamics-based models of well-defined damages and critical mass imbalances in rotating disks and rotors.

Baaklini, George Y.↗

Solar Sail Topology Variations Due to On-Orbit Thermal Effects

The objective of this research was to predict the influence of non-uniform temperature distribution on solar sail topology and the effect of such topology variations on sail performance (thrust, torque). Specifically considered were the thermal effects due to on orbit attitude control maneuvers. Such maneuvers are expected to advance the sail to a position off-normal to the sun by as much as 35 degrees; a solar sail initially deformed by typical pre-tension and solar pressure loads may suffer significant thermally induced strains due to the non-uniform heating caused by these maneuvers. This on-orbit scenario was investigated through development of an automated analytical shape model that iterates many times between sail shape and sail temperature distribution before converging on a final coupled thermal structural affected sail topology. This model utilizes a validated geometrically non-linear finite element model and a thermal radiation subroutine. It was discovered that temperature gradients were deterministic for the off-normal solar angle cases as were thermally induced strains. Performance effects were found to be moderately significant but not as large as initially suspected. A roll torque was detected, and the sail center of pressure shifted by a distance that may influence on-orbit sail control stability.

Banik, Jeremy A.↗

Aeroelastic Analysis of SUGAR Truss-Braced Wing Wind-Tunnel Model Using FUN3D and a Nonlinear Structural Model

Considerable attention has been given in recent years to the design of highly flexible aircraft. The results of numerous studies demonstrate the significant performance benefits of strut-braced wing (SBW) and trussbraced wing (TBW) configurations. Critical aspects of the TBW configuration are its larger aspect ratio, wing span and thinner wings. These aspects increase the importance of considering fluid/structure and control system coupling. This paper presents high-fidelity Navier-Stokes simulations of the dynamic response of the flexible Boeing Subsonic Ultra Green Aircraft Research (SUGAR) truss-braced wing wind-tunnel model. The latest version of the SUGAR TBW finite element model (FEM), v.20, is used in the present simulations. Limit cycle oscillations (LCOs) of the TBW wing/strut/nacelle are simulated at angle-of-attack (AoA) values of -1, 0 and +1 degree. The modal data derived from nonlinear static aeroelastic MSC.Nastran solutions are used at AoAs of -1 and +1 degrees. The LCO amplitude is observed to be dependent on AoA. LCO amplitudes at -1 degree are larger than those at +1 degree. The LCO amplitude at zero degrees is larger than either -1 or +1 degrees. These results correlate well with both wind-tunnel data and the behavior observed in previous studies using linear aerodynamics. The LCO onset at zero degrees AoA has also been computed using unloaded v.20 FEM modes. While the v.20 model increases the dynamic pressure at which LCO onset is observed, it is found that the LCO onset at and above Mach 0.82 is much different than that produced by an earlier version of the FEM, v. 19.

Bartels, Robert E.↗

A Model-Based, Bayesian Solution for Characterization of Complex Damage Scenarios in Aerospace Composite Structures

Ultrasonic damage detection and characterization is commonly used in nondestructive evaluation (NDE) of aerospace composite components. In recent years there has been an increased development of guided wave based methods. In real materials and structures, these dispersive waves result in complicated behavior in the presence of complex damage scenarios. Model-based characterization methods utilize accurate three dimensional finite element models (FEMs) of guided wave interaction with realistic damage scenarios to aid in defect identification and classification. This work describes an inverse solution for realistic composite damage characterization by comparing the wavenumber-frequency spectra of experimental and simulated ultrasonic inspections. The composite laminate material properties are first verified through a Bayesian solution (Markov chain Monte Carlo), enabling uncertainty quantification surrounding the characterization. A study is undertaken to assess the efficacy of the proposed damage model and comparative metrics between the experimental and simulated output. The FEM is then parameterized with a damage model capable of describing the typical complex damage created by impact events in composites. The damage is characterized through a transdimensional Markov chain Monte Carlo solution, enabling a flexible damage model capable of adapting to the complex damage geometry investigated here. The posterior probability distributions of the individual delamination petals as well as the overall envelope of the damage site are determined.

H. Reed↗

The Middeck Active Control Experiment (MACE)

The Middeck Active Control Experiment (MACE) is a NASA In-Step and Control Structure Interaction (CSI) Office funded Shuttle middeck experiment. The objective is to investigate the extent to which closed-loop behavior of flexible spacecraft in zero-gravity (0-g) can be predicted. This prediction becomes particularly difficult when dynamic behavior during ground testing exhibits extensive suspension and direct gravity coupling. On-orbit system identification and control reconfiguration is investigated to improve performance which would otherwise be limited due to errors in prediction. The program is presently in its preliminary design phase with launch expected in the summer of 1994. The MACE test article consists of three attitude control torque wheels, a two axis gimballing payload, inertial sensors and a flexible support structure. With the acquisition of a second payload, this will represent a multiple payload platform with significant structural flexibility. This paper presents on-going work in the areas of modelling and control of the MACE test article in the zero and one-gravity environments. Finite element models, which include suspension and gravity effects, and measurement models, derived from experimental data, are used as the basis for Linear Quadratic Gaussian controller designs. Finite element based controllers are analytically used to study the differences in closed-loop performance as the test article transitions between the 0-g and 1-g environments. Measurement based controllers are experimentally applied to the MACE test article in the 1-g environment and achieve over an order of magnitude improvement in payload pointing accuracy when disturbed by a broadband torque disturbance. The various aspects of the flight portion of the experiment are also discussed.

Miller, David W.↗

Flutter Analysis of the Transonic Truss-Braced Wing Aircraft Using Transonic Correction

This paper describes a flutter analysis method for the Transonic Truss-Braced Wing aircraft using a vortex lattice method coupled to an unsteady transonic correction method to account for unsteady aerodynamics in transonic flow. A steady-state vortex-lattice model of the Truss-Braced Wing aircraft is developed using vortex lattice code VSPAERO. A transonic and viscous flow correction method is implemented in the VSPAERO model to account for steady-state transonic and viscous flow effects using transonic small disturbance code TSFOIL coupled to an in-house integral boundary layer code. In addition, a wing-strut interference correction method is developed to account for the transonic interference aerodynamics in the strut juncture region using high-fidelity CFD code FUN3D. A structural dynamic finite-element model of the Truss-Braced Wing aircraft is developed using BEAM3D in-house finite-element code and is coupled to the VSPAERO. The BEAM3D model includes a geometric nonlinearity due to the tension in the strut which causes a deflection-dependent nonlinear stiffness. An unsteady transonic correction method is developed to better capture the unsteady aerodynamics in transonic flow. The unsteady transonic correction method makes use of the Theodorsen’s theory to account for the amplitude and phase shift of the unsteady lift coefficient in transonic flow. A preliminary flutter analysis of the Truss-Braced Wing aircraft is conduct to illustrate the unsteady transonic correction approach.

Nhan Nguyen↗

Computer Aided Modeling and Post Processing with NASTRAN Analysis

Computer aided engineering systems are invaluable tools in performing NASTRAN finite element analysis. These techniques are implemented in both the pre-processing and post-processing phases of the NASTRAN analysis. The finite element model development, or pre-processing phase, was automated with a computer aided modeling program called Supertabl, and the review and interpretation of the results of the NASTRAN analysis, or post-processing phase, was automated with a computer aided plotting program called Output Display. An intermediate program, Nasplot, which was developed in-house, has also helped to cut down on the model checkout time and reduce errors in the model. An interface has been established between the finite element computer aided engineering system and the Learjet computer aided design system whereby data can be transferred back and forth between the two. These systems have significantly improved productivity and the ability to perform NASTRAN analysis in response to product development requests.

Boroughs, R. R.↗

Mapping methods for computationally efficient and accurate structural reliability

Mapping methods are developed to improve the accuracy and efficiency of probabilistic structural analyses with coarse finite element meshes. The mapping methods consist of the following: (1) deterministic structural analyses with fine (convergent) finite element meshes; (2) probabilistic structural analyses with coarse finite element meshes; (3) the relationship between the probabilistic structural responses from the coarse and fine finite element meshes; and (4) a probabilistic mapping. The results show that the scatter in the probabilistic structural responses and structural reliability can be efficiently predicted using a coarse finite element model and proper mapping methods with good accuracy. Therefore, large structures can be efficiently analyzed probabilistically using finite element methods.

Shiao, Michael C.↗

Experimental and Analytical Evaluation of a Composite Honeycomb Deployable Energy Absorber

In 2006, the NASA Subsonic Rotary Wing Aeronautics Program sponsored the experimental and analytical evaluation of an externally deployable composite honeycomb structure that is designed to attenuate impact energy during helicopter crashes. The concept, which is designated the Deployable Energy Absorber (DEA), utilizes an expandable Kevlar honeycomb structure to dissipate kinetic energy through crushing. The DEA incorporates a unique flexible hinge design that allows the honeycomb to be packaged and stowed flat until needed for deployment. A variety of deployment options such as linear, radial, and/or hybrid methods can be used. Experimental evaluation of the DEA utilized a building block approach that included material characterization testing of its constituent, Kevlar -129 fabric/epoxy, and flexural testing of single hexagonal cells. In addition, the energy attenuation capabilities of the DEA were demonstrated through multi-cell component dynamic crush tests, and vertical drop tests of a composite fuselage section, retrofitted with DEA blocks, onto concrete, water, and soft soil. During each stage of the DEA evaluation process, finite element models of the test articles were developed and simulations were performed using the explicit, nonlinear transient dynamic finite element code, LS-DYNA. This report documents the results of the experimental evaluation that was conducted to assess the energy absorption capabilities of the DEA.

Jackson, Karen E.↗

Machine learning–assisted prediction of heat fluxes through thermally anisotropic building envelopes

Thermally anisotropic building envelope (TABE) is a novel active building envelope that can save energy use to maintain thermal comfort in buildings by redirecting heat and coolness from building envelopes to thermal loops. Finite element models (FEMs) can be used to compute the heat fluxes through TABEs, but the high computational cost of finite element simulations has prevented parametric studies and design optimizations. This paper proposes a domain knowledge–informed, finite element–based machine learning framework to reduce the computation cost for the energy management of buildings installed with TABE that uses a ground thermal loop. First, the training heat flux data set was generated by FEM simulations with different thermal loop schedules. Then, both shallow learning models (i.e., multivariate linear regression and eXtreme Gradient Boost, or XGBoost) and a deep learning model (i.e., deep neural network, or DNN) were trained to predict the heat fluxes. Domain knowledge was used for data preprocessing and feature selection. Finally, the suitability of the selected machine learning model was tested under different thermal loop schedules. Herein, the case study results showed that: (1) XGBoost can be as accurate as DNN (coefficient of determination equal to 0.81) with much less training time; (2) the annual energy cost savings for different thermal loop schedules obtained by the XGBoost-predicted and FEM-calculated heat fluxes are consistent, having a difference of only 4%; and (3) XGBoost can reduce the computation time for the annual energy analysis of the case study building with a given thermal loop schedule from around 12 h by using FEM to less than 1 min.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Modeling aluminum-lithium alloy welding characteristics

The purpose of this project was to develop a finite element model of the heat-affected zone in the vicinity of a weld line on a plate in order to determine an accurate plastic strain history. The resulting plastic strain increments calculated by the finite element program were then to be used to calculate the measure of damage D. It was hoped to determine the effects of varying welding parameters, such as beam power, efficiency, and weld speed, and the effect of different material properties on the occurrence of microfissuring. The results were to be compared first to the previous analysis of Inconel 718, and then extended to aluminum 2195.

Bernstein, Edward L.↗

Hierarchic plate and shell models based on p-extension

Formulation of hierarchic sequences of finite element models for beams, arches, plates and shells based on the principle of virtual work is described. The exact solutions corresponding to models in the hierarchic sequence converge to the exact solution of the fully three-dimensional problem of linear elasticity. The stopping criterion is that the functionals of interest must be substantially independent of the choice of the model. This process is closely related to p-extensions. Aspects of implementation are discussed in connection with axisymmetric shells and an example is presented. An application of superconvergent extraction methods for the computation of stress resultants is demonstrated.

Szabo, B. A.↗

A Flexible Method for Producing F.E.M. Analysis of Bone Using Open-Source Software

Individuals who experience decreases in load-bearing bone densities can be subject to a higher risk of bone fracture during daily activity. Astronauts may lose up to nine percent of their load-bearing bone density for every month they spend in space [1]. Because of this, specialized countermeasures reduce percent loss in bone density and reduce fracture risk upon returning to Earth. Astronauts will typically not be at risk for fracture during spaceflight, because of the lesser loads experienced in microgravity conditions. However, once back on Earth, astronauts have an increased risk for bone fracture as a result of weakened bone and return to 1G conditions [2]. It is therefore important to understand the significance of any bone density loss in addition to developing exercises in an attempt to limit losses in bone strength. NASA seeks to develop a deeper understanding of fracture risk through the development of a computational bone strength model to assess the bone fracture risk of astronauts pre-flight and post-flight. This study addresses the several key processes needed to develop such strength analyses using medical image processing and finite element modeling.

biomechanics↗

Quantifying spectral albedo effects on bifacial photovoltaic module measurements and system model predictions

We provide a comprehensive analysis of the effect of spectral albedo on photovoltaic (PV) module measurements and system model predictions. We demonstrate how to account for albedo in indoor bifacial device measurements by adjusting the applied irradiance using the scaled rear irradiance method, exemplified on fabricated silicon heterojunction (SHJ) modules. System model performance is studied using a detailed 3D finite-element model, DUET, for fixed-tilt and horizontal single-axis tracked (SAT) arrays between 15 and 75°N. Spectral effects cause variations in measured SHJ module short-circuit current up to 2% and efficiency variation up to 0.3% abs. We further demonstrate that rear-side spectral mismatch factors (SMMs) resulting from including or omitting spectral albedo in PV system modeling vary between ±13%, while total (front+rear) SMMs vary up to 3%, depending on the deployment configuration and latitude. SAT array SMMs are weakly correlated with latitude, while fixed-tilt array SMMs increase with latitude, driven by an increasing proportion of ground-reflected light on the front-side of modules. Ground-reflections can constitute between 2% and 32% of total incident module irradiance, with notably high (>10%) contributions for fixed-tilt arrays at high latitude. Effects of spectral albedo are most significant for: (1) fixed-tilt deployments at high latitudes, (2) wide bandgap technologies such as perovskite and cadmium telluride cells, (3) albedos which vary steeply over the technology's absorption range, and (4) high albedo ground covers. Overall, we demonstrate that omitting spectral albedo effects can result in PV measurement and system-level modeling uncertainties on the order of several percent in these cases.

14 SOLAR ENERGY↗

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗

NASTRAN thermal analyzer: A general purpose finite element heat transfer computer program

The program not only can render temperature distributions in solids subjected to various thermal boundary conditions, including effects of diffuse-gray thermal radiation, but is fully compatible in capacity and in the finite-element model representation with that of its structural counterpart in the NASTRAN system. The development history of the finite-element approach for determining temperatures is summarized. The scope of analysis capability, program structure, features, and limitations are given with the objective of providing NASTRAN users with an overall veiw of the NASTRAN thermal analyzer.

Lee, H.↗

The LDCM actuator for vibration suppression

A linear dc motor (LDCM) has been proposed as an actuator for the COFS I mast and the COFS program ground test Mini-Mast. The basic principles of operation of the LDCM as an actuator for vibration suppression in large flexible structures are reviewed. Because of force and stroke limitations, control loops are required to stabilize the actuator, which results in a non-standard actuator-plant configuration. A simulation model that includes LDCM actuator control loops and a finite element model of the Mast is described, with simulation results showing the excitation capability of the actuator.

Ide, Eric N.↗

A study of fiber volume fraction effects in notched unidirectional SCS-6/Ti-15V-3Cr-3Al-3Sn composite

Notched unidirectional SCS-6/Ti-15-3 composite of three different fiber volume fractions (vf = 0.15, 0.37, and 0.41) was investigated for various room temperature microstructural and material properties including: fatigue crack initiation, fatigue crack growth, and fracture toughness. While the matrix hardness is similar for all fiber volume fractions, the fiber/matrix interfacial shear strength and matrix residual stress increases with fiber volume fraction. The composite fatigue crack initiation stress is shown to be matrix controlled and occurs when the net maximum matrix stress approaches the endurance limit stress of the matrix. A model is presented which includes residual stresses and presents the composite initiation stress as a function of fiber volume fraction. This model predicts a maximum composite initiation stress at vf approximately 0.15 which agrees with the experimental data. The applied composite stress levels were increased as necessary for continued crack growth. The applied Delta(K) values at crack arrest increase with fiber volume fraction by an amount better approximated using an energy based formulation rather than when scaled linear with modulus. After crack arrest, the crack growth rate exponents for vf37 and vf41 were much lower and toughness much higher, when compared to the unreinforced matrix, because of the bridged region which parades with the propagating fatigue crack. However, the vf15 material exhibited a higher crack growth rate exponent and lower toughness than the unreinforced matrix because once the bridged fibers nearest the crack mouth broke, the stress redistribution broke all bridged fibers, leaving an unbridged crack. Degraded, unbridged behavior is modeled using the residual stress state in the matrix ahead of the crack tip. Plastic zone sizes were directly measured using a metallographic technique and allow prediction of an effective matrix stress intensity which agrees with the fiber pressure model if residual stresses are considered. The sophisticated macro/micro finite element models of the 0.15 and 0.37 fiber volume fractions presented show good agreement with experimental data and the fiber pressure model when an estimated effective fiber/matrix debond length is used.

Covey, Steven J.↗