Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

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 433 records · Page 24

Prediction of thermal and mechanical stress-strain responses of TMC's subjected to complex TMF histories

This paper presents an experimental and analytical evaluation of cross-plied laminates of Ti-15V-3Cr-3Al-3Sn (Ti-15-3) matrix reinforced with continuous silicon-carbide fibers (SCS-6) subjected to a complex TMF loading profile. Thermomechanical fatigue test techniques were developed to conduct a simulation of a generic hypersonic flight profile. A micromechanical analysis was used. The analysis predicts the stress-strain response of the laminate and of the constituents in each ply during thermal and mechanical cycling by using only constituent properties as input. The fiber was modeled as elastic with transverse orthotropic and temperature-dependent properties. The matrix was modeled using a thermoviscoplastic constitutive relation. The fiber transverse modulus was reduced in the analysis to simulate the fiber-matrix interface failures. Excellent correlation was found between measured and predicted laminate stress-strain response due to generic hypersonic flight profile when fiber debonding was modeled.

Johnson, W. S.↗

Controller design via structural reduced modeling by FETM

The Finite Element - Transfer Matrix (FETM) method has been developed to reduce the computations involved in analysis of structures. This widely accepted method, however, has certain limitations, and does not directly produce reduced models for control design. To overcome these shortcomings, a modification of FETM method has been developed. The modified FETM method easily produces reduced models that are tailored toward subsequent control design. Other features of this method are its ability to: (1) extract open loop frequencies and mode shapes with less computations, (2) overcome limitations of the original FETM method, and (3) simplify the procedures for output feedback, constrained compensation, and decentralized control. This semi annual report presents the development of the modified FETM, and through an example, illustrates its applicability to an output feedback and a decentralized control design.

Yousuff, A.↗

Band Anticrossing in Highly Mismatched Compound Semiconductor Alloys

Compound semiconductor alloys in which metallic anions are partially replaced with more electronegative isoelectronic atoms have recently attracted significant attention. Group IIIN(x)V(1-x), alloys with a small amount of the electronegative N substituting more metallic column V elements has been the most extensively studied class of such Highly Mismatched Alloys (HMAs). We have shown that many of the unusual properties of the IIIN(x),V(1-x) alloys can be well explained by the Band Anticrossing (BAC) model that describes the electronic structure in terms of an interaction between highly localized levels of substitutional N and the extended states of the host semiconductor matrix. Most recently the BAC model has been also used to explain similar modifications of the electronic band structure observed in Te-rich ZnS(x)Te(l-x) and ZnSe(Y)Te(1-y) alloys. To date studies of HMAs have been limited to materials with relatively small concentrations of highly electronegative atoms. Here we report investigations of the electronic structure of ZnSe(y)Te(1-y) alloys in the entire composition range, 0 less than or equal to y less than or equal to 1. The samples used in this study are bulk ZnSe(y)Te(1-y) crystals grown by either a modified Bridgman method or by physical vapor transport. Photomodulated reflection (PR) spectroscopy was used to measure the composition dependence of optical transitions from the valence band edge and from the spin-orbit split off band to the conduction band. The pressure dependence of the band gap was measured using optical absorption in a diamond anvil cell. We find that the energy of the spin-orbit split off valence band edge does not depend on composition and is located at about 3 eV below the conduction band edge of ZnSe. On the Te-rich side the pressure and the composition dependence of the optical transitions are well explained by the BAC model which describes the downward shift of the conduction band edge in terms of the interaction between localized Se states and the conduction band. On the other hand we show that the large band gap reduction observed on the Se-rich side of the alloy system is a result of an interaction between the localized Te level and the valence bands. This interaction leads to the formation of a Te-like valence band edge that strongly interacts with the light hole valence band. Calculations based on a modified k p model account for the reduction of the band gap and the large increase of the spin-orbit splitting observed in Se-rich ZnSe(y)Te(l-y) alloys. We will also discuss the importance of these new results for understanding of the electronic structure and band offsets in other highly mismatched alloy systems.

Yu, Kin Man↗

Band Anticrossing in Highly Mismatched Compound Semiconductor Alloys

