Results for “finite element modeling (FEM)”

NASA NTRS2019

Large Male Anthropomorphic Test Device (ATD) Finite Element Model (FEM) Correlation Improvement

As part of an activity for the Commercial Crew Program (CCP), a NASA Engineering and Safety Center (NESC) team conducted tests of small female, medium male, and large male Hybrid III anthropomorphic test device (ATD) head-neck complexes. The results were compared to simulation results from ATD finite element model (FEM) predictions to determine the validity and fidelity of ATD models for use in certifying CCP vehicles. The NESC was requested to help improve the large male ATD FEM to better represent physical responses. This report contains the outcome of the NESC assessment.

Kelly, Michael J.↗

NASA NTRS1996

A Thermostructural Analysis of a Diboride Composite Leading Edge

In an effort to support the design of zirconium diboride composite leading edges for hypersonic vehicles, a finite element model (FEM) of a prototype leading edge was created and finite element analysis (FEA) was employed to assess its thermal and structural response to aerothermal boundary conditions. Unidirectional material properties for the structural components of the leading edge, a continuous fiber reinforced diboride composite, were computed with COSTAR. These properties agree well with those experimentally measured. To verify the analytical approach taken with COSMOS/M, an independent FEA of one of the leading edge assembly components was also done with COSTAR. Good agreement was obtained between the two codes. Both showed that a unidirectional lay-up had the best margin of safety for a simple loading case. Both located the maximum stress in the same region and ply. The magnitudes agreed within 4 percent. Trajectory based aerothermal heating was then applied to the leading edge assembly FEM created with COSMOS/M to determine steady state temperature response, displacement, stresses, and contact forces due to thermal expansion and thermal strains. Results show that the leading edge stagnation line temperature reached 4700 F. The maximum computed failure index for the laminated composite components peaks at 4.2, and is located at the bolt flange in layer 2 of the side bracket. The temperature gradient in the tip causes a compressive stress of 279 ksi along its width and substantial tensile stresses within its depth.

Kowalski, Tom↗

NASA NTRSDate not supplied

Modal Correlation of Complex Aerospace Joints Using Automated Variable Substitution

A critical task involved with being able to predict flight loads accurately in aerospace finite element models (FEMs) is the prior verification of the FEMs by conducting modal survey testing (MST). Experience comparing dynamic response of initial FEMs to MST data tends to demonstrate that FEMs can have unacceptable accuracy even when best modeling practices are followed. One inherent source of inaccuracy in linear dynamic FEMs is the modeling of nonlinear joints with mechanisms such as spherical bearings. These joints are usually designed to freely translate or rotate under the high levels of loading experienced in flight. Engineers who create linear FEMs conventionally model these joints without any stiffness in the mechanism degrees of freedom to meet this design intent. However, inaccuracy is observed during test validation of these FEMs, which usually relies on low-level force excitation orders of magnitude below flight load levels. This low-level modal test rarely overcomes the joint friction that is present, and thus the mechanism joints are able to react loads. This divide between the test results and the FEM creates a significant challenge to the engineer who is performing the correlation, in that the engineer has no basis for what stiffness value should be used to make the FEM match the test results. A compounding challenge is that complicated built-up aerospace structures commonly have multiple joints through a load path where each joint will “stick” and “slip” at different levels of force input. Explicitly matching the dynamics of a system containing these nonlinear mechanisms would require a nonlinear FEM, which is prohibitively costly for dynamic simulations of most aerospace systems. The objective of this paper is to present a workflow that can efficiently cycle through many iterations of a FEM, allowing a Monte Carlo style examination of the design space to identify candidate stiffness values for nonlinear mechanism joints. The outlined approach is specific to MSC Nastran and utilizes MSC Nastran’s symbolic substitution capabilities, coupled with the IMAT™ and Attune™ software packages developed by ATA Engineering, Inc. (ATA). The workflow is demonstrated with a case study from the correlation effort for The Boeing Company’s Crew Space Transportation (CST)-100 Starliner FEM. Keywords: Correlation, MSC Nastran, IMAT™, Attune™, Finite Element Model, Modal Testing

Correlation↗

NASA NTRSDate not supplied

Modal Test and Model Correlation of NASA Plum Brook Station Mechanical Vibration Facility Head Expander –Lessons Learned from the Perspective of an Early-Career Engineer

In preparation for the Sierra Nevada Corporation’s (SNC) Dream Chaser spacecraft vibration test campaign at the Mechanical Vibration Facility (MVF) at NASA Plum Brook Station (PBS) in Sandusky, Ohio, a test-verified model of MVF is needed in order to be able to perform accurate pretest analysis used for determining response limits and abort levels. MVF was designed to vibration test MPCV Orion and was used to perform the system level vibration test of the European Service Module Structural Test Article (E-STA) in 2016. MVF is comprised of an 18 ft diameter annulus table that is driven with sixteen hydraulic vertical actuator assemblies and four hydraulic horizontal actuator assemblies, which allow it to perform single axis vibration testing in the vertical axis and in each of the two orthogonal horizontal axes without the need for reconfiguring the test article. A head expander for the MVF Table has been designed and built that fills in the center opening providing a continuous flat mounting surface with a maximum diameter of 16.25feet that expands the vibration testing capabilities of MVF. The MVF Table with this head expander will be used during the SNC Dream Chaser spacecraft vibration test campaign. Therefore, a critical element in a test-verified model of the MVF will be a test correlated finite element model (FEM) of the head expander. To obtain this, engineers from the Structural Dynamics Lab (SDL) at NASA Glenn Research Center (GRC) in Cleveland, Ohio performed a modal pretest analysis, conducted a modal test in July 2019, and most recently correlated the head expander finite element model to the modal test data up to 300 Hz. From the initial test preparations to the final delivery of a correlated finite element model, all efforts mentioned were led by the same early-career engineers at NASA GRC. From the viewpoint of an early-career engineer, lessons learned about modal pretest analysis, modal testing, and finite element model correlation of the MVF Table expander head will be presented and discussed. This will include the importance of understanding the limitations of using uncorrelated finite element models in the modal pretest analysis and planning, the importance of orthogonality metrics in judging adequacy and accuracy of test mode shapes, and the importance of having the FEM match the as built hardware in the model correlation effort.

Emma L Pierson↗

NASA NTRS2017

Simulation of Guided-Wave Ultrasound Propagation in Composite Laminates: Benchmark Comparisons of Numerical Codes and Experiment

Ultrasonic wave methods constitute the leading physical mechanism for nondestructive evaluation (NDE) and structural health monitoring (SHM) of solid composite materials, such as carbon fiber reinforced polymer (CFRP) laminates. Computational models of ultrasonic wave excitation, propagation, and scattering in CFRP composites can be extremely valuable in designing practicable NDE and SHM hardware, software, and methodologies that accomplish the desired accuracy, reliability, efficiency, and coverage. The development and application of ultrasonic simulation approaches for composite materials is an active area of research in the field of NDE. This paper presents comparisons of guided wave simulations for CFRP composites implemented using four different simulation codes: the commercial finite element modeling (FEM) packages ABAQUS, ANSYS, and COMSOL, and a custom code executing the Elastodynamic Finite Integration Technique (EFIT). Benchmark comparisons are made between the simulation tools and both experimental laser Doppler vibrometry data and theoretical dispersion curves. A pristine and a delamination type case (Teflon insert in the experimental specimen) is studied. A summary is given of the accuracy of simulation results and the respective computational performance of the four different simulation tools.

Composite↗

NASA NTRS2015

Integration of a Finite Element Model with the DAP Bone Remodeling Model to Characterize Bone Response to Skeletal Loading

NASA's Digital Astronaut Project (DAP) has developed a bone remodeling model that has been validated for predicting volumetric bone mineral density (vBMD) changes of trabecular and cortical bone in the absence of mechanical loading. The model was recently updated to include skeletal loading from exercise and free living activities to maintain healthy bone using a new daily load stimulus (DLS). This new formula was developed based on an extensive review of existing DLS formulas, as discussed in the abstract by Pennline et al. The DLS formula incorporated into the bone remodeling model utilizes strains and stress calculated from finite element model (FEM) of the bone region of interest. The proximal femur was selected for the initial application of the DLS formula, with a specific focus on the femoral neck. METHODS: The FEM was generated from CAD geometry of a femur using de-identified CT data. The femur was meshed using linear tetrahedral elements Figure (1) with higher mesh densities in the femoral neck region, which is the primary region of interest for the initial application of the DLS formula in concert with the DAP bone remodeling model. Nodal loads were applied to the femoral head and the greater trochanter and the base of the femur was held fixed. An L2 norm study was conducted to reduce the length of the femoral shaft without significantly impacting the stresses in the femoral neck. The material properties of the FEM of the proximal femur were separated between cortical and trabecular regions to work with the bone remodeling model. Determining the elements with cortical material properties in the FEM was based off of publicly available CT hip scans [4] that were segmented, cleaned, and overlaid onto the FEM.

