Search NASA⌕ Search

SEARCH · Search NASA

Results for “stiffness 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Node-recording method for stiffness matrix wavefront reduction in structural analysis

Method provides approach to automatic node relabeling that is consistent with requirements of wavefront concept. Specific applications are in analysis of aircraft, building structures, radar and surveillance structures, bridges, etc., or any other structure that is studied with aid of large and complex analytical model. Minimum growth sequencing is effective, rapid, and capable of producing economies.

Levy, R.↗
NASA NTRS1976

Symmetric stiffness matrix for incompressible hyperelastic materials

Symmetric structure matrices are derived for solving plane strain and axisymmetric problems involving incompressible hyperelastic materials. An infinite hollow cylinder subjected to internal pressure is considered as an example. Displacement and hydrostatic pressure profiles are calculated using the Newton-Raphson iteration technique. The results are in good agreement with the exact curves.

Takamatsu, T.↗
NASA NTRS1986

Ultrasonic determination of the elastic constants of the stiffness matrix for unidirectional fiberglass epoxy composites

The elastic constants of a fiberglass epoxy unidirectional composite are determined by measuring the phase velocities of longitudinal and shear stress waves via the through transmission ultrasonic technique. The waves introduced into the composite specimens were generated by piezoceramic transducers. Geometric lengths and the times required to travel those lengths were used to calculate the phase velocities. The model of the transversely isotropic medium was adopted to relate the velocities and elastic constants.

Marques, E. R. C.↗
NASA NTRS1993

Identification of a second-order mechanical system model from a state space realization

Many system identification algorithms produce models in modern state variable form. When a mechanical system is identified one knows that there must be a model in second order form. In structural dynamics one is interested in either a modal model of the system, or in a model including a mass matrix, stiffness matrix and a damping matrix. In this paper, algorithms are developed that convert a modern state space realization into the above two representations of interest in mechanical systems. The algorithm can be used with any identification methods that produce a modern state variable representation such as the Eigensystem Realization Algorithm (ERA), its modified version by data correlation (ERA/DC), and combinations of these with Observer/Kalman Filter Identification (OKID). An algorithm is developed that allows one to identify the damping matrix in a model representation, so that one determines how near to modal the damping is in the system, and can therefore understand the degree to which the modes are coupled by the damping in the system. A second algorithm produced the mass, damping, and stiffness matrices, given the input and output matrices. It is assumed that the number of sensors (or the number of actuators) is greater than or equal to the number of modes in the system, and it is proved that this is a necessary condition to be able to uniquely identify these matrices. Experience with examples suggests that a singular value truncation involved in these algorithms can help one determine the true system order. This truncation has the benefit of knowledge of the correct form for a mechanical system model, which is not present in the modern control identification algorithm.

Tseng, Dong-Huei↗
NASA NTRS1971

A Stability Analysis of Cylindrical Panels Using a Finite Element Formulation

A cylindrical finite element suitable for the linear stability analysis of cylindrical shells is developed. Energy principles and variational methods lead to a problem formulation which lends itself to physical interpretations of the governing matrices of the finite element. By properly grouping the terms which result from taking the second variation of the potential energy of the element, it is possible to identify three distinct types of matrices. The first matrix is the conventional stiffness matrix; the second is an initial stress stiffness matrix; and the third is an initial displacement stiffness matrix. With the assumption of linearity, the buckling problem is stated in terms of the classical linear real eigenvalue equation. This problem formulation was programmed on the CDC 6600 series computer. The computer program is used to analyze the buckling of a variety of structures. Columns, arches, flat plates and curved panels with and without cutouts are considered.

Snyder, R. E.↗
NASA NTRS1973

A NASTRAN DMAP alter for determining a local stiffness modification to obtain a specified eigenvalue

A technique is described which has been programmed as a DMAP Alter to Rigid Format 3, for determining a stiffness matrix modification to obtain a specified eigenvalue for a structure. The stiffness matrix modifications allowable are those that can be described as the product of a single scalar variable and a matrix of constant coefficients input by the user. The program solves for the scalar variable multiplier which will yield a specified eigenvalue for the complete structure (provided it exists), makes the modification to the stiffness matrix, and proceeds in Rigid Format 3 to obtain the eigenvalues and eigenvectors of the modified structure.

Case, W. R., Jr.↗
NASA NTRS1996

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↗
NASA NTRS1997

Frictionless Contact of Multilayered Composite Half Planes Containing Layers With Complex Eigenvalues

A previously developed local-global stiffness matrix methodology for the response of a composite half plane, arbitrarily layered with isotropic, orthotropic or monoclinic plies, to indentation by a rigid parabolic punch is further extended to accommodate the presence of layers with complex eigenvalues (e.g., honeycomb or piezoelectric layers). First, a generalized plane deformation solution for the displacement field in an orthotropic layer or half plane characterized by complex eigenvalues is obtained using Fourier transforms. A local stiffness matrix in the transform domain is subsequently constructed for this class of layers and half planes, which is then assembled into a global stiffness matrix for the entire multilayered half plane by enforcing continuity conditions along the interfaces. Application of the mixed boundary condition on the top surface of the half plane indented by a rigid punch results in an integral equation for the unknown pressure in the contact region. The integral possesses a divergent kernel which is decomposed into Cauchy-type and regular parts using the asymptotic properties of the local stiffness matrix and a relationship between Fourier and finite Hilbert transform of the contact pressure. The solution of the resulting singular integral equation is obtained using a collocation technique based on the properties of orthogonal polynomials developed by Erdogan and Gupta. Examples are presented that illustrate the important influence of low transverse properties of layers with complex eigenvalues, such as those exhibited by honeycomb, on the load versus contact length response and contact pressure distributions for half planes containing typical composite materials.

Zhang, Wang↗
NASA NTRS1976

Stiffness coupling application to modal synthesis program, users guide

A FORTRAN IV computer program used to perform modal synthesis of structures by stiffness coupling, using the dynamic transformation method is described. The program was named SCAMP (Stiffness Coupling Approach Modal-Synthesis Program). The program begins with the entry of a substructure's physical mode shapes and eigenvalues or a substructure's mass and stiffness matrix. If the mass and stiffness matrices are entered, the eigen problem for the individual substructure is solved. Provisions are included for a maximum of 20 substructures which are coupled by stiffness matrix springs. Each substructure has a number degrees of freedom (DOF), except that for DOF greater than 100; vector sets having maximum row and column size of 100 were generated prior to entering SCAMP. The substructures are then coupled together via coupling springs, and the dynamic transformation is used to reduce the size of the eigen problem.

Kuhar, E. J.↗
NASA NTRS1990

Vibration transmission through rolling element bearings in geared rotor systems

A new mathematical model is proposed to examine the vibration transmission through rolling element bearings in geared rotor systems. Current bearing models, based on either ideal boundary conditions for the shaft or purely translational stiffness element description, cannot explain how the vibratory motion may be transmitted from the rotating shaft to the casing. Experimental results have shown that the casing plate motion is primarily flexural. Here, this issue is clarified qualitatively and quantitatively by developing a comprehensive bearing stiffness matrix of dimension 6 to model precision rolling element bearings using basic principles. The proposed bearing stiffness matrix is partially verified using available analytical and experimental data, and is completely characterized. The study extends the proposed bearing formulation to analyze the overall geared rotor system dynamics including casing and mounts. The bearing stiffness matrix is included in discrete system models using lumped parameter and/or dynamic finite element techniques.

Singh, Rajendra↗
NASA NTRS1982

Calculation of exact vibration modes for plane grillages by the dynamic stiffness method

A dynamic stiffness method is developed for the calculation of the exact modal parameters for plane grillages which consist of straight and uniform beams with coincident elastic and inertial axes. Elementary bending-torsion beam theory is utilized, and bending translation is restricted to one direction. The exact bending-torsion dynamic stiffness matrix is obtained for a straight and uniform beam element with coincident elastic and inertial axes. The element stiffness matrices are assembled using the standard procedure of the static stiffness method to form the dynamic stiffness matrix of the complete grillage. The exact natural frequencies, mode shapes, and generalized masses of the grillage are then calculated by solving a nonlinear eigenvalue problem based on the dynamic stiffness matrix. The exact modal solutions for an example grillage are calculated and compared with the approximate solutions obtained by using the finite element method.

Hallauer, W. L., Jr.↗
NASA NTRS1973

Placing three-dimensional isoparametric elements into NASTRAN

