1 files changed, 31 insertions, 0 deletions
diff --git a/3705/CH14/EX14.5/Ex14_5.sce b/3705/CH14/EX14.5/Ex14_5.sce
new file mode 100644
index 000000000..799cc3e53
--- /dev/null
+++ b/3705/CH14/EX14.5/Ex14_5.sce
@@ -0,0 +1,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)