SEARCH · Search NASA

Results for “LAGRANGE MULTIPLIER”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Staggered solution procedures for multibody dynamics simulation

A computational procedure is presented for the direct integration of the multibody dynamical (MBD) equations with constraints. Numerical experiments conducted with this MBD procedure indicate that it yields robust solutions in cases where the step size furnishes more than 20 samples for the period of highest apparent response frequency in a given multibody system. The MBD solution procedure is implemented in two separate modules: the generalized coordinated solver, 'CINT', and the constraint Lagrange multiplier solver, 'LINT'.

Park, K. C.↗

NASA NTRS1990

Explicit-implicit staggered procedure for multibody dynamics analysis

An explicit-implicit staggered time-integration procedure is presented for the solution of multibody dynamical equations involving large rotations and constraints. The algorithms adopts a two-stage modification of the central difference algorithm for integrating the translational coordinates and the angular velocity vector, and the midpoint implicit algorithm to solve the kinematical relation in terms of the Euler parameters for updating the angular orientations. The Lagrange multipliers to enforce the system constraints are obtained by implicitly integrating a parabolically regularized differential equation for the multipliers. The performance of the present procedure has been evaluated to applying the procedure to solve several sample problems. The results indicate that the procedure is robust in dealing with a variety of constraints and spatial kinematic motions, hence it is recommended for applications to general multibody dynamics analyses.

Park, K. C.↗

NASA NTRS1991

An active set algorithm for tracing parametrized optima

Optimization problems often depend on parameters that define constraints or objective functions. It is often necessary to know the effect of a change in a parameter on the optimum solution. An algorithm is presented here for tracking paths of optimal solutions of inequality constrained nonlinear programming problems as a function of a parameter. The proposed algorithm employs homotopy zero-curve tracing techniques to track segments where the set of active constraints is unchanged. The transition between segments is handled by considering all possible sets of active constraints and eliminating nonoptimal ones based on the signs of the Lagrange multipliers and the derivatives of the optimal solutions with respect to the parameter. A spring-mass problem is used to illustrate all possible kinds of transition events, and the algorithm is applied to a well-known ten-bar truss structural optimization problem.

Rakowska, J.↗

NASA NTRS1990

Electromagnetic axisymmetric finite elements based on a gauged four-potential variational principle

Electromagnetic finite elements are derived based on a variational principle that uses the electromagnetic four-potential as a primary variable. The Lorentz gage normalization is incorporated as a constraint condition through a Lagrange multiplier field to construct elements suitable for downstream coupling with mechanical and thermal finite elements for the analysis of high-temperature superconductor devices with aerospace applications. The main advantages are: jump discontinuities on interfaces are naturally handled; no a priori approximations are invoked; and the number of degrees of freedom per node remains modest as the problem dimensionality increases. The new elements are tested on two magnetostatic axisymmetric problems. The results are in excellent agreement with analytical solutions and previous solutions for the 1D problem of a conducting infinite wire, in which case the multiplier field has no effect. For materials of widely different permeability, jump conditions are naturally accommodated by the present formulation.

Schuler, J.↗

NASA NTRS1992

Structural optimization of large structural systems by optimality criteria methods

The fundamental concepts of the optimality criteria method of structural optimization are presented. The effect of the separability properties of the objective and constraint functions on the optimality criteria expressions is emphasized. The single constraint case is treated first, followed by the multiple constraint case with a more complex evaluation of the Lagrange multipliers. Examples illustrate the efficiency of the method.

Berke, Laszlo↗

NASA NTRS1992

Aerodynamic design and optimization in one shot

This paper describes an efficient numerical approach for the design and optimization of aerodynamic bodies. As in classical optimal control methods, the present approach introduces a cost function and a costate variable (Lagrange multiplier) in order to achieve a minimum. High efficiency is achieved by using a multigrid technique to solve for all the unknowns simultaneously, but restricting work on a design variable only to grids on which their changes produce nonsmooth perturbations. Thus, the effort required to evaluate design variables that have nonlocal effects on the solution is confined to the coarse grids. However, if a variable has a nonsmooth local effect on the solution in some neighborhood, it is relaxed in that neighborhood on finer grids. The cost of solving the optimal control problem is shown to be approximately two to three times the cost of the equivalent analysis problem. Examples are presented to illustrate the application of the method to aerodynamic design and constraint optimization.

Ta'asan, Shlomo↗

NASA NTRS1984

Equations of motion for control of the SCOLE laboratory experiment

The objectives of this study are listed as follows: (1) to develop Lagrange's equations of motion for the shuttle antenna configuration in orbit; (2) to modify equations using the Lagrange multiplier method to develop equations of motion for the laboratory experiment; and (3) to discuss methods for simulation and control. The equations are presented in graph form.

Meirovitch, L.↗

NASA NTRS1993

Coupled structural, thermal, phase-change and electromagnetic analysis for superconductors, volume 2

Two families of parametrized mixed variational principles for linear electromagnetodynamics are constructed. The first family is applicable when the current density distribution is known a priori. Its six independent fields are magnetic intensity and flux density, magnetic potential, electric intensity and flux density and electric potential. Through appropriate specialization of parameters the first principle reduces to more conventional principles proposed in the literature. The second family is appropriate when the current density distribution and a conjugate Lagrange multiplier field are adjoined, giving a total of eight independently varied fields. In this case it is shown that a conventional variational principle exists only in the time-independent (static) case. Several static functionals with reduced number of varied fields are presented. The application of one of these principles to construct finite elements with current prediction capabilities is illustrated with a numerical example.

Felippa, Carlos A.↗

NASA NTRS1992

Two neural network algorithms for designing optimal terminal controllers with open final time

Multilayer neural networks, trained by the backpropagation through time algorithm (BPTT), have been used successfully as state-feedback controllers for nonlinear terminal control problems. Current BPTT techniques, however, are not able to deal systematically with open final-time situations such as minimum-time problems. Two approaches which extend BPTT to open final-time problems are presented. In the first, a neural network learns a mapping from initial-state to time-to-go. In the second, the optimal number of steps for each trial run is found using a line-search. Both methods are derived using Lagrange multiplier techniques. This theoretical framework is used to demonstrate that the derived algorithms are direct extensions of forward/backward sweep methods used in N-stage optimal control. The two algorithms are tested on a Zermelo problem and the resulting trajectories compare favorably to optimal control results.

Plumer, Edward S.↗

NASA NTRS1992

A robust direct-integration method for rotorcraft maneuver and periodic response

The Newmark-Beta method and the Newton-Raphson iteration scheme are combined to develop a direct-integration method for evaluating the maneuver and periodic-response expressions for rotorcraft. The method requires the generation of Jacobians and includes higher derivatives in the formulation of the geometric stiffness matrix to enhance the convergence of the system. The method leads to effective convergence with nonlinear structural dynamics and aerodynamic terms. Singularities in the matrices can be addressed with the method as they arise from a Lagrange multiplier approach for coupling equations with nonlinear constraints. The method is also shown to be general enough to handle singularities from quasisteady control-system models. The method is shown to be more general and robust than the similar 2GCHAS method for analyzing rotorcraft dynamics.

Panda, Brahmananda↗

NASA NTRS1993

On solving the irregularly-shaped plate buckling problem

This paper develops methods for finding critical buckling loads of isotropic plates. In particular, this paper is motivated by the lack of tabular data for non-rectangular plates such as general trapezoids. High order displacement functions are used, and boundary conditions are imposed by utilizing the Lagrange multiplier method, which allows for a variety of constraint conditions. The element is derived with the use of minimum potential energy and subparametric mappings. An example is given, and several test cases are presented to show the validity of the approach in the one element case.

Smith, James P.↗

NASA NTRS1993

Shape optimization governed by the Euler equations using an adjoint method

A numerical approach for the treatment of optimal shape problems governed by the Euler equations is discussed. Focus is on flows with embedded shocks. A very simple problem is considered: the design of a quasi-one-dimensional Laval nozzle. A cost function and a set of Lagrange multipliers are introduced to achieve the minimum. The nature of the resulting costate equations is discussed. A theoretical difficulty that arises for cases with embedded shocks is pointed out and solved. Finally, some results are given to illustrate the effectiveness of the method.

Iollo, Angelo↗

NASA NTRS1993

Merits and limitations of optimality criteria method for structural optimization

The merits and limitations of the optimality criteria (OC) method for the minimum weight design of structures subjected to multiple load conditions under stress, displacement, and frequency constraints were investigated by examining several numerical examples. The examples were solved utilizing the Optimality Criteria Design Code that was developed for this purpose at NASA Lewis Research Center. This OC code incorporates OC methods available in the literature with generalizations for stress constraints, fully utilized design concepts, and hybrid methods that combine both techniques. Salient features of the code include multiple choices for Lagrange multiplier and design variable update methods, design strategies for several constraint types, variable linking, displacement and integrated force method analyzers, and analytical and numerical sensitivities. The performance of the OC method, on the basis of the examples solved, was found to be satisfactory for problems with few active constraints or with small numbers of design variables. For problems with large numbers of behavior constraints and design variables, the OC method appears to follow a subset of active constraints that can result in a heavier design. The computational efficiency of OC methods appears to be similar to some mathematical programming techniques.

Patnaik, Surya N.↗

NASA NTRS1994

Optimization for minimum sensitivity to uncertain parameters

A procedure to design a structure for minimum sensitivity to uncertainties in problem parameters is described. The approach is to minimize directly the sensitivity derivatives of the optimum design with respect to fixed design parameters using a nested optimization procedure. The procedure is demonstrated for the design of a bimetallic beam for minimum weight with insensitivity to uncertainties in structural properties. The beam is modeled with finite elements based on two dimensional beam analysis. A sequential quadratic programming procedure used as the optimizer supplies the Lagrange multipliers that are used to calculate the optimum sensitivity derivatives. The method was perceived to be successful from comparisons of the optimization results with parametric studies.

Pritchard, Jocelyn I.↗

NASA NTRS1994

Finite element approximation of an optimal control problem for the von Karman equations

This paper is concerned with optimal control problems for the von Karman equations with distributed controls. We first show that optimal solutions exist. We then show that Lagrange multipliers may be used to enforce the constraints and derive an optimality system from which optimal states and controls may be deduced. Finally we define finite element approximations of solutions for the optimality system and derive error estimates for the approximations.

Hou, L. Steven↗

NASA NTRS1994

Multigrid one shot methods for optimal control problems: Infinite dimensional control

The multigrid one shot method for optimal control problems, governed by elliptic systems, is introduced for the infinite dimensional control space. ln this case, the control variable is a function whose discrete representation involves_an increasing number of variables with grid refinement. The minimization algorithm uses Lagrange multipliers to calculate sensitivity gradients. A preconditioned gradient descent algorithm is accelerated by a set of coarse grids. It optimizes for different scales in the representation of the control variable on different discretization levels. An analysis which reduces the problem to the boundary is introduced. It is used to approximate the two level asymptotic convergence rate, to determine the amplitude of the minimization steps, and the choice of a high pass filter to be used when necessary. The effectiveness of the method is demonstrated on a series of test problems. The new method enables the solutions of optimal control problems at the same cost of solving the corresponding analysis problems just a few times.

Arian, Eyal↗

NASA NTRS1994

Comparison of transform coding methods with an optimal predictor for the data compression of digital elevation models

Statistical encoding techniques enable the reduction of the number of bits required to encode a set of symbols, and are derived from their probabilities. Huffman encoding is an example of statistical encoding that has been used for error-free data compression. The degree of compression given by Huffman encoding in this application can be improved by the use of prediction methods. These replace the set of elevations by a set of corrections that have a more advantageous probability distribution. In particular, the method of Lagrange Multipliers for minimization of the mean square error has been applied to local geometrical predictors. Using this technique, an 8-point predictor achieved about a 7 percent improvement over an existing simple triangular predictor.

Lewis, Michael↗

NASA NTRS1994

Compression strength of composite primary structural components

The linear elastic response is determined for an internally pressurized, long circular cylindrical shell stiffened on the inside by a regular arrangement of identical stringers and identical rings. Periodicity of this configuration permits the analysis of a portion of the shell wall centered over a generic stringer-ring joint; i.e., a unit cell model. The stiffeners are modeled as discrete beams, and the stringer is assumed to have a symmetrical cross section and the ring an asymmetrical section. Asymmetery causes out-of-plane bending and torsion of the ring. Displacements are assumed as truncated double Fourier series plus simple terms in the axial coordinate to account for the closed and pressure vessel effect (a non-periodic effect). The interacting line loads between the stiffeners and the inside shell wall are Lagrange multipliers in the formulation, and they are also assumed as truncated Fourier series. Displacement continuity constraints between the stiffeners and shell along the contact lines are satisfied point-wise. Equilibrium is imposed by the principle of virtual work. A composite material crown panel from the fuselage of a large transport aircraft is the numerical example. The distributions of the interacting line loads, and the out-of-plane bending moment and torque in the ring, are strongly dependent on modeling the deformations due to transverse shear and cross-sectional warping of the ring in torsion. This paper contains the results from the semiannual report on research on 'Pressure Pillowing of an Orthogonally Stiffened Cylindrical Shell'. The results of the new work are illustrated in the included appendix.

Johnson, Eric R.↗