Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Method (FEM)”

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 271 records · Page 15

Analysis of hourglass instabilities and control in underintegrated finite element methods

Belytschko et al. (1981, 1984) has developed stabilization methods for the treatment of underintegrated FEM problems; these methods involve the computation of an underintegrated stiffness matrix, which is rank-deficient, and the addition of a stabilization matrix which effectively eliminates the spurious modes. An attempt is presently made to give this a priori stabilization method a mathematical means of support. Attention is also given to an a posteriority stabilization method for hourglass control, in which an approximate solution of the underintegrated system is obtained and then subjected to a special projection in order to eliminate the hourglass modes. A proof is obtained for the convergence of this stabilized underintegrated approximation to the exact solution of a model problem at almost the same rate (as the mesh is refined) as the fully integrated solutions.

Jacquotte, O.-P.↗

Large Area Detection of Microstructural Defects with Multi-Mode Ultrasonic Signals

Cyclic loading or other stresses can lead to development of cracks and crack growth in mechanical structures, leading to eventual failure. While ultrasound imaging can be used for non-destructive testing of such structures, conventional ultrasound techniques are often limited by crack size, density, and areal coverage. An effective characterization of real-world, large-area structures is required at an early damage stage to prevent catastrophic failure and predict remaining life. In this study, a new nonlinear ultrasonic testing (NUT) method is proposed for large-area monitoring of practical structures with arbitrary complexity by using multiple-mode guided-wave ultrasonic signals. The proposed guided-wave NUT technique requires single-element transducers, simple electronics, and a mixed time-frequency domain signal processing. As a proof-of-concept demonstration, numerical simulations and experiments are performed on an A36 carbon steel beam assembly with previously formed microstructural defects that cause nonlinearities in ultrasonic response. The quadratic dependence of the nonlinear wave excitation on the input ultrasonic signal amplitude is shown by numerical simulations, and such a nonlinear ultrasonic response is experimentally observed in the zone with a high density of microstructural defects.

36 MATERIALS SCIENCE↗

Sizing and shape design sensitivity analysis using a hybrid finite element code

The practical (software) implementation of the theoretical FEM-based design sensitivity analysis method of Haug et al. (1986) is described and demonstrated. The fundamental principles of the method are reviewed; the problems involved in applying it to the Engineering Analysis Language (EAL) hybrid FEM code (Whetstone, 1983) are discussed; the formulations for beams, plate bending, and built-up structures are given in detail; programming aspects are considered; and numerical results for sample problems are presented in extensive tables. It is shown that the use of a data-base management system with an FEM such as EAL makes it possible to obtain accurate sizing and shape design sensitivity information using only one program with one data base.

Dopker, Bernhard↗

Improving Computational Efficiency of Mechanical Finite Element Method Simulations for PV Modules: Preprint

In this work we have elucidated the trade off between structural mechanics FEM model accuracy and computation time by employing lower fidelity viscoelastic models. Results indicate that computation time can easily be cut in half while only expecting a potential maximum error of 10 % by considering lower fidelity models. A novel approach to produce the Prony Series fit for viscoelastic characterization has also been presented. Employing this approach, we were able to achieve a further 10 % reduction in computation time without further sacrificing simulation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Augmenting the Structural Dynamics Workflow to Reduce Development Times [Slides]

The structural dynamics workflow at LANL stands to be improved by implementing methods and techniques related to reduced-order finite-element modeling (FEM), pretest predictions, and FEM-test correlation and updating. The benefits of these industry-standard practices are demonstrated in this work. Software used for this project included Ansys 2023 R2 (Mechanical and Spaceclaim), Creo Parametric 8.0, MathCAD Prime 8.0, Matlab R2023b.

42 ENGINEERING↗

An Innovative Structural Mode Selection Methodology: Application for the X-33 Launch Vehicle Finite Element Model

An innovative methodology for determining structural target mode selection and mode selection based on a specific criterion is presented. An effective approach to single out modes which interact with specific locations on a structure has been developed for the X-33 Launch Vehicle Finite Element Model (FEM). We presented Root-Sum-Square (RSS) displacement method computes resultant modal displacement for each mode at selected degrees of freedom (DOF) and sorts to locate modes with highest values. This method was used to determine modes, which most influenced specific locations/points on the X-33 flight vehicle such as avionics control components, aero-surface control actuators, propellant valve and engine points for use in flight control stability analysis and for flight POGO stability analysis. Additionally, the modal RSS method allows for primary or global target vehicle modes to also be identified in an accurate and efficient manner.

