Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary finite element”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Why does near ridge extensional seismicity occur primarily in the Indian Ocean?

It is argued that though thermoelastic stresses provide a low level background in all plates, the data favoring their contributing significantly to the stress field and seismicity in the young oceanic lithosphere may be interpreted in terms of stresses resulting from individual plate geometry and local boundary effects. The dramatic concentration of extensional seismicity in the Central Indian Ocean region is shown to be consistent with finite element results for the intraplate stress incorporating the effects of the Himalayan collision and the various subduction zones. Most of the data for both ridge-parallel extension and depth stratification are provided by earthquakes in this area, and it is suggested that these effects may be due more to the regional stress.

Stein, Seth↗
NASA NTRS1987

Design modeling for shape optimization

Some important aspects of design modeling for shape optimization are discussed for both stamped sheet metal components and cast solid components. For stamped components the basis for the modeling approach is a boundary design function. Design parameters control the shape of 2-D regions. For more complex, folded plate components, the 2-D regions can be assembled using translation and rotation operations. The analysis model is automatically created using a mesh generation procedure requiring only boundary data. For less complex solid components, it was found that this approach is not suitable. For these structures, the finite element models are typically created using very sophisticated graphical modeling systems. A new approach which overlays a parameterized surface design model on an existing analysis model is described. To summarize, the future needs for solid shape design is described in terms of an extension of the previously described 2-D capability.

Botkin, M. E.↗
NASA NTRS1988

Thermo-viscoelastic analysis of composite materials, volume 1

Advanced composite materials, especially graphite/epoxy, are being applied to aircraft structures in order to improve performance and save weight. An important consideration in composite design is the residual strength of a structure containing holes, delaminations, or interlaminar damage when subjected to compressive loads. Recent studies have revealed the importance of viscoelastic effects in polymer-based composites. The viscoelastic effect is particularly significant at elevated temperature/moisture conditions since the matrix material is strongly affected by the environment. The solution of viscoelastic problems in composites was limited to special cases which can be solved by classical lamination theory. A finite element procedure is presented for calculating time-dependent stresses and strains in composite structures with general configurations and complicated boundary conditions. Using this procedure the in-plane and interlaminar stress distributions and histories in notched and unnotched composites were obtained for mechanical and thermal loads. Both two-dimensional and three-dimensional viscoelastic problems are analyzed. The effects of layup orientation and load spectrum on creep response and stress relaxation were also studied.

Lin, K. Y.↗
NASA NTRS1992

Numerical simulation of aerothermal loads in hypersonic engine inlets due to shock impingement

The effect of shock impingement on an axial corner simulating the inlet of a hypersonic vehicle engine is modeled using a finite-difference procedure. A three-dimensional dynamic grid adaptation procedure is utilized to move the grids to regions with strong flow gradients. The adaptation procedure uses a grid relocation stencil that is valid at both the interior and boundary points of the finite-difference grid. A linear combination of spatial derivatives of specific flow variables, calculated with finite-element interpolation functions, are used as adaptation measures. This computational procedure is used to study laminar and turbulent Mach 6 flows in the axial corner. The description of flow physics and qualitative measures of heat transfer distributions on cowl and strut surfaces obtained from the analysis are compared with experimental observations. Conclusions are drawn regarding the capability of the numerical scheme for enhanced modeling of high-speed compressible flows.

Ramakrishnan, R.↗
NASA NTRS1993

FEM/SINDA: Combining the strengths of NASTRAN, SINDA, I-DEAS, and PATRAN for thermal and structural analysis

This paper describes the interface/integration between FEM/SINDA, a general purpose geometry driven thermal analysis code, and the FEM software: I-DEAS, PATRAN, and NASTRAN. FEM/SINDA brings together the advantages of the finite element method to model arbitrary geometry and anisotropic materials and SINDA's finite difference capability to model thermal properties, loads, and boundary conditions that vary with time or temperature. I-DEAS and PATRAN thermal entities are directly supported since FEM/SINDA uses the nodes of the FEM model as the point at which the temperature is determined. Output from FEM/SINDA (as well as the FEM/SINDA input deck) can be used directly by NASTRAN for structural analysis.

