From 7bc77cb1ed33745c720952c92b3b2747c5cbf2df Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Sat, 3 Feb 2018 11:01:52 +0530
Subject: Added new code

---
 3863/CH10/EX10.19/Ex10_19.sce | 31 +++++++++++++++++++++++++++++++
 1 file changed, 31 insertions(+)
 create mode 100644 3863/CH10/EX10.19/Ex10_19.sce

(limited to '3863/CH10/EX10.19')

diff --git a/3863/CH10/EX10.19/Ex10_19.sce b/3863/CH10/EX10.19/Ex10_19.sce
new file mode 100644
index 000000000..2c4b532e3
--- /dev/null
+++ b/3863/CH10/EX10.19/Ex10_19.sce
@@ -0,0 +1,31 @@
+clear
+//
+//
+//Given
+//Variable declaration
+h=20              //height in m
+D=4               //External diameter in m
+d=2               //Internal diameter in m
+p=1               //Horizontal wind pressure in kN/sq.m
+w=22              //specific weight in kN/m^3
+K=2/3             //Co-efficient of wind resistance
+
+//Calculation
+A1=(%pi/4)*(D**2-d**2)           //Area of cross-section
+W=w*A1*h                             //Weigth of the chimney in kN
+sigma0=W/A1                          //Direct stress in kN/sq.mm
+A2=D*h                               //Projected area of the surface exposed to wind in sq.m
+F=K*p*A2                             //Wind Force in kN
+M=F*h/2                              //Bending moment in kNm
+I=(%pi/64)*(D**4-d**4)           //Moment of inertia
+y=D/2                                //Distance between C.G. of the base section and extreme edge of the base
+Z=I/y                                //Section modulus
+sigmab=M/Z                           //Bending stress
+sigma_max=(sigma0+sigmab)     //Maximum stress in kN/sq.m
+
+sigma_min=(sigma0-sigmab)     //Minimum stress in kN/sq.m
+
+
+//Result
+printf("\n Maximum stress = %0.3f  kN/m^2",sigma_max)
+printf("\n Minimum stress = %0.3f  kN/m^2",sigma_min)
-- 
cgit