Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Modelling”

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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Magnetic Materials Characterization and Modeling for the Enhanced Design of Magnetic Shielding of Cryomodules in Particle Accelerators

Particle accelerators produce beams of high-energy particles, which are used for both fundamental and applied scientific research and are critical to the development of accelerator driven sub-critical reactor systems. An effective magnetic shield is very important to achieve higher quality factor (Qo) of the cryomodule of a particle accelerator. The allowed value of field inside the cavity due to all external fields (particularly the Earth’s magnetic field) is ~15 mG or less. The goal of this PhD dissertation is to comprehensively study the magnetic properties of commonly used magnetic shielding materials at both cryogenic and room temperatures. This knowledge can be used for the enhanced design of magnetic shields of cryomodes (CM) in particle accelerators. To this end, we first studied the temperature dependent magnetization behavior (M-H curves) of Amumetal and A4K under different annealing and deformation conditions. This characterized the effect of stress or deformation induced during the manufacturing processes and subsequent restoration of high permeability with appropriate heat treatment. Next, an energy based stochastic model for temperature dependent anhysteretic magnetization behavior of ferromagnetic materials was proposed and benchmarked against experimental data. We show that this model is able to simulate and explain the magnetic behavior of as rolled, deformed and annealed amumetal and A4K over a large range of temperatures. The experimental results for permeability are then used in a finite element model (FEM) in COMSOL to evaluate the shielding effectiveness of multiple shield designs at room temperature as well as cryogenic temperature. This work could serve as a guideline for future design, development and fabrication of magnetic shields of CMs.

Sah, Sanjay↗
NASA NTRS2020

Investigating Local Buckling and Plastic Collapse in Multifunctional Sandwich Composite Cores using a Unit Cell Method

The primary goal of developing analytical models of composite sandwich cores has been to obtain simple expressions for elastic and failure response to improve the rapid analysis of composite structures. An earlier investigation by the authors was made to verify the accuracy of widely-used analytical models that calculate equivalent elastic moduli of a uniform homogeneous single layer of material in place of the complex honeycomb core properties in sandwich composites. This transformation of elastic properties permits a coarsening of the finite element model (FEM) of the core to a fraction of the computational cost in simulating the original complex core architecture. In addition, the core was assumed constructed of multi-layered isotropic walls to support potential additional functionality beyond simple load-carrying capability.

Erik Saether↗
NASA NTRS1989

Fatigue delamination onset prediction in tapered composite laminates

Tapered (0 deg) laminates of S2/CE9000 and S2/SP250 glass/epoxies, and IM6/1827I graphite/epoxy were tested in cyclic tension. The specimens usually showed some initial stable delaminations in the tapered region, but these did not affect the stiffness of the specimens, and loading was continued until the specimens either delaminated unstably, or reached 10(exp 6) to 2 x 10(exp 7) million cycles with no unstable delamination. The final unstable delamination originated at the junction of the thin and tapered regions. A finite-element model was developed for the tapered laminate with and without the initial stable delaminations observed in the tests. The analysis showed that for both cases the most likely place for an opening (Mode 1) delamination to originate is at the junction of the taper and thin regions. For each material type, the models were used to calculate the strain energy release rate, G, associated with delaminations originating at that junction and growing either into the thin region or tapered region. For the materials tested, cyclic G(sub Imax) values from DCB tests were used with the maximum strain energy release rates calculated from the finite-element analysis to predict the onset of unstable delamination at the junction as a function of fatigue cycles. The predictions were compared to experimental values of maximum cyclic load as a function of cycles to unstable delamination from fatigue tests in tapered laminates. For the IM6/1827I and S2/SP250 laminates, the predictions agreed very well with the test data. Predicted values for the S2/CE9000 were conservative compared to the test data.

Murri, Gretchen Bostaph↗
NASA NTRS1991

Fatigue delamination onset prediction in unidirectional tapered laminates

Tapered (0 deg) laminates of S2/CE9000 and S2/SP250 glass/epoxies, and IM6/1827I graphite/epoxy were tested in cyclic tension. The specimens usually showed some initial stable delaminations in the tapered region, but these did not affect the stiffness of the specimens, and loading was continued until the specimens either delaminated unstably, or reached 10(exp 6) to 2 x 10 (exp 7) million cycles with no unstable delamination. The final unstable delamination originated at the junction of the thin and tapered regions. A finite-element model was developed for the tapered laminate with and without the initial stable delaminations observed in the tests. The analysis showed that for both cases the most likely place for an opening (Mode I) delamination to originate is at the junction of the taper and thin regions. For each material type, the models were used to calculate the strain energy release rate, G, associated with delaminations originating at that junction and growing either into the thin region or tapered region. For the materials tested, cyclic G(sub Imax) values from DCB tests were used with the maximum strain energy release rates calculated from the finite-element analysis to predict the onset of unstable delamination at the junction as a function of fatigue cycles. The predictions were compared to experimental values of maximum cyclic load as a function of cycles to unstable delamination from fatigue tests in tapered laminates. For the IM6/1827I and S2/SP250 laminates, the predictions agreed very well with the test data. Predicted values for the S2/CE9000 were conservative compared to the test data.

Murri, Gretchen B.↗
NASA NTRS1996

Neutron diffraction measurements and modeling of residual strains in metal matrix composites

Neutron diffraction measurements at room temperature are used to characterize the residual strains in tungsten fiber-reinforced copper matrix, tungsten fiber-reinforced Kanthal matrix, and diamond particulate-reinforced copper matrix composites. Results of finite element modeling are compared with the neutron diffraction data. In tungsten/Kanthal composites, the fibers are in compression, the matrix is in tension, and the thermal residual strains are a strong function of the volume fraction of fibers. In copper matrix composites, the matrix is in tension and the stresses are independent of the volume fraction of tungsten fibers or diamond particles and the assumed stress free temperature because of the low yield strength of the matrix phase.

Saigal, A.↗
NASA NTRS2017

Augmented Method to Improve Thermal Data for the Figure Drift Thermal Distortion Predictions of the JWST OTIS Cryogenic Vacuum Test

The JWST Optical Telescope Element (OTE) assembly is the largest optically stable infrared-optimized telescope currently being manufactured and assembled, and is scheduled for launch in 2018. The JWST OTE, including the 18 segment primary mirror, secondary mirror, and the Aft Optics Subsystem (AOS) are designed to be passively cooled and operate near 45K. These optical elements are supported by a complex composite backplane structure. As a part of the structural distortion model validation efforts, a series of tests are planned during the cryogenic vacuum test of the fully integrated flight hardware at NASA JSC Chamber A. The successful ends to the thermal-distortion phases are heavily dependent on the accurate temperature knowledge of the OTE structural members. However, the current temperature sensor allocations during the cryo-vac test may not have sufficient fidelity to provide accurate knowledge of the temperature distributions within the composite structure. A method based on an inverse distance relationship among the sensors and thermal model nodes was developed to improve the thermal data provided for the nanometer scale WaveFront Error (WFE) predictions. The Linear Distance Weighted Interpolation (LDWI) method was developed to augment the thermal model predictions based on the sparse sensor information. This paper will encompass the development of the LDWI method using the test data from the earlier pathfinder cryo-vac tests, and the results of the notional and as tested WFE predictions from the structural finite element model cases to characterize the accuracies of this LDWI method.

Thermal Modeling↗
DOE OSTI2014

SU-E-T-58: Calculation of Dose Distribution of Accuboost Brachytherapy in Deformable Polyvinil Alcohol Breast Phantom Using Biomechanical Modeling and Monte Carlo Simulation

