# The Resonance Reliability and Global Sensitivity Analysis of Curved Pipe Conveying Fluid Based on TIS Method

Authors: Zhao Yuzhen; Liu Yongshou; Guo Qing; Li Baohui
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Based on the Flügge curved beam theory and total inextensible assumption, the dynamic equations of curved pipe’s in-plane vibration are established using the Newton method. The wave propagation method is proposed for calculating the natural frequency of curved pipes with clamped-clamped supported at both ends. Then, the performance function of the resonance reliability of curved pipe conveying fluid is established. Main and total effect indices of global sensitivity analysis (GSA) are introduced. The truncated importance sampling (TIS) method is used for calculating these indices. In the example, the natural frequency and critical velocity of a semicircular pipe are calculated. The importance ranking of input variables is obtained at different working conditions. The method proposed in this paper is valuable and leads to reliability estimation and antiresonance design of curved pipe conveying fluid.


---

## Body

## 1. Introduction

Curved pipes conveying fluid are widely used in aerospace, aircraft, nuclear, petroleum, and ocean engineering. All kinds of uncertainties in fluid velocity, pressure, pipe geometry, and materials will lead to a large change on the natural frequency of the curved pipe conveying fluid through the fluid-structure interaction. When the natural frequency is close to the excitation frequency, this will lead to the resonance failure of pipes. The resonance will have a large damage to the curved pipe. So it is essential and important to study the resonance reliability analysis of curved pipe conveying fluid.Up to now, there have been a lot of researches in solving the natural frequency of curved pipes conveying fluid [1–10]. Chen S.S. [2] put forward the semicircular pipe’s vibration dynamic model. In this model, the total axial extension is none, and the Newton and Hamilton method is used to deduce the motional differential equation of curved pipes conveying fluid. Ni Q. et al. [3] have a lot of work in the dynamic researching of curved pipes. They used the differential quadrature method (DQM) to study the semicircular pipe’s vibration and stability analysis. Misra A.K. et al. [4, 5] used the finite element method (FEM) to study the dynamics and stability of curved pipes conveying fluid. Li B.H. et al. [6] put forward the wave propagation method to solve the natural frequency of curved pipes, including the total extensible assumption and total inextensible assumption. Based on the Hamilton principle, Jung D.H. et al. [7] put forward the novel expression of fluid acceleration of semicircular pipe. Hu J.Y. et al. [8] studied the fluid-conveying curved pipe with an arbitrary undeformed configuration.Although researchers have made a lot of progress in vibration analysis of pipes conveying fluid, the study of resonance failure of pipes conveying fluid is scarce [9–12]. Nokland T.E. et al. [9] provided a guideline for selecting the appropriate importance measure for different application of risk and reliability. Kvassay M. et al. [10] proposed a new approach for importance analysis of multistate Systems (MSSs). Zhai H.B. et al. [13] studied the dynamic reliability of straight pipe conveying fluid based on a refined response surface method, but the pipe is uniformly straight. Alizadeh A.-A. et al. [14] used the Monte Carlo method to study the self-vibration and stability analysis of fluid-conveying pipe. Li Y.L. et al. [15] computed the reliability of subsea under spatially varying ground motions by using subset simulation. Among these researches, only the single straight pipe was considered for calculating the resonance failure and local sensitivity analysis (LSA). Then random parameters such as pipe length and diameter and fluid density are considered and the reliability of fuel pipelines in different flight states using different gradual enlarged pipes is analyzed. Sensitivity of the fuel pipeline’s functional reliability is explored, and pipe length is chief factor with influences on functional reliability. The parameter sensitivity has been also analyzed by Monte Carlo method. They also studied the parameter sensitivity of the metal bellows based on Monte Carlo method. But the fluid-structure interaction was not considered. Ritto T.G. et al. [16] proposed a probabilistic model for the fluid-structure interaction considering modeling errors and analyzed the stability of the stochastic system. Yun W.Y. et al. [17] proposed a modification to the importance sampling in reliability analysis.The resonance failure is getting more and more serious in engineering especially in the pipe systems. But there are few papers studying the resonance failure of the pipes conveying fluid especially for the curved pipe. Because the curved pipes are widely used in engineering, studying the resonance failure of curved pipes conveying fluid is essential and important. There are many input variables affecting the failure probabilities and the uncertainty of each variable has different effects on the resonance failure probability. The sensitivity analysis aims to measure the uncertainty of each input variable on the failure probabilities including local sensitivities analysis (LSA) and global sensitivity analysis (GSA). Compared to the LSA, the GSA reflects the effect of the variable’s distribution on the failure probabilities globally. Many global sensitivity analyses such as the variance-based GSA and moment independent GSA are available. This work is proposed for reducing the failure probability of curved pipes conveying fluid, so the global reliability sensitivity indices defined by Cui L.J. [18] and Li L.Y. [19] can be used in this paper. The main and total effect indices are introduced. The failure probability can be reduced mostly by decreasing the uncertainty of input variables with high main effect indices, and the failure probabilities cannot be reduced significantly by decreasing the uncertainty of the input variables with small total effect induces close to zero [20]. This can help designers reduce the pipe’s resonance failure probabilities. This work aims at studying the effect of the uncertainty of each input variable on the resonance failure probability and deciding how to reduce the failure probability by decreasing the uncertainty of the input variables.This paper is organized as follows: Section2 introduces the wave propagation method for calculating the natural frequency of curved pipes conveying fluid. Section 3 establishes the performance function of resonance reliability of curved pipes conveying fluid. The resonance reliability and Sobol’s variance-based GSA based on the TIS method are introduced. Section 4 calculates the resonance reliability and Sobol’s variance-based main and total effect indices of the clamped-clamped supported curved pipe conveying fluid. Section 5 gives the conclusion to this work.


## 2. Motion Equations

There are two assumptions of approximating curved pipe, i.e., the total extensible assumption and total inextensible assumption. In this paper, on the base of total inextensible assumption, the Flügge curved beam was applied, and the wave propagation method is proposed for calculating the natural frequency of in-plane vibration.


### 2.1. Dynamic Equations of In-Plane Vibration

The force diagram of the pipe is sketched in Figure1. The N, Q, and M are the axial force, shear force, and bending moment along the pipe, respectively. N- and f denote normal and shear force of the unit fluid acting on the pipe, as shown in Figure 1.Figure 1
The force sketch of curved pipe.The equilibrium equation along the tangent line is as follows:(1)N+∂N∂θdθcos⁡dθ2-Ncos⁡dθ2+Q+∂Q∂θdθsin⁡dθ2+Qsin⁡dθ2-fRdθ=ρpApRdθ∂2u∂t2where ds=Rdθ, ∂/∂θ=R∂/∂s, cosdθ/2≐1, and sindθ/2≐dθ/2. Equation (1) can be simplified as(2)∂N∂s+QR-f=mp∂2u∂t2where mp=ρpAp is the unit mass of the curved pipe.The equilibrium equation along the normal line is(3)NR-∂Q∂s-N-=mp∂2w∂t2The moment equilibrium equation is(4)∂M∂θdθ+QRdθ=0Equation (4) also can be written as (5)∂M∂s+Q=0The force diagram of fluid is sketched as Figure2. In this paper, the pressure dropping along the pipe is neglected. P is the pressure on the cross section, and P=PAf. The pressure is along the direction of minus tangent line of the pipe axial, and P=-Pτ.Figure 2
The force sketch of the fluid.The relationship of coordinate systems before and after deformation is(6)τn=1uR+∂w∂s-uR-∂w∂s1τ0n0CombiningP and (6), the pressure P can be written as(7)P=Pτ0+PuR+∂w∂sn0Then, the equilibrium equation along the tangent line can be written as(8)Pcos⁡dθ2+PuR+∂w∂ssin⁡dθ2-Pcos⁡dθ2+PuR+∂w∂ssin⁡dθ2+fRdθ=ρfAfRdθaτ0Equation (8) can be simplified as(9)PRuR+∂w∂s+f=mfaτ0where mf=ρfAf is the fluid mass of unit length.The normal equilibrium equation is(10)N--PR=mfan0Combining (3), (4), (9), and (10), the fluid and pipe coupled dynamic equations are obtained as follows:(11)Tangent:∂N∂s+QR-mfaτ0+PRuR+∂w∂s=mp∂2u∂t2(12)Normal:NR-∂Q∂s-mfan0-PR=mp∂2w∂t2In this paper, the Flügge beam model is used to approximate the curved pipe. The axial forceN and the shear force Q are [21, 22](13)N=EAp∂u∂s-wR+EIR∂2w∂s2+wR2=EIR∂2w∂s2+wR2Q=EI∂3w∂s3+1R2∂w∂sBecause the axial strainε=0, the first item of axial force N is zero. E is Young’s modulus, and I is the inertia moment across the section.The fluid acceleration in the pipe can be described as(14)af=aτ0τ0+an0n0=∂2u∂t2+V∂2u∂t∂s-1R∂w∂t-V2RuR+∂w∂sτ0+∂2w∂t2+2V1R∂u∂t+∂2w∂t∂s+V21R+1R∂u∂s+∂2w∂s2n0Combining (13), (14), (11), and (12), we can get the tangent and normal dynamic equations as follows:Tangent dynamic equation:(15)2EIR∂3w∂s3+1R2∂w∂s+1RP+mfV2uR+∂w∂s-mfV∂2u∂t∂s-1R∂w∂t=mf+mp∂2u∂t2Normal dynamic equation:(16)EI∂4w∂s4-wR4+PR+2mfV1R∂u∂t+∂2w∂t∂s+mfV21R+1R∂u∂s+∂2w∂s2=-mf+mp∂2w∂t2


### 2.2. Solution of the Dynamic Equation

When the motion is simple harmonic vibration, tangential displacement and radial displacement can be written as [21, 22](17)us,t=Cuexp⁡iks+ωt,ws,t=Cwexp⁡iks+ωtwhere i=-1 is the imaginary number unit, k and ω are the wave number and frequency, and Cu and Cw are the undetermined coefficient.Substituting (17) into (15), (16) and neglecting the same factor exp⁡i(ks+ωt), we can get the matrix as follows:(18)-2EIik3R+2EIR3+P+mfV2RikP+mfV2R2+mp+mfω2EIk4-mfV2k2-2mfVωk-EIR4+mfV2R2-mp+mfω22mfViωRCwCu=0To solve the equation above, the determinant on the left side of (18) must be nonzero. And a fourth-order polynomial equation is obtained which is called to be the frequency dispersion equation.(19)k4-2βωRαk3-mfV2EIk2+2βωαR3+βωP+mfV2αREI-2mfVωEIk-1R4+mfV2R2EI-mf+mpω2EI=0where α=P+mfV2/R2+(mp+mf)ω2, β=2mfV/R.Equation (19) has four roots. Two of them are real numbers and are opposite to each other. The others are conjugate numbers for each other. They represent the wave positive and negative propagating along the pipe, respectively. In this paper, k1 and k2 are regarded as the negative propagation waves along the curved pipe. k3 and k4 are regarded as the positive propagation waves along the curved pipe.As shown in Figure3, the positive propagation wave and negative propagation wave along the semicircular pipe with clamped support are described as W+, W-. The subscripts l and r are left and right sides of the curved pipe, respectively.Figure 3
The wave mode in semicircular pipe.whereTl and Tr represent the left wave propagation and right wave propagation and Rl and Rr are the reflect matrix of left and right wave, it is easy to know that(20)Wr+=TrWl+,Wl-=TlWr-The transformation matrices of left and right wave propagation are(21)Tl=exp⁡-ik1Rπ00exp⁡-ik2Rπ,Tr=exp⁡ik3Rπ00exp⁡ik4RπThe reflection relation on the left side is(22)Wl+=RlWl-Combined with (14), (22) could be expanded as (23)Cu1+Cu2+Cu3+Cu4=0Cw1+Cw2+Cw3+Cw4=0According to the total inextensible assumption∂u/∂s-w/R=0, substituting (23) to (17), we can get(24)Cw=RikCuSubstituting (24) into (23), we can get(25)Cu1+Cu2+Cu3+Cu4=0k1Cu1+k2Cu2+k3Cu3+k4Cu4=0(26)11k1k2Wl-+11k3k4Wl+=0Substituting (22) into (26), the reflect matrix on the left side of the pipe is (27)Rl=-11k3k4-111k1k2And the reflect matrix on the right side of the pipe is(28)Rr=-11k1k2-111k3k4The wave propagation and reflection in the pipe can be expressed as follows [23]:(29)Wr+=TrWl+,Wr-=RrWr+,Wl-=TlWr-,Wl+=RlWl-The equation above can be written as(30)I-RlTlRrTrWl+=0where I is the second-order matrix. This equation which has solutions must satisfy(31)det⁡I-RlTlRrTr=0It can be seen that (31) includes ω, so the natural frequency can be solved.


## 2.1. Dynamic Equations of In-Plane Vibration

The force diagram of the pipe is sketched in Figure1. The N, Q, and M are the axial force, shear force, and bending moment along the pipe, respectively. N- and f denote normal and shear force of the unit fluid acting on the pipe, as shown in Figure 1.Figure 1
The force sketch of curved pipe.The equilibrium equation along the tangent line is as follows:(1)N+∂N∂θdθcos⁡dθ2-Ncos⁡dθ2+Q+∂Q∂θdθsin⁡dθ2+Qsin⁡dθ2-fRdθ=ρpApRdθ∂2u∂t2where ds=Rdθ, ∂/∂θ=R∂/∂s, cosdθ/2≐1, and sindθ/2≐dθ/2. Equation (1) can be simplified as(2)∂N∂s+QR-f=mp∂2u∂t2where mp=ρpAp is the unit mass of the curved pipe.The equilibrium equation along the normal line is(3)NR-∂Q∂s-N-=mp∂2w∂t2The moment equilibrium equation is(4)∂M∂θdθ+QRdθ=0Equation (4) also can be written as (5)∂M∂s+Q=0The force diagram of fluid is sketched as Figure2. In this paper, the pressure dropping along the pipe is neglected. P is the pressure on the cross section, and P=PAf. The pressure is along the direction of minus tangent line of the pipe axial, and P=-Pτ.Figure 2
The force sketch of the fluid.The relationship of coordinate systems before and after deformation is(6)τn=1uR+∂w∂s-uR-∂w∂s1τ0n0CombiningP and (6), the pressure P can be written as(7)P=Pτ0+PuR+∂w∂sn0Then, the equilibrium equation along the tangent line can be written as(8)Pcos⁡dθ2+PuR+∂w∂ssin⁡dθ2-Pcos⁡dθ2+PuR+∂w∂ssin⁡dθ2+fRdθ=ρfAfRdθaτ0Equation (8) can be simplified as(9)PRuR+∂w∂s+f=mfaτ0where mf=ρfAf is the fluid mass of unit length.The normal equilibrium equation is(10)N--PR=mfan0Combining (3), (4), (9), and (10), the fluid and pipe coupled dynamic equations are obtained as follows:(11)Tangent:∂N∂s+QR-mfaτ0+PRuR+∂w∂s=mp∂2u∂t2(12)Normal:NR-∂Q∂s-mfan0-PR=mp∂2w∂t2In this paper, the Flügge beam model is used to approximate the curved pipe. The axial forceN and the shear force Q are [21, 22](13)N=EAp∂u∂s-wR+EIR∂2w∂s2+wR2=EIR∂2w∂s2+wR2Q=EI∂3w∂s3+1R2∂w∂sBecause the axial strainε=0, the first item of axial force N is zero. E is Young’s modulus, and I is the inertia moment across the section.The fluid acceleration in the pipe can be described as(14)af=aτ0τ0+an0n0=∂2u∂t2+V∂2u∂t∂s-1R∂w∂t-V2RuR+∂w∂sτ0+∂2w∂t2+2V1R∂u∂t+∂2w∂t∂s+V21R+1R∂u∂s+∂2w∂s2n0Combining (13), (14), (11), and (12), we can get the tangent and normal dynamic equations as follows:Tangent dynamic equation:(15)2EIR∂3w∂s3+1R2∂w∂s+1RP+mfV2uR+∂w∂s-mfV∂2u∂t∂s-1R∂w∂t=mf+mp∂2u∂t2Normal dynamic equation:(16)EI∂4w∂s4-wR4+PR+2mfV1R∂u∂t+∂2w∂t∂s+mfV21R+1R∂u∂s+∂2w∂s2=-mf+mp∂2w∂t2


## 2.2. Solution of the Dynamic Equation

When the motion is simple harmonic vibration, tangential displacement and radial displacement can be written as [21, 22](17)us,t=Cuexp⁡iks+ωt,ws,t=Cwexp⁡iks+ωtwhere i=-1 is the imaginary number unit, k and ω are the wave number and frequency, and Cu and Cw are the undetermined coefficient.Substituting (17) into (15), (16) and neglecting the same factor exp⁡i(ks+ωt), we can get the matrix as follows:(18)-2EIik3R+2EIR3+P+mfV2RikP+mfV2R2+mp+mfω2EIk4-mfV2k2-2mfVωk-EIR4+mfV2R2-mp+mfω22mfViωRCwCu=0To solve the equation above, the determinant on the left side of (18) must be nonzero. And a fourth-order polynomial equation is obtained which is called to be the frequency dispersion equation.(19)k4-2βωRαk3-mfV2EIk2+2βωαR3+βωP+mfV2αREI-2mfVωEIk-1R4+mfV2R2EI-mf+mpω2EI=0where α=P+mfV2/R2+(mp+mf)ω2, β=2mfV/R.Equation (19) has four roots. Two of them are real numbers and are opposite to each other. The others are conjugate numbers for each other. They represent the wave positive and negative propagating along the pipe, respectively. In this paper, k1 and k2 are regarded as the negative propagation waves along the curved pipe. k3 and k4 are regarded as the positive propagation waves along the curved pipe.As shown in Figure3, the positive propagation wave and negative propagation wave along the semicircular pipe with clamped support are described as W+, W-. The subscripts l and r are left and right sides of the curved pipe, respectively.Figure 3
The wave mode in semicircular pipe.whereTl and Tr represent the left wave propagation and right wave propagation and Rl and Rr are the reflect matrix of left and right wave, it is easy to know that(20)Wr+=TrWl+,Wl-=TlWr-The transformation matrices of left and right wave propagation are(21)Tl=exp⁡-ik1Rπ00exp⁡-ik2Rπ,Tr=exp⁡ik3Rπ00exp⁡ik4RπThe reflection relation on the left side is(22)Wl+=RlWl-Combined with (14), (22) could be expanded as (23)Cu1+Cu2+Cu3+Cu4=0Cw1+Cw2+Cw3+Cw4=0According to the total inextensible assumption∂u/∂s-w/R=0, substituting (23) to (17), we can get(24)Cw=RikCuSubstituting (24) into (23), we can get(25)Cu1+Cu2+Cu3+Cu4=0k1Cu1+k2Cu2+k3Cu3+k4Cu4=0(26)11k1k2Wl-+11k3k4Wl+=0Substituting (22) into (26), the reflect matrix on the left side of the pipe is (27)Rl=-11k3k4-111k1k2And the reflect matrix on the right side of the pipe is(28)Rr=-11k1k2-111k3k4The wave propagation and reflection in the pipe can be expressed as follows [23]:(29)Wr+=TrWl+,Wr-=RrWr+,Wl-=TlWr-,Wl+=RlWl-The equation above can be written as(30)I-RlTlRrTrWl+=0where I is the second-order matrix. This equation which has solutions must satisfy(31)det⁡I-RlTlRrTr=0It can be seen that (31) includes ω, so the natural frequency can be solved.


## 3. Resonance Reliability and GSA

### 3.1. Performance Function

In order to avoid resonance behavior, the excitation frequency must be far away from the natural frequency of the curved pipe conveying fluid [24, 25]. The resonance failure probability analysis aims to estimate the failure probability of the structure.In the reality engineering, due to the uncertainties among materials, manufacturing, installing, and servicing, the excited frequency, natural frequency, and vibration response are random variables. According to the traditional design criterion, there exist risks in resonance design if1-k1<RR/S<1+k2, where RR is the natural frequency, S is the excited frequency, and k1 and k2 are the constants based on the materials. For different materials, k1 and k2 are different. Now it is assumed that S1=(1-k1)S, S2=(1+k2)S, where S1 and S2 are the minimum and maximum of excitation circular frequency. The performance function of pipe resonance failure [13]: (32)zjjx=Rjj-S1·Rjj-S2,jj=1,2,⋯,mzx=z1x,z2x,⋯zmx.where zjj(x) is thejj-th performance function, x is the basic variables, Rjj is thejj-th natural frequency, and m is the number of orders and is also the failure mode number.The resonance reliability can also be expressed as follows:(33)Pr=Pzx>0=∫zx>0fXx,θxdx.where fX is the probability density function (PDF) of x and θx are the distribution indices of variable x.Equation (33) is a nonlinear invisible performance function, so it could not be solved analytically. The equations could be solved by numerical or approximated approach.


### 3.2. Global Sensitivity Indices

Suppose the limit state function is given asY=g(x), where x=x1,x2,⋯,xn is the random input variable vector with joint PDF.The failure domain of this structure system is defined as(34)F=x:gx<0Suppose the indicator function of this failure domain is given asIFx(35)IFx=1x∈F0x∉FThen the failure probability can be expressed as(36)Pf=EIF=∫RnIFxfxdx=∫FfxdxA modified version of global reliability sensitivity index defined by Cui in [18] is given as(37)δiP=EPf-Pf∣Xi2=∫-∞+∞PfY-Pf∣Xi2fixidxiwhere Pf∣Xi is the failure probability conditional on Xi.Further, Li proved that(38)δiP=EPf-Pf∣Xi2=VEIF∣XiIn GSA, the global sensitivity indices of a group of input variables are defined as main and total effect indices.The main effect index of the single input variable is defined as(39)Si=VEIF∣XiVIFThe total effect index of the single input variables is defined as(40)STi=1-VEIF∣X~iVIF=EVIF∣X~iVIFwhere vector X~i is the input variables excluding Xi.The Sobol indices satisfied 0≤Si≤STi≤1.It is easy to know that the larger the main effect indexSi of one input variable, the higher the contribution of input variable Xi on the failure probability Pf. Reducing the uncertainty of Xi, the reduction of failure probability is much more. Obviously, the smaller the total effect index STi of one input variable, the lower the contribution of input variables Xi on the failure probability Pf. If the total effect indices of one input variable are small enough or close to zero, the influence of this input variable can be ignored to reduce the dimensionality and save the calculating cost.Due to the well-known low total expectation,(41)Si=VEIF∣xiVIF=EE2IF∣xi-E2EIF∣xiPf-Pf2=EE2IF∣xi-Pf2Pf-Pf2=Di-Pf2D(42)STi=1-VEIF∣x~iVIF=1-EE2IF∣x~i-E2EIF∣x~iPf-Pf2=1-EE2IF∣x~i-Pf2Pf-Pf2=1-D~i-Pf2Dwhere D=VIF=Pf-Pf2, Di=EE2IFx~i, D~i=EE2IFx~i. Thus for computing Si and STi, we need to calculate D, Pf, Pf2, Di, and D~i.


### 3.3. TIS for ComputingSi and STi

In this paper, the TIS procedure is used to calculate the resonance failure probability. Compared with the other procedure, the TIS procedure can further improve the efficiency of solving.The basic idea of TIS is as follows: in the standard normal space, the most probable point (MPP)xMPP is computed by the advanced first-order second-moment (AFOSM) procedure; the reliability index is the distance from the MPP to the coordinate origin. Then a hypersphere with radius of β can be obtained, denoted as β-sphere. The indicator function of the outer space of the β-sphere is defined as follows [26]:(43)Iβ=0,x2<β21,x2≥β2The PDF is truncated by theβ-sphere, and the truncated importance sampling density function can be written as(44)htrx=0,x2<β21Phtrhxx,x2≥β2We can rewrite (33) by introducing importance sampling PDF as follows:(45)Pf=∫⋯∫RnIFxfxxdx=∫⋯∫RnIFxfxxhxxhxxdxFurthermore, by inducing the indicator functionIβx, the failure probability can be rewritten as(46)Pf=∫⋯∫RnIFxfxxdx=∫⋯∫RnIFxIβxfxxhxxhxxdx=EIFxIβxfxxhxxwhere ω2(x)=IF(x)Iβ(x)fx/hx and Pf=Eω2(x). So the term Pf2 could be rewritten as(47)Pf2=∫R2nIFxIβxfxhxhxIFx′Iβx′fx′hx′hx′dxdx′=EIFxIβxfxhxIFx′Iβx′fx′hx′Generating N samples according to importance sampling PDF, the mean value of samples can be regarded as the expectation value.P^f which is the estimated value of Pf could be computed by(48)P^f=1N∑k=1NIFxkIβxkfxxkhxxkThe TIS method just needs to calculate the samples which are in the outer space of the hypersphere. And the traditional method needs to calculate the limit state function values of all sampling points. So the TIS method is more efficient than the traditional method.Above all, the following five steps based on TIS can be implemented for computing theSi and STi.Step 1.
Generate2N samples of the input variables by h(x), and distribute them to two N×n sample matrices A2 and B2.Step 2.
Generate anotherN×n sample matrix C2(i). Each column is equal to the corresponding column of B2 except thei-th column, which is equal toi-th column of A2.Step 3.
Compute the values of functionsω2(x), F3, and F4 for each sample in A2, B2, and C2(i); then we can obtain four N-dimensional vectors, i.e.,(49)ω2A=ω2A2,ω2B=ω2B2,F3A=F3A2,C2i,F4A=F4B2,C2iStep 4.
ComputeSi and STi as follows:(50)S^i=D^i-P^f2D^,S^Ti=1-D^~i-P^f2D^where P^f=1/2N∑j=1Nω2Aj+ω2Bj, P^f2=1/N∑j=1Nω2Ajω2Bj, D^=P^f-P^f2, and D^i=1/N∑j=1NF3A(j), D^~i=1/N∑j=1NF4B(j). ω2Aj, ω2Bj, F3A(j), and F4B(j) are thej-th element of ω2A, ω2B, F3A, and F4B, respectively. In this step, we first judge if the sample is contained in the β-sphere and if so, then we do not need to compute the limit state function value for this sample.Step 5.
Compute the probable errors ofS^i and S^Ti:(51)δS^i≈D^δD^i+D^iδD^D^2,δS^Ti≈D^δD^~i+D^~iδD^D^2where(52)δD^≈0.6745N1N∑j=1Nω2Aj+ω2Bj22-1N∑j=1Nω2Aj+ω2Bj22(53)δD^i≈0.6745N1N∑j=1NF3Aj2-1N∑j=1NF3Aj2(54)δD^~i≈0.6745N1N∑j=1NF4Bj2-1N∑j=1NF4Bj2


## 3.1. Performance Function

In order to avoid resonance behavior, the excitation frequency must be far away from the natural frequency of the curved pipe conveying fluid [24, 25]. The resonance failure probability analysis aims to estimate the failure probability of the structure.In the reality engineering, due to the uncertainties among materials, manufacturing, installing, and servicing, the excited frequency, natural frequency, and vibration response are random variables. According to the traditional design criterion, there exist risks in resonance design if1-k1<RR/S<1+k2, where RR is the natural frequency, S is the excited frequency, and k1 and k2 are the constants based on the materials. For different materials, k1 and k2 are different. Now it is assumed that S1=(1-k1)S, S2=(1+k2)S, where S1 and S2 are the minimum and maximum of excitation circular frequency. The performance function of pipe resonance failure [13]: (32)zjjx=Rjj-S1·Rjj-S2,jj=1,2,⋯,mzx=z1x,z2x,⋯zmx.where zjj(x) is thejj-th performance function, x is the basic variables, Rjj is thejj-th natural frequency, and m is the number of orders and is also the failure mode number.The resonance reliability can also be expressed as follows:(33)Pr=Pzx>0=∫zx>0fXx,θxdx.where fX is the probability density function (PDF) of x and θx are the distribution indices of variable x.Equation (33) is a nonlinear invisible performance function, so it could not be solved analytically. The equations could be solved by numerical or approximated approach.


## 3.2. Global Sensitivity Indices

Suppose the limit state function is given asY=g(x), where x=x1,x2,⋯,xn is the random input variable vector with joint PDF.The failure domain of this structure system is defined as(34)F=x:gx<0Suppose the indicator function of this failure domain is given asIFx(35)IFx=1x∈F0x∉FThen the failure probability can be expressed as(36)Pf=EIF=∫RnIFxfxdx=∫FfxdxA modified version of global reliability sensitivity index defined by Cui in [18] is given as(37)δiP=EPf-Pf∣Xi2=∫-∞+∞PfY-Pf∣Xi2fixidxiwhere Pf∣Xi is the failure probability conditional on Xi.Further, Li proved that(38)δiP=EPf-Pf∣Xi2=VEIF∣XiIn GSA, the global sensitivity indices of a group of input variables are defined as main and total effect indices.The main effect index of the single input variable is defined as(39)Si=VEIF∣XiVIFThe total effect index of the single input variables is defined as(40)STi=1-VEIF∣X~iVIF=EVIF∣X~iVIFwhere vector X~i is the input variables excluding Xi.The Sobol indices satisfied 0≤Si≤STi≤1.It is easy to know that the larger the main effect indexSi of one input variable, the higher the contribution of input variable Xi on the failure probability Pf. Reducing the uncertainty of Xi, the reduction of failure probability is much more. Obviously, the smaller the total effect index STi of one input variable, the lower the contribution of input variables Xi on the failure probability Pf. If the total effect indices of one input variable are small enough or close to zero, the influence of this input variable can be ignored to reduce the dimensionality and save the calculating cost.Due to the well-known low total expectation,(41)Si=VEIF∣xiVIF=EE2IF∣xi-E2EIF∣xiPf-Pf2=EE2IF∣xi-Pf2Pf-Pf2=Di-Pf2D(42)STi=1-VEIF∣x~iVIF=1-EE2IF∣x~i-E2EIF∣x~iPf-Pf2=1-EE2IF∣x~i-Pf2Pf-Pf2=1-D~i-Pf2Dwhere D=VIF=Pf-Pf2, Di=EE2IFx~i, D~i=EE2IFx~i. Thus for computing Si and STi, we need to calculate D, Pf, Pf2, Di, and D~i.


## 3.3. TIS for ComputingSi and STi

In this paper, the TIS procedure is used to calculate the resonance failure probability. Compared with the other procedure, the TIS procedure can further improve the efficiency of solving.The basic idea of TIS is as follows: in the standard normal space, the most probable point (MPP)xMPP is computed by the advanced first-order second-moment (AFOSM) procedure; the reliability index is the distance from the MPP to the coordinate origin. Then a hypersphere with radius of β can be obtained, denoted as β-sphere. The indicator function of the outer space of the β-sphere is defined as follows [26]:(43)Iβ=0,x2<β21,x2≥β2The PDF is truncated by theβ-sphere, and the truncated importance sampling density function can be written as(44)htrx=0,x2<β21Phtrhxx,x2≥β2We can rewrite (33) by introducing importance sampling PDF as follows:(45)Pf=∫⋯∫RnIFxfxxdx=∫⋯∫RnIFxfxxhxxhxxdxFurthermore, by inducing the indicator functionIβx, the failure probability can be rewritten as(46)Pf=∫⋯∫RnIFxfxxdx=∫⋯∫RnIFxIβxfxxhxxhxxdx=EIFxIβxfxxhxxwhere ω2(x)=IF(x)Iβ(x)fx/hx and Pf=Eω2(x). So the term Pf2 could be rewritten as(47)Pf2=∫R2nIFxIβxfxhxhxIFx′Iβx′fx′hx′hx′dxdx′=EIFxIβxfxhxIFx′Iβx′fx′hx′Generating N samples according to importance sampling PDF, the mean value of samples can be regarded as the expectation value.P^f which is the estimated value of Pf could be computed by(48)P^f=1N∑k=1NIFxkIβxkfxxkhxxkThe TIS method just needs to calculate the samples which are in the outer space of the hypersphere. And the traditional method needs to calculate the limit state function values of all sampling points. So the TIS method is more efficient than the traditional method.Above all, the following five steps based on TIS can be implemented for computing theSi and STi.Step 1.
Generate2N samples of the input variables by h(x), and distribute them to two N×n sample matrices A2 and B2.Step 2.
Generate anotherN×n sample matrix C2(i). Each column is equal to the corresponding column of B2 except thei-th column, which is equal toi-th column of A2.Step 3.
Compute the values of functionsω2(x), F3, and F4 for each sample in A2, B2, and C2(i); then we can obtain four N-dimensional vectors, i.e.,(49)ω2A=ω2A2,ω2B=ω2B2,F3A=F3A2,C2i,F4A=F4B2,C2iStep 4.
ComputeSi and STi as follows:(50)S^i=D^i-P^f2D^,S^Ti=1-D^~i-P^f2D^where P^f=1/2N∑j=1Nω2Aj+ω2Bj, P^f2=1/N∑j=1Nω2Ajω2Bj, D^=P^f-P^f2, and D^i=1/N∑j=1NF3A(j), D^~i=1/N∑j=1NF4B(j). ω2Aj, ω2Bj, F3A(j), and F4B(j) are thej-th element of ω2A, ω2B, F3A, and F4B, respectively. In this step, we first judge if the sample is contained in the β-sphere and if so, then we do not need to compute the limit state function value for this sample.Step 5.
Compute the probable errors ofS^i and S^Ti:(51)δS^i≈D^δD^i+D^iδD^D^2,δS^Ti≈D^δD^~i+D^~iδD^D^2where(52)δD^≈0.6745N1N∑j=1Nω2Aj+ω2Bj22-1N∑j=1Nω2Aj+ω2Bj22(53)δD^i≈0.6745N1N∑j=1NF3Aj2-1N∑j=1NF3Aj2(54)δD^~i≈0.6745N1N∑j=1NF4Bj2-1N∑j=1NF4Bj2


## 4. Example

The curved pipe is sketched in Figure4, and the pipe is clamped-clamped supported at both ends. The materials and geometry are listed in Table 1. Each input variable is independent of others and has the same coefficient of variance.Table 1
The distribution parameter of the input variables of example.
Input variables/unit Distribution Mean value variance Pipe Yang’s modulus(GPa) Gauss 68.9 3.43 Pipe outer diameter(m) Gauss 0.0232 0.00116 Pipe thickness(m) Gauss 0.0041 0.000205 Pipe density(kg/m3) Gauss 7197 359.85 Fluid density(kg/m3) Gauss 1000 50 Velocity(m/s) Gauss 0 0.05 Maximum excitation frequency(rad/s) Gauss 410 20.5 Minimum excitation frequency(rad/s) Gauss 210 10.5 Radius of curvature(m) Gauss 0.5 0.025 Pressure(MPa) Gauss 10 0.5Figure 4
The curved pipe illustration of supported.


### 4.1. Natural Frequency and Critical Velocity

The first four natural frequencies are listed in Table2. It is shown from Figure 5 that with increasing of velocity, the natural frequency decreases slowly until the point about V=135.1m/s, where the natural frequency falls to zero suddenly. It means that the critical velocity is Vcr=135.1m/s.Table 2
The relationship of four natural frequencies and velocities.
V (m/s) First (rad/s) Second (rad/s) Third (rad/s) Fourth (rad/s) 0 165.46 501.67 992.99 1644.95 10 165.41 501.53 992.83 1644.78 30 165.01 500.42 991.51 1643.36 50 164.22 498.19 988.88 1640.52 100 160.65 487.64 976.43 1627.14Figure 5
The first natural frequency and velocity relationship.


### 4.2. The Resonance Reliability Analysis

The input variables are listed in Table1, and all the input variables follow the normal distribution. The coefficient of variance ζ is 0.05, the pressure p is 10MPa, the maximum excitation circular frequency s2 is 410rad/s, and the minimum excitation circular frequency s1 is 210rad/s. Here, different working conditions are considered, v=0m/s, v=10m/s, and v=20m/s. The resonance failure probabilities under different conditions are listed in Table 3.Table 3
The resonance failure probability of curved pipe conveying fluid at different velocity conditions.
Working condition First failure probability Second failure probability v=0 m/s 0.0071463 0.0231528 v=10 m/s 0.0070448 0.0308182 v=20 m/s 0.0062525 0.0373823It can be seen from Table3 that, with increasing of velocity, the first-order resonance failure probability decreases, and the second-order resonance failure probability becomes larger. This is because the natural frequency is decreasing along with the increasing of velocity as shown in Table 2. So the first-order failure area is decreasing, and the second-order failure area is increasing. The results in this example coincide with the theoretical analysis.


### 4.3. GSA

Table4 shows the Sobol indices with different coefficients of variance 0.05 and 0.1. Under the working condition, fluid velocity v=0m/s, pressure p=10MPa, the range of excitation circular frequency s1-s2:210-410rad/s, and the first natural frequency resonance failure probabilities are calculated with the coefficients of variance 0.05 and 0.1. The results are listed in Table 4, where the resonance failure probabilities are Pf,ζ=0.05=7.0448E-4 and Pf,ζ=0.1=0.1445, respectively. The results demonstrate that, with the increasing of variance, the failure probability is increasing sharply. The main effect indices are getting larger with the increasing of variance from 0.05 to 0.1. The total effect indices are getting smaller with the increasing of variance. This is because the failure area is getting larger with the increasing of variance, and the resonance failure probability is getting bigger. According to the definition of total effect indices, the total effect indices are inversely proportional to the main effect indices, as demonstrated in Figure 6.Table 4
The sobol indices with different coefficients of variance.
ζ index E w d h d p d f d v s 1 s 2 R p 0.05 S ^ i(E-5) 0.6528 8.4508 6.6361 477.4792 6.1049 0.2442 8414.9 9.0357 475.8557 0.3755 S ^ T i 0.4032 0.5625 0.2930 0.2784 0.3987 0.4590 0.9991 0.4626 0.8455 0.2553 0.1 S ^ i 0.0796 0.1262 0.00085 0.09031 0.00043 0 0.1038 0 0.0697 0 S ^ T i 0.09796 0.38659 0.0319 0.20620 0.08094 0.0983 0.80845 0.0890 0.58185 0.09839Figure 6
The Sobol indices of different coefficient of variance: main effect indices (a) and total effect indices (b).
(a) (b)From Figure6, it can be seen that importance ranking of the main effect indices is S1>pd>R>S2>wd>hd>fd>E>p>v. It indicates that, by decreasing the uncertainty of s1, the failure probability can be reduced most, which is followed by pipe density pd, while the last one is radius of curvature R. The other main effect indices are nearly zero. The failure probability will not change as we reduce the variance of input variable which has small main effect indices close to zero. In the reliability design of curved pipe, we can omit the uncertainty of the input variable which has small main effect indices. In this way we can reduce the dimension of model inputs. The importance ranking based on total effect indices is obtained as follows: S1>R>wd>S2>v>E>fd>hd>pd>p, and in the total effect indices, none of the indices is close to zero. Obviously, the total indices are larger than main indices, which validates that the interactions along the variables affect the failure probability greatly.When the coefficient of variance is 0.1, the main effect indices change. Particularly the main index of out pipe diameterwd changes quite a lot. It is demonstrated that the uncertainty of wd will have much more influence on the failure probability. This is because the change of wd will have a lot of influence on the moment of inertia of the cross section. And the changes of moment of inertia will make the mass term change a lot in the motion equations.Figure7 shows the Sobol indices with different velocities v=0m/s, v=10m/s, and v=30m/s. From Figure 7, it is shown that the velocity has little effect on the main and total effect indices, because the natural frequency is not obviously relevant to the velocity.Figure 7
The Sobol indices of different velocity: main effect indices (a) and total effect indices (b).
(a) (b)When the velocity increases from 0 m/s to 30 m/s, the main and total effect indices remain unchanged almost. The importance ranking under three velocity conditions still remains asS1>pd>R>S2>wd>hd>fd>E>p>v. The other main effect indices are close to zero. It indicates that, by decreasing the uncertainty of s1, the failure probability can be reduced most, which is followed by pipe density pd, while the last one is radius of curvature R.From Table5, the total effect indices stay unchanged with the increasing of velocity. In the total effect indices, none of the indices is close to zero. Obviously, the total indices are larger than main indices, which validates that the interactions along the variables affect the failure probability greatly.Table 5
The resonance reliability GSA with different velocitiesv.
v index E w d h d p d f d v s 1 s 2 R p 0m/s S ^ i 6.651E-6 8.603E-5 6.752 E-5 4.834E-3 6.211 E-5 2.494E-6 8.442E-2 9.195 E-5 4.8215E-3 3.817E-6 S ^ T i 0.4028 0.5619 0.2926 0.2781 0.3981 0.4582 0.9991 0.4619 0.8452 0.2549 10m/s S ^ i 6.528e-6 8.450E-5 6.636e-5 4.774E-3 6.105E-5 2.442E-6 8.415E-2 9.035E-5 4.7586E-3 3.755E-6 S ^ T i 0.40322 0.56258 0.29308 0.27849 0.39877 0.459 0.99913 0.46267 0.8456 0.25537 30m/s S ^ i 5.6376e-6 7.3446e-5 5.795e-5 0.0043393 5.3316e-5 2.0692e-6 0.078957 7.8782e-5 0.0043244 3.301E-6 S ^ T i 0.40583 0.58354 0.29737 0.27726 0.4016 0.46314 0.99923 0.47664 0.84785 0.25557Table6 shows the main and total effect indices with different two pressures: 0MPa and 10MPa. When the pressure is p=0MPa, the importance ranking of main effect indices is S1>R>pd>fd>hd>wd>v>S2>E>p. It indicates that, by decreasing the uncertainty of s1, the failure probability can be reduced most, which is followed by radius of curvature R, while the last one is pipe density pd. The importance ranking of total effect indices is S1>R>v>hd>fd>E>S2>p>wd>pd, and none of the input variables is close to zero. Obviously, the total indices are larger than main indices, which validates that the interactions along the variables affect the failure probability greatly.Table 6
The main and total effect indices with different pressures.
p index E w d h d p d f d v s 1 s 2 R p 0 Mpa S ^ i 0 1.935e-85 1.955E-85 2.385E-85 2.35E-85 1.35E-87 0.34995 5.88E-100 4.34E-85 0 S ^ T i 0.1110 0.06435 0.18587 0.05004 0.111135 0.26033 1 0.11100 0.26248 0.11100 10Mpa S ^ i 6.528e-6 8.4508e-5 6.6361e-5 0.0047747 6.1049e-5 2.4422e-6 0.08414 9.0357e-5 0.004758 3.7551e-6 S ^ T i 0.4032 0.5625 0.29307 0.2784 0.3987 0.4590 0.9991 0.4626 0.84559 0.25537From Figure8, it could be concluded that, with the increasing of pressure, the importance ranking will have a big change. Because the pressure has much influence on the pipe’s dynamic equation, when the pressure fluctuates, the natural frequency and the resonance failure probability will change a lot.Figure 8
The Sobol indices with different pressure: main effect indices (a) and total effect indices (b).
(a) (b)


## 4.1. Natural Frequency and Critical Velocity

The first four natural frequencies are listed in Table2. It is shown from Figure 5 that with increasing of velocity, the natural frequency decreases slowly until the point about V=135.1m/s, where the natural frequency falls to zero suddenly. It means that the critical velocity is Vcr=135.1m/s.Table 2
The relationship of four natural frequencies and velocities.
V (m/s) First (rad/s) Second (rad/s) Third (rad/s) Fourth (rad/s) 0 165.46 501.67 992.99 1644.95 10 165.41 501.53 992.83 1644.78 30 165.01 500.42 991.51 1643.36 50 164.22 498.19 988.88 1640.52 100 160.65 487.64 976.43 1627.14Figure 5
The first natural frequency and velocity relationship.


## 4.2. The Resonance Reliability Analysis

The input variables are listed in Table1, and all the input variables follow the normal distribution. The coefficient of variance ζ is 0.05, the pressure p is 10MPa, the maximum excitation circular frequency s2 is 410rad/s, and the minimum excitation circular frequency s1 is 210rad/s. Here, different working conditions are considered, v=0m/s, v=10m/s, and v=20m/s. The resonance failure probabilities under different conditions are listed in Table 3.Table 3
The resonance failure probability of curved pipe conveying fluid at different velocity conditions.
Working condition First failure probability Second failure probability v=0 m/s 0.0071463 0.0231528 v=10 m/s 0.0070448 0.0308182 v=20 m/s 0.0062525 0.0373823It can be seen from Table3 that, with increasing of velocity, the first-order resonance failure probability decreases, and the second-order resonance failure probability becomes larger. This is because the natural frequency is decreasing along with the increasing of velocity as shown in Table 2. So the first-order failure area is decreasing, and the second-order failure area is increasing. The results in this example coincide with the theoretical analysis.


## 4.3. GSA

Table4 shows the Sobol indices with different coefficients of variance 0.05 and 0.1. Under the working condition, fluid velocity v=0m/s, pressure p=10MPa, the range of excitation circular frequency s1-s2:210-410rad/s, and the first natural frequency resonance failure probabilities are calculated with the coefficients of variance 0.05 and 0.1. The results are listed in Table 4, where the resonance failure probabilities are Pf,ζ=0.05=7.0448E-4 and Pf,ζ=0.1=0.1445, respectively. The results demonstrate that, with the increasing of variance, the failure probability is increasing sharply. The main effect indices are getting larger with the increasing of variance from 0.05 to 0.1. The total effect indices are getting smaller with the increasing of variance. This is because the failure area is getting larger with the increasing of variance, and the resonance failure probability is getting bigger. According to the definition of total effect indices, the total effect indices are inversely proportional to the main effect indices, as demonstrated in Figure 6.Table 4
The sobol indices with different coefficients of variance.
ζ index E w d h d p d f d v s 1 s 2 R p 0.05 S ^ i(E-5) 0.6528 8.4508 6.6361 477.4792 6.1049 0.2442 8414.9 9.0357 475.8557 0.3755 S ^ T i 0.4032 0.5625 0.2930 0.2784 0.3987 0.4590 0.9991 0.4626 0.8455 0.2553 0.1 S ^ i 0.0796 0.1262 0.00085 0.09031 0.00043 0 0.1038 0 0.0697 0 S ^ T i 0.09796 0.38659 0.0319 0.20620 0.08094 0.0983 0.80845 0.0890 0.58185 0.09839Figure 6
The Sobol indices of different coefficient of variance: main effect indices (a) and total effect indices (b).
(a) (b)From Figure6, it can be seen that importance ranking of the main effect indices is S1>pd>R>S2>wd>hd>fd>E>p>v. It indicates that, by decreasing the uncertainty of s1, the failure probability can be reduced most, which is followed by pipe density pd, while the last one is radius of curvature R. The other main effect indices are nearly zero. The failure probability will not change as we reduce the variance of input variable which has small main effect indices close to zero. In the reliability design of curved pipe, we can omit the uncertainty of the input variable which has small main effect indices. In this way we can reduce the dimension of model inputs. The importance ranking based on total effect indices is obtained as follows: S1>R>wd>S2>v>E>fd>hd>pd>p, and in the total effect indices, none of the indices is close to zero. Obviously, the total indices are larger than main indices, which validates that the interactions along the variables affect the failure probability greatly.When the coefficient of variance is 0.1, the main effect indices change. Particularly the main index of out pipe diameterwd changes quite a lot. It is demonstrated that the uncertainty of wd will have much more influence on the failure probability. This is because the change of wd will have a lot of influence on the moment of inertia of the cross section. And the changes of moment of inertia will make the mass term change a lot in the motion equations.Figure7 shows the Sobol indices with different velocities v=0m/s, v=10m/s, and v=30m/s. From Figure 7, it is shown that the velocity has little effect on the main and total effect indices, because the natural frequency is not obviously relevant to the velocity.Figure 7
The Sobol indices of different velocity: main effect indices (a) and total effect indices (b).
(a) (b)When the velocity increases from 0 m/s to 30 m/s, the main and total effect indices remain unchanged almost. The importance ranking under three velocity conditions still remains asS1>pd>R>S2>wd>hd>fd>E>p>v. The other main effect indices are close to zero. It indicates that, by decreasing the uncertainty of s1, the failure probability can be reduced most, which is followed by pipe density pd, while the last one is radius of curvature R.From Table5, the total effect indices stay unchanged with the increasing of velocity. In the total effect indices, none of the indices is close to zero. Obviously, the total indices are larger than main indices, which validates that the interactions along the variables affect the failure probability greatly.Table 5
The resonance reliability GSA with different velocitiesv.
v index E w d h d p d f d v s 1 s 2 R p 0m/s S ^ i 6.651E-6 8.603E-5 6.752 E-5 4.834E-3 6.211 E-5 2.494E-6 8.442E-2 9.195 E-5 4.8215E-3 3.817E-6 S ^ T i 0.4028 0.5619 0.2926 0.2781 0.3981 0.4582 0.9991 0.4619 0.8452 0.2549 10m/s S ^ i 6.528e-6 8.450E-5 6.636e-5 4.774E-3 6.105E-5 2.442E-6 8.415E-2 9.035E-5 4.7586E-3 3.755E-6 S ^ T i 0.40322 0.56258 0.29308 0.27849 0.39877 0.459 0.99913 0.46267 0.8456 0.25537 30m/s S ^ i 5.6376e-6 7.3446e-5 5.795e-5 0.0043393 5.3316e-5 2.0692e-6 0.078957 7.8782e-5 0.0043244 3.301E-6 S ^ T i 0.40583 0.58354 0.29737 0.27726 0.4016 0.46314 0.99923 0.47664 0.84785 0.25557Table6 shows the main and total effect indices with different two pressures: 0MPa and 10MPa. When the pressure is p=0MPa, the importance ranking of main effect indices is S1>R>pd>fd>hd>wd>v>S2>E>p. It indicates that, by decreasing the uncertainty of s1, the failure probability can be reduced most, which is followed by radius of curvature R, while the last one is pipe density pd. The importance ranking of total effect indices is S1>R>v>hd>fd>E>S2>p>wd>pd, and none of the input variables is close to zero. Obviously, the total indices are larger than main indices, which validates that the interactions along the variables affect the failure probability greatly.Table 6
The main and total effect indices with different pressures.
p index E w d h d p d f d v s 1 s 2 R p 0 Mpa S ^ i 0 1.935e-85 1.955E-85 2.385E-85 2.35E-85 1.35E-87 0.34995 5.88E-100 4.34E-85 0 S ^ T i 0.1110 0.06435 0.18587 0.05004 0.111135 0.26033 1 0.11100 0.26248 0.11100 10Mpa S ^ i 6.528e-6 8.4508e-5 6.6361e-5 0.0047747 6.1049e-5 2.4422e-6 0.08414 9.0357e-5 0.004758 3.7551e-6 S ^ T i 0.4032 0.5625 0.29307 0.2784 0.3987 0.4590 0.9991 0.4626 0.84559 0.25537From Figure8, it could be concluded that, with the increasing of pressure, the importance ranking will have a big change. Because the pressure has much influence on the pipe’s dynamic equation, when the pressure fluctuates, the natural frequency and the resonance failure probability will change a lot.Figure 8
The Sobol indices with different pressure: main effect indices (a) and total effect indices (b).
(a) (b)


## 5. Conclusion

In this paper, the uncertainty of input variables on the natural frequency of the curved pipe conveying fluid is considered. The resonance failure probability and Sobol’s GSA are calculated based on TIS method. By decreasing the uncertainty of input variables with high main effect indices, the most reduction of failure probability can be obtained. By decreasing the uncertainty of the input variables with small total effect indices close to zero, the failure probability will not be reduced significantly. Different velocity and pressure working conditions are calculated, respectively.According to the results of natural frequency, with the increasing of velocity, the natural frequency decreases slowly until the velocityV=135.1m/s, where the natural frequency falls to zero suddenly. So the critical velocity of this example is Vcr=135.1m/s. From the main and total effect indices, it can be seen that, reducing the uncertainty of minimum excitation frequency s1, the radius of curvature R, the outer diameter wd, and the resonance failure probability can be reduced most. The method proposed in this paper will help the designers to reduce the resonance failure probability by decreasing the uncertainty of the input variables and will also be guiding for antiresonance designing of curved pipe conveying fluid.



---

*Source: 3980790-2018-11-28.xml*