Bone Remodeling↗

NASA NTRSDate not supplied

Improving Kinetic Consistency Across Elastic Model Transitions With Quadratic Inequality Constrained Weighted Least Squares (WLSQI)

Flexible body modeling presents a large challenge in the development of simulations to aid in design of flight control systems for launch and landing vehicles. Typically, the flexible body model is not a single continuous model but are rather discrete sets of Linear Time Invariant (LTI) Finite Element Models (FEM) incremented by propellant levels. This introduces the problem of smoothly transitioning modal and physical states of the vehicle when switching from one FEM to the next. This paper introduces a new approach to optimally transition flexible body states with weighted least squares, building off previous methods.

Convex Optimization↗

NASA NTRSDate not supplied

Improving Kinetic Consistency Across Elastic Model Transitions With Quadratic Inequality Constrained Weighted Least Squares (WLSQI)

Flexible body modeling presents a large challenge in the development of simulations to aid in design of flight control systems for launch and landing vehicles. Typically, the flexible body model is not a single continuous model but are rather discrete sets of Linear Time Invariant (LTI) Finite Element Models (FEM) incremented by propellant levels. This introduces the problem of smoothly transitioning modal and physical states of the vehicle when switching from one FEM to the next. This paper introduces a new approach to optimally transition flexible body states with weighted least squares, building off previous methods.

GN&C↗

NASA NTRSDate not supplied

Automated Probabilistic Finite Element Model Calibration Tool Based on Uncertainty Quantification and Machine Learning

Qualification and certification of safety critical parts is a hurdle to the adoption of metallic additively manufactured components for aerospace vehicle applications. Challenges include variability in part properties due to inconsistent defect distribution and microstructure. Understanding of the process through finite element modeling (FEM), and process control through in-situ monitoring, may result in significant improvements; however, solutions useful to manufacturers will require large volumes of data and automated data utilization. Toward this end, a generalizable automated FEM calibration paradigm is developed. This paradigm leverages existing and novel tools from machine learning and uncertainty quantification to enable the automatic calibration of FEMs without requiring prior knowledge of the model performance across input parameter space, including meshing and solver settings, which can require time consuming manual model probing or cause noisy and inconsistent predictions. The result is a probabilistic distribution of calibrated and validated FEM input parameters targeting measured data.

Additive manufacturing model calibration finite el↗

NASA NTRS1994

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗

NASA NTRS1985

A high-precision finite element method for shock-tube calculations

A two-pass explicit scheme is developed in order to exploit some of the capabilities of finite difference modeling (FDM) for finite element modeling (FEM), which offers the opportunity to account for any type of geometry in fluid flow modeling. Features of the first-order upwind and the Lax-Wendroff high precision explicit finite difference algorithms are reviewed. A flux limiter is developed for FEM to serve as an analog for the single limiter function which has been defined for the various FDMs. It is shown that an antidiffusive limiter must be introduced into the weighting function which normally multiplies the time-derivative term in the variational equation. The two-pass scheme which results is demonstrated to be the equivalent of FDMs with five-point support. However, the present scheme is valid only for one-dimensional calculations and linear shape functions for shock tube flow phenomena. Further work is required for its use with nonlinear hyperbolic systems.

Hughes, T. J. R.↗

NASA NTRS2008

Optimized Non-Obstructive Particle Damping (NOPD) Treatment for Composite Honeycomb Structures

