Results for “Finite Element Modeling” — Page 71

NASA NTRS2001

Modeling and Analysis of Structural Dynamics for a One-Tenth Scale Model NGST Sunshield

New modeling and analysis techniques have been developed for predicting the dynamic behavior of the Next Generation Space Telescope (NGST) sunshield. The sunshield consists of multiple layers of pretensioned, thin-film membranes supported by deployable booms. Modeling the structural dynamic behavior of the sunshield is a challenging aspect of the problem due to the effects of membrane wrinkling. A finite element model of the sunshield was developed using an approximate engineering approach, the cable network method, to account for membrane wrinkling effects. Ground testing of a one-tenth scale model of the NGST sunshield were carried out to provide data for validating the analytical model. A series of analyses were performed to predict the behavior of the sunshield under the ground test conditions. Modal analyses were performed to predict the frequencies and mode shapes of the test article and transient response analyses were completed to simulate impulse excitation tests. Comparison was made between analytical predictions and test measurements for the dynamic behavior of the sunshield. In general, the results show good agreement with the analytical model correctly predicting the approximate frequency and mode shapes for the significant structural modes.

Johnston, John↗

DOE OSTI2019

HAMT extension for EnergyPlus encompasses moisture sources due to air leakage

As building envelopes become more energy efficient, new buildings are more prone to hygrothermal failure. The combined Heat and Moisture Transfer (HAMT) model in EnergyPlus is able to provide moisture profiles through multiple material layers of an envelope assembly using a coupled, one-dimensional, finite element model simulating the movement and storage of heat and moisture. However, HAMT is currently unable to properly assess the effects of moisture sources inside an assembly to determine impacts on durability, indoor air quality, and energy performance.This paper reviews existing tools for modeling moisture sources due to air leakage and the HAMT model equations in EnergyPlus. A sub-model is then developed that enables air leakage to be applied to a moisture sink/source at an arbitrary location by extending the existing mass transfer equation with a source term that includes water vapor deposit from air flow through an assembly. We demonstrate use of this model, discuss limitations, and make the source code publicly available for external use.

Antretter, Florian↗

DOE OSTI2014

Quench protection analysis of the Mu2e production solenoid

The Muon-to-Electron conversion experiment (Mu2e), under development at Fermilab, seeks to detect direct muon to electron conversion to provide evidence for a process violating muon and electron lepton number conservation that cannot be explained by the Standard Model of particle physics. The Mu2e magnet system consists of three large superconducting solenoids. In case of a quench, the stored magnetic energy is extracted to an external dump circuit. However, because of the fast current decay, a significant fraction of the energy dissipates inside of the cryostat in the coil support shells made of structural aluminum, and in the radiation shield. A 3D finite-element model of the complete cold-mass was created in order to simulate the quench development and understand the role of the quench-back. The simulation results are reported at the normal and non-standard operating conditions.

43 PARTICLE ACCELERATORS↗

NASA NTRS1990

General purpose program to generate compatibility matrix for the integrated force method

An efficient procedure for obtaining the compatibility conditions of finite-element models involves the generation of both field and compatibility conditions from deformation-displacement relations, using (1) the compatibility bandwidth, and (2) the node-determinacy concept. A computer program thus structured will generate sparse and banded compatibility conditions for a structure that is idealized by the finite elements.

Nagabhusanam, J.↗

NASA NTRS2001

Structural Modeling of a Five-Meter Thin Film Inflatable Antenna/Concentrator With Rigidized Support Struts

