Study on Establish a Brittle Fracture Prediction Considering Different Crack Opening Modes Using Mixed-Mode Ratio

: In this study, we propose a method for predicting the occurrence of brittle fractures in the beam-to-column joints of steel structures, considering different crack opening modes. We conducted experiments on beam-to-diaphragm joint specimens with varying plastically constrained cracks to reproduce brittle fractures. The experiments’ results demonstrated the effectiveness of the toughness scale model and the Weibull stress approach. In addition, we propose the mixed-mode ratio, which is a quantitative index of the mode difference, and we applied it to the ﬁnite element models of the specimens. In this study, we evaluate the validity of the mixed-mode ratio and explore the differences in crack opening modes, as they pertain to the occurrence of brittle fractures.


Introduction
Brittle fractures in the welded joints of building structures were first observed in the devastating 1994 Northridge earthquake [1]. In Japan, the 1995 Kobe earthquake, which registered a magnitude of 7.2 on the Richter scale, resulted in significant damage to buildings due to brittle fracture. These fractures resulted from welding defects or cope holes in the webs on the bottom flange side of the welds [2]. After the earthquakes, various experimental studies, such as those conducted by the Architectural Institute of Japan [3], proved that brittle fractures occur owing to ductile cracks propagation by singularities in the shapes around the bottoms of cope holes and end tabs.
Several indices, such as J-value and crack tip opening displacement (CTOD) [4], which are used for the evaluation of brittle fractures, allow an overly safe assessment of building structures. This is because brittle fractures occur as a result of significant deformations of building structures, the stress-strain state at the crack tip constantly changes owing to ductile crack propagation, whereas the plastic constraint state at the crack tip of the material specimen used to obtain the fracture toughness significantly differs from that of the actual building structure. There currently exist various methods focused on the differences in plastic constraints, which are useful for the prediction and evaluation of brittle fractures. For example, the Japan Welding Engineering Society [5] proposed WES2808, which is a method for correcting the fracture driving force using the plastic constraint correction factor and the temperature shift concept. The applicability of their proposed method also expands to the CTOD design curve. Meanwhile, we demonstrated the effectiveness of the toughness scale model (TSM) [6] and Weibull stress approach [7], which is an extension of the Ritchie-Knott-Rice model [8], for the effective prediction of the failure of buildings based on previous studies on the occurrence of brittle fractures in beam-to-diaphragm joint specimens [9]. A study on the estimation of plastic deformation capacity [10] showed that the Weibull stress approach could be used to account for the effects of plastic constraints and cyclic loading. There is interest in the reliability of cracks, and not just in building structures. Coro et al. proposed a new model for structures made by additive manufacturing by extending the method developed for welded components [11]. Regarding fatigue crack propagation, a unique code that includes crack propagation was proposed for gas turbine welded structures [12].
In addition, other factors, such as different crack opening modes, affect the onset of brittle fractures. The defects at the ends of welds in building structures are affected by complex stresses, and the modes mentioned above act in a combinatorial manner [13]. Crack deformations can be classified into three modes, as follows: opening mode (Mode I), in-plane shear (Mode II), and out-of-plane shear (Mode III). As the current Weibull stress approaches assume the dominance of Mode I, they cannot be used to determine the influence of Modes II and III [14]. Thus, the development of an approach that ensures the accurate evaluation of fracture initiation considering different crack deformation modes would play a significant role in improving our understanding of the risk of brittle fracture occurrence in building structures.
In this study, we focused on the different modes acting on weld defects and we propose a mixed-mode ratio as a quantitative index for the ratios of Mode I to Mode II. The mixed-mode ratio may be inaccurate when a ductile crack propagates significantly. We targeted ductile crack growth of up to 20 mm in this study. Furthermore, we classified the brittle fracture results of several experiments using the mixed-mode ratio as the index and subsequently investigated its relationship with plastic deformation.

Fracture Tests
We conducted fracture tests on beam-to-diaphragm joint specimens that simulate full-penetration weld joints to reproduce the brittle fracture of a building structure after significant deformation. The experiments were intended to replicate the phenomenon of ductile crack growth from the end of the full-penetration weld between the flange and the diaphragm, followed by a brittle fracture. This study aimed to reproduce the experimental results through finite element analysis, achieve an effective prediction of the occurrence of brittle fractures through the TSM and Weibull stress approaches, and investigate the effect of different modes on the prediction accuracy. The variables used in the experiments included the location, shape, depth, and inserted tip radius of the crack. Among the several specimens used in the experiments, six specimens showed the occurrence of brittle fractures.

Specimens and Experimental Variables
The specimens used in the experiments included an I-section beam connected to rectangular hollow section (RHS) column joints through a diaphragm, as shown in Figure 1.
Appl. Mech. 2021, 2, FOR PEER REVIEW [10] showed that the Weibull stress approach could be used to account for the effec plastic constraints and cyclic loading. There is interest in the reliability of cracks, an just in building structures. Coro et al. proposed a new model for structures made b ditive manufacturing by extending the method developed for welded components Regarding fatigue crack propagation, a unique code that includes crack propagation proposed for gas turbine welded structures [12].
In addition, other factors, such as different crack opening modes, affect the ons brittle fractures. The defects at the ends of welds in building structures are affecte complex stresses, and the modes mentioned above act in a combinatorial manner Crack deformations can be classified into three modes, as follows: opening mode (M I), in-plane shear (Mode II), and out-of-plane shear (Mode III). As the current We stress approaches assume the dominance of Mode I, they cannot be used to determin influence of Modes II and III [14]. Thus, the development of an approach that ensure accurate evaluation of fracture initiation considering different crack deformation m would play a significant role in improving our understanding of the risk of brittle fra occurrence in building structures.
In this study, we focused on the different modes acting on weld defects and we pose a mixed-mode ratio as a quantitative index for the ratios of Mode I to Mode II mixed-mode ratio may be inaccurate when a ductile crack propagates significantly targeted ductile crack growth of up to 20 mm in this study. Furthermore, we classifie brittle fracture results of several experiments using the mixed-mode ratio as the index subsequently investigated its relationship with plastic deformation.

Fracture Tests
We conducted fracture tests on beam-to-diaphragm joint specimens that sim full-penetration weld joints to reproduce the brittle fracture of a building structure significant deformation. The experiments were intended to replicate the phenomen ductile crack growth from the end of the full-penetration weld between the flange an diaphragm, followed by a brittle fracture. This study aimed to reproduce the experim results through finite element analysis, achieve an effective prediction of the occurr of brittle fractures through the TSM and Weibull stress approaches, and investigat effect of different modes on the prediction accuracy. The variables used in the experim included the location, shape, depth, and inserted tip radius of the crack. Among the eral specimens used in the experiments, six specimens showed the occurrence of b fractures.

Specimens and Experimental Variables
The specimens used in the experiments included an I-section beam connected to tangular hollow section (RHS) column joints through a diaphragm, as shown in Figu  The bottom flange, diaphragm, and end of the full-penetration weld, which is shown in red in Figure 1, were simulated to create a beam-to-diaphragm joint specimen, which is shown in Figure 2.
The bottom flange, diaphragm, and end of the full-penetration weld, which is sh in red in Figure 1, were simulated to create a beam-to-diaphragm joint specimen, w is shown in Figure 2. The flange was 25 mm thick and 100 mm wide, the diaphragm was 32 mm thick 225 mm wide, and the web was 16 mm thick and 150 mm long. All the plates were m ufactured using SN490B, in accordance with the specifications provided by the Japan Industrial Standards (JIS). All the welding assemblies were completed in the factory. bottom flange was welded to the diaphragm through full-penetration welding u YGW18 in accordance with JIS. The welding process comprised seven layers and 10 pa The heat input of each layer was below 40 kJ/cm, and the temperature between passes under 350 °C.
The experimental variables included insertion location, shape, depth, and tip ra of the defect inserted at the beginning and at the end of the weld. Figure 3 shows the de types. Table 1 shows a list of the experimental variables.  We considered two defect insertion locations: the diaphragm side and the flange s We considered two defect shapes: through crack (TC) and surface crack (SC). The dista from the face of the flange to the crack tip represented the defect depth. The different c The flange was 25 mm thick and 100 mm wide, the diaphragm was 32 mm thick and 225 mm wide, and the web was 16 mm thick and 150 mm long. All the plates were manufactured using SN490B, in accordance with the specifications provided by the Japanese Industrial Standards (JIS). All the welding assemblies were completed in the factory. The bottom flange was welded to the diaphragm through full-penetration welding using YGW18 in accordance with JIS. The welding process comprised seven layers and 10 passes. The heat input of each layer was below 40 kJ/cm, and the temperature between passes was under 350 • C.
The experimental variables included insertion location, shape, depth, and tip radius of the defect inserted at the beginning and at the end of the weld. Figure 3 shows the defect types. Table 1 shows a list of the experimental variables.
Appl. Mech. 2021, 2, FOR PEER REVIEW 3 The bottom flange, diaphragm, and end of the full-penetration weld, which is shown in red in Figure 1, were simulated to create a beam-to-diaphragm joint specimen, which is shown in Figure 2. The flange was 25 mm thick and 100 mm wide, the diaphragm was 32 mm thick and 225 mm wide, and the web was 16 mm thick and 150 mm long. All the plates were manufactured using SN490B, in accordance with the specifications provided by the Japanese Industrial Standards (JIS). All the welding assemblies were completed in the factory. The bottom flange was welded to the diaphragm through full-penetration welding using YGW18 in accordance with JIS. The welding process comprised seven layers and 10 passes. The heat input of each layer was below 40 kJ/cm, and the temperature between passes was under 350 °C.
The experimental variables included insertion location, shape, depth, and tip radius of the defect inserted at the beginning and at the end of the weld. Figure 3 shows the defect types. Table 1 shows a list of the experimental variables.  We considered two defect insertion locations: the diaphragm side and the flange side. We considered two defect shapes: through crack (TC) and surface crack (SC). The distance from the face of the flange to the crack tip represented the defect depth. The different crack  We considered two defect insertion locations: the diaphragm side and the flange side. We considered two defect shapes: through crack (TC) and surface crack (SC). The distance from the face of the flange to the crack tip represented the defect depth. The different crack tip radii of the defects were produced using two methods: copperplate insertion and mechanical notching. Based on the result after the brittle fracture, the crack tip radius obtained through a copper plate exhibited a sharp tip radius of approximately 0 mm. Figure 4 shows the copper plate insertion of the specimen. tip radii of the defects were produced using two methods: copperplate insertion and mechanical notching. Based on the result after the brittle fracture, the crack tip radius obtained through a copper plate exhibited a sharp tip radius of approximately 0 mm. Figure  4 shows the copper plate insertion of the specimen. The microscopic measurement of the mechanical notch defect confirmed the bluntness of the crack tip with a radius of approximately 0.15 mm. Figure 5 shows a magnified image of the notch tip.

Material Properties
In this study, we focused on ductile cracks propagating from the tip of the defect to the brittle fracture after significant deformation. To reproduce the brittle fracture initiation after ductile crack propagation, the fracture toughness of each specimen was determined through Charpy impact tests. The experimental temperature of the beam-to-diaphragm joint specimens was investigated based on the energy transition temperature. Table 2 shows the results of the Charpy impact test. The positions of the test pieces are the base metal of the flange plate, diaphragm plate, heat-affected zone (HAZ), and deposited weld metal (DEPO). The microscopic measurement of the mechanical notch defect confirmed the bluntness of the crack tip with a radius of approximately 0.15 mm. Figure 5 shows a magnified image of the notch tip.
tip radii of the defects were produced using two methods: copperplate insertion and mechanical notching. Based on the result after the brittle fracture, the crack tip radius obtained through a copper plate exhibited a sharp tip radius of approximately 0 mm. Figure  4 shows the copper plate insertion of the specimen. The microscopic measurement of the mechanical notch defect confirmed the bluntness of the crack tip with a radius of approximately 0.15 mm. Figure 5 shows a magnified image of the notch tip.

Material Properties
In this study, we focused on ductile cracks propagating from the tip of the defect to the brittle fracture after significant deformation. To reproduce the brittle fracture initiation after ductile crack propagation, the fracture toughness of each specimen was determined through Charpy impact tests. The experimental temperature of the beam-to-diaphragm joint specimens was investigated based on the energy transition temperature. Table 2 shows the results of the Charpy impact test. The positions of the test pieces are the base metal of the flange plate, diaphragm plate, heat-affected zone (HAZ), and deposited weld metal (DEPO).

Material Properties
In this study, we focused on ductile cracks propagating from the tip of the defect to the brittle fracture after significant deformation. To reproduce the brittle fracture initiation after ductile crack propagation, the fracture toughness of each specimen was determined through Charpy impact tests. The experimental temperature of the beam-to-diaphragm joint specimens was investigated based on the energy transition temperature. Table 2 shows the results of the Charpy impact test. The positions of the test pieces are the base metal of the flange plate, diaphragm plate, heat-affected zone (HAZ), and deposited weld metal (DEPO).
For the FTC and FSC specimens, the experimental temperature was set to −20 • C, which is slightly above the transition temperature of the welded metal, because the transition temperature of the flange was high. Although high transition temperatures were recorded at the diaphragms of the DTC and DSC specimens, the brittle fracture was expected to occur in the welded metal or in the HAZ. Therefore, the experimental temperature was set to −45 • C, which was between that of the two specimens.

Specimen
Section The tensile coupon tests were conducted using the test pieces cooled to a temperature similar to that used in the experiments. Table 3 shows the results of the tensile coupon tests. Jc, which is the critical J value, was determined to be 152.46 N/mm through the single-edge notched bending (SENB) tests of the FTC and FSC specimens.

Testing Procedures
The diaphragm was fixed to the base using high-tension bolts and was loaded quasistatically using hydraulic jacks. Figure 6 shows the installation and loading procedures of the test specimens. The bending moment M and the relative rotation angle θ were calculated by: The bending moment M and the relative rotation angle θ were calculated by: where L is the distance from the fixed end to the loading position, and u 1 and u 2 are the values measured using a high-sensitivity displacement meter at the position shown in Figure 6.
The tests were conducted at a ram head displacement rate of 0.05 mm/s. Cyclic loading comprised a few cycles in the elastic range. Subsequently, cantilever rotations of θ p , 2θ p , and 3θ p with positive (tension) and negative (compression) displacement at each amplitude were conducted until failure. The full plastic moment M p was calculated using the measured yield strengths of the materials. The rotations at the full plastic moment θ p were calculated by dividing M p by the elastic stiffness of the cantilever, which was determined using the slopes at the unloading portions of the hysteresis loops.

Test Results
Brittle fractures occurred in all specimens. Table 4 shows the experimental results of the strength and deformation of the specimens. We used the cumulative deformation factor η to indicate the plastic deformation capacity (see Appendix A for additional information on η).

Failure Moment Prediction
To numerically reproduce the experimental results, we performed finite element (FE) analysis using the general-purpose analysis software called ABAQUS (Ver. 2021, Dassault Systèmes, Paris, French). The elements were eight-node reduced integration elements following the von Mises yield condition. The hardening law follows isotropic hardening, and the Poisson's ratio was set to 0.3. For the material definition, the tensile test results were

Failure Moment Prediction
To numerically reproduce the experimental results, we performed finite element (FE) analysis using the general-purpose analysis software called ABAQUS (Ver. 2021, Dassault Systèmes, Paris, French). The elements were eight-node reduced integration elements following the von Mises yield condition. The hardening law follows isotropic hardening, and the Poisson's ratio was set to 0.3. For the material definition, the tensile test results were converted to true stress-logarithmic strains, and 13 points were used to approximate the data. The model was fabricated as a half model, and the minimum element size was set to 0.05 mm at the defect tip. The defect shape was determined by evaluating the fracture surface after the experiment. Monotonic loading was followed. This analysis did not consider the propagation of the ductile crack. An example of the analytical model is shown in Figure 9.
Appl. Mech. 2021, 2, FOR PEER REVIEW 8 converted to true stress-logarithmic strains, and 13 points were used to approximate the data. The model was fabricated as a half model, and the minimum element size was set to 0.05 mm at the defect tip. The defect shape was determined by evaluating the fracture surface after the experiment. Monotonic loading was followed. This analysis did not consider the propagation of the ductile crack. An example of the analytical model is shown in Figure 9.  Figure 10 shows the comparison between the experimental and the analysis results As the analysis was based on one-sided loading, the experimental results involved skeleta curves, with the history curve replaced by monotonic loading. The experimental and analytical results corresponded well. Therefore, the analysis captured the macroscopic be-  Figure 10 shows the comparison between the experimental and the analysis results. As the analysis was based on one-sided loading, the experimental results involved skeletal curves, with the history curve replaced by monotonic loading. The experimental and analytical results corresponded well. Therefore, the analysis captured the macroscopic behavior of the specimens during the experiments.  Figure 10 shows the comparison between the experimental and the analysis results. As the analysis was based on one-sided loading, the experimental results involved skeletal curves, with the history curve replaced by monotonic loading. The experimental and analytical results corresponded well. Therefore, the analysis captured the macroscopic behavior of the specimens during the experiments.

TSM Approach
Under small-yield conditions, the fracture toughness at the crack tip can be determined univocally using the J-integral. However, as the deformation of the member increases and the plasticity increases, the fracture limit of the material cannot be determined

TSM Approach
Under small-yield conditions, the fracture toughness at the crack tip can be determined univocally using the J-integral. However, as the deformation of the member increases and the plasticity increases, the fracture limit of the material cannot be determined using the J-integral due to a plastic constraint. The SENB specimens used to obtain J were subjected to a more significant plastic constraint at the crack tip. Anderson and Dodds [6] focused on the extent of this high-stress region. They proposed that for specimens under different plastic constraints, the probability of the occurrence of brittle fractures is equal when the volumes of the high-stress regions at the crack tips are equal.
In this study, the apparent fracture toughness value app J c was determined by increasing the J value through the high-stress region. The flat area of the crack tip bounded by a constant maximum principal stress contour replaced the volume of the highly stressed region. The maximum principal stress was assumed to be σ 1 = 3σ Y . Figure 11 shows an example of the determination of app J c . The calculated fracture moment for each specimen was the bending moment when J reached app J c , not J c .
Appl. Mech. 2021, 2, FOR PEER REVIEW 9 using the J-integral due to a plastic constraint. The SENB specimens used to obtain J were subjected to a more significant plastic constraint at the crack tip. Anderson and Dodds [6] focused on the extent of this high-stress region. They proposed that for specimens under different plastic constraints, the probability of the occurrence of brittle fractures is equal when the volumes of the high-stress regions at the crack tips are equal. In this study, the apparent fracture toughness value appJc was determined by increasing the J value through the high-stress region. The flat area of the crack tip bounded by a constant maximum principal stress contour replaced the volume of the highly stressed region. The maximum principal stress was assumed to be σ1 = 3σY. Figure 11 shows an example of the determination of appJc. The calculated fracture moment for each specimen was the bending moment when J reached appJc, not Jc.

Weibull Stress Approach
A significant variation in fracture toughness was observed, especially, when transitioning from ductile cracks to brittle fractures. This variation can be attributed to the factors that cause brittle fractures at the crack tip. Beremin [7] proposed a local approach for

Weibull Stress Approach
A significant variation in fracture toughness was observed, especially, when transitioning from ductile cracks to brittle fractures. This variation can be attributed to the factors that cause brittle fractures at the crack tip. Beremin [7] proposed a local approach for evaluating the variations in fracture toughness using a statistical method. Based on this approach, brittle fractures occur when a microcrack is subjected to a specific effective stress, in which the limit of the weakest microcrack is regarded as the limit of the entire material. In other words, brittle fractures are considered to follow the weakest link model. Because the failure probability of the weakest link model follows a Weibull distribution, the probability of brittle fracture occurrence is expressed as: where σ W denotes the Weibull stress, σ u denotes the scale parameter corresponding to the Weibull stress with a cumulative probability of failure that is equal to 63.2%, and m is the Weibull slope which is a shape parameter that characterizes the scatter of the Weibull stress.
For structural steel, m ranges from 10 to 50 [15]. Typically, m has a low value in the case of low-level, broadly distributed fracture toughness. In this study, the fracture process zone was the region with a maximum principal stress at the crack tip of 3σ Y or greater. This parameter can be expressed as follows: where σ eff denotes the effective stress, V p denotes the fracture process zone, and V 0 is a unit volume, set to 1 mm 3 . The maximum likelihood estimation proposed by Riesch-Oppermann and Diegele [16] was used to calculate m. In this study, we used the maximum principal stress as the effective stress. The fracture process zone was determined as the region in which the maximum principal stress exceeded 3σ Y or greater. Based on the results obtained through several SENB tests, m was determined to be 32. The critical value of the Weibull stress, σ W,cr = 1324 N/mm 2 , was determined as the Weibull stress at J c = 152.46 N/mm in the analysis. Figure 12 shows the relationship between the Weibull stress and the bending moment for each specimen. The calculated failure moment was achieved when the Weibull stress reached the limit determined through elastic-plastic analysis.

Appl. Mech. 2021, 2, FOR PEER REVIEW
The maximum likelihood estimation proposed by Riesch-Oppermann and Die [16] was used to calculate m. In this study, we used the maximum principal stress as effective stress. The fracture process zone was determined as the region in which the m imum principal stress exceeded 3σY or greater. Based on the results obtained through eral SENB tests, m was determined to be 32. The critical value of the Weibull stress, σW 1324 N/mm 2 , was determined as the Weibull stress at Jc = 152.46 N/mm in the analysi Figure 12 shows the relationship between the Weibull stress and the bending ment for each specimen. The calculated failure moment was achieved when the Wei stress reached the limit determined through elastic-plastic analysis.

Prediction Results
The predicted failure moments of the beam-to-diaphragm joint specimens were culated according to the procedure described in Sections 3.2 and 3.3. The results are s marized in Table 5, where TSMMpre and WMpre represent the predicted bending momen failure calculated using the TSM approach and the Weibull approach, respectively.

Prediction Results
The predicted failure moments of the beam-to-diaphragm joint specimens were calculated according to the procedure described in Sections 3.2 and 3.3. The results are summarized in Table 5, where TSM M pre and W M pre represent the predicted bending moments of failure calculated using the TSM approach and the Weibull approach, respectively. Based on the similar magnitudes of the predicted bending moments, the results of the fracture prediction methods were verified. For the FTC-8C, FSC-7C, and FSC-14C defects with sharp crack tips, both approaches accurately predicted the occurrence of fractures. Both approaches also demonstrated high accuracy in predicting the fracture behavior of FSC-7C, in which brittle fractures occurred after significantly high plastic deformations (η = 24), suggesting that both approaches can be used to effectively predict brittle fractures after a significant deformation of building structures. In contrast, low prediction accuracy was noted for the DTC-8W, DTC-12W, and DSC-21W specimens. This can be attributed to the crack tips and relaxed plastic constraint of these specimens. Thus, this contradicts the idea that the Weibull stress approach is unaffected by plastic constraints. In addition to the differences in the crack tip radii, prediction accuracy could also be affected by the defect location. The location and shape of the defect also affect the plastic constraint; however, the crack opening mode is also different. The following section will focus on differences in the crack opening mode and discuss their impact on prediction accuracy.

Effect of the Crack Opening Modes
Different crack opening modes affect the crack tip depending on the type, depth, and location of a crack. Although various modes act in combination in a building structure, only Mode I is involved in the SENB test when obtaining the fracture toughness value. As Mode III, which is an out-of-plane shear type fracture, has minimal effect of the occurrence of brittle fractures, this study focused on Modes I and II only and proposes a mixedmode ratio R I as the first step to capture the effect of different modes on the initiation of brittle fractures.

Mixed-Mode Ratio Methodology
In their study on the effect of different crack opening modes through material testing, Shimizu et al. [14] defined the ratio of Mode II to Mode I using the energy release rate in the initial deformation stage. In this study, we considered the ratio of Mode II to Mode I using the stress intensity factor. The mixed-mode ratio R I is defined as follows: where K I and K II represent the stress intensity factors for Mode I and Mode II, respectively, and R I represents the slope of K I and K II , which has a maximum value of 1.0 and decreases with the increase in Mode II.

Elastic Analysis of the Stress Intensity Factor
The stress intensity factor is the maximum value on the circumferential integral path at 0.2 mm from the crack tip. The stress intensity factor was defined for small-scale yielding conditions and is calculated for linear isotropic materials using ABAQUS. The direction of crack propagation is considered as the direction of the maximum energy release rate. The analytical model and other conditions for the elastic analysis of the stress intensity factor were similar to those of the elastic-plastic analysis. Figure 13 shows an example of the comparison between the elastic and the elasticplastic analyses. Although we observed a gap between the results after yielding, the transitions of the elastic and elastic-plastic analyses appeared almost identical in the initial stage of the deformation. Therefore, the K I values and K II at θ = 0.01 rad, which were linearly correlated, were adopted to calculate R I .

Relationship between Mixed-Modes and Fracture Toughness
The effect of the mixed-mode ratio, which quantifies the difference in the crack opening modes, on fracture toughness was investigated. The mixed-mode ratios for each specimen are listed in Table 6. To compare the crack tip and mode differences, the results of FTC-8W, FSC-7W, and FSC-14W are also listed in Table 6, in addition to the six tested specimens. Comparing FTC-8C and FTC-8W, FSC-7C and FSC-7W, and FSC-14C and FSC-14W, it appears that RI was not affected by the differences in the crack tip radius. Meanwhile, the defect depth caused a slight difference, such as that between FSC-7C and FSC-14C. As for the defect position, RI decreased for defects on the diaphragm side, as observed for FTC-8W and DTC-8W. This result is consistent with the assumption that the defects on the diaphragm side, which are horizontal to the loading direction, were subjected to stronger shear deformation than those on the flange side. For FTC-8C and FSC-7C, the difference in the defect shape had a relatively considerable effect on RI. In both cases, the

Relationship between Mixed-Modes and Fracture Toughness
The effect of the mixed-mode ratio, which quantifies the difference in the crack opening modes, on fracture toughness was investigated. The mixed-mode ratios for each specimen are listed in Table 6. To compare the crack tip and mode differences, the results of FTC-8W, FSC-7W, and FSC-14W are also listed in Table 6, in addition to the six tested specimens. Comparing FTC-8C and FTC-8W, FSC-7C and FSC-7W, and FSC-14C and FSC-14W, it appears that R I was not affected by the differences in the crack tip radius. Meanwhile, the defect depth caused a slight difference, such as that between FSC-7C and FSC-14C. As for the defect position, R I decreased for defects on the diaphragm side, as observed for FTC-8W and DTC-8W. This result is consistent with the assumption that the defects on the diaphragm side, which are horizontal to the loading direction, were subjected to stronger shear deformation than those on the flange side. For FTC-8C and FSC-7C, the difference in the defect shape had a relatively considerable effect on R I . In both cases, the defects demonstrated higher R I values than the surface defects compared to the flange-or diaphragm-side defects.
For the defects, the crack tip line was parallel to the axis of symmetry of the specimen and was under a crack opening mode similar to that of the notched tensile test at both ends, resulting in different R I values. Considering the experimental results of η as a measure of plastic deformation, the plastic deformation capacity tended to increase as R I decreased, i.e., as the percentage of Mode II crack increased. However, these variations also resulted in brittle fractures, as shown for DTC-8W and DTC-12W, where the trends of the defect depth and plastic deformation capacity were reversed, whereas that of other conditions remained the same. Thus, further research is required to determine the factors that affect fracture toughness. Figure 14 shows the relationship between R I and fracture prediction accuracy, as shown in Table 5. High prediction accuracy was obtained for the three specimens with R I above 0.90 and with Mode I as the dominant mode, whereas the prediction accuracy for the three specimens with lower values of R I varied. No clear correlation was obtained. This suggests that the increased influence of Mode II resulted in an increased difficulty of failure prediction.
Appl. Mech. 2021, 2, FOR PEER REVIEW 13 conditions remained the same. Thus, further research is required to determine the factors that affect fracture toughness. Figure 14 shows the relationship between RI and fracture prediction accuracy, as shown in Table 5. High prediction accuracy was obtained for the three specimens with RI above 0.90 and with Mode I as the dominant mode, whereas the prediction accuracy for the three specimens with lower values of RI varied. No clear correlation was obtained. This suggests that the increased influence of Mode II resulted in an increased difficulty of failure prediction.

Conclusions
In this study, the TSM and Weibull stress approaches were used to predict the occurrence of brittle fractures in building structures using beam-to-diaphragm joint specimens. In addition, we attempted to quantify the mixed-mode ratio of the crack opening modes as an initial approach to investigate the prediction of brittle fractures under different crack opening modes. The findings of this study are as follows: 1. The TSM and Weibull stress approaches generated similar predictions. Further, both methods can be used for the prediction of brittle fractures. The prediction results obtained through both methods were similar. 2. We calculated the mixed-mode ratio RI using the stress intensity factor at the initial stage of deformation. Up to θ = 0.01 rad, the results obtained through the elastic and the elastic-plastic analyses were similar, and RI followed a linear trend. 3. RI demonstrated differences in the crack opening modes owing to deformations in crack tips, regardless of the crack tip radii. Based on RI, crack opening mode differences resulting from crack location and shape can be appropriately classified. The relationship between fracture prediction and RI suggests that different crack opening modes affect the brittle fractures of building structures.
Ductile crack growth was not considered in the analysis of this study. Therefore, if a substantial ductile crack propagates, the mixed-mode ratio may not capture the difference in modes for the actual specimen. Further study is needed to clarify the mixed-mode ratio's range of application and evaluate ductile crack growth.

Conclusions
In this study, the TSM and Weibull stress approaches were used to predict the occurrence of brittle fractures in building structures using beam-to-diaphragm joint specimens. In addition, we attempted to quantify the mixed-mode ratio of the crack opening modes as an initial approach to investigate the prediction of brittle fractures under different crack opening modes. The findings of this study are as follows: 1.
The TSM and Weibull stress approaches generated similar predictions. Further, both methods can be used for the prediction of brittle fractures. The prediction results obtained through both methods were similar.

2.
We calculated the mixed-mode ratio R I using the stress intensity factor at the initial stage of deformation. Up to θ = 0.01 rad, the results obtained through the elastic and the elastic-plastic analyses were similar, and R I followed a linear trend.

3.
R I demonstrated differences in the crack opening modes owing to deformations in crack tips, regardless of the crack tip radii. Based on R I , crack opening mode differences resulting from crack location and shape can be appropriately classified. The relationship between fracture prediction and R I suggests that different crack opening modes affect the brittle fractures of building structures.
Ductile crack growth was not considered in the analysis of this study. Therefore, if a substantial ductile crack propagates, the mixed-mode ratio may not capture the difference in modes for the actual specimen. Further study is needed to clarify the mixed-mode ratio's range of application and evaluate ductile crack growth.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study or in the collection, analyses, or interpretation of data as well as in the writing of the manuscript or in the decision to publish the results.