# Stress Intensity Factors for Corner Cracks under Arbitrary Stress Fields in the Finite Plates Based on the Point Weight Function Method

Authors: Guo Li; Huimin Zhou; Junbo Liu; Tongge Xu; Shuiting Ding
Journal: International Journal of Aerospace Engineering
 (2022)

---

## Abstract

Probabilistic damage tolerance assessment is an essential method to evaluate the safety of aeroengine rotors. The stress intensity factors (SIFs) are the core parameters. The weight function method can calculate SIF efficiently. The available weight functions for corner cracks are suitable for cracks with universal stresses. However, when cracks are under bivariant stress distributions, the lack of the weight function database makes the damage tolerance assessment impractical. Therefore, this paper derives the point weight functions for corner cracks, which are suited for cracks in real-world components with two-directional stress distributions. The response surface method is used to build the surrogate model of the point weight functions. By integrating the weight functions with the stress distributions, SIFs can be calculated with high accuracy. In sum, 81% of the differences between the point weight function method and finite element results are less than 10%.


---

## Body

## 1. Introduction

Aeroengine rotors are typical life-limited parts [1] that operate in complex environments with remarkably high temperatures, pressures, and rotational speeds [2, 3]. The probabilistic damage tolerance assessment methodology is applied to address the threat of defects that may occur anywhere in the components [4–7]. The probability of failure of these rare defects is predicted based on a generalized stress-strength interference model, which requires the simulation of the crack propagation process under complex loads [8, 9]. The stress intensity factor (SIF) is the core parameter assessed during a crack propagation simulation, as it reflects the severity and stress concentration around a crack [10].SIF solution research has been a central topic in fracture mechanics in recent decades. Various analytical and numerical methods for determining SIFs have been proposed. Representative SIF solutions of through cracks in two-dimensional geometries have been generated and collected in several well-known handbooks [11]. As to cracks in three-dimensional geometries, Newman and Raju [12, 13] have proposed SIF solutions for the corner, surface, and embedded cracks in three-dimensional geometries. Unfortunately, these solutions are expressed in simple load cases, including uniform remote tension and bending stress. Instead, finite element analysis or boundary element method is developed to obtain the SIFs of the cracks with arbitrary stresses. Nevertheless, the detailed three-dimensional numerical approach only calculates the singular SIF of the crack with specific loading and geometry with a complex preprocess to model and mesh the geometry. Thus, these numerical methods hardly meet the requirements of the highly efficient probabilistic damage tolerance assessment because millions of crack propagation processes need to be calculated [14–16].To efficiently obtain the SIFs for any specific cracked body geometry that is applied with an arbitrary stress distribution, the weight function method was envisioned by Bueckner [17] and studied further by Rice [18]. The weight function, also called Green’s function of a crack surface, is determined by the geometry of the cracks in the plate. The SIFs are calculated efficiently by integrating Green’s function and the stress distribution on the crack surface. The standardized analytical weight function method was developed mainly concentrating on the through-thickness cracks in two-dimensional geometries [19, 20]. For elliptical embedded cracks and semielliptical surface cracks in three-dimensional geometries, the universal weight function was proposed [21–23]. These two weight functions are widely used to calculate SIFs with univariate stress distributions [24, 25].Nevertheless, these weight functions are limited when the cracks are under complex bivariant stress distributions [26, 27]. The corner cracks are more likely to bear complex bivariant loadings because of the stress concentration, thermal loading, and machining or strengthening residual stress. Oore and Burns [28] proposed the point weight function method to calculate SIFs under bivariant stress distributions. Orynyak et al. then developed point weight functions for elliptical cracks [29, 30]. Finally, the point weight function method was applied on the corner cracks at holes for calculating SIFs under bivariant stress gradients and implemented in production software for design and life management [31–33]. However, few efforts have been devoted to explaining the details for establishing point weight functions for corner cracks in three-dimensional finite rectangular plates. The lack of the weight function database makes the probabilistic damage tolerance assessment or reliability analysis of corner cracks with arbitrary complex loadings impractical. Moreover, previous research ignored analyzing the accuracy of the point weight function method. To ensure the reliability of the derived point weight function database, the quantified error of the point weight function is supposed to be considered when calculating crack propagation process.This study derives the point weight functions for corner cracks to calculate SIFs of corner defects under arbitrary stress distributions and help conduct the probabilistic damage tolerance assessment [25]. The challenge of the point weight function method is to derive the coefficients of the weight functions determined by the geometric parameters of a structure with cracks. The response surface method [34, 35] is an effective way to derive coefficients in the point weight functions. The errors induced by the fitting method are less than 1%. Then, the SIFs for corner cracks calculated by the point weight function method are validated. When bivariant stress distributions are applied onto the crack surface, the differences between the point weight function method and finite element results are less than 5MPam when 0<a/L≤0.91 and 0<a/T≤0.91. Meanwhile, 81% of the relative errors are less than 10%. The errors of the SIFs calculated by the point weight function method may be induced by the response surface method and the accuracy of the reference SIFs. The errors induced by the response surface method are less than 1%. Moreover, the 1% error of the reference SIFs can cause a 1.25% deviation of the calculated SIFs.This paper is divided into five sections. The theoretical formula of the point weight function is analyzed in Section2. Section 3 introduces the derivation of the coefficients in point weight functions determined by multiple geometric parameters, including the length and thickness of a flat plate. In Section 4, the accuracy of the point weight function is quantified. Finally, the conclusions are summarized in Section 5.


## 2. Point Weight Function Theory for Bivariant Stress Field

Stress distributions of structural components are complicated such as aeroengines. These complex stress gradients are mainly induced by local stress concentrations or multiple loads, including centrifugal, thermal, and residual stresses. When the stresses vary in bivariant directions (vertical and along with the crack depth), the point weight function can help obtain accurate SIFs of cracks. The SIF considered in this paper is a mode I corner crack in rectangular geometry with finite length and width, as shown in Figure1. The SIFs for two points along the corner cracks, A and B, are critical for the crack propagation calculation [31].Figure 1
Geometry and coordinate system.When the point weight function is known, the SIF, the core parameter used in probabilistic damage tolerance analysis, can be obtained by integrating the product of the load with the weight function for the crack surface. This method considers two variations compared with the universal weight function method, which only considers stress distribution changes in thex direction.The weight function is called Green’s function of a crack surface. The SIF is given by [27]
(1)K=∫0c∫0a1−x2/c2σx,yWQQ′dydx,where σx,y is the stress perpendicular to the crack plane in the corresponding uncracked body and WQQ′ is the point weight function. The integration is carried out for the crack dimensions, characterized by a in the x direction and c in the y direction. Both a and c are the radii of the corner crack. Only circular cracks are considered in this paper; thus, a=c.The point weight functionsWQQ′ for quarter-circular corner cracks are, respectively, given by [26]
(2)WQQ′=R2−r2ππR1lQQ′21+lQQ′2lQxQ′2+lQQ′2lQyQ′21+π11−rR+π21−yy′+π31−xx′,where Qx,y denotes the location where the point load is applied. The two points Qx and Qy do not exist on the crack surface. The locations of Qx and Qy are those symmetrical to the location of point load Q with respect to the x- and y-axes. Q′ is the point along the crack perimeter. x and y are the values of the x-axes and y-axes. R represents the distance between the point along with the crack front Q′ to the ellipse center O, r represents the distance between the point in the crack surface Q to the ellipse center O,lQQ′ is the distance between Q and Q′,lQxQ′ is the distance between Qx and Q′, and lQyQ′ is the distance between Qy and Q′. The length parameters, lQQ′, lQxQ′, and lQyQ′, correct the free surface effects for a surface crack by assuming a symmetrical stress distribution for an imaginary prolonged crack extending into the other semicircles. Π1, Π2, and Π3 are the critical coefficients of the point weight function, as determined by the crack size and geometry parameters.ΠEquation (2) can be expressed in cylindrical coordinates
(3)WQQ′=R2−r2π1R2+r2−2Rrcosφ+θ1+R2+r2−2Rrcosφ−θR2+r2−2Rrcosφ+θ+R2+r2−2Rrcosφ−θR2+r2+2Rrcosφ+θ1+∏11−rR+∏21−rsinφR2−r2cosφ2+∏31−rcosφR2−r2sinφ2,where φ and θ are the angles of Q and Q′.With the integration of stress distributionσx,y and point weight function, the SIF of Q′ along the surface crack can be obtained by
(4)K=∫0π/2∫0Rσr,θR2−r2ππR1R2+r2−2Rrcosφ−θ1+R2+r2−2Rrcosφ−θR2+r2−2Rrcosφ+θ+R2+r2−2Rrcosφ−θR2+r2+2Rrcosφ+θ1+∏11−rR+∏21−rsinφR2−r2cosφ2+∏31−rcosφR2−r2sinφ2rdrdθ,The SIFs for pointsA and B are concerned Q′=A,B. For each point, by using three known SIFs Krlk,l=1,2,3,k=A,B and the associated three reference stress distributions (σrlr,θ,i=1,2,3), six simultaneous equations are given as
(5)Krlk=∫0c∫0a1−x2/c2σrlr,θR2/r2ππR1R2+r2−2Rrcosφ−θ1+R2+r2−2Rrcosφ−θR2+r2−2Rrcosφ+θ+R2+r2−2Rrcosφ−θR2+r2+2Rrcosφ+θ1+Π1k1−rR+Π2k1−rsinφR2−r2cosφ2+Π3k1−rcosφR2−r2sinφ2rdrdθl=1,2,3,k=A,B.The coefficientsΠ1k, Π2k, and Π3k can be determined by solving Equation (5). Finally, the point weight functions can be derived, and the SIFs of corner cracks under unknown stress distributions can be calculated by integrating the weight functions and stress distributions as Equation (4).


## 3. Derivation of the Point Weight Functions for Corner Cracks

### 3.1. Reference Stress Intensity Factors

Π1A, Π2A, Π3A, Π1B, Π2B, and Π3B are the critical coefficients of the point weight functions obtained by solving Equation (5). The construction of the equations leads to the acquisition of the reference SIFs under reference stress distributions. Three-dimensional finite elements are used to model the symmetric structure containing a quarter-circular corner crack to obtain the reference SIFs. Figure 1 shows the geometry and the coordinate system of the corner cracks. T and L are the thickness and length of the finite plate, respectively.Finite element analysis is performed using ABAQUS version 6.14 with 20-node isoperimetric three-dimensional solid elements and reduced integration. As shown in Figure2, three-dimensional prism elements with four midside nodes at the quarter points are used to model the square root singularity at the crack tip, and the separate crack tip nodal points are constrained and assumed to have the same displacement.Figure 2
Typical finite element mesh.The SIFs for two points of corner cracks,A and B, shown in Figure 1, are critical for the crack propagation calculation [31]. It is worth mentioning that the reference solutions are determined as the SIFs taken at 2.5° inside the two points A and B. This approach aims to avoid the numerical ambiguities resulting from the corner point crack singularity associated with the calculation of the SIF directly at the surface.The SIF is calculated from theJ-integral by using the virtual crack extension method. The analysis is conducted with a linear classic elastic material model with Young’s modulus E (210 GPa) and Poisson’s ratio ν (0.361). The relationship between J and SIF can be used to calculate the reference SIF
(6)Kr=JE′where E′=E is the equation for the plane stress and E′=E/1−ν2 is the formula for the plane strain.The reference loads are applied directly onto the crack surface. Three types of loading are applied to each crack geometry, with the stress distribution given by(7)σr1x,y=σ0,σr2x,y=σ01−,σr3x,y=σ01−,where σ0 is a constant, in which the value of 400 MPa is assigned and a is the crack radius whose value is 0.003 m in this paper.


### 3.2. Validation of the Finite Element Method

The accuracy of the finite element results needs to be validated because the reference SIFs are critical inputs for the point weight function method. Therefore, the SIFs calculated by ABAQUS are compared with FADD3D results [32] to validate the accuracy of the reference SIFs. The univariant stress distribution is given by
(8)σx,y=σ01−ya1−ya1−ya.The bivariant stress distribution is given by(9)σx,y=σ01−ya1−ya1−ya1−xa1−xa1−xa.The bivariant stress distribution decreases to zero at both surface tips, while the univariant stress distribution is uniform in they-direction and applies a much larger resultant force on the crack plane than the bivariant stress field.Figure3 shows the SIF results along the crack tip calculated by FADD3D and ABAQUS. The circular crack is analyzed, and a/c=1. The geometry of the finite plate is set as a/T=0.5 and c/L=0.5. Q is the shape factor associated with the a/c ratio, defined by
(10)Q=1+1.464ac1.65ac≤1,1+1.464ca1.65ac>1.Figure 3
SIF results when cracks under univariant and bivariant stress distributions calculated by FADD3D and ABAQUS (c/a=1,a/T=0.5,c/L=0.5).Notably, the predicted results exist sharp changes at the intersection of the crack perimeter with a free surface (A and B tips) because of the complex nature of the singularity at that point. Therefore, the solution values for A and B tips are generally taken at the first node inside the free surface, 2.5° from the free surface.MAPE is the average of the absolute relative errors of the actual and predicted values. The smaller value of MAPE indicates effective model prediction.(11)MAPE=100%Nall∑i=1NallKWFM−KFEMKFEM,where KWFM is the SIF results calculated by the weight function method, KFEM is the SIF results calculated by finite element method, and Nall represents the number of the SIF values.The maximum error between these two numerical methods is less than 3.12%. The average absolute relative error MAPE of universal stress distribution between ABAQUS and FADD3D is 1.37%, while the average absolute relative error MAPE of bivariant stress distribution is 1.80%. The validation ensures the accuracy of the reference SIFs in this paper.


### 3.3. Weight Function Derivation

The SIFs are determined by the geometry size parametersa,T, and L. Supposing H is infinite, and the other three variables are then regarded as critical parameters and normalized as a/Landa/T. The coefficients of weight functions Πlkl=1,2,3,k=A,B are also determined by the geometry size parameters a,T,andL, given by
(12)Πlk=FaT,aL,where a/L=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06. a/T is the same as a/L.The discrete data are calculated by finite element analysis. The response surface method [34, 35] is applied to obtain the relationship between the coefficients and the normalized size parameters. The theoretical formula is given by
(13)Πlk=p0+∑i1=12pi1ei1+∑i1=12∑i2=12pi1i2ei1ei2+∑i3=12∑i1=12∑i2=12pi1i2i3ei1ei2ei3+∑i3=12⋯∑ij=12pi1⋯ijei1⋯eij+ε,where j is the highest order of polynomial; p0,pi1⋯ij represent the undetermined coefficients; ei1,⋯,eij represent the normalized size parameters; e1 and e2 represent a/T and a/L, respectively; and ε is the error term. The undetermined coefficient is supposed to be obtained by the regression algorithm, which is denoted as p^0,p^i1⋯ij. Thus, Equation (14) can be expressed by
(14)Πlk=p^0+∑i1=12p^i1ei1+∑i1=12∑i2=12p^i1i2ei1ei2+∑i3=12∑i1=12∑i2=12p^i1i2i3ei1ei2ei3+∑i3=12⋯∑ij=12p^i1⋯ijei1⋯eij+ε.


### 3.4. Stress Fit Formula

The SIF denoted byK can be expressed by integrating the stress and the weight function as Equation (4). In this case, the bivariate stress function is another critical input of the point weight function method. After a finite element analysis of rotor disks, the stress distribution is extracted. Then, the stress fit process can help obtain the stress distribution function, given by
(15)σx,y=∑κ=1m∑η=1nβκηxκ−1yη−1,where x and y represent the coordinate value on the crack surface, m and n are the degrees of a polynomial function, and β is the coefficient. The Cartesian coordinates are converted to cylindrical coordinates as
(16)σφ,r=∑κ=1m∑η=1nβκηrcosφκ−1rsinφη−1.


## 3.1. Reference Stress Intensity Factors

Π1A, Π2A, Π3A, Π1B, Π2B, and Π3B are the critical coefficients of the point weight functions obtained by solving Equation (5). The construction of the equations leads to the acquisition of the reference SIFs under reference stress distributions. Three-dimensional finite elements are used to model the symmetric structure containing a quarter-circular corner crack to obtain the reference SIFs. Figure 1 shows the geometry and the coordinate system of the corner cracks. T and L are the thickness and length of the finite plate, respectively.Finite element analysis is performed using ABAQUS version 6.14 with 20-node isoperimetric three-dimensional solid elements and reduced integration. As shown in Figure2, three-dimensional prism elements with four midside nodes at the quarter points are used to model the square root singularity at the crack tip, and the separate crack tip nodal points are constrained and assumed to have the same displacement.Figure 2
Typical finite element mesh.The SIFs for two points of corner cracks,A and B, shown in Figure 1, are critical for the crack propagation calculation [31]. It is worth mentioning that the reference solutions are determined as the SIFs taken at 2.5° inside the two points A and B. This approach aims to avoid the numerical ambiguities resulting from the corner point crack singularity associated with the calculation of the SIF directly at the surface.The SIF is calculated from theJ-integral by using the virtual crack extension method. The analysis is conducted with a linear classic elastic material model with Young’s modulus E (210 GPa) and Poisson’s ratio ν (0.361). The relationship between J and SIF can be used to calculate the reference SIF
(6)Kr=JE′where E′=E is the equation for the plane stress and E′=E/1−ν2 is the formula for the plane strain.The reference loads are applied directly onto the crack surface. Three types of loading are applied to each crack geometry, with the stress distribution given by(7)σr1x,y=σ0,σr2x,y=σ01−,σr3x,y=σ01−,where σ0 is a constant, in which the value of 400 MPa is assigned and a is the crack radius whose value is 0.003 m in this paper.


## 3.2. Validation of the Finite Element Method

The accuracy of the finite element results needs to be validated because the reference SIFs are critical inputs for the point weight function method. Therefore, the SIFs calculated by ABAQUS are compared with FADD3D results [32] to validate the accuracy of the reference SIFs. The univariant stress distribution is given by
(8)σx,y=σ01−ya1−ya1−ya.The bivariant stress distribution is given by(9)σx,y=σ01−ya1−ya1−ya1−xa1−xa1−xa.The bivariant stress distribution decreases to zero at both surface tips, while the univariant stress distribution is uniform in they-direction and applies a much larger resultant force on the crack plane than the bivariant stress field.Figure3 shows the SIF results along the crack tip calculated by FADD3D and ABAQUS. The circular crack is analyzed, and a/c=1. The geometry of the finite plate is set as a/T=0.5 and c/L=0.5. Q is the shape factor associated with the a/c ratio, defined by
(10)Q=1+1.464ac1.65ac≤1,1+1.464ca1.65ac>1.Figure 3
SIF results when cracks under univariant and bivariant stress distributions calculated by FADD3D and ABAQUS (c/a=1,a/T=0.5,c/L=0.5).Notably, the predicted results exist sharp changes at the intersection of the crack perimeter with a free surface (A and B tips) because of the complex nature of the singularity at that point. Therefore, the solution values for A and B tips are generally taken at the first node inside the free surface, 2.5° from the free surface.MAPE is the average of the absolute relative errors of the actual and predicted values. The smaller value of MAPE indicates effective model prediction.(11)MAPE=100%Nall∑i=1NallKWFM−KFEMKFEM,where KWFM is the SIF results calculated by the weight function method, KFEM is the SIF results calculated by finite element method, and Nall represents the number of the SIF values.The maximum error between these two numerical methods is less than 3.12%. The average absolute relative error MAPE of universal stress distribution between ABAQUS and FADD3D is 1.37%, while the average absolute relative error MAPE of bivariant stress distribution is 1.80%. The validation ensures the accuracy of the reference SIFs in this paper.


## 3.3. Weight Function Derivation

The SIFs are determined by the geometry size parametersa,T, and L. Supposing H is infinite, and the other three variables are then regarded as critical parameters and normalized as a/Landa/T. The coefficients of weight functions Πlkl=1,2,3,k=A,B are also determined by the geometry size parameters a,T,andL, given by
(12)Πlk=FaT,aL,where a/L=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06. a/T is the same as a/L.The discrete data are calculated by finite element analysis. The response surface method [34, 35] is applied to obtain the relationship between the coefficients and the normalized size parameters. The theoretical formula is given by
(13)Πlk=p0+∑i1=12pi1ei1+∑i1=12∑i2=12pi1i2ei1ei2+∑i3=12∑i1=12∑i2=12pi1i2i3ei1ei2ei3+∑i3=12⋯∑ij=12pi1⋯ijei1⋯eij+ε,where j is the highest order of polynomial; p0,pi1⋯ij represent the undetermined coefficients; ei1,⋯,eij represent the normalized size parameters; e1 and e2 represent a/T and a/L, respectively; and ε is the error term. The undetermined coefficient is supposed to be obtained by the regression algorithm, which is denoted as p^0,p^i1⋯ij. Thus, Equation (14) can be expressed by
(14)Πlk=p^0+∑i1=12p^i1ei1+∑i1=12∑i2=12p^i1i2ei1ei2+∑i3=12∑i1=12∑i2=12p^i1i2i3ei1ei2ei3+∑i3=12⋯∑ij=12p^i1⋯ijei1⋯eij+ε.


## 3.4. Stress Fit Formula

The SIF denoted byK can be expressed by integrating the stress and the weight function as Equation (4). In this case, the bivariate stress function is another critical input of the point weight function method. After a finite element analysis of rotor disks, the stress distribution is extracted. Then, the stress fit process can help obtain the stress distribution function, given by
(15)σx,y=∑κ=1m∑η=1nβκηxκ−1yη−1,where x and y represent the coordinate value on the crack surface, m and n are the degrees of a polynomial function, and β is the coefficient. The Cartesian coordinates are converted to cylindrical coordinates as
(16)σφ,r=∑κ=1m∑η=1nβκηrcosφκ−1rsinφη−1.


## 4. Validation of the Point Weight Functions

### 4.1. Comparison with the Universal Weight Function Method

The universal weight function derived by Glinka is suitable to calculate the stress intensity factor (SIF) for quarter-elliptical corner crack [23] in a finite thickness plate. As shown in Figure 4, T1,L1,andH1 are the plate’s thickness, width, and height. The universal weight function derived by Glinka is suitable when L1andH1 are infinite.Figure 4
Quarter-elliptical corner crack in finite thickness plate [23].The universal weight function method takes univariant stress distribution into consideration. Thus, two types of stresses are applied to the crack surface, given by(17)σx,y=σ0=1−xa1−xa1−xa,(18)σx,y=σ01−xa1+xa.Then, the SIFs are calculated using the point weight function method and universal weight function method [22, 23]. Meanwhile, the finite element method results calculated by ABAQUS are compared with these two weight function methods, as shown in Figure 5.Figure 5
Comparison of SIF results based on universal weight function method, point weight function method, and finite element method when cracks under univariant stress distributions (L1⟶∞, H1⟶∞).
(a)
Comparison of SIF results when cracks under stress distribution Equation (17)(b)
Comparison of SIF results when cracks under stress distribution Equation (18)Figure5 illustrates that the point weight function method and the universal weight function method fit the finite element method results well. The MAPE between the point weight function and the universal weight function method is less than 0.64% for point A and 1.46% for point B. Therefore, the accuracy of the point weight function method is verified under univariant stress distributions based on the comparison of the point weight function method results and the universal weight function method results.


### 4.2. Comparison with the Finite Element Results

#### 4.2.1. Stress Distributions

To validate the SIF results calculated by the point weight function derived in Section3, arbitrary stress distributions are applied onto the crack surface to replace the reference stresses (σrlx,y,l=1,2,3). Fourteen stress distributions along the x and y directions are applied onto the crack surface, as shown in Figure 6. The stress gradients are listed in Table 1. σ0 is assigned 400 MPa, and a is assigned 0.003 m in this paper. Stresses 1-5 change in the x direction as shown in Figure 6(a). Only these stress distributions are concerned by the universal weight function method. Stresses 6-9 change in the y direction as shown in Figure 6(b), while stresses 10-14 are the bivariate gradients shown in Figures 6(c)–6(g). The point weight function method can calculate the SIFs when there are cracks under these three types of stresses. The Πlk varies associated with two dimensionless geometric ratios a/L and a/T. To obtain the weight function coefficients when 0<a/T≤0.91 and 0<a/L≤0.91, 49 geometries are considered when a/T=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06 and a/L=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06.Figure 6
Stress distributions on the corner crack.
(a)
Stresses 1-5(b)
Stresses 6-9(c)
Stress 10(d)
Stress 11(e)
Stress 12(f)
Stress13(g)
Stress 14Table 1
Stress distributions applied onto the crack surface.
Stress numberxyStress distributions1x2σx,y=σ01−2x/a1−x/a2x3σx,y=σ01−x/a1−x/a1−x/a3x4σx,y=−σ01−4x/a1−6x/a1−x/a1−2x/a4x5σx,y=−σ01−4x/a1−2x/a1−1.05x/a1−0.9x/a1−0.8x/a5x6σx,y=−σ01−8x/a1−4x/a1−2x/a1−x/a1−0.9x/a1−0.8x/a1−0.6x/a+5/4σ06y1σx,y=σ0/41−6y/a7y2σx,y=σ01−y/a1+y/a8y3σx,y=σ0/801−6y/a1−6y/a1−6y/a9y4σx,y=σ01−y/a1−y/a1−2y/a1−3y/a10xyσx,y=σ01−2x/a1−6y/a11xy3σx,y=σ01−6y/a1−3y/a1−2x/a1−y/a12xy7σx,y=σ0x/ay/a713x2yσx,y=−σ01−8x/a1−4x/a1−2y/a+5/4σ014x2y3σx,y=−σ01−x/a1−0.5x/a1−4y/a1−2y/a1−y/a


#### 4.2.2. Validation Results

TheΠlk varies associated with two dimensionless geometric ratios a/L and a/T. To obtain the weight function coefficients when 0<a/T≤0.91 and 0<a/L≤0.91, the surrogate model is established as Equation (10), and the surrogate model of the point weight function coefficients is shown in the appendix. The SIFs calculated by the obtained point weight functions with Πlk (surrogate model) are compared with finite element results to validate the accuracy of the weight function method. The SIF values are normalized by
(19)NormalizedSIF=Kσ0πa.These 49 geometries with different sizes are numbered as shown in Figure7(a). The SIF results under stresses 1-14 are calculated, and Figures 7(a)–7(l) show the results of stresses 1, 5, 6, 7, 10, and 14. The normalized SIFs calculated by the finite element method for point B and point A of the 49 geometries are marked with hollow and solid points, respectively. With the surrogate model, the SIFs are calculated and shown as three-dimensional surfaces (0<a/T≤0.91 and 0<a/L≤0.91) in Figures 7(a), 7(c), 7(e), 7(g), 7(i), and 7(k). Figures 7(b), 7(d), 7(f), 7(h), 7(j), and 7(l) show the relationship between the finite element results and the weight function method results of the 49 geometries.Figure 7
Comparison of SIFs on corner cracks for pointA and point B with different stress distributions.
(a)
3D SIF results under stress 1(b)
2D SIF results under stress 1(c)
3D SIF results under stress 5(d)
2D SIF results under stress 5(e)
3D SIF results under stress 6(f)
2D SIF results under stress 6(g)
3D SIF results under stress 7(h)
2D SIF results under stress 7(i)
3D SIF results under stress 10(j)
2D SIF results under stress 10(k)
3D SIF results under stress 14(l)
2D SIF results under stress 14Results show that the SIFs are determined by the geometries with the same loading. Whena/T and a/L are close to 1, the SIFs are 3-6 times larger than SIFs when a/T and a/L are close to 0. The validation result also shows that the SIFs calculated by the two methods are consistent when a/L and a/T are small.The relative errorCWFM between the weight function method and finite element method is defined as
(20)CWFM=KWFM−KFEMKFEM,where KWFM is the SIF results calculated by the weight function method and KFEM is the SIF results calculated by a finite element method. Similarly, the absolute error ΔKWFM between the weight function method and finite element method is defined as
(21)ΔKWFM=KWFM−KFEM.The specific errors between the weight function method and finite element method results of these 14 stresses are shown in Table2. The relative errors are larger than 100% when the absolute SIF value is small, less than 1 MPam. Most of the differences between the weight function method and the finite element results are less than 5 MPam when 0<a/L≤0.91 and 0<a/T≤0.91 except for stress 3, 10, 11, and 13 cases. The errors are related to the monotonicity of the stress distribution. Stress 1 can represent most stress concentration conditions that occur in rotor disks.Table 2
Errors between point weight function method and finite element results.
Stress typeMax errorAverage error (MAPE)ABABABABKUWF/KFEM−1KUWF−KFEMKUWF/KFEM−1KUWF−KFEM1127.27%6.92%6.2972.70137.92%1.84%1.2680.507226.33%3.44%3.9661.35710.59%1.05%0.8050.273357.53%30.86%10.9043.34915.87%6.44%2.2330.584456.23%27.46%2.3016.43721.00%7.55%0.4991.19358.06%10.17%8.5029.5902.64%2.95%1.6561.708625.26%1.44%0.0470.4001.90%0.91%0.0080.24571.74%5.59%1.2943.2950.62%2.54%0.2660.623858.93%31.64%2.00910.82022.86%7.60%0.4152.001920.56%45.80%5.7582.7466.26%17.74%1.1400.58710105.03%16.86%19.73620.49133.50%4.79%4.8034.2841187.32%59.31%15.2342.92719.47%13.92%3.1120.5321235.36%9.98%0.1250.16514.05%3.58%0.0220.0331352.64%33.00%14.6122.39612.20%11.98%2.9360.541145.32%6.75%0.6320.5424.23%2.15%0.1330.103The cumulative distributions ofCWFM and ΔKWFM are shown in Figures 8(a) and 8(b). Results show that 81% CWFM of point B and 54% of point A are less than 10%. As shown in Figure 8(b), for point B, 97% ΔKWFM are less than 5 MPam while 93% for point A. The absolute values of point B are usually larger than point A; thus, the relative error of point B is more minor.Figure 8
The errors of SIFs calculated by weight function method and finite element method.
(a)
Relative errors(b)
Absolute errorIn sum, the derived point weight functions can help obtain SIFs of corner cracks in disks under bivariant stress distributions. The relative error may increase when the absolute SIFs are less than 1MPam, and 81% of the relative errors for point B and 54% for point A are less than 10%. Meanwhile, 90% of the absolute errors between the point weight function and finite element method results are less than 5 MPam.


