Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRIX STRESS CALCULATION”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Micromechanical Modeling of the Thermal Expansion of Graphite/copper Composites with Nonuniform Microstructure

Two micromechanical models were developed to investigate the thermal expansion of graphite/copper (Gr/Cu) composites. The models incorporate the effects of temperature-dependent material properties, matrix inelasticity, initial residual stresses due to processing history, and nonuniform fiber distribution. The first model is based on the multiple concentric cylinder geometry, with each cylinder treated as a two-phase composite with a characteristic fiber volume fractions. By altering the fiber volume fraction of the individual cylinders, unidirectional composites with radially nonuniform fiber distributions can be investigated using this model. The second model is based on the inelastic lamination theory. By varying the fiber content in the individual laminae, composites with nonuniform fiber distribution in the thickness direction can be investigated. In both models, the properties of the individual regions (cylinders or laminae) are calculated using the method of cells micromechanical model. Classical incremental plasticity theory is used to model the inelastic response of the copper matrix at the microlevel. The models were used to characterize the effects of nonuniform fiber distribution on the thermal expansion of Gr/Cu. These effects were compared to the effects of matrix plasticity, choice of stress-free temperature, and slight fiber misalignment. It was found that the radially nonuniform fiber distribution has little effect on the thermal expansion of Gr/Cu but could become significant for composites with large fiber-matrix transverse CTE and Young's modulus mismatch. The effect of nonuniform fiber distribution in the through-thickness direction of a laminate was more significant, but only approached that of the stress-free temperature for the most extreme cases that include large amounts of bending. Subsequent comparison with experimental thermal expansion data indicated the need for more accurate characterization of the graphite fiber thermomechanical properties. Correlation with cyclic data revealed the presence of a mechanism not considered in the developed models. The predicted response did, however, exhibit ratcheting behavior that has been observed experimentally in Gr/Cu. Finally, simulation of the actual fiber distribution of particular specimens had little effect on the predicted thermal expansion.

Bednarcyk, Brett A.↗
NASA NTRS1994

Experimental verification of a progressive damage model for composite laminates based on continuum damage mechanics

Progressive failure is a crucial concern when using laminated composites in structural design. Therefore the ability to model damage and predict the life of laminated composites is vital. The purpose of this research was to experimentally verify the application of the continuum damage model, a progressive failure theory utilizing continuum damage mechanics, to a toughened material system. Damage due to tension-tension fatigue was documented for the IM7/5260 composite laminates. Crack density and delamination surface area were used to calculate matrix cracking and delamination internal state variables, respectively, to predict stiffness loss. A damage dependent finite element code qualitatively predicted trends in transverse matrix cracking, axial splits and local stress-strain distributions for notched quasi-isotropic laminates. The predictions were similar to the experimental data and it was concluded that the continuum damage model provided a good prediction of stiffness loss while qualitatively predicting damage growth in notched laminates.

Coats, Timothy William↗
NASA NTRS2001

Notch Sensitivity of Woven Ceramic Matrix Composites Under Tensile Loading: An Experimental, Analytical, and Finite Element Study

The stress concentrations associated with circular notches and subjected to uniform tensile loading in woven ceramic matrix composites (CMCs) have been investigated for high-efficient turbine engine applications. The CMC's were composed of Nicalon silicon carbide woven fabric in SiNC matrix manufactured through polymer impregnation process (PIP). Several combinations of hole diameter/plate width ratios and ply orientations were considered in this study. In the first part, the stress concentrations were calculated measuring strain distributions surrounding the hole using strain gages at different locations of the specimens during the initial portion of the stress-strain curve before any microdamage developed. The stress concentration was also calculated analytically using Lekhnitskii's solution for orthotropic plates. A finite-width correction factor for anisotropic and orthotropic composite plate was considered. The stress distributions surrounding the circular hole of a CMC's plate were further studied using finite element analysis. Both solid and shell elements were considered. The experimental results were compared with both the analytical and finite element solutions. Extensive optical and scanning electron microscopic examinations were carried out for identifying the fracture behavior and failure mechanisms of both the notched and notched specimens. The stress concentration factors (SCF) determined by analytical method overpredicted the experimental results. But the numerical solution underpredicted the experimental SCF. Stress concentration factors are shown to increase with enlarged hole size and the effects of ply orientations on stress concentration factors are observed to be negligible. In all the cases, the crack initiated at the notch edge and propagated along the width towards the edge of the specimens.

Haque, A.↗
NASA NTRS1985

Lubrication of nonconformal contacts

Minimum film thickness results for piezoviscous-rigid regime of lubrication are developed for a compressible Newtonian fluid with Roelands viscosity. The results provide a basis for the analysis and design of a wide range of machine elements operating in the piezoviscous-rigid regime of lubrication. A new numerical method of calculating elastic deformation in contact stresses is developed using a biquadratic polynomial to approximate the pressure distribution on the whole domain analyzed. The deformation of every node is expressed as a linear combination of the nodal pressures whose coefficients can be combined into an influence coefficient matrix. This approach has the advantages of improved numerical accuracy, less computing time and smaller storage size required for influence matrix. The ideal elastohydrodynamic lubrication is extended to real bearing systems in order to gain an understanding of failure mechanisms in machine elements. The improved elastic deformation calculation is successfully incorporated into the EHL numerical scheme. Using this revised numerical technique and the flow factor model developed by Patir and Cheng (1978) the surface roughness effects on the elastohydrodynamic lubrication of point contact is considered. Conditions typical of an EHL contact in the piezoviscous-elastic regime entrained in pure rolling are investigated. Results are compared with the smooth surface solutions. Experiments are conducted to study the transient EHL effects in instrument ball bearings.

Jeng, Y. R.↗
NASA NTRS1992

Stiffnesses by (TtN) ensemble molecular dynamics

Calculation of elastic constants requires initializing a lattice of atoms, an h matrix defining the computational cell boundaries, and the parameters defining the interatomic potentials, then numerically integrating the equations of motion for the atoms and the computational cell boundaries, calculating strains and stresses at each timestep, and updating at each timestep the average quantities required for calculation of fluctuations. Enforcements of the DYNAMO FORTRAN code included addition of documentation, improved memory management and data flow, and enabling use of several interatomic potentials. Documentation included headers for each subroutine, which document the parameters input to the routine, algorithms employed, and output.

Tibbits, Patrick↗
NASA NTRS1991

Analysis of delamination in cross ply laminates initiating from impact induced matrix cracking

Several two dimensional finite element analyses of (0 sub 2/90 sub 8/0 sub 2) glass/epoxy and graphite-epoxy composite laminates were performed to study some of the characteristics of damage development due to an impact load. A cross section through the thickness of the laminate with fixed ends, and carrying a transverse load in the center was analyzed. Inclined matrix cracks such as those produced by low velocity impact were modeled in the 90 deg ply group. The introduction of the matrix cracks caused large interlaminar tension and shear stresses in the vicinity of both crack tips in the 0/90 and 90/0 interfaces. The large interlaminar stresses at the ends of the matrix cracks indicate that matrix cracking may give rise to delamination. The ratio of mode I to total strain energy release rate at the beginning of delamination calculated at the two matrix crack tips was 60 and 28 pct., respectively, in the glass/epoxy laminate. The corresponding ratio was 97 and 77 pct. in the graphite-epoxy laminate. Thus, a significant mode I component of strain energy release rate may be present at the delamination initiation due to an impact load.

Salpekar, S. A.↗
NASA NTRS1990

Thermal expansion of laminated, woven, continuous ceramic fiber/chemical-vapor-infiltrated silicon carbide matrix composites

Thermal expansions of three two-dimensional laminate, continuous fiber/chemical-vapor-infiltrated silicon carbide matrix composites reinforced with either FP-Alumina (alumina), Nextel (mullite), or Nicalon (Si-C-O-N) fibers are reported. Experimental thermal expansion coefficients parallel to a primary fiber orientation were comparable to values calculated by the conventional rule-of-mixtures formula, except for the alumina fiber composite. Hysteresis effects were also observed during repeated thermal cycling of that composite. Those features were attributed to reoccurring fiber/matrix separation related to the micromechanical stresses generated during temperature changes and caused by the large thermal expansion mismatch between the alumina fibers and the silicon carbide matrix.

Eckel, Andrew J.↗
NASA NTRS1978

Predicted inlet gas temperatures for tungsten fiber reinforced superalloy turbine blades

Tungsten fiber reinforced superalloy composite (TFRS) impingement cooled turbine blade inlet gas temperatures were calculated taking into account material spanwise strength, thermal conductivity, material oxidation resistance, fiber-matrix interaction, and coolant flow. Measured values of TFRS thermal conductivities are presented. Calculations indicate that blades made of 30 volume percent fiber content TFRS having a 12,000 N-m/kg stress-to-density ratio while operating at 40 atmospheres and a 0.06 coolant flow ratio could permit a turbine blade inlet gas temperature of over 1900K. This is more than 150K greater than similar superalloy blades.

Winsa, E. A.↗
NASA NTRS1978

Predicted inlet gas temperatures for tungsten fiber reinforced superalloy turbine blades

Tungsten fiber-reinforced superalloy composite (TFRS) impingement-cooled turbine blade inlet gas temperatures were calculated taking into account material spanwise strength, thermal conductivity, material oxidation resistance, fiber-matrix interaction, and coolant flow. Measured values of TFRS thermal conductivities are presented. Calculations indicate that blades made of 30 volume percent fiber content TFRS having a 12,000 N-m/kg stress-to-density ratio while operating at 40 atm and a 0.06 coolant flow ratio could permit a turbine blade inlet gas temperature of over 1900 K. This is more than 150 K greater than similar superalloy blades.

Winsa, E. A.↗
NASA NTRS1992

The optimal fiber volume fraction and fiber-matrix property compatibility in fiber reinforced composites

Although the question of minimum or critical fiber volume fraction beyond which a composite can then be strengthened due to addition of fibers has been dealt with by several investigators for both continuous and short fiber composites, a study of maximum or optimal fiber volume fraction at which the composite reaches its highest strength has not been reported yet. The present analysis has investigated this issue for short fiber case based on the well-known shear lag (the elastic stress transfer) theory as the first step. Using the relationships obtained, the minimum spacing between fibers is determined upon which the maximum fiber volume fraction can be calculated, depending on the fiber packing forms within the composites. The effects on the value of this maximum fiber volume fraction due to such factors as fiber and matrix properties, fiber aspect ratio and fiber packing forms are discussed. Furthermore, combined with the previous analysis on the minimum fiber volume fraction, this maximum fiber volume fraction can be used to examine the property compatibility of fiber and matrix in forming a composite. This is deemed to be useful for composite design. Finally some examples are provided to illustrate the results.

Pan, Ning↗
NASA NTRS1989

A Monte Carlo-finite element model for strain energy controlled microstructural evolution - 'Rafting' in superalloys

This paper presents a specialized microstructural lattice model, MCFET (Monte Carlo finite element technique), which simulates microstructural evolution in materials in which strain energy has an important role in determining morphology. The model is capable of accounting for externally applied stress, surface tension, misfit, elastic inhomogeneity, elastic anisotropy, and arbitrary temperatures. The MCFET analysis was found to compare well with the results of analytical calculations of the equilibrium morphologies of isolated particles in an infinite matrix.

Gayda, J.↗
NASA NTRS1991

Stability analysis of bridged cracks in brittle matrix composites

The bridging of matrix cracks by fibers is an important toughening mechanism in fiber reinforced brittle matrix composites. This paper presents the results of a nonlinear finite element analysis of the Mode-I propagation of a bridged matrix crack in a finite size specimen. The composite is modeled as an orthotropic continuum and the bridging due to the fibers is modeled as a distribution of tractions which resist crack opening. A critical stress intensity factor criterion is employed for matrix crack propagation while a critical crack opening condition is used for fiber failure. The structural response of the specimen (load-deflection curves) as well as the stress intensity factor of the propagating crack are calculated for various constituent properties and specimen configurations for both tensile and bending loading. By controlling the length of the bridged crack results are obtained which highlight the transition from stable to unstable behavior of the propagating crack.

Ballarini, Roberto↗
NASA NTRS1992

Application of neural networks to prediction of advanced composite structures mechanical response and behavior

Two types of neural networks were used to evaluate acousto-ultrasonic (AU) data for material characterization and mechanical reponse prediction. The neural networks included a simple feedforward network (backpropagation) and a radial basis functions network. Comparisons of results in terms of accuracy and training time are given. Acousto-ultrasonic (AU) measurements were performed on a series of tensile specimens composed of eight laminated layers of continuous, SiC fiber reinforced Ti-15-3 matrix. The frequency spectrum was dominated by frequencies of longitudinal wave resonance through the thickness of the specimen at the sending transducer. The magnitude of the frequency spectrum of the AU signal was used for calculating a stress-wave factor based on integrating the spectral distribution function and used for comparison with neural networks results.

Cios, K. J.↗
NASA NTRS1992

The effect of fiber microstructure on evolution of residual stresses in silicon carbide/titanium aluminide composites