Zarda, P. Richard↗
NASA NTRS1994

Radiation and scattering from cylindrically conformal printed antennas

Microstrip patch antennas offer considerable advantages in terms of weight, aerodynamic drag, cost, flexibility, and observables over more conventional protruding antennas. These flat patch antennas were first proposed over thirty years ago by Deschamps in the United States and Gutton and Baisinot in France. Such antennas have been analyzed and developed for planar as well as curved platforms. However, the methods used in these designs employ gross approximations, suffer from extreme computational burden, or require expensive physical experiments. The goal of this thesis is to develop accurate and efficient numerical modeling techniques which represent actual antenna structures mounted on curved surfaces with a high degree of fidelity. In this thesis, the finite element method is extended to cavity-backed conformal antenna arrays embedded in a circular, metallic, infinite cylinder. Both the boundary integral and absorbing boundary mesh closure conditions will be used for terminating the mesh. These two approaches will be contrasted and used to study the scattering and radiation behavior of several useful antenna configurations. An important feature of this study will be to examine the effect of curvature and cavity size on the scattering and radiation properties of wraparound conformal antenna arrays.

Kempel, Leo C.↗
NASA NTRSDate not supplied

Test and Analysis Correlation of Standard and Hybrid Standard/Thin-ply Composite Notched Test Specimens

Hybrid composite laminates, a combination of standard-and thin-ply layers, are structurally efficient and can offer opportunities for weight savings over composite laminates composed of only standard-thickness plies. Three hybrid laminates and three standard laminates were tested and analyzed to provide a baseline comparison for not only weight savings and strength benefits of the hybrid stacking sequences, but also to provide fidelity for progressive failure modeling techniques using the two different laminate designs. The standard and hybrid composite specimens were designed to have the same thickness and in-plane stiffness properties, as closely as possible, and were loaded with axial tension and axial compression. Small notched specimens were tested under compression and tension loading, and larger intermediate specimens were tested under tension loading, only. The finite element modeling represented the as-designed panel dimensions as well as the real-life testing conditions through various boundary conditions. The progressive failure modeling and analysis technique described, exhibited very good correlation with the hybrid test specimens, and therefore is appropriate for use in future investigations into more appealing stacking sequences. The progressive failure modeling can also be used to determine the proper usage of thin-ply materials in order to determine the most ideal strength benefits over their standard-ply counterparts.

Alana Zahn↗
NASA NTRS2008

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗
NASA NTRS1984

Finite-analytic numerical method for unsteady two-dimensional Navier-Stokes equations

A finite analytic (FA) numerical solution is developed for unsteady two-dimensional Navier-Stokes equations. The FA method utilizes the analytic solution in a small local element to formulate the algebraic representation of partial differential equations. The combination of linear and exponential functions that satisfy the governing equation is adopted as the boundary function, thereby improving the accuracy of the finite analytic solution. Two flows, one a starting cavity flow and the other a vortex shedding flow behind a rectangular block, are solved by the FA method. The starting square cavity flow is solved for Reynolds number of 400, 1000, and 2000 to show the accuracy and stability of the FA solution. The FA solution for flow over a rectangular block (H x H/4) predicts the Strouhal number for Reynolds numbers of 100 and 500 to be 0.156 and 0.125. Details of the flow patterns are given. In addition to streamlines and vorticity distribution, rest-streamlines are given to illustrate the vortex motion downstream of the block.

Chen, C.-J.↗
NASA NTRS1991

Lithospheric extension near Lake Mead, Nevada - A model for ductile flow in the lower crust

Small variations in gravity anomalies and topographic elevation observed in areas that have undergone highly variable amounts of upper crustal thinning can be satisfactorily explained by ductile flow of lower crustal material under the proper conditions. The boundary between the unextended Colorado Plateau and a strongly extended domain in the Basin and Range Province in the Lake Mead (Nevada) region is examined. Finite element modeling of Newtonian flow and power law creep shows that flow over the length scale of the eastern Basin and Range (500) km or more) corresponding to upper crustal extension by a factor of 1.4-3 over 10 million years requires effective viscosities less than 10 to the 18th - 10 to the 20th Pa s for ductile channels 10-25 km thick. Modeling suggests that these effective viscosities may be sustained by lower crustal material deforming at laboratory-derived power law creep rates. The longer-scale flow may require elevated crustal temperatures (more than 700 C), depending on the composition and material properties assumed. Under the boundary conditions assumed in this study the linear viscous flow models yield a satisfactory approximation to deformation by power law creep. This work suggests that flow in the lower crust may be a viable mechanism for producing small variations in total crustal thickness between strongly extended and less extended regions, and thereby explaining the relative uniformity in gravity and topography between such regions.

