Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite strains”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Alloying effects on deformation induced microstructure evolution in copper

In this work, we investigated the effects of alloying elements on plastic deformation and microstructure evolution in polycrystalline copper (Cu) and Cu alloyed with 1 wt. % lead (Cu-1 %Pb). These materials were selected due to the size mismatch between Cu and Pb, with the latter forming precipitates at grain boundaries. Multi-modal characterization techniques, including neutron diffraction, electron backscatter diffraction (EBSD), and transmission electron microscopy (TEM), along with finite element simulations were employed to study the deformation behavior across multiple length scales. While both Cu and Cu-1%Pb exhibited similar macroscale response and final deformation textures, both dislocation line profile analysis and TEM revealed increased dislocation density in deformed Cu-1%Pb specimens. The presence of lead precipitates also significantly affected local plastic deformation during compression, with their influence diminishing with increasing strain. These results demonstrate the complex relationships between alloying elements, plastic deformation, microstructural evolution, and material behavior under load. The insights gained from this multi-scale and multi-technique approach contribute to the fundamental understanding of microstructural evolution in immiscible alloys and are valuable for tailoring the properties of structural materials for specific engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Characterization of the Fusion Zone in Laser Beam Welds on SS304L

Designing a welded joint requires consideration of the welding process, component/assembly structure and weld geometry, material composition, properties, and desired joint performance. The designer must also specify a set of acceptance criteria used to evaluate the welds. This study was conducted to support the design team responsible for specifying the spot weld profile by investigating the relationship between weld process parameters and the weld geometry and properties to further inform design optimization. The subjects of this study are two thin stainless steel 304L (SS304L) components joined by Laser Beam Welding (LBW). The components are the Housing, which contains flexible circuitry, and the Cover that protects the circuitry during handling, installation, and under service loads. Figure 1 illustrates coupon versions of the Cover and Housing used to simulate the interface between these components for weld testing. These coupon representations of assembly components were designed to match the material and relevant geometry related to the welded joint while other features not relevant to the welded joint are excluded.The objective of this study was to characterize the sensitivity of the weld geometry and mechanical properties to the input weld parameters to inform further Finite Element Analysis (FEA) efforts in balancing joint strength against assembly deformation caused by the welding process. Characterization was accomplished by varying laser welding parameters around the current design values and measuring weld geometry on the cross-sections. Nanoindentation with a spherical and flat-ended tip was used on polished section cuts of the weld to obtain direct estimates of the stress-strain response of the Fusion Zone (FZ) for comparison to the Base Metal (BM).

42 ENGINEERING↗
DOE OSTI2024

Validation of Hole-Drilling Residual Stress Measurements in Workpieces of Various Thickness

A recent revision to the ASTM E837 standard for near-surface residual stress measurement by the hole-drilling method describes a new thickness-dependent stress calculation procedure applicable to “thin” and “intermediate” workpieces for which strain versus depth response depends on workpiece thickness. This new calculation procedure differs from that of the prior standard, which applies only to thick workpieces with strain versus depth response independent of thickness. Herein we assess the new calculation procedures by performing hole-drilling residual stress measurements in samples with a range of thickness. Near-surface residual stress is measured in a thick aluminum plate containing near-surface residual stress from a uniform shot peening treatment, and in samples of different thickness removed from the plate at the peened surface. A finite element (FE) model is used to assess consistency between measured residual stress across the range of sample thickness. Measured residual stress varies with sample thickness, with thinner samples exhibiting smaller near-surface compressive stress and a larger gradient of subsurface stress. These trends are consistent with both observed bending (curvature) of the removed samples and the trend in FE-calculated expected residual stress. The measured and expected residual stresses are in good agreement for samples of intermediate thickness, but the agreement decreases with sample thickness. Measured residual stress is invariant with gage circle diameter. The new thickness-dependent stress calculation procedure for hole-drilling provides meaningful improvement compared to thick-workpiece calculations.

42 ENGINEERING↗
DOE OSTI2020

Toward large-scale simulation of residual stress and distortion in wire and arc additive manufacturing

In this study, we aim to advance the structural analysis of wire and arc additive manufacturing (WAAM) by considering the thermomechanical features inherent in direct energy deposition. Simulation approaches including the iterative substructure method (ISM), dynamic mesh refining method (DMRM), and graphics processing unit (GPU) based explicit finite element method (FEM) were developed for accelerating additive manufacturing stress analysis that is very time consuming by conventional numerical methods. ISM and DMRM take advantage of the strong nonlinearity phenomenon near the moving heat source by reducing the global iterations and refining the local mesh, respectively. In addition to the use of GPUs, the explicit FEM is accelerated by a time scaling technique based on the inherent strain concept. The residual stress and distortion of two large builds were analyzed, showing very consistent numerical results and good agreement with experiments. Compared with the commercial software Abaqus, the novel approaches reduced the computational cost substantially without compromising accuracy. The GPU code showed the highest computational efficiency (30~70×), while DMRM and ISM had acceleration factors of 9× and 3×, respectively. Finally, such high-fidelity modeling approaches will be very useful for building up a digital twin of WAAM to reduce development time and cost.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimentally Validated Computational Modeling of Creep and Creep-Cracking for Nuclear Concrete Structures

In a Nuclear Power Plant, one of the most important components is the concrete nuclear reactor cavity, which serves both a structural and protective function as the biological radiation shield. Given that creep has been identified as a major knowledge gap in the assessment of nuclear structures (NUREG/CR-7153), this work helps to further the understanding of creep behavior of massive concrete containment structures for decades to enable safe and long-term operation of these facilities. This project has developed a robust, experimentally validated model to predict creep in nuclear concrete structures for up 60 years using short-term creep data thereby enabling a longer service life of critical facilities and early detection of structural failure. The work presented in this report is a pairing of computational and experimental methods. For the first time, the time temperature superposition (TTS) principle was successfully used to generate a uniaxial creep compliance master curve to predict mortar creep response for up to 22,500 days (nearly 60 years) at a reference temperature of 20°C. These data were used as input into finite element analysis (FEA) codes that use highly realistic random, 3D concrete microstructures from reconstructed coarse limestone aggregates. Finite element analysis performed provides the ability to quickly upscale mortar viscoelastic behavior to long-term concrete creep/relaxation data. A master creep compliance curve, constructed from the TTS principle, spanning 27 years, was used to validate two and a half decades of simulated concrete creep. Concurrently, three different simulated wall specimens were designed to mimic the behavior of post-tensioned concrete nuclear containment facility vessel walls over time as a result of concrete creep. The specimens were designed with different thicknesses, transverse and longitudinal reinforcement ratios, and level of post-tensioning stress. Each specimen contained various instrumentation to measure internal concrete temperature, concrete strain, and post-tensioning strain hourly for over 3 years. The concrete creep model developed in this project, based on the FEA concrete simulations, was applied to simulate the structural-scale experiments of prestressed concrete walls conducted in this project using the Grizzly code. These models can represent the effects of reinforcing and prestressing. Although there are some discrepancies with the experimental data, the model can predict the overall trends of the creep response in these experiments. One of these experimental models was also applied to an extended time to demonstrate how the findings from this study can be used to predict the behavior of actual structures of interest that have been in service for extended periods of time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Grain2Mesh: Mesh Generation for Grain-Scale Nonlinear Elasticity Modeling

The nonlinear hysteretic behavior of rocks under cyclic loading is a crucial area of study in geomechanics. The macroscopic response of a variety of materials has been found to be contingent upon the behavior of the micro-scale structure. This project aims to develop a functional and maintainable software package for generating a multi-phase numerical mesh and accompanying simulation files for finite element modeling used in computational mechanics solvers. Meshes generated from images often lack key preprocessing that reduces noise and prevents mesh element distortion that can increase computational cost. By incorporating user feedback throughout, grain2mesh ensures a high-fidelity mesh that can be used to model grain-scale interactions such as shearing, crack propagation, and interfacial material contrast. Scientific applications of this software include material fracturing, stress-strain analysis for natural and engineered materials, and nonlinear meso-scale analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Role of crystallographic orientation on intragranular void growth in polycrystalline FCC materials

In this study, we study the effect of crystallographic orientation and applied triaxiality on the growth of intragranular voids. Two 3D full-field micromechanics methods are used, the dilatational visco-plastic fast-Fourier transform (DVP-FFT) and the crystal plasticity Finite Elements (CP-FE), both of which incorporate a combination of crystalline plasticity and dilatational plasticity. We demonstrate with several select cases that predictions of void growth from both formulations agree qualitatively. With the more computationally efficient DVP-FFT, additional effects of polycrystalline microstructure and the influence of nearest neighborhood are investigated. Crystals bearing a single intracrystalline void are studied in three types of 3D microstructural environments: isolated single crystals, individual equal-sized grains within a regular polycrystal, and individual variable sized grains within a polycrystal with grains and voids randomly located. We show that loading type plays a significant role. In strain-rate controlled conditions, voids in the hardest [111]-crystals grow the fastest in time, whereas in stress-controlled conditions, voids in the softest [100]-crystal grow the fastest in time. The analysis reveals that on average void growth is slower for the same starting orientation in the polycrystal than in the single crystal. We find that at the highest triaxiality tested that the correlation between crystal orientation and void growth rate in the polycrystal strengthens, drawing closer to that seen in the isolated single crystals. These results and model can help guide the microstructural design of polycrystalline materials with high strength and damage-tolerance in high-rate deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Comparative Study of Layer Heating and Continuous Heating Methods on Prediction Accuracy of Residual Stresses in Selective Laser Melted Tube Samples

Thermal distortion and residual stresses are two important factors that affect the quality and reliability of steel parts manufactured by laser powder bed fusion (LPBF) processes. A cost-effective model for evaluation of those heat effects is needed to refine the manufacturing process and provides insights into the product design and heat treatment. In this study, the layer heating method and sophisticated track-layer scanning method were applied to simulate the thermo-mechanical response of IN625 tube parts built by LPBF. Based on the similarity of temperature field in each layer deposit, a swept mesh was constructed to perform the thermal analysis for top layer, with the rest of layers referring to the temperature by node number offsetting. A novel explicit finite element analysis code accelerated by graphics processing unit was used for the massive-element numerical analysis. The computational accuracy and efficiency of the layer heating and track-layer scanning methods were compared in detail. It is shown that layer heating method can efficiently capture the pattern of stress distribution with reasonable accuracy in stress magnitude. The grouped track-layer scanning method can predict the residual stress and strain more accurately at a higher cost (5 ~ 10×). The elastic strain distribution was compared with the measurement by X-ray diffraction, confirming the accuracy of residual stress prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗
DOE OSTI2026

Computationally inexpensive part-scale thermal history of additive friction-stir deposition

This study presents an analytical model for steady-state power generation and tool heat loss in additive friction-stir deposition (AFSD), developed to enable part-scale thermal simulation while remaining computationally inexpensive. The model predicts total generated power, yielding 3.7–4.7 kW across deposition temperature setpoints of 400–460 °C for the deposition of AA6061 with a Be-Cu tool. This corresponds to 90–95% of the reported spindle power. Tool heat loss is experimentally determined by calibrating a steady-state energy balance between the generated power, the substrate-deposition thermal gradient, and a temperature dependent tool heat loss term: q tool (T) = a + b (T - 400°C) with a = 2.7 x 10 6 Wm -2 and b = 9.5 x 10 3 Wm -2 K -1 . The calibration indicates that about 69% of the generated heat is conducted into the tool for this configuration, which is much higher than previously reported. The calibrated heat-source is implemented in finite element software (Adamantine) to simulate the transient thermal history of a 100 cm 3 representative build in 8 min on a standard desktop (at 0.635 mm build-height resolution). For the first three layers, the substrate temperatures between simulation and experiment are within 10% mean absolute percentage error. Sensitivity analysis indicates that uncertainties in average deposition temperature and deformation localization (stir-zone geometry, depth, and spatial dependance of strain-rate and flow stress) dominate model variance, motivating additional experimental verification.

Additive Friction-Stir Deposition↗
DOE OSTI2024

Strength-ductility synergy through microstructural and compositional heterogeneity in directed energy deposition additive manufacturing of face-centered cubic materials

Directed energy deposition (DED) is an additive manufacturing (AM) process based on welding technology and offers the advantages of large build volume, high deposition rate, and ability to fabricate multi-material parts. Epitaxial continuous columnar grain growth is a characteristic microstructural feature of DED processed alloys. In this study, a bamboo-like microstructure (periodic alternation of equiaxed and columnar structure) was produced by adopting an intermittent deposition strategy in 316L stainless steel and Inconel 625. The formation of a bamboo-like alternating microstructure was confirmed through electron backscattered diffraction (EBSD) analysis. Hardness mapping showed that the columnar to equiaxed transition (CET) occurred at the region right below the fusion line. A finite element (FE) model was used to investigate the relationship between the temperature gradient (G) and the solidification rate (R). The FE model showed a low G/R ratio at the region right below the interface promoting the CET. The grain size and material-dependent deformation behaviors are analyzed using digital image correlation (DIC). The lower deformation on the fine-grain regions observed in DIC analysis is attributed to a higher strain hardening rate, which is confirmed through dislocation density analysis on a tensile-interrupted specimen. The periodically alternating grain size coupled with the microstructural changes caused by intermittent deposition strategy result in a better strength-ductility synergy in both single-material and bimetallic specimens.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Symmetry-driven persistent spin texture for the two-dimensional nonsymmorphic CdTe and ZnTe crystal structures

In this paper, two nonsymmorphic two-dimensional structures of CdTe and ZnTe are modeled, and using state-of-the-art density functional theory with the group theory of solids, their symmetry-enforced electronic properties are studied. The in-plane ferroelectricity coupled with strong spin-orbit coupling induces a unidirectional out-of-plane Rashba spin-orbit field that can host a momentum-independent uniform spin configuration known as persistent spin texture (PST) at the Brillouin zone center. PST in these structures is found to be robust against external perturbations such as strain, structural distortion, and independent of layer thickness. These unprecedented intrinsic spin transport properties hold utmost importance in spintronics, as the experimental stringent condition of equal Rashba and Dresselhaus constants [Phys. Rev. Lett. 90, 146801 (2003)] is eliminated. The calculated persistent spin helix wavelength of <~5 nm paves the way for developing next-generation nanosized nonballistic spin field-effect transistors compared with micrometer-sized GaAs/AlGaAs quantum wells. Further, these materials exhibit finite spin Hall conductivity at the band edges and hence can be used in ferromagnet-free spin Hall transistors. Although CdTe and ZnTe systems have been widely studied for photocatalysis and solar cell applications over the past few decades, their potential application in spintronic devices has not been explored. Mono/few layers of CdTe and ZnTe synthesized from (110) facets of bulk zinc-blende crystals [Nat. Commun. 3, 1057 (2012)] satisfy all symmetry operations of the nonsymmorphic space group and hence can be considered ideal materials to verify our theoretical results experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A crystal plasticity finite element method modeling of Zircaloy with hydride phases based on scanning electron micrographs

This work developed a microstructure-based finite element model to predict the stress state of alloys with second phase inclusions. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. Generation of digital microstructure was achieved through two steps of tessellations using software Neper and Matlab. The process is demonstrated using an example of Zircaloy material with secondary phases of hydrides dispersed within and stress-strain response of Zircaloy containing hydrides was predicted. The constitutive material model for Zircaloy in this study was based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as isotropic elasto-plastic material. A parametric study had been conducted to understand the effect of volume fraction, and lamellae thickness of the hydride phase on the mechanical properties of the overall material. Results can help designers to alter the manufacturing process to obtain the enhanced mechanical properties for components used in nuclear applications made by Zircaloy material.

Kulkarni, Shank S.↗
DOE OSTI2024

Digital image correlation and infrared thermography data for seven unique geometries of 304L stainless steel

Material Testing 2.0 (MT2.0) is a paradigm that advocates for the use of rich, full-field data, such as from digital image correlation and infrared thermography, for material identification. By employing heterogeneous, multi-axial data in conjunction with sophisticated inverse calibration techniques such as finite element model updating and the virtual fields method, MT2.0 aims to reduce the number of specimens needed for material identification and to increase confidence in the calibration results. To support continued development, improvement, and validation of such inverse methods—specifically for rate-dependent, temperature-dependent, and anisotropic metal plasticity models—we provide here a thorough experimental data set for 304L stainless steel sheet metal. The data set includes full-field displacement, strain, and temperature data for seven unique specimen geometries tested at different strain rates and in different material orientations. Commensurate extensometer strain data from tensile dog bones is provided as well for comparison. We believe this complete data set will be a valuable contribution to the experimental and computational mechanics communities, supporting continued advances in material identification methods.

36 MATERIALS SCIENCE↗
DOE OSTI2023

From Well Log to Formation Model: A Novel Laboratory Calibrated Methodology with Demonstration

This work demonstrates how the characterization and modeling of both elastic and creep properties are essential to describe zones in layered rock formations as either low-stress targets for stimulation or high-stress barriers to fracture growth. Prediction of fracture height is critical for designing stimulation operations in oil and gas wells. Ideally, fractures are placed in target zones which will produce hydrocarbons and should not propagate into zones expected to be unproductive or to produce unwanted fluids such as water which in turn must be treated and/or disposed. The essential task in designing stimulation plans is predicting which zones have low horizontal stresses and which will be high-stress barriers to fracture growth. Despite this importance, there are gaps in current knowledge and a complete workflow from laboratory characterization to a finite element model which includes time dependent rock deformation is required. While the research and methodology presented here also have application to CO 2or hydrogen storage, wastewater injection, and geothermal applications, the focus will be on hydrocarbon extraction. This thesis presents the results of a characterization-to-prediction workflow for the Caney shale, which is an emerging hydrocarbon resource in Oklahoma, USA. It begins with an investigation to enable critical evaluation of the Caney zonation into nominally “brittle” and “ductile” zones based on properties observed from well logs. It shows none of the zones are consistently “brittle” or “ductile” mechanical behavior based on the variety of definitions of these terms. However, the nominally ductile zones are weaker and more prone to creep. A laboratory investigation of samples including strength, elastic, and creep properties, is then used in a finite element model of stress evolution. The model includes both elastic deformation and viscoplastic creep. Results predict the least creep-prone layers to have the lowest horizontal stresses, therefore comprising hydraulic fracturing targets. The most creep-prone layers attain a horizontal stress similar to the vertical stress and therefore are predicted to be high stress barriers to hydraulic fracture stimulation. In addition to defining stimulation target intervals, the model shows how as tectonic strain rate increases, there is a transition from creep-dominated stresses to stresses dominated by elasticity.

Benge, Margaret↗
DOE OSTI2024

Agrivoltaic Racking Design Optimization Based on Wind and Snow Loading Finite Element Analysis

The racking structure of photovoltaic (PV) systems plays a critical role in ensuring the PV panels generate power properly as it provides integral structural support and sometimes even solar tracking ability. Distinct applications of PV systems require variations of racking structure designs, which yield different mechanical performances under external loading conditions from the environment, such as heavy wind and snow. Past works, such as Reddy et al. [1], have analyzed the pressure effects of wind and snow loading on the racking structure of conventional rooftop and utility PV systems, but there is little knowledge of how various agrivoltaic racking systems perform, in terms of stress and strain, under the same loading conditions. There is also a knowledge gap in the industry on a set of optimized agrivoltaic racking design standards. This study investigates the mechanical performance of various existing racking systems and proposes novel designs that are optimized for agrivoltaics applications under wind and snow loading.

Liao, Quanhuan↗
DOE OSTI2022

Generalized thermo-mechanical framework for heterogeneous materials through asymptotic homogenization

Abstract A fundamental understanding of the interaction between microstructure and underlying physical mechanisms is essential, especially for developing more accurate multi-physics models for heterogeneous materials. Effects of microstructure on the material response at the macroscale are modeled by using the generalized thermomechanics. In this study, strain gradient theory is employed as a higher-order theory on the macroscale with thermodynamics modeled as a first-order theory on the microscale. Hence, energy depends only on the temperature such that we circumvent an extension of Fourier’s law and analyze the “simplest” thermo-mechanical model in strain gradient elasticity. Developing multiphysics models for heterogeneous materials is indeed a challenge and even this “simplest” model in generalized thermomechanics creates dozens of parameters to be determined. We develop a thermo-mechanical framework, in which microstructure is modeled as a periodic structure and through asymptotic homogenization approach, higher-order parameters at macroscopic scale are calculated. To illustrate the importance of higher-order parameters in overall thermo-mechanical response of a heterogeneous materials, finite element method (FEM) is employed with the aid of open-source codes (FEniCS). Verification example of a bulk system and several case studies of porous structures demonstrate how such numerical framework can be beneficial in the design of materials with tailored microstructures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

In-situ effect in cross-ply laminates under various loading conditions analyzed with hybrid macro/micro-scale computational models

In this article, multi-scale finite element analyses based on three-dimensional (3D) hybrid macro/micro-scale computational models subjected to various loading conditions are carried out to examine the in-situ effect imposed by the neighboring plies on the failure initiation and propagation of cross-ply laminates. A detailed comparative study on crack suppression mechanisms due to the effect of embedded laminar thickness and adjacent ply orientation is presented. Furthermore, we compare the results of in-situ transverse failure strain and strength between the computational models and analytical predictions. Good agreements are generally observed, indicating the constructed computational models are highly accurate to quantify the in-situ effect. Subsequently, empirical formulas for calculating the in-situ strengths as a function of embedded ply thickness and different ply angle between embedded and adjacent plies are developed, during which several material parameters are obtained using a reverse fitting method. Finally, a new set of failure criteria for σ 22 -τ 12 , σ 22 -τ 23 , and σ 11 -τ 12 accounting for the in-situ strengths are proposed to predict laminated composites failure under multi-axial stress states. This study demonstrates an effective and efficient computational technique towards the accurate prediction of the failure behaviors and strengths of cross-ply laminates by including the in-situ effects.

36 MATERIALS SCIENCE↗