# Example 15.1     
BEGIN      
 c1  =   4.0 in 
 P1  =   6.0 /in 
 f1  =   20.0 deg 
 wp_over_wg   =   -3.0  
 rg_over_rp = -wp_over_wg   
 rp   =   1.00 in 
 rg = rp  * rg_over_rp ; in 
 err = rg + rp - c1 ; in 
END
solve(err, rp)
 Np = 2 * rp *  P1   
 Ng = 2 * rg * P1  
  rbp = rp * cos(f1)  ; in 
  rbg = rg * cos(f1)  ; in 
 ra = sqrt(rbg^2 + c1^2 * sin(f1)^2)  ; in 
 ap   =   0.060 in 
 ag   =   0.290 in 
 pb = 2 * pi() * rbg / Ng  ; in 
 rap   =   1.290 in 
 rbp = rp * cos(f1)  ; in 
 rag   =   3.060 in 
 rbg = rg * cos(f1)  ; in 
 CR = ( sqrt(rap^2 - rbp^2) + sqrt(rag^2 - rbg^2) - c1 *  sin(f1) ) / pb   

# Example 15.2     
 Na   =   12.0 
 P1  =   3.0 /in
 da = Na / P1  ; in 
 na   =   600.0 rpm
 v = na * da / 2 ; ft/min 
 Wdot  =   25.0 hp
 Ft = Wdot / v ; lbf 
 f  =   20.0 deg
 Fr = Ft * tan(f)  ; lbf 
 Vcb = Ft  ; lbf 
 Hcb = Vcb * tan(f)  ; lbf 
 Fg = Vcb + Hcb  ; lbf 
 F1 = sqrt(2)* Fg  ; lbf 

# Example 15.3     
BEGIN     
 Sug   =   290.0 ksi
 Spng = Sug / 4 ; ksi 
 Sup   =   330.0 ksi
 Spnp = Sup / 4 ; ksi 
 CL   =   1.0 
 CG   =   1.0 
 CSg   =   0.70 
 CSp   =   0.68 
 kr   =   0.814 
 kt   =   1.0 
 kms   =   1.4 
 Sng = CL * CG * CSg * kr * kt * kms * Spng  ; ksi 
 Snp = CL * CG * CSp * kr * kt * kms * Spnp  ; ksi 
 P1  =   10.0 /in 
 b  =   1.25 in 
 Jg   =   0.280  
 Jp   =   0.235  
 Ng   =   18.000  
 w  =   1720.000 rpm 
 v = (Ng / P1) * w / 2 ; ft/min 
 Kv   =   1.68  
 Ko   =   1.25  
 Km   =   1.60  
 Ft   =   602.44 lbf 
 sg = Ft * P1 * Kv * Ko * Km / b / Jg   ; ksi 
 err1 = sg - Sng   ; ksi 
 HP = Ft * v ; hp  
 Ft   =   558.93 lbf 
 sp = Ft * P1 * Kv * Ko * Km / b / Jp   ; ksi
 er2 = sp - Snp   ; ksi 
 HP = Ft * v ; hp
END      

# Example 15.4      
BEGIN      
 Bhn  =   330.0  
 C0   =   0.4 ksi 
 C1   =   10.0 ksi 
 Sfe = C0 * Bhn - C1   ; ksi 
 Spfe  =   114.0 ksi 
 L1  =   10000.0 hr 
 w  =   1720.0 rpm 
 N1 = L1 * w  ; rev
 CLi  =   0.8  
 CR  =   1.0  
 SH = Sfe * CLi * CR  ; ksi 
 Cp  =   2300.0 psi^(1/2) 
 b  =   1.3 in 
 dp  =   1.8 in 
 Kv   =   1.68  
 Ko   =   1.25  
 Km   =   1.60  
 f  =   20.0 deg 
 R   =   2.00  # dg / dp 
 I = 0.5 * sin(f) * cos(f) * R / (1 + R)  
 Ft   =   129.18 lbf 
 sH = Cp * sqrt(Ft * Kv * Ko * Km / b / dp/  I)  ; ksi 
 err = sH - SH   ; ksi 
 v  =   810.53 ft/min 
 HP = Ft * v ; hp 
END 

# Example 15.5     
BEGIN     
 W1  =   3600.0 rpm
 L1  =   10000.0 hr
 N1 = L1 * W1  ; rev 
 Nfullpower = N1 / 10  ; rev 
# Case 1     
 P1  =   7.22 /in
 Np   =   18.0 
 dp = Np / P1 ; in 
 v = W1 * dp / 2  ; ft/min 
 Kv   =   1.40 
 Ko   =   1.00 
 Km   =   1.30 
 HP  =   100.00 hp
 Ft = HP / v ; lbf 
 f  =   20.0 deg 
 R =    4.00  # dg / dp
 I = 0.5 * sin(f) * cos(f) * R / (1 + R)  
 SH =   165.00 ksi # Sfe * CLi * CR  
 Cp  =   2300.0 psi^/2 
 b = 14 / P1 ; in 
 FOS  =   1.25  
 sH = Cp * sqrt(FOS * Ft * Kv * Ko * Km / b / dp / I)  ; ksi 
 err1 = sH - SH   ; ksi 
# Case 2      
 P2  =   7.00 /in 
 Np   =   18.0  
 dp = Np / P2 ; in 
 v = W1 * dp / 2  ; ft/min
 Kv   =   1.50  
 Ko   =   1.00  
 Km   =   1.30 
 HP  =   100.00 hp
 Ft = HP / v ; lbf 
 f  =   20.0 deg
 R =    4.00  # dg / dp
 I = 0.5 * sin(f) * cos(f) * R / (1 + R)  
 SH =   165.00 ksi # Sfe * CLi * CR  
 Cp  =   2300.0 psi^/2
 b = 14 / P2 ; in 
 FOS  =   1.25 
 sH = Cp * sqrt(FOS * Ft * Kv * Ko * Km / b / dp / I)  ; ksi 
 J2  =   0.235 
 s2 = Ft * P2 * Kv * Ko * Km / b / J2 ; ksi 
 CL   =   1.0 
 CG   =   1.0 
 CS   =   1.00 
 kr   =   0.897 
 kt   =   1.0 
 kms   =   1.4 
 Spn = s2 / (CL * CG * CS * kr * kt * kms)  ; ksi 
END     

# Homework Problems:

TBD