Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “composition and microstructure”

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 919 records · Page 51

Heat affected zone liquation cracking evaluation on FeMnAl alloys

FeMnAl steels are currently generating a lot of interest with potential applications for structural parts in armored vehicles due to their lower density and outstanding mechanical properties. Despite the extensive mechanical performance and heat treatment exploration of this alloy class, further weldability investigation is required for future large-scale deployment. In the present study, the liquation cracking of four heats of cast FeMnAl alloys was investigated by the spot-Varestraint technique. The study focuses primarily on the effect of the major elements of the FeMnAl system: C, Mn and Al. Optical and electron microscopy were employed to investigate the microstructural features, and CALPHAD was employed to aid the discussion regarding the alloy’s composition differences and their effect on the liquation cracking susceptibility. The study was able to identify that compositions with the higher Mn, C, and lower Al presented the highest liquation cracking susceptibility. Conversely, composition presenting lower Mn, C, and Al showed the most resistant behavior. Furthermore, lower Al content promoted a fully-γ microstructure at low temperatures, which encouraged the appearance of longer cracks as a γ-matrix is more susceptible to HAZ cracking than a fully ferritic (α) or duplex (α + γ) microstructure.

36 MATERIALS SCIENCE↗

A Novel Ultrasonic Method for Characterizing Microstructural Gradients in Tubular Structures

Ultrasonic velocity and time-of-flight (TOF) imaging that uses back surface reflections to gage volumetric material quality is highly suited for quantitative characterization of microstructural gradients including those due to pore fraction, density, fiber fraction, and chemical composition variations. However, a weakness of conventional pulse echo ultrasonic velocity and TOF imaging is that the image shows the effects of thickness as well as microstructural variations, unless the part is uniformly thick. This limits this imaging method's usefulness in practical applications. Prior studies have described a pulse echo TOF based ultrasonic imaging method that requires using a single transducer in combination with a reflector plate placed behind samples which, eliminates the effect of thickness variation in the image. In those studies, this method was successful at isolating ultrasonic variations due to material in plate like samples of silicon nitride, metal matrix composite, and polymer matrix composite. In this study, the method is engineered for inspection of more complex shaped structures- those having (hollow) tubular or curved geometry. The experimental inspection technique and results are described as applied to a polymer matrix composite "proof of concept" tube that contains machined patches of various depths and an as manufactured monolithic silicon nitride ceramic tube that might be used in "real world" applications.

Roth, Don J.↗

Microstructural Evolution and Mechanical Properties of LP-DED NASA HR-1 – A Hydrogen Resistant AM Superalloy for Space Propulsion Applications

