Next Article in Journal
Continuous Field Determination and Ecological Risk Assessment of Pb in the Yellow Sea of China
Previous Article in Journal
Ship Maneuvering in Shallow and Narrow Waters: Predictive Methods and Model Development Review
Previous Article in Special Issue
Hydrodynamic Interactions and Enhanced Energy Harnessing amongst Many WEC Units in Large-Size Wave Parks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multivariate Data Analysis of Maximum Stress Concentration Factors in FRP-Retrofitted Two-Planar KT-Joints under Axial Loads for Offshore Renewables

by
Esmaeil Zavvar
1,2,*,
Fernanda Sousa
1,3,
Francisco Taveira-Pinto
1,2 and
Paulo Rosa Santos
1,2,*
1
Department of Civil Engineering, Faculty of Engineering of the University of Porto (FEUP), 4200-465 Porto, Portugal
2
Interdisciplinary Centre of Marine and Environmental Research of the University of Porto (CIIMAR), 4450-208 Matosinhos, Portugal
3
Research Centre for Territory, Transports and Environment (CITTA), 4200-465 Porto, Portugal
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(8), 1451; https://doi.org/10.3390/jmse12081451
Submission received: 12 July 2024 / Revised: 16 August 2024 / Accepted: 18 August 2024 / Published: 22 August 2024
(This article belongs to the Special Issue Development and Utilization of Offshore Renewable Energy)

Abstract

:
With growing concerns about the danger of global climate change and worldwide demand for energy, the interest in the investigation and construction of renewable energy technologies has increased. Fixed platforms are a type of support structure for wind turbines composed of different types of tubular joints. These structures are under different kinds of cyclic loadings in ocean environmental conditions, which must be designed and reinforced against fatigue. In the present paper, the relationships between the parameters in DKT-joints reinforced with FRP under axial loads are investigated using several models, under 16 axial loading cases, with different nondimensional parameters and different FRP materials, and orientations were generated in ANSYS (total 5184) and analyzed. The four loading conditions that cause the maximum stress concentration factors were selected. After analyzing the 1296 reinforced models, relevant data were extracted, and possible samples were created. The extracted data were used in a multivariate data analysis of maximum stress concentration factors. The Pearson correlation coefficient is utilized to study the relationship between parameters and subsequently to make predictions. To reduce the number of variables and to group the data points into clusters based on certain similarities, hierarchical and non-hierarchical classifications are used, respectively.

1. Introduction

Tubular steel joints (TSJs), as seen in Figure 1, are common and usually the main components of the steel structures used in offshore structures, bridges, buildings, and other applications [1,2,3,4]. These TSJs, whether applied offshore or in onshore areas, are connected by welding, which is sensitive to the type of loads, especially cyclic loadings, like wind and waves [5,6,7], which can cause fatigue damage. Thus, to ensure structural integrity in these types of structures, the joints must be designed considering fatigue. The stress method [4,7,8,9] is mostly applied to calculate the fatigue resistance of TSJs. For this purpose, the SCF [10,11,12,13,14,15] is a key parameter for studying fatigue. The common way to predict the fatigue of joints is through empirical and numerical studies [16,17,18]. As an example, Bao et al. [19] worked on Y-three-planner-joints under axial loading and proposed equations to calculate SCFs. Rahmanli and Becque [20] worked on two-planner KT-joints with the help of SolidWorks under balanced axial loading. Their results revealed that among the different non-dimensional parameters, γ and τ are critical parameters for SCFs, so increasing them results in an increase in SCFs. In addition, the maximum SCF on the chord member took place in the sectors around 180° and 225° < φ < 270°. Ahmadi et al. and Zavvar et al. [21,22] worked on multiplanar KT-joints in ANSYS APDL 2024 to derive equations to calculate SCF values. They proposed a set of equations with higher R2, which indicates their accuracy. Kuang et al. [23] developed formulae to calculate SCFs in some types of joins, including the T/Y-, K-, and KT-joints. The authors did not model the weld profile thanks to the capability of the shell element used in their approach.
Some methods are used to reinforce steel structures and joints, such as internal and external steel ring plates [15,24,25,26,27,28,29,30,31,32], fiber-reinforced polymer (FRP) [33,34,35,36,37,38,39], doubler plates [40,41,42,43,44,45,46,47], and concrete and grout [48,49,50]. In FRP methods, Zavvar et al. and Hosseini et al. [51,52] performed an investigation on KT-joints against axial loadings and bending moments. More than 3000 FE models were created and analyzed in ABAQUS to derive formulae to calculate the value of SCFs. Zavvar et al. [53] conducted an investigation on uniplanar DKT-joints to calculate the maximum SCF using FE modeling. The solid element in ANSYS was used to create and analyze the models. SPSS was used to derive the formulae. The proposed formulae (with a high R2, i.e., >0.9) were considered accurate enough to calculate the maximum SCF. Zhao et al. [54] conducted a study on rectangular section joints reinforced with CFRP. The authors concluded that the crippling web capacity can be increased with CFRP.
In this research, finite element (FE) models of uniplanar KT-joints reinforced with FRP (which were previously verified in [16,55,56,57]) are used to extract data to calculate the maximum SCF. ANSYS v24 was utilized to generate the FE models and perform analysis under 16 axial loading conditions with different non-dimensional parameters (Table 1). Later, among all loading conditions, the four parameters leading to the maximum SCF were selected (details in Section 2). In the models with four loading conditions, FRP combined with five types of material was used. A collection of sample databases was generated in accordance with the findings of the FE investigation. SPSS and MATLAB were used to apply the statistical methods to the obtained data.

2. Models Characterization

In this section, the properties of the models and the modeling process are described. The models were created with different dimensions (e.g., diameters, length, and thickness). The combination of those dimensions created some well-known dimensionless parameters (Table 1). The selected specifications include a diverse range of tubular connectors that are often used in marine structures [58]. SOLID 186 and SHELL 281 were used to model the DKT-joints in ANSYS; the solid elements were used for modeling the weld profile and members whereas SHELL 281 was used for FRP. The main purpose of the FE modeling in ANSYS is to extract the HSS from the weld toe. Therefore, there is no need to have a fine mesh in all parts of the members, which would increase analysis time significantly. It is sufficient to have a fine mesh in the HSS region (International Institute of Welding (IIW) [59] recommendations) and weld profile. To this end, the sub-zone method is utilized to control the meshing quality (Figure 2).
In total, 16 possible loading conditions were applied to uniplanar KT-joints. The four loading conditions that create the maximum SCFs are presented in Table 2. A convergence investigation was performed on the meshing through the chord thickness and two elements were selected. Additional information is included in the references [12,13,53,60].
To model the weld profile, AWS [61] recommendations are used. An important parameter for the weld design is the dihedral angle (Ψ). It can be calculated as indicated in Figure 3 [22,50,62].
FRP is like a sheet of paper. Hence, the best option to model it is using shell elements. Members are modeled with solid elements, and FRP is modeled with shell elements. In order to simulate the interaction between the solid members and the FRP, an accurate method is using the contact capability of ANSYS. In this technique, the contact is applied to the joint outer surfaces (element contact 174), and the target is applied to the FRP (target 170). Regarding the length of the selected FRP, previous research [63] indicated that having an FRP longer than 6 D T for the chord members and 6 d t for the brace members is not necessary (Figure 4).
Based on previous studies [51,52], fiber orientations of 90° and 0° have been used. Table 3 provides information on the FRP materials. The third scheme was selected [51,52]. Additional information is provided in the references [12,13,53].
SCFs are calculate as follows [9,12,64]:
S C F = σ H S S / σ n
where σ n is the nominal stress. The hot-spot stress (HSS) is defined as the largest value of stress around the weld toe. The HSS is determined according to IIW [59]. Hence, HSS is achieved from two points, 0.4 T and 1.4 T (Figure 5). The nominal stress for these types of loading patterns is calculated as follows [12]:
σ n = F A
where σ n represents the nominal stress, F is the axial load, and A is the cross-section area of the loaded brace. The stress σ is determined as follows:
σ h s s = 1.4 σ 1 0.4 σ 2
where the 1st point is represented by σ 1 , and the 2nd point is indicated by σ 2 .

3. Statistical Procedures

The statistical analysis of the FRP-retrofitted two-planar KT-Joint under AX loading is described in this section. The models were generated considering seven parameters, namely the joint traditional non-dimensional and FRP parameters (β, γ, τ, θ, η, and ξ) and the resulting SCFs. The processing of sample data involves a quantitative statistical analysis aimed at determining global indicators across six dimensions. This procedure allowed for the creation of five sample databases that present the values of SCF for this type of joint. Then, using software such as SPSS v29 and MATLAB v24, the statistical analysis was conducted. Finally, the samples were submitted to detailed statistical procedures and methods.

3.1. Applied Methods

The statistical methods applied to the created datasets focus on six quantitative indicators. In the first step, with the help of univariate analysis, each parameter is described individually using some graphical representations, such as box plots. In the second step, the bivariate analysis, Pearson’s correlation coefficient is applied to analyze the linear relationships between the variables. Finally, in the multivariate analysis, two approaches, including multiple regression [65,66,67] and hierarchical classification (HC) [68,69] are used to create classes of data and classes of indicators.
HC, often referred to as Agglomerative Hierarchical Clustering (AHC), is a method used to classify a set of elements by grouping them into a smaller number of classes. This classification is based on the principle that elements within the same class share similarities with each other and exhibit differences from elements in other classes. Unlike some other clustering methods, the number of classes in AHC is not predetermined but is determined as part of the clustering process itself. This iterative approach builds the classes gradually, merging similar objects into clusters until a stopping criterion is met, resulting in a hierarchical structure of clusters.
The hierarchical classification (HC) method involves two critical decisions. The first one being the determination of the measure of comparison between the pairs of elements to be classified. This choice depends on whether the elements being classified are variables (similarity type) or individuals (dissimilarities). The second consists of selecting the criterion to measure the proximity between two classes. It is essential to recognize that these two choices are pivotal factors as they can significantly influence the outcome of the clustering process. Careful consideration of these choices is crucial since they play a vital role in shaping the final clustering results.

3.2. Univariate Analysis

The models were analyzed to obtain SCFs, which were then categorized into five distinct samples for further analysis. The 1st sample comprised the SCFs from the 1st loading condition, while the 2nd sample included SCFs from the 2nd loading condition. Similarly, the third and fourth samples consisted of SCFs from the 3rd and 4th loading conditions, respectively. The fifth sample contained the maximum SCFs of all the 16 loading conditions. Table 4 displays, for each sample, the SCFs’ descriptive measures. A positive skewness indicates that the distribution tail is skewed towards higher values; therefore, all samples show positive asymmetry. Moreover, the kurtosis value in the third sample surpasses that of the other samples, suggesting a probability distribution with a pronounced peak. In contrast, the kurtosis value in the fourth sample is lower than that of the other four, indicating a probability distribution with a less pronounced peak for this sample. In this research, the number of the considered parameters is 6 (Table 1).
Thus, looking at the 5 samples in Table 4, sample 1 and sample 5 have higher values: sample 1 has an average of 20.51, a variance of 110.1, a mode of 13.42, a standard deviation of 10.50, and a maximum of 56.76; sample 5 has a median of 18.02. Sample 2 (the second loading condition) has a skewness of 0.85 and the maximum is 34.49.
Figure 6 displays box plots illustrating the distributions of sample 1 through sample 5. Each box plot is composed of the minimum, first quartile, median, 3rd quartile, and maximum, dividing the distribution into four intervals, each showing 25 percent of the respective distribution. This set of box plots enables comparisons between sample 1 and sample 5 within each dimension, as well as comparisons of indicators across dimensions. The samples demonstrate significant concentrations in the lower range of the distribution. (positive asymmetry). The sample 1 and sample 5 box plots reveal an accentuated negative asymmetry of the distribution. Sample 1 and sample 5 have a significant concentration between the first quartile and the median. Other samples (2, 3, and 4) are symmetrical. Upon comparing the samples, it can be inferred that sample 1 and 5 exhibit the widest range of values, while sample 4 has the narrowest range. The sample 4 box plot displays a few lower–moderate outliers. Figure 6 also shows the behavior of the samples. Sample 5 has the widest range of values, and sample 4 has the lowest. Furthermore, sample 5 and sample 1 have values that are greater than the values of the other samples.

3.3. The Histograms

In order to construct a density histogram, it is necessary to partition the range (R) into many classes. There are several suggestions to calculate the number of classes, such as the Sturges [70] and Freedman–Diaconis (FD) [71,72,73,74] method.
In this study, the FD [71] rule was used to calculate the number of classes. It is determined as follows:
n c = 2 ( i n t e r q u a r t i l e   r a n g e ) n 3
Table 5 indicates the FD values for each sample. Figure 7 indicates the histograms with normal distribution. According to the skewness and kurtosis values (Table 5), it can be inferred that the histograms exhibit a larger right tail compared to the left tail. Figure 7 indicates that the values of the kurtosis are higher than 3 for four samples, including 1, 2, 3, and 5 (Table 4), hence suggesting a leptokurtic distribution thinner than the standard normal distribution. Only for sample 4 is the kurtosis value smaller than 3, indicating a platykurtic distribution.

3.4. Bivariate Analysis

This section addresses the bivariate analysis, which involves examining the relationship between two variables. Measures of association are very common for this analysis, and one of the most popular methods is the Pearson correlation coefficient (PCC), which can involve two or more variables. Correlation simply measures the linear association between variables without any implication of cause and effect. The sign of the correlation coefficient indicates the direction of the relationship, a positive correlation signifies that the variables fluctuate in the same direction, while a negative correlation signifies that they vary in opposite directions. The Pearson coefficient is defined as follows:
R ( a , b ) = c o v ( a , b ) v a r a × v a r ( b )
Therefore, the Pearson correlation coefficient is a standardized covariance between −1 and +1. When it is close to 1 it indicates the best positive correlation. In this section, instead of presenting all sample results, just the results of sample 5 are presented for the sake of brevity.
In the bivariate analysis of the samples, Pearson’s correlation determines correlation values for samples and parameters (Table 6 and Table 7). This analysis indicates that γ values are correlated with the η values in each dimension and also between them. The γ-SCF and τ-SCF correlation values are higher than for the others. The correlation between the β-SCF, θ-SCF, and ξ-SCF is weak between all parameters, while the η with other indicators is the lowest.
Table 6 and Table 7 indicate that parameters γ and τ have the best correlation with the SCF and that parameters γ and ξ have a high correlation together. Furthermore, it shows that, among other parameters, these three are important and have significant roles in predicting SCF. The lowest correlation is related to θ; hence, it means that this parameter has no relevant effect on the structure’s behavior against fatigue. Therefore, the thickness of the braces and the chord are much more important than the other variables in protecting tubular joints against fatigue.

4. Scatter and P-P Plots

Figure 8 indicates P-P and scatter plots of the samples. P-P plots (Figure 8 left) are typically used in regression analysis to assess the assumption of normality of the residuals. The X-axis represents the cumulative probability of the observed residuals, and the Y-axis represents the cumulative probability that would be expected. Deviations from the line indicate departures from normality. Large deviations can suggest problems with the model or the need for transformation of the dependent variable. The scatter plots (Figure 8, right) provide the relationship between the regression standardized residuals and the SCF (for all samples) values. The X-axis indicates the regression standardized residuals, whereas the Y-axis indicates the SCF values. Scatter plots with independent and dependent variables are presented. The joint parameters (β, γ, τ, θ, η, and ξ) are the independent variables, whereas SCF is defined as the dependent variable. In all samples, it can be seen from P-P plots (Figure 8 left) that the points closely follow the diagonal line, which presents that the data are approximately normally distributed. Any major deviations from the line might indicate a problem with normality, but in plots, the points seem to be well-aligned with the line, suggesting that the normality assumption is acceptably met.
The scatter plots (Figure 8, right) show that, for example, in sample 5, distribution seems to be a fan-shaped pattern, suggesting heteroscedasticity (variance of residuals increases with the predicted values). In the lower range of the SCF (up to about 20), the data are more tightly clustered around zero, but as the SCF increases, the spread of residuals increases, forming a distinct pattern. There are clusters of points, particularly in the mid-range of SCF values, indicating that certain ranges of SCF values have more residual variation than others. The presence of outliers can be seen at both ends of the SCF values, with some points far from the main cluster of data.

5. Multiple Regression

In this section, the multiple regression analysis is explained with the joint geometry parameters, FRP, and SCF. Several key parameters are important in the multiple regression analysis, such as R2 and adjusted R2 which can be calculated as follows:
R 2 = 1 i = 1 n ( b i b ^ i ) 2 i = 1 n ( b i b ¯ i ) 2
A d j u s t e d   R 2 = 1 ( 1 R 2 ) ( n 1 ) n k 1
where b i represents the actual values of the dependent variable, b ^ i represents the predicted values from the model, and b ¯ i represents the mean of the actual values; n indicates the sample size, and k presents the number of predictors.
Table 8 indicates the summary information of the samples. Samples 1 and 5 both have the highest R2 values of 0.903, indicating that these models explain 90.3% of the SCF. Sample 4 has the lowest standard error of the estimate (2.03349), suggesting the predictions for this sample are the closest to the actual values. All samples show high R and R2 values, suggesting that the samples fit the data well. Sample 5 is the most important sample, presenting the SCFmax in FRP-reinforced two-planar KT-joints subjected to AX loads. The results indicate that the model performs well across all samples, with strong correlations and high R2 values, though the accuracy of predictions (as measured by the standard error) varies somewhat between samples. It shows that samples 1 and 5 demonstrate the best overall performance with the highest R and R2 values.
Table 9 provides the variance analysis (ANOVA). The ANOVA table details the analysis of variance for five different regression models (Sample 1 to Sample 5). The sum of squares (SS) (regression) represents the variability explained by the regression model, and a higher value presents that the model explains a significant portion of the total variability. The sum of squares (residual) represents the variability not explained by the model, and a lower value shows that the model fits the data well. The degrees of freedom (DoFs) for the regression are 6, and for the residuals, there are 1289. The mean square (MS) is the sum of squares (SS) divided by the respective degrees of freedom. It is used in the calculation of the F-statistic. A higher F-value proves that the model is significantly more appropriate at determining the outcome than a model with no predictors. A value of <0.001 demonstrates that the regression model is significant.
Table 10 proves that all models have a p-value (sig.) less than 0.001 and that the predictors collectively have a considerable influence on the SCF. The F-statistics are very high for all samples, further confirming the models’ overall significance. Higher F-values indicate that the model explains a substantial proportion of the variability in the SCF. For example, samples 1, 2, and 3 have a substantial proportion of the variability (F = 1993.656, p < 0.001 and F = 1573.234, p < 0.001). Samples 3 and 4 are significant with an F-value of 1080.158 (p < 0.001) and F = 752.596 (p < 0.001). The consistently high F-values and low p-values across all samples indicate that the regression models are effective in explaining the variability in the dependent variable for each sample.
The standardized coefficients are shown in Table 10. Standardized coefficients help to evaluate the effects of the other materials on the SCF. All the p-values are less than 0.001, indicating that all predictors are statistically significant at the 0.05 level. The variable τ has the highest standardized coefficient, suggesting it has the strongest effect on the dependent variable. The variable β has a negative influence on the dependent variable. For example, in sample 1, β, η , and ξ are negative, indicating a negative relationship with the dependent variable. τ has a high value, indicating a strong positive relationship. γ and θ are positive, showing a positive relationship. Based on the Beta values in Table 10, it can be seen that two parameters, including γ and τ, have considerable effects on the SCF. The lowest effects for variables belong to the inclination angle (θ). For example, the value of the θ in samples is close to zero, especially in samples 2, 3, and 4.

6. Hierarchical and Non-Hierarchical Classification

The importance of clustering in multivariate statistical analysis is that it allows for the identification of patterns in the data that may not be readily apparent from the raw data and to group similar data points together, reducing the dimensionality of the data and making it easier to interpret and visualize. In this case, considering the huge dataset available, with more than 1296 observations, this classification is essential for the multivariate analysis of these data.
The classification of observations in each of the clusters is, in general, more rigorous in non-hierarchical methods; however, it is advisable in a cluster analysis problem to start with hierarchical methods for the purpose of exploration and then proceed with the non-hierarchical, in the case of k-means, to refine and interpret the cluster solution. In this work, the hierarchical method was not fully successful due to the size of the dataset, making it impossible to apply this method to the entire set, but instead only to a minimal part, for the learning and understanding criteria of the method. The k-means method was able to be applied to the entire dataset, resulting in a very useful analysis.
A good distance method is one that results in a clear separation of groups in the dendrogram, with similar samples forming a single branch and different samples forming separate branches. The Euclidean distance (it is known that the comparison criterion has more effect on the classification than the distance measure used) was used, and the results are presented in the dendrograms below (Figure 9 and Figure A1). This type of graph makes it possible to easily understand the connections made between individuals and identify those who are more similar, which are used as the final result obtained via this hierarchical classification system.
In this case, the first observations show a peculiar behavior, which is the clusters of housing pairs with indices very close to each other. For the sake of brevity, only data from sample 5 are represented. When comparing the dendrograms resulting from the four comparison criteria mentioned, it can be observed that only the two formed via ward and single linkage present a tree that is quite different from the others. Therefore, the method of average distance between clusters, or average between groups, was used with the Euclidean distance measure. The proximity matrix, Table 11, reveals the dissimilarity between the cases, while the agglomeration scheme indicates the order of aggregation of individuals in the respective clusters. The first cluster to be formed contains individuals with indexes 1153 and 1156, which have a smaller average distance between them. In the second step, the cluster is formed with subjects 901 and 904, and so on.
The dendrogram was omitted due to its long length. In the initial stages, clusters are combined at very low coefficients (0.1), indicating they are quite similar or close together. The “Coefficients” column contains the proximities shown in Table 11. As the process continues, these coefficients increase, reflecting the growing dissimilarity between the clusters being combined. As we move down the table to higher stages, the coefficients increase, indicating that clusters being combined are becoming less similar. For example, by stage 1294, the coefficient is 17.291, and by stage 1295, it increases significantly to 29.821. The last stages (1293–1295) involve the combination of larger clusters, resulting in significantly higher coefficients. This suggests that in the final stages, the clusters being merged are quite distinct from one another. Table 11 provides a detailed view of the hierarchical agglomerative clustering process, showcasing how clusters are incrementally combined from initial stages with low dissimilarity to final stages where the clusters are significantly different. This progressive increase in coefficients reflects the growing differences between clusters as the clustering process consolidates more diverse groups.
It can be seen in Table 11 that most of the combined clusters were the result of close index pairs, with a few exceptions. This is due to the previously mentioned fact that data are very similar from one time to the next, and therefore tend to be grouped together.
(a) 
Non-hierarchical cluster (NHC) grouping
NHC grouping approaches are intended to exclusively group individuals, or observations, into a set of clusters where the number is defined by the analyst. This is very useful when working with very large datasets. The method used for this classification is the k-means method, which starts from a previously defined number of classes (in this case, the number of classes used is that resulting from the application of the previously presented hierarchical classification, and the respective cutoff of the dendrogram in k classes, that is, 16 classes), and each of the classes is assigned an individual who will function as the center of the respective class. Table A1 presents the average of each variable in each of the 10 clusters.
The analysis carried out with SPSS, using 10 clusters with 10 iterations, proved to be reasonably appropriate for the purpose of this work. However, for the purpose of comparison, a second analysis was carried out using the MATLAB programming language, which allows for more iterations without compromising the computer’s memory, like SPSS. Therefore, an analysis was carried out with the same clusters but with 300 iterations, as seen in Table 12.
Table 13 shows the variation in the center of the clusters at each step of the iteration. It would be possible to increase the number of iterations, but this requires a high-performance computer. The parameters like β, τ, and γ have more uniform values across clusters, suggesting these are more consistent features. Other parameters like η, ξ, and SCF show more variability, indicating these are more distinctive features that differentiate the initial clusters. The values show the initial conditions or centers from which the clustering process starts. Some variables have a wide range (e.g., SCF from 5.12 to 56.76), indicating heterogeneity among initial clusters. Others are more constrained (e.g., β between 0.4 and 0.7), suggesting some homogeneity.
Table 14 shows the number of cases assigned to each of the 16 clusters after the initial clustering process. The total of 1296 cases being valid and with none missing indicates the completeness and integrity of the dataset used for clustering. The number of cases per cluster ranges from as few as 14 cases (Cluster 5) to as many as 147 cases (Cluster 10). There is considerable variability in cluster sizes, indicating a diverse spread of data across clusters. The largest cluster, Cluster 10, has 147 cases, making it the most populous, but the smallest cluster, Cluster 5, has 14 cases, making it the least populous. Clusters with a large number of cases (e.g., Cluster 10 and Cluster 9) may represent more common patterns or groupings within the data. Clusters with a small number of cases (e.g., Cluster 5 and Cluster 6) might represent outliers or less common groupings, which could be significant depending on the context. The presence of clusters with moderate sizes (e.g., Clusters 1, 2, 3, 4, 7, 8, 11, 12, 13, 14, and 15) suggests a balance between common and less common patterns.
It can be seen from Table 14 that θ, τ, γ, β, η, and ξ represent the average of the respective features within each cluster. For example, for cluster 0, the average θ is 0.5236, τ is 0.4, etc. SCF shows the average SCF (presumably some important metric or feature) for each cluster. For instance, cluster 0 has an average SCF of 7.9351. The SCF values show significant variation across clusters, indicating this feature plays a crucial role in differentiating the clusters. For instance, cluster 4 has the highest SCF value (52.8629), while cluster 0 has the lowest (7.9351). The centroids provide a summary of each cluster’s central tendency, highlighting the average values for each feature within the clusters. This information helps to understand the characteristics and differences between clusters.

7. Conclusions

A total of 5184 analyses were carried out on two-planar KT-joints under 16 AX loadings, and 1296 models were selected and analyzed under four loading conditions which created maximum stress. In the first stage, an FE analysis was conducted, and then a univariate analysis was carried out, dealing with the descriptive statistics of each variable, making it possible to identify important information such as measures of dispersion and location, such as mean, median, and standard deviation, as well as histograms to verify the distribution of their values and box plots to identify possible outliers and distribution of quartiles. A bivariate analysis was then carried out, where the correlations between the variables were identified using the Pearson coefficient and also multiple forms of regression to verify the degree of linearity between pairs of variables.
Among three angles, including 0°, ±45°, and 90°, for FRP orientation, the effective fiber orientations are 90° and 0° in the chord, while FRP orientations on the brace have no effects on chord SCFs. The highest and lowest SCFs were 56.76 and 0.012, respectively. The SCFmax was located at the saddle point of the central brace under the 1st loading condition.
In the multivariate analysis stage, the analysis of the main components was first carried out, where it was possible to decrease the size of the dataset into two variables that contained more than 80% of the initial information, proving to be an efficient approach to reducing size in this case, a variable housed the direction variables, and another housed the height and period variables.
Bivariate analysis shows that θ has the lowest correlation and τ and γ have the highest correlation among others. Hence, it means that θ has no effects on the structure’s behavior against fatigue. However, the thickness of the braces and the chord are much more important in tubular joints against fatigue.
Hierarchical and non-hierarchical classification analyses of the observations were carried out. The first proved to be incapable of being carried out on such an extensive set of data, compromising the computer’s memory, but based on knowledge of the methodology, dendrograms and tables were created using a smaller subset of data. The non-hierarchical approach, using the k-means clustering technique, proved capable of being applied to a large dataset and to be quite efficient in grouping representative SCF samples.
Multivariate data analysis provides a theoretical and practical framework for studying complex data by examining relationships between multiple variables simultaneously, by reducing the number of variables

Author Contributions

Conceptualization, E.Z. and F.S.; methodology, E.Z. and F.S.; software, E.Z.; validation, E.Z.; formal analysis, E.Z.; investigation, E.Z.; resources, E.Z.; data curation, E.Z.; writing—original draft preparation, E.Z.; writing—review and editing, E.Z., F.S. and P.R.S.; visualization, E.Z. and F.S.; supervision, P.R.S. and F.T.-P.; project administration, P.R.S. and F.T.-P.; funding acquisition, P.R.S. and F.T.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Portuguese Foundation for Science and Technology (Fundação para a Ciência e Tecnologia-FCT), grant number 2023.01289.BD.

Institutional Review Board Statement

Not applicable for studies not involving humans or animals.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in this study are included in this article, further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

θAngle between outer braces and the chord
AWSAmerican Welding Society
AXAxial load
dBrace diameter
lBrace length
tBrace thickness
τBrace thickness to the chord thickness t/T
DChord diameter
γChord diameter to twice the thickness of chord D/(2T)
LChord length
TChord thickness
βDiameter of brace to diameter of chord d/D
EcElastic modulus of concrete
FRPFiber-reinforced polymer
FEFinite element
FEMFinite element model
HSSHot spot stress
IIWInternational Institute of Welding
IPBIn-plane bending
OPBOut-of-plane bending
SCFStress concentration factor
τThe brace thickness to the chord thickness t/T
ξThe ratio of Es to EFRP
EsYoung’s modulus for steel

Appendix A

Figure A1. Hierarchical cluster grouping types.
Figure A1. Hierarchical cluster grouping types.
Jmse 12 01451 g0a1
Table A1. Change in cluster centers.
Table A1. Change in cluster centers.
12345678910111213141516
11.8791.0360.5511.8732.2420.4500.9500.4530.5861.9943.0390.8090.7900.9260.9200.768
20.7000.0000.5020.3601.0981.0880.1500.6870.6150.4301.1800.0800.1230.2990.0850.514
30.3100.0000.5760.2100.5610.8620.0970.2610.1590.1100.7210.0540.0790.2770.0280.596
40.1750.0400.3990.1350.0000.0710.3420.0890.0940.0490.2890.0550.0770.1410.1420.457
50.0990.0810.3460.1110.0000.1690.1350.0000.0690.0000.2830.0560.0750.1200.1120.200
60.0960.0800.2680.0650.0000.0870.0890.0180.0860.0000.0690.0290.0370.0760.1740.066
70.0620.0400.2070.0000.0000.0830.0440.0360.0810.0000.2050.0860.1060.0000.0570.068
80.0310.0000.1630.0000.0000.0000.0890.0550.0740.0000.1290.1530.1630.0000.0580.134
90.0600.0000.0450.0000.0000.0000.0450.0580.0960.0000.0640.0960.0940.0000.0000.065
100.0300.0800.0000.0000.0000.0000.0880.0190.0390.0000.0000.0650.0620.0000.0000.000
110.0000.0410.0000.0000.0000.0000.0870.0380.0340.0000.0000.0660.0600.0000.0000.065
120.0310.0000.0000.0000.0000.0000.0450.0450.0390.0160.0000.0330.0300.0000.0000.063
130.0000.0000.0000.0000.0000.0000.0000.0630.0330.0170.0000.0000.0000.0000.0000.000
140.0000.0000.0000.0000.0000.0000.0000.0850.0490.0160.0000.0000.0000.0000.0000.000
150.0000.0000.0000.0000.0000.0000.0000.0410.0320.0000.0000.0000.0000.0000.0000.000
160.0310.0000.0000.0000.0000.0000.0000.0000.0220.0000.0000.0000.0000.0000.0000.000
170.0000.0000.0000.0000.0000.0000.0000.0420.0320.0000.0000.0000.0000.0000.0000.000
180.0310.0000.0000.0000.0000.0000.0000.0210.0370.0000.0000.0000.0000.0000.0000.000
190.0000.0000.0000.0000.0000.0000.0000.0210.0160.0000.0000.0000.0000.0000.0000.000
200.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.000

References

  1. Ziaeinejad, A.; Tousignant, K.; Sun, M.; Yang, Z. SCFs for CHS-to-CHS T- and X-type moment-loaded end connections: Proposal for formulae. J. Constr. Steel Res. 2024, 218, 108734. [Google Scholar] [CrossRef]
  2. Nejad, A.Z.; Sun, M.; Tousignant, K. Circular hollow section X-connections near an open chord end: Stress concentration factors. J. Constr. Steel Res. 2021, 177, 106454. [Google Scholar] [CrossRef]
  3. Ahmadi, H.; Ziaei Nejad, A. Geometrical effects on the local joint flexibility of two-planar tubular DK-joints in jacket substructure of offshore wind turbines under OPB loading. Thin-Walled Struct. 2017, 114, 122–133. [Google Scholar] [CrossRef]
  4. Ziaeinejad, A.; Tousignant, K.; Sun, M.; Yang, Z. Stress concentration factors for RHS-to-RHS T-type moment-loaded end connections. Thin-Walled Struct. 2024, 203, 112179. [Google Scholar] [CrossRef]
  5. Zavvar, E.; Chen, B.-Q.; Uzunoglu, E.; Guedes Soares, C. Stress distribution on the CENTEC-TLP in still water and rated wind speed. In Trends in Maritime Technology and Engineering; Taylor & Francis Group: London, UK, 2022; pp. 519–526. [Google Scholar]
  6. Zavvar, E.; Abdelwahab, H.S.; Uzunoglu, E.; Chen, B.Q.; Soares, C.G. Numerical study of the wave induced motions and loads on the CENTEC-TLP floating wind turbine. In Trends in Renewable Energies Offshore; CRC Press: Boca Raton, FL, USA, 2022; pp. 567–573. [Google Scholar]
  7. Zavvar, E.; Abdelwahab, H.S.; Uzunoglu, E.; Chen, B.-Q.; Guedes Soares, C. Stress Distribution on the Preliminary Structural Design of the CENTEC-TLP under Still Water and Wave-Induced Loads. J. Mar. Sci. Eng. 2023, 11, 951. [Google Scholar] [CrossRef]
  8. Chang, E.; Dover, W.D. Prediction of degree of bending in tubular X and DT joints. Int. J. Fatigue 1999, 21, 147–161. [Google Scholar] [CrossRef]
  9. Zavvar, E.; Ahmadi, H. Parametric study of stress concentration factors in three-planar tubular KT-joints subjected to out-of-plane bending loads. In Proceedings of the International Conference on Coasts, Ports and Marine Structures (ICOPMAS 2016), Tehran, Iran, 31 October–2 November 2016; Ports & Maritime Organization: Tehran, Iran, 2016. [Google Scholar]
  10. Recho, N.; Brozzetti, J. Stress concentration at tubular V-joints. In Welding of Tubular Structures, Proceedings of the Second International Conference, Boston, MA, USA, 16–17 July 1984; Pergamon: Oxford, UK, 1984; pp. 517–524. [Google Scholar] [CrossRef]
  11. Chang, E.; Dover, W.D. Stress concentration factor parametric equations for tubular X and DT joints. Int. J. Fatigue 1996, 18, 363–387. [Google Scholar] [CrossRef]
  12. Zavvar, E.; Guedes Soares, C. Effects of fibre reinforced polymer on stress concentration factors in uniplanar DKT-joints subjected to the compression loading. In Trends in Maritime Technology and Engineering; Guedes Soares, C., Santos, T.A., Eds.; Taylor and Francis: London, UK, 2022. [Google Scholar]
  13. Zavvar, E.; Guedes Soares, C. Stress distribution in uniplanar KT joints reinforced with fibre reinforced polymer subjected to the axial loadings. In Advances in the Analysis and Design of Marine Structures; CRC Press: Boca Raton, FL, USA, 2023; pp. 565–575. [Google Scholar]
  14. Ahmadi, H.; Lotfollahi-Yaghin, M.A. Geometrically parametric study of central brace SCFs in offshore three-planar tubular KT-joints. J. Constr. Steel Res. 2012, 71, 149–161. [Google Scholar] [CrossRef]
  15. Ahmadi, H.; Ali Lotfollahi-Yaghin, M.; Yong-Bo, S.; Aminfar, M.H. Parametric study and formulation of outer-brace geometric stress concentration factors in internally ring-stiffened tubular KT-joints of offshore structures. Appl. Ocean Res. 2012, 38, 74–91. [Google Scholar] [CrossRef]
  16. Chiew, S.P.; Soh, C.K.; Wu, N.W. Experimental and numerical stress analyses of tubular XT-joint. J. Struct. Eng. 1999, 125, 1239–1248. [Google Scholar] [CrossRef]
  17. Chiew, S.P.; Soh, C.K.; Wu, N.W. Experimental and numerical SCF studies of multiplanar tubular XX-joint. J. Struct. Eng. 2000, 126, 1331–1338. [Google Scholar] [CrossRef]
  18. Nassiraei, H.; Mojtahedi, A.; Lotfollahi-Yaghin, M.A. Static strength of X-joints reinforced with collar plates subjected to brace tensile loading. Ocean Eng. 2018, 161, 227–241. [Google Scholar] [CrossRef]
  19. Bao, S.; Wang, W.; Li, X.; Zhou, J. Stress concentration factor distribution formulas for three-planar tubular Y-joints under axial force. Ocean Eng. 2022, 265, 112687. [Google Scholar] [CrossRef]
  20. Rahmanli, U.; Becque, J. Stress concentration factors for multi-planar DKT tubular truss joints under axially balanced loads. J. Constr. Steel Res. 2021, 184, 106781. [Google Scholar] [CrossRef]
  21. Ahmadi, H.; Zavvar, E. The effect of multi-planarity on the SCFs in offshore tubular KT-joints subjected to in-plane and out-of-plane bending loads. Thin-Walled Struct. 2016, 106, 148–165. [Google Scholar] [CrossRef]
  22. Zavvar, E.; Hectors, K.; De Waele, W. Stress concentration factors of multi-planar tubular KT-joints subjected to in-plane bending moments. Mar. Struct. 2021, 78, 103000. [Google Scholar] [CrossRef]
  23. Kuang, J.G.; Potvin, A.B.; Leick, R.D. Stress concentration in tubular joints. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 5–8 May 1975; p. OTC-2205. [Google Scholar]
  24. Ahmadi, H.; Lotfollahi-Yaghin, M.A. Stress concentration due to in-plane bending (IPB) loads in ring-stiffened tubular KT-joints of offshore structures: Parametric study and design formulation. Appl. Ocean Res. 2015, 51, 54–66. [Google Scholar] [CrossRef]
  25. Nassiraei, H.; Rezadoost, P. Stress concentration factors in tubular T-joints reinforced with external ring under in-plane bending moment. Ocean Eng. 2022, 266, 112551. [Google Scholar] [CrossRef]
  26. Ahmadi, H.; Lotfollahi-Yaghin, M.A.; Yong-Bo, S. Chord-side SCF distribution of central brace in internally ring-stiffened tubular KT-joints: A geometrically parametric study. Thin-Walled Struct. 2013, 70, 93–105. [Google Scholar] [CrossRef]
  27. Ahmadi, H.; Alinezhad, R.; Alizadeh Atalo, A. Stress concentration analysis of internally ring-stiffened two-planar tubular KK-joints. Ships Offshore Struct. 2022, 17, 2203–2217. [Google Scholar] [CrossRef]
  28. Nassiraei, H. Geometrical effects on the LJF of tubular T/Y-joints with doubler plate in offshore wind turbines. Ships Offshore Struct. 2022, 17, 481–491. [Google Scholar] [CrossRef]
  29. Ahmadi, H.; Yeganeh, A.; Mohammadi, A.H.; Zavvar, E. Probabilistic analysis of stress concentration factors in tubular KT-joints reinforced with internal ring stiffeners under in-plane bending loads. Thin-Walled Struct. 2016, 99, 58–75. [Google Scholar] [CrossRef]
  30. Nassiraei, H.; Rezadoost, P. Local joint flexibility of tubular X-joints stiffened with external ring or external plates. Mar. Struct. 2021, 80, 103085. [Google Scholar] [CrossRef]
  31. Nassiraei, H.; Yara, A. Local joint flexibility of tubular K-joints reinforced with external plates under IPB loads. Mar. Struct. 2022, 84, 103199. [Google Scholar] [CrossRef]
  32. Nassiraei, H.; Yara, A. Numerical Analysis of Local Joint Flexibility of K-joints with External Plates under Axial Loads in Offshore Tubular Structures. J. Mar. Sci. Appl. 2022, 21, 134–144. [Google Scholar] [CrossRef]
  33. Nassiraei, H.; Rezadoost, P. Stress concentration factors in tubular T/Y-joints strengthened with FRP subjected to compressive load in offshore structures. Int. J. Fatigue 2020, 140, 105719. [Google Scholar] [CrossRef]
  34. Nassiraei, H.; Rezadoost, P. Parametric study and formula for SCFs of FRP-strengthened CHS T/Y-joints under out-of-plane bending load. Ocean Eng. 2021, 221, 108313. [Google Scholar] [CrossRef]
  35. Mohamed, H.S.; Shao, Y.; Chen, C.; Shi, M. Static strength of CFRP-strengthened tubular TT-joints containing initial local corrosion defect. Ocean Eng. 2021, 236, 109484. [Google Scholar] [CrossRef]
  36. Ramezani, F.; Nunes, P.D.P.; Carbas, R.J.C.; Marques, E.A.S.; da Silva, L.F.M. The joint strength of hybrid composite joints reinforced with different laminates materials. J. Adv. Join. Process. 2022, 5, 100103. [Google Scholar] [CrossRef]
  37. Nassiraei, H. Local joint flexibility of CHS X-joints reinforced with collar plates in jacket structures subjected to axial load. Appl. Ocean Res. 2019, 93, 101961. [Google Scholar] [CrossRef]
  38. Nassiraei, H.; Rezadoost, P. Probabilistic analysis of the SCFs in tubular T/Y-joints reinforced with FRP under axial, in-plane bending, and out-of-plane bending loads. Structures 2022, 35, 1078–1097. [Google Scholar] [CrossRef]
  39. Ramezani, F.; Simões, B.D.; Carbas, R.J.C.; Marques, E.A.S.; da Silva, L.F.M. Developments in Laminate Modification of Adhesively Bonded Composite Joints. Materials 2023, 16, 568. [Google Scholar] [CrossRef]
  40. Soh, A.K.; Soh, C.K. Stress analysis of axially loaded T tubular joints reinforced with doubler plates. Comput. Struct. 1995, 55, 141–149. [Google Scholar] [CrossRef]
  41. Nassiraei, H. Probability distribution models for the ultimate strength of tubular T/Y-joints reinforced with collar plates at room and different fire conditions. Ocean Eng. 2023, 270, 113557. [Google Scholar] [CrossRef]
  42. Nassiraei, H. Probability distribution functions for the ultimate strength of X-joints with collar plates in compressive load at room and fire conditions. Structures 2024, 59, 105703. [Google Scholar] [CrossRef]
  43. Nassiraei, H. Probabilistic Analysis of Strength in Retrofitted X-Joints under Tensile Loading and Fire Conditions. Buildings 2024, 14, 2105. [Google Scholar] [CrossRef]
  44. Nassiraei, H. Static strength of tubular T/Y-joints reinforced with collar plates at fire induced elevated temperature. Mar. Struct. 2019, 67, 102635. [Google Scholar] [CrossRef]
  45. Nassiraei, H.; Chavoshi, H.R. Collar plate effect on local joint flexibility of tubular K-joints under IPB moments: Parametric study and predictive equations. Mar. Struct. 2024, 98, 103676. [Google Scholar] [CrossRef]
  46. Nassiraei, H.; Zhu, L.; Gu, C. Static capacity of collar plate reinforced tubular X-connections subjected to compressive loading: Study of geometrical effects and parametric formulation. Ships Offshore Struct. 2021, 16, 54–69. [Google Scholar] [CrossRef]
  47. Nassiraei, H. Local joint flexibility of CHS T/Y-connections strengthened with collar plate under in-plane bending load: Parametric study of geometrical effects and design formulation. Ocean Eng. 2020, 202, 107054. [Google Scholar] [CrossRef]
  48. Sundaravadivelu, R.; Nandakumar, C.G.; Srivastava, S.K.; Ganapathy, C. Experimental studies to determine strain concentration factors for space tubular joints. J. Strain Anal. Eng. Des. 1987, 22, 237–245. [Google Scholar] [CrossRef]
  49. Morahan, D.J.; van Foeken, R.J.; Lalani, M.; Wardenier, J. New data on the fatigue and ultimate limit state of grouted tubular joints. In Proceedings of the 1996 Offshore Technology Conference, Houston, TX, USA, 6–9 May 1996; MSL Engineering Ltd.: Cork, Ireland, 1996. Part 3 (of 4). Volume 2, pp. 467–471. [Google Scholar]
  50. Zavvar, E.; Giannini, G.; Taveira-Pinto, F.; Rosa Santos, P. Parametric study of stress concentration factors in KT connections reinforced with concrete subjected to axial loads. Ocean Eng. 2024, 298, 117209. [Google Scholar] [CrossRef]
  51. Zavvar, E.; Sadat Hosseini, A.; Lotfollahi-Yaghin, M.A. Stress concentration factors in steel tubular KT-connections with FRP-Wrapping under bending moments. Structures 2021, 33, 4743–4765. [Google Scholar] [CrossRef]
  52. Sadat Hosseini, A.; Zavvar, E.; Ahmadi, H. Stress concentration factors in FRP-strengthened steel tubular KT-joints. Appl. Ocean Res. 2021, 108, 102525. [Google Scholar] [CrossRef]
  53. Zavvar, E.; Henneberg, J.; Guedes Soares, C. Stress concentration factors in FRP-reinforced tubular DKT joints under axial loads. Mar. Struct. 2023, 90, 103429. [Google Scholar] [CrossRef]
  54. Zhao, X.-L.; Fernando, D.; Al-Mahaidi, R. CFRP strengthened RHS subjected to transverse end bearing force. Eng. Struct. 2006, 28, 1555–1565. [Google Scholar] [CrossRef]
  55. UK Health and Safety Executive. Stress Concentration Factors for Tubular Complex Joints; OTH 91 353; Lloyd’s Register of Shipping: London, UK, 1992. [Google Scholar]
  56. Sadat Hosseini, A.; Bahaari, M.R.; Lesani, M. Experimental and parametric studies of SCFs in FRP strengthened tubular T-joints under axially loaded brace. Eng. Struct. 2020, 213, 110548. [Google Scholar] [CrossRef]
  57. Tong, L.; Xu, G.; Zhao, X.L.; Zhou, H.; Xu, F. Experimental and theoretical studies on reducing hot spot stress on CHS gap K-joints with CFRP strengthening. Eng. Struct. 2019, 201, 109827. [Google Scholar] [CrossRef]
  58. Chang, E.; Dover, W.D. Prediction of stress distributions along the intersection of tubular Y and T-joints. Int. J. Fatigue 1999, 21, 361–381. [Google Scholar] [CrossRef]
  59. Zhao, X.-L.; Packer, J.A. Recommended Fatigue Design Procedure for Welded Hollow Section Joints; International Institute of Welding: Paris, France, 1999. [Google Scholar]
  60. Zavvar, E.; Sousa, F.; Giannini, G.; Taveira-Pinto, F.; Santos, P.R. Probability of maximum values of stress concentration factors in tubular DKT-joints reinforced with FRP under axial loads. Structures 2024, 66, 106809. [Google Scholar] [CrossRef]
  61. AWS D1.1/D1.1M:2015; Structural Welding Code. American National Standard: New York, NY, USA, 2015.
  62. Ahmadi, H.; Zavvar, E. Degree of bending (DoB) in offshore tubular KT-joints under the axial, in-plane bending (IPB), and out-of-plane bending (OPB) loads. Appl. Ocean Res. 2020, 95, 102015. [Google Scholar] [CrossRef]
  63. Sadat Hosseini, A.; Bahaari, M.R.; Lesani, M. Parametric study of FRP strengthening on stress concentration factors in an offshore tubular T-joint subjected to in-plane and out-of-plane bending moments. Int. J. Steel Struct. 2019, 19, 1755–1766. [Google Scholar] [CrossRef]
  64. Ahmadi, H.; Zavvar, E. Stress concentration factors induced by out-of-plane bending loads in ring-stiffened tubular KT-joints of jacket structures. Thin-Walled Struct. 2015, 91, 82–95. [Google Scholar] [CrossRef]
  65. Jolliffe, I.T.; Cadima, J. Principal component analysis: A review and recent developments. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2016, 374, 20150202. [Google Scholar] [CrossRef]
  66. Johnson, D.E. Applied Multivariate Methods for Data Analysts; Duxbury Press: Pacific Grove, CA, USA, 1998. [Google Scholar]
  67. Jobson, J.D. Applied Multivariate Data Analysis: Volume II: Categorical and Multivariate Methods; Springer Science & Business Media: New York, NY, USA, 2012. [Google Scholar]
  68. Gordon, A.D. A review of hierarchical classification. J. R. Stat. Soc. Ser. A (Gen.) 1987, 150, 119–137. [Google Scholar] [CrossRef]
  69. Gordon, A.D. Classification; CRC Press: Boca Raton, FL, USA, 1999. [Google Scholar]
  70. Sturges, H.A. The choice of a class interval. J. Am. Stat. Assoc. 1926, 21, 65–66. [Google Scholar] [CrossRef]
  71. Freedman, D.; Diaconis, P. On the histogram as a density estimator: L 2 theory. Z. Wahrscheinlichkeitstheor. Verwandte Geb. 1981, 57, 453–476. [Google Scholar] [CrossRef]
  72. Mann, H.B.; Wald, A. On the choice of the number of class intervals in the application of the chi square test. Ann. Math. Stat. 1942, 13, 306–317. [Google Scholar] [CrossRef]
  73. Williams, C.A., Jr. The choice of the number and width of classes for the chi-square test of goodness of fit. J. Am. Stat. Assoc. 1950, 45, 77–86. [Google Scholar]
  74. Stuart, A.; Ord, K. Kendall’s Advanced Theory of Statistics, Distribution Theory; John Wiley & Sons: Chichester, UK, 2010; Volume 1. [Google Scholar]
Figure 1. Uniplanar tubular DKT-joints in a fixed platform with related non-dimensional parameters.
Figure 1. Uniplanar tubular DKT-joints in a fixed platform with related non-dimensional parameters.
Jmse 12 01451 g001
Figure 2. Meshing the uniplanar KT-joints with the help of the sub-zone approach.
Figure 2. Meshing the uniplanar KT-joints with the help of the sub-zone approach.
Jmse 12 01451 g002
Figure 3. Dihedral angle Ψ (American Welding Society (AWS) [61]).
Figure 3. Dihedral angle Ψ (American Welding Society (AWS) [61]).
Jmse 12 01451 g003
Figure 4. The FRP length for a DKT-joint.
Figure 4. The FRP length for a DKT-joint.
Jmse 12 01451 g004
Figure 5. The proposed extrapolation points via IIW [59].
Figure 5. The proposed extrapolation points via IIW [59].
Jmse 12 01451 g005
Figure 6. SCF box plots for the five samples.
Figure 6. SCF box plots for the five samples.
Jmse 12 01451 g006
Figure 7. Density histograms of the SCF sample.
Figure 7. Density histograms of the SCF sample.
Jmse 12 01451 g007
Figure 8. P-P plots (left) and scatter diagrams (right).
Figure 8. P-P plots (left) and scatter diagrams (right).
Jmse 12 01451 g008aJmse 12 01451 g008b
Figure 9. Hierarchical cluster grouping types, cutting from 5.
Figure 9. Hierarchical cluster grouping types, cutting from 5.
Jmse 12 01451 g009
Table 1. The non-dimensional parameters.
Table 1. The non-dimensional parameters.
ParameterDefinitionValue(s)
βd/D0.4, 0.5, 0.6, 0.7
γD/2T12, 18, 24
τt/T0.4, 0.7, 1.0
θ-30°, 45°, 60°
NNumber of FRP sheets4, 8, 12, 16
ξEFRP/ESteel0.14–0.87
ηTFRP/Tchord0.3–1.20
βd/D0.4, 0.5, 0.6, 0.7
γD/2T12, 18, 24
Table 2. Loading conditions.
Table 2. Loading conditions.
No LoadingBrace
ABCDEF
1 111111
2 −11−1111
3 −1−11111
4 111−1−1−1
1 = Compression, −1 = Tension
Table 3. Materials of FRP.
Table 3. Materials of FRP.
NameE1 (MPa)υ32G23 (MPa)
Glass/vinyl ester28,0000.382540
Kevlar 9/Epoxy76,0000.372010
Carbon T300/5208181,0000.593240
Table 4. SCF descriptive statistics for samples.
Table 4. SCF descriptive statistics for samples.
12345
Mean20.5113.5511.8510.4520.56
Median17.9212.1411.0810.1918.02
Mode13.429.8912.245.399.05
Skewness0.890.850.970.410.89
Kurtosis 3.233.153.952.633.25
Standard deviation10.506.625.624.3110.47
Variance110.1243.8231.6018.52109.51
Max56.7635.0134.4925.4556.76
Min5.124.033.682.985.12
Count12961296129612961296
1st quartile12.678.647.896.4912.73
3rd quartile26.9417.3114.7513.7226.94
Table 5. Values of FD.
Table 5. Values of FD.
Sample12345
FD (nc)2019241720
Significance level0.050.050.050.050.05
Table 6. Pearson’s correlation coefficient of sample 5.
Table 6. Pearson’s correlation coefficient of sample 5.
Confidence Intervals
Pearson CorrelationSig. (2-Tailed)95% Confidence Intervals (2-Tailed)
LowerUpper
β-SCF−0.1110.058−0.0020.107
τ-SCF0.6930.0000.5110.587
θ-SCF0.1340.0000.0820.189
η-SCF0.0820.000−0.152−0.044
γ-η0.5060.0000.4650.546
ξ-SCF−0.1610.082−0.1030.006
γ-SCF0.5590.000−0.153−0.046
Table 7. Pearson’s correlation matrix of sample 5.
Table 7. Pearson’s correlation matrix of sample 5.
βτθηξγSCF
β10.0000.0000.0000.0000.000−0.111
τ0.00010.0000.0000.0000.0000.693
θ0.0000.00010.0000.0000.0000.134
η0.0000.0000.00010.000.5060.082
ξ0.0000.0000.0000.00010.000−0.161
γ0.0000.0000.0000.5060.00010.559
SCF−0.1110.6930.1340.082−0.1610.5591
Table 8. Model summary.
Table 8. Model summary.
SamplesRR2Adjusted R2Std. Error of the Estimate
10.9500.9030.9023.27270
20.9380.8800.8792.30080
30.9130.8340.8332.29590
40.8820.7780.7772.03349
50.9500.9030.9023.27276
Table 9. Coefficients.
Table 9. Coefficients.
ModelUnstandardized CoefficientsStandardized CoefficientstSig.95% Confidence Interval for B
BStd. ErrorBetaLower BoundUpper Bound
Sample 1(Constant)−18.6210.719 −25.893<0.001−20.032−17.210
β−10.3530.813−0.111−12.733<0.001−11.948−8.758
τ29.5910.3710.69379.730<0.00128.86230.319
θ6.5720.4250.13415.454<0.0015.7387.407
η −14.6110.545−0.270−26.816<0.001−15.679−13.542
ξ−5.4530.295−0.161−18.513<0.001−6.031−4.875
γ1.4870.0220.69669.078<0.0011.4441.529
Sample 2(Constant)−4.5480.506 −8.996<0.001−5.540−3.556
β−10.0850.572−0.170−17.642<0.001−11.206−8.963
τ18.2170.2610.67469.818<0.00117.70518.729
θ−1.1450.299−0.037−3.830<0.001−1.732−0.559
η −7.8540.383−0.230−20.504<0.001−8.605−7.102
ξ−4.1820.207−0.195−20.197<0.001−4.588−3.776
γ0.9190.0150.68060.742<0.0010.8890.949
Sample 3(Constant)−4.4310.504 −8.783<0.001−5.421−3.441
β−8.8590.570−0.176−15.530<0.001−9.978−7.740
τ15.5960.2600.68059.901<0.00115.08516.107
θ2.1840.2980.0837.319<0.0011.5982.769
η −8.0000.382−0.275−20.931<0.001−8.750−7.250
ξ−3.5190.207−0.193−17.030<0.001−3.924−3.114
γ0.7230.0150.63047.888<0.0010.6930.753
Sample 4(Constant)−2.6800.447 −5.999<0.001−3.557−1.804
β1.2550.5050.0332.4840.0130.2642.246
τ13.9950.2310.79760.687<0.00113.54214.447
θ.0910.2640.0050.3440.731−0.4270.609
η −5.8020.339−0.261−17.139<0.001−6.466−5.138
ξ−2.6090.183−0.187−14.257<0.001−2.968−2.250
γ0.3250.0130.37024.327<0.0010.2990.352
Sample 5(Constant)−18.6200.719 −25.892<0.001−20.031−17.210
β−10.3530.813−0.111−12.732<0.001−11.948−8.758
τ29.5910.3710.69379.730<0.00128.86330.319
θ6.5720.4250.13415.453<0.0015.7387.406
η −14.6100.545−0.270−26.815<0.001−15.679−13.541
ξ−5.4520.295−0.161−18.512<0.001−6.030−4.875
γ1.4860.0220.69669.075<0.0011.4441.529
Table 10. ANOVA.
Table 10. ANOVA.
ModelSSDf (Degree of Freedom)MSFSignificance
Sample 1Regression128,119.296621,353.2161993.656<0.001
Residual13,805.937128910.711
Total141,925.2331295
Sample 2Regression49,969.29568328.2161573.234<0.001
Residual6823.56812895.294
Total56,792.8631295
Sample 3Regression34,162.05165693.6751080.158<0.001
Residual6794.51512895.271
Total40,956.5661295
Sample 4Regression18,672.31063112.052752.596<0.001
Residual5330.12712894.135
Total24,002.4371295
Sample 5Regression128,118.641621,353.1071993.576<0.001
Residual13,806.424128910.711
Total141,925.0651295
Table 11. Proximity matrix.
Table 11. Proximity matrix.
Agglomeration Schedule
StageCluster CombinedCoefficientsStage Cluster First AppearsNext Stage
Cluster 1Cluster 2Cluster 1Cluster 2
1115311560.100145
29019040.100615
31601630.10073
49759780.100429
58938960.100716
68308330.100231
1290128.443125512851292
129114268.460128412871293
129211310.484129012881294
1293141512.317129112891294
129411417.291129212931295
129512729.821129412860
Table 12. Centroids of k-means cluster analysis with MATLAB, sample 5.
Table 12. Centroids of k-means cluster analysis with MATLAB, sample 5.
Cluster Numberθτγβ η ξSCF
00.52360.4120.40.1920.36717.9351
11.04720.7240.70.5760.874412.456
20.78541180.50.0960.135316.2223
31.04720.7120.60.3840.874429.9557
41.04720.7120.70.7680.874452.8629
51.04720.7120.70.2880.874439.703
61.04720.7120.70.4320.874446.0828
71.04720.7120.70.1440.874425.1819
81.04720.7120.70.5760.874420.6465
91.04720.7120.70.5760.874434.5585
Table 13. Initial Cluster Centers.
Table 13. Initial Cluster Centers.
12345678910111213141516
β0.500000.400000.500000.700000.500000.400000.500000.400000.600000.700000.500000.400000.500000.700000.400000.60000
τ1.000001.000001.000000.400001.000001.000001.000000.700001.000000.400001.000000.700000.400000.700001.000001.00000
θ1.047201.047200.785400.523601.047200.785401.047200.785400.785400.523601.047200.523600.523600.523600.785400.78540
η 0.09600.76800.28800.57600.19200.19200.76800.38400.38400.38400.14400.57600.57600.57600.57600.3840
ξ0.135270.874400.874400.367150.135270.367150.874400.874400.367150.874400.135270.874400.874400.367150.874400.87440
γ12.0024.0018.0018.0024.0024.0024.0012.0012.0012.0018.0024.0024.0018.0018.0024.00
SCF24.1426.2230.928.3156.7648.8335.1011.5517.915.1241.2718.6511.7415.1221.7642.51
Table 14. Number of Cases in each Cluster.
Table 14. Number of Cases in each Cluster.
Cluster180
282
382
496
514
636
774
888
9117
10147
1154
1284
1391
1486
15114
1651
Valid1296
Missing0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zavvar, E.; Sousa, F.; Taveira-Pinto, F.; Rosa Santos, P. Multivariate Data Analysis of Maximum Stress Concentration Factors in FRP-Retrofitted Two-Planar KT-Joints under Axial Loads for Offshore Renewables. J. Mar. Sci. Eng. 2024, 12, 1451. https://doi.org/10.3390/jmse12081451

AMA Style

Zavvar E, Sousa F, Taveira-Pinto F, Rosa Santos P. Multivariate Data Analysis of Maximum Stress Concentration Factors in FRP-Retrofitted Two-Planar KT-Joints under Axial Loads for Offshore Renewables. Journal of Marine Science and Engineering. 2024; 12(8):1451. https://doi.org/10.3390/jmse12081451

Chicago/Turabian Style

Zavvar, Esmaeil, Fernanda Sousa, Francisco Taveira-Pinto, and Paulo Rosa Santos. 2024. "Multivariate Data Analysis of Maximum Stress Concentration Factors in FRP-Retrofitted Two-Planar KT-Joints under Axial Loads for Offshore Renewables" Journal of Marine Science and Engineering 12, no. 8: 1451. https://doi.org/10.3390/jmse12081451

APA Style

Zavvar, E., Sousa, F., Taveira-Pinto, F., & Rosa Santos, P. (2024). Multivariate Data Analysis of Maximum Stress Concentration Factors in FRP-Retrofitted Two-Planar KT-Joints under Axial Loads for Offshore Renewables. Journal of Marine Science and Engineering, 12(8), 1451. https://doi.org/10.3390/jmse12081451

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop