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 1,081 records · Page 60

Multi-Scale Modeling of Cross-Linked Nanotube Materials

The effect of cross-linking single-walled carbon nanotubes on the Young's modulus of a nanotube-reinforced composite is modeled with a multi-scale method. The Young's modulus is predicted as a function of nanotube volume fraction and cross-link density. In this method, the constitutive properties of molecular representative volume elements are determined using molecular dynamics simulation and equivalent-continuum modeling. The Young's modulus is subsequently calculated for cross-linked nanotubes in a matrix which consists of the unreacted cross-linking agent. Two different cross-linking agents are used in this study, one that is short and rigid (Molecule A), and one that is long and flexible (Molecule B). Direct comparisons between the predicted elastic constants are made for the models in which the nanotubes are either covalently bonded or not chemically bonded to the cross-linking agent. At a nanotube volume fraction of 10%, the Young's modulus of Material A is not affected by nanotube crosslinking, while the Young's modulus of Material B is reduced by 64% when the nanotubes are cross-linked relative to the non-cross-linked material with the same matrix.

Frankland, S. J. V.↗

Incorporation of a progressive failure analysis method in the CSM testbed software system

Analysis of the postbuckling behavior of composite shell structures pose many difficult and challenging problems in the field of structural mechanics. Current analysis methods perform well for most cases in predicting the postbuckling response of undamaged components. To predict component behavior accurately at higher load levels, the analysis must include the effects of local material failures. The CSM testbed software system is a highly modular structural analysis system currently under development at Langley Research Center. One of the primary goals of the CSM testbed is to provide a software environment for the development of advanced structural analysis methods and modern numerical methods which will exploit advanced computer architecture such as parallel-vector processors. Development of a progressive failure analysis method consists of the design and implementation of a processor which will perform the ply-level progressive failure analysis and the development of a geometrically nonlinear analysis procedure which incorporates the progressive failure processor. Regarding the development of the progressive failure processor, two components are required: failure criteria and a degradation model. For the initial implementation, the failure criteria of Hashin will be used. For a matrix failure which typically indicates the development of transverse matrix cracks, the ply properties will be degraded. Work to date includes the design of the progressive failure analysis processor and initial plans for the controlling geometrically nonlinear analysis procedure. The implementation of the progressive failure analysis has begun. Access to the model database and the Hashin failure criteria are completed. Work is in progress on the input/output operations for the processor related data and the finite element model updating procedures. In total the progressive failure processor is approximately one-third complete.

Arenburg, Robert T.↗

Free-edge stress analysis of glass-epoxy laminates with matrix cracks

The effect of matrix cracks on the composite delamination and interlaminar stresses is investigated in (+15/90n/-15)s glass-epoxy laminates (with values of n equal to 0, 1, 2, or 3) subjected to monotonically increasing tension loads. Three-dimensional (3D) and quasi-3D (Q3D) finite-element analyses are used to model the free-edge stress states in the laminates with and without a matrix crack, respectively. The Q3D results show that in-plane transverse tensile stresses exist in the +15 deg plies near the free edges of all of the laminates used and that only the interlaminar shear stress is high at the +15/theta interface. The results of 3D analysis indicate that large tensile interlaminar normal as well as shear stresses develop at the intersection of the matrix crack and the free edge. This suggests that the interlaminar normal stress plays a significant role in the failure of these laminates.

Fish, John C.↗

Introduction to the analysis of delamination related fracture processes in composites

This research concerns the analysis and prediction of delamination damage that occurs in composite structures on the sublaminate scale - that is, the scale of individual plies or groups of plies. The objective was to develop analytical models for fixed-mode delamination in composites. These include: (1) the influence of residual thermal and moisture strains; (2) local or transverse crack tip delamination originating at the tip of transverse matrix cracks; and (3) delamination in tapered composite under tensile loading. Computer codes based on the analytical models were developed and comparisons of predictions with available experimental and analytical results in the literature were performed.

Armanios, Erian A.↗

The reduced order model problem in distributed parameter systems adaptive identification and control

The basic assumption that a large space structure can be decoupled preceding the application of reduced order active control was considered and alternative solutions to the control of such structures (in contrast to the strict modal control) were investigated. The transfer function matrix from the actuators to the sensors was deemed to be a reasonable candidate. More refined models from multivariable systems theory were studied and recent results in the multivariable control field were compared with respect to theoretical deficiencies and likely problems in application to large space structures.

Johnson, C. R., Jr.↗

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

Effects of noise variance model on optimal feedback design and actuator placement

In optimal placement of actuators for stochastic systems, it is commonly assumed that the actuator noise variances are not related to the feedback matrix and the actuator locations. In this paper, we will discuss the limitation of that assumption and develop a more practical noise variance model. Various properties associated with optimal actuator placement under the assumption of this noise variance model are discovered through the analytical study of a second order system.

Ruan, Mifang↗

Vertical Transport Processes for Inert and Scavenged Species: TRACE-A Measurements

The TRACE-A mission of the NASA DC-8 aircraft made a large-scale survey of the tropical and subtropical atmosphere in September and October of 1992. Both In-situ measurements of CO (G. Sachsen NASA Langley) and aerosol size (J. Browell group, NASA Langley) provide excellent data sets with which to constrain vertical transport by planetary boundary layer mixing and deep-cloud cumulus convection. Lidar profiles of aerosol-induced scattering and ozone (also by Bremen) are somewhat require more subtle interpretation as tracers, but the vertical information on layering largely compensates for these complexities. The reason this DC-8 dataset is so useful is that very large areas of biomass burning over Africa and South America provide surface sources of appropriate sizes with which to characterize vertical and horizontal motions; the major limitation of our source description is that biomass burning patterns move considerably every few days, and daily burning inventories are a matter of concurrent, intensive research. We use the Penn State / NCAR MM5 model in an assimilation mode on the synoptic and intercontinental scale, and assess the success it shows in vertical transport descriptions. We find that the general level of emissions suggested by the climatological approach (Will. Has, U. of Montana) appears to be approximately correct, possibly a bit low, for this October, 1992, time period. Vertical transport in planetary boundary layer mixing to 5.5 kin was observed and reproduced in our simulations. Furthermore we find evidence that Blackader "transilient" or matrix-transport scheme is needed, but may require some adaptation in our tracer model: CO seems to exhibit very high values at the top of the planetary boundary layer, a process that stretches the eddy-diffusion parameterization. We will report on progress in improving the deep convective transport of carbon monoxide: the Grail scheme as we used it at 100 kin resolution did not transport enough material to the upper troposphere. We expect to be able to attribute this to either parameterization reasons (inadequacy of this parameterization at the large 100km scale) or other reasons. Nevertheless, the qualitative nature of deep transport by clouds shows up well in the simulations. As for scavengable species, the simulations predict tens of micrograms per standard cubic meter of smoke aerosol in the boundary layer. In a straightforward illustration of our simple bulk-mass scavenging parameterization, to one or two micrograms per standard cubic meter of smoke aerosol in the free troposphere just above the source regions: very high concentrations for the free troposphere. We expect to report on comparisons of these predictions to a variety of observations.

Chatfield, Robert B.↗

Advanced composites for large Navy spacecraft

An overview is given of work conducted on contract for the Naval Sea Systems Command. The objective of this contract was to provide direction for the development of high modulus graphite reinforced metal matrix composites. These advanced materials can have a significant effect on the performance of a spacecraft before, during and after an evasive maneuver. The work conducted on this program was organized into seven technical tasks. Task 1 was development of a generic Navy spacecraft model. Finite element models of candidate structural designs were developed. In Task 2, the finite-element model(s) of the structure were used to conduct analytical assessments involving conventional materials, resin matrix composites and metal matrix composites (MMC). In Task 3 and 4, MMC material design, fabrication and evaluation was conducted. This consisted of generating material designs and developing a data base for a broad range of graphite reinforced MMC materials. All material was procured according to specifications which set material quality and material property standards. In Task 5, a set of evasive maneuvering requirements were derived and used in Task 6 to conduct analytical simulations. These analytical simulations used current SOA material properties and projected material properties to provide an indication of key payoffs for material development. In Task 7, a set of material development recommendations was generated.

Davis, William E.↗

The Effect of Fiber Architecture on Matrix Cracking in Sic/sic Cmc's

Applications incorporating silicon carbide fiber reinforced silicon carbide matrix composites (CMC's) will require a wide range of fiber architectures in order to fabricate complex shape. The stress-strain response of a given SiC/SiC system for different architectures and orientations will be required in order to design and effectively life-model future components. The mechanism for non-linear stress-strain behavior in CMC's is the formation and propagation of bridged-matrix cracks throughout the composite. A considerable amount of understanding has been achieved for the stress-dependent matrix cracking behavior of SiC fiber reinforced SiC matrix systems containing melt-infiltrated Si. This presentation will outline the effect of 2D and 3D architectures and orientation on stress-dependent matrix-cracking and how this information can be used to model material behavior and serve as the starting point foe mechanistic-based life-models.

Morscher, Gregory N.↗

A fractal model for crustal deformation

It is hypothesized that crustal deformation occurs on a scale-invariant matrix of faults. For simplicity, a two-dimensional pattern of hexagons on which strike-slip faulting occurs is considered. The behavior of the system is controlled by a single parameter, the fractal dimension. Deformation occurs on all scales of faults. The fractal dimension determines the fraction of the total displacement that occurs on the first-order or primary faults. The value of the fractal dimension can be obtained from the frequency-magnitude relation for earthquakes. The results are applied to the San Andreas fault system in central California. Earthquake studies give D = 1.90. The main strand of the San Andreas fault is associated with the primary faults of the fractal system. It is predicted that the relative velocity across the main strand is 2.93 cm/yr. The remainder of the relative velocity of 5.5 cm/yr between the Pacific and North American plates occurs on higher-order faults. The predicted value is in reasonably good agreement with the value 3.39 + or - 0.29 cm/yr obtained from geological studies.

Turcotte, D. L.↗

Modeling brittle high-temperature composite structures

A major inhibitor to the projected application of brittle matrix composite structures, such as ceramic matrix composites (CMC), to advanced high temperature airframes is premature matrix microcracking. However, the phenomenon has been experienced with polymer matrix composite (PMC) systems that are now in an advanced stage of development. Consequently, the subsequent mechanical behavior of these composite systems must be predictable, which demands a more complete understanding of the mechanisms such as nonlinear features, localized breakdown of the matrix, interface, and fiber phases. To establish the need for such development we will commence with a review of the key differences between polymer, and ceramic as well as metal matrix composites (MMC). Next some lessons learned from the development of the most mature PMC structures will be discussed to indicate how this experience may be carried over to assist the development on high temperature composites subjected to monotonic and cyclic loads as discussed and selected areas for future technical development are indicated.

Kedward, K. T.↗

An investigation of the thermoviscoplastic behavior of a metal matrix composite at elevated temperatures

This research investigates the response of a fiberless 13 layer hot isostatically pressed Ti-15-3 laminate to creep, constant strain rate, and cyclic constant strain rate loading at temperatures ranging from 482C to 649C. Creep stresses from 48 to 260 MPa and strain rates of .0001 to .01 m/m/sec were used. Material parameters for three unified constitutive models (Bodner-Partom, Miller, and Walker models) were determined for Ti-15-3 from the experimental data. Each of the three models was subsequently incorporated into a rule of mixtures and evaluated for accuracy and ease of use in predicting the thermoviscoplastic response of unidirectional metal matrix composite laminates (both 0 and 90). The laminates were comprised of a Ti-15-3 matrix with 29 volume percent SCS6 fibers. The predicted values were compared to experimentally determined creep and constant strain rate data. It was found that all three models predicted the viscoplastic response of the 0 specimens reasonably well, but seriously underestimated the viscoplastic response of the 90 specimens. It is believed that this discrepancy is due to compliant and/or weak fiber-matrix interphase. In general, it was found that of the three models studied, the Bodner-Partom model was easiest to implement, primarily because this model does not require the use of cyclic constant strain rate tests to determine the material parameters involved. However, the version of the Bodner-Partom model used in this study does not include back stress as an internal state variable, and hence may not be suitable for use with materials which exhibit a pronounced Baushinger effect. The back stress is accounted for in both the Walker and Miller models; determination of the material parameters associated with the Walker model was somewhat easier than in the Miller model.

Rogacki, John R.↗

An independent refinement and integration procedure in multiregion finite element analysis

An independent refinement and integration procedure is developed to couple together independently modeled (global and local) regions in a single analysis. The models can have different levels of refinement and along the interface between them the finite element nodes need not coincide with one another. In the local model all the nodes except the nodes at the interface are statically condensed and the reduced stiffness matrix is obtained. For this static condensation a modified frontal solution technique is employed. A spline interpolation function that satisfies the linear isotropic plate bending differential equation is used to relate the local model interface nodal displacements to the global model interface displacements. The proposed independent refinement and integration procedure is evaluated by applying it to two- and three-dimensional cases involving inplane and out-of-plane deformation. The procedure yielded very accurate results for all the examples studied.

Krishnamurthy, T.↗

Remote sensing of Earth terrain

The layered random medium model is used to investigate the fully polarimetric scattering of electromagnetic waves from vegetation. The vegetation canopy is modeled as an anisotropic random medium containing nonspherical scatterers with preferred alignment. The underlying medium is considered as a homogeneous half space. The scattering effect of the vegetation canopy are characterized by 3-D correlation functions with variances and correlation lengths respectively corresponding to the fluctuation strengths and the physical geometries of the scatterers. The strong fluctuation theory is used to calculate the anisotropic effective permittivity tensor of the random medium and the distorted Born approximation is then applied to obtain the covariance matrix which describes the fully polarimetric scattering properties of the vegetation field. This model accounts for all the interaction processes between the boundaries and the scatterers and includes all the coherent effects due to wave propagation in different directions such as the constructive and destructive interferences. For a vegetation canopy with low attenuation, the boundary between the vegetation and the underlying medium can give rise to significant coherent effects.

Kong, Jin AU↗

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

Aerothermal modeling

The objectives, approach, and status of a program to develop the computational fluid dynamics tools needed to improve combustor design and analysis are outlined. The calculation procedure selected consists of a finite difference solution of the time averaged, steady state, primitive variable, elliptic form of the Reynolds equations. Standard TEACH type numerics are used to solve the resulting equations. These include hybrid differencing, SIMPLE algorithm for the pressure field, line by line iterative solution using the ADI method and the tridiagonal matrix algorithm (TDMA). Convergence is facilitated by using under relaxation. The physical processes are modeled by a two equation eddy viscosity model for turbulence; combustion is represented by a simple, irreversible, one step chemical reaction whose rate is influenced only by the time scale of the turbulence. The model evaluation procedure is also described.

Sturgess, G. J.↗

C-H Hot Bands in the Near-IR Emission Spectra of Leonids

The reported infrared (IR) emission spectra from 1999 Leonid fireballs show a 3.4 micron C-H emission band and unidentified bands at longer wavelengths. Upon atmospheric entry, the Leonid meteorites were flash-heated to temperatures around 2400K, which would destroy any organics on the surface of the meteorite grains. We propose that the nu(sub )CH emission band in the Leonid emission spectra arises from matrix-embedded C(sub n)-H-O entities that are protected from instant pyrolysis. Our model is based on IR absorption nu(sub )CH bands, which we observed in laboratory-grown MgO and natural olivine single crystals, where they arise from C(sub n)-H-O units imbedded in the mineral matrix, indicative of aliphatic -CH2- and -CH3 organics. Instead of being pyrolyzed, the C(sub n)-H-O entities in the Leonid trails become vibrationally excited to higher levels n = 1, 2, 3 etc. During de-excitation they emit at 3.4 microns, due to the (0 => 1) transition, and at longer wavelengths, due to hot bands. As a first step toward verifying this hypothesis we measured the C-H vibrational manifold of hexane (C6H14). The calculated positions of the (2 => l ) , (3 => 2), and possibly (4 => 3) hot bands agree with the Leonid emission bands at 3.5, 3.8 and 4.l microns.

Freund, F. T.↗