Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite expression method”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Significance of Strain in Formulation in Theory of Solid Mechanics

The basic theory of solid mechanics was deemed complete circa 1860 when St. Venant provided the strain formulation or the field compatibility condition. The strain formulation was incomplete. The missing portion has been formulated and identified as the boundary compatibility condition (BCC). The BCC, derived through a variational formulation, has been verified through integral theorem and solution of problems. The BCC, unlike the field counterpart, do not trivialize when expressed in displacements. Navier s method and the stiffness formulation have to account for the extra conditions especially at the inter-element boundaries in a finite element model. Completion of the strain formulation has led to the revival of the direct force calculation methods: the Integrated Force Method (IFM) and its dual (IFMD) for finite element analysis, and the completed Beltrami-Michell formulation (CBMF) in elasticity. The benefits from the new methods in elasticity, in finite element analysis, and in design optimization are discussed. Existing solutions and computer codes may have to be adjusted for the compliance of the new conditions. Complacency because the discipline is over a century old and computer codes have been developed for half a century can lead to stagnation of the discipline.

Patnaik, Surya N.↗
NASA NTRS1975

Analytical theory for artificial satellites

A theory for generating segmented ephemerides is discussed as a means for fast generation and simple retrieval of nominal orbit data. Over a succession of finite intervals of time, the orbit is represented by a best approximation expressed by Chebyshev polynomials. Storage of coefficients tables for Chebyshev polynomials is seen as a method to reduce data and decrease transmission costs. A general algorithm was constructed and computer programs were designed. The possibility of storing an ephemeris for a few days in the on-board computer, or in microprocessors attached to the data collectors is suggested.

Deprit, A.↗
NASA NTRS1976

Calculation of subsonic and supersonic steady and unsteady aerodynamic forces using velocity potential aerodynamic elements

Expressions for calculation of subsonic and supersonic, steady and unsteady aerodynamic forces are derived, using the concept of aerodynamic elements applied to the downwash velocity potential method. Aerodynamic elements can be of arbitrary out of plane polygon shape, although numerical calculations are restricted to rectangular elements, and to the steady state case in the supersonic examples. It is suggested that the use of conforming, in place of rectangular elements, would give better results. Agreement with results for subsonic oscillating T tails is fair, but results do not converge as the number of collocation points is increased. This appears to be due to the form of expression used in the calculations. The methods derived are expected to facilitate automated flutter analysis on the computer. In particular, the aerodynamic element concept is consistent with finite element methods already used for structural analysis. The method is universal for the complete Mach number range, and, finally, the calculations can be arranged so that they do not have to be repeated completely for every reduced frequency.

Haviland, J. K.↗
NASA NTRS1981

Effect of Surface Roughness on Hydrodynamic Bearings

A theoretical analysis on the performance of hydrodynamic oil bearings is made considering surface roughness effect. The hydrodynamic as well as asperity contact load is found. The contact pressure was calculated with the assumption that the surface height distribution was Gaussian. The average Reynolds equation of partially lubricated surface was used to calculate hydrodynamic load. An analytical expression for average gap was found and was introduced to modify the average Reynolds equation. The resulting boundary value problem was then solved numerically by finite difference methods using the method of successive over relaxation. The pressure distribution and hydrodynamic load capacity of plane slider and journal bearings were calculated for various design data. The effects of attitude and roughness of surface on the bearing performance were shown. The results are compared with similar available solution of rough surface bearings. It is shown that: (1) the contribution of contact load is not significant; and (2) the hydrodynamic and contact load increase with surface roughness.

Majumdar, B. C.↗
NASA NTRS1973

Analysis of three-component aeromagnetic data

The numerical method of obtaining three field components from total field measurements, using double Fourier series expansion, is presented. The expressions for moments of the anomalous field components over a finite area are given. The magnitude and direction of the magnetization vector indicate that the vertical component of the magnetic field calculated from total field observations is more accurate at higher geomagnetic latitudes than at lower latitudes. The opposite is true for the horizontal components. The error in determining the magnetization vector directions are significantly large over most of the range of variation of declination and inclination of the vector, demonstrating the practical limitations of computing field components from total field data even under the best of conditions.

Bhattacharyya, B. K.↗
NASA NTRS1975

Fully unsteady subsonic and supersonic potential aerodynamics for complex aircraft configurations for flutter applications

A general theory for study, oscillatory or fully unsteady potential compressible aerodynamics around complex configurations is presented. Using the finite-element method to discretize the space problem, one obtains a set of differential-delay equations in time relating the potential to its normal derivative which is expressed in terms of the generalized coordinates of the structure. For oscillatory flow, the motion consists of sinusoidal oscillations around a steady, subsonic or supersonic flow. For fully unsteady flow, the motion is assumed to consist of constant subsonic or supersonic speed for time t or = 0 and of small perturbations around the steady state for time t 0.

Tseng, K.↗
NASA NTRS1976

Dynamic ductile fracture of a central crack

A central crack, symmetrically growing at a constant speed in a two dimensional ductile material subject to uniform tension at infinity, is investigated using the integral transform methods. The crack is assumed to be the Dugdale crack, and the finite stress condition at the crack tip is satisfied during the propagation of the crack. Exact expressions of solution are obtained for the finite stress condition at the crack tip, the crack shape, the crack opening displacement, and the energy release rate. All those expressions are written as the product of explicit dimensional quantities and a nondimensional dynamic correction function. The expressions reduce to the associated static results when the crack speed tends to zero, and the nondimensional dynamic correction functions were calculated for various values of the parameter involved.

Tsai, Y. M.↗
NASA NTRS1988

Kuhn-Tucker optimization based reliability analysis for probabilistic finite elements

The fusion of probability finite element method (PFEM) and reliability analysis for fracture mechanics is considered. Reliability analysis with specific application to fracture mechanics is presented, and computational procedures are discussed. Explicit expressions for the optimization procedure with regard to fracture mechanics are given. The results show the PFEM is a very powerful tool in determining the second-moment statistics. The method can determine the probability of failure or fracture subject to randomness in load, material properties and crack length, orientation, and location.

Liu, W. K.↗
NASA NTRS1990

Numerical derivative techniques for trajectory optimization

The adoption of robust numerical optimization techniques in trajectory simulation programs has resulted in powerful design and analysis tools. These trajectory simulation/optimization programs are widely used, and a representative list includes the GTS system, the POST program, and newer collocation methods such as OTIS and FONPAC. All of these programs rely on optimization algorithms which require objective function and constraint gradient data during the iteration process. However, most trajectory optimization problems lack simple analytical expressions for these derivatives. In the general case a function evaluation involves integrating aerodynamic, propulsive, and gravity forces over multiple trajectory phases with complex control models. With the newer collocation methods, the integration is replaced by defect constraints and cubic approximations for the state. While analytic gradient expressions can sometimes be derived for trajectory optimization problems, the derivation is cumbersome, time consuming, and prone to mistakes. Fortunately, an alternate method exists for the gradient evaluation, namely finite difference approximations. In this paper some finite difference gradient techniques developed for use with the GTS system are presented. These techniques include methods for computing first and second partial derivatives of single and multiple sets of functions. A key feature of these methods is an error control mechanism which automatically adjusts the perturbation size to obtain accurate derivative values.

Hallman, Wayne P.↗
NASA NTRS1991

A recursive approach to the equations of motion for the maneuvering and control of flexible multi-body systems

Interest lies in a mathematical formulation capable of accommodating the problem of maneuvering a space structure consisting of a chain of articulated flexible substructures. Simultaneously, any perturbations from the 'rigid body' maneuvering and any elastic vibration must be suppressed. The equations of motion for flexible bodies undergoing rigid body motions and elastic vibrations can be obtained conveniently by means of Lagrange's equations in terms of quasi-coordinates. The advantage of this approach is that it yields equations in terms of body axes, which are the same axes that are used to express the control forces and torques. The equations of motion are nonlinear hybrid differential quations. The partial differential equations can be discretized (in space) by means of the finite element method or the classical Rayleigh-Ritz method. The result is a set of nonlinear ordinary differential equations of high order. The nonlinearity can be traced to the rigid body motions and the high order to the elastic vibration. Elastic motions tend to be small when compared with rigid body motions.

Kwak, Moon K.↗
NASA NTRS1962

BUCKLING OF INITIALLY IMPERFECT AXIALLY COMPRESSED CYLINDRICAL SHELLS

The effects of imperfections on the buckling strength of pressurized cylinders under axial compression are studied. The imperfection factor is considered in the finite-deflection compatibility equation, as well as in the equilibrium equation. A new relation to express the imperfection as an explicit exponential function of pressure and radius thickness ratio is proposed. A method for finding the exponential function is described, and the solution for the critical stress is found in a fairly simple form. The relation of decrease of stability with respect to magnitude of imperfections if found.

CYLINDRICAL SHELL↗
NASA NTRS1984

Automatic finite element generators

The design and implementation of a software system for generating finite elements and related computations are described. Exact symbolic computational techniques are employed to derive strain-displacement matrices and element stiffness matrices. Methods for dealing with the excessive growth of symbolic expressions are discussed. Automatic FORTRAN code generation is described with emphasis on improving the efficiency of the resultant code.

Wang, P. S.↗
NASA NTRS2010

Measuring Crack Length in Coarse Grain Ceramics

Due to a coarse grain structure, crack lengths in precracked spinel specimens could not be measured optically, so the crack lengths and fracture toughness were estimated by strain gage measurements. An expression was developed via finite element analysis to correlate the measured strain with crack length in four-point flexure. The fracture toughness estimated by the strain gaged samples and another standardized method were in agreement.

Salem, Jonathan A.↗
NASA NTRS1986

Sublaminate analysis of interlaminar fracture in composites

A simple analysis method based upon a transverse shear deformation theory and a sublaminate approach is utilized to analyze a mixed-mode edge delamination specimen. The analysis provides closed form expressions for the interlaminar shear stresses ahead of the crack, the total energy release rate, and the energy release rate components. The parameters controlling the behavior are identified. The effect of specimen stacking sequence and delamination interface on the strain energy release rate components is investigated. Results are compared with a finite element simulation for reference. The simple nature of the method makes it suitable for preliminary design analyses which require a large number of configurations to be evaluated quickly and economically.

Armanios, E. A.↗
NASA NTRS1986

Design component method for sensitivity analysis of built-up structures

A 'design component method' that provides a unified and systematic organization of design sensitivity analysis for built-up structures is developed and implemented. Both conventional design variables, such as thickness and cross-sectional area, and shape design variables of components of built-up structures are considered. It is shown that design of components of built-up structures can be characterized and system design sensitivity expressions obtained by simply adding contributions from each component. The method leads to a systematic organization of computations for design sensitivity analysis that is similar to the way in which computations are organized within a finite element code.

Choi, Kyung K.↗
NASA NTRS1996

Comparison of Finite-Difference and Analytical Microwave Calculation Methods

Normal modes and power absorption distribution in microwave cavities containing lossy dielectric samples were calculated for problems of interest in materials processing. The calculations were performed both using a commercially available finite-difference electromagnetic solver and by numerical evaluation of exact analytic expressions.

power↗
NASA NTRS1973

High speed flow past wings

The analytical solution to the transonic small perturbation equation which describes steady compressible flow past finite wings at subsonic speeds can be expressed as a nonlinear integral equation with the perturbation velocity potential as the unknown function. This known formulation is substituted by a system of nonlinear algebraic equations to which various methods are applicable for its solution. Due to the presence of mathematical discontinuities in the flow solutions, however, a main computational difficulty was to ensure uniqueness of the solutions when local velocities on the wing exceeded the speed of sound. For continuous solutions this was achieved by embedding the algebraic system in an one-parameter operator homotopy in order to apply the method of parametric differentiation. The solution to the initial system of equations appears then as a solution to a Cauchy problem where the initial condition is related to the accompanying incompressible flow solution. In using this technique, however, a continuous dependence of the solution development on the initial data is lost when the solution reaches the minimum bifurcation point. A steepest descent iteration technique was therefore, added to the computational scheme for the calculation of discontinuous flow solutions. Results for purely subsonic flows and supersonic flows with and without compression shocks are given and compared with other available theoretical solutions.

Norstrud, H.↗
NASA NTRS1995

Analytic Sensitivity and Approximation of Skin Buckling Constraints in Wing-Shape Synthesis

Explicit expressions for terms of the stiffness and geometric stiffness matrices are derived for the buckling analysis of trapezoidal fiber composite wing skin panels. The formulation is based on Ritz analysis using simple polynomials, and leads to explicit expressions for the analytic sensitivities of the stiffness and geometric stiffness matrices with respect to layer thickness, fiber directions, and panel shape. Integration with wing box analysis using either the equivalent plate approach or the finite element method, makes it possible to obtain sensitivities of panel buckling constraints with respect to wing planform shape or locations of internal ribs and spars. The analytic sensitivities are used to construct approximations of panel buckling constraints for integrated wing/panel design synthesis.

Livne, Eli↗