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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1990

A weak Hamiltonian finite element method for optimal control problems

A temporal finite element method based on a mixed form of the Hamiltonian weak principle is developed for dynamics and optimal control problems. The mixed form of Hamilton's weak principle contains both displacements and momenta as primary variables that are expanded in terms of nodal values and simple polynomial shape functions. Unlike other forms of Hamilton's principle, however, time derivatives of the momenta and displacements do not appear therein; instead, only the virtual momenta and virtual displacements are differentiated with respect to time. Based on the duality that is observed to exist between the mixed form of Hamilton's weak principle and variational principles governing classical optimal control problems, a temporal finite element formulation of the latter can be developed in a rather straightforward manner. Several well-known problems in dynamics and optimal control are illustrated. The example dynamics problem involves a time-marching problem. As optimal control examples, elementary trajectory optimization problems are treated.

Hodges, Dewey H.↗

NASA NTRS1990

Finite element solution of optimal control problems with inequality constraints

A finite-element method based on a weak Hamiltonian form of the necessary conditions is summarized for optimal control problems. Very crude shape functions (so simple that element numerical quadrature is not necessary) can be used to develop an efficient procedure for obtaining candidate solutions (i.e., those which satisfy all the necessary conditions) even for highly nonlinear problems. An extension of the formulation allowing for discontinuities in the states and derivatives of the states is given. A theory that includes control inequality constraints is fully developed. An advanced launch vehicle (ALV) model is presented. The model involves staging and control constraints, thus demonstrating the full power of the weak formulation to date. Numerical results are presented along with total elapsed computer time required to obtain the results. The speed and accuracy in obtaining the results make this method a strong candidate for a real-time guidance algorithm.

Bless, Robert R.↗

NASA NTRS1991

Advanced development of the boundary element method for elastic and inelastic thermal stress analysis

The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.

Henry, Donald P., Jr.↗

NASA NTRS1991

Opposed Flow Flame Spread in Normal, Enhanced and Reduced Gravity

Experimental and theoretical aspects of opposed-flow flame spread over solid fuels are presented with emphasis on the microgravity environments of spacecraft. For high opposing flow velocities, spread rate decreases with increasing velocity eventually leading to flame blowoff due to kinetic effects. At low opposing flow velocities, where diffusional effects are slowed and radiation becomes important, flame spread rate increases with increasing flow velocity. Extinction at low velocities is due to radiative effects. Modeling efforts that include radiation, both solid surface radiation and gas-phase radiation, predict qualitatively the experimental trends observed. Computationally, gas-phase radiation is conveniently included in solution of the conservation equations by employing a Plank mean absorption coefficient, a fraction of radiation that is fed back to the surface, and a shape function that describes the radiative flux distribution along the surface.

Altenkirch, Robert A.↗

NASA NTRS1991

An integrated airloads-inflow model for use in rotor aeroelasticity and control analysis

The airloads-inflow model presented is ideally suited for rotor aeroelasticity and control analysis, especially where eigenvalues are sought. The model, which is derived from thin-airfoil theory with an extension that allows for other airfoil effects, is of hierarchical-lift type and may be simplified to any degree desired; the inflow is also expanded in shape functions which may be truncated to any desired texture. Both models are written in closed-form, state-variable equations without hidden states and without integrals over the wake.

Peters, David A.↗

NASA NTRS1992

ANS shell elements with improved transverse shear accuracy

A method of forming assumed natural coordinate strain (ANS) plate and shell elements is presented. The ANS method uses equilibrium based constraints and kinematic constraints to eliminate hierarchical degrees of freedom which results in lower order elements with improved stress recovery and displacement convergence. These techniques make it possible to easily implement the element into the standard finite element software structure, and a modified shape function matrix can be used to create consistent nodal loads.

Jensen, Daniel D.↗

NASA NTRS1992

A curvilinear, anisotropic, p-version, brick finite element based on geometric entities

A 'brick' solid finite element is presently developed on the basis of the p-version analysis, and used to demonstrate the FEM concept of 'geometric entities'. This method eliminates interelement discontinuities between low- and high-order elements, allowing very fine control over the shape-function order in various parts of the model. Attention is given to the illustrative cases of a one-element model of an elliptic pipe, and a square cross-section cantilevered beam.

Hinnant, Howard E.↗

NASA NTRS1992

A far wing line shape theory and its application to the foreign-broadened water continuum absorption. III

The far wing line shape theory developed previously and applied to the calculation of the continuum absorption of pure water vapor is extended to foreign-broadened continua. Explicit results are presented for H2O-N2 and H2O-CO2 in the frequency range from 0 to 10,000/cm. For H2O-N2 the positive and negative resonant frequency average line shape functions and absorption coefficients are computed for a number of temperatures between 296 and 430 K for comparison with available laboratory data. In general the agreement is very good.

Ma, Q.↗

NASA NTRS1992

A finite element-boundary integral method for cavities in a circular cylinder

Conformal antenna arrays offer many cost and weight advantages over conventional antenna systems. However, due to a lack of rigorous mathematical models for conformal antenna arrays, antenna designers resort to measurement and planar antenna concepts for designing non-planar conformal antennas. Recently, we have found the finite element-boundary integral method to be very successful in modeling large planar arrays of arbitrary composition in a metallic plane. We extend this formulation to conformal arrays on large metallic cylinders. In this report, we develop the mathematical formulation. In particular, we discuss the shape functions, the resulting finite elements and the boundary integral equations, and the solution of the conformal finite element-boundary integral system. Some validation results are presented and we further show how this formulation can be applied with minimal computational and memory resources.

Kempel, Leo C.↗

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