TRANSIENT ANALYSIS OF PENSTOCK

The structural design of penstock mainly refers to the thickness calculation for the desired hoop stress in case of exposed pipe and hoop and external stresses for buried pipe.

The thickness of pipe shall be checked by either Indian design code of IS 11639 (Part 2): 1995 or other similar codes.

Allowable stress

a. According to ASME code, a penstock may be designed under the following conditions

(i) Normal condition- Maximum static head plus pressure rise due to normal operation
Allowable stress = ultimate tensile strength/3 or 2/3 (0.5) * minimum yield stress, whichever is lesser.

(ii) Intermittent Condition-Condition during filling & draining the penstock
Allowable stress = Ultimate tensile strength/2.25 (2.5) or 0.8 (2/3) * yield strength, whichever is smaller

(iii) Emergency condition- Gate closure
Allowable stress = Ultimate tensile strength /1.5 (2) or 0.8* yield strength, whichever is smaller

(iv) Exceptional condition
Allowable stress < Minimum yield stress.

b. According to Indian Standard (IS 11639: 1995)

(i) Normal condition- Maximum static head plus pressure rise due to normal operation
Allowable stress = ultimate tensile strength/3 or 0.6* minimum yield stress, whichever is lesser.

(ii) Intermittent Condition-Condition during filling & draining the penstock
Allowable stress = 0.4* Ultimate tensile strength or 2/3 of yield strength, whichever is smaller

(iii) Emergency condition- Gate closure
Allowable stress = 2/3* Ultimate tensile strength or 0.9* yield strength, whichever is smaller

(iv) Exceptional condition
Allowable stress < Minimum yield stress

c. According to American Water Works Association (AWWA)

(i) Allowable working pressure, allowable stress = 50% of Yield strength

(ii) Allowable transient, allowable stress = 75% of Yield strength

Loading condition

As per IS 11639 (Part 2): 1995

a) Normal condition: Maximum static head plus pressure rise due to normal operation or head at transient maximum surge whichever is higher.

b) Intermittent Condition: Condition during filling & draining the penstocks and maximum surge in combination with pressure rise during normal operation.

c) Emergency condition- Condition includes partial gate closure in critical time of penstock (2L/a seconds) at maximum rate, and the cushioning stroke being inoperative in one unit

d) Exceptional condition-Condition includes slam shut, malfunctioning of control equipment in the most adverse manner resulting in odd situation of extreme loading. This should not be taken as a design criterion.

WATER HAMMER

Propagation of wave is given by transient formula

\(\alpha=\frac{1}{\sqrt\frac{γ}{g} (\frac{1}{k}+\frac{1}{E}×\frac{D}{t}) }\)

Where,

\(γ\) = Specific weight of water = 9.8 kN/m3

\(k\) = Bulk modulus of water = 2.04 x 106 kN/m2

\(E\) = Elastic material of pipe material

\(D\) = Diameter of pipe

\(g\) = Acceleration due to gravity = 9.8 m2/s

As per IS 11639 (Part 2): 1995, the pressure wave velocity in a steel penstock carrying water may be computed as

\(\alpha=\frac{1}{1+\frac{D}{100t}}\)

Where,

\(D\) = Diameter of pipe in m

\(t\) = Thickness of pipe in mm

Thickness of pipe to calculate rise in level.

Allievi Method

Design process

1) Calculate the individual wave velocities for different length and pipe diameter from Equations above

2) Calculate average wave velocity from following formula

\(\alpha_m=\frac{ΣL}{Σ(\frac{L}{\alpha}})\)

3) Calculate average flow velocity

\(V_m=\frac{Σ(LV)}{ΣL}\)

4) Calculate travelling time of wave

\(T=\frac{2ΣL}{\alpha_m}\)

5) Calculate increase of height of water due to wave

\(\frac{h}{H}=(\frac{0.75}{\theta\sqrt(\theta)}+1.25)×n\)

Where,

\(\rho=\frac{\alpha_m V_m}{2gH}\)


\(\theta=\frac{\alpha_m T}{2L}\)


\(n=\frac{\rho}{\theta}\)

Where,

\(\rho\) = Pipe coefficient of Allievi

\(\theta\) = Closing time of closure

\(h\) = Rise of water level due to hammer

\(H\) = Head at that point (Closure) = Net head = Total Head-Head loss

\(V_m\) = Average velocity of pipe

\(T\) = Closing time of closure

\(\alpha_m\) = Average wave velocity in pipes

6) Calculation of surge at each point (Bends, reducer, valve etc.)

\(h_i=\frac{ΣL_i}{ΣL}×h\)

Surge level at a point = Normal water level + hi

Surge height at a point = Surge level – Pipe center line level

7) Calculation of thickness

i) Minimum thickness of pipe for handling

\(t=\frac{D+800}{400}\)

Where,

\(t\) = Thickness in mm

\(D\) = Diameter of pipe in mm

As per Indian standard (Clause 7.2 of IS 11639: 1995)

\(t_0=\frac{R+0.25}{200}\)

Where,

\(R\) = Radius of pie in m

ii) Assume pipe thickness

iii) Calculation of circumferential stress

\(\sigma_y=\frac{P(D+2\epsilon_1)}{2(t-\epsilon_1-\epsilon_2)η}\)

Where,

\(P\) = Water pressure =\( γ H\)

\(D\) = Diameter of pipe

\(\epsilon_1\) = Internal corrosion thickness

\(\epsilon_2\) = External corrosion thickness

\(t\) = Thickness of pipe

\(η\) = Welding efficiency = 0.9 to 1.0

iv) Calculation of axial stress

\(\sigma_x=ν\sigma_y\)

Where,

\(ν\) = Poisson's ratio

v) Calculation of combined circumferential stress

\(\sigma_g=\sqrt(σ_x^2+σ_y^2-σ_x σ_y ) \)

If \(\sigma_a>sigma_y\) and \(\sigma_a>\sigma_g\), The assumed thickness of pipe is safe.

Design for external pressure (As per IS 11639 Part 2: 1995)

A) Vaughan’s formula

\(\frac{13K^2}{4E'}P_{cr}^2+2P_{cr}(1+3K\frac{Y_o}{R}-\frac{F_yK}{2E'})-(\frac{4F_y}{K}-\frac{(F_y)^2}{E'})=0\)

Where,

\(K\) = ratio of pipe diameter to plate thickness

\(E'=\frac{E_s}{1-(\mu_s)^2}\)

\(E_s\) = modulus of elasticity of steel on N/m2

\(mu_s\) = Poisson’s ratio of steel

\(P_{cr}\) = critical external pressure at buckling in N/m2

\(Y_o\) = initial gap between steel lining and concrete in m

\(R\) = radius of steel liner in m

\(F_y\) = yield point stress in steel in N/m2


B) Amstutz’s formula

\((\frac{f_n}{E'}+\frac{Y_o}{R})[1+\frac{3K^2f_n}{E'}]^2=1.68K\frac{(f_y'-f_n)}{E'}×[1-\frac{K}{4}×\frac{f_y'-f_n}{E'}]\)

And

\(1-\frac{P_{cr}K}{2f_n}=0.175\frac{K}{E'}(f_y'-f_n)\)

Where,

\(K=\frac{2R}{t}\)

\(E'=\frac{E_s}{1-mu_s^2}\)

\(f_y'=\frac{f_y}{\sqrt(1-μ_s+μ_s^2) }\)

\(f_y\) = yield stress in steel in N/m2

\(E_s\) = modulus of elasticity of steel in N/m2

\(mu_s\) = Poisson’s ratio of steel

\(R\) = Radius of liner in m

\(Y_o\) = initial gap between steel lining and concrete in m

\(f_n\) = allowable stress in steel in N/m2

\(P_{cr}\) = critical external pressure at buckling in N/m2

\(t\) = thickness of steel liner in m.

The maximum external pressure on the penstock should not exceed two-thirds times the critical external pressure calculated according to Vaughan’s equation or Amstutz’s equation. If the maximum pressure exceeds the value equal to two-thirds times the critical external pressure for unstiffened shell, stiffener should be provided to prevent buckling.