### 4.3. Error Analysis

#### 4.3.1. Errors Induced by the Surrogate Model

Based on the finite element analysis mentioned in Section3.1, Forty-nine geometries with different sizes (a/T=a/L=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06) are established, and the associated discrete weight function coefficients Πlk are obtained by solving Equations (5). But to obtain the point weight function coefficients when 0<a/T≤0.91 and 0<a/L≤0.91, the surrogate model is established by the response surface method. The surrogate model will induce errors in the SIFs.The hexagonal polynomials (j=6) fit well, and the coefficient of determination (R-square) value reaches 0.9995. The stress distributions 1-14 mentioned in Table 1 are applied on the crack surface, and the SIFs of the 49 geometries are calculated. The relative error of the SIFs calculated by discrete Πlk and surrogate model is given as
(22)Cfit=Krps−KdisKdis,where Krps is the SIF results calculated by the surrogate model and Kdis is the SIF results calculated by discrete Πlk. The cumulative distribution of Cfit is shown in Figure 9. Results show that 93% Cfit is less than 1%. The response surface method helps obtain the surrogate model of weight function coefficients Πlk with relatively small errors.Figure 9
The relative errors of SIF calculated by discrete coefficients and surrogate model.


#### 4.3.2. Errors Induced by the Reference Stress Intensity Factors

The accuracy of the reference SIFs directly ensures the effectiveness of theΠlk as shown in Equations (5). It is assumed that the three reference solutions for point AKr1A, Kr2A, and Kr3A deviate to the actual finite element results ranging from -20% to 20%. Thus, the three reference solutions are given as
(23)KrlAd=KrlA+αKrlA,where α is the offset scale and ranges from -20% to 20%.The proportion of the new coefficients derived by the offset reference solutions are defined as(24)γlAd=ΠlAdΠlA,where ΠlAd are the coefficients derived by the deviated reference solutions.Similarly, the proportion of the SIFs calculated by the new coefficients are given as(25)δlAd=,where KlAd are the SIFs calculated by the weight function method. Equations (23)–(25) also apply to point B.TheΠlk show a multivariate linear correlation of the offset scale in Figure 10 because the Πlk are solved by Equation (5). When the reference SIFs have a 1% error, the Πlk will change to 5 times of the original value. The SIFs calculated by the point weight function also change, as shown in Figure 11. When the reference SIFs have a 1% error, the SIFs for points A and B have an error up to 1.25%. The coefficients of the point weight function are highly dependent on the reference SIFs. Thus, the precise reference SIFs are significant to the point weight function method.Figure 10
The proportion of the coefficientsγlAd and γlBd with different offset scale.
(a)
PointA(b)
PointBFigure 11
The proportion of the coefficientsδlAd with different offset scale.The errors of the SIFs calculated by the point weight function method may be induced by the response surface method and the accuracy of the reference SIFs. First, the response surface fitting process of the discrete weight function coefficients causes the deviation of the SIFs. By fitting the coefficients with high order polynomial, the errors are less than 1%. Second, the accuracy of the reference SIFs calculated by the finite element method will influence the results of the weight function coefficients up to 5 times since these coefficients are calculated by solving simultaneous equations. A 1% error of the reference SIFs can cause a 1.25% deviation of the calculated SIFs.In addition to these two factors, the selection of the types of the reference stress distributions influences the accuracy of the point weight function. Since the uniform and linear reference stresses are used in this paper, the SIF results of linear distributions are accurate. Further research is necessary to study the influence of the selection of the reference loadings. Moreover, the point weight function equations can approximate Green’s function instead of equaling to Green’s function completely. The SIF errors are inevitably induced because the weight functions deviate from Green’s function, which is the limitation of the weight function method itself.In this case, during the crack propagation calculation and damage tolerance assessment, the SIF errors are supposed to be considered. Meanwhile, the development of the accurate finite element solutions is highly essential, which can provide the reference solutions to the point weight function method. Apart from the weight function method, a new surrogate model between the SIFs and the stress distributions of different geometries with cracks may be developed by other methods such as the machine learning method with enough SIF data.


## 4.1. Comparison with the Universal Weight Function Method

The universal weight function derived by Glinka is suitable to calculate the stress intensity factor (SIF) for quarter-elliptical corner crack [23] in a finite thickness plate. As shown in Figure 4, T1,L1,andH1 are the plate’s thickness, width, and height. The universal weight function derived by Glinka is suitable when L1andH1 are infinite.Figure 4
Quarter-elliptical corner crack in finite thickness plate [23].The universal weight function method takes univariant stress distribution into consideration. Thus, two types of stresses are applied to the crack surface, given by(17)σx,y=σ0=1−xa1−xa1−xa,(18)σx,y=σ01−xa1+xa.Then, the SIFs are calculated using the point weight function method and universal weight function method [22, 23]. Meanwhile, the finite element method results calculated by ABAQUS are compared with these two weight function methods, as shown in Figure 5.Figure 5
Comparison of SIF results based on universal weight function method, point weight function method, and finite element method when cracks under univariant stress distributions (L1⟶∞, H1⟶∞).
(a)
Comparison of SIF results when cracks under stress distribution Equation (17)(b)
Comparison of SIF results when cracks under stress distribution Equation (18)Figure5 illustrates that the point weight function method and the universal weight function method fit the finite element method results well. The MAPE between the point weight function and the universal weight function method is less than 0.64% for point A and 1.46% for point B. Therefore, the accuracy of the point weight function method is verified under univariant stress distributions based on the comparison of the point weight function method results and the universal weight function method results.


## 4.2. Comparison with the Finite Element Results

### 4.2.1. Stress Distributions

To validate the SIF results calculated by the point weight function derived in Section3, arbitrary stress distributions are applied onto the crack surface to replace the reference stresses (σrlx,y,l=1,2,3). Fourteen stress distributions along the x and y directions are applied onto the crack surface, as shown in Figure 6. The stress gradients are listed in Table 1. σ0 is assigned 400 MPa, and a is assigned 0.003 m in this paper. Stresses 1-5 change in the x direction as shown in Figure 6(a). Only these stress distributions are concerned by the universal weight function method. Stresses 6-9 change in the y direction as shown in Figure 6(b), while stresses 10-14 are the bivariate gradients shown in Figures 6(c)–6(g). The point weight function method can calculate the SIFs when there are cracks under these three types of stresses. The Πlk varies associated with two dimensionless geometric ratios a/L and a/T. To obtain the weight function coefficients when 0<a/T≤0.91 and 0<a/L≤0.91, 49 geometries are considered when a/T=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06 and a/L=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06.Figure 6
Stress distributions on the corner crack.
(a)
Stresses 1-5(b)
Stresses 6-9(c)
Stress 10(d)
Stress 11(e)
Stress 12(f)
Stress13(g)
Stress 14Table 1
Stress distributions applied onto the crack surface.
Stress numberxyStress distributions1x2σx,y=σ01−2x/a1−x/a2x3σx,y=σ01−x/a1−x/a1−x/a3x4σx,y=−σ01−4x/a1−6x/a1−x/a1−2x/a4x5σx,y=−σ01−4x/a1−2x/a1−1.05x/a1−0.9x/a1−0.8x/a5x6σx,y=−σ01−8x/a1−4x/a1−2x/a1−x/a1−0.9x/a1−0.8x/a1−0.6x/a+5/4σ06y1σx,y=σ0/41−6y/a7y2σx,y=σ01−y/a1+y/a8y3σx,y=σ0/801−6y/a1−6y/a1−6y/a9y4σx,y=σ01−y/a1−y/a1−2y/a1−3y/a10xyσx,y=σ01−2x/a1−6y/a11xy3σx,y=σ01−6y/a1−3y/a1−2x/a1−y/a12xy7σx,y=σ0x/ay/a713x2yσx,y=−σ01−8x/a1−4x/a1−2y/a+5/4σ014x2y3σx,y=−σ01−x/a1−0.5x/a1−4y/a1−2y/a1−y/a


