Search NASA⌕ Search

SEARCH · Search NASA

Results for “force matrix”

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 217 records · Page 12

Experimental study of unsteady hydrodynamic force matrices on whirling centrifugal pump impellers

An experimental facility was constructed and instrumented. A set of centrifugal flow pumps whose impellers were made to follow a controlled circular whirl motion were studied. The aim was to characterize the steady and unsteady fluid forces measured on the impeller under various pump operating conditions. The postulation was that the unsteady lateral forces result from interactions between the impeller and the surrounding diffuser and/or volute (via the working fluid), and that under certain flow regimes these forces can drive unstable lateral motions of the pump rotor. The lateral hydrodynamic forces were decomposed into their steady and unsteady parts, the latter being further expressed in terms of a generalized fluid stiffness matrix. Conclusions regarding the effect of impeller geometry could not be reached given the similarity of the tested designs. However, other results on phenomena such as skin friction and leakage flow are presented. Some of the findings are compared to experimental and theoretical data from other sources. Finally, the rotordynamic consequences of the results are discussed as the present data were applied to the case of the Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump (HPOTP).

Belgacem, Jery↗

Tensegrity: the architectural basis of cellular mechanotransduction

Physical forces of gravity, hemodynamic stresses, and movement play a critical role in tissue development. Yet, little is known about how cells convert these mechanical signals into a chemical response. This review attempts to place the potential molecular mediators of mechanotransduction (e.g. stretch-sensitive ion channels, signaling molecules, cytoskeleton, integrins) within the context of the structural complexity of living cells. The model presented relies on recent experimental findings, which suggests that cells use tensegrity architecture for their organization. Tensegrity predicts that cells are hard-wired to respond immediately to mechanical stresses transmitted over cell surface receptors that physically couple the cytoskeleton to extracellular matrix (e.g. integrins) or to other cells (cadherins, selectins, CAMs). Many signal transducing molecules that are activated by cell binding to growth factors and extracellular matrix associate with cytoskeletal scaffolds within focal adhesion complexes. Mechanical signals, therefore, may be integrated with other environmental signals and transduced into a biochemical response through force-dependent changes in scaffold geometry or molecular mechanics. Tensegrity also provides a mechanism to focus mechanical energy on molecular transducers and to orchestrate and tune the cellular response.

Review↗

Numerical studies of major and minor stratospheric warmings caused by orographic forcing

A primitive equation spectral model using spherical harmonics is formulated to study dynamic interactions between the troposphere and stratosphere in association with sudden stratospheric warmings. Using sigma coordinates for five tropospheric layers and log-pressure coordinates for 26 stratospheric and mesospheric layers, separate model equations for each system are combined to form single matrix governing equations. The gradual introduction of large scale topography to balanced initial states representative of observed mean winter conditions in the Northern Hemisphere is used for the generation of planetary waves during 40-day time integrations. Results of these integrations indicate that stratospheric warmings can be simulated by this orographic forcing and that mean momentum flux divergence due to zonal mean motion appears to be an essential mechanism of these simulated sudden warmings. It was found that the strength of the polar night jet can be a determining factor whether a warming becomes 'major' or 'minor'.

Koermer, J. P.↗

Design Spectrum Analysis in NASTRAN

The utility of Design Spectrum Analysis is to give a mode by mode characterization of the behavior of a design under a given loading. The theory of design spectrum is discussed after operations are explained. User instructions are taken up here in three parts: Transient Preface, Maximum Envelope Spectrum, and RMS Average Spectrum followed by a Summary Table. A single DMAP ALTER packet will provide for all parts of the design spectrum operations. The starting point for getting a modal break-down of the response to acceleration loading is the Modal Transient rigid format. After eigenvalue extraction, modal vectors need to be isolated in the full set of physical coordinates (P-sized as opposed to the D-sized vectors in RF 12). After integration for transient response the results are scanned over the solution time interval for the peak values and for the times that they occur. A module called SCAN was written to do this job, that organizes these maxima into a diagonal output matrix. The maximum amplifier in each mode is applied to the eigenvector of each mode which then reveals the maximum displacements, stresses, forces and boundary reactions that the structure will experience for a load history, mode by mode. The standard NASTRAN output processors have been modified for this task. It is required that modes be normalized to mass.

Butler, T. G.↗

Computational Simplification of Robot-Arm Dynamics

Report presents two general methods for reducing mathematical complexities of state equations of robot-arm dynamics. Robot arms contain both rotary and linear joints. Both methods start with homogenous coordinates and Lagrangian formulation of mechanics briefly summarized in report. First method uses matrix-analysis techniques; second, vector-analysis techniques. Vector-analysis method includes new differential-vector representation of centripetal and Coriolis forces.

