summaryrefslogtreecommitdiff
path: root/3886/CH9
diff options
context:
space:
mode:
Diffstat (limited to '3886/CH9')
-rw-r--r--3886/CH9/EX9.1/9_1.sce6
-rw-r--r--3886/CH9/EX9.1/9_1.txt3
-rw-r--r--3886/CH9/EX9.10/9_10.sce13
-rw-r--r--3886/CH9/EX9.10/9_10.txt6
-rw-r--r--3886/CH9/EX9.11/9_11.sce11
-rw-r--r--3886/CH9/EX9.11/9_11.txt6
-rw-r--r--3886/CH9/EX9.12/9_12.sce9
-rw-r--r--3886/CH9/EX9.12/9_12.txt5
-rw-r--r--3886/CH9/EX9.13/9_13.sce19
-rw-r--r--3886/CH9/EX9.13/9_13.txt6
-rw-r--r--3886/CH9/EX9.14/9_14.sce13
-rw-r--r--3886/CH9/EX9.14/9_14.txt5
-rw-r--r--3886/CH9/EX9.15/9_15.sce4
-rw-r--r--3886/CH9/EX9.15/9_15.txt3
-rw-r--r--3886/CH9/EX9.16/9_16.sce7
-rw-r--r--3886/CH9/EX9.16/9_16.txt3
-rw-r--r--3886/CH9/EX9.17/9_17.sce9
-rw-r--r--3886/CH9/EX9.17/9_17.txt3
-rw-r--r--3886/CH9/EX9.18/9_18.sce4
-rw-r--r--3886/CH9/EX9.18/9_18.txt3
-rw-r--r--3886/CH9/EX9.19/9_19.sce5
-rw-r--r--3886/CH9/EX9.19/9_19.txt3
-rw-r--r--3886/CH9/EX9.2/9_2.sce11
-rw-r--r--3886/CH9/EX9.2/9_2.txt5
-rw-r--r--3886/CH9/EX9.20/9_20.sce8
-rw-r--r--3886/CH9/EX9.20/9_20.txt5
-rw-r--r--3886/CH9/EX9.3/9_3.sce16
-rw-r--r--3886/CH9/EX9.3/9_3.txt3
-rw-r--r--3886/CH9/EX9.4/9_4.sce20
-rw-r--r--3886/CH9/EX9.4/9_4.txt5
-rw-r--r--3886/CH9/EX9.5/9_5.sce21
-rw-r--r--3886/CH9/EX9.5/9_5.txt5
-rw-r--r--3886/CH9/EX9.6/9_6.sce43
-rw-r--r--3886/CH9/EX9.6/9_6.txt4
-rw-r--r--3886/CH9/EX9.7/9_7.sce8
-rw-r--r--3886/CH9/EX9.7/9_7.txt4
-rw-r--r--3886/CH9/EX9.8/9_8.sce25
-rw-r--r--3886/CH9/EX9.8/9_8.txt3
-rw-r--r--3886/CH9/EX9.9/9_9.sce13
-rw-r--r--3886/CH9/EX9.9/9_9.txt7
40 files changed, 352 insertions, 0 deletions
diff --git a/3886/CH9/EX9.1/9_1.sce b/3886/CH9/EX9.1/9_1.sce
new file mode 100644
index 000000000..81138e071
--- /dev/null
+++ b/3886/CH9/EX9.1/9_1.sce
@@ -0,0 +1,6 @@
+//Locating the centroid of T section
+//Refer fig. 9.10
+//due to symmetry the centroid lies on y axis
+//distance of centroid from top is
+ybar=(100*20*10+20*100*70)/(100*20+20*100)
+printf("Centroid of T-section is on the symmetric axis at a distance of %.2f mm from the top",ybar)
diff --git a/3886/CH9/EX9.1/9_1.txt b/3886/CH9/EX9.1/9_1.txt
new file mode 100644
index 000000000..b227287c1
--- /dev/null
+++ b/3886/CH9/EX9.1/9_1.txt
@@ -0,0 +1,3 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.1.sce', -1)
+Centroid of T-section is on the symmetric axis at a distance of 40.00 mm from the top \ No newline at end of file
diff --git a/3886/CH9/EX9.10/9_10.sce b/3886/CH9/EX9.10/9_10.sce
new file mode 100644
index 000000000..af8e232a5
--- /dev/null
+++ b/3886/CH9/EX9.10/9_10.sce
@@ -0,0 +1,13 @@
+//Moment of inertia of L-section
+//Divide the section into two rectangles A1 and A2
+A1=125*10 //mm^2
+A2=75*10 //mm^2
+A=A1+A2 //mm^2
+xbar=((1250*5)+750*(10+75/2))/A //mm
+ybar=((1250*125/2)+(750*5))/A //mm
+Ixx=((10*125^3)/12)+(1250*21.56^2)+((75/12)*10^3)+(750*39.94^2) //mm^4
+Iyy=((125*10^3)/12)+(1250*15.94^2)+((10*75^3)/12)+(750*26.56^2) //mm^4
+Izz=Ixx+Iyy //mm^4
+printf("\nIxx=%.2f mm^4\nIyy=%.2f mm^4\nIzz=%.2f mm^4",Ixx,Iyy,Izz)
+
+
diff --git a/3886/CH9/EX9.10/9_10.txt b/3886/CH9/EX9.10/9_10.txt
new file mode 100644
index 000000000..78517996d
--- /dev/null
+++ b/3886/CH9/EX9.10/9_10.txt
@@ -0,0 +1,6 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.10.sce', -1)
+
+Ixx=3411298.87 mm^4
+Iyy=1208658.87 mm^4
+Izz=4619957.73 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.11/9_11.sce b/3886/CH9/EX9.11/9_11.sce
new file mode 100644
index 000000000..7ecc83b92
--- /dev/null
+++ b/3886/CH9/EX9.11/9_11.sce
@@ -0,0 +1,11 @@
+//Moment of inertia of I section
+//Refer fig. 9.36
+A1=200*9 //mm^2
+A2=(250-9*2)*6.7 //mm^2
+A3=200*9 //mm^2
+A=A1+A2+A3 //mm^2
+Ixx=((200*9^3)/12)+(1800*120.5^2)+((6.7*232^3)/12)+(0)+((200*9^3)/12)+(1800*120.5^2) //mm^4
+Iyy=((9*200^3)/12)+((232*6.7^3)/12)+((9*200^3)/12) //mm^4
+Izz=Ixx+Iyy //mm^4
+printf("\nIxx=%.2f mm^4\nIyy=%.2f mm^4\nIzz=%.2f mm^4",Ixx,Iyy,Izz)
+//The answers vary due to round off error
diff --git a/3886/CH9/EX9.11/9_11.txt b/3886/CH9/EX9.11/9_11.txt
new file mode 100644
index 000000000..7646a75cc
--- /dev/null
+++ b/3886/CH9/EX9.11/9_11.txt
@@ -0,0 +1,6 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.11.sce', -1)
+
+Ixx=59269202.13 mm^4
+Iyy=12005814.75 mm^4
+Izz=71275016.88 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.12/9_12.sce b/3886/CH9/EX9.12/9_12.sce
new file mode 100644
index 000000000..a4fa20ad3
--- /dev/null
+++ b/3886/CH9/EX9.12/9_12.sce
@@ -0,0 +1,9 @@
+//Second moment of Inertia
+//refer fig. 9.37
+A1=100*13.5 //mm^2
+A2=(400-27)*8.1 //mm^2
+A3=100*13.5 //mm^2
+A=A1+A2+A3 //mm^2
+Ixx=((100*13.5^3)/12)+(1350*193.25^2)+((8.1*373^3)/12)+((100*13.5^3)/12)+(1350*193.25^2) //mm^4
+Iyy=((13.5*100^3)/12)+(1350*24.27^2)+((373*8.1^3)/12)+(3021.3*21.68^2)+((13.5*100^3)/12)+(1350*24.27^2)
+printf("\nIxx=%.2d mm^4\nIyy=%.2d mm^4",Ixx,Iyy)
diff --git a/3886/CH9/EX9.12/9_12.txt b/3886/CH9/EX9.12/9_12.txt
new file mode 100644
index 000000000..b4cad0683
--- /dev/null
+++ b/3886/CH9/EX9.12/9_12.txt
@@ -0,0 +1,5 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.12.sce', -1)
+
+Ixx=135903228 mm^4
+Iyy=5276986 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.13/9_13.sce b/3886/CH9/EX9.13/9_13.sce
new file mode 100644
index 000000000..ad31e5007
--- /dev/null
+++ b/3886/CH9/EX9.13/9_13.sce
@@ -0,0 +1,19 @@
+//Polar moment of inertia
+//Refer fig.9.38
+//section can be divided into three triangles
+A1=80*12 //mm^2
+A2=(150-22)*12 //mm^2
+A3=120*10 //mm^2
+A=A1+A2+A3 //mm^2
+Ixx=((80*12^3)/12)+(960*74.22^2)+((12*128^3)/12)+(1536*4.22^2)+((120*10^3)/12)+(1200*64.78^2) //mm^4
+Iyy=((12*80^3)/12)+((128*12^3)/12)+((10*120^3)/12) //mm^4
+//Polar moment of Inertia (Izz)
+Izz=Ixx+Iyy //mm^4
+kxx=sqrt(Ixx/A) //mm
+Kyy=sqrt(Iyy/A) //mm The answer provided in the textbook is wrong
+printf("\nIzz=%.2d mm^4\nkxx=%.2f mm\nkyy=%.2f mm",Izz,kxx,kyy)
+
+
+
+
+
diff --git a/3886/CH9/EX9.13/9_13.txt b/3886/CH9/EX9.13/9_13.txt
new file mode 100644
index 000000000..553fb7d87
--- /dev/null
+++ b/3886/CH9/EX9.13/9_13.txt
@@ -0,0 +1,6 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.13.sce', -1)
+
+Izz=14440459 mm^4
+kxx=58.09 mm
+kyy=31.21 mm \ No newline at end of file
diff --git a/3886/CH9/EX9.14/9_14.sce b/3886/CH9/EX9.14/9_14.sce
new file mode 100644
index 000000000..f1a6022b9
--- /dev/null
+++ b/3886/CH9/EX9.14/9_14.sce
@@ -0,0 +1,13 @@
+//Determine moment of inertia
+//refer fig.9.39
+A1=100*30 //mm^2
+A2=100*25 //mm^2
+A3=200*20 //mm^2
+A4=87.5*20/2 //mm^2
+A5=87.5*20/2 //mm^2
+A=A1+A2+A3+A4+A5 //mm^2
+ybar=(3000*135+2500*70+4000*10+875*(20/3+20)*2)/A //mm
+Ixx=((100*30*30*30)/(12))+(3000*(75.74^2))+((25*(100^3))/(12))+(2500*(10.74^2))+((200*(20^3))/(12))+(4000*(49.26^2))+((87.5*(20^3))/(36))+(875*(32.59^2))+((87.5*(20^3))/(36))+(875*(32.59^2)) //mm^4
+Iyy=((30*(100^3))/(12))+((100*(25^3))/(12))+((20*(200^3))/(12))+((20*(87.5^3))/(36))+(875*(41.66^2))+((20*(87.5^3))/(36))+(875*(41.66^2)) //mm^4
+printf("\nIxx=%.2d mm^4\nIyy=%.2d mm^4",Ixx,Iyy)
+//The answers vary due to round off error
diff --git a/3886/CH9/EX9.14/9_14.txt b/3886/CH9/EX9.14/9_14.txt
new file mode 100644
index 000000000..8a4ef1863
--- /dev/null
+++ b/3886/CH9/EX9.14/9_14.txt
@@ -0,0 +1,5 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.14.sce', -1)
+
+Ixx=31543446 mm^4
+Iyy=19745121 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.15/9_15.sce b/3886/CH9/EX9.15/9_15.sce
new file mode 100644
index 000000000..c4a5b88a8
--- /dev/null
+++ b/3886/CH9/EX9.15/9_15.sce
@@ -0,0 +1,4 @@
+//Moment of inertia of section
+//refer fig. 9.40
+IAB=((400*20^3)/12)+(400*20*10^2)+(((100*10^3)/12)+100*10*(20+5)^2)*2+(((10*380^3)/12)+10*380*(30+190)^2)*2+(((100*10^3)/12)+100*10*(20+10+380+5)^2)*2 //mm^4
+printf("IAB=%.2d mm^4",IAB)
diff --git a/3886/CH9/EX9.15/9_15.txt b/3886/CH9/EX9.15/9_15.txt
new file mode 100644
index 000000000..3aa3569e1
--- /dev/null
+++ b/3886/CH9/EX9.15/9_15.txt
@@ -0,0 +1,3 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.15.sce', -1)
+IAB=806093333 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.16/9_16.sce b/3886/CH9/EX9.16/9_16.sce
new file mode 100644
index 000000000..eefe5eedd
--- /dev/null
+++ b/3886/CH9/EX9.16/9_16.sce
@@ -0,0 +1,7 @@
+//Moment of inertia
+//refer fig.9.41
+sumAiyi=250*10*5+2*40*10*(10+20)+40*10*(10+5)+40*10*255+250*10*(10+125) //mm^3
+A=2*250*10+40*10*4 //mm^2
+ybar=sumAiyi/A
+Ixx=((250*10^3)/12)+(250*10*(73.03-5)^2)+(((10*40^3)/12)+40*10*(73.03-30)^2)*2+((40*10^3)/12)+40*10*(73.03-15)^2+((10*250^3)/12)+250*10*(73.03-135)^2+((40*10^3)/12)+40*10*(73.03-255)^2 //mm^4
+printf("Ixx=%.2d mm^4",Ixx)
diff --git a/3886/CH9/EX9.16/9_16.txt b/3886/CH9/EX9.16/9_16.txt
new file mode 100644
index 000000000..60165db26
--- /dev/null
+++ b/3886/CH9/EX9.16/9_16.txt
@@ -0,0 +1,3 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.16.sce', -1)
+Ixx=50399393 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.17/9_17.sce b/3886/CH9/EX9.17/9_17.sce
new file mode 100644
index 000000000..c396c0e98
--- /dev/null
+++ b/3886/CH9/EX9.17/9_17.sce
@@ -0,0 +1,9 @@
+//moment of inertia
+//refer fig.9.42
+sumAiyi=(600*15)*((600/2)+20)+140*10*2*(70+30)+150*10*2*(5+20)+400*20*10 //mm^3
+A=600*15+140*10*2+150*10*2+400*20 //mm^2
+ybar=sumAiyi/A //mm
+Ixx=((15*(600^3))/(12))+(600*15*((145.39-320)^2))+((10*2*(140^3))/(12))+(1400*2*((145.39-100)^2))+((150*2*(10^3))/(12))+(1500*2*((145.39-15)^2))+((400*(20^3))/(12))+(400*20*((145.39-10)^2))
+printf("Ixx=%.2f mm^4",Ixx)
+//The answer provided in the textbook is wrong
+
diff --git a/3886/CH9/EX9.17/9_17.txt b/3886/CH9/EX9.17/9_17.txt
new file mode 100644
index 000000000..b9b667f95
--- /dev/null
+++ b/3886/CH9/EX9.17/9_17.txt
@@ -0,0 +1,3 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.17.sce', -1)
+Ixx=752679847.88 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.18/9_18.sce b/3886/CH9/EX9.18/9_18.sce
new file mode 100644
index 000000000..e4c6b47e3
--- /dev/null
+++ b/3886/CH9/EX9.18/9_18.sce
@@ -0,0 +1,4 @@
+//Compute moment of inertia
+//refer fig. 9.43
+Ixx=((125*60^3)/36)+(125*(60/2)*(60+60/3)^2)+((125*60^3)/36)+(125*(60/2)*(2*60/3)^2)+((125*60^3)/36)+(125*(60/2)*(60/3)^2)+((125*60^3)/36)+(125*(60/2)*(60/3)^2) //mm^4
+printf("Ixx=%.2d mm^4",Ixx)
diff --git a/3886/CH9/EX9.18/9_18.txt b/3886/CH9/EX9.18/9_18.txt
new file mode 100644
index 000000000..89bd66c2c
--- /dev/null
+++ b/3886/CH9/EX9.18/9_18.txt
@@ -0,0 +1,3 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.18.sce', -1)
+Ixx=36000000 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.19/9_19.sce b/3886/CH9/EX9.19/9_19.sce
new file mode 100644
index 000000000..04ea5aa68
--- /dev/null
+++ b/3886/CH9/EX9.19/9_19.sce
@@ -0,0 +1,5 @@
+//moment of inertia of shaded region
+//refer fig. 9.44
+//The figure is divided into simple geometry
+IAB=((80*80^3)/12)+((%pi*80^4)/128)-((%pi*40^4)/64) //mm^4
+printf("IAB=%.2d mm^4",IAB)
diff --git a/3886/CH9/EX9.19/9_19.txt b/3886/CH9/EX9.19/9_19.txt
new file mode 100644
index 000000000..1e51c7115
--- /dev/null
+++ b/3886/CH9/EX9.19/9_19.txt
@@ -0,0 +1,3 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.19.sce', -1)
+IAB=4292979 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.2/9_2.sce b/3886/CH9/EX9.2/9_2.sce
new file mode 100644
index 000000000..9f5a5dd70
--- /dev/null
+++ b/3886/CH9/EX9.2/9_2.sce
@@ -0,0 +1,11 @@
+//Centroid of angle
+//refer fig. 9.11
+//the given figure can be divided into two rectangles
+A1=150*12 //mm^2
+A2=(200-12)*12 //mm^2
+//total area
+A=A1+A2 //mm^2
+xbar=(1800*75+2256*6)/(4056) //mm
+ybar=(1800*6+2256*106)/(4056) //mm
+printf("The centroid is at \nxbar=%.2f mm\nybar=%.2f mm ",xbar,ybar)
+
diff --git a/3886/CH9/EX9.2/9_2.txt b/3886/CH9/EX9.2/9_2.txt
new file mode 100644
index 000000000..e16ddcce8
--- /dev/null
+++ b/3886/CH9/EX9.2/9_2.txt
@@ -0,0 +1,5 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.2.sce', -1)
+The centroid is at
+xbar=36.62 mm
+ybar=61.62 mm \ No newline at end of file
diff --git a/3886/CH9/EX9.20/9_20.sce b/3886/CH9/EX9.20/9_20.sce
new file mode 100644
index 000000000..8acc440c9
--- /dev/null
+++ b/3886/CH9/EX9.20/9_20.sce
@@ -0,0 +1,8 @@
+//Find second moment of inertia
+ybar=28.47 //mm
+xbar=39.21 //mm
+Ixx=((80*20^3)/36)+(80*(20/2)*(60-(2*20/3)-28.47)^2)+((80*40^3)/12)+(80*40*(28.47-20)^2)-((0.0068598*20^4)+(20^2)*(%pi/2)*(28.47-((4*20)/(3*%pi)))^2) //mm^4
+Iyy=((20*30^3)/36)+((20/2)*30*(39.21-(2*30/3))^2)+((20*50^3)/36)+(20/2)*50*(39.21-(30+50/3))^2+((40*80^3)/12)+(40*80*(39.21-40)^2)-(%pi*40^4)/(2*64)-((%pi)*(40^2)*(40-39.21)^2)/(4*2) //mm^4
+printf("\nIxx=%.2d mm^4\nIyy=%.2d mm^4",Ixx,Iyy)
+
+
diff --git a/3886/CH9/EX9.20/9_20.txt b/3886/CH9/EX9.20/9_20.txt
new file mode 100644
index 000000000..ba2e7ee4b
--- /dev/null
+++ b/3886/CH9/EX9.20/9_20.txt
@@ -0,0 +1,5 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.20.sce', -1)
+
+Ixx=686944 mm^4
+Iyy=1868392 mm^4 \ No newline at end of file
diff --git a/3886/CH9/EX9.3/9_3.sce b/3886/CH9/EX9.3/9_3.sce
new file mode 100644
index 000000000..928930bcc
--- /dev/null
+++ b/3886/CH9/EX9.3/9_3.sce
@@ -0,0 +1,16 @@
+//locating centroid
+//refer fig. 9.12
+//due to symmetry centroid must lie on y-axis
+xbar=0
+A1=100*20 //mm^2
+//for A1
+y1=30+100+20/2 //mm
+//for A2
+A2=100*20 //mm^2
+y2=30+100/2 //mm
+A3=150*30 //mm^2
+y3=30/2 //mm
+A=2000+2000+4500
+ybar=(A1*y1+A2*y2+A3*y3)/A //mm
+printf("The centroid is on the symmetric axis at a distance of %.2f mm from the bottom as shown in figure 9.12",ybar)
+
diff --git a/3886/CH9/EX9.3/9_3.txt b/3886/CH9/EX9.3/9_3.txt
new file mode 100644
index 000000000..6088ee43f
--- /dev/null
+++ b/3886/CH9/EX9.3/9_3.txt
@@ -0,0 +1,3 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.3.sce', -1)
+The centroid is on the symmetric axis at a distance of 59.71 mm from the bottom as shown in figure 9.12 \ No newline at end of file
diff --git a/3886/CH9/EX9.4/9_4.sce b/3886/CH9/EX9.4/9_4.sce
new file mode 100644
index 000000000..114872319
--- /dev/null
+++ b/3886/CH9/EX9.4/9_4.sce
@@ -0,0 +1,20 @@
+//Centroid of dam
+//refer fig. 9.13 and select axis accordingly
+A1=2*6/2 //m^2
+A2=2*7.5 //m^2
+A3=3*5/2 //m^2
+A4=1*4 //m^2
+A=A1+A2+A3+A4 //m^2
+//centroid of simple geometries are
+x1=2*2/3 //m
+y1=6/3 //m
+x2=2+1 //m
+y2=7.5/2 //m
+x3=2+2+3/3 //m
+y3=1+5/3 //m
+x4=4+4/2 //m
+y4=0.5 //m
+xbar=(A1*x1+A2*x2+A3*x3+A4*x4)/(A) //m
+ybar=(A1*y1+A2*y2+A3*y3+A4*y4)/(A) //m
+printf("centroid is at\nxbar=%.3f mm\nybar=%.3f mm",xbar,ybar)
+
diff --git a/3886/CH9/EX9.4/9_4.txt b/3886/CH9/EX9.4/9_4.txt
new file mode 100644
index 000000000..16dfa7458
--- /dev/null
+++ b/3886/CH9/EX9.4/9_4.txt
@@ -0,0 +1,5 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.4.sce', -1)
+centroid is at
+xbar=3.523 mm
+ybar=2.777 mm \ No newline at end of file
diff --git a/3886/CH9/EX9.5/9_5.sce b/3886/CH9/EX9.5/9_5.sce
new file mode 100644
index 000000000..c341eed72
--- /dev/null
+++ b/3886/CH9/EX9.5/9_5.sce
@@ -0,0 +1,21 @@
+//Determine centroid
+//refer fig. 9.14
+//This figure is divided into three simple figures
+A1=3*4/2 //m^2
+A2=6*4 //m^2
+A3=%pi*(1/2)*2^2 //m^2
+A=A1+A2+A3 //m^2
+//Co-ordinates of centroid
+x1=6+3/3 //m
+ya=4/3 //m
+x2=3 //m
+y2=2 //m
+R=2
+x3=(-4*R)/(3*%pi) //m
+y3=2 //m
+xbar=(A1*x1+A2*x2+A3*x3)/(A) //m
+ybar=(A1*ya+A2*y2+A3*y3)/(A) //m
+printf("\nxbar=%.3f m\nybar=%.3f m",xbar,ybar)
+
+
+
diff --git a/3886/CH9/EX9.5/9_5.txt b/3886/CH9/EX9.5/9_5.txt
new file mode 100644
index 000000000..b987a37fd
--- /dev/null
+++ b/3886/CH9/EX9.5/9_5.txt
@@ -0,0 +1,5 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.5.sce', -1)
+
+xbar=2.995 m
+ybar=1.890 m \ No newline at end of file
diff --git a/3886/CH9/EX9.6/9_6.sce b/3886/CH9/EX9.6/9_6.sce
new file mode 100644
index 000000000..d097410de
--- /dev/null
+++ b/3886/CH9/EX9.6/9_6.sce
@@ -0,0 +1,43 @@
+//Centroid of the gusset plate
+//refer fig. 9.15
+//The composite area is divided into algebraic sum and differences of simple geometries
+//for rectangle
+A1=160*280 //mm^2
+x1=140 //mm
+y1=80 //mm
+//for triangle
+A2=280*40/2 //mm^2
+x2=2*280/3 //mm
+y2=160+40/3 //mm
+//1st hole
+A3=(-%pi*21.5^2)/(4) //mm^2
+x3=70 //mm
+y3=50 //mm
+//second hole
+A4=-363.05 //mm^2
+x4=140 //mm
+y4=50 //mm
+//third hole
+A5=-363.05 //mm^2
+x5=210 //mm
+y5=50 //mm
+//fourth hole
+A6=-363.05 //mm^2
+x6=70 //mm
+y6=120 //mm
+//fifth hole
+A7=-363.05 //mm^2
+x7=140 //mm
+y7=130 //mm
+//sixth hole
+A8=-363.05 //mm^2
+x8=210 //mm
+y8=140 //mm
+A=A1+A2+A3+A4+A5+A6+A7+A8 //mm^2
+sumAixi=A1*x1+A2*x2+A3*x3+A4*x4+A5*x5+A6*x6+A7*x7+A8*x8 //mm^3
+xbar=sumAixi/A //mm
+sumAiyi=A1*y1+A2*y2+A3*y3+A4*y4+A5*y5+A6*y6+A7*y7+A8*y8
+ybar=sumAiyi/A //mm
+printf("\xbar=%.3f mm \nybar=%.3f mm",xbar,ybar)
+
+
diff --git a/3886/CH9/EX9.6/9_6.txt b/3886/CH9/EX9.6/9_6.txt
new file mode 100644
index 000000000..9fad73cab
--- /dev/null
+++ b/3886/CH9/EX9.6/9_6.txt
@@ -0,0 +1,4 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.6.sce', -1)
+xbar=145.419 mm
+ybar=90.387 mm \ No newline at end of file
diff --git a/3886/CH9/EX9.7/9_7.sce b/3886/CH9/EX9.7/9_7.sce
new file mode 100644
index 000000000..47a4f987a
--- /dev/null
+++ b/3886/CH9/EX9.7/9_7.sce
@@ -0,0 +1,8 @@
+//Determine co-ordinates
+//total area
+A=200*150-(100*75/2)-((%pi*100^2)/(4)) //mm^2
+xc=2375000/26250 //mm
+yc=1781250/26250 //mm
+printf("\xc=%.3f mm\nyc=%.3f mm",xc,yc)
+
+
diff --git a/3886/CH9/EX9.7/9_7.txt b/3886/CH9/EX9.7/9_7.txt
new file mode 100644
index 000000000..c4bac107a
--- /dev/null
+++ b/3886/CH9/EX9.7/9_7.txt
@@ -0,0 +1,4 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.7.sce', -1)
+xc=90.476 mm
+yc=67.857 mm \ No newline at end of file
diff --git a/3886/CH9/EX9.8/9_8.sce b/3886/CH9/EX9.8/9_8.sce
new file mode 100644
index 000000000..f8406efd5
--- /dev/null
+++ b/3886/CH9/EX9.8/9_8.sce
@@ -0,0 +1,25 @@
+//Locate centroid
+//refer fig.9.17
+x=40 //mm
+A1=168*x^2
+A2=12*x^2
+A3=-16*x^2
+A4=-8*%pi*x^2
+A5=-4*%pi*x^2
+x1=7*x
+x2=16*x
+x3=2*x
+x4=6*x
+x5=12.3023*x
+y1=6*x
+y2=4*x/3
+y3=10*x
+y4=(16*x/(3*%pi))
+y5=10.3023*x
+A=126.3009*x^2
+sumAixi=1030.6083*x^3
+sumAiyi=691.8708*x^3
+xbar=1030.6083*x/126.3009 //mm
+ybar=691.8708*x/126.3009 //mm
+printf("centroid is at (%.2f, %.2f)",xbar,ybar)
+
diff --git a/3886/CH9/EX9.8/9_8.txt b/3886/CH9/EX9.8/9_8.txt
new file mode 100644
index 000000000..b795417e5
--- /dev/null
+++ b/3886/CH9/EX9.8/9_8.txt
@@ -0,0 +1,3 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.8.sce', -1)
+centroid is at (326.40, 219.12) \ No newline at end of file
diff --git a/3886/CH9/EX9.9/9_9.sce b/3886/CH9/EX9.9/9_9.sce
new file mode 100644
index 000000000..b50f36305
--- /dev/null
+++ b/3886/CH9/EX9.9/9_9.sce
@@ -0,0 +1,13 @@
+//Determine moment of inertia
+//refer fig.9.34
+//composite section can be divided into simple ones
+A1=150*10 //mm^2
+A2=140*10 //mm^2
+A=A1+A2 //mm^2
+//due to symmetry centroid lies on y-axis
+ybar=(1500*5+1400*(10+70))/(2900) //mm
+Ixx=((150/12)*10^3)+(1500*36.21^2)+((10/12)*140^3)+(1400*38.79^2) //mm^4
+Iyy=((10*150^3)/(12))+((140*10^3)/(12)) //mm^4
+kxx=sqrt(Ixx/A) //mm
+kyy=sqrt(Iyy/A) //mm
+printf("\nIxx=%.2f mm^4\nIyy=%.2f mm^4\nkxx=%.2f mm\nkyy=%.2f mm",Ixx,Iyy,kxx,kyy)
diff --git a/3886/CH9/EX9.9/9_9.txt b/3886/CH9/EX9.9/9_9.txt
new file mode 100644
index 000000000..fbb8c2e0d
--- /dev/null
+++ b/3886/CH9/EX9.9/9_9.txt
@@ -0,0 +1,7 @@
+
+--> exec('E:\My program EM\9. Centroid and moment of inertia\9.9.sce', -1)
+
+Ixx=6372442.56 mm^4
+Iyy=2824166.67 mm^4
+kxx=46.88 mm
+kyy=31.21 mm \ No newline at end of file