Calcvessel
Calcvessel
Calcvessel
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Introduction:
This Mathcad template is made for the purpose of aiding, automating and simplifiying the calculations of various parameters
concerning unfired pressure vessels, according to the standard NF EN 13445-3. Please note that the worksheet is intended as a
supplementary tool to the standard, and should therefore always be used alongside the standard, and never alone. To ensure
minimal risk of error, it is recommended (although not necessary) to open a new version of the worksheet for each new calculation.
When printing the worksheet, it is recommended to save as pdf and print in A4 size. Due to the nature of hidden areas in Mathcad, it
may be necessary to add a few lines to ensure no text or calcultations gets split in two between pages.
To use the worksheet, all relevant input parameters should be entered first. The worksheet uses two sets of input parameters:
Global input parameters and local input parameters. The global input variables are values concerning the materials of the vessel and
the inteded pressure, and are found under this introduction. The local input variables concern the specific geometry of the vessel.
Each section contains its own relevant input parameters at the top of its hidden area. Every input parameter is coloured blue . After
all relevant input parameters have been entered, the worksheet will automatically calculate the output values further down. The
most relevant output parameters are green . An orange value might be of interest to the reader, but it is not the end result.
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
The value of fdn above, where n is the desired steel designation, should be entered here as f f ≔ 154 [MPa]
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Calculation values:
Analysis thickness:
1
ea ≔ ―⋅ ⎛⎝De - Di⎞⎠ = 1.82
2
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Output parameters:
P ⋅ Di
e ≔ ―――― = 0.329 7.4-1
2 f⋅z-P
P ⋅ De
OR e ≔ ―――― = 0.367 7.4-2
2 f⋅z+P
2 f ⋅ z ⋅ ea
Pmax ≔ ―――= 20.901 7.4-3
Dm
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Conditions of applicability:
‖ e || ‖ e ||
status_R1 ≔ ‖ if ― ≤ 0.16 || if only the internal diameter is known: status_R1i ≔ ‖ if ――― ≤ 0.16| |
‖ De || ‖ Di + 2 e ||
‖ ‖ | ‖ ‖ |
‖ ‖ return “OK” | | ‖ ‖ return “OK” | |
‖
‖ “Not Valid” | ‖
‖ “Not Valid” |
Requirements Results
e
R1: ― ≤ 0.16 status_R1 = “OK”
De
status_R1i = “OK”
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Calculation values:
Analysis thickness:
1
ea ≔ ―⋅ ⎛⎝De - Di⎞⎠ = 0.445
2
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Output parameters:
P ⋅ Di
e ≔ ―――― = 0.294 7.4-4
4 f⋅z-P
P ⋅ De
OR e ≔ ―――― = 0.296 7.4-5
4 f⋅z+P
4 f ⋅ z ⋅ ea
Pmax ≔ ―――= 6.032 7.4-6
Dm
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Conditions of applicability:
‖ e || ‖ e ||
status_R1 ≔ ‖ if ― ≤ 0.16 || if only the internal diameter is known: status_R1i ≔ ‖ if ――― ≤ 0.16| |
‖ De || ‖ Di + 2 e ||
‖ ‖ | ‖ ‖ |
‖ ‖ return “OK” | | ‖ ‖ return “OK” | |
‖
‖ “Not Valid” | ‖
‖ “Not Valid” |
Requirements Results
e
R1: ― ≤ 0.16 status_R1 = “OK” status_R1i = “OK”
De
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
NF EN 13445-3 Worksheet template prepared by Jelena Jovanovic and Eskil Vik
Output parameters:
P ⋅ Di 1
e ≔ ―――― ⋅ ――― = 0.304 7.6-2
2 f ⋅ z - P cos ((α))
P ⋅ De 1
OR e ≔ ―――― ⋅ ――― = 0.326 7.6-3
2 f ⋅ z + P cos ((α))
NOTE: both the input and output parameters in this section only considers one specific point on the cone, so
multiple calculations may be required
Conditions of applicability:
‖ e ( ) ||
con.a ⋅ cos (α)
status_R1 ≔ ‖ if α ≤ 75 || status_R2 ≔ ‖ if ――――― > 0.001| |
‖ ‖ | ‖ Dc ||
‖ ‖ return “OK” | | ‖ ‖ |
‖ “Not Valid” | ‖ ‖ return “OK” ||
‖ | ‖ “Not Valid” |
‖ | ‖
‖ |
Requirements Results
Input parameters:
7.5.3.2 Design
r
―= 0.093
Di
⎛ Di ⎞ P
⎜0.75 + 0.2 ―⎟ ⋅ ―= 0.023
⎝ R⎠ f
The value of βa can be read from the chart above or calculated through the iterative βa ≔ [unitless]
process shown in 3.5.3.5, a)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Calculation values:
Rp0.2;T
fb1 ≔ ――― = 133.333
1.5
7.5-4
OR, for cold spun seamless austentic stainless steel:
1.6 Rp0.2;T
fb2 ≔ ―――― = 213.333 7.5-5
1.5
The correct value for fb should be entered with regard to the material used: fb ≔ 133 [MPa]
(at test conditions the value 1.5 in the fb equations shall be replaced by 1.05)
The required tickness e shall be the greatest of es, ey and eb, where:
P⋅R
es ≔ ――――― = 7.843
2 f ⋅ z - 0.5 P
7.5-1
⎛⎝0.75 R + 0.2 Di⎞⎠
ey ≔ βa ⋅ P ⋅ ―――――― =?
f
⎛ 1 ⎞
⎜―― ⎟ 7.5-2
⎝ 1.5 ⎠
⎛ P ⎛ Di ⎞ 0.825⎞
eb ≔ ⎛⎝0.75 R + 0.2 Di⎞⎠ ⋅ ⎜――― ⋅ ⎜―⎟ ⎟ = 8.286
⎜⎝ 111 fb ⎝ r ⎠ ⎟⎠ 7.5-3
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Output parameters:
NOTE 1 For stainless steel ends that are not cold spun, fb will be less than f.
NOTE 2 The 1.6 factor for cold spun ends takes account of strain hardening.
NF EN 13445-3 Worksheet template prepared by Jelena Jovanovic and Eskil Vik
‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
⎛ Di ⎞ ⎛ Di ⎞
hi ≔ R - ⎜R - ―⎟ ⋅ ⎜R + ―- 2 r⎟ = 68.046
⎝ 2 ⎠ ⎝ 2 ⎠
7.5.3.3 Rating
ea
―= 0.003
R
r
―= 0.093
Di
The value of βb can be read from the chart above or calculated through the process βb ≔ [unitless]
shown in 3.5.3.5, b)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Calculation values:
Rp0.2;T
fb1 ≔ ――― = 133.333
1.5
7.5-4
OR, for cold spun seamless austentic stainless steel:
1.6 Rp0.2;T
fb2 ≔ ―――― = 213.333 7.5-5
1.5
The correct value for fb should be entered with regard to the material used: fb ≔ 133 [MPa]
(at test conditions the value 1.5 in the fb equations shall be replaced by 1.05)
For a given geometry Pmax shall be the least of Ps, Py and Pb, where:
⎛⎝2 f ⋅ z ⋅ ea⎞⎠
Ps ≔ ―――― = 1.025 7.5-6
R + 0.5 ea
f ⋅ ea
Py ≔ ――――――― =? 7.5-7
βb ⋅ ⎛⎝0.75 R + 0.2 Di⎞⎠
⎛ ea ⎞ 1.5 ⎛ r ⎞ 0.825
Pb ≔ 111 fb ⋅ ⎜―――――⎟ ⎜―⎟ = 0.474 7.5-8
⎝ 0.75 R + 0.2 Di ⎠ ⎝ Di ⎠
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Output parameters:
7.5.3.4 Exceptions
Functions
⎛ ey ⎞ ⎛ 1 ⎞ 1
Y ⎛⎝ey⎞⎠ ≔ min ⎜―, 0.04⎟ Z ⎛⎝ey⎞⎠ ≔ log ⎜――⎟ N ((Y)) ≔ 1.006 - ―――――
⎝ Y ⎛⎝ey⎞⎠ ⎠ ⎛ 4⎞
⎝R ⎠
⎝6.2 + ((90 Y)) ⎠
Calculation quantities
r
X ≔ ―= 0.093
Di
Initial beta ( β ): β0 ≔ 1 0.4 ≤ β ≤ 1.8 (The initial value does not have any impact on the result)
⎛ ea ⎞
Y ≔ min ⎜―, 0.04⎟ = 0.003 7.5-9
⎝R ⎠
⎛1⎞
Z ≔ log ⎜―⎟ = 2.477 7.5-10
⎝Y⎠
r
X ≔ ―= 0.093 7.5-11
Di
1
N ≔ 1.006 - ――――― = 0.845 7.5-12
⎛ 4⎞
⎝6.2 + ((90 Y)) ⎠
βb ≔ ‖ if X = 0.06 | |
‖ ‖ | |
‖ ‖ β ← β0.06 | |
‖ | |
‖ if 0.06 < X < 0.1| |
‖ ‖ β ← βA |
‖ ‖ ||
‖ if X = 0.1 | |
‖ ‖β←β | |
0.1 |
‖ ‖ |
‖ if 0.1 < X < 0.2| |
‖ ‖β←β | |
‖ ‖ B | |
‖ if X = 0.2 | |
‖ ‖ | |
‖ ‖ β ← β0.2 |
‖ |
f ⋅ ea
βb = 1.12 => Py ≔ ――――――― = 0.513
βb ⋅ ⎛⎝0.75 R + 0.2 Di⎞⎠
Requirements Results
Input parameters:
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Calculation values:
At C:
P ⋅ ((a + L))
σm.C ≔ ―――― = 556 15.5.1.2-1
e
At D:
At B:
P ⋅ ⎛⎝a + l1⎞⎠
σm.B ≔ ―――― = 600 15.5.1.2-2
e
At A:
At a corner:
P
σm.BC ≔ ―⋅ ⎛⎝a + ‾‾‾‾‾‾
L 2 + l1 2 ⎞⎠ = 790.47 15.5.1.2-3
e
3
I2 ≔ I1 = 2.25
The following equations apply to calculate the bending stresses (eq. 15.5.1.2-5 - eq. 15.5.1.2-9)
⎛ l1 ⎞ 15.5.1.2-10
θ ≔ atan ⎜―⎟ = 0.828
⎝L⎠
NOTE: the value of θ varies thorughout the corner, but the value given here will give the highest bending stress
a
ϕ ≔ ―= 0.125 15.5.1.2-15
l1
L
α3 ≔ ―= 0.918 15.5.1.2-14
l1
⋅
MA ≔ P ⋅ ⎛⎝-K3⎞⎠ = ⎛⎝3.526 ⋅ 10 4 ⎞⎠ ――― 15.5.1.2-10
At C:
e ⎛
± σb.C ≔ ―― ⎛ 2 ⎞⎞
⎝2 MA + P ⋅ ⎝2 a ⋅ L - 2 a ⋅ l1 + L ⎠⎠ = ⎛⎝1.987 ⋅ 10 ⎞⎠
5
15.5.1.2-5
4 I1
At D:
⎛ e ⎞
± σb.D ≔ abs ⎜―― ⎛ 2 M + P ⋅ ⎛ 2 a ⋅ L - 2 a ⋅ l + L 2 - l 2 ⎞⎞⎟ = ?
⎝ A ⎝ 1 1 ⎠⎠
⎝ 4 I1 ⎠ 15.5.1.2-6
At A:
⎛ MA ⋅ e ⎞
± σb.A ≔ abs ⎜――⎟ = ?
⎝ 2 I1 ⎠ 15.5.1.2-7
At B:
⎛ e ⎞
± σb.B ≔ abs ⎜―― ⋅ ⎛⎝2 MA + P ⋅ L 2 ⎞⎠⎟ = ⎛⎝2.031 ⋅ 10 5 ⎞⎠
⎝ 4 I1 ⎠ 15.5.1.2-8
At a corner:
⎛ e ⎞
± σb.BC ≔ abs ⎜―― ⋅ ⎛⎝2 MA + P ⋅ ⎛⎝2 a ⋅ ⎛⎝L ⋅ cos ((θ)) - l1 ⋅ ((1 - sin ((θ))))⎞⎠ + L 2 ⎞⎠⎞⎠⎟ = ⎛⎝2.221 ⋅ 10 5 ⎞⎠ 15.5.1.2-9
⎝ 4 I1 ⎠
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Output parameters::
The maximum stress at any location is equal to the sum of the bending and membrane stresses at that location (as stated in
15.4). The maximum stress is therefore:
At C:
At D:
At A:
At B:
At the corners:
σm ≤ f ⋅ z 15.5.3-1
This gives:
σm + σb ≤ 1.5 ⋅ f ⋅ z 15.5.3-2
This gives: