Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape function”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Automating FEA programming

In this paper we briefly describe a combined symbolic and numeric approach for solving mathematical models on parallel computers. An experimental software system, PIER, is being developed in Common Lisp to synthesize computationally intensive and domain formulation dependent phases of finite element analysis (FEA) solution methods. Quantities for domain formulation like shape functions, element stiffness matrices, etc., are automatically derived using symbolic mathematical computations. The problem specific information and derived formulae are then used to generate (parallel) numerical code for FEA solution steps. A constructive approach to specify a numerical program design is taken. The code generator compiles application oriented input specifications into (parallel) FORTRAN77 routines with the help of built-in knowledge of the particular problem, numerical solution methods and the target computer.

Sharma, Naveen↗
NASA NTRS1993

Development of an hp-version finite element method for computational optimal control

The purpose of this research effort was to begin the study of the application of hp-version finite elements to the numerical solution of optimal control problems. Under NAG-939, the hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element. One possible drawback is the increased computational effort within each element required in implementing hp-version finite elements. We are trying to determine whether this computational effort is sufficiently offset by the reduction in the number of time elements used and improved Newton-Raphson convergence so as to be useful in solving optimal control problems in real time. Because certain of the element interior unknowns can be eliminated at the element level by solving a small set of nonlinear algebraic equations in which the nodal values are taken as given, the scheme may turn out to be especially powerful in a parallel computing environment. A different processor could be assigned to each element. The number of processors, strictly speaking, is not required to be any larger than the number of sub-regions which are free of discontinuities of any kind.

Hodges, Dewey H.↗
NASA NTRS1992

A piecewise continuous Timoshenko beam model for the dynamic analysis of tapered beam-like structures

Distributed parameter modeling offers a viable alternative to the finite element approach for modeling large flexible space structures. The introduction of the transfer matrix method into the continuum modeling process provides a very useful tool to facilitate the distributed parameter model applied to some more complex configurations. A uniform Timoshenko beam model for the estimation of the dynamic properties of beam-like structures has given comparable results. But many aeronautical and aerospace structures are comprised of non-uniform sections or sectional properties, such as aircraft wings and satellite antennas. This paper proposes a piecewise continuous Timoshenko beam model which is used for the dynamic analysis of tapered beam-like structures. A tapered beam is divided into several segments of uniform beam elements. Instead of arbitrarily assumed shape functions used in finite element analysis, the closed-form solution of the Timoshenko beam equation is used. Application of the transfer matrix method relates all the elements as a whole. By corresponding boundary conditions and compatible conditions a characteristic equation for the global tapered beam has been developed, from which natural frequencies can be derived. A computer simulation is shown in this paper, and compared with the results obtained from the finite element analysis. While piecewise continuous Timoshenko beam model decreases the number of elements significantly; comparable results to the finite element method are obtained.

Shen, Ji Yao↗
NASA NTRS1992

Wall shear stress estimates in coronary artery constrictions

Wall shear stress estimates from laminar boundary layer theory were found to agree fairly well with the magnitude of shear stress levels along coronary artery constrictions obtained from solutions of the Navier Stokes equations for both steady and pulsatile flow. The relatively simple method can be used for in vivo estimates of wall shear stress in constrictions by using a vessel shape function determined from a coronary angiogram, along with a knowledge of the flow rate.

Back, L. H.↗
NASA NTRS1992

Flap-lag damping in hover and forward flight with a three-dimensional wake

Prediction of lag damping is difficult owing to the delicate balance of drag, induced drag and Coriolis forces in the in-plane direction. Moreover, induced drag is sensitive to dynamic wake, both shed and trailing components, and thus its prediction requires adequate unsteady-wake representation. Accordingly, rigid-blade flap-lag equations are coupled with a three-dimensional finite-state wake model; three isolated rotor configurations with three, four and five blades are treated over a range of thrust levels, Lock numbers, lag frequencies and advance ratios. The investigation includes convergence characteristics of damping with respect to the number of shape functions and harmonics of the wake model for multiblade modes of low frequency (less than 1/rev.) to high frequency (greater than l/rev.). Predicted flap and lag damping levels are then compared with similar predictions with (1) rigid wake (no unsteady induced flow), (2) Loewy lift deficiency, and (3) dynamic inflow.

Manjunath, A. R.↗
NASA NTRS1993

Compton scattering and the gamma-ray power-law spectrum in Markarian 421

The nearest BL Lac object, Mrk 421, has a gamma-ray spectrum which is approximately flat in EF-sub E from E less than about 50 MeV to E greater than about 1 TeV. Inverse Compton scattering can explain this smooth spectrum, despite the structure in the Klein-Nishina cross section, if the injected electron distribution function is proportional to gamma exp -2, where gamma is the electron Lorentz factor. When this is the case, the structure imprinted on the steady state electron distribution function by the structure in the Klein-Nishina cross section is almost exactly compensated in the radiated spectrum. Because particle acceleration in strong shocks injects particles with this distribution function, this shape injection function is in fact quite plausible. Other blazars may be explained by the same model if the cutoff below TeV energies observed in other objects is due to pair production on the cosmological IR background, as suggested by Stecker et al. (1992).

