The differences of damage initiation and accumulation of DP steels: a numerical and experimental analysis

Many studies have examined the damage behaviour of dual-phase steels already. It is a topic of high interest, since understanding the mechanisms of damage during forming processes enables the production of steels with improved properties and damage tolerance. However, the focus was rarely on the comparison between representatives of this steel class, and the numerical simulation for the quantification of damage states was not thoroughly used. Therefore, this study compares the damage initiation and accumulation of two dual-phase steels (DP800 and DP1000), which are used in the automotive industry. Additionally, parameter sets of a phenomenological damage mechanics model with coupled damage evolution are calibrated for each material. The combined analysis reveals an earlier initiation of damage for the DP800, where the damage accumulation phase is prolonged. For DP1000 the damage nucleates only shortly before material failure. The material model is able to correctly predict the behaviour, while experimental analysis confirms the prediction via light optical and SEM metallography.


Introduction
The usage of dual-phase (DP) steels has been on the rise in recent years. Especially the automotive industry shows high interest in the development of these advanced high strength steels (AHSS), since DP steels show high strength values while still maintaining good formability. Thus, a lightweight component design can be achieved by reducing component thickness while still keeping the identical safety conditions (Davies and Magee, 1979). These specific properties result from a distinct microstructure, that is composed of a soft ferritic phase with hard martensite islands on the grain boundaries and triple points of ferrite grains. Due to the difference in mechanical properties of the two phases, plastic behaviour of DP steels shows a relatively low yield to tensile ratio, pronounced strain hardening and excellent global formability. The reason for this extraordinary property profile lies in the partitioning of stress and strain between the involved phases, allowing for multiple degrees of freedom for microstructural design (Bieler et al., 2009) .
The strain partitioning between ferrite and martensite depends vastly on the specific microstructure. Marteau et al. reported that the local microstructural neighbourhood is the critical factor for strain heterogeneity (Marteau et al., 2013). Strain accumulates mostly in the ferrite forming localized bands with an angle of 45-50 • with respect to the loading direction (Ghadbeigi et al., 2010;Tasan et al., 2014a), whereas the martensite carries the majority of the applied stress (Tasan et al., 2014b). Therefore, martensite is elastically deformed for materials with low martensite content, while its deformation behaviour is plastic for high contents (Shen et al., 1986). The local microstructure especially is determining the strain distribution, e.g. average size of martensite islands and global distribution (Park et al., 2014;Saai et al., 2014).
Due to this inhomogeneity in the material constituents' behaviour, the microscopic damage modes of dual-phase steels differ quite significantly to those of common structural steels. Where for structural steels the inclusions play the major role for void initiation, in DP steels damage incidents occur in relation to the two phases, martensite and ferrite (Tasan et al., 2010). Mechanisms for damage initiation in dual-phase steels are mostly decohesion of the martensite/ferrite interface, cracking of the martensite phase, or a localization of plastic strain in the ferrite phase, which results in debonding of the ferrite grain boundaries (Ahmad et al., 2000). The mode for the damage initiation depends on the microstructure and the resulting strain heterogeneity (Kadkhodapour et al., 2011). Therefore, grain size and martensite content do play an important role (Maire et al., 2008;Ramazani et al., 2013;Tasan et al., 2015). Additionally, martensite morphology influences the early damage nucleation (Ghadbeigi et al., 2013;He et al., 1984). Besides, observations have shown, that for banded martensite cracking is far more likely than a decohesion of the interface boundary of ferrite and martensite (Avramovic-Cingara et al., 2009).
To assess the material's properties and predict the load bearing capabilities of structures, damage mechanics models are widely used for DP steels, e.g. in the automotive industry. In the field of damage mechanics, two different model types exist: Coupled and uncoupled models (Besson, 2010). For the coupled damage mechanics models usually a damage variable is employed to reduce the yield potential according to the softening resulting from ductile damage in the material during deformation. In case of the coupled models, micromechanical models are very popular, for instance the Gurson-Tvergaard-Needleman (GTN) model (Gurson, 1977;Tvergaard, 1981;Tvergaard and Needleman, 1984). Micromechanical models are characterized by the depiction of physical phenomena like void nucleation, growth and coalescence through sets of parameters. Therefore, the parameters are interdependent and thus, an extensive iteration process is necessary for the parameter calibration (West et al., 2012). Alternatively, phenomenological, coupled models are used to describe the damage in materials numerically. In contrast to the micromechanical models, damage evolution is treated in a macroscopic way, where a number of effects are described by a mathematical formulation. A good example for this type of model is the Lemaitre model (Lemaitre, 1985(Lemaitre, , 1992, which describes damage as an irreversible process. Contrary to that, uncoupled models describe the material behaviour including fracture without taking damage into account. Both the Johnson -Cook (Johnson and Cook, 1985), as well as the Bai-Wierzbicki (BW) model are good examples for this type of model (Bai and Wierzbicki, 2008). Further development has been applied by Lian et al., who combined the advantages of uncoupled and coupled models into a hybrid formulation, making it the modified Bai-Wierzbicki model (MBW) (Lian et al., 2013). The model therefore holds an easy formulation and combines it with the influence of damage onto material behaviour. The model has been developed further since its inception. For the first version a locus for the damage initiation point, which was dependent on both stress triaxiality and Lode angle was utilized. Additionally, a set of critical values for the damage variable was applied, at which material fracture was assumed in the numerical simulation. Wu et al. changed that considerably by implementing a locus for the fracture, as well as considering non proportional loading paths until the inception of ductile damage (Wu et al., 2017). A further development of the MBW model was made by Shen et al. to characterise the influence of loading orientation, which was used to describe the anisotropic ductile damage and fracture behaviour of pipeline steels . Since the MBW damage mechanics model is easy to use and calibrate, while also depicting the damage behaviour accurately, it is applied here for the characterisation of damage behaviour in DP steels.
While many studies focused on the damage in dual-phase steels from an experimental standpoint, it is hard to experimentally determine the evolution of damage during the tests. Therefore, this study aims to enhance the experimental investigation by performing finite element (FE) based numerical simulations that are utilized to quantify the damage in the material during forming processes. Thus, in this study, two dual-phase steels, DP800 and DP1000 were compared. Their damage and fracture properties are distinctly different, while the strength is not very far apart. To compare the materials behaviour, a damage mechanics model has been used that can describe both, damage initiation as well as ductile material fracture, while also taking the changes of the stress state during deformation, due to non-proportional loading effects, into account . This allows a comparison of the damage initiation for different stress states between the materials. Additionally, by means of a calibrated fracture locus, the damage accumulation phase can be analysed and compared. Thus, tensile tests were conducted on flat specimens of different geometries to gather information about materials deformations and damage properties under different stress states. On that basis, the material parameters of the modified Bai Wierzbicki model were calibrated. For the validation of the numerical results regarding the damage initiation and accumulation of the investigated material, interrupted tensile tests were conducted and a metallography analysis was performed by using the light optical microscopy and scanning electron microscopy (SEM).

Materials Characterization
In the present study, two dual-phase steels were evaluated for comparison purpose. Even though both materials are dualphase steels, vastly different properties are observable. These varying characteristics are obtained by distinct alloying concepts as well as heat treatment processes.   Fig. 1 shows a comparison of the respective microstructures at a magnification of 1000. It is very well observable that the average grain size of DP1000 is significantly smaller than that of DP800. Additionally, DP1000 has increased martensite contents of approximately 38% while DP800 contains about 32%. For DP800 a pronounced banding of the martensite in the microstructure is noticeable. Since the martensite bands run parallel to the rolling direction, there will be a significant influence on the mechanical properties. A certain extent of failure anisotropy is expected due to the banded microstructure, however, the anisotropic fracture properties are beyond the scope of this study and all tensile specimens were manufactured perpendicular to the rolling direction of both DP steels. Both steels were delivered with a thickness of 1.5mm; their respective chemical compositions are given in Table 1. While the alloying concepts show noticeable similarities, some minor differences are present.
The carbon content for DP1000 is decreased compared to DP800, thus leading to higher carbon concentration in the martensite phase for the DP800, since the phase fraction of martensite is higher for DP1000. Thus, it is to be expected that the strength of the martensite is reduced for DP1000 due to the decreased carbon content, therefore leading to a bigger contrast of properties between ferrite and martensite in the DP800. Furthermore, manganese and chromium contents are different, which leads to slight disparities due to solid solution hardening. Additionally, Si as well as Mn and Cr reduce the critical cooling rate needed for forming martensite, thus influencing the respective time -temperature -transformation graphs. On top of that the solubility of carbon in ferrite is reduced by silicon. Therefore, both materials will have very distinct processing routes tailored to the respective production process. For the characterization of the resulting mechanical properties, isothermal, uniaxial quasistatic tensile tests were carried out on flat specimens without a notch. To ensure a proper depiction of the material's properties, three tests were carried out per DP steel. A video extensometer was used to capture the elongation of the specimen during deformation, where the starting length of 40mm was used. The necking took place inside the area tracked by the extensometer for all 6 tensile tests, ensuring a good comparability.
The results for both materials are shown in Fig. 2. The scatter shown is the difference in fracture, resulting from the three tensile tests per material mentioned before. From this figure it is obvious, that DP1000 shows higher strength, while the strain at fracture of DP800 is considerably higher. The higher strength is a result of both; the higher martensite phase fraction as well as significantly refined grains in DP1000. In addition, there is a clear scatter in the elongation values, with elongation at fracture between 18 and 23 percent for DP800, while the scatter for the DP1000 is about one percent. This variation in elongation at fracture can be explained by the clearly pronounced band structures, which can lead to significant deviations depending on the position of the bands in the specimen. Since DP1000 fractures shortly after uniform strain, the necking is far less pronounced than in DP800. During necking the stress state in the sample can change quite significantly leading to non-proportional loading paths during the deformation of the sample. Thus, to describe the materials behaviour after the uniform strain, it is necessary for the material model to consider the effects of the changes of stress state during deformation. Therefore, Fig. 2 Engineering stress-strain curve of uniaxial tensile tests of DP800 and DP1000 a development of the existing MBW model was required to determine the material behaviour and damage accumulation more realistically.

Extension of MBW model for non-proportional loading paths (npMBW-19)
In the framework of continuum damage mechanics, the modified Bai Wierzbicki (MBW) model has been proposed by Lian et al. (2013) and widely applied to describe the damage and fracture behaviours of various grades of steels (Lian et al., 2013;Münstermann et al., 2017;Wu et al., 2017;Novokshanov et al., 2015;Shen et al., 2020;Liu et al., 2020). Like in other damage mechanics models, the significant influences of stress state on the ductile fracture are considered through defining a strain based criterion which is usually a weighted function of two particularly important variables, the stress triaxiality η and the Lode-angle parameter θ that are related to the three stress invariants. (1) where σ 1 , σ 2 and σ 3 are principal stresses and σ is the von Mises equivalent stress. For the material model, the Lodeangle parameter θ was used, which has a linear relationship with the Lode-angle θ .
The effects of stress state on plasticity in some metallic materials have been reported, while steels typically show a negligible pressure sensitivity, therefore, only the effects of Lode-angle parameter are considered in the yield criterion of the MBW model.
Where D is a scalar variable to quantify the damage effects, yield stress σ y is determined by the equivalent plastic strain ε p and Lode-angle parameter θ , and σ y (ε p ) corresponds to the flow stress at given equivalent plastic strain ε p under the reference stress state. c s θ , c t θ , c c θ are the normalised strength under shear, tension and compression state and m is a material parameter with positive integral values that describes the Lode-angle sensitivity. γ is another stress state parameter with unique correlation to the Lode-angle parameter θ . According to the derivation of Lian et al., (Lian et al., 2013), the yield locus of MBW model is convex if the material parameters are located within the specific range of In the coupled damage mechanics model, two individual criteria have been defined to identify the ductile damage initiation (DDI) and ductile fracture (DF), which corresponds to the initiation of degradation on microscopic scale in the material and the loss of load carrying capacity on the macroscopic scale. Numerically, damage initiation, in this study, is defined as the onset of macroscopic softening due to damage, which must be taken into account by the numerical representation of the material behaviour. In order to consider the change of stress state during plastic deformation, the average values of the stress triaxiality η avg and the Lode-angle parameter θ avg have been used to describe the stress state for non-proportional loading paths (Wu et al., 2017;Mu et al., 2020).
Since the damage is dependent on stress state it is necessary to define equations for the initiation of damage, as well as the fracture, that represent this dependency. Therefore, the damage initiation locus (DIL) and ductile fracture locus (DFL) have been defined as two individual equations f di and f d f with the stress triaxiality and the Lode-angle parameter as independent variables. The instantaneous and average values of the independent stress state variables have been used in the damage and fracture criteria under proportional and non-proportional loading conditions, respectively. Under non-proportional loading conditions, these two equations describe the critical equivalent plastic strains at the moment of damage initiation and ductile fracture, respectively. (15) where D 1 -D 6 and F 1 -F 6 are material parameters used to define the damage initiation locus and ductile fracture locus. Under the condition that D 1 = D 5 , D 2 = D 6 and F 1 = F 5 , F 2 = F 6 , the DIL and DFL are symmetric with respect to the Lode-angle parameter and four independent parameters are enough to define the corresponding loci. Based on previous experimental observations, a cut-off value of the stress triaxiality η c exists, below which the initiation and evolution of ductile damage cannot be triggered due to pressure effects. η c = −1 3 as a reasonable estimation has been adopted in the MBW model (Wu et al., 2017). Therefore, when the stress triaxiality is lower than η c , the equations f di and f d f are set to be infinite. The damage initiation specified by this model is unrelated to the materials mechanisms of damage initiation, e.g. micro crack formation, void formation. Instead, it aims to describe the aggregative accumulation of the defects and their influence on the load bearing capabilities. For this step a plasticity model is no longer able to describe the materials mechanical behaviour (Keim et al., 2019). For the non-proportional loading, two indicators have been applied to describe the ductile damage initiation I dd and ductile fracture I d f respectively to consider the effects of stress state evolution.
The values of equivalent plastic strain and equivalent stress at the moment of damage initiation (I dd = 1) are defined as two characteristic variables ε p,c di and σ c di , respectively: After the damage initiation criterion is fulfilled, damage evolution is controlled according to the energy dissipation theory. Depending on the shape of damage initiation locus and ductile fracture locus, when the indicator of the ductile fracture I d f reaches unity, the damage variable D does not necessarily reach unity. Therefore, a critical value of the damage variable D cr exists, at which the material point will fail regardless of the value of the D variable: Where G f is a material parameter which controls the damage evolution rate. Linear damage evolution is assumed in the MBW model, which is expressed as: In summary, the damage evolution is determined by the two independent damage initiation and ductile fracture criteria. After a certain damage nucleation period, which is controlled by plastic deformation, damage evolution takes place. Once the indicator of the ductile fracture I d f reaches unity, the final crack propagation is triggered and failure occurs. Therefore, the model, hereafter called npMBW-19, is capable of representing the influence of the necking, and thus the change of stress state, during deformation.

Calibration of the new model for materials DP800 and DP1000
The calibration approach for the material models for both steels follows roughly the approach of Lian et al. (Lian et al., 2013(Lian et al., , 2014. Since the calibrated npMBW-19 model needs to be able to account for various stress states, the calibration of the material model is carried out on a variety of sample geometries. By varying the sample geometries in tensile tests, different stress states can be accomplished. In this study three differently notched specimen types were applied for the calibration of the material model in addition to the uniaxial tensile test. Used specimen types were: Notched dogbone samples (varying notch geometries at the edge of the sample), central hole samples (round, as well as elliptical holes in the center of the specimen) and plane strain samples (notch with different radii over the thickness of the sample). The applied specimens for each material can be seen in Fig.  3 and Fig. 4 The type of notch of the sample is abbreviated with an r continuing with the radius, for the notched dog bone samples.
The corresponding stress states, characterised by the Lodeangle parameter and the stress triaxiality in the applied samples are delineated in Table 2. To achieve multiple stress states, notches were modified with various radii to gain geometries of different stress states within one sample type. Per specimen geometry, three tensile tests were performed in For the determination of the base flow curve, the uniaxial tensile tests (T), presented in the previous chapter, were utilised. From the determined engineering stress-strain curve, the true stress-true strain curve was calculated until the uniform elongation point. This data was then used to fit the Hollomon -Voce hardening model to the material's flow curve via the Matlab curve fitting tool.
This specific hardening model was chosen, since it shows a good compromise between accurate representation at low plastic strains and realistic hardening behaviour for higher strains. In Table 3 the parameters for the Hollomon-Voce models are given for both, DP800 and DP1000.
After the calibration of the flow curve the basic parameters of the MBW model were determined (c s θ , c t θ , c c θ , m). This was done by iterating over multiple simulations using a range of different sample geometries.
Subsequently the damage and fracture parameters of the npMBW-19 model were determined. Damage and fracture criteria in this material model are described by equations 14 and 15. Therefore, the specified locus needs to be calibrated for both events, damage initiation and fracture (Lian et al., 2013;Wu et al., 2017). For the damage initiation locus, a comparison of force and displacement curve between simulation and experimental results was used. Since the damage described in this model is related to the accumulated damage incidents, a threshold method has been utilised to find the numerical damage initiation. For that reason, the numerical onset of damage was determined as the point where the deviation between simulated and experimental force and displacement curves was apparent. Similar methods have been used by other authors within the damage mechanics field (Børvik et al., 2001;Bouchard et al., 2011). At this step, the Lode-angle parameter and stress triaxiality as well as the equivalent plastic strain (PEEQ) are taken from the simulation. Since these parameters may vary locally, the element is chosen that shows the most critical state of stress and thus is most likely to encounter damage first. By extracting the Lode-angle parameter, stress triaxiality and equivalent plastic strain for a multitude of different tensile geometries, data points are gathered in the space defined by these three variables. Applying the curve fitting tool of Matlab, a function can be defined that describes the desired surface while using the obtained results as supporting points. For non-proportional loading paths it is necessary to average the stress state of the critical element, where damage happens first, over the simulated steps (Wu et al., 2017).
After determining the locus for the onset of damage (DIL), the effect that damage has on the component needs to be adjusted. In the npMBW-19 model parameter G f is calibrated to adjust the speed at which damage accumulates in the simulated material. G f is defined as the energy dissipation between damage initiation and complete fracture. When the softening is specified the fracture locus can be determined. The approach used for this determination follows the one from the damage initiation locus closely. This time the point for the experimental fracture is compared to the simulation. The step where the fracture should occur is identified and Lode-angle parameter, stress triaxiality and equivalent plastic strain are extracted for the critical elements. Again, the stress states are averaged from the point of damage initiation to the presumed fracture of the sample. After gathering the data for all sample geometries the locus is fitted in regards to the obtained points using the Matlab curve fitting tool. In Fig. 5 the final results of this calibration process are depicted for steel DP800. From this figure it is obvious, that a good match between experimental data and simulations was been obtained. The scatter for the experimental testing can be seen in the shaded areas. Likewise, the material model for DP1000 was calibrated. The applied flow curve can be seen in Table 3. Additionally, a damage initiation, as well as a ductile fracture locus were calibrated using the same approach as described above for the DP800. By duplicating the approach stated above, a good agreement with the experimental data could be reached (Fig. 6). Contrary to the DP800 almost no scatter could be found during the tests of the DP1000 material, which also shows no significant banding in its microstructure. The applied set of parameters can be found in Table  4. Interestingly, the calibrated G f parameter for DP1000 is higher which results in a slower development of the damage variable. This results in a fairly slow accumulation of damage after the initiation.

Damage behaviour prediciton
For the scope of this study, it is important to differentiate between failure and damage of a material or component. Be- Fig. 6 Comparison between experimentally obtained data(background and lighter color) and simulation results for DP1000 cause damage is the deterioration of materials properties before failure, especially the load bearing capacity (Lemaitre, 1992), it is not to be equated with component failure. Damage occurs on a microscale and is usually described as the development of voids inside the microstructure, while on a macroscale damage usually equates to cracks in the component and therefore can be seen as component failure. It is therefore highly relevant to differentiate between micro and macroscale (Tekkaya et al., 2017). For numerical analysis, damage is defined as the macroscopic reduction of the stress during loading, that cannot be described by basic plasticity modelling. Thus, Lemaitre introduced a factor for damage in a microstructure, which results in a reduction of the flow potential by the term (1 − D), where D is the damage variable (Lemaitre, 1985). The damage variable adopted by Lian et al. shows some differences to the one postulated by Lemaitre. While Lemaitre's damage variable is calculated based on the area fraction of defects, Lian et al. refer to the stress at damage initiation, divided by the energy required to create new surfaces in a volume of the material (see G f ), an adaptation of the damage evolution law used by ABAQUS finite element code (Lian et al., 2013). Accordingly, both damage variables are scalar, but there are quite pronounced differences between both numerical damage rules. These differences between the damage models must be distinguished, as well as the differences between micro-and macroscopic damage phenomena.
Because the damage law used in the MBW model does not refer to a physical material characteristic, like the area fraction of voids, except for the energy for cracks, it is a rather more phenomenological approach to model the influence of damage on the materials flow potential. The damage evolution, as discussed before, starts when a specific equivalent plastic strain (PEEQ) locally exceeds a certain threshold, which changes with stress states. The respective value for PEEQ is determined by the damage initiation locus. Af- Fig. 7 Comparison of the ductile damage initiation locus and the ductile fracture locus for DP800 (left) and DP1000 (right) ter this point, softening occurs in the simulation, which leads to a direct reduction of the resulting stress compared to virgin materials. The length of this following phase where damage accumulates depends on the stress state which is considered in the ductile fracture locus. The comparisons of the ductile damage initiation locus and the ductile fracture locus for each respective material are shown in Fig. 7. The shape of the loci for DP800 and DP1000 are different, as was to be expected. The distance between the plots is higher for the DP 800 material which leads to a longer damage accumulation phase. Merely for higher triaxialities and Lodeangle parameters around zero, the differences between the loci of DP800 and DP1000 is minimal. Some research suggests a different shape for the ductile damage initiation locus and the ductile fracture loci, especially for the area around a stress triaxiality of 0, namely shear stress state (Papasidero et al., 2015). Nevertheless, based on the experimental and numerical results in this study, the loci in Fig. 7 constructed for both steels using the corresponding calibrated damage and fracture parameters is validated within the range of investigated stress states. In the case of an application of the calibrated material model for even lower or higher stress triaxialities, the loci would have to be revisited to confirm or adapt their shape.
Due to the differences in the distance of the DIL and the DFL the damage accumulation phase is significantly different between the two steels. Fig. 8 shows this difference utilising the flow curves obtained from the uniaxial tensile test of both materials, as well as calculating the points for damage initiation and fracture under uniaxial tension condition (η = 1 3 , θ = 1) based on calibrated material parameters. The point for the damage initiation takes place at roughly the same strain for both materials, while fracture is delayed significantly for DP800.
To further examine precision of the numerical results, interrupted tensile tests were conducted for both materials. For each material a sample was therefore first tested until failure and subsequent specimens of identical geometry were stopped after a distinct strain was reached. The lowest elongation used in this investigation was the uniform elongation, as no or little damage is expected below this. This way a metallographic damage analysis could be carried out to investigate the average amount of damage that could be observed in a sample. For both materials unnotched dog bone specimens were utilised to ensure a good comparability. For the analysis of the damage in the material, light optical microscopy was chosen, since a bigger area can be investigated by light optical analysis, where scanning electro microscopy (SEM) pictures resolve only smaller areas of the samples. However, it is not easily possible to differentiate between voids and inclusions in the material. Thus the area fraction that is detected is not quantitatively representative of the actual void fraction. To find out about the area fraction for each sample, multiple pictures were taken to gather information about the scatter band where the actual values lie. For this analysis, the light optical pictures were converted to greyscale images, which were subsequently evaluated by a threshold method, with which a differentiation between matrix material and voids/inclusion could be made. For these steps Fiji was used as image analysis software (Rueden et al., 2017;Schindelin et al., 2012).
The results of this analysis are depicted in Fig. 9. To better compare the values for both materials, a normalisa- Fig. 9 Comparison of light optically detected area fraction for DP800 and DP1000 for different strains tion was carried out, where the current strain was divided by the respective fracture strain. A comparison of the values for the detected area fraction reveals a gradual increase for DP800, while for DP1000 no significant rise in fraction can be observed until just before fracture of the sample. The large scatter, especially at the beginning, can be explained by the lack of necking, which means that the region of interest cannot be identified accurately.
Thus, the damage accumulation phase for DP800 starts at lower strains relative to the fracture strain of the material. By contrast, the damage accumulation phase for DP1000 starts very late and just before fracture. Therefore, the damage in the material behaves exactly as predicted using the npMBW-19 model. To assess the damage state in the microstructure, pictures were taken in the SEM. Especially for DP1000 an analysis for higher magnifications was necessary to reveal if damage forms earlier than shortly before fracture. For 70% of the fracture strain, only very few events of damage initiation could be found under high magnification (Fig. 10). While the amount of these initiation locations increases with the strain, growth is very limited (Fig. 11).
A comparison of the samples that are at 95% of fracture strain reveals, that the voids in DP1000 (Fig. 11) are significantly smaller than in DP800 (Fig. 12). It is therefore concluded, that the damage accumulation phase for DP1000 is indeed significantly shortened. In particular, it is noticeable, that no void in DP1000 exceeds a length of 1µm , while the Fig. 10 Evolution of damage in DP1000. Single martensite crack with magnification of 5000 for 0.7 of fracture strain DP800 features multiple larger voids. Additionally, voids for DP800 are more circular, while they are shaped like cracks for DP1000 again leading to the conclusion, that there has been no time for growth after initiation. This is in line with the results demonstrated in Fig. 8, where a shorter damage accumulation phase is present in DP1000 and thus a lower decrease of load bearing capabilities is to be expected.
Besides, the damage initiation modes were investigated. For both materials, the prevalent modus for damage initia-  tion was the cracking of martensite islands. For DP800 the martensite bands especially were sites for damage initiation. Furthermore, decohesion of ferrite and martensite islands was found in the DP800 after about 80% of fracture strain.

Conclusions
This study showed significant differences between two industrially produced dual-phase steels. Starting with the experimental results, the difference in fracture strain was found to be significant with large scatter for the DP800. This scatter was attributed to the pronounced banding found within the material. For the numerical analysis it was found, that the change in the stress state during necking needs to be considered by the material model for proper simulation results. Therefore, the MBW model was extended to account for non-proportional loading paths. The parameters of the model were then fitted for both materials to reveal the disparities in the material behaviour numerically. Subsequently, the damage initiation and fracture loci were calibrated. The comparison of simulation results to the experimentally obtained force-displacement curves reveals a high agreement for both materials. Especially the differences in damage behaviour were modelled precisely.
The found differences during the experimental testing and analysis can be attributed to the differences in the microstructure. Especially grain size and martensite content, but also the pronounced banding in the DP800 play an important role for mechanical properties, as well as damage initiation and accumulation. In this study, it was shown, that the damage in both dual-phase steels initiates at similar equivalent plastic strains. Oppositely, the fracture happens at vastly different equivalent plastic, as well as global strains. This leads to completely different damage accumulation phases in the material. The numerical simulations showed an exceedingly different length of the damage accumulation phase for the two steels. This difference was subsequently verified by experimental tests, where the amount of damage in the material after an interrupted tensile test was examined. For these tests it could be shown, that DP800 exhibits a pronounced damage accumulation phase, while DP1000 fractures shortly after a critical amount of voids forms in the material. Thus fracture occurs with almost no damage accumulation and the void growth phase is nearly skipped. Therefore, the different microstructures lead to specific damage characteristics, which in turn influence and change the specific properties of the material. Additionally, the contrast in the mechanical properties between the two phases for both materials reinforce this effect. Since the carbon content in martensite is relatively higher in DP800 than in DP1000, the martensite fractures earlier leading to a relatively early damage initiation and longer damage accumulation phase. For the DP1000 the the contrasts are not so distinct, which leads to a comparatively late initiation of damage and failure shortly afterwards.
During the comparison of experimental and numerical results, it was quite obvious, that the presented material model is able to accurately represent the experimental tensile test curves, both uniaxial and notched specimens. The stress states do, however show only minor variance. Thus, for higher deviations, an adjustment of the averaging scheme for determining stress triaxiality and Lode angle parameter for the ductile damage and fracture loci might be appropriate. Furthermore, the model shows a discontinuity around the value of η c = −1 3 , since a fixed value at which no damage is applied in the model, will be difficult to deal with, when getting close to it (e.g. η = −0.32). However, since the loading paths in this study are exclusively above this value, this is a promising and important concern for future development.
The comparison of numerical and experimental ductile damage showed, that the presented material model is able to accurately predict the damage initiation, damage accumulation and fracture of both materials. Nevertheless, the accuracy of the damage initiation point in the material model is still an important topic for further investigation. Since the location of the fracture locus strongly depends on the damage initiation locus, a high precision for the DIL is desirable. However, the commonly used method of direct current potential drop (DCPD) is not feasible for DP steel, since its void volume is comparably low. Therefore, an improvement of the method to determine aforementioned damage initiation locus is necessary and currently examined. One possibility is to measure the density of the material to determine the time of damage initiation Schowtjak et al., 2019;Meya et al., 2019).
The analysis of damage initiation point and damage accumulation by light optical microscopy is rather qualitative than quantitative. Since statistical representativeness and accuracy have to be balanced for this type of examination, micro voids are not detected in the pictures. Thus, the values received are not the void area fraction. On top of that, light optical pictures show inclusions in a similar colour to voids, isolation of voids for analysis purpose is rather difficult for light optical microscopy. For a more quantitative result of void area fraction SEM pictures with high resolution over a big area seem to be more promising.
The presented results suggest that the damage, and therefore the materials mechanical properties depend on the microstructure of the respective steel. It is therefore of high interest to investigate the influence of each microstructural parameter on the damage characteristics as well as the mechanical properties of the material. Three-dimensional representative volume elements seem to be a promising approach to investigate the influence of different microstructural characteristics, like martensite volume content, martensite morphology and grain size.

Acknowledgement
The funding of this research by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) -Projektnummer 278868966 -TRR 188, is gratefully acknowledged.