Next Article in Journal
Water-Ionic Liquid Binary Mixture Tailored Resorcinol-Formaldehyde Carbon Aerogels without Added Catalyst
Previous Article in Journal
Physical Characteristics of Coupled Plasma and Its Influence on Weld Formation in Hybrid Laser-Double-Arc Welding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Cement Composites Fracture Parameters Using Deformation Criterion

by
Marta Kosior-Kazberuk
1,
Andrzej Kazberuk
2 and
Anna Bernatowicz
1,*
1
Faculty of Civil and Environmental Engineering, Bialystok University of Technology, 15-351 Bialystok, Poland
2
Faculty of Mechanical Engineering, Bialystok University of Technology, 15-351 Bialystok, Poland
*
Author to whom correspondence should be addressed.
Materials 2019, 12(24), 4206; https://doi.org/10.3390/ma12244206
Submission received: 17 November 2019 / Revised: 9 December 2019 / Accepted: 11 December 2019 / Published: 14 December 2019
(This article belongs to the Section Construction and Building Materials)

Abstract

:
A simple deformation criterion based on Dugdale’s cohesive zone model is presented. The criterion can be used for both the experimental determination of the critical stress intensity factor, KIc, and the critical tip opening displacement, CTODc. It can also be applied for the evaluation of the load capacity of structural elements. The criterion is presented in explicit and compact form, which allows straightforward calculations to be performed for the estimation of KIc and CTODc values from the experimental data obtained from samples with a U-shaped notch, rounded with an arbitrary radius. Thanks to the simple form of the approximate relationship between the maximal load level and the dimensionless notch tip opening displacement, the reverse procedure was obtained, i.e., the estimation of the value of the maximal force loading the structural element as a function of the known critical stress intensity factor.

1. Introduction

The process of cracking an element of quasi-brittle material with a heterogeneous microstructure (e.g., cement composite) is associated with the formation and spreading of a major crack called a critical defect. The methods of fracture mechanics make it possible to describe the process of the destruction of such materials [1,2,3,4,5], as well as the degradation resulting from the impact of operational factors [6,7]. The parameters of material fracture mechanics, such as critical stress intensity factor, crack energy, and crack opening displacement, are determined experimentally using samples with cracks [8]. For metals, the procedure for determining these parameters has been standardized and is widely used. The cracking process begins with a fatigue-initiated crack in the notch tip.
In the case of cement composites and other brittle materials, the assessment of fracture mechanics parameters requires experimental tests using samples with a pre-formed primary crack or notch. Obtaining an initial crack of strictly defined geometric parameters is difficult. In standard quasi-brittle materials tests, an elongated U-shaped notch with a width of 1 to 10 mm is used, ignoring the effect of the rounding of its apex. Typically, a notch is produced at the stage of sample formation or by cutting with a diamond saw [9]. This approach is easy to implement, but rather narrow notches are obtained. Notches rounded with a greater radius are manufactured at the forming stage using reusable or disposable inserts. Notches of this kind are useful to analyze the influence of the notch geometry on the fracture process and especially to explain to what extent the geometry of the primary defect influences this. The crack propagation process should be analyzed based on the strength criteria of materials and on the basis of the criteria of fracture mechanics [10,11].
Most of the research on the influence of the shape of the notch on the fracture process has been focused on theoretical considerations rather than test methods and has been conducted using materials such as PMMA (poly(methyl methacrylate)) [10], metal alloys [11,12], and graphite [13]. The research results indicate a significant relationship between the notch rounding radius, the value of the damage force, and the cracking process evolution. Research conducted using concrete samples has mainly concerned the influence of the notch depth on fracture characteristics [6]. In the case of samples with notches where vertices are rounded with a large rounding radius (relative to the notch depth), the stress singularity disappears, and consequently, crack criteria based on linear elastic fracture mechanics cannot be used, even for brittle materials [12].
This work presents the application of a simple deformation criterion based on Dugdale’s cohesive model of the process zone [14], both for the experimental determination of the critical stress intensity factor and crack tip opening displacement as well as for the assessment of the loading capacity of the structural element.

