SEARCH · Engineering Papers

Results for “finite element methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Numerical Implementation of a Multiple-ISV Thermodynamically-Based Work Potential Theory for Modeling Progressive Damage and Failure in Fiber-Reinforced Laminates

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Damage is considered to be the effect of any structural changes in a material that manifest as pre-peak non-linearity in the stress versus strain response. Conversely, failure is taken to be the effect of the evolution of any mechanisms that results in post-peak strain softening. It is assumed that matrix microdamage is the dominant damage mechanism in continuous fiber-reinforced polymer matrix laminates, and its evolution is controlled with a single ISV. Three additional ISVs are introduced to account for failure due to mode I transverse cracking, mode II transverse cracking, and mode I axial failure. Typically, failure evolution (i.e., post-peak strain softening) results in pathologically mesh dependent solutions within a finite element method (FEM) setting. Therefore, consistent character element lengths are introduced into the formulation of the evolution of the three failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs is derived. The theory is implemented into commercial FEM software. Objectivity of total energy dissipated during the failure process, with regards to refinements in the FEM mesh, is demonstrated. The model is also verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared to the experiments.

Pineda, Evan J.↗

NASA NTRS2012

Computational Implementation of a Thermodynamically Based Work Potential Model For Progressive Microdamage and Transverse Cracking in Fiber-Reinforced Laminates

A continuum-level, dual internal state variable, thermodynamically based, work potential model, Schapery Theory, is used capture the effects of two matrix damage mechanisms in a fiber-reinforced laminated composite: microdamage and transverse cracking. Matrix microdamage accrues primarily in the form of shear microcracks between the fibers of the composite. Whereas, larger transverse matrix cracks typically span the thickness of a lamina and run parallel to the fibers. Schapery Theory uses the energy potential required to advance structural changes, associated with the damage mechanisms, to govern damage growth through a set of internal state variables. These state variables are used to quantify the stiffness degradation resulting from damage growth. The transverse and shear stiffness of the lamina are related to the internal state variables through a set of measurable damage functions. Additionally, the damage variables for a given strain state can be calculated from a set of evolution equations. These evolution equations and damage functions are implemented into the finite element method and used to govern the constitutive response of the material points in the model. Additionally, an axial failure criterion is included in the model. The response of a center-notched, buffer strip-stiffened panel subjected to uniaxial tension is investigated and results are compared to experiment.

Pineda, Evan J.↗

NASA NTRS2012

A Thermodynamically-Based Mesh Objective Work Potential Theory for Predicting Intralaminar Progressive Damage and Failure in Fiber-Reinforced Laminates

Pineda, Evan J.↗

NASA NTRS2012

Iterative Nonlinear Tikhonov Algorithm with Constraints for Electromagnetic Tomography

Low frequency electromagnetic tomography such as the capacitance tomography (ECT) has been proposed for monitoring and mass-gauging of gas-liquid two-phase system under microgravity condition in NASA's future long-term space missions. Due to the ill-posed inverse problem of ECT, images reconstructed using conventional linear algorithms often suffer from limitations such as low resolution and blurred edges. Hence, new efficient high resolution nonlinear imaging algorithms are needed for accurate two-phase imaging. The proposed Iterative Nonlinear Tikhonov Regularized Algorithm with Constraints (INTAC) is based on an efficient finite element method (FEM) forward model of quasi-static electromagnetic problem. It iteratively minimizes the discrepancy between FEM simulated and actual measured capacitances by adjusting the reconstructed image using the Tikhonov regularized method. More importantly, it enforces the known permittivity of two phases to the unknown pixels which exceed the reasonable range of permittivity in each iteration. This strategy does not only stabilize the converging process, but also produces sharper images. Simulations show that resolution improvement of over 2 times can be achieved by INTAC with respect to conventional approaches. Strategies to further improve spatial imaging resolution are suggested, as well as techniques to accelerate nonlinear forward model and thus increase the temporal resolution.

Xu, Feng↗

NASA NTRS2012

Multiscale Modeling of Grain Boundaries in ZrB2: Structure, Energetics, and Thermal Resistance

A combination of ab initio, atomistic and finite element methods (FEM) were used to investigate the structures, energetics and lattice thermal conductance of grain boundaries for the ultra high temperature ceramic ZrB2. Atomic models of idealized boundaries were relaxed using density functional theory. Information about bonding across the interfaces was determined from the electron localization function. The Kapitza conductance of larger scale versions of the boundary models were computed using non-equilibrium molecular dynamics. The interfacial thermal parameters together with single crystal thermal conductivities were used as parameters in microstructural computations. FEM meshes were constructed on top of microstructural images. From these computations, the effective thermal conductivity of the polycrystalline structure was determined.

Lawson, John W.↗

NASA NTRS2004

Fracture Analysis of Semi-Elliptical Surface Cracks in Ductile Materials

Accurate life assessment of structural components may require advanced life prediction criteria and methodologies. Structural components often exhibit several different types of defects, among the most prevalent being surface cracks. A semi-elliptical surface crack subjected to monotonic loading will exhibit stable crack growth until the crack has reached a critical size, at which the crack loses stability and fracture ensues (Newman, 2000). The shape and geometry of the flaw are among the most influential factors. When considering simpler crack configurations, such as a through-the-thickness crack, a three-dimensional (3D) geometry may be modeled under the approximation of two-dimensional (2D) plane stress or plane strain. The more complex surface crack is typically modeled numerically with the Finite Element Method (FEM). A semi-elliptical surface crack is illustrated in Figure 1-1.

Daniewicz, S. R.↗

NASA NTRS2014

Calibration of International Space Station (ISS) Node 1 Vibro-Acoustic Model

Reported here is the ability of utilizing the Energy Finite Element Method (E-FEM) to predict the vibro-acoustic sound fields within the International Space Station (ISS) Node 1 and to compare the results with actual measurements of leak sounds made by a one atmosphere to vacuum leak through a small hole in the pressure wall of the Node 1 STA module during its period of storage at Stennis Space Center (SSC). While the E-FEM method represents a reverberant sound field calculation, of importance to this application is the requirement to also handle the direct field effect of the sound generation. It was also important to be able to compute the sound fields in the ultrasonic frequency range. This report demonstrates the capability of this technology as applied to this type of application.

Zhang, Weiguo↗

NASA NTRS2014

Calibration of International Space Station (ISS) Node 1 Vibro-Acoustic Model-Report 2

Reported here is the capability of the Energy Finite Element Method (E-FEM) to predict the vibro-acoustic sound fields within the International Space Station (ISS) Node 1 and to compare the results with simulated leak sounds. A series of electronically generated structural ultrasonic noise sources were created in the pressure wall to emulate leak signals at different locations of the Node 1 STA module during its period of storage at Stennis Space Center (SSC). The exact sound source profiles created within the pressure wall at the source were unknown, but were estimated from the closest sensor measurement. The E-FEM method represents a reverberant sound field calculation, and of importance to this application is the requirement to correctly handle the direct field effect of the sound generation. It was also important to be able to compute the sound energy fields in the ultrasonic frequency range. This report demonstrates the capability of this technology as applied to this type of application.

Zhang, Weiguo↗

NASA NTRS2014

Simulating Damage Due to a Lightning Strike Event: Effects of Temperature Dependent Properties on Interlaminar Damage

A multidirectional, carbon fiber-epoxy, composite panel is subjected to a simulated lightning strike, within a finite element method framework, and the effect of material properties on the failure (delamination) response is investigated through a detailed numerical study. The numerical model of the composite panel consists of individual homogenized plies with user-defined, cohesive interface elements between them. Lightning strikes are simulated as an assumed combination of excessive heat and high pressure loadings. It is observed that the initiation and propagation of lightning-induced delamination is a significant function of the temperature dependency of interfacial fracture toughness. This dependency must be defined properly in order to achieve reliable predictions of the present lightning-induced delamination in the composite panel.

computational mechanics↗

NASA NTRS2015

A Numerical Study on the Effect of Facesheet-Core Disbonds on the Buckling Load of Curved Honeycomb Sandwich Panels

A numerical study on the effect of facesheet-core disbonds on the post-buckling response of curved honeycomb sandwich panels is presented herein. This work was conducted as part of the development of a damage tolerance approach for the next-generation Space Launch System heavy lift vehicle payload fairing. As such, the study utilized full-scale fairing barrel segments as the structure of interest. The panels were composed of carbon fiber reinforced polymer facesheets and aluminum honeycomb core. The panels were analyzed numerically using the finite element method. Facesheet and core nodes in a predetermined circular region were detached to simulate a disbond induced via low-speed impact between the outer mold line facesheet and honeycomb core. Surface-to-surface contact in the disbonded region was invoked to prevent interpenetration of the facesheet and core elements. The diameter of this disbonded region was varied and the effect of the size of the disbond on the post-buckling response was observed. A significant change in the slope of the edge load-deflection response was used to determine the onset of global buckling and corresponding buckling load.

structural failure↗

NASA NTRS2016

Fatigue Analysis of Notched Laminates: A Time-Efficient Macro-Mechanical Approach

A coupled transversely isotropic deformation and damage fatigue model is implemented within the finite element method and was utilized along with a static progressive damage model to predict the fatigue life, stiffness degradation as a function of number of cycles, and post-fatigue tension and compression response of notched, multidirectional laminates. Initially, the material parameters for the fatigue model were obtained utilizing micromechanics simulations and the provided [0], [90] and [plus or minus 45] experimental composite laminate S-N (stress-cycle) data. Within the fatigue damage model, the transverse and shear properties of the plies were degraded with an isotropic scalar damage variable. The damage in the longitudinal (fiber) ply direction was suppressed, and only the strength of the fiber was degraded as a function of fatigue cycles. A maximum strain criterion was used to capture the failure in each element, and once this criterion was satisfied, the longitudinal stiffness of the element was decreased by a factor of 10 (sup 4). The resulting, degraded properties were then used to calculate the new stress state. This procedure was repeated until final failure of the composite laminate was achieved or a specified number of cycles reached. For post-fatigue tension and compression behavior, four internal state variables were used to control the damage and failure. The predictive capability of the above-mentioned approach was assessed by performing blind predictions of the notched multidirectional IM7/977-3 composite laminates response under fatigue and post-fatigue tensile and compressive loading, followed by a recalibration phase. Although three different multidirectional laminates were analyzed in the course of this study, only detailed results (i.e., stiffness degradation and post-fatigue stress-train curves as well as damage evolution states for a single laminate ([30/60/90/minus 30/minus 60] (sub 2s)) are discussed in detail here.

Composite Fatigue analysis↗

NASA NTRS2017

Structural Design Exploration of an Electric Powered Multi-Propulsor Wing Configuration

Advancements in aircraft electric propulsion may enable an expanded operational envelope for electrically powered vehicles compared to their internal combustion engine counterparts. High aspect ratio wings provide additional lift and drag reduction for a proposed multi-propulsor design, however, the challenge is to reduce the weight of wing structures while maintaining adequate structural and aeroelastic margins. Design exploration using a conventional design-and-build philosophy coupled with a finite element method (FEM)-based design of experiments (DOE) strategy are presented to examine high aspect ratio wing structures that have spanwise distributed electric motors. Multiple leading-edge-mounted engine masses presented a challenge to design a wing within acceptable limits for dynamic and aeroelastic stability. Because the first four primary bending eigenmodes of the proposed wing structure are very sensitive to outboard motor placement, safety-of-flight requirements drove the need for multiple spars, rib attachments, and outboard structural reinforcements in the design. Global aeroelasticity became an increasingly important design constraint during the on-going design process, with outboard motor pod flutter ultimately becoming a primary design constraint. Designers successively generated models to examine stress, dynamics, and aeroelasticity concurrently. This research specifically addressed satisfying multi-disciplinary design criteria to generate fluid-structure interaction solution sets, and produced high aspect ratio primary structure designs for the NASA Scalable Convergent Electric Propulsion Technology and Operations Research (SCEPTOR) project in the Aeronautic Research Mission Directorate at NASA. In this paper, a dynamics-driven, quasi-inverse design methodology is presented to address aerodynamic performance goals and structural challenges encountered for the SCEPTOR demonstrator vehicle. These results are compared with a traditional computer aided design based approach.

Moore, James B.↗

NASA NTRS2018

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗

NASA NTRS2018

Analysis of the Influence of Layer Shifting on the Elastic Response and Damage Nucleation and Growth in Woven Composite Laminates

The influence of relative layer shifting on the elastic and damage response of plain weave composite laminates is analyzed using a continuum damage mechanics approach in combination with the finite element method. First, the homogenized properties of the woven composite as a function of the number of layers and of layer shifting are presented. Next, the damage development in various shifting configurations is studied using different damage constitutive models for the matrix and the fiber bundles. It is shown that the impact of layer shifting on both the elastic response and the nonlinear damage response is significant. Most notably, the model captures changes in the damage mechanisms within the woven composite that occur due to layer shifting, resulting in stiffer, more brittle behavior, which has been shown experimentally in the literature. Model results in the linear and nonlinear regimes are shown to be consistent with both an independent analytical model and reported experiments.

layer shifting↗

NASA NTRS2018

Sensitivity Analysis for Multidisciplinary Systems (SAMS)

This report describes the research conducted under an interagency collaboration agreement between the Aerospace Systems Directorate of the Air Force Research Laboratory (AFRL/RQ) and the Computational AeroSciences Branch of NASA Langley (NASA LaRC). Both organizations have a long-term goal of developing a modular computational system for coupling fluids and structures to enable both analysis and optimization of aerospace vehicles. Ultimately, the system should support multiple solvers within the fluid and structure domains to allow the best combination for the task at hand, as well as to allow for institutional preferences of specific software components. Towards this goal, the current research was focused on enhancing the existing modal aeroelastic analysis in the NASA FUN3D (Fully-UNstructured three-dimensional CFD (Computational Fluid Dynamics) code) software (Biedron et al. 2018), as well as developing new aeroelastic analysis and optimization capabilities based on a non-linear finite-element method. The methods and enhancements described in this document pertain to FUN3D Version 13.4.

Biedron, Robert T.↗

NASA NTRS2018

Impedance Eduction in a Duct Using the Linearized Euler Equations

An impedance eduction technique that is based upon a numerical solution to the linearized Euler equations coupled with the minimization of an objective function is presented. The linearized Euler equations are solved numerically using a cubic finite element method and the minimum of the objective function is obtained via a gradient based optimizer. The mean flow is allowed to have a gradient in two coordinate directions and the impedance of the liner may vary arbitrarily over the liner surface. A new regularized wall impedance boundary condition presented by Rienstra which contains the effects of the mean boundary layer thickness is implemented. Impedances are educed for a perforate over honeycomb liner using test data acquired in the Langley Grazing Flow Impedance Tube. The impedance of the liner is educed both with an upstream source and a downstream source. The method is validated by comparing the educed impedance to a benchmark method. The primary conclusion of the study is that the impedances educed for upstream and downstream sources are not well matched and that the primary effect of the boundary thickness is to reduce the resistance of the liner.

Watson, Willie R.↗

NASA NTRS2019

X-57 60kW Permanent Magnet Synchronous Cruise Motor Finite Element Electromagnetic Modeling

The X-57 60kW Permanent Magnet Synchronous Motor for cruise applications was modeled utilizing a two-dimensional electromagnetics simulation software called Finite Element Method Magnets (FEMM, D. Meeker). Through FEMM, the simulated induction and torque characteristics of the X-57 PMSM were obtained. These parameters and other values were compared to actual static laboratory measurements. A three-dimensional electromagnetic model of the X-57 cruise motor was created utilizing OperaFEA (Dassault Systemes SE, Velizy-Villacoublay, France). Torque, RPM, power, resistance, and inductance characteristics were examined along with establishing work to begin examining heat flow and heat dissipation for efficiency purposes.

Marrufo, Marco↗

NASA NTRSDate not supplied

Optimization of a Solver for Computational Materials and Structures Problems on NVIDIA Volta and AMD Instinct GPUs

The Scalable Implementation of Finite Elements by NASA (ScIFEN) is a software package developed to solve complex computational materials and structures problems using the finite element method (FEM). In this paper, we describe optimization techniques to speed up the linear solver computation that occurs within the ScIFEN application. We consider GPUs from two different vendors, NVIDIA and AMD as our target platforms for optimization and highlight differences in performance and optimization techniques. The NVIDIA GPU Volta V100 is used in the Summit system deployed at Oak Ridge National Laboratory, and the new exascale system, Frontier, will be using AMD Radeon Instinct GPU. We evaluated the performance of various optimization techniques on test matrices, ranging in size from100K to 4M, that are representative of ScIFEN applications. The linear solver computation is memory-bound on both GPUs. Our experiments show that on the NVIDIA GPU we obtained up to79%of the theoretical peak bandwidth, while the AMD GPU achieved 59%. Overall, the NVIDIA V100 GPU outperforms the AMD MI 25 GPU1. We observed an overall speedup of up to37X on an NVIDIA V100 compared to an Intel Skylake 12-coremachine. The solver for a 4M degree of freedom system took under 2.5 seconds.

Mohammad Zubair↗