### 4.2.2. Validation Results

TheΠlk varies associated with two dimensionless geometric ratios a/L and a/T. To obtain the weight function coefficients when 0<a/T≤0.91 and 0<a/L≤0.91, the surrogate model is established as Equation (10), and the surrogate model of the point weight function coefficients is shown in the appendix. The SIFs calculated by the obtained point weight functions with Πlk (surrogate model) are compared with finite element results to validate the accuracy of the weight function method. The SIF values are normalized by
(19)NormalizedSIF=Kσ0πa.These 49 geometries with different sizes are numbered as shown in Figure7(a). The SIF results under stresses 1-14 are calculated, and Figures 7(a)–7(l) show the results of stresses 1, 5, 6, 7, 10, and 14. The normalized SIFs calculated by the finite element method for point B and point A of the 49 geometries are marked with hollow and solid points, respectively. With the surrogate model, the SIFs are calculated and shown as three-dimensional surfaces (0<a/T≤0.91 and 0<a/L≤0.91) in Figures 7(a), 7(c), 7(e), 7(g), 7(i), and 7(k). Figures 7(b), 7(d), 7(f), 7(h), 7(j), and 7(l) show the relationship between the finite element results and the weight function method results of the 49 geometries.Figure 7
Comparison of SIFs on corner cracks for pointA and point B with different stress distributions.
(a)
3D SIF results under stress 1(b)
2D SIF results under stress 1(c)
3D SIF results under stress 5(d)
2D SIF results under stress 5(e)
3D SIF results under stress 6(f)
2D SIF results under stress 6(g)
3D SIF results under stress 7(h)
2D SIF results under stress 7(i)
3D SIF results under stress 10(j)
2D SIF results under stress 10(k)
3D SIF results under stress 14(l)
2D SIF results under stress 14Results show that the SIFs are determined by the geometries with the same loading. Whena/T and a/L are close to 1, the SIFs are 3-6 times larger than SIFs when a/T and a/L are close to 0. The validation result also shows that the SIFs calculated by the two methods are consistent when a/L and a/T are small.The relative errorCWFM between the weight function method and finite element method is defined as
(20)CWFM=KWFM−KFEMKFEM,where KWFM is the SIF results calculated by the weight function method and KFEM is the SIF results calculated by a finite element method. Similarly, the absolute error ΔKWFM between the weight function method and finite element method is defined as
(21)ΔKWFM=KWFM−KFEM.The specific errors between the weight function method and finite element method results of these 14 stresses are shown in Table2. The relative errors are larger than 100% when the absolute SIF value is small, less than 1 MPam. Most of the differences between the weight function method and the finite element results are less than 5 MPam when 0<a/L≤0.91 and 0<a/T≤0.91 except for stress 3, 10, 11, and 13 cases. The errors are related to the monotonicity of the stress distribution. Stress 1 can represent most stress concentration conditions that occur in rotor disks.Table 2
Errors between point weight function method and finite element results.
Stress typeMax errorAverage error (MAPE)ABABABABKUWF/KFEM−1KUWF−KFEMKUWF/KFEM−1KUWF−KFEM1127.27%6.92%6.2972.70137.92%1.84%1.2680.507226.33%3.44%3.9661.35710.59%1.05%0.8050.273357.53%30.86%10.9043.34915.87%6.44%2.2330.584456.23%27.46%2.3016.43721.00%7.55%0.4991.19358.06%10.17%8.5029.5902.64%2.95%1.6561.708625.26%1.44%0.0470.4001.90%0.91%0.0080.24571.74%5.59%1.2943.2950.62%2.54%0.2660.623858.93%31.64%2.00910.82022.86%7.60%0.4152.001920.56%45.80%5.7582.7466.26%17.74%1.1400.58710105.03%16.86%19.73620.49133.50%4.79%4.8034.2841187.32%59.31%15.2342.92719.47%13.92%3.1120.5321235.36%9.98%0.1250.16514.05%3.58%0.0220.0331352.64%33.00%14.6122.39612.20%11.98%2.9360.541145.32%6.75%0.6320.5424.23%2.15%0.1330.103The cumulative distributions ofCWFM and ΔKWFM are shown in Figures 8(a) and 8(b). Results show that 81% CWFM of point B and 54% of point A are less than 10%. As shown in Figure 8(b), for point B, 97% ΔKWFM are less than 5 MPam while 93% for point A. The absolute values of point B are usually larger than point A; thus, the relative error of point B is more minor.Figure 8
The errors of SIFs calculated by weight function method and finite element method.
(a)
Relative errors(b)
Absolute errorIn sum, the derived point weight functions can help obtain SIFs of corner cracks in disks under bivariant stress distributions. The relative error may increase when the absolute SIFs are less than 1MPam, and 81% of the relative errors for point B and 54% for point A are less than 10%. Meanwhile, 90% of the absolute errors between the point weight function and finite element method results are less than 5 MPam.


## 4.2.1. Stress Distributions

To validate the SIF results calculated by the point weight function derived in Section3, arbitrary stress distributions are applied onto the crack surface to replace the reference stresses (σrlx,y,l=1,2,3). Fourteen stress distributions along the x and y directions are applied onto the crack surface, as shown in Figure 6. The stress gradients are listed in Table 1. σ0 is assigned 400 MPa, and a is assigned 0.003 m in this paper. Stresses 1-5 change in the x direction as shown in Figure 6(a). Only these stress distributions are concerned by the universal weight function method. Stresses 6-9 change in the y direction as shown in Figure 6(b), while stresses 10-14 are the bivariate gradients shown in Figures 6(c)–6(g). The point weight function method can calculate the SIFs when there are cracks under these three types of stresses. The Πlk varies associated with two dimensionless geometric ratios a/L and a/T. To obtain the weight function coefficients when 0<a/T≤0.91 and 0<a/L≤0.91, 49 geometries are considered when a/T=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06 and a/L=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06.Figure 6
Stress distributions on the corner crack.
(a)
Stresses 1-5(b)
Stresses 6-9(c)
Stress 10(d)
Stress 11(e)
Stress 12(f)
Stress13(g)
Stress 14Table 1
Stress distributions applied onto the crack surface.
Stress numberxyStress distributions1x2σx,y=σ01−2x/a1−x/a2x3σx,y=σ01−x/a1−x/a1−x/a3x4σx,y=−σ01−4x/a1−6x/a1−x/a1−2x/a4x5σx,y=−σ01−4x/a1−2x/a1−1.05x/a1−0.9x/a1−0.8x/a5x6σx,y=−σ01−8x/a1−4x/a1−2x/a1−x/a1−0.9x/a1−0.8x/a1−0.6x/a+5/4σ06y1σx,y=σ0/41−6y/a7y2σx,y=σ01−y/a1+y/a8y3σx,y=σ0/801−6y/a1−6y/a1−6y/a9y4σx,y=σ01−y/a1−y/a1−2y/a1−3y/a10xyσx,y=σ01−2x/a1−6y/a11xy3σx,y=σ01−6y/a1−3y/a1−2x/a1−y/a12xy7σx,y=σ0x/ay/a713x2yσx,y=−σ01−8x/a1−4x/a1−2y/a+5/4σ014x2y3σx,y=−σ01−x/a1−0.5x/a1−4y/a1−2y/a1−y/a


