CC BY
29
УДК 537.37
Crystal plasticity study on stress and strain partitioning in a measured 3D dual phase steel microstructure
M. Diehl, D. An, P. Shanthraj, S. Zaefferer, F. Roters, and D. Raabe
Dual phase steels are advanced high strength alloys typically used for structural parts and reinforcements in car bodies. Their good combination of strength and ductility and their lean composition render them an economically competitive option for realizing multiple lightweight design options in automotive engineering. The mechanical response of dual phase steels is the result of the strain and stress partitioning among the ferritic and martensitic phases and the individual crystallographic grains and subgrains of these phases. Therefore, understanding how these microstructural features influence the global and local mechanical properties is of utmost importance for the design of improved dual phase steel grades. While multiple corresponding simulation studies have been dedicated to the investigation of dual phase steel micromechanics, numerical tools and experiment techniques for characterizing and simulating real 3D microstructures of such complex materials have been emerged only recently. Here we present a crystal plasticity simulation study based on a 3D dual phase microstructure which is obtained by EBSD tomography, also referred to as 3D EBSD (EBSD — electron backscatter diffraction). In the present case we utilized a 3D EBSD serial approach based on mechanical polishing. Moreover, sections of the 3D microstructure are used as 2D models to study the effect of this simplification on the stress and strain distribution. The simulations are conducted using a phenomenological crystal plasticity model and a spectral method approach implemented in the Düsseldorf Advanced Material Simulation Kit (DAMASK).
Keywords: crystal plasticity, dual phase steel, microstructure, DAMASK, spectral method, 3D EBSD
Исследование распределения напряжений и деформаций в трехмерной микроструктуре двухфазной стали в рамках теории пластичности кристаллов
M. Diehl, D. An, P. Shanthraj, S. Zaefferer, F. Roters, D. Raabe
Двухфазные стали — это современные высокопрочные сплавы, применяемые для изготовления конструкций и силовых элементов автомобильных кузовов. Сочетание прочности и пластичности, а также обедненный легирующими элементами состав делают их экономически выгодным материалом для производства облегченных конструкций, в частности в автомобилестроении. Механический отклик двухфазных сталей обусловлен распределением деформаций и напряжений между ферритной и мартенситной фазами, а также отдельными зернами и субзернами этих фаз. Таким образом, изучение влияния микроструктурных особенностей на глобальные и локальные механические свойства имеет первостепенное значение для разработки новых марок двухфазных сталей. Микромеханика двухфазных сталей довольно подробно изучена методами компьютерного моделирования, в то время как численные инструменты и экспериментальные методы для описания и моделирования реальной трехмерной микроструктуры таких сложных материалов появились совсем недавно. В статье проведено компьютерное моделирование в рамках теории пластичности кристаллов для трехмерной двухфазной микроструктуры, полученной методом EBSD-томографии (3D electron backscatter diffraction). Использован метод послойного EBSD-анализа. Срезы трехмерной микроструктуры выступают в качестве двумерных моделей для изучения влияния такого упрощения на распределение напряжений и деформаций. Моделирование проводили с использованием феноменологической модели пластичности кристаллов и спектрального метода, реализованного в программном комплексе DAMASK для моделирования современных материалов.
Ключевые слова: пластичность кристаллов, двухфазная сталь, микроструктура, программный комплекс DAMASK, спектральный метод, трехмерный EBSD-анализ
Max-Planck-Institut für Eisenforschung GmbH, Düsseldorf, 40237, Germany
Институт исследования железа им. Макса Планка, Дюссельдорф, 40237, Германия
1. Introduction
bodies. Their good combination of strength and ductility is a result of the strain and stress partitioning among these two different phases and the individual crystallographic grains and subgrains [1]. Therefore, understanding how
Dual phase steels, consisting of a relatively soft ferritic matrix reinforced by strong martensitic islands, are frequently used for structural parts and reinforcements in car
© Diehl M., An D., Shanthraj P., Zaefferer S., Roters F., Raabe D., 2017
these microstructural features influence the global mechanical properties is of very high importance for the design of improved dual phase steel grades. Hence, a large number of simulation studies — starting with Karlsson and Sundström [2] in 1974 — has been dedicated to the investigation of dual phase steel micromechanics. However, even recent (crystal plasticity) simulations aiming at a better understanding of dual phase steel micromechanics are often performed on models that are significantly simplified representations of the microstructure: Either (a) 2D microstructures obtained from measurements of grain morphology and crystallo-graphic orientation [3-5] or (b) idealized 3D microstructures without crystallographic orientation information [610] and/or simplified grain morphology, phase distribution, and orientation scatter [7, 8, 11-13] are used. Both approaches have their drawbacks: Simulations using 2D microstructures reduce real materials to artificially extruded, hence columnar materials and thus leave out a significant part of the grain and phase neighborhood, imposing local boundary conditions that are quantitatively and qualitatively different than in a 3D situation [12, 14, 15]. Simplified structures obviously miss details associated with grain morphology, martensite distribution, phase and texture percolation effects, and in-grain crystallographic orientation scatter present in real dual phase microstructures.1 In the present study we, therefore, investigate the micromechanical behavior of a measured 3D dual phase microstructure to overcome these limitations. Moreover, sections of the 3D microstructure are used in the form of extruded 2D models to study how the reduction of a 3D microstructure into a 2D one influences the stress and strain distribution.
The study is structured as follows: Fist, the crystal plasticity model and the numerical method for solving the mechanical boundary value problem are presented. Afterwards, details of the microstructure model — starting from the experimental characterization — are given, followed by the presentation and discussion of the results. The study finishes with concluding remarks and an outlook on further work.
2. Simulation framework
The simulations presented in this work are performed by using the Düsseldorf Advanced Material Simulation Kit (DAMASK) available as free and open source software on https://damask.mpie.de [16]. In DAMASK, crystal plasticity is modeled within a continuum mechanics approach by taking into account the distinct deformation modes of crystalline matter. The details of the underlying continuum mechanics description, the crystal plasticity formulation, and the numerical technique to solve the boundary value problem for static mechanical equilibrium are presented in the following.
2.1. Continuum mechanical description
B denotes a body that occupies the region B0 in the reference configuration and B t in the current configuration. The location of the material points in the reference state is given by x, x e B0 and in a deformed configuration
by y, y e B t.
A deformation map /(x) : x e ®0 ^ y e Bt maps points x in the reference configuration to points y in the current configuration. The displacement u of a material point is the difference vector between these configurations:
u(x, t) = /(x, t) - x. (1)
Focusing on a point in time, i.e. a fixed deformation state, allows simplifying the notation to u (x ) = / (x) - x.
A line segment dx in an infinitesimal neighborhood of a material point x is pushed forward by
fy , /O
y + dy = y + — dx + (9(dx2). dx
(2)
Neglecting terms of higher order, dy can be expressed as
(3)
dy =-dy dx = M^dx = Vx dx, dx dx ^
=:F (x)
where F(x) is the deformation gradient and V is the "del" operator. The deformation gradient maps the vector dx at x in the reference configuration to the vector dy at y in the current configuration.
For a moving body, the position of the material points varies with time. The material velocity field is defined as2
du(x) . . v = —— = ii = x. dt
(4)
ii = x holds because the points in the reference configuration do not change their position, i.e. dx/ dt = 0. The spatial gradient of the velocity field is
L = (5)
where L is called the velocity gradient. Using the chain rule, it can be expressed as
L = F • F-1. (6)
2.2. Crystal plasticity formulation
As common in large-strain crystal plasticity modeling, the deformation gradient F is multiplicatively split up [17] into an elastic part Fe and a plastic part Fp as
F = Fe Fp (7)
to model the combined elastic-plastic response usually seen in metals.
For small strains (and short loading times), the behaviour is usually purely elastic, i.e. reversible and linear (Hookean) elasticity can be assumed. Therefore, the second Piola-Kirchhoff stress S depends on the elastic Green-Lagrange strain via the anisotropic elastic stiffness C. This relation is
1 A recent example that shows how most of these features can be included into artificial microstructures is given by Pagenkopf et al. [13].
2 To simplify the notation, in the following, the argument x is dropped whenever it is possible, i.e. F(x) is denoted as F only.
written in an artificial, only elastically deformed "intermediate configuration"1 [17]:
S = C(FeTFp -1)/2. (8)
Since C relates two symmetric tensors, it is itself symmetric and can be written as a 6x6 matrix.
The stress S is acting as a driving force for the plastic velocity gradient Lp. Lp depends on the underlying microstructure represented by a state variable vector 4 of the plasticity model and possibly other variables:
Lp = f (S, $,...) (9)
with f depending on the details of the used plasticity model which might by sensitive to strain rate, temperature, etc. as indicated by "...".
From Eq. (6) follows
Fp = Lp Fp (10)
for the evolution of the plastic deformation gradient.
The set of nonlinear Eqs. (7) to (10) needs to be solved iteratively. In DAMASK, a Newton-Raphson scheme is employed for that. More details about the implementation are given by Kords [19].
The core of each crystal plasticity model is the formulation shown in Eq. (9). In this equation, the plastic deformation of the material and its dependence to the aforementioned influences are incorporated. For the evolution of the state 4 a crystal plasticity model needs to be completed by a second equation:
4 = g(s, £,...). (11)
The two integration schemes are performed staggered, i.e. Eqs. (7) to (10) are solved at a fixed plastic state, followed by a state update, the solution of Eqs. (7) to (10), and so forth until a converged solution is achieved within the tolerance limits specified.
Using P = F S for conversion from S to P, a crystal plasticity model (i.e. Eqs. (7) to (10)) can be summarized as
P(x) = f (x, F, 4,...) (12)
which is a stress-strain relation with dependence on the material state and model dependent parameters (denoted by "..." in Eq. (12)) such as loading rate and temperature.
The crystal plasticity model used in this study is based on a phenomenological description proposed by Hutchinson [20] for face-centered cubic (fcc) crystals. Here, it has been adopted for describing the behavior of body-centered cubic (bcc) crystallites. This model is based on the assumption that plastic deformation occurs on a slip system when the resolved shear stress exceeds a critical value.
The microstructural state is parameterized in terms of resistances £ on Ns = 24 slip systems of which 12 have a
1 This is formally equivalent to the stress-strain relation in the reference configuration. The given relation is therefore a valid approximation for small elastic strains, where the difference between reference and intermediate configuration is negligible.
{110} and another 12 have a {112} plane. On both families of slip planes (111) directions are densely packed and serve as slip directions. The resistances on the k = 1, ..., Ns slip systems evolve from their initial value 40 asymptotically to a system-dependent saturation value ^ with shear ym on all slip systems according to the relationship
4* = ho Elïm
«=1
1—2-
(
sign
1—
km '
(13)
where hlm is an interaction matrix and h0 and a are system dependent fitting parameters.
Given a set of current slip resistances, a slip system shears with a rate
Y * = Y o
sign T
(14)
with Tk = S • (sk <8> n ) being (according to Schmids law [21]) the resolved shear stress for applied stress S and the unit vectors along the slip direction sk and along the slip plane normal nk. Non-Schmid contributions [22] are not taken into account.
The sum of shear rates on all systems determines finally the plastic velocity gradient
Lp = EYk sk(
(15)
k=1
2.3 Boundary value solver
For general cases, the solution for (static) equilibrium and strain compatibility under the given boundary conditions has to be sought numerically. A variety of numerical techniques exist for solving the differential equations. Most often, the finite element method is used. As an alternative solver, a spectral method using fast Fourier transforms was introduced by Moulinec and Suquet [23] in the field of material micromechanics and has gained significant attention during the last years [24-26]. Limited to periodic boundary conditions, this technique usually exceeds the finite element method in terms of solution quality because it uses trigonometric polynomials as ansatz functions. Moreover, since it operates in Fourier space, the use of fast Fourier transforms allows for a very time- and memory-efficient iterative solution algorithm. Various improvements and extensions were reported in the recent years: Lebensohn [27] extended it to the context of crystal viscoplasticity and showed the capabilities of this approach as well as several applications in a number of studies [28-30]. Crystal plasticity-based constitutive laws have also been successfully employed by, e.g., Suquet et al. [31] and Grennerat et al. [32].
Simulations of heterogeneous materials, however, are limited by the slow convergence of the original fixed-point iterative method when it is applied to materials with a large contrast in the local stiffness [33]. Several approaches have been proposed to overcome this limitation. Accelerated
72
schemes have been introduced by Eyre and Milton [34] and Monchiet and Bonnet [35] for materials with large property contrasts. Michel et al. [33] suggested a method based on augmented Lagrangians that also works in the case of materials with infinite property contrast. Using the original formulation and substituting the fixed-point method by advanced solution methods is another option to improve convergence as shown by Zeman et al. [36] and Brisard and Dormieux [37].
The spectral formulation used in this study is the one presented by Eisenlohr et al. [38] and its extension by Shanthraj et al. [39]. Since the reader might not be familiar with this method, we outline it briefly in the following.
2.4. Formulation
The deformation map x(x) is expressed as a sum of a homogeneous deformation, characterized by a constant deformation gradient F, and a superimposed deformation fluctuation field w:
X(x) = Fx + w (x), (16)
for which periodicity conditions hold, i.e. w- = w+ on corresponding surfaces dB- and dB + on B.
Equation (16) allows writing the deformation gradient F as the sum of a spatially homogeneous deformation part F and a locally fluctuating displacement part F :
F = F + F. (17)
The material response, Eq. (12), is formally written as a relation between the deformation gradient F and the first Piola-Kirchhoff stress P through a strain energy density functional W:
8W
P(x) =1W = f (x, F, $,...). (18)
oF( x)
2.5. Direct variational formulation
The direct variational formulation is closely related to the original spectral method suggested by Moulinec and Suquet [23] and the corresponding large-strain formulation introduced by Lahellec et al. [40] and Eisenlohr et al. [38]. Here it is written in a general form that allows replacing the original fixed-point approach by more powerful numerical solvers [39].
The equilibrated deformation field is obtained by minimizing W over all deformation fields that fulfill Eq. (16) for an externally prescribed average deformation. Static equilibrium expressed in real and Fourier1 space follows as
min W ^ V-P(x) = F_1[P(k)ik] = 0, (19)
x
which is equivalent to finding the root of the residual body force field
1 Quantities in real space and Fourier space are distinguished by notation Q(x) and Q(k), respectively, with x the position in real space, k the frequency vector in Fourier space, and i2 = -1. F_1[.] denotes the inverse Fourier transform.
F[x(k )]:= P (k )ik = 0. (20)
The differential Eq. (20) in Fourier space is numerically difficult to solve because of its high condition number. Introducing, in the spirit of Eshelby and Mura [41], a linear comparison material of stiffness D allows reformulation ofEq. (20) into an equivalent problem P(x) = DF(x) = = DVx with better numerical properties, i.e. a lower condition number. Equilibrium in this reference material is fulfilled if, for a given deformation map x, the residual body force field vanishes
P[x(k)] := D[x(k) ® ik ]ik = A (k )x (k) = 0. (21)
The acoustic tensor A(k) is a shorthand notation for A (k) a (k) := D[a (k) ® ik ]ik for any given vector field a(k). It corresponds to an operator on a deformation map producing the body forces resulting in the reference material. The inverse A-1 therefore gives the deformation map that would result from a known body force field in the reference material. This deformation map vanishes if the body force field vanishes, i.e., in static equilibrium for a positive-definite D. Next, an operator that results in the deformation map causing the same body force field in the reference material as a given deformation map in the original material is defined. This corresponds to a preconditioning operation of P-1 on the nonlinear operator F. P is straightforward to invert since it is local in k, with P-1 = = A(k)-1. The preconditioned system thus reads (Vk ^ 0):
P _1F(x(k)) = A(k )-1P (k )ik = 0. (22)
The deformation gradient field corresponding to this deformation map is obtained from the gradient in real space of Eq. (22)
P _1F[x(k)] <8> i k = [A(k )-1 P (k )i k ] <g> i k = 0. (23)
This is equivalent to Eq. (22) except for a constant residual field, i.e. at k = 0 where the prescribed average deformation gradient is known to hold. Expressed in terms of the deformation gradient field, Eq. (23) reads
Fdir[F (k )]:=r(k )P(k) = 0, (24)
where the Gamma operator r(k) is defined as a shorthand notation to r(k )T(k ):= [A(k )-1T(k )ik ] ® i k for a tensor field T(k).
2.6. Mixed variational formulation
The idea of using a mixed variational formulation was firstly outlined by Michel et al. [42] to overcome convergence problems with high phase contrast. Here, this idea is presented in a general form given by Shanthraj et al. [39]. In this formulation the equilibrium deformation field is obtained by minimizing W over all deformation gradient fields fulfilling Eq. (16) for an externally prescribed average deformation. Unlike in the direct variational formulation, the compatibility of the deformation gradient VxF = 0, field is imposed as an auxiliary constraint
minW (F) subject to F = V/.
(25)
This can be reformulated as an unconstrained optimization problem by introducing a Lagrange multiplier field A (x), and a penalty term. Conveniently, the reference stiffness D as defined above is used as the penalty parameter. The resulting augmented Lagrange multiplier functional reads L[F, /, A] = W + J A(x) • [V/(x) - F(x)] dx + B>
+ j [Vx(x) - F(x)] • D[Vx(x) - F(x)] dx (26)
®0 2 The equilibrium condition is equivalent to the saddle point of Eq. (26) (for details see [43]). This results in the following three stationary conditions:
8L - = P(x) - A(x) + D{F(x) - V/(x)} = 0,
D =
argmax||dP/ dF (x)||F + argmin ||dP/dF (x)||F
. (30)
8F (x)
8L
§/(x) 8£
= V[A (x)-D{F (x) - V/ (x)}] = 0,
= V/ (x) -F(x) = 0.
(27)
8A (x)
After the application of the Fourier transform, it reads as P(k) - Ao(x) + D{F(k) - /(k) ® Zk} = 0,
A-1(k){DF(k) - Ao(x)}zk = /(k), (28)
/(k) ® ik - F(k) = 0 with the Lagrange multiplier field relative to the reference configuration as A 0(x) = F(x)A (x). Eliminating / from Eqs. (28)13 using Eq. (28)2 gives the equilibrium deformation gradient field as the solution to
^mix[F(k), Ao(k)]:=
= |P(k) - Ao (k) + D(F(k) - r(k)(DF(k) - Ao (k))) ]
: |F(k) -r(k)(DF(k) -Ao(k))
= 0. (29)
The condition F(k) - r(k)(DF(k) - A0 (k)) = 0 in Eq. (29) is necessary and sufficient to enforce both the compatibility of the deformation gradient field, given by F (k) -- r(k )DF (k) = 0, and equilibrium of the Lagrange multiplier field, given by T(k) A 0(k) = 0 (for details see [35]). The condition P(k) - A0(k) = 0 enforces equilibrium of the stress field corresponding to F.
2. 7. Implementation
The hexahedral domain B0 with side lengths dx, dy, dz is discretized into a regular grid of Nx x Ny x Nz = N points with unit spacing. The solution field is approximated in the discrete Fourier space associated with this real space grid. For the transformation between real and Fourier space, the Fastest Fourier Transform in the West [44] developed by Frigo and Johnson [45] is used.
The choice of the reference stiffness D has a strong influence on stability and convergence rate as shown by Michel et al. [33]. In absence of an analytic expression for the large strain formulation the reference stiffness D is selected as
- F -
PF (x) - F
-1
(31)
A collocation-based discretization approach at the grid points in real space is used to discretize Eq. (29). The static equilibrium condition is expressed in terms of the deformation gradient field and a rescaled polarization field, FT(x) := D-1(F-1(x)A0(x)) + F(x). To get a well-conditioned system of equations, the inverse of the reference stiffness D-1 is used to scale the stress term in Eq. (29) resulting in the following expression:
F [F (x), FT (x)]:=
iD-1(P(x) - A0(x)) + pF(x) -
jT(k)(PDF(k)-aA0(k),if k * 0," [pFBC,if k = 0
" r(k)(pDF(k)-aA0(k)," - if k * 0, pFBc,if k = 0
The coefficients a and P with default value a = P = 1.0 enable to weight the conditions for static equilibrium and strain compatibility. Their optimal choice depends on the problem type, e.g. the ratio of increased or decreased stiffness compared to the average stiffness. A detailed study of their influence on the convergence rate is given by Moulinec and Silva [46] for a similar formulation in a small strain framework.
The desired deformation boundary condition FBC of the volume element is conveniently prescribed by setting AFBC = F - FBC. This allows adjusting the deformation gradient such that a stress boundary condition is component-wise as outlined by Eisenlohr et al. [38].
The nonlinear solution methods to solve the resulting systems of discretized equations are taken from PETSc [47] developed by Balay et al. [48]. The residual results from evaluating Eq. (31) at the grid points. Shanthraj et al. [39] compared three solution methods implemented in PETSc. The method showing the best performance in this benchmark, namely the nonlinear GMRES method [49], is used in this study.
2.8. Convergence criteria
A solution is accepted if the stress field is equilibrated, the deformation gradient field is compatible and the macroscopic applied boundary conditions are fulfilled.
Using the differentiation property F[df (x)/dx] = i2nk x xf (k), the deviation from static equilibrium V • P and compatibility VxF can be calculated in Fourier space. Moreover, the root mean square (RMS) as a volume average is easily accessible as discussed in Eisenlohr et al. [38] in Fourier space. The corresponding equilibrium criterion reads:
max(£e,rel
P £
I II max' e'abs
> RMS(||V- P(x)||2), (32)
Polishing number, times
Fig. 1. The surface removal per polishing step is measured on markers perpendicular to the region of interest. (a) Markers on one of the parallel sides. Dashed line 1 indicates region of interest. (b) Surface removal per polishing step and linear fit
where e rel and e bs are the relative and absolute equilibrium tolerances. The expression to test for a compatible solution is
max(ecrel|If - ill , ecabs)
--llmax c,abs > RMS( VxF(x)_), (33)
m
where e rel and e bs are the relative and absolute compatibility tolerances.
The fulfillment of complementary macroscopic deformation gradient and stress boundary conditions is determined by
max(eBC,rel ||P||max, eBC,abs) > H^BclLx (34)
with APbcij = Ajki(Fbc - F)ki, if fbcij is prescribed and APbc j = (PBC - P)j, if PBC j prescribed and the relative and absolute tolerances e bcrel and e bcabs.
3. Microstructure characterization and model creation
The investigated material is a commercial dual phase steel grade containing the main alloy elements C (0.14 wt. %), Mn (1.7 wt. %), Si (0.24 wt. %), Al (0.04 wt. %), and Cr + Mo (0.56 wt. %). The average grain size of ferrite is 5.0 ^m and the martensite islands (each consisting of several laths with different crystallographic orientations) have on average a diameter of 0.5 ^m.
In the following it is outlined how the 2D and 3D microstructural models are created from this material. Additionally, the parameters for the constitutive model and the loading boundary conditions are provided.
3.1. Sample preparation and microstructural characterization
For obtaining a three-dimensional microstructure data set of the investigated material, a serial-sectioning approach that consists of multiple cycles of mechanical polishing and electron backscatter diffraction imaging [50] was employed. To this end, first a sample with two parallel sides of high surface quality was prepared by mechanical grinding and polishing. On each of those parallel surfaces, four marker lines were engraved by focused ion beam milling using a
Zeiss Cross-Beam instrument (Fig. l, a). Two of these lines are parallel to each other and serve as a reference for magnification and positioning. The two other inclined lines have including angles of 45° with the reference lines and serve as markers to determine the polishing depth. In fact, after mounting the material and polishing the surface of the region of interest perpendicular to the marked surfaces, the marker lines appear as sharply peaked spikes at the edge of the sample. The change of distance between the depth markers and the reference markers then corresponds directly to the depth of polishing. In the present case the sample was warm mounted into conductive resin (Polyfast by Struers) into a cylindric mount of 25mm diameter.
For each measurement cycle, the sample was polished with silicon-oxide particle suspension produced by Struers at a pressure of 60 kPa for 4 min. This led to a removal of, in average, 0.13 ^m per slice as displayed in Fig. 1, b. Figure 1, b moreover shows that the removal rate was approximately constant and varied by approximately 50 nm per slice. After every polishing step the sample was manually placed into a Zeiss XB 1540 high resolution scanning electron microscope. The sample was tilted to 70° and then manually moved to the original measurement position using the reference markers. Finally, an orientation map was taken using EBSD (OIM 6.x by EDAX/TSL using a Hikari camera) operating at 15 kV electron beam acceleration voltage. One of these acquired orientation maps is shown in Fig. 2.
In total, 22 slices, each sized 20x70 ^m2 (transverse x rolling direction) at an in-plane EBSD step size of 0.2 ^m
Fig. 2. Orientation map (colors according to inverse pole figure parallel to the normal direction) overlaid with the image quality (darker values indicate lower image quality, i.e. typically martensite). Rolling direction is horizontal
Fig. 3. 3D model constructed from 22 sections where the markers can clearly be seen. Color code displays image quality. A lower image quality (dark) indicates typically martensite. Arrow 1 — rolling direction, arrow 2 — transverse direction, arrow 3 — normal direction
have been acquired. The CUBE software from Bruker [51] has been used to create a 3D model from all sections. The sample 20.0x70.0x2.86 ^m3 discretized by 351x100x22 = = 772200 points is shown in Fig. 3.
3.2. Model creation
The model required for the crystal plasticity simulations to be conducted using the spectral solver can be directly created from the measured orientation data on the regular 3D voxel grid. However, due to the similar lattice constants, it is not directly possible to distinguish martensite and ferrite. The conventional approach of using the image quality as an indicator [5, 52] (which is lower for the heavily distorted rnartensitic crystals) is not applicable here as polishing and sample preparation introduced a substantial variation in image quality between the sections and even within one section. Therefore, the grain size was used to indicate whether a measured point pertains to the martensitic or to the ferritic phase. This approach is based on the ratio that (i) ferritic grains are 10 times larger than martesitic islands and (ii) in each martensite grain, several lath of distinct orientation can be identified as individual grains [52-55]. To this end, first each orientation map was cleaned using the features available within the TSL OIM software. After cleaning, a 3D model was created on which individual grains (i.e., in the case of martensite, laths) have been identified using a tight tolerance angle of 3°. Now, all grains of a size smaller than 30 voxels have been considered as martensite. As a final step, all grains identified as martensite but con-
sisting of only one voxel are considered and treated as a measurement artifact and are hence assigned with ferrite properties to avoid associated numerical artifacts. This leads to a volume fraction of 5.8 % of martensite. A visual inspection of the segmentation achieved by this approach indicates a reasonable phase separation. However, few ferrite grains are erroneously detected to be martensite because they are smaller then the size threshold set for the phase separation. This issue on grain segmentation in three dimensions is directly related to the serial sectioning approach: The crystallographic orientations obtained from the individual EBSD measurements are mapped in a reference system relative to the assumed position of the sample on the microscope stage. Despite careful sample preparation and alignment, the position of the sample with respect to the EBSD detector differs between the individual measurements. Therefore, orientation differences significantly larger than the EBSD angular resolution of less than 2° (which is achieved in each individual map) appear between points belonging to the same grain but acquired in different sections. Moreover, at the surfaces of the measurement area where grains are cut off, detected grain sizes are smaller than in the volume and the wrong assignment occurs slightly more frequent than inside of the measurement block.
The three-dimensional model constructed from all 772200 measured points has been mirrored in all three directions to avoid artifacts at the sample boundaries introduced by the periodic repetition due to the approximation of the deformation gradient and stress fields by a Fourier series. Additionally, from each of the original 22 slices, a 2D model (still mirrored in both in-plane directions) is created to investigate the influence of a simplified columnar grain structure on the stress and strain partitioning in the spirit of Diehl et al. [12]. The martensite volume fraction ranges from 1.8 to 8.9 % in the 2D models.
Material parameters that were before fitted to a similar material by Tasan et al. [5] are used in this study (Table 1).
Uniaxial tensile loading along the rolling direction at an engineering strain rate of 1 • 10-3 s-1 was applied for 200 s in 400 increments, i.e. a final true strain of e = 0.182 was reached. This load case reflects the situation in a tensile test and has little deformation constraints compared to, e.g., plane stress loading. It should be noted that the values are prescribed as volume average for a periodically repeated body (compare Sect. 2.3) and hence differ from typical (non-periodic) boundary conditions used in finite element analysis.
Table 1
Constitutive parameters of ferrite and martensite used for the simulations taken from [5]
Property Cii, GPa C12' GPa C44, GPa Y 0' s-1 £0'{110}' MPa ^{110}' MPa ^0,{112}' MPa ^,{112}' MPa h0' GPa hkm n a
Ferrite 233.3 235.5 128.0 0.001 95 222 96 412 1 1 20 2.25
Martensite 417.4 242.4 211.1 406 873 457 971 563
Fig. 4. Equivalent strain EvM (a) and stress ovM (b) at increasing strain level e from top to bottom mapped onto the deformed configuration for the case of the 3D simulation. The surface shown is parallel to the normal plane. The black frame indicates the shape in the undeformed configuration. A logarithmic scale is used to for mapping the stress and strain values to the color bar in order to resolve details in both phases with their vastly differing mechanical properties
4. Results
4.1. Three-dimensional model
The spatial distribution of the equivalent (von Mises) strain evM and stress ctvM is shown in Fig. 4 for three different loading states, i.e. at e = 0.020, 0.095, and 0.182, where e denotes the logarithmic (true) strain along loading direction. A significant stress and strain partitioning among the phases can be seen at each load level. The spatial distribution over time, i.e. from top to bottom in Fig. 4, is largely constant as the same locations at e = 0.020, 0.095, and 0.182 can be identified as being low respectively high stressed/ strained. Among different ferrite grains and within the same grain, a large variation in stress can be observed at all three load levels (Fig. 4, b). The extreme values of stress and
strain at e = 0.182 in both phases are given together with the corresponding mean value in Table 2. As expected, the stress in martensite is much higher than in ferrite while the opposite holds for the strain. The strain in martensite between different regions differs by more than one order of magnitude and the stress by almost three orders of magnitude. In ferrite, the differences are slightly lower, i.e. by almost one order of magnitude in strain and by more than two orders of magnitude in stress.
4.2. Comparison between the three-dimensional model and the two-dimensional models
The global stress-strain curves in terms of the average Cauchy stress and logarithmic strain in loading direction
Table 2
Minimum, mean, and maximum values of the equivalent stress and strain found in the 3D simulation at a global strain in loading direction of e = 0.182
Phase Equivalent stress ovM, GPa Equivalent strain evM
Minimum Mean Maximum Minimum Mean Maximum
Ferrite 0.17 0.67 1.32 0.0026 0.2150 0.8648
Martensite 0.32 1.91 7.03 0.0008 0.0454 0.6905
Table 3
Minimum, mean, and maximum values of the equivalent stress and strain found in the 2D simulations at a global strain in loading direction of e = 0.182
Phase Equivalent stress ovM, GPa Equivalent strain evM
Minimum Mean Maximum Minimum Mean Maximum
Ferrite 0.21 0.67 1.29 0.0322 0.2195 1.2503
Martensite 1.39 2.07 11.40 0.0291 0.0465 0.4851
ct, MPa -6004002000 0.00 0.04 0.08 0.12 0.16 e
Fig. 5. Stress-strain response in loading direction obtained by the different simulations. The values of the 3D simulation is given by the solid line, the average value of the 2D simulations is given by the dashed line and the grey area indicates minimum and maximum values among all 2D simulations. Cauchy stress and true strain are computed from F and P
for the 3D simulation together with the average and the extrema of the 2D individual simulations are given in Fig. 5. The average response of all 2D simulations, i.e. using all points in the 3D volume and doing individual simulations for each section, is almost the same as the response for the full model. A difference between the softest and the stiffest section of approximately 70 MPa can be observed.
Table 3 gives the minimum, mean, and maximum values of the equivalent strain and stress obtained by the 2D simulations for a global strain of e = 0.182. A comparison of the stress values (Table 3) to the results obtained from the 3D simulation (given in Table 2) reveals that the 2D simulation predicts a more homogeneous distribution around the same average in ferrite while minimum, average, and maximum values in the martensitic phase are substantially higher in the 2D situation than in their 3D counterpart. The mean values of the equivalent strain show no significant difference between the 2D assumption and the 3D model (compare Table 3 with Table 2). A wider distribution of the strain values in martensite results from the 3D simulation while in ferrite the minimum and maximum values are significantly higher in the 2D case.
In Fig. 6, probability density maps ("heat maps") are given for comparison of the stress-strain partitioning in the 3D simulation (left) to the sum of all 2D simulations (right). More specifically, Fig. 6, a displays the probability density of equivalent strain and stress for ferrite and martensite together. Fig. 6, b shows it for ferrite only, and Fig. 6, c for martensite only. A bipartitioned distribution can be seen in Fig. 6, a: At the bottom (low stress) the values belonging to the ferrite are located and at low strains (horizontal center) the values of the martensitic phase. This representation confirms the findings that the 2D situation leads to higher strain and stress values in martensite while the ferrite is more strained, but stays at a similar stress level. Figures 6, b and 6, c showing the heat maps for both phases individually allow to investigate behavior of martensite and ferrite more closely. From Fig. 6, b the significant different stress levels
in ferrite can be accessed. Figure 6, c reveals that in the 2D simulations a substantial share of martensite shows stress values above the maximum value found in the 3D simulation.
A point-to-point comparison of the equivalent strain and stress maps for one section in the center of the 3D microstructure is given in Fig. 7. It can be observed that the stress in martensite is higher in the 2D situation, a trend perceived also from the comparison of Table 3 with Table 2. To visualize the differences between the results obtained by both simulations, first the ratio between the shown values is computed followed by taking the natural logarithm of it: ln( P2d/ Рзб) for P e CTvM > evM- This measure is motivated from the true (logarithmic) strain measure. Its desired property is the asymptotic behavior, lim ln(p2D/p3D) = ±«>,
in contrast to the relative difference which has lower bound of -1 and upper bound of+«>. Hence, this measure is better suited then a "standard" relative difference when large deviations to both, lower and higher values, occur. Maps showing the obtained values for stress and strain are given in the bottom row of Fig. 7. It can be seen that the relative differences in strain are especially large in or in the vicinity of martensite, while the largest relative stress difference is often found in ferrite.
OvM, GPa ctvM, GPa
0.5 0.3 0.1 c 0.1 0.3 0.5
WM
Fig. 6. Probability density maps showing the correlation between equivalent stress and strain within the 3D simulation (left) and all 2D simulations (right). Darker value indicates higher probability to find the stress-strain correlation. Ferrite and martensite (a), ferrite only (b) and martensite only (c)
.00 0.04 0.08 0.12 0.16
0.0009
P = ^vM -2.5 +2.5
OvM, GPa 0.25 2.5
b
? s - t > * f -.*'>»
k r i ' "N i : - ,r * r u i- " C ».« 's . ■ J, ' V-'ï + r • - ■ . - .-
P = °vM -1.0 +1.0
v. £ :
Fig. 7. Equivalent strain evM (a) and stress ovM (b) in a section of the 3D microstructure (top) and in the corresponding 2D microstructure (middle) and the logarithmic of the ratio between them (bottom) mapped onto the deformed configuration of the 3D simulation (e = = 0.182). The section shown is parallel to the normal plane. A logarithmic scale is used for mapping the stress and strain values in order to resolve details in both phases with their vastly differing mechanical properties
5. Discussion
The results presented above clearly emphasize the severe strain and stress partitioning occurring in dual phase steels caused by the very different mechanical behavior of martensite and ferrite. The magnitude and pattern of the partitioning depends on details such as martensite shape and distribution, crystallographic orientation and misori-entation that have been included into the simulation by using an experimentally characterized microstructure as input data.
The global stress-strain response of a 3D model and the 2D microstructures built from this 3D models differs only slightly. This is especially astonishing as the martensite volume fraction of the individual sections differs significantly from their average given by the 3D model. In contrast, the local response is vastly changed when a 2D section is used alone compared to embedding it into the experimentally obtained neighborhood. In fact, the 2D simulations show strain hot spots which are up to 11 times higher in strain and up to 3 times higher in stress than the same positions in the 3D simulations. The reason for this is that the columnar structure of the 2D simulation does not allow for relaxation or support through different deformation paths above or below the observed structure. As a consequence all strain paths follow a 2D-pattern resulting in the observed intense shear patterns in the 2D simulations. In contrast, the 3D simulations allow strain paths in 3 dimensions, thus resulting in lower stress and strain concentrations. These observations show that the analysis of local strain features, for example for the understanding of crack nucleation, imperatively need 3D calculations with realistic microstructure input.
6. Conclusion and outlook
A dual phase steel microstructure model based directly on a data set obtained by a 3D EBSD characterization approach has been utilized for a simulation study using a crystal plasticity approach in combination with a fast spectral solver included in the free and open source software DAMASK. The microstructural model, by including details such as martensite shape and distribution, crystallographic orientation and misorientation, allows to investigate stress and strain partitioning without significant simplifications. The proposed modeling method and the approach for obtaining the microstructure by serial sectioning are therefore a promising route for the efficient investigation of stress and strain partitioning in dual phase steels and other materials with complex microstructures.
2D models, assuming a columnar grain structure, have been used to investigate the influence of this simplification. In agreement with Diehl et al. [12] it is clearly shown how important full 3D microstructures are for capturing the details of stress and strain partitioning in dual phase steels for the investigated loading conditions. While no large effect on the global response could be found, the local stress and strain values differ quantitatively and qualitatively when a 2D model is used instead of a 3D model. Especially for simulations tackling the question of damage initiation in dual phase steels (for damage models implemented in DAMASK see [57]), 2D modelling is therefore an unsuitable and erroneous approach: As these damage models predict damage based on the local stress or strain response, relying on values that are off by a factor of up to 3 for the stress and by a factor of up to 11 for the strain when employing a 2D model would in most cases result in entirely incorrect simulation results and conclusions.
For building the model of dual phase steel, the distinct separation of both phases remains a challenge. It can be expected that replacing the manual serial sectioning approach presented here by an automated system will help to solve this problem. First, a more reproducible surface preparation will restore the possibility to use the image quality as a phase indicator. Second, a more exact placement of the sample in the EBSD system will reduce the deviations in measured crystallographic orientation between different acquisition runs. Combining both indicators, i.e. grain size and (grain average) image quality into one algorithm should then allow to reliably separate martensite and ferrite. Moreover, an automated system will reduce experimental efforts and enable to investigate larger volumes and improve statistical confidence.
Acknowledgments
M. Diehl acknowledges the funding of the TCMPrecip-Steel project in the framework of the SPP 1713 (Strong coupling of thermo-chemical and thermo-mechanical states in applied materials) by the Deutsche Forschungsgemeinschaft (DFG).
References
1. Tasan C.C., Diehl M., Yan D., Bechtold M, Roters F., Schemmann L., Zheng C., Peranio N., Ponge D., Koyama M., Tsuzaki K., Raabe D. An overview of dual-phase steels: Advances in microstructure-ori-ented processing and micromechanically guided design // Annu. Rev. Mater. Res. - 2015. - V. 45. - P. 391-431.
2. Karlsson B., Sundström B. O. Inhomogeneity in plastic deformation of two-phase steels // Mater. Sci. Eng. - 1974. - V. 16(1-2). - P. 161168.
3. Kadkhodapour J., Butz A., Ziaei-Rad S., Schmauder S. A microme-chanical study on failure initiation of dual phase steels under tension using single crystal plasticity model // Int. J. Plasticity. - 2011. -V. 27(7). - P. 1103-1125.
4. Tasan C.C., Hoefnagels J.P.M., Diehl M., Yan D, Roters F., Raabe D. Strain localization and damage in dual phase steels investigated by coupled in-situ deformation experiments-crystal plasticity simulations // Int. J. Plasticity. - 2014. - V 63. - P. 198-210.
5. Tasan C.C., Diehl M., Yan D., Zambaldi C., Shanthraj P., Roters F., Raabe D. Integrated experimental-numerical analysis of stress and strain partitioning in multi-phase alloys // Acta Mater. - 2014. - V. 81.-P. 386-400.
6. Scherff F., Goldschmidt F., Scholl S., Diebels S. High-resolution simulation of microstructures in dual-phase steel // PAMM. - 2016. -V. 16(1). - P. 391-392.
7. Uthaisangsuk V., Prahl U., Bleck W. Stretch-flangeability characterization of multiphase steel using a microstructure based failure modeling // Comp. Mater. Sci. - 2009. - V. 45(3). - P. 617-623.
8. Fillafer A., Krempaszky C., Werner E. On strain partitioning and micro-damage behavior of dual-phase steels // Mater. Sci. Eng. A. - 2014. -V. 614. - P. 180-192.
9. Scheunemann L., Balzani D., Brands D., Schröder J., Raabe D. Statistically similar RVE construction based on 3D dual-phase steel microstructures // Proc. 5th Int. Conf. Structural Engineering, Mechanics and Computation. - Boca Raton: CRC Press, 2013. - P. 411-416.
10. Brands D., Balzani D., Scheunemann L., Schröder J., Richter H., Raabe D. Computational modeling of dual-phase steels based on representative three-dimensional microstructures obtained from EBSD data // Arch. Appl. Mech. - 2016. - V 86(3). - P. 575-598.
11. Woo W, Em V.T., Kim E-Y, Han S.H., Han Y.S., Choi S-H Stressstrain relationship between ferrite and martensite in a dual-phase steel studied by in situ neutron diffraction and crystal plasticity theories // Acta Mater. - 2012. - V. 60(20). - P. 6972-6981.
12. Diehl M., Shanthraj P., Eisenlohr P., Roters F. Neighborhood influences on stress and strain partitioning in dual-phase microstructures. An investigation on synthetic polycrystals with a robust spectral-based numerical method // Meccanica. - 2016. - V. 51(2). - P. 429-441.
13. Pagenkopf J., Butz A., Wenk M., Helm D. Virtual testing of dualphase steels: Effect of martensite morphology on plastic flow behavior // Mater. Sci. Eng. A. - 2016. - V. 674. - P. 672-686.
14. Zeghadi A., N'guyen F., Forest S., Gourgues A-F, Bouaziz O. Ensemble averaging stress-strain fields in polycrystalline aggregates with a constrained surface microstructure. Part 1: Anisotropic elastic behavior // Philos. Mag. - 2007. - V 87(8-9). - P. 1401-1424.
15. Zeghadi A., Forest S., Gourgues A-F, Bouaziz O. Ensemble averaging stress-strain fields in polycrystalline aggregates with a constrained surface microstructure. Part 2: Crystal plasticity // Philos. Mag. -V. 87(8-9).
16. Roters F., Eisenlohr P., Kords C., Tjahjanto D.D., DiehlM., RaabeD. DAMASK: the Dusseldorf Advanced Material Simulation Kit for Studying Crystal Plasticity Using an FE Based or a Spectral Numerical Solver // Proc. IUTAM: From Microstructure to Macroscale Properties / Ed. by O. Cazacu. - Amsterdam: Elsevier, 2012. - V. 3. -P. 3-10.
17. Reina C., Conti S. Kinematic description of crystal plasticity in the finite kinematic framework: A micromechanical understanding of F = FeFp // J. Mech. Phys. Solids. - 2014. - V. 67. - P. 40-61.
18. Roters F., Eisenlohr P., Bieler T.R., Raabe D. Crystal Plasticity Finite Element Methods in Materials Science and Engineering. - Weinheim: Wiley-VCH, 2010.
19. Kords C. On the Role of Dislocation Transport in the Constitutive Description of Crystal Plasticity: PhD Thesis. - Berlin: RWTH Aachen, 2013. - http://darwin.bth.rwth-aachen.de/opus3/volltexte/2014/4862.
20. Hutchinson J.W. Bounds and self-consistent estimates for creep of polycrystalline materials // Proc. R. Soc. A. - 1976. - V. 348. - P.101-127.
21. Gottstein G. Physical Foundations of Material Science. - Berlin: Springer, 2004.
22. Gröger R., Vitek V. Breakdown of the Schmid law in bcc molybdenum related to the effect of shear stress perpendicular to the slip direction // Mat. Sci. Forum. - 2005. - V. 482. - P. 123-126.
23. Moulinec H., Suquet P. A fast numerical method for computing the linear and nonlinear properties of composites // Comptes Rendus de l'Academie des Sciences. Ser. II. Mecanique, Physique, Chimie, Astronomie. - 1994. - V. 318. - P. 1417-1423.
24. Schneider M., Ospald F., Kabel M. Computational homogenization of elasticity on a staggered grid // Int. J. Numer. Meth. Eng. - 2015. -V. 105(9). - P. 693-720.
25. Kabel M., Fliegener S., Schneider M. Mixed boundary conditions for FFT-based homogenization at finite strains // Comp. Mater. Sci. -2016. - V. 57(2). - P. 193-210.
26. Zeman J., de Geus T.W.J., Vondrejc J., Peerlings R.H.J., Geers M.G.D. A finite element perspective on non-linear FFT-based micromechanical simulations // Int. J. Numer. Meth. Eng. - 2017. - doi 10.1002/ nme.5481.
27. Lebensohn R.A. N-site modeling of a 3D viscoplastic polycrystal using fast Fourier transform // Acta Mater. - 2001. - V. 49(14). - P. 27232737.
28. Lebensohn R.A., Castelnau O., Brenner R., Gilormini P. Study of the antiplane deformation of linear 2D polycrystals with different microstructures // Int. J. Solids. Struct. - 2005. - V 42(20). - P. 5441-5459.
29. Lefebvre G., Sinclair C.W., Lebensohn R.A., Mithieux J.-D. Accounting for local interactions in the prediction of roping of ferritic stainless steel sheets // Model. Simul. Mater. Sci. Eng. - 2012. - V. 20(2). -P. 024008.
30. Lebensohn R.A., Kanjarla A.K., Eisenlohr P. An elasto-viscoplastic formulation based on fast Fourier transforms for the prediction of
micromechanical fields in polycrystalline materials // Int. J. Plasticity. - 2012. - V. 32-33. - P. 59-69.
31. Suquet P., Moulinec H., Castelnau O., Montagnat M., Lahellec N., Grennerat F., Duval P., Brenner R. Mult-Scale Modeling of the Mechanical Behavior of Polycrystalline Ice under Transient Creep // Procedia IUTAM. - 2012. - V.3. - P. 76-90.
32. Grennerat F., Montagnat M., Castelnau O., Vacher P., Moulinec H., Suquet P., Duval P. Experimental characterization of the intragranular strain field in columnar ice during transient creep // Acta Mater. -2012. - V. 60(8). - P. 3655-3666.
33. Michel J.C., Moulinec H., Suquet P. A computational scheme for linear and non-linear composites with arbitrary phase contrast // Int. J. Numer. Meth. Eng. - 2001. - V. 52(12). - P. 139-160.
34. Eyre D.J., Milton G.W. A fast numerical scheme for computing the response of composites using grid refinement // Eur. Phys. J. Appl. Phys. - 1999. - V. 6(1). - P. 41-47.
35. Monchiet V., Bonnet G. A polarization-based FFT iterative scheme for computing the effective properties of elastic composites with arbitrary contrast // Int. J. Numer. Meth. Eng. - 2012. - V. 89(11). -P. 1419-1436.
36. Zeman J., Vondrejc J., Novak J., Marek I. Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients // J. Comput. Phys. - 2010. - V. 229(21). - P. 8065-8071.
37. Brisard S., Dormieux L. FFT-based methods for the mechanics of composites: A general variational framework // Comp. Mater. Sci. -2010. - V. 49(3). - P. 663-671.
38. Eisenlohr P., Diehl M., Lebensohn R.A., Roters F. A spectral method solution to crystal elasto-viscoplasticity at finite strains // Int. J. Plasticity. - 2013. - V. 46. - P. 37-53.
39. Shanthraj P., Eisenlohr P., Diehl M., Roters F. Numerically robust spectral methods for crystal plasticity simulations of heterogeneous materials // Int. J. Plasticity. - 2015. - V. 66. - P. 31-45.
40. Lahellec N., Michel J.C., Moulinec H., Suquet P. Analysis ofInhomo-geneous Materials at Large Strains Using Fast Fourier Transforms // Solid Mechanics and Its Applications: IUTAM Symposium on Computational Mechanics of Solid Materials at Large Strains / Ed. by C. Miehe. - Dordrecht: Kluwer Academic Publishers, 2001. - V. 108. -P. 247-258.
41. Mura T. Micromechanics of Defects in Solids. - Dordrecht: Martinus Nijhoff Publishers, 1987.
42. Michel J.C., Moulinec H., Suquet P. A computational method based on augmented Lagrangians and fast Fourier transforms for composites with high contrast // Comput. Model. Eng. Sci. - 2000. - V. 1(2). -P. 79-88.
43. Augmented Lagrangian Methods: Applications to the Numerical Solution of Boundary-Value Problems / Ed. by M. Fortin, R. Glowinski. -1983.
44. Frigo M., Johnson S.G. FFTW, 2014. - http://www.fftw.org.
45. Frigo M., Johnson S.G. The design and implementation of FFTW3 // Proc. IEEE. - 2005. - V. 93(2). - P. 216-231.
46. Moulinec H., Silva F. Comparison of three accelerated FFT-based schemes for computing the mechanical response of composite materials // Int. J. Numer. Meth. Eng. - 2014. - V. 97(13). - P. 960-985.
47. PETSc Team. - 2015. - http://www.mcs.anl.gov/petsc.
48. Balay S., Brown J., Buschelman K., Eijkhout V., Gropp W.D., Kaus-hikD., Knepley M.G., Mclnnes L.C., Smith B.F., Zhang H. PETSc Users Manual: Technical Report. - 2013.
49. Oosterlee C.W., Washio T. Krylov subspace acceleration of nonlinear multigrid with application to recirculating flows // SIAM J. Sci. Comput. - 2000. - V. 21(5). - P. 1670-1690.
50. Zaefferer S., Wright S.I., Raabe D. Three-dimensional orientation microscopy in a focused ion beam-scanning electron microscope: A new dimension of microstructure characterization // Metall. Mater. Trans. A. - 2008. - V. 39(2). - P. 374-389.
51. Konijnenberg P.J., Zaefferer S., Lee S.B., Rollett A.D., Rohrer G.S., Raabe D. Advanced methods and tools for reconstruction and analysis of grain boundaries from 3D-EBSD data sets // Mat. Sci. Forum. -2011. - V. 702-703. - P. 475-478.
52. Schemmann L., Zaefferer S., Raabe D., Friedel F., Mattissen D. Alloying effects on microstructure formation of dual phase steels // Acta Mater. - 2015. - V. 95. - P. 386-398.
53. Morito S., Tanaka H., Konishi R., Furuhara T., Maki T. The morphology and crystallography of lath martensite in Fe-C alloys // Acta Mater. - 2003. - V. 51(6). - P. 1789-1799.
54. Calcagnotto M., Adachi Y., Ponge D., Raabe D. Deformation and fracture mechanisms in fine- and ultrafinegrained ferrite/martensite dual-phase steels and the effect of aging // Acta Mater. - 2011. -V. 59(2). - P. 658-670.
55. Yan D., Tasan C.C., Raabe D. High resolution in situ mapping of microstrain and microstructure evolution reveals damage resistance criteria in dual phase steels // Acta Mater. - 2015. - V. 96. - P. 399409.
56. Shanthraj P., Sharma L., Svendsen B., Roters F., Raabe D. A phase field model for damage in elasto-viscoplastic materials // Comput. Meth. Appl. Mech. Eng. - 2016. - V. 312. - P. 167-185.
57. Shanthraj P., Svendsen B., Sharma L., Roters F., Raabe D. Elasto-viscoplastic phase field modeling of anisotropic cleavage fracture // J. Mech. Phys. Solids. - 2017. - V. 99. - P. 19-34.
nocTynma b peaaKUHro 24.01.2017 r.
CeedeHun 06 aemopax
Martin Diehl, Dr.-Ing., PhD Res., Max-Planck-Institut fur Eisenforschung GmbH, Germany, [email protected] Dayong An, MSc, PhD Stud., Max-Planck-Institut fur Eisenforschung GmbH, Germany, [email protected] Pratheek Shanthraj, Dr., PhD Res., Max-Planck-Institut fur Eisenforschung GmbH, Germany, [email protected] Stefan Zaefferer, PD Dr., Group Leader, Max-Planck-Institut fur Eisenforschung GmbH, Germany, [email protected] Franz Roters, PD Dr., Group Leader, Max-Planck-Institut fur Eisenforschung GmbH, Germany, [email protected] Dierk Raabe, Prof. Dr.-Ing. habil., Director, Max-Planck-Institut für Eisenforschung GmbH, Germany, [email protected]
