Damage Evolution in Nakajima Tests of DP800 Dual P
Damage Evolution in Nakajima Tests of DP800 Dual P
Damage Evolution in Nakajima Tests of DP800 Dual P
Abstract. In order to extend the understanding of damage evolution in sheet metal forming,
standardized Nakajima tests are carried out on a DP800 dual phase steel. Sample geometries for
characteristic stress-strain states are drawn in incrementing stages and their damage states are
analyzed using light and scanning electron microscopy as well as micro hardness measurements.
Numerically analyzed load paths are correlated with the respective damage states to allow pre-
diction of damage evolution in deep drawing processes. The influence of anisotropy is investi-
gated by testing samples cut at various angles to the rolling direction of the sheet material. The
result of the conducted research is the understanding of interactions between load paths and
damage evolution in sheet forming. These results will later be used to optimize load paths in a
deep drawing process, taking into account Lode parameter and stress triaxiality to produce dam-
age controlled parts.
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution
of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Published under licence by IOP Publishing Ltd 1
International Deep Drawing Research Group 37th Annual Conference IOP Publishing
IOP Conf. Series: Materials Science and Engineering 418 (2018) 012048 doi:10.1088/1757-899X/418/1/012048
1234567890‘’“”
tensor [7]. The value of these invariants can be identified with certain stress states. The von Mises
equivalent describes the magnitude of the stress state. The Lode parameter 𝐿 is a normalized value and
takes values of 𝐿 = −1 … 1. A Lode parameter of 𝐿 = −1 describes biaxial tension with uniaxial com-
pression, 𝐿 = 0 describes pure shear stress and 𝐿 = 1 describes uniaxial tension. Stress triaxiality 𝜂
takes values of 𝜂 = −∞ … ∞ with 𝜂 = −2/3 describing plane compressive stress, 𝜂 = −1/3 uniaxial
compressive stress, 𝜂 = 0 pure shear stress, 𝜂 = 1/3 uniaxial tensile stress and 𝜂 = 2/3 plane tensile
stress [8]. Yin has identified positive values of stress triaxiality 𝜂 as an indication of damage critical
stress states in in-plane torsion tests of DP600 dual phase steel [9].
Figure 1. Nakajima test principle: (a) specimen geometry; (b) test setup.
Evolving stress states are described using load paths. These show the evolution of loads, in this case
of the three invariants, over time [10]. To achieve different load paths, Nakajima tests are executed.
These were originally developed to record forming limit diagrams, which help in predicting failure in
sheet metal forming. Nakajima tests use specimens with varying cutout widths (Figure 1 (a)) to achieve
different stress-strain states. These specimens are clamped using a circular blank holder and then drawn
until fracture through a drawing die using a hemispherical punch (Figure 1 (b)). Major and minor strains
on the specimen surface are recorded when fracture occurs and plotted in the forming limit diagram,
resulting in a curve representing the limit of formability of the sheet material [11].
Pathak [12] investigated damage evolution in DP780 dual phase steel during Edge Stretching. Holes
were alternatively reamed and sheared into a tensile test specimen and the specimen was then drawn.
Experiments were interrupted at different strains and the damage states were analyzed. Void nucleation,
void growth and void coalescence were investigated as damage mechanisms. Void density, as an indi-
cator for void nucleation, showed a strong increase with increasing strain. Both average void diameter
and void orientation did not change significantly with increased strain [12].
2
International Deep Drawing Research Group 37th Annual Conference IOP Publishing
IOP Conf. Series: Materials Science and Engineering 418 (2018) 012048 doi:10.1088/1757-899X/418/1/012048
1234567890‘’“”
C3D8R). Convergence analysis has shown element lengths of 𝑙p = 0.3 mm in the sheet plane and 𝑙s =
0.5 mm in sheet thickness direction to give consistent results. The simulations have yielded load paths
for various cutout widths [12]. The simulations were validated with strain fields obtained from physical
Nakajima tests. While a significant drop-off in strain towards the specimen center could be observed in
the simulations that did not occur in physical experiments, strain results starting at 20 mm from the
center were consistent with the fields observed in experiment. Load paths were therefore extracted for
elements at the center line of the specimen, at the end of the cutout (see Figure 2 (a)). The extracted load
paths are shown in Figure 3. Specimens with cutout width 𝑤C = 47 mm were selected to use in exper-
iments since their load paths show gradually increasing values of the Lode parameters 𝐿 at the top sur-
face and a steep increase to the maximum value at the bottom surface.
Figure 2. Finite Element model of a 𝑤C = 95 mm Nakajima test: (a) location of investigated ele-
ments at 𝑡 = 10 s; (b) load paths at the upper specimen surface.
Figure 3. Load paths with respect to Lode parameter (a) at top surface, (b) at bottom surface for vari-
ous cutout widths.
3
International Deep Drawing Research Group 37th Annual Conference IOP Publishing
IOP Conf. Series: Materials Science and Engineering 418 (2018) 012048 doi:10.1088/1757-899X/418/1/012048
1234567890‘’“”
Figure 4. BUP 1000 Sheet Metal Testing Machine: (a) machine setup; (b) lubrication detail.
4
International Deep Drawing Research Group 37th Annual Conference IOP Publishing
IOP Conf. Series: Materials Science and Engineering 418 (2018) 012048 doi:10.1088/1757-899X/418/1/012048
1234567890‘’“”
bottom surface, average void percentage increases slightly from 𝑣 = 0.058 % at 𝑠d = 7.5 mm to
𝑣 = 0.060 % at 𝑠d = 12.5 mm and then decreases to 𝑣 = 0.053 % when the drawing depth is in-
creased to 𝑠d = 17.5 mm. Again, regarding standard deviation, no effect of larger drawing depths can
be seen.
Figure 5. Void percentage in specimens: (a) void identification in micrograph; (b) void percentage
averages for different drawing depths.
Possible reasons for the lack of effect on measured damage are issues with the evaluation procedure,
a lack of resolution with the light microscope, or an actual lack of detectable effects. Damage initiation
in the form of void generation happens through brittle fracture of the martensite phase, through ductile
deformation of the ferrite phase or through stress concentration at inclusions or precipitations [13].
These effects occur at a scale that cannot be investigated using light microscopy. Continued deformation
of the specimen does eventually lead to voids that can be detected with a light microscope, but these
deformations were not achieved in the experiments.
Due to symmetry, the direction of the maximum principal stress at the investigated location is ex-
pected to be perpendicular to the plane in which the micrograph was taken. Since voids will grow pri-
marily in the direction of the maximum principal stress, detection of growth is difficult.
5
International Deep Drawing Research Group 37th Annual Conference IOP Publishing
IOP Conf. Series: Materials Science and Engineering 418 (2018) 012048 doi:10.1088/1757-899X/418/1/012048
1234567890‘’“”
The results of the presented work and further investigations will be used in subprojects within the
Collaborative Research Center CRC/Transregio 188 Damage Controlled Forming Processes that are
specifically concerned with damage modelling to extend and calibrate various damage models for the
specific dual-phase steel used within the CRC.
Acknowledgements
The project presented in this paper is funded by the German Research Foundation (DFG) in the Collab-
orative Research Center CRC/Transregio 188 Damage Controlled Forming Processes, subproject Influ-
encing Damage in Deep Drawing, and in the Individual Grant Program Investigation of the Influence of
Reversing Strokes on Work Piece Properties in Deep Drawing with Servo-Press Technology. The au-
thors would like to also thank Prof. Dr.-Ing. A. Erman Tekkaya and Alexander Schowtjak, M.Sc. of the
Institute of Forming Technology and Lightweight Construction at TU Dortmund for support in the ex-
ecution of the experimental work.
References
[1] J. Rogelj et al.: Paris Agreement climate proposals need a boost to keep warming well below 2 C.
Nature 534.7609 (2016), pp. 631-639.
[2] H.J. Kim, C. McMillan, G.A. Keolian and S.J. Skerlos: Greenhouse Gas emissions payback for
lightweighted vehicles using aluminum and high-strength steel. Journal of Industrial Ecology
14 (2010) No. 6, pp. 929-946
[3] F. Klocke and D. Breuer: Using the finite element method and artificial neural networks to predict
ductile fracture in cold forming processes. Materials Processing and Design: Modeling,
Simulation and Applications, Columbus 2004, pp. 1944-1949
[4] M. Cockcroft, D. Latham: A simple Criterion Of Fracture For Ductile Metals. National
Engineering Laboratory, 1966
[5] Y. Bao, T. Wierzbicki: On fracture locus in the equivalent strain and stress triaxiality space.
International Journal of Mechanical Sciences 46.1 (2004): 81-98.
[6] J. Lemaitre: A Continous Damage Mechanics Model for Ductile Fracture. Journal of Engineering
Material and Technology, 107 (1985), pp. 83-89
[7] Y. Li; M. Luo, J. Gerlach; T. Wierzbicki: Prediction of shear-induced fracture in sheet metal
forming. Journal of Materials Processing Technology 210 (2010) No. 14, 1858-1869
[8] M. Dunant, D. Mohr: Effect of Lode parameter on plastic flow localization after proportional
loading at low stress triaxialities. Journal of the Mechanics and Physics of Solids 66 (2014),
pp. 133-153
[9] Q. Yin: Verfestigung- und Schädigungsverhalten von Blechwerkstoffen im ebenen
Torsionsversuch. Dissertation TU Dortmund, 2014
[10] N. Shamsaei, M. Gladskyi, K. Panasovskyi, S. Shukaev, A. Fatemi: Multiaxial fatigue of titanium
including step loading and load path alteration and sequence effects. International Journal of
Fatigue 32 (2010), No. 11, pp. 1862-1874
[11] ISO 12004-2:2008: Metallic materials – Sheet and strip – Determination of forming-limit curves
– Part 2: Determination of forming-limit curves in the laboratory
[12] N. Pathak, C. Butcher, M.J. Worswick, E. Bellhouse, J. Gao: Damage Evolution In Complex-
Phase and Dual-Phase Steels during Edge Stretching. Materials 10 (2017), No. 4, pp. 1-29
[13] T. Bergs, M. Nick, A. Feuerhack, D. Trauth, F. Klocke: Numerical-experimental Investigation Of
Load Paths in DP800 Dual Phase Steel During Nakajima Test. ESAFORM 2018, April 2018,
Palermo. Unpublished conference paper. AIP Publishing, 2018
[14] M. Calcagnotto, Y. Adachi, D. Ponge, D. Raabe: Deformation and fracture mechanisms in fine-
and ultrafine-grained ferrite/martensite dual-phase steels and the effect of aging. Acta
Materialia 59 (2011), pp. 658–670