Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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,099 records · Page 61

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗

HydroChrono: An Open-Source Hydrodynamics Package for Project Chrono

In this paper we present the development and verification of HydroChrono, a hydrodynamics package for the Project Chrono physics engine. This package includes the implementation of hydrodynamics equations, the added mass for multibody systems, the development of I/O functions as well as a Python API, and comparison against standard reference cases and other existing tools. HydroChrono provides a flexible, fully open-source solution for simulating wave energy converters (WECs), floating offshore wind turbines (FOWTs) platforms, and other hydrodynamic systems. Here we show, via comparisons with existing tools for benchmark verification cases, that HydroChrono accurately models hydrodynamic forces - making it a useful tool for the design and optimization of these systems. Additionally, the integration of HydroChrono with Project Chrono offers access to finite element modeling capabilities and high-fidelity modelling - with Chrono's existing coupling to CFD and SPH codes. This provides numerical modelers with a multifidelity simulation framework for designing and validating these systems. The development of HydroChrono provides a new, open-source solution for simulating hydrodynamic systems. Its compatibility with other simulation tools enables a more streamlined and efficient design process, advancing the field and providing new opportunities for innovation in this area.

BEM↗

Static shape determination and control for a large space antenna

An integral operator approach is used to derive solutions to static shape determination and control problems associated with large space structures. Problem assumptions include a linear self-adjoint system model, observations and control forces at discrete points, and quadratic performance criteria for the comparison of estimates or control forces. Results are illustrated by simulations with a finite element model of a large space antenna. Modal expansions for terms in the solution algorithms are presented, using modes from the static or associated dynamic model. These expansions provide approximate solutions in the event that a closed form analytical solution to the system boundary value problem is not available.

Weeks, C.↗

Influence of local fiber undulations on the global buckling behavior of filament-wound cylinders

The predicted compressive stiffness and buckling strength of filament-wound cylinders using classical lamination theory is significantly higher than those observed experimentally. This discrepancy is partially influenced by the variation of mechanical properties in the region of fiber undulations. These regions are localized geometric defects intrinsic to the filament-winding, weaving, and braiding processes. In the present work, the average mechanical properties of the fiber undulation region are quantified using modified models of woven-fabric composites to account for the 3-dimensional effects. The mechanical properties thus determined can be incorporated as local element properties into global finite-element models. Preliminary results from large-displacement analyses of filament-wound cylinders are relatively more accurate when fiber undulations are accounted for.

Jensen, D. W.↗

A software tool to design thermal barrier coatings

This paper summarizes work completed for a NASA Phase 1 SBIR program which demonstrated the feasibility of developing a software tool to aid in the design of thermal barrier coating (TBC) systems. Toward this goal, three tasks were undertaken and completed. Task 1 involved the development of a database containing the pertinent thermal and mechanical property data for the top coat, bond coat and substrate materials that comprise a TBC system. Task 2 involved the development of an automated set-up program for generating two dimensional (2D) finite element models of TBC systems. Most importantly, Task 3 involved the generation of a rule base to aid in the design of a TBC system. These rules were based on a factorial design of experiments involving FEM results, and were generated using a Yates analysis. A previous study has indicated the suitability and benefit of applying finite element analysis to perform computer based experiments to decrease but not eliminate physical experiments on TBC's. This program proved feasibility by expanding on these findings by developing a larger knowledge base and developing a procedure to extract rules to aid in TBC design.

Petrus, G.↗

A software tool to design thermal barrier coatings

This paper summarizes work completed for a NASA Phase 1 SBIR program which demonstrated the feasibility of developing a software tool to aid in the design of thermal barrier coating (TBC) systems. Toward this goal, three tasks were undertaken and completed. Task 1 involved the development of a database containing the pertinent thermal and mechanical property data for the top coat, bond coat and substrate materials that comprise a TBC system. Task 2 involved the development of an automated set-up program for generating two dimensional (2D) finite element models of TBC systems. Most importantly, task 3 involved the generation of a rule base to aid in the design of a TBC system. These rules were based on a factorial design of experiments involving FEM results and were generated using a Yates analysis. A previous study had indicated the suitability and benefit of applying finite element analysis to perform computer based experiments to decrease but not eliminate physical experiments on TBC's. This program proved feasibility by expanding on these findings by developing a larger knowledgebase and developing a procedure to extract rules to aid in TBC design.

Petrus, Gregory↗

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗

Fiber shape effects on metal matrix composite behavior

The effects of different fiber shapes on the behavior of metal matrix composites is computationally simulated. A three-dimensional finite element model consisting of a group of nine unidirectional fibers in a three by three unit cell array of a SiC/Ti-15-3 metal matrix composite is used in the analysis. The model is employed to represent five fiber shapes that include a circle, an ellipse, a kidney, and two different cross shapes. The distribution of stresses and the composite material properties, such as moduli, coefficients of thermal expansion, and Poisson's ratios, are obtained from the finite element analysis using the various fiber shapes. Comparisons of these results are used to determine the sensitivity of the composite behavior to the different fiber shapes. In general, fiber dominated properties are not affected by fiber geometry and matrix dominated properties are only moderately affected.

Brown, H. C.↗

Structural dynamics modeling of spent nuclear fuel during hypothetical package drop events

The response of spent nuclear fuel (SNF) to hypothetical package drop events is of particular interest in the scope of spent fuel storage and transportation because of the mechanical shock encountered in such scenarios. Previous testing and modeling by the U.S. Department of Energy has demonstrated that the shock and vibration environment of normal shipping and handling conditions (excluding package drop events) is relatively benign and does not challenge the integrity of spent nuclear fuel. Cask drop events are worth considering because SNF packages are required to withstand free drops onto unyielding surfaces as part of their licensing basis. The acceleration experienced during drop events can be orders of magnitude higher, and thus more advanced models are needed to encompass potential nonlinear behavior of the fuel, such as spacer grid buckling and rod-to-rod impact. This work describes a number of finite element models developed to calculate the response of spent nuclear fuel to various hypothetical drop events that have been validated by package and fuel assembly drop tests conducted in the last decade. Sensitivity of the model response to factors such as package drop orientation, secondary impacts, and irradiated material properties as well as their potential impacts to fuel cladding integrity, was also investigated. Cask drops are not expected as a regular occurrence during SNF transportation, but this work helps raise the understanding of SNF mechanical loads to the point of consistency with the package design requirements.

Kadooka, Kevin↗

Weight estimation of unconventional structures by structural optimization

Automated techniques are presented that are used in structural optimization technology, with emphasis on modifications of finite element models to obtain an optimal material distribution for minimum weight while satisfying the prescribed design requirements. It is anticipated that the future development of computer aided engineering (CAE) system will provide environments where structural analysis, a design optimization, and weight evaluation modules are integrated, sharing a common data base. Structural optimization capabilities obtained by integrating a finite element structural analysis program and a numerical optimization code are developed and applied to two illustrative examples: marine gear housing structural weight minimization and joined wing structures.

Miura, Hirokazu↗

Structural Analysis of Wind Tunnel Blades for the National Transonic Facility

NASA has a significant maintenance challenge with the agency’s aging wind tunnel infrastructure. Wind tunnel blades, such as the ones in the National Transonic Facility at NASA Langley Research Center, have been subjected to damage from impacts and fatigue. Accurate analytical models of these blades are necessary to either recertify replacement blades to be built or for redesign of these blades. In this paper a process for generating structural finite element models of these wind tunnel blades is presented. The model developed in this paper was analyzed under loading conditions from tests conducted on the original blades and the results are compared with legacy strain gage data.

Brian H Mason↗

Measurement of helicon waves with phase contrast imaging on DIII-D – A theoretical feasibility study

A DIII-D high-beta H-mode discharge, with I p = 850 kA, B t = 2.1 T n e = 4 10 19 m –3 , has been designed to validate full wave modeling of helicon waves by optimizing the expected response of the Phase Contrast Imaging diagnostic. Helicon waves have been predicted to have high current drive efficiency off-axis without facing the accessibility issues of lower-hybrid waves. To test these predictions experimentally, DIII-D has recently commissioned a high-power helicon antenna. To confidently predict the behavior of helicon waves in future devices, measurements of their fundamental properties and validation against models will be essential. Phase contrast imaging (PCI) is an absolutely calibrated internal reference interferometer able to measure density fluctuations with radial wavenumbers k R between 1.5 cm –1 and 20 cm –1 . For helicon waves 2 cm –1 < k R < 10 cm –1 is expected, allowing PCI to measure their envelope and wavenumber spectrum. This makes PCI a powerful tool for the validation of state-of-the-art models, like the AORSA full wave code. AORSA is used to compute the density perturbations measured by the PCI with 2D calculations corresponding to 11 different toroidal mode numbers combined to resolve the trajectory of the helicon wave in 3D. This is necessary because the waves travel 120 degrees toroidally from the antenna to the PCI. Here, a cold plasma finite element model (CPFEM) [4] is used to predict propagation through the scrape-off layer. The result of the CPFEM model is connected to AORSA by creating an artificial Gaussian antenna on the last closed flux surface. PCI shows best results for waves with small vertical wavenumbers k z . Modeling the helicon waves for several past DIII-D experiments shows that k z is minimized if the intersection of the helicon and the PCI laser beams occurs in the midplane. For such an optimized scenario the predicted signal level is two orders of magnitude larger than the background density fluctuations arising from broadband turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Crystal plasticity model of residual stress in additive manufacturing using the element elimination and reactivation method

Selective laser melting is receiving increasing interest as an additive manufacturing technique. Residual stresses induced by the large temperature gradients and inhomogeneous cooling process can favour the generation of cracks. In this work, a crystal plasticity finite element model is developed to simulate the formation of residual stresses and to understand the correlation between plastic deformation, grain orientation and residual stresses in the additive manufacturing process. The temperature profile and grain structure from thermal-fluid flow and grain growth simulations are implemented into the crystal plasticity model. An element elimination and reactivation method is proposed to model the melting and solidification and to reinitialize state variables, such as the plastic deformation, in the reactivated elements. The accuracy of this method is judged against previous method based on the stiffness degradation of liquid regions by comparing the plastic deformation as a function of time induced by thermal stresses. The method is used to investigate residual stresses parallel and perpendicular to the laser scan direction, and the correlation with the maximum Schmid factor of the grains along those directions. The magnitude of the residual stress can be predicted as a function of the depth, grain orientation and position with respect to the molten pool. The simulation results are directly comparable to X-ray diffraction experiments and stress–strain curves.

316 stainless steel↗

Thermomechanical design criteria for ceramic-coated surfaces

Some early history of ceramic applications is presented. Finite element modeling of components to determine service and fabrication loads found inelastic behavior and residual stresses to be significant to component life. Inelastic behavior mitigates peak strains but enhances residual strains. Results of furnace, Mach 0.3 burner, and engine tests are discussed and categorized into design criteria (loading, geometry, fabrication, materials, analysis, and testing). These design rules and finite element analyses are brought to bear on two test cases: turboshaft engine seals, and rocket thrust chambers.

Mullen, R. L.↗

Compression failure mechanisms of uni-ply composite plates with a circular cutout

The effect of circular-hole size on the failure mode of uniply graphite-epoxy composite plates is investigated experimentally and analytically for uniaxial compressive loading. The test specimens are sandwiched between polyetherimide plastic for nondestructive evaluations of the uniply failure mechanisms associated with a range of hole sizes. Finite-element modeling based on classical lamination theory is conducted for the corresponding materials and geometries to reproduce the experimental results analytically. The type of compressive failure is found to be a function of hole size, with fiber buckling/kinking at the hole being the dominant failure mechanism for hole diam/plate width ratios exceeding 0.062. The results of the finite-element analysis supported the experimental data for these failure mechanisms and for those corresponding to smaller hole sizes.

Khamseh, A. R.↗

NASTRAN analysis of Tokamak vacuum vessel using interactive graphics

Isoparametric quadrilateral and triangular elements were used to represent the vacuum vessel shell structure. For toroidally symmetric loadings, MPCs were employed across model boundaries and rigid format 24 was invoked. Nonsymmetric loadings required the use of the cyclic symmetry analysis available with rigid format 49. NASTRAN served as an important analysis tool in the Tokamak design effort by providing a reliable means for assessing structural integrity. Interactive graphics were employed in the finite element model generation and in the post-processing of results. It was felt that model generation and checkout with interactive graphics reduced the modelling effort and debugging man-hours significantly.

Miller, A.↗

Dynamic Characterization of an Inflatable Concentrator for Solar Thermal Propulsion

An inflatable structural system that is a technology demonstrator for solar thermal propulsion and other applications is characterized for structural dynamic behavior both experimentally and computationally. The inflatable structure is a pressurized assembly developed for use in orbit to support a Fresnel lens or inflatable lenticular element for focusing sunlight into a solar thermal rocket engine. When the engine temperature reaches a pre-set level, the propellant is injected into the engine, absorbs heat from an exchanger, and is expanded through the nozzle to produce thrust. The inflatable structure is a passively adaptive system in that a regulator and relief valve are utilized to maintain pressure within design limits during the full range of orbital conditions. Modeling and test activities are complicated by the fact that the polyimide film material used for construction of the inflatable is nonlinear, with modulus varying as a function of frequency, temperature, and level of excitation. Modal vibration testing and finite element modeling are described in detail in this paper. The test database is used for validation and modification of the model. This work is highly significant because of the current interest in inflatable structures for space application, and because of the difficulty in accurately modeling such systems.

Leigh, Larry↗

Application of Benchmark Examples to Assess the Single and Mixed-Mode Static Delamination Propagation Capabilities in ANSYS

The application of benchmark examples for the assessment of quasi-static delamination propagation capabilities is demonstrated for ANSYS. The examples are independent of the analysis software used and allow the assessment of the automated delamination propagation in commercial finite element codes based on the virtual crack closure technique (VCCT). The examples selected are based on two-dimensional finite element models of Double Cantilever Beam (DCB), End-Notched Flexure (ENF), Mixed-Mode Bending (MMB) and Single Leg Bending (SLB) specimens. First, the quasi-static benchmark examples were recreated for each specimen using the current implementation of VCCT in ANSYS . Second, the delamination was allowed to propagate under quasi-static loading from its initial location using the automated procedure implemented in the finite element software. Third, the load-displacement relationship from a propagation analysis and the benchmark results were compared, and good agreement could be achieved by selecting the appropriate input parameters. The benchmarking procedure proved valuable by highlighting the issues associated with choosing the input parameters of the particular implementation. Overall the results are encouraging, but further assessment for three-dimensional solid models is required.

Krueger, Ronald↗