Hidalgo, Homero, Jr.↗

Compression After Impact on Honeycomb Core Sandwich Panels with Thin Facesheets, Part 2: Analysis

A two part research study has been completed on the topic of compression after impact (CAI) of thin facesheet honeycomb core sandwich panels. The research has focused on both experiments and analysis in an effort to establish and validate a new understanding of the damage tolerance of these materials. Part 2, the subject of the current paper, is focused on the analysis, which corresponds to the CAI testings described in Part 1. Of interest, are sandwich panels, with aerospace applications, which consist of very thin, woven S2-fiberglass (with MTM45-1 epoxy) facesheets adhered to a Nomex honeycomb core. Two sets of materials, which were identical with the exception of the density of the honeycomb core, were tested in Part 1. The results highlighted the need for analysis methods which taken into account multiple failure modes. A finite element model (FEM) is developed here, in Part 2. A commercial implementation of the Multicontinuum Failure Theory (MCT) for progressive failure analysis (PFA) in composite laminates, Helius:MCT, is included in this model. The inclusion of PFA in the present model provided a new, unique ability to account for multiple failure modes. In addition, significant impact damage detail is included in the model. A sensitivity study, used to assess the effect of each damage parameter on overall analysis results, is included in an appendix. Analysis results are compared to the experimental results for each of the 32 CAI sandwich panel specimens tested to failure. The failure of each specimen is predicted using the high-fidelity, physicsbased analysis model developed here, and the results highlight key improvements in the understanding of honeycomb core sandwich panel CAI failure. Finally, a parametric study highlights the strength benefits compared to mass penalty for various core densities.

Mcquigg, Thomas D.↗

A finite element method for nonlinear forced vibrations of beams

Techniques for defining a finite element model (FEM) for analysis of nonlinear vibrations in beam structures subjected to harmonic excitation are presented. The resulting model covers longitudinal deformation and inertial effects. The nonlinear oscillations of a beam element under forced excitation are modeled by a harmonic force matrix based on first order approximations of the Jacobian elliptic forcing function. Harmonic force and nonlinear stiffness matrices are derived and the nonlinear forced responses of beams are calculated under various boundary conditions. The results of FEM computations for simply-supported and clamped beams show that midplane stretching caused by large deflections increases the nonlinearity. Axially-restrained beams experience only hardening nonlinearity, while axially-free beams have reduced nonlinearity in deformation and inertia and an increase in linearity due to large deflection.

Mei, C.↗

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification↗

A finite-element method for large-amplitude, two-dimensional panel flutter at hypersonic speeds

The nonlinear flutter behavior of a two-dimensional panel in hypersonic flow is investigated analytically. An FEM formulation based unsteady third-order piston theory (Ashley and Zartarian, 1956; McIntosh, 1970) and taking nonlinear structural and aerodynamic phenomena into account is derived; the solution procedure is outlined; and typical results are presented in extensive tables and graphs. A 12-element finite-element solution obtained using an alternative method for linearizing the assumed limit-cycle time function is shown to give predictions in good agreement with classical analytical results for large-amplitude vibration in a vacuum and large-amplitude panel flutter, using linear aerodynamics.

Mei, Chuh↗

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗

Finite Element Flux-Corrected Transport (FEM-FCT) for the Euler and Navier-Stokes equations

A high resolution finite element method for the solution of problems involving high speed compressible flows is presented. The method uses the concepts of flux-corrected transport and is presented in a form which is suitable for implementation on completely unstructured triangular or tetrahedral meshes. Transient and steady state examples are solved to illustrate the performance of the algorithm.

Loehner, Rainald↗

Recent developments in FEM-CFD

The current status of CFD with regard to unstructured grids employing finite element methods and Eulerian frames is reviewed. Algorithms suitable for the computation of large three-dimensional problems involving flow past arbitrary geometries are developed. Adaptive mesh refinement strategy is reviewed, and domain splitting or local time-stepping are briefly addressed. The development of search algorithms of optimal order, variable time-stepping Jacobi smoothers for elliptic problems, and transport concepts for hyperbolics to help achieve good performance for unstructured multigrid processes is discussed. As examples, transient supersonic flow in a channel, regular shock reflection of a wall, viscous flow past a protruberance, potential flow past a cylinder, and Burgers equation are considered.

