Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 559 records · Page 31

A frequency domain analysis for damped space structures

A method is presented for the analysis of damped structural systems in which the structural components are represented by impedance models and analyzed in the frequency domain. Methods are presented to assemble and condense system impedance matrices, and then to identify approximate mass, stiffness, and damping matrices for systems whose impedances are complicated functions of frequency. Formulas are derived for determination of approximate values for system natural frequencies and damping using frequency domain quantities. The sensitivities of these approximate values to system parameter changes are analyzed. The implementation of these analysis tools is discussed and applied to a simple mechanical system.

Hagood, Nesbitt W.↗

Investigation of scene identification algorithms for radiation budget measurements

The computation of earth radiation budget from satellite measurements requires the identification of the scene in order to select spectral factors and bidirectional models. A scene identification procedure is developed for AVHRR SW and LW data by using two radiative transfer models. These AVHRR GAC pixels are then attached to corresponding ERBE pixels and the results are sorted into scene identification probability matrices. These scene intercomparisons show that there generally is a higher tendency for underestimation of cloudiness over ocean at high cloud amounts, e.g., mostly cloudy instead of overcast, partly cloudy instead of mostly cloudy, for the ERBE relative to the AVHRR results. Reasons for this are explained. Preliminary estimates of the errors of exitances due to scene misidentification demonstrates the high dependency on the probability matrices. While the longwave error can generally be neglected the shortwave deviations have reached maximum values of more than 12 percent of the respective exitances.

Diekmann, F. J.↗

Single-mode projection filters for modal parameter identification for flexible structures

Single-mode projection filters are developed for eigensystem parameter identification from both analytical results and test data. Explicit formulations of these projection filters are derived using the orthogonal matrices of the controllability and observability matrices in the general sense. A global minimum optimization algorithm is applied to update the filter parameters by using the interval analysis method. The updated modal parameters represent the characteristics of the test data. For illustration of this new approach, a numerical simulation for the MAST beam structure is shown by using a one-dimensional global optimization algorithm to identify modal frequencies and damping. The projection filters are practical for parallel processing implementation.

Huang, Jen-Kuang↗

A petrographic thin sectioning technique for evaluating composite materials

Petrographic thin sectioning by a low-speed diamond saw has been used in conjunction with transmission polarized light microscopy for the characterization of the microstructure and deformation mechanisms of a variety of polymer systems. It has proven possible by these means to study three types of thermoplastic matrices for composite applications: PEEK, BPA-based polycarbonate (PC), and a rubber-modified PC. The reinforcing fibers for these matrices were in all cases AS4 carbon fibers, unidirectionally arrayed. Superior analyzability of matrix morphology and subsurface fracture processes is achieved by thin sectioning.

Parker, D. S.↗

Factors influencing the Mode I interlaminar fracture toughness of a rubber toughened thermoplastic matrix composite

The use of a rubber modified thermoplastic resin has been investigated as a method to improve the Mode I interlaminar fracture toughness of a unidirectional continuous carbon fiber composite. Test results show that the improvement in the fracture toughness is less than expected due to rubber particle agglomeration, solvent and molding induced crystallization of the matrix and poor fiber/matrix adhesion. The plastic zone in composites utilizing tough matrices can extend well beyond a single interfibrillar spacing. However, the development of the plastic zone is limited due to the failure of the fiber/matrix interface. In order to fully evaluate the potential of tough composites using toughened matrices, any improvement made in the matrix toughness must be coupled with improvements in the fiber/matrix adhesion.

Parker, D. S.↗

Modal techniques for analyzing airplane dynamics

A series of techniques are presented for analyzing airplane control and output characteristics. State-space matrix equations describing the linear perturbation dynamics are transformed from physical coordinates into scaled, modal coordinates. Techniques for analyzing the impacts of system inputs on the fundamental modes of motion and the appearance of these modes in the system outputs are explained. Scaled modal matrices are used to compute steady-state control inputs which optimize the steady-state response of selected system outputs. Graphics which promote quick understanding of the analysis are presented to display the resulting vectors and matrices. The defined analysis techniques are applied to an example airplane model, illustrating the insight which can be acquired using the described modal techniques.

Arbuckle, P. Douglas↗

Coordination of two redundant robots

A solution is presented for the optimal force distributions of two redundant coordinating robots. Two robots whose end-effectors firmly grasp a common inertial load form a closed-chain structure, and thus the robots must observe a set of constraint conditions on (1) the end-effector positions, orientations, linear and angular velocities and (2) the forces and torques inserted onto the end-effectors by the load. To achieve the maximum load capability, the optimal force distribution under the known posture is developed. If the robots are redundant, joint solutions may be obtained if some optimization criterion such as manipulability is imposed. The singularity, manipulability, and force ellipsoids are investigated in terms of the two manipulator Jacobian matrices. The manipulability is shown to correspond to the cascade structure, while the force ellipsoid corresponds to the parallel structure in terms of the manipulator Jacobian matrices.

Tao, Jian M.↗

Robust adaptive regulation without persistent excitation

A globally convergent adaptive regulator for minimum- or nonminimum-phase systems subject to bounded disturbances and unmodeled dynamics is presented. The control strategy is designed for a particular input-output representation obtained from the state space representation of the system. The leading coefficient of the representation is the product of the observability and controllability matrices of the system. The controller scheme uses a Least-Squares identification algorithm a with dead zone. The dead zone is chosen to obtain convergence properties on the estimates and on the covariance matrix as well. This allows the definition of modified estimates which secure well-conditioned matrices in the adaptive control law. Explicit bounds on the plant output are given.

Lozano-Leal, Rogelio↗

Dynamic substructuring by the boundary flexibility vector method of component mode synthesis

Component mode synthesis (CMS) is a method of dynamic analysis, for structures having a large number of degrees of freedom (DOF). These structures often required lengthy computer CPU time and large computer memory resources, if solved directly by the finite-element method (FEM). In CMS, the structure is divided into independent components in which the DOF are defined by a set of generalized coordinates defined by displacement shapes. The number of the generalized coordinates are much less than the original number of physical DOF, in the component. The displacement shapes are used to transform the component property matrices and any applied external loads, to a reduced system of coordinates. Reduced system property matrices are assembled, and any type of dynamic analysis is carried out in the reduced coordinate system. Any obtained results are back transformed to the original component coordinate systems. In all conventional methods of CMS, the mode shapes used for components are dynamic mode shapes, supplemented by static deflected shapes. Historically, all the dynamic mode shapes used in conventional CMS are the natural modes (eigenvectors) of components. A new method of CMS, namely the boundary flexibility vector method of CMS, is presented. The method provides for the incorporation of a set of static Ritz vectors, referred to as boundary flexibility vectors, as a replacement and/or supplement to conventional eigenvectors, as displacement shapes for components. The generation of these vectors does not require the solution of a costly eigenvalue problem, as in the case of natural modes in conventional CMS, and hence a substantial saving in CPU time can be achieved. The boundary flexibility vectors are generated from flexibility (or stiffness) properties of components. The formulation presented is for both free and fixed-interface components, and for both the free and forced vibration problems. Free and forced vibration numerical examples are presented to verify the accuracy of the method and the saving in CPU time. Compared to conventional methods of CMS, the results indicate that by using the new method, more accurate results can be obtained with a substantial saving in CPU time.

Abdallah, Ayman Ahmed↗

An order (n) algorithm for the dynamics simulation of robotic systems

The formulation of an Order (n) algorithm for DISCOS (Dynamics Interaction Simulation of Controls and Structures), which is an industry-standard software package for simulation and analysis of flexible multibody systems is presented. For systems involving many bodies, the new Order (n) version of DISCOS is much faster than the current version. Results of the experimental validation of the dynamics software are also presented. The experiment is carried out on a seven-joint robot arm at NASA's Goddard Space Flight Center. The algorithm used in the current version of DISCOS requires the inverse of a matrix whose dimension is equal to the number of constraints in the system. Generally, the number of constraints in a system is roughly proportional to the number of bodies in the system, and matrix inversion requires O(p exp 3) operations, where p is the dimension of the matrix. The current version of DISCOS is therefore considered an Order (n exp 3) algorithm. In contrast, the Order (n) algorithm requires inversion of matrices which are small, and the number of matrices to be inverted increases only linearly with the number of bodies. The newly-developed Order (n) DISCOS is currently capable of handling chain and tree topologies as well as multiple closed loops. Continuing development will extend the capability of the software to deal with typical robotics applications such as put-and-place, multi-arm hand-off and surface sliding.

Chun, H. M.↗

Development and testing of methodology for evaluating the performance of multi-input/multi-output digital control systems

A Controller Performance Evaluation (CPE) methodology for multi-input/multi-output digital control systems was developed and tested on an aeroelastic wind-tunnel model. Modern signal processing methods were used to implement control laws and to acquire time domain data of the whole system (controller and plant) from which appropriate transfer matrices of the control system could be generated. Matrix computational procedures were used to calculate singular values of return-difference matrices at the plant input and output points to evaluate the performance of the control system. The CPE procedures effectively identified potentially destabilizing controllers and confirmed the satisfactory performance of stabilizing ones.

Pototzky, Anthony S.↗

Sparse Gaussian elimination with controlled fill-in on a shared memory multiprocessor

It is shown that in sparse matrices arising from electronic circuits, it is possible to do computations on many diagonal elements simultaneously. A technique for obtaining an ordered compatible set directly from the ordered incompatible table is given. The ordering is based on the Markowitz number of the pivot candidates. This technique generates a set of compatible pivots with the property of generating few fills. A novel heuristic algorithm is presented that combines the idea of an order-compatible set with a limited binary tree search to generate several sets of compatible pivots in linear time. An elimination set for reducing the matrix is generated and selected on the basis of a minimum Markowitz sum number. The parallel pivoting technique presented is a stepwise algorithm and can be applied to any submatrix of the original matrix. Thus, it is not a preordering of the sparse matrix and is applied dynamically as the decomposition proceeds. Parameters are suggested to obtain a balance between parallelism and fill-ins. Results of applying the proposed algorithms on several large application matrices using the HEP multiprocessor (Kowalik, 1985) are presented and analyzed.

Alaghband, Gita↗

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

Park, K. C.↗

Development and testing of methodology for evaluating the performance of multi-input/multi-output digital control systems

The development of a controller performance evaluation (CPE) methodology for multiinput/multioutput digital control systems is described. The equations used to obtain the open-loop plant, controller transfer matrices, and return-difference matrices are given. Results of applying the CPE methodology to evaluate MIMO digital flutter suppression systems being tested on an active flexible wing wind-tunnel model are presented to demonstrate the CPE capability.

Polotzky, Anthony S.↗

A linear quadratic regulator approach to the stabilization of uncertain linear systems

This paper presents a linear quadratic regulator approach to the stabilization of uncertain linear systems. The uncertain systems under consideration are described by state equations with the presence of time-varying unknown-but-bounded uncertainty matrices. The method is based on linear quadratic regulator (LQR) theory and Liapunov stability theory. The robust stabilizing control law for a given uncertain system can be easily constructed from the symmetric positive-definite solution of the associated augmented Riccati equation. The proposed approach can be applied to matched and/or mismatched systems with uncertainty matrices in which only their matrix norms are bounded by some prescribed values and/or their entries are bounded by some prescribed constraint sets. Several numerical examples are presented to illustrate the results.

Shieh, L. S.↗

Ada Namelist Package

Ada Namelist Package, developed for Ada programming language, enables calling program to read and write FORTRAN-style namelist files. Features are: handling of any combination of types defined by user; ability to read vectors, matrices, and slices of vectors and matrices; handling of mismatches between variables in namelist file and those in programmed list of namelist variables; and ability to avoid searching entire input file for each variable. Principle benefits derived by user: ability to read and write namelist-readable files, ability to detect most file errors in initialization phase, and organization keeping number of instantiated units to few packages rather than to many subprograms.

Klumpp, Allan R.↗

Polynomial Transformations For Discrete-Time Linear Systems

Transformations based on polynomial matrices of finite degree developed for use in computing functions for compensation, inversion, and approximation of discrete-time, multivariable, linear systems. Method derived from z-transform transfer-function form of matrices. Applicable to cascade-compensation problems in design of control systems.

Baram, Yoram↗

Physically weighted approximations of unsteady aerodynamic forces using the minimum-state method

The Minimum-State Method for rational approximation of unsteady aerodynamic force coefficient matrices, modified to allow physical weighting of the tabulated aerodynamic data, is presented. The approximation formula and the associated time-domain, state-space, open-loop equations of motion are given, and the numerical procedure for calculating the approximation matrices, with weighted data and with various equality constraints are described. Two data weighting options are presented. The first weighting is for normalizing the aerodynamic data to maximum unit value of each aerodynamic coefficient. The second weighting is one in which each tabulated coefficient, at each reduced frequency value, is weighted according to the effect of an incremental error of this coefficient on aeroelastic characteristics of the system. This weighting yields a better fit of the more important terms, at the expense of less important ones. The resulting approximate yields a relatively low number of aerodynamic lag states in the subsequent state-space model. The formulation forms the basis of the MIST computer program which is written in FORTRAN for use on the MicroVAX computer and interfaces with NASA's Interaction of Structures, Aerodynamics and Controls (ISAC) computer program. The program structure, capabilities and interfaces are outlined in the appendices, and a numerical example which utilizes Rockwell's Active Flexible Wing (AFW) model is given and discussed.

Karpel, Mordechay↗