Kruse, Sarah↗
NASA NTRS1998

Model Analysis of an Aircraft Fueslage Panel using Experimental and Finite-Element Techniques

The application of Electro-Optic Holography (EOH) for measuring the center bay vibration modes of an aircraft fuselage panel under forced excitation is presented. The requirement of free-free panel boundary conditions made the acquisition of quantitative EOH data challenging since large scale rigid body motions corrupted measurements of the high frequency vibrations of interest. Image processing routines designed to minimize effects of large scale motions were applied to successfully resurrect quantitative EOH vibrational amplitude measurements

Fleming, Gary A.↗
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 NTRS1993

An analytical/numerical correlation study of the multiple concentric cylinder model for the thermoplastic response of metal matrix composites

The utility of a recently developed analytical micromechanics model for the response of metal matrix composites under thermal loading is illustrated by comparison with the results generated using the finite-element approach. The model is based on the concentric cylinder assemblage consisting of an arbitrary number of elastic or elastoplastic sublayers with isotropic or orthotropic, temperature-dependent properties. The elastoplastic boundary-value problem of an arbitrarily layered concentric cylinder is solved using the local/global stiffness matrix formulation (originally developed for elastic layered media) and Mendelson's iterative technique of successive elastic solutions. These features of the model facilitate efficient investigation of the effects of various microstructural details, such as functionally graded architectures of interfacial layers, on the evolution of residual stresses during cool down. The available closed-form expressions for the field variables can readily be incorporated into an optimization algorithm in order to efficiently identify optimal configurations of graded interfaces for given applications. Comparison of residual stress distributions after cool down generated using finite-element analysis and the present micromechanics model for four composite systems with substantially different temperature-dependent elastic, plastic, and thermal properties illustrates the efficacy of the developed analytical scheme.

Pindera, Marek-Jerzy↗
NASA NTRS1985

Thermal stress analysis of a new turbine shroud seal concept

The thermal stress field of a two piece turbine shroud seal concept was analyzed and results compared to one piece designs by finite element analysis. The two piece seal has independently formed structure (substrate) and ceramic components that are assembled at ambient conditions. The boundary conditions used for analysis were hot gas surface temperatures of 1370 and 1650 C (2500 and 3000 F) and cooled surface temperature of 700 C (1285 F). The resulting thermal stress field, of the two piece seal when compared to the one piece seals in the region of all ceramic material, was reduced substantially.

Handschuh, R. F.↗
NASA NTRS1986

Preconditioned conjugate residual methods for the solution of spectral equations

Conjugate residual methods for the solution of spectral equations are described. An inexact finite-difference operator is introduced as a preconditioner in the iterative procedures. Application of these techniques is limited to problems for which the symmetric part of the coefficient matrix is positive definite. Although the spectral equation is a very ill-conditioned and full matrix problem, the computational effort of the present iterative methods for solving such a system is comparable to that for the sparse matrix equations obtained from the application of either finite-difference or finite-element methods to the same problems. Numerical experiments are shown for a self-adjoint elliptic partial differential equation with Dirichlet boundary conditions, and comparison with other solution procedures for spectral equations is presented.

Wong, Y. S.↗
NASA NTRS1986

Mathematical aspects of finite element methods for incompressible viscous flows

Mathematical aspects of finite element methods are surveyed for incompressible viscous flows, concentrating on the steady primitive variable formulation. The discretization of a weak formulation of the Navier-Stokes equations are addressed, then the stability condition is considered, the satisfaction of which insures the stability of the approximation. Specific choices of finite element spaces for the velocity and pressure are then discussed. Finally, the connection between different weak formulations and a variety of boundary conditions is explored.

Gunzburger, M. D.↗