summaryrefslogtreecommitdiff
path: root/3862/CH4
diff options
context:
space:
mode:
Diffstat (limited to '3862/CH4')
-rw-r--r--3862/CH4/EX4.1/Ex4_1.sce27
-rw-r--r--3862/CH4/EX4.10/Ex4_10.sce30
-rw-r--r--3862/CH4/EX4.11/Ex4_11.sce33
-rw-r--r--3862/CH4/EX4.12/Ex4_12.sce38
-rw-r--r--3862/CH4/EX4.13/Ex4_13.sce41
-rw-r--r--3862/CH4/EX4.15/Ex4_15.sce34
-rw-r--r--3862/CH4/EX4.16/Ex4_16.sce47
-rw-r--r--3862/CH4/EX4.17/Ex4_17.sce43
-rw-r--r--3862/CH4/EX4.18/Ex4_18.sce15
-rw-r--r--3862/CH4/EX4.19/Ex4_19.sce25
-rw-r--r--3862/CH4/EX4.2/Ex4_2.sce30
-rw-r--r--3862/CH4/EX4.20/Ex4_20.sce26
-rw-r--r--3862/CH4/EX4.21/Ex4_21.sce24
-rw-r--r--3862/CH4/EX4.22/Ex4_22.sce10
-rw-r--r--3862/CH4/EX4.23/Ex4_23.sce41
-rw-r--r--3862/CH4/EX4.27/Ex4_27.sce39
-rw-r--r--3862/CH4/EX4.3/Ex4_3.sce38
-rw-r--r--3862/CH4/EX4.5/Ex4_5.sce24
-rw-r--r--3862/CH4/EX4.6/Ex4_6.sce25
-rw-r--r--3862/CH4/EX4.7/Ex4_7.sce32
-rw-r--r--3862/CH4/EX4.8/Ex4_8.sce30
-rw-r--r--3862/CH4/EX4.9/Ex4_9.sce41
22 files changed, 693 insertions, 0 deletions
diff --git a/3862/CH4/EX4.1/Ex4_1.sce b/3862/CH4/EX4.1/Ex4_1.sce
new file mode 100644
index 000000000..265a46521
--- /dev/null
+++ b/3862/CH4/EX4.1/Ex4_1.sce
@@ -0,0 +1,27 @@
+clear
+//
+
+//variable declaration
+
+L1=600.0 //length of wire AB,mm
+L2=200.0 //length of wire BC,mm
+L3=300.0 //length of wire CD,mm
+theta=45*%pi/180
+//The wire is divided into three segments AB, BC and CD. Taking A as origin the coordinates of the centroids of AB, BC and CD are (X1,Y1),(X2,Y2),(X3,Y3)
+
+X1=300.0
+X2=600.0
+X3=600.0-150*cos(theta)
+Y1=0
+Y2=100
+Y3=200+150*sin(theta)
+L=L1+L2+L3 //Total length,mm
+
+xc=(L1*X1+L2*X2+L3*X3)/L
+
+printf("\n xc= %0.2f mm",xc)
+
+
+yc=(L1*Y1+L2*Y2+L3*Y3)/L
+
+printf("\n yc= %0.2f mm",yc)
diff --git a/3862/CH4/EX4.10/Ex4_10.sce b/3862/CH4/EX4.10/Ex4_10.sce
new file mode 100644
index 000000000..bff9ed699
--- /dev/null
+++ b/3862/CH4/EX4.10/Ex4_10.sce
@@ -0,0 +1,30 @@
+clear
+// If xc and yc are the coordinates of the centre of the circle, centroid also must have the coordinates xc and yc as per the condition laid down in the problem. The shaded area may be considered as a rectangle of size 200 mm × 150 mm minus a triangle of sides 100 mm × 75 mm and a circle of diameter 100 mm.
+
+//
+//variable declaration
+
+Ap=200.0*150.0 //Area of plate,mm^2
+At=100.0*75.0/2.0 //Area of triangle,mm^2
+Ah=%pi*(100**2)/4.0 //Area of hole **mm^2
+
+A=Ap-At-Ah
+
+
+X1=100.0
+X2=200.0-100.0/3.0
+//X3=Xc
+
+Y1=75.0
+Y2=150.0-25.0
+//Y3=Yc
+
+A=Ap-At-Ah
+
+xc=(Ap*X1-At*X2)/(Ah+A)
+
+printf("\n xc= %0.2f mm",xc)
+
+yc=(Ap*Y1-At*Y2)/(Ah+A)
+
+printf("\n yc= %0.2f mm",yc)
diff --git a/3862/CH4/EX4.11/Ex4_11.sce b/3862/CH4/EX4.11/Ex4_11.sce
new file mode 100644
index 000000000..c39288a06
--- /dev/null
+++ b/3862/CH4/EX4.11/Ex4_11.sce
@@ -0,0 +1,33 @@
+clear
+//
+
+//variable declaration
+X=40.0
+A1=14.0*12.0*(X**2) //Area of rectangle**mm^2
+A2=6.0*4.0*(X**2)/2.0 //Area of triangle**mm^2
+A3=-4*4*(X**2) //Area of removed subtracted**mm^2
+A4=-%pi*(4*X**2)/2.0 //Area of semicircle to be subtracted**mm^2
+A5=-%pi*(4*X**2)/4.0 //Area of quarter of circle to be subtracted**mm^2
+
+X1=7.0*X
+X2=14*X+2*X
+X3=2*X
+X4=6.0*X
+X5=14.0*X-(16*X/(3*%pi))
+
+Y1=6.0*X
+Y2=4.0*X/3.0
+Y3=8.0*X+2.0*X
+Y4=(16.0*X)/(3*%pi)
+Y5=12*X-4*(4*X/(3*%pi))
+
+A=A1+A2+A3+A4+A5
+
+xc=(A1*X1+A2*X2+A3*X3+A4*X4+A5*X5)/A
+
+printf("\n xc= %0.2f m",xc)
+
+
+yc=(A1*Y1+A2*Y2+A3*Y3+A4*Y4+A5*Y5)/A
+
+printf("\n yc= %0.2f m",yc)
diff --git a/3862/CH4/EX4.12/Ex4_12.sce b/3862/CH4/EX4.12/Ex4_12.sce
new file mode 100644
index 000000000..571e8eddb
--- /dev/null
+++ b/3862/CH4/EX4.12/Ex4_12.sce
@@ -0,0 +1,38 @@
+clear
+//The given composite section can be divided into two rectangles
+
+//
+//variable declaration
+
+
+A1=150.0*10.0 //Area of 1,mm^2
+A2=140.0*10.0 //Area of 2,mm^2
+A=A1+A2 //Total area,mm^2
+//Due to symmetry, centroid lies on the symmetric axis y-y. The distance of the centroid from the top most fibre is given by:
+
+Y1=5.0
+Y2=10.0+70.0
+
+yc=(A1*Y1+A2*Y2)/A
+
+//Referring to the centroidal axis x-x and y-y, the centroid of A1 is g1 (0.0, yc-5) and that of A2 is g2 (0.0, 80-yc)
+
+//Moment of inertia of the section about x-x axis Ixx = moment of inertia of A1 about x-x axis + moment of inertia of A2 about x-x axis.
+
+
+Ixx=(150*(10**3)/12)+(A1*((yc-5)**2))+(10*(140**3)/12)+(A2*((80-yc)**2))
+
+printf("\n Ixx= %0.1f mm^4",Ixx)
+
+Iyy=(10*(150**3)/12)+(140*(10**3)/12)
+
+printf("\n Iyy= %0.1f mm^4",Iyy)
+
+//Hence, the moment of inertia of the section about an axis passing through the centroid and parallel to the top most fibre is Ixxmm^4 and moment of inertia of the section about the axis of symmetry is Iyy mm^4.
+//The radius of gyration is given by
+
+kxx=sqrt(Ixx/A)
+printf("\n kxx= %0.2f mm",kxx)
+
+kyy=sqrt(Iyy/A)
+printf("\n kyy= %0.2f mm",kyy)
diff --git a/3862/CH4/EX4.13/Ex4_13.sce b/3862/CH4/EX4.13/Ex4_13.sce
new file mode 100644
index 000000000..10766b63c
--- /dev/null
+++ b/3862/CH4/EX4.13/Ex4_13.sce
@@ -0,0 +1,41 @@
+clear
+//The given composite section can be divided into two rectangles
+
+
+//variable declaration
+
+
+A1=125.0*10.0 //Area of 1,mm^2
+A2=75.0*10.0 //Area of 2,mm^2
+A=A1+A2 //Total area,mm^2
+
+//First, the centroid of the given section is to be located. Two reference axis (1)–(1) and (2)–(2)
+
+//The distance of centroid from the axis (1)–(1)
+
+X1=5.0
+X2=10.0+75.0/2
+
+xc=(A1*X1+A2*X2)/A
+
+//Similarly, the distance of the centroid from the axis (2)–(2)
+
+Y1=125.0/2
+Y2=5.0
+
+yc=(A1*Y1+A2*Y2)/A
+
+//With respect to the centroidal axis x-x and y-y, the centroid of A1 is g1 (xc-5, (85/2)-xc) and that of A2 is g2 ((135/2)-yc, yc-5).
+Ixx=(10*(125**3)/12)+(A1*(21.56**2))+(75.0*(10.0**3.0)/12)+(A2*((39.94)**2))
+
+printf("\n Ixx= %0.1f mm^4",Ixx)
+
+Iyy=(125*(10**3)/12)+(A1*(15.94**2))+(10*(75**3)/12)+(A2*(26.56**2))
+
+printf("\n Iyy= %0.1f mm^4",Iyy)
+
+//Izz=Polar moment of inertia
+
+Izz=Ixx+Iyy
+
+printf("\n Izz= %0.1f mm^4",Izz) \ No newline at end of file
diff --git a/3862/CH4/EX4.15/Ex4_15.sce b/3862/CH4/EX4.15/Ex4_15.sce
new file mode 100644
index 000000000..405e99812
--- /dev/null
+++ b/3862/CH4/EX4.15/Ex4_15.sce
@@ -0,0 +1,34 @@
+clear
+//
+
+//variable declaration
+
+
+A1=100.0*13.5 //Area of 1,mm^2
+A2=(400.0-27.0)*8.1 //Area of 2,mm^2
+A3=100.0*13.5 //Area of 3,mm^2
+
+A=A1+A2+A3 //Total area,mm^2
+
+//The given section is symmetric about horizontal axis passing through the centroid g2 of the rectangle A2.
+
+X1=50.0
+X2=8.1/2.0
+X3=50.0
+
+xc=(A1*X1+A2*X2+A3*X3)/A
+
+Y1=386.5+13.5/2.0
+Y2=200.0
+Y3=13.5/2
+
+yc=(A1*Y1+A2*Y2+A3*Y3)/A
+
+//With reference to the centroidal axis x-x and y-y
+
+Ixx=(100.0*(13.5**3)/12.0)+(A1*((200-(13.5/2))**2))+(8.1*(373**3.0)/12.0)+0+(100*(13.5**3)/12.0)+(A3*((200-(13.5/2))**2))
+printf("\n Ixx= %0.1f mm^4",Ixx)
+
+Iyy=(13.5*(100.0**3)/12.0)+(A1*((50-xc)**2))+(373.0*(8.1**3.0)/12.0)+A2*(21.68**2)+(13.5*(100**3)/12.0)+(A3*((50-xc)**2))
+
+printf("\n Iyy= %0.1f mm^4",Iyy)
diff --git a/3862/CH4/EX4.16/Ex4_16.sce b/3862/CH4/EX4.16/Ex4_16.sce
new file mode 100644
index 000000000..ce6ff8492
--- /dev/null
+++ b/3862/CH4/EX4.16/Ex4_16.sce
@@ -0,0 +1,47 @@
+clear
+// The section is divided into three rectangles A1, A2 and A3
+
+//
+
+//variable declaration
+
+
+A1=80.0*12.0 //Area of 1,mm^2
+A2=(150.0-22.0)*12.0 //Area of 2,mm^2
+A3=120.0*10.0 //Area of 3,mm^2
+
+A=A1+A2+A3 //Total area,mm^2
+
+//Due to symmetry, centroid lies on axis y-y. The bottom fibre (1)–(1) is chosen as reference axis to locate the centroid
+
+Y1=150-6
+Y2=(128/2) +10
+Y3=5
+X1=60.0
+X2=60.0
+X3=60.0
+
+yc=(A1*X1+A2*X2+A3*X3)/A
+
+
+
+xc=(A1*Y1+A2*Y2+A3*Y3)/A
+
+//With reference to the centroidal axis x-x and y-y, the centroid of the rectangles A1 is g1 (0.0, 150-6-yc), that of A2 is g2 (0.0, 75-yc) and that of A3 is g3 (0.0, yc-5 ).
+
+Iyy=(12*((80**3))/12)+(128*((12**3))/12)+(10*((120**3))/12)
+
+Ixx=(80.0*(12.0**3)/12.0)+(A1*((150-6-yc)**2))+(12*(128**3.0)/12.0)+(A2*((75-yc)**2))+(120*(10**3)/12.0)+(A3*((150-10-6-yc)**2))
+
+
+
+PolarmomentofInertia=Ixx+Iyy
+
+printf("\n Polar moment of Inertia= %0.0f mm^4",PolarmomentofInertia)
+
+kxx=sqrt(Ixx/A)
+printf("\n kxx= %0.2f mm",kxx)
+
+
+kyy=sqrt(Iyy/A)
+printf("\n kyy= %0.2f mm",kyy)
diff --git a/3862/CH4/EX4.17/Ex4_17.sce b/3862/CH4/EX4.17/Ex4_17.sce
new file mode 100644
index 000000000..ac0c13a12
--- /dev/null
+++ b/3862/CH4/EX4.17/Ex4_17.sce
@@ -0,0 +1,43 @@
+clear
+//
+
+//The given composite section may be divided into simple rectangles and triangle
+
+//variable declaration
+
+
+A1=100.0*30.0 //Area of 1,mm^2
+A2=100.0*25.0 //Area of 2,mm^2
+A3=200.0*20.0 //Area of 3,mm^2
+A4=87.5*20.0/2.0 //Area of 4,mm^2
+A5=87.5*20.0/2.0 //Area of 5,mm^2
+
+A=A1+A2+A3+A4+A5 //Total area,mm^2
+
+//Due to symmetry, centroid lies on the axis y-y. A reference axis (1)–(1) is choosen as shown in the figure. The distance of the centroidal axis from (1)–(1)
+
+X1=100.0
+X2=100.0
+X3=100.0
+X4=2.0*87.5/3.0
+X5=200-X4
+xc=(A1*X1+A2*X2+A3*X3+A4*X4+A5*X5)/A
+
+Y1=135.0
+Y2=70.0
+Y3=10.0
+Y4=(20.0/3.0)+20.0
+Y5=Y4
+
+yc=(A1*Y1+A2*Y2+A3*Y3+A4*Y4+A5*Y5)/A
+
+//With reference to the centroidal axis x-x and y-y, the centroid of the rectangle A1 is g1 (0.0,135.0-yc), that of A2 is g2(0.0,70.00-yc), that of A3 is g3 (0.0, yc-10.0), the centroid of triangle A4 is g4 (41.66,yc-20.0-(20.0/3.0) ) and that of A5 is g5 (41.66,yc-20.0-(20.0/3.0)).
+
+
+Ixx=(100.0*(30**3)/12.0)+(A1*((135.0-yc)**2))+(25.0*(100**3.0)/12.0)+(A2*((70.0-yc)**2))+(200*(20**3)/12.0)+(A3*((yc-10.0)**2))+((87.5*(20**3)/36.0)+(A4*((yc-20.0-(20.0/3.0))**2)))*2
+
+printf("\n Ixx= %0.1f mm^4",Ixx)
+
+Iyy=(30.0*(100**3)/12.0)+(100.0*(25**3.0)/12.0)+(20*(200**3)/12.0)+((20.0*(87.5**3)/36.0)+(A4*((41.66)**2)))*2
+
+printf("\n Iyy= %0.1f mm^4",Iyy)
diff --git a/3862/CH4/EX4.18/Ex4_18.sce b/3862/CH4/EX4.18/Ex4_18.sce
new file mode 100644
index 000000000..a9d56d4cd
--- /dev/null
+++ b/3862/CH4/EX4.18/Ex4_18.sce
@@ -0,0 +1,15 @@
+clear
+//In this problem, it is required to find out the moment of inertia of the section about an axis AB. So there is no need to find out the position of the centroid.
+//The given section is split up into simple rectangles
+//Moment of inertia about AB = Sum of moments of inertia of the rectangle about AB
+
+//variable declaration
+
+A1=400*20.0
+A2=100*10
+A3=10*380.0
+A4=100*10.0
+
+IAB=(400.0*(20**3)/12)+(A1*(10**2))+((100*(10**3)/12)+(A2*(25**2)))*2+((10*(380**3)/12)+(A3*(220**2)))*2+((100*(10**3)/12)+(A4*(415**2)))*2
+
+printf("\n IAB= %0.0f mm^4",IAB)
diff --git a/3862/CH4/EX4.19/Ex4_19.sce b/3862/CH4/EX4.19/Ex4_19.sce
new file mode 100644
index 000000000..e55f95842
--- /dev/null
+++ b/3862/CH4/EX4.19/Ex4_19.sce
@@ -0,0 +1,25 @@
+clear
+// The built-up section is divided into six simple rectangles
+
+//variable declaration
+
+
+A1=250.0*10.0 //Area of 1,mm^2
+A2=40.0*10.0 //Area of 2,mm^2
+
+A=A1*2+A2*4 //Total area,mm^2
+
+
+Y1=5.0
+Y2=30.0
+Y3=15.0
+Y4=255.0
+Y5=135.0
+
+yc=(A1*Y1+2*A2*Y2+A2*Y3+A2*Y4+A1*Y5)/A
+
+//Now, Moment of inertia about the centroidalaxis=Sum of the moment of inertia of the individual rectangles
+
+Ixx=(250.0*(10**3)/12.0)+(A1*((yc-5)**2))+((10.0*(40**3.0)/12.0)+(A2*((yc-30.0)**2)))*2+(40*(10**3)/12.0)+(A2*((yc-15.0)**2))+(10.0*(250.0**3.0)/12.0)+(A1*((yc-135.0)**2))+(40.0*(10.0**3)/12)+(A2*((yc-255)**2))
+
+printf("\n Ixx= %0.1f mm^4",Ixx)
diff --git a/3862/CH4/EX4.2/Ex4_2.sce b/3862/CH4/EX4.2/Ex4_2.sce
new file mode 100644
index 000000000..fa6b88563
--- /dev/null
+++ b/3862/CH4/EX4.2/Ex4_2.sce
@@ -0,0 +1,30 @@
+clear
+//
+
+//The composite figure is divided into three simple figures and taking A as origin coordinates of their centroids
+
+//variable declaration
+
+L1=400.0 //length of wire AB,mm
+L2=150.0*%pi //length of wire BC,mm
+L3=250.0 //length of wire CD,mm
+theta=30*%pi/180
+
+//The wire is divided into three segments AB, BC and CD. Taking A as origin the coordinates of the centroids of AB, BC and CD are (X1,Y1),(X2,Y2),(X3,Y3)
+X1=200.0
+X2=475.0
+X3=400+300.0+250*cos(theta)/2
+
+Y1=0
+Y2=2*150/%pi
+Y3=125*sin(theta)
+L=L1+L2+L3 //Total length,mm
+
+xc=(L1*X1+L2*X2+L3*X3)/L
+
+printf("\n xc= %0.2f mm",xc)
+
+
+yc=(L1*Y1+L2*Y2+L3*Y3)/L
+
+printf("\n yc= %0.2f mm",yc)
diff --git a/3862/CH4/EX4.20/Ex4_20.sce b/3862/CH4/EX4.20/Ex4_20.sce
new file mode 100644
index 000000000..6d9ef22f4
--- /dev/null
+++ b/3862/CH4/EX4.20/Ex4_20.sce
@@ -0,0 +1,26 @@
+clear
+//Each angle is divided into two rectangles
+
+//variable declaration
+
+A1=600.0*15.0 //Area of 1,mm^2
+A2=140.0*10.0 //Area of 2,mm^2
+A3=150.0*10.0
+A4=400.0*20.0
+A=A1+A2*2+A3*2+A4 //Total area,mm^2
+
+//The distance of the centroidal axis from the bottom fibres of section
+
+Y1=320.0
+Y2=100.0
+Y3=25.0
+Y4=10.0
+
+yc=(A1*Y1+2*A2*Y2+A3*Y3*2+A4*Y4)/A
+printf("\n yc")
+//Now, Moment of inertia about the centroidalaxis=Sum of the moment of inertia of the individual rectangles
+
+Ixx=(15.0*(600**3)/12.0)+(A1*((yc-320)**2))+((10.0*(140**3.0)/12.0)+(A2*((yc-100.0)**2)))*2+((150*(10**3)/12.0)+(A3*((yc-15.0)**2)))*2+(400.0*(20.0**3.0)/12.0)+(A4*((yc-10.0)**2))
+
+
+printf("\n Ixx= %0.1f mm^4",Ixx)
diff --git a/3862/CH4/EX4.21/Ex4_21.sce b/3862/CH4/EX4.21/Ex4_21.sce
new file mode 100644
index 000000000..202ef8b1e
--- /dev/null
+++ b/3862/CH4/EX4.21/Ex4_21.sce
@@ -0,0 +1,24 @@
+clear
+//
+
+//The rectangle is divided into four triangles
+//The lines AE and FC are parallel to x-axis
+
+//variable declaration
+
+theta=asin(4.0/5.0)
+
+AB=100.0
+BK=AB*sin((90*%pi/180)-theta)
+ND=BK
+FD=60.0/sin(theta)
+AF=150.0-FD
+AE=AB/cos((90*%pi/180)-theta)
+FC=AE
+A=125.0*60.0/2.0
+
+//Moment of inertia of the section about axis x-x=Sum of the momentsof inertia of individual triangular areasabout axis
+
+Ixx=(125*(60**3)/36)+(A*((ND*4.0/3.0)**2))+(125*(60**3)/36)+(A*((ND*2.0/3.0)**2))+(125*(60**3)/36)+(A*((ND*1.0/3.0)**2))+(125*(60**3)/36)+(A*((ND*1.0/3.0)**2))
+
+printf("\n Ixx= %0.0f mm^4",Ixx)
diff --git a/3862/CH4/EX4.22/Ex4_22.sce b/3862/CH4/EX4.22/Ex4_22.sce
new file mode 100644
index 000000000..7a011265b
--- /dev/null
+++ b/3862/CH4/EX4.22/Ex4_22.sce
@@ -0,0 +1,10 @@
+clear
+//
+
+//The section is divided into a triangle PQR, a semicircle PSQ having base on axis AB and a circle having its centre on axis AB
+
+//variable declaration
+//Now,Moment of inertia of the section about axis AB
+IAB=(80*(80**3)/12)+(%pi*(80**4)/128)-(%pi*(40**4)/64)
+
+printf("\n IAB= %0.0f mm^4",IAB)
diff --git a/3862/CH4/EX4.23/Ex4_23.sce b/3862/CH4/EX4.23/Ex4_23.sce
new file mode 100644
index 000000000..4ffbca00e
--- /dev/null
+++ b/3862/CH4/EX4.23/Ex4_23.sce
@@ -0,0 +1,41 @@
+clear
+//
+
+//The section is divided into three simple figures viz., a triangle ABC, a rectangle ACDE and a semicircle.
+
+//variable declaration
+
+r=20.0 //radius of semicircle
+A1=80.0*20.0/2 //Area of triangle ABC
+A3=40.0*80.0 //Area of rectangle ACDE
+A4=%pi*(r**2)/2 //Area of semicircle
+At1=30.0*20.0/2.0
+At2=50.0*20.0/2.0
+A=A1+A3-A4 //Total area
+
+X1=2.0*30.0/3.0
+X2=50.0*30.0/3.0
+X3=40.0
+X4=40.0
+
+xc=(At1*X1+At2*X2+A3*X3-A4*X4)/A
+//mistake in book
+
+Y1=(20.0/3.0)+40.0
+Y3=20.0
+Y4=(4.0*20.0)/(3.0*%pi)
+
+yc=(A1*Y1+A3*Y3-A4*Y4)/A
+printf("\n %0.3f %0.3f ",xc,yc)
+
+//
+//Moment of inertia of the section about axis x-x=Sum of the momentsof inertia of individual triangular areasabout axis
+
+Ixx=(80.0*(20.0**3)/36) +A1*((60.0-(2*20.0/3.0)-yc)**2)+(80*(40**3)/12)+(A3*((yc-20.0)**2))-((0.0068598*(20**4))+(A4*((yc-Y4)**2)))
+
+printf("\n Ixx= %0.0f mm^4",Ixx)
+
+
+Iyy=(20.0*(30.0**3)/36) +At1*((xc-(2*30.0/3.0))**2)+(20*(50**3)/36)+(At2*((xc-(30.0+(50/3)))**2))+((40*(80**3)/12)+(A3*((xc-40)**2)))-((%pi*(40**4))/(2*64))-(A4*((40-xc)**2))
+
+printf("\n Iyy= %0.0f mm^4",Iyy)
diff --git a/3862/CH4/EX4.27/Ex4_27.sce b/3862/CH4/EX4.27/Ex4_27.sce
new file mode 100644
index 000000000..615d9a78d
--- /dev/null
+++ b/3862/CH4/EX4.27/Ex4_27.sce
@@ -0,0 +1,39 @@
+clear
+//A concrete block of size 0.60 m × 0.75 m × 0.5 m is cast with a hole of diameter 0.2 m and depth 0.3 m
+//The hole is completely filled with steel balls weighing 2500 N. Locate the centre of gravity of the body.
+
+//
+
+//variable declaration
+
+W=25000.0 // weight of concrete=25000, N/m^3
+W1=0.6*0.75*0.5*W //Weight of solid concrete block
+W2=%pi*(0.2**2)*0.3*W/4 //Weight of concrete (W2) removed for making hole:
+W3=2500
+
+//Taking origin as shown in the figure, the centre of gravity of solid block is (0.375, 0.3, 0.25) and that of hollow portion is (0.5, 0.4, 0.15)
+
+X1=0.375
+X2=0.5
+X3=0.5
+
+Y1=0.3
+Y2=0.4
+Y3=0.4
+
+Z1=0.25
+Z2=0.15
+Z3=0.15
+
+Wt=W3+W1-W2
+printf("\n %0.3f %0.3f %0.3f %0.3f ",W,W1,W2,Wt)
+
+xc=(W1*X1-W2*X2+W3*X3)/Wt
+
+yc=(W1*Y1-W2*Y2+W3*Y3)/Wt
+
+zc=(W1*Z1-W2*Z2+W3*Z3)/Wt
+
+printf("\n xc= %0.3f m",xc)
+printf("\n yc= %0.3f m",yc)
+printf("\n zc= %0.3f m",zc)
diff --git a/3862/CH4/EX4.3/Ex4_3.sce b/3862/CH4/EX4.3/Ex4_3.sce
new file mode 100644
index 000000000..1bac15a13
--- /dev/null
+++ b/3862/CH4/EX4.3/Ex4_3.sce
@@ -0,0 +1,38 @@
+clear
+//
+// The length and the centroid of portions AB, BC and CD
+// portion AB is in x-z plane, BC in y-z plane and CD in x-y plane. AB and BC are semi circular in shape
+
+//variable declaration
+
+L1=100.0*%pi //length of wire AB,mm
+L2=140.0*%pi //length of wire BC,mm
+L3=300.0 //length of wire CD,mm
+theta=45*%pi/180
+
+//The wire is divided into three segments AB, BC and CD. Taking A as origin the coordinates of the centroids of AB, BC and CD are (X1,Y1),(X2,Y2),(X3,Y3)
+X1=100.0
+X2=0
+X3=300*sin(theta)
+
+Y1=0
+Y2=140
+Y3=280+300*cos(theta)
+Z1=2*100/%pi
+Z2=2*140/%pi
+Z3=0
+
+L=L1+L2+L3 //Total length,mm
+
+xc=(L1*X1+L2*X2+L3*X3)/L
+
+printf("\n xc= %0.2f mm",xc)
+
+
+yc=(L1*Y1+L2*Y2+L3*Y3)/L
+
+printf("\n yc= %0.2f mm",yc)
+
+zc=(L1*Z1+L2*Z2+L3*Z3)/L
+
+printf("\n zc= %0.2f mm",zc)
diff --git a/3862/CH4/EX4.5/Ex4_5.sce b/3862/CH4/EX4.5/Ex4_5.sce
new file mode 100644
index 000000000..60a719ee0
--- /dev/null
+++ b/3862/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,24 @@
+clear
+//variable declaration
+
+A1=150.0*12.0 //Area of 1 ,mm^2
+A2=(200.0-12.0)*12.0 //Area of 2,mm^2
+
+X1=75
+X2=6
+
+Y1=6
+Y2=12+(200-12)/2
+
+A=A1+A2
+
+xc=(A1*X1+A2*X2)/A
+
+printf("\n xc= %0.2f ",xc)
+
+yc=(A1*Y1+A2*Y2)/A
+
+printf("\n yc= %0.2f mm",yc)
+
+printf("\nThus, the centroid is at x = 36.62 mm and y = 61.62 mm ")
+
diff --git a/3862/CH4/EX4.6/Ex4_6.sce b/3862/CH4/EX4.6/Ex4_6.sce
new file mode 100644
index 000000000..c52f5dddc
--- /dev/null
+++ b/3862/CH4/EX4.6/Ex4_6.sce
@@ -0,0 +1,25 @@
+clear
+//variable declaration
+
+A1=100.0*20 //Area of 1 ,mm^2
+A2=100.0*20.0 //Area of 2,mm^2
+A3=150.0*30.0 //Area of 3,mm^2
+
+//Selecting the coordinate system, due to symmetry centroid must lie on y axis,
+
+X1=0
+X2=0
+
+Y1=30+100+20/2
+Y2=30+100/2
+Y3=30/2
+
+A=A1+A2+A3
+
+
+yc=(A1*Y1+A2*Y2+A3*Y3)/A
+
+printf("\n yc= %0.2f mm",yc)
+
+printf("\n Thus, the centroid is on the symmetric axis at a distance 59.71 mm from the bottom")
+
diff --git a/3862/CH4/EX4.7/Ex4_7.sce b/3862/CH4/EX4.7/Ex4_7.sce
new file mode 100644
index 000000000..bdd81f46d
--- /dev/null
+++ b/3862/CH4/EX4.7/Ex4_7.sce
@@ -0,0 +1,32 @@
+clear
+// Note that it is convenient to take axis in such a way that the centroids of all simple figures are having positive coordinates. If coordinate of any simple figure comes out to be negative, one should be careful in assigning the sign of moment of area
+
+//variable declaration
+
+A1=2.0*6.0*1.0/2.0 //Area of 1,m^2
+A2=2.0*7.5 //Area of 2,m^2
+A3=3.0*5.0*1.0/2 //Area of 3,m^2
+A4=1.0*4.0 //Area of 4,m^2
+
+//The composite figure can be conveniently divided into two triangles and two rectangle
+
+X1=2.0*2.0/3.0
+X2=2.0+1.0
+X3=2.0+2.0+(1.0*3.0/3.0)
+X4=4.0+4.0/2.0
+
+Y1=6.0/3.0
+Y2=7.5/2.0
+Y3=1.0+5.0/3.0
+Y4=1/2.0
+
+A=A1+A2+A3+A4
+
+xc=(A1*X1+A2*X2+A3*X3+A4*X4)/A
+
+printf("\n xc= %0.3f m",xc)
+
+
+yc=(A1*Y1+A2*Y2+A3*Y3+A4*Y4)/A
+
+printf("\n yc= %0.3f m",yc)
diff --git a/3862/CH4/EX4.8/Ex4_8.sce b/3862/CH4/EX4.8/Ex4_8.sce
new file mode 100644
index 000000000..2d226987d
--- /dev/null
+++ b/3862/CH4/EX4.8/Ex4_8.sce
@@ -0,0 +1,30 @@
+clear
+// The composite section is divided into three simple figures, a triangle, a rectangle and a semicircle
+
+//
+//variable declaration
+
+A1=1.0*3.0*4.0/2.0 //Area of 1,m^2
+A2=6.0*4.0 //Area of 2,m^2
+A3=1.0*%pi*(2**2)/2 //Area of 3**m^2
+
+//The coordinates of centroids of these three simple figures are:
+
+X1=6.0+3.0/3.0
+X2=3.0
+X3=-(4*2)/(3.0*%pi)
+
+Y1=4.0/3.0
+Y2=2.0
+Y3=2.0
+
+A=A1+A2+A3
+
+xc=(A1*X1+A2*X2+A3*X3)/A
+
+printf("\n xc= %0.4f m",xc)
+
+
+yc=(A1*Y1+A2*Y2+A3*Y3)/A
+
+printf("\n yc= %0.3f m",yc)
diff --git a/3862/CH4/EX4.9/Ex4_9.sce b/3862/CH4/EX4.9/Ex4_9.sce
new file mode 100644
index 000000000..0177d4636
--- /dev/null
+++ b/3862/CH4/EX4.9/Ex4_9.sce
@@ -0,0 +1,41 @@
+clear
+//The composite area is equal to a rectangle of size 160 × 280 mm plus a triangle of size 280 mm base width and 40 mm height and minus areas of six holes. In this case also the can be used for locating centroid by treating area of holes as negative. The area of simple figures and their centroids are
+
+//
+
+//variable declaration
+
+Ar=160.0*280.0 //Area of rectangle,mm^2
+At=280.0*40.0/2.0 //Area of triangle,mm^2
+d=21.5 //diameter of hole,mm
+Ah=-%pi*(d**2)/4 //Area of hole**mm^2
+
+A=Ar+At+Ah*6
+
+
+Xr=140.0
+Xt=560/3.0
+Xh1=70.0
+Xh2=140.0
+Xh3=210.0
+Xh4=70.0
+Xh5=140.0
+Xh6=210.0
+
+Yr=80.0
+Yt=160.0+40.0/3.0
+Yh1=50.0
+Yh2=50.0
+Yh3=50.0
+Yh4=120.0
+Yh5=130.0
+Yh6=140.0
+
+xc=(Ar*Xr+At*Xt+Ah*(Xh1+Xh2+Xh3+Xh4+Xh5+Xh6))/A
+
+printf("\n xc= %0.2f m",xc)
+
+
+yc=(Ar*Yr+At*Yt+Ah*(Yh1+Yh2+Yh3+Yh4+Yh5+Yh6))/A
+
+printf("\n yc= %0.2f m",yc)
generated by cgit (git 2.34.1) at 2024-12-19 12:27:21 +0000