Abstract
Superior mechanical and thermal properties, high wear resistance and a competitive price of compacted graphite iron (CGI) have made it an integral part of industry worldwide. In its applications in automotive engines, high-temperature environments cause thermal expansion that can result in emergence of interfacial damage in CGI. Although graphite-matrix interfacial damage is considered the main damage mechanism that can lead to total fracture of CGI, extensive research on CGI has not yet fully investigated this phenomenon at the microscale, especially under pure thermal loading. This paper focuses on the high-temperature performance of CGI and the onset of damage in graphite in thermal cycles. Three-dimensional numerical models are developed, with a single graphite inclusion embedded in a unit cell of the metallic matrix. Elastoplastic behaviour is considered for both phases in simulations. The effects of morphology and orientation of graphite inclusions on a response of an entire unit cell to thermal loading are investigated. Also, the influence of periodic and fully-fixed boundary conditions on the damage behaviour of CGI is discussed. The results can give a better understanding of the fracture mechanisms of CGI exposed to elevated temperatures.
Similar content being viewed by others
1 Introduction
Compacted graphite iron was first used in 1948 (Dawson 2008) and is considered an important metal matrix composite in the industry thanks to its beneficial physical and mechanical properties, thermal conductivity, and competitive price. CGI has wide applications in pipes and machinery, as well as in the automotive industry for disc brakes (Behera 2012), cylinder heads and blocks (Dawson 2008).
As a typical in-situ metal matrix composite, cast iron is composed of graphite, pearlite and ferrite (Chen and Chen 2011). The microstructure of CGI includes graphite inclusions and a metallic matrix, two phases with considerably different thermomechanical properties. Based on their circularity, graphite particles are classified as flake, vermicular and spherical (Fig. 1).
Graphite-matrix debonding is experimentally observed as the main damage mechanism of compacted graphite iron (Di Cocco et al. 2010; Qiu et al. 2016a). Graphite is brittle and soft, generally separating from the metallic matrix at low stress levels (Di Cocco et al. 2014; Norman and Calmunger 2019) and resulting in plastic deformation of the matrix (Di Cocco et al. 2010). This interfacial-debonding (decohesion) phenomenon initiates microcracks in the matrix; the process of coalescence of such microcracks may cause main cracks or macrocracks (Di Cocco et al. 2010). Furthermore, under cyclic mechanical loads (Di Cocco et al. 2013), interfacial debonding can cause the initiation of fatigue cracks (Kohout 2001; Endo and Yanase 2014). Also, microcracks between the graphite particles and the matrix can result in yielding in cast irons (Hornbogen 1985; Wang et al. 2007; Qiu et al. 2016a).
In diesel engines that are manufactured in many cases of cast irons, the inner wall of the combustion chamber can reach 400–600 °C (Tang et al. 2022). Due to such high-temperature service conditions of automotive engines, CGI experiences intense thermal loadings. The mismatch in coefficients of thermal expansion between the graphite particles and the matrix can lead to interfacial debonding between these two phases, resulting in microcracks even without mechanical loads. The propagating microcracks can form a network and result in a main crack under high temperatures (Qiu et al. 2016a). Also, CGI softens when exposed to severe thermal loadings in engines enhancing debonding at lower stress. At room temperature, debonding occurs above 495 MPa (Wu et al. 2019). On the contrary, it was observed at stress levels below 50 MPa at 723 K (450 °C) (Qiu et al. 2016b). At the macroscale, the tensile strength of CGI decreases when the temperature surpasses 300 °C (Selin 2010). Therefore, the mechanical and thermal properties of graphite and matrix are significantly affected by thermal loads, which requires further investigation.
Morphology, size and distribution of the graphite phase are important factors in determining the damage mechanisms of cast irons (Di Cocco et al. 2014; Zhang et al. 2018; Salomonsson and Jarfors 2018; Zhan et al. 2022). The higher content and size of vermicular graphite particles cause easier initiation of cracks (Qiu et al. 2016a). Also, it was found that longer graphite inclusions reduced the fatigue life (Kihlberg et al. 2021). Both stiffness and toughness of composites changes with increasing volume fraction and aspect ratio of filler (Safaei et al. 2015). But in contrast to polymer-based composites, in cast irons the reinforcement is a weaker and softer constituent. Under increasing environmental temperature, the Young’s modulus of nodular cast iron remains constant whereas its ultimate tensile strength and yield stress decrease (Šamec et al. 2011).
Thanks to their various applications in industry, cast irons were researched for many decades with the development of science and technology. Advanced equipment for observation such as microscopy provided opportunities for researchers to examine and explore the surfaces of cast irons at microscopic scale in 1863 (Stefanescu 2019). The first macroscopic constitutive model for grey cast iron was developed in 1976, providing an accurate load to predict its failure (Frishmuth and McLaughlin 1976). In modern research, the thermomechanical behaviour of cast iron is numerically analysed mainly with phenomenological and micromechanical modelling schemes. Based on phenomenological descriptions, the yield surface and hardening parameters were modified to accommodate microstructural features (Frishmuth and McLaughlin 1976). The micromechanical methodology, on the other hand, focuses on predicting the effect of microstructure of a material directly from experimental observations (Andriollo et al. 2019; Yang et al. 2021). In addition, the constituents were usually assumed to be anisotropic (Andriollo et al. 2016) or elastoplastic (Andriollo et al. 2015).
So far, interfacial debonding was studied mainly under mechanical loads and there is not enough information in the literature about the thermal damage mechanism of CGI at the microscale under pure thermal loading (Palkanoglou et al. 2022). Due to differences in the coefficients of thermal expansion between its main constituents, CGI is vulnerable to high temperatures. However, because of its complex microstructure, it is difficult to study thermal debonding thoroughly only with experimental or microstructural studies. Therefore, this work assesses the effects of graphite morphology on high-temperature behaviours of CGI under thermal loading with three-dimensional numerical models. The inputs of into these models were obtained from the statistical characterisation of SEM images and mechanical testing.
2 Methodology
2.1 Microstructural characterisation of CGI
Statistical analysis of microstructure was performed for CGI specimens (EN-GJV-450). A set of 20 images was acquired using scanning electron microscopy (SEM) and analysed with ImageJ software. The image processes in ImageJ include scale setting, threshold adjustment, and outline shape adjustment of graphite (Fig. 2). The obtained results present the microstructure characterisation of CGI in terms of the perimeter, area, volume fraction, circularity, lengths of major and minor axes of graphite inclusions.
Based on the obtained 2D micrographs, the following assumptions about the geometry of both graphite and matrix in 3D numerical simulations were adopted:
-
nodular (spherical) graphite in 2D was simulated as a sphere in 3D;
-
vermicular graphite particles in 2D corresponded to regular ellipsoids in 3D;
-
the major axis of the vermicular graphite was considered equal to the diameter of a nodular graphite;
-
the rectangular matrix domain in 2D was assumed as a cube with equal length, width, and height in 3D.
A set of unit cells in 3D was created in Abaqus (Version 2021) using the methodology of representative volume element (RVE), where a graphite inclusion was embedded into a cubic metallic matrix domain (Fig. 3). An RVE with equivalent morphology could represent a constituent at the macroscale using a micro-unit cell (Zhan et al. 2022).
2.2 Numerical models
2.2.1 Geometry
Based on the assumptions for the microstructure characterisation results, the dimension of the spherical graphite was calculated with the following relations:
The 2D parameters of nodular graphite particle was obtained as follows:
The 3D dimensions of the numerical model follow the expressions:
Based on Eqs. (1)–(3) and (4)–(6), the dimension of metallic matrix domain in the 3D numerical model was acquired with the following equations:
The aspect ratio of vermicular graphite particle \({A}_{\text{R}}\) was acquired as follows:
The numerical models developed comprised a metallic matrix cube containing a nodular or a vermicular graphite inclusion. The dimensions of developed FEA models were established according to the range of spherical graphite diameters, as measured from the CGI micrographs. For nodular graphite, the diameter varied between 0.29–42.53 μm; hence, 15 μm was selected. For vermicular graphite particles, a major axis of 15 μm length was used; the maximum aspect ratio in the statistical analysis was 5, thus, the corresponding minor-axis length was 3 μm. Based on the 2D analysis, the dimension of cubic matrix in 3D was in the range of 25.74–32.23 μm; therefore, 30 μm was chosen in this research for comparability. The volume fraction of nodular graphite was 6.5%, in the range of statistical results (5.2–11.37%) (Palkanoglou et al. 2020). The dimensions of spherical and vermicular particles in the matrix are given in Table 1 for all the studied models. To assess the effects of the orientation of graphite inclusions on the mechanical behaviour of CGI, different orientations of graphite (at 0° and 45° with respect to the X axis) in the central XY plane of the cubic matrix were selected (Fig. 3). The length of diameter for the nodular graphite inclusion and the major axis for the vermicular graphite inclusion was 15 μm, and the distance between boundary of particle and the RVE face was 7.5 μm. Following the path AB of measurement for presentation of results, the region from 0 to 7.5 μm is inside the graphite domain, while the matrix region stretches from 7.5 to 15 μm in all diagrams (Fig. 3).
2.2.2 Constitutive behaviour
As a soft and brittle material, graphite was considered to display a limited plastic behaviour (Seldin 1966; Greenstreet et al. 1973; Andriollo et al. 2015). The graphite and matrix domain of the model had the properties of CGI used in the previous 2D study in Table 2 (Palkanoglou et al. 2020).
The mechanical behaviours of both graphite and the metallic matrix in this research were described with the J2-flow theory (Palkanoglou et al. 2020). For the graphite phase, beyond the elastic region, a damage model was employed. The initiation of damage was considered when a plastic-strain-based criterion was fulfilled. When the plastic displacement approaches the displacement at failure, the damage initiation was considered to make sure the debonding of graphite happens at the determined temperature matching the SEM experiments. After the onset of damage, the stiffness of the material degraded gradually to 0. The integral ductile damage criterion is given by the following equation (Hooputra et al. 2004)
where \({\omega }_{\text{D}}\) is the state variable increasing monotonically. At each increment, the state variable increment:
where \({\overline{\varepsilon }}_{\text{D}}^{\text{pl}}\) is the plastic strain at the onset of damage.
Damage was considered to evolve linearly with deformation after its initiation. Element deletion occurred when the value of damage exceeded a set threshold at all integration points of a given element. This damage criterion was only used in graphite inclusion and applied to every element in the graphite. The four neighbouring elements surrounding the node with maximum damage were selected. In simulations, the damage parameter was calculated as the average for the four neighbouring elements with maximum damage.
2.2.3 Boundary and loading conditions
Eight-node 3D stress elements (C3D8) were selected to mesh the models. After conducting a mesh convergence study, 1 μm was selected as the size of elements (Fig. 4). 26,400 elements were used leading to 28,957 nodes in model A.
Periodic boundary conditions (PBCs) and fully-fixed boundary conditions (FFBC) were considered in this study to represent the extreme cases of the in-service constraints. PBCs control the pairs of nodes on the corresponding surfaces with the same level of displacement in opposite directions. These conditions are widely applied in finite-elements analysis to simulate the microscopic and mesoscopic behaviours of the constituents with RVE (Garoz et al. 2019) and allow simulations of deformation of the RVE model using a small domain that represents the corresponding infinitely large system permitting the distortion of boundary surfaces (Omairey et al. 2019). Under PBCs, a pair of any two points \(x\) and \(x+d\) on the two corresponding boundary surfaces with distance \(d\) should meet the condition
where \({\mathbf{u}}\) is the displacement and \({\mathbf{t}}\) is the traction at \(x\), respectively; \(\overline{\varepsilon }\) is the average infinitesimal strain over the volume (Drago and Pindera 2007). In some simulations, fully-fixed boundary conditions (FFBCs) were applied at all the nodes along the boundary surfaces of the cubic metallic matrix, constraining all six degrees of freedom of these nodes. For each morphology of graphite unit cell, PBCs or FFBCs were applied for comparability. Fully fixed and periodic boundary conditions can give the upper and lower bounds of the realistic case in practise. The realistic case should be somewhere in between. Thermal loading was used for all unit cells with a linear increase from 25 to 500 °C. Specially, Model A2 was also followed by cooling down to 25 °C and then under further multiple thermal loading cycles (single cycle of 25–500–25 °C).
In summary, the details of FEA models are shown in Table 3:
Thus, models A1 to A2 described nodular graphite particles, with all other models dealing with vermicular ones.
Complex 3D shapes of vermicular and flake graphite particles in CGI—as opposed to inclusions in nodular cast irons—prevent the establishment of general features of its thermomechanical behaviour at microscale. In contrast to direct introduction of such shapes (obtained, e.g., with computed tomography) into numerical simulations, this study focuses on the simplified–but controlled–morphology of the inclusions. Together with the concept of RVE, this allows the assessment of the geometrical features on evolution of thermal stresses and strains at microscale and distribution of the resultant damage.
2.2.4 Validation of numerical models
A tensile test was implemented on a CGI (EN-GJV-450) specimen at room temperature and the measured mechanical response was compared with results obtained with the developed FEA models for validation in Fig. 5. The mechanical data for the matrix domain in simulations were based on the literature (Niu et al. 2021). The tensile loading and PBCs were applied to the 3D numerical model at room temperature to emulate the experiments. The validation of the proposed models was implemented by the comparison between numerical simulations and tensile-test experiments since the interaction in microstructure was difficult to measure directly. Monitoring of three-dimensional particle debonding with SEM is not possible. Some in-situ SEM studies of microspecimens loaded with miniature load cells allowed 2d observations of separation between the nodular graphite particles and the matrix in ferritic–pearlitic ductile cast only on the specimen’s surface (Iacoviello et al. 2008).
3 Results and discussions
3.1 Effect of particle shape orientation
Results for the effect of graphite orientation were obtained with the respective models; the temperature of damage initiation in graphite in models with FFBCs is shown in Fig. 6. The graphite with a larger aspect ratio exhibited damage at lower temperatures since it is much sharper than the spherical particle (an aspect ratio of 1). These nodular particles are more difficult to damage. Damage appeared earlier in the thermal cycle as the inclusion aspect ratio increased, in a monotonic but non-linear way. The trend was similar for both graphite orientations, 0° and 45°. It can be concluded that damage in graphite initiates at lower temperatures in less spherical particles while graphite orientation does not affect this process significantly when FFBCs are used, to represent behaviour of constrained parts of the components exposed to the thermal load.
3.2 Effect of graphite morphology
The evolution of the maximum damage in graphite particles with different aspect ratios is shown in Fig. 7 for FFBCs. As the aspect ratio of the inclusion increased, and the damage due to the application of thermal load appeared earlier. The cases with spherical graphite demonstrated higher strength, with the damage occurring at higher temperatures. The evolution of damage was affected by the inclusion shape with the sharper-shaped particles (i.e., higher aspect ratios) exhibiting a more rapid damage evolution. The trends of damage evolution differed significantly even for small variations in the inclusion’s aspect ratio.
Under FFBCs, the unit cell was fully constrained; thus, at heating, the graphite and matrix were compressed by each other due to the mismatch in the coefficients of thermal expansion. The graphite with a smaller aspect ratio has higher circularity and larger volume so that it can withstand a higher compressive load. With an increase in the aspect ratio, the stress–strain behaviour of graphite particles became nonlinear due to the initiation of graphite-matrix decohesion. In contrast, the decreased aspect ratio reduces of the level of compressive stresses at FFBCs and made their mechanical response almost linear. Kinematically constrained graphite particles continue resisting compression even after local failures, under nearly hydrostatic loading; the current model does not account for the local cracking in compressed graphite.
Spatial distribution of damage in graphite inclusions with different morphologies is shown in Fig. 8 after heating to 500 °C. It was found that damage was mainly localised near the interface boundary between the particle and the matrix. For a higher aspect ratio, the volume of graphite inclusion in the direction of the Z axis in model E1 was smaller than that in model A1 (for the spherical particle) leading to a higher volume of the metallic matrix (Fig. 8). As a result, the higher compressive load from the matrix under FFBCs caused a more significant damage in the graphite inclusion in model E1 in Z direction.
3.3 Effect of boundary conditions
The effect of boundary conditions was studied by comparing the distributions of normal stress along the major axis of the graphite inclusion, from the graphite centre to the edge of the analysis domain (see path AB in Fig. 3a). In FFBCs, the mismatch caused compressive normal stresses in both constituents. The increased aspect ratio resulted in significantly higher stress levels. Transition to the PBCs that allow the thermal expansion of the unit cell fully changed the character of the normal stresses, generating tensile stresses both in the graphite particles and the matrix. Also, the charter of their distribution depending on the aspect ratio changed: stresses in the nodular particle were lower. With the increasing aspect ratio, the stresses grew, while the situation in the matrix was reversed: for the low aspect ratio, the normal stresses were smaller. For the nodular particle, for the PBCs, they were the highest in the matrix.
The calculated levels of the normal stresses in the graphite particles were rather high, potentially in excess of its strength and close to the maximum values observed in the literature (AZOM 2022). Hence, the future models should incorporate the account for the potential fracture of graphite.
The distribution of von Mises stress along the path AB in the models with FFBCs is presented in Fig. 9a. Under thermal loading, the thermal expansion of both phases is restricted by the fixed boundaries resulting in compressive stress. This interaction is also affected by the mismatch of thermomechanical parameters of the two constituents and the particle morphology. When the inclusion shape was closer to spherical, the stress state was effectively hydrostatic, i.e., with vanishing effective stresses. As a result, the von Mises stress in spherical graphite was rather low in all cases. In contrast, its magnitude in the matrix increased at the periphery of the unit cell with the aspect ratio.
The distribution of von Mises stresses in CGI changes for PBCs are in Fig. 9b; being able to expand, the unit cell did not develop excessive thermal stresses. As a result, the developed stresses in the matrix at the maximum temperature were far below the yield point even in the case of vermicular graphite inclusion. The trend of all curves for different aspect ratios was similar, characterised by considerable gradients in parts of the matrix close to the interface.
In FFBCs, the thermal stresses were still present after the cooling-down stage. The highest level of von Mises stress was 460 MPa at 500 °C and 553 MPa at 25 °C surrounding graphite particle (Fig. 10); this results in the plastic strain of the matrix.
In case of PBCs, the mismatch of thermal expansion between the two constituents led to stress up to 160 MPa after the increase in the temperature (Fig. 11a). Residual stresses were observed in the model after each cycle; however, the level under PBC was negligibly small (up to 6.6 MPa) because of the relatively unconstrained thermal expansion after cooling down to 25 °C (see the difference in the scales at Fig. 11). The stresses in the models with PBC were lower than the yield point due to non-restricted boundary conditions. Thus, graphite damage was not observed in these PBC cases.
3.4 Outlook
The current model did not consider interfacial damage processes directly, although some analysis for formulations has been done by authors previously. This is related to the high intensity of the three-dimensional simulation using cohesive zone elements. The next step will be the introduction of these types of elements into simulations and to consider the interface between graphite and matrix.
The present paper approached 3D simulations and used the same damage criterion of graphite without interface. The simplified model without interface would be the basic stage of the series of model. The interface will be considered in the future models. However, direct measurement of interfacial damage properties is not possible and needs development of some additional experimental methods.
From the simplified models, damage was assumed to appear on the outer layer of graphite because of its brittle properties and inability to accommodate plastic deformation. Metallic matrix was generally considered isotropic and ductile. The ductile damage criterion has been used in ferrite matrix phase of cast iron because of the high ductility and moderate yield strength of matrix (Andriollo and Hattel 2016). However, the limitations of ductile damage criterion include its ability to properly describe the plastic volumetric strain (Murakami 2012).
Pure thermal load results in high levels of hydrostatic stresses, the J2 based plastic model is not fully adequate. Thus, the authors consider looking into other plastic potentials (such as Drucker-Prager) for their future simulations.
4 Conclusions
In this work, the effects of morphology and orientation of graphite inclusions as well as the boundary conditions of the studied domain on damage and thermomechanical behaviour of CGI were investigated by developing and employing a set of three-dimensional numerical models. The main conclusions are as follows:
-
The temperature of damage initiation in graphite were similar for both orientations (0° and 45°). Hence, orientation of graphite particles did not affect the damage onset significantly.
-
In contrast, the morphology of graphite influenced the onset of damage and plasticisation. When the aspect ratio of graphite inclusion increased, the initiation of damage in graphite started earlier and develops with increasing temperature.
-
Two types of used boundary conditions—PBCs and FFBCs—provide the lower and upper boundaries for graphite damage mechanisms, respectively. The results of models with FFBCs are more conservative than PBCs in FEA models (models A to E).
This work demonstrated a consistent trend for the graphite damage mechanism for various three-dimensional shapes with the change in the aspect ratios. The next step would be to consider the effect of deviation of real-life shapes of graphite particles from ellipsoids used in this study on thermomechanical behaviour of modern cast irons.
Abbreviations
- \(A\) :
-
Area of graphite (2D)
- \({A}_{\text{R}}\) :
-
Aspect ratio of vermicular graphite
- \(D\) :
-
Diameter of spherical graphite
- \(\text{FFBC}\) :
-
Fully-fixed boundary condition
- \({L}_{2\text{D}}\) :
-
Length of square matrix in 2D
- \({L}_{3\text{D}}\) :
-
Length of cubic matrix in 3D
- \({L}_{\text{major}}\) :
-
Major axis of vermicular graphite
- \({L}_{\text{minor}}\) :
-
Minor axis of vermicular graphite
- \(\text{PBC}\) :
-
Periodic boundary condition
- \(p\) :
-
Pressure
- \(q\) :
-
Von Mises equivalent stress
- \(S\) :
-
Area of square matrix (2D)
- \(V\) :
-
Volume of cubic matrix in 3D
- \({V}_{\text{f}}\) :
-
Volume fraction of graphite
- \({V}_{\text{g}}\) :
-
Volume of spherical graphite in 3D
- \({\overline{\varepsilon }}_{\text{D}}^{\text{pl}}\) :
-
Plastic strain at the onset of damage
- \({\omega }_{\text{D}}\) :
-
State variable increasing monotonically
References
Andriollo T, Hattel J (2016) On the isotropic elastic constants of graphite nodules in ductile cast iron: analytical and numerical micromechanical investigations. Mech Mater 96:138–150. https://doi.org/10.1016/j.mechmat.2016.02.007
Andriollo T, Thorborg J, Tiedje NS, Hattel J (2015) Modeling of damage in ductile cast iron—the effect of including plasticity in the graphite nodules. In: IOP conference series: materials science and engineering, vol 84, p 12027. https://doi.org/10.1088/1757-899x/84/1/012027
Andriollo T, Thorborg J, Tiedje N, Hattel J (2016) A micro-mechanical analysis of thermo-elastic properties and local residual stresses in ductile iron based on a new anisotropic model for the graphite nodules. Model Simul Mater Sci Eng 24:55012. https://doi.org/10.1088/0965-0393/24/5/055012
Andriollo T, Zhang Y, Fæster S, Thorborg J, Hattel J (2019) Impact of micro-scale residual stress on in-situ tensile testing of ductile cast iron: digital volume correlation vs. model with fully resolved microstructure vs. periodic unit cell. J Mech Phys Solids 125:714–735. https://doi.org/10.1016/j.jmps.2019.01.021
AZOM (2022) Graphite (C)—classifications, properties & applications. https://www.azom.com/
Behera A (2012) A novel material used in automotive industry: compacted graphite iron. Mater Res 1:1–4. https://doi.org/10.1680/emr.12.00002
Chen JK, Chen SF (2011) Thermal Conductivity of an in-situ metal matrix Composite—cast Iron. In: Cuppoletti J (ed) Metal, ceramic and polymeric composites for various uses. IntechOpen, Taiwan, pp 211–224
Dawson S (2008) Compacted graphite iron—a material solution for modern diesel engine cylinder blocks and heads. In: 68th WFC—World Foundry Congress, Chennai, India, pp 93–99
Di Cocco V, Iacoviello F, Cavallini M (2010) Damaging micromechanisms characterization of a ferritic ductile cast iron. Eng Fract Mech 77:2016–2023. https://doi.org/10.1016/j.engfracmech.2010.03.037
Di Cocco V, Iacoviello F, Rossi A, Cavallini M, Natali S (2013) Graphite nodules and fatigue crack propagation micromechanisms in a ferritic ductile cast iron. Fatigue Fract Eng Mater Struct 36:893–902. https://doi.org/10.1111/FFE.12056
Di Cocco V, Iacoviello F, Rossi A, Iacoviello D (2014) Macro and microscopical approach to the damaging micromechanisms analysis in a ferritic ductile cast iron. Theor Appl Fract Mech 69:26–33. https://doi.org/10.1016/j.tafmec.2013.11.003
Drago A, Pindera MJ (2007) Micro-macromechanical analysis of heterogeneous materials: macroscopically homogeneous vs periodic microstructures. Compos Sci Technol 67:1243–1263. https://doi.org/10.1016/j.compscitech.2006.02.031
Endo M, Yanase K (2014) Effects of small defects, matrix structures and loading conditions on the fatigue strength of ductile cast irons. Theor Appl Fract Mech 69:34–43. https://doi.org/10.1016/j.tafmec.2013.12.005
Frishmuth RE, McLaughlin PV (1976) Failure analysis of cast irons under general three-dimensional stress states. J Eng Mater Technol 98:69–75. https://doi.org/10.1115/1.3443339
Garoz D, Gilabert FA, Sevenois RDB, Spronk SWF, Van Paepegem W (2019) Consistent application of periodic boundary conditions in implicit and explicit finite element simulations of damage in composites. Compos Part B Eng 168:254–266. https://doi.org/10.1016/j.compositesb.2018.12.023
Greenstreet WL, Yahr GT, Valachovic RS (1973) The behavior of graphite under biaxial tension. Carbon N Y 11:43–57. https://doi.org/10.1016/0008-6223(73)90007-9
Hooputra H, Gese H, Dell H, Werner H (2004) A comprehensive failure model for crashworthiness simulation of aluminium extrusions. Int J Crashworthiness 9:449–464. https://doi.org/10.1533/ijcr.2004.0289
Hornbogen E (1985) Fracture toughness and fatigue crack growth of grey cast irons. J Mater Sci 20:3897–3905. https://doi.org/10.1007/BF00552378
Iacoviello F, Di Bartolomeo O, Cocco Di, Piacente V (2008) Damaging micromechanisms in ferritic–pearlitic ductile cast irons. Mater Sci Eng A 478:181–186. https://doi.org/10.1016/j.msea.2007.05.110
Kihlberg E, Norman V, Skoglund P, Schmidt P, Moverare J (2021) On the correlation between microstructural parameters and the thermo-mechanical fatigue performance of cast iron. Int J Fatigue 145:106112. https://doi.org/10.1016/j.ijfatigue.2020.106112
Kohout J (2001) A simple relation for deviation of grey and nodular cast irons from Hooke’s law. Mater Sci Eng A 313:16–23. https://doi.org/10.1016/S0921-5093(01)01145-5
Murakami S (2012) Continuum damage mechanics: a continuum mechanics approach to the analysis of damage and fracture. Springer, Berlin
Niu J, Huang C, Su R, Jiang G, Liu H (2021) Constitutive equation of compacted graphite iron (GJV450) at high temperature and high strain rate. Int J Adv Manuf Technol 113:2163–2174. https://doi.org/10.1007/s00170-020-06513-1
Norman V, Calmunger M (2019) On the micro- and macroscopic elastoplastic deformation behaviour of cast iron when subjected to cyclic loading. Int J Plast 115:200–215. https://doi.org/10.1016/j.ijplas.2018.11.019
Omairey SL, Dunning PD, Sriramula S (2019) Development of an ABAQUS plugin tool for periodic RVE homogenisation. Eng Comput 35:567–577. https://doi.org/10.1007/s00366-018-0616-4
Palkanoglou EN, Baxevanakis KP, Silberschmidt V V (2020) Interfacial debonding in compacted graphite iron: effect of thermal loading. Procedia Struct Integr 28:1286–1294. https://doi.org/10.1016/j.prostr.2020.11.110
Palkanoglou EN, Baxevanakis KP, Silberschmidt VV (2022) Thermal debonding of inclusions in compacted graphite iron: effect of matrix phases. Eng Fail Anal. https://doi.org/10.1016/j.engfailanal.2022.106476
Qiu Y, Pang JC, Li SX, Yang EN, Fu WQ et al (2016a) Influence of thermal exposure on microstructure evolution and tensile fracture behaviors of compacted graphite iron. Mater Sci Eng A 664:75–85. https://doi.org/10.1016/j.msea.2016.03.118
Qiu Y, Pang JC, Yang EN, Li SX, Zhang ZF (2016b) Transition of tensile strength and damaging mechanisms of compacted graphite iron with temperature. Mater Sci Eng A 677:290–301. https://doi.org/10.1016/j.msea.2016.09.035
Safaei M, Sheidaei A, Baniassadi M, Ahzi S, Mosavi Mashhadi M et al (2015) An interfacial debonding-induced damage model for graphite nanoplatelet polymer composites. Comput Mater Sci 96:191–199. https://doi.org/10.1016/j.commatsci.2014.08.036
Salomonsson K, Jarfors AEW (2018) Three-dimensional microstructural characterization of cast iron alloys for numerical analyses. Mater Sci Forum 925:427–435. https://doi.org/10.4028/www.scientific.net/MSF.925.427
Šamec B, Potrč I, Šraml M (2011) Low cycle fatigue of nodular cast iron used for railway brake discs. Eng Fail Anal 18:1424–1434. https://doi.org/10.1016/j.engfailanal.2011.04.002
Seldin EJ (1966) Stress-strain properties of polycrystalline graphites in tension and compression at room temperature. Carbon N Y 4:177–191. https://doi.org/10.1016/0008-6223(66)90079-0
Selin M (2010) Tensile and thermal properties in compacted graphite irons at elevated temperatures. Metall Mater Trans A Phys Metall Mater Sci 41:3100–3109. https://doi.org/10.1007/s11661-010-0385-1
Stefanescu DM (2019) The meritocratic ascendance of cast iron: from magic to virtual cast iron. Int J Met 13:726–752. https://doi.org/10.1007/s40962-019-00347-8
Tang C, Liu L, Yang Z, Tao D, Li J et al (2022) Surface evolution of vermicular cast iron in ultra-high temperature combustion with different single-pulsing duration. Eng Fail Anal 141:106679. https://doi.org/10.1016/j.engfailanal.2022.106679
Wang W, Jing T, Gao Y, Qiao G, Zhao X (2007) Properties of a gray cast iron with oriented graphite flakes. J Mater Process Technol 182:593–597. https://doi.org/10.1016/j.jmatprotec.2006.09.028
Wu Y, Li J, Yang Z, Guo Y, Ma Z et al (2019) Thermal conductivity analysis of compacted graphite cast iron after a creep test. Metall Mater Trans A 50:3697–3704. https://doi.org/10.1007/s11661-019-05278-x
Yang W, Pang J, Wang L, Wang S, Liu Y et al (2021) Tensile properties and damage mechanisms of compacted graphite iron based on microstructural simulation. Mater Sci Eng A 814:141244. https://doi.org/10.1016/j.msea.2021.141244
Zhan Y, Kaddouri W, Kanit T, Jiang Q, Liu L et al (2022) From unit inclusion cell to large Representative Volume Element: comparison of effective elastic properties. Eur J Mech A/Solids 92:104490. https://doi.org/10.1016/j.euromechsol.2021.104490
Zhang Y, Pang J, Shen R, Qiu Y, Li S et al (2018) Investigation on tensile deformation behavior of compacted graphite iron based on cohesive damage model. Mater Sci Eng A 713:260–268. https://doi.org/10.1016/j.msea.2017.12.077
Author information
Authors and Affiliations
Contributions
VVS and KPB conceptualized the study. MC implemented numerical simulations and prepared the draft together with graphical results. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Cao, M., Baxevanakis, K.P. & Silberschmidt, V.V. High-temperature behaviour and interfacial damage of CGI: 3D numerical modelling. Multiscale and Multidiscip. Model. Exp. and Des. (2023). https://doi.org/10.1007/s41939-023-00188-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s41939-023-00188-3