Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient stiffness”

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 91 records · Page 5

A survey of the core-congruential formulation for geometrically nonlinear TL finite elements

This article presents a survey of the core-congruential formulation (CCF) for geometrically nonlinear mechanical finite elements based on the total Lagrangian (TL) kinematic description. Although the key ideas behind the CCF can be traced back to Rajasekaran and Murray in 1973, it has not subsequently received serious attention. The CCF is distinguished by a two-phase development of the finite element stiffness equations. The initial phase developed equations for individual particles. These equations are expressed in terms of displacement gradients as degrees of freedom. The second phase involves congruential-type transformations that eventually binds the element particles of an individual element in terms of its node-displacement degrees of freedom. Two versions of the CCF, labeled direct and generalized, are distinguished. The direct CCF (DCCF) is first described in general form and then applied to the derivation of geometrically nonlinear bar, and plane stress elements using the Green-Lagrange strain measure. The more complex generalized CCF (GCCF) is described and applied to the derivation of 2D and 3D Timoshenko beam elements. Several advantages of the CCF, notably the physically clean separation of material and geometric stiffnesses, and its independence with respect to the ultimate choice of shape functions and element degrees of freedom, are noted. Application examples involving very large motions solved with the 3D beam element display the range of applicability of this formulation, which transcends the kinematic limitations commonly attributed to the TL description.

Felippa, Carlos A.↗

Structural optimization of framed structures using generalized optimality criteria

The application of a generalized optimality criteria to framed structures is presented. The optimality conditions, Lagrangian multipliers, resizing algorithm, and scaling procedures are all represented as a function of the objective and constraint functions along with their respective gradients. The optimization of two plane frames under multiple loading conditions subject to stress, displacement, generalized stiffness, and side constraints is presented. These results are compared to those found by optimizing the frames using a nonlinear mathematical programming technique.

Kolonay, R. M.↗

Toughened High-Temperature Thermoset Composites

Laminated composites of carbon fibers and either PMR-15 or LARC(TM)RP-46 polyimides toughened selectively by incorporating thin layers of Matrimid(TM) 5218 thermoplastic polyimide at interfaces between plies to form gradient semi-interpenetrating microstuctures. Toughening of composites accompanied by acceptably small decreases in strength and stiffness. Thick laminates of toughened composites fabricated easily at 500 psi. Materials useful in components of aircraft engines, where exposed to high temperatures for short times.

Johnston, Norman J.↗

Control Force Compensation in Ground-Based Flight Simulators

This paper presents the results of a study that investigated if controller force compensations accounting for the inertial force and moment due to the aircraft motion during flight have a significant effect on pilot control behavior and performance. Seven rotorcraft pilots performed a side-step and precision hovering task in light turbulence in the Vertical Motion Simulator. The effects of force compensation were examined for two different simulated rotorcraft: linear and UH-60 dynamics with two different force gradient of the lateral stick control. Four motion configurations were used: large motion, hexapod motion, fixed-base motion, and fixed-base motion with compensation. Control-input variables and task performance such as the time to translate to the designated hover position, station-keeping position errors, and handling qualities ratings were used as measures. Control force compensation enabled pilot control behavior and performance more similar to that under high- or medium-fidelity motion to some extent only. Control force compensation did not improve overall task performance considering both rotorcraft models at the same time. The control force compensation had effects on the linear model with lighter force gradient, but only a minimal effect on pilots? control behavior and task performance for the UH-60 model, which had a higher force gradient. This suggests that the control force compensation has limited benefits for controllers that have higher stiffness.

Chung, William W.↗

Control Force Compensation in Ground-Based Flight Simulators

This paper presents the results of a study that investigated if controller force compensations accounting for the inertial force and moment due to the aircraft motion during flight have a significant effect on pilot control behavior and performance. Seven rotorcraft pilots performed a side-step and precision hovering task in light turbulence in the Vertical Motion Simulator. The effects of force compensation were examined for two different simulated rotorcraft: linear and UH-60 dynamics with two different force gradient of the lateral stick control. Four motion configurations were used: large motion, hexapod motion, fixed-base motion, and fixed-base motion with compensation. Control-input variables and task performance such as the time to translate to the designated hover position, station-keeping position errors, and handling qualities ratings were used as measures. Control force compensation enabled pilot control behavior and performance more similar to that under high- or medium-fidelity motion to some extent only. Control force compensation did not improve overall task performance considering both rotorcraft models at the same time. The control force compensation had effects on the linear model with lighter force gradient, but only a minimal effect on pilots? control behavior and task performance for the UH-60 model, which had a higher force gradient. This suggests that the control force compensation has limited benefits for controllers that have higher stiffness.

Chung, William↗

Physics constrained learning for data-driven inverse modeling from sparse observations

Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Furthermore, our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.

97 MATHEMATICS AND COMPUTING↗

The effect of nonsymmetric pressure stiffness on the dynamic characteristics of Solid Rocket Motors

This paper discusses the effect of pressure on the dynamics of pre-stiffened structures such as the Advanced Solid Rocket Motor (ASRM). Previous work in which the stiffness terms resulting from constant pressure were derived has been extended to enable modeling of nonconstant pressure applied over nonenclosed volumes. These conditions will result in nonsymmetric terms in the global stiffness matrix which will not cancel out. Three new pressure stiffness elements incorporating these nonsymmetric terms have been implemented as dummy elements in COSMIC NASTRAN and have been tested on various simple examples as well as an existing ASRM NASTRAN finite element model. The results indicate that for all load cases of practical interest to the ASRM program, the nonsymmetric terms have very little effect on the dynamic characteristics. In addition, the pressure stiffness elements developed in the previous work which assumed constant pressure gave virtually the same results as the new elements even for problems in which the pressures are not constant. The original elements appear to work well as long as the pressure gradient across any individual element is no larger than about 0.75 psi/inch. The new elements are therefore most useful for determining the conditions under which the original pressure stiffness elements can be used.

Christensen, Eric R.↗

Kinematics and control algorithm development and simulation for a redundant two-arm robotic manipulator system

An efficient approach to cartesian motion and force control of a 7 degree of freedom (DOF) manipulator is presented. It is based on extending the active stiffness controller to the 7 DOF case in general and use of an efficient version of the gradient projection technique for solving the inverse kinematics problem. Cooperative control is achieved through appropriate configuration of individual manipulator controllers. In addition, other aspects of trajectory generation using standard techniques are integrated into the controller. The method is then applied to a specific manipulator of interest (Robotics Research T-710). Simulation of the kinematics, dynamics, and control are provided in the context of several scenarios: one pertaining to a noncontact pick and place operation; one relating to contour following where contact is made between the manipulator and environment; and one pertaining to cooperative control.

Hennessey, Michael P.↗

Parallel Aeroelastic Analysis Using ENSAERO and NASTRAN

A high fidelity parallel static structural analysis capability is created and interfaced to the multidisciplinary analysis package ENSAERO-MPI of Ames Research Center. This new module replaced ENSAERO's lower fidelity simple finite element and modal modules. Full aircraft structures may be more accurately modeled using the new finite element capability. Parallel computation is performed by breaking the full structure into multiple substructures. This approach is conceptually similar to ENSAERO's multi-zonal fluid analysis capability. The new substructure code is used to solve the structural finite element equations for each substructure in parallel. NASTRAN/COSMIC is utilized as a front end for this code. Its full library of elements can be used to create an accurate and realistic aircraft mode. It is used to create the stiffness matrices for each sub-structure. The new parallel code then uses an iterative preconditioned conjugate gradient method to solve the global structural equations for the sub-structure boundary nodes. Results are presented for a wing-body configuration.

Eldred, Lloyd B.↗

A comparison of experimental and theoretical results for leakage, pressure gradients, and rotordynamic coefficients for tapered annular gas seal

A brief review of current annular seal theory and a discussion of the predicted effect on stiffness of tapering the seal stator are presented. An outline of Nelson's analytical-computational method for determining rotordynamic coefficients for annular compressible-flow seals is included. Modifications to increase the maximum rotor speed of an existing air-seal test apparatus at Texas A&M University are described. Experimental results, including leakage, entrance-loss coefficients, pressure distributions, and normalized rotordynamic coefficients, are presented for four convergent-tapered, smooth-rotor, smooth-stator seals. A comparison of the test results shows that an inlet-to-exit clearance ratio of 1.5 to 2.0 provides the maximum direct stiffness, a clearance ratio of 2.5 provides the greatest stability, and a clearance ratio of 1.0 provides the least stability. The experimental results are compared to theoretical results from Nelson's analysis with good agreement. Test results for cross-coupled stiffness show less sensitivity of fluid prerotation than predicted.

Elrod, D. A.↗

Isotope identity experiments in JET-ILW with H and D L-mode plasmas

NBI-heated L-mode plasmas have been obtained in JET with the Be/W ITER-like wall (JET-ILW) in H and D, with matched profiles of the dimensionless plasma parameters, ρ*, ν*, β and q in the plasma core confinement region and same T i /T e and Z eff . The achieved isotope identity indicates that the confinement scale invariance principle is satisfied in the core confinement region of these plasmas, where the dominant instabilities are Ion Temperature Gradient (ITG) modes. The dimensionless thermal energy confinement time, Ω i τ E ,t h , and the scaled core plasma heat diffusivity, A χ eff /B T , are identical in H and D within error bars, indicating lack of isotope mass dependence of the dimensionless L-mode thermal energy confinement time in JET-ILW. Predictive flux driven simulations with JETTO-TGLF of the H and D identity pair is in very good agreement with experiment for both isotopes: the stiff core heat transport, typical of JET-ILW NBI heated L-modes, overcomes the local gyro-Bohm scaling of gradient-driven TGLF, explaining the lack of isotope mass dependence in the confinement region of these plasmas. The effect of E × B shearing on the predicted heat and particle transport channels is found to be negligible for these low beta and low momentum input plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mechanisms for bacterial gliding motility on soft substrates

The motility mechanism of certain prokaryotes has long been a mystery, since their motion, known as gliding, involves no external appendages. The physical principles behind gliding still remain poorly understood. Using myxobacteria as an example of such organisms, we identify here the physical principles behind gliding motility and develop a theoretical model that predicts a 2-regime behavior of the gliding speed as a function of the substrate stiffness. Our theory describes the elasto-capillary–hydrodynamic interactions between the membrane of the bacteria, the slime it secretes, and the soft substrate underneath. Defining gliding as the horizontal translation under zero net force, we find the 2-regime behavior is due to 2 distinct mechanisms of motility thrust. On mildly soft substrates, the thrust arises from bacterial shape deformations creating a flow of slime that exerts a pressure along the bacterial length. This pressure in conjunction with the bacterial shape provides the necessary thrust for propulsion. On very soft substrates, however, we show that capillary effects must be considered that lead to the formation of a ridge at the slime–substrate–air interface, thereby creating a thrust in the form of a localized pressure gradient at the bacterial leading edge. To test our theory, we perform experiments with isolated cells on agar substrates of varying stiffness and find the measured gliding speeds in good agreement with the predictions from our elasto-capillary–hydrodynamic model. Finally, the mechanisms reported here serve as an important step toward an accurate theory of friction and substrate-mediated interactions between bacteria proliferating in soft media.

59 BASIC BIOLOGICAL SCIENCES↗

Static and Dynamic Model Update of an Inflatable/Rigidizable Torus Structure

The present work addresses the development of an experimental and computational procedure for validating finite element models. A torus structure, part of an inflatable/rigidizable Hexapod, is used to demonstrate the approach. Because of fabrication, materials, and geometric uncertainties, a statistical approach combined with optimization is used to modify key model parameters. Static test results are used to update stiffness parameters and dynamic test results are used to update the mass distribution. Updated parameters are computed using gradient and non-gradient based optimization algorithms. Results show significant improvements in model predictions after parameters are updated. Lessons learned in the areas of test procedures, modeling approaches, and uncertainties quantification are presented.

Horta, Lucas G.↗

Thickness and orientational design for a maximum stiff membrane

Recent results from sensitivity analysis for strain energy with anisotropic elasticity are applied to thickness and orientational design of laminated membranes. Primarily, the first order gradients of the total elastic energy are used in an optimality criteria based method. This traditional method is shown to give slow convergence with respect to design parameters, although the convergence of strain energy is very good. To get a deeper insight into this rather general characteristic, second order derivatives are included and it is shown how they can be obtained by first order sensitivity analysis. Examples of only thickness design, only orientational design, and combined thickness--orientational design are presented.

Pedersen, Pauli↗

Influence of freestream values on k-omega turbulence model predictions

The Wilcox (1988, 1991) k-omega model for eddy-viscosity turbulence does not require damping functions in the viscous sublayer, and its equations are less stiff near the wall. It has been designed to predict the requisite wake length in equilibrium, adverse pressure-gradient boundary-layer flows. When applied to free shear layers, however, a strong dependency of its results on the freestream value of omega has been noted. This feature is presently investigated via the self-similar equations for incompressible equilibrium boundary layers.

Menter, F. R.↗

Stiffness, thermal expansion, and thermal bending formulation of stiffened, fiber-reinforced composite panels

A method is presented for formulating stiffness terms and thermal coefficients of stiffened, fiber-reinforced composite panels. The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. Nonlinear, temperature and load dependent constitutive material data of each laminate are used to 'build-up' the stiffened panel membrane, bending, and membrane-bending coupling stiffness terms and thermal coefficients. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN analysis program (FEA) is revealed that provides accurate solutions of entire airframes or engines with coarsely meshed models. An example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. The erroneous results and the correct ones produced from this method are compared to analysis from discretely meshed three-dimensional FEA.

Collier, Craig S.↗

Quasi-Optimal Schwarz Methods for the Conforming Spectral Element Discretization

Fast methods are proposed for solving the system K(sub N)x = b resulting from the discretization of self-adjoint elliptic equations in three dimensional domains by the spectral element method. The domain is decomposed into hexahedral elements, and in each of these elements the discretization space is formed by polynomials of degree N in each variable. Gauss-Lobatto-Legendre (GLL) quadrature rules replace the integrals in the Galerkin formulation. This system is solved by the preconditioned conjugate gradients method. The conforming finite element space on the GLL mesh consisting of piecewise Q(sub 1) elements produces a stiffness matrix K(sub h) that is spectrally equivalent to the spectral element stiffness matrix K(sub N). The action of the inverse of K(sub h) is expensive for large problems, and is therefore replaced by a Schwarz preconditioner B(sub h) of this finite element stiffness matrix. The preconditioned operator then becomes B(sub h)(exp -l)K(sub N). The technical difficulties stem from the nonregularity of the mesh. Tools to estimate the convergence of a large class of new iterative substructuring and overlapping Schwarz preconditioners are developed. This technique also provides a new analysis for an iterative substructuring method proposed by Pavarino and Widlund for the spectral element discretization.

Casarin, Mario↗

On the non-stiffness of edge transport in L-mode tokamak plasmas

Transport analyses using first-principle turbulence codes and 11/2 -D transport codes usually study radial transport properties between the tokamak plasma magnetic axis and a normalized minor radius around 0.8. In this region, heat transport shows significantly stiff properties resulting in temperature scalelength values (R∕L{sub T}) that are relatively independent of the level of the radial heat flux. We have studied experimentally in the tokamak à configuration variable [F. Hofmann et al., Plasma Phys. Controlled Fusion 36, B277 (1994)] the radial electron transport properties of the edge region, close to the last closed flux surface, namely, between ρ{sub V}=√(V/V{sub edge})=0.8 to 1. It is shown that electron transport is not stiff in this region and high R∕L{sub Te} values (∼20–40) can be attained even for L-mode confinement. We can define a “pedestal” location, already in L-mode regimes, where the transport characteristics change from constant logarithmic gradient, inside ρ{sub V} = 0.8, to constant gradient between 0.8 and 1.0. In particular, we demonstrate, with well resolved T{sub e} and n{sub e} profiles, that the confinement improvement with plasma current I{sub p}, with or without auxiliary heating, is due to this non-stiff edge region. This new result is used to explain the significant confinement improvement observed with negative triangularity, which could not be explained by theory to date. Preliminary local gyrokinetic simulations are now consistent with an edge, less stiff, region that is more sensitive to triangularity than further inside. We also show that increasing the electron cyclotron heating power increases the edge temperature inverse scalelength, in contrast to the value in the main plasma region. The dependence of confinement on density in ohmic plasmas is also studied and brings new insight in the understanding of the transition between linear and saturated confinement regimes, as well as of the density limit and appearance of a 2/1 tearing mode. The results presented in this paper provide an important new perspective with regards to radial transport in tokamak plasmas which goes beyond L-mode plasmas and explains some previous puzzling results. It is proposed that understanding the transport properties in this edge non-stiff region will also help in understanding the improved and high confinement edge properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