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 1,099 records · Page 61

Tungsten fiber reinforced copper matrix composites: A review

Tungsten fiber reinforced copper matrix (W/Cu) composites have served as an ideal model system with which to analyze the properties of metal matrix composites. A series of research programs were conducted to investigate the stress-strain behavior of W/Cu composites; the effect of fiber content on the strength, modulus, and conductivity of W/Cu composites; and the effect of alloying elements on the behavior of tungsten wire and of W/Cu composites. Later programs investigated the stress-rupture, creep, and impact behavior of these composites at elevated temperatures. Analysis of the results of these programs as allows prediction of the effects of fiber properties, matrix properties, and fiber content on the properties of W/Cu composites. These analyses form the basis for the rule-of-mixtures prediction of composite properties which was universally adopted as the criteria for measuring composite efficiency. In addition, the analyses allows extrapolation of potential properties of other metal matrix composites and are used to select candidate fibers and matrices for development of tungsten fiber reinforced superalloy composite materials for high temperature aircraft and rocket engine turbine applications. The W/Cu composite efforts are summarized, some of the results obtained are described, and an update is provided on more recent work using W/Cu composites as high strength, high thermal conductivity composite materials for high heat flux, elevated temperature applications.

Mcdanels, David L.↗

Chemical fractionations in meteorites. XI - C2 chondrites

Measurements of the compositions of 20 trace elements in the representative C2 chondrites Boriskino, Cold Bokkeveld, Erakot, Essebi, Haripura, Santa Cruz and Al Rais are reported. The contents of Ag, Au, Bi, Cd, Cs, Ge, In, Ir, Ni, Os, Pd, Rb, Re, Sb, Se, Sn, Te, Tl, U, and Zn were determined by radiochemical neutron activation analysis. The siderophile abundances of the C2 chondrites are found to be less uniform than those of other carbonaceous chondrites, while the C2R chondrite Al Rais is systematically lower in 12 volatiles than the C2M chondrites. Enrichment of Bi and Tl found in Erakot and Haripura indicate the possible presence of the late condensate mysterite. Volatile abundances are shown to agree with matrix contents for meteorites that have suffered little aqueous alteration, however to be 20-30% lower for the more altered meteorites. Finally, the decline of element abundance with volatility is shown to be consistent with the sigmoid curve explained by the two-component model.

Wolf, R.↗

Elemental Spectra from the First ATIC Flight

The Advanced Thin Ionization Calorimeter (ATIC) instrument is a balloon-borne experiment designed to measure the composition and energy spectra of Z = l to 26 cosmic rays over the energy range from approx. 10(exp 11) to approx. 10(exp 14) eV. The instrument consists of a silicon matrix charge detector, plastic scintillator strip hodoscopes interleaved with graphite interaction targets, and a fully active Bismuth Germanate (BGO) calorimeter. ATIC had two successful Long Duration Balloon flights launched from McMurdo Station, Antarctica in 2000 and 2002. In this paper, spectra of various elements measured during the first 16 day flight are presented.

Ahn, H. S.↗

Origins of Solar Systems Program

The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments.

O'D. Alexander, Conel↗

Final Summary of Research Report to the National Aeronautics and Space Administration Origins of Solar Systems Program

The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments. Cosmic spherules present an opportunity to test our evaporation models under flash heating conditions that would be difficult to simulate experimentally. However, there is surprisingly little known about the isotopic compositions of silicate cosmic spherules, and a number of questions need to be addressed. Is the range of compositions they exhibit due to evaporation? If they are, are the relative volatilities consistent with the models/experiments and are the isotopic fractionations consistent with Rayleigh conditions? For instance, do the alkalis and S evaporate prior to significant melting so that conditions did not meet the Rayleigh criteria of rapid diffusion? If so, their isotopic fractionation might be considerably suppressed. Could this mechanism of K loss apply to chondrule formation? The Fe isotopic fractionation during evaporation of silicates has not been measured, so cosmic spherules might provide a clue to whether FeO diffusion is fast enough to maintain Rayleigh conditions during evaporation. And so on.

O'D. Alexander, Conel↗

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

An accurate mathematical model for sound propagation in axisymmetric aircraft engine ducts with compressible mean flow is reported. The model is based on the usual perturbation of the basic fluid mechanics equations for small motions. Mean flow parameters are derived in the absence of fluctuating quantities and are then substituted into the equations for the acoustic quantities which were linearized by eliminating higher order terms. Mean swirl is assumed to be zero from the restriction of axisymmetry. A linear rectangular serendipity element is formulated from these equations using a Galerkin procedure and assembled in a special purpose computer program in which the matrix map for a rectangular mesh was specifically coded. Representations of the fluctuating quantities, mean quantities and coordinate transformations are isoparametric. The global matrix is solved by foreward and back substitution following an L-U decomposition with pivoting restricted internally to the blocks. Results from the model were compared with results from several alternative analyses and yielded satisfactory agreement.

Abrahamson, A. L.↗

The Apollo 17 'melt sheet' - Chemistry, age and Rb/Sr systematics

Major, minor, and trace-element compositions, age data, and Rb/Sr systematics of Apollo 17 boulders have been compiled, and additional analyses performed on a norite breccia clast (77215) included in the Apollo 17, Station 7 boulder. The Apollo 17 boulders are found to be identical or nearly so in major, minor, and trace-element composition, suggesting that they all originated as an impact melt analogous to melt sheets found in larger terrestrial craters. The matrix dates (Ar-40/Ar-39) and Rb/Sr systematics available suggest that this impact melt formed by a single impact about 4 billion years ago. This impact excavated, shocked, brecciated, and melted norites, norite cumulates, and possibly anorthositic gabbros and dunites about 4.4 billion years old. The impact was likely a major one, possibly the Serenitatis basin-forming event.

Winzer, S. R.↗

Accuracy of an approximate static structural analysis technique based on stiffness matrix eigenmodes

Use of the stiffness matrix eigenmodes, instead of the vibration eigenmodes, as generalized coordinates is proposed for condensation of static load deflection equations in finite element stiffness method. The modes are selected by strain energy criteria and the resulting fast, approximate analysis technique is evaluated by applications to idealized built-up wings and a fuselage segment. The best results obtained are a two-order of magnitude reduction of the number of degrees of freedom in a high aspect ratio wing associated with less than one percent error in prediction of the largest displacement.

Sobieszczanski-Sobieski, J.↗

Selecting band combinations with thematic mapper data

A problem arises in making color composite images because there are 210 different possible color presentations of TM three-band images. A method is given for reducing that 210 to a single choice, decided by the statistics of a scene or subscene, and taking into full account any correlations that exist between different bands. Instead of using total variance as the measure for information content of the band triplets, the ellipsoid of maximum volume is selected which discourages selection of bands with high correlation. The band triplet is obtained by computing and ranking in order the determinants of each 3 x 3 principal submatrix of the original matrix M. After selection of the best triplet, the assignment of colors is made by using the actual variances (the diagonal elements of M): green (maximum variance), red (second largest variance), blue (smallest variance).

Sheffield, C. A.↗

Dipole and slot elements and arrays on semi-infinite substrates

The printed dipole or slot antenna on a semi-infinite substrate and infinite phased arrays of these elements are investigated. The solution is based on the moment method in the Fourier transform domain. The generalized impedance or admittance matrix can be expressed in rapidly converging infinite-integral or infinite-summation forms, allowing the accurate determination of the current distributions. Using the present formulation, the input impedance, resonant length, and radiation pattern for the isolated antennas, and the reflection coefficient for infinite phased arrays, are calculated.

Kominami, M.↗

A mathematical model for simulating noise suppression of lined ejectors

A mathematical model containing the essential features embodied in the noise suppression of lined ejectors is presented. Although some simplification of the physics is necessary to render the model mathematically tractable, the current model is the most versatile and technologically advanced at the current time. A system of linearized equations and the boundary conditions governing the sound field are derived starting from the equations of fluid dynamics. A nonreflecting boundary condition is developed. In view of the complex nature of the equations, a parametric study requires the use of numerical techniques and modern computers. A finite element algorithm that solves the differential equations coupled with the boundary condition is then introduced. The numerical method results in a matrix equation with several hundred thousand degrees of freedom that is solved efficiently on a supercomputer. The model is validated by comparing results either with exact solutions or with approximate solutions from other works. In each case, excellent correlations are obtained. The usefulness of the model as an optimization tool and the importance of variable impedance liners as a mechanism for achieving broadband suppression within a lined ejector are demonstrated.

Watson, Willie R.↗

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗

Finite Element Analysis of Geodesically Stiffened Cylindrical Composite Shells Using a Layerwise Theory

Layerwise finite element analyses of geodesically stiffened cylindrical shells are presented. The layerwise laminate theory of Reddy (LWTR) is developed and adapted to circular cylindrical shells. The Ritz variational method is used to develop an analytical approach for studying the buckling of simply supported geodesically stiffened shells with discrete stiffeners. This method utilizes a Lagrange multiplier technique to attach the stiffeners to the shell. The development of the layerwise shells couples a one-dimensional finite element through the thickness with a Navier solution that satisfies the boundary conditions. The buckling results from the Ritz discrete analytical method are compared with smeared buckling results and with NASA Testbed finite element results. The development of layerwise shell and beam finite elements is presented and these elements are used to perform the displacement field, stress, and first-ply failure analyses. The layerwise shell elements are used to model the shell skin and the layerwise beam elements are used to model the stiffeners. This arrangement allows the beam stiffeners to be assembled directly into the global stiffness matrix. A series of analytical studies are made to compare the response of geodesically stiffened shells as a function of loading, shell geometry, shell radii, shell laminate thickness, stiffener height, and geometric nonlinearity. Comparisons of the structural response of geodesically stiffened shells, axial and ring stiffened shells, and unstiffened shells are provided. In addition, interlaminar stress results near the stiffener intersection are presented. First-ply failure analyses for geodesically stiffened shells utilizing the Tsai-Wu failure criterion are presented for a few selected cases.

Gerhard, Craig Steven↗

Progressive Damage Analysis of Bonded Composite Joints

The present work is related to the development and application of progressive damage modeling techniques to bonded joint technology. The joint designs studied in this work include a conventional composite splice joint and a NASA-patented durable redundant joint. Both designs involve honeycomb sandwich structures with carbon/epoxy facesheets joined using adhesively bonded doublers.Progressive damage modeling allows for the prediction of the initiation and evolution of damage within a structure. For structures that include multiple material systems, such as the joint designs under consideration, the number of potential failure mechanisms that must be accounted for drastically increases the complexity of the analyses. Potential failure mechanisms include fiber fracture, intraply matrix cracking, delamination, core crushing, adhesive failure, and their interactions. The bonded joints were modeled using highly parametric, explicitly solved finite element models, with damage modeling implemented via custom user-written subroutines. Each ply was discretely meshed using three-dimensional solid elements. Layers of cohesive elements were included between each ply to account for the possibility of delaminations and were used to model the adhesive layers forming the joint. Good correlation with experimental results was achieved both in terms of load-displacement history and the predicted failure mechanism(s).

Leone, Frank A., 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.↗

A Three-Dimensional Thermodynamically Based Function for the Progressive Failure of Unidirectional Composites

A thermodynamically-based work potential theory for modeling progressive damage for laminated, unidirectional composites assuming plane stress (2D Schaperystheory) is extended to three dimensioanl (3D). An internal state variable, S, is defined to account for the dissipated energy due to damage evolution in the form of microstructure changes in the matrix. With the stationary of the total work potential with respect to the internal state variable, a thermodynamically-consistent set of evolution equationsis derived. The internal state variable is related to the transverse and shear modulithrough microdamage functions. In the first part of this work, coupon specimensare prepared to conduct experiments to characterize the relations between the internalstate variable and the transverse modulus as well as shear modulus. The informationis subsequent used for the prediction of three point bending test. In the second partof this work, objectivity is studied. Three separate methods utilizing different definitions of a reduced internal state variable or of the order of the polynomials are used to represent the matrix micro-damage functions are employed. The three methods are implemented in a user defined subroutine within a commercial finite element method software package. Results from numerical simulations of a center-notched composites panel are compared. The agreement in the maximum stress predictions among the three methods indicates that objectivity, with respect to the functional form of themicrodamage functions, is satisfied.

carbon fiber reinforced polymer composites↗

An Improved Plastically Dilatant Unified Viscoplastic Constitutive Formulation for Multiscale Analysis of Polymer Matrix Composites Under High Strain Rate Loading

Polymer matrix composites are commonly used to fabricate energy-absorbing structures expected to experience impact loading. As such, a detailed understanding of the dynamic response of the constituent materials is necessary. Since the rate, temperature, and pressure dependence of carbon fiber reinforced polymer matrix composites are primarily manifestations of the rate, temperature, and pressure dependence of the polymer matrix, it is crucial that the constitutive behavior of the matrix be accurately characterized. In this work, an existing unified viscoplastic constitutive formulation is extended to ensure thermodynamic consistency and to more accurately account for the tension-compression asymmetry observed in the response of polymeric materials. A new plastic potential function is proposed, and elementary loading conditions are utilized to determine relations between model constants to ensure nonnegative plastic dissipation, a necessary thermodynamic requirement. Expressions for plastic Poisson’s ratios are derived and are bounded by enforcing nonnegative plastic dissipation. The model is calibrated against available experimental data from tests conducted over a range of strain rates, temperatures, and loading cases on a representative thermoset epoxy; good correlation between simulations and experimental data is obtained. Temperature rises due to the conversion of plastic work to heat are computed via the adiabatic heat energy equation. The viscoplastic polymer model is then used as a constitutive model in the generalized method of cells micromechanics theory to investigate the effects of matrix adiabatic heating on the high strain rate response of a unidirectional composite. The thermodynamic consistency of the model ensures plastic dissipation can only cause an increase in temperature. Simulation results indicate that significant thermal softening due to the conversion of plastic work to heat is observed in the composite for matrix dominated deformation modes.

Impact loading↗

A viscoplastic constitutive theory for metal matrix composites at high temperature

A viscoplastic constitutive theory is presented for representing the high-temperature deformation behavior of metal matrix composites. The point of view taken is a continuum one where the composite is considered a material in its own right, with its own properties that can be determined for the composite as a whole. It is assumed that a single preferential (fiber) direction is identifiable at each material point (continuum element) admitting the idealization of local transverse isotropy. A key ingredient in this work is the specification of an experimental program for the complete determination of the material functions and parameters for characterizing a particular metal matrix composite. The parameters relating to the strength of anisotropy can be determined through tension/torsion tests on longitudinally and circumferentially reinforced thin-walled tubes. Fundamental aspects of the theory are explored through a geometric interpretation of some basic features analogous to those of the classical theory of plasticity.

Robinson, D. N.↗