Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element problems”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

An inverse problem approach for structural damage detection - Finite element model refinement

In this work, a methodology for incorporating measured modal data into an existing refined finite element model is examined with the objective of detecting and locating structural damage. This same algorithm is also useful in terms of finite element model refinement. The algorithm is based on the partial inverse problem, in that only partial spectral information is required. The technique utilizes a symmetric eigenstructure assignment algorithm to perform the partial spectral assignment. Algorithms to enhance mode shape assignability and to preserve sparsity in the updated model are developed. The sparsity preservation is of particular importance when considering damage detection in truss-like structures. Several examples are presented which highlight the key points made within the paper.

Zimmerman, D. C.↗
NASA NTRS1993

Massively parallel computation of RCS with finite elements

One of the promising combinations of finite element approaches for scattering problems uses Whitney edge elements, spherical vector wave-absorbing boundary conditions, and bi-conjugate gradient solution for the frequency-domain near field. Each of these approaches may be criticized. Low-order elements require high mesh density, but also result in fast, reliable iterative convergence. Spherical wave-absorbing boundary conditions require additional space to be meshed beyond the most minimal near-space region, but result in fully sparse, symmetric matrices which keep storage and solution times low. Iterative solution is somewhat unpredictable and unfriendly to multiple right-hand sides, yet we find it to be uniformly fast on large problems to date, given the other two approaches. Implementation of these approaches on a distributed memory, message passing machine yields huge dividends, as full scalability to the largest machines appears assured and iterative solution times are well-behaved for large problems. We present times and solutions for computed RCS for a conducting cube and composite permeability/conducting sphere on the Intel ipsc860 with up to 16 processors solving over 200,000 unknowns. We estimate problems of approximately 10 million unknowns, encompassing 1000 cubic wavelengths, may be attempted on a currently available 512 processor machine, but would be exceedingly tedious to prepare. The most severe bottlenecks are due to the slow rate of mesh generation on non-parallel machines and the large transfer time from such a machine to the parallel processor. One solution, in progress, is to create and then distribute a coarse mesh among the processors, followed by systematic refinement within each processor. Elimination of redundant node definitions at the mesh-partition surfaces, snap-to-surface post processing of the resulting mesh for good modelling of curved surfaces, and load-balancing redistribution of new elements after the refinement are auxiliary steps expected to result in a robust low i/o system for very large finite element problems.

Parker, Jay↗
NASA NTRS2018

A Higher-Order Space-Time Finite-Element Method for Moving-Body and Fluid-Structure Interaction Problems

We present a high-order finite-element method for moving body and fluid/structure interaction problems. Our solution strategy is based on a space-time discontinuous Galerkin (DG) spectral-element discretization which extends to arbitrary order of accuracy. The space-time DG discretization is a natural choice for moving body and fluid-structure interaction problems as moving surfaces are incorporated simply by considering curved space-time elements whose space-time faces align with the moving body. We present a discontinuous-Galerkin in time discretization for six-degree of motion modeling of rigid bodies, and a continuous-Galerkin discretization for equations of linear elasticity to generate curved space-time meshes. Numerical results for several simple 2D test cases are presented in order to verify the implementation of the different models. Finally we present a preliminary dynamic simulation of a parachute.

High-Order↗
NASA NTRS1984

Algorithms for elasto-plastic-creep postbuckling

This paper considers the development of an improved constrained time stepping scheme which can efficiently and stably handle the pre-post-buckling behavior of general structure subject to high temperature environments. Due to the generality of the scheme, the combined influence of elastic-plastic behavior can be handled in addition to time dependent creep effects. This includes structural problems exhibiting indefinite tangent properties. To illustrate the capability of the procedure, several benchmark problems employing finite element analyses are presented. These demonstrate the numerical efficiency and stability of the scheme. Additionally, the potential influence of complex creep histories on the buckling characteristics is considered.

Padovan, J.↗
NASA NTRS1985

Automatic crack propagation tracking

A finite element based approach to fully automatic crack propagation tracking is presented. The procedure presented combines fully automatic mesh generation with linear fracture mechanics techniques in a geometrically based finite element code capable of automatically tracking cracks in two-dimensional domains. The automatic mesh generator employs the modified-quadtree technique. Crack propagation increment and direction are predicted using a modified maximum dilatational strain energy density criterion employing the numerical results obtained by meshes of quadratic displacement and singular crack tip finite elements. Example problems are included to demonstrate the procedure.

Shephard, M. S.↗
NASA NTRS1986

Adaptive finite element methods for compressible flow problems

Some recent work on adaptive FEMs for solving transient Euler equations in two-dimensional domains is summarized. The formulation of an FEM model of the Euler equations is shown, and the application of the adaptive strategies to data management schemes is addressed. Sample numerical results from the application of the model and strategies to the flow over a step and to transient cases are given.

Oden, J. T.↗
NASA NTRS1990

Accuracy of the QUAD4 thick shell element

The accuracy of the relatively new QUAD4 thick shell element is assessed via comparison with a theoretical solution for thick homogeneous and honeycomb flat simply supported plates under the action of a uniform pressure load. The theoretical thick plate solution is based on the theory developed by Reissner and includes the effects of transverse shear flexibility which are not included in the thin plate solutions based on Kirchoff plate theory. In addition, the QUAD4 is assessed using a set of finite element test problems developed by the MacNeal-Schwendler Corp. (MSC). Comparison of the COSMIC QUAD4 element as well as those from MSC and Universal Analytics, Inc. (UAI) for these test problems is presented. The current COSMIC QUAD4 element is shown to have excellent comparison with both the theoretical solutions and also those from the two commercial versions of NASTRAN that it was compared to.

Case, William R.↗
NASA NTRS1992

Accuracy of the TRIA3 thick shell element

The accuracy of the new TRIA3 thick shell element is assessed via comparison with a theoretical solution for thick homogeneous and honeycomb flat simply supported plates under the action of a uniform pressure load. The theoretical thick plate solution is based on the theory developed by Reissner and includes the effects of transverse shear flexibility which are not included in the thin plate solutions based on Kirchoff plate theory. In addition, the TRIA3 is assessed using a set of finite element test problems developed by the MacNeal-Schwendler Corp. (MSC). Comparison of the COSMIC TRIA3 element as well as those from MSC and Universal Analytics Inc. (UAI) for these problems is presented. The current COSMIC TRIA3 element is shown to have excellent comparison with both the theoretical solutions and also those from the two commercial versions of NASTRAN with which it is compared.

Case, William R.↗
NASA NTRS2004

A Discontinuous Galerkin Method for Parabolic Problems with Modified hp-Finite Element Approximation Technique

A recent paper is generalized to a case where the spatial region is taken in R(sup 3). The region is assumed to be a thin body, such as a panel on the wing or fuselage of an aerospace vehicle. The traditional h- as well as hp-finite element methods are applied to the surface defined in the x - y variables, while, through the thickness, the technique of the p-element is employed. Time and spatial discretization scheme based upon an assumption of certain weak singularity of double vertical line u(sub t) double vertical line 2, is used to derive an optimal a priori error estimate for the current method.

Kaneko, Hideaki↗
NASA NTRS2016

Tangle-Free Finite Element Mesh Motion for Ablation Problems

Mesh motion is the process by which a computational domain is updated in time to reflect physical changes in the material the domain represents. Such a technique is needed in the study of the thermal response of ablative materials, which erode when strong heating is applied to the boundary. Traditionally, the thermal solver is coupled with a linear elastic or biharmonic system whose sole purpose is to update mesh node locations in response to altering boundary heating. Simple mesh motion algorithms rely on boundary surface normals. In such schemes, evolution in time will eventually cause the mesh to intersect and "tangle" with itself, causing failure. Furthermore, such schemes are greatly limited in the problems geometries on which they will be successful. This paper presents a comprehensive and sophisticated scheme that tailors the directions of motion based on context. By choosing directions for each node smartly, the inevitable tangle can be completely avoided and mesh motion on complex geometries can be modeled accurately.

Droba, Justin↗
NASA NTRS2016

Tangle-Free Finite Element Mesh Motion for Ablation Problems

In numerical simulations involving boundaries that evolve in time, the primary challenge is updating the computational mesh to reflect the physical changes in the domain. In particular, the fundamental objective for any such \mesh motion" scheme is to maintain mesh quality and suppress unphysical geometric anamolies and artifacts. External to a physical process of interest, mesh motion is an added component that determines the specifics of how to move the mesh given certain limited information from the main system. This paper develops a set of boundary conditions designed to eliminate tangling and internal collision within the context of PDE-based mesh motion (linear elasticity). These boundary conditions are developed for two- and three-dimensional meshes. The paper presents detailed algorithms for commonly occuring topological scenarios and explains how to apply them appropriately. Notably, the techniques discussed herein make use of none of the specifics of any particular formulation of mesh motion and thus are more broadly applicable. The two-dimensional algorithms are validated by an extensive verification procedure. Finally, many examples of diverse geometries in both two- and three-dimensions are shown to showcase the capabilities of the tangle-free boundary conditions.

Droba, Justin↗
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRS1990

Finite element analysis of structural components using viscoplastic models with application to a cowl lip problem

The viability of advanced viscoplastic models for nonlinear finite element analyses of structural components is investigated. Several uniaxial and a multiaxial problem are analyzed using the finite element implementation of Freed's viscoplastic model. Good agreement between the experimental and calculated uniaxial results validates the finite element implementation and gives confidence to apply it to more complex multiaxial problems. A comparison of results for a sample structural component (the cowl lip of a hypersonic engine inlet) with the earlier elastic, elastic-plastic, and elastic-plastic-creep analyses available in the literature shows that the elastic-viscoplastic analyses yield more reasonable stress and strain distributions. Finally, the versatility of the finite-element-based solution technology presented herein is demonstrated by applying it to another viscoplastic model.

Arya, V. K.↗