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 253 records · Page 14

An analysis of cross-coupling of a multicomponent jet engine test stand using finite element modeling techniques

A two axis thrust measuring system was analyzed by using a finite a element computer program to determine the sensitivities of the thrust vectoring nozzle system to misalignment of the load cells and applied loads, and the stiffness of the structural members. Three models were evaluated: (1) the basic measuring element and its internal calibration load cells; (2) the basic measuring element and its external load calibration equipment; and (3) the basic measuring element, external calibration load frame and the altitude facility support structure. Alignment of calibration loads was the greatest source of error for multiaxis thrust measuring systems. Uniform increases or decreases in stiffness of the members, which might be caused by the selection of the materials, have little effect on the accuracy of the measurements. It is found that the POLO-FINITE program is a viable tool for designing and analyzing multiaxis thrust measurement systems. The response of the test stand to step inputs that might be encountered with thrust vectoring tests was determined. The dynamic analysis show a potential problem for measuring the dynamic response characteristics of thrust vectoring systems because of the inherently light damping of the test stand.

Schweikhard, W. G.↗

Finite-element modeling of layered, anisotropic composite plates and shells: A review of recent research

Finite element papers published in the open literature on the static bending and free vibration of layered, anisotropic, and composite plates and shells are reviewed. A literature review of large-deflection bending and large-amplitude free oscillations of layered composite plates and shells is also presented. Non-finite element literature is cited for continuity of the discussion.

Reddy, J. N.↗

Finite element modeling of cracked bodies using the Bodner-Partom flow law

The Bodner-Partom flow law which models viscoplastic material behavior has been used to represent two nickel-base superalloys, Gatorized IN100 and Inconel 718 at elevated temperature. Procedures for the determination of the material parameters are presented along with a discussion of the physical significance of each parameter. The material model is then used in finite element computations to evaluate the response of cracked bodies to monotonic, sustained, or cyclic loading. Geometries investigated include the center cracked panel, the compact tension specimen, and the single cracked ring under tension. A Hybrid Experimental Numerical (HEN) procedure has been used to deduce crack growth rates from experimental displacement measurements which are input into finite element computations. The results of several studies conducted over the last several years are summarized.

Nicholas, T.↗

Stress Distribution During Deformation of Polycrystalline Aluminum by Molecular-Dynamics and Finite-Element Modeling

In this paper, a multiscale modelling strategy is used to study the effect of grain-boundary sliding on stress localization in a polycrystalline microstructure with an uneven distribution of grain size. The development of the molecular dynamics (MD) analysis used to interrogate idealized grain microstructures with various types of grain boundaries and the multiscale modelling strategies for modelling large systems of grains is discussed. Both molecular-dynamics and finite-element (FE) simulations for idealized polycrystalline models of identical geometry are presented with the purpose of demonstrating the effectiveness of the adapted finite-element method using cohesive zone models to reproduce grain-boundary sliding and its effect on the stress distribution in a polycrystalline metal. The yield properties of the grain-boundary interface, used in the FE simulations, are extracted from a MD simulation on a bicrystal. The models allow for the study of the load transfer between adjacent grains of very different size through grain-boundary sliding during deformation. A large-scale FE simulation of 100 grains of a typical microstructure is then presented to reveal that the stress distribution due to grain-boundary sliding during uniform tensile strain can lead to stress localization of two to three times the background stress, thus suggesting a significant effect on the failure properties of the metal.

Yamakov, V.↗

Finite-element modelling of multilayer X-ray optics

Multilayer optical elements for hard X-rays are an attractive alternative to crystals whenever high photon flux and moderate energy resolution are required. Prediction of the temperature, strain and stress distribution in the multilayer optics is essential in designing the cooling scheme and optimizing geometrical parameters for multilayer optics. The finite-element analysis (FEA) model of the multilayer optics is a well established tool for doing so. Multilayers used in X-ray optics typically consist of hundreds of periods of two types of materials. The thickness of one period is a few nanometers. Most multilayers are coated on silicon substrates of typical size 60 mm × 60 mm × 100–300 mm. The high aspect ratio between the size of the optics and the thickness of the multilayer (10 7 ) can lead to a huge number of elements for the finite-element model. For instance, meshing by the size of the layers will require more than 10 16 elements, which is an impossible task for present-day computers. Conversely, meshing by the size of the substrate will produce a too high element shape ratio (element geometry width/height > 10 6 ), which causes low solution accuracy; and the number of elements is still very large (10 6 ). In this work, by use of ANSYS layer-functioned elements, a thermal-structural FEA model has been implemented for multilayer X-ray optics. The possible number of layers that can be computed by presently available computers is increased considerably.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A finite element model investigation of ultrasonic array performance for inspecting polycrystalline materials

Microstructural noise has long hindered ultrasonic NDE of polycrystalline materials. In recent years however, arrays have enabled new possibilities to advance ultrasonic inspection of these materials. A Finite Element (FE) model is used to explore the different phenomena caused by grain scattering which may hinder detection of defects by an array. These include multiple scattering and beam deviation due to anisotropy; two aspects of the physics which are often required to be ignored in analytical models due to theoretical assumptions or computational limitations. We rely on a GPU based FE solver, Pogo, to provide fast computation and thereby enable parametric studies. The impact on array detection performance when varying center-frequency and aperture size is investigated. Preliminary results show that there exists an optimum for both the array aperture and frequency for inspection of these materials.

36 MATERIALS SCIENCE↗

Evaluation of Human Spaceflight-Related Tissue Weight Relief Using Whole Body Finite Element Model Simulations

Tissue Weight Relief (TWR) is a physiological condition observed in human spaceflight. It not only impacts the injury biomechanics of soft tissue but also the physiological responses of the cardiovascular system both due to fluid redistribution and the effect of tissue-related transmural pressure on the large venous blood vessels. Understanding the effects of tissue weight relief is especially important because of the role it may play in understanding the cause of Space Associated Neuro-Ocular Syndrome (SANS). SANS can be characterized by a number of ocular changes which reduce visual acuity and SANS related symptoms occur in up to 51% of astronauts. A prevailing theory for the causation of SANS is that of headward (cephalad) fluid shift and a prolonged increase of Intracranial Pressure (ICP) similar to intracranial hypertension, which is not fully supported by the experimental data or astronaut symptom reporting. However, it is still believed that SANS is caused by a pressure change in the eye and the surrounding tissues. It has been proposed that TWR plays a substantial role in affecting internal pressures and fluid shifts in microgravity. In this effort, two whole-body Finite Element (FE) models – Elemance and THUMS – are used to ascertain the microgravity-associated TWR of the musculature surrounding the lower body veins. Elemance and THUMS are physics-based computational models that have been validated and verified for several automotive and domestic applications, and as such, can simulate the relief of soft tissue weight due to changes in the gravitation vector. Specifically, the current effort modeled the transition of the gravitational vector from 1G to 0G, applied across the whole-body model in a supine position. For each 1G to 0G transition simulation, the lower body vein’s transmural pressure-time profile was extracted and averaged around the anterior portion of the thigh muscle. The ascertained transmural pressure changes from 1G to 0G transition are given in Figure 1 for the Elemance and the THUMS FE models. The transmural pressure changes of 10 mmHg and 21 mmHg are in the same order of magnitude as Lu’s value of 44 mmHg. It is to be noted that Lu implemented a 0D to 1D lumped parameter model and the Elemance and THUMS are 3D higher order computational models. This proof-of-concept approach demonstrates that TWR pressure can be adequately estimated with in silico techniques however, further in silico investigations need to be conducted to address the unique contributions to the transmural pressure from each of the computational models.

Finite element modeling↗

Finite element model for brittle fracture and fragmentation

A new computational model for brittle fracture and fragmentation has been developed based on finite element analysis of non-linear elasticity equations. The proposed model propagates the cracks by splitting the mesh nodes alongside the most over-strained edges based on the principal direction of strain tensor. To prevent elements from overlapping and folding under large deformations, robust geometrical constraints using the method of Lagrange multipliers have been incorporated. In conclusion, the model has been applied to 2D simulations of the formation and propagation of cracks in brittle materials, and the fracture and fragmentation of stretched and compressed materials.

97 MATHEMATICS AND COMPUTING↗

Finite Element Modeling of Plastic Deformation During Spin Forming of Aluminum 6061-O

Spin and flow forming are metal deformation techniques in which a disk or tube of material is radially thinned and axially lengthened over a rotating mandrel. A finite-element continuum model (FEM) based on commercial software (DEFORM®) was developed to determine the stress/strain distribution and damage accumulation in a flow-formed, near-net-shape part. Computational validation and experimental verification of the model was leveraged to assess the effectiveness of computer simulations. Quantitative comparisons between experimental and computational results show promise for applying a FEM to the flow forming process. Experimental failure locations correlated well with computed damage gradients, and there was good agreement between the measured and predicted roller forces during forming.

Elizabeth Urig↗

Thermal Analysis of a Finite Element Model in a Radiation Dominated Environment

This paper presents a brief overview of thermal analysis, evaluating the University of Arizona mirror design, for the Next Generation Space Telescope (NGST) Pre-Phase A vehicle concept. Model building begins using Thermal Desktop(Tm), by Cullimore and Ring Technologies, to import a NASTRAN bulk data file from the structural model of the mirror assembly. Using AutoCAD(R) capabilities, additional surfaces are added to simulate the thermal aspects of the problem which, for due reason, are not part of the structural model. Surfaces are then available to accept thermophysical and thermo-optical properties. Thermal Desktop(Tm) calculates radiation conductors using Monte Carlo simulations. Then Thermal Desktop(Tm) generates the SINDA/Fluint input file having a one-to-one correspondence with the NASTRAN node and element definitions. A model is now available to evaluate the mirror design in the radiation dominated environment conduct parametric trade studies of the thermal design, and provide temperatures to the finite element structural model.

Page, Arhur T.↗

Thermal Analysis of a Finite Element Model in a Radiation Dominated Environment

This paper presents a brief overview of thermal analysis, evaluating the University of Arizona mirror design, for the Next Generation Space Telescope (NGST) Pre-Phase A vehicle concept. Model building begins using Thermal Desktop(TM), by Cullimore and Ring Technologies, to import a NASTRAN bulk data file from the structural model of the mirror assembly. Using AutoCAD(R) capabilities, additional surfaces are added to simulate the thermal aspects of the problem which, for due reason, are not part of the structural model. Surfaces are then available to accept thermophysical and thermo-optical properties. Thermal Desktop(TM) calculates radiation conductors using Monte Carlo simulations. Then Thermal Desktop(TM) generates the SINDA input file having a one-to-one correspondence with the NASTRAN node and element definitions. A model is now available to evaluate the mirror design in the radiation dominated environment, conduct parametric trade studies of the thermal design, and provide temperatures to the finite element structural model.

Page, Arthur T.↗

An enriched finite element model with q-refinement for radiative boundary layers in glass cooling

Radiative cooling in glass manufacturing is simulated using the partition of unity finite element method. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary simplified P{sub 1} approximation for the radiation in non-grey semitransparent media. To integrate the coupled equations in time we consider a linearly implicit scheme in the finite element framework. A class of hyperbolic enrichment functions is proposed to resolve boundary layers near the enclosure walls. Using an industrial electromagnetic spectrum, the proposed method shows an immense reduction in the number of degrees of freedom required to achieve a certain accuracy compared to the conventional h-version finite element method. Furthermore the method shows a stable behaviour in treating the boundary layers which is shown by studying the solution close to the domain boundaries. The time integration choice is essential to implement a q-refinement procedure introduced in the current study. The enrichment is refined with respect to the steepness of the solution gradient near the domain boundary in the first few time steps and is shown to lead to a further significant reduction on top of what is already achieved with the enrichment. The performance of the proposed method is analysed for glass annealing in two enclosures where the simplified P{sub 1} approximation solution with the partition of unity method, the conventional finite element method and the finite difference method are compared to each other and to the full radiative heat transfer as well as the canonical Rosseland model.

97 MATHEMATICS AND COMPUTING↗

Broadband Calibration of Acoustic Emission and Ultrasonic Sensors from Generalized Ray Theory and Finite Element Models

This paper describes a calibration method for acoustic emission and ultrasonic sensors that is effective from 1 kHz to 1 MHz. The method combines generalized ray theory and finite element analysis to model wave propagation at higher and lower frequencies, respectively. A ball impact is used as a calibration source, a thick aluminum plate is used as the test block, and hot glue is used as the couplant. We demonstrate this method on five commercial piezoelectric sensors: Physical Acoustics (PAC) R15a, PAC WSa, Panametrics V101, Panametrics V103, and Valpey-Fisher Pinducer. Our calibration results show that reflections and other wave phases can be more clearly identified with the less-resonant Panametrics sensors. The PAC sensors have the greatest sensitivity and are able to detect surface normal displacements at least down to 1 pm amplitude in the 100s of kHz frequency band. Aperture effect is minimized by the small size of the Pinducer. Our method focuses on the amplitude response of the sensors (phase is ignored) and extends the calibration to a frequency band that is lower than typical analyses. Low frequency information is useful for determining the seismic moment of a seismic source (analogous to the magnitude of an earthquake) and can increase the amount of information acquired in a single recording.

42 ENGINEERING↗

FINITE ELEMENT MODEL MESH REFINEMENT EFFECTS ON QUALIFICATION OF NUCLEAR GRADE GRAPHITE COMPONENTS

The American Society of Mechanical Engineers (ASME) provides the full and simplified design-by-analysis probabilistic assessments for determining acceptance of nuclear grade graphite core components. The assessments can be characterized by three parts: (1) a component stress distribution, often determined by a finite element (FE) model; (2) a Weibull probability density function (pdf) that characterizes the experimental tensile strength distribution; and (3) the post-processor, which combines the FE model and the Weibull strength distribution in accordance with the full and simplified assessments to determine component acceptance. It is known that the level of mesh refinement in FE models can affect the modeled component’s calculated stress distribution. Depending on the component geometry, the stress distribution may converge with sufficient refinement. It was previously unknown whether the acceptance decision resulting from the full and simplified assessments might change even with sufficient mesh refinement. This study explores that question using experimental strength results for a dog-bone geometry for two graphite grades, IG-110 and PCEA. The simplified assessment has two criteria that must be met, the first limits the combined membrane stress by the allowable stress and the second limits the peak equivalent stress by the allowable stress scaled by the ratio of flexural to tensile strength. In the application of the simplified assessment, convergence of the peak equivalent stress required extreme mesh refinement, however, the acceptance decision was not affected. It is hypothesized that more complex geometries with stress concentrations may present mesh refinement effects on the simplified assessment acceptance decision. Mesh refinement did affect the acceptance decision in the full assessment for the applied pressure loadings in this study. This work suggests component stress distribution convergence is not a sufficient criteria for POF convergence in the full assessment and that mesh refinement should continue until the POF has converged, especially where the resulting POF is bordering the SRC acceptable POF limit.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A 2-D Interface Element for Coupled Analysis of Independently Modeled 3-D Finite Element Subdomains

Over the past few years, the development of the interface technology has provided an analysis framework for embedding detailed finite element models within finite element models which are less refined. This development has enabled the use of cascading substructure domains without the constraint of coincident nodes along substructure boundaries. The approach used for the interface element is based on an alternate variational principle often used in deriving hybrid finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise to a non-positive definite system which causes difficulties with many of the equation solvers in general-purpose finite element codes. Hence the global system of equations is generally solved using, a decomposition procedure with pivoting. The research reported to-date for the interface element includes the one-dimensional line interface element and two-dimensional surface interface element. Several large-scale simulations, including geometrically nonlinear problems, have been reported using the one-dimensional interface element technology; however, only limited applications are available for the surface interface element. In the applications reported to-date, the geometry of the interfaced domains exactly match each other even though the spatial discretization within each domain may be different. As such, the spatial modeling of each domain, the interface elements and the assembled system is still laborious. The present research is focused on developing a rapid modeling procedure based on a parametric interface representation of independently defined subdomains which are also independently discretized.

Kandil, Osama A.↗