Multiscale modeling of the mechanical behavior of brazed Ni-based superalloy sheet metals

The elastic–plastic behavior of brazed Ni-based superalloys used in abradable turbomachinery sealing systems is analyzed by means of numerical simulations. A sequential multiscale modeling approach is employed to analyze the relevant effects on the mechanical behavior of a layered composite consisting of the braze metal and the joining partners. The focus of the investigations lies on the role of the multiphase microstructure within the brazing layer and the significance of microscopic parameters such as, for example, the volume fraction of the phases compared to macroscopic parameters (e.g., the brazing layer thickness). A representative volume element is employed on the microscale, and a layered composite is modeled on the macroscale to capture all relevant effects on both length scales during mechanical loading of the layered material composite. Virtual tensile tests at different temperatures and strain rates are chosen as a controlled testing environment that captures the characteristic loading conditions during a rubbing event in a turbomachinery. The parameter on the microscale with the highest influence on the mechanical behavior is the volume ratio of the brittle and ductile phases. It also has a significant effect on the macroscopic mechanical behavior. To reduce the risk of damage, it is advised to minimize the fraction of brittle phases in the brazing layer. According to the simulation results, this can be even more effective than reducing the overall thickness of the brazing layer. Additionally, an improvement of an existing analytical model for the estimation of the flow stress in the ductile phase of a dual-phase microstructure is proposed. By increasing the order of the root function in the analytic model, the effect of stress distribution between brittle and ductile phases can be incorporated in an empirical manner. This reduces the deviation between the analytical and the numerical approaches significantly.

formed in the joining process can also deteriorate the mechanical properties of the material composite. An example where the brazing process has a particularly detrimental effect is abradable sealing systems in a turbomachinery. Honeycomb liners are typically brazed to a carrier plate and form the stator of the sealing system. Sealing fins attached to the rotor rub into the honeycomb liners either in order to achieve a minimal leakage air flow or due to a change of the operation conditions [1]. Due to capillary action, the brazing material can rise from the carrier plate up to the rubbing surface and form a layered composite consisting of the base metal and the braze filler metal (see Fig. 1, left). Brittle phases that form during brazing pose a high risk of rotor damage as they show strong mechanical resistance up to high temperatures. Additionally, brittle phases can act as an accelerator for critical stator damage such as cracking and large breakouts in the honeycomb liners, which can significantly increase the leakage flow and can also damage other parts of the turbomachinery.
The microstructure of the brazing layer between the joining partners shows different characteristics depending on the duration of brazing and the braze gap size [2][3][4]. A common feature is the coexistence of brittle and ductile phases in the brazing zone. Ulan kyzy et al. [5] identified three major microstructural constituents in the double-foiled segments of a Hastelloy X honeycomb liner brazed with BNi-5 filler metal. Those were a matrix of γ -nickel solid solution (1), NiSi (2) and several silicides of Ni y (Cr,Mo) x Si type (3). A similar microstructure could be identified in Haynes 214 sheet metals brazed with BNi-5 (see contrast of back-scattered electrons in Fig. 1, right). Previous studies compared Hastelloy X with Haynes 214 as honeycomb liner materials [6,7]. It was shown that the latter alloy leads to higher mechanical and thermal loads on the rotor during a rubbing event. For this reason, we focus on Haynes 214 instead of Hastelloy X in this study.
According to Hawk et al. [3], the matrix (1) has a microstructure similar to the base metal Haynes 214 consisting of γ -nickel and nanometer-sized γ -precipitates. Measurements of microhardness showed a strong local increase in hardness in the Si-rich phase (2) and the silicides (3) [3,5]. Table 1 shows the chemical composition of two base metals typically used for the honeycomb foils as well as the braze metal  Since in the literature mainly macroscopic mechanical testing has been performed as, e.g., by Hawk et al. [3] as well as by Li and Schulson [8], the effect of the microstructure on the mechanical behavior of brazed joints is not fully understood. Isolated mechanical testing of the identified phases is cumbersome due to the difficulty of capturing the conditions in the brazing layer that stem from the process of brazing. Numerical simulations can help to overcome this challenge. However, current modeling approaches treat the brazing Fig. 1 Left: schematic top view of a honeycomb double-foiled segment with capillary braze, right: exemplary microstructure (contrast of back-scattered electrons) of the brazing layer in a honeycomb double-foiled segment with phases similar to those identified by Ulan kyzy et al. [5] for brazed Hastelloy X sheets Table 1 Chemical composition of two common base metals used for honeycomb foils as well as the braze metal BNi-5 in mass-%, balance is nickel [5] Alloy  1 9 . 0 0 -----1 0 .00 --zone as either purely ductile or purely brittle depicted by an elastic material behavior while neglecting the multiphase microstructure in the brazing layer [7,9].

Scope and goals of this work
In this work, numerical simulations are conducted to identify the most influential effects on the elastic-plastic behavior of the layered composite on the macroscale as well as on the microscale. The Ni-based superalloy Haynes 214 is chosen as base metal and BNi-5 as braze filler metal. Virtual high-temperature tensile tests at different strain rates will serve as a model basis for the simulative studies. They depict the characteristics of the loading conditions during a rubbing event in a controlled environment that will allow a quantitative comparison between the simulations and future experiments. The goal of this work is twofold. First, the necessary depth of modeling and crucial parameters for experimental characterization of brazed Ni-based sheet metals are identified and validated. Second, guidelines for a beneficial microstructure as well as a beneficial macroscopic geometry of base metal-braze composites in honeycomb liners will be derived. In particular, this means 1. lowering the resistance against the rubbing partner by reducing the yield limit as well as strain hardening, and 2. decreasing the risk of critical damage in the honeycomb liner by lowering stress peaks on both the microscale and the macroscale considering the type and location of damage. This work is organized as follows: In Sect. 2, the constitutive relations employed on each scale are described. A sequential multiscale simulation setup is presented in Sect. 3. The employed micromechanical model is validated with experimental results from the literature in Sect. 4. Sensitivity studies on both scales are discussed in Sect. 5. The work is closed in Sect. 6 with the most important conclusions as well as an outlook on further studies.

Constitutive models used on different scales
On each scale, the relevant mechanisms influencing the mechanical material behavior should be incorporated. In order to capture the effective material behavior and to assess potential micro-damage, the brazing layer is microstructurally resolved. The base metal is only considered on the macroscale as an isotropic elastic-plastic material with its uniaxial flow curve taken directly from experiment [10].

Microscale
A phenomenological elastic-viscoplastic material law considering the face-centered cubic (fcc) crystal structure is used for the ductile nickel solid solution phase (1). The generalized Hooke's law is used to calculate the second Piola-Kirchhoff stress tensor in the plastic configuration S p caused by elastic deformation F e depicted by the Green-Lagrange strain tensor E = 1/2(F T e F e − I) with I being the second rank identity tensor, In Eq. (1), C is the fourth-order elastic stiffness tensor inheriting the symmetries of the underlying crystal lattice. For the cubic crystal of the nickel solid solution phase, three independent elastic constants are required (C 11 , C 12 and C 44 ) [11]. In the framework of large deformations, the total deformation gradient F is multiplicatively split into a lattice-preserving plastic deformation step followed by a lattice distorting elastic deformation step, F = F e F p . Plastic deformation is depicted by the plastic deformation gradient F p . Its evolution is given by the single dot product of the plastic velocity gradient L p and the plastic deformation gradienṫ For crystal plasticity, it is assumed that the plastic deformation follows the slip systems of the underlying crystal lattice. A suitable flow rule for a fcc crystal lattice can be formulated as: with the slip direction m α and the normal of the slip plane n α defining the slip system α [12]. The shear rate of plastic deformationγ α in the slip system α caused by the Schmid shear stress τ α = S p : m α ⊗ n α is calculated using a phenomenological power law by Rice, Hutchinson and Peirce [13] γ α =γ α Therein a reference shear rateγ α 0 acts as a proportional factor on the ratio of the Schmid shear stress acting on the slip system α and the corresponding critical resolved shear stress τ α c . The evolution of the critical resolved shear stress is given as:τ where h αα is a hardening matrix accounting for active and latent hardening In the equation above h 0 , a and τ ∞ are hardening parameters and q αα assigns a factor of 1.0 for coplanar slip systems and 1.4 for systems with distinct slip planes. For the fcc crystal structure, the hardening parameters are assumed to be the same for all 12 slip systems [12]. The corresponding material parameters for Haynes 214 have been identified by Fischer et al. [6] by comparing experimental tensile tests at different strain rates with virtual tensile tests on a polycrystalline representative volume element (RVE). According to Jamaloei et al. [4], Ni 3 Si precipitates can be found in the ductile phase of the brazing layer. These precipitates show the same crystal structure as the γ phase of Haynes 214 and have a similar size as the γ precipitates. Therefore, the already identified material parameters of Haynes 214 should be suitable for the ductile phase in the brazing layer. The crystal elastic constants were chosen to match not only the effective Young's modulus E and Poisson's ratio ν but also to achieve a single crystal anisotropy factor typical for a Ni-based alloy (Zener ratio of A = 2.4) [14]. This is an enhancement to [6] in which the elastic anisotropy of the grains is neglected.
The remaining phases (phases 2 and 3 in Fig. 1) are assumed as purely brittle. As a first approach, a purely isotropic linear elastic material behavior that matches the effective elastic properties of polycrystalline Haynes 214 [15] is assigned to the brittle phases. Table 2 summarizes the material parameters used for the ductile and brittle phases at room temperature (T = 20 • C) as well as at T = 1100 • C which is close to the melting point of Haynes 214 (T melt = 1355 • C).

