Abstract
A reliable determination of the onset of void coalescence is critical to the modelling of ductile fracture. Numerical models have been developed but rely mostly on analyses on single defect cells, thus underestimating the interaction between voids. This study aims to provide the first extensive analysis of the response of microstructures with random distributions of voids to various loading conditions and to characterize the dispersion of the results as a consequence of the randomness of the void distribution. Cells embedding a random distribution of identical spherical voids are generated within an elastoplastic matrix and subjected to a macroscopic loading with constant stress triaxiality and Lode parameter under periodic boundary conditions in finite element simulations. The failure of the cell is determined by a new indicator based on the loss of full rankedness of the average deformation gradient rate. It is shown that the strain field developing in random microstructures and the one in unit cells feature different dependencies on the Lode parameter L owing to different failure modes. Depending on L, the cell may fail in extension (coalescence) or in shear. Moreover the random void populations lead to a significant dispersion of failure strain, which is present even in simulations with high numbers of voids.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.References
Bandstra JP, Koss DA (2008) On the influence of void clusters on void growth and coalescence during ductile fracture. Acta Mater 56:4429–4439. https://doi.org/10.1016/j.actamat.2008.05.009
Bao Y, Wierzbicki T (2004) On fracture locus in the equivalent strain and stress triaxiality space. Int J Mech Sci 46:81–98. https://doi.org/10.1016/j.ijmecsci.2004.02.006
Barsoum I, Faleskog J (2007) Rupture mechanisms in combined tension and shear-micromechanics. Int J Solids Struct 44:5481–5498. https://doi.org/10.1016/j.ijsolstr.2007.01.010
Barsoum I, Faleskog J (2011) Micromechanical analysis on the influence of the Lode parameter on void growth and coalescence. Int J Solids Struct 48:925–938. https://doi.org/10.1016/j.ijsolstr.2010.11.028
Becker R (1987) The effect of porosity distribution on ductile failure. J Mech Phys Solids 35:577–599. https://doi.org/10.1016/0022-5096(87)90018-4
Benzerga AA, Leblond JB (2014) Effective yield criterion accounting for microvoid coalescence. J Appl Mech 81:031009. https://doi.org/10.1115/1.4024908
Besson J (ed) (2004) Local Approach to Fracture. Les Presses de l’École des Mines, Paris
Besson J (2010) Continuum models of ductile fracture?: a review. Int J Damage Mech 19:3–52. https://doi.org/10.1177/1056789509103482
Besson J, Cailletaud G, Chaboche JL, Forest S, Blétry M (2009) Non-linear mechanics of materials. Solid mechanics and its applications. Springer, New York
Besson J, Foerch R (1998) Object-Oriented Programming Applied to the Finite Element Method Part I. General Concepts. Revue Européenne des Éléments Finis 7. Publisher: Taylor & Francis, pp 535–566. _eprint: https://doi.org/10.1080/12506559.1998.10511321
Bilger N, Auslender F, Bornert M, Michel JC, Moulinec H, Suquet P, Zaoui A (2005) Effect of a nonuniform distribution of voids on the plastic response of voided materials: a computational and statistical analysis. Int J Solids Struct 42:517–538. https://doi.org/10.1016/j.ijsolstr.2004.06.048
Bilger N, Auslender F, Bornert M, Moulinec H, Zaoui A (2007) Bounds and estimates for the effective yield surface of porous media with a uniform or a nonuniform distribution of voids. Eur J Mech A 26:810–836. https://doi.org/10.1016/j.euromechsol.2007.01.004
Boyce BL, Kramer SLB, Bosiljevac TR et al (2016) The second Sandia Fracture Challenge: predictions of ductile failure under quasi-static and moderate-rate dynamic loading. Int J Fract 198:5–100. https://doi.org/10.1007/s10704-016-0089-7
Boyce BL, Kramer SLB, Fang HE et al (2014) The Sandia fracture challenge: blind round robin predictions of ductile tearing. Int J Fract 186:5–68. https://doi.org/10.1007/s10704-013-9904-6
Brown L, Embury J (1973) Initiation and growth of voids at second phase particles. Inst Met (London) Monogr Rep Ser 1: 164–169
Brown MB, Forsythe AB (1974) Robust tests for the equality of variances. J Am Stat Assoc 69:364–367. https://doi.org/10.1080/01621459.1974.10482955
Cawthorne C, Fulton EJ (1967) Voids in irradiated stainless steel, vol 5115. Nature Publishing Group, Berlin, pp 575–576. https://doi.org/10.1038/216575a0
Coenen EWC, Kouznetsova VG, Geers MGD (2012) Novel boundary conditions for strain localization analyses in microstructural volume elements. In J Numer Methods Eng. https://doi.org/10.1002/nme.3298
Dæhli LEB, Morin D, Børvik T, Benallal A, Hopperstad OS (2020) A numerical study on ductile failure of porous ductile solids with rate-dependent matrix behavior. J Appl Mech 87:031014. https://doi.org/10.1115/1.4045524
Danas K, Ponte Castañeda P (2009) A finite-strain model for anisotropic viscoplastic porous media: I-Theory. Eur J Mech A 28:387–401. https://doi.org/10.1016/j.euromechsol.2008.11.002
Danas K, Ponte Castañeda P (2012) Influence of the Lode parameter and the stress triaxiality on the failure of elasto-plastic porous materials. Int J Solids Struct 49:1325–1342. https://doi.org/10.1016/j.ijsolstr.2012.02.006
Dunand M, Mohr D (2011) On the predictive capabilities of the shear modified Gurson and the modified Mohr–Coulomb fracture models over a wide range of stress triaxialities and Lode angles. J Mech Phys Solids 59:1374–1394. https://doi.org/10.1016/j.jmps.2011.04.006
Dunand M, Mohr D (2014) Effect of Lode parameter on plastic flow localization after proportional loading at low stress triaxialities. J Mech Phys Solids 66:133–153. https://doi.org/10.1016/j.jmps.2014.01.008
Fritzen F, Forest S, Böhlke T, Kondo D, Kanit T (2012) Computational homogenization of elastoplastic porous metals. Int J Plast 29:102–119. https://doi.org/10.1016/j.ijplas.2011.08.005
Fritzen F, Forest S, Kondo D, Böhlke T (2013) Computational homogenization of porous materials of Green type. Comput Mech 52:121–134. https://doi.org/10.1007/s00466-012-0801-z
Gamito M, Maddock S (2009) Accurate multidimensional poisson-disk sampling. ACM Trans Graph 29:8. https://doi.org/10.1145/1640443.1640451
Gao X, Zhang T, Hayden M, Roe C (2009) Effects of the stress state on plasticity and ductile failure of an aluminum 5083 alloy. Int J Plast 25:2366–2382. https://doi.org/10.1016/j.ijplas.2009.03.006
Gilioli A, Manes A, Giglio M, Allahverdizadeh N (2013) Effect of triaxiality and lode angle on the plasticity behaviour of a Ti-6Al-4V Titanium alloy. Key Eng Mater 577–578:413–416
Gologanu M, Leblond JB, Devaux J (2001) Theoretical models for void coalescence in porous ductile solids. II. Coalescence in columns. Int J Solids Struct 38:5595–5604. https://doi.org/10.1016/S0020-7683(00)00355-3
Green RJ (1972) A plasticity theory for porous solids. Int J Mech Sci 14:215–224. https://doi.org/10.1016/0020-7403(72)90063-X
Guo HJ, Ling C, Busso EP, Zhong Z, Li DF (2020) Crystal plasticity based investigation of micro-void evolution under multi-axial loading conditions. Int J Plast 129:102673. https://doi.org/10.1016/j.ijplas.2020.102673
Guo TF, Wong WH (2018) Void-sheet analysis on macroscopic strain localization and void coalescence. J Mech Phys Solids 118:172–203. https://doi.org/10.1016/j.jmps.2018.05.002
Gurson AL (1977) Continuum theory of ductile rupture by void nucleation and Growth: part I-yield criteria and flow rules for porous ductile media. J Eng Mater Technol 99:2–15. https://doi.org/10.1115/1.3443401
Heiman GW (2001) Understanding research methods and statistics: an integrated introduction for psychology. Houghton, Mifflin and Company, Boston
Helbert AL, Feaugas X, Clavel M (1996) The influence of stress triaxiality on the damage mechanisms in an equiaxed Ti-6AI-4V alloy. Metall Mater Trans A 27:3043–3058. https://doi.org/10.1007/BF02663853
Hure J (2021) Yield criterion and finite strain behavior of random porous isotropic materials. Eur J Mech A 85:104143. https://doi.org/10.1016/j.euromechsol.2020.104143
Hutchinson JW, Tvergaard V (2012) Comment on Influence of the Lode parameter and the stress triaxiality on the failure of elasto-plastic porous materials by K. Danas and P. Ponte Castañeda. Int J Solids Struct 49:3484–3485. https://doi.org/10.1016/j.ijsolstr.2012.07.009
Kanit T, Forest S, Galliet I, Mounoury V, Jeulin D (2003) Determination of the size of the representative volume element for random composites: statistical and numerical approach. Int J Solids Struct 40:3647–3679. https://doi.org/10.1016/S0020-7683(03)00143-4
Keralavarma SM, Reddi D, Benzerga AA (2020) Ductile failure as a constitutive instability in porous plastic solids. J Mech Phys Solids 139:103917. https://doi.org/10.1016/j.jmps.2020.103917
Khan IA, Bhasin V (2017) On the role of secondary voids and their distribution in the mechanism of void growth and coalescence in porous plastic solids. Int J Solids Struct 108:203–215. https://doi.org/10.1016/j.ijsolstr.2016.12.016
Khdir YK, Kanit T, Zaïri F, Naït-Abdelaziz M (2014) Computational homogenization of plastic porous media with two populations of voids. Mater Sci Eng A 597:324–330. https://doi.org/10.1016/j.msea.2013.12.095
Khdir YK, Kanit T, Zaïri F, Naït-Abdelaziz M (2015) A computational homogenization of random porous media: effect of void shape and void content on the overall yield surface. Eur J Mech A 49:137–145. https://doi.org/10.1016/j.euromechsol.2014.07.001
Koplik J, Needleman A (1988) Void growth and coalescence in porous plastic solids. Int J Solids Struct 24:835–853. https://doi.org/10.1016/0020-7683(88)90051-0
Kramer SLB et al (2019) The third Sandia fracture challenge: predictions of ductile fracture in additively manufactured metal. Int J Fract 218:5–61. https://doi.org/10.1007/s10704-019-00361-1
Lacourt L (2019) Étude numérique de la nocivité des défauts. Thèse de doctorat. Paris Sciences & Lettres, Paris
Lacourt L, Ryckelynck D, Forest S, Rancourt Vd, Flouriot S (2020) Hyper-reduced direct numerical simulation of voids in welded joints via image-based modeling. Int J Numer Methods Eng 121:2581–2599. https://doi.org/10.1002/nme.6320
Leblond JB, Mottet G (2008) A theoretical approach of strain localization within thin planar bands in porous ductile materials. Comptes Rendus Mécanique 336:176–189. https://doi.org/10.1016/j.crme.2007.11.008
Li XY, Hao Q, Shi YW, Lei YP, Marquis G (2003) Influence of mechanical mismatching on the failure of welded joints by void nucleation and coalescence. Int J Press Vess Pip 80:647–654. https://doi.org/10.1016/S0308-0161(03)00094-2
Ling C, Besson J, Forest S, Tanguy B, Latourte F, Bosso E (2016) An elastoviscoplastic model for porous single crystals at finite strains and its assessment based on unit cell simulations. Int J Plast 84:58–87. https://doi.org/10.1016/j.ijplas.2016.05.001
Liu ZG, Wong WH, Guo TF (2016) Void behaviors from low to high triaxialities: transition from void collapse to void coalescence. Int J Plast 84:183–202. https://doi.org/10.1016/j.ijplas.2016.05.008
Luo T, Gao X (2018) On the prediction of ductile fracture by void coalescence and strain localization. J Mech Phys Solids 113:82–104. https://doi.org/10.1016/j.jmps.2018.02.002
Matern B (1986) Spatial variation, 2nd edn. Lecture notes in statistics. Springer, New York. https://doi.org/10.1007/978-1-4615-7892-5
McClintock FA (1968) A criterion for ductile fracture by the growth of holes. J Appl Mech 35:363–371. https://doi.org/10.1115/1.3601204
Morin D, Blystad Dæhli LE, Børvik T, Benallal A, Hopperstad OS (2019) Numerical study of ductile failure under non-proportional loading. Eur J Mech A 74:221–241. https://doi.org/10.1016/j.euromechsol.2018.11.001
Morin L, Leblond JB, Benzerga AA (2015) Coalescence of voids by internal necking: theoretical estimates and numerical results. J Mech Phys Solids 75:140–158. https://doi.org/10.1016/j.jmps.2014.11.009
Morin L, Leblond JB, Benzerga AA, Kondo D (2016) A unified criterion for the growth and coalescence of microvoids. J Mech Phys Solids 97:19–36. https://doi.org/10.1016/j.jmps.2016.01.013
Needleman A, Tvergaard V (1992) Analyses of plastic flow localization in metals. Appl Mech Rev 45:3–18. https://doi.org/10.1115/1.3121390
Nguyen VD, Pardoen T, Noels L (2020) A nonlocal approach of ductile failure incorporating void growth, internal necking, and shear dominated coalescence mechanisms. J Mech Phys Solids 137:103891. https://doi.org/10.1016/j.jmps.2020.103891
Pedregosa F et al (2011) Scikit-learn: machine learning in python. J Mach Learn Res 12:2825–2830
Ponte Castañeda P, Willis JR (1995) The effect of spatial distribution on the effective behavior of composite materials and cracked media. J Mech Phys Solids 43:1919–1951. https://doi.org/10.1016/0022-5096(95)00058-Q
Rice JR, Tracey DM (1969) On the ductile enlargement of voids in triaxial stress fields. J Mech Physics of Solids 17:201–217
Rice JR (1976) The localization of deformation. In: Theoretical and Applied Mechanics, Proceedings of the 14th IUTAM congress, Delft. Ed. by WT Koiter. Delft: North-Holland Publishing Company: 207–220
Sarre B (2018) Influence du soudage laser Nd:YAG sur les propriétés métallurgiques et mécaniques de l’alliage de titane TA6V. Thèse de doctorat. Troyes: Université de Technologie de Troyes
Schöberl J (1997) NETGEN: an advancing front 2D/3D-mesh generator based on abstract rules. Comput Visual Sci 1:41–52
Shakoor M, Bernacki M, Bouchard PO (2015) A new body-fitted immersed volume method for the modeling of ductile fracture at the microscale: analysis of void clusters and stress state effects on coalescence. Eng Fract Mech 147:398–417. https://doi.org/10.1016/j.engfracmech.2015.06.057
Shakoor M, Bernacki M, Bouchard PO (2018) Ductile fracture of a metal matrix composite studied using 3D numerical modeling of void nucleation and coalescence. Eng Fract Mech 189:110–132. https://doi.org/10.1016/j.engfracmech.2017.10.027
Tekoğlu C, Hutchinson JW, Pardoen T (2015) On localization and void coalescence as a precursor to ductile fracture. Philos Trans R Soc A 373:20140121. https://doi.org/10.1098/rsta.2014.0121
Thomason PF (1985) Three-dimensional models for the plastic limit-loads at incipient failure of the intervoid matrix in ductile porous solids. Acta Metall 33:1079–1085. https://doi.org/10.1016/0001-6160(85)90201-9
Torki ME (2019) A unified criterion for void growth and coalescence under combined tension and shear. Int J Plast 119:57–84. https://doi.org/10.1016/j.ijplas.2019.02.002
Trejo Navas VM, Bernacki M, Bouchard PO (2018) Void growth and coalescence in a three dimensional non-periodic void cluster. Int J Solids Struct 139–140:65–78. https://doi.org/10.1016/j.ijsolstr.2018.01.024
Tvergaard V (2016) Effect of void cluster on ductile failure evolution. Meccanica 51:3097–3105. https://doi.org/10.1007/s11012-016-0537-5
Tvergaard V (2017) Nucleation from a cluster of inclusions, leading to void coalescense. Int J Mech Sci 133:631–638. https://doi.org/10.1016/j.ijmecsci.2017.09.027
Tvergaard V, Needleman A (1984) Analysis of the cup-cone fracture in a round tensile bar. Acta Metall 32:157–169. https://doi.org/10.1016/0001-6160(84)90213-X
Vincent PG, Monerie Y, Suquet P (2009) Porous materials with two populations of voids under internal pressure: I. Instantaneous constitutive relations. Int J Solids Struct 46:480–506. https://doi.org/10.1016/j.ijsolstr.2008.09.003
Virtanen P et al (2020) SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat Method 17:261–272. https://doi.org/10.1038/s41592-019-0686-2
Vishwakarma V, Keralavarma SM (2019) Micromechanical modeling and simulation of the loading path dependence of ductile failure by void growth to coalescence. Int J Solids Struct 166:135–153. https://doi.org/10.1016/j.ijsolstr.2019.02.015
Wong WH, Guo TF (2015) On the energetics of tensile and shear void coalescences. J Mech Phys Solids 82(259–286):2020. https://doi.org/10.1016/j.jmps.2015.05.013
www.zset-software.com (2020)
Xiao X, Mu Z, Pan H, Lou Y (2018) Effect of the Lode parameter in predicting shear cracking of 2024-T351 aluminum alloy Taylor rods. Int J Impact Eng 120:185–201. https://doi.org/10.1016/j.ijimpeng.2018.06.008
Zhai J, Luo T, Gao X, Graham SM, Baral M, Korkolis YP, Knudsen E (2016) Modeling the ductile damage process in commercially pure titanium. Int J Solids Struct 91:26–45. https://doi.org/10.1016/j.ijsolstr.2016.04.031
Zhang KS, Bai JB, François D (1999) Ductile fracture of materials with high void volume fraction. Int J Solids Struct 36:3407–3425. https://doi.org/10.1016/S0020-7683(98)00157-7
Zhang K, Badreddine H, Saanouni K (2020) Ductile fracture prediction using enhanced CDM model with Lode angle-dependency for titanium alloy Ti-6Al-4V at room temperature. J Mater Proces Technol 277:116462. https://doi.org/10.1016/j.jmatprotec.2019.116462
Zhu JV, Ben BM, Abed-Meraim F (2020) Investigation of the competition between void coalescence and macroscopic strain localization using the periodic homogenization multiscale scheme. J Mech Phys Solids 143:104042. https://doi.org/10.1016/j.jmps.2020.104042
Zhu Y, Engelhardt MD, Kiran R (2018) Combined effects of triaxiality, Lode parameter and shear stress on void growth and coalescence. Eng Fract Mech 199:410–437. https://doi.org/10.1016/j.engfracmech.2018.06.008
Acknowledgements
The authors would like to thank L. Lacourt for valuable discussion and his help in the analysis of localization modes.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
About the \(\delta \) indicator
This appendix provides several complements about the \(\delta \) indicator. Its expression is first derived by computing the homogeneous deformation of Green matrix. An example of application is then presented. Finally a sensitivity analysis regarding the threshold coefficients is carried out.
1.1 Derivation of the expression for \(\delta \)
Consider a perfectly plastic volume element (neglecting here the elasticity) which deforms homogeneously when subjected to the loading conditions (13). In order to simply represent the porous nature of the cell, the material behavior will obey Green’s (1972) isotropic yield criterion (also used by Fritzen et al. (2013)):
with C a constant (\(C=0\) corresponds to a von Mises material, and for \(C=1/2\), there is no lateral contraction of the cube in tension). The other equations in Eq. (4) are unchanged, but they are applied here to macroscopic quantities.
As the material behavior is isotropic, stays diagonal in the diagonalizing basis of . Then can be written as:
where \(b_2\) and \(b_3\) are functions to be determined. As , and is diagonal, can be written as:
For a perfectly plastic Green material, the behavior law in (4) reads:
is diagonal so there are three constants \(\alpha _1\), \(\alpha _2\) and \(\alpha _3\) such that:
Combining (25) and (28) yields the system:
The plastic multiplier is then \({\dot{p}}=\frac{{\dot{\epsilon }}/\alpha _1}{1+{\dot{\epsilon }}t}\) and the differential equations can be solved with the initial conditions \(b_2(0)=1\), \(b_3(0)=1\):
Finally,
The function comparing the behavior of and the homogeneous plastic deformation case is then:
The \(\delta \) criterion used throughout the article is recovered by setting \(C=0\), which corresponds to the simplified case of a von Mises material. In this case, the criterion depends no more on the applied .
The evolution of \(\delta \) for a simulation with \(T=1\) and \(L=-1\) (coalescence in uniaxial strain state) is shown in Fig. 15, for two values of C: 0 and 1/2. For both values of C, the vanishing of \(\delta _C\) is simultaneous with the stabilization of transverse displacement. However the sharp drop of \(\delta _C\) allows a more precise numerical determination of the onset of coalescence than the more progressive stabilization of the transverse strain. For \(C=1/2\), \(\delta _{1/2}\) is approximately constant at the beginning of the simulation, so that the hypothesis of homogeneous flow in a Green volume element (taking into account the porosity) well represents the overall behavior of the cell with a von Mises matrix. However, with \(C=0\), \(\delta _0\) does not depend anymore on the stress state, while still keeping the sudden drop of \(\delta _C\) necessary for the determination of the coalescence onset.
1.2 Sensibility analysis regarding the threshold coefficients
Finally we verify that the \(\delta \) indicator is a reliable indicator of failure by assessing its sensitivity to the choice of the empirically chosen threshold values. As the Eq. (16) shows, the determination of the onset of coalescence relies on two thresholds: a relative one A, which compares the current value of \(\delta \) to its maximum, and an absolute one B mostly active in shear-like conditions. The values for those were chosen as \(A=0.05\) and \(B=0.005\) but a robust indicator should not be too sensitive to these values.
Figure 16 compares the effect of different A and B values on the \(E_c-L\) curve (common triaxiality \(T=1\), microstructure R1). At constant B, the effect of A is only visible in the HLEMZ and the LLEMZ, and generally negligible. At constant A, B only affects the coalescence strain values in the SMZ. Although a change in B can modify the strain by 0.05, the global aspect of the curve is preserved. The determination of failure by the indicator therefore appears robust with respect to changes in the coefficients.
Effect of different meshing parameters and boundary conditions
In this section, we review the simulation hypotheses and assess their influence on the results presented up to now, showing therefore how representative the results are and how far they can be generalized. First we verify that finite element discretization effects can be neglected, and investigate the effect of different boundary conditions.
1.1 Effect of the meshing parameters
All the simulations described up to now were carried out on meshes of cells with the same meshing parameter. To determine the influence of mesh size on coalescence results, the same microstructure R1 was meshed with different meshing parameters \(h_{cell}/r_0\in \{1.25, 1, 0.875, 0.625\}\) (with the notation of Sect. 2.1). The maximum element size near the voids is also adapted to keep the ratio \(h_{cell}/h_{void}=5\) constant. The same loading condition \(T=1,\,L=-1\) is applied to the four meshes. Figure 17 shows that stress values during the simulations differ between the meshes, but the relative difference between the finest and coarsest meshes is about \(5\%\), which remains acceptable. The onset of coalescence \(E_c\) which is our main quantity of interest, is almost identical between the meshes, at \(E_c=0.33\pm 1\%\). Therefore the influence of mesh refinement for random microstructure cells appear limited (although there was only a ratio of 2 between the element sizes of the coarsest and the finest mesh), which justifies the value \(h_{cell}=0.08\) adopted throughout this study.
1.2 Effect of the boundary conditions
We here investigate the influence of boundary conditions. The results from Sect. 4 are first compared to those obtained with different boundary conditions. Namely we investigate the influence of conditions on the average gradient, and of planar faces conditions. The consistency of results at \(L=-1\) is also checked by a comparison with simulations on axisymmetric cells.
The conditions imposed on the average gradient to prevent rigid body motion are first investigated. In Sect. 2.3 we imposed symmetric, as for Ling et al. (2016). However another reasonable choice would be to fix some degrees of freedom at the vertices of the cubic cell, as depicted in Fig. 18, which is the standard method for boundary value problems. A vertex is already fixed in order to prevent translations, but by fixing two degree of freedom on a second one, and a last one on a third vertex, all rotations are fixed. This can be reformulated as:
i.e. is an upper triangular matrix. Due to the mixed conditions imposed by the macroscopic spring element, the results from the symmetric case cannot be easily transposed to the triangular case. These two choices lead to distinct proportional loading path classes and should therefore be compared.
On the microstructure R1, at fixed triaxiality \(T=1\), simulations were performed for several Lode parameters to compare the two sets of conditions on (Fig. 19a). The evolution of \(E_c\) is close between the two types of conditions, and the same ductility zones can be identified for the triangular gradient condition. However, in that case, cusps seem to be less pronounced than for a symmetric gradient; this may be due to the different treatment of shear components by the two conditions. Therefore the influence of the conditions on remains limited.
We then compare the effects of periodic and parallel faces boundary conditions. Parallel faces conditions mean that the cubic cell retains parallel flat faces throughout the computation (for instance all the points on the \(x_0=0\) face have the same x-displacement). This condition is more constraining than periodic boundary conditions. As the comparison in Fig. 19b shows, the two types of conditions lead to qualitatively different responses. For the parallel faces, no separation between three ductility zones can be seen (except near \(L=0\)) and the response of the random microstructure is closer to that typical of the unit cell. Moreover no decrease of ductility near \(L=0\) is observed for the unit cell. Results for the unit cell differ between the parallel faces and periodic boundary conditions, because in the periodic case, faces are allowed not to remain strictly parallel and planar. On the contrary, boundary conditions made of parallel sides strongly hinder the shear mode failure and only the extension mode remains possible. The competition between these two modes tends to postpone failure (see the cusps on Fig. 4). Therefore, the reduced competition between modes may explain an earlier coalescence for parallel unit cells. The preceding results show that boundary conditions exert a strong influence on the response of the cell.
Finally, the consistency of results obtained at \(L=-1\) is checked. As this type of loading is axisymmetric, a computation with a 2D axisymmetric unit cell was also performed for comparison. Such unit cells are frequent in ductile fracture studies (Morin et al. (2015) for instance). The diameter and the height of the cylinder were chosen equal to \(L_{cube}\). The porosity is still \(6\%\), so the radius of the void was modified to \(0.22 L_{cube}\). The boundary conditions for this cell differ slightly from those described in Sect. 2.3: they are no more periodic and are replaced by straight edges conditions. Besides the virtual constant triaxiality element is not linked to the average deformation gradient but to the displacement of the top left node.
Figure 20 compare results for the unit cell and the 2D axisymmetric cell at varying stress triaxiality for \(L=-1\). In this type of loading, the unit cell was shown in Sect. 4 to exhibit the same behavior as random microstructures, compatible with a Rice-Tracey evolution. At initial porosity \(f_0=6\%\), the axisymmetric cell presents however a significantly higher exponent (in absolute value) for the evolution of \(E_c\) with respect to T. This effect seems due to the relatively high porosity in the axisymmetric cell: as depicted in Fig. 20, the evolution of \(E_c\) with respect to T for the low porosity \(f_0=1\%\) axisymmetric cell is much closer to the one predicted by Rice and Tracey.
Rights and permissions
About this article
Cite this article
Cadet, C., Besson, J., Flouriot, S. et al. Ductile fracture of materials with randomly distributed voids. Int J Fract 230, 193–223 (2021). https://doi.org/10.1007/s10704-021-00562-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10704-021-00562-7