Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress strain relations”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Performance of Chevron-notch short bar specimen in determining the fracture toughness of silicon nitride and aluminum oxide

Ease of preparation and testing are advantages unique to the chevron-notch specimen used for the determination of the plane strain fracture toughness of extremely brittle materials. During testing, a crack develops at the notch tip and extends stably as the load is increased. For a given specimen and notch configuration, maximum load always occurs at the same relative crack length independent of the material. Fracture toughness is determined from the maximum load with no need for crack length measurement. Chevron notch acuity is relatively unimportant since a crack is produced during specimen loading. In this paper, the authors use their previously determined stress intensity factor relationship for the chevron-notch short bar specimen to examine the performance of that specimen in determining the plane strain fracture toughness of silicon nitride and aluminum oxide.

Munz, D.↗
DOE OSTI2025

Enrichment of root-associated Streptomyces strains in response to drought is driven by diverse functional traits and does not predict beneficial effects on plant growth

The genus Streptomyces has consistently been found enriched in drought-stressed plant root microbiomes, yet the ecological basis and functional variation underlying this enrichment at the strain and isolate level remain unclear. Using two 16S rRNA sequencing methods with different levels of taxonomic resolution, we confirmed drought-associated enrichment (DE) of Streptomyces in field-grown sorghum roots and identified five closely related but distinct amplicon sequence variants (ASVs) belonging to the genus with variable drought enrichment patterns. From a culture collection of sorghum root endophytes, we selected 12 Streptomyces isolates representing these ASVs for phenotypic and genomic characterization. Whole-genome sequencing revealed substantial variation in gene content, even among closely related isolates, and exometabolomic profiling showed distinct metabolic responses to media supplemented with drought- versus well-watered root tissue. Traits linked to drought survival, including osmotic stress tolerance, siderophore production, and carbon utilization, varied widely among isolates and were not phylogenetically conserved. Using a broader panel of 48 Streptomyces, we demonstrate that DE scores, determined through mono-association experiments in gnotobiotic sorghum systems, showed high variability and lacked correlation with plant growth promotion. Pangenome-wide association identified orthogroups involved in osmolyte transport (e.g., proP) and membrane biosynthesis (e.g., fabG) as positively associated with DE, though most associations lacked phylogenetic signal. Collectively, these results demonstrate that Streptomyces DE is not a conserved genus-level trait but is instead strain-specific and functionally heterogeneous. Furthermore, DE in the root microbiome was shown not to predict beneficial effects on plant growth. This work underscores the need to resolve functional traits at the strain level and highlights the complexity of microbe-host-environment interactions under abiotic stress.

Fonseca-Garcia, Citlali↗
NASA NTRS1985

Strain energy release rate determination of stress intensity factors by finite element methods

The stiffness derivative finite element technique is used to determine the Mode I stress intensity factors for three-crack configurations. The geometries examined include the double edge notch, single edge notch, and the center crack. The results indicate that when the specified guidelines of the Stiffness Derivative Method are used, a high degree of accuracy can be achieved with an optimized, relatively coarse finite element mesh composed of standard, four-node, plane strain, quadrilateral elements. The numerically generated solutions, when compared with analytical ones, yield results within 0.001 percent of each other for the double edge crack, 0.858 percent for the single edge crack, and 2.021 percent for the center crack.

Walsh, R. M., Jr.↗
DOE OSTI2023

The deformation and burst behavior of Zircaloy-4 cladding tubes with hydride rim features subject to internal pressure loads

This study is motivated by the desire to conduct a separate effects test replicating specific aspects of an in-pile pulse reactor test on high exposure light water reactor (LWR) fuel rods. The principal purpose of the tests was to determine whether temperature and or hydrogen concentration thresholds exist in assessing the vulnerability of a stress relieved zirconium alloy cladding tubes to low strain ruptures. A secondary objective was to determine whether observations of the cladding stress during deformation could aid in the understanding of the crack formation and propagation in cladding tubes with hydride rim features. Here, the study utilized artificially hydrided Zircaloy-4 cladding tubes and an isothermal, internal pressure test that was controlled with an in-situ strain measurement. The test results show a clear transition in failure vulnerability between room temperature and 150 °C. Below 150 °C, cladding tubes rupture at very low uniform plastic strain—even at relatively modest cladding hydrogen levels. At and above 150 °C, cladding tubes only rupture at low uniform plastic strain if the cladding hydrogen content is greater than 500 ppm and usually experience at least some modest amount of plastic deformation prior to rupture. Post-test metallographs from non-ruptured cladding tubes show numerous primary cracks forming in the hydride rim and arresting at the metal/hydride interface. This observation combined with the decrease in bulk yielding of hydrided cladding tubes when compared to fresh non-hydrided tubes leads to the conclusion that the fracture stress of the hydride rims is related to and bounded by the yield stress of the zirconium metal. Some evidence is presented for a change in crack propagation mechanism from one of void formation on hydride plates and linkage at temperatures below 150 °C to more classical ductile shearing beneath the primary crack tip at temperatures above 150 °C. This apparent change in crack propagation mechanism may help explain the observed temperature and hydrogen concentration thresholds.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1972

Thin-film light-intensity measurement strain-analysis technique.

The optical response to loading of a thin metallic film deposited on a low-modulus structural substrate is studied theoretically and experimentally. Two types of optical properties called total and central-image transmittance (or reflectance) are shown to be related to the mechanical state of the substrate. Empirical optical-mechanical relationships are proposed between these optical properties and the substrate strain field of a general plane-stress problem. A technique based on wrinkle and microfracture patterns is described for determining principal directions of strain. Experimental results for uniaxially loaded specimens show that it is possible to obtain a nearly linear relationship between transmittance and strain for certain materials combinations.

Williams, J. G.↗
NASA NTRS2005

Nonlinear Visco-Elastic Response of Composites via Micro-Mechanical Models

Micro-mechanical models for a study of nonlinear visco-elastic response of composite laminae are developed and their performance compared. A single integral constitutive law proposed by Schapery and subsequently generalized to multi-axial states of stress is utilized in the study for the matrix material. This is used in conjunction with a computationally facile scheme in which hereditary strains are computed using a recursive relation suggested by Henriksen. Composite response is studied using two competing micro-models, viz. a simplified Square Cell Model (SSCM) and a Finite Element based self-consistent Cylindrical Model (FECM). The algorithm is developed assuming that the material response computations are carried out in a module attached to a general purpose finite element program used for composite structural analysis. It is shown that the SSCM as used in investigations of material nonlinearity can involve significant errors in the prediction of transverse Young's modulus and shear modulus. The errors in the elastic strains thus predicted are of the same order of magnitude as the creep strains accruing due to visco-elasticity. The FECM on the other hand does appear to perform better both in the prediction of elastic constants and the study of creep response.

Gates, Thomas S.↗
NASA NTRS1981

Strain buildup and release, earthquake prediction and selection of VBL sites for margins of the north Pacific

During the past year, the grant supported research on several aspects of crustal deformation. The relation between earthquake displacements and fault dimensions was studied in an effort to find scaling laws that relate static parameters such as slip and stress drop to the dimensions of the rupture. Several implications of the static relations for the dynamic properties of earthquakes such as rupture velocity and dynamic stress drop were proposed. A theoretical basis for earthquake related phenomena associated with slow rupture growth or propagation, such as delayed multiple events, was developed using the stress intensity factor defined in fracture mechanics and experimental evidence from studies of crack growth by stress corrosion. Finally, extensive studies by Japanese geologists have established the offset across numerous faults in Japan over the last one hundred thousand years. These observations of intraplate faulting are being used to establish the spatial variations of the average strain rate of subregions in southern Japan.

Scholz, C. H.↗
NASA NTRS1999

Composite Fuselage Technology

Work was conducted over a ten-year period to address the central issue of damage in primary load-bearing aircraft composite structure, specifically fuselage structure. This included the three facets of damage resistance, damage tolerance, and damage arrest. Experimental, analytical, and numerical work was conducted in order to identify and better understand the mechanisms that control the structural behavior of fuselage structures in their response to the three aspects of damage. Furthermore, work was done to develop straightforward design methodologies that can be employed by structural designers in preliminary design stages to make intelligent choices concerning the material, layup, and structural configurations so that a more efficient structure with structural integrity can be designed and built. Considerable progress was made towards achieving these goals via this work. In regard to damage tolerance considerations, the following were identified as important effects: composite layup and associated orthotropy/structural anisotropy, specifics of initial local damage mechanisms, role of longitudinal versus hoop stress, and large deformation and associated geometric nonlinearity. Means were established to account for effects of radius and for the nonlinear response. In particular, nondimensional parameters were identified to characterize the importance of nonlinearity in the response of pressurized cylinders. This led to the establishment of a iso-nonlinear-error plot for reference in structural design. Finally, in the case of damage tolerance, the general approach of the original methodology to predict the failure pressure involving extending basic plate failure data by accounting for the local stress intensification was accomplished for the general case by accounting for the mechanisms noted by utilizing the capability of the STAGS finite element code and numerically calculating the local stress intensification for the particular configuration to be considered. For the issue of damage arrest, placement of and configuration of stiffeners (including stiffener curvature), and magnitude and orientation of principal strains due to local bending were found to be key considerations. Means were established to account for stiffener effectiveness quantitatively based on radius, slit size, stiffener curvature' and relative bending stifffiesses involved. Geometric nonlinearity was also found to play an - 24 - important role here. Furthermore, it was determined that damage propagation is controlled by different mechanisms (hoop stress versus flapping stress and the associated factors involved in each) depending upon the direction of damage propagation. This latter item results in an inability to scale these phenomena in one test due to the different factors involved. Finally, the importance of shell curvature and associated instability in response to transverse loading including impact were found to be important considerations in damage resistance. A technique, involving asymmetric meshing of a finite element mesh, was developed to predict this behavior and showed excellent correlation with experimental results. Further details of these ten years of work are presented herein with references made to the fourteen documents produced during this work where full details can be found. Implications of this work are discussed and recommendations made. Although it is clear that there is more work to be done to fully understand composite fuselage technology and specifically the overall issue of damage in primary load-bearing composite structures, important understanding and capability has been extended via this work.

Lagace, Paul A.↗
DOE OSTI2021

Measuring residual stress in Ti-6Al-4V with HR-EBSD, using reference patterns from annealed material

In this study, elastic strain and residual stress values in an engineering material were determined by high-resolution EBSD (HR-EBSD) using reference patterns from annealed, recrystallized material. Specifically, rolled Ti-6Al-4V sheet was analyzed using HR EBSD in both an as-annealed and a strained state, the latter established via a static two-point bend fixture. The commercial software CrossCourt4 (CC4) was used to evaluate the relative intragranular strain tensor components of the alpha grains mapped by HR-EBSD. The intragranular reference patterns identified in each grain by the software were then replaced with similarly-oriented patterns from scans of the as-annealed, recrystallized material. The CC4 analysis was then repeated to provide calibrated strain and residual stress values. Elastic strains on the order of 10–3 were detected, with corresponding residual stress values on the order of 1 GPa. Subsequent analyses using high energy synchrotron x-ray diffraction showed good correlation with the elastic strains measured by HR-EBSD using the annealed reference patterns. Additionally, the maximum macroscopic stress analytically predicted for two-point bending was in agreement with the residual stress observed. Good repeatability of the methodology and insight into the state of individual grains were also observed. The findings provide encouraging evidence that accurate intragranular and long-range residual stress measurements of engineering significance can be made with HR-EBSD using a physical reference pattern library.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Modeling and life prediction methodology for Titanium Matrix Composites subjected to mission profiles

Titanium matrix composites (TMC) are being evaluated as structural materials for elevated temperature applications in future generation hypersonic vehicles. In such applications, TMC components are subjected to complex thermomechanical loading profiles at various elevated temperatures. Therefore, thermomechanical fatigue (TMF) testing, using a simulated mission profile, is essential for evaluation and development of life prediction methodologies. The objective of the research presented in this paper was to evaluate the TMF response of the (0/90)2s SCS-6/Timetal-21S subjected to a generic hypersonic flight profile and its portions with a temperature ranging from -130 C to 816 C. It was found that the composite modulus, prior to rapid degradation, had consistent values for all the profiles tested. A micromechanics based analysis was used to predict the stress-strain response of the laminate and of the constituents in each ply during thermomechanical loading conditions by using only constituent properties as input. The fiber was modeled as elastic with transverse orthotropic and temperature dependent properties. The matrix was modeled using a thermoviscoplastic constitutive relation. In the analysis, the composite modulus degradation was assumed to result from matrix cracking and was modeled by reducing the matrix modulus. Fatigue lives of the composite subjected to the complex generic hypersonic flight profile were well correlated using the predicted stress in 0 degree fibers.

Mirdamadi, M.↗
DOE OSTI2021

Development of Digital Twin Predictive Model for PWR Components: Updates on Multi Times Series Temperature Prediction Using Recurrent Neural Network, DMW Fatigue Tests, System Level Thermal-Mechanical-Stress Analysis

The long-term operation (LTO) of nuclear power plant (NPP) beyond their original design life of 40 years, can lead to more material damage associated with cyclic fatigue under thermal-mechanical loading cycles and associated long-term exposure of reactor material to the deleterious reactor-coolant environments. However, under this LTO condition the reactor components can still safely operate but may require more frequent Nondestructive Evaluation (NDE) of reactor components. Frequent NDE requirement may lead to frequent shutdown of the NPP. This in turn can lead to power outage and additional NDE-inspection-cost related economic loss. The economic loss can be minimized by reducing uncertainty in life estimation of safety-critical pressure boundary components and by implementing more digital approach such as by using upcoming digital-twin (DT) technology for predicting the structural states (e.g., time and location dependent inside/outside thickness temperature, stress, strain, plastic deformation, etc.) and associated fatigue life of a component in real time. Towards this goal Argonne National Laboratory (ANL) with the sponsorship of DOE Light Water Reactor Sustainability (LWRS) program is working on the development of a DT framework that can be used for real time environmental fatigue prediction of reactor components. The DT framework is based on limited experiment-data, Artificial-intelligence (AI) – Machine-Learning (ML) - Deep-Learning (DL) based techniques and Multiphysics-computational-mechanics such as finite element (FE) based modeling tools. Towards this overall goal, following are some of the major contributions made during the FY21: 1) Multiple 82/182 dissimilar metal weld (DMW) specimens (both solid-weld and joint-weld representing the actual reactor multi-metal nozzles) were fatigue tested. The resulting fatigue lives were compared to the NUREG-6909 based best-fit and design fatigue curves. Additionally, the results of 52/152 DMW fatigue specimens (which were recently tested at Republic of Korea under the sponsorship of International Nuclear Energy Research Initiative - INERI program) were compared to the NUREG-6909 based best-fit and design fatigue curves. From the comparison of 82/182 and 52/152 DMW test data with NUREG-6909 best-fit curve, most of the reported test data fall way away from the NUREG-6909 suggested best-fit or mean curve. The NUREG-6909 suggested best-fit curve is the best-fit curve of austenitic stainless steel and due to lack of enough data on Nickel-based welds, this is currently being used for predicting the life of Nickel-alloy-based welded components. However, the above observation may require higher scaling factor (e.g., ASME suggested factor of 20 on cycles rather than the current NUREG-6909 suggested factor of 12 on cycles) for scaling the austenitic-stainless-steel best-fit-curve for estimating the design or safe-life of a welded component. Accordingly, for example, if a DMW component experience a strain amplitude of 0.6% the PWR-water life of the component would be 52 cycles instead of 85 cycles. However, more DMW tests are required to further ascertain the above-mentioned observations. 2) A system level CAD and finite element model were developed which consists of reactor pressure vessel (RPV), part of steam generator (SG), part of pressurizer (PRZ), hot leg (HL), and surge line (SL). This is with detailed nozzle geometry and thermal-mechanical material properties of different metals to simulate realistic thermal-mechanical stress under connected system global thermal-mechanical boundary conditions. 3) Different system level heat transfer analyses were performed with estimation of relevant heat transfer coefficients. The resulting data were used in subsequent system level thermal-mechanical stress analysis and for generating spatial-temporal training and validation data for a system level digital-twin based temperature predictor. Transient heat transfer analyses were performed considering thermal boundary condition under design-basis (DB) loading and EDF (Électricité de France) data-based grid-load-following (EDF-GLF) loading cycles. 4) System level thermal-mechanical stress analysis was performed for identifying damage-prone hotspots and for future extension of the model for cyclic state prediction. From the system-level model simulation under DB loading cycle it is found that HL and the SL nozzle that connect to the HL can experience significant stress and strain and could be one of the weakest links in the overall reactor coolant system (RCS). 5) An AI/ML based DT model was developed for multi-time-series temperature prediction at any inside/outside thickness locations of PWR pressure boundary components. This is by using Recurrent-neural-network (RNN) and keras machine learning libraries. The RNN model was validated against two laboratory test-based data sets with one obtained through ANL’s in-air fatigue test system and other through PWR-water test loop. The experimentally validated DT model further validated against FE model results to predict thermal scarification related spatialtemporal temperatures at random locations of a component. The well validated DT model was then used for demonstrating spatial-temporal temperature prediction under 100+ years of reactor operation subjected to combined DB, EDF-GLF and randomized grid-load-following (RANDOMGLF) loading Cycles. The expert-elicitation DT model framework was developed assuming field/input/process measurements can be available from a few existing plant sensors and can readily be used by the NPP operators. The above temperature prediction model will feed to the next-step stress analysis model based on which the life of a component can be predicted in realtime, which is one of our future works.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1992

1400 and 1500 K Compressive Creep Properties of an NiAl/AlN Composite

Compressive creep properties of an NiAl/AlN(p) composite produced by a reaction milling process were investigated at 1400 and 1500 K and at slow strain rates, to investigate the relative strength of this composite at high temperatures, and to determine if the consolidation technique affects the 1400 K creep properties. Results indicate that the stress exponent of the NiAl/AlN(p) composite was similar to that for unreinforced NiAl. However, the activation energy for the composite was found to be more than twice that measured in the unreinforced matrix. Oxidation did not affect the composite at 1400 K, but a significant attack was observed in a sample subjected to fast deformation at 1500 K.

Whittenberger, J. D.↗
DOE OSTI2024

Sister Rod Destructive Examinations (FY23) Appendix F: Cyclic Integrated Reversible-Bending Fatigue Tests

Spent nuclear fuel (SNF) assemblies must be shipped to other sites for processing and disposal. During shipment, the fuel is typically oriented horizontally, and the fuel rods are subject to periodic alternating loads related to the movement of the vehicle that results in alternating bending of the SNF fuel rods. The number of bending cycles is related to the length of the shipping route, with longer routes producing more cycles. Because it is well-known that cyclic loads can produce failures, even when the stress and strain imposed are below the material’s yield point, investigation of the SNF’s fatigue behavior is prudent. This report discusses the results of fatigue testing conducted at Oak Ridge National Laboratory (ORNL) using the Cyclic Integrated Reversible-Bending Fatigue Tester (CIRFT) created by the US Department of Energy Office of Nuclear Energy for the High Burnup Spent Fuel Data Project and its sister rods. The results of the tests are compared with results obtained over the last decade using the same CIRFT for the US Nuclear Regulatory Commission (NRC).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Understanding Uncertainty in Microstructure Evolution and Constitutive Properties in Additive Process Modeling

Coupled process–microstructure–property modeling, and understanding the sources of uncertainty and their propagation toward error in part property prediction, are key steps toward full utilization of additive manufacturing (AM) for predictable quality part development. The OpenFOAM model for process conditions, the ExaCA model for as-solidified grain structure, and the ExaConstit model for constitutive mechanical properties are used as part of the ExaAM modeling framework to examine a few of the various sources of uncertainty in the modeling workflow. In addition to “random” uncertainty (due to random number generation in the orientations and locations of grains present), the heterogeneous nucleation density N 0 and the mean substrate grain spacing S 0 are varied to examine their impact of grain area development as a function of build height in the simulated microstructure. While mean grain area after 1 mm of build is found to be sensitive to N 0 and S 0 , particularly at small N 0 and large S 0 (despite some convergence toward similar values), the resulting grain shapes and overall textures develop in a reasonably similar manner. As a result of these similar textures, ExaConstit simulation using ExaCA representative volume elements (RVEs) from various permutations of N 0 , S 0 , and location within the build resulted in similar yield stress, stress–strain curve shape, and stress triaxiality distributions. It is concluded that for this particular material and scan pattern, 15 layers is sufficient for ExaCA texture and ExaConstit predicted properties to become relatively independent of additional layer simulation, provided that reasonable estimates for N 0 and S 0 are used. However, additional layers of ExaCA will need to be run to obtain mean grain areas independent of build height and baseplate structure.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Shear deformation plate continua of large double-layered space structures

The energy equivalence to construct equivalent continua for the actual lattice structure is described. An energy equivalent continuum is defined as that which has the same amount of strain and kinetic energies stored in it as the original lattice structure when both are subjected to the same loading conditions. The equivalent continuum is characterized by its strain and kinetic energies from which the constitutive relations and the equations of motion can be derived. A simple method to model large rigid jointed lattices as continuous media with couple stresses is presented. The transition from the discrete system to the continuous medium is achieved by expanding the displacements and the rotations of the nodal points in a Taylor series about a suitable chosen origin. Basic kinematic assumptions are introduced to insure that the assumptions used in deriving the governing equations of the modeled continuum are satisfied. The number of terms retained in the Taylor series expansion will depend upon the properties to be evaluated. This implies that the kind of continuum needed to model from the discrete lattice, before the actual properties are derived was to predetermined.

Hefzy, M. S.↗
NASA NTRS1981

A Single-Expression Formula for Inverting Strain-Life and Stress-Strain Relationships

Starting with the basic fatigue lift formula, an inversion formula is derived. The inversion formula is valid over the entire life range of engineering interest for all materials examined. Conformity between the two equations is extremely close, suitable for all engineering problems. The approach used to invert the life relation is also suitable for the inversion of other formulas involving the sum of two power-law terms.

Manson, S. S.↗
DOE OSTI2024

Neutron and photon in-vivo materials characterization at the evolving plasma-material interface in plasma-burning fusion environments

This report includes the first set of MOSS data collected in the context of fusion-relevant plasma-material interactions. The data collected show that the surface curvature and stress irreversibly change during low-energy deuterium ion irradiation. This holds for both tungsten thin films, bulk tungsten, and single-crystal silicon. Meanwhile, argon irradiation did not produce this same effect. The surface morphology of the exposed samples was virtually identical regardless of gas species. There were, however, some surface chemistry differences between argon-exposed and deuterium-exposed samples, mainly concerning the oxide states of the W4f XPS signal. The changes in surface chemistry could be due to hydrogen dynamics on the tungsten surfaces, though further investigation will be required with in-situ irradiation and analysis. The changes in surface stress could also not be explained by temperature alone, as the temperature reached during the ion irradiations was too low to induce non-elastic strain. With this in mind, deuterium retention in the materials could change the surface curvature and stress. To confirm this, further work will be needed to quantify the change in curvature and relate it to ion flux, fluence, and deuterium retention. Unfortunately, when this work was conducted, the experimental setup for thermal desorption spectroscopy was not operating correctly, and that work could not be included in this dissertation. Despite experimental difficulties, this work shows the first results of a promising technique to study deuterium dynamics in tungsten-based reconstituted thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Predicting long-term stress relaxation on Alloy 709 using the crystal plasticity finite element method

This report summarizes the results of physics-based, crystal plasticity simulations for the long-term stress relaxation behavior of Alloy 709. The purpose of the study was to provide insight into five key questions related to long-term behavior in high temperatures materials which are difficult or impossible to answer experimentally: (1) is there a threshold stress for long-term relaxation? (2) is there strain threshold for relaxation damage, below which significant damage does not accumulate? (3) does damage continue to accumulate as the material relaxes or will damage accumulation plateau under some loading conditions? (4) does stress relaxation loading inevitably lead to failure? and (5) which, if any, engineering models for relaxation damage accumulation reasonably match the simulation results? The report summarizes the numerical simulations used to address these five questions and provides at least partial answers to each question.

36 MATERIALS SCIENCE↗