Analysis model
To comprehensively assess the performance of the test joint, a corresponding finite element model was constructed. Initially, as depicted in Fig. 7(a), the joint model was established using the default ANSYS preprocessor. However, challenges arose in achieving a seamless transition between the main pipe and the branches. Consequently, the joint was accurately modeled using 3D modeling software SolidWorks to ensure consistency with the actual joint configuration, as illustrated in Fig. 7 (b).
The established joint model was imported into the finite element software ANSYS for analysis. The material properties used in the model were obtained from the material test. The elastic modulus of the material (E) was determined to be 2.0×10^5 N/mm², the yield strength (fy) was 235 MPa, and the Poisson's ratio (µ) was 0.3. The constitutive behavior of the material was selected to be the ideal elastic-plastic model, utilizing the Von-Mises yield criterion and associated flow rule for the elastic-plastic analysis [20]. To accurately simulate the real-world scenario, the boundary conditions of the joint were set as follows: the end part of the main pipe was fixed, while the ends of the branch pipes were fixed vertically. Additionally, the load was applied to the ends of the branch pipes in the form of surface pressure.
After extensive research and comparison, the three-dimensional solid element Solid65 was selected from the ANSYS element type library for modeling the joint. This element type features quadratic displacement and is well-suited for irregular grid division [21]. The finite element mesh of the joint, depicted in Fig. 8, was meticulously crafted to ensure accurate representation and analysis of the structural behavior.
Analysis results of the testing joint
The finite element analysis was conducted for the test joint under varying load conditions. Figure 9 illustrates the stress contours of the joint corresponding to load levels of 1000 kN, 2000 kN, 3000 kN, 4000 kN, and 5000 kN, respectively. From Fig. 9, it is evident that the overall stress level of the joint under a 1000 kN load is relatively low. The maximum stress value is predominantly concentrated in the vicinity of the joint core area. Furthermore, the stress observed in the main pipe and branch pipes is notably lower compared to that in the core area of the joint. Specifically, the stress in the former accounts for only 11.5% of the stress observed in the latter.
As the load increases to 2000 kN, the stress level of the joint exhibits a gradual linear increase, reaching a maximum stress value of 219 MPa. Upon reaching a load of 3000 kN, signs of yielding in the steel become apparent. However, the yield region is primarily concentrated at three points in the chamfer between the branch pipes. Subsequently, as the load further increases to 4000 kN, the yield region expands outward. This expansion manifests in two ways: firstly, the plastic region enlarges, and secondly, the chamfer between the main pipe and branch pipes also enters the plastic zone. Upon reaching the maximum load of 5000 kN, the expansion of the plastic zone within the joint intensifies, although it remains primarily concentrated in the vicinity of the joint core area. At this stage, nearly the entire core area of the joint enters the yield state, indicating the formation of a plastic hinge and marking the entry of the load-displacement curve into stage III (as illustrated in Fig. 6). Despite this, the stresses observed in the main pipe and branch pipes remain relatively low, approximately 106 MPa. This underscores the concentration of stress within the core region of the cast-steel joint with branches, which significantly influences its ultimate load-carrying capacity.
Figure 10 illustrates the vertical displacement of the joint under the maximum load of 5000 kN. It is observed that the maximum vertical displacement obtained from the finite element analysis (4.366 mm) is smaller than that obtained from the test (8.064 mm). This disparity can be primarily attributed to the non-tight contact between the joint and the test equipment piston during the experimental testing process.
Verification of the finite element model through the experiment results
To validate the numerical model of the cast-steel joint with branches, a comparison is made between the results of finite element analysis and those of the verification experiment. The stress values obtained from representative measuring points in the finite element model are compared with the experimentally derived stress values, and their relative differences are listed in Table 3.
Table 3
The joints analysis results with different parameters under axial loading
Joint
number
|
θ
(o)
|
L
(mm)
|
β
|
γ
|
R1
(mm)
|
R2
(mm)
|
R3
(mm)
|
Ultimate load-carrying capacity (kN)
|
J1
|
20
|
800
|
0.7
|
20
|
1000
|
0
|
0
|
2000.24
|
J2
|
30
|
800
|
0.7
|
20
|
1000
|
0
|
0
|
2301.79
|
J3
|
40
|
800
|
0.7
|
20
|
1000
|
0
|
0
|
2317.40
|
J4
|
50
|
800
|
0.7
|
20
|
1000
|
0
|
0
|
1660.00
|
J5
|
30
|
800
|
0.7
|
10
|
1000
|
0
|
0
|
5720.70
|
J6
|
30
|
800
|
0.7
|
15.2
|
1000
|
0
|
0
|
3344.50
|
J7
|
30
|
800
|
0.7
|
25
|
1000
|
0
|
0
|
1705.06
|
J8
|
30
|
800
|
0.7
|
29.9
|
1000
|
0
|
0
|
1334.73
|
J9
|
30
|
800
|
0.7
|
20
|
500
|
0
|
0
|
2056.27
|
J10
|
30
|
800
|
0.7
|
20
|
1500
|
0
|
0
|
2344.98
|
J11
|
30
|
800
|
0.7
|
20
|
2000
|
0
|
0
|
2233.56
|
J12
|
30
|
800
|
0.7
|
20
|
1000
|
10
|
0
|
2298.95
|
J13
|
30
|
800
|
0.7
|
20
|
1000
|
20
|
0
|
2268.51
|
J14
|
30
|
800
|
0.7
|
20
|
1000
|
30
|
0
|
2208.79
|
J15
|
30
|
800
|
0.7
|
20
|
1000
|
0
|
50
|
2824.84
|
J16
|
30
|
800
|
0.7
|
20
|
1000
|
0
|
100
|
3275.00
|
J17
|
30
|
800
|
0.7
|
20
|
1000
|
0
|
150
|
3296.32
|
J18
|
30
|
800
|
0.6
|
20
|
1000
|
0
|
0
|
2223.88
|
J19
|
30
|
800
|
0.66
|
20
|
1000
|
0
|
0
|
2281.90
|
J20
|
30
|
800
|
0.74
|
20
|
1000
|
0
|
0
|
2352.50
|
J21
|
30
|
800
|
0.8
|
20
|
1000
|
0
|
0
|
2418.69
|
J22
|
30
|
500
|
0.7
|
20
|
1000
|
0
|
0
|
2310.59
|
J23
|
30
|
600
|
0.7
|
20
|
1000
|
0
|
0
|
2312.61
|
J24
|
30
|
700
|
0.7
|
20
|
1000
|
0
|
0
|
2314.97
|
Analyzing the data in Table 3 reveals that:
1) The stress distribution of the cast-steel joint with branches obtained from the experiment aligns closely with that calculated from the finite element model. The calculated stresses from the finite element analysis and the experimental values of measuring points exhibit consistency. The maximum error between the calculated stresses and the experimental results is 9.02%, substantiating the validity of the finite element model utilized in this study.
2) Both the finite element modeling and the experiment confirm that the area of large stress is concentrated in the core area of the joint. The stresses in the main pipe and the branch pipes are comparatively small, approximately half of the largest stress observed in the core area of the joint.
3) The casting precision of cast steel joints presents challenges in control. In this study, it was observed that the chamfer between the main pipe and the branch pipes was slightly larger than the design value. Additionally, the wall thickness exceeded the design specifications, resulting in a smaller diameter thickness ratio. Consequently, the stresses predicted by the finite element model tend to be generally higher than those observed in the experiment.
4) In summary, the results obtained from the finite element model align closely with those from the experiment. The numerical model effectively captures the actual stress and deformation states of the cast-steel joint with branches. Therefore, it can serve as a reliable tool for investigating the load-carrying capacity of such joints in further research.
Effect of joint parameters on compression behavior of the cast-steel joint
To investigate the influence of different parameters on the ultimate load-carrying capacity of the joint, a parametric study is conducted wherein individual variables are varied while keeping other parameters constant. This approach allows for a systematic analysis of how each variable impacts the joint's performance independently.
The modeling results are shown in Table 4, which could be summarized as:
1) During compression testing, the fifth joint (J5) demonstrates the highest load-bearing capacity, reaching 5720.7 kN. Conversely, the eighth joint (J8) exhibits the lowest load-bearing capacity, registering only 1334.73 kN. Although both J8 and J5 share similar geometric characteristics, J8 boasts the greatest diameter thickness ratio (γ), while J5 possesses the smallest. Thus, the diameter thickness ratio (γ) significantly impacts the load-bearing capacity of joints.
2) Increasing only θ while holding other variables constant substantially decreases the joint's ultimate load-bearing capacity. This observation underscores θ's substantial influence on the joint's load-carrying capability.
3) Gradual increments in β and R3, while keeping other factors constant, substantially enhance the joint's ultimate load-carrying capacity. This finding highlights the significant impact of β and R3 on the joint's load-bearing capability.
4) When dimensions L, R1, and R2 undergo gradual increments while all other variables remain constant, the ultimate load-carrying capacity of the joint exhibits minimal variation, with the largest observed change being less than 5%. This indicates that dimensions L, R1, and R2 exert negligible influence on the ultimate load-carrying capacity of the joint.
5) Based on the findings of finite element modeling, it is deduced that thorough consideration of geometric parameters is imperative when analyzing the load-carrying capacity of cast-steel joints with branches. Careful selection of dimensional parameters for the joint is essential to ensure the structural safety and reliability.
Load-carrying capacity estimation of the three-branch cast-steel joint
In existing literature [19–20], load-carrying capacity formulas for welded tubular T-joints, steel tubular XK-joints, and multi-planar KX and KT-joints under axial loads are consistently represented as the product of the material yield strength and the square of the pipe wall thickness. Accordingly, the estimation of load-carrying capacity for cast-steel joints with branches can be similarly expressed as:
where Fu is load-carrying capacity of the joint; K is a parameter that contains the geometric parameters such as θ, γ and β of the joint; T is the pipe wall thickness; and fy is the material yield strength of the joint.
In Eq. (1), the expression of parameter K serves as the primary research focus across various types of joints. Given that K encompasses a range of geometric parameters affecting the load-carrying capacity of the joint, the focus has shifted from solely examining the relationship between K and material yield strength (Fu) to conducting multiple studies on the correlation between each individual parameter and K.
The finite element analysis results indicate that dimensions L, R1, and R2 of the joint exert minimal influence on the ultimate load-carrying capacity. Consequently, these parameters are disregarded during the analysis of the comprehensive index K. Utilizing line charts depicting the relationships between θ, γ, β, and R3 with K, as illustrated in Fig. 11, a regression analysis is performed.
Following the regression analysis, the relationship between θ and Kθ is initially examined. Through this analysis, the relationship between the sine value of θ and Kθ can be expressed as:
$${K_\theta }=0.60022 - 2.5311\sin \theta +6.59681{\sin ^2}\theta - 5.07388{\sin ^3}\theta$$
2
For the relationship between γ and Kγ, it could be expressed as a power function:
$${K_\gamma }=4.3725{\gamma ^{0.66242}}$$
3
From the observations in Fig. 10, it is apparent that a linear relationship exists between β and Kβ. This relationship can be expressed as:
$${K_\beta }=1+0.58856\beta$$
4
Finally, the regression analysis is performed between R3 and KR3. According to the principle of dimensional analysis, it is essential that the parameter R3 in the formula is dimensionless. To account for the influence of R3, a dimensionless chamfer coefficient ρ is defined as:
$$\rho =\frac{{{R_3}}}{{\sqrt {dt} }}$$
5
where d is the outer diameter of the branch pipe; t is the wall thickness of the branch pipe. The finite element model shows that the joint ultimate load-carrying ca-pacity is very small when R3 is greater than or equal to 100mm. So R3 is limited less than or equal to 100mm on the ultimate load-carrying capacity calculation formula for the cast-steel joint with three branches. Through regression analysis, the relationship between the chamfer coefficient ρ and KR3 follows a linear relation, which is expressed:
$${K_{{R_3}}}=1+0.33738\frac{{{R_3}}}{{\sqrt {dt} }}$$
6
Because these four parameters are independent of each other, the overall formula for the joint load-carrying capacity can be obtained by multiplying them, following the method of establishing the load-carrying capacity of joints in the existing standards, which is expressed as:
$$\begin{gathered} F=4.37251{\gamma ^{0.66242}}(1+0.58856\beta )(1+0.33738\frac{{{R_3}}}{{dt}})(0.60022 - 2.5311\sin \theta +6.59681{\sin ^2}\theta - \hfill \\ 5.07388{\sin ^3}\theta ){f_y}{T^2} \hfill \\ \end{gathered}$$
7
To validate the accuracy of Eq. (7), a comparison between the results obtained from finite element modeling and those derived from the regression formula is conducted. The comparative results are presented in Table 5. Notably, the disparity between the calculated values obtained from the formula and those from finite element analysis is minimal, with the maximum error amounting to only 1.9%. Consequently, it is deduced that the proposed formula effectively predicts the ultimate load-carrying capacity of the cast-steel joint with three branches with a high level of accuracy.
Table 5
Calculation formula error table
Joint number
|
θ
(o)
|
γ
|
β
|
R3
(mm)
|
F/
(fyo*T²)
|
Regression formula results
|
Difference percentage(%)
|
J5
|
30
|
10
|
0.7
|
0
|
9.74
|
9.92
|
1.90
|
J6
|
30
|
15.2
|
0.7
|
0
|
13.07
|
13.09
|
0.19
|
J2
|
30
|
20
|
0.7
|
0
|
15.67
|
15.70
|
0.21
|
J7
|
30
|
25
|
0.7
|
0
|
18.14
|
18.21
|
0.37
|
J8
|
30
|
29.9
|
0.7
|
0
|
20.37
|
20.50
|
0.65
|
J1
|
20
|
20
|
0.7
|
0
|
13.62
|
13.62
|
0.00
|
J3
|
40
|
20
|
0.7
|
0
|
15.78
|
15.78
|
0.00
|
J4
|
50
|
20
|
0.7
|
0
|
11.30
|
11.30
|
0.00
|
J18
|
30
|
20
|
0.6
|
0
|
15.14
|
15.05
|
-0.61
|
J19
|
30
|
20
|
0.66
|
0
|
15.54
|
15.44
|
-0.61
|
J20
|
30
|
20
|
0.74
|
0
|
16.02
|
15.97
|
-0.32
|
J21
|
30
|
20
|
0.8
|
0
|
16.47
|
16.36
|
-0.66
|
J15
|
30
|
20
|
0.7
|
50
|
19.23
|
19.09
|
-0.75
|
J16
|
30
|
20
|
0.7
|
100
|
22.40
|
22.47
|
0.32
|