Simplified Fracture Mechanics Analysis at the Zinc–Adhesive Interface in Galvanized Steel–CFRP Single-Lap Joints

TL;DR

This study simplifies the analysis of how cracks form at the zinc-adhesive interface in joints between galvanized steel and carbon fiber composites.

Contribution

A simplified analytical framework is introduced to qualitatively assess geometric effects on fracture behavior in galvanized steel–CFRP joints.

Findings

Increasing overlap length reduces crack driving force (Gint), while increasing adhesive thickness raises it.

Thicker adhesive layers shift fracture mode from shear- to opening-dominated.

The simplified model captures qualitative trends despite limited quantitative accuracy.

Abstract

Adhesively bonded joints between galvanized steel and carbon fiber-reinforced polymers (CFRPs) are critical in modern lightweight structures, but their performance is often limited by failure at the zinc–adhesive interface. This study presents a parametric analysis to investigate the influence of key geometric parameters on interfacial cracking in a single-lap joint (SLJ) configuration, employing a simplified analytical methodology based on Interface Fracture Mechanics (IFM). The model combines the Goland–Reissner approach for estimating crack-tip loads with highly simplified, constant shape functions to calculate the energy release rate (Gint) and phase angle (ψ). To provide a practical reference, experimental data from shear tests on S350 GD galvanized steel bonded to CFRP were used to estimate the range of interfacial fracture toughness for this material system. The parametric results demonstrate that, for a constant load, increasing the overlap length reduces the crack driving force (Gint), while increasing the adhesive thickness raises it. Crucially, the model indicates that a thicker adhesive layer shifts the fracture mode from shear- to opening-dominated, a trend consistent with the established mechanics of SLJs, where increased joint rotation amplifies peel stresses. The study concludes that while the use of constant shape functions limits the model’s quantitative accuracy, this simplified analytical framework effectively captures the qualitative influence of key geometric parameters on the joint’s fracture behavior. It serves as a valuable and resource-efficient tool for preliminary design explorations and for interpreting experimentally observed failure trends in galvanized steel–CFRP joints.

Figures

Full text

1. Introduction

Adhesive bonding is increasingly recognized as an effective method for joining diverse structural materials, including steel with carbon fiber-reinforced polymers (CFRPs) [1,2,3], especially in Rzeszut et al. [4]. The application of CFRP to reinforce or connect steel elements offers numerous advantages, such as mass reduction, increased load-bearing capacity, and improved fatigue and corrosion resistance [5,6]. However, designing durable and reliable joints of this type, particularly involving galvanized steel, presents several engineering challenges [7]. The zinc layer, applied to protect steel from corrosion [8,9], introduces an additional material interface (zinc–adhesive) that can become a site for damage initiation [10,11,12,13]. Understanding the failure mechanisms at this specific interface and the influence of various structural and material parameters on its behavior is crucial for ensuring the integrity of the entire joint [14,15].

Interface Fracture Mechanics (IFM) provides theoretical tools for analyzing cracks located at the boundary between two dissimilar materials [16,17,18]. This approach allows for the quantification of conditions for crack initiation and propagation through parameters such as stress intensity factors ( $[eqn]$ ), energy release rate ( $[eqn]$ ), and phase angle ( $[eqn]$ ), which describes the mixity of fracture modes [19,20]. Accurate IFM analyses are often complex and require advanced numerical methods, such as the finite element method employing special techniques for extracting fracture parameters [21]. However, there is a need to develop simplified analytical models that, despite certain approximations, would allow for effective parametric studies and preliminary assessment of joint behavior.

A critical challenge in the application of adhesively bonded steel–CFRP systems is the prevention of premature failure, with debonding at the material interface being recognized as a primary and critical failure mode [22]. The analysis of such failures is complex, as crack propagation along the interface typically occurs under mixed-mode loading conditions, involving a combination of normal (Mode I) and shear (Mode II) stresses [23]. Significant research has been dedicated to understanding this phenomenon. For instance, Zeng et al. (2018) [24] provided detailed investigations into the interfacial behavior and debonding mechanisms in full-scale CFRP-strengthened steel beams, highlighting the importance of fracture mechanics in predicting joint performance [24]. The field continues to evolve, with recent studies exploring the fracture characteristics of various bi-material combinations, including the Mode I tensile fracture behavior of metal–CFRP joints [25]. This body of work underscores the necessity of applying fracture mechanics principles to accurately characterize the interfacial behavior and ensure the long-term reliability of these hybrid structures.

Experimental studies on adhesively bonded galvanized steel–CFRP joints, such as those conducted by Dybizbański et al. (2024) [26], provide valuable data on the failure mechanisms and load-bearing capacities of these systems. In the cited work, various failure modes were observed, including failure at the steel–adhesive interface (interpreted here as zinc–adhesive), fiber rupture, and mixed-mode failures. It was also noted that an increase in adhesive thickness led to a decrease in joint strength, while an increase in overlap length generally improved the load-bearing capacity.

The motivation for the present work is the need for a systematic analysis of the influence of key geometric and material parameters on the conditions for crack initiation at the zinc–adhesive interface, using an analytical tool that is less resource-intensive than full numerical modelling. Such a parametric study, despite its inherent simplifications, can provide important insights into the sensitivity of the joint to design changes and aid in the interpretation of experimental results.

The main aim of this work is to conduct a parametric study of the influence of selected geometric parameters (overlap length $[eqn]$ ; adhesive layer thickness $[eqn]$ ; initial crack length $[eqn]$ ) on the energy release rate $[eqn]$ and phase angle $[eqn]$ at the zinc–adhesive interface in steel–CFRP single-lap joints (SLJs), using a simplified analytical methodology based on IFM.

2. Theoretical Background

2.1. Interface Fracture Mechanics (IFM)–Key Parameters

Interface Fracture Mechanics (IFM) describes the behavior of cracks located at the boundary of two dissimilar materials. In the context of this work, material 1 is the zinc layer on the steel surface, and material 2 is the epoxy adhesive.

The stress field near an interfacial crack tip is characterized by a complex stress intensity factor $[eqn]$ [16,27]. The components $[eqn]$ and $[eqn]$ correspond to the conventional Mode I (opening) and Mode II (in-plane shear).

The mismatch in elastic properties between material 1 (zinc) and material 2 (adhesive) is described by the Dundurs parameters $[eqn]$ i $[eqn]$ [28,29] and the bimaterial constant $[eqn]$ (oscillation index) [16]:

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

where $[eqn]$ for plane strain (PS), $[eqn]$ is the shear modulus, and $[eqn]$ for PS.

The energy release rate for an interfacial crack $[eqn]$ is related to $[eqn]$ i $[eqn]$ by $[eqn]$ [16,20], where the effective Young’s modulus for the interface is calculated as follows [17]:

$$[eqn]$$

The relative proportion of fracture modes is characterized by the phase angle $[eqn]$ , defined with respect to a characteristic length $[eqn]$ [16,17]:

$$[eqn]$$

In this work $[eqn]$ is taken as the adhesive layer thickness $[eqn]$ .

Interfacial Fracture Toughness $[eqn]$ is the critical value of $[eqn]$ required for crack propagation, which generally depends on the phase angle $[eqn]$ [11,17]. Figure 1 presents those key concepts of IFM.

2.2. Goland–Reissner Model for Single-Lap Joints

The Goland–Reissner model [9,31] is a classical analytical model used to determine the stress distribution in the adhesive layer and the forces and moments in the adherends of lap joints, accounting for the eccentricity of the applied load. This leads to rotation of the overlap and induces bending moments in the adherends, as well as peel (normal) and shear stresses in the adhesive.

In this work, to estimate the bending moment $[eqn]$ (per unit width) in the steel adherend at the edge of the overlap, the following simplified formula is adopted:

$$[eqn]$$

where $[eqn]$ is the force per unit width of the joint, and $[eqn]$ is the steel adherend thickness. The coefficient $[eqn]$ accounts for the adhesive flexibility and joint geometry; for the analysis, an approximate form is calculated as follows:

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

where $[eqn]$ is the overlap length, $[eqn]$ is Young’s modulus of steel, and $[eqn]$ is the shear modulus of the adhesive.

2.3. Simplified Shape Functions for Stress Intensity Factors

To calculate the stress intensity factors $[eqn]$ and $[eqn]$ for an edge crack at an interface, subjected to stresses $[eqn]$ (average axial stress in the adherend) and $[eqn]$ (bending stress in the adherend), the following general form is used [32,33]:

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ is the crack length. In this work, to achieve maximum simplification of the analytical methodology for the parametric study, the following very general, constant values for the shape functions are adopted:

• $[eqn]$ : This value is a classical solution for the Mode I stress intensity factor for an edge crack in a semi-infinite homogeneous plate subjected to uniform tension, for $[eqn]$ [32,33].
• $[eqn]$ : This assumption arises from the fact that, in a homogeneous material, a symmetric tensile load on an edge crack does not induce Mode II.
• $[eqn]$ : For an edge crack in a homogeneous material subjected to bending, the Mode I shape function depends on $[eqn]$ [32]. The value of 1.0 is adopted as a further simplification for short cracks.
• $[eqn]$ : This value is taken as an estimate to account for a small potential Mode II component induced by the bending of an edge crack, particularly at an interface [17,34].

It must be strongly emphasized that these constant shape function values are a significant simplification and do not reflect the complex dependence on the relative crack length $[eqn]$ or the interface material mismatch parameters $[eqn]$ . They represent a compromise between accuracy and the goal of creating a maximally simple analytical methodology suitable for a parametric study. Their impact on the results and limitations will be discussed later in the paper.

It should be emphasized that this work focuses on the energy release rate ( $[eqn]$ ) and the phase angle ( $[eqn]$ ), rather than on the individual stress intensity factors ( $[eqn]$ , $[eqn]$ ). For cracks at bimaterial interfaces, the near-tip stress fields exhibit an oscillatory character, which makes the separation of $[eqn]$ and $[eqn]$ non-unique and dependent on an arbitrarily chosen length scale. The analytical framework based on the ( $[eqn]$ , $[eqn]$ ) pair was developed to overcome this ambiguity, providing a more physically robust and unambiguous description of the crack driving force ( $[eqn]$ ) and the mode mixity ( $[eqn]$ ). It is therefore the preferred, modern approach for analyzing interfacial crack problems.

3. Methodology

3.1. Analyzed Joint Description and Model Idealizations

The analysis focuses on a single-lap joint (SLJ) with geometry and material characteristics corresponding to the specimens tested by Dybizbański et al. (2024) [26]. SLJ geometry is presented in Figure 2.

The joint consists of a steel adherend (S350 GD, thickness $[eqn]$ ), the surface of which is coated with a zinc layer (thickness $[eqn]$ ), an epoxy adhesive layer (SikaDur 330, thickness $[eqn]$ ), and a carbon fiber-reinforced polymer (CFRP SikaWrap 230 C) adherend. Both SikaWrap 230 C and SikaDur 330 were provided by Sika Poland sp. z o.o. (Warsaw, Poland).

An initial short crack of length a is assumed to be located at the interface between the zinc layer and the adhesive, at one end of the overlap. The crack length a is one of the parameters in the analysis. All materials (steel, zinc, adhesive) are treated as linear elastic and isotropic. The mechanical properties of zinc $[eqn]$ and the adhesive $[eqn]$ are crucial for determining the interface parameters. The overlap length is denoted as $[eqn]$ , and the joint width as $[eqn]$ .

Consistent with previous considerations, in the global SLJ analysis using the Goland–Reissner model, the influence of the zinc layer thickness on the global stiffness of the steel adherend is neglected, and $[eqn]$ is used as the steel thickness. For the local IFM analysis at the zinc–adhesive interface, the mechanical properties of zinc are used to define the Dundurs parameters. In the arguments of the shape functions $[eqn]$ , $[eqn]$ is used as the dimension characterizing the adherend whose stresses generate the loading at the crack tip located on the very thin zinc layer.

3.2. Calculation Steps of the Analytical IFM Model

The calculation procedure for IFM parameters for a given external force $[eqn]$ (over the entire joint width) is performed in four steps. Step 1 includes the calculation of elastic parameters for the zinc–adhesive interface. Based on the properties of zinc (material 1) and adhesive (material 2), for an assumed plane strain (PS) condition, the following are calculated: $[eqn]$ and $[eqn]$ , according to the formulas presented in point 2.1. In step 2 the estimation of loads $[eqn]$ using the Goland–Reissner model are carried out for the force per unit width $[eqn]$ , the bending moment $[eqn]$ at the overlap edge in the steel adherend is calculated using the formula $[eqn]$ with the adopted form for k′. Then, applying simplifications, $[eqn]$ and $[eqn]$ are assumed as the axial force and bending moment (per unit width) in the steel adherend cross-section at the crack tip. The calculation of $[eqn]$ and $[eqn]$ , using simplified shape functions, are performed in step 3. Then, based on $[eqn]$ and $[eqn]$ , the $[eqn]$ and $[eqn]$ are calculated. Next, using the adopted constant shape function values ( $[eqn]$ , $[eqn]$ , $[eqn]$ , $[eqn]$ ), $[eqn]$ and $[eqn]$ are calculated according to the formulas in Section 2.3. Finally in step 4 the calculation of $[eqn]$ are performed. The energy release rate $[eqn]$ and the phase angle $[eqn]$ (with $[eqn]$ ) are calculated based on $[eqn]$ , $[eqn]$ , $[eqn]$ and $[eqn]$ according to the formulas in Section 2.1.

3.3. Parametric Study Plan

A parametric study was conducted to investigate the influence of key geometric parameters on the calculated $[eqn]$ and $[eqn]$ values. The parameters as reference values for material and geometric properties were adopted based on Dybizbański et al. (2024) [26] and general material data. The joint width w was taken from the experimental specimens. The varied parameters that were taken into account include the following: overlap length $[eqn]$ , adhesive layer thickness $[eqn]$ , and initial crack length $[eqn]$ . For each parametric analysis (varying one parameter while keeping others at baseline values), calculations of $[eqn]$ and $[eqn]$ were performed for a constant, reference external force $[eqn]$ . The choice of this force (e.g., corresponding to the average failure load from one of the experimental series [26] or a typical design load) allows for direct comparison of the parameter’s influence on the crack driving force. The primary output quantities analyzed were $[eqn]$ and $[eqn]$ .

3.4. Processing of Experimental Data for Gint,cψ Estimation

Independently of the parametric study, experimental data concerning failure loads $[eqn]$ from Dybizbański et al. (2024) [26] were processed using the described analytical methodology (Steps 1–4 from Section 3.2). For each specimen that underwent failure interpreted as adhesive at the zinc–adhesive interface, a pair ( $[eqn]$ ) was calculated. The resulting set of points was used to attempt an estimation of a segment of the interfacial fracture toughness curve $[eqn]$ for the studied material system.

4. Results

4.1. Estimation of the Interfacial Fracture Toughness Curve Segment, Gint,cψ, from Experimental Data

Experimental failure load data ( $[eqn]$ ) for single-lap joints (SLJs) composed of S350 GD galvanized steel, SikaDur 330 adhesive, and SikaWrap 230 C carbon fiber textile, as reported by Dybizbański et al. (2024) [26], were processed using the analytical IFM model (Section 3.2). For each specimen that exhibited failure at the steel–adhesive interface (interpreted as the zinc–adhesive interface), a pair of ( $[eqn]$ ) values was calculated. The resulting data points are plotted in Figure 3.

As observed in Figure 3, the calculated phase angles ( $[eqn]$ ) for the experimental failure points fall within a range of approximately 10° to 14°. While the individual data points show considerable scatter, varying from approximately 0.1 N/mm to over 0.5 N/mm, the addition of a polynomial trendline (red dotted line) reveals an underlying non-linear relationship. The trend suggests that the interfacial fracture toughness ( $[eqn]$ ) is dependent on the phase angle, reaching a maximum value of approximately 0.4 N/mm at an optimal phase angle of around $[eqn]$ . For phase angles deviating from this optimum, the toughness appears to decrease. This parabolic trend is consistent with the theoretical framework of Interface Fracture Mechanics, although it is partially obscured by the significant experimental variability, which can be attributed to factors discussed in Section 5.2.

4.2. Influence of Overlap Length (L) on IFM Parameters

To investigate the effect of overlap length ( $[eqn]$ ) on the crack driving force, calculations of $[eqn]$ and $[eqn]$ were performed for a continuous range of $[eqn]$ values, from 0 mm to 100 mm. For this analysis, other parameters were kept constant at their baseline values, and a constant reference applied load $[eqn]$ was used to isolate the effect of the overlap length. The results for the energy release rate, $[eqn]$ , for different thickness of the steel adherend $[eqn]$ are presented in Figure 4.

One can notice that for increased thickness of the steel adherend, the energy release rate is minimized for the whole considered range of overlap length. Another interesting observation is that the energy release rate becomes constant for overlap lengths greater than 20 mm.

4.3. Influence of Adhesive Layer Thickness (tk) on IFM Parameters

The thickness of the adhesive layer is a critical parameter known to significantly influence the mechanical performance and failure behavior of adhesively bonded joints. To investigate this effect within the framework of the developed analytical model, a parametric study was conducted where the adhesive thickness ( $[eqn]$ ) was varied from 0 to 1.0 mm. For this analysis, other parameters were held constant at their baseline values ( $[eqn]$ , $[eqn]$ ), and a constant reference applied load $[eqn]$ was used to isolate the influence of $[eqn]$ on the Interfacial Fracture Mechanics parameters.

The results of this analysis are presented in Figure 5 ( $[eqn]$ vs. $[eqn]$ ).

In the case of the investigation of the influence of adhesive thickness on energy release rate, the increased thickness of the steel adherend caused the reduction in energy release rate for the whole considered range of overlap length. Moreover, the energy release rate becomes almost constant for a large thickness of the steel adherend, and for the small, the result are of a nonliteral form.

Figure 6 presents phase angle–adhesive thickness relationships ( $[eqn]$ vs. $[eqn]$ ) for different adherend thicknesses.

It suggests that the increase in the adhesive layer thickness changes the fracture mode from shear-dominated to opening-dominated. Moreover, an increase in adherend thickness escalates this effect.

4.4. Influence of Assumed Initial Crack Length (a) on IFM Parameters

In any analysis based on Linear Elastic Fracture Mechanics (LEFM), it is necessary to assume the presence of an initial crack or flaw. The size of this assumed crack, $[eqn]$ , can significantly influence the calculated fracture mechanics parameters. To assess the sensitivity of the developed analytical model to this parameter, a study was conducted where $[eqn]$ was varied from a very small value (0.01 mm) up to 1.0 mm, while other parameters were held constant ( $[eqn]$ , $[eqn]$ , and a reference applied load $[eqn]$ ).

The results of this analysis are presented in Figure 7 ( $[eqn]$ vs. $[eqn]$ ).

Based on the conducted analysis, it can be noticed that for all thicknesses of the steel adherend, the energy release rate remains proportional to the crack size. Furthermore, the minimum value of the energy release rate is obtained for the large thickness of the steel adherend.

5. Discussion

The parametric study, conducted using the simplified analytical Interface Fracture Mechanics (IFM) methodology, provided key insights into how geometric variables influence the crack driving force ( $[eqn]$ ) and fracture mode mixity ( $[eqn]$ ) at the zinc–adhesive interface in single-lap joints. This analysis allows for a deeper interpretation of the experimental results reported by Dybizbański et al. (2024) [26].

5.1. Model Simplifications and Limitations

The analytical methodology employed in this study was intentionally simplified to serve as a tool for rapid, qualitative assessment of parametric trends. It is crucial, however, to acknowledge the limitations stemming from these simplifications. The most significant approximation is the use of constant, general-purpose shape functions to calculate the stress intensity factors. In reality, these functions are complex and depend on the relative crack length (a/t_a) and the elastic mismatch between the materials, as described by the Dundurs parameters. The adopted constant values are classical solutions for edge cracks in homogeneous materials under specific conditions and do not capture the full complexity of the interfacial stress field [22]. This choice inherently limits the quantitative accuracy of the calculated energy release rate values. Therefore, the results of this model should be interpreted as indicative of qualitative trends rather than as precise numerical predictions. For enhanced quantitative accuracy, a finite element-based calibration of the shape functions for the specific joint geometry and material combination would be a necessary and valuable next step.

A second important simplification is the treatment of the thin (≈9 µm) zinc layer as a discrete, isotropic elastic material. Given that its thickness is orders of magnitude smaller than the adherend and adhesive dimensions, a more physically representative approach for a high-fidelity analysis would be to model it as a compliant interphase or using a gradient property model. Such advanced models, often employing cohesive zone elements, can better capture the specific compliance and potential failure mechanisms within the coating itself. The first-order approximation used in this study was chosen for its consistency with the overall simplified analytical framework, which assumes that failure occurs perfectly at the zinc–adhesive interface. Investigating more sophisticated interphase models remains a compelling and important avenue for future research.

5.2. Analysis of Experimental Fracture Toughness Data

The significant scatter observed in the experimental interfacial fracture toughness data, as plotted in Figure 3, is a key finding that warrants further discussion. This variability, while substantial, is not unexpected and can be attributed to a combination of factors related to the material’s microstructure, the adhesive bonding process, and the idealizations within the analytical model. It is important to note that the data points correspond exclusively to specimens that exhibited failure at the zinc–adhesive interface, as verified by post-failure visual inspection [26].

A primary source of this scatter is the complex, non-homogeneous nature of the hot-dip galvanized zinc layer itself. The coating is not a monolithic material but a multi-phase system comprising several brittle zinc–iron intermetallic phases (such as Γ and δ phases) beneath the outer pure zinc layer [35]. The fracture behavior and damage mechanisms within the zinc coating are highly sensitive to these microstructural features [36,37]. For instance, research has demonstrated a direct relationship between the fracture toughness of zinc coatings and their grain boundary distribution, with different boundary types exhibiting varied cracking resistance [38]. Consequently, the initiation of a crack at the zinc–adhesive interface is not governed by a single material property but is highly dependent on these local microstructural variations, which are a significant source of the scatter observed in the macroscopic joint strength.

Furthermore, such variability is a well-documented characteristic of adhesive joint testing in general. Extensive reviews of adhesion science have shown that high coefficients of variation, often between 20% and 50%, are common due to a multitude of factors, including microscopic voids, non-uniformities in the bondline thickness, and subtle variations in the curing process [39]. The apparent strength of a joint is also highly dependent on the bending stiffness of the adherends, which directly affects the peel stress distribution and can introduce further variability [40]. Finally, the simplifications inherent in the analytical model, particularly the use of constant shape functions as discussed in Section 5.1, contribute to the calculated scatter by not capturing the full complexity of the stress state at the crack tip for each specific specimen. Considering these factors, the plot should not be interpreted as a precise, single-valued failure locus but rather as an estimation of the range of interfacial fracture toughness for this specific material system under the given test conditions.

The relationship presented in Figure 3, with the addition of a trendline, allows for a more nuanced interpretation of these results. The underlying parabolic trend suggests a fundamental, mechanics-driven relationship between fracture toughness and mode mixity, which is consistent with IFM theory. This indicates that while local material inconsistencies and experimental variability dictate the precise failure load of any single specimen, the average behavior of the interface is governed by the phase angle, with a clear peak in toughness under a specific mixed-mode condition. Therefore, the plot should be interpreted as showing two coexisting phenomena: a deterministic, mechanics-based trend (the parabola) superimposed with stochastic, material- and process-induced variability (the scatter).

5.3. Effect of Adhesive Thickness on Fracture Mode

The parametric study indicates a clear trend where an increase in adhesive thickness shifts the fracture mode towards being more opening-dominated, as reflected by the rising phase angle $[eqn]$ . While this result is derived from a simplified analytical model, the underlying physical mechanism is well-established in the literature on single-lap joints (SLJs). The eccentric load path inherent in the SLJ configuration induces joint rotation, which in turn generates significant peel stresses (the Mode I component) at the edges of the overlap—a phenomenon extensively analyzed since the work of Hart-Smith [41].

A thicker adhesive layer results in a more compliant bondline. This increased compliance allows for greater joint rotation under load, which amplifies the bending moment in the adherends and, consequently, the magnitude of the peel stresses at the bond ends [42]. This effect often leads to a reduction in the overall strength of the joint as the adhesive thickness increases [43]. This shift in the local stress state, with a higher contribution from Mode I, naturally favors a more opening-dominated fracture. Indeed, several studies have confirmed that as the bond thickness increases, the entire fracture mechanism within the adhesive joint can change completely [44]. Therefore, the qualitative trend captured by the presented model is a correct reflection of the fundamental mechanics of single-lap joints, corroborated by extensive experimental and numerical research.

5.4. Comparison of Parametric Trends with the Existing Literature

The parametric trends identified in this study, despite being derived from a simplified model, show strong agreement with the findings reported in more complex experimental and numerical investigations of CFRP-steel bonded joints. A significant body of research, including the work of Zeng et al. (2018) [24] and Teng et al. (2012) [22], has focused on debonding failures in large-scale CFRP-strengthened steel structures, establishing fracture mechanics as the essential framework for analyzing interfacial behavior. While these studies often involve more complex geometries, the fundamental mechanical principles governing the influence of joint parameters are consistent with the trends observed in the presented analysis.

Presented finding that increasing the overlap length reduces the energy release rate ( $[eqn]$ ) for a constant load is consistent with the well-known principle that stress concentrations in single-lap joints are localized at the ends of the overlap. For longer overlaps, a central portion of the bondline remains practically unloaded, which effectively distributes the load and reduces the crack driving force at the joint’s edge [45]. This explains the observed plateau in $[eqn]$ for longer overlap lengths in the presented analysis and aligns with findings that joint strength does not increase proportionally with overlap length, especially for brittle adhesives.

Similarly, the model’s prediction that a thicker adhesive layer increases $[eqn]$ and shifts the mode mixity towards opening-dominated fracture (a higher phase angle $[eqn]$ ) is strongly supported by the literature. As detailed in Section 5.3, this phenomenon is attributed to increased joint rotation and amplified peel stresses, a trend widely reported in studies on metal–composite joints [41,42]. The consistency of presented results with the broader literature findings confirms that the simplified analytical model, while not quantitatively precise, successfully captures the correct qualitative relationships between key geometric parameters and the governing fracture mechanics parameters. This validates its utility as a tool for preliminary research and design exploration and for interpreting experimentally observed failure trends.

Furthermore, the parametric trends predicted by the model provide a direct, mechanical explanation for the key experimental observations reported in. The experimental study found that joint strength decreased with increasing adhesive thickness and generally increased with longer overlaps. The model aligns perfectly with these findings on a qualitative basis: a thicker adhesive increases the crack driving force ( $[eqn]$ ) for a given load (Figure 5), which means the critical fracture toughness is reached at a lower external load, thus reducing strength. Conversely, a longer overlap decreases $[eqn]$ (Figure 4), requiring a higher applied force to initiate failure. This strong consistency confirms the model’s utility in capturing the fundamental mechanisms governing the joint’s behavior and enhances the interpretation of the experimental results.

6. Conclusions

This paper presented a parametric analysis of interfacial cracking in steel–CFRP single-lap joints using a simplified IFM methodology, contextualized with experimental data. The analysis leads to the following key conclusions:

• The experimental data, when analyzed with a trendline, suggests a non-linear, parabolic relationship between the interfacial fracture toughness ( $[eqn]$ ) and the phase angle ( $[eqn]$ ). The interface exhibits a peak fracture toughness at a phase angle of approximately $[eqn]$ , indicating that the joint’s resistance to cracking is highest under a specific mixed-mode condition. The significant scatter around this trend highlights the influence of local microstructural variations and process-induced variability.
• The parametric model demonstrates that for a constant load, increasing the overlap length reduces the crack driving force ( $[eqn]$ ), while increasing the adhesive thickness raises it. These trends provide a direct mechanical explanation for experimental observations where longer overlaps increased joint strength and thicker adhesives decreased it, thus validating the model’s qualitative predictive capabilities.
• A key finding, strongly supported by established joint mechanics, is that a thicker adhesive layer shifts the fracture mode from shear-dominated to opening-dominated (a higher phase angle $[eqn]$ ). This is attributed to the increased joint rotation that amplifies peel stresses at the bond ends.
• The study confirmed the model’s sensitivity to the assumed initial crack length, with $[eqn]$ being directly proportional to it, which is characteristic of LEFM-based analyses. Furthermore, increasing the thickness of the steel adherend was found to consistently reduce the crack driving force across all analyzed geometric parameters.
• In summary, the simplified analytical model, despite its acknowledged limitations—primarily the use of constant, general-purpose shape functions—effectively captures the qualitative influence of key geometric parameters on the joint’s fracture behavior. It provides a sound mechanical basis for interpreting experimentally observed failure trends and serves as a useful, resource-efficient tool for preliminary research and design explorations.