SEARCH · Search NASA

Results for “Variational formulation”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Geometrically nonlinear analysis of layered composite shells

Two kinds of finite-element analyses are developed for the geometrically nonlinear study of the large deformations in laminated composite structures, especially shells. The first kind of finite-element analysis utilizes the general incremental variational formulation as well as the total Lagrangian description of motion, and a three-dimensional degenerate element is adopted. The second kind of analysis employs a formulation based on deformable shell theory, and the plate-bending element is used. Numerical results for bending are presented for five plate and shell structures of isotropic as well as orthotropic composition, including an isotropic cylindrical panel with uniform loading and a laminated cylindrical panel with uniform loading. The results obtained using these analyses are found to be in good agreement with those available in the literature.

Chao, W. C.↗

NASA NTRS1986

Transformation of two and three-dimensional regions by elliptic systems

Grid smoothing and orthogonalization procedures were developed and implemented in the construction of two and three dimensional grids. The procedures are based on the variational methods of grid generation. The two-dimensional examples were computed using the MSU IRIS Graphics Workstation. It was demonstrated that the elliptic grid generation equations, with arbitrary forcing functions, can be solved, in their variational formulation, using a gradient method. Since gradient methods have a global convergence property, the divergence problems often encountered when using SOR iterative methods can be avoided. It is not to be concluded, however, that SOR methods should be abandoned, since gradient methods tend to converge very slowly. In fact, slow convergence was the major problem encountered in the three-dimensional grids. Further progress was made on the continuing effort to develop conservative interpolation formulas for overlapping grids.

Mastin, C. Wayne↗

NASA NTRS1988

Reduced order models for nonlinear aerodynamics

Reduced order models are needed for reliable, efficient and accurate prediction of aerodynamic forces to analyze fluid structure interaction problems in turbomachinery, including propfans. Here, a finite difference, time marching Navier-Stokes code is validated for unsteady airfoil motion by comparing results with those from classical potential flow. The Navier-Stokes code is then analyzed for calculation of primitive and exact estimates of eigenvalues and eigenvectors associated with fluid-airfoil interaction. A variational formulation for the Euler equations and Navier-Stokes equations will be the basis for reduction of order through an eigenvector transformation.

Mahajan, Aparajit J.↗

NASA NTRS1989

Vector potential methods

Vector potential and related methods, for the simulation of both inviscid and viscous flows over aerodynamic configurations, are briefly reviewed. The advantages and disadvantages of several formulations are discussed and alternate strategies are recommended. Scalar potential, modified potential, alternate formulations of Euler equations, least-squares formulation, variational principles, iterative techniques and related methods, and viscous flow simulation are discussed.

Hafez, M.↗

NASA NTRS1992

Boundary-conforming mapping applied to computations of highly deformed solidification interfaces

A new boundary-conforming mapping is developed for the calculation of highly deformed cellular solidification interfaces in a model of directional solidification of a binary alloy. The mapping is derived through a variational formulation that is designed so that the grid penetrates the grooves between cells along the interface without causing a loss of ellipticity of the mapping equations. A finite element/Newton method is presented for simultaneous solution of the free boundary problem described by the solutal model of directional solidification and the mapping equations. Results are compared to previous calculations and demonstrate the importance of accurate representation of the interface shape for understanding the solution structure.

Tsiveriotis, K.↗

NASA NTRS1992

Identification and feedback control in structures with piezoceramic actuators

In this lecture we give fundamental well-posedness results for a variational formulation of a class of damped second order partial differential equations with unbounded input or control coefficients. Included as special cases in this class are structures with piezoceramic actuators. We consider approximation techniques leading to computational methods in the context of both parameter estimation and feedback control problems for these systems. Rigorous convergence results for parameter estimates and feedback gains are discussed.

Banks, H. T.↗

NASA NTRS1993

A new interface element for connecting independently modeled substructures

A new interface element based on the hybrid variational formulation is presented and demonstrated. The element provides a means of connecting independently modeled substructures whose nodes along the common boundary need not be coincident. The interface element extends previous work to include connecting an arbitrary number of substructures, the use of closed and generally curved interfaces, and the use of multiple, possibly nested, interfaces. Several applications of the element are presented and aspects of the implementation are discussed.

Ransom, Jonathan B.↗

NASA NTRS1993

A nonlinear theory for spinning anisotropic beams using restrained warping functions

A geometrically nonlinear theory is developed for spinning anisotropic beams having arbitrary cross sections. An assumed displacement field is developed using the standard 3D kinematics relations to describe the global beam behavior supplemented with an additional field that represents the local deformation within the cross section and warping out of the cross section plane. It is assumed that the magnitude of this additional field is directly proportional to the local stress resultants. In order to take into account the effects of boundary conditions, a restraining function is introduced. This function plays the role of reducing the amount of free warping deformation throughout the field due to the restraint of the cross section(s) at the end(s) of the beam, e.g., in the case of a cantilever beam. Using a developed ordering scheme, the nonlinear strains are calculated to the third order. The FEM is developed using the weak form variational formulation. Preliminary interesting numerical results have been obtained that indicate the role of the restraining function in the case of a cantilever beam with circular cross section. These results are for the cases of a tip displacement (static) and free vibration studies for both isotropic and anisotropic materials with varied fiber orientations.

Ie, C. A.↗

NASA NTRS1988

On the inclusion principle and spillover effect in the identification of distributed structures

Distributed structures are Distributed-Parameter Systems (DPS) and must be modeled by large-dimensional systems. Practical limitations require that only a subset of the modes be actively controlled. Control of the residual (uncontrolled) modes is known as control spillover, and it can degrade the system performance. The contamination of the sensor signals by the residual modes is known as observation spillover and it can destabilize the system. An analogous effect exists in the identification problem. Distributed structures must be represented by discretized (finite-dimensional) models or lumped-parameter systems. The contamination of the measurements by the residual (unmodeled) modes is referred to as identification spillover, and it is the object of this study. To study the effects of identification spillover, we use a variational formulation to obtain the natural frequencies and mode shapes. This formulation has direct application to DPS. The identification technique consists of finding the stationary values of a pseudo-Rayleigh quotient.

Norris, Mark. A.↗

NASA NTRS1993

An adaptive grid with directional control

An adaptive grid generator for adaptive node movement is here derived by combining a variational formulation of Winslow's (1981) variable-diffusion method with a directional control functional. By applying harmonic-function theory, it becomes possible to define conditions under which there exist unique solutions of the resulting elliptic equations. The results obtained for the grid generator's application to the complex problem posed by the fluid instability-driven magnetic field reconnection demonstrate one-tenth the computational cost of either a Eulerian grid or an adaptive grid without directional control.

Brackbill, J. U.↗

NASA NTRS1994

Cross-surface interface element for coupling built-up structural subdomains

A new finite element for coupling built-up shell substructures is presented. The present work extends the hybrid variational formulation of the interface element developed by Aminpour and Ransom to permit coupling between two intersecting substructures. Designed for the assembly of independently built-up finite element models, this technique provides a level of modeling flexibility previously unavailable.

Davila, C. G.↗

NASA NTRS1995

Multiple burn fuel-optimal orbit transfers: Numerical trajectory computation and neighboring optimal feedback guidance

This report describes current work in the numerical computation of multiple burn, fuel-optimal orbit transfers and presents an analysis of the second variation for extremal multiple burn orbital transfers as well as a discussion of a guidance scheme which may be implemented for such transfers. The discussion of numerical computation focuses on the use of multivariate interpolation to aid the computation in the numerical optimization. The second variation analysis includes the development of the conditions for the examination of both fixed and free final time transfers. Evaluations for fixed final time are presented for extremal one, two, and three burn solutions of the first variation. The free final time problem is considered for an extremal two burn solution. In addition, corresponding changes of the second variation formulation over thrust arcs and coast arcs are included. The guidance scheme discussed is an implicit scheme which implements a neighboring optimal feedback guidance strategy to calculate both thrust direction and thrust on-off times.

Chuang, C.-H.↗

NASA NTRS1993

A global/local analysis method for treating details in structural design

A method for analyzing global/local behavior of plate and shell structures is described. In this approach, a detailed finite element model of the local region is incorporated within a coarser global finite element model. The local model need not be nodally compatible (i.e., need not have a one-to-one nodal correspondence) with the global model at their common boundary; therefore, the two models may be constructed independently. The nodal incompatibility of the models is accounted for by introducing appropriate constraint conditions into the potential energy in a hybrid variational formulation. The primary advantage of this method is that the need for transition modeling between global and local models is eliminated. Eliminating transition modeling has two benefits. First, modeling efforts are reduced since tedious and complex transitioning need not be performed. Second, errors due to the mesh distortion, often unavoidable in mesh transitioning, are minimized by avoiding distorted elements beyond what is needed to represent the geometry of the component. The method is applied reduced to a plate loaded in tension and transverse bending. The plate has a central hole, and various hole sixes and shapes are studied. The method is also applied to a composite laminated fuselage panel with a crack emanating from a window in the panel. While this method is applied herein to global/local problems, it is also applicable to the coupled analysis of independently modeled components as well as adaptive refinement.

Aminpour, Mohammad A.↗

NASA NTRS1995

Characteristics of Three-Node Smoothing Element under Penalty Constraints

The project is based upon research of Tessler et al. (1994) on an improved variational formulation for post-processing stress predictions in Finite Element Analysis. The methodology, called Smoothing Element Analysis (SEA), employs a three-node smoothing finite element. The present effort focused on verifying the basic constant strain criterion for the three-node smoothing element subject to a set of internal penalty constraints. The convergence characteristics of the element are assessed by first deriving the constrained form of the assumed element stress and stress gradient fields, and then by verifying the validity of the constant strain criterion once the element penalty constraints are explicitly imposed. The analytical investigation is carried out with the use of the symbolic manipulation code 'Mathematica'.

Sobel, Seth S.↗

NASA NTRS1997

Interface Technology for Geometrically Nonlinear Analysis of Multiple Connected Subdomains

Interface technology for geometrically nonlinear analysis is presented and demonstrated. This technology is based on an interface element which makes use of a hybrid variational formulation to provide for compatibility between independently modeled connected subdomains. The interface element developed herein extends previous work to include geometric nonlinearity and to use standard linear and nonlinear solution procedures. Several benchmark nonlinear applications of the interface technology are presented and aspects of the implementation are discussed.

Ransom, Jonathan B.↗

NASA NTRS1999

Development of Analysis Methods for Designing with Composites

The project involved the development of new analysis methods to achieve efficient design of composite structures. We developed a complex variational formulation to analyze the in-plane and bending coupling response of an unsymmetrically laminated plate with an elliptical cutout subjected to arbitrary edge loading as shown in Figure 1. This formulation utilizes four independent complex potentials that satisfy the coupled in-plane and bending equilibrium equations, thus eliminating the area integrals from the strain energy expression. The solution to a finite geometry laminate under arbitrary loading is obtained by minimizing the total potential energy function and solving for the unknown coefficients of the complex potentials. The validity of this approach is demonstrated by comparison with finite element analysis predictions for a laminate with an inclined elliptical cutout under bi-axial loading.The geometry and loading of this laminate with a lay-up of [-45/45] are shown in Figure 2. The deformed configuration shown in Figure 3 reflects the presence of bending-stretching coupling. The validity of the present method is established by comparing the out-of-plane deflections along the boundary of the elliptical cutout from the present approach with those of the finite element method. The comparison shown in Figure 4 indicates remarkable agreement. The details of this method are described in a manuscript by Madenci et al. (1998).

Madenci, E.↗

NASA NTRS2000

Bolted Double-Lap Composite Joints Under Mechanical and Thermal Loading

This study concerns the determination of the contact stresses and contact region around bolt holes and the bolt load distribution in single- and double-lap joints of composite laminates with arbitrarily located bolts under general mechanical loading conditions and uniform temperature change. The unknown contact stress distribution and contact region between the bolt and laminates and the interaction among the bolts require the bolt load distribution, as well as the contact stresses, to be as part of the solution. The present method is based on the complex potential theory and the variational formulation in order to account for bolt stiffness, bolt-hole clearance, and finite geometry of the composite laminates.

Kradinov, V.↗

NASA NTRS2000

Examples of Unsteady Combustion in Non-Metallized Propellants

For over thirty years the T-Burner has been the de facto standard for determining the unstable response of solid propellants. This paper reviews some of the different approaches that have been used, and briefly summarizes advantages and disadvantages of the different techniques. Data from a variety of sources have been reviewed and summarized showing the effect of various test conditions and propellant formulation variations. Much of the data taken in the 60's and early 70's were with so-called research propellants, with relatively low solids loading and usually monomodal AP as the oxidizer. Later data were obtained with propellants containing solids loading more typical and practical propellants. Several distinct differences were observed in the data for the different types of propellants. Data varying AP particle size distributions indicate that either additional fine AP or finer AP usually increases the propellant response. A significant reduction in the response can be achieved by reducing the amount of fine AP, while holding burning rate constant either by adding a catalyst or by varying the overall AP size distribution. A significant amount of data has been obtained for reduced smoke propellants containing stability additives. The data indicate that the additives provide significant particle damping, and most often, also reduce the propellant response. Data from double base propellants indicate a much broader response over a wide frequency range. This implies that double base propellants should be more unstable in tangential modes than composite propellants. This trend seems to be consistent with observed motor data. Measuring the response of metallized propellants has been a challenge due to the large damping due to metal oxide particles in the gas. Data from metallized propellants are not presented here, but will be summarized in a future paper.

Beckstead, M. W.↗