1. Introduction
The design of the monopile is a pivotal factor in determining the maximum height of an offshore wind turbine. This design is critical for ensuring the safe, efficient, and economical operation of offshore wind turbines, significantly as they increase in height and are subjected to more demanding environmental conditions. It requires a balance of engineering, environmental science, and economics to achieve the best outcome for sustainable energy production [
1]. In the monopile design, tubular steel structures are essential due to their robustness and ability to withstand dynamic loads and vibrations. They offer advantages such as prefabrication off-site, requiring less on-site labor, and easier mounting of turbine components. Additionally, these structures allow wind turbines to be positioned at optimal heights, maximizing energy production from faster winds at higher altitudes [
2]. Overall, tubular structures play a crucial role in ensuring wind turbines’ safety, longevity, and energy efficiency for sustainable energy production. Innovation in tubular sections analysis is not only vital in the tower in the tower structure, but also in the support structure system like the monopiles (as shown in
Figure 1), jacket structures, or in the legs of jack-up barge arrangements. Hence, an optimum and simplified installation process of these support systems can also reduce the construction hassles and costs, ensure structural integrity, and enable innovation in wind energy projects [
3]. Because of its significance, “DNV-ST-0126 Support structures for wind turbines” provides general principles and guidelines for the structural design, installation, and in-service inspection of wind turbine support structure systems [
4].
When making decisions about energy system investments, it is essential to consider all aspects of planning, development, and construction. In general, it was found that the initial cost of a wind turbine system covered one-fourth of the total cost. Therefore, designers should prioritize building a robust and cost-effective structure in the planning phase [
5,
6]. A detachable tubular structure for a wind turbine tower offers several advantages in a support structure system. These benefits include faster and easier installation, reduced construction costs, and minimized environmental impact. Additionally, it enables improved transportation and assembly, increased robustness, better interaction with the wind, and cost savings that align with environmental considerations by reducing energy consumption and promoting sustainable practices. Furthermore, in times of emergency, the mobility of the structure gives engineers an advantage in preventing potential damage during natural calamities. However, due to limited studies, the potential of detachable joints in circular hollow section (CHS) structures still needs to be explored. Connections for hollow sections in steel construction are achieved mainly through welding or bolting. Welded circular hollow section joints are used for longitudinal joints and connecting braces to continuous chords, as Ummenhofer et al. (2002) [
7] used in their study.
On detachable tubular structures the most commonly used method is to use prestressed bolt in the flange of the CHS structure. The creation of L-type flanges is believed to be connected to the standardization initiatives for pressure vessels and piping systems during the Industrial Revolution. In the early 20th century, the American Society of Mechanical Engineers (ASME) [
8] spearheaded establishing and documenting flange standards, including L-type flanges, laying the groundwork for safety and quality in engineering practices. In a recent study, Tran et al. (2023) [
9] also developed a finite element model of a 5 MW L-type flange joint and verified its accuracy through experimental test results. They investigated the influences of geometrical imperfections (such as flange and tower-sided gaps) on the structural response using the same loading conditions and material properties. Another significant probabilistic study was conducted on the two-planar local joint flexibility by Ahmadi and Mayeli in 2018 [
10]. In their study, nine theoretical probability density function curves (PDFs) were prepared and fitted to the developed histograms. The maximum likelihood (ML) method was then applied to evaluate the parameters of the fitted PDFs. However, these connections are not easily detachable, so it is imperative to discover new hybrid connections that offer easier installation and mobility when needed. As a result, in a previous study, Islam et al. (2023) [
11] tried to develop a detachable joint named Pile-in-Pile (PIP) slip joint and explored the structural behavior under pure bending. For further development in this research, the response of a PIP slip joint subjected to bending and axial compression, considering perturbation, has been analyzed.
PIP slip joints are a one-of-a-kind new structural innovation that allows the contact surfaces of the CHSs to slide apart from each other, as demonstrated in
Figure 2. Some prior studies found that moment capacity, self-centering ability, and rotational stiffness are fundamental for analyzing the support structure system [
11]. Therefore, in the first step, the behavior of the PIP slip joint was investigated for a pure bending case to determine the optimum overlapping length for this CHS structure arrangement. For the next step, the combined loading of axial force and bending in wind turbine structures is crucial for several reasons. The combined loading of axial force (due to gravity and wind loads) and bending (due to wind-induced moments) affects the overall efficiency of the tubular structures. According to NORSOK N-004, tubular members subjected solely to axial tension, axial compression, bending, shear, or hydrostatic pressure should be designed to meet the strength and stability requirements. Tubular members subjected to combined loads without hydrostatic pressure should also be designed to meet the strength and stability requirements [
12]. By considering both axial and bending loads, engineers can also optimize the design to minimize material usage while ensuring structural integrity. Likewise, buckling is another critical failure mode for slender structures especially for CHS arrangements. Combined axial and bending loads influence buckling behavior and hence proper design accounts for buckling resistance need to be considered to prevent the possibility of failure. Wind-induced bending moments can also cause sliding between tower segments (e.g., flanges or slip joints), which can impact the overall tower stability. Analyzing these sliding distances ensures a safe operation and prevents excessive wear. Ideally, in the current study, the overlap will resist external bending, and two ring stoppers will be used inside the joint area to resist axial compression. Three different overlapping arrangements have been assessed thoroughly, while the dimensions of the ring stoppers remain unchanged. Local buckling near the overlapping area was also extensively studied as that part was found to be most interesting. After conducting a thorough analysis, we developed design equations for the entire system and its components. These will be compared with the well-established existing code.
3. Result and Discussion
After analysis, it has been determined that the overlapped area experiences significantly higher stress compared to other regions of the arrangement.
Figure 10 illustrates the equivalent stress for the entire arrangement, revealing a consistent stress distribution. However, localized buckling is evident in the overlapped area, with the maximum deformation occurring at the free end of the arrangement. Notably, a visible gap between two stoppers in the ring stopper area is causing a substantial concentration of stress in that specific area. Furthermore, the buckling effect has led to a noticeable change in the shape of the inner pile, as depicted in
Figure 11b.
The compression zone exhibited the highest stress, resulting in buckling in those regions. The maximum stress occurred at the interface of the inner and outer pile regions. When examining deformation, it was observed that the deflection was consistently varied from the support to the free end, as shown in
Figure 11a. However, a notable change in deformation was observed after the overlapping area. Notably, a dent was also present in the ring stopper area, where it connects with the shell of the slip joint arrangement.
A more accurate condition assessment was obtained through meticulous analysis of the in the overlapped area, where stress concentration is known to occur. Additionally, contact pressure has also become a concern with the deformation in the ring stopper. However, introducing separate stiffeners can effectively mitigate these issues and make this topic a compelling one to address in the near future.
Valuable insights into pile-to-pile interaction were gained through a thorough examination of the frictional stress at the contact surface. These findings offer significant insights into the impact of frictional stress on load-bearing structures and suggest potential strategies for optimizing the overall stability and load-bearing capacity of such arrangements [
28]. Additionally, they serve as a crucial indicator of the load-transfer mechanism from the outer to inner piles. The potential to substantially enhance lateral resistance exists by strategically managing this frictional stress in the future. In
Figure 12, the distribution of frictional stress is visually depicted, showing a higher concentration at the outer pile’s edge (labeled as “zone B”) where it intersects with the inner pile. A visible concentration of frictional stress at the edge of the inner pile, denoted as “zone A”, was also observed. Notably,
Figure 12b demonstrates a distinct behavior in the 2D overlapping length compared to the other two overlapping lengths in relation to frictional stress. The 2.5D and 1.5D overlapping lengths exhibit maximum frictional stress at 120 MPa and 116 MPa, respectively, while the 2D overlapping length peaks at 175 MPa. After a total displacement of 110 mm to 115 mm, a plateau in frictional stress was observed, followed by an increase up to 180 mm. This hike was caused by the very local deformation in the ring.
While examining a sliding ring joint, the sliding distance study is obviously crucial. Longer sliding distances between surfaces lead to increased wear and material loss. The resulting sliding-induced fatigue cycles could significantly impact the components’ lifespans in this slip joint configuration. This study plots the sliding distance between the ring stoppers of varying overlapping lengths against the total displacement. It is evident from
Figure 13b that the sliding pattern is almost identical for 2D and 2.5D overlapping lengths, while the behavior is notably different for 1.5D overlapping lengths. The increase in sliding for 2D and 2.5D overlapping lengths compared to 1.5D overlapping lengths is 15% and 25%, respectively. This variance may be attributed to the rotation in the overlapped area, consequently causing lateral pressure on the ring and potential deformation.
Consequently, the distortion in the ring joint was also checked. The deformed ring stopper is shown in
Figure 14a, and the pattern of deformation for different overlapping lengths is shown in
Figure 14b. From the figures, it can be observed that although the 2.5D and 1.5D overlapping lengths have the same pattern for ring deformation, the maximum value of deformation is 59 mm and 48 mm, respectively. The 2D overlapping length, however, showcases interesting behavior. It peaked at 62 mm with a total deformation of 120 mm. Subsequently, it decreased with a negative slope until a total deformation of 150 mm. Beyond this point, it ascended to 180 mm of total arrangement, with the maximum deformation in the ring reaching 78 mm. The observation suggests that after 150 mm of total displacement, the ring experienced local buckling, resulting in a loss of integrity. Therefore, the need to investigate this local buckling further in future studies remains essential.
In order to evaluate the axial load and moment capacity in the ring stopper area, this study considered the area where both of these parameters are crucial. After analyzing, it was found that the area lies in the region with overlapping lengths. To improve the understanding of the scenario, a model of a pile without any joints was also analyzed and compared with the result of pile slip joint arrangements.
Figure 15a,b present the axial load and moment capacity with respect to the normalized rotation.
Figure 16 shows the location of moment and resisting axial load in the PIP slip joint arrangement. To understand the behavior of the PIP slip joint based on bending and compression, this location is optimum. The axial compression and moment capacities were plotted against the rotational ratio to normalize it with respect to a dimensionless parameter. Surprisingly, in terms of axial load-carrying capability, the single pile section exhibited better results than the PIP slip joint arrangements. Among the slip joint arrangements, the 2.5D configuration demonstrated the maximum axial load-carrying capacity, with the axial compression capacity for 2D and 1.5D being almost identical. The axial load carrying capacity of 2.5D is nearly 5% more than the other two configurations, while in the case of the only pile, it exhibited a 9% resistance against the axial load. This might be attributed to stress concentration in a smaller region of the ring stopper area.
On the other hand, for moment capacity, the PIP slip joint has better moment capacity than the single pile section. Among the three different overlapping lengths, the 2D overlapping length provides better rotational resistance than the other two configurations. Following that, the 2.5D and 1.5D overlapping length arrangements have a greater capability of resisting rotation in chronological order. The rotational capacity for 1.5D, 2.5D, and 2D overlapping lengths is 4.77%, 36.02%, and 67.14% more than the single pipe case. The moment capacity increased significantly due to the change in stiffness around the overlapping area.
Based on this analysis the normalized compression-moment envelopes are presented in
Figure 17. Here, N
y and M
y represent the maximum axial load and the yield moment of the section. All the curves for the slip joint arrangements are nearly straight lines. To compare with the NORSOK code, a shaded area is also marked in graph. Within that area every line satisfies the design criteria proposed by the code. In all the cases, the lines lie within that area.
From the aforementioned analysis and discussion, it was found that the slip joint arrangement has to be designed to satisfy the following estimated equation.
In the “NORSOK N-004” code (Section 7.3) the equations are presented as follows [
12]:
Here,
- -
NSD and MSD are the design axial compression and design bending moment.
- -
Ny and My are the axial load capacity and yield moment, respectively.
- -
A and B are two dimensionless constants.
- -
D and t are the diameter and the thickness of the section.
- -
fy and E are the yield strength and the Young’s modulus.
- -
Ce and are the elastic buckling coefficient (0.3) and material factor, respectively.
- -
fCLE is the characteristic elastic local buckling.
Both equations are applicable for the following conditions:
In Equation (4) for different values of A and B the general design equation can be found for three different overlapping lengths. For different D/t ratios, these values will be different, and the P-M curves will be shifted parallel to the top and bottom of the current curves. In
Table 7 the values of these constants have been given for different overlapping lengths.
The members that are subjected to combined bending and axial compression need to satisfy the following condition as per Eurocode 3 (Section 6.3, Part 1-1) for a class-4 section [
20]:
Here,
- -
NED and My.ED are the design axial compression and design bending moment.
- -
NRk and My.Rk are the design value of the resistance depending upon the section.
- -
is reduction factor due to flexural buckling.
- -
is the interaction factor.
Based on the findings, it is evident that PIP slip joints outperform tubular sections without a joint in flexure. However, it is possible to develop the PIP slip joint using alternative methods for axial compression. Previous research has shown that in the case of pure bending, the 2D overlapping length outperforms the other two slip joint arrangements. This finding remains consistent when considering combined loading.
4. Conclusions
The current study explores structural analysis of the PIP slip joints, particularly focusing on stress distribution, deformation, and load-bearing capacities under axial compression and bending moments. The study observes that the overlapped area of the joint arrangement is prone to higher stress levels, leading to localized buckling and deformation.
Frictional stress at the contact surface is also being studied to understand how it affects the structural behavior and load-transfer mechanism. The distribution of frictional stress indicates higher concentrations at specific zones where the outer and inner piles intersect. The 2D overlapping length stands out for reaching a peak frictional stress of 175 MPa, and a plateau is followed by an increase in stress after certain displacement thresholds.
The study also analyzed the sliding distance between ring stoppers and found that longer overlapping is associated with higher wear and potential deformation. The sliding patterns of 2D and 2.5D overlapping lengths are similar, while the 1.5D length behaves differently. This difference is due to rotation in the overlapped area, which applies lateral pressure on the ring.
The deformation of the ring joint was also investigated and found that the 2D overlapping length exhibited a distinct pattern of deformation, reaching a peak of 62 mm then decreasing before rising again. This indicates the occurrence of local buckling and gap formation after a specific displacement point, resulting in a loss of structural integrity. Additionally, a general design equation has been estimated and compared with NORSOK and Eurocode for the same condition. The findings are significant for designing and analyzing PIP slip joints, offering a comprehensive understanding of their structural behavior under various loading conditions.
For further expansion of the current study, it is deemed to have an experimental analysis for validating the current study of PIP slip joint designs. Moreover, varying D/t ratios, load distances, and eccentricities are other focal points to be developed.