Macroscale
The most important effects on the macroscale are expected to be the temperature and rate dependencies of the elastic-plastic material behavior. A preliminary study considering fully resolved grains assigned with crystal elastic-plastic material behavior has been conducted to analyze the influence of varying crystal orientations on the base metal layers. It was shown that the stress heterogeneity due to the contrast between ductile and brittle layers is a magnitude of order larger than the stress heterogeneity due to the contrast between differently orientated grains in the base metal. Therefore, the crystalline structure of the base metal can be regarded as a second-order effect on the fluctuations of mechanical fields in the layered composite. Consequently, the base metal can be regarded as homogeneous and isotropic. A further study on the propagation of microstructural variation to the macroscopic material behavior showed that also the brazing layer can be regarded as homogeneous and isotropic (see Sect. 3.1). These considerations lead to the macroscopic material behavior of the base metal-braze layered composite being modeled as follows: • Homogeneous material in each layer, • Isotropic elastic-plastic material with von Mises yield criterion and isotropic hardening, and • Strain rate dependency of the yield surface.
All necessary material parameters for the base metal are taken from tensile tests conducted at different strain rates by Ulan kyzy [10]. For the brazing layer, a homogenization approach using a dual-phase RVE is employed (see Sect. 3.1). The uniaxial yield limit at different levels of plastic equivalent strain and strain rate is given as tabular data to the FE-solver Abaqus [16].

Representative volume element of the microstructure in a brazing layer
On the microscale the focus lies on capturing the relevant effects of the multiphase microstructure consisting of brittle and ductile phases. To model the effect of the anisotropic flow behavior of the grains in the ductile phase, the above-described crystal elastic-plastic constitutive law is employed. Compared to modeling isotropic phases (as done by, for example, [17,18]), additional contributions to the stress and strain fluctuations due to crystallographic misorientation can be depicted. This is potentially important for assessing the potential risk of micro-damage.
A fully periodic RVE is used to model the microstructure of the brazing layer in a highly resolved manner. A necessary condition in the concept of the RVE is the separation of scales [19], which can be formulated as In the microstructure of the brazing layer, the characteristic length of the heterogeneity L h is the diameter of the smallest silicides (3 in Fig. 1). The macroscopic length scale L macro is set to be the brazing layer thickness t braze , which lies in the range of 50 µm to 150 µm. For the generation of the geometry of the RVE, the software Neper is employed [20]. The implemented grain growth algorithm is used which is based on a Laguerre-Voronoï tessellation and lets the user introduce statistics for the grains' equivalent diameter (diameter d of a sphere with same volume as the grain) and shape. In order to account for differently sized grains, the grain diameter is set to lie in the range L h = [0.5 − −2.0] µm. The grains' equivalent diameter follows a logarithmic normal distribution (mean value of 0.87 µm and standard deviation of 0.3 µm). This logarithmic normal distribution was chosen as it is generally able to describe measured grain size distributions as shown by, for example, [21] and allows for a few larger precipitates that were also observed, even though those were of small quantity. Through a "coloring" process in which one out of two phases is assigned to each grain the RVE is divided into a brittle and a ductile phase. Thereby, the grains can be grouped to form a bigger microstructural cluster. The coloring process is random with the fraction of the ductile phase being the only condition that has to be satisfied. In combination with the varying grain diameter, a wide range of microstructural cluster sizes and shapes can be achieved. Compared to a multilevel tessellation [22], this approach has the advantage that the internal boundaries in a single phase are conserved. This is relevant for the ductile phase wherein varying crystallographic orientations depict a possible source of heterogeneity inside a single phase. In this approach, the grain size dependency of the flow stress (Hall-Petch relation) is not transferred to the grains individually as suggested by Meier et al. [23] or Fischer et al. [6]. Instead, the critical resolved shear stress τ c and the saturation stress τ ∞ are chosen to be the same for all grains in the ductile phase, and their values are set according to the mean equivalent grain diameter using the Hall-Petch constants k c and k ∞ from [7]. It is assumed that the geometrically necessary plastic deformation at the grain boundaries depicted by a spatial jump in the crystallographic orientation in combination with the hardening behavior following the equations (5) and (6) already accounts for a varying resistance against dislocation slip for grains with distinct sizes. Smaller grains have a higher proportion of areas of geometrically necessary deformation than bigger grains and therefore exhibit a higher mean flow stress. A shortcoming of this approach is that the contrast of the flow stress between two differently sized grains only builds up gradually with progressing deformation but is not present at the initial plastic deformation. As the mechanical contrast between the brittle and ductile phases is expected to be much higher than the contrast between grains of the ductile phase due to distinct grain sizes, the above-mentioned shortcoming is accepted in order to reduce the implementation effort. At this point, it should also be noted that the grain size hardening is only applicable to polycrystals up the recrystallization temperature (T /T melt < 0.3-0.5) [24]. For the elevated temperatures (T = 1100 • C) considered in this work, the grain size effect is assumed to be insignificant and is therefore neglected (see Table 2). The geometric tessellation consists of 2000 grains. This is a rather large number of grains compared to studies by other authors [9,17] who employed RVEs consisting of only 150 grains. However, regarding the complex microstructural clusters showing features of different size this is necessary to ensure the representativeness of the volume element. Figure 2 shows the Laguerre-Voronoï tessellation with varying grain sizes as well as the colored RVE showing the dual-phase microstructure with a fraction of ductile phases of P ductile = 0.50.
Due to its superior computational performance for fully periodic solution fields, a spectral method for solving mechanical equilibrium is employed on the microscale. The RVE is discretized by 128 voxels in each spatial direction. A convergence study that justifies this numerical resolution can be found in "Appendix B". The spectral method implemented in the simulation toolkit DAMASK [12] was verified against a FEM simulation. This was done for the case of the high contrast in the numerical tangential stiffness between the purely elastic brittle and the yielding ductile phases. This characteristic of the microstructure in the brazing layer is challenging for the spectral solver due to jumps in the stress field at phase boundaries. Both solution methods showed good agreement, and the spectral solver is considered to be suitable for the microstructures that are observed in the honeycomb brazing layer. The micromechanical model itself was validated with an experimental analysis of a similar microstructure taken from the literature (see Sect. 4).
In order to assess the significance of microstructural variation influencing the mechanical behavior on the macroscale, virtual uniaxial tensile tests were performed on different microstructures. Concerning the mechanical model on the macroscale (see Sect. 3.2), it is important to evaluate whether the effective material parameters show significant uncertainties leading to a potentially heterogeneous and anisotropic material behavior in the brazing layer on the macroscale. The tests were carried out at room temperature and under quasistatic conditions. As a measure for the uncertainties on the macroscale due to variations in the microstructure, the coefficient of variation CV is introduced as the ratio of the standard deviation to the mean value of the mechanical response of the RVE Two types of microstructural variation were analyzed. The first one was the variation due to different phase morphology. A random grain-to-phase-assignment process was applied to generate 10 different microstructures For the two types of microstructural variation, the CV of the mechanical response at 5% volume-averaged strain is shown in Fig. 3.
The stress measures evaluated are the stress component in the load direction σ xx as well as the 90% quantiles of the first principal stress in the brittle phases P 90 {σ I,brittle } and the hydrostatic stress in the ductile phase P 90 {σ m,ductile }. The CV is small for both types of microstructural variation. For a varying phase morphology, the mechanical response of the RVE varies stronger than for the varied crystallographic texture, but the CV is still less than 7%. The chosen fraction of the ductile phase of P ductile = 0.50 is expected to form the most complex grain clusters and subsequently the strongest possible morphological variation. Regarding the crystallographic variation, it is expected that an increasing fraction of the ductile phase could increase the microstructural variation further to a certain degree. However, since the observed effect of crystallographic variation on the stress measures of interest is significantly smaller than the effect of the morphological variation, any potential increase in the CV is expected to be negligible. Therefore, an assumption of isotropic and homogeneous material behavior in the brazing layer on the macroscale is reasonable as long as the volume fraction of the ductile phase is constant.

Macroscopic model of the layered composite
The macroscopic tensile sample is modeled as a layered composite of thickness t comp with a brazing layer that shows spatial thickness fluctuations. This is motivated by the observation that the depth of interdiffusion varies along the in-plane directions forming a non-flat interface between the brazing layer and the base sheet metal. The layers are generated by coloring a Voronoï tessellation with 2000 equally sized cells on the macroscale according to the cell centroid positions, which leads to an interface roughness in the order of magnitude of the Voronoï cell diameter. Figure 4 shows the procedure schematically. As a model input parameter, the mean thickness of the brazing layer t braze is used.
The interface roughness is smaller than the sample's in-plane dimensions by a factor of approximately 12. Therefore, the condition of the separation of scales is satisfied, and a semi-periodic macroscopic RVE with periodic boundary conditions applied in the x-and z-directions is employed to reduce the computational time significantly. For the y-direction force free surface boundary conditions are applied. The uniaxial tensile load  is introduced over a multipoint constraint that prescribes a volume-averaged deformation in the x-direction. Details of the implementation of periodic boundary conditions can be found in, for example, [25]. The cubeshaped computational domain is discretized by 65 3 brick finite elements with linear shape functions (C3D8R). An example of the layered composite model with a mean thickness of the brazing layer of t braze = 0.5 t comp is provided in Fig. 5 together with the sources of the material behavior for each layer that are described in Sect. 2. To be able to deal with the free surfaces in the y-direction of the tensile samples, the finite element method is chosen as the solution method for mechanical equilibrium on the macroscale. This modeling approach still accounts for an accurate representation of the interface between the base metal and the brazing layer while neglecting second-order effects caused by anisotropic material behavior of the oligocrystalline base sheet metals. A summary of the reference configuration in terms of geometry and loading for the simulations on the microscale as well as the macroscale is given in "Appendix C". The values reported there are used if not stated differently. Fig. 6 Comparison of simulation and experiment [8] in terms of engineering stress-strain curves for uniaxial compression. The segments on the strain axis indicate regimes of distinct degrees of accordance between simulation and experiment

Validation of the microscopic model
Li et al. [8] experimentally analyzed the mechanical properties of a nickel silicide similar to the silicides in the brazing layer of the honeycombs (phases 2 and 3 in Fig. 1). The brittle silicide matrix was enhanced with ductile particles. Those particles consisted of nickel solid solution with Ni 3 Si precipitates. The precipitates have a similar size and the same crystallographic structure [26] as the γ precipitates in Haynes 214. A very similar microstructure is also found in the brazing layer of honeycomb double-foiled segments. Therefore, the mechanical tests conducted by Li et al. [8] are suitable for validating the micromechanical model. In the experiment, the ductile particles were encased by Ni 3 Si rims. It is not clarified if these rims show a brittle or ductile behavior. Their hardness lies between the hardness of the brittle matrix and the ductile particles but is closer to the hardness of the brittle matrix. Therefore, in the simulation model the bounding cases of purely brittle and purely ductile rims are both considered, which leads to a fraction of the ductile phase of P ductile = 0.25 and P ductile = 0.5, respectively. The comparison of the volume-averaged stress-strain curves between simulation and experiment (Fig. 6) shows a good agreement regarding the bounds in the fraction of the ductile phase. Three segments on the volume-averaged strain axis can be distinguished. For elastic deformation and small plastic strains (I) the experimental results almost coincide with the simulation results assuming purely brittle rims (P ductile = 0.25). With increasing plastic deformation, the stress response is overestimated by the simulation, but the stress value as well as the hardening slope remains in the bounds given by the simulation (II). At approximately 2.5% volume-averaged strain, the hardening slope falls below the lower bound, and at approximately 5% volume-averaged strain also the stress response falls out of the given bounds (III). This deviation could partly be explained by the increased hardness of the Ni 3 Si rims. The onset of plastic deformation in the rims would give additional contributions to the volume-averaged deformation at higher stresses. However, this effect can explain the behavior not further than up to segment II. Regarding the drastic decrease of the hardening slope in segment III, there has to be an additional source of inelastic deformation which is possibly micro-damage. This hypothesis is in accordance with the observation of micro-cracks in the brittle matrix at rather small strainsε = 2.3% reported by Li et al. [8].

Sensitivity studies on the micro-and the macroscale
In order to obtain a better understanding of the constitutive behavior of the brazing layer, sensitivity studies are conducted on the microscale. The flow curves extracted from the micromechanical simulations will be incorporated in the material behavior used for the macroscopic studies.

Flow behavior in dependency of the volume fraction of the ductile phase
According to Hawk et al. [3], the fractions of brittle and ductile phases depend significantly on brazing parameters such as brazing temperature and duration. Therefore, a wide range of possible phase fractions has to be considered. Figure 7a shows the volume-averaged stress-strain curves for five different fractions of ductile and brittle phases. A strong influence of the phase fraction is apparent. For a decreasing fraction of the ductile phase, the hardening of the dual-phase microstructure increases following a nonlinear rule of mixture. As we do not yet have a quantitative criterion for damage initiation, the validity limit of the micromechanical model is drawn into the stress-strain curves in Fig. 7a. The limit is an estimate of the maximal stress and strain based on the experimental findings of Li et al. [8] who observed a fracture stress of approximately σ max = 2400 MPa and a fracture strain of approximately ε max = 0.025. Taking the first principal stress in the brittle phase as a failure criterion, one obtains an estimate of the maximum volume-averaged stresses and strains that are possible. In this estimate of the validity limit, it is assumed that the micro-damage starts in the brittle phases at the micro-stress level that leads to the failure in the microstructures used for validation (P ductile = 0.375). For an increasing fraction of the ductile phase, the failure criterion has to be extended in a way that it incorporates the failure of the ductile phase as well.
Looking at the yield limit in dependency of the fraction of the ductile phase (Fig. 7b), one can observe that for microstructures dominated by the brittle phases (P ductile < 0.5) the hardening slope increases progressively following a strongly nonlinear rule of mixture. However, the initial yield point is hardly dependent on the fraction of the ductile phase. This behavior was already predicted by Fischmeister and Karlsson [27] for dualphase microstructures that show negligible phase contrast in the purely elastic regime and can be verified when reducing the amount of plastic strain used for the definition of the yield limit R p successively (see Fig. 7b).   [28] for which phase-specific deformations are taken from the RVE simulations. Volume-averaged von Mises stress in the ductile phase containing a hard second phase plotted over plastic strain of whole aggregate. For better readability, a successive offset of 250 MPa in the stress axis is employed for each phase fraction starting from P ductile = 0.75 At this point, it should be remembered that there is some heterogeneity in the ductile phase due to the crystal elastic anisotropy of the nickel solid solution phase. Therefore, there is also a small contrast between the phases in terms of their elastic behavior. However, according to the simulation results this contrast is not sufficient for a noticeable dependency of the initial yield point on the fraction of the ductile phase.
Fischmeister and Karlsson [27] suggested a model for the estimation of the flow stress of a soft matrix including a hard second phase. Therefore, this model can be used to verify the inverse case to a brittle matrix containing ductile particles that is analyzed in Sect. 4. The estimation is based on the hardening of the material by means of geometrically necessary dislocations (GNDs) at grain boundaries in the matrix as well as at phase boundaries. In the presented case, the soft phase is the ductile nickel solid solution and the hard phase are the brittle silicides. The flow stress of the soft phase containing a hard second phase can then be derived as: In this equation, the friction stress σ r,ductile , the collecting constant C (for details, see [28]), the shear modulus G and the magnitude of the Burgers vector b are constants. As we do not alter the average grain size in the ductile phase when varying the fraction of the ductile phase, the mean free path for dislocation slip λ ductile is also constant. This means that the flow stress in the ductile phase is a function of only the volume-averaged strain of the whole aggregateε, the averaged strain of the brittle phasesε brittle and the ratio of the phase fractions P brittle P ductile . As Fischmeister considers GNDs inside the soft phase while assuming a purely elastic deformation in the hard phase his model is well applicable to the presented microstructure. Figure 8 compares the flow stress in the ductile phase extracted directly from the RVE simulations with the flow stress in the ductile phase calculated by Eq. (9). The difference in averaged strain between the aggregate and the brittle phases is taken from the RVE simulations. Both approaches show a similar trend of the flow stress of the ductile phase with increasing deformation as well as varying fractions of the ductile phase. However, the initial hardening of the ductile phase is underestimated by Eq. (9) for small strains. In contrast, for larger strains the semi-analytic approach overestimates the amount of hardening. A possible mechanism that most likely should be incorporated to explain this mismatch is the redistribution of stress from the ductile to the brittle phases. One can indeed verify this by plotting the ratio of the volumeaveraged von Mises stress in the brittle and ductile phases over the averaged plastic strain (see Fig. 9). For a decreasing fraction of the ductile phase, the accordance deteriorates. This is to be expected as the assumption of a continuous ductile matrix becomes increasingly violated which is also clear when inspecting Fig. 9. The effect of stress redistribution from the ductile to the brittle phases becomes significantly stronger for a decreasing fraction of the ductile phase. In the limit case, a brittle matrix impedes the plastic deformation of the ductile phase. An empirical modification of Eq. (9) to reduce the mismatch between the two approaches is given in "Appendix D".

Strain rate and temperature dependency of the flow behavior
As two important material properties for the analysis of rubbing events, the strain rate as well as the temperature dependency of the initial flow behavior are investigated for different fractions of the ductile phase. The positive relative difference in the yield limit for the extremes in the strain rate or the temperature is evaluated as: where R max p0.2 corresponds to the maximum value of the strain rate or the minimum value of the temperature and R min p0.2 corresponds to the minimum value of the strain rate or the maximum value of the temperature. As can be seen in Fig. 10, the temperature dependency of the initial yield limit is hardly influenced by the fraction of the ductile phase. The reason for this observation could be found in the temperature-dependent material behavior of the single phases. The ductile phase shows a temperature dependency in the elastic constants as well as in the yield limit in the form of thermal softening. Therefore, it might be expected that the temperature dependency of the dual-phase microstructure rises in its magnitude with an increasing fraction of the ductile phase. However, this is not the case as the brittle phases also show a temperature dependency in their elastic constants.

Assessment of the risk of micro-damage
According to experimental observations, micro-damage is assumed to occur first in the brittle phases [8]. However, the stress state in the ductile phase should also be monitored since the experimental observations also show small areas of the fracture surface that indicate ductile material failure. The risk of damage in the brittle phases can be assessed in terms of the first principal stress according to the Rankine criterion [29]. In the ductile phase, the hydrostatic stress can lead to the formation and growth of micro-pores [29] (an alternative approach would involve the plastic deformation as driving force for ductile damage). Figure 11 shows the spatial distribution of the first principal stress as well as the hydrostatic stress for uniaxial tension. The volume-averaged strain is chosen asε = 2.5%, which corresponds approximately to the strain level at which micro-damage occurs in the validation example (transition between segments II and III in Fig. 6). Stress peaks occur mainly in regions where narrow segments of the brittle phases are surrounded by a dominant ductile phase (see marks 1 in the left part of Fig. 11). However, these characteristics of the microstructure are only a necessary but not a sufficient criterion for stress peaks, as the orientation of the narrow segments with respect to the macroscopic loading direction is also of importance. If the flux of the load is not aligned with the narrow segments of the brittle phases meaning that the brittle segment and the embedding ductile phase are locally aligned in series, no stress peaks can be observed (see marks 2 in the left part of Fig. 11). Regarding the hydrostatic stress in the ductile phase, one can observe a similar behavior with the exception that highly stressed narrow segments of the ductile phase embedded in a dominant brittle phase are oriented in series with respect to the load. For the brittle as well as the ductile phases, the stress peaks can be found predominantly near the phase boundaries. At this point, it should be mentioned that for a quantitative assessment the third dimension needs to be incorporated when comparing stress peaks in the microstructure. However, the observations based on two-dimensional cross sections do give first insights in the microstructural characteristics that are predestined for damage initiation. When looking at the same level of deformation but now considering a uniaxial compression test, one can observe approximately an inverse of the distribution of stresses in both the brittle and ductile phases compared to the tensile loading case (see Fig. 12). For the brittle phases, this means that damage will most likely occur in narrow segments that are aligned in series with an embedding dominant ductile phase with respect to the external load. The stress levels in brittle phases are similar for tension and compression. In the ductile phase, the maximal hydrostatic stress is significantly reduced in compression compared to tension. This indicates that damage will most probably occur in the brittle phases in the case of a compressive load.
The stress exceeding a critical value only at a single location in the microstructure does not necessarily indicate failure on the macroscale. Therefore, a quantile-based assessment of the relevant stress measures is more reasonable than monitoring the maximum stress value with respect to the whole microstructure. Fillafer [17] suggests the 90% quantile for the assessment of damage in steels with a dual-phase microstructure and a strong phase contrast. As the analyzed microstructure shares these characteristics, the 90% quantile is also chosen for the evaluation of the relevant stress measures (see Fig. 13).
In tension and for both phases, the stress measure assumed as most relevant for damage initiation increases for decreasing fraction of the ductile phase. The first principal stress in the brittle phases shows critically In compression, the level of the stress measure assumed as most relevant for damage initiation show a different trend with a varying fraction of the ductile phase. While the 90% quantile of the hydrostatic stress stays in the compressive regime for all analyzed phase fractions, the 90% quantile of the first principal stress stays on a high level hardly influenced by the fraction of the ductile phase. Still the two stress measures approach each other which makes a shift from brittle to ductile damage likely to occur for very high fractions of the ductile phase also in this load case.

Thickness of the brazing layer
First, a sensitivity study on the macroscale with respect to the influence of the braze thickness on the mechanical behavior of the composite is conducted. Figure 14 shows the volume-averaged stress-strain curves (left) as well as the first principal stress averaged over 20 different paths through the layers of the composite (right). Due to similar elastic properties of the layers, the yield point of the composite is not influenced by the thickness of the brazing layer. A linear rule of mixture can be recognized from the average stress-strain curves. Except for superimposed stress peaks of second order at the interfaces between the layers the case of ideal parallel loading is recovered for t braze > 0.2 t comp . At this point, it should be noted that the small peaks observed are artifacts from the modeling approach that can be reduced by further increasing the semi-periodic volume element. For braze thicknesses greater than this threshold, the stress in the brazing layer decreases continuously as the layer becomes increasingly incomplete (see Fig. 14, right). This decrease continues until no brazing layer is present, in which case the stress is homogeneously distributed.

Fraction of the ductile phase in the brazing layer
To assess the influence of the fraction of the ductile phase in the brazing layer on the mechanical response of the composite, the effective material behavior for three different phase fractions is transferred from the microscale to the macroscale as described in Sect. 2.2. Figure 15 shows the volume-averaged stress-strain curves (left) as well as the first principal stress averaged over 20 different paths through the layers of the composite (right). In contrast to the study of the thickness of the brazing layer, a nonlinear rule of mixture is observed. Also, the influence on the stress increase in the brazing layer is significantly higher when reducing the fraction of the ductile phase than for an increasing thickness of the brazing layer.

Conclusions and outlook
The comparison between the simulation results and experimental findings in the literature shows that the simulation model is capable of predicting the flow behavior of the brittle-ductile microstructure at small plastic strains reasonably well. To reduce the mismatch between simulation and experiment for advancing deformation, the introduction of micro-damage in the brittle phases appears to be necessary. A suitable microdamage model must incorporate local stress peaks and should be able to depict stress redistribution in the microstructure during damage evolution. The established phase field approach for damage applied to the microscale might be able to serve these demands.
The behavior of the layered composite on the macroscale approximately follows an iso-strain behavior. This makes the calibration of the micromechanical model with uniaxial tests on a layered specimen easier. Tensile tests at different temperatures and strain rates should serve as a calibration and validation basis for further model enhancements.
Furthermore, based on the simulation results the fraction of the ductile phase could be identified as a central microstructural parameter in the brazing layer, whereas variations in the morphology of the phases and in the crystallographic orientation of the grains in the ductile phase appeared as second-order effects. Compared to the thickness of the brazing layer, small fractions of the ductile phase in the brazing layer can cause a significantly higher maximum stress in the layered composite. Therefore, an increase of the fraction of the ductile phase is considered as most important for a reduction of the loads during rubbing events in turbomachinery.
Using the results from the fully numerical RVE approach, the analytical model proposed by Fischmeister and Karlsson [27] can be modified regarding the employed root function to account for the effect of stress redistribution between the ductile and the brittle phases. However, when doing so, the physical basis for such a modification has yet to be found.

B Convergence of the spectral solver
The convergence of the spectral solver with respect to the number of voxels per edge of the RVE is assessed based on the 90% quantile of the stress measures that are relevant for damage analysis. For all considered stress measures, good convergence is achieved at rather small numbers of voxels (see Fig. 17). To ensure a good approximation of the local stress fields, 128 voxels per edge of the cubic RVE are chosen.

D Empirical modification of semi-analytical approach
The mismatch between the results of the fully numeric RVE simulations and the semi-analytic approach by Fischmeister motivates a modification of Eq. (9). Instead of the square root the fourth root is proposed to fit the RVE simulations over the whole range of plastic strain that is considered. The modified equation reads: Figure 18 shows a significantly improved accordance between the fully numerical and the semi-analytical approaches. However, this modification is only based on an empirical analysis of the RVE simulations. Fig. 18 Improved accordance between the fully numeric and semi-analytic approaches. Volume-averaged von Mises stress in the ductile phase containing a hard second phase plotted over plastic strain of whole aggregate. For a better readability, a successive offset of 250 MPa in the stress axis is employed for each phase fraction starting from P ductile = 0.75