blob: 799cc3e53765fa7ea6610e9f57b12bd836fafea8 (
plain)
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
|
clear//
//Variable Declaration
Ix_bar=37.37*10**6 //Moment of inertia in mm^4
Iy_bar=21.07*10**6 //Moment of inertia in mm^4
Ixy_bar=-16.073*10**6 //Moment of inertia in mm^4
//Calculations
b=(Ix_bar+Iy_bar)*0.5 //Parameter for the circle in mm^4
R=sqrt(((Ix_bar-Iy_bar)*0.5)**2+Ixy_bar**2) //Radius of the Mohr's Circle in mm^4
//Part 1
I1=b+R //MI in mm^4
I2=b-R //MI in mm^4
theta1=asin(abs(Ixy_bar)/R)*180*%pi**-1*0.5 //Angle in degrees
theta2=theta1+90 //Angle in degrees
//Part 2
alpha=(100-theta1*2)*0.5 //Angle in degrees
Iu=(b)+R*(cos(alpha*%pi*180**-1)) //MI in mm^4
Iv=(b)-R*(cos(alpha*%pi*180**-1)) //MI in mm^4
Iuv=R*sin(2*alpha*%pi*180**-1) //MI in mm^4
//Result
printf("\n The Principal Moment of inertias are as follows")
printf("\n I1= %0.0f mm^4 and I2= %0.0f mm^4",I1,I2)
printf("\n Princial direction are theta1= %0.1f degrees theta2= %0.1f degrees" ,theta1,theta2)
printf("\n The moment of inertia along the uv-axis is %0.0f mm^4" ,Iuv)
|
