Search NASASearch

Engineering topics

Watson, L. T.

Publications and source records attributed to Watson, L. T..

Tracing the efficient curve for multi-objective control-structure optimization

A recently developed active set algorithm for tracing parameterized optima is adapted to multiobjective optimization. The algorithm traces a path of Kuhn-Tucker points using homotopy curve tracking techniques, and is based on identifying and maintaining the set of active constraints. Second order necessary optimality conditions are used to determine nonoptimal stationary points on the path. In the bi-objective optimization case the algorithm is used to trace the curve of efficient solution (Pareto optima). As an example, the algorithm is applied to the simultaneous minimization of the weight and control force of a ten-bar truss with two collocated sensors and actuators, with some interesting results.

Rakowska, J.

Note on unit tangent vector computation for homotopy curve tracking on a hypercube

Probability-one homotopy methods are a class of methods for solving nonlinear systems of equations that are globally convergent from an arbitrary starting point. The essence of all such algorithms is the construction of an appropriate homotopy map and subsequent tracking of some smooth curve in the zero set of the homotopy map. Tracking a homotopy curve involves finding the unit tangent vector at different points along the zero curve, which amounts to calculating the kernel of the n x (n + 1) Jacobian matrix. While computing the tangent vector is just one part of the curve tracking algorithm, it can require a significant percentage of the total tracking time. This note presents computational results showing the performance of several different parallel orthogonal factorization/triangular system solving algorithms for the tangent vector computation on a hypercube.

Chakraborty, A.

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.

Parallel homotopy curve tracking on a hypercube

An investigation is conducted to find good parallel algorithms for solving systems of nonlinear equations using probability-one homotopy methods. Particular attention is paid to algorithms for the hypercube. Methods for one of the most computationally expensive steps of the homotopy approach, the computation of the kernel of the Jacobian matrix of the homotopy map, are studied. General nonlinear systems of equations with small and dense Jacobian matrices are considered, however, polynomial systems are not, since their structure leads to different strategies for parallelism. The mathematics behind the homotopy algorithm is summarized and the use of orthogonal factorizations is discussed. Parallel algorithms for orthogonal factorizations and triangular system solving are described. Computational results are presented and discussed.

Chakraborty, A.

Efficient identification of critical stresses in structures subject to dynamic loads

Optimum structural design problems generally employ constraints which are parametric in terms of space and time variables. A parametric constraint may be replaced by equivalent critical point constraints at its local minima for optimization applications. In complex structures, accurate identification of such critical points is computationally expensive due to the cost of finite element analyses. Three techniques are described for efficiently and accurately identifying critical points for space- and time-dependent parametric constraints. An adaptive search technique and a spline interpolation technique are developed for exactly known response. A least squares spline approximation is suggested for noisy behavior. A helicopter tail-boom structure subjected to transient loading is used as an example to demonstrate the techniques described. All three techniques are shown to be computationally efficient for critical point identification and the least squares approximation also removes noise from the data. The case of multiple constraints per element is shown to be particularly suited to the use of spline techniques.

Haftka, R. T.

Design oriented identification of critical times in transient response

Two techniques are presented for reducing the computational effort in identifying the critical time points. One approach is an adaptive search technique, well suited for the case where the response is exactly known. The other technique, useful for noisy response, is based on a least-squares spline approximation of the response. The possibility of grouping several closely spaced local peaks to identify a single super peak from each group is also investigated. The computational efficiency of the techniques proposed here is illustrated by two examples.

Haftka, R. T.

A quasi-Newton versus a homotopy method for nonlinear structural analysis

The globally convergent quasi-Newton minimization algorithm and the homotopy algorithms are discussed in detail and their effectiveness in solving certain classes of highly nonlinear problems of structural analysis is demonstrated. The application of the double dogleg strategy controls the directions and step-lengths of the quasi-Newtonian algorithm and overcomes the problem of nonpositive definite Hessians being produced during the iteration process. The algorithms are applied to a centrally loaded clamped beam, the snap-through of a shallow arch, and a shallow reticulated dome.

Kamat, M. P.

Minimization versus homotopy algorithms

The relative merits and demerits of the minimization techniques are assessed using globally convergent quasi-Newton algorithms on the one hand and the homotopy algorithms on the other hand for the solution of problems of nonlinear structural analysis. Like the homotopy algorithms, the globally convergent quasi-Newton algorithms are equally suited for the solution of the nonlinear equations of structural analysis directly without having to pose the problem as an equivalent minimization problem. In the close neighborhood of the limit and bifurcation points quasi-Newton algorithms experience difficulties. Homotopy algorithms are robust for practically all types of nonlinear problems but are computationally not as cost effective since they provide an extremely accurate prediction of the response by calculating it as a large number of points. Globally convergent algorithms can perform well with very approximate Hessians, while homotopy algorithms require extremely accurate Hessians. While quasi-Newton algorithms can be very easily structured to exploit sparsity and symmetry, homotopy algorithms are not presently so structured and would require special modifications for exploitation of such features without sacrificing robustness and global convergence.

Kamat, M. P.

Acoustic propagation in partially choked converging-diverging ducts

A computer model based on the wave-envelope technique is used to study acoustic propagation in converging-diverging hard walled and lined circular ducts carrying near sonic mean flows. The influences of the liner admittance, boundary layer thickness, spinning mode number, and mean Mach number are considered. The numerical results indicate that the diverging portion of the duct can have a strong reflective effect for partially choked flows.

Kelly, J. J.

Propagation of spinning acoustic modes in partially choked converging ducts

A computer model based on the wave-envelope technique is used to study the propagation of spinning acoustic modes in converging hard-walled and lined circular ducts carrying near sonic mean flows. The results show that with increasing spinning mode number the intensification of the acoustic signal at the throat decreases for upstream propagation. The influence of the throat Mach number, frequency, boundary-layer thickness, and liner admittance on the propagation of spinning modes is considered.

Nayfeh, A. H.

Non-linear propagation in near sonic flows

A nonlinear analysis is developed for sound propagation in a variable-area duct in which the mean flow approaches choking conditions. A quasi-one-dimensional model is used and the nonlinear analysis represents the acoustic disturbance as a sum of interacting harmonics. The numerical procedure is stable for cases of strong interaction and is able to integrate through the throat region without any numerical instability.

Nayfeh, A. H.