Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Composite 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Oxidation of Boron Nitride in Composites

Boron nitride (BN) is a prime candidate for fiber coatings in silicon carbide (SiC) fiber-reinforced SiC matrix composites. The properties of BN allow the fiber to impart beneficial composite properties to the matrix, even at elevated temperatures. The problem with BN is that it is readily attacked by oxygen. Although BN is an internal component of the composite, a matrix crack or pore can create a path for hot oxygen to attack the BN. This type of attack is not well understood. A variety of phenomena have been observed. These include borosilicate glass formation, volatilization of the BN, and under some conditions, preservation of the BN. In this study at the NASA Lewis Research Center, a series of BN materials and BN-containing model composites were methodically examined to understand the various issues dealing with the oxidation of BN in composites. Initial studies were done with a series of monolithic BN materials prepared by hot pressing and chemical vapor deposition (CVD). From these studies, we found that BN showed a strong orientation effect in oxidation and was extremely sensitive to the presence of water vapor in the environment. In addition, CVD material deposited at a high temperature showed much better oxidation behavior than CVD material deposited at a lower temperature.

Jacobson, Nathan S.↗
NASA NTRSDate not supplied

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar M Rao↗
NASA NTRSDate not supplied

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar Rao↗
NASA NTRS1995

Micromechanical models for textile structural composites

The objective is to develop micromechanical models for predicting the stiffness and strength properties of textile composite materials. Two models are presented to predict the homogeneous elastic constants and coefficients of thermal expansion of a textile composite. The first model is based on rigorous finite element analysis of the textile composite unit-cell. Periodic boundary conditions are enforced between opposite faces of the unit-cell to simulate deformations accurately. The second model implements the selective averaging method (SAM), which is based on a judicious combination of stiffness and compliance averaging. For thin textile composites, both models can predict the plate stiffness coefficients and plate thermal coefficients. The finite element procedure is extended to compute the thermal residual microstresses, and to estimate the initial failure envelope for textile composites.

Marrey, Ramesh V.↗
NASA NTRS1988

Modeling creep behavior of fiber composites

A micromechanical model for the creep behavior of fiber composites is developed based on a typical cell consisting of a fiber and the surrounding matrix. The fiber is assumed to be linearly elastic and the matrix nonlinearly viscous. The creep strain rate in the matrix is assumed to be a function of stress. The nominal stress-strain relations are derived in the form of differential equations which are solved numerically for off-axis specimens under uniaxial loading. A potential function and the associated effective stress and effective creep strain rates are introduced to simplify the orthotropic relations.

Chen, J. L.↗
NASA NTRS1990

Nonmethane hydrocarbon chemistry in the remote marine boundary layer

A photochemical model of the remote marine boundary layer (MBL) is presented, with focus placed on the role of reactive nonmethane hydrocarbons (NMHC). A wide range of NMHC air-sea fluxes with various relative distributions of NMHC regions are considered. In particular, the flux magnitude at which NMHC emissions become significant, and then dominant, players in MBL chemistry is identified. Emphasis is placed on diurnal variability, diurnal ozone variations and sensitivity to NMHC emission fluxes, to CO, O3, H2O, and UV light, and to kinetics and isometric composition. Model runs indicate that, in the range consistent with current observations, the NMHCs may either dominate MBL chemistry, or simply be contributors at the 10-percent level. These model runs also show that existing observations of NMHCs in ocean water find them to scarce for fluxes from bulk-flux air-sea gas exchange models to be consistent with the fluxes needed in the proposed model to maintain the lowest observed MBL NMHC.

Donahue, Neil M.↗
NASA NTRS1999

A Shell/3D Modeling Technique for Delaminations in Composite Laminates

A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provide a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with plate or shell finite elements. For simple double cantilever beam (DCB), end notched flexure (ENF), and single leg bending (SLB) specimens, mixed mode energy release rate distributions were computed across the width from nonlinear finite element analyses using the virtual crack closure technique. The analyses served to test the accuracy of the shell/3D technique for the pure mode I case (DCB), mode II case (ENF) and a mixed mode I/II case (SLB). Specimens with a unidirectional layup where the delamination is located between two 0 plies, as well as a multidirectional layup where the delamination is located between two non-zero degree plies, were simulated. For a local 3D model extending to a minimum of about three specimen thicknesses in front of and behind the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗
NASA NTRS1989

Composite fracture using the normal stress ratio theory

Experiments were performed on center-notched specimens of unidirectional AS4/3501-6 graphite-epoxy subjected to varied far-field loading conditions. Specimens tested include on-axis and off-axis tensile coupons and Iosipescu shear specimens, allowing the study of normal, shear, and combined far-field loading conditions. From the tests, the direction of initial crack extension and the far-field stresses at crack initiation and fracture are obtained. The normal stress ratio theory is subsequently applied within a macroscopic-level homogeneous anisotropic stress solution to predict the direction of crack extension, critical stresses, and other crack growth behavior for the experiments. Comparison between experimental results and theoretical predictions indicates that application of the normal stress ratio theory within a macroscopic-level stress analysis has significant potential as an approach to modeling composite crack growth.

Beuth, J. L., Jr.↗
NASA NTRS1997

Global Geodetic Signatures of the Antarctic Ice Sheet

Four scenarios of present day Antarctic ice sheet mass change are developed from comprehensive reviews of the available glaciological and oceanographic evidence. The gridded scenarios predict widely varying contributions to secular sea level change xi ranging from -1.1 to 0.45 mm/yr, and predict polar motion m and time-varying low-degree gravitational coefficients J(sub l) that differ significantly from earlier estimates. A reasonably linear relationship between the rate of sea level change from Antarctica xi(sub A) and the predicted Antarctic J(sub l) is found for the four scenarios. This linearity permits a series of forward models to be constructed that incorporate the effects of ice mass changes in Antarctica, Greenland, and distributed smaller glaciers, as well as postglacial rebound (assuming the ICE-3G deglaciation history), with the goal of obtaining optimum reconciliation between observed constraints on J(sub l) and sea level rise xi. Numerous viable combinations of lower mantle viscosity and hydrologic sources are found that safely "observed" in the range of 1 to 2-2.5 mm/yr and observed J(sub l) for degrees 2, 3, and 4. In contrast, rates of global sea level rise above 2.5 mm/yr are inconsistent with available J(sub l) observations. The successful composite models feature a pair of lower mantle viscosity solutions arising from the sensitivity of J(sub l) to glacial rebound. The paired values are well separated at xi = 1 mm/yr, but move closer together as xi is increased, and, in fact, merge around xi = 2 - 2.5 mm/yr, revealing an intimate relation between xi and preferred lower mantle viscosity. This general pattern is quite robust and persists for different J(sub l) solutions, for variations in source assumptions, and for different styles of lower mantle viscosity stratification. Tighter J(sub l) constraints for l greater than 2 may allow some viscosity stratification schemes and source assumptions to be excluded in the future. For a given total observed xi the sea level rise from Antarctica xi(sub A) is tightly constrained and ranges from 0 to + 1 mr/yr (corresponding to an ablating ice sheet) as estimates of xi are raised from 1 to 2.5 mm/yr. However, when the degree 3 zonal harmonic constraint is removed, the solutions show little sensitivity to Antarctic mass balance, emphasizing the need for a well determined odd-degree secular zonal harmonic for determining polar ice mass balance.

James, Thomas S.↗
NASA NTRS1978

A composition and thickness model for lunar impact crater and basin deposits

A model was developed for determining the composition and thicknesses of continuous deposits of large lunar impact craters and basins. Results of a photogeologic study and topographic analysis of continuous deposits of the lunar crater Delisle, which support the model, show that numerous secondary craters and crater chains with concentric dunes on their uprange rims occur well within the mapped deposits. At any given radius from Delisle, the secondary craters are of equal freshness and not vastly different in size, but with increasing radius they become better defined. The upper surface of the continuous deposits reaches the level of the pre-existing mare surface well within the mapped boundary of the deposits.

Morrison, R. H.↗
NASA NTRS2015

Multiscale Modeling of Ceramic Matrix Composites

Results of multiscale modeling simulations of the nonlinear response of SiC/SiC ceramic matrix composites are reported, wherein the microstructure of the ceramic matrix is captured. This micro scale architecture, which contains free Si material as well as the SiC ceramic, is responsible for residual stresses that play an important role in the subsequent thermo-mechanical behavior of the SiC/SiC composite. Using the novel Multiscale Generalized Method of Cells recursive micromechanics theory, the microstructure of the matrix, as well as the microstructure of the composite (fiber and matrix) can be captured.

ceramic matrix composites↗
NASA NTRS1986

On fracture of fibrous composites

The problem of an existing through-the-thickness center notch in a unidirectional fibrous composite is investigated using a homogeneous anisotropic elasticity analysis in conjunction with various crack growth criteria. Tensor polynomial, strain energy density, and normal stress ratio criteria are applied as crack growth direction criteria and compared to experimental observations. For the cases and criteria studied, only the normal stress ratio criterion exhibits an ability to correctly predict crack growth direction. The normal stress ratio criterion is subsequently applied to the same cases as a failure criterion. Comparison with experiment indicates that these predictions exhibit reasonable correlation with observed critical stresses. Finally, the macroscopic approach to the problem used in this study is compared to a microscopic shear lag analysis of experimental results in the literature. The comparison indicates that both approaches yield essentially identical predictions. The comparisons made with other criteria and approaches indicate that application of the normal stress ratio within a macroscopic level analysis has great potential as an approach to modeling composite crack growth.

Beuth, J. L., Jr.↗
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRS2019

Micromechanical Modeling of High-Strain Thin-Ply Composites

This paper presents a micromechanical analysis for the elastic and viscoelastic behavior of high-strain thin-ply composites. The modeling approach is based on unit cell homogenization. The geometry of the internal weave architecture and ply configuration is characterized via micrographic analysis and explicitly modeled in the unit cell. The composites are modeled as Kirchhoff plates and the homogenization analysis computes the effective relaxation ABD matrix represented by Prony series using the elastic and viscoelastic properties of the constituent fiber and matrix. The formulation of the micromechanical model and numerical implementation d. Composite laminates with 3-ply and 4-ply configurations are studied.

Hamillage, Milinda Yapa↗