2. Deformation Criterion

It is known that the limit equilibrium in the elastic–plastic body occurs when the opening displacement in the rounded notch tip is equal to its critical value:
δ I = δ c
where δc is the critical value of the crack tip opening displacement (CTODc) (for a sample in plane stress). It is assumed that δc is constant for a particular material and that it depends on neither the notch opening angle (it is the same for a crack and for a V-shaped notch) nor the notch tip rounding radius, ρ.
Assuming small inelastic deformations, i.e., when the length of the inelastic band is small in comparison to the crack length and the size of the body (e.g., [15,16]), the opening displacement at the vertex of the notch is related to the critical stress intensity factor for a sample with a crack [17,18]. In plane stress conditions this relationship takes the form
δ c = K c 2 E σ Y
where E is the modulus of elasticity and σY is the stress limit. Criterial Equation (1), for small inelastic deformations, can also be written in the following way [19]:
δ ˜ I ( K I ) 2 = K c 2
where δ ˜ I is the dimensionless opening displacement at the tip of a rounded U-shaped notch and KI is the stress intensity factor at the sharp vertex of the corresponding crack.
Equation (3) can be used to determine the unknown failure load level, by which cracking starts, and the related stress intensity factor, KI, while the other values are regarded as known. Additionally, when the load that causes cracking is known and consequently the value of KI is known, Equation (3) can be used for the experimental determination of fracture toughness, Kc, with the use of samples with rounded U-shaped notches.
The asymptotic solution to the elastic−plastic problem of an infinite rounded V-notch with a plastic band at the vertex was shown in [19,20]. The value of the dimensionless opening displacement at the notch vertex and the range of the inelastic zone are dependent on parameter γY, which is the dimensionless load level referenced to a material’s yield point [19,21], according to the following equation:
γ Y = 1 2 π ρ K I σ Y .
The dimensionless values of the length of the inelastic zone, l ˜ Y , arising from an infinite U-shaped notch and notch opening displacement, δ ˜ I , in function of parameter γY are shown in Figure 1 [19].
The graphs in Figure 1a,b, both for the length of the inelastic zone and the notch tip opening displacement, start from certain non-zero values of load level, because traction stress at the notch tip must achieve the yield point of the material:
( σ s ) m a x = K I 2 π ρ R I = σ Y ( γ Y ) m i n = 1 R I = 0.3341
where RI = 2.993 is the stress rounding factor [22,23]. For the U-shaped notch and parameter γY → ∞, the dimensionless inelastic band length approached l ˜ Y = π / 8 and the dimensionless notch opening displacement achieved δ ˜ I = 1 —the theoretical value of the crack opening displacement for a half-infinite crack with an inelastic band at the tip [24,25] (see also [26,27]).
The relationship between the dimensionless notch tip opening displacement, δ ˜ I , at the vertex of a U-shaped notch and the load level parameter, γY, can be approximated by
δ ˜ I = ( 1 a γ Y b ) 2 ,   1 R I < γ Y <
where parameter a = 0.12124 was fitted by the least-squares method and the value of parameter b,
b = ln a ln R I = 1.9247 ,
results from condition δ ˜ I ( 1 / R ) = 0 . The relative error of Equation (6) is less than 0.3% in the whole interval, 1 / R I < γ Y < .
Allowing a less precise estimation of δ ˜ I (with an error not exceeding 1%), the function δ ˜ I ( γ Y ) can be written down as
δ ˜ I = ( 1 ( R I γ Y ) 2 ) 2 ,   1 R I < γ Y < .
The estimated relationship between the dimensionless crack length, l ˜ Y , and parameter γY is described by the following equation:
l ˜ Y = π 8 c γ Y d ,   1 R I < γ Y <
where parameter c = 0.0088351 was fitted by the least-squares method and parameter d,
d = ln ( 8 a / π ) ln R I = 3.4611 ,
was calculated from the condition l ˜ Y ( 1 / R ) = 0 . The relative error of the estimation of l ˜ Y from Equation (9) is less than 0.5% in a whole interval, 1 / R I < γ Y < .
Assuming that a material’s parameters Kc and σY and the geometrical characteristics of a sample are known, the fracture load can be determined by the equation
K I δ ˜ I = K c .
The approximate relationship in Equation (10) between the critical displacement in the notch tip and the dimensionless load level parameter, γY, can be written as
δ I = 1 1 R I 2 γ Y 2
where
γ Y = K I σ Y 2 π ρ
By linking Equations (11) and (12) to Equation (10) the following quadratic equation is obtained:
K I [ 1 ( σ Y 2 π ρ R I K I ) 2 ] = K c K I 2 K c K I = 2 π ρ σ Y R I 2 = 0 .
By inserting the values of the stress intensity factor, which depend on the shape of the sample and the loading conditions,
K I = F I σ n π l ,
Equation (13) can be rewritten in the following form:
σ p 2 K c σ n F I π l 2 σ Y 2 ε F I 2 R I = 0 .
Thus, the positive root of Equation (15) is equal to
σ n σ Y = 1 2 F I [ 8 ε R I 2 + ( K c σ Y π l ) 2 + K c σ Y π l ]
where ε = ρ / l is the relative notch radius at the vertex. For ε = 0 (sharp crack), the minimal value of nominal stress is obtained:
σ n σ Y | m i n = 1 F I [ K c σ Y π l ] .

3. Experimental Tests

For the experimental verification of the proposed deformation criterion, the results presented in [28], which analyzed the impact of the U-shaped notch rounding radius on the parameters of fracture mechanics calculated in accordance with the two-parameter Shah criterion [2,29] (recommended by RILEM [30] and similar to JCI standard [31]), were used. The most important elements of the description and selected results are listed below.
The samples were made of mortar based on Portland cement. A mixture of natural aggregate with a grain size not exceeding 2 mm was used. The reduction of grain size was aimed at eliminating the possible impact of coarse aggregate on cracking characteristics. The water−cement ratio (w/c) of the mortar was 0.35. The selected material made it possible to make custom test elements with different notch geometry. The average compressive strength of concrete after 28 days of curing was 51.7 MPa, the average tensile strength at bending was 4.01 MPa, and the tangential modulus of elasticity determined during the measurements of fracture mechanics parameters reached the value of 22.3 GPa. In beam samples with dimensions of 100 × 100 × 100 mm, notches 30 mm deep (ratio a0/d = 0.30) were produced during formation. The notch rounding radii were 0.15, 0.35, 0.65, 1.35, 3.15, 8.0, 15.0, and 30.0 mm. The tests were performed after 28 days of curing the samples in water at 20 ± 2 °C.
U-notched beam elements were loaded with a concentrated force under three-point bending conditions. The loading rate was chosen so that the maximum value was reached within about 5 min. The dimensions and the loading scheme of the sample element are shown in Figure 2.
During the tests, the diagram of the loading force, P, as a function of the notch mouth opening displacement was recorded. The design of the test stand provided the conditions for stabilized damage of the samples. For the measurements, a universal hydraulic testing machine MTS 322 (manufactured by MTS Systems Corporation in Eden Prairie in the USA) was used, enabling fast deformation control and load regulation. The results were recorded automatically. A view of the test stand with the sample element is shown in Figure 3.

4. Results Analysis

