NON-LINEAR SHALLOW SHELL ANALYSIS BY THE MATRIX FORCE METHOD
Nonlinear shallow shell analysis by the matrix force method
SEARCH · Search NASA
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.
Nonlinear shallow shell analysis by the matrix force method
Microscopic coupling between the electron spin and the lattice vibration is responsible for an array of exotic properties from morphic effects in simple non magnets to magnetodielectric coupling in multiferroic spinels and hematites. Traditionally, a single spin–phonon coupling constant is used to characterize how effectively the lattice can affect the spin, but it is hardly enough to capture novel electromagnetic behaviors to the full extent. Here, we introduce a concept of spin–phonon dispersion to project the spin moment change along the phonon crystal momentum direction, so the entire spin change can be mapped out. Different from the phonon dispersion, the spin–phonon dispersion has both positive and negative frequency branches even in the equilibrium ground state, which correspond to the spin enhancement and spin reduction, respectively. Our study of bcc Fe and hcp Co reveals that the spin force matrix, that is, the second-order spatial derivative of spin moment, is similar to the vibrational force matrix, but its diagonal elements are smaller than the off-diagonal ones. This leads to the distinctive spin–phonon dispersion. Furthermore, the concept of spin–phonon dispersion expands the traditional Elliott-Yafet theory in nonmagnetic materials to the entire Brillouin zone in magnetic materials, thus opening the door to excited states in systems such as CoF 2 and NiO, where a strong spin-lattice coupling is detected in the THz regime.
Geometric nonlinearities for large amplitude free and forced vibrations of circular plates are investigated. Inplane displacement and inertia are included in the formulation. Finite element method is used. Harmonic force matrix for nonlinear forced vibration analysis is introduced and derived. Various out-of-plane and inplane boundary conditions are considered. The relations of amplitude - frequency ratio for different boundary conditions and various loads conditions are presented.
Multiple-mode nonlinear free and forced vibration of a beam is analyzed by the finite element method. The geometric nonlinearity is investigated. Inplane displacement and inertia (IDI) are also considered in the formulation. Harmonic force matrix is derived and explained. Nonlinear free vibration can be simply treated as a special case of the general forced vibration by setting the harmonic force matrix equal to zero. The effect of the higher modes is more pronouced for the clamped supported beam than the simply supported one. Beams without IDI yield more effect of the higher modes than the one with IDI. The effects of IDI are to reduce nonlinearity. For beams with end supports restrained from axial movement (immovable cases), only the hardening type nonlinearity is observed. However, beams of small slenderness ratio (L/R = 20) with movable end supports, the softening type nonlinearity is found. The concentrated force case yields a more severe response than the uniformly distributed force case. Finite element results are in good agreement with the solution of simple elliptic response, harmonic balance method, and Runge-Kutte method and experiment.
Effects of large deflection geometric nonlinearity and multiple modes for free and forced vibrations of beams are investigated. Inplane displacement and inplane inertia are included in the formulation. The finite element method is employed. A harmonic force matrix is derived for forced vibration analysis. The relations of amplitude, frequency ratio and forcing intensity for beams of different boundary conditions and various load conditions are presented.
Geometric non-linearities for large amplitude free and forced vibrations of circular plates are investigated. In-plane displacement and in-plane inertia are included in the formulation. The finite element method is used. An harmonic force matrix for non-linear forced vibration analysis is introduced and derived. Various out-of-plane and in-plane boundary conditions are considered. The relations of amplitude and frequency ratio for different boundary conditions and various load conditions are presented.
Techniques for defining a finite element model (FEM) for analysis of nonlinear vibrations in beam structures subjected to harmonic excitation are presented. The resulting model covers longitudinal deformation and inertial effects. The nonlinear oscillations of a beam element under forced excitation are modeled by a harmonic force matrix based on first order approximations of the Jacobian elliptic forcing function. Harmonic force and nonlinear stiffness matrices are derived and the nonlinear forced responses of beams are calculated under various boundary conditions. The results of FEM computations for simply-supported and clamped beams show that midplane stretching caused by large deflections increases the nonlinearity. Axially-restrained beams experience only hardening nonlinearity, while axially-free beams have reduced nonlinearity in deformation and inertia and an increase in linearity due to large deflection.
The finite element method has been extended to determine the response of large amplitude forced vibrations of thin plates. A harmonic force matrix of a rectangular element under uniform harmonic excitation is developed for nonlinear forced vibration analysis. Inplane deformation and inertia are both considered in the formulation. Results obtained are compared with simple elliptic response, perturbation and other approximation solutions.
A partition procedure for forced-convection conduction transient problems is presented. Mixed time partitions are defined wherein coupled conduction force-matrix equations are discretized using an implicit integration method, followed by derivation of a mixed time integration technique. Explicit-implicit and explicit-explicit partitions are performed for a stability analysis for transient conditions, e.g., those found in an actively air-cooled engine and airframe structure.
A close link is established between open loop optimal control theory and optimal design by noting certain similarities in the gradient calculations. The resulting benefits include a unified approach, together with physical insights in design sensitivity analysis, and an efficient approach for simultaneous optimal control and design. Both matrix displacement and matrix force methods are considered, and results are presented for dynamic systems, structures, and elasticity problems.
The RWDMES is a tool for modeling the disturbances imparted on spacecraft by spinning reaction wheels. Reaction wheels are usually the largest disturbance source on a precision pointing spacecraft, and can be the dominating source of pointing error. Accurate knowledge of the disturbance environment is critical to accurate prediction of the pointing performance. In the past, it has been difficult to extract an accurate wheel disturbance model since the forcing mechanisms are difficult to model physically, and the forcing amplitudes are filtered by the dynamics of the reaction wheel. RWDMES captures the wheel-induced disturbances using a hybrid physical/empirical model that is extracted directly from measured forcing data. The empirical models capture the tonal forces that occur at harmonics of the spin rate, and the broadband forces that arise from random effects. The empirical forcing functions are filtered by a physical model of the wheel structure that includes spin-rate-dependent moments (gyroscopic terms). The resulting hybrid model creates a highly accurate prediction of wheel-induced forces. It accounts for variation in disturbance frequency, as well as the shifts in structural amplification by the whirl modes, as the spin rate changes. This software provides a point-and-click environment for producing accurate models with minimal user effort. Where conventional approaches may take weeks to produce a model of variable quality, RWDMES can create a demonstrably high accuracy model in two hours. The software consists of a graphical user interface (GUI) that enables the user to specify all analysis parameters, to evaluate analysis results and to iteratively refine the model. Underlying algorithms automatically extract disturbance harmonics, initialize and tune harmonic models, and initialize and tune broadband noise models. The component steps are described in the RWDMES user s guide and include: converting time domain data to waterfall PSDs (power spectral densities); converting PSDs to order analysis data; extracting harmonics; initializing and simultaneously tuning a harmonic model and a wheel structural model; initializing and tuning a broadband model; and verifying the harmonic/broadband/structural model against the measurement data. Functional operation is through a MATLAB GUI that loads test data, performs the various analyses, plots evaluation data for assessment and refinement of analysis parameters, and exports the data to documentation or downstream analysis code. The harmonic models are defined as specified functions of frequency, typically speed-squared. The reaction wheel structural model is realized as mass, damping, and stiffness matrices (typically from a finite element analysis package) with the addition of a gyroscopic forcing matrix. The broadband noise model is realized as a set of speed-dependent filters. The tuning of the combined model is performed using nonlinear least squares techniques. RWDMES is implemented as a MATLAB toolbox comprising the Fit Manager for performing the model extraction, Data Manager for managing input data and output models, the Gyro Manager for modifying wheel structural models, and the Harmonic Editor for evaluating and tuning harmonic models. This software was validated using data from Goodrich E wheels, and from GSFC Lunar Reconnaissance Orbiter (LRO) wheels. The validation testing proved that RWDMES has the capability to extract accurate disturbance models from flight reaction wheels with minimal user effort.
The force exerted on a targeting sequence by the electrical potential across the inner mitochondrial membrane is calculated on the basis of continuum electrostatics. The force is found to vary from 3.0 pN to 2.2 pN (per unit elementary charge) as the radius of the inner membrane pore (assumed aqueous) is varied from 6.5 to 12 A, its measured range. In the present model, the decrease in force with increasing pore width arises from the shielding effect of water. Since the pore is not very much wider than the distance between water molecules, the full shielding effect of water may not be present; the extreme case of a purely membranous pore without water gives a force of 3.2 pN per unit charge, which should represent an upper limit. When applied to mitochondrial import experiments on the protein barnase, these results imply that forces between 11 +/- 2 pN and 13.5 +/- 2.5 pN catalyze the unfolding of barnase in those experiments. A comparison of these results with unfolding forces measured using atomic force microscopy is made.
Thrust vector control (TVC) systems for rocket engine propulsion traditionally use a simple linear relationship to convert between actuator forces and torques about the engine gimbal’s center-of-rotation (COR). As shown in Equation (1), the torque about the gimbal COR is proportional to the applied actuator force and the TVC moment arm (MA)—the perpendicular distance between the TVC actuator’s line-of-action (LOA) and the engine gimbal’s COR. While this fundamental relationship remains valid and accurate in a two-dimensional (2D), one-degree-of-freedom (1-DOF) context—particularly in its non-linear formulation as described in the ER63 Technical Bulletin TB-02 (Derivation of Thrust Vector Control (TVC) Engine-Gimbal / Actuator Moment-Arm Geometry)—it becomes limited when extended to three-dimensional (3D), two-degree-of-freedom (2-DOF) analyses. In 3D, out-of-plane angular displacements can arise, causing the gimbaled engine plane-of-motion to become non-coplanar with the respective actuator plane-of-motion. Such misalignments occur due to the rod-end (RE) and/or tail-stock (TE) mounting geometry of individual TVC actuators. These geometric complexities lead to inaccuracies in calculating engine torques and corresponding actuator forces if only the traditional TVC MA relationship is employed. For higher gimbal angular displacements (e.g., greater than 5 degrees), these inaccuracies become more pronounced, necessitating a more robust mathematical framework. This document presents a comprehensive 3D (2-DOF) derivation of a transformation matrix that accurately converts between gimbaled engine torques and TVC actuator forces. By incorporating the necessary geometric and rotational considerations, this new approach corrects the limitations of the traditional TVC MA method. Subsequent sections compare the newly formulated approach to the traditional method, demonstrating its enhanced accuracy and reliability for a broad range of gimbaled engine conditions.
The electrostatic force exerted on a targeting sequence by the electrical potential across the inner mitochondrial membrane is calculated and found to vary from 1.4 pN to 2.2 pN (per unit elementary charge) as the radius of the inner membrane pore (assumed aqueous) is varied from 12 to 6.5 Angstroms, its measured range. Since the pore is not very much wider than the distance between water molecules, the full shielding effect of water may not be present; the extreme case of a nonaqueous pore gives a force of 3.1 pN per unit charge, which represents an upper limit. When applied to mitochondrial import experiments on the protein harness, these results imply that a force of 11 plus or minus 4 pN is sufficient to catalyze the unfolding of harness during import. Comparison of these results with unfolding forces measured using atomic force microscopy suggests that the two are not inconsistent.
The governing equations for the analysis of open branch-chain mechanical systems are developed in a form suitable for implementation in a general purpose finite element computer program. Lagrange's form of d'Alembert's principle is used to derive the system mass matrix and force vector. The generalized coordinates are selected as the unconstrained relative degrees of freedom giving the position and orientation of each slave link with respect to their master link. Each slave link may have from zero to six degrees of freedom relative to the reference frames of its master link. A strategy for automatic generation of the system mass matrix and force vector is described.
An efficient procedure for obtaining the compatibility conditions of finite-element models involves the generation of both field and compatibility conditions from deformation-displacement relations, using (1) the compatibility bandwidth, and (2) the node-determinacy concept. A computer program thus structured will generate sparse and banded compatibility conditions for a structure that is idealized by the finite elements.
Explore the source record for details and available documents.