Zdziarski, Andrzej A.↗
NASA NTRS1993

Development of an hp-version finite element method for computational optimal control

The purpose of this research effort is to develop a means to use, and to ultimately implement, hp-version finite elements in the numerical solution of optimal control problems. The hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element.

Hodges, Dewey H.↗
NASA NTRS1992

Improved assumed-stress hybrid shell element with drilling degrees of freedom for linear stress, buckling, and free vibration analyses

An improved four-node quadrilateral assumed-stress hybrid shell element with drilling degrees of freedom is presented. The formulation is based on Hellinger-Reissner variational principle and the shape functions are formulated directly for the four-node element. The element has 12 membrane degrees of freedom and 12 bending degrees of freedom. It has nine independent stress parameters to describe the membrane stress resultant field and 13 independent stress parameters to describe the moment and transverse shear stress resultant field. The formulation encompasses linear stress, linear buckling, and linear free vibration problems. The element is validated with standard tests cases and is shown to be robust. Numerical results are presented for linear stress, buckling, and free vibration analyses.

Rengarajan, Govind↗
NASA NTRS1994

Development of an adaptive hp-version finite element method for computational optimal control

In this research effort, the usefulness of hp-version finite elements and adaptive solution-refinement techniques in generating numerical solutions to optimal control problems has been investigated. Under NAG-939, a general FORTRAN code was developed which approximated solutions to optimal control problems with control constraints and state constraints. Within that methodology, to get high-order accuracy in solutions, the finite element mesh would have to be refined repeatedly through bisection of the entire mesh in a given phase. In the current research effort, the order of the shape functions in each element has been made a variable, giving more flexibility in error reduction and smoothing. Similarly, individual elements can each be subdivided into many pieces, depending on the local error indicator, while other parts of the mesh remain coarsely discretized. The problem remains to reduce and smooth the error while still keeping computational effort reasonable enough to calculate time histories in a short enough time for on-board applications.

Hodges, Dewey H.↗
NASA NTRS1995

Sensitivity analysis of the non-linear dynamic viscoplastic response of 2-d structures with respect to material parameters

A computational procedure is presented for evaluating the sensitivity coefficients of the viscoplastic response of structures subjected to dynamic loading. A state of plane stress is assumed to exist in the structure, a velocity strain-Cauchy stress formulation is used, and the geometric non-linearities arising from large strains are incorporated. The Jaumann rate is used as a frame indifferent stress rate. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. The equations of motion emanating from a finite element semi-discretization are integrated using an explicit central difference scheme with an implicit stress update. The sensitivity coefficients are evaluated using a direct differentiation approach. Since the domain of integration is the current configuration, the sensitivity coefficients of the spatial derivatives of the shape functions must be included. Numerical results are presented for a thin plate with a central cutout subjected to an in-plane compressive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients, and spatial distributions at selected times are presented.

Kulkarni, Makarand↗
NASA NTRS1994

A survey of the core-congruential formulation for geometrically nonlinear TL finite elements

This article presents a survey of the core-congruential formulation (CCF) for geometrically nonlinear mechanical finite elements based on the total Lagrangian (TL) kinematic description. Although the key ideas behind the CCF can be traced back to Rajasekaran and Murray in 1973, it has not subsequently received serious attention. The CCF is distinguished by a two-phase development of the finite element stiffness equations. The initial phase developed equations for individual particles. These equations are expressed in terms of displacement gradients as degrees of freedom. The second phase involves congruential-type transformations that eventually binds the element particles of an individual element in terms of its node-displacement degrees of freedom. Two versions of the CCF, labeled direct and generalized, are distinguished. The direct CCF (DCCF) is first described in general form and then applied to the derivation of geometrically nonlinear bar, and plane stress elements using the Green-Lagrange strain measure. The more complex generalized CCF (GCCF) is described and applied to the derivation of 2D and 3D Timoshenko beam elements. Several advantages of the CCF, notably the physically clean separation of material and geometric stiffnesses, and its independence with respect to the ultimate choice of shape functions and element degrees of freedom, are noted. Application examples involving very large motions solved with the 3D beam element display the range of applicability of this formulation, which transcends the kinematic limitations commonly attributed to the TL description.

Felippa, Carlos A.↗
NASA NTRS1994

Advance finite element modeling of rotor blade aeroelasticity

An advanced beam finite element has been developed for modeling rotor blade dynamics and aeroelasticity. This element is part of the Element Library of the Second Generation Comprehensive Helicopter Analysis System (2GCHAS). The element allows modeling of arbitrary rotor systems, including bearingless rotors. It accounts for moderately large elastic deflections, anisotropic properties, large frame motion for maneuver simulation, and allows for variable order shape functions. The effects of gravity, mechanically applied and aerodynamic loads are included. All kinematic quantities required to compute airloads are provided. In this paper, the fundamental assumptions and derivation of the element matrices are presented. Numerical results are shown to verify the formulation and illustrate several features of the element.

Straub, F. K.↗
NASA NTRS1993

Flap-lag damping in hover and forward flight with a three-dimensional wake

Prediction of lag damping is difficult owing to the delicate balance of drag, induced drag and Coriolis forces in the in-plane direction. Moreover, induced drag is sensitive to dynamic wake, both shed and trailing components, and thus its prediction requires adequate unsteady-wake representation. Accordingly, rigid-blade flap-lag equations are coupled with a three-dimensional finite-state wake model; three isolated rotor configurations with three, four and five blades are treated over a range of thrust levels, Lock numbers, lag frequencies and advance ratios. The investigation includes convergence characteristics of damping with respect to the number of radial shape functions and harmonics of the wake model for multiblade modes of low frequency (less than 1/rev.) to high frequency (greater than 1/rev.). Predicted flap and lag damping levels are then compared with similar predictions with (1) rigid wake (no unsteady induced flow, (2) Loewy lift defficiency and (3) dynamic inflow. The coverage also includes correlations with the measured lag regressive-mode damping in hover and forward flight and comparisons with similar correlations with dynamic wake model are consistently higher than the predictions with the dynamic inflow model; even for the low frequency lag regressive mode, the number of wake harmonics should at least be equal to twice the number of blades.

Manjunath, A. R.↗
NASA NTRS1994

The detailed balance requirement and general empirical formalisms for continuum absorption

Two general empirical formalisms are presented for the spectral density which take into account the deviations from the Lorentz line shape in the wing regions of resonance lines. These formalisms satisfy the detailed balance requirement. Empirical line shape functions, which are essential to provide the continuum absorption at different temperatures in various frequency regions for atmospheric transmission codes, can be obtained by fitting to experimental data.

Ma, Q.↗
NASA NTRS1994

An improved quasistatic line-shape theory: The effects of molecular motion on the line wings

A theory is presented for the modification of the line-shape functions and absorption coefficient due to the breakdown of the quasistatic approximation. This breakdown arises from the effects of molecular motion and increases the absorption in the near wings. Numerical calculations for the high-frequency wing of the nu(sub 3) band of CO2 broadened by Ar are reported and it is shown that these effects are significant near the bandhead. The importance of such corrections in other spectral regions and for other systems is discussed briefly.

Ma, Q.↗
NASA NTRS1995

Application of the boundary element method to the micromechanical analysis of composite materials

A new boundary element formulation for the micromechanical analysis of composite materials is presented in this study. A unique feature of the formulation is the use of circular shape functions to convert the two-dimensional integrations of the composite fibers to one-dimensional integrations. To demonstrate the applicability of the formulations, several example problems including elastic and thermal analysis of laminated composites and elastic analyses of woven composites are presented and the boundary element results compared to experimental observations and/or results obtained through alternate analytical procedures. While several issues remain to be addressed in order to make the methodology more robust, the formulations presented here show the potential in providing an alternative to traditional finite element methods, particularly for complex composite architectures.

Goldberg, R. K.↗
NASA NTRS1995

Vacuum Ultraviolet Absorption Measurements of Atomic Oxygen in a Shock Tube

The absorption of vacuum ultraviolet light by atomic oxygen has been measured in the Electric Arc-driven Shock Tube (EAST) Facility at NASA-Ames Research Center. This investigation demonstrates the instrumentation required to determine atomic oxygen concentrations from absorption measurements in impulse facilities. A shock wave dissociates molecular oxygen, producing a high temperature sample of atomic oxygen in the shock tube. A probe beam is generated with a Raman-shifted ArF excimer laser. By suitable tuning of the laser, absorption is measured over a range of wavelengths in the region of the atomic line at 130.49 nm. The line shape function is determined from measurements at atomic oxygen densities of 3x10(exp 17) and 9x10(exp 17) cm(exp -3). The broadening coefficient for resonance interactions is deduced from this data, and this value is in accord with available theoretical models.

Meyer, Scott Andrew↗
NASA NTRS1995

Vacuum Ultraviolet Absorption Measurements of Atomic Oxygen in a Shock Tube

The absorption of vacuum ultraviolet light by atomic oxygen has been measured in the Electric Arc-driven Shock Tube (EAST) Facility at NASA-Ames Research Center. This investigation demonstrates the instrumentation required to determine atomic oxygen concentrations from absorption measurements in impulse facilities. A shock wave dissociates molecular oxygen, producing a high temperature sample of atomic oxygen in the shock tube. A probe beam is generated with a Raman-shifted ArF excimer laser. By suitable tuning of the laser, absorption is measured over a range of wavelengths in the region of the atomic line at 130.49 nm. The line shape function is determined from measurements at atomic oxygen densities of 3 x 10(exp 17) and 9 x 10(exp 17) cm(exp -3). The broadening coefficient for resonance interactions is deduced from this data, and this value is in accord with available theoretical models.

Meyer, Scott Andrew↗