Bejczy, A. K.↗

A Transient Response Method for Linear Coupled Substructures

A method is presented for determining the transient response of a discrete coordinate model of a linear structural system composed of substructures. The method is applicable to systems consisting of any number of substructures, both determinate and indeterminate interface boundaries, and any topological arrangement of the substructures. The method is simple to implement from a computational point of view because the equations of motion of each of the substructures are solved independently, and the interface boundary compatibility conditions are enforced at each integration time step by a matrix multiplication. The method is demonstrated for a structural system consisting of two beam segments and acted upon by a time dependent force. The numerical results from the demonstration problem validates the accuracy of the method. The application of this method to structural systems with changing interface boundary conditions between substructures is discussed.

Admire, J. R.↗

Static aeroelastic analysis of a three-dimensional generic wing

A continuation of research on the static aeroelastic analysis of a generic wing configuration is presented. Results of the study of the asymmetric oblique wing model developed by Rockwell International, in conjunction with the NASA Oblique Wing Research Aircraft Program, are reported. The capability to perform static aeroelastic analyses of an oblique wing at arbitrary skew positions is demonstrated by applying the MSC/NASTRAN static analysis scheme modified by the aerodynamic influence coefficient matrix created by the NASA Ames aerodynamic panel codes. The oblique wing is studied at two skew angles, and in particular, the capability to calculate 3-D thickness effects on the aerodynamic properties of the wing is investigated. The ability to model asymmetric wings in both subsonic and supersonic Mach numbers is shown. The aerodynamic influence coefficient matrix computed by the external programs is inserted in MSC/NASTRAN static aeroelasticity analysis run stream to compute the aeroelastic deformation and internal forces. Various aerodynamic coefficients of the oblique wing were computed for two Mach numbers, 0.7 and 1.4, and the angle of attach -5 through 15 deg.

Green, John A.↗

Random field estimation approach to robot dynamics

The difference equations of Kalman filtering and smoothing recursively factor and invert the covariance of the output of a linear state-space system driven by a white-noise process. Here it is shown that similar recursive techniques factor and invert the inertia matrix of a multibody robot system. The random field models are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. They are easier to describe than the models based on classical mechanics, which typically require extensive derivation and manipulation of equations of motion for complex mechanical systems. With the spatially random models, more primitive locally specified computations result in a global collective system behavior equivalent to that obtained with deterministic models. The primary goal of applying random field estimation is to provide a concise analytical foundation for solving robot control and motion planning problems.

Rodriguez, Guillermo↗

Thermoviscoplastic analysis of fibrous periodic composites using triangular subvolumes

The nonlinear viscoplastic behavior of fibrous periodic composites is analyzed by discretizing the unit cell into triangular subvolumes. A set of these subvolumes can be configured by the analyst to construct a representation for the unit cell of a periodic composite. In each step of the loading history, the total strain increment at any point is governed by an integral equation which applies to the entire composite. A Fourier series approximation allows the incremental stresses and strains to be determined within a unit cell of the periodic lattice. The nonlinearity arising from the viscoplastic behavior of the constituent materials comprising the composite is treated as fictitious body force in the governing integral equation. Specific numerical examples showing the stress distributions in the unit cell of a fibrous tungsten/copper metal matrix composite under viscoplastic loading conditions are given. The stress distribution resulting in the unit cell when the composite material is subjected to an overall transverse stress loading history perpendicular to the fibers is found to be highly heterogeneous, and typical homogenization techniques based on treating the stress and strain distributions within the constituent phases as homogeneous result in large errors under inelastic loading conditions.

Walker, Kevin P.↗

Adjoints and Low-rank Covariance Representation

Quantitative measures of the uncertainty of Earth System estimates can be as important as the estimates themselves. Second moments of estimation errors are described by the covariance matrix, whose direct calculation is impractical when the number of degrees of freedom of the system state is large. Ensemble and reduced-state approaches to prediction and data assimilation replace full estimation error covariance matrices by low-rank approximations. The appropriateness of such approximations depends on the spectrum of the full error covariance matrix, whose calculation is also often impractical. Here we examine the situation where the error covariance is a linear transformation of a forcing error covariance. We use operator norms and adjoints to relate the appropriateness of low-rank representations to the conditioning of this transformation. The analysis is used to investigate low-rank representations of the steady-state response to random forcing of an idealized discrete-time dynamical system.

Tippett, Michael K.↗

Effects of Solute Concentrations on Kinetic Pathways in Ni-Al-Cr Alloys

The kinetic pathways resulting from the formation of coherent gamma'-precipitates from the gamma-matrix are studied for two Ni-Al-Cr alloys with similar gamma'-precipitate volume fractions at 873 K. The details of the phase decompositions of Ni-7.5Al-8.5Cr at.% and Ni-5.2Al-14.2Cr at.% for aging times from 1/6 to 1024 h are investigated by atom-probe tomography, and are found to differ significantly from a mean-field description of coarsening. The morphologies of the gamma'-precipitates of the alloys are similar, though the degrees of gamma'-precipitate coagulation and coalescence differ. Quantification within the framework of classical nucleation theory reveals that differences in the chemical driving forces for phase decomposition result in differences in the nucleation behavior of the two alloys. The temporal evolution of the gamma'-precipitate average radii and the gamma-matrix supersaturations follow the predictions of classical coarsening models. The compositional trajectories of the gamma-matrix phases of the alloys are found to follow approximately the equilibrium tie-lines, while the trajectories of the gamma'-precipitates do not, resulting in significant differences in the partitioning ratios of the solute elements.

Booth-Morrison, Christopher↗

Working Toward Nanotube Composites

One of the most attractive applications of single-wall carbon nanotubes (SWNT) is found in the area of structural materials. Nanotubes have a unique combination of high strength, modulus, and elongation to failure, and therefore have potential to significantly enhance the mechanical properties of today's composites. This is especially attractive for the aerospace industry looking for any chance to save weight. This is why NASA has chosen to tackle this difficult application of SWNT. Nanotube properties differ significantly from that of conventional carbon fibers, and a whole new set of problems, including adhesion and dispersion in the adhesive polymer matrix, must be resolved in order to engineer superior composite materials. From recent work on a variety of applications it is obvious that the wide range of research in nanotubes will lead to advances in physics, chemistry, and engineering. However, the possibility of ultralightweight structures is what causes dreamers to really get excited. One of the important issues in composite engineering is aspect ratio of the fibers, since it affects load transfer in composites. Nanotube length was a gray area for years, since they are formed in bundles, making it impossible to monitor individual nanotube length. Even though bundles are observed to be tens and hundreds of microns long, they can be built of relatively short tubes weakly bound by Van der Waals forces. Nanotube length can be affected by subsequent purification and ultrasound processing, which has been necessary in order to disperse nanotubes and introduce them into a polymer matrix. Some calculations show that nanotubes with 10(exp 5) aspect ratio may be necessary to achieve good load transfer. We show here that nanotubes produced in our laser system are as much as tens of microns long and get cut into lengths of hundreds of nanometers during ultrasound processing. Nanotube length was measured by AFM on pristine nanotube specimens as well, as after sonication. In each case great care was taken to measure individual nanotubes, rather than bundles. Pristine nanotubes were collected on quartz substrates placed directly in the laser oven and exposed for 0.5 s. This results in an equal mix of bundles and individual nanotubes. Nanotube length measurements were limited practically by the lateral span of the AFM scanner, but nanotube length is certainly in excess of 20 micrometers.

Arepalli, Sivaram↗

A new look at the simultaneous analysis and design of structures

The minimum weight optimization of structural systems, subject to strength and displacement constraints as well as size side constraints, was investigated by the Simultaneous ANalysis and Design (SAND) approach. As an optimizer, the code NPSOL was used which is based on a sequential quadratic programming (SQP) algorithm. The structures were modeled by the finite element method. The finite element related input to NPSOL was automatically generated from the input decks of such standard FEM/optimization codes as NASTRAN or ASTROS, with the stiffness matrices, at present, extracted from the FEM code ANALYZE. In order to avoid ill-conditioned matrices that can be encountered when the global stiffness equations are used as additional nonlinear equality constraints in the SAND approach (with the displacements as additional variables), the matrix displacement method was applied. In this approach, the element stiffness equations are used as constraints instead of the global stiffness equations, in conjunction with the nodal force equilibrium equations. This approach adds the element forces as variables to the system. Since, for complex structures and the associated large and very sparce matrices, the execution times of the optimization code became excessive due to the large number of required constraint gradient evaluations, the Kreisselmeier-Steinhauser function approach was used to decrease the computational effort by reducing the nonlinear equality constraint system to essentially a single combined constraint equation. As the linear equality and inequality constraints require much less computational effort to evaluate, they were kept in their previous form to limit the complexity of the KS function evaluation. To date, the standard three-bar, ten-bar, and 72-bar trusses have been tested. For the standard SAND approach, correct results were obtained for all three trusses although convergence became slower for the 72-bar truss. When the matrix displacement method was used, correct results were still obtained, but the execution times became excessive due to the large number of constraint gradient evaluations required. Using the KS function, the computational effort dropped, but the optimization seemed to become less robust. The investigation of this phenomenon is continuing. As an alternate approach, the code MINOS for the optimization of sparse matrices can be applied to the problem in lieu of the Kreisselmeier-Steinhauser function. This investigation is underway.

Striz, Alfred G.↗

The Turning Point of the Aerosol Era

Over the CMIP6 historical period (1850–2014), aerosols provided the largest negative forcing compared to all other climate forcings via their ability to absorb or scatter solar radiation and alter clouds. Aerosols played an important role in counterbalancing warming by greenhouse gases (GHGs). Here we study aerosol forcing trends in the CMIP6 simulations of the NASA Goddard Institute for Space Studies (GISS) ocean-atmosphere ModelE version 2.1 (GISS-E2.1-G) using a fully coupled atmospheric composition configuration, including interactive gas-phase chemistry, and either an aerosol microphysical (MATRIX) or a mass-based aerosol (OMA) module. Simulations of the CMIP6 historical period are analyzed as well as four Shared Socioeconomic Pathway (SSP) future scenarios for 2015–2100: SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5. The main conclusion of this study is that aerosol forcing in the GISS model has reached its turning point, switching from globally increasing to a decreasing trend in the first decade of the 21st century. This result is robust, independent of which aerosol module or SSP scenario is used. Non-linear aerosol-cloud interactions dominate as a forcing agent over aerosol-radiation interactions. Aerosols' ability to counterbalance GHG forcing on the global scale is today at a level comparable to that at the beginning of the last century. In the 1980s, the decade of largest global aerosol loads, aerosols balanced up to 80% of GHG forcing. As a consequence, global warming of the last decades, which is primarily driven by greenhouse gases, has been augmented by the effect of decreasing aerosol cooling in our model. By the end of this century, following the SSP scenarios, aerosols will only counterbalance 0%–20% of GHG forcing, depending on model and on scenario.

Aerosol forcing↗

Dynamics of a spinning space station with a counterweight connected by multiple cables.

A unique and relatively simple approach is presented for obtaining the linearized equations of motion. A conservative system consisting of two rigid bodies connected by any number of massless cables with linear axial stiffness is assumed. The cable forces and torques are expanded in a Taylor series about the equilibrium values of the system coordinates which results in a cable stiffness matrix. The method of obtaining the equilibrium values of the coordinates is discussed and results are presented. The range of validity of the linear model is determined by comparing results with a digital simulation of the nonlinear system.

Nixon, D. D.↗

Flapping response characteristics of hingeless rotor blades by a gereralized harmonic balance method

Linearized equations of motion for the flapping response of flexible rotor blades in forward flight are derived in terms of generalized coordinates. The equations are solved using a matrix form of the method of linear harmonic balance, yielding response derivatives for each harmonic of the blade deformations and of the hub forces and moments. Numerical results and approximate closed-form expressions for rotor derivatives are used to illustrate the relationships between rotor parameters, modeling assumptions, and rotor response characteristics. Finally, basic hingeless rotor response derivatives are presented in tabular and graphical form for a wide range of configuration parameters and operating conditions.

Peters, D. A.↗

The role of rotor impedance in the vibration analysis of rotorcraft, part 4

A method for a strongly idealized case of vertical excitation and for rolling and pitching moment excitation of a four bladed hingeless rotor on an up-focussing flexible mount is developed. The aeroelastic rotor impedances are computed directly with a finite blade element method that includes aerodynamics. The rotor impedance matrix for three or more blades is determined from the root moment impedance for a single blade by a simple multiblade transformation rule. Force and moment amplitudes transferred from the rotor to support are found to be critically dependent on the support dynamics.

Hohenemser, K. H.↗

The role of rotor impedance in the vibration analysis of rotorcraft

In an improved method which retains the advantage of separate treatment of rotor and airframe, the rotor impedance is used to correct the input to the airframe. This improved method is illustrated for a strongly idealized case of vertical excitation and then for rolling and pitching moment excitation of a four bladed hingeless rotor on an up-focussing flexible mount. Contrary to the usual approach that represents aeroelastic blade motions by a series of normal blade modes in vacuum, the aeroelastic rotor impedances are computed directly with a finite blade element method that includes aerodynamics. The rotor impedance matrix for three or more blades is determined from the root moment impedance for a single blade by a simple multiblade transformation rule. Force and moment amplitudes transferred from the rotor to the support are found to be critically dependent on the support dynamics.

Hohenemser, K. H.↗