Mathcad - Example 6
Mathcad - Example 6
Darren Ting
Q_s 34690lbf
M 1312600in lbf
L_l 30ft
Depth of Head
L_h 15in
L_a 4ft
E_c 0.85
V_l 0lbf
M_t 0lbf in
F_axial 0lbf
Assuming
RSF_a 0.9
0
Refer to 6.2 of API 579 for and limitations of procedure. All these limitations are assumed
for this analysis.
April 21 2013
1
2
3
4
5
6
7
8
9
Pit_Couple 0 P_k
11
12
13
14
15
16
17
18
19
0.6
1.8
.9
1.2
1.2
2.2
.5
1
.8
d_jk .3 in
.8
.5
1.6
.8
.8
.2
1.2
.6
.5
3.5
4.2
2.7
2.1
4.6
3.1
2.9
3.1
2.6
2.2 in
1.8
2.5
3.8
1.9
1.8
1
2.5
1.5
1.3
10
15
22
30
5
15
20
45
60
_k 0
10
20
35
90
0
22
45
67
90
.4
.65
.75
.6
.5
.45
.6
.75
.2
w_jk .75 in
.5
.7
.75
.5
.5
.7
.4
.7
.7
0.5
1.6
0.9
1
.7
1.1
.8
.5
1.3
d_ik .4 in
1.5
.6
2.4
.4
1
.6
.9
.6
.8
.5
.6
.5
.7
.6
.5
.65
.4
.5
w_ik .55 in
.4
.75
.5
.25
.7
.75
.3
.5
.4
1
2
0.45
0.625
W_avg_1
( w_ik w_jk)
0.625
0.65
0.55
0.475
0.625
0.575
0.35
W_avg_1
0.65
0.45
0.725
0.625
0.375
0.6
0.725
0.35
0.6
0.55
in
R_c D_o
R_m
30.6025 in
2
From Annex A .7.3 API 579(2009)-Note that 2007 equation is slighty different
2( R_m) 2 L_h2
( L_l)
( 3 Q_s L_l)
4 L_a 0.3633 in
T_sl
2
( 4 L_h)
L_l
SA E_l ( R_m)
1
3 L_l
Check
L_l 30 ft
R_m 30.6025 in
L_h 15 in
_1
T_c
E_c
_2
R_c
P 0.4
1.2861 104 psi
T_c T_sl
2E_l
MAWP_c
MAWP_l
( SA E_c T_c)
R_c 0.6T_c
506.216 psi
680.3423 psi
Governing MAWP
MAWP min( MAWP_c MAWP_l ) 506.216 psi
Calculate avg1
0.55
0.8429
1.7
0.5952
0.9
0.6667
1.1
0.4762
0.95
0.7935
1.65
0.4677
0.65
0.7759
0.75
D_avg_k
( d_ik d_jk)
2
0.7581
1.05
0.35
in
_avg1
( P_k D_avg_k)
P_k
1.15
0.5962
0.8409
0.3611
0.55
0.78
0.4737
0.6
0.6842
0.9
0.5
0.4
0.6
1.05
0.58
0.6
0.6
0.65
0.5
Calculate _21
Calculate _11
0
0
_11
_1
_avg1
2.0507759104
2.90389868104
2.59276668104
3.62987335104
2.17839757104
3.69543756104
2.22785877104
2.2801636104 psi
2.899438104
2.05552673104
10
4.78664617104
11
2.21603989104
12
3.64907903104
13
2.52628548104
14
3.45702223104
15
...
_21
_2
_avg1
1.52590197104
2.16067719104
1.92917606104
2.70084649104
1.62086025104
2.74963025104
6
7
1.65766239104
2.15735818104
1.52943688104
10
3.56155581104
11
1.64886843104
12
2.71513668104
13
1.87971001104
14
2.57223475104
15
...
1.69658037104 psi
2
4
2
3 ( sin( 2_k) ) ( _11 _21) ( sin( _k) ) 4 ( sin( 2_k) ) 2 ( _21) 2
_1 _11 ( cos( _k) ) ( sin( 2_k) )
Note typo in api 2007 example for _1. There should not be
numerical values in given formula
0
0
4.17324402108
8.28405144108
6.46071703108
1.2191117109
4.73635557108
1.34156465109
4.80493929108
_1 7
Check
5.9020898108
4.22519015108
10
2.27352649109
11
4.75409373108
12
1.19379849109
13
3.53330972108
14
1.19510027109
15
...
4.29419254108 psi
5.4908 10 psi
Note -unable to directly use _avg1 max( _11 _21 _11 _21 ) due to inner product being taken.
require element by element operation. To combine all three vectors and filter for maximum value as shown
below:
2.0508104
1.5259104
5.2487103
2.9039104
2.1607104
7.4322103
2.5928104
1.9292104
6.6359103
3.6299104
2.7008104
9.2903103
2.1784104
1.6209104
5.5754103
3.6954104
2.7496104
9.4581103
6
augment( _11 _21 _11 _21) 7
2.2279104
1.6577104
5.702103
2.2802104
1.6966104
5.8358103 psi
2.8994104
2.1574104
7.4208103
2.0555104
1.5294104
5.2609103
10
4.7866104
3.5616104
1.2251104
11
2.216104
1.6489104
5.6717103
12
3.6491104
2.7151104
9.3394103
13
2.5263104
1.8797104
6.4658103
14
3.457104
2.5722104
8.8479103
15
2.8809104
2.1435104
...
Combine three vectors into matrix for _11 _21 _11 _21 and choose maximum value and output as vector
i 0 rows( augment( _11 _21 _11 _21) ) 1
j 0 cols( augment( _11 _21 _11 _21) ) 1
Max_val max( submatrix( augment( _11 _21 _11 _21) i i 0 cols( augment( _11 _21 _11 _21) ) 1 ) )
i
0
0
2.0508104
2.9039104
2.5928104
3.6299104
2.1784104
3.6954104
6
Max_val 7
2.2279104
2.8994104
2.0555104
10
4.7866104
11
2.216104
12
3.6491104
13
2.5263104
14
3.457104
15
...
2.2802104 psi
Calculate _1
0
0
1.7285104
1.7285104
1.7285104
1.7285104
1.7285104
1.7285104
_1 ( _avg1 Max_val) 7
1.7285104
1.7285104
1.7285104
10
1.7285104
11
1.7285104
12
1.7285104
13
1.7285104
14
1.7285104
15
...
1.7285104 psi
Calculate E_avg1
0
E_avg_1
0.8461
0.8461
0.6006
0.6006
0.68
0.68
0.4951
0.4951
0.7942
0.7942
0.4719
0.4719
0.7885
0.7885
0.8341
0.7115
0.8409
_1
7
0.5
( _1)
8
Check
_1
0.8341
0.7115
( _1)
0.5
0.8409
10
0.3625
10
0.3625
11
0.7928
11
0.7928
12
0.5003
12
0.5003
13
0.9196
13
0.9196
14
0.5
14
0.5
15
0.612
15
...
16
0.6382
17
0.7458
18
0.672
Calculate RSF_1
0.9337
0.7611
0.8086
0.6859
0.8917
0.76
0.8735
0.9087
RSF_1 1
[ W_avg_1 ( 1 E_avg_1 ) ]
T_c
0.9034
0.901
0.7255
0.8562
0.7011
0.9711
0.7129
0.7308
0.8788
0.8541
0.8274
Check
W_avg_1 0.45 in
0
T_c 1.045 in
E_avg_1 0.8461
0
RSF_1i
RSF_pit
0.8256
length ( RSF_1)
Check
length ( RSF_1) 19
RSF_1i 15.6857
i
Note that PSF_pit value is different from example due to one entry being different
RSF_pit 0.8256
RSF_a 0.9
Check longitudinal stress condition
Statement 3
result_3 if ( RSF_pit RSF_a true false)
result_3 "true"
Hence determine MAWP_r using equation in Part2 paragraph 2.4.2.2
MAWP_r MAWP
RSF_pit
464.3475 psi
RSF_a
Calculate R_t
0
R_t
0.5694
0.4737
0.5694
0.378
0.4737
0.5694
6
T_c FCA w_ik
7
T_c
0.4258
0.5694
0.5215
10
0.6651
11
0.3301
12
0.5694
13
0.8086
14
0.378
15
...
0.6651
Calculate Q
From table 4.5
Q 1.123
( 1 R_t) 2
R_t
RSF_a
0.5
0.6869
0.5439
0.6869
0.435
0.5439
0.6869
0.4865
0.9028
0.6869
0.6095
10
0.9028
Check
11
0.3878
R_t 0.5694
12
0.6869
RSF_a 0.9
13
1.7942
14
0.435
15
...
Determine if
5.4388
4.3067
5.4388
3.4447
4.3067
5.4388
3.8523
Q ( D T_c) 0.5
7
7.1489
0.5
1.6
0.9
0.7
1.1
6
d_ik 7
0.8
Check
D 60 in
T_c 1.045 in
Q 0.6869
0
in
d values
0.5
5.4388
1.3
4.8264
0.4
10
7.1489
10
1.5
11
3.0704
11
0.6
12
5.4388
12
2.4
13
14.2069
13
0.4
14
3.4447
14
15
...
15
...
in
d_ik Q ( D T_c)
0.5
10
11
12
13
14
15
if mind_ik Q ( D T_c)
16
17
18
0.5
T_c 1.045 in
RSF_pit 0.8256
MAWP_r 464.3475 psi
RSF_a 0.9
T_sl 0.3633 in
Weight Case
4
M 1.3126 10 in lbf
No Thermal Load assumed
Calculate B
B min 1
RSF_pit
RSF_a
0.9173
Calculate T_eq
T_eq B T_c 0.9586 in
M 1.3126 10 in lbf
MAWP_r
_cm
R_c
T_eq
0.6
RSF_pit cos( )
Check
D_o 62.25 in
A_m
T_eq 0.9586 in
D_f 60.3329 in
D_o 62.25 in
D_o
I_x
31.125 in
2
3 2
A_a D_f 2.8589 10 in
4
A_t
Calculate
M_t
a A_t T_eq
V_l
A_m
0 psi
Check
2
A_m 184.5751 in
T_eq 0.9586 in
A_t 2.9505 10 in
Longitudinal Membrane Stress
No applied section axial force for weight or weight plus thermal load
Calculate
Tensile
_lmt
1
A_a MAWP_r F_axial M a
9.0159 103 psi
E_c
cos
(
)
A_m
A_m
I_x
A_m 1.2818 ft
Check
2
3
A_a 2.8589 10 in
2
A_m 184.5751 in
F_axial 0 lbf
A_m 184.5751 in
6
M 1.3126 10 lbf in
4
I_x 8.6694 10 in
a 31.125 in
Compressive
_lmc
1
A_a MAWP_r F_axial M a
7.9071 103 psi
I_x
E_c
cos
(
)
A_m
A_m
Weight case-tensile
0.5
1.5578 10 psi
Check
0
3
No Thermal Loads
0.5
1.5616 10 psi
Check
4
SA 1.75 10 psi
1.9444 10 psi
RSF_a
SA
RSF_a 0.9
Statement 4
true false
RSF_a
SA
result_4 "true"
Conclusion from API 579 2007
If maximum longitudinal stress in step 7 is compressive, this stress should be less than or equal to allowable
compressive stress computed using paragrpah A.4.4 or allowable tensile stress, whichever smaller.When using
this methodology to establish an allowable compressive stress, an average thickness representative of the region
of pitting damage in the compressive stress zone should be used in the calculations.
The maximum longitudinal stress in STEP 7 is NOT compressive.
SUMMARY
MAWP_r 464.3475 psi
The longitudinal stress is acceptable. Equipment fir for service for MAWP_r calculated.