Non-Obstructive Particle Damping (NOPD) technology is a passive vibration damping approach whereby metallic or non-metallic particles in spherical or irregular shapes, of heavy or light consistency, and even liquid particles are placed inside cavities or attached to structures by an appropriate means at strategic locations, to absorb vibration energy. The objective of the work described herein is the development of a design optimization procedure and discussion of test results for such a NOPD treatment on honeycomb (HC) composite structures, based on finite element modeling (FEM) analyses, optimization and tests. Modeling and predictions were performed and tests were carried out to correlate the test data with the FEM. The optimization procedure consisted of defining a global objective function, using finite difference methods, to determine the optimal values of the design variables through quadratic linear programming. The optimization process was carried out by targeting the highest dynamic displacements of several vibration modes of the structure and finding an optimal treatment configuration that will minimize them. An optimal design was thus derived and laboratory tests were conducted to evaluate its performance under different vibration environments. Three honeycomb composite beams, with Nomex core and aluminum face sheets, empty (untreated), uniformly treated with NOPD, and optimally treated with NOPD, according to the analytically predicted optimal design configuration, were tested in the laboratory. It is shown that the beam with optimal treatment has the lowest response amplitude. Described below are results of modal vibration tests and FEM analyses from predictions of the modal characteristics of honeycomb beams under zero, 50% uniform treatment and an optimal NOPD treatment design configuration and verification with test data.

Panossian, H.↗

NASA NTRS2022

Accelerance Decoupling: An Approach for Removing the Influence of the Test Stand from the Integrated Modal Test

The main objective for launch vehicle (LV) modal testing is to quantify the LV’s modal properties in the free-free state (post pad separation). However, given the size of most LV systems, free-free testing is a challenge and often not feasible. With this, a test stand, typically the launch pad itself, is introduced as the means of support. This shifts the challenge to developing robust numerical methods for removing the influence of the launch pad from the integrated system modal test. The Space Launch System (SLS) is no exception where the mobile launcher (ML) is used to support the vehicle for the integrated modal test (IMT). For the IMT, it is well understood from pre-test analysis with finite element models (FEMs) of the SLS and SLS coupled to ML that the ML has a significant influence on the SLS modal properties especially in the lower frequency range where the primary SLS bending modes exist. An accelerance decoupling (AD) method has been formulated for the purpose of “subtracting out” the influence of the ML from the IMT results. With AD, the SLS decoupled frequency response functions (FRFs) are directly extracted from the IMT FRFs. The subject approach is aimed to utilize measured data only and achieve a robust FRF decoupling scheme. AD is derived from a widely used coupling technique called “Receptance Coupling” (RC). The AD core equation reverses the RC process and utilizes a pair of auxiliary equations that enable the core equation to be resolved based on measured data only. In AD, the decoupled component FRFs are extracted from the coupled system FRFs with a transformation to remove the contribution of the “subtractive component”. This paper addresses the AD’s operational flexibility to resolve SLS free-free modal properties from coupled system measured data but also the possibility to include data from FEM if there is enough confidence in the FEM or if it is asserted that the effect to the final outcome is reasonable.

Accelerance decoupling↗

NASA NTRS2022

Accelerance Decoupling: An Approach for Removing the Influence of the Test Stand from the Integrated Modal Test

The main objective for launch vehicle (LV) modal testing is to quantify the LV’s modal properties in the free-free state (post pad separation). However, given the size of most LV systems, free-free testing is a challenge and often not feasible. With this, a test stand, typically the launch pad itself, is introduced as the means of support. This shifts the challenge to developing robust numerical methods for removing the influence of the launch pad from the integrated system modal test. The Space Launch System (SLS) is no exception where the mobile launcher (ML) is used to support the vehicle for the integrated modal test (IMT). For the IMT, it is well understood from pre-test analysis with finite element models (FEMs) of the SLS and SLS coupled to ML that the ML has a significant influence on the SLS modal properties especially in the lower frequency range where the primary SLS bending modes exist. An accelerance decoupling (AD) method has been formulated for the purpose of “subtracting out” the influence of the ML from the IMT results. With AD, the SLS decoupled frequency response functions (FRFs) are directly extracted from the IMT FRFs. The subject approach is aimed to utilize measured data only and achieve a robust FRF decoupling scheme. AD is derived from a widely used coupling technique called “Receptance Coupling” (RC). The AD core equation reverses the RC process and utilizes a pair of auxiliary equations that enable the core equation to be resolved based on measured data only. In AD, the decoupled component FRFs are extracted from the coupled system FRFs with a transformation to remove the contribution of the “subtractive component”. This paper addresses the AD’s operational flexibility to resolve SLS free-free modal properties from coupled system measured data but also the possibility to include data from FEM if there is enough confidence in the FEM or if it is asserted that the effect to the final outcome is reasonable.

Accelerance decoupling↗