Compound semiconductor alloys in which metallic anions are partially replaced with more electronegative isoelectronic atoms have recently attracted significant attention. Group IIIN(sub x)V(sub 1-x) alloys with a small amount of the electronegative N substituting more metallic column V elements has been the most extensively studied class of such Highly Mismatched Alloys (HMAs). We have shown that many of the unusual properties of the IIIN(sub x)V(sub 1-x) alloys can be well explained by the Band Anticrossing (BAC) model that describes the electronic structure in terms of an interaction between highly localized levels of substitutional N and the extended states of the host semiconductor matrix. Most recently the BAC model has been also used to explain similar modifications of the electronic band structure observed in Te-rich ZnS(sub x)Te(sub 1-x) and ZnSe(sub y)Te(sub 1-y) alloys. To date studies of HMAs have been limited to materials with relatively small concentrations of highly electronegative atoms. Here we report investigations of the electronic structure of ZnSe(sub y)Te(sub 1-y) alloys in the entire composition range, y between 0 and 1. The samples used in this study are bulk ZnSe(sub y)Te(sub 1-y) crystals grown by either a modified Bridgman method or by physical vapor transport. Photomodulated reflection (PR) spectroscopy was used to measure the composition dependence of optical transitions from the valence band edge and from the spin-orbit split off band to the conduction band. The pressure dependence of the band gap was measured using optical absorption in a diamond anvil cell. We find that the energy of the spin-orbit split off valence band edge does not depend on composition and is located at about 3 eV below the conduction band edge of ZnSe. On the Te-rich side the pressure and the composition dependence of the optical transitions are well explained by the BAC model which describes the downward shift of the conduction band edge in terms of the interaction between localized Se states and the conduction band. On the other hand we show that the large band gap reduction observed on the Se-rich side of the alloy system is a result of an interaction between the localized Te level and the valence bands. This interaction leads to the formation of a Te-like valence band edge that strongly interacts with the light hole valence band. Calculations based on a modified k(sup dot)p model account for the reduction of the band gap and the large increase of the spin-orbit splitting observed in Se-rich ZnSe(sub y)Te(sub 1-y) alloys. We will also discuss the importance of these new results for understanding of the electronic structure and band offsets in other highly mismatched alloy systems.

Yu, Kin Man↗

Modeling stiffness loss in boron/aluminum laminates 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. Previously announced in STAR as N82-24298

Johnson, W. S.↗

Low-order design and high-order simulation of active closed-loop control for aerospace structures under construction

Partially constructed/assembled structures in space are complicated enough but their dynamics will also be operating in closed-loop with feedback controllers. The dynamics of such structures are modeled by large-scale finite element models. The model dimension L is extremely large (approximately 10,000) while the numbers of actuators (M) and sensors (P) are small. The model parameters M(sub m) mass matrix, D(sub o) damping matrix, and K(sub o) stiffness matrix, are all symmetric and sparse (banded). Thus simulation of open-loop structure models of very large dimension can be accomplished by special integration techniques for sparse matrices. The problem of simulation of closed-loop control of such structures is complicated by the addition of controllers. Simulation of closed-loop controlled structures is an essential part of the controller design and evaluation process. Current research in the following areas is presented: high-order simulation of actively controlled aerospace structures; low-order controller design and SCI compensation for unmodeled dynamics; prediction of closed-loop stability using asymptotic eigenvalue series; and flexible robot manipulator control experiment.

Balas, Mark J.↗

Oxidation Kinetics and Strength Degradation of Carbon Fibers in a Cracked Ceramic Matrix Composite

Experimental results and oxidation modeling will be presented to discuss carbon fiber susceptibility to oxidation, the oxidation kinetics regimes and composite strength degradation and failure due to oxidation. Thermogravimetric Analysis (TGA) was used to study the oxidation rates of carbon fiber and of a pyro-carbon interphase. The analysis was used to separately obtain activation energies for the carbon constituents within a C/SiC composite. TGA was also conducted on C/SiC composite material to study carbon oxidation and crack closure as a function of temperature. In order to more closely match applications conditions C/SiC tensile coupons were also tested under stressed oxidation conditions. The stressed oxidation tests show that C/SiC is much more susceptible to oxidation when the material is under an applied load where the cracks are open and allow for oxygen ingress. The results help correlate carbon oxidation with composite strength reduction and failure.

Halbig, Michael C.↗

Effect of microstructural damage on ply stresses in laminated composites

The mechanisms involved in damage and failure of laminated orthotropic composites are investigated theoretically. The continuum model developed accounts for both matrix cracks and interply delamination using second-order tensor-valued internal-state variables based on the locally averaged microcrack dynamics. The derivation of the model is given in detail, and numerical results for sample problems are presented in extensive graphs and tables. The model is shown to be effective in predicting stresses at the ply level, and significant damage-induced decreases in laminate stress states are found.

Allen, D. H.↗

A methodology to predict damage initiation, damage growth and residual strength in titanium matrix composites

In this research, a methodology to predict damage initiation, damage growth, fatigue life, and residual strength in titanium matrix composites (TMC) is outlined. Emphasis was placed on micromechanics-based engineering approaches. Damage initiation was predicted using a local effective strain approach. A finite element analysis verified the prevailing assumptions made in the formulation of this model. Damage growth, namely, fiber-bridged matrix crack growth, was evaluated using a fiber bridging (FB) model which accounts for thermal residual stresses. This model combines continuum fracture mechanics and micromechanics analyses yielding stress-intensity factor solutions for fiber-bridged matrix cracks. It is assumed in the FB model that fibers in the wake of the matrix crack are idealized as a closure pressure, and an unknown constant frictional shear stress is assumed to act along the debond length of the bridging fibers. This frictional shear stress was used as a curve fitting parameter to the available experimental data. Fatigue life and post-fatigue residual strength were predicted based on the axial stress in the first intact 0 degree fiber calculated using the FB model and a three-dimensional finite element analysis.

Bakuckas, J. G., Jr.↗

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗

Outgassing and dimensional changes of polymer matrix composites in space

A thermal-vacuum outgassing model and test protocol for predicting outgassing times and dimensional changes for polymer matrix composites is described. Experimental results derived from a 'control' sample are used to provide the basis for analytical predictions to compare with the outgassing response of Long Duration Exposure Facility (LDEF) flight samples.

Tennyson, R. C.↗

Experimental Space Shuttle Orbiter Studies to Acquire Data for Code and Flight Heating Model Validation

In an experimental study to obtain detailed heating data over the Space Shuttle Orbiter, CUBRC has completed an extensive matrix of experiments using three distinct models and two unique hypervelocity wind tunnel facilities. This detailed data will be employed to assess heating augmentation due to boundary layer transition on the Orbiter wing leading edge and wind side acreage with comparisons to computational methods and flight data obtained during the Orbiter Entry Boundary Layer Flight Experiment and HYTHIRM during STS-119 reentry. These comparisons will facilitate critical updates to be made to the engineering tools employed to make assessments about natural and tripped boundary layer transition during Orbiter reentry. To achieve the goals of this study data was obtained over a range of Mach numbers from 10 to 18, with flight scaled Reynolds numbers and model attitudes representing key points on the Orbiter reentry trajectory. The first of these studies were performed as an integral part of Return to Flight activities following the accident that occurred during the reentry of the Space Shuttle Columbia (STS-107) in February of 2003. This accident was caused by debris, which originated from the foam covering the external tank bipod fitting ramps, striking and damaging critical wing leading edge heating tiles that reside in the Orbiter bow shock/wing interaction region. During investigation of the accident aeroheating team members discovered that only a limited amount of experimental wing leading edge data existed in this critical peak heating area and a need arose to acquire a detailed dataset of heating in this region. This new dataset was acquired in three phases consisting of a risk mitigation phase employing a 1.8% scale Orbiter model with special temperature sensitive paint covering the wing leading edge, a 0.9% scale Orbiter model with high resolution thin-film instrumentation in the span direction, and the primary 1.8% scale Orbiter model with detailed thin-film resolution in both the span and chord direction in the area of peak heating. Additional objectives of this first study included: obtaining natural or tripped turbulent wing leading edge heating levels, assessing the effectiveness of protuberances and cavities placed at specified locations on the orbiter over a range of Mach numbers and Reynolds numbers to evaluate and compare to existing engineering and computational tools, obtaining cavity floor heating to aid in the verification of cavity heating correlations, acquiring control surface deflection heating data on both the main body flap and elevons, and obtain high speed schlieren videos of the interaction of the orbiter nose bow shock with the wing leading edge. To support these objectives, the stainless steel 1.8% scale orbiter model in addition to the sensors on the wing leading edge was instrumented down the windward centerline, over the wing acreage on the port side, and painted with temperature sensitive paint on the starboard side wing acreage. In all, the stainless steel 1.8% scale Orbiter model was instrumented with over three-hundred highly sensitive thin-film heating sensors, two-hundred of which were located in the wing leading edge shock interaction region. Further experimental studies will also be performed following the successful acquisition of flight data during the Orbiter Entry Boundary Layer Flight Experiment and HYTHIRM on STS-119 at specific data points simulating flight conditions and geometries. Additional instrumentation and a protuberance matching the layout present during the STS-119 boundary layer transition flight experiment were added with testing performed at Mach number and Reynolds number conditions simulating conditions experienced in flight. In addition to the experimental studies, CUBRC also performed a large amount of CFD analysis to confirm and validate not only the tunnel freestream conditions, but also 3D flows over the orbiter acreage, wing leading edge, and controlurfaces to assess data quality, shock interaction locations, and control surface separation regions. This analysis is a standard part of any experimental program at CUBRC, and this information was of key importance for post-test data quality analysis and understanding particular phenomena seen in the data. All work during this effort was sponsored and paid for by the NASA Space Shuttle Program Office at the Johnson Space Center in Houston, Texas.

Wadhams, T. P.↗

Modeling the Elastic Modulus of 2D Woven CVI SiC Composites

The use of fiber, interphase, CVI SiC minicomposites as structural elements for 2D-woven SiC fiber reinforced chemically vapor infiltrated (CVI) SiC matrix composites is demonstrated to be a viable approach to model the elastic modulus of these composite systems when tensile loaded in an orthogonal direction. The 0deg (loading direction) and 90deg (perpendicular to loading direction) oriented minicomposites as well as the open porosity and excess SiC associated with CVI SiC composites were all modeled as parallel elements using simple Rule of Mixtures techniques. Excellent agreement for a variety of 2D woven Hi-Nicalon(TradeMark) fiber-reinforced and Sylramic-iBN reinforced CVI SiC matrix composites that differed in numbers of plies, constituent content, thickness, density, and number of woven tows in either direction (i.e, balanced weaves versus unbalanced weaves) was achieved. It was found that elastic modulus was not only dependent on constituent content, but also the degree to which 90deg minicomposites carried load. This depended on the degree of interaction between 90deg and 0deg minicomposites which was quantified to some extent by composite density. The relationships developed here for elastic modulus only necessitated the knowledge of the fractional contents of fiber, interphase and CVI SiC as well as the tow size and shape. It was concluded that such relationships are fairly robust for orthogonally loaded 2D woven CVI SiC composite system and can be implemented by ceramic matrix composite component modelers and designers for modeling the local stiffness in simple or complex parts fabricated with variable constituent contents.

Morscher, Gregory N.↗

A scatter model for vegetation up to Ku-band

A scatter model is developed based on the matrix doubling method for volume scattering and the Kirchhoff method in rough surface scattering. Scattering from vegetation is assumed to be dominated by leaves and a single leaf is modeled by a thin dielectric disk. In developing the phase matrix for the disk, field within the disk is taken to be constant over the disk thickness, but phase changes across the surface of the disk are accounted for. Comparisons of this scatter model with radar measurements indicate good agreements in polarization, angular trends, and frequency up to Ku-band. This represents a considerable improvement over low-frequency scatter models which are valid up to S-band.

Eom, H. J.↗

Associative Flow Rule Used to Include Hydrostatic Stress Effects in Analysis of Strain-Rate-Dependent Deformation of Polymer Matrix Composites

designing reliable composite engine cases that are lighter than the metal cases in current use. The types of polymer matrix composites that are likely to be used in such an application have a deformation response that is nonlinear and that varies with strain rate. The nonlinearity and the strain-rate dependence of the composite response are due primarily to the matrix constituent. Therefore, in developing material models to be used in the design of impact-resistant composite engine cases, the deformation of the polymer matrix must be correctly analyzed. However, unlike in metals, the nonlinear response of polymers depends on the hydrostatic stresses, which must be accounted for within an analytical model. By applying micromechanics techniques along with given fiber properties, one can also determine the effects of the hydrostatic stresses in the polymer on the overall composite deformation response. First efforts to account for the hydrostatic stress effects in the composite deformation applied purely empirical methods that relied on composite-level data. In later efforts, to allow polymer properties to be characterized solely on the basis of polymer data, researchers at the NASA Glenn Research Center developed equations to model the polymers that were based on a non-associative flow rule, and efforts to use these equations to simulate the deformation of representative polymer materials were reasonably successful. However, these equations were found to have difficulty in correctly analyzing the multiaxial stress states found in the polymer matrix constituent of a composite material. To correct these difficulties, and to allow for the accurate simulation of the nonlinear strain-rate-dependent deformation analysis of polymer matrix composites, in the efforts reported here Glenn researchers reformulated the polymer constitutive equations from basic principles using the concept of an associative flow rule. These revised equations were characterized and validated in an experimental program carried out through a university grant with the Ohio State University, wherein tensile and shear deformation data were obtained for a representative polymer for strain rates ranging from quasi-static to high rates of several hundred per second. Tensile deformation data also were obtained over a variety of strain rates and fiber orientation angles for a representative polymer matrix composite composed using the polymer.

Goldberg, Robert K.↗

A Study of Influencing Factors on the Tensile Response of a Titanium Matrix Composite With Weak Interfacial Bonding

The generalized method of cells micromechanics model is utilized to analyze the tensile stress-strain response of a representative titanium matrix composite with weak interfacial bonding. The fiber/matrix interface is modeled through application of a displacement discontinuity between the fiber and matrix once a critical debonding stress has been exceeded. Unidirectional composites with loading parallel and perpendicular to the fibers are examined, as well as a cross-ply laminate. For each of the laminates studied, analytically obtained results are compared to experimental data. The application of residual stresses through a cool-down process was found to have a significant effect on the tensile response. For the unidirectional laminate with loading applied perpendicular to the fibers, fiber packing and fiber shape were shown to have a significant effect on the predicted tensile response. Furthermore, the interface was characterized through the use of semi-emperical parameters including an interfacial compliance and a "debond stress;" defined as the stress level across the interface which activates fiber/matrix debonding. The results in this paper demonstrate that if architectural factors are correctly accounted for and the interface is appropriately characterized, the macro-level composite behavior can be correctly predicted without modifying any of the fiber or matrix constituent properties.

Goldberg, Robert K.↗

Spatial orientation of the vestibular system

1. A simplified three-dimensional state space model of visual vestibular interaction was formulated. Matrix and dynamical system operators representing coupling from the semicircular canals and the visual system to the velocity storage integrator were incorporated into the model. 2. It was postulated that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix while another rotates the pitch and roll eigenvectors with the head, while maintaining the yaw axis eigenvector approximately spatially invariant. Using this representation, the response characteristics of the pitch, roll, and yaw eye velocity were obtained in terms of the eigenvalues and associated eigenvectors. 3. Using OKAN data obtained from monkeys and comparing to the model predictions, the eigenvalues and eigenvectors of the system matrix were identified as a function of tilt to the side or of tilt to the prone positions, using a modification of the Marquardt algorithm. The yaw eigenvector for right-side-down tilt and for downward pitch cross-coupling was approximately 30 degrees from the spatial vertical. For the prone position, the eigenvector was computed to be approximately 20 degrees relative to the spatial vertical. For both side-down and prone positions, oblique OKN induced along eigenvector directions generated OKAN which decayed to zero along a straight line with approximately a single time constant. This was verified by a spectral analysis of the residual sequence about the straight line fit to the decaying data. The residual sequence was associated with a narrow autocorrelation function and a wide power spectrum. 4. Parameters found using the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced to represent the direct pathway and the coupling of the visual system to the integrator. Model simulations predicted the behavior of yaw and pitch OKN and OKAN when the animal was upright, as well as the cross-coupling in the tilted position. The trajectories in velocity space were also accurately simulated. 5. There were similarities between the monkey eigenvectors and human perception of the spatial vertical. For side-down tilts and downward eye velocity cross-coupling, there was only an Aubert (A) effect. For upward eye velocity cross-coupling there were both Muller (E) and Aubert (A) effects. The mean of the eigenvectors for upward and downward eye velocities overlay human 1 x g perceptual data.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Neuroscience↗

Wave normal directions of chorus near the equatorial source region

OGO 5 triaxial search coil magnetometer data are used to determine the wave propagation directions of postmidnight chorus in the near-equatorial region at L shells of 6 to 7. The methods used to estimate the wave normal directions involved minimum variance, the imaginary part of the cross-spectral matrix, the eigenvector of the Hermitian cross-spectral matrix, and a fitting of dispersion relations for one-wave and two-wave models to the cross-spectral matrix. Wave propagation at all frequencies within chorus tones is found to occur most frequently along the magnetic field with median and average cone angles of 9.1 deg and 12.2 deg, respectively. It is concluded that for waves propagating parallel to B, wave growth is maximum.

Goldstein, B. E.↗