1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
|
clc;
n=4; // no of component
d=25;//mm, diameter of holes
t=12.5;//mm Thickness
l=100;//mm, length in z direction
r=50;//mm , radius of quarter circle
V=[l*r*t,1/4*%pi*(r^2)*t,-%pi/4*d^2*t,-%pi/4*d^2*t];//mm^3, Volumes of part I,II , III and IV
x=[t/2,(r+t)/2,t/2,t/2];//mm, x components of centroids of part I,II , III and IV respectively
y=[-r/2,-4/3*r/%pi,-r/2,-r/2,r-r/2];//mm, y components of centroids of part I,II , III and IV respectively
z=[56.25,6.25,87.5,37.5];//mm, z components of centroids of part I,II , III and IV respectively
sumV=0;
sumxV=0;
sumyV=0;
sumzV=0;
for(i=1:n)
sumV=sumV+V(i);
sumxV=sumxV+x(i)*V(i);
sumyV=sumyV+y(i)*V(i);
sumzV=sumzV+z(i)*V(i);
end
//Location of centre of gravity
X=sumxV/sumV;// X co-ordinate
Y=sumyV/sumV;// Y co-ordinate
Z=sumzV/sumV;// Z co-ordinate
printf("Co-ordinates of centroid are X= %.2f mm, Y= %.1f mm and Z= %.2f mm \n",X,Y,Z);
|
