Mathcad - Hyrotestv11
Mathcad - Hyrotestv11
Mathcad - Hyrotestv11
1. INPUT PARAMETER
−3 C≡ K
pcf := lb ⋅ ft
Equivalent condition
Phase :Hydrotest
Wave/Current Data : 1 yr return period wave + 1 yr return period current
Soil Type : Clay
−5 2 −1
• Kinematic Viscosity of Seawater ν := 1.076 ⋅ 10 ft s
2. CALCULATIONS
Elastic Modulus EI := E ⋅ ⎢ (
⎡ π ⋅ Ds4 − ID4 ⎤
⎥ )
⎣ 64 ⎦
9 2
EI = 1.58 × 10 ft ⋅ lbf
Inertia I :=
π
64
( 4
⋅ Ds − ID
4 )
3 4
I = 7.583 × 10 in
Weight of Pipe
lb
Wst = 209.652
ft
π
⋅ ⎡⎣( Ds + 2 ⋅ tcorr) − Ds ⎤⎦ ⋅ ρcorr
2 2
Corrosion Coating Weight Wcorr :=
4
lb
Wcorr = 4.857
ft
π
⋅ ⎡⎣ Dtot − ( Ds + 2 ⋅ tcorr)
2 2⎤
Concrete Coating Weight Wcc := ⎦ ⋅ ρcc
4
lb
Wcc = 462.058
ft
π 2
Content Weight Wcont := ⋅ ID ⋅ ρcont
4
lb
Wcont = 193.174
ft
π 2
Added Mass/Displaced water Wadd := ⋅ Dtot ⋅ ρsw
4
lb
Wadd = 518.356
ft
3 lb
Weff = 1.388 × 10
ft
Pe := ρsw ⋅ g ⋅ d
External Pressure
Pe = 36.454 psi
Pressure difference ΔP := Po − Pe
ΔP = 307.296 psi
2.2 WAVE INDUCED VELOCITY CALC.
L = 317.561 ft
Wave Length
d 1 shallow approximation
λ ( L , d) := Tp ⋅ g ⋅ d if <
L 25
2
g ⋅ Tp d 1 deep approximation
if >
2 ⋅π L 2
2
g ⋅ Tp ⎛ 2π ⋅ d ⎞ otherwise intermediate depth approximation
tanh ⎜ ⎟
2 ⋅π ⎝ L ⎠
Given
λ ( L , d) = L
L ( d) := Find ( L)
2π
k ( d) :=
L ( d)
⎡⎢ Hs g d 1 ⎤⎥ shallow
u ( d , θ ) := ⋅ if < ⋅ cos ( θ )
⎢ 2 d L ( d) 25 ⎥ approximation
⎢ π ⋅ Hs k( d) ⋅ ( Dtot−d) d 1 ⎥
⎢ ⋅e if > ⎥ deep
⎢ Tp L ( d) 2 ⎥ approximation
⎢ Hs g ⋅ Tp cosh ( k ( d) ⋅ Dtot) ⎥
⎢ 2 ⋅ ⋅ otherwise ⎥ intermediate depth
⎣ L ( d) cosh ( k ( d) ⋅ d) ⎦ approximation
uw := max ( u ( d , θ ) )
Reduction factor due to wave directionality Rwave := sin ( φwave) Rwave = 1
ft
Uw ( d , θ ) = 1.202
sec
⎡ Hs ⋅ π g d 1 ⎤ shallow
Aw ( d , θ ) := ⎢ ⋅ if < ⎥ ⋅ sin ( θ ) approximation
⎢ Tp d L ( d) 25
⎥
⎢ 2 ⎥
⎢ 2 ⋅ Hs ⋅ ⎛⎜ π ⎞⎟ ek( d) ⋅ ( Dtot−d) if d > 1
deep
⎥ approximation
⎢ ⎝ Tp ⎠ L ( d) 2 ⎥
⎢ g ⋅ π cosh [ k ( d) ⋅ ( Dtot − d) ] ⎥ intermediate depth
⎢ Hs ⋅ ⋅ otherwise ⎥ approximation
⎣ L ( d) cosh ( k ( d) ⋅ d) ⎦
Aw := max ( Aw ( d , θ ) )
Significant wave acceleration perpendicular to pipe
Aw ( d , θ ) := Aw ⋅ Rwave
ft
Aw ( d , θ ) = 2.61
2
sec
2 ⋅ Weff ⋅ δ
Stability Number Ks := Ks = 0.53
2
ρsw ⋅ Dtot
Note :
con1 := "In-line Oscillation"
con2 := "Cross-flow Oscillation"
5
Re ( tcc) = 7.505 × 10
9.87 EI Vr
Critical Pipespan Lcr := ⋅ ⋅ Dtot ⋅
2π Weff − Wadd Uw ( d , θ ) + Ur
Lcr = 159.067 ft
Yield Requirement
i := 0 .. 100 Li := i ⋅ m
⎡ ⎛ Wst ⋅ L2 ⋅ 0.5Ds ⋅ g ⎞ ⎡ ( Po − Pe) ⋅ Ds⎤ ⎤
σx ( L) := ⎢ ⎜ ⎟+⎢ ⎥⎥
Longitudinal Stress ⎣⎝ 10 ⋅ I ⎠ ⎣ 4 ⋅ ts ⎦⎦
5 .10
8
4 .10
8
3 .10
8
σx( L )
2 .10
8
1 .10
8
0
0 20 40 60 80 100
L
4
Limiting longitudinal pressure σxa ( L) := 0.9 ⋅ SMYS σxa ( L) = 4.68 × 10 psi
Lcrit ( L) := L ← 1m
while σx ( L − 1m) < σxa ( L)
Lcrit ← L − 1m
L ← L + 1m
Lcrit
Lcrit ( L) = 282.152 ft
( Po − Pe) ⋅ Ds
σh :=
Hoop stress 2 ⋅ ts
3
σh = 7.867 × 10 psi
2 2
σe := σx ( Lcrit ( L) ) + σh
Von mises
4
σe = 4.686 × 10 psi
4
0.9 ⋅ SMYS = 4.68 × 10 psi
3. SUMMARY RESULT