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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1985

High temperature metal and ceramic composites

The Materials Division at NASA Lewis is engaged in research and development efforts on behalf of fiber-reinforced composite materials that are lighter, stiffer, and more structurally reliable than conventional monolithic alloys and ceramics in applications that range from the cryogenic to the refractory. Attention is presently given to metal matrix composites, in which high performance depends on stiff, strong and thermally stable large diameter fibers, with chemically stable interfacial bonding and good coefficient of thermal expansion matching between fibers and matrices, and to ceramic matrix composites, in which intermediate strength interfacial bonds must allow cracks to propagate through the matrix only, while retaining good load transfer characteristics between fiber and matrix.

Signorelli, R. A.↗

NASA NTRS1988

A proposed method for enhanced eigen-pair extraction using finite element methods: Theory and application

The paper covers two distinct parts: theory and application. The goal of this work was the reduction of model size with an increase in eigenvalue/vector accuracy. This method is ideal for the condensation of large truss- or beam-type structures. The theoretical approach involves the conversion of a continuum transfer matrix beam element into an 'Exact' dynamic stiffness element. This formulation is implemented in a finite element environment. This results in the need to solve a transcendental eigenvalue problem. Once the eigenvalue is determined the eigenvectors can be reconstructed with any desired spatial precision. No discretization limitations are imposed on the reconstruction. The results of such a combined finite element and transfer matrix formulation is a much smaller FEM eigenvalue problem. This formulation has the ability to extract higher eigenvalues as easily and as accurately as lower eigenvalues. Moreover, one can extract many more eigenvalues/vectors from the model than the number of degrees of freedom in the FEM formulation. Typically, the number of eigenvalues accurately extractable via the 'Exact' element method are at least 8 times the number of degrees of freedom. In contrast, the FEM usually extracts one accurate (within 5 percent) eigenvalue for each 3-4 degrees of freedom. The 'Exact' element results in a 20-30 improvement in the number of accurately extractable eigenvalues and eigenvectors.

Jara-Almonte, J.↗

NASA NTRS1982

Modeling stiffness loss in boron/aluminum below the fatigue limit

Boron/aluminum can develop significant internal matrix cracking when fatigued. These matrix cracks can result in a 40 percent secant modulus loss in some laminates, even when fatigued below the fatigue limit. It is shown that the same amount of fatigue damage will develop during stress or strain-controlled tests. Stacking sequence has little influence on secant modulus loss. The secant modulus loss in unidirectional composites is small, whereas the losses are substantial in laminates containing off-axis plies. A simple analysis is presented that predicts unnotched laminate secant modulus loss due to fatigue. The analysis is based upon the elastic modulus and Poisson's ratio of the fiber and matrix, fiber volume fraction, fiber orientations, and the cyclic-hardened yield stress of the matrix material. Excellent agreement was achieved between model predictions and experimental results. With this model, designers can project the material stiffness loss for design load or strain levels and assess the feasibility of its use in stiffness-critical parts.

Johnson, W. S.↗

NASA NTRS1995

Experimental Verification of a Progressive Damage Model for IM7/5260 Laminates Subjected to Tension-Tension Fatigue

The durability and damage tolerance of laminated composites are critical design considerations for airframe composite structures. Therefore, the ability to model damage initiation and growth and predict the life of laminated composites is necessary to achieve structurally efficient and economical designs. The purpose of this research is to experimentally verify the application of a continuum damage model to predict progressive damage development in a toughened material system. Damage due to monotonic and tension-tension fatigue was documented for IM7/5260 graphite/bismaleimide laminates. Crack density and delamination surface area were used to calculate matrix cracking and delamination internal state variables to predict stiffness loss in unnotched laminates. A damage dependent finite element code predicted the stiffness loss for notched laminates with good agreement to experimental data. It was concluded that the continuum damage model can adequately predict matrix damage progression in notched and unnotched laminates as a function of loading history and laminate stacking sequence.

Coats, Timothy W.↗

NASA NTRS2013