This paper examines the effect of the morphology of the SCS6 silicon carbide fiber on the evolution of residual stresses in SiC/Ti composites. A micromechanics model based on the concentric cylinder concept is presented which is used to calculate residual stresses in a SiC/Ti composite during axisymmetric cooling by a spatially uniform temperature change. The silicon carbide fiber is modeled as a layered material with five distinct transversely isotropic and orthotropic, elastic layers, whereas the titanium matrix is taken to be isotropic, with temperature-dependent elastoplastic properties. The results arc compared with those obtained based on the assumption that the silicon carbide fiber is isotropic and homogeneous.

Pindera, Marek-Jerzy↗
NASA NTRS2012

Calculation of the Frequency Distribution of the Energy Deposition in DNA Volumes by Heavy Ions

Radiation quality effects are largely determined by energy deposition in small volumes of characteristic sizes less than 10 nm representative of short-segments of DNA, the DNA nucleosome, or molecules initiating oxidative stress in the nucleus, mitochondria, or extra-cellular matrix. On this scale, qualitatively distinct types of molecular damage are possible for high linear energy transfer (LET) radiation such as heavy ions compared to low LET radiation. Unique types of DNA lesions or oxidative damages are the likely outcome of the energy deposition. The frequency distribution for energy imparted to 1-20 nm targets per unit dose or particle fluence is a useful descriptor and can be evaluated as a function of impact parameter from an ions track. In this work, the simulation of 1-Gy irradiation of a cubic volume of 5 micron by: 1) 450 (1)H(+) ions, 300 MeV; 2) 10 (12)C(6+) ions, 290 MeV/amu and 3) (56)Fe(26+) ions, 1000 MeV/amu was done with the Monte-Carlo simulation code RITRACKS. Cylindrical targets are generated in the irradiated volume, with random orientation. The frequency distribution curves of the energy deposited in the targets is obtained. For small targets (i.e. <25 nm size), the probability of an ion to hit a target is very small; therefore a large number of tracks and targets as well as a large number of histories are necessary to obtain statistically significant results. This simulation is very time-consuming and is difficult to perform by using the original version of RITRACKS. Consequently, the code RITRACKS was adapted to use multiple CPU on a workstation or on a computer cluster. To validate the simulation results, similar calculations were performed using targets with fixed position and orientation, for which experimental data are available [5]. Since the probability of single- and double-strand breaks in DNA as function of energy deposited is well know, the results that were obtained can be used to estimate the yield of DSB, and can be extended to include other targeted or non-target effects.

Plante, Ianik↗
NASA NTRS1988

Thermo-viscoelastic analysis of composite materials, volume 1

Advanced composite materials, especially graphite/epoxy, are being applied to aircraft structures in order to improve performance and save weight. An important consideration in composite design is the residual strength of a structure containing holes, delaminations, or interlaminar damage when subjected to compressive loads. Recent studies have revealed the importance of viscoelastic effects in polymer-based composites. The viscoelastic effect is particularly significant at elevated temperature/moisture conditions since the matrix material is strongly affected by the environment. The solution of viscoelastic problems in composites was limited to special cases which can be solved by classical lamination theory. A finite element procedure is presented for calculating time-dependent stresses and strains in composite structures with general configurations and complicated boundary conditions. Using this procedure the in-plane and interlaminar stress distributions and histories in notched and unnotched composites were obtained for mechanical and thermal loads. Both two-dimensional and three-dimensional viscoelastic problems are analyzed. The effects of layup orientation and load spectrum on creep response and stress relaxation were also studied.

Lin, K. Y.↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Shuttle Leading Edge Materials with Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites, such as the Reinforced Carbon Carbon (RCC) material used on the leading edges of the Space Shuttle. In the developed model, the differences in the tension and compression deformation behaviors have also been accounted for. State variable viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear response is independent of the stiffness in the normal directions. The developed equations have been implemented into LS-DYNA through the use of user defined subroutines (UMATs). Several sample qualitative calculations have been conducted, which demonstrate the ability of the model to qualitatively capture the features of the deformation response present in woven ceramic matrix composites.

Goldberg, Robert K.↗
NASA NTRS1977

Comparison of linear inversion methods by examination of the duality between iterative and inverse matrix methods

Linear numerical inversion methods applied to atmospheric remote sounding generally can be categorized in two ways: (1) iterative, and (2) inverse matrix methods. However, these two categories are not unrelated; a duality exists between them. In other words, given an iterative scheme, a corresponding inverse matrix method exists, and conversely. This duality concept is developed for the more familiar linear methods. The iterative duals are compared with the classical linear iterative approaches and their differences analyzed. The importance of the initial profile in all methods is stressed. Calculations using simulated data are made to compare accuracies and to examine the dependence of the solution on the initial profile.

Fleming, H. E.↗