Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATRIX STRESS CALCULATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Thermal/Mechanical Durability of Polymer-Matrix Composites in Cryogenic Environments

In order to increase the reliability of the next generation of space transportation systems, the mechanical behavior of polymeric-matrix composite (PMC) materials at cryogenic temperatures must be investigated. This paper presents experimental data on the residual mechanical properties of a carbon fiber polymeric composite, IM7/PETI-5 as a function of temperature and aging. Tension modulus and strength were measured at room temperature, -196 C, and -269 C on five different specimens ply lay-ups. Specimens were preconditioned with one set of coupons being isothermally aged for 576 hours at -184 C, in an unloaded state. Another set of corresponding coupons were mounted in constant strain fixtures such that a constant uniaxial strain was applied to the specimens for 576 hours at -184 C. A third set was mechanically cycled in tension at -184 C. The measured properties indicated that temperature, aging, and loading mode can all have significant influence on performance. Moreover, this influence is a strong function of laminate stacking sequence. Thermal-stress calculations based on lamination theory predicted that the transverse tensile ply stresses could be quite high for cryogenic test temperatures. Microscopic examination of the surface morphology showed evidence of degradation along the exposed edges of the material because of aging at cryogenic temperatures. ________________

Gates, Thomas S.↗
NASA NTRS1984

Advanced composite combustor structural concepts program

An analytical study was conducted to assess the feasibility of and benefits derived from the use of high temperature composite materials in aircraft turbine engine combustor liners. The study included a survey and screening of the properties of three candidate composite materials including tungsten reinforced superalloys, carbon-carbon and silicon carbide (SiC) fibers reinforcing a ceramic matrix of lithium aluminosilicate (LAS). The SiC-LAS material was selected as offering the greatest near term potential primarily on the basis of high temperature capability. A limited experimental investigation was conducted to quantify some of the more critical mechanical properties of the SiC-LAS composite having a multidirection 0/45/-45/90 deg fiber orientation favored for the combustor linear application. Rigorous cyclic thermal tests demonstrated that SiC-LAS was extremely resistant to the thermal fatigue mechanisms that usually limit the life of metallic combustor liners. A thermal design study led to the definition of a composite liner concept that incorporated film cooled SiC-LAS shingles mounted on a Hastelloy X shell. With coolant fluxes consistent with the most advanced metallic liner technology, the calculated hot surface temperatures of the shingles were within the apparent near term capability of the material. Structural analyses indicated that the stresses in the composite panels were low, primarily because of the low coefficient of expansion of the material and it was concluded that the dominant failure mode of the liner would be an as yet unidentified deterioration of the composite from prolonged exposure to high temperature. An economic study, based on a medium thrust size commercial aircraft engine, indicated that the SiC-LAS combustor liner would weigh 22.8N (11.27 lb) less and cost less to manufacture than advanced metallic liner concepts intended for use in the late 1980's.

Sattar, M. A.↗
NASA NTRS2002

Thermal/Mechanical Response and Damage Growth in Polymeric Composites at Cryogenic Temperatures

In order to increase the reliability of the next generation of space transportation systems, the mechanical behavior of polymeric matrix composite (PMC) materials at cryogenic temperatures must be investigated. This paper presents experimental data on the residual mechanical properties of a carbon fiber polymeric composite, IM7/PETI-5 both before and after aging at cryogenic temperatures. Tension modulus and strength were measured at room temperature, -196 C, and -269 C on five different specimen ply lay-ups, [0](sub 12), [90](sub 12), [+/-45](sub 3S), [+/-25](sub 3s) and [45,90(sub 3),-45,0(sub 3),-45,90(sub 3),45]. Specimens were preconditioned with one set of coupons being isothermally aged for 555 hours at -184 C in an unloaded state. Another set of corresponding coupons were mounted in constant displacement fixtures such that a constant uniaxial strain was applied to the specimens for 555 hours at -184 C. The measured lamina level properties indicated that cryogenic temperatures have an appreciable influence on behavior, and residual stress calculations based on lamination theory showed that the transverse tensile ply stresses could be quite high for cryogenic test temperatures. Microscopic examination of the surface morphology showed evidence of degradation along the exposed edges of the material due to aging at cryogenic temperatures.

Whitley, Karen S.↗
DOE OSTI2018

A quantifiable and automated volume fraction characterization technique for secondary and tertiary γ′ precipitates in Ni-based superalloys

Highlights: • A new, automated γ′ volume fraction characterization technique has been introduced for Ni-base superalloys. • New image processing techniques can automatically and accurately segment γ′ precipitates from the surrounding γ matrix. • CALPHAD and γ′ phase extraction measurements validate the new characterization technique. - Abstract: Ni-based superalloys are extensively used in high temperature, high stress applications. The presence of coherent γ′-Ni{sub 3}Al precipitates provides the requisite elevated temperature strength and creep resistance in these alloys. Therefore, quantification of the precipitate size distribution and volume fraction is essential for accurate model predictions and future alloy design. Current γ′ volume fraction characterization techniques suffer from a lack of precision and repeatability. Without a reproducible procedure, unacceptable variations and errors have plagued the past 6 decades of superalloys research and have inhibited the implementation of accurate precipitation models for these alloys. Through the use of high resolution, low kV, scanning electron microscopy and state-of-the-art image post-processing, a new, accurate γ′ characterization technique is presented. Using the new technique, secondary γ′ volume fractions in superalloy ME3 were calculated, and confirmed through phase extraction, thus opening the possibility for validation of presently used superalloy precipitation databases.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Characterization of elastic mechanical properties of Tuscaloosa Marine Shale from well logs using the vertical transversely isotropic model

To avoid steep declines in the Tuscaloosa Marine Shale (TMS) production, wells are fracture-stimulated to release the hydrocarbons trapped in the matrix of the formation. An accurate estimation of Young’s modulus and Poisson’s ratio is essential for hydraulic fracture propagation. In addition, ignoring the highly heterogeneous and anisotropic character of TMS can lead to erroneous stress values, which subsequently affect hydraulic fracture width estimates and the overall hydraulic fracturing process. We have developed an empirical 1D geomechanical model that takes into account VTI anisotropy, and it is used to characterize the elastic mechanical properties of TMS in two wells. In the analyzed formation, the vertical Poisson’s ratio is less than the horizontal Poisson’s ratio, which suggests the necessity of an alternative to the ANNIE equations. The stiffness coefficients [Formula: see text] and [Formula: see text] were estimated using the relationships developed from the ultrasonic core data available for the two TMS. Further, correlations between the static and dynamic properties from laboratory tests were used to improve the minimum horizontal stress calculation. We compare VTI Young’s moduli, Poisson’s ratios, and minimum horizontal stress with the isotropic solution. VTI modeling improves the estimation of the elastic mechanical properties. The isotropic solution underestimates the minimum horizontal stress in the formation. Moreover, it was shown that the 20 ft shale interval below the TMS base is characterized by a low Young’s modulus (the vertical Young’s modulus is equal to 20 GPa, whereas the horizontal Young’s modulus is equal to 40 GPa) and may be a frac barrier.

Geochemistry & Geophysics↗
NASA NTRS1989

The interface in tungsten fiber reinforced niobium metal-matrix composites

The creep resistance of tungsten fiber reinforced niobium metal-matrix composites was evaluated. The interface region between the fiber and matrix was characterized by microhardness and electron probe microanalysis measurements which indicated that its properties were between those of fiber and matrix. However, the measured properties of the composite exceeded those calculated by the rule of mixtures even when the interface zone was assumed to retain all the strength of the fiber. The composite structure appeared to enhance the strengths of both the fibers and the matrix above what they exhibited in stand-alone tests. The effect of fiber orientation and matrix alloy composition on the fiber/matrix interface were also evaluated. Small alloying additions of zirconium and tungsten to the niobium matrix affected the creep resistance of the composites only slightly. A decrease in the creep resistance of the composite with increasing zirconium content in the matrix was ascribed to an increase in the diffusion rate of the fiber/matrix interdiffusion reaction, and a slight increase in the creep resistance of the composite was observed with an addition of 9 w percent tungsten to the matrix. In addition, Kirkendall void formation was observed at the fiber/matrix interface; the void distribution differed depending on the fiber orientation relative to the stress axis.

Grobstein, Toni L.↗
NASA NTRS1993

Thermoelastic response of metal matrix composites with large-diameter fibers subjected to thermal gradients

A new micromechanical theory is presented for the response of heterogeneous metal matrix composites subjected to thermal gradients. In contrast to existing micromechanical theories that utilize classical homogenization schemes in the course of calculating microscopic and macroscopic field quantities, in the present approach the actual microstructural details are explicitly coupled with the macrostructure of the composite. Examples are offered that illustrate limitations of the classical homogenization approach in predicting the response of thin-walled metal matrix composites with large-diameter fibers when subjected to thermal gradients. These examples include composites with a finite number of fibers in the thickness direction that may be uniformly or nonuniformly spaced, thus admitting so-called functionally gradient composites. The results illustrate that the classical approach of decoupling micromechanical and macromechanical analyses in the presence of a finite number of large-diameter fibers, finite dimensions of the composite, and temperature gradient may produce excessively conservative estimates for macroscopic field quantities, while both underestimating and overestimating the local fluctuations of the microscopic quantities in different regions of the composite. Also demonstrated is the usefulness of the present approach in generating favorable stress distributions in the presence of thermal gradients by appropriately tailoring the internal microstructure details of the composite.

Aboudi, Jacob↗
NASA NTRS2006

Nonlinear Spring Finite Elements for Predicting Mode I-Dominated Delamination Growth in Laminated Structure with Through-Thickness reinforcement

One particular concern of polymer matrix composite laminates is the relatively low resistance to delamination cracking, in particular when the dominant type of failure is mode I opening. One method proposed for alleviating this problem involves the insertion pultruded carbon pins through the laminate thickness. The pins, known as z-pins, are inserted into the prepreg laminate using an ultrasonic hammer prior to the curing process, resulting in a field of pins embedded normal to the laminate plane as illustrated in Figure. 1. Pin diameters range between 0.28-mm to 0.5-mm and standard areal densities range from 0.5% to 4%. The z-pins are provided by the manufacturer, Aztex(Registered TradeMark) , in a low-density foam preform, which acts to stabilize orientation of the pins during the insertion process [1-3]. Typical pin materials include boron and carbon fibers embedded in a polymer matrix. A number of methods have been developed for predicting delamination growth in laminates reinforced with z-pins. During a study on the effect of z-pin reinforcement on mode I delamination resistance, finite element analyses of z-pin reinforced double cantilever beam (DCB) specimens were performed by Cartie and Partridge [4]. The z-pin bridging stresses were modeled by applying equivalent forces at the pin locations. Single z-pin pull-out tests were performed to characterize the traction law of the pins under mode I loading conditions. Analytical solutions for delamination growth in z-pin reinforced DCB specimens were independently derived by Robinson and Das [5] and Ratcliffe and O'Brien [6]. In the former case, pin bridging stresses were modeled using a distributed load and in the latter example the bridging stresses were discretely modeled by way of grounded springs. Additionally, Robinson and Das developed a data reduction strategy for calculating mode I fracture toughness, G(sub Ic), from a z-pin reinforced DCB specimen test [5]. In both cases a traction law similar to that adopted by Cartie and Partridge was used to represent z-pin failure under mode I loading conditions. In the current work spring elements available in most commercial finite element codes were used to model z-pins. The traction law used in previous analyses [4-6] was employed to represent z-pin damage. This method is intended for and is limited to simulating z-pins in composite laminate structure containing mode I-dominated delamination cracking. The current technique differs from previous analyses in that spring finite elements (available in commercial codes) are employed for simulating zpins, reducing the complexity of the analysis construction process. Furthermore, the analysis method can be applied to general structure that experiences mode I-dominated delamination cracking, in contrast to existing analytical solutions that are only applicable to coupon DCB specimens.

Ratcliffe, James G.↗
NASA NTRS1990

Materials with periodic internal structure: Computation based on homogenization and comparison with experiment

The combination of thermal and mechanical loading expected in practice means that constitutive equations of metal matrix composites must be developed which deal with time-independent and time-dependent irreversible deformation. Also, the internal state of composites is extremely complicated which underlines the need to formulate macroscopic constitutive equations with a limited number of state variables which represent the internal state at the micro level. One available method for calculating the macro properties of composites in terms of the distribution and properties of the constituent materials is the method of homogenization whose formulation is based on the periodicity of the substructure of the composite. A homogenization procedure was developed which lends itself to the use of the finite element procedure. The efficiency of these procedures, to determine the macroscopic properties of a composite system from its constituent properties, was demonstrated utilizing an aluminum plate perforated by directionally oriented slits. The selection of this problem is based on the fact that, extensive experimental results exist, the macroscopic response is highly anisotropic, and that the slits provide very high stress gradients which severely test the effectiveness of the computational procedures. Furthermore, both elastic and plastic properties were investigated so that the application to practical systems with inelastic deformation should be able to proceed without difficulty. The effectiveness of the procedures was rigorously checked against experimental results and with the predictions of approximate calculations. Using the computational results it is illustrated how macroscopic constitutive equations can be expressed in forms of the elastic and limit load behavior.

Jansson, S.↗
DOE OSTI2017

Influence of Alloying Elements and Effect of Stress on Anisotropic Hydrogen Diffusion in Zr-Based Alloys Predicted by Accelerated Kinetic Monte Carlo Simulations

The formation of hydrides degrades the toughness of Zr-based claddings, and thus is one of the primary concerns regarding cladding integrity during fuel operation and used fuel storage. However, the formation mechanism of hydrides, particularly the nucleation mechanism, is yet to be fully understood. To explore the nucleation process, atomic scale studies including density functional theory, molecular dynamics and kinetic Monte Carlo have been used to investigate H diffusion, segregation, and the possible nucleation paths of hydrides. Here, density functional theory calculations are used to parameterize a kinetic Monte Carlo model, which can accurately predict the effective H diffusivity in various Zr-based claddings including Zircaloy, Zirlo and Zr-2.5Nb. Assisted by density functional theory calculations, a recently developed charge-optimized-many-body potential is further improved and used to investigate the possible nucleation paths of hydrides. It’s shown that the formation of face-centered hydrides starts with coherent zeta-phase, shearing of which in the basal plane changes the hcp stacking into fcc and leads to the formation of gamma-hydride nuclei. The shearing involves a negligible barrier and negligible strain accumulation in the matrix by adopting three equivalent shear partials. The results obtained at the atomic scale are consistent with previous theory and experiments, and improves the current understanding on how and how fast hydride forms in Zr based cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2016

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

On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH

We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1997

Determination of Residual Stress in Composite Materials Using Ultrasonic Waves

The performance of high temperature composites can be significantly affected by the presence of residual stresses. These stresses arise during cooling processes from fabrication to room temperature due to mismatch of thermal expansion coefficients between matrix and fiber materials. This effect is especially pronounced in metal matrix and intermetallic composites. It can lead to plastic deformations, matrix cracking and fiber/matrix interface debonding. In this work the feasibility of ultrasonic techniques for residual stress assessment in composites is addressed. A novel technique for absolute stress determination in orthotropic materials from angular dependencies of ultrasonic velocities is described. The technique is applicable for determination of both applied and residual stresses and does not require calibration measurements on a reference sample. The important advantage of this method is that stress is determined simultaneously with stress-dependent elastic constants and is thus decoupled from the material texture. It is demonstrated that when the principal plane stress directions coincide with acoustical axes, the angular velocity data in the plane perpendicular to the stress plane may be used to determine both stress components. When the stress is off the acoustical axes, the shear and the difference of the normal stress components may be determined from the angular dependence of group velocities in the plane of stresses. Synthetic sets of experimental data corresponding to materials with different anisotropy and stress levels are used to check the applicability of the technique. The method is also verified experimentally. A high precision ultrasonic wave transmission technique is developed to measure angular dependence of ultrasonic velocities. Examples of stress determination from experimental velocity data are given. A method is presented for determination of velocities of ultrasonic waves propagating through the composite material with residual stresses. It is based on the generalized self-consistent multiple scattering model. Calculation results for longitudinal and shear ultrasonic wave velocities propagating perpendicular to the fibers direction in SCS-6/Ti composite with and without residual stresses are presented. They show that velocity changes due to presence of stresses are of order 1%.

Rokhlin, S. I.↗
DOE OSTI2017

A Comparison of Monte Carlo and Deterministic Solvers for keff and Sensitivity Calculations

Verification and validation of our solutions for calculating the neutron reactivity for nuclear materials is a key issue to address for many applications, including criticality safety, research reactors, power reactors, and nuclear security. Neutronics codes solve variations of the Boltzmann transport equation. The two main variants are Monte Carlo versus deterministic solutions, e.g. the MCNP [1] versus PARTISN [2] codes, respectively. There have been many studies over the decades that examined the accuracy of such solvers and the general conclusion is that when the problems are well-posed, either solver can produce accurate results. However, the devil is always in the details. The current study examines the issue of self-shielding and the stress it puts on deterministic solvers. Most Monte Carlo neutronics codes use continuous-energy descriptions of the neutron interaction data that are not subject to this effect. The issue of self-shielding occurs because of the discretisation of data used by the deterministic solutions. Multigroup data used in these solvers are the average cross section and scattering parameters over an energy range. Resonances in cross sections can occur that change the likelihood of interaction by one to three orders of magnitude over a small energy range. Self-shielding is the numerical effect that the average cross section in groups with strong resonances can be strongly affected as neutrons within that material are preferentially absorbed or scattered out of the resonance energies. This affects both the average cross section and the scattering matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2006

Anisotropic Laminar Piezocomposite Actuator Incorporating Machined PMN-PT Single Crystal Fibers

The design, fabrication, and testing of a flexible, laminar, anisotropic piezoelectric composite actuator utilizing machined PMN-32%PT single crystal fibers is presented. The device consists of a layer of rectangular single crystal piezoelectric fibers in an epoxy matrix, packaged between interdigitated electrode polyimide films. Quasistatic free-strain measurements of the single crystal device are compared with measurements from geometrically identical specimens incorporating polycrystalline PZT-5A and PZT-5H piezoceramic fibers. Free-strain actuation of the single crystal actuator at low bipolar electric fields (+/- 250 V/mm) is approximately 400% greater than that of the baseline PZT-5A piezoceramic device, and 200% greater than that of the PZT-5H device. Free-strain actuation under high unipolar electric fields (0-4kV/mm) is approximately 200% of the PZT-5A baseline device, and 150% of the PZT-5H alternate piezoceramic device. Performance increases at low field are qualitatively consistent with predicted increases based on scaling the low-field d33 piezoelectric constants of the respective piezoelectric materials. High-field increases are much less than scaled d33 estimates, but appear consistent with high-field freestrain measurements reported for similar bulk single-crystal and piezoceramic compositions. Measurements of single crystal actuator capacitance and coupling coefficient are also provided. These properties were poorly predicted using scaled bulk material dielectric and coupling coefficient data. Rules-of-mixtures calculations of the effective elastic properties of the single crystal device and estimated actuation work energy densities are also presented. Results indicate longitudinal stiffnesses significantly lower (50% less) than either piezoceramic device. This suggests that single-crystal piezocomposite actuators will be best suited to low induced-stress, high strain and deflection applications.

Wilkie, W. Keats↗
NASA NTRS2004

Anisotropic Piezocomposite Actuator Incorporating Machined PMN-PT Single Crystal Fibers

The design, fabrication, and testing of a flexible, planar, anisotropic piezoelectric composite actuator utilizing machined PMN-32%PT single crystal fibers is presented. The device consists of a layer of rectangular single crystal piezoelectric fibers in an epoxy matrix, packaged between interdigitated electrode polyimide films. Quasistatic free-strain measurements of the single crystal device are compared with measurements from geometrically identical specimens incorporating polycrystalline PZT-5A and PZT-5H piezoceramic fibers. Free-strain actuation of the single crystal actuator at low bipolar electric fields (+/- 250 V/mm) is approximately 400% greater than that of the baseline PZT-5A piezoceramic device, and 200% greater than that of the PZT-5H device. Free-strain actuation under high unipolar electric fields (0-4kV/mm) is approximately 200% of the PZT-5A baseline device, and 150% of the PZT-5H alternate piezoceramic device. Performance increases at low field are qualitatively consistent with predicted increases based on scaling the low-field d(sub 33) piezoelectric constants of the respective piezoelectric materials. High-field increases are much less than scaled d(sub 33) estimates, but appear consistent with high-field freestrain measurements reported for similar bulk single-crystal and piezoceramic compositions. Measurements of single crystal actuator capacitance and coupling coefficient are also provided. These properties were poorly predicted using scaled bulk material dielectric and coupling coefficient data. Rules-of-mixtures calculations of the effective elastic properties of the single crystal device and estimated actuation work energy densities are also presented. Results indicate longitudinal stiffnesses significantly lower (50% less) than either piezoceramic device. This suggests that single-crystal piezocomposite actuators will be best suited to low induced-stress, high strain and deflection applications.

Wilkie, W. Keats↗
NASA NTRSDate not supplied

Multiscale Prediction of Yarn Pullout Failure Mode in Unreinforced Textile Fabrics

Unreinforced woven fabrics have been implemented in a variety of high performance applications, including body armor, deployable structures, and as the reinforcement material in composites. Multiscale modeling techniques have significantly improved the capabilities of simulation-based tools to capture fabric mechanics efficiently and accurately, but often lack in their prediction of failure and require pairing with finite element analysis (FEA) software, limiting their application to the design of ‘fit-for-purpose’ materials. NASA’s Multiscale Analysis Tool (NASMAT) is a standalone multiscale program that has been traditionally used in the analysis of reinforced composites materials. More recently, it has been amended to simulate unreinforced fabric behavior by allowing the geometric state of the tows to change with applied loading due to the lack of a reinforcement material, such as the matrix seen in composites. Previous work has shown the ability of NASMAT to capture nonlinear macroscale behavior by predicting geometric changes in the state of each subcell as a function of the applied loading and allowing each subcell in the analysis to rotate according to these predicted changes, as well as predict nonlinear behavior due to the fiber breakage failure mode. In this work, the capability of predicting the onset and propagation of failure in plain woven fabrics in NASMAT is presented for the yarn pullout failure mode, which occurs when a fabric is loaded at an off-axis angle relative to the warp of weft tow direction. Yarn pullout behavior is initiated by determining the applied load in which the shear resistance of the contact area between yarn families is overcome. When failure is initiated, yarn pullout is determined to have occurred when the applied displacement, calculated from global strain, exceeds the deformed position of a given contact points between yarn families, determined from pin-joint kinematics. Contact points where pullout has occurred contribute to a global damage parameter used to modify the homogenized stiffness of the fabric, resulting in nonlinear behavior observed at the macroscale. The off-axis loading behavior and yarn pullout failure theory have been developed and implemented into NASMAT such that users can simulate off-axis tensile behavior of fabrics in a single, standalone multiscale tool. Simulations are compared to uniaxial tensile tests at various off-axis angles to demonstrate the capability of the tool in its prediction of both the onset of failure at each off-axis angles and the stress-strain behavior as failure progresses.

Materials↗
NASA NTRS2000

Assessment of Internal Fabric of Particulate Materials

Particle arrangement and distribution within a soil matrix has long been recognized as having significant influence on the mechanical behavior of cohesionless soils. It is well known that two soil specimens having the same grain type (e.g., quartz, feldspar, etc.), same grain size distribution and relative density (or void ratio) can display completely different mechanical behavior. Because of the different fabric configurations in the otherwise similar specimens, they are likely to have different mechanical properties such as stress-strain response, anisotropy, dilatancy, etc. Soil Fabric is defined as the arrangement of particles, particle groups and associated pore space. In the literature, fabric analysis techniques are mainly classified as destructive (e.g., specimen stabilization, thin-sectioning, and microscopy), and nondestructive techniques (e.g., magnetic resonance imaging, ultrasonic testing, x-ray radiography, and computed tomography). Quantifying the void ratio and its distribution is the main parameter used to describe the fabric of particulate materials. This paper presents a comprehensive literature review of fabric analysis techniques applied to particulate materials. In addition, the results of a comprehensive investigation to quantify void ratio of sand specimens will be presented and discussed. The sand used in the experiments in a natural, uniform rounded to sub-rounded silica sand known as F-75 banding sand with mean particle size of 0.22 mm. Uniform specimens and specimens subjected to different axial-strain levels tested under triaxial and biaxial conditions are examined to evaluate void ratio evolution and its distribution using destructive thin-sectioning and nondestructive Computed Tomography (CT) techniques. Details of a new innovative polygon generation technique called Voronoi tessellation used to quantify void ratio of microscopic images of sand grains will be presented and compared to classical Oda's method. Finally, frequency distribution of void ratio calculated using CT analysis will be presented and discussed.

Alshibi, Khalid A.↗