Purpose: the accuboost is the most modern method of breast brachytherapy that is a boost method in compressed tissue by a mammography unit. the dose distribution in uncompressed tissue, as compressed tissue is important that should be characterized. Methods: In this study, the mechanical behavior of breast in mammography loading, the displacement of breast tissue and the dose distribution in compressed and uncompressed tissue, are investigated. Dosimetry was performed by two dosimeter methods of Monte Carlo simulations using MCNP5 code and thermoluminescence dosimeters. For Monte Carlo simulations, the dose values in cubical lattice were calculated using tally F6. The displacement of the breast elements was simulated by Finite element model and calculated using ABAQUS software, from which the 3D dose distribution in uncompressed tissue was determined. The geometry of the model is constructed from MR images of 6 volunteers. Experimental dosimetery was performed by placing the thermoluminescence dosimeters into the polyvinyl alcohol breast equivalent phantom and on the proximal edge of compression plates to the chest. Results: The results indicate that using the cone applicators would deliver more than 95% of dose to the depth of 5 to 17mm, while round applicator will increase the skin dose. Nodal displacement, in presence of gravity and 60N forces, i.e. in mammography compression, was determined with 43% contraction in the loading direction and 37% expansion in orthogonal orientation. Finally, in comparison of the acquired from thermoluminescence dosimeters with MCNP5, they are consistent with each other in breast phantom and in chest's skin with average different percentage of 13.7±5.7 and 7.7±2.3, respectively. Conclusion: The major advantage of this kind of dosimetry is the ability of 3D dose calculation by FE Modeling. Finally, polyvinyl alcohol is a reliable material as a breast tissue equivalent dosimetric phantom that provides the ability of TLD dosimetry for validation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

PACKAGE IMPACT MODELS AS A PRECURSOR TO CLADDING ANALYSIS

The evaluation of spent nuclear fuel storage casks and transportation packages under impact loading is an important safety topic that is reviewed as part of cask and package certification by the United States Nuclear Regulatory Commission. Explicit dynamic finite element models of full systems are increasingly common in industry for determining structural integrity during hypothetical drop accidents. Full cask and package model results are also used as the loading basis for single fuel pin impact models, which evaluate the response of fuel cladding under drop conditions. In this paper, a simplified package system is evaluated to illustrate several important structural dynamic phenomena, including the effect of gaps between components, the difference in local response at various points on a package during impact, and the effect of modeling various simplified representations of the basket and fuel assemblies. This paper focuses on the package impact analysis, and how loading conditions for a subsequent fuel assembly or fuel cladding analysis can be extracted.

Klymyshyn, Nicholas A.↗
DOE OSTI2025

Micromechanical Surrogate Machine Learning Model for Creep Deformation Modeling

Process variability during the manufacture of gas turbine engine hot section components can significantly affect the material’s resulting microstructure. In casting, for instance, geometric variation within a component (thin sections versus thick sections, radial location) influences cooling rates and the resulting grain size. The high temperature creep response is known to be sensitive to grain size owing to a diffusional creep mechanism which occurs more readily along grain boundaries. Microstructural variation correspondingly drives mechanical behavior which propagates into component scale performance uncertainty. These factors are essential when planning inspection, maintenance, and repair strategies within a reliability framework. These benefits provide opportunities to increase overall energy efficiency through refined margins. Critically, there is an opportunity to bolster existing data-driven reliability models using physics-driven process-structure-property relations. Here we present recent work establishing a framework for evaluating the probabilistic creep performance of high-temperature materials. A novel microstructure-sensitive crystal plasticity finite element model is established that captures both grain boundary and crystallographic deformation effects. The computationally expensive physics model is calibrated using a statistical approach and this high-fidelity model is subsequently used to train a computationally efficient machine learning surrogate model. The surrogate model is essential for sampling a large ensemble of simulated structure-property pair results. The ensemble data are then mined to extract salient trends to be incorporated into a microstructure-sensitive reliability model. The proposed approach represents a novel way to capture microstructure-sensitive trends from physics-based models within a modern reliability framework.

Fernandez-Zelaia, Patxi [ORNL]↗
NASA NTRS1989

Thermal analysis of a hypersonic wing test structure

The three-dimensional finite element modeling techniques developed for the thermal analysis of a hypersonic wing test structure (HWTS) are described. The computed results are compared to measured test data. In addition, the results of a NASA two-dimensional parameter finite difference local thermal model and the results of a contractor two-dimensional lumped parameter finite difference local thermal model will be presented.

Sandlin, Doral R.↗
NASA NTRS1989

Space Shuttle Main Engine structural analysis and data reduction/evaluation. Volume 3A: High pressure oxidizer turbo-pump preburner pump housing stress analysis report

The model generation and structural analysis performed for the High Pressure Oxidizer Turbopump (HPOTP) preburner pump volute housing located on the main pump end of the HPOTP in the space shuttle main engine are summarized. An ANSYS finite element model of the volute housing was built and executed. A static structural analysis was performed on the Engineering Analysis and Data System (EADS) Cray-XMP supercomputer