NASA NTRS2022

Accelerance Decoupling: An Approach for Removing the Influence of the Test Stand from the Integrated Modal Test

The main objective for launch vehicle (LV) modal testing is to quantify the LV’s modal properties in the free-free state (post pad separation). However, given the size of most LV systems, free-free testing is a challenge and often not feasible. With this, a test stand, typically the launch pad itself, is introduced as the means of support. This shifts the challenge to developing robust numerical methods for removing the influence of the launch pad from the integrated system modal test. The Space Launch System (SLS) is no exception where the mobile launcher (ML) is used to support the vehicle for the integrated modal test (IMT). For the IMT, it is well understood from pre-test analysis with finite element models (FEMs) of the SLS and SLS coupled to ML that the ML has a significant influence on the SLS modal properties especially in the lower frequency range where the primary SLS bending modes exist. An accelerance decoupling (AD) method has been formulated for the purpose of “subtracting out” the influence of the ML from the IMT results. With AD, the SLS decoupled frequency response functions (FRFs) are directly extracted from the IMT FRFs. The subject approach is aimed to utilize measured data only and achieve a robust FRF decoupling scheme. AD is derived from a widely used coupling technique called “Receptance Coupling” (RC). The AD core equation reverses the RC process and utilizes a pair of auxiliary equations that enable the core equation to be resolved based on measured data only. In AD, the decoupled component FRFs are extracted from the coupled system FRFs with a transformation to remove the contribution of the “subtractive component”. This paper addresses the AD’s operational flexibility to resolve SLS free-free modal properties from coupled system measured data but also the possibility to include data from FEM if there is enough confidence in the FEM or if it is asserted that the effect to the final outcome is reasonable.

Accelerance decoupling↗

NASA NTRS2007

Build-Up Approach to Updating the Mock Quiet Spike Beam Model

When a new aircraft is designed or a modification is done to an existing aircraft, the aeroelastic properties of the aircraft should be examined to ensure the aircraft is flight worthy. Evaluating the aeroelastic properties of a new or modified aircraft can include performing a variety of analyses, such as modal and flutter analyses. In order to produce accurate results from these analyses, it is imperative to work with finite element models (FEM) that have been validated by or correlated to ground vibration test (GVT) data, Updating an analytical model using measured data is a challenge in the area of structural dynamics. The analytical model update process encompasses a series of optimizations that match analytical frequencies and mode shapes to the measured modal characteristics of structure. In the past, the method used to update a model to test data was "trial and error." This is an inefficient method - running a modal analysis, comparing the analytical results to the GVT data, manually modifying one or more structural parameters (mass, CG, inertia, area, etc.), rerunning the analysis, and comparing the new analytical modal characteristics to the GVT modal data. If the match is close enough (close enough defined by analyst's updating requirements), then the updating process is completed. If the match does not meet updating-requirements, then the parameters are changed again and the process is repeated. Clearly, this manual optimization process is highly inefficient for large FEM's and/or a large number of structural parameters. NASA Dryden Flight Research Center (DFRC) has developed, in-house, a Mode Matching Code that automates the above-mentioned optimization process, DFRC's in-house Mode Matching Code reads mode shapes and frequencies acquired from GVT to create the target model. It also reads the current analytical model, as we11 as the design variables and their upper and lower limits. It performs a modal analysis on this model and modifies it to create an updated model that has similar mode shapes and frequencies as those of the target model. The Mode Matching Code output frequencies and modal assurance criteria (MAC) values that allow for the quantified comparison of the updated model versus the target model. A recent application of this code is the F453 supersonic flight testing platform, NASA DFRC possesses a modified F-15B that is used as a test bed aircraft for supersonic flight experiments. Traditionally, the finite element model of the test article is generated. A GVT is done on the test article ta validate and update its FEM. This FEM is then mated to the F-15B model, which was correlated to GVT data in fall of 2004, A GVT is conducted with the test article mated to the aircraft, and this mated F-15B/ test article FEM is correlated to this final GVT.

Herrera, Claudia Y.↗

NASA NTRS1992

Modal Test/Analysis Correlation of Space Station Structures Using Nonlinear Sensitivity

The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlation. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

Gupta, Viney K.↗

NASA NTRS1992

Modal test/analysis correlation of Space Station structures using nonlinear sensitivity

The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlations. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

Gupta, Viney K.↗