## 4.2.2. Validation Results

TheΠlk varies associated with two dimensionless geometric ratios a/L and a/T. To obtain the weight function coefficients when 0<a/T≤0.91 and 0<a/L≤0.91, the surrogate model is established as Equation (10), and the surrogate model of the point weight function coefficients is shown in the appendix. The SIFs calculated by the obtained point weight functions with Πlk (surrogate model) are compared with finite element results to validate the accuracy of the weight function method. The SIF values are normalized by
(19)NormalizedSIF=Kσ0πa.These 49 geometries with different sizes are numbered as shown in Figure7(a). The SIF results under stresses 1-14 are calculated, and Figures 7(a)–7(l) show the results of stresses 1, 5, 6, 7, 10, and 14. The normalized SIFs calculated by the finite element method for point B and point A of the 49 geometries are marked with hollow and solid points, respectively. With the surrogate model, the SIFs are calculated and shown as three-dimensional surfaces (0<a/T≤0.91 and 0<a/L≤0.91) in Figures 7(a), 7(c), 7(e), 7(g), 7(i), and 7(k). Figures 7(b), 7(d), 7(f), 7(h), 7(j), and 7(l) show the relationship between the finite element results and the weight function method results of the 49 geometries.Figure 7
Comparison of SIFs on corner cracks for pointA and point B with different stress distributions.
(a)
3D SIF results under stress 1(b)
2D SIF results under stress 1(c)
3D SIF results under stress 5(d)
2D SIF results under stress 5(e)
3D SIF results under stress 6(f)
2D SIF results under stress 6(g)
3D SIF results under stress 7(h)
2D SIF results under stress 7(i)
3D SIF results under stress 10(j)
2D SIF results under stress 10(k)
3D SIF results under stress 14(l)
2D SIF results under stress 14Results show that the SIFs are determined by the geometries with the same loading. Whena/T and a/L are close to 1, the SIFs are 3-6 times larger than SIFs when a/T and a/L are close to 0. The validation result also shows that the SIFs calculated by the two methods are consistent when a/L and a/T are small.The relative errorCWFM between the weight function method and finite element method is defined as
(20)CWFM=KWFM−KFEMKFEM,where KWFM is the SIF results calculated by the weight function method and KFEM is the SIF results calculated by a finite element method. Similarly, the absolute error ΔKWFM between the weight function method and finite element method is defined as
(21)ΔKWFM=KWFM−KFEM.The specific errors between the weight function method and finite element method results of these 14 stresses are shown in Table2. The relative errors are larger than 100% when the absolute SIF value is small, less than 1 MPam. Most of the differences between the weight function method and the finite element results are less than 5 MPam when 0<a/L≤0.91 and 0<a/T≤0.91 except for stress 3, 10, 11, and 13 cases. The errors are related to the monotonicity of the stress distribution. Stress 1 can represent most stress concentration conditions that occur in rotor disks.Table 2
Errors between point weight function method and finite element results.
Stress typeMax errorAverage error (MAPE)ABABABABKUWF/KFEM−1KUWF−KFEMKUWF/KFEM−1KUWF−KFEM1127.27%6.92%6.2972.70137.92%1.84%1.2680.507226.33%3.44%3.9661.35710.59%1.05%0.8050.273357.53%30.86%10.9043.34915.87%6.44%2.2330.584456.23%27.46%2.3016.43721.00%7.55%0.4991.19358.06%10.17%8.5029.5902.64%2.95%1.6561.708625.26%1.44%0.0470.4001.90%0.91%0.0080.24571.74%5.59%1.2943.2950.62%2.54%0.2660.623858.93%31.64%2.00910.82022.86%7.60%0.4152.001920.56%45.80%5.7582.7466.26%17.74%1.1400.58710105.03%16.86%19.73620.49133.50%4.79%4.8034.2841187.32%59.31%15.2342.92719.47%13.92%3.1120.5321235.36%9.98%0.1250.16514.05%3.58%0.0220.0331352.64%33.00%14.6122.39612.20%11.98%2.9360.541145.32%6.75%0.6320.5424.23%2.15%0.1330.103The cumulative distributions ofCWFM and ΔKWFM are shown in Figures 8(a) and 8(b). Results show that 81% CWFM of point B and 54% of point A are less than 10%. As shown in Figure 8(b), for point B, 97% ΔKWFM are less than 5 MPam while 93% for point A. The absolute values of point B are usually larger than point A; thus, the relative error of point B is more minor.Figure 8
The errors of SIFs calculated by weight function method and finite element method.
(a)
Relative errors(b)
Absolute errorIn sum, the derived point weight functions can help obtain SIFs of corner cracks in disks under bivariant stress distributions. The relative error may increase when the absolute SIFs are less than 1MPam, and 81% of the relative errors for point B and 54% for point A are less than 10%. Meanwhile, 90% of the absolute errors between the point weight function and finite element method results are less than 5 MPam.


## 4.3. Error Analysis

### 4.3.1. Errors Induced by the Surrogate Model

Based on the finite element analysis mentioned in Section3.1, Forty-nine geometries with different sizes (a/T=a/L=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06) are established, and the associated discrete weight function coefficients Πlk are obtained by solving Equations (5). But to obtain the point weight function coefficients when 0<a/T≤0.91 and 0<a/L≤0.91, the surrogate model is established by the response surface method. The surrogate model will induce errors in the SIFs.The hexagonal polynomials (j=6) fit well, and the coefficient of determination (R-square) value reaches 0.9995. The stress distributions 1-14 mentioned in Table 1 are applied on the crack surface, and the SIFs of the 49 geometries are calculated. The relative error of the SIFs calculated by discrete Πlk and surrogate model is given as
(22)Cfit=Krps−KdisKdis,where Krps is the SIF results calculated by the surrogate model and Kdis is the SIF results calculated by discrete Πlk. The cumulative distribution of Cfit is shown in Figure 9. Results show that 93% Cfit is less than 1%. The response surface method helps obtain the surrogate model of weight function coefficients Πlk with relatively small errors.Figure 9
The relative errors of SIF calculated by discrete coefficients and surrogate model.


### 4.3.2. Errors Induced by the Reference Stress Intensity Factors

The accuracy of the reference SIFs directly ensures the effectiveness of theΠlk as shown in Equations (5). It is assumed that the three reference solutions for point AKr1A, Kr2A, and Kr3A deviate to the actual finite element results ranging from -20% to 20%. Thus, the three reference solutions are given as
(23)KrlAd=KrlA+αKrlA,where α is the offset scale and ranges from -20% to 20%.The proportion of the new coefficients derived by the offset reference solutions are defined as(24)γlAd=ΠlAdΠlA,where ΠlAd are the coefficients derived by the deviated reference solutions.Similarly, the proportion of the SIFs calculated by the new coefficients are given as(25)δlAd=,where KlAd are the SIFs calculated by the weight function method. Equations (23)–(25) also apply to point B.TheΠlk show a multivariate linear correlation of the offset scale in Figure 10 because the Πlk are solved by Equation (5). When the reference SIFs have a 1% error, the Πlk will change to 5 times of the original value. The SIFs calculated by the point weight function also change, as shown in Figure 11. When the reference SIFs have a 1% error, the SIFs for points A and B have an error up to 1.25%. The coefficients of the point weight function are highly dependent on the reference SIFs. Thus, the precise reference SIFs are significant to the point weight function method.Figure 10
The proportion of the coefficientsγlAd and γlBd with different offset scale.
(a)
PointA(b)
PointBFigure 11
The proportion of the coefficientsδlAd with different offset scale.The errors of the SIFs calculated by the point weight function method may be induced by the response surface method and the accuracy of the reference SIFs. First, the response surface fitting process of the discrete weight function coefficients causes the deviation of the SIFs. By fitting the coefficients with high order polynomial, the errors are less than 1%. Second, the accuracy of the reference SIFs calculated by the finite element method will influence the results of the weight function coefficients up to 5 times since these coefficients are calculated by solving simultaneous equations. A 1% error of the reference SIFs can cause a 1.25% deviation of the calculated SIFs.In addition to these two factors, the selection of the types of the reference stress distributions influences the accuracy of the point weight function. Since the uniform and linear reference stresses are used in this paper, the SIF results of linear distributions are accurate. Further research is necessary to study the influence of the selection of the reference loadings. Moreover, the point weight function equations can approximate Green’s function instead of equaling to Green’s function completely. The SIF errors are inevitably induced because the weight functions deviate from Green’s function, which is the limitation of the weight function method itself.In this case, during the crack propagation calculation and damage tolerance assessment, the SIF errors are supposed to be considered. Meanwhile, the development of the accurate finite element solutions is highly essential, which can provide the reference solutions to the point weight function method. Apart from the weight function method, a new surrogate model between the SIFs and the stress distributions of different geometries with cracks may be developed by other methods such as the machine learning method with enough SIF data.


## 4.3.1. Errors Induced by the Surrogate Model

Based on the finite element analysis mentioned in Section3.1, Forty-nine geometries with different sizes (a/T=a/L=0.91, 0.85, 0.75, 0.6, 0.3, 0.12, and 0.06) are established, and the associated discrete weight function coefficients Πlk are obtained by solving Equations (5). But to obtain the point weight function coefficients when 0<a/T≤0.91 and 0<a/L≤0.91, the surrogate model is established by the response surface method. The surrogate model will induce errors in the SIFs.The hexagonal polynomials (j=6) fit well, and the coefficient of determination (R-square) value reaches 0.9995. The stress distributions 1-14 mentioned in Table 1 are applied on the crack surface, and the SIFs of the 49 geometries are calculated. The relative error of the SIFs calculated by discrete Πlk and surrogate model is given as
(22)Cfit=Krps−KdisKdis,where Krps is the SIF results calculated by the surrogate model and Kdis is the SIF results calculated by discrete Πlk. The cumulative distribution of Cfit is shown in Figure 9. Results show that 93% Cfit is less than 1%. The response surface method helps obtain the surrogate model of weight function coefficients Πlk with relatively small errors.Figure 9
The relative errors of SIF calculated by discrete coefficients and surrogate model.


## 4.3.2. Errors Induced by the Reference Stress Intensity Factors

The accuracy of the reference SIFs directly ensures the effectiveness of theΠlk as shown in Equations (5). It is assumed that the three reference solutions for point AKr1A, Kr2A, and Kr3A deviate to the actual finite element results ranging from -20% to 20%. Thus, the three reference solutions are given as
(23)KrlAd=KrlA+αKrlA,where α is the offset scale and ranges from -20% to 20%.The proportion of the new coefficients derived by the offset reference solutions are defined as(24)γlAd=ΠlAdΠlA,where ΠlAd are the coefficients derived by the deviated reference solutions.Similarly, the proportion of the SIFs calculated by the new coefficients are given as(25)δlAd=,where KlAd are the SIFs calculated by the weight function method. Equations (23)–(25) also apply to point B.TheΠlk show a multivariate linear correlation of the offset scale in Figure 10 because the Πlk are solved by Equation (5). When the reference SIFs have a 1% error, the Πlk will change to 5 times of the original value. The SIFs calculated by the point weight function also change, as shown in Figure 11. When the reference SIFs have a 1% error, the SIFs for points A and B have an error up to 1.25%. The coefficients of the point weight function are highly dependent on the reference SIFs. Thus, the precise reference SIFs are significant to the point weight function method.Figure 10
The proportion of the coefficientsγlAd and γlBd with different offset scale.
(a)
PointA(b)
PointBFigure 11
The proportion of the coefficientsδlAd with different offset scale.The errors of the SIFs calculated by the point weight function method may be induced by the response surface method and the accuracy of the reference SIFs. First, the response surface fitting process of the discrete weight function coefficients causes the deviation of the SIFs. By fitting the coefficients with high order polynomial, the errors are less than 1%. Second, the accuracy of the reference SIFs calculated by the finite element method will influence the results of the weight function coefficients up to 5 times since these coefficients are calculated by solving simultaneous equations. A 1% error of the reference SIFs can cause a 1.25% deviation of the calculated SIFs.In addition to these two factors, the selection of the types of the reference stress distributions influences the accuracy of the point weight function. Since the uniform and linear reference stresses are used in this paper, the SIF results of linear distributions are accurate. Further research is necessary to study the influence of the selection of the reference loadings. Moreover, the point weight function equations can approximate Green’s function instead of equaling to Green’s function completely. The SIF errors are inevitably induced because the weight functions deviate from Green’s function, which is the limitation of the weight function method itself.In this case, during the crack propagation calculation and damage tolerance assessment, the SIF errors are supposed to be considered. Meanwhile, the development of the accurate finite element solutions is highly essential, which can provide the reference solutions to the point weight function method. Apart from the weight function method, a new surrogate model between the SIFs and the stress distributions of different geometries with cracks may be developed by other methods such as the machine learning method with enough SIF data.


## 5. Conclusions

This paper derives the point weight functions for corner cracks in finite plates to efficiently calculate SIFs when cracks are subjected to complex stress distributions induced by thermal, residual, and other types of stresses and conduct the probabilistic damage tolerance assessment. The conclusions are as follows:(1)
Based on the point weight function for the cracked structure, the SIFs for corner cracks can be obtained by integrating the two-directional stress in the uncracked body. In most cases (more than 90%), when the stress distributions are close to that of the rotor disks, the differences between the point weight function method and finite element results are less than 5MPam when 0<a/L≤0.91 and 0<a/T≤0.91. The relative error may increase when the absolute SIFs are less than 1 MPam, and 81% of the relative errors for point B and 54% for point A are less than 10%(2)
Obtaining the coefficients in the weight function, which are determined by the geometric parameters of the plates, is the critical process for the point weight function method. The response surface method is efficient for deriving the point weight function to reflect the relationship between geometric parameters and SIFs. The errors induced by the fitting method are less than 1%(3)
The reference stress intensity factors under associated loadings are essential input for the point weight function method because the unknown point weight function coefficients are calculated by solving simultaneous equations. The finite element method can help obtain the reference stress intensity factors. The accuracy of the reference SIFs will influence the results of the weight function coefficients. A 1% error of the reference SIFs can cause five times of deviation of the coefficients and a 1.25% deviation of the calculated SIFs. Thus, the accurate SIFs are the basis of the point weight function



---

*Source: 7575667-2022-07-01.xml*