SEARCH · Search NASA
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.
Closed form expressions are derived for the position and velocity of a spacecraft during a finite burn using the method involving the theory of asymptotic expansion of the optimal impulsive solution. The small parameter is given by the reciprocal of the mass flow rate. The expansion is given in terms of the impulsive solution and holds up through third-order for the velocity and fourth-order for the position.
Grid quality solutions of a wide variety of equations can be generated to fit reasonable functional boundary conditions in two dimensions using very coarse rectangular finite elements. Some wavy operators (with no natural variational expression) were tried in order to demonstrate the method's versatility.
This paper summarizes ongoing research on the sensitivity analysis of dynamic aeroelastic response of wings. Two approaches are being used to express the unsteady aerodynamic loads: (1) the frequency-domain approach, and (2) the state-space approach. The frequency-domain approach is demonstrated on a three-dimensional box wing and the state-space domain approach is demonstrated on a simple two-dimensional sectional model. Three different methods are used to find the sensitivities: (1) a purely finite difference approach, (2) a semi-analytical approach in which an analytical expression is used for calculating the sensitivity of an eigenvalue of the complex valued aeroelastic matrix, however, the derivatives of the components of the matrix are obtained using finite difference, and (3) a semi-analytic approach that differs from (2) in the sense that the sensitivity of the aerodynamic matrix is now obtained analytically. A good agreement is seen between the three sets of results. For the two-dimensional sectional model the results for the sensitivities of the flutter speed with respect to the various parameters are obtained.
Method of estimating design survivability of structural component incorporates finite-element and probabilistic properties of materials. Involves evaluation of design parameters through direct comparisons of survivability of component expressed in terms of percentages of like components that survive at various lifetimes. Probabilistic properties of materials, given in terms of Weibull parameters, coupled with stress field computed by finite-element analysis to determine fatigue life based on initiation of cracks. Method applied to rotating disk containing bolt holes, representative of disks used in aerospace propulsion turbines. Also used in early stages of design process to optimize life-based designs, reducing testing of full-sized components needed to validate designs.
Numerical methods for mixed boundary value problem of axisymmetric shells of revolution, using truncated series expansion and finite difference expressions
A mixed finite-difference scheme is presented for the free-vibration analysis of simply supported closed noncircular cylindrical shells. The problem is formulated in terms of eight first-order differential equations in the circumferential coordinate which possess a symmetric coefficient matrix and are free of the derivatives of the elastic and geometric characteristics of the shell. In the finite-difference discretization, two interlacing grids are used for the different fundamental unknowns in such a way as to avoid averaging in the difference-quotient expressions used for the first derivative. The resulting finite-difference equations are symmetric. The inverse-power method is used for obtaining the eigenvalues and eigenvectors.
The hydrodynamic impact of a falling body upon a viscous incompressible fluid was investigated by numerically solving the equations of motion. Initially the mathematical model simulated the axisymmetric impact of a rigid right circular cylinder upon the initially quiescent free surface of a fluid. A compressible air layer exists between the falling cylinder and the liquid free surface. The mathematical model was developed by applying the Navier-Stokes equations to the incompressible air layer and the incompressible fluid. Assuming the flow to be one dimensional within the air layer, the average velocity, pressure and density distributions were calculated. The liquid free surface was allowed to deform as the air pressure acting on it increases. For the liquid the normalized equations were expressed in two-dimensional cylindrical coordinates. The governing equations for the air layer and the liquid were expressed in finite difference form and solved numerically. For the liquid a modified version of the Marker-and-Cell method was used. The mathematical model has been reexamined and a new approach has recently been initiated. Essentially this consists of examining the impact of an inclined plate onto a quiesent water surface with the equations now formulated in cartesian coordinates.
A method for calculating a component's design survivability by incorporating finite element analysis and probabilistic material properties was developed. The method evaluates design parameters through direct comparisons of component survivability expressed in terms of Weibull parameters. The analysis was applied to a rotating disk with mounting bolt holes. The highest probability of failure occurred at, or near, the maximum shear stress region of the bolt holes. Distribution of failure as a function of Weibull slope affects the probability of survival. Where Weibull parameters are unknown for a rotating disk, it may be permissible to assume Weibull parameters, as well as the stress-life exponent, in order to determine the disk speed where the probability of survival is highest.
A method for estimating a component's design survivability by incorporating finite element analysis and probabilistic material properties was developed. The method evaluates design parameters through direct comparisons of component survivability expressed in terms of Weibull parameters. The analysis was applied to a rotating disk with mounting bolt holes. The highest probability of failure occurred at, or near, the maximum shear stress region of the bolt holes. Distribution of material failure as a function of Weibull slope affects the probability of survival. Where Weibull parameters are unknown for a rotating disk, it may be permissible to assume Weibull parameters, as well as the stress-life exponent, in order to determine the qualitative effect of disk speed on the probability of survival.
The ability to accurately predict the thermomechanical deformation response of advanced composite materials continues to play an important role in the development of these strategic materials. Analytical models that predict the effective behavior of composites are used not only by engineers in performing structural analysis of large-scale composite components but also by material scientists in developing new material systems. For an analytical model to fulfill these two distinct functions, it must be based on a micromechanics approach that uses physically based deformation and life constitutive models, and it must allow one to generate the average (macro) response of a composite material given the properties of the individual constituents and their geometric arrangement. Only then can such a model be used by a material scientist to investigate the effect of different deformation mechanisms on the overall response of the composite and, thereby, identify the appropriate constituents for a given application. However, if a micromechanical model is to be used in a large-scale structural analysis it must be (1) computationally efficient, (2) able to generate accurate displacement and stress fields at both the macro and micro level, and (3) compatible with the finite element method. In addition, new advancements in processing and fabrication techniques now make it possible to engineer the architectures of these advanced composite systems. Full utilization of these emerging manufacturing capabilities require the development of a computationally efficient micromechanics analysis tool that can accurately predict the effect of microstructural details on the internal and macroscopic behavior of composites. Computational efficiency is required because (1) a large number of parameters must be varied in the course of engineering (or designing) composite materials and (2) the optimization of a material's microstructure requires that the micromechanics model be integrated with optimization algorithms. From this perspective, analytical approaches that produce closed-form expressions which describe the effect of a material's internal architecture on the overall material behavior are preferable to numerical methods such as the finite element or finite difference schemes.
In the dynamics modeling of a flexible structure, finite element analysis employs reduction techniques, such as Guyan's reduction, to remove some of the insignificant physical coordinates, thus producing a dynamics model that has smaller mass and stiffness matrices. But this reduction is limited in the sense that it removes certain degrees of freedom at a node points themselves in the model. From the standpoint of linear control design, the resultant model is still too large despite the reduction. Thus, some form of the model reduction is frequently used in control design by approximating a large dynamical system with a fewer number of state variables. However, a problem arises from the placement of sensors and actuators in the reduced model, because a model usually undergoes, before being reduced, some form of coordinate transformations that do not preserve the physical meanings of the states. To correct such a problem, a method is developed that expresses a reduced model in terms of a subset of the original states. The proposed method starts with a dynamic model that is originated and reduced in finite element analysis. Then the model is converted to the state space form, and reduced again by the internal balancing method. At this point, being in the balanced coordinate system, the states in the reduced model have no apparent resemblance to those of the original model. Through another coordinate transformation that is developed, however, this reduced model is expressed by a subset of the original states.
An important consideration in the global local finite-element method (GLFEM) is the availability of global functions for the given problem. The role and mathematical requirements of these global functions in a GLFEM analysis of localized stress states in prismatic structures are discussed. A method is described for determining these global functions. Underlying this method are theorems due to Toupin and Knowles on strain energy decay rates, which are related to a quantitative expression of Saint-Venant's principle. It is mentioned that a mathematically complete set of global functions can be generated, so that any arbitrary interface condition between the finite element and global subregions can be represented. Convergence to the true behavior can be achieved with increasing global functions and finite-element degrees of freedom. Specific attention is devoted to mathematically two-dimensional and three-dimensional prismatic structures. Comments are offered on the GLFEM analysis of NASA flat panel with a discontinuous stiffener. Methods for determining global functions for other effects are also indicated, such as steady-state dynamics and bodies under initial stress.
We introduce an exact-two-component complete active space self-consistent-field (X2C-CASSCF) method formulated under the restricted-magnetic-balance condition. This framework allows for the nonperturbative treatment of static magnetic fields using gauge-including atomic orbitals (GIAOs). The GIAO-X2C-CASSCF methodology effectively captures all microstates within the same 2J + 1-degenerate manifold and their splitting in a static magnetic field, which are not accessible through single-reference-based methods. We also present mathematical recursive expressions for evaluating one-electron relativistic integrals by using GIAOs in the presence of a finite magnetic field. Benchmark studies include oxygen and nitrogen K-edge X-ray magnetic circular dichroism spectroscopy (XMCD) for closed-shell organic compounds, as well as L-edge XMCD spectroscopy for the high-spin open-shell transition metal ion Mn 2+ and the tetrahedral Mn(II)O 4 6– complex.
To eliminate the need to use finite-element modeling for structure shape predictions, a new method was invented. This method is to use the Displacement Transfer Functions to transform the measured surface strains into deflections for mapping out overall structural deformed shapes. The Displacement Transfer Functions are expressed in terms of rectilinearly distributed surface strains, and contain no material properties. This report is to apply the patented method to the shape predictions of non-symmetrically loaded slender curved structures with different curvatures up to a full circle. Because the measured surface strains are not available, finite-element analysis had to be used to analytically generate the surface strains. Previously formulated straight-beam Displacement Transfer Functions were modified by introducing the curvature-effect correction terms. Through single-point or dual-point collocations with finite-elementgenerated deflection curves, functional forms of the curvature-effect correction terms were empirically established. The resulting modified Displacement Transfer Functions can then provide quite accurate shape predictions. Also, the uniform straight-beam Displacement Transfer Function was applied to the shape predictions of a section-cut of a generic capsule (GC) outer curved sandwich wall. The resulting GC shape predictions are quite accurate in partial regions where the radius of curvature does not change sharply.
A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.
Through the method of truncation functions, the oceanic geoid undulation is divided into two constituents: an inner zone contribution expressed as an integral of surface gravity disturbances over a spherical cap; and an outer zone contribution derived from a finite set of potential harmonic coefficients. Global, average error estimates are formulated for undulation differences, thereby providing accuracies for a relative geoid. The error analysis focuses on the outer zone contribution for which the potential coefficient errors are modeled. The method of computing undulations based on gravity disturbance data for the inner zone is compared to the similar, conventional method which presupposes gravity anomaly data within this zone.
A computer program was developed to account approximately for the effects of finite wing thickness in the transonic potential flow over an oscillating wing of finite span. The program is based on the original sonic-box program for planar wing which was previously extended to include the effects of the swept trailing edge and the thickness of the wing. Account for the nonuniform flow caused by finite thickness is made by application of the local linearization concept. The thickness effect, expressed in terms of the local Mach number, is included in the basic solution to replace the coordinate transformation method used in the earlier work. Calculations were made for a delta wing and a rectangular wing performing plunge and pitch oscillations, and the results were compared with those obtained from other methods. An input quide and a complete listing of the computer code are presented.
The algorithms that comprise the Variational Trajectory Optimization Tool Set (VTOTS) package are briefly described. The VTOTS is a software package for solving nonlinear constrained optimal control problems from a wide range of engineering and scientific disciplines. The VTOTS package was specifically designed to minimize the amount of user programming; in fact, for problems that may be expressed in terms of analytical functions, the user needs only to define the problem in terms of symbolic variables. This version of the VTOTS does not support tabular data; thus, problems must be expressed in terms of analytical functions. The VTOTS package consists of two methods for solving nonlinear optimal control problems: a time-domain finite-element algorithm and a multiple shooting algorithm. These two algorithms, under the VTOTS package, may be run independently or jointly. The finite-element algorithm generates approximate solutions, whereas the shooting algorithm provides a more accurate solution to the optimization problem. A user's manual, some examples with results, and a brief description of the individual subroutines are included.