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

The NESSUS finite element code

The objective of this development is to provide a new analysis tool which integrates the structural modeling versatility of a modern finite element code with the latest advances in the area of probabilistic modeling and structural reliability. Version 2.0 of the NESSUS finite element code was released last February, and is currently being exercised on a set of problems which are representative of typical Space Shuttle Main Engine (SSME) applications. NESSUS 2.0 allows linear elastostatic and eigenvalue analysis of structures with uncertain geometry, material properties and boundary conditions, which are subjected to a random mechanical and thermal loading environment. The NESSUS finite element code is a key component in a broader software system consisting of five major modules. NESSUS/EXPERT is an expert system under development at Southwest Research Institute, with the objective of centralizing all component-specific knowledge useful for conducting probabilistic analysis of typical Space Shuttle Main Engine (SSME) components. NESSUS/FEM contains the finite element code used for the structural analysis and parameter sensitivity evaluation of these components. The task of parametrizing a finite element mesh in terms of the random variables present is facilitated with the use of the probabilistic data preprocessor in NESSUS/PRE. An external database file is used for managing the bulk of the data generated by NESSUS/FEM.

Dias, J. B.↗
NASA NTRS1989

Finite element large-amplitude free and forced vibrations of rectangular thin composite plates

A finite element formulation is presented for determining the large-amplitude free and steady-state forced vibration response of arbitrarily laminated anisotropic composite rectangular thin plates. The nonlinear stiffness and harmonic force matrices of an arbitrarily laminated composite rectangular plate element are developed for nonlinear free and forced vibration analyses. The linearized updated-mode method with nonlinear time function approximation is employed for the solution of the system nonlinear eigenvalue equations. The amplitude-frequency relations for convergence with gridwork refinement, different boundary conditions, aspect ratios, lamination angles and number of plies are presented. The finite element results are compared with available approximate continuum solutions.

Chiang, C. K.↗
NASA NTRS1985

On the existence and stability conditions for mixed-hybrid finite element solutions based on Reissner's variational principle

The extensions of Reissner's two-field (stress and displacement) principle to the cases wherein the displacement field is discontinuous and/or the stress field results in unreciprocated tractions, at a finite number of surfaces ('interelement boundaries') in a domain (as, for instance, when the domain is discretized into finite elements), is considered. The conditions for the existence, uniqueness, and stability of mixed-hybrid finite element solutions based on such discontinuous fields, are summarized. The reduction of these global conditions to local ('element') level, and the attendant conditions on the ranks of element matrices, are discussed. Two examples of stable, invariant, least-order elements - a four-node square planar element and an eight-node cubic element - are discussed in detail.

Karlovitz, L. A.↗
NASA NTRS1987

Integration of geometric modeling and advanced finite element preprocessing

The structure to a geometry based finite element preprocessing system is presented. The key features of the system are the use of geometric operators to support all geometric calculations required for analysis model generation, and the use of a hierarchic boundary based data structure for the major data sets within the system. The approach presented can support the finite element modeling procedures used today as well as the fully automated procedures under development.

Shephard, Mark S.↗
NASA NTRS2015

Investigation of Micro-Scale Architectural Effects on Damage of Composites

This paper presents a three-dimensional, energy based, anisotropic, stiffness reduction, progressive damage model for composite materials and composite material constituents. The model has been implemented as a user-defined constitutive model within the Abaqus finite element software package and applied to simulate the nonlinear behavior of a damaging epoxy matrix within a unidirectional composite material. Three different composite microstructures were considered as finite element repeating unit cells, with appropriate periodicity conditions applied at the boundaries. Results representing predicted transverse tensile, longitudinal shear, and transverse shear stress-strain curves are presented, along with plots of the local fields indicating the damage progression within the microstructure. It is demonstrated that the damage model functions appropriately at the matrix scale, enabling localization of the damage to simulate failure of the composite material. The influence of the repeating unit cell geometry and the effect of the directionality of the applied loading are investigated and discussed.

damage↗
NASA NTRS2004

Investigation of Three-Dimensional Stress Fields and Slip Systems for FCC Single Crystal Superalloy Notched Specimens

Metals and their alloys, except for a few intermetallics, are inherently ductile, i.e. plastic deformation precedes fracture in these materials. Therefore, resistance to fracture is directly related to the development of the plastic zone at the crack tip. Recent studies indicate that the fracture toughness of single crystals depends on the crystallographic orientation of the notch as well as the loading direction. In general, the dependence of crack propagation resistance on crystallographic orientation arises from the anisotropy of (i) elastic constants, (ii) plastic deformation (or slip), and (iii) the weakest fracture planes (e.g. cleavage planes). Because of the triaxial stress state at the notch tips, many slip systems that otherwise would not be activated during uniaxial testing, become operational. The plastic zone formation in single crystals has been tackled theoretically by Rice and his co-workers and only limited experimental work has been conducted in this area. The study of the stresses and strains in the vicinity of a FCC single crystal notch tip is of relatively recent origin. We present experimental and numerical investigation of 3D stress fields and evolution of slip sector boundaries near notches in FCC single crystal tension test specimens, and demonstrate that a 3D linear elastic finite element model that includes the effect of material anisotropy is shown to predict active slip planes and sectors accurately. The slip sector boundaries are shown to have complex curved shapes with several slip systems active simultaneously near the notch. Results are presented for surface and mid-plane of the specimens. The results demonstrate that accounting for 3D elastic anisotropy is very important for accurate prediction of slip activation near FCC single crystal notches loaded in tension. Results from the study will help establish guidelines for fatigue damage near single crystal notches.

Arakere, Nagaraj K.↗
NASA NTRS1999

Development of Analysis Methods for Designing with Composites

The project involved the development of new analysis methods to achieve efficient design of composite structures. We developed a complex variational formulation to analyze the in-plane and bending coupling response of an unsymmetrically laminated plate with an elliptical cutout subjected to arbitrary edge loading as shown in Figure 1. This formulation utilizes four independent complex potentials that satisfy the coupled in-plane and bending equilibrium equations, thus eliminating the area integrals from the strain energy expression. The solution to a finite geometry laminate under arbitrary loading is obtained by minimizing the total potential energy function and solving for the unknown coefficients of the complex potentials. The validity of this approach is demonstrated by comparison with finite element analysis predictions for a laminate with an inclined elliptical cutout under bi-axial loading.The geometry and loading of this laminate with a lay-up of [-45/45] are shown in Figure 2. The deformed configuration shown in Figure 3 reflects the presence of bending-stretching coupling. The validity of the present method is established by comparing the out-of-plane deflections along the boundary of the elliptical cutout from the present approach with those of the finite element method. The comparison shown in Figure 4 indicates remarkable agreement. The details of this method are described in a manuscript by Madenci et al. (1998).

Madenci, E.↗
NASA NTRS1982

Three-dimensional analysis of 0/90s and 90/0s laminates with a central circular hole

Stress distributions were calculated near a circular hole in laminates, using a three dimensional finite element analysis. These stress distributions were presented three ways: through the thickness at the hole boundary, along radial lines at the 0/90 and 90/0 interfaces, and around the hole at these interfaces. The interlaminar normal stress, and the shear stress, distributions had very steep gradients near the hole boundary, suggesting interlaminar stress singularities. The largest compressive stress occurred at about 60 deg from the load axis. A simple procedure was introduced to calculate interlaminar stresses near the hole boundary. It used stresses calculated by an exact two dimensional analysis of a laminate with a hole as input to a quasi three dimensional model. It produced stresses that agreed closely with those from the three dimensional finite element model.

Raju, I. S.↗
NASA NTRS1989

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Chung, T. J.↗
NASA NTRS2001

A Hybrid Numerical Analysis Method for Structural Health Monitoring

A new hybrid surface-integral-finite-element numerical scheme has been developed to model a three-dimensional crack propagating through a thin, multi-layered coating. The finite element method was used to model the physical state of the coating (far field), and the surface integral method was used to model the fatigue crack growth. The two formulations are coupled through the need to satisfy boundary conditions on the crack surface and the external boundary. The coupling is sufficiently weak that the surface integral mesh of the crack surface and the finite element mesh of the uncracked volume can be set up independently. Thus when modeling crack growth, the finite element mesh can remain fixed for the duration of the simulation as the crack mesh is advanced. This method was implemented to evaluate the feasibility of fabricating a structural health monitoring system for real-time detection of surface cracks propagating in engine components. In this work, the authors formulate the hybrid surface-integral-finite-element method and discuss the mechanical issues of implementing a structural health monitoring system in an aircraft engine environment.

Forth, Scott C.↗
NASA NTRSDate not supplied

Blade Mistuning Analysis for Boundary Layer Inlet Flow Distortion-Tolerant Fan Blades

NASA Glenn Research Center has tested the BLI2DTF propulsor including the fan rotor. BLI2DTF stands for Boundary Layer Ingesting Inlet/Distortion-Tolerant Fan. A comparison of the vibratory response between analysis results and experimental wind tunnel test results yielded reasonable agreement, yet some differences were noted. Blade mistuning effects have been identified as a potential cause for these differences. A tuned forced response analysis in a turbomachinery engine rotor system assumes that all of the blades are identical, so that the natural frequencies of individual blades are the same. In contrast to this assumption, the wind tunnel tests of the BLI2DTF rotor system revealed that the individual blade response frequencies varied. These variations in individual blade frequencies may be attributed to blade mistuning due to potential differences in the structural and material properties of individual blades. For an understanding of these variations, it is necessary to analyze the influence of blade mistuning on the BLI2DTF rotor system, which will provide a more precise assessment of its vibratory stresses and high cycle fatigue life. The purpose of the study was to investigate the influence of blade mistuning on the dynamic characteristics of this rotor system. At the same fan operation condition, a comparative forced response analysis was performed for both tuned and mistuned BLI2DTF rotor systems. A comparison of the results from these analyses found that the mistuning effects on the present BLI2DTF rotor system are notable. As a result, fan blade mistuning is to be considered a contributor to the levels of vibratory dynamic stresses observed during the wind tunnel experiment.

Blade mistuning↗
DOE OSTI2020

Arctic Coastal Erosion: Modeling and Experimentation

Increasing Arctic coastal erosion rates have put critical infrastructure and native communities at risk while also mobilizing ancient organic carbon into modern carbon cycles. Although the Arctic comprises one-third of the global coastline and has some of the fastest eroding coasts, current tools for quantifying permafrost erosion are unable to explain the episodic, storm-driven erosion events. Our approach, mechanistically coupling oceanographic predictions with a terrestrial model to capture the thermo-mechanical dynamics of erosion, enables this much needed treatment of transient erosion events. The Arctic Coastal Erosion Model consists of oceanographic and atmospheric boundary conditions that force a coastal terrestrial permafrost environment in Albany (a multi-physics based finite element model). An oceanographic modeling suite (consisting of WAVEWATCH III, Delft3D-FLOW, and Delft3D-WAVE) produced time-dependent surge and run-up boundary conditions for the terrestrial model. In the terrestrial model, a coupling framework unites the mechanical and thermal aspects of erosion. 3D stress/strain fields develop in response to a plasticity model of the permafrost that is controlled by the frozen water content determined by modeling 3D heat conduction and solid-liquid phase change. This modeling approach enables failure from any allowable deformation (block failure, slumping, etc.). Extensive experimental work has underpinned the ACE Model development including field campaigns to measure in situ ocean and erosion processes, strength properties derived from thermally driven geomechanical experiments, as well as extensive physical composition and geochemical analyses. Combined, this work offers the most comprehensive and physically grounded treatment of Arctic coastal erosion available in the literature. The ACE model and experimental results can be used to inform scientific understanding of coastal erosion processes, contribute to estimates of geochemical and sediment land-to-ocean fluxes, and facilitate infrastructure susceptibility assessments.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

COMGEN: A computer program for generating finite element models of composite materials at the micro level

COMGEN (Composite Model Generator) is an interactive FORTRAN program which can be used to create a wide variety of finite element models of continuous fiber composite materials at the micro level. It quickly generates batch or session files to be submitted to the finite element pre- and postprocessor PATRAN based on a few simple user inputs such as fiber diameter and percent fiber volume fraction of the composite to be analyzed. In addition, various mesh densities, boundary conditions, and loads can be assigned easily to the models within COMGEN. PATRAN uses a session file to generate finite element models and their associated loads which can then be translated to virtually any finite element analysis code such as NASTRAN or MARC.

Melis, Matthew E.↗
DOE OSTI2020

A continuum framework for coupled solid deformation–fluid flow through anisotropic elastoplastic porous media

We present a continuum framework for coupled solid deformation–fluid flow in anisotropic elastoplastic porous media. A thermodynamic formulation of the coupled processes gives rise to an anisotropic Biot tensor that is a function of drained elastic tangent moduli tensor of the solid skeleton and the intrinsic bulk modulus of the solid constituent. Two effective stress measures emerge from the formulation, namely, σ', which is energy-conjugate to the elastic strain, and σ'', which is energy-conjugate to the plastic strain. For the special case of transverse isotropy that is commonly encountered in natural rocks, the Biot tensor can be expressed in terms of its normal and tangential components to the bedding plane, along with a microstructure tensor. Apart from its thermodynamic consistency, an advantage of this new formulation is that standard mixed finite element formulation can be employed to discretize the domain and solve initial boundary-value problems. We conduct plane strain simulations of coupled solid deformation–fluid flow in a transversely isotropic porous medium to demonstrate the impacts of material anisotropy, stress history, and the Biot tensor on the system response.

42 ENGINEERING↗
NASA NTRS1950

Linear theory of boundary effects in open wind tunnels with finite jet lengths

In the first part, the boundary conditions for an open wind tunnel (incompressible flow) are examined with special reference to the effects of the closed entrance and exit sections. Basic conditions are that the velocity must be continuous at the entrance lip and that the velocities in the upstream and downstream closed portions must be equal. In the second part, solutions are derived for four types of two-dimensional open tunnels, including one in which the pressures on the two free surfaces are not equal. Numerical results are given for every case. In general, if the lifting element is more than half the tunnel height from the inlet, the boundary effect at the lifting element is the same as for an infinitely long open tunnel. In the third part, a general method is given for calculating the boundary effect in an open circular wind tunnel of finite jet length. Numerical results are given for a lifting element concentrate at a point on the axis.

Katzoff, S↗
NASA NTRS1973

Finite element analysis of unsteady incompressible flow around an oscillating obstacle of arbitrary shape.

An analytical procedure based on Navier-Stokes equations was developed for representing unsteady flow patterns around oscillating obstacles. A variational formulation of the Helmholtz vorticity equation was discretized in finite element form and integrated numerically. At each step of the numerical integration the velocity field around the obstacle was determined from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints at each time step of the numerical integration. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