Fatigue life of tubular joints in offshore structures is significantly influenced by the degree of bending (DoB). The DoB exhibits considerable scatter calling for greater emphasis in accurate determination of its governing probability distribution which is a key input for the fatigue reliability analysis of a tubular joint. Although the tubular X-joints are commonly found in offshore jacket structures, as far as the authors are aware, no comprehensive research has been carried out on the probability distribution of the DoB in tubular X-joints. In the present paper, results of parametric equations available for the calculation of the DoB have been used to develop probability distribution models for the DoB in the chord member of tubular X-joints subjected to four types of bending loads. Based on a parametric study, a set of samples was prepared and density histograms were generated for these samples using Freedman-Diaconis method. Twelve different probability density functions (PDFs) were fitted to these histograms. In each case, Kolmogorov-Smirnov test was used to evaluate the goodness of fit. Finally, after substituting the values of estimated parameters for each distribution, a set of fully defined PDFs have been proposed for the DoB in tubular X-joints subjected to bending loads.
1. Introduction
Offshore jacket-type platforms are mainly fabricated with circular hollow section (CHS) members. The intersection between CHS members is called a tubular joint. Figure 1 shows a tubular X-joint along with the three commonly named positions along the brace/chord intersection: saddle, crown toe, and crown heel. Nondimensional geometrical parameters including α, β, γ, τ, and αB which are used to feasibly relate the behavior of a tubular joint to its geometrical characteristics are defined in Figure 1.
Geometrical notation for a tubular X-joint.
Tubular joints are subjected to wave induced cyclic loads and thus are susceptible to fatigue damage. The stress-life (S-N) approach, based on the hot-spot stress (HSS), is widely used to estimate the fatigue life of the joint. The HSS can be calculated through the multiplication of nominal stress by the stress concentration factor (SCF). However, the investigation of a large number of fatigue test results has shown that tubular joints with different geometry or loading type but with similar HSSs often can endure significantly different numbers of cycles before failure [1]. These differences are thought to be attributable to changes in crack growth rate which is dependent on the through-the-thickness stress distribution as well as the HSS. The stress distribution across the wall thickness is assumed to be a linear combination of membrane and bending stresses. It can be characterized by the degree of bending (DoB), that is, the ratio of bending stress to total stress.
As mentioned before, it has become evident that the HSS is not enough to characterize all aspects of fatigue failure. Therefore, the standard stress-life approach may be unconservative for the joints with low DoB. Hence, the current standard HSS-based S-N approach can be modified to include the effect of the DoB representing the through-the-thickness stress distribution in the tubular joint in order to reduce the scatter in the S-N curve and to obtain more accurate fatigue life prediction. The other shortcoming of the S-N approach is that this method gives only the total life and cannot be used to predict fatigue crack growth and the remaining life of cracked joints. For the fatigue analysis of cracked joints, fracture mechanics (FM) should be used. The accurate determination of a stress intensity factor (SIF) is the key for FM calculations. It is well known that it is necessary to take the complex stress field in tubular welded joints into account to have accurate SIF data. Owing to the complexities introduced by the structural geometry and the nature of the local stress fields, it is impossible to calculate the SIFs analytically. This problem is often tackled by using the simplified models, such as the flat plate solution or T-Butt weight function based method, with an appropriate load shedding model. In order to use these simplified SIF models to calculate the remaining fatigue life of tubular joints, the information is required again on the distribution of through-the-thickness stress acting in the anticipated crack path, which can be characterized by the DoB. Thus, the DoB is an important input parameter for the calculation of fatigue crack growth in tubular welded joints.
Deterministic fatigue analyses typically produce conservative results, since limiting assumptions are to be made on key input parameters. Some of the key parameters of the problem can exhibit stochastic behavior. This highlights the necessity of conducting a reliability analysis in which these key parameters can be modeled as random quantities. The fundamentals of reliability assessment, if properly applied, can provide immense insight into the performance and safety of the structural system. Under any specific loading condition, the DoB value in a tubular joint is mainly determined by the joint geometry and exhibits considerable scatter calling for greater emphasis in accurate determination of its governing probability distribution which is an essential input for the fatigue reliability analysis of a tubular joint. As far as the authors are aware, despite the considerable research work accomplished on the study of SCFs and SIFs in tubular joints and a few projects defined about the deterministic analysis of the DoB (see the next paragraph), no comprehensive research has been carried out on the probability distribution of the DoB in tubular joints. What has been used so far as the probability distribution of the DoB in the fatigue reliability analysis of offshore structures is mainly based on assumptions and limited observations, especially in terms of distribution parameters.
Bowness and Lee [2] investigated the fatigue crack curvature under the weld toe in an offshore tubular joint. Lee et al. [3] numerically studied the cracked tubular T-, Y-, and K-joints under combined loads. Shao [4] analyzed the stress intensity factor (SIF) for cracked tubular K-joints subjected to balanced axial load. Wordsworth and Smedley [5] studied stress concentrations at unstiffened tubular joints. Efthymiou [6] developed a set of SCF formulae and generalized influence functions for use in fatigue analysis. Chang and Dover [7] proposed parametric equations to predict stress distributions along the intersection of tubular X- and DT-joints. Lotfollahi-Yaghin and Ahmadi [8] investigated geometric stress distribution along the weld toe of the outer brace in two-planar tubular DKT-joints. Ahmadi and Lotfollahi-Yaghin [9] performed a geometrically parametric study on central brace SCFs in offshore three-planar tubular KT-joints. Ahmadi et al. [10] studied chord-side SCF distribution of central brace in internally ring-stiffened tubular KT-joints. A series of systematic thin shell FE analyses were carried out by Chang and Dover [11] for 330 tubular X- and DT-joints typical of those found in offshore structures, under six different types of loading. Mean and design equations for DoBs at critical positions in axially loaded tubular K-joints were derived by Morgan and Lee [12] from a previously established FE database of 254 joints. Design equations met all the acceptance criteria recommended by the UK DoE [13]. Lee and Bowness [14] proposed an engineering methodology for estimating SIF solutions for semielliptical weld-toe cracks in tubular joints. The SIFs for a grouted tubular joint were determined both numerically and empirically by Shen and Choo [15].
In the present paper, results of parametric equations available for the calculation of the DoB have been used to propose probability distribution models for the DoB in the chord member of tubular X-joints subjected to four different types of bending loads including single and double in-plane bending (IPB) and out-of-plane bending (OPB) loadings (Figure 2). Based on a parametric study, a set of samples was prepared and density histograms were generated for these samples using Freedman-Diaconis method. Twelve different probability density functions (PDFs) were fitted to these histograms. The maximum likelihood (ML) method was used to determine the parameters of fitted distributions. In each case, Kolmogorov-Smirnov test was used to evaluate the goodness of fit. Finally, after substituting the values of estimated parameters in distribution models, a set of fully defined PDFs have been proposed for the DoB in tubular X-joints under bending loads.
Considered loading conditions.
2. DoB in Tubular X-Joints Subjected to Bending Loads
As mentioned earlier, the degree of bending (DoB) is the ratio of bending stress over total stress expressed as(1)DoB=σBσT=σBσB+σM,where σB and σM are the bending and membrane stress components and σT is the total stress on the outer surface of the chord (Figure 3).
Chang and Dover [11] proposed a set of equations for the calculation of the DoB of the chord member in tubular X-joints subjected to bending loads (2). In these equations, the DoB is corresponding to the position of the HSS, that is, the position in which the maximum SCF occurs. DoBS-IPB, DoBS-OPB, DoBD-IPB, and DoBD-OPB denote the DoB under the single IPB, single OPB, double IPB, and double OPB loadings, respectively (Figure 2), and DoB+ and DoB− stand for the DoB at the positive and negative HSS positions. It should be noted that single IPB and single OPB loadings can actually occur in a tubular joint of an offshore jacket structure depending on the wave incident angle, location of the joint, relative position of the wave crest, and design load combination. In such loading cases, the single IPB/OPB moment is balanced by the internal forces of the chord member instead of the other brace member(2)DoB-S-IPB=exp-0.7153+0.654sinθ-1.491β5γ2βsedxps+0.00129αθ+0.00832βγτ2-0.123βθsesxps+0.00316γθ-0.000004γ4-0.2sin2θsesxps+1.19β4+0.01261β2sesxps-0.0009γβ2+0.0001γ2β,DoB+S-IPB=0.7924-0.0661lnβθ+0.00963βγτ2-0.1192lnγθ-0.0428τ4-0.209β5-0.536lnsinθ-0.0126γ+0.0068β2γ-0.056lnβ-0.000004γ4+0.00071γ2-0.2lnθ+0.262sin2θ+0.00082α,DoBS-OPB=exp0.092972-0.1492sinθ-0.1593βγθsexsps+0.00269βγτ2+0.0392βθsesxps-0.0136τθ-0.00104γθ+0.00002γ2,DoB-D-IPB=exp-1.2643-0.0184γτ+1.67sinθγθsexsps+0.0149βγτ2+0.00264γθ-0.0655β2γτsexsps-0.885sin2θ+0.0057β2+0.057βγτsesxps+0.091τθ+0.127θ+0.00024αγsexsps+0.0018αθ-0.00012α2,DoB+D-IPB=1.521+0.00063αγ+0.011392β2γτ-1.116sinθ-0.079βθ-0.0053γθ+0.00551β2+0.00011γ2-0.0129β3γ-0.00005γ2β+0.575sin2θ-0.083θ-0.00012α2+25γ3-0.0078α,DoBD-OPB=0.8388-0.0183lnβθ+0.00215βγτ2-0.2123β5+0.00475β2γ-0.00003γ2-0.0371lnβ-0.0433lnsinθ-0.0052lnτθ.
The validity ranges for the application of (2) are as follows: (3)6.0≤α≤40.0,0.2≤β≤0.8,7.6≤γ≤32.0,0.2≤τ≤1.0,35∘≤θ≤90∘.
3. Preparation of the DoB Samples
A MATLAB code was developed to generate six samples for the DoB based on (2). These equations have five variables including α, β, γ, τ, and θ. Developed MATLAB code divided the validity range for each parameter (3) into equal intervals and calculated the DoB for all possible combinations of the boundary values. For example, to generate the DoB-S-IPB sample, developed code calculated the DoB for all possible combinations of five values of α (6, 14.5, 23, 31.5, and 40), five values of β (0.2, 0.35, 0.5, 0.65, and 0.8), five values of γ (7.6, 13.7, 19.8, 25.9, and 32), five values of τ (0.2, 0.4, 0.6, 0.8, and 1.0), and five values of θ (35°, 48.75°, 62.5°, 76.25°, and 90°) which led to 3125 (55) data points for this sample.
Values of the size (n), mean (μ), standard deviation (σ), coefficient of skewness (α3), and coefficient of kurtosis (α4) for generated samples are listed in Table 1.
Values of statistical measures for the DoB samples.
Statistical measure
Sample
DoB-S-IPB
DoB+S-IPB
DoBS-OPB
DoB-D-IPB
DoB+D-IPB
DoBD-OPB
n
3125
3125
625
3125
243
625
μ
0.4220
−0.6714
0.9054
0.8171
0.5890
0.8983
σ
0.3050
1.65
0.0274
0.0747
0.2002
0.0337
α3
−0.0633
−0.9183
−0.1602
−0.2505
0.0183
−0.2918
α4
1.4590
2.3547
2.3787
2.9113
2.2840
3.2364
The value of α3 for DoB-S-IPB, DoB+S-IPB, DoBS-OPB, DoB-D-IPB, and DoBD-OPB samples is negative meaning that in these cases; the distribution is expected to have a longer tail on the left, which is toward decreasing values, than on the right. However, the DoB+D-IPB sample has a positive α3 value which means that its distribution is expected to have a longer tail on the right. Moreover, in DoB-S-IPB, DoB+S-IPB, DoBS-OPB, DoB-D-IPB, and DoB+D-IPB samples, the value of α4 is smaller than three which means that, in these cases, the probability distribution is expected to be mild-peak (platykurtic). On the contrary, in DoBD-OPB sample, the value of α4 is greater than three meaning that, in this case, a sharp-peak (leptokurtic) probability distribution is to be expected.
4. Generation of the Density Histograms
For generating a density histogram, the range (R) should be divided into a number of classes/cells/bins. The number of occurrences in each class is counted and tabulated. These are called frequencies. Then, the relative frequency of each class can be obtained through dividing its frequency by the sample size. Afterwards, the density is calculated for each class through dividing the relative frequency by the class width. The width of classes is usually made equal to facilitate interpretation.
Care should be exercised in the choice of the number of classes (nc). Too few will cause an omission of some important features of the data; too many will not give a clear overall picture because there may be high fluctuations in the frequencies. In the present research, Freedman-Diaconis rule was adapted to determine the number of classes:(4)nc=Rn1/32IQR,where R is the range of sample data, n is the sample size, and IQR is the interquartile range calculated as follows:(5)IQR=Q3-Q1,where Q1 is the lower quartile which is the median of the lower half of the data and, likewise, Q3 is the upper quartile that is the median of the upper half of the data.
For example, density histograms of DoBS-OPB, DoB-D-IPB, DoB+D-IPB, and DoBD-OPB samples are shown in Figure 4. As it was expected from values of α3 and α4 (Table 1), histograms of (a), (b), and (d) have a longer tail on the left than on the right, while the histogram of (c) has a longer tail on the right. It can also be seen that histograms of (a), (b), and (c) are platykurtic; while the histogram of (d) is leptokurtic.
Generated histograms for DoB samples: (a) DoBS-OPB, (b) DoB-D-IPB, (c) DoB+D-IPB, and (d) DoBD-OPB.
5. Application of Maximum Likelihood Method for PDF Fitting
In order to investigate the degree of fitting of various distributions to the sample data, twelve different PDFs were fitted to the generated histograms. For example, PDFs fitted to density histograms of DoBS-OPB, DoB-D-IPB, DoB+D-IPB, and DoBD-OPB samples are shown in Figure 5.
PDFs fitted to the generated histograms of DoB samples: (a) DoBS-OPB, (b) DoB-D-IPB, (c) DoB+D-IPB, and (d) DoBD-OPB.
In each case, distribution parameters were estimated using the maximum likelihood (ML) method. Results are given in Table 2. It should be noted that none of the considered distributions was acceptably fitted to the DoB-S-IPB and DoB+S-IPB samples. Hence, no data is provided for these two samples in Table 2.
Estimated parameters of PDFs fitted to the density histograms of DoB samples.
Fitted PDF
Parameters
Estimated values
DoB-S-IPB
DoB+S-IPB
DoBS-OPB
DoB-D-IPB
DoB+D-IPB
DoBD-OPB
Birnbaum-Saunders
β0γ0
0.90498 0.0303296
0.813539 0.093508
—
0.897649 0.037795
Extreme value
μσ
0.91884 0.0251115
0.853587 0.0693142
0.688765 0.190828
0.91473 0.0318889
Gamma
ab
1089.95 0.000830676
116.665 0.0070038
—
704.135 0.00127574
Generalized extreme value
kσμ
−0.342762 0.0280815 0.896484
−0.331497 0.0772041 0.79235
−0.286603 0.19848 0.518683
−0.325209 0.0350571 0.887031
Log-logistic
μs
−0.0988711 0.0177489
−0.201561 0.05286
−0.56315 0.219722
−0.106302 0.0209965
Logistic
μs
—
—
0.906074 0.0160357
0.819049 0.0427215
0.590133 0.118319
0.899405 0.0187924
Lognormal
μσ
−0.0998412 0.0303511
−0.206291 0.0934124
−0.596375 0.386297
−0.107973 0.0378179
Nakagami
μΩ
273.263 0.820492
—
2.2093 0.386864
177.064 0.80806
Normal (Gaussian)
μσ
0.905396 0.0273892
0.817095 0.074663
0.589015 0.200225
0.89829 0.0337215
Beta
ab
97.8619 10.2226
—
—
64.4534 7.29193
Rician
sσ
—
0.813649 0.0748107
0.546732 0.209701
—
Weibull
λk
—
—
0.657831 3.29088
—
The ML procedure is an alternative to the method of moments. As a means of finding an estimator, statisticians often give it preference. For a random variable X with a known PDF, fX(x), and observed values x1,x2,…,xn, in a random sample of size n, the likelihood function of θ, where θ represents the vector of unknown parameters, is defined as(6)Lθ=∏i=1nfXxi∣θ.
The objective is to maximize L(θ) for the given data set. This is easily done by taking m partial derivatives of L(θ), where m is the number of parameters, and equating them tozero. Then the maximum likelihood estimators (MLEs) of the parameter set θ are found by solving the equations. In this way, the greatest probability is given to the observedset of events, provided that the true form of the probability distribution is known.
6. Using Kolmogorov-Smirnov Test to Evaluate the Goodness of Fit
The Kolmogorov-Smirnov goodness-of-fit test is a nonparametric test based on the cumulative distribution function (CDF) of a continuous variable. It is not applicable to discrete variables. The test statistic, in a two-sided test, is the maximum absolute difference (i.e., usually the vertical distance) between the empirical and hypothetical CDFs. For a continuous variate X, let x1,x(1),…,x(n) represent the order statistics of a sample of the size n, that is, the values arranged in increasing order. The empirical or sample distribution function Fn(x) is a step function. This gives the proportion of values not exceeding x and is defined as(7)Fnx=0,Forx<x(1)kn,Forxk≤x<xk+1k=1,2,…,n-11,Forx≥xn.
For example, empirical CDFs for DoBS-OPB, DoB-D-IPB, DoB+D-IPB, and DoBD-OPB samples have been shown in Figure 6.
Empirical CDFs for generated DoB samples: (a) DoBS-OPB, (b) DoB-D-IPB, (c) DoB+D-IPB, and (d) DoBD-OPB.
Let F0(x) denote a completely specified theoretical continuous CDF. The null hypothesis H0 is that the true CDF of X is the same as F0(x). That is, under the null hypothesis,(8)limn→∞PrFn(x)=F0(x)=1.
The test criterion is the maximum absolute difference between Fn(x) and F0(x), formally defined as(9)Dn=supxFn(x)-F0(x).
Theoretical continuous CDFs fitted to the empirical distribution functions of DoBS-OPB, DoB-D-IPB, DoB+D-IPB, and DoBD-OPB samples have been shown in Figure 7.
Theoretical CDFs fitted to the empirical CDFs of DoB samples: (a) DoBS-OPB, (b) DoB-D-IPB, (c) DoB+D-IPB, and (d) DoBD-OPB.
A large value of this statistic (Dn) indicates a poor fit. Hence, the acceptable values must be known. The critical values Dn,ξ for large samples, say n>35, are (1.3581/n) and (1.6276/n) for ξ=0.05 and 0.01, respectively [16] where ξ is the level of significance in hypothesis testing.
Results of Kolmogorov-Smirnov test for DoBS-OPB, DoB-D-IPB, DoB+D-IPB, and DoBD-OPB samples are given in Tables 3−6, respectively. It should be noted that, according to the results of Kolmogorov-Smirnov test, none of considered continuous CDFs was acceptably fitted to the DoB-S-IPB and DoB+S-IPB samples. Hence, no table is provided here for these two samples.
Results of Kolmogorov-Smirnov goodness-of-fit test for DoBS-OPB sample.
Fitted distribution
Test statistic
Critical value
Test result
ξ=0.05
ξ=0.01
ξ=0.05
ξ=0.01
Birnbaum-Saunders
0.0613
0.0541
0.0648
Reject
Accept
Extreme value
0.0515
Accept
Accept
Gamma
0.0596
Reject
Accept
Generalized extreme value
0.0425
Accept
Accept
Log-logistic
0.0547
Reject
Reject
Logistic
0.0521
Accept
Accept
Lognormal
0.0612
Reject
Accept
Nakagami
0.0723
Reject
Accept
Normal (Gaussian)
0.0578
Reject
Reject
Beta
0.0404
Accept
Accept
Results of Kolmogorov-Smirnov goodness-of-fit test for DoB-D-IPB sample.
Fitted distribution
Test statistic
Critical value
Test result
ξ=0.05
ξ=0.01
ξ=0.05
ξ=0.01
Birnbaum-Saunders
0.0411
0.0242
0.0291
Reject
Reject
Extreme value
0.0435
Reject
Reject
Gamma
0.0344
Reject
Reject
Generalized extreme value
0.0185
Accept
Accept
Log-logistic
0.0278
Reject
Accept
Logistic
0.0235
Accept
Accept
Lognormal
0.0408
Reject
Reject
Normal (Gaussian)
0.0221
Accept
Accept
Rician
0.0221
Accept
Accept
Results of Kolmogorov-Smirnov goodness-of-fit test for DoB+D-IPB sample.
Fitted distribution
Test statistic
Critical value
Test result
ξ=0.05
ξ=0.01
ξ=0.05
ξ=0.01
Extreme value
0.0746
0.0864
0.1037
Accept
Accept
Generalized extreme value
0.0666
Accept
Accept
Log-logistic
0.0811
Accept
Accept
Logistic
0.0788
Accept
Accept
Lognormal
0.1062
Reject
Reject
Nakagami
0.0721
Accept
Accept
Normal (Gaussian)
0.0710
Accept
Accept
Rician
0.0671
Accept
Accept
Weibull
0.0660
Accept
Accept
Results of Kolmogorov-Smirnov goodness-of-fit test for DoBD-OPB sample.
Fitted distribution
Test statistic
Critical value
Test result
ξ=0.05
ξ=0.01
ξ=0.05
ξ=0.01
Birnbaum-Saunders
0.0650
0.0541
0.0648
Reject
Reject
Extreme value
0.0721
Reject
Reject
Gamma
0.0628
Reject
Accept
Generalized extreme value
0.0611
Reject
Accept
Log-logistic
0.0315
Accept
Accept
Logistic
0.0284
Accept
Accept
Lognormal
0.0651
Reject
Reject
Nakagami
0.0608
Reject
Accept
Normal (Gaussian)
0.0588
Reject
Accept
Beta
0.0587
Reject
Accept
It is evident in Tables 3−6 that Beta, Generalized Extreme Value, Weibull, and Logistic distributions have the smallest values of test statistic for DoBS-OPB, DoB-D-IPB, DoB+D-IPB, and DoBD-OPB samples, respectively.
7. Developed Probability Models for DoBs
Based on the results of Kolmogorov-Smirnov goodness-of-fit test (Tables 3−6), it can be concluded that Beta, Generalized Extreme Value, Weibull, and Logistic distributions are the best probability models for DoBS-OPB, DoB-D-IPB, DoB+D-IPB, and DoBD-OPB in tubular X-joints under bending loads, respectively. The PDFs of these distributions are given by the following equations:(10)fXx=Γa+bΓaΓbxa-11-xb-1Betadistribution,fX(x)=1σexp-1+kx-μσ-1/k1+kx-μσ-1-1/kGeneralizedExtremeValuedistribution,fXx=kλxλk-1exp-xλk;x≥0Weibulldistribution,fXx=exp-x-μ/ss1+exp-x-μ/s2Logisticdistribution,where Γ(a) is the Gamma function defined as follows:(11)Γa=∫0∞e-rra-1dr.
After substituting the values of estimated parameters from Table 2, the following probability density functions are proposed for the DoBS-OPB, DoB-D-IPB, DoB+D-IPB, and DoBD-OPB in tubular X-joints subjected to bending loads:(12)fXx=5.9328×1014x96.86191-x9.2226DoBS-OPB,fXx=12.9527exp4.2938x-4.40223.0166·4.4022-4.2938x2.01660DoB-D-IPB,fXx=5.0026x0.6578312.29088·exp-x0.6578313.29088;x≥0ssssssssssssssssssssssssss(DoB+D-IPB),fXx=exp-x-0.899405/0.01879240.01879241+exp-x-0.899405/0.01879242dddddddddddddddddddddddddddddd(DoBD-OPB).
Developed PDFs, shown in Figure 8, can be adapted in the fatigue reliability analysis of tubular X-joints commonly found in offshore jacket structures subjected to bending loads.
Proposed PDFs for the DoB values in tubular X-joints subjected to the bending loads: (a) DoBS-OPB—Beta distribution, (b) DoB-D-IPB—Generalized Extreme Value distribution, (c) DoB+D-IPB—Weibull distribution, and (d) DoBD-OPB—Logistic distribution.
8. Conclusions
In the present paper, results of parametric equations available for the calculation of the DoB were used to propose probability distribution models for the DoB in the chord member of tubular X-joints subjected to four different types of bending loads including single and double IPB and OPB loadings. Based on a parametric study, a set of samples was prepared and density histograms were generated for these samples using Freedman-Diaconis method. Twelve different PDFs were fitted to these histograms. The ML method was used to determine the parameters of fitted distributions. In each case, Kolmogorov-Smirnov test was used to evaluate the goodness of fit. Based on the results of this test, it was concluded that Beta, Generalized Extreme Value, Weibull, and Logistic distributions are the best probability models for the DoBS-OPB, DoB-D-IPB, DoB+D-IPB, and DoBD-OPB in tubular X-joints under bending loads, respectively. Finally, after substituting the values of estimated parameters in distribution models, a set of fully defined PDFs were proposed for the DoB in tubular X-joints subjected to the bending loads.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
ConnollyM. P. M.1986London, UKUniversity College LondonBownessD.LeeM. M. K.Fatigue crack curvature under the weld toe in an offshore tubular joint199820648149010.1016/s0142-1123(98)00017-62-s2.0-0032120185LeeC. K.TjhenL. S.ChiewS. P.YongboS.Numerical models verification of cracked tubular T, Y and K-joints under combined loads2005727983100910.1016/j.engfracmech.2004.07.0062-s2.0-12444252160ShaoY.-B.Analysis of stress intensity factor (SIF) for cracked tubular K-joints subjected to balanced axial load2006131446410.1016/j.engfailanal.2004.12.0312-s2.0-27244454342WordsworthA. C.SmedleyG. P.Stress concentrations at unstiffened tubular jointsProceedings of the European Offshore Steels Research Seminar, Paper 311978Cambridge, UKEfthymiouM.Development of SCF formulae and generalized influence functions for use in fatigue analysisProceedings of the Offshore Tubular Joints Conference (OTJ '88)1988Surrey, UKChangE.DoverW. D.Parametric equations to predict stress distributions along the intersection of tubular X and DT-joints19992166196352-s2.0-003266041610.1016/S0142-1123(99)00018-3Lotfollahi-YaghinM. A.AhmadiH.Geometric stress distribution along the weld toe of the outer brace in two-planar tubular DKT-joints: parametric study and deriving the SCF design equations201124323926010.1016/j.marstruc.2011.02.0052-s2.0-79958072834AhmadiH.Lotfollahi-YaghinM. A.Geometrically parametric study of central brace SCFs in offshore three-planar tubular KT-joints20127114916110.1016/j.jcsr.2011.10.0242-s2.0-84855997463AhmadiH.Lotfollahi-YaghinM. A.Yong-BoS.Chord-side SCF distribution of central brace in internally ring-stiffened tubular KT-joints: a geometrically parametric study2013709310510.1016/j.tws.2013.04.0112-s2.0-84878351821ChangE.DoverW. D.Prediction of degree of bending in tubular X and DT joints19992121471612-s2.0-0033078409MorganM. R.LeeM. M. K.Prediction of stress concentrations and degrees of bending in axially loaded tubular K-joints1998451679710.1016/S0143-974X(97)00059-X2-s2.0-0031599119UK Department of Energy (DoE)1995London, UKUK Department of Energy (DoE)LeeM. M. K.BownessD.Estimation of stress intensity factor solutions for weld toe cracks in offshore tubular joints200224886187510.1016/S0142-1123(01)00209-22-s2.0-0036680718ShenW.ChooY. S.Stress intensity factor for a tubular T-joint with grouted chord201235374710.1016/j.engstruct.2011.10.0142-s2.0-84855461407KottegodaN. T.RossoR.20082ndSouthampton, UKBlackwell Publishing