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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS1974

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

NASA NTRS1987

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

NASA NTRS1987

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

NASA NTRS1987

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↗

NASA NTRS1990

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

NASA NTRS1996

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↗

NASA NTRS1999

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

NASA NTRS2003

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

NASA NTRS2007

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

NASA NTRS2018

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

NASA NTRSDate not supplied

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↗

NASA NTRS1986

Computation of turbulent boundary layer flows with an algebraic stress turbulence model

An algebraic stress turbulence model is presented, characterized by the following: (1) the eddy viscosity expression is derived from the Reynolds stress turbulence model; (2) the turbulent kinetic energy dissipation rate equation is improved by including a production range time scale; and (3) the diffusion coefficients for turbulence equations are adjusted so that the kinetic energy profile extends further into the free stream region found in most experimental data. The turbulent flow equations were solved using a finite element method. Examples include: fully developed channel flow, fully developed pipe flow, flat plate boundary layer flow, plane jet exhausting into a moving stream, circular jet exhausting into a moving stream, and wall jet flow. Computational results compare favorably with experimental data for most of the examples considered. Significantly improved results were obtained for the plane jet flow, the circular jet flow, and the wall jet flow; whereas the remainder are comparable to those obtained by finite difference methods using the standard kappa-epsilon turbulence model. The latter seems to be promising with further improvement of the expression for the eddy viscosity coefficient.

Kim, Sang-Wook↗

NASA NTRS1980

Stress analysis and buckling of J-stiffened graphite-epoxy panel

A graphite epoxy shear panel with bonded on J stiffeners was investigated. The panel, loaded to buckling in a picture frame shear test is described. Two finite element models, each of which included the doubler material bonded to the panel skin under the stiffeners and at the panel edges, were used to make a stress analysis of the panel. The shear load distributions in the panel from two commonly used boundary conditions, applied shear load and applied displacement, were compared with the results from one of the finite element models that included the picture frame test fixture.

Davis, R. C.↗

NASA NTRS1994

Dynamic analysis of Space Shuttle/RMS configuration using continuum approach

The initial assembly of Space Station Freedom involves the Space Shuttle, its Remote Manipulation System (RMS) and the evolving Space Station Freedom. The dynamics of this coupled system involves both the structural and the control system dynamics of each of these components. The modeling and analysis of such an assembly is made even more formidable by kinematic and joint nonlinearities. The current practice of modeling such flexible structures is to use finite element modeling in which the mass and interior dynamics is ignored between thousands of nodes, for each major component. The model characteristics of only tens of modes are kept out of thousands which are calculated. The components are then connected by approximating the boundary conditions and inserting the control system dynamics. In this paper continuum models are used instead of finite element models because of the improved accuracy, reduced number of model parameters, the avoidance of model order reduction, and the ability to represent the structural and control system dynamics in the same system of equations. Dynamic analysis of linear versions of the model is performed and compared with finite element model results. Additionally, the transfer matrix to continuum modeling is presented.

Ramakrishnan, Jayant↗

NASA NTRS1980

Application of a Galerkin finite element method to atmospheric transport problems

Numerical simulation of the movement of a contaminant within the atmosphere presents difficulties due to the multidimensionality of the problem and the fact that the horizontal transport is usually convection dominated, that the boundary conditions are mixed, and that both slow and fast atmospheric chemical reactions can be important. In this study, numerical experiments using a Crank-Nicolson Galerkin finite element method to solve the time-dependent partial differential equations demonstrate the applicability and accuracy of this method for the variety of conditions encountered in atmospheric pollutant modeling. The Crank-Nicolson Galerkin method using piecewise linear, piecewise cubic Hermite polynomials, and upwind finite elements is shown to accurately model the pure convection of initial wave forms. Numerical results studying the interactions of convection, diffusion, chemical reaction, pollutant removal, and the effects of contaminant emission source strength, source location and multiple sources are also presented.

Carmichael, G. R.↗

NASA NTRS2006

A Numerical and Experimental Study of Compression-Loaded Composite Panels With Cutouts

Results from a numerical and experimental study on the effects of laminate orthotropy and circular cutout size on the response of compression-loaded composite curved panels are presented. Several 60-in-radius composite panels with four different laminate configurations were tested with cutout diameters that range from 10% to 60% of the panel width. Finite-element analyses were performed for each panel in order to identify the effects boundary conditions, measured initial geometric imperfections and thickness variations had on the nonlinear and buckling behavior of the panels. The compression-loaded panels considered herein exhibited two separate types of behavior depending on the laminate stacking sequence and cutout size. More specifically, some of the panels exhibited the classical snap-through type buckling response; however, some of the panels exhibited a monotonically increasing stable response and achieved compressive loads in excess of twice the predicted linear bifurcation buckling load. In general, the finite-element analyses were able to predict accurately the nonlinear response and buckling loads of the panels and the prebuckling and postbuckling out-of-plane deformations and strains.

Thornburgh, Robert P.↗

NASA NTRS1989

Thermo-viscoelastic analysis of composite materials

The thermo-viscoelastic boundary value problem for anisotropic materials is formulated and a numerical procedure is developed for the efficient analysis of stress and deformation histories in composites. The procedure is based on the finite element method and therefore it is applicable to composite laminates containing geometric discontinuities and complicated boundary conditions. Using the present formulation, the time-dependent stress and strain distributions in both notched and unnotched graphite/epoxy composites have been obtained. The effect of temperature and ply orientation on the creep and relaxation response is also studied.

Lin, Kuen Y.↗

NASA NTRS2022

A Benchmark Example for Delamination Propagation Predictions Based on the Calibrated End-Loaded Split Specimen

A benchmark example based on the Calibrated End-Loaded Split (C-ELS) specimen is developed and used to assess the performance of recently developed delamination propagation capabilities in the Abaqus/Standard finite element code. The C-ELS specimen has the advantage of a longer region of stable delamination propagation compared to the existing mode II benchmark case. The new benchmark example may therefore provide a better assessment tool by enabling more stable crack growth in regions further away from the boundary conditions or load application. First, a benchmark result is created manually using two-dimensional finite element models of the C-ELS specimen with different delamination lengths. Second, the performance of the delamination propagation capabilities in the Abaqus/Standard finite element code are assessed by comparing the results to the benchmark case. Two examples with different starter delamination lengths are studied. A shorter starter length is chosen to create a scenario with unstable delamination propagation. A longer delamination causes stable delamination propagation. Detailed results from three-dimensional analyses with aligned and mis-aligned meshes and two levels of mesh refinement are provided for several permutations of numerical input parameters. In general, good agreement can be achieved between the results obtained from the quasi-static propagation analysis and the benchmark analysis. However, particular non-default settings are found to be most reliable, accurate, and numerically efficient. Numerical artifacts including anomalous unreleased nodes in the crack wake and zig-zag crack fronts occur, and further development of the Abaqus/Standard VCCT propagation may be required. Use of the benchmark case to assess the continuous improvements in one finite element code illustrates the value of establishing benchmark solutions.

delamination↗