diff options
Diffstat (limited to '564/CH1/EX1.2')
-rwxr-xr-x | 564/CH1/EX1.2/1_02.sce | 18 |
1 files changed, 18 insertions, 0 deletions
diff --git a/564/CH1/EX1.2/1_02.sce b/564/CH1/EX1.2/1_02.sce new file mode 100755 index 000000000..a4bf98e69 --- /dev/null +++ b/564/CH1/EX1.2/1_02.sce @@ -0,0 +1,18 @@ +pathname=get_absolute_file_path('1_02.sce')
+filename=pathname+filesep()+'1_02data.sci'
+exec(filename)
+printf("\nbending moment due to direct loading in a vertical plane: %f N/mm^2",Load*t);
+Sx1= Load/(%pi*(d^2)/4);//σx (axial load)
+printf("\nσx (axial load) is: %f N/mm^2",Sx1);
+I= %pi*(d^4)/64;// moment of Inertia
+Sx2= t*Load*(d/2)*(1/I);//σx (bending moment)
+printf("\nσx (bending moment) is: %f N/mm^2",Sx2)
+Sx=Sx1+Sx2;//σx
+J= %pi*(d^4)/32;//torsion constant
+Txy=T*(d/2)*(1/J);//τxy
+printf("\nSince the element is positioned at the bottom of the beam\nτxy:%fN/mm^2",-Txy);
+printf("\nσx: %f N/mm^2",-Sx);
+Sn=(-Sx)*(cos(theta))^2 +(-Txy)*sin(2*theta);//σn
+T= (-Sx-0)*sin(2*theta)/2 -(-Txy)*cos(2*theta);//τ
+printf("\nσn: %f N/mm^2",Sn);
+printf("\n τ: %f N/mm^2)",T);
\ No newline at end of file |