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 649 records · Page 36

A finite element computational method for high Reynolds number laminar flows

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables are interpolated using complete quadratic shape functions, and the pressure is interpolated using linear shape functions which are defined on a triangular element for the two-dimensional case and on a tetrahedral element for the three-dimensional case. The triangular element and the tetrahedral element are contained inside the complete bi- and tri-quadratic elements for velocity variables for two and three dimensional cases, respectively, so that the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow of Reynolds numbers 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorably with the finite difference computational results and/or experimental data available. It was found that the present method can capture the delicate pressure driven recirculation zones, that the method did not yield any spurious pressure modes, and that the method requires fewer grid points than the finite difference methods to obtain comparable computational results.

Kim, Sang-Wook↗

Structural analysis of ultra-high speed aircraft structural components

The buckling characteristics of a hypersonic beaded skin panel were investigated under pure compression with boundary conditions similar to those found in a wing mounted condition. The primary phases of analysis reported include: (1) experimental testing of the panel to failure; (2) finite element structural analysis of the beaded panel with the computer program NASTRAN; and (3) summary of the semiclassical buckling equations for the beaded panel under purely compressive loads. A comparison of each of the analysis methods is also included.

Lenzen, K. H.↗

Stress analysis of multilayered plates around circular holes

The elasticity problem for a thick plate consisting of two bonded dissimilar homogeneous layers is considered. It is assumed that the plate is infinite, contains a circular hole, and is subjected to axisymmetric external loads. A technique suitable to the formulation of three-dimensional boundary value problems is introduced and applied to the three-dimensional elasticity solution of an infinite plate with a hole. The effect of Poisson's ratio is studied in some detail, and the axisymmetric problem of a laminated plate is solved. The stresses at the interface are calculated and are observed to become unbounded at the hole boundary. An approximate treatment of the singularity is presented and the stress intensity factors are computed. Some of the results are compared with those found from the finite element method. The agreement between the two sets with regard to trend and magnitude is found to be very satisfactory.

Delale, F.↗

On the penetration of a hot diapir through a strongly temperature-dependent viscosity medium

The ascent of a hot spherical body through a fluid with a strongly temperature-dependent viscosity has been studied using an axisymmetric finite element method. Numerical solutions range over Peclet numbers of 0.1 - 1000 from constant viscosity up to viscosity variations of 100,000. Both rigid and stress-free boundary conditions were applied at the surface of the sphere. The dependence of drag on viscosity variation was shown to have no dependence on the stress boundary condition except for a Stokes flow scaling factor. A Nusselt number parameterization based on the stress-free constant viscosity functional dependence on the Peclet number scaled by a parameter depending on the viscosity structure fits both stress-free and rigid boundary condition data above viscosity variations of 100. The temperature scale height was determined as a function of sphere radius. For the simple physical model studied in this paper pre-heating is required to reduce the ambient viscosity of the country rock to less than 10 to the 22nd sq cm/s in order for a 10 km diapir to penetrate a distance of several radii.

Daly, S. F.↗

Component mode synthesis and large deflection vibration of complex structures. Volume 2: Single-mode large deflection vibrations of beams and plates using finite element method

A finite element method is presented for the large amplitude vibrations of complex structures that can be modelled with beam and rectangular plate elements subjected to harmonic excitation. Both inplane deformation and inertia are considered in the formulation. Derivation of the harmonic force and nonlinear stiffness matrices for a beam and a rectangular plate element are presented. Solution procedures and convergence characteristics of the finite element method are described. Nonlinear response to uniform and concentrated harmonic loadings and improved nonlinear free vibration results are presented for beams and rectangular plates of various boundary conditions.

Mei, Chuh↗

Multiple-mode nonlinear free and forced vibrations of beams using finite element method

Effects of large deflection geometric nonlinearity and multiple modes for free and forced vibrations of beams are investigated. Inplane displacement and inplane inertia are included in the formulation. The finite element method is employed. A harmonic force matrix is derived for forced vibration analysis. The relations of amplitude, frequency ratio and forcing intensity for beams of different boundary conditions and various load conditions are presented.

Mei, Chuh↗

A global-local finite element method suitable for parallel computations