To determine the critical stress intensity factor using the deformation criterion, the only significant value is the maximum value of the force destroying the sample. The initial fragments of the P-CMOD relationship (shown in Figure 4 for all notch rounding radii) were numerically smoothed and piecewise approximated with the polynomial of third degree to remove measuring noise and obtain a repeatable determination of the maximum force value.
The nominal stress value and the corresponding dimensionless load level factor were calculated for each sample:
σ n = 3 2 P S 2 t d 2 , γ Y = F I σ n σ Y 2 ε
where the dimensionless stress intensity factor value is equal to FI = 1.042 (for a0/d = 0.3 [32]).
The dimensionless notch tip opening displacement corresponding to the value of γY was calculated according to Equation (6):
δ ˜ I = 1 0.12124 γ Y 1.9247 .
The resulting values of the critical stress intensity factor were calculated using the relationship in Equation (3) in the form of
K c = 1.0426 σ n δ I π a 0 .
The obtained average values of the critical stress intensity factor for each analyzed relative notch rounding radius are shown in Figure 5. The value of δc = CMODc was calculated from Equation (2).
By using the obtained average critical value of the stress intensity factor, Kc = 0.46 MPa√m, the criterion values of nominal stresses destroying the tested element were calculated from Equation (16).
In dimensionless form, the criterion values are presented in Figure 6 with the experimental results.

5. Conclusions

The paper presents the application of the deformation criterion for determining the parameters of fracture mechanics of quasi-brittle materials using U-notched samples. Due to the ease of making this type of edge notch or hole with parallel edges and rounded vertices in samples at the forming stage, the criterion is particularly well suited for determining the parameters of concrete and cement composites.
This criterion has been presented in short form, which allows the simple calculation of the critical stress intensity factor, Kc, and of the critical tip opening displacement, CTODc, for any sample geometry with a rounded U-shaped notch and an arbitrary relative rounding radius. Due to the simple form of the approximating formulas used to calculate the dependence of the load level on the dimensionless notch opening displacement, formulas have also been obtained that enable the inverse procedure, i.e., the estimation of the destructive force of the structural element based on the known stress intensity factor.

Author Contributions

Conceptualization, M.K.-K., A.K., and A.B.; methodology, M.K.-K., A.K., and A.B.; software, M.K.-K., A.K., and A.B.; validation, M.K.-K., A.K., and A.B.; formal analysis, M.K.-K., A.K., and A.B.; investigation, M.K.-K., A.K., and A.B.; resources, M.K.-K., A.K., and A.B.; data preparation, M.K.-K., A.K., and A.B.; writing—original draft preparation, M.K.-K., A.K., and A.B.; writing—review and editing, M.K.-K., A.K., and A.B.; visualization, M.K.-K., A.K., and A.B.; project administration, M.K.-K., A.K., and A.B.; funding acquisition, M.K.-K. and A.B.

Funding

This research was funded by the Ministry of Science and Higher Education, grant numbers S/WBiIŚ/2/2017 and W/WBiIŚ/10/2019.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bažant, Z.P. Concrete fracture models: Testing and practice. Eng. Fract. Mech. 2002, 69, 165–205. [Google Scholar] [CrossRef]
  2. Shah, S.P.; Swartz, S.E.; Ouyang, C. Fracture Mechanics of Concrete: Applications of Fracture Mechanics to Concrete, Rock and other Quasi-Brittle Materials; John Wiley & Sons, Inc.: New York, NY, USA, 1995. [Google Scholar]
  3. Bažant, Z.P.; Gettu, R.; Jirásek, M.; Barr, B.I.G.; Carol, I.; Carpinteri, A.; Elices, M.; Huet, C.; Mihashi, H.; Nemati, K.M.; et al. RILEM TC QFS Quasibrittle fracture scaling and size effect- Final report. Mater. Struct. 2004, 37, 547–568. [Google Scholar]
  4. Karihaloo, B.L.; Nallathambi, P. An improved effective crack model for the determination of fracture toughness of concrete. Cem. Concr. Res. 1989, 19, 603–610. [Google Scholar] [CrossRef]
  5. Spagnoli, A.; Carpinteri, A.; Ferretti, D.; Vantadori, S. An experimental investigation on the quasi-brittle fracture of marble rocks. Fatigue Fract. Eng. Mater. Struct. 2016, 39, 956–968. [Google Scholar] [CrossRef]
  6. Jenq, Y.S.; Shah, S.P. Features of mechanics of quasi-brittle crack propagation in concrete. Int. J. Fract. 1991, 51, 103–120. [Google Scholar]
  7. Astm, D-7313. Standard Test Method for Determination Fracture Energy of Asphalt-Aggregate Mixtures Using the Disc-Shaped Compact Tension Geometry; ASTM International: West Conshohocken, PA, USA, 2007. [Google Scholar]
  8. ASTM. Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness KIC of Metallic Materials; ASTM International: West Conshohocken, PA, USA, 2009. [Google Scholar]
  9. Kosior-Kazberuk, M.; Kazberuk, A. Determination of the fracture toughness parameters of quasi-brittle materials using cylindrical samples. Ceram. Mater. 2010, 62, 244–248. [Google Scholar]
  10. Berto, F.; Campagnolo, A.; Elices, M.; Lazzarin, P. A synthesis of Polymethylmethacrylate data from U-notched specimens and V-notches with end holes by means of local energy. Mater. Des. 2013, 49, 826–833. [Google Scholar] [CrossRef]
  11. Berto, F.; Barati, E. Fracture assessment of U-notches under three point bending by means of local energy density. Mater. Des. 2011, 32, 822–830. [Google Scholar] [CrossRef]
  12. Gomez, F.J.; Elices, M.; Berto, F.; Lazzarin, P. A generalized notch stress intensity factor for U-notched components loaded under mixed mode. Eng. Fract. Mech. 2008, 75, 4819–4833. [Google Scholar] [CrossRef]
  13. Torabi, A.R. Sudden fracture from U-notches in fine-grained isostatic graphite under mixed mode I/II loading. Int. J. Fract. 2013, 181, 309–316. [Google Scholar] [CrossRef]
  14. Dugdale, D. Yielding of steel sheets containing slits. J. Mech. Phys. Solids 1960, 8, 100–104. [Google Scholar] [CrossRef]
  15. Rice, J.R. Limitations to the small scale yielding approximation for crack tip plasticity. J. Mech. Phys. Solids 1974, 22, 17–26. [Google Scholar] [CrossRef]
  16. Rice, J.R. The localization of plastic deformation. In Theoretical Applied Mechanics; North-Holland Publishing Company: Amsterdam, The Netherlands, 1976; Volume 1, pp. 207–220. [Google Scholar]
  17. Panasyuk, V.V.; Andreykiv, A.E.; Parton, V.Z. The Bases of Materials Fracture Mechanics; Naukova Dumka: Kiev, Ukraine, 1988. (In Russian) [Google Scholar]
  18. Parton, V.Z. Fracture Mechanics: From Theory to Practice; CRC Press: Boca Raton, FL, USA, 1992. [Google Scholar]
  19. Savruk, M.P.; Kazberuk, A. Two-dimensional fracture mechanics problems for solids with sharp and rounded V-notches. Int. J. Fract. 2010, 161, 79–95. [Google Scholar] [CrossRef]
  20. Savruk, M.P.; Kazberuk, A. Stress concentration near sharp and rounded V-notches in orthotropic and quasi-orthotropic bodies. Theor. Appl. Fract. Mech. 2016, 84, 166–176. [Google Scholar] [CrossRef]
  21. Savruk, M.P.; Kazberuk, A. Problems of fracture mechanics of solid bodies with V-shaped notches. Mater. Sci. 2009, 45, 162–180. [Google Scholar] [CrossRef]
  22. Savruk, M.P.; Kazberuk, A. Relationship between the stress intensity and stress concentration factors for sharp and rounded notches. Mater. Sci. 2006, 42, 725–738. [Google Scholar] [CrossRef]
  23. Savruk, M.P.; Kazberuk, A. Stress Concentration at Notches; Springer: Cham, Switzerland, 2017. [Google Scholar]
  24. Rice, J.R. Plastic yielding at a crack tip. In Proceedings of the 1st International Conference on Fracture, Sendai, Japan, 17 October 2012; pp. 283–308. [Google Scholar]
  25. Cherepanov, G.P. Mechanics of Brittle Fracture; McGraw-Hill: New York, NY, USA, 1979. [Google Scholar]
  26. Panasyuk, V.V.; Savruk, M.P. Model for plasticity bands in elastoplastic failure mechanics. Mater. Sci. 1992, 28, 41–57. [Google Scholar] [CrossRef]
  27. Panasyuk, V.V.; Yarema, S.Y. On the Origin of the k—Model and the Model of Plastic Strips. Mater. Sci. 2001, 37, 346–353. [Google Scholar] [CrossRef]
  28. Kosior-Kazberuk, M.; Kazberuk, A. Application of deformation criterion to assess fracture properties of concrete. In Proceedings of the ICF 2017—14th International Conference on Fracture, Fairfax, VA, USA, 18–23 June 2017; pp. 348–349. [Google Scholar]
  29. Shah, S.; Carpinteri, A. Fracture Mechanics Test Methods for Concrete; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  30. Shah, S.P. Determination of fracture parameters (KIC and CTODC) of plain concrete using three-point bend tests, Materials and Structures. Kluwer Acad. Publ. 1990, 23, 457–460. [Google Scholar]
  31. Kitsutaka, Y.; Kanakubo, T. Outline of JCI standard (JCI-S-001-003). Concr. J. 2006, 44, 10–15. [Google Scholar] [CrossRef] [Green Version]
  32. Savruk, M.P. Stress Intensity Factors in Cracked Bodies; Naukova Dumka: Kiev, Ukraine, 1988. (In Russian) [Google Scholar]
Figure 1. Relationship between the dimensionless length of an inelastic zone l ˜ Y (a) and the dimensionless opening displacement δ ˜ Y (b) at the vertex of a U-notch and the dimensionless load level parameter γY.
Figure 1. Relationship between the dimensionless length of an inelastic zone l ˜ Y (a) and the dimensionless opening displacement δ ˜ Y (b) at the vertex of a U-notch and the dimensionless load level parameter γY.
Materials 12 04206 g001
Figure 2. The loading scheme and dimensions of a sample with a U-shaped notch.
Figure 2. The loading scheme and dimensions of a sample with a U-shaped notch.
Materials 12 04206 g002
Figure 3. An image of the test stand.
Figure 3. An image of the test stand.
Materials 12 04206 g003
Figure 4. The initial phase of the P-CMOD graph for single-edge, U-notched samples with various relative notch rounding radii. ε = ρ/a0.
Figure 4. The initial phase of the P-CMOD graph for single-edge, U-notched samples with various relative notch rounding radii. ε = ρ/a0.
Materials 12 04206 g004
Figure 5. Average values of critical stress intensity factor, Kc, as a function of relative U-notch rounding radius.
Figure 5. Average values of critical stress intensity factor, Kc, as a function of relative U-notch rounding radius.
Materials 12 04206 g005
Figure 6. The critical values of nominal stresses for beams with U-shaped notches as a function of relative U-notch rounding radius.
Figure 6. The critical values of nominal stresses for beams with U-shaped notches as a function of relative U-notch rounding radius.
Materials 12 04206 g006

Share and Cite

MDPI and ACS Style

Kosior-Kazberuk, M.; Kazberuk, A.; Bernatowicz, A. Estimation of Cement Composites Fracture Parameters Using Deformation Criterion. Materials 2019, 12, 4206. https://doi.org/10.3390/ma12244206

AMA Style

Kosior-Kazberuk M, Kazberuk A, Bernatowicz A. Estimation of Cement Composites Fracture Parameters Using Deformation Criterion. Materials. 2019; 12(24):4206. https://doi.org/10.3390/ma12244206

Chicago/Turabian Style

Kosior-Kazberuk, Marta, Andrzej Kazberuk, and Anna Bernatowicz. 2019. "Estimation of Cement Composites Fracture Parameters Using Deformation Criterion" Materials 12, no. 24: 4206. https://doi.org/10.3390/ma12244206

APA Style

Kosior-Kazberuk, M., Kazberuk, A., & Bernatowicz, A. (2019). Estimation of Cement Composites Fracture Parameters Using Deformation Criterion. Materials, 12(24), 4206. https://doi.org/10.3390/ma12244206

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop