{'fraction_non_alphanumeric': 0.033083034062465995, 'fraction_numerical': 0.012950266623136359, 'mean_word_length': 5.298377884395705, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Integrated Physical ‐ Constitutive Computational Framework for Plastic Deformation Modeling

: An integrated framework for deformation modeling has been developed, which combines a physical state parameter ‐ based formulation for microstructure evolution during plastic deformation processes with constitutive creep models of polycrystalline materials. The implementations of power law, Coble, Nabarro–Herring and Harper–Dorn creep and grain boundary sliding are described and their contributions to the entire stress response at a virtual applied strain rate are discussed. The present framework simultaneously allows calculating the plastic deformation under prescribed strain rate or constant stress, as well as stress relaxation after preceding stress or strain loading. The framework is successfully applied for the construction of deformation mechanism maps.


Introduction
The physical understanding of creep mechanisms has been of significant academic interest within the last decades. Especially, the creep response of pure Al has been investigated over a wide temperature and stress range. The physical process of plastic flow is still the subject of controversial discussions for some deformation regimes [1][2][3][4], but semi-empirical equations in literature satisfactorily reproduce the experimentally obtained stress-strain relations. The power law regime is associated with vacancy-assisted climb of dislocations, while diffusional flow of atoms is attributed to Coble creep [5] and Nabarro-Herring creep [6]. At low stresses and large grain sizes, Harper-Dorn creep dominates [7][8][9][10][11]. In fine grain structures, grain boundary sliding can control the deformation process [12][13][14][15][16].
In the present work, characteristic creep relations are brought into a homogeneous form and used to calculate the stress response of polycrystalline materials to a constant applied strain rate or the strain rate resulting from an applied external stress load, respectively. Since this framework is based on fundamental plastic deformation mechanisms, thermomechanical simulations are possible for a broad spectrum of metallic materials. Therefore, only characteristic material parameters, which are described in the following discussion, need to be calibrated for each material. Within this work, the model is exemplarily applied to Al. An inverse summation rule of the individual stress contributions ensures that the dominant deformation mechanism mainly defines the overall thermally activated stress response. Microstructural changes during the deformation process are represented by the average dislocation density evolution. The utilized Kocks-Mecking approach considers dislocation generation as well as dynamic and static recovery for dislocation annihilation. The combination of thermally activated stress contributions and the athermal stress contribution from dislocation hardening allows simulating the stress-strain relations over a wide range of temperatures and strain rates. This framework allows for modeling of plastic deformation by dislocation glide as well as of primary creep and secondary creep. The strain rate acceleration and damage accumulation due to cavitation, cracks or both is not considered, which is typical for tertiary creep. Consequently, this framework is valid for experiments, where creep damage is effectively suppressed, or for materials which show a pronounced secondary creep stage. As an additional feature, stress relaxation tests from literature are successfully simulated within this work to compare the simulation framework with experimental data.
Although not elaborated in the present work, the framework is fully capable of considering microstructural features, such as fine grain hardening, dislocation locking by Cottrell atmospheres or precipitates. Application of the latter two mechanisms has been demonstrated recently by Soliman et al. [17]. These features provide further barriers for ongoing dislocation glide, leading to a higher stress response to an applied strain rate. Diffusion creep, such as Coble creep and Nabarro-Herring creep, are not affected by those microstructural barriers.

The Model
The following two subchapters describe the implementation of both aspects of this framework, which are the constitutive creep relations on the one hand, and physically-based microstructure evolution on the other hand. A combination of the prevailing stress contribution is presented in section 3, Results and Discussion.

Constitutive Modeling of Creep
The deformation characteristics of crystalline materials are mainly controlled by the applied strain rate / external stress loading as well as temperature. Dependent on the particular combination of these parameters, the deformation process is made possible by different physical mechanisms. Generally, at high temperatures, the deformation characteristics show considerable dependency on the applied stress in combination with a low amount of strain hardening. At high temperatures, various creep mechanisms can be operative simultaneously, whereas mostly dislocation glide facilitates plastic deformation at low temperature, high stress or both. Interactions with obstacles, such as grain boundaries, solute atoms, precipitates or other dislocations cause a barrier for an ongoing movement of the mobile dislocations that carry the plastic deformation process. Consequently, the various obstacle strengths determine the rate of dislocation motion and, thus, the strain rate, , which is typically represented by an Arrhenius equation [18][19][20][21][22][23] as is a constant [24], the Boltzmann constant, T is the temperature and ∆ is the shear stressdependent free activation enthalpy, which can be expressed by [20] ∆ ∆ · 1 . (2) ∆ is the is the total free energy to overcome the obstacle barrier, the applied shear stress, ̂ the mechanical threshold, which represents the activation energy barrier in the absence of thermal activation, and p and q define the shape of the energy barrier. In the following, a box shaped barrier is assumed by p = q = 1. A combination of Equations (1) and (2) leads to At lower stresses, higher temperatures or both, (T 0.3· ), diffusion-controlled deformation regimes, instead of dislocation glide, become dominant.
is the melting temperature. Although the following phenomenological equation does not describe the creep process on a physical basis, it is often successfully applied in literature [20,25,26] and it is also included within this framework as is the tracer diffusion coefficient of the matrix-forming element, the shear modulus and b the Burgers vector. Brown and Ashby [26] empirically correlated the coefficients, A and n, based on detailed evaluations of creep data of bcc, fcc, hpc metals (including their alloys), dispersion hardened materials, cubic elements, alkali halides and oxides. Missing explanations for the observed values of n and a huge variety in A of many orders of magnitude motivated Frost and Ashby to include dislocation core diffusion additionally to the bulk matrix diffusion coefficient [20]. Consequently, two rate equations are obtained: At low temperatures or high stresses, dislocation core diffusion is dominant, while lattice diffusion is the dominant mechanism at high temperatures and low stresses. In the following, both mechanisms, which are termed as low-temperature creep, LT creep, and hightemperature creep, HT creep, are dealt with separately. The weighted combination of creep mechanisms is discussed in section 3. The strain rate, , for HT creep can be expressed as is a coefficient and is microstructure-dependent within the present framework (see section 3).
is the effective diffusion coefficient that includes trapping of vacancies at solute atoms [27], excess vacancies [28] and dislocation pipe diffusion [29,30] as inherent diffusion mechanisms. The migration of vacancies predominantly determines diffusion in polycrystalline solids. Here, the creation and annihilation of vacancies at grain boundaries, dislocation jogs or Frank loops are considered within the FSAK model [28]. Additionally, an increasing dislocation density during deformation causes a higher effective diffusion coefficient based on a higher pipe diffusion contribution. The strain rate for the LT creep contribution is calculated as is a constant, is a state-dependent parameter and is the dislocation core diffusion coefficient, which, for Al-based alloys, is given by a correlation with the tracer diffusion coefficient [31]: and are the activation energies for bulk diffusion and dislocation core diffusion. R is the universal gas constant. Above a certain stress of approximately 10 -3 , the stress-strain relation typically changes, which is referred to as power law breakdown (PLB). This phenomenon is intensively discussed in literature [32,33]. Within the present framework, the following empirical equation is utilized , A , D and refer to the prevailing deformation mechanisms and α′ specifies the stress level at which the power law break down starts.
At low applied stresses ( / 10 −5 ) and small grain sizes, diffusional flow determines the deformation process. In contrast to power law creep, strain is caused by the diffusion of vacancies, rather than dislocation mechanisms. This deformation mode has first been described by Nabarro and Herring [6]. The corresponding creep rate , as implemented within the present framework, is given by [20] where d is the grain size. Since lattice diffusion is the predominate mechanism, the effective matrix diffusion coefficient, , is utilized. Vacancy diffusion along grain boundaries instead of lattice diffusion leads to Coble creep [5]. The strain rate, , for Coble creep is formulated as is the effective boundary thickness and is the grain boundary diffusion coefficient, which, for Al-based alloys is given [31] as is the activation energy for grain boundary diffusion. At low stresses and large grain sizes, another deformation mechanism occurs, which was first reported by Harper et al. in 1957 [7]. In the following fifty years, the Harper-Dorn regime has been identified in many studies [7][8][9][10][11]. A phenomenological relation between the applied stress, , and the creep rate, , reads where is a coefficient. Equal to the Nabarro-Herring creep, increases linearly with the applied stress, , (Newtonian nature) and the activation energy is of the magnitude of self-diffusion. However, measured steady-state creep rates of Harper-Dorn creep are often higher by a factor of approximately 1400. Low initial dislocation densities, as well as high purity of the material, favor Harper-Dorn creep, as reported by [4,11]. The low steady-state dislocation density typical for this deformation mechanism remains constant and is independent of [8]. A possible explanation is given by the Ardell-Przystupa-Lee (APL) dislocation network theory [1]. The mechanism of the deformation process during Harper-Dorn creep has been discussed intensively in literature [4,[7][8][9][10] and its analysis is beyond the of scope of this work. Another deformation mechanism, which can control the plastic behavior of fine-grained polycrystalline materials, is based on grain boundary sliding (GBS). Many observations confirm an inverse dependency of the strain rate, , on squared or cubed grain size [12][13][14][15][16]. If the activation energy is associated with grain boundary diffusion, , the creep rate is [12]: where the coefficient 2 10 . If effective lattice diffusion is the rate controlling mechanism, the following equation is utilized [12]: , (14) where the coefficient 8 10 . The impact of grain boundary sliding processes on conventional Ashby-Frost [20] and Langdon-Mohamed [34,35] deformation mechanism maps are described by Lüthy et al. [12]. It is shown for Ni that -and -controlled GBS strongly reduces Coble creep and suppresses Nabarro-Herring creep to insignificant normalized stresses within the Ashby maps. Although Coble creep and Nabarro-Herring creep are overtaken by grain boundary sliding in Al at all temperatures at an applied normalized stress of / 10 , these mechanisms are incorporated into the framework to achieve a generalized form of all creep mechanisms. Except for the dislocation glide mechanism, all thermally activated deformation processes introduced so far are diffusion-controlled, either by lattice diffusion, pipe diffusion or grain boundary diffusion. The present models describe the activation of dislocation glide, dislocation climb and vacancy diffusion; however, under the assumption of a constant microstructure and no accumulation of creep damage, which is the condition of steady state creep. In real materials, the dislocation density changes during deformation before the steady-state dislocation density is reached. This microstructure evolution is characteristic for the transient range of creep and is discussed in the following.

Physical Modeling of Microstructure Evolution
For the simulation of the stress evolution during thermomechanical treatments, the knowledge of the current microstructure is mandatory. Therefore, physically based models are utilized, which can simulate the evolution of grain sizes [36,37], nucleation and growth of precipitates, or the current dislocation density. All these models are already implemented within the thermo-kinetic software package MatCalc (http://matcalc.at). To incorporate the dislocation strengthening contribution accompanying the transient deformation region, following extended one-parameter model of Kocks and Mecking is used to calculate the temperature and strain rate dependent dislocation density evolution [24]: The Taylor factor, M, relates the macroscopically observed flow stress to the critical resolved shear stress acting on a slip plane within a polycrystalline material. A and B are material-specific coefficients and is the critical annihilation distance between two dislocations [24]. C is a constant and is the equilibrium dislocation density. The generation of dislocations is caused by dislocation storage in the crystal, being inversely proportional to the current mean distance between dislocations and the dislocation density, , respectively. Annihilation of dislocations occurs by cross slip processes at low and intermediate temperatures and by vacancy-assisted climb at high temperatures. The latter annihilation process represents static recovery, which is marked by the index s. The Taylor equation correlates the athermal, plastic stress contribution with the dislocation density [38,39] as The MatCalc software package includes models for the growth and shrinkage of subgrains during plastic deformation, which are out of scope in the present work. The calculation of the average dislocation density for the athermal stress contribution has often been validated [24] and possible industrial applications, such as hot rolling processes, which additionally include recrystallization, have been demonstrated recently [40].
In terms of the possible evolution of the precipitate microstructure, the SFFK model [41][42][43] is utilized as implemented in the MatCalc software package. Within this framework, a thermodynamic system out of equilibrium state is described by the Gibbs energy, G, given by The index i represents the component i in the matrix and the index k represents a specific precipitate with the radius .
is the number of moles, the chemical potential, c the concentration and accounts for the elastic energy, which is associated with the elastic stress field around the precipitate. is the interface free energy density.
The steady-state nucleation rate, , which is defined as the number of formed nuclei per unit volume and unit time, is given by [44,45] where is the Zeldovich factor, * the atomic attachment rate, the number of available nucleation sites and ∆ * is the critical nucleation energy. Equations (17) and (18) provide the basis for rather complex precipitation evolution models, which cause strengthening contributions by shearing-and non-shearing mechanisms. A detailed discussion of precipitation strengthening is beyond the scope of this work and the authors refer to Ahmadi et al. [46]. Solid solution strengthening models as well as the cross core diffusion effect, which can lead to negative strain-rate sensitivity, are implemented in MatCalc and successfully applied, e.g., in Kreyca and Kozeschnik [24]. These thermally activated stress contributions are linked to the present framework by the mechanical threshold ̂ (Equation (3)). In the following, ̂ is assumed to remain constant for the benefit of clear discussion.

Results and Discussion
In classical creep tests, a constant stress is applied and the resulting creep rate is measured. In an inverse methodology, the creep relations described in the previous section can be utilized to calculate stresses at an applied constant true strain rate. The stress response of pure dislocation glide is calculated by Equation (3) and leads to Since Equation (19) delivers negative stress values above a certain critical temperature, , the following alternative equation is widely used in literature [24]: The following equations summarize the present implementation of the LT creep, the HT creep and the Harper-Dorn creep in the integrated framework of deformation modeling. The indices display the different associated pre-factors, A, the diffusion coefficients, D, as well as the exponents, n, with The power law breakdown regime is implemented with The stress response of Nabarro-Herring creep, , and Coble creep, , are calculated by Equations (9) and (10) and -and -controlled GBS are calculated by Equations (13) and (14). Kreyca [47] has correlated the thermally activated dynamic yield stress contribution at high temperatures with the HT power law creep relation, Equation (5), leading to the relation where ∆F is the activation energy for lattice diffusion and is the effective diffusion coefficient. Although it is not very common to connect the mechanical threshold concept ̂ with the creep behavior, this approach takes the influence of microstructure features, such as precipitates on the creep rate, into account. It should be emphasized that the particular microstructural state also influences the diffusion coefficient, D , of Equation (23) by trapping effects, pipe diffusion or by excess vacancies. For correlating of Equation (6) with the mechanical threshold concept, the core diffusion z as well as the activation energy for pipe diffusion are utilized. This concept is used for calibrating and in the following. Figure 1 shows the Harper-Dorn creep (HD) and power law creep (PL) regime, including the power law breakdown (PLB) at higher stresses. The symbols indicate the normalized strain rate, · · / · · , as a function of the normalized stress / from Straub and Blum [48] for pure 99.999 Al.  [48]. The normalized strain rate is defined by · · / · · .
All simulation input parameters are listed in Table 1, but instead of αʹ = 1000, which is a typical value in literature [20], better agreement with experimental results are achieved with αʹ = 1300 in this case. In determination of the overall creep behavior, the assumption is made that the creep mechanism with the lowest resulting stress and, therefore, lowest deformation resistance, is dominant. The following summation rule is utilized = , (24) where is the coupling coefficient. The total stress contribution at a certain temperature and strain rate is the sum of the athermal stress contribution, , of Equation (16) and of Equation (24). Except for the dislocation glide mechanism, each creep regime is characterized by a linear relation in a double-logarithmic plot until the power law breakdown is reached at high stresses. Exemplary stress-strain rate relations are illustrated in Figure 2 for all discussed deformation regimes at different temperatures from 100 to 500 °C, whereas of Equation (24) is represented by the bold line, labeled "DYN". Table 1 lists all simulation input parameters. Activation energy for grain boundary J/mol 60,200 11 [31] The exponent n of the power law mechanism (Equations (5) and (6)) determines the slope of the LT creep and the HT creep range, until the PLB becomes dominant. The Nabarro-Herring creep, the Coble creep and the Harper-Dorn creep deliver parallel shifted results based on the linear relation of stress and strain rate. The slope of the GBS regime is steeper and crosses the Coble creep mechanism at a specific stress value, which depends on the initial grain size, which is 5 μm in this case, and the deformation temperature. At 100 °C, 200 °C and 300 °C, grain boundary diffusion-controlled GBS is dominant at intermediate stress values, which is finally overtaken by dislocation glide (at 100 °C) or high temperature creep (at 200 °C and 300 °C) at high stresses. Above 400 °C, -controlled GBS replaces -controlled GBS as dominant deformation mode. Figure 2 illustrates that GBS controls plastic flow in the intermediate stress range, separating diffusional flow and power law creep range, in good agreement with literature [12]. The transition of these deformation regimes is governed by the exponent in Equation (24). The higher the coupling coefficient , the closer is the resulting stress, , to the prevailing mechanism. A value of = 10 is used in the present analysis, which assures that the resulting curves are rather close to the strain rates predicted by the dominating mechanism and numerical artifacts are minimized. The mechanical threshold, , which represents the activation energy barrier in the absence of thermal activation, needs to be defined to calculate the dislocation glide contribution, especially at low temperatures (see Equation (20)). Since ̂ comprises solid solution hardening, precipitation hardening and cross core diffusion hardening (see section 2.2, Physical modeling of microstructure evolution), the evolution of the microstructure during the deformation processes is inherently included in the integrated model. In the present analysis, no precipitation processes occur in this setup, leading to a constant mechanical threshold ̂ = 145 MPa and consequently, to constants and . The present framework allows calculating both, stresses during constant strain rates as well as time-dependent stress relaxation. The stress reduction, , during relaxation is assumed to occur along the elastic region with where E is the Youngs modulus. Hence, the stress reduction, , is defined by the current dynamic strain rate, . The new stress value, ∆ , in each iteration step, Δt, of the numerical integration procedure is given by For a given materials' microstructure, each stress value corresponds to a unique strain rate, which determines the next stress relaxation step in Equation (25). In this setup, simulations of relaxation tests are carried out. Figure 3a shows relaxation curves of an AlMg4.5Mn alloy, after compressing cylindrical specimens with a constant strain rate of 0.01 (1/s) at different temperatures (350 °C, 400 °C, 450 °C, 500 °C). The experimental points are adopted from the work of Falkinger and Simon [50]. Figure 3b depicts the decreasing strain rate during the relaxation process, starting from 0.01 (1/s). The higher the deformation temperature in the pre-step, the lower the initial stress level for relaxation is, resulting in lower strain rates. The stress reduction of the HT creep (HT), the LT creep (LT), the Coble creep, lattice diffusion controlled GBS as well as grain boundary diffusion-controlled GBS are illustrated for 350 °C in Figure 3c and 500 °C in Figure 3d. According to Equation (24), Nabarro-Herring and Harper-Dorn creep are negligible in this example setup, due to their exceeding stress contributions. The calibration settings are summarized in Table 1 Figure 3 is to illustrate the possibility of calculating stress relaxations within this framework for a specific set of parameters. When the stress relaxation starts after compressing the cylinder, HT creep controls the plastic flow at 350 °C and 500 °C, as shown in Figure 3c and Figure 3d at 0 s. As the strain rate rapidly decreases in both cases, lattice diffusioncontrolled GBS becomes the dominant mechanism. At very low stresses and strain rates, Coble creep becomes the dominant mechanism at 500 °C. The total thermal stress contribution, , of Equation (24) is named "DYN" in Figure 3c and Figure 3d. The athermal stress contributions at the high temperatures and small strain rates (Equation (16)) are not significant due to a negligible amount of dislocation generation. An additional possible application of this framework is the construction of deformation maps, which provide information about the dominant deformation mechanism at specific temperatures and applied shear stresses. Figure 4 shows the calculated normalized stresses ( / ) for different strain rates and temperatures in dependence of the grain diameter ((a) 5 μm, (b) 50 μm) for pure Al. Since grain rotation is not considered in the present work, application of the framework to sub-micrometer grain sizes might deliver inaccurate results (ultra-fine grained materials). The numbers in Figure 4 denote the corresponding deformation mechanisms and the black, bold lines indicate the transitions between the different deformation regimes.   Figure 4 shows representative deformation maps that have been generated without explicit calibration and thus they allow only for qualitative and exemplary analysis. In this sense, Figure 4a illustrates that GBS becomes the deformation-controlling mechanism at small grain sizes for a wide temperature range, which is in good agreement with literature [12]. Increasing GBS-fields are accompanied by shrinking LT creep and HT creep fields within the deformation maps. The dislocation density rapidly increases at high strain rates, leading to an increasing pipe diffusion coefficient, , and a larger LT creep field (number 2). A detailed description of deformation maps and a catalogue of maps for various groups of materials are provided by Frost and Ashby [20]. The effects of GBS on Frost-Ashby deformation maps are discussed in [12].

Conclusions
The implementation of creep mechanisms and their contributions to the stress response to an applied strain rate within an integrated framework for deformation modeling is discussed. Although, steady state is usually presumed for the secondary creep relations under no creep damage condition, athermal dislocation strengthening contribution and the mechanical threshold concept are considered. In the present framework, an integration of the influence of precipitates on the creep regimes is possible but was beyond the scope of this work. It is demonstrated that plastic flow is either governed by dislocation glide, diffusion of either vacancies or dislocations or by grain boundary sliding. The effective diffusion coefficient, which is influenced by pipe diffusion enhancement, excess vacancies or trapping effects, is accounted for. This new framework allows calculating the creep or plastic deformation behavior, while applying either a constant stress, stress relaxation after previous deformation or constant strain rate loading. Exemplary deformation maps for Al are calculated with the present model. Considering the microstructure evolution within a wide range of temperature and strain rate makes this framework suitable for many industrial applications, such as hot forming processes or strain relaxation during stress relieve annealing.