Dynamic characterization of a non-rigidized thin film inflatable antenna/solar concentrator structure with rigidized composite support struts is described in detail. A two-step finite element modeling approach in MSC/NASTRAN is utilized, consisting of: (1) a nonlinear static pressurization procedure used to obtain the updated stiffness matrix, and (2) a modal "restart" eigen solution that uses the modified stiffness matrix. Unique problems encountered in modeling of this large 5-m lightweight inflatable are identified, including considerable difficulty in obtaining convergence in the nonlinear pressurization solution. It was found that the extremely thin polyimide film material (.001 in or I mil) presents tremendous problems in obtaining a converged solution when internal pressure loading is applied. It was concluded that the ratios of film thickness to other geometric dimensions such as torus cross-sectional and ring diameter and lenticular diameter are the critical parameters for convergence of the pressurization procedure. Comparison of finite element predictions for frequency and mode shapes with experimental results indicated reasonable agreement considering the complexity of the structure, the film-to-air interaction, and the nonlinear material properties of the film. It was also concluded that analysis should be done using different finite element to codes to determine if a more robust and stable solution can be obtained.

Smalley, Kurt B.↗

NASA NTRS2010

Predicting the Dynamic Crushing Response of a Composite Honeycomb Energy Absorber Using Solid-Element-Based Models in LS-DYNA

This paper describes an analytical study that was performed as part of the development of an externally deployable energy absorber (DEA) concept. The concept consists of a composite honeycomb structure that can be stowed until needed to provide energy attenuation during a crash event, much like an external airbag system. One goal of the DEA development project was to generate a robust and reliable Finite Element Model (FEM) of the DEA that could be used to accurately predict its crush response under dynamic loading. The results of dynamic crush tests of 50-, 104-, and 68-cell DEA components are presented, and compared with simulation results from a solid-element FEM. Simulations of the FEM were performed in LS-DYNA(Registered TradeMark) to compare the capabilities of three different material models: MAT 63 (crushable foam), MAT 26 (honeycomb), and MAT 126 (modified honeycomb). These material models are evaluated to determine if they can be used to accurately predict both the uniform crushing and final compaction phases of the DEA for normal and off-axis loading conditions

Jackson, Karen E.↗

DOE OSTI2024

Observation and Modeling of Dynamic Fracture Behaviors of Battery Cell Under Impact Loading Using Enhanced Representative Volume Element Concept

The burgeoning electric automobile industry has increased interest in battery safety. Battery cells experience significant mechanical stress during operation, including the impact of accidents and vibrations from driving. The potential for thermal runaway reactions in battery cells raises safety concerns. Although numerous researchers have defined the dynamic behavior of battery cells and proposed numerical models to describe it, few studies have focused on the high-strain rate mechanical impact phase correlated with the onset of fracture. In this study, we describe the dynamic behavior of pouch battery cells and propose a modeling method to study their mechanical failure under impact situations. Impact tests are conducted at various velocities and heights. To overcome numerical issues commonly encountered under rapid deformation scenarios, a new finite element model is developed based on the representative volume element model. The proposed approach efficiently simulates continuous crack propagation and brittleness behavior during impact by permitting the individual behavior of the cell components. Therefore, engineers can reliably design safer electric vehicle battery cells by measuring the properties of the cell components.

ENERGY STORAGE↗

DOE OSTI2024

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗

NASA NTRS1984

A critical assessment of viscous models of trench topography and corner flow

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

Zhang, J.↗

NASA NTRS1985

A critical assessment of viscous models of trench topography and corner flow

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

Zhang, J.↗

NASA NTRS1989

Incorporation of a progressive failure analysis method in the CSM testbed software system

Analysis of the postbuckling behavior of composite shell structures pose many difficult and challenging problems in the field of structural mechanics. Current analysis methods perform well for most cases in predicting the postbuckling response of undamaged components. To predict component behavior accurately at higher load levels, the analysis must include the effects of local material failures. The CSM testbed software system is a highly modular structural analysis system currently under development at Langley Research Center. One of the primary goals of the CSM testbed is to provide a software environment for the development of advanced structural analysis methods and modern numerical methods which will exploit advanced computer architecture such as parallel-vector processors. Development of a progressive failure analysis method consists of the design and implementation of a processor which will perform the ply-level progressive failure analysis and the development of a geometrically nonlinear analysis procedure which incorporates the progressive failure processor. Regarding the development of the progressive failure processor, two components are required: failure criteria and a degradation model. For the initial implementation, the failure criteria of Hashin will be used. For a matrix failure which typically indicates the development of transverse matrix cracks, the ply properties will be degraded. Work to date includes the design of the progressive failure analysis processor and initial plans for the controlling geometrically nonlinear analysis procedure. The implementation of the progressive failure analysis has begun. Access to the model database and the Hashin failure criteria are completed. Work is in progress on the input/output operations for the processor related data and the finite element model updating procedures. In total the progressive failure processor is approximately one-third complete.

Arenburg, Robert T.↗

NASA NTRS1992

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for a basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness, and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few apriori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, Mordechay↗

NASA NTRS1992

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few a priori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, M.↗

NASA NTRS2001

Development of Reliability Based Life Prediction Methods for Thermal and Environmental Barrier Coatings in Ceramic Matrix Composites

Literature survey related to the EBC/TBC (environmental barrier coating/thermal barrier coating) fife models, failure mechanisms in EBC/TBC and the initial work plan for the proposed EBC/TBC life prediction methods development was developed as well as the finite element model for the thermal/stress analysis of the GRC-developed EBC system was prepared. Technical report for these activities is given in the subsequent sections.

Shah, Ashwin↗

NASA NTRSDate not supplied

Component Characterization of an eVTOL Reference Model for Crashworthiness Studies

Researchers at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) have conducted a series of structural component and seat level tests to improve finite element model (FEM) characterization of a representative vertical take-off and landing (eVTOL) test article developed by NASA. A full-scale dynamic test was conducted on the representative eVTOL test article in November of 2022. The test article represented a high wing, six passenger eVTOL design concept and is referred to as the lift plus cruise (LPC) test article. The full-scale test identified limitations in the analytical models used to predict aircraft structural response, in particular the composite material models did not effectively capture brittle failure of the structure which were measured during dynamic loading. To better understand the mechanism behind the composite material failure mechanisms observed and to improve the FEM, intact sample specimens of the composite airframe structure were recovered from the test article post-test and used in material characterization testing. In addition, the seat configurations used in the LPC test article were further studied using isolated seat and anthropomorphic test device (ATD) drop tower testing. Dynamic compression tests and three-point bend tests, conducted at varied impact speeds, were performed on the recovered frame section specimens. Additional testing was conducted to characterize the material properties of the forming foam, which remained in the frames after fabrication. These tests were used to improve characterization of the damage and failure parameters of the composite material model used in the FE model of the LPC test article. Seat level tests were conducted on the seats used in the LPC test article using acceleration pulses inclusive of current general aviation and rotorcraft certification load levels as well as conditions representative of those measured at the seat base during the LPC test. The structural material models and seat environment models of the LPC test article FEM were calibrated using the generated component test data. The updates made to these models were then integrated into the LPC FEM and simulated in the full-scale test condition. Results demonstrated the effectiveness of component testing to improve predictive capability of composite aerospace structural models within the crash and dynamic loading environments. Demonstration of the LPC FEM response across an accumulation of coupon, component, seat environment, and full-scale test levels provides confidence in the predictive capability of this model for future use in the study of occupant safety within eVTOL relevant crash environments.

Craswhorthiness↗

NASA NTRSDate not supplied

Component Characterization of an eVTOL Reference Model for Crashworthiness Studies

Researchers at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) have conducted a series of structural component and seat level tests to improve finite element model (FEM) characterization of a representative vertical take-off and landing (eVTOL) test article developed by NASA. A full-scale dynamic test was conducted on the representative eVTOL test article in November of 2022. The test article represented a high wing, six passenger eVTOL design concept and is referred to as the lift plus cruise (LPC) test article. The full-scale test identified limitations in the analytical models used to predict aircraft structural response, in particular the composite material models did not effectively capture brittle failure of the structure which were measured during dynamic loading. To better understand the mechanism behind the composite material failure mechanisms observed and to improve the FEM, intact sample specimens of the composite airframe structure were recovered from the test article post-test and used in material characterization testing. In addition, the seat configurations used in the LPC test article were further studied using isolated seat and anthropomorphic test device (ATD) drop tower testing. Dynamic compression tests and three-point bend tests, conducted at varied impact speeds, were performed on the recovered frame section specimens. Additional testing was conducted to characterize the material properties of the forming foam, which remained in the frames after fabrication. These tests were used to improve characterization of the damage and failure parameters of the composite material model used in the FE model of the LPC test article. Seat level tests were conducted on the seats used in the LPC test article using acceleration pulses inclusive of current general aviation and rotorcraft certification load levels as well as conditions representative of those measured at the seat base during the LPC test. The structural material models and seat environment models of the LPC test article FEM were calibrated using the generated component test data. The updates made to these models were then integrated into the LPC FEM and simulated in the full-scale test condition. Results demonstrated the effectiveness of component testing to improve predictive capability of composite aerospace structural models within the crash and dynamic loading environments. Demonstration of the LPC FEM response across an accumulation of coupon, component, seat environment, and full-scale test levels provides confidence in the predictive capability of this model for future use in the study of occupant safety within eVTOL relevant crash environments.

Craswhorthiness↗

DOE OSTI2019

Multilayered triangular and quadrilateral flat shell elements based on the Refined Zigzag Theory

The paper presents a class of C 0 -continuous, flat shell elements based on the Refined Zigzag Theory (RZT) for the analysis of multilayered and curved composite and sandwich structures. In this study, the use of the interdependent interpolation strategy allows eliminating the shear-locking phenomenon and introducing the drilling rotation necessary to complete the set of classical nodal degrees of freedom (three displacements and three rotations). Additional kinematic variables are present in RZT, the zigzag rotations around the in-plane axes that measure the normal distortion typical of multilayered structures. An additional “drilling” zigzag rotation is therefore included among the nodal degrees of freedom in order to properly model curved and built-up structures. A stabilization procedure is adopted to suppress spurious zero-energy modes. A three-node triangular and a four-node quadrilateral flat shell element are formulated with 9 degrees of freedom per node. Example problems involving flat and curved multilayered structures are presented and discussed in order to assess the accuracy and convergence properties of the presented elements. Both static response predictions and free vibrations analyses are considered and the comparison is made with analytic RZT solutions, high-fidelity 3D finite element models and FSDT-based flat shell elements.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Uncertainty Quantification and Comparison of Weld Residual Stress Measurements and Predictions.

In pressurized water reactors, the prevention, detection, and repair of cracks within dissimilar metal welds is essential to ensure proper plant functionality and safety. Weld residual stresses, which are difficult to model and cannot be directly measured, contribute to the formation and growth of cracks due to primary water stress corrosion cracking. Additionally, the uncertainty in weld residual stress measurements and modeling predictions is not well understood, further complicating the prediction of crack evolution. The purpose of this document is to develop methodology to quantify the uncertainty associated with weld residual stress that can be applied to modeling predictions and experimental measurements. Ultimately, the results can be used to assess the current state of uncertainty and to build confidence in both modeling and experimental procedures. The methodology consists of statistically modeling the variation in the weld residual stress profiles using functional data analysis techniques. Uncertainty is quantified using statistical bounds (e.g. confidence and tolerance bounds) constructed with a semi-parametric bootstrap procedure. Such bounds describe the range in which quantities of interest, such as means, are expected to lie as evidenced by the data. The methodology is extended to provide direct comparisons between experimental measurements and modeling predictions by constructing statistical confidence bounds for the average difference between the two quantities. The statistical bounds on the average difference can be used to assess the level of agreement between measurements and predictions. The methodology is applied to experimental measurements of residual stress obtained using two strain relief measurement methods and predictions from seven finite element models developed by different organizations during a round robin study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