Loehner, R.↗

Finite element flux-corrected transport (FEM-FCT) for the Euler and Navier-Stokes equations

A high resolution finite element method for the solution of problems involving high speed compressible flows is presented. The method uses the concepts of flux-corrected transport and is presented in a form which is suitable for implementation on completely unstructured triangular or tetrahedral meshes. Transient and steady-state examples are solved to illustrate the performance of the algorithm.

Lohner, Rainald↗

FEM/SINDA: Combining the strengths of NASTRAN, SINDA, I-DEAS, and PATRAN for thermal and structural analysis

This paper describes the interface/integration between FEM/SINDA, a general purpose geometry driven thermal analysis code, and the FEM software: I-DEAS, PATRAN, and NASTRAN. FEM/SINDA brings together the advantages of the finite element method to model arbitrary geometry and anisotropic materials and SINDA's finite difference capability to model thermal properties, loads, and boundary conditions that vary with time or temperature. I-DEAS and PATRAN thermal entities are directly supported since FEM/SINDA uses the nodes of the FEM model as the point at which the temperature is determined. Output from FEM/SINDA (as well as the FEM/SINDA input deck) can be used directly by NASTRAN for structural analysis.

Zarda, P. Richard↗

Solar Cruiser TRAC Boom Development

The TRAC boom studied is a 30 meter, high-aspect ratio and highly nonlinear structural element planned for the NASA Solar Cruiser spacecraft where it would serve as the solar sail’s skeletal system designed to deploy the packaged reflective sail membrane as well as support operational sail tension loads for flight. Key design drivers for this boom are understanding and characterizing load-deformation behavior and quantifying buckling performance–under flight conditions (zero gravity). Knowing that the full-scale solar sail system cannot be feasibly tested for flight conditions under Earth’s gravity, a building block approach is methodically followed to verify finite element analysis (FEA) predictions and to evolve the finite element model (FEM) accordingly, beginning with the lowest level complexity to match global elastic behavior to the highest level complexity to bound theoretical responses under conservative expected flight conditions. Although the testing campaign for the TRAC boom has not developed yet to produce a full-scale test data point that loads the boom to failure, the fidelity of simulations have evolved with success and warrant further interest and work on this high-aspect ratio, but highly stable, deployable structural boom.

Lee Nguyen↗

Arbitrarily shaped dual-stacked patch antennas: A hybrid FEM simulation

A dual-stacked patch antenna is analyzed using a hybrid finite element - boundary integral (FE-BI) method. The metallic patches of the antenna are modeled as perfectly electric conducting (PEC) plates stacked on top of two different dielectric layers. The antenna patches may be of any shape and the lower patch is fed by a coaxial cable from underneath the ground plane or by an aperture coupled microstrip line. The ability of the hybrid FEM technique for the stacked patch antenna characterization will be stressed, and the EM coupling mechanism is also discussed with the aid of the computed near field patterns around the patches.

Gong, Jian↗

Investigation of difficult component effects on finite element model vibration prediction for the Bell AG-1G helicopter. Volume 2: Correlation results

Under the NASA-sponsored DAMVIBS (Design Analysis Methods for VIBrationS) program, a series of ground vibration tests and NASTRAN finite element model (FEM) correlations were conducted on the Bell AH-1G helicopter gunship to investigate the effects of difficult components on the vibration response of the airframe. Previous correlations of the AG-1G showed good agreement between NASTRAN and tests through 15 to 20 Hz, but poor agreement in the higher frequency range of 20 to 30 Hz. Thus, this effort emphasized the higher frequency airframe vibration response correlations and identified areas that need further R and T work. To conduct the investigations, selected difficult components (main rotor pylon, secondary structure, nonstructural doors/panels, landing gear, engine, furl, etc.) were systematically removed to quantify their effects on overall vibratory response of the airframe. The entire effort was planned and documented, and the results reviewed by NASA and industry experts in order to ensure scientific control of the testing, analysis, and correlation exercise. In particular, secondary structure and damping had significant effects on the frequency response of the airframe above 15 Hz. Also, the nonlinear effects of thrust stiffening and elastomer mounts were significant on the low frequency pylon modes below main rotor 1p (5.4 Hz). The results of the NASTRAN FEM correlations are given.

Dompka, R. V.↗