Improved Subcell Model for the Prediction of Braided Composite Response

In this work, the modeling of triaxially braided composites was explored through a semi-analytical discretization. Four unique subcells, each approximated by a "mosaic" stacking of unidirectional composite plies, were modeled through the use of layered-shell elements within the explicit finite element code LS-DYNA. Two subcell discretizations were investigated: a model explicitly capturing pure matrix regions, and a novel model which absorbed pure matrix pockets into neighboring tow plies. The in-plane stiffness properties of both models, computed using bottom-up micromechanics, correlated well to experimental data. The absorbed matrix model, however, was found to best capture out-of- plane flexural properties by comparing numerical simulations of the out-of-plane displacements from single-ply tension tests to experimental full field data. This strong correlation of out-of-plane characteristics supports the current modeling approach as a viable candidate for future work involving impact simulations.

Cater, Christopher R.↗

NASA NTRS1988

Shape sensitivity analysis of flutter response of a laminated wing

A method is presented for calculating the shape sensitivity of a wing aeroelastic response with respect to changes in geometric shape. Yates' modified strip method is used in conjunction with Giles' equivalent plate analysis to predict the flutter speed, frequency, and reduced frequency of the wing. Three methods are used to calculate the sensitivity of the eigenvalue. The first method is purely a finite difference calculation of the eigenvalue derivative directly from the solution of the flutter problem corresponding to the two different values of the shape parameters. The second method uses an analytic expression for the eigenvalue sensitivities of a general complex matrix, where the derivatives of the aerodynamic, mass, and stiffness matrices are computed using a finite difference approximation. The third method also uses an analytic expression for the eigenvalue sensitivities, but the aerodynamic matrix is computed analytically. All three methods are found to be in good agreement with each other. The sensitivities of the eigenvalues were used to predict the flutter speed, frequency, and reduced frequency. These approximations were found to be in good agreement with those obtained using a complete reanalysis.

Bergen, Fred D.↗

NASA NTRS1989

Shape sensitivity analysis of flutter response of a laminated wing

Kapania, Rakesh K.↗

NASA NTRS1983

Modeling stiffness loss in boron/aluminum laminates below the fatigue limit

Johnson, W. S.↗

NASA NTRS1994

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.↗

NASA NTRS1988

Fatigue damage criteria - Matrix, fibers and interfaces of continuous fiber reinforced metal matrix composites

Continuous fiber reinforced metal matrix composites (MMC) are projected for use in high temperature, stiffness critical parts that will be subjected to cyclic loadings. Depending on the relative fatigue behavior of the fiber and matrix, and the interface properties, the failure modes of MMC can be grouped into four catagories: (1) matrix dominated, (2) fiber dominated, (3) self-similar damage growth, and (4) fiber/matrix interfacial failures. These four types of damage are discussed and illustrated by examples. The emphasis is on the fatigue of unnotched laminates.

Johnson, W. S.↗

NASA NTRS1976

On a numerical solution of the supersonic panel flutter eigenproblem

An automated digital computer procedure is presented in this paper which enables efficient solution of the eigenvalue problem associated with the supersonic panel flutter phenomena. The step-by-step incremental solution procedure is based on an inverse iteration technique which effectively utilizes solution results from the previous step in determining such results during the current solution step. Also, the computations are limited to the determination of a few specific roots only, which are expected to contain the flutter mode, and this is achieved at each step without having to compute any other root. The structural discretization achieved by the finite-element method yields highly banded stiffness, mass, and aerodynamic matrices; the aerodynamic matrix evaluated by the linearized piston theory is real but unsymmetric in nature. The solution algorithm presented in this paper fully exploits the banded form of the associated matrices, and the resulting computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most efficient and economical when compared to existing procedures of such analysis. Numerical results are presented for a two-dimensional panel flutter problem.

Gupta, K. K.↗

NASA NTRS1977

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗

NASA NTRS1983

Optimal trajectories for aeroassisted orbital transfer

Consideration is given to classical and minimax problems involved in aeroassisted transfer from high earth orbit (HEO) to low earth orbit (LEO). The transfer is restricted to coplanar operation, with trajectory control effected by means of lift modulation. The performance of the maneuver is indexed to the energy expenditure or, alternatively, the time integral of the heating rate. Firist-order optimality conditions are defined for the classical approach, as are a sequential gradient-restoration algorithm and a combined gradient-restoration algorithm. Minimization techniques are presented for the aeroassisted transfer energy consumption and time-delay integral of the heating rate, as well as minimization of the pressure. It is shown that the eigenvalues of the Jacobian matrix of the differential system is both stiff and unstable, implying that the sequential gradient restoration algorithm in its present version is unsuitable. A new method, involving a multipoint approach to the two-poing boundary value problem, is recommended.

Miele, A.↗

NASA NTRS1984

Rheological characterization of addition polyimide matrix resins and prepregs

Although graphite-reinforced polyimide matrix composites offer outstanding specific strength and stiffness, together with high thermal oxidative stability, processing problems connected with their rheological behavior remain to be addressed. The present rheological studies on neat polyimide resin systems encountered outgassing during cure. A staging technique has been developed which can successfully handle polyimide samples, and novel methods were applied to generate rheological curves for graphite-reinforced prepregs. The commercial graphite/polyimide systems studied were PRM 15, LARC 160, and V378A.

Maximovich, M. G.↗

NASA NTRS1988

A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach

A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.

Koenig, Herbert A.↗

NASA NTRS1990

Environment enhanced fatigue of advanced aluminum alloys and composites

The objective is to characterize and understand the environmental fatigue crack propagation behavior of advanced, high stiffness and strength, aluminum alloys and metal matrix composites. Those gases and aqueous electrolytes which are capable of producing atomic hydrogen by reactions on clean crack surfaces are emphasized. Characterizations of the behavior of new materials are sought to provide data for damage tolerant component life prediction. Mechanistic models are sought for crack tip damage processes which are generally applicable to structural aluminum alloys. Such models will enable predictions of cracking behavior outside of the data, metallurgical improvements in material cracking resistance, and insight on hydrogen compatibility.

Slavik, Donald C.↗

NASA NTRSDate not supplied

Metallic Powder Core Tubular Wire (PCTW) Development for Additive Manufacturing

NASA Langley Research Center (LaRC) in Hampton, VA, and the Center for Welding, Joining, and Coating Research (CWJCR) in the Department of Metallurgical and Materials Engineering at the Colorado School of Mines (CSM) in Golden, CO, have collaborated for nearly fifteen years to develop metallic powder core tubular wire (PCTW) feedstock for additive manufacturing (AM) using electron beam deposition. The motivation for this development was to offset dealloying observed during deposition of titanium and aluminum alloy solid wire feedstock in the LaRC electron beam freeform fabrication (EBF3) system and to explore novel alloy compositions by tailoring the powder fill. EBF3 is a fusion-based AM process that is performed under vacuum and is consequently prone to vaporization loss of low vapor pressure alloying elements. Research spanning four graduate student research projects succeeded in developing the basic methodology to fabricate PCTW, the mass balance equations used to define the powder fill, microstructure control through particulate inoculation, and the feasibility of producing metal matrix composite materials. Metallic PCTW is comprised of a thin foil sheath wrapped around a powder blend. Fabrication of PCTW at CSM was accomplished using a tubular wire mill system to first bend a flat strip into a u-shaped cross section, deposit a powder fill, create a lap joint to seal the wire, and draw to final diameter. Powder management systems were developed to ensure uniform mixing of powder constituents and identify feed rates for optimum fill ratios. Processing parameters were successfully developed for titanium and aluminum alloy sheaths to account for material work hardening during drawing. Mass balance equations were developed to account for compositional contributions of the sheath and powder fill. PCTW composition can be adjusted through selection of the sheath or powder materials. PCTWs were successfully developed to offset Al loss in deposits of Ti-6-4 and Mg loss in Al 6061, resulting in deposits that were within the composition limits for each material. Ti-6-4 deposits were fabricated using PCTW comprised of commercially pure (CP) titanium sheath and a fill of pre-alloyed Ti-6-4 powder with additions of elemental Al and V powders. Additions of both elements were required to account for the CP Ti from the sheath, and Al was further enhanced to offset documented losses in EBF3 deposits. For Al 6061 deposits, composition was modulated through the choice of sheath material. The PCTW was made using an Al 5052 sheath to supplement Mg content and a fill composition of pre-alloyed 6061 powder with additions of Si and Cu.The slower solidification rate of the EBF3 process compared with laser powder bed fusion processes can lead to coarse-grained microstructures and mechanical property anisotropy. EBF3deposits of Ti-6-4 exhibited large epitaxially grown β grains, which result in degraded mechanical properties in the deposition direction. Microstructure refinement was achieved through a combination of electron beam modulation and the use of a PCTW designed for Ti-6-4 deposits and modified by the addition of powder made from Fe and B. Resulting deposits exhibited reductions in β grain size, α phase colony intercept length and lath width, and an increased volume fraction of β phase. The mechanism of microstructure refinement was related to precipitation of boride particles. Hardness testing confirmed that mechanical property anisotropy was reduced and an overall increase in hardness was achieved.Metal matrix composite (MMC) materials offer strength and stiffness benefits over unreinforced materials. Two methods were investigated using PCTW technology for producing Al 6061 basedMMC materials; the direct addition of SiC reinforcing particles and the in-situ generation of reinforcements through reaction synthesis of precursor powders. Al 6061 based MMCs reinforced with SiC exhibited particle clustering and formation of aluminum carbides during deposition. The application of a Ni coating to the SiC particles prior to incorporation into the powder fill eliminated aluminum carbide formation and aided particle distribution during deposition. Modulation of beam focus and power further reduced aluminum carbide formation and aided particle distribution. In reaction synthesis, precursor additions undergo an exothermic reaction during deposition to form new ceramic products, generally of finer scale than the initial precursor components. The product phases provide reinforcement in the MMC and act as nucleation sites during solidification, resulting in grain refinement, improved mechanical properties, and reduced susceptibility to solidification cracking. Higher tensile strength and stiffness were demonstrated in Al 6061 based MMCs with 2 vol% reinforcement. These combined research efforts demonstrate the potential of PCTW technology for customizing of alloy composition and achieving microstructure control. A manufacturing process for titanium- and aluminum-based PCTWs was developed based on iterative design, theoretical calculations, and empirical data. PCTW technology successfully mitigated Al loss from Ti-6-4 and Mg loss in Al 6061, resulted in grain refinement in Ti-6-4, and demonstrated fabrication of Al based MMCs

Marcia S Domack↗

NASA NTRSDate not supplied

Stiffness and Fatigue Life Estimator for Polymer Composite Laminates Using Machine Learning

Machine learning (ML) models are increasingly being used in many engineering fields due to the advancements in ML algorithms and availability of high-speed computing power. One of the most popular ML class of models is artificial neural networks (ANN). ML is increasingly being used in the design and analysis of composite materials and structures, specifically in the constitutive modeling of composite materials with the focus on greatly accelerating multiscale analyses of composite materials and structures through development of surrogate models. Towards that end, Python-based neural nets have been developed to predict initial stiffness and fatigue life of an eight-ply symmetric polymer matrix composite laminate. Two types of neural networks, a Multilayer Perceptron (MLP) and a Recurrent Neural Network (RNN), have been established. Results show that both neural net type algorithms can provide an excellent estimate of initial laminate stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminates (PMCs). RNNs are better able to capture the shape of the fatigue curve of a laminate. The resulting tool and GUI can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. Further, the associated surrogate models can also be used in composite multiscale analyses to replace the actual physics-based calculations at lower scales and thereby significantly increase the computational efficiency of such analyses and thus make micromechanics-based multiscale analyses a viable industrial tool for large scale structural problems.

multiscale analysis↗