A global-local finite element analysis procedure is developed based on the fast convergent nature of FEM in displacement. A special scheme is used to utilize the global displacement solution as boundary conditions for local regions of interest. In the local region, a refined mesh is used for further stress analysis. This global-local procedure can be easily programmed in parallel on MIMD multi-processor computers for significant time savings. A Sequent Balance 21000 system is used for demonstrating the parallel programming.

Sun, C. T.↗

Free vibrations of thin-walled semicircular graphite-epoxy composite frames

A detailed study is made of the effects of variations in lamination and material parameters of thin walled composite frames on their vibrational characteristics. The structures considered are semicircular thin walled frames with I and J sections. The flanges and webs of the frames are modeled by using 2-D shell and plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frame. The frequencies and modes predicted by the 2-D finite element model are compared with those obtained from experiments, as well as with the predictions of a non-dimensional thin walled beam finite element model. A detailed study is made of the sensitivity of the vibrational response to variations in the fiber orientation, material properties of the individual layers, and boundary conditions.

Carden, Huey D.↗

Free vibrations of thin-walled semicircular graphite-epoxy composite frames

A detailed study is made of the effects of variations in lamination and material parameters of thin walled composite frames on their vibrational characteristics. The structures considered are semicircular thin walled frames with I and J sections. The flanges and webs of the frames are modelled by using 2-D shell and plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frames. The frequencies and modes predicted by the 2-D finite element model are compared with those obtained from experiments, as well as with the predictions of a 1-D thin walled beam finite element model. A detailed study is made of the sensitivity of the vibrational response to variations in the fiber orientation, material properties of the individual layers, and boundary conditions.

Noor, Ahmed K.↗

Free vibrations of thin-walled semicircular graphite-epoxy composite frames

A detailed study is made of the effects of variations in lamination and material parameters of thin walled composite frames on their vibrational characteristics. The structures considered are semicircular thin walled frames with I and J sections. The flanges and webs of the frames are modeled by using 2-D shell and plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frame. The frequencies and modes predicted by the 2-D finite element model are compared with those obtained from experiments, as well as with the predictions of a non-dimensional thin walled beam finite element model. A detailed study is made of the sensitivity of the vibrational response to variations in the fiber orientation, material properties of the individual layers, and boundary conditions.

Carden, Huey D.↗

Free vibrations of thin-walled semicircular graphite-epoxy composite frames

A detailed study is made of the effects of variations in lamination and material parameters of thin-walled composite frames on their vibrational characteristics. The structures considered are semicircular thin-walled frames with I and J sections. The flanges and webs of the frames are modeled by using two-dimensional shell and plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frame. The frequencies and modes predicted by the two-dimensional finite-element model are compared with those obtained from experiments, as well as with the predictions of a one-dimensional, thin-walled-beam, finite-element model. A detailed study is made of the sensitivity of the vibrational response to variations in the fiber orientation, material properties of the individual layers, and boundary conditions.

Noor, Ahmed K.↗

Probabilistic finite elements for fatigue and fracture analysis

Attenuation is focused on the development of Probabilistic Finite Element Method (PFEM), which combines the finite element method with statistics and reliability methods, and its application to linear, nonlinear structural mechanics problems and fracture mechanics problems. The computational tool based on the Stochastic Boundary Element Method is also given for the reliability analysis of a curvilinear fatigue crack growth. The existing PFEM's have been applied to solve for two types of problems: (1) determination of the response uncertainty in terms of the means, variance and correlation coefficients; and (2) determination the probability of failure associated with prescribed limit states.

Belytschko, Ted↗

Drop Ejection From an Oscillating Rod

The dynamics of a drop of a Newtonian liquid that is pendant from or sessile on a solid rod that is forced to undergo time-periodic oscillations along its axis is studied theoretically. The free boundary problem governing the time evolution of the shape of the drop and the flow field inside it is solved by a method of lines using a finite element algorithm incorporating an adaptive mesh. When the forcing amplitude is small, the drop approaches a limit cycle at large times and undergoes steady oscillations thereafter. However, drop breakup is the consequence if the forcing amplitude exceeds a critical value. Over a wide range of amplitudes above this critical value, drop ejection from the rod occurs during the second oscillation period from the commencement of rod motion. Remarkably, the shape of the interface at breakup and the volume of the primary drop formed are insensitive to changes in forcing amplitude. The interface shape at times close to and at breakup is a multi-valued function of distance measured along the rod axis and hence cannot be described by recently popularized one-dimensional approximations. The computations show that drop ejection occurs without the formation of a long neck. Therefore, this method of drop formation holds promise of preventing formation of undesirable satellite droplets.