The National Aeronautics and Space Administration (NASA) has actively pursued metal additive manufacturing (AM) technologies for spaceflight applications since the late 2000s. AM offers transformative advantages in cost, schedule, part consolidation, and design flexibility. Among the various AM techniques, laser powder directed energy deposition (LP-DED) is particularly well suited for fabricating complex geometries with fine feature resolution. In propulsion systems that utilize high-pressure gaseous hydrogen—such as liquid hydrogen rocket engines—hydrogen environment embrittlement (HEE) presents a serious threat to material performance 1,2 . Mechanical property degradation under these conditions can compromise component reliability, especially under cyclic loading. To address this challenge, NASA developed NASA HR-1 (Hydrogen Resistant-1) as a solution for liquid rocket engine components operating in hydrogen-rich environments, using the LP-DED technique 3-9 . A key component in a liquid rocket engine is the exhaust nozzle, which is typically regeneratively cooled (regen) due to the high heat flux. NASA HR-1 was specifically developed for regen nozzle applications using hydrogen as a propellant, providing resistance to HEE, a critical issue for many materials. The AM version of NASA HR-1 was also formulated to achieve high ultimate tensile strength, along with high yield strength and ductility in this environment 5,6 . Low-cycle fatigue (LCF) is another important consideration in nozzle design, as components are expected to endure multiple starts and missions. Additionally, the LP-DED version of the alloy exhibits improved thermal conductivity compared to its wrought counterpart, which benefits nozzle cooling. Overall, NASA HR-1 offers an excellent balance of high strength, HEE resistance, LCF performance, thermal conductivity, and ductility to meet the demanding requirements of channel-cooled nozzles and other components used with hydrogen and other propellants. The LP-DED–processed NASA HR-1 requires several post-processing heat treatment steps to achieve the material properties desirable for its intended application 6 . These steps include stress relief, homogenization, solution annealing, and aging for precipitation hardening. The stress relief treatment mitigates residual stresses accumulated during the LP-DED process and minimizes the potential for distortion. Homogenization, a common step for AM materials, reduces elemental segregation and promotes recrystallization to develop a more equiaxed grain structure. The subsequent solution annealing treatment heats the part to a solid solution temperature to dissolve the undesirable η-phase that forms during cooling from homogenization, followed by rapid cooling to retain an η-phase–free microstructure. Finally, aging promotes precipitation of the strengthening γ′ phase in the alloy. The integration of compositional design and optimized thermal processing enables high-quality LP-DED NASA HR-1 components with excellent microstructural and mechanical stability. Improved chemical and microstructure homogeneity enhances ductility and fatigue resistance—both critical for safe and reliable operation in high-pressure hydrogen environments. NASA has successfully fabricated and hot-fire tested multiple subscale and full-scale channel wall nozzles using LP-DED NASA HR-1 5,6, 9-14 . These efforts included process refinements to support thin-wall construction and various channel geometries. Throughout development, several key observations emerged. After homogenization, the as-built columnar grain structure transforms into a fully equiaxed microstructure. However, subsequent treatments—such as solution annealing and aging—result in changes that are more difficult to track. The grain structure remains largely unchanged, and the γ′ precipitates, typically 5–10 nm in diameter, are beyond the resolution of scanning electron microscopy (SEM). While transmission electron microscopy (TEM) can resolve these fine precipitates, TEM sample preparation is time-consuming and difficult for LP-DED material. As an alternative, differential scanning calorimetry (DSC) offers a useful, qualitative approach to monitor precipitate evolution throughout different stages of heat treatment. The overall goal is to improve the understanding of how heat treatment affects the microstructure and mechanical performance of LP-DED NASA HR-1. This paper presents heat treatment design considerations, microstructural characterization, mechanical testing – including tensile and LCF testing in both air and hydrogen environments.

Superalloy↗

Macromolecular gelatin properties affect fibrin microarchitecture and tumor spheroid behavior in fibrin-gelatin gels

The biophysical properties of extracellular matrices (ECM) are known to regulate cell behavior, however decoupling cell behavior changes due to the relative contributions of material microstructure versus biomechanics or nutrient permeability remains challenging, especially within complex, multi-material matrices. We developed four gelatin-fibrin interpenetrating network (IPN) formulations which are identical in composition but possess variable gelatin molecular weight distributions, and display differences in microstructure, biomechanics, and diffusivity. In this work we interrogate the response of multicellular tumor spheroids to these IPN formulations and found that a high stiffness, gelatin-network dominated IPNs impeded remodeling and invasion of multicellular tumor spheroids; whereas relatively lower stiffness, fibrin-network dominated IPNs permitted protease-dependent remodeling and spheroid invasion. Cell proliferation correlated to nutrient diffusivity across tested IPN formulations. These findings demonstrate the complexity of ECM IPNs, relative to single polymer matrices, and highlight that cell response does not derive from a single aspect of the ECM, but rather from the interplay of multiple biomechanical properties. The methodology developed in this study represents a framework for future studies which aim to characterize cellular phenotypic responses to biophysical cues present within complex, multi-material matrices.

36 MATERIALS SCIENCE↗

High-temperature ordered intermetallic alloys III; Proceedings of the Third Symposium, Boston, MA, Nov. 29-Dec. 1, 1988

The present conference on high-temperature ordered intermetallic alloys discusses alloy theory and phase stability, defects and microstructures, mechanical behavior, alloy design and material processing, multiphase and composite materials, nickel aluminides, titanium aluminides, and other ordered intermetallics. Attention is given to compositional effects on processing and properties of nickel aluminides, dynamic compaction and hot-isostatic-pressing of nickel aluminides, improvement of elevated-temperature material properties in Ni3Al, the effect of microalloying B on the high-temperature mechanical properties of Ti3Al, the effects of structure on creep of Ti-53.4 mol pct Al intermetallics, the deformation microstructure in Ni3Si polycrystals strained over the range of temperature of flow stress anomaly, and the microstructure and mechanical properties of dual phase alloys consisting of the intermetallic phases.

Liu, C. T.↗

Fatigue crack propagation behavior of the Ni-Ni3Cb eutectic composite

The experimental procedures used in the investigation are considered, taking into account solidification and heat treatment, fatigue testing, questions of test specimen and loading fixture design verification, metallographic techniques, and electron microscopy techniques. Attention is given to the characteristics of the microstructure, aspects of fatigue crack propagation (FCP) in the Ni-Ni3Nb composite, the effect of test variables on FCP, the effect of metallurgical variables on FCP, and fatigue fracture mechanisms. Fractographic and metallographic observations are discussed.

Mills, W. J.↗

Surface Finishing and Coating Parameters Impact on Additively Manufactured Binder-Jetted Steel–Bronze Composites

In this paper, electroless nickel plating is explored for the protection of binder-jetting-based additively manufactured (AM) composite materials. Electroless nickel plating was attempted on binder-jetted composites composed of stainless steel and bronze, resulting in differences in the physicochemical properties. We investigated the impact of surface finishing, plating solution chemistry, and plating parameters to attain a wide range of surface morphologies and roughness levels. We employed the Keyence microscope to quantitatively evaluate dramatically different surface properties before and after the coating of AM composites. Scanning electron microscopy revealed a wide range of microstructural properties in relation to each combination of surface finishing and coating parameters. We studied chempolishing, plasma cleaning, and organic cleaning as the surface preparation methods prior to coating. We found that surface preparation dictated the surface roughness. Taguchi statistical analysis was performed to investigate the relative strength of experimental factors and interconnectedness among process parameters to attain optimum coating qualities. The quantitative impacts of phosphorous level, temperature, surface preparation, and time factor on the roughness of the nickel-plated surface were 17.95%, 8.2%, 50.02%, and 13.21%, respectively. On the other hand, the quantitative impacts of phosphorous level, temperature, surface preparation, and time factor on the thickness of nickel plating were 35.12%, 41.40%, 3.87%, and 18.24%, respectively. The optimum combination of the factors’ level projected the lowest roughness of Ra at 7.76 µm. The optimum combination of the factors’ level projected the maximum achievable thickness of ~149 µm. This paper provides insights into coating process for overcoming the sensitivity of AM composites in hazardous application spaces via robust coating.

36 MATERIALS SCIENCE↗

Thermoelectric Properties in the TiO2/SnO2 System

Nanotechnology has provided a new interest in thermoelectric technology. A thermodynamically driven process is one approach in achieving nanostructures in bulk materials. TiO2/SnO2 system exhibits a large spinodal region with exceptional stable phase separated microstructures up to 1400 C. Fabricated TiO2/SnO2 nanocomposites exhibit n-type behavior with Seebeck coefficients greater than -300 .V/K. Composites exhibit good thermal conductance in the range of 7 to 1 W/mK. Dopant additions have not achieved high electrical conductivity (<1000 S/m). Formation of oxygen deficient composites, TixSn1-xO2-y, can change the electrical conductivity by four orders of magnitude. Achieving higher thermoelectric ZT by oxygen deficiency is being explored. Seebeck coeffcient, thermal conductivity, electrical conductance and microstructure will be discussed in relation to composition and doping.

Dynys, F.↗

Synthesis of a 316L stainless steel-copper composite by laser melting

This paper demonstrates a laser melting method to produce a metal matrix composite with a higher thermal conductivity compared to 316 L stainless steel using a Jettable Cu ink. A novel canister was used to control the oxygen level during laser processing to emulate a hybrid inkjet-powder bed fusion technique for doping the stainless steel powder bed with Cu precursor ink. Here, the thermal conductivity of the composite was found to be 187% higher than that of 316 L stainless steel, while the microhardness decreased by 39%. Microstructural results show opportunities to further enhance the thermal conductivity and mechanical properties of the composite by eliminating interfacial cracking.

36 MATERIALS SCIENCE↗

Diversifying Composition Leads to Hierarchical Composites with Design Flexibility and Structural Fidelity

Although significant progress has been made in the self-assembly of nanostructures, present successes heavily rely on precision in building block design, composition, and pair interactions. These requirements fundamentally limit our ability to synthesize macroscopic materials where the likelihood of impurity inclusion escalates and, more importantly, to access molecular-to-nanoscopic-to-microscopic-to-macroscopic hierarchies, since the types and compositions of building blocks vary at each stage. Inspired by biological blends and high-entropy alloys, we hypothesize that diversifying the blend’s composition can overcome these limitations. Increasing the number of components increases mixing entropy, leading to the dispersion of different components and, as a result, enhances interphase miscibility, weakens the dependence on specific pair interactions, and enables long-range cooperativity. This hypothesis is validated in complex blends containing small molecules, block copolymer-based supramolecules, and nanoparticles/colloidal particles. Hierarchically structured composites can be obtained with formulation flexibility in the filler selection and blend composition. It is worth noting that, by adding small molecules, we can solve the size constraint that plagues traditional block copolymer/nanoparticle blends. Detailed characterization and simulation further confirm that each component is distributed to locally mediate unfavorable interactions, cooperatively mitigate composition fluctuations, and retain structural fidelity. Furthermore, the blends have sufficient mobility to access tunable microstructures without compromising the order of the nanostructure. Besides establishing a kinetically viable pathway to release current constraints in the composite design and to navigate uncertainties during structure formation over multiple length scales, the present study demonstrates that entropy-driven behaviors can be realized in systems beyond high-entropy alloys despite inherent differences between metal alloys and organic/inorganic hybrids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MS212--A Homogeneous Sputtered Solid Lubricant Coating for Use to 800 C

Composite coatings containing chromium carbide, stable fluorides and silver were prepared by magnetron sputtering. The microstructure of the coatings is very homogeneous compared to that of plasma sprayed and sintered versions of the same chemical composition. Friction and wear of MS212-coated and baseline uncoated aluminum and Inconel X-750 are compared. At room temperature, the friction and wear of coated aluminum is dramatically better compared to the baseline. The acceptable load is limited by deformation of the soft aluminum substrate. In the case of the nickel alloy, lower friction is observed for the coated alloy at all temperatures up to the maximum test temperature of 800 C. Pin wear factors for sliding against the coated alloy are lower than the baseline at room temperature and 350 C, and comparable to baseline wear at higher test temperatures. Low baseline wear at high temperatures is due to the lubricious nature of the natural oxides formed on nickel-chromium alloys in a hot, oxidizing atmosphere. No load limit was found for coated Inconel X-750 at loads up to five times the load limit for coated aluminum.

Sliney, Harold E.↗

Microscopic Scale Simulation of the Ablation of Fibrous Materials

Ablation by oxidation of carbon-fiber preforms impregnated in carbonized phenolic matrix is modeled at microscopic scale. Direct numerical simulations show that the carbonized phenolic matrix ablates in volume leaving the carbon fibers exposed. This is due to the fact that the reactivity of carbonized phenolic is higher than the reactivity of carbon fibers. After the matrix is depleted, the fibers ablate showing progressive reduction of their diameter. The overall material recession occurs when the fibers are consumed. Two materials with the same carbon-fiber preform, density and chemical composition, but with different matrix distributions are studied. These studies show that at moderate temperatures (< 1000 K) the microstructure of the material influences its recession rate; a fact that is not captured by current models that are based on chemical composition only. Surprisingly, the response of these impregnated-fiber materials is weakly dependent on the microstructure at very high temperatures (e.g., Stardust peak heating conditions: 3360K).

Lachaud, Jean Romain↗

Mesoscale modeling of microstructure-dependent thermal conductivity in U-Zr fuels

In uranium-zirconium (U-Zr) based metallic fuels, different phases can form at different compositions and temperatures. Typically, lamellar δ-UZr 2 and α-U phases are the dominant microstructures in U-rich U-Zr alloys at temperatures below 880 K. In this work, a finite element method based mesoscale modeling technique is used to calculate the effective thermal conductivities of such heterogeneous microstructures, using the thermal conductivities of two individual phases and their interphase thermal resistance (Kapitza resistance) as input parameters. The Kapitza resistance between δ-UZr 2 and α-U is determined at different temperatures, which shows an approximately T 3 dependence in the temperature range between 300 and 800 K. In addition, the Kapitza resistance exhibits a strong dependence on the aspect ratio of the δ-UZr 2 phase. Further, an analytical model is therefore developed to quantify the effects of both temperature and δ-UZr 2 aspect ratio on the Kapitza resistance. Using this newly developed Kapitza resistance model, the effective thermal conductivities of a number of δ-UZr 2 + α-U heterogeneous microstructures in U-Zr alloys, including non-lamellar microstructures, can be estimated accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Cooling rate and dendrite spacing control in direct metal deposition printed Cu-Fe alloys

Over the last several decades, additive manufacturing (AM) has been primarily used for rapid prototyping or to create novel geometries that would be difficult or impossible to create by normal manufacturing methods. More recently, research has been focused on expanding the list of materials that can be made through additive manufacturing, opening a greater range of material properties for this manufacturing method. Due to the unusual conditions during AM, including the high cooling rates and voxel by voxel method of production, AM parts often have anisotropic material microstructures and properties. In this investigation, the laser power, composition, and nozzle head speed during direct metal deposition of copper-iron alloys was varied to understand how the grain structure within the printed parts could be changed and controlled. The resulting dendrite spacing was measured and compared to calculated cooling rates from Gaussian beams on flat plates under similar material and laser properties, which resulted in a cooling rate to dendrite spacing relationship following an inverse square root, as is found in other dendritic systems [Young and Kerkwood, Metall. Trans. A 6, 197–205 (1975)]. Thus, it is demonstrated that in the Cu-Fe system, dendrite spacing can be controlled through manipulation of printing parameters.

36 MATERIALS SCIENCE↗

Characterization of Ternary NiTiPt High-Temperature Shape Memory Alloys

Pt additions substituted for Ni in NiTi alloys are known to increase the transformation temperature of the alloy but only at fairly high Pt levels. However, until now only ternary compositions with a very specific stoichiometry, Ni50-xPtxTi50, have been investigated and then only to very limited extent. In order to learn about this potential high-temperature shape memory alloy system, a series of over twenty alloys along and on either side of a line of constant stoichiometry between NiTi and TiPt were arc melted, homogenized, and characterized in terms of their microstructure, transformation temperatures, and hardness. The resulting microstructures were examined by scanning electron microscopy and the phase compositions quantified by energy dispersive spectroscopy."Stoichiometric" compositions along a line of constant stoichiometry between NiTi to TiPt were essentially single phase but by any deviations from a stoichiometry of (Ni,Pt)50Ti50 resulted in the presence of at least two different intermetallic phases, depending on the overall composition of the alloy. Essentially all alloys, whether single or two-phase, still under went a martensitic transformation. It was found that the transformation temperatures were depressed with initial Pt additions but at levels greater than 10 at.% the transformation temperature increased linearly with Pt content. Also, the transformation temperatures were relatively insensitive to alloy stoichiometry within the range of alloys examined. Finally, the dependence of hardness on Pt content for a series of Ni50-xPtxTi50 alloys showed solution softening at low Pt levels, while hardening was observed in ternary alloys containing more than about 10 at.% Pt. On either side of these "stoichiometric" compositions, hardness was also found to increase significantly.