Shannon, Robert V., Jr.↗
NASA NTRS1996

A Frequency-Domain Substructure System Identification Algorithm

A new frequency-domain system identification algorithm is presented for system identification of substructures, such as payloads to be flown aboard the Space Shuttle. In the vibration test, all interface degrees of freedom where the substructure is connected to the carrier structure are either subjected to active excitation or are supported by a test stand with the reaction forces measured. The measured frequency-response data is used to obtain a linear, viscous-damped model with all interface-degree of freedom entries included. This model can then be used to validate analytical substructure models. This procedure makes it possible to obtain not only the fixed-interface modal data associated with a Craig-Bampton substructure model, but also the data associated with constraint modes. With this proposed algorithm, multiple-boundary-condition tests are not required, and test-stand dynamics is accounted for without requiring a separate modal test or finite element modeling of the test stand. Numerical simulations are used in examining the algorithm's ability to estimate valid reduced-order structural models. The algorithm's performance when frequency-response data covering narrow and broad frequency bandwidths is used as input is explored. Its performance when noise is added to the frequency-response data and the use of different least squares solution techniques are also examined. The identified reduced-order models are also compared for accuracy with other test-analysis models and a formulation for a Craig-Bampton test-analysis model is also presented.

Blades, Eric L.↗
NASA NTRS1998

Analysis and Testing of Plates with Piezoelectric Sensors and Actuators

Piezoelectric material inherently possesses coupling between electrostatics and structural dynamics. Utilizing linear piezoelectric theory results in an intrinsically coupled pair of piezoelectric constitutive equations. One equation describes the direct piezoelectric effect where strains produce an electric field and the other describes the converse effect where an applied electrical field produces strain. The purpose of this study is to compare finite element analysis and experiments of a thin plate with bonded piezoelectric material. Since an isotropic plate in combination with a thin piezoelectric layer constitutes a special case of a laminated composite, the classical laminated plate theory is used in the formulation to accommodated generic laminated composite panels with multiple bonded and embedded piezoelectric layers. Additionally, the von Karman large deflection plate theory is incorporated. The formulation results in laminate constitutive equations that are amiable to the inclusion of the piezoelectric constitutive equations yielding in a fully electro-mechanically coupled composite laminate. Using the finite element formulation, the governing differential equations of motion of a composite laminate with embedded piezoelectric layers are derived. The finite element model not only considers structural degrees of freedom (d.o.f.) but an additional electrical d.o.f. for each piezoelectric layer. Comparison between experiment and numerical prediction is performed by first treating the piezoelectric as a sensor and then again treating it as an actuator. To assess the piezoelectric layer as a sensor, various uniformly distributed pressure loads were simulated in the analysis and the corresponding generated voltages were calculated using both linear and nonlinear finite element analyses. Experiments were carried out by applying the same uniformly distributed loads and measuring the resulting generated voltages and corresponding maximum plate deflections. It is found that a highly nonlinear relationship exists between maximum deflection and voltage versus pressure loading. In order to assess comparisons of predicted and measured piezoelectric actuation, sinusoidal excitation voltages are simulated/applied and maximum deflections are calculated/measured. The maximum deflection as a function of time was determined using the linear finite elements analysis. Good correlation between prediction and measurement was achieved in all cases.

Bevan, Jeffrey S.↗
NASA NTRS2000

Calcaneal loading during walking and running

PURPOSE: This study of the foot uses experimentally measured kinematic and kinetic data with a numerical model to evaluate in vivo calcaneal stresses during walking and running. METHODS: External ground reaction forces (GRF) and kinematic data were measured during walking and running using cineradiography and force plate measurements. A contact-coupled finite element model of the foot was developed to assess the forces acting on the calcaneus during gait. RESULTS: We found that the calculated force-time profiles of the joint contact, ligament, and Achilles tendon forces varied with the time-history curve of the moment about the ankle joint. The model predicted peak talocalcaneal and calcaneocuboid joint loads of 5.4 and 4.2 body weights (BW) during walking and 11.1 and 7.9 BW during running. The maximum predicted Achilles tendon forces were 3.9 and 7.7 BW for walking and running. CONCLUSIONS: Large magnitude forces and calcaneal stresses are generated late in the stance phase, with maximum loads occurring at approximately 70% of the stance phase during walking and at approximately 60% of the stance phase during running, for the gait velocities analyzed. The trajectories of the principal stresses, during both walking and running, corresponded to each other and qualitatively to the calcaneal trabecular architecture.

NASA Center ARC↗
DOE OSTI2021

Modeling lattice rotation fields from discrete crystallographic slip bands in superalloys

Here in this work, we investigate the relationship between an intense slip band (ISB) and the zone of large lattice rotations that forms ahead of the tip of the ISB. We develop a crystal plasticity finite element model of a discrete ISB lying within an oligocrystalline assembly and calculate the local crystalline stress and lattice rotation fields generated by the ISB. The calculations demonstrate that, first, a region of severe lattice rotations, commonly referred to as a microvolume, does not form without the ISB, and second, large amounts of accumulated slip in the ISB are required to enlarge the microvolume to sizes and rotation magnitudes observed experimentally. Ahead of the ISB tip, the quintessential plastic zone always forms, but the atypical microvolume forms when non-concentrated and spatially diffuse slip is activated by the ISB-induced stress field. This result suggests that the detrimental ISB/microvolume pair will likely appear in pairs of crystals in which transmission of the slip from the ISB is severely blocked by the grain boundary, a hypothesis that we verify with a few target cases.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Aeroelastic Optimization Study Based on the X-56A Model

One way to increase the aircraft fuel efficiency is to reduce structural weight while maintaining adequate structural airworthiness, both statically and aeroelastically. A design process which incorporates the object-oriented multidisciplinary design, analysis, and optimization (MDAO) tool and the aeroelastic effects of high fidelity finite element models to characterize the design space was successfully developed and established. This paper presents two multidisciplinary design optimization studies using an object-oriented MDAO tool developed at NASA Armstrong Flight Research Center. The first study demonstrates the use of aeroelastic tailoring concepts to minimize the structural weight while meeting the design requirements including strength, buckling, and flutter. Such an approach exploits the anisotropic capabilities of the fiber composite materials chosen for this analytical exercise with ply stacking sequence. A hybrid and discretization optimization approach improves accuracy and computational efficiency of a global optimization algorithm. The second study presents a flutter mass balancing optimization study for the fabricated flexible wing of the X-56A model since a desired flutter speed band is required for the active flutter suppression demonstration during flight testing. The results of the second study provide guidance to modify the wing design and move the design flutter speeds back into the flight envelope so that the original objective of X-56A flight test can be accomplished successfully. The second case also demonstrates that the object-oriented MDAO tool can handle multiple analytical configurations in a single optimization run.

flutter constraints↗
NASA NTRSDate not supplied

Real-Time State Estimation of Structural Modes for an Aeroelastic Wind Tunnel Model

A method is presented for estimating displacements, velocities, and accelerations of structural modes in generalized coordinates from measured sensor data in real time. Specifically, strain data from conventional strain gauges and fiber optic strain sensors (FOSS) were combined with the strain modes (obtained from a finite element model) in a least-squares estimator to produce structural mode displacement estimates and uncertainties. Likewise, accelerometer data were combined with displacement mode shapes in a second least-squares estimator to produce structural mode acceleration estimates. A Kalman filter was then used to refine the displacement estimates and obtain velocity estimates. The method was applied to the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. The technique was found to be useful for real-time control and system identification applications.

Jared A. Grauer↗
NASA NTRS1989

Scaling effects in the static large deflection response of graphite-epoxy composite beams

Scaling effects in the large deflection response of graphite-epoxy composite beams was investigated. Eight different scale model beams ranging from 1/6 to full-scale were subjected to an eccentric axial compressive load to promote large bending deformations and failures. Beams having laminate stacking sequences including unidirectional, angle ply, cross ply, and quasi-isotropic were tested to examine a wide variety of composite response and failure modes. The model beams were loaded under scaled test conditions until catastrophic failure. Data acquired included load, end displacement, and strain measurements, and qualitative failure measurements. The experimental data is compared to a large rotation beam analysis and a finite element model analysis. Results from the tests indicate that the beam response becomes nonlinear. Failure modes are consistent between scale models within a laminate family, however, a significant scale effect is observed in strength of the scaled beams.

Jackson, Karen E.↗