Wilkes, E. D.↗

Simulating the Structural Response of a Preloaded Bolted Joint

The present paper describes the structural analyses performed on a preloaded bolted-joint configuration. The joint modeled was comprised of two L-shaped structures connected together using a single bolt. Each L-shaped structure involved a vertical flat segment (or shell wall) welded to a horizontal segment (or flange). Parametric studies were performed using elasto-plastic, large-deformation nonlinear finite element analyses to determine the influence of several factors on the bolted-joint response. The factors considered included bolt preload, washer-surface-bearing size, edge boundary conditions, joint segment length, and loading history. Joint response is reported in terms of displacements, gap opening, and surface strains. Most of the factors studied were determined to have minimal effect on the bolted-joint response; however, the washer-bearing-surface size affected the response significantly.

Knight, Norman F., Jr.↗

Experimental investigation of the buckling characteristics of a beaded skin panel for a hypersonic aircraft - Including comparisons with finite element and classical analyses

Results of a compression test of a beaded panel intended for a proposed hypersonic aircraft are presented. The panel was tested to failure at room temperature to determine its buckling characteristics, in particular, to study the buckling caused by pure compression. The boundary conditions of the panel simulated as nearly as possible a wing mounted condition. Strain, out-of-plane deflection, and load data were measured, and elastic buckling strength as well as mode shapes of the panel were determined. Application of the moire technique is described.

Siegel, W. H.↗

Application of NASTRAN to large deflection supersonic flutter of panels

Flat panel flutter at high supersonic Mach number is analyzed using NASTRAN Level 16.0 by means of modifications to the code. Two-dimensional plate theory and quasi-steady aerodynamic theory are employed. The finite element formulation and solution procedure are presented. Modifications to the NASTRAN code are discussed. Convergence characteristics of the iteration processes are also briefly discussed. Effects of aerodynamic damping, boundary support condition and applied in-plane loading are included. Comparison of nonlinear vibration and linear flutter results with analytical solutions demonstrate that excellent accuracy is obtained with NASTRAN.

Mei, C.↗

Aerodynamic optimization, comparison, and trim design of canard and conventional high performance general aviation configurations

A design study has been conducted to optimize trim cruise flight of high performance general aviation canard aircraft which achieve minimum drag. In order to investigate the advantages and disadvantages of canard configured aircraft, corresponding conventional tail-aft 'baseline' aircraft were designed and used for comparison. Two-dimensional predictions were obtained by coupling inviscid results from a vortex panel multi-element program to a momentum integral boundary layer analysis. Using the results of the two-dimensional vortex panel analysis, a vortex lattice method was employed to predict the finite wing results. The analysis utilized a turbulent airfoil and a natural laminar airfoil which are two NASA state-of-the-art airfoil sections. The canard aircraft designs give quantitative results of wing and canard loadings, wing-to-canard moment arm ratios, and aspect ratio effects for trim cruise flight for a wide range of wing-to-canard area ratios. Both canard and baseline aircraft achieved a 25 to 30 percent improvement in performance over typical current technology aircraft, but high canard loading necessary for trim resulted in slightly poorer performance of the canard aircraft as compared to the baseline designs.

Keith, M. W.↗

Application of High Order Acoustic Finite Elements to Transmission Losses and Enclosure Problems

A family of acoustic finite elements was developed based on C continuity (acoustic pressure being the nodal variable) and the no-flow condition. The family include triangular, quadrilateral and hexahedral isoparametric elements with linear quadratic and cubic variation in modelling and distortion. Of greatest use in problems with irregular boundaries are the cubic isoparametric elements: the 32 node hexahedral element for three-dimensional systems; and the twelve node quadrilateral and ten node triangular elements for two-dimensional/axisymmetric applications. These elements were applied to problems involving cavity resonances, transmission loss in silencers and the study of end effects, using a Floating Point Systems 164 attached array processor accessed through an Amdahl 5860 mainframe. The elements are presently being used to study the end effects associated with duct terminations within finite enclosures. The transmission losses with various silencers and sidebranches in ducts is also being studied using the same elements.

Craggs, A.↗