Loading [MathJax]/jax/output/SVG/jax.js
Research article

Hypertemperature effects in heterogeneous media and thermal flux at small-length scales

  • We propose an enriched microscopic heat conduction model that can account for size effects in heterogeneous media. Benefiting from physically relevant scaling arguments, we improve the regularity of the corrector in the classical problem of periodic homogenization of linear elliptic equations in the three-dimensional setting and, while doing so, we clarify the intimate role that correctors play in measuring the difference between the heterogeneous solution (microscopic) and the homogenized solution (macroscopic). Moreover, if the data are of form f=divF with FL3(Ω,R3), then we recover the classical corrector convergence theorem.

    Citation: Grigor Nika, Adrian Muntean. Hypertemperature effects in heterogeneous media and thermal flux at small-length scales[J]. Networks and Heterogeneous Media, 2023, 18(3): 1207-1225. doi: 10.3934/nhm.2023052

    Related Papers:

    [1] Ayman Elsharkawy, Hoda Elsayied, Abdelrhman Tawfiq, Fatimah Alghamdi . Geometric analysis of the pseudo-projective curvature tensor in doubly and twisted warped product manifolds. AIMS Mathematics, 2025, 10(1): 56-71. doi: 10.3934/math.2025004
    [2] Zihan Yue, Wei Jiang, Boying Wu, Biao Zhang . A meshless method based on the Laplace transform for multi-term time-space fractional diffusion equation. AIMS Mathematics, 2024, 9(3): 7040-7062. doi: 10.3934/math.2024343
    [3] Zhe Su, Yiying Tong, Guo-Wei Wei . Persistent de Rham-Hodge Laplacians in Eulerian representation for manifold topological learning. AIMS Mathematics, 2024, 9(10): 27438-27470. doi: 10.3934/math.20241333
    [4] Rabia Cakan Akpınar, Esen Kemer Kansu . Metallic deformation on para-Sasaki-like para-Norden manifold. AIMS Mathematics, 2024, 9(7): 19125-19136. doi: 10.3934/math.2024932
    [5] Chao Wang, Fajie Wang, Yanpeng Gong . Analysis of 2D heat conduction in nonlinear functionally graded materials using a local semi-analytical meshless method. AIMS Mathematics, 2021, 6(11): 12599-12618. doi: 10.3934/math.2021726
    [6] H. Thameem Basha, R. Sivaraj, A. Subramanyam Reddy, Ali J. Chamkha, H. M. Baskonus . A numerical study of the ferromagnetic flow of Carreau nanofluid over a wedge, plate and stagnation point with a magnetic dipole. AIMS Mathematics, 2020, 5(5): 4197-4219. doi: 10.3934/math.2020268
    [7] Wei Sun, Yan Wang, Shiyong Li . An optimal resource allocation scheme for virtual machine placement of deploying enterprise applications into the cloud. AIMS Mathematics, 2020, 5(4): 3966-3989. doi: 10.3934/math.2020256
    [8] Abdon Atangana, Ali Akgül . Analysis of a derivative with two variable orders. AIMS Mathematics, 2022, 7(5): 7274-7293. doi: 10.3934/math.2022406
    [9] Fahad Sikander, Tanveer Fatima, Sharief Deshmukh, Ayman Elsharkawy . Curvature analysis of concircular trajectories in doubly warped product manifolds. AIMS Mathematics, 2024, 9(8): 21940-21951. doi: 10.3934/math.20241066
    [10] Jun Ma, Junjie Li, Jiachen Sun . A novel adaptive safe semi-supervised learning framework for pattern extraction and classification. AIMS Mathematics, 2024, 9(11): 31444-31469. doi: 10.3934/math.20241514
  • We propose an enriched microscopic heat conduction model that can account for size effects in heterogeneous media. Benefiting from physically relevant scaling arguments, we improve the regularity of the corrector in the classical problem of periodic homogenization of linear elliptic equations in the three-dimensional setting and, while doing so, we clarify the intimate role that correctors play in measuring the difference between the heterogeneous solution (microscopic) and the homogenized solution (macroscopic). Moreover, if the data are of form f=divF with FL3(Ω,R3), then we recover the classical corrector convergence theorem.



    Abbreviations: Wf : Total strain energy attributed to fracture; GIc : Surface release energy or critical mode I fracture toughness; JIc : Critical mode I fracture toughness; KIc : Fracture toughness; PQ : Applied load; we : Essential work of fracture in the elastic zone; βp : Geometric shape factors related to the plastic zone during tearing after necking; βy : Geometric shape factors related to the plastic zone during ligament yielding; : Pre-crack length; L: Ligament length; t: Thickness of the specimen; w: Sample width; wp : Non-essential work of the fracture; Wp : Work of the tearing and necking of the plastic zone; Wpp : Relative plastic energy in tearing and necking; Wpy: Relative energy in plastic and yielding of ligament length; Wy : Elastic energy of the elastic and yielding ligament length; wy : Essential work of fracture in the elastic zone; β: Plasticity shape factor; βwp: Slope of linear fitting regression; δ˳: 0.2% offset displacement; f(aw) : Correction factor

    In recent years, aluminium alloys have been sought in the ground vehicle industry because of their special properties, namely, a combination of strength, fatigue resistance, formability, and corrosion resistance. This helped classify these as special metals. Magnesium yields 5754 H-111 aluminium alloys. Their strength originates from the solid-solution strengthening. Chromium, manganese, and zirconium are added to control the grain and subgrain structures and increase the strength of the alloys. The H111 designation for the condition differs from the O condition in that semi-finished products supplied in this condition are rolled after annealing to improve their dimensional properties such as flatness [1,2].

    It is important to understand the fracture mechanisms to prevent in-service failures owing to design and usage problems. Fracture mechanics are dependent on crack-like material defects regardless of their origin. It is a gradual as well as temporally and spatially inhomogeneous process that begins with the initiation of a crack, propagates, and eventually causes the fracture or coalescence of cracks. This occurs when stresses do not exceed the yield stress. The critical stress concentration for fracture is still in the range of linear elastic fracture mechanics and is an inherent material property, KIC [3]. Toughness is an important parameter for measuring the capability of a material to absorb energy before failure. This energy is considered to represent a large plastic deformation before separation, which is the main difference between brittle and ductile fractures.

    The essential work of fracture (EWF) approach [4] is effective for dividing the total energy expended in the FPZ into two components: the EWF we and non-EWF wp. The crack propagation resistance is represented by the EWF, we. The critical fracture energy (GIC) differs from the J-integral and represents the actual crack-propagation resistance after the initiation of a fracture.

    Shinde et al. [5] used a modified one-sided notched specimen design to measure thin 6061-T6 alloy aluminium plates. The modified design reduced buckling of the specimen. It was observed that the fracture toughness of the thin 6061-T6 aluminium alloy plates was higher than the plane load toughness. However, a linear finite element model was used to calculate the stress field distribution.

    Tippeswamy [6] conducted an experiment using aluminium 6082-T6 in accordance with ASTM standards. The specimens were tested by varying the notch length within a specified constraint. The experimental results showed that the notch length significantly affected the size of the plastic zone at the end of the cut. Furthermore, the rates of fracture formation and propagation varied depending on the notch length of each specimen.

    Pardoen et al. [7] investigated the influence of sheet thickness on the fracture strength of a thin sheet of the aluminium alloy 6082-o with a thickness of 1–6 mm for the final separation using the EWF and J-integral methods in a tensile test (double-edge notched tensile specimen (DENT)). Their study aimed to determine the relationship between these two strategies. The EWF is a measure of the fracture initiation resistance. Go when there are few variations in GC.

    Using numerical simulations, Moiseenko et al. [8] analysed hybrid discrete-continuous cellular automata that enable the calculation of the local moments of force and the simulation of the dissipation of torsional energy. This resulted in the development of new defect structures. The simulations included the calculation of thermal expansion with the simultaneous build-up of thermal stresses and initiation of microrotation. This allowed for the assessment of local entropy and tracking of crystal defect evolution from the onset to storage. The mechanical behaviours of materials subjected to thermally-induced twinning or phase transitions were simulated, and the algorithms were explained.

    Balokhonov et al. [9] presented a numerical analysis of the microstructural influences on the deformation and fracture of friction stir welds in aluminium alloys. The mechanical response of individual grains was simulated using an elastic–plastic formulation of the problem that considers isotropic strain hardening including the Hall–Petch effect, and a fracture model that permits crack initiation and growth in the regions of the maximum equivalent plastic strain. The strength of welded materials has been shown to influence the plastic strain and fracture localisation.

    Balokhonov et al. [10] investigated the localisation of plastic deformation and fractures in a porous coated material. The numerical simulations were performed using the finite-difference method. The presence of localised tensile zones around the pores and at the interface between the coating and substrate (both in tension and compression of the coated material) was demonstrated to be the cause of the unique deformation and fracture characteristics of the composite studied. The relationship between the crack propagation in the coating and the inhomogeneous plastic flow in the steel substrate was investigated.

    A thin aluminium strip was measured by Abdellah [11] using the EWF. A simple finite-element numerical model was constructed to predict the EWF parameters. Although the sample was measured in millimetres, the experimental results and proposed model agreed considerably. The concept of different energy zones around the fracture tip was first introduced by Broberg [12]. The fracture process zone (FPZ) and outer plastic zone are two subdivisions of the crack tip generated by Broberg. Two new surfaces were generated using energy from the FPZ. The FPZ is not affected by the loading conditions or stress at the fracture tip. The fracture length, load, shape, and stress state affect the outer plastic zone. Cotterell and Reddel first proposed the central thesis on the fracture concept in 1976 [4]. Cotterell and Reddel used the Broberg notion of energy expended in an independent autonomous end zone (sheet). In linear elastic fracture mechanics, the plastic zone around the crack tip is assumed to be marginal.

    Korsunsky and Kyungmok [13] investigated the feasibility of quantifying the crucial work in ductile cracking using single tensile tests on dog-bone specimens without notches. They used different types of heat-treatable aluminium alloys. The obtained results were compared with those of a typical DENT. The authors concluded that a laser scanner could be used.

    Most conventional methods [14,15,16,17] for determining fracture parameters such as GIC or the stress intensity factor (SIF), KIC, were developed considering isotropic materials when linear elastic fracture mechanics (LEFM) was originally developed.

    Masuda et al. [17] used A7075-T6 compact tensile (CT) specimens of various thicknesses (1–21 mm) in fatigue crack growth studies. As the thickness decreases, the radius of curvature of the leading edge of the fatigue crack increases. The extent of the lateral plastic contraction (indentation depth d) at the crack tip during fatigue loading was estimated quantitatively using the three-dimensional elastoplastic finite element method. The results of the experiment are as follows. The rate of fatigue crack growth da/dN at a constant K increased with an increase in t from 1 to 11 mm in the range where K ≥ 5 MPam. Between t = 11 and 21 mm, da/dN remained constant. Meanwhile, in the range where K is less than 5 MPam. A J-R curve (or D–R curve) is a plot of the resistance to stable crack extension measured as J (or R-d) plotted against ductile crack extension (generally considered as Dap, the measured physical crack extension). However, in many cases, only estimates of this crack extension are available. ASTM E1820 [18] was developed for measuring the elastic–plastic initiation toughness, JIC, and J-R curves or the corresponding dIc and d-R curves under plane-strain conditions. However, in [19], it was developed for evaluating the CTOD at the onset of cleavage, dc. In thin-walled materials with low constraints, the CTOA parameter is used to describe stable crack extension.

    Complex crack configurations in finite plates pose a significant challenge for SIF calculations. Byskov [20] proposed a thorough numerical FEM to solve the above problems by focusing on specific cracked elements, Here, the stiffness matrix is connected to the crack elements. To determine the K-R curves according to ASTM E1820, Lu and Wang [21] tested C(T) and M(T) specimens over a range of thicknesses (plane stress and plane strain). A novel method was presented for determining the CTOA using the K-R test data. Considering the effective performance of the ASTM E1820 test, Di et al. [22] recommended the use of a modified C(T) specimen to evaluate the CTOA of X80 pipeline steel. Finite element analysis (FEA) was used to predict the failure of the C(T) and M(T) specimens based on the CTOA fracture criteria.

    As mentioned in the previous paragraph, many researchers have studied the fracture toughness of the aluminium alloy 5754-H111 (also known as Al-Mg) under ductile damage using linear elastic fracture mechanics. Therefore, this study focuses on using other indirect methods with a simple technique to complete the description of the plastic behaviour of materials. The objectives of this work are as follows: (1) to establish that the EWF method is reasonable for measuring the fracture toughness in ductile damage energi JIC, (2) to compare the results of standard compact tensile specimens with EWF, and (3) to investigate the fracture topography and failure modes in ductile damage.

    The paper is structured as follows: In the first section, the EWF concepts are outlined. Then, the material and method are explained. In the third section, the results of the EWF and CT are illustrated and discussed. Finally, the comparison is limited to the applicability of the method.

    The EWF method can be used to evaluate the fracture toughness of thin layers of materials. It has gained considerable popularity in recent years [5]. According to the EWF method, the total energy (Wf) required to fracture a notched specimen is divided into two parts: the essential work (We) used to generate new surfaces in the so-called fracture process zone and the non-essential work (Wp) used to plastically deform the area around the process zone. Therefore the specific work, Wf, can be expressed as the combination of the following two terms [11] (see Eqs 1–3):

    Wf=δ0pdδ (1)
    Wf=We+Wp (2)
    Wf=We+βWp (3)

    Figure 1 shows the EWF technique used to break the surfaces of typical DENT specimens. When the two plastic zones are formed at the crack touch, the ligament L yields completely under the maximum force applied to the specimen [23]. The curve between the load and displacement when the ligament fails completely owing to ductile cracks is shown in Figure 2 [11,24].

    Figure 1.  Double-edge notched specimen with a plastic zone [24,25].
    Figure 2.  EWF data (a) load–displacement curve and (b) EWF fitting [25].

    Here, Wp is the plastic deformation zone behind the fracture process zone, we is the instability of the crack tip, and we is the surface release work in the crack process zone. Furthermore, the failure displacement δ. For a given specimen thickness, the surface release depends on the ligament length L. The volume energy (L 2t) is proportional to the volume. It describes the plastic work (Wp). Equation 2 is divided by the ligament area, Lt, to obtain the following expression for energy:

    wf=WfLt=we+βwpL (4)

    where β is the plastic work per unit volume of the plastic deformation zone in front of the crack tip as wp, the specific non-essential work of the fracture, and, the shape factor of the plastic deformation. The surface release energy is also necessary for the formation of the cracked surface. The relationship between wf and the ligament length L is shown in Eq 4. It is a linear regression. The range of the ligament lengths for the effective EWF tests was the plastic work per unit volume of the plastic deformation zone before the crack tip. In addition, we considers the surface release energy required to initiate the crack surface formation. The relationship of wf shown in Eq 4 is a linear regression relating it to the ligament length L. According to the "rule of thumb" of Cotterell and Reddel [4], the range of ligament lengths for the valid EWF tests is shown in Eq 5:

    Lmin=(35)t<L<Lmax=W3 (5)

    The positive intercept at L = 0 corresponds to a specific EWF. The slope of the regression line was determined by linearly fitting the data with the non-EWF, wp (see Figure 2b). After a load is applied and the ligament yields completely, Eq 2 for a DENT specimen could be rewritten as follows:

    Wf=Wy+Wpp (6)

    where Wy is the mechanical energy in the elastic zone and Wpp is the plastic energy used for constriction or subsequent tearing in the plastic zone (see Figure 2a). Using Eq 4, we can be divided into two zones: the elastic zone of EWF Wey associated with crack initiation and plastic zone of EWF Wepp associated with cracking prior to necking before the crack tip (see Eq 7):

    we=Wey+Wepp (7)

    Eq 8 expresses the divisions of the slope:

    βWp=βywpy+βpwpp (8)

    where βy and βp are the geometric slope parameters associated with the plastic zone during ligament yielding and tearing after necking, respectively.

    The commercial aluminium 5754-H111 alloy specimens were obtained from Egypt Alum. Co. The chemical composition of the A5754 -T alloy is presented in Table 1.

    Table 1.  Chemical composition of aluminium 5754-H111 [26].
    Zn Cu Mn Si Fe Mg Cr Ti Al
    0.2 0.1 0.5 0.4 0.4 2.6–3.6 0.3 0.15 Bal.

     | Show Table
    DownLoad: CSV

    These substances significantly affect aluminium alloy 5754. These provide a remarkable corrosion resistance, particularly in the presence of seawater and industrially polluted atmospheres. It is a medium-strength alloy. In accordance with the symbol H111, the alloy was subjected to formative work-hardening. The aluminium alloy 5754-H111 is characterised by a high fatigue strength, good machinability, and good cold formability. The alloy belongs to the high-strength 5xxx series. Therefore, 5754-H111 is highly suitable for use in flooring, shipbuilding, and structures addressing chemicals and nuclear energy. It is also widely used in structural components and interior trims of automobiles [27,28]. The chemical components were determined using XRD (see Figure 3). It showed that the maximum dislocation density was achieved by the following factors:

    Figure 3.  XRF for aluminium sheet A 5754.

    Peak Angle: 78.24°, FWHM: 0.3380°, Scherrer Crystallite Size: 53.86 nm

    Peak Angle: 65.1°, FWHM: 0.2860°, Scherrer Crystallite Size: 68.53 nm

    Peak Angle: 44.74°, FWHM: 0.2460°, Scherrer Crystallite Size: 88.49 nm

    Start: 5°, End: 80°, Step: 0.02°, Speed: 2.5°/min, Time/Step: 7.999999E-03°, Wavelength: 1.54056 (Cu).

    It was observed that this aluminium alloy has a larger amount of Mg (2.6–3.6). This makes it more brittle and vulnerable to environmental effects. Meanwhile, Mn and Si have relatively higher values and can precipitate during cracking action.

    A standard tensile test was performed according to ASTM E399-81 [29] to determine the following mechanical properties: tensile strength, elongation strength, Young's modulus, and percentage elongation. A standard rectangular tensile specimen with dimensions 90 × 15 × 10 mm was used (Figure 4). The tests were conducted using a computer-controlled electromechanical universal testing machine (model WDW-100-Jinan Victory Instrument Co. Ltd., China) [30] with a loading capacity of 100 kN and controlled speed of 2 mm/min. To understand the ductile damage, the topography was investigated. The aim was to conduct these tests via SEM because the tensile test is a basic and standard test for investigating failure and damage.

    Figure 4.  Tensile test sample.

    The specimens were subjected to an ASTM Standard D 5045 compliant fracture toughness test [31]. To obtain accurate results for the fracture toughness and plane strain fracture toughness KIC (which indicates the resistance of a material to fracture), the crack resistance measurement should be stopped at the point of degradation. According to the measurements provided in ASTM D5045 [31], two compact tensile specimens were fabricated using a milling machine (see Figure 5). The specimens were clamped between two sacrificial plates made of an identical material to prevent damage. Load holes were drilled using carbide drills. The initial crack was generated using a diamond blade with a thickness of 1 mm. The CT geometry is recommended because it enables planar loading with smaller specimen sizes than the other designs. The widths (W), thicknesses (t), and crack lengths (a) of the specimens are shown in Figure 4. The crack length a is 0.45–0.55 times W. The nominal value of the W / B ratio is two.

    0.45<a/W<0.55 (9)
    Figure 5.  CT test specimen geometry [33].

    Two thicknesses (t) of 1.8 mm and 5 mm were used. Both the specimens had a width of 60 mm and length of 65 mm, according to [32] (see Figure 5). These thicknesses were selected for porosities to study the effect of the Mg reaction during the test when oxygen penetrated the crack surface, from an available thin plate with a thickness of 1.8 mm and a larger thickness of 5 mm. It is also established that the fracture toughness depends on the thickness according to LEFM. Therefore, it would be effective to review this concept for the EWF technique [30]. The load corresponding to an apparent increase in the crack extension of 2.5% was determined by a certain deviation from the linear part of the record. KIC was calculated from this load using equations established based on elastic stress analyses of specimens of the type described in the test methods. The generation of a sharp crack at the crack tip in a specimen sufficiently large to exhibit linear elastic behaviour is necessary to validate the evaluation of the KIC value determined using these test methods.

    A DENT was tested using a universal testing machine at a transverse speed of 2 mm/min at room temperature. The specimens had a length of five ligaments, as shown in Figure 1. The ligaments had lengths of 4, 6, 10, 12, and 14 mm. The specimen was cut with a CNC milling machine. The cracks were cut with a sharp blade of 1 mm thickness according to [34,35,36,37]. The load was applied on both sides of the specimen to complement failure, and the load and displacement were recorded. The number of specimens for each ligament was three. After complete failure attained the mechanical energies, Wy and Wpp were calculated using Eqs 10 and 11 by integrating the resulting load–displacement curves as follows:

    Wy=δ00pdδ (10)
    Wpp=δδ0pdδ (11)

    where δ0 and δ are the displacements at 0.02% offset and failure, respectively. p is the applied load. The obtained total energy Wf (measured using Eq 6) was plotted against the ligament length L. The displacement at fracture (δ) could be plotted against L according to [38].

    Figure 6 shows the stress–strain relationship obtained by the tensile test of the aluminium alloy 5754. It is observed that the stress was uniform, yield strength was 153.9, tensile strength was 265 MPa, and Young's modulus was 68 GPa. The percent elongation was 1.2, and the ductile fracture behaviour was characterised by a large plastic zone. The fracture investigation was performed by SEM image analysis. It was observed that in the case of the aluminium alloy before fracture, the alumina was distributed uniformly over the aluminium matrix, and no voids were observed (see Figure 7a). Although the case occurred after the fracture and deformation, it was observed that the deep pits were distributed over the matrix. This was owing to void coalescence. Moreover, MgO formed during the test. As the oxygen insert induced by the microcrach, these hard particles of MgO participated at the grain boundary of the lighter lines (see Figure 7b). The formation of these oxides rendered the alloy brittle, thereby resulting in microcracks. It was also observed that the higher strength was owing to the strong bonding between the alloying elements [39,40]. The results are summarised in Table 2.

    Figure 6.  Stress and strain relationship for aluminium alloy.
    Figure 7.  SEM image of fracture surface in tension.
    Table 2.  Mechanical properties of Al 5754-H111 alloy.
    Yield strength
    (MPa)
    Ultimate tensile strength
    (MPa)
    Young's modulus
    (GPa)
    % elongation
    153.9 265 68 1.2

     | Show Table
    DownLoad: CSV

    Figure 8a, b shows the force–displacement curves obtained from the EWF measurement of the 5 mm and 1.8 mm aluminium DENT specimens at room temperature. The curves show a linear increase up to a certain point where an oscillation occurs indicating the flow range. These continue to increase until these attain the peak value (full flow). Then, ductile cracking and eventual failure occur. The load-displacement curves as a function of ligament length show similar geometries. The total energy Wf stored in the cracked specimen is represented by the area under the force-displacement curve. This area is measured using Eq 1. The internal work increases with an increase in the ligament length (middle part). This is owing to the increase in the amount of material to be deformed. The crack then propagates through the material as the resistance increases (the length of the fracture zone increases as the ligament length increases). A decrease in the ligament length implies an increase in the crack initiation length a0. To obtain the elastic work of fracture Wy, the area under the linear intercept of the force–displacement curve was plotted. The total work Wf under the curve was then divided by the area of the ligament portion (Lt). Thereby, the relationship between the total work Wf with each ligament was poled as shown in Figure 9. The linear regression of the total work done per ligament area is shown in Figure 9a, b for the DENTs with thicknesses of 5 mm and 1.8 mm, respectively. It was observed that the we as the intersection point between the extension of the linear regression data with the yaxis was separated. It was 273 kJ/m2 and 63 kJ/m2 for the plates with thicknesses of 5 mm and 1.8 mm, respectively. This demonstrated that the EWF sensitively tested the thickness. The higher value is owing to the MgO formation caused by oxygen penetration into the crack during the crack propagation. The MgO particles were distributed over the crack surfaces for the 5 mm-thick plates [41]. These solid microscopic particles bridge the crack surfaces, thereby reducing the crack propagation or advancement. At a low thickness, the fracture toughness was relatively low at 63 kJ/m2. This was owing to the thin crack surface having only a marginal amount of bridging by the crack. In addition, the Si would crack and thereby contribute to the fracture [42]. The elastic EWF Wey was shown in Figure 10a, b. It was 36 kJ/m2 and 11.5 kJ/m2 for the thicknesses of 5 mm and 1.8 mm, respectively. The value for 5 mm is near the standard range of fracture toughness in the ASM standard [43], whereas the value for 1.8 mm is below the range. This can be attributed to the fact that the fracture toughness is sensitive to the thickness. Even at a small thickness, there was a larger amount of localised stress. The modes of failure were net tension as shown in Figure 11a, b for all the specimen with thicknesses of 5 mm and 1.8 mm, respectively.

    Figure 8.  Load-displacement curve for DENT for (a) 5 mm thickness, (b) 1.8 mm thickness.
    Figure 9.  EWF fitting for DENT for (a) 5 mm thickness, (b) 1.8 mm thickness.
    Figure 10.  EWF data of elastic energy Wy for DENT for (a) 5 mm thickness, (b) 1.8 mm thickness.
    Figure 11.  Modes of failure for DENT for (a) 5 mm thickness, (b) 1.8 mm thickness.

    Figure 12a, b shows the curve of loading and displacement of the crack mouth. It was observed that for a 5 mm-thick slab (Figure 12a), the curve attains a mixing maximum and then, gradually decreases as the crack propagates. Meanwhile, for a thin slab with a thickness of 1.8 mm, it yields a flat plateau after attaining the maximum value. Figure 13a, b shows the crack propagation proceeding in a straight line with an undulating pattern owing to the plastic deformation by the crack surfaces.

    Figure 12.  Load-displacement curve of CT specimens of (a) 5 mm thickness (b) 1.8 mm thickness.
    Figure 13.  Modes of failure of CT test for thin-plate aluminium specimens: (a) 5 mm thickness (b) 1.8 mm thickness.

    The fracture toughness values (KIC) were calculated from the peak load values (MPam) using Eq 12. According to the ASTM E399 standard, the critical stress intensity factor for the breaking load (pQ) is given by [44]

    KIC=pQtWf(aW) (12)

    where (t) is the specimen thickness, mm; (W) is the specimen width, mm; (a) is the crack length, mm; (pQ) is the load at 5% secant; and f(aW) is the shape correction factor (see Eq 13):

    f(a/w)=2+a/w(1a/w)1.5[0.886+4.64(a/w)13.32(a/w)2+14.72(a/w)35.6(a/w)4 (13)

    The critical strain energy release rate JIC can be measured using Eq 14:

    JIC=KIC²E (14)

    where E is the elastic modulus. Value of the energy release rate, JIC, when the crack propagates (referred to as unstable crack growth). For materials and loaded configurations with an increasing resistance curve (or R-curve), JIC is not unique. Examples of increasing and flat R curves are shown in Figure 12. With an increasing R-curve, the JIC value is generally defined for crack initiation (assuming that a marginal amount of stable cracking occurs before unstable crack growth). This is analogous to the 0.2% proof stress derived from tensile tests on materials that do not have a distinct yield point (see Table 3).

    Table 3.  Fracture toughness corresponding to thickness.
    Thickness
    (t, mm)
    Crack length
    (a, mm)
    pQ
    (N)
    KIC
    (MPam)
    JIC
    (kJ/m2)
    5 34 3000 48.3 34.5
    1.8 34 600 26.88 10.6

     | Show Table
    DownLoad: CSV

    The shape of the R-curve was influenced significantly by the size and geometry of the component. The R-curve resulting from a crack in a thin sheet is generally steeper than that resulting from a crack in a thick sheet. This is because a thin sheet is generally loaded in the plane, whereas the material at the crack tip of a thick sheet is under stress and is simultaneously loaded in the plane. The JIC values were 34.5 kJ/m2 and 10.6 kJ/m2 for the specimens with thickness of 5 mm and 1.8 mm, respectively. The failure modes are illustrated in Figure 13. It was evidently a simple tension mode for the 5 mm-thick plate as shown in Figure 13-a. Here, there was a small bending through the specimen with the smaller thickness (1.8 mm) (see Figure 13-b).

    It was observed that the values of the elastic EWF wy were close to the release energy of the surface CT (see Table 4). The percentage variation was 5.8% in the case of the specimens with a thickness of 5 mm, whereas it was 8.4% for the 1.8 mm-thick slabs. This was because the CT indicates the fracture toughness at the beginning of the failure or at the beginning of the peak load based on a small plastic zone before the crack tip. This concept was developed using an elastic EWF. The elastic blow behaviour is characterised by a larger plastic deformation in front of the crack tip. Consequently, a linear fraction of elastic deformation exists in addition to the energy stored by the plastic deformation. Kobayshi and Yamada [45] determined the fracture toughness of a ductile damaged metal by correlating the average values of the initial and failure fracture toughness values over the entire R-curve. Test methods with different thicknesses yield different values. This is a generally debatable concept of fracture toughness measurement recommended by the ASTM standards [46]. The ASTM provides many methods for testing the fracture toughness of metals. Linear-elastic fracture mechanics refers to the fracture toughness of linear or elastic materials, whereas elastic–plastic fracture toughness (EPFM) refers to non-linear materials or materials with large plastic deformation (such as metals). The percentage error produced by the various methods was within a reasonable range from an engineer's perspective. In addition, the EWF method, which is based on the energy stored throughout the entire body, uses a different measurement criterion that does not depend on the LEFM and J-integral concepts. It rather depends on the strip length and the total work and energy stored in a cracked specimen [47,48] according to the EPFM [49]. It should be acknowledged that the EWF (we) considers the crack initiation resistance (GIC). If the variation in (GIC) is marginal [50], the relationship between the EWF (we) and J-integral (Ji) is valid according to the EPFM [51]. However, EWF is characterised by its simplicity of data reduction, sample preparation, and evaluation. This makes it attractive for measuring the fracture toughness of ductile thin films [52,54].

    Table 4.  Comparison between EWF and CT.
    Specimen Surface release energy JIC, kJ/m2 % of variation
    EWF CT
    5 mm 36 34.5 5.8
    1.8 mm 11.5 10.6 8.4

     | Show Table
    DownLoad: CSV

    The mechanical and fracture properties of a structural material (the aluminium alloy 5754-H11) were determined. It was observed that the concept of essential work of fracture (EWF) was reasonable for measuring the surface energy in the case of elastic–plastic fracture mechanics JIC. The EWF was affected by the sheet thickness. It was 273 kJ/m2 for the specimens with a thickness of 5 mm, and 63 kJ/m2 for the specimens with a thickness of 1.8 mm. These values represent the total fracture toughness of the plastic flow of ductile fracture materials (essential (elastic) and non-essential (plastic) fracture work). Therefore, the elastic work of fracture Wy could be determined as the area under the elastic zone of the load-displacement curves. Thus, the elastic EWF wy was 36 kJ/m² and 11.5 kJ/m2 for the specimens with thicknesses of 5 mm and 1.8 mm, respectively. These two values of elastic work of fracture are close to the values obtained with the conventional standard specimen CT. Here, the surface release energies JIC were 34.5 kJ/m2 and 10.6 kJ/m2 for the specimens with thicknesses of 5 mm and 1.8 mm, respectively. These values correspond to percentage errors of 5.8% and 8.4%, respectively, compared with the elastic EWF data. The large amount of Mg in the aluminium alloy 5754 was the main reason for the material embrittlement. It reduced the fracture toughness in the lower thickness. Meanwhile, the elements with Si formed a solid particle through the crack surfaces when these reacted with oxygen. It functioned as a bridge and thereby, reduced the crack propagation or crack propagation and then reinforced the fracture.

    The authors thank the members of Egyptian Aluminium Company (Egyptalum) in Nag Hammadi, Egypt, for providing the raw material.

    The authors declare no conflict of interest.



    [1] E. C. Aifantis, Internal length gradient (ILG) material mechanics across scales and disciplines, Adv. Appl. Mech., 49 (2016), 1–110. https://doi.org/10.1016/bs.aams.2016.08.001 doi: 10.1016/bs.aams.2016.08.001
    [2] G. Allaire, Homogenization and two-scale convergence, SIAM J. Math. Anal., 23 (1992), 1482–1518. https://doi.org/10.1137/0523084 doi: 10.1137/0523084
    [3] G. Allaire, M. Amar, Boundary layer tails in periodic homogenization, ESAIM Contr Optim CA, 4 (1999), 209–243. https://doi.org/10.1051/cocv:1999110 doi: 10.1051/cocv:1999110
    [4] L. Ambrosio, A. Carlotto, A. Massaccesi, Lectures on Elliptic Partial Differential Equations, Berlin: Springer, 2019.
    [5] A. Bensoussan, J. L. Lions, G. Papanicolaou, Asymptotic Analysis for Periodic Structures. Netherlands: Elsevier, 1978.
    [6] A. Berezovski, Internal variables representation of generalized heat equations, Contin. Mech. Thermodyn., 31 (2019), 1733–1741. https://doi.org/10.1007/s00161-018-0729-4 doi: 10.1007/s00161-018-0729-4
    [7] F. Bonetto, J. Lebowitz, L. Rey-Bellet, Fourier's law: a challenge for theorists, Mathematical Physics 2000, (2000), 128–150.
    [8] M. Cherdantsev, K. D. Cherednichenko, Two-Scale Γ–Convergence of Integral Functionals and its Application to Homogenisation of Nonlinear High-Contrast Periodic composites, Arch. Rational Mech. Anal., 204 (2012), 445–478. https://doi.org/10.1007/s00205-011-0481-4 doi: 10.1007/s00205-011-0481-4
    [9] C. I. Christov, On a higher-gradient generalization of Fourier's law of heat conduction, AIP Conference Proceedings, (2007), 11–22.
    [10] C. Ciorănescu, P. Donato, An Introduction to Homogenization, Oxford: Oxford University Press, 2000.
    [11] D. Ciorănescu, A. Damlamian, P. Donato, G. Griso, R. Zaki, The periodic unfolding method in domains with holes, SIAM J. Math. Anal., 44 (2012), 718–760. https://doi.org/10.1137/100817942 doi: 10.1137/100817942
    [12] D. Ciorănescu, A. Damlamian, G. Griso, Éclatement périodique et homogénéisation, Comptes Rendus Math., 335 (2022), 99–104. https://doi.org/10.1016/S1631-073X(02)02429-9 doi: 10.1016/S1631-073X(02)02429-9
    [13] D. Ciorănescu, A. Damlamian, G. Griso, The periodic unfolding method in homogenization, SIAM J. Math. Anal., 40 (2008), 1585–1620.
    [14] D. Ciorănescu, A. Damlamian, G. Griso, The periodic unfolding method, Theory and Applications to Partial Differential Problems, Singapore: Springer Singapore, 2018.
    [15] B. D. Coleman, V. J. Mizel, Thermodynamics and departures from Fourier's law of heat conduction, Arch. Ration. Mech. Anal., 13 (1963), 245–261. https://doi.org/10.1007/BF01262695 doi: 10.1007/BF01262695
    [16] A. Damlamian, An elementary introduction to periodic unfolding, Gakuto Int. Series, Math. Sci. Appl., 24 (2005), 1651–1684. https://doi.org/10.1377/hlthaff.24.6.1684 doi: 10.1377/hlthaff.24.6.1684
    [17] F. Demengel, G. Demengel, Functional Spaces for the Theory of Elliptic Partial Differential Equations, London: Springer-Verlag London, 2012.
    [18] A. Fehér, R. Kovács, On the evaluation of non-Fourier effects in heat pulse experiments, Int. J. Eng. Sci., 169 (2021), 103577. https://doi.org/10.1016/j.ijengsci.2021.103577 doi: 10.1016/j.ijengsci.2021.103577
    [19] A. Fehér, N. Lukács, L. Somlai, T. Fodor, M. Szücs, T. Fülöp, et al., Size effects and beyond-Fourier heat conduction in room-temperature experiments, J. Non-Equil. Thermody., 46 (2021), 403–411.
    [20] S. Forest, Milieux continus généralisés et matériaux hétérogènes, Paris: Presses des MINES, 2006.
    [21] S. Forest, E. C. Aifantis, Some links between recent gradient thermo-elasto-plasticity theories and the thermomechanics of generalized continua, Int. J. Solids Structures, 47 (2010), 3367–3376.
    [22] S. Forest, M. Amestoy, Hypertemperature in thermoelastic solids, CR MECANIQUE, 336 (2008), 347–353.
    [23] S. Forest, F. Pradel, K. Sab, Asymptotic analysis of heterogeneous Cosserat media, Int J Solids Struct, 38 (2001), 4585–4608.
    [24] P. Germain, La méthode des puissances virtuelles en mécanique des milieux continus, I: Théorie du second gradient, J. Mécanique, 12 (1973), 235–274.
    [25] P. Germain, The method of virtual power in continuum mechanics. Part 2: Microstructure, SIAM J. Appl. Math., 25 (1973), 556–575. https://doi.org/10.1016/0026-2714(73)90243-6 doi: 10.1016/0026-2714(73)90243-6
    [26] M. Giaquinta, L. Martinazzi, An Introduction to the Regularity Theory for Elliptic Systems, Harmonic Maps and Minimal Graphs, Berlin: Springer Science & Business Media, 2013.
    [27] G. Griso, Error estimates and unfolding for periodic homogenization, Asym. Anal., 40 (2004), 269–286.
    [28] D. Lukkassen, G. Nguetseng, P. Wall, Two-scale convergence, Int. J. Pure Appl. Math., 2 (2002), 35–86.
    [29] A. Mielke, On evolutionary Γ–convergence, In A. Muntean, J. Rademacher, A. Zagaris, editors, Macroscopic and Large Scale Phenomena: Coarse Graining, Mean Field Limits and Ergodicity, Cham: Springer, 2016,187–249.
    [30] R. D. Mindlin, N. N. Eshel, On first strain-gradient theories in linear elasticity, Int J Solids Struct, 4 (1968), 109–124. https://doi.org/10.1177/002188636800400105 doi: 10.1177/002188636800400105
    [31] S. Moskow, M. Vogelius, First-order corrections to the homogenised eigenvalues of a periodic composite medium. a convergence proof, P ROY SOC EDINB A, 127 (1997), 1263–1299. https://doi.org/10.1017/S0308210500027050 doi: 10.1017/S0308210500027050
    [32] G. Nguetseng, A general convergence result for a functional related to the theory of homogenization, SIAM J. Math. Anal., 20 (1989), 608–623. https://doi.org/10.1137/0520043 doi: 10.1137/0520043
    [33] G. Nika, A gradient system for a higher-gradient generalization of Fourier's law of heat conduction, Mod. Phys. Lett. B, (2023), 1–11. https://doi.org/10.1142/S0217984923500112 doi: 10.1142/S0217984923500112
    [34] G. Nika, A. Muntean, Effective medium theory for second-gradient nonlinear elasticity with chirality, arXiv: 2202.00644, [Preprint], (2022) [cited 2023 April 12]. Available from: https://doi.org/10.48550/arXiv.2202.00644
    [35] G. Nika, B. Vernescu, Rate of convergence for a multiscale model of dilute emulsions with non-uniform surface tension, Discrete Continuous Dyn Syst Ser B, 9 (2016), 1553–1564. https://doi.org/10.3934/dcdss.2016062 doi: 10.3934/dcdss.2016062
    [36] D. Onofrei, B. Vernescu, Error estimates for periodic homogenization with non-smooth coefficients, Asym. Anal., 54 (2007), 103–123. https://doi.org/10.1353/tcj.2007.0008 doi: 10.1353/tcj.2007.0008
    [37] L. S. Pan, D. Xu, J. Lou, Q. Yao, A generalized heat conduction model in rarefied gas, EPL, 73 (2006), 846–850. https://doi.org/10.1209/epl/i2005-10473-7 doi: 10.1209/epl/i2005-10473-7
    [38] D. Ruelle, A mechanical model for Fourier's law of heat conduction, Commun. Math. Phys, 311 (2012), 755–768.
    [39] E. Sanchez-Palencia, Non-homogeneous media and vibration theory, Lect. Notes Phys., 320 (1980), 57–65.
    [40] G. Stampacchia, The spaces L(p,λ), N(p,λ) and interpolation, Ann. Scuola Norm-Sci., 19 (1965), 443–462.
    [41] A. S. J. Suiker, C. S. Chang, Application of higher-order tensor theory for formulating enhanced continuum models, Acta Mech., 142 (2000), 223–234. https://doi.org/10.1007/BF01190020 doi: 10.1007/BF01190020
  • This article has been cited by:

    1. Cuili Li, Wan Zhou, Angel Gutierrez Ruiz, Yasaman Mohammadi, Qingning Li, Shuting Zhang, XiuJun Li, Guanglei Fu, Microfluidic immunoassays for point-of-care testing of SARS-CoV-2 antigens and antibodies, 2024, 177, 01659936, 117809, 10.1016/j.trac.2024.117809
    2. Michael Mallory, Jennifer E. Munt, Tara M. Narowski, Izabella Castillo, Edwing Cuadra, Nora Pisanic, Paul Fields, John M. Powers, Alexandria Dickson, Rohan Harris, Richard Wargowsky, Seamus Moran, Ahmed Allabban, Kristin Raphel, Timothy A. McCaffrey, James D. Brien, Christopher D. Heaney, John E. Lafleur, Ralph S. Baric, Lakshmanane Premkumar, COVID-19 point-of-care tests can identify low-antibody individuals: In-depth immunoanalysis of boosting benefits in a healthy cohort, 2024, 10, 2375-2548, 10.1126/sciadv.adi1379
    3. Simone Cavalera, Fabio Di Nardo, Thea Serra, Valentina Testa, Claudio Baggiani, Sergio Rosati, Barbara Colitti, Ludovica Brienza, Irene Colasanto, Chiara Nogarol, Domenico Cosseddu, Cristina Guiotto, Laura Anfossi, A semi-quantitative visual lateral flow immunoassay for SARS-CoV-2 antibody detection for the follow-up of immune response to vaccination or recovery, 2024, 12, 2050-750X, 2139, 10.1039/D3TB02895J
    4. Shanil Gandhi, Xhorxhina Shaulli, Jeppe Fock, Frank Scheffold, Rodolphe Marie, IgG and IgM differentiation in a particle-based agglutination assay by control over antigen surface density, 2024, 8, 2473-2877, 10.1063/5.0196224
    5. Seyyed Mojtaba Mousavi, Masoomeh Yari Kalashgrani, Ahmad Gholami, Navid Omidifar, Mojtaba Binazadeh, Wei-Hung Chiang, Recent Advances in Quantum Dot-Based Lateral Flow Immunoassays for the Rapid, Point-of-Care Diagnosis of COVID-19, 2023, 13, 2079-6374, 786, 10.3390/bios13080786
    6. Tadej Slatinek, Janez Slapnik, Upcycling of SARS-CoV-2 Rapid Antigen Test Cassettes into Flame Retardant Plastics, 2024, 17, 1996-1944, 2384, 10.3390/ma17102384
    7. Bartolomeo Della Ventura, Rosalba Fittipaldi, Margherita Borriello, Ivan Gentile, Giuseppe Portella, Daniela Terracciano, Antonio Vecchione, Diego Ingrosso, Raffaele Velotta, 2024, Rapid Colorimetric Assay for Salivary Detection of Sars-CoV-2, 979-8-3503-6925-0, 1, 10.1109/SAS60918.2024.10636688
    8. Stella Schuster da Silva, Mylena Lemes Cunha, Lucas B. Ayres, Carlos D. Garcia, Lucas Blanes, Advancements and future directions in cardiac biomarker detection using lateral flow assays, 2023, 15, 1759-9660, 3610, 10.1039/D3AY01081C
    9. Rushi Liang, Aiping Fan, Feiqian Wang, Yajing Niu, Optical lateral flow assays in early diagnosis of SARS-CoV-2 infection, 2024, 40, 0910-6340, 1571, 10.1007/s44211-024-00596-6
    10. Le Zhang, Shihui Wen, Jawairia Umar Khan, Yuan Liu, Mahnaz Maddahfar, Jiajia Zhou, Dayong Jin, Ultrasensitive Rapid Antigen Test by Geometric Lateral Flow Assays and Highly Doped Upconversion Nanoparticles, 2024, 96, 0003-2700, 16581, 10.1021/acs.analchem.4c02379
    11. Di Jiao, Fan Jiao, Zhen-Jie Qian, Lin Luo, Yu Wang, Yu-Dong Shen, Hong-Tao Lei, Zhen-Lin Xu, Formation and Detection of Gizzerosine in Animal Feed Matrices: Progress and Perspectives, 2024, 72, 0021-8561, 3247, 10.1021/acs.jafc.3c05973
    12. Valeria Secchi, Alice Armanni, Linda Barbieri, Antonia Bruno, Alessandro Colombo, Sara Fumagalli, Ekaterina A. Kukushkina, Roberto Lorenzi, Letizia Marchesi, Hind Moukham, Alberto Paleari, Alessandra Ronchi, Giulia Tomaino, Farida Tripodi, Miriam Colombo, Laura Sironi, Angelo Monguzzi, Advanced techniques and nanotechnologies for point-of-care testing, 2025, 6, 2673-3013, 10.3389/fnano.2024.1465429
    13. Shuma Yazaki, Misaki Komatsu, Jinhua Dong, Hiroshi Ueda, Ryoichi Arai, Crystal Structures of Antigen-Binding Fragment of Anti-Osteocalcin Antibody KTM219, 2025, 26, 1422-0067, 648, 10.3390/ijms26020648
    14. Jianqiang Ma, Scott Kaniper, Yuliya Vabishchevich, Nana Nyantakyi, Dorret Lynch, Fan Chun, Hai-Lung Dai, Glenn S. Gerhard, SARS-CoV-2 neutralizing antibody determination after vaccination using spectrophotometric measurement of lateral flow immunochromatography, 2025, 15, 2045-2322, 10.1038/s41598-025-90730-9
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1690) PDF downloads(47) Cited by(1)

Figures and Tables

Figures(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog