1. Introduction
The demand for adhesive bonding in structural applications is growing in various industries. Adhesives can vary in terms of their ductility, brittleness, temperature resistance, water absorption, and whether they are made of thermoset or thermoplastic polymers [
1]. Epoxy adhesives are widely utilized for structural purposes due to their excellent bonding capabilities to various materials and ease of processing [
2]. The market value of epoxy was USD 5.9 billion in 2019 and is projected to reach USD 10.3 billion by 2027 [
3]. Despite their widespread use and economic importance, the non-biodegradable nature of epoxy adhesives raises environmental concerns especially when adhesively bonded structures reach their end-of-life (ELV) [
3,
4]. Typically, polymers and composite materials end up in landfills or are incinerated, which are the least favored methods of disposal according to the waste management hierarchy [
5]. By 2050, it is estimated that the world will face up to 43 million tons of waste from decommissioned wind turbine blades which are predominantly made of glass fiber-reinforced polymers (GFRPs) [
6]. With the electric vehicle market growing, structural adhesives play a significant role in assembling battery packs. This growth could lead to an increase in battery waste. For example, 1 million electric vehicles could generate 250,000 tons of used batteries [
7]. With over 16 million electric vehicles as of 2022, it is important to adopt a circular economy now to manage waste and recover valuable battery materials before it becomes unmanageable. These statistics and projections underscore the urgency for the development of modern adhesives exhibiting switchable (i.e., reversible) adhesion characteristics while maintaining high performance and strength properties under service conditions, as depicted in
Figure 1 [
8]. This approach helps prevent damage and support the reuse, recycling, and repurposing of materials, aligning with the principles of the waste management hierarchy.
The linear economic model, as illustrated in
Figure 2a, relies on fundamentally unsustainable production processes due to their lack of material recycling and reuse. The future development and testing of products will focus heavily on the disassembly processes of adhesively bonded items, as depicted in
Figure 2, with a particular focus on enabling raw material recovery and facilitating repairs [
9]. Upcoming product designs should include explicit disassembly requirements in their specifications to align with the principles of the circular economy and contribute positively to the ecological balance. Thus, the selection of adhesives should not only be based on their mechanical and thermal properties, but also on their reusability, recyclability, and degradability without causing harm to the environment. Incorporating such considerations early in the design phase could lead to the development of innovative adhesive materials that are both robust for use and sustainable for end-of-life processing.
Numerous technologies are currently being developed to facilitate debonding-on-demand, employing various types of triggers (
Figure 3). In the Circular Bonding initiative, various debonding methods, namely, convection, induction, thermally expandable particles (TEPs), electricity, microwave, and ultraviolet treatments were examined for their application at a laboratory scale for smartphones [
10]. The use of electricity was by far the most cost-effective method [
11]. This cost-effectiveness is a critical factor in the commercial viability and potential widespread adoption of such technologies. Based on their ranking for each debonding technique across various criteria, induction, TEPs, electricity, and convection were identified as the only debonding technologies with market potential [
12]. These four techniques require heat generation, which is necessary for the thermal debonding process.
When incorporating thermal stimuli techniques into composite materials like GFRP, it is vital to thoroughly understand the thermal and structural performance of the materials. This ensures that the debonding process does not cause thermal and/or mechanical deformation, which could affect the structural integrity and durability of the GFRP joints under service conditions. Apart from the geometrical and mechanical properties of the adherend materials, the mechanical properties of the adhesive, such as ductility, stiffness, and toughness, play a significant role in stress distribution along the bondline. One of the key issues here is the occurrence of peak stresses at the overlap ends, particularly for brittle adhesives [
13]. These peak stresses can lead to premature joint failures. To address this issue, bi-adhesive (functionally graded) and flat-joggle-flat (FJF) joints have been identified as effective methods [
14,
15]. The common method of making functionally graded adhesive (FGA) joints is using two adhesives with different stiffness such as Hysol
® EA 9696 and 3M™ DP490 epoxy adhesives which have gained widespread acceptance by various industries in structural applications, with the former being particularly renowned for its exceptional toughness, an important characteristic for the adhesion of different materials, as indicated in [
16]. Conversely, DP490 has been characterized as a more brittle structural adhesive, as detailed in [
17]. Both adhesives exhibit high temperature resistance with distinct stiffness at both room and elevated temperatures. Silva et al. [
18] found that graded adhesive joints using both a low-temperature adhesive (LTA) and a high-temperature adhesive (HTA) exhibited improved performance, including increased load capacity, in comparison to using HTA at low temperatures or LTA at high temperatures individually. However, the mixed joints did not surpass the performance of LTA alone at low temperatures or HTA alone at high temperatures. Stein et al. [
19] introduced a novel computational method to optimize the distribution of thermal and mechanical stresses in FGA joints. The new analytical approach presented by these researchers is significant for evaluating stress in composite adhesive lap joints under different loading conditions. Their results, confirmed by FEA, demonstrate the framework’s ability to predict stress distributions and identify the most effective grading functions for lightweight FGA joint design.
To improve bonding and debonding characteristics of adhesive joints, fillers are commonly employed to functionalize the adhesive bondline. According to Gupta et al. [
20], neat epoxy typically exhibits explicit brittleness and poor fracture toughness, leading to catastrophic failure. The researchers go on to state that the incorporation of nanoparticles, also called fillers, can enhance both properties. In our previous studies, different approaches were employed, such as interleaving three carbon fiber veils to facilitate the Joule heating of the EA9696 structural film adhesive [
21]. In another study, thermally expandable particles (TEPs) with varying weight content were mixed with the DP490 two-part epoxy adhesive to induce stress concentration at high temperatures [
22]. Iron oxide (Fe
3O
4) particles were incorporated into the DP490 adhesive to enable the electromagnetic induction heating of dissimilar joints [
23].
An alternative approach to the fabrication of FGA joints involves the incorporation of filler materials within the bondline. This method allows for the tailoring of thermal and mechanical properties to meet the specific requirements of the joint. Incorporating fillers into the joints can enable them to outperform both the use of LTAs at low temperatures and that of HTAs at elevated temperatures. However, the challenge lies in the uniform distribution of these filler materials to avoid weak spots that could compromise the integrity of the joint. To address this challenge, recent studies have focused on the development of advanced mixing techniques that can ensure a more homogeneous dispersion of fillers within the adhesive. Jia et al. [
24] introduced graphene nanoplatelets (GNPs) into epoxy resin to create composite adhesives with increased modulus. They formed a bondline with a gradient in the modulus by placing high-modulus adhesives containing GNPs at the center and medium-modulus neat epoxy at the edges. FGA joints with a lateral size configuration of high-modulus edge region–medium-modulus center region–high-modulus edge region at a ratio of 1:3:1 demonstrate remarkable enhancements of 210.1%, 350%, and 1118.58% in failure load, elongation at break, and toughness, respectively. In their study, Kumar et al. [
25] observed that by adjusting the stiffness at the center and compliance at the edges of the adhesive joints, they achieved significant improvements in strength and toughness without compromising joint stiffness. Both the center stiffness- and edge compliance-tailored adhesive designs exhibited over 100% increase in strength and a 150% increase in toughness compared to non-tailored (with added stiffness/compliance) and constant-modulus counterparts. These improvements were attributed to a reduction in peak peel and shear strains at the ends of the adhesive. Considering these studies, stress analysis at elevated temperatures is as crucial as room temperature analysis for evaluating FRP joints bonded with various epoxy adhesives and fillers. Understanding the behavior of these materials under thermal stress is key to ensuring the reliability and safety of structures with these bonded joints. This approach to thermal stimuli debonding, combined with stress analysis across temperatures, marks an important step forward in adhesive technology and the use of composite materials.
Bending moment plays an important role in the analysis of adhesive joints, particularly in understanding stress distribution and failure mechanisms. An early study by Volkersen [
26] introduced the shear lag theory account for only shear stresses, prompting the development of more comprehensive stress analysis methods by Goland and Reissner [
27], which incorporate bending moments at the ends of the overlap region as boundary conditions. Redmann et al. [
28] demonstrated the significance of testing both single-lap joint samples and block shear samples. Block shear samples exhibited more than 100% strength compared to single-lap samples. This was because bending in the overlap region was eliminated during the tests, leading to pure shear. To address the issue of insufficient determination of the bending moment, Hart-Smith [
29] developed a novel bending moment factor. Zhao et al. [
30] proposed a new method for accurately determining the bending moments of joints of similar adherends, surpassing Goland and Reissner’s method and Hart-Smith’s method for overlap lengths up to 25 mm.
Timoshenko [
31] predicted that stress concentrations occur at the free edge of a bi-metal interface. Subsequent studies have similarly indicated that failure initiation tends to occur at the bi-material corner within the adhesive bondline, where singularities manifest as a locus of failure at room temperature [
32,
33] and under thermal loading [
34,
35]. Adhesive joints typically involve multiple materials, resulting in nonlinear behavior under applied loads. Finite element analysis (FEA) is a valuable numerical technique for predicting stress–strain distribution including the singularity at the adhesive/adherend interface, particularly in the presence of nonlinearities such as geometry, material properties, and boundary conditions.
This comprehensive study presents a critical examination of the behavior of adhesive joints under thermal and mechanical loading. Understanding the underlying mechanism is essential for the success of debonding-on-demand applications. This study considers the concept of on-demand debonding facilitated by thermomechanical stimuli, a technique that holds promise for applications requiring reversible adhesion. This aspect of the research explores a new frontier in adhesive technology, potentially paving the way for advancements in material design and structural engineering. By investigating how different thermal and mechanical loadings can cause different stress distributions, this study aims to provide insights into the performance and disassembly of bonded GFRP joints under a variety of temperature regimes, thereby contributing to the broader understanding of adhesive joint behavior in practical applications. Peel stress at the mid-section is more critical than shear stress as well as peel stress at the interface in determining SLJ shear strength, emphasizing the need to carefully analyze peel stress and the bending moment factor in SLJ design. This study also explores the analytical solutions proposed by Goland and Reissner, with modifications by Hart-Smith (HS) and Zhao, to understand the limitations and discrepancies when compared to FEA results. Through finite element analysis (FEA), we were able to obtain detailed mappings of the stress concentrations that occur both at the adhesive/adherend interface, as well as at the mid-section of the bondline thickness of the adhesive. This analysis enabled a clearer visualization of how both shear and peel stresses are distributed across these critical points in the bonding area. By mapping out stress distributions, the study emphasizes the importance of adhesive selection based on stress type, temperature, and solution methods in optimizing adhesive bonding applications. The insights deduced from this research aim to guide the development of modern adhesives, balancing high performance with environmental considerations, and contributing to the advancement of sustainable adhesive technologies.
Lastly, this study explores the concept of FGA joints and assess their effects on stress distribution. By adjusting the stiffness of these FGA joints using the experimental data from two structural epoxy adhesives, we computationally considered the incorporation of fillers at different weight percentages (wt.%). This comprehensive exploration of the FGA joints’ potential can reveal innovative strategies to significantly enhance the performance and extend the longevity in real-world applications. The exploration of FGA joints leads us to consider the adaptability of adhesive technologies in response to changing environmental conditions, potentially revolutionizing the way we approach the design and reusing of composite materials.
3. Results and Discussion
Figure 5,
Figure 6,
Figure 7 and
Figure 8 show the peel stress and shear stress computed by finite element analysis at two locations: the adhesive/adherend interface and the mid-section of the bondline thickness. The load applied to various joint configurations ranged from 40 N/mm to 300 N/mm; EA9696 and its combinations with CN34, CN10, and N34 were simulated at room temperature under a 300 N/mm load, while the same configurations at 100 °C were modeled at reduced loads of 150 N/mm for EA9696 and EA9696 + CN10, and 40 N/mm for EA9696 + CN34 and EA9696 + N34. These adjustments in the loading conditions were made in light of the findings presented in [
21], which detailed the performance outcomes of single-lap joint experiments.
The joint configurations using the DP490 adhesive and its variants with 5, 10, and 15 wt.% TEPs at room temperature conditions were modeled under a load of 240 N/mm, while the same configurations at 145 °C were simulated at a reduced load of 60 N/mm. This reduction in the applied load during high-temperature testing was based on the empirical results documented in [
22].
The presence of load eccentricity in the joint results in substantial rotations, leading to significant peak stresses at specific locations. In other words, the uneven distribution of the load causes pronounced rotations within the joint, causing high stress concentrations at those specific locations. The stress distributions are asymmetric at the interface, with a singularity at the square edge corner. At the mid-section, the stress distributions are relatively symmetric, with decreasing peak stresses. The conventional analytical models assume a symmetrical distribution of adhesive displacement with respect to the central axis. However, empirical observations reveal that displacement is greater on the load application side (actuator side) compared to the stationary grip side [
43]. Additionally, in high-temperature configurations, it is anticipated that the thermal mismatch between different materials will introduce additional stress singularities at the interface [
31].
It has been observed that shear stresses at the overlap edge are higher than at the center due to differential shearing. Peel stresses exhibit similar behavior. This phenomenon is a primary contributor to the delamination in FRP adherends. As reported in [
14,
15], reinforcing the adhesive layer leads to a reduction in the percentage of delamination of composites and the percentage of cohesive failure. The stress distribution for both neat DP490 and EA9696 joints at room temperature, as shown in
Figure 5 and
Figure 7, respectively, indicates that the peak peel stresses for both neat DP490 and EA9696 were higher than the remaining configurations. This observation further supports the benefits of reinforcing the adhesive layer in SLJs. Consistently, when we reinforced the adhesive layer with either TEPs or by interleaving carbon fiber veil, we observed an increase in the average percentage of adhesive residue and a decrease in the percentage of light-fiber-tearing (i.e., delamination) [
21,
22]. This highlights the benefits of reinforcing adhesive layers in terms of reducing delamination and improving overall joint performance, specifically improving static strength and extending fatigue life by delaying the onset of crack initiation associated with stress singularities [
44].
The analysis of
Figure 6 and
Figure 8 reveals a significant decrease in peak peel stresses at elevated temperatures across all joint configurations. For DP490 joint configurations, the peak peel stresses at the adhesive/adherend interface were reduced by 85.7% (0 wt.% TEPs), 75.5% (5 wt.% TEPs), 81.9% (10 wt.% TEPs), and 83.5% (15 wt.% TEPs). This reduction correlates with the absence of fiber-tearing observed in experimental tests, as reported in [
22]. In the case of EA9696 configurations, the results were more varied. Neat EA9696 and EA9696 + CN10 joints exhibited a reduction in peak peel stresses at the adhesive/adherend interface of 48.9% and 36.1%, respectively, which corresponded to a noticeable decrease in light-fiber-tearing, as reported in [
21]. More pronounced effects were observed in EA9696 + CN34 and EA9696 + N34 joints, where peak peel stresses at the adhesive/adherend interface dropped by 86.6% and 84.6%, respectively, resulting in no observable fiber-tearing.
The evaluation of
Figure 5 and
Figure 6 unveiled that the shear stress of the joints initially decreases and subsequently increases with the increment in TEPs wt.%. This behavior suggests that the shear stress is related to the shear modulus and the Young’s modulus of the adhesive, as detailed in
Table 4.
The ratios of peak peel stresses at the interface,
, to experimental shear strength,
[
22], for 0, 5, 10, and 15 wt.% TEPs at room temperature are 8.17, 5.79, 6.45, and 7.06, respectively. The ratios of peak peel stress at the bondline mid-section,
, to
[
22] for the same weight percentages of TEPs are consistently at 2.42, regardless of the amount of TEPs added. The ratios of peak shear stresses at the interface,
, to
[
22] for 0, 5, 10, and 15 wt.% TEPs at room temperature are 5.69, 5.74, 5.53, and 6.09, respectively. The ratios of peak shear stress at the bondline mid-section,
, to
[
22] for the same wt.% TEPs at room temperature are 3.06, 3.31, 3.22, and 3.18, respectively.
The ratios of peak peel stresses at the interface,
, to
[
22] for 0, 5, 10 and 15 wt.% TEPs at 145 °C are 8.71, 6.06, 8.71 and 8.71, respectively. The ratios of peak peel stress at the bondline mid-section,
, to
[
22] remain unchanged at 3.52 across all configurations for the different percentages of TEPs. The ratios of peak shear stresses at the interface,
, to
[
22] for 0, 5, 10, and 15 wt.% TEPs at 145 °C are 5.00, 4.96, 4.85, and 4.85, respectively. The ratios of peak shear stress at the bondline mid-section,
, to
[
22] for 0, 5, 10, and 15 wt.% TEPs at 145 °C are 2.95, 2.92, 2.99, and 2.92, respectively.
The observed variations in stress-to-strength ratios clearly indicate that the concentration of TEPs and temperature variations have a significant impact on the magnitude of stresses experienced. This intricate interplay suggests that adjusting TEP levels could be a key factor in tailoring the performance of adhesive joints. Notably, the stress-to-strength ratios measured at the bondline mid-section exhibited an upward trend with rising temperatures. This implies that the bondline mid-section has a more significant effect on the strength of the adhesive at higher temperatures. Compared to room temperature (
Figure 9), this observation also confirms that the rotation angle was lowered at 145 °C (
Figure 10) This evidence emphasizes the importance of monitoring temperature variations when determining the debonding temperature of the adhesive joints.
As the stiffness of the adherend increases, there is a corresponding decrease in the ratio of peak peel stress to peak shear stress,
[
45]. This trend indicates that stiffer adherends tend to redistribute stresses in a way that reduces the relative magnitude of peel stresses compared to shear stresses within the joint. In neat DP490 joints, the ratio of
was 2.06 at room temperature, decreasing to 1.74 at 145 °C. This decrease indicates that as temperature rises above the T
g of GFRP (135 °C), GFRP loses its stiffness significantly, as detailed in
Table 1. Conversely, neat EA9696 joints demonstrated a
ratio of 1.48 at room temperature, which slightly increased to 1.52 at 100 °C—below the T
g of GFRP, thus maintaining its stiffness and showing no significant change in the ratio. Although the adherends were 32% thicker in the case of neat DP490 compared to the neat EA9696 joints, and hence stiffer, the DP490 joints exhibited a higher ratio of
because DP490 is stiffer than EA9696 under room temperature conditions (see
Table 4 and
Table 5). Compared to the EA9696 adhesive joints, the DP490 adhesive joints exhibited lower peel and shear stresses because of the lower Young’s modulus of DP490 and using thicker adherends which improves the load-bearing capacity of the SLJ [
46]. However, using thicker adherends increases both cost and the weight of structures. A comparison of tests conducted at room temperature to those performed at elevated temperatures revealed a substantial decrease in both peel and shear stresses experienced by the joints. When operating at temperatures above T
g, such as when using the DP490 adhesive, the adhesive properties underwent a significant decrease, leading to a heightened sensitivity of the joint’s behavior to even minor temperature variations within this temperature range. This significant reduction in stress levels can be attributed to several influential factors. Firstly, the lower Young’s modulus of the materials employed in the joints plays a crucial role in reducing the stresses. Moreover, the diminished levels of joint rotation observed under the different conditions also help reduce stress. Additionally, it is worth emphasizing that the applied load during these tests was notably lower (see
Table 4 and
Table 5), which further contributes to the overall decrease in stress levels within the joints. An increase in temperature has two effects: it reduces the yield stress of the adherend material and leads to a reduction in the failure envelope. Consequently, as the temperature rises, the stress required for failure decreases due to the decreased yield stress and the resulting shrinkage in the failure envelope [
47].
The aforementioned findings indicate that peel stress plays a dominant role in accurately determining the experimental strength of SLJ. The analytical prediction of peel stress depends on the precise calculation of the bending moment factor. Thus, three different methods for determining bending moment factors were considered to identify the most accurate one.
Figure 9,
Figure 10,
Figure 11 and
Figure 12 illustrate comparisons of the peel stress and shear stress distributions in identical joints using both analytical and FEA solutions. This study evaluates three analytical methods—Goland and Reissner’s, Hart-Smith’s, and Zhao’s—against the FEA results for different adhesive types and temperature conditions. When comparing the outcomes for 10 wt.% TEP–epoxy adhesive joints analyzed through Goland and Reissner’s technique, Hart-Smith’s procedure resulted in an 8.5% increase in peak peel stress and a 2.9% rise in peak shear stress, whereas Zhao’s method led to a decrease of 0.9% in peak peel stress and 0.3% in shear stress under room temperature conditions, as shown in
Figure 9a,b. Furthermore, in comparison to the findings for 10 wt.% TEP–epoxy adhesive joints using Goland and Reissner’s method, Hart-Smith’s approach demonstrated an 11.1% increase in peak peel stress and a 0.6% rise in peak shear stress, while Zhao’s method exhibited a reduction of 0.4% in peak peel stress and 0.02% in peak shear stress at a temperature of 145 °C (
Figure 10a,b).
Compared to the findings for CN34 epoxy adhesive joints using Goland and Reissner’s method, Hart-Smith’s approach led to a 4.2% increase in peak peel stress and a 1.8% increase in peak shear stress, while Zhao’s approach resulted in a decrease of 3.3% in peak peel stress and 1.5% in peak shear stress under room temperature conditions, as illustrated in
Figure 11a,b. At an elevated temperature of 100 °C, Hart-Smith’s approach led to a 16.9% increase in peak peel stress and a 9.4% increase in peak shear stress, while Zhao’s approach resulted in a decrease of 0.3% in peak peel stress and 0.1% in peak shear stress, as illustrated in
Figure 12a,b.
The alternative methods proposed by Hart-Smith and Zhao to determine the bending moment factor showed some improvements compared to the original Goland and Reissner model. However, these modifications had a marginal impact on stress due to the neglect of transverse deflection in the overlap region, which led to an underestimation of the bending behavior [
41]. The analytical solutions generally aligned well with the FEA results at the adhesive joint’s mid-section. However, the FEA revealed significant stress differences between the interface and mid-section of the adhesive bondline, particularly at the corner interface where stress singularity was observed. It should be noted that this is not typically taken into account in the design of adhesive joints [
32]. The stress concentration at the corners can be attributed to the assumption of perfect elasticity in the solutions. The results from the Goland and Reissner solutions, including the modified approaches by Hart-Smith and Zhao, showed good agreement with peel stress around the central overlap length obtained from the FEA. However, there was a significant difference in the maximum shear stress. When subjected to different loads at a temperature of 100 °C, the disparity between the peak shear stress predicted by the Goland and Reissner solutions and the FEA models increased. It is important to note that the analyses conducted in this study did not consider bond damage behavior, such as the use of the cohesive zone modeling (CZM) approach. The thickness of the GFRP adherend in a single-lap joint significantly influences stress distribution and the behavior of the joint under bending moment and load eccentricity. When the thickness of the adherend is increased, the bending stress decreases due to a reduction in the applied bending moment. Zhao’s [
30] study indicated that thicker adherends (up to 6.35 mm) exhibit higher stiffness, resulting in reduced bending compared to thinner adherends, leading to more accurate results for thicker adherends. The elastic solutions proposed by Goland and Reissner, as well as by Hart-Smith, exhibit a considerable degree of overlap for both adhesives. However, it is important to note that these models tend to underestimate the results from the FEA across different adhesives and joint configurations. Still, the results using the methods of Goland and Reissner, Hart-Smith, and Zhao will be invaluable for the analysis of failure criteria, especially since the thickness of the adhesive layer (
ta) is 0.5 mm and the ratio of overlap length to adhesive layer thickness (c/t
a) ≥ 20, in accordance with the findings from both studies [
48,
49].
As demonstrated in
Figure 13a, for DP490 adhesive joints at room temperature (RT) and under peel stress, the parabolic solution yielded a stress of 45.2 MPa, while the biquadratic solution resulted in a slightly higher value of 52.0 MPa, indicating a difference of −15.38%. Similarly, for the shear stress type at room temperature, the parabolic solution produced a stress of 30.8 MPa, and the biquadratic solution yielded 36.1 MPa, showing a difference of −14.15%. This trend continued to be observed for the cases with an elevated temperature of 145 °C, as indicated in
Figure 13b. In the peel stress type, the parabolic solution yielded a stress of 2.86 MPa, whereas the biquadratic solution provided a higher stress of 4.46 MPa, resulting in a difference of −55.15%. In the shear stress type at the same elevated temperature, the parabolic solution yielded a stress of 2.49 MPa, while the biquadratic solution resulted in a slightly greater stress of 3.47 MPa, demonstrating a difference of −28.17%. Comparing these results, it is evident that at both room temperature and elevated temperature, the biquadratic solution consistently provided higher stress values compared to the parabolic solution for both peel and shear stress types. The percentage differences between the two solutions indicate the extent to which these stress values diverged, with the largest difference observed at 145 °C in the peel stress type.
As shown in
Figure 14a, for EA9696 adhesive joints at RT and under peel stress, the parabolic solution exhibited a stress measurement of 53.3 MPa, whereas the biquadratic solution yielded a notably higher value of 73.9 MPa, indicating a substantial variance of −38.91%. Likewise, in the shear stress category at room temperature, the parabolic solution resulted in a stress reading of 35.5 MPa, while the biquadratic solution yielded a significantly higher value of 51.3 MPa, illustrating a considerable difference of −44.45%. Upon increasing the temperature to 100 °C, as shown in
Figure 14b, the comparative trends persisted. For the peel stress type, the parabolic solution yielded a stress level of 13.6 MPa, while the biquadratic solution resulted in a significantly higher value of 38.9 MPa, leading to a notable discrepancy of −64.91%. Similarly, within the shear stress category at the same temperature, the parabolic solution yielded a stress measurement of 7.70 MPa, whereas the biquadratic solution resulted in a significantly higher value of 26.3 MPa, indicating a substantial variance of −70.68%. Comparing the results at both room temperature and elevated temperature, the biquadratic solution consistently provided higher stress values compared to the parabolic solution for both peel and shear stress types.
Overall, the comparison indicates that both adhesives, DP490 and EA9696, displayed different behaviors in response to stress and temperature, as presented in
Figure 13 and
Figure 14. While EA9696 consistently exhibited higher stress values, particularly under the biquadratic solution, DP490 generally showed smaller differences between solution methods. The results suggest that the parabolic solution is beneficial in reducing peak stresses at room and elevated temperatures because the parabolic solution creates a relatively more gradual change in Young’s modulus distribution near the overlap [
42]. This smoother transition in material properties mitigates the stress concentrations in these critical areas. Conversely, the biquadratic solution led to higher peak stresses. While this may initially appear to be a drawback, it presents a notable advantage in easier crack initiations, leading to easier debonding at elevated temperature. These observations underscore the importance of adhesive selection based on the specific stress type, temperature conditions, and solution method to ensure accurate and reliable results in adhesive bonding but also to optimize the debonding process.