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 361 records · Page 20

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.↗

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↗

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.↗

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↗

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.↗

Numerical calculations of velocity and pressure distribution around oscillating airfoils

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

Bratanow, T.↗

Significance of foundation-soil separation in dynamic soil-structure interaction

THe dynamic response of flexible surface strip-foundations allowed to uplift is numerically obtained for externally applied forces of a transient time variation. The soil medium is represented by an isotropic, homogeneous and linear half-space. The soil is treated by a time domain boundary element method, while the flexible foundation is treated by the finite element method. It was concluded that intermediate relative stiffness leads to moderate deformations when uplift is permitted. Very flexible footings produce higher deformations in unilateral contact compared to bilateral contact, and thus should be considered in their design. Unilateral contact does not significantly increase deformations for stiff footings subjected to concentrated central loading. However, relatively large deformation differences occur when the loading is eccentric, necessitating consideration of uplift in their design.

Spyrakos, Constantine C.↗

Analytical approach to peel stresses in bonded composite stiffened panels

A closed-form solution was obtained for the stresses and displacements of two bonded beams. A system of two fourth-order and two second-order differential equations with the associated boundary equations was determined using a variational work approach. A FORTRAN computer program was devised to solve for the eigenvalues and eigenvectors of this system and to calculate the coefficients from the boundary conditions. The results were then compared with NASTRAN finite-element solutions and shown to agree closely.

Barkey, Derek A.↗

Calculation of skin-stiffener interface stresses in stiffened composite panels

A method for computing the skin-stiffener interface stresses in stiffened composite panels is developed. Both geometrically linear and nonlinear analyses are considered. Particular attention is given to the flange termination region where stresses are expected to exhibit unbounded characteristics. The method is based on a finite-element analysis and an elasticity solution. The finite-element analysis is standard, while the elasticity solution is based on an eigenvalue expansion of the stress functions. The eigenvalue expansion is assumed to be valid in the local flange termination region and is coupled with the finite-element analysis using collocation of stresses on the local region boundaries. Accuracy and convergence of the local elasticity solution are assessed using a geometrically linear analysis. Using this analysis procedure, the influence of geometric nonlinearities and stiffener parameters on the skin-stiffener interface stresses is evaluated.

Cohen, David↗

Effects of turbulence mixing, variable properties, and vaporization on spray droplet combustion

Combustion of liquid fuels in the form of spray droplets is simulated numerically. Various vaporization models are examined as to their performance in finite element calculations involving a turbulent flow field. The Eulerian coordinate for the gas and Lagrangian coordinate for the liquid spray droplets are coupled through source terms being updated in the equations of continuity, momentum, and energy. The k-epsilon and modified eddy breakup models are used for simulating turbulent spray combustion flow field. Numerical results for the droplet trajectories, droplet heating, recirculation characteristics, and effects of evaporation models are evaluated. It is also shown that the finite element method is advantageous in dealing with complex geometries, complex boundary conditions, adaptive unstructured grids.

Kim, Y. M.↗

A Conforming Multigrid Method for the Pure Traction Problem of Linear Elasticity: Mixed Formulation

A multigrid method using conforming P-1 finite element is developed for the two-dimensional pure traction boundary value problem of linear elasticity. The convergence is uniform even as the material becomes nearly incompressible. A heuristic argument for acceleration of the multigrid method is discussed as well. Numerical results with and without this acceleration as well as performance estimates on a parallel computer are included.

Lee, Chang-Ock↗

Influence of Ply Waviness on Fatigue Life of Tapered Composite Flexbeam Laminates

Nonlinear tapered flexbeam laminates, with significant ply waviness, were cut from a full-size composite rotor hub flexbeam. The specimens were tested under combined axial tension and cyclic bending loads. All of the specimens had wavy plies through the center and near the surfaces (termed marcelled areas), although for some of the specimens the surface marcels were very obvious, and for others they were much smaller. The specimens failed by first developing cracks through the marcels at the surfaces, and then delaminations grew from those cracks, in both directions. Delamination failure occurred in these specimens at significantly shorter fatigue lives than similar specimens without waviness, tested in ref. 2. A 2D finite element model was developed which closely approximated the flexbeam geometry, boundary conditions, and loading. In addition, the FE model duplicated the waviness observed in one of the test specimens. The model was analyzed using a geometrically nonlinear FE code. Modifications were made to the original model to reduce the amplitude of the marcels near the surfaces. The analysis was repeated for each modification. Comparisons of the interlaminar normal stresses, sigma(sub n), in the various models showed that under combined axial-tension and cyclic-bending loading, for marcels of the same aspect ratio, sigma(sub n) stresses increased as the distance along the taper, from thick to thin end, increased. For marcels of the same aspect ratio and at the same X-location along the taper, sigma(sub n) stresses decreased as the distance from the surface into the flexbeam interior increased. A technique was presented for determining the smallest acceptable marcel aspect ratio at various locations in the flexbeam.

Murri, Gretchen B.↗

Aeroheating Thermal Model Correlation for Mars Global Surveyor (MGS) Solar Array

The Mars Global Surveyor (MGS) Spacecraft made use of aerobraking to gradually reduce its orbit period from a highly elliptical insertion orbit to its final science orbit. Aerobraking produces a high heat load on the solar arrays, which have a large surface area exposed to the airflow and relatively low mass. To accurately model the complex behavior during aerobraking, the thermal analysis needed to be tightly coupled to the spatially varying, time dependent aerodynamic heating. Also, the thermal model itself needed to accurately capture the behavior of the solar array and its response to changing heat load conditions. The correlation of the thermal model to flight data allowed a validation of the modeling process, as well as information on what processes dominate the thermal behavior. Correlation in this case primarily involved detailing the thermal sensor nodes, using as-built mass to modify material property estimates, refining solar cell assembly properties, and adding detail to radiation and heat flux boundary conditions. This paper describes the methods used to develop finite element thermal models of the MGS solar array and the correlation of the thermal model to flight data from the spacecraft drag passes. Correlation was made to data from four flight thermal sensors over three of the early drag passes. Good correlation of the model was achieved, with a maximum difference between the predicted model maximum and the observed flight maximum temperature of less than 5%. Lessons learned in the correlation of this model assisted in validating a similar model and method used for the Mars Odyssey solar array aeroheating analysis, which were used during onorbit operations.

Amundsen, Ruth M.↗

Computational Fluid Dynamics-Based Aeroservoelastic Analysis with Hyper-X Applications

A finite element computational fluids dynamics-based aeroservoelastic analysis methodology is presented in this paper, in which both structural and fluids discretization are achieved by the finite element method, and their interaction is modeled by the transpiration boundary condition technique. In the fluids discipline either inviscid or viscous flow may be accounted for, usually employing unstructured grids.Adescription of a novel viscous flow solver employing unstructured grids is given in detail. Provisions are made for digital as well as analog controllers. These new aeroservoelastic analysis techniques are next applied for the solution of a number of example problems including the novel Hyper-X launch vehicle. Experimental and actual flight test data are also compared with analysis results that signify to the efficacy and accuracy of the newly developed solution procedures.

Gupta, K. K.↗

Use of Simple Continuum Solutions in Finite Element Alternating Method for Fracture Problems

The performance of the finite element alternating (FEAM) method for two-dimensional crack problems is studied with respect to a polynomial pressure distribution fitted to the crack face stresses. The FEAM alternates between the analytical solution of crack in an infinite plate subjected to arbitrary polynomial distribution and a finite element solution of an uncracked body to satisfy the required boundary conditions in the crack problem. In this paper, the FEAM is applied to embedded crack and edge crack problems. For embedded crack problems, all of the constant, linear, and quadratic ( N=0,1, or 2, respectively) pressure distributions yield very accurate results with this algorithm with 4 to 5 iterations. The edge crack problems, on the other hand, require much higher order polynomials distributions (N=5 to 6) to yield accurate solutions. For slant edge crack problems, the mode-I stress-intensity factors have better accuracy than the mode-II stress-intensity factors for the same convergence tolerance.

Krishnamurthy, T.↗

TPSAS-NF1676L-17800-DND

Currently, there are two national challenge problems that guide much of the research in durability and damage tolerance at NASA Langley. The first, Airframe Digital Twin, is a concept that combines as-built vehicle components, as-experienced loads and environments, and other vehicle-specific characteristics to enable ultrahigh fidelity modeling of aircraft and spacecraft throughout their service lives. The second, Materials Genome Initiative, is an analog to the Human Genome Project, and is intended to improve the rate at which materials scientists can discover, understand fundamental physics, and improve material systems. This presentation will highlight several research projects ongoing at NASA Langley that are in support of the above challenge problems. Two of those topics will be the subject of detailed discussion. First, investigations of microstructurally-small fatigue cracking (MSFC) in Al-2Cu and Al-4Cu, fabricated in-house, will be presented. Single- and oligo-crystals of Al-Cu specimens were loaded in uniaxial fatigue, while high-resolution in-situ measurements of deformation were made using image correlation (IC) in a scanning-electron microscope (SEM). The Al-Cu specimens were then replicated as crystal plasticity finite element models (CPFEM), where evolution of slip localization near grain boundaries was computed. Comparison among experiment and CPFEM is made. In addition, XRay diffraction measurements of the as-fabricated specimens were carried out, where direct measurements of the embedded copper precipitates were made, and their influence on growing MSFCs were directly observed. The second main topic will illustrate ongoing work in the area of so-called damage-sensing particles. In this work, shape-memory alloys are embedded in an aluminum alloy matrix. Upon the propagation of a fatigue crack, these particles undergo a strain-induced phase transformation which is detected using an acoustic sensor, providing real-time information on propagating cracks. Experiments and simulations regarding the development of this system will also be detailed.

Jacob Hochhalter↗