Rios, Orlando↗

Quantifying electron transport in aggregated colloidal suspensions in the strong flow regime

Electron transport in complex fluids, biology, and soft matter is a valuable characteristic in processes ranging from redox reactions to electrochemical energy storage. These processes often employ conductor–insulator composites in which electron transport properties are fundamentally linked to the microstructure and dynamics of the conductive phase. While microstructure and dynamics are well recognized as key determinants of the electrical properties, a unified description of their effect has yet to be determined, especially under flowing conditions. In this work, the conductivity and shear viscosity are measured for conductive colloidal suspensions to build a unified description by exploiting both recent quantification of the effect of flow-induced dynamics on electron transport and well-established relationships between electrical properties, microstructure, and flow. These model suspensions consist of conductive carbon black (CB) particles dispersed in fluids of varying viscosities and dielectric constants. In a stable, well-characterized shear rate regime where all suspensions undergo self-similar agglomerate breakup, competing relationships between conductivity and shear rate were observed. To account for the role of variable agglomerate size, equivalent microstructural states were identified using a dimensionless fluid Mason number, Mn f , which allowed for isolation of the role of dynamics on the flow-induced electron transport rate. At equivalent microstructural states, shear-enhanced particle–particle collisions are found to dominate the electron transport rate. This work rationalizes seemingly contradictory experimental observations in literature concerning the shear-dependent electrical properties of CB suspensions and can be extended to other flowing composite systems.

Science & Technology - Other Topics↗

High-Performance Scanning Acousto-Ultrasonic System

A high-performance scanning acousto-ultrasonic system, now undergoing development, is designed to afford enhanced capabilities for imaging microstructural features, including flaws, inside plate specimens of materials. The system is expected to be especially helpful in analyzing defects that contribute to failures in polymer- and ceramic-matrix composite materials, which are difficult to characterize by conventional scanning ultrasonic techniques and other conventional nondestructive testing techniques. Selected aspects of the acousto-ultrasonic method have been described in several NASA Tech Briefs articles in recent years. Summarizing briefly: The acousto-ultrasonic method involves the use of an apparatus like the one depicted in the figure (or an apparatus of similar functionality). Pulses are excited at one location on a surface of a plate specimen by use of a broadband transmitting ultrasonic transducer. The stress waves associated with these pulses propagate along the specimen to a receiving transducer at a different location on the same surface. Along the way, the stress waves interact with the microstructure and flaws present between the transducers. The received signal is analyzed to evaluate the microstructure and flaws. The specific variant of the acousto-ultrasonic method implemented in the present developmental system goes beyond the basic principle described above to include the following major additional features: Computer-controlled motorized translation stages are used to automatically position the transducers at specified locations. Scanning is performed in the sense that the measurement, data-acquisition, and data-analysis processes are repeated at different specified transducer locations in an array that spans the specimen surface (or a specified portion of the surface). A pneumatic actuator with a load cell is used to apply a controlled contact force. In analyzing the measurement data for each pair of transducer locations in the scan, the total (multimode) acousto-ultrasonic response of the specimen is utilized. The analysis is performed by custom software that extracts parameters of signals in the time and frequency domains. The computer hardware and software provide both real-time and postscan processing and display options. For example, oscilloscope displays of waveforms and power spectral densities are available in real time. Images can be computed while scanning continues. Signals can be digitally preprocessed and/or post-processed by filtering, windowing, time-segmenting, and running-waveform-averaging algorithms. In addition, the software affords options for off-line simulation of the waveform-data-acquisition and scanning processes. In tests, the system has been shown to be capable of characterizing microstructural changes and defects in SiC/SiC and C/SiC ceramic-matrix composites. Delaminations, variations in density, microstructural changes attributable to infiltration by silicon, and crack-space indications (defined in the next sentence) have been revealed in images formed from several time- and frequency-domain parameters of scanning acousto-ultrasonic signals. The crack-space indications were image features that were not revealed by other nondestructive testing methods and are so named because they turned out to mark locations where cracking eventually occurred.

Roth, Don↗

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