Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Reliability analysis of laminated CMC components through shell subelement techniques

An updated version of the integrated design program Composite Ceramics Analysis and Reliability Evaluation of Structures (C/CARES) was developed for the reliability evaluation of ceramic matrix composites (CMC) laminated shell components. The algorithm is now split into two modules: a finite-element data interface program and a reliability evaluation algorithm. More flexibility is achieved, allowing for easy implementation with various finite-element programs. The interface program creates a neutral data base which is then read by the reliability module. This neutral data base concept allows easy data transfer between different computer systems. The new interface program from the finite-element code Matrix Automated Reduction and Coupling (MARC) also includes the option of using hybrid laminates (a combination of plies of different materials or different layups) and allows for variations in temperature fields throughout the component. In the current version of C/CARES, a subelement technique was implemented, enabling stress gradients within an element to be taken into account. The noninteractive reliability function is now evaluated at each Gaussian integration point instead of using averaging techniques. As a result of the increased number of stress evaluation points, considerable improvements in the accuracy of reliability analyses were realized.

Starlinger, Alois↗
NASA NTRS1978

NASTRAN implementation of an isoparametric doubly-curved quadrilateral shell element

A quadrilateral shell element, CQUAD4, was added to level 15.5 and subsequently to level 16.0 of NASTRAN. The element exhibited doubly curved surfaces and used biquadratic interpolation functions. Reduced integration techniques were used to improve the performance of the element in thin shell problems. The creation of several new bulk data items is discussed, along with a special module, GPNORM, to process SHLNORM bulk data cards. In addition to the theoretical basis for the element stiffness matrix, consistent mass and load matrices are presented. Several potential sources of degenerate behavior of the element were investigated. Guidelines for proper use of the element were suggested. Performance of the element on several widely published classical examples was demonstrated. The results showed a significant improvement over presently available NASTRAN shell elements for even the coarsest meshes. Potential applications to two classes of practical problems are discussed.

Potvin, A. B.↗
NASA NTRS1982

Alternative ways for formulation of hybrid stress elements

An element stiffness matrix can be derived by the conventional potential energy principle and, indirectly, also by generalized variational principles, such as the Hu-Washizu principle and the Hellinger-Reissner principle. The present investigation has the objective to show an approach which is concerned with the formulation of incompatible elements for solid continuum and for plate bending problems by the Hellinger-Reissner principle. It is found that the resulting scheme is equivalent to that considered by Tong (1982) for the construction of hybrid stress elements. In Tong's scheme the inversion of a large flexibility matrix can be avoided. It is concluded that the introduction of additional internal displacement modes in mixed finite element formulations by the Hellinger-Reissner principle and the Hu-Washizu principle can lead to element stiffness matrices which are equivalent to the assumed stress hybrid method.

Pian, T. H. H.↗
NASA NTRS1987

An efficient quadrilateral element for plate bending analysis

A simple, shear flexible, quadrilateral plate element is developed based on the Hellinger/Reissner mixed variational principle with independently assumed displacement and stress fields. The crucial point of the selection of appropriate stress parameters is emphasized in the formulation. For this purpose, a set of guidelines is formulated based on the follwoing considerations: (1) suppression of all kinematic deformation modes, (2) the element has a favorable value for the constraint index in the thin plate limit, and (3) element properties are frame-invariant. For computer implementation the components of the element stiffness matrix are evaluated analytically using the symbolic manipulation package MACSYMA. The effectiveness and practical usefulness of the proposed element are demonstrated by the numerical results of a variety of problems involving thin and moderately thick plates under different loading and support conditions.

Saleeb, A. F.↗
DOE OSTI2022

An energy-based study of the embedded element method for explicit dynamics

Abstract The embedded finite element technique provides a unique approach for modeling of fiber-reinforced composites. Meshing fibers as distinct bundles represented by truss elements embedded in a matrix material mesh allows for the assignment of more specific material properties for each component rather than homogenization of all of the properties. However, the implementations of the embedded element technique available in commercial software do not replace the material of the matrix elements with the material of the embedded elements. This causes a redundancy in the volume calculation of the overlapping meshes leading to artificially increased stiffness and mass. This paper investigates the consequences in the energy calculations of an explicit dynamic model due to this redundancy. A method for the correction of the edundancy within a finite element code is suggested which removes extra energy and is shown to be effective at correcting the energy calculations for large amounts of redundant volume.

Martin, Valerie A.↗
NASA NTRS1988

Improved finite strip Mindlin plate bending element using assumed shear strain distributions