Linear (8 node), parabolic (20 node), cubic (32 node) and mixed (some edges linear, some parabolic and some cubic) have been inserted into NASTRAN, level 15.1. First the dummy element feature was used to check out the stiffness matrix generation routines for the linear element in NASTRAN. Then, the necessary modules of NASTRAN were modified to include the new family of elements. The matrix assembly was changed so that the stiffness matrix of each isoparametric element is only generated once as the time to generate these higher order elements tends to be much longer than the other elements in NASTRAN. This paper presents some of the experiences and difficulties of inserting a new element or family of elements into NASTRAN.

Newman, M. B.↗
NASA NTRS1984

Parallel triangularization of substructured finite element problems

Much of the computational effort of the finite element process involves the solution of a system of linear equations. The coefficient matrix of this system, known as the global stiffness matrix, is symmetric, positive definite, and generally sparse. An important technique for reducing the time required to solve this system is substructuring or matrix partitioning. Substructuring is based on the idea of dividing a structure into pieces, each of which can then be analyzed relatively indepenently. As a result of this division, each point in the finite element discretization is either interior to a substructure or on a boundary between substructures. Contributions to the global stiffness matrix from connections between boundary points from the K(bb) matrix are reported. The triangularization of a general K(bb) matrix on a parallel machine is specifically discussed.

Leuze, M. R.↗
NASA NTRS1985

Improved accuracy in solid isoparametric elements using selectively reduced integration

Currently COSMIC NASTRAN uses the same Gaussian integration order for every term in the stiffness matrix of the solid isoparametric elements. This results in overly stiff elements under some conditions. A modification to the CIHEXI, CIHEX2 and CIHEX3 elements is described which provides the user with the capability of controlling the order of integration for selected terms in the stiffness matrix. Examples of improved accuracy using selectively reduced integration be provided are provided.

Case, W. R.↗
NASA NTRS1993

Vibration of thermally buckled composite plates with initial deflections using triangular elements

A consistent finite element formulation is presented for the analysis of thermal postbuckling and free vibration of thermally buckled thin, laminated composite plates subjected to large temperature change. The influence of moderately large initial imperfections in deflection on the thermal postbuckling deflection and the vibration characteristics of the buckled plate is also investigated. The finite element equations of motion are derived from the principle of virtual work. These equations can be mathematically separated into two sets and solved in sequence. The first set of equations yields the particular solution of static thermal postbuckling deflection, and the second set of equations gives the homogeneous solution of vibration characteristics on the buckled plate. The first set of static equations is solved by using Newton-Raphson iteration method. The tangent stiffness matrix in the final iteration is equal to the total stiffness matrix of the second set of dynamic equations. This feature saves tremendous computation time in comparing with using the conventional approach. The influence of lamination angle, temperature distribution, plate planform of arbitrary shape, and boundary support conditions on postbuckling and vibration behavior are investigated.

Zhou, R. C.↗
NASA NTRS1983

A refined finite element for vibration analysis of twisted blades based on beam theory

A finite element method of discretizing beam segments of pretwisted rotating blades is presented. Employing the matrix displacement method, stiffness and mass properties are developed from basic mechanics of a pretwisted beam theory. By introducing the proper displacement functions, the effect of rotor blade rotational motion on the stiffness matrix is obtained systematically from the kinetic energy expression. Comparing with other beam elements the derivation of this element is more fundamental. This allows one to apply the same approach to more complicated problems including nonlinear effects or complex dynamic motions. Illustrative examples are given comparing numerical results with available data and other numerical solutions from rotating and nonrotating force fields. These examples show that accurate prediction of vibration frequencies for pretwisted blades can be obtained by employing a quite modest number of degrees of freedom.

Sisto, F.↗
NASA NTRS1984

System identification of analytical models of damped structures

A procedure is presented for identifying linear nonproportionally damped system. The system damping is assumed to be representable by a real symmetric matrix. Analytical mass, stiffness and damping matrices which constitute an approximate representation of the system are assumed to be available. Given also are an incomplete set of measured natural frequencies, damping ratios and complex mode shapes of the structure, normally obtained from test data. A method is developed to find the smallest changes in the analytical model so that the improved model can exactly predict the measured modal parameters. The present method uses the orthogonality relationship to improve mass and damping matrices and the dynamic equation to find the improved stiffness matrix.

Fuh, J.-S.↗