summaryrefslogtreecommitdiff
path: root/3588/CH9/EX9.5/EX9_5.sce
diff options
context:
space:
mode:
Diffstat (limited to '3588/CH9/EX9.5/EX9_5.sce')
-rw-r--r--3588/CH9/EX9.5/EX9_5.sce24
1 files changed, 24 insertions, 0 deletions
diff --git a/3588/CH9/EX9.5/EX9_5.sce b/3588/CH9/EX9.5/EX9_5.sce
new file mode 100644
index 000000000..ada6a46e9
--- /dev/null
+++ b/3588/CH9/EX9.5/EX9_5.sce
@@ -0,0 +1,24 @@
+//Clearing console
+clc
+clear
+
+//Intializing Variables
+P = 10
+r = 3
+z0 = 0
+z1 = 1
+
+//As we have pressure on one face only and no axial pressure, we immediately observe that
+fr2 = 0
+fz1 = 0
+fz2 = 0
+fz3 = 0
+
+//Calculating fr1 and fr3
+fr1 = 2*%pi*P*r*[integrate('z','z',z0,z1)]
+fr3 = 2*%pi*P*r*[integrate('1-z','z',z0,z1)]
+
+printf('\nResults\n')
+printf('\nNodal Forces')
+printf('\nfr1 =%flb fr2 =%flb fr3 =%flb',fr1,fr2,fr3)
+printf('\nfz1 =%flb fz2 =%flb fz3 =%flb',fz1,fz2,fz3)