A linear finite strip plate element based on Mindlin/Reissner plate theory is developed. The analysis is suitable for both thin and thick plates. In the formulation new transverse shear strains are introduced and assumed constant in each two-code linear strip. The element stiffness matrix is explicitly formulated for efficient computation and computer implementation. Numerical results showing the efficiency and predictive capability of the element for the analysis of plates are presented for different support and loading conditions and a wide range of thicknesses. No sign of shear locking phenomenon was observed with the newly developed element.

Chulya, Abhisak↗
NASA NTRS1989

Assumed strain distributions for a finite strip plate bending element using Mindlin-Reissner plate theory

A linear finite strip plate element based on Mindlin-Reissner plate theory is developed. The analysis is suitable for both thin and thick plates. In the formulation, new transverse shear strains are introduced and assumed constant in each two-node linear strip. The element stiffness matrix is explicitly formulated for efficient computation and computer implementation. Numerical results showing the efficiency and predictive capability of the element for the analysis of plates are presented for different support and loading conditions and a wide range of thicknesses. No sign of shear locking is observed with the newly developed element.

Chulya, Abhisak↗
NASA NTRS2005

Alloy NASA-HR-1

NASA-HR-1 is a high-strength Fe-Ni-base superalloy that resists high-pressure hydrogen environment embrittlement (HEE), oxidation, and corrosion. Originally derived from JBK-75, NASA-HR-1 has exceptional HEE resistance that can be attributed to its gamma-matrix and eta-free (Ni3Ti) grain boundaries. The chemistry was formulated using a design approach capable of accounting for the simultaneous effects of several alloy additions. This approach included: (1) Systematically modifying gamma-matrix compositions based on JBK-75; (2) Increasing gamma (Ni3(Al,Ti)) volume fraction and adding gamma-matrix strengthening elements to obtain higher strength; and (3) Obtaining precipitate-free grain boundaries. The most outstanding attribute of NASA-HR-1 is its ability to resist HEE while showing much improved strength. NASA-HR-1 has approximately 25% higher yield strength than JXK-75 and exhibits tensile elongation of more than 20% with no ductility loss in a hydrogen environment at 5 ksi, an achievement unparalleled by any other commercially available alloy. Its Cr and Ni contents provide exceptional resistance to environments that promote oxidation and corrosion. Microstructural stability was maintained by improved solid solubility of the gamma-matrix, along with the addition of alloying elements to retard eta (Ni3Ti) precipitation. NASA-HR-1 represents a new system that greatly extends the compositional ranges of existing HEE-resistant Fe-Ni-base superalloys.

Chen, Po-Shou↗
NASA NTRS1988

An efficient Mindlin finite strip plate element based on assumed strain distribution

A simple two node, linear, finite strip plate bending element based on Mindlin-Reissner plate theory for the analysis of very thin to thick bridges, plates, and axisymmetric shells is presented. The new transverse shear strains are assumed for constant distribution in the two node linear strip. The important aspect is the choice of the points that relate the nodal displacements and rotations through the locking transverse shear strains. The element stiffness matrix is explicitly formulated for efficient computation and ease in computer implementation. Numerical results showing the efficiency and predictive capability of the element for analyzing plates with different supports, loading conditions, and a wide range of thicknesses are given. The results show no sign of the shear locking phenomenon.

Chulya, Abhisak↗
NASA NTRS1986

Reduced modeling of flexible structures for decentralized control

Based upon the modified finite element-transfer matrix method, this paper presents a technique for reduced modeling of flexible structures for decentralized control. The modeling decisions are carried out at (finite-) element level, and are dictated by control objectives. A simply supported beam with two sets of actuators and sensors (linear force actuator and linear position and velocity sensors) is considered for illustration. In this case, it is conjectured that the decentrally controlled closed loop system is guaranteed to be at least marginally stable.

Yousuff, A.↗
NASA NTRS1986

Triangular element for analysis of perforated plates under inplane and transverse loads

A C0-type triangular element formulation in orthogonal curvilinear coordinates has been developed, based on assumptions of transverse inextensibility and constant shear angle through thickness, for analysis of perforated plates subjected to in-plane and transverse loads. The assumed quadratic-displacement potential-energy approach is utilized in obtaining an element stiffness matrix and consistent load vector, which are numerically integrated. Numerical results have been obtained using a straight-sided triangular version, which behaves like a subparametric element, for stretching and bending analyses of perforated plates.

Chaudhuri, R. A.↗
DOE OSTI2024

Graph network heterogeneity predicts interplant wake losses

Wind plants generate large-scale wakes, which can affect the performance of neighboring installations. Such wakes are challenging to estimate due to the inherent complexity in modeling wake interactions between large quantities of turbines at various distances. Weighted directed graph networks can inform complex models by linking turbine pairs into chains of upstream and downstream neighbors for a given wind direction. A novel interpretation of the graph network adjacency matrix is proposed where each element of the matrix represents the cumulative impact of upstream turbines on an individual. In this study, wake losses were estimated with an engineering wake model across a range of inflow conditions for nine parametric variations of a system containing two neighboring wind plants. The parametric nature of the study isolates turbine spacing within the plant, separation distance between plants, and wind direction as the main drivers of wake losses. Spatial heterogeneity is computed from the weighted average adjacency matrix of each plant arrangement. The proposed method is orders of magnitude faster than wake modeling and does not require detailed turbine information or atmospheric conditions. Furthermore, the weighted average adjacency matrix provides insight on the spatial organization of wake losses at various scales. Plant heterogeneity is correlated with wake losses within and among plants. Framing wind plant wake interaction in terms of graph network spatial heterogeneity provides an efficient approach for predicting wake losses within and among neighboring wind plants with applications to other complex systems where wake interactions are key factor.

17 WIND ENERGY↗
NASA NTRS1976

Multivariate normal integration

Monte Carlo program evaluates integrals over rectangular regions for dimensions less than six and over elliptical regions in bivariate case. Program gives positive definite symmetric variance/covariance matrix factorization and calculates reciprocal of lower triangular matrix and product of diagonal elements of triangular matrix.

Falls, L. W.↗
NASA NTRS1986

Triangular element for analysis of a stretched plate weakened by a part-through hole

A C0-type triangular element formulation in orthogonal curvilinear coordinates has been developed, based on assumptions of transverse inextensibility and layerwise constant shear angle, for analysis of a stretched homogeneous plate weakened by a part-through hole. The element stiffness matrix and consistent load vector have been derived using an assumed-quadratic-displacement (in the curvilinear coordinate plane) potential-energy approach. Numerical results obtained using the straight-sided triangular element version indicate the presence of transverse shear deformation in a stretched homogeneous plate weakened by a symmetrically located concentric hole.

Chaudhuri, R. A.↗
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 NTRS2001

Free Vibration of Uncertain Unsymmetrically Laminated Beams

Monte Carlo Simulation and Stochastic FEA are used to predict randomness in the free vibration response of thin unsymmetrically laminated beams. For the present study, it is assumed that randomness in the response is only caused by uncertainties in the ply orientations. The ply orientations may become random or uncertain during the manufacturing process. A new 16-dof beam element, based on the first-order shear deformation beam theory, is used to study the stochastic nature of the natural frequencies. Using variational principles, the element stiffness matrix and mass matrix are obtained through analytical integration. Using a random sequence a large data set is generated, containing possible random ply-orientations. This data is assumed to be symmetric. The stochastic-based finite element model for free vibrations predicts the relation between the randomness in fundamental natural frequencies and the randomness in ply-orientation. The sensitivity derivatives are calculated numerically through an exact formulation. The squared fundamental natural frequencies are expressed in terms of deterministic and probabilistic quantities, allowing to determine how sensitive they are to variations in ply angles. The predicted mean-valued fundamental natural frequency squared and the variance of the present model are in good agreement with Monte Carlo Simulation. Results, also, show that variations between plus or minus 5 degrees in ply-angles can affect free vibration response of unsymmetrically and symmetrically laminated beams.

Kapania, Rakesh K.↗
NASA NTRS1989

Part B: Pattern control of horn antennas

During this period, the computations of the impedance elements were completed. These include interactions between the two electric current modes, the elecric current mode and the magnetic current mode, and the two magnetic current modes. An accurate and efficient formulation of computing interactions between electric current mode and magnetic current mode was accomplished. This, together with other subroutines allows for the fill-in of all the elements in the matrix. After the fill-in of the impedance elements in the matrix, the forward problem is accomplished. That is, given the specification of the horn and the excitating waveguide mode, the radiation pattern of the antenna based on the integral equation can be obtained. An example case was run for a standard X-band gain-horn (DBG-520). The H- and E-plane patterns of this horn antenna with perfectly conducting walls are compared with the gain pattern available from the manufacturer for up to the first side lobe. Good agreements are obtained although the cross polarization has not yet been accounted for. The effect of the lossy coating on the radiation pattern was also investigated. The resulting E-plane pattern shows about 3-dB improvement in the first sidelobe and 4-dB improvement in the second sidelobe.

Balanis, Constantine A.↗