This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. The ...model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.
Display omitted
This paper presents a microstructure sensitive model for predicting mechanical response and texture evolution of metals in the dynamic recrystallization regime. A recently proposed viscoplastic ...self-consistent (VPSC) formulation for the prediction of recrystallization driven by strain energy and intragranular misorientation is extended to hexagonal close-packed (hcp) metals. The model is applied to the dynamic recrystallization of magnesium alloy WE43 at different temperatures and strain rates. Model predictions in terms of stress-strain response and texture evolution are compared to the experimental measurements and acceptable agreement is achieved. According to the model predictions, superplastic behavior of nuclei was found to be the dominant softening mechanism at high temperatures and low strain rates. High concentration of precipitates at the grain boundaries and presence of alloying elements are the likely causes of low boundary mobility, resulting in nucleation dominated dynamic recrystallization. Relatively strong basal compression textures indicate dominant activity of basal slip, which can be achieved only through large difference in slip resistance between soft basal and hard prismatic and pyramidal modes.
•Dynamic recrystallization nuclei remain small due to low boundary mobility.•Superplastic behavior of recrystallization nuclei causes strong softening.•Pyramidal and prismatic slips are considerably harder than basal slip.•Sharp increase of dislocation removal observed at elevated temperatures.
Gels are comprised of polymer networks swelled by some interstitial solvent. They are under wide investigation by material scientists and engineers for their broad applicability in fields ranging ...from adhesives to tissue engineering. Gels’ mechanical properties greatly influence their efficacy in such applications and are largely dictated by their underlying microstructures and constituent-scale properties. Yet predictively mapping the local-to-global property functions of gels remains difficult due - in part - to the complexity introduced by solute-solvent interactions. We here introduce a novel, discrete mesoscale modeling method that preserves local solute concentration-dependent gradients in osmotic pressure through the Flory-Huggins mixing parameter, χ. The iteration of the model used here replicates gels fabricated from telechelically crosslinked star-shaped polymers and intakes χ, macromer molecular weight (Mw), crosslink functionality (f), and as-prepared solute concentration (ϕ*) as its inputs, all of which are analogues to the control parameters of experimentalists. Here we demonstrate how this method captures solvent-dependent homogenization (χ≤0.5) or phase separation (χ>0.5) of polymer suspensions in the absence of phenomenological pairwise potentials. We then demonstrate its accurate, ab initio prediction of gel topology, isotropic swelling mechanics, and uniaxial tensile stress for a 10k tetra-PEG gel. Finally, we use the model to predict trends in the mechanical response and failure of multi-functional PEG-based gels over a range of Mw and f, while investigating said trends’ micromechanical origins. The model predicts that increased crosslink functionality results in higher initial chain stretch (as measured at the equilibrated swollen state) for gels of the same underlying chain length, which improves modulus and failure stress but decreases failure strain and toughness.
Both the mechanical properties of oxide nuclear fuel and the high burn-up grain restructuring are known to be related to dislocation motion. Previous studies show that at temperatures above 1000 K ...slip on {001} crystallographic planes in UO2 is dominant, and 1/2〈110〉{001} dislocations should have the highest mobility. However, no information is available on the temperature and stress dependence of dislocation velocities in UO2. This information is required as input in higher scale simulation techniques, such as dislocation dynamics, to gain an insight into the unusual plastic behaviour of this oxide. In this study we employed molecular dynamics to gather a representative set of data on the 1/2〈110〉{001} dislocation motion in UO2 at temperatures from 1000 to 2000 K and shear stress from 25 to 1000 MPa. We found that dislocation motion was thermally activated, with two distinct modes of motion. Based on this knowledge, we selected an analytic expression for dislocation velocity and performed fitting to MD data. This expression described the simulation data well, and was also used to make a physically justified prediction for dislocation velocities at typical operating conditions of the oxide nuclear fuel.
Display omitted
•Motion of the ½〈110〉{001} dislocations in UO2 is investigated by molecular dynamics.•This motion is thermally activated at T ≤ 2000 K and σ ≤ 1000 MPa.•At T < 1750 K dislocations move by propagation of double kinks.•At higher temperatures the dislocation line is a vibrating string.•Analytic function for dislocation glide mobility in the kink chain mode is proposed.
•This study investigates the cyclic plasticity behavior of AM 316 L stainless steel.•Hill's anisotropic yield criterion is utilized in a 3D return mapping algorithm.•Kinematic and isotropic hardening ...is integrated together with a damage model.•The model accurately captures full cyclic hysteresis loops with historical effects.•An energy-based fatigue life model is derived from the algorithm and fatigue tests.
This study addresses the critical need for a constitutive model to analyze the cyclic plasticity of additively manufactured 316L stainless steel. The anisotropic behavior at both room temperature and 300 °C is investigated experimentally based on cyclic hysteresis loops performed in different orientations with respect to the build direction. A comprehensive constitutive model is proposed, that integrates the Armstrong-Frederick nonlinear kinematic hardening, Voce nonlinear isotropic hardening and Hill's anisotropic yield criterion within a 3D return mapping algorithm. The model was calibrated to specimens in the 0° and 90° orientations and validated with specimens in the 45° orientation. A single set of hardening parameters successfully represented the elastoplastic response for all orientations at room temperature. The algorithm effectively captured the full cyclic hysteresis loops, including historical effects observed in experimental tests. A consistent trend of reduced hardening was observed at elevated temperature, while the 45° specimen orientation consistently exhibited the highest degree of strain hardening. The applicability of the model was demonstrated by computing energy dissipation for stabilized hysteresis loops, which was combined with fatigue tests to propose an energy-based fatigue life prediction model.
This paper presents a recrystallization model driven by intragranular orientation gradients and strain energy fields calculated by means of the viscoplastic self-consistent (VPSC) formulation. The ...VPSC model is extended for calculation of the coupling between intragranular stress fluctuations with corresponding second moments of lattice spin and misorientation fields in the grains. Access to these quantities allows modelling of transition bands and nucleation kinetics. In the proposed recrystallization model, grain growth is assumed to be proportional to the difference between the stored energy of each grain and that of the effective medium. Recrystallization textures for several cubic metals are simulated, showing good agreement with corresponding experiments. The model reveals the importance of considering appropriate, microstructurally-based and orientation-dependent recrystallization nucleation mechanisms. The recrystallization texture of heavily rolled copper with a strong cube texture component is found to be a consequence of nucleation at transition bands, which is also the cause of the recrystallization textures in compressed iron and drawn copper wire. In contrast, the recrystallization texture of rolled interstitial-free steel is found to be caused by grain boundary nucleation occurring in grains with the highest strain energy.
Display omitted
Exascale applications: skin in the game Alexander, Francis; Almgren, Ann; Bell, John ...
Philosophical transactions of the Royal Society of London. Series A: Mathematical, physical, and engineering sciences,
03/2020, Letnik:
378, Številka:
2166
Journal Article
Recenzirano
Odprti dostop
As noted in Wikipedia,
refers to having 'incurred risk by being involved in achieving a goal', where '
is a synecdoche for the person involved, and
is the metaphor for actions on the field of play ...under discussion'. For exascale applications under development in the US Department of Energy Exascale Computing Project, nothing could be more apt, with the
being exascale applications and the
being delivering comprehensive science-based computational applications that effectively exploit exascale high-performance computing technologies to provide breakthrough modelling and simulation and data science solutions. These solutions will yield high-confidence insights and answers to the most critical problems and challenges for the USA in scientific discovery, national security, energy assurance, economic competitiveness and advanced healthcare. This article is part of a discussion meeting issue 'Numerical algorithms for high-performance computational science'.
Using neural networks to handle intractability problems and solve complex computation equations is becoming common practices in academia and industry. It has been shown that, although complicated, ...these problems can be formulated as a set of equations and the key is to find the zeros of them. Zeroing neural networks (ZNN), as a class of neural networks particularly dedicated to find zeros of equations, have played an indispensable role in the online solution of time-varying problem in the past years and many fruitful research outcomes have been reported in the literatures. The aim of this paper is to provide a comprehensive survey of the research on ZNNs, including continuous-time and discrete-time ZNN models for various problems solving as well as their applications in motion planning and control of redundant manipulators, tracking control of chaotic systems, or even populations control in mathematical biosciences. By considering the fact that real-time performance is highly demanded for time-varying problems in practice, stability and convergence analyses of different continuous-time ZNN models are reviewed in detail in a unified way. For the case of discrete-time problems solving, the procedures on how to discretize a continuous-time ZNN model and the techniques on how to obtain an accuracy solution are summarized. Concluding remarks and future directions of ZNN are pointed out and discussed.
Strain-gradient (SG) plasticity refers to a class of non-local theories in which gradients of plastic slip determine the storage of geometrically necessary dislocations, introducing a length-scale ...dependence in the mechanical behavior of crystalline materials, which is otherwise lacking in local theories. In this work, we incorporate lower-order (LO) and higher-order energetic (HOE) strain-gradient effects into a crystal plasticity fast Fourier transform (FFT)-based formulation to investigate the interplay of the length scale that each strain-gradient term introduces at the microscale, and the mechanical properties that result at the macroscale. For an applicable range of length scales, we consider two systems: a 1-D two-phase face centered cubic (FCC) laminate and a 3-D FCC polycrystal, and two uniaxial deformation modes: monotonic tension and cyclic tension–compression. We show that increases in the individual LO and HOE length scales increase the hardening rate and strength of the material, respectively. When combined, the strong LO hardening is less pronounced than the effect alone due to the lowering of the gradients due to the HOE microstress. We demonstrate that the LO and HOE hardening manifest as “isotropic” (yield surface expansion) and “kinematic” (yield surface shift) effects, respectively, consistent with their theoretical origins. We show that in cyclic loading, the Bauschinger effect emerges in both local and non-local calculations and link its origins and severity to the behavior in the strain field, slip-system rates, and the HOE microforce.
•The interplay between strain-gradient length-scale parameters is studied.•Increasing the higher-order energetic parameter increases polycrystal yield strength.•Increasing the lower-order (LO) length scale increases the strain-hardening rate.•The higher-order effect reduces gradients, limiting the influence of the LO effect.•The higher-order energetic effect enhances the local Bauschinger effect.