diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /564/CH1/EX1.3 | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '564/CH1/EX1.3')
-rwxr-xr-x | 564/CH1/EX1.3/1_03.sce | 27 | ||||
-rwxr-xr-x | 564/CH1/EX1.3/1_03MohrCircle.JPG | bin | 0 -> 72564 bytes |
2 files changed, 27 insertions, 0 deletions
diff --git a/564/CH1/EX1.3/1_03.sce b/564/CH1/EX1.3/1_03.sce new file mode 100755 index 000000000..3f913c5f6 --- /dev/null +++ b/564/CH1/EX1.3/1_03.sce @@ -0,0 +1,27 @@ +pathname=get_absolute_file_path('1_03.sce')
+filename=pathname+filesep()+'1_03data.sci'
+exec(filename)
+Txy=((lLoad-Sx)*(lLoad-Sy))^0.5;
+printf("\nτxy:%f N/mm^2",Txy);
+printf(",%f N/mm^2",-Txy);
+S=poly(0,'S');
+y=(S^2)-(S*(Sx+Sy))+(Sx*Sy)-(Txy^2);
+m=roots(y);
+printf("\nσ1: %f N/mm^2",m(2,:));
+printf("\nσ2:%f N/mm^2",m(1,:));
+Tmax=abs((m(2,:)-m(1,:))/2);
+printf("\nτmax: %f N/mm^2",Tmax);
+//plotting mohr circle
+x=abs((m(2,:)+m(1,:))/2);//Centre of circle
+plot2d(0,0,-1,"031"," ",[x-Tmax-50,-Tmax-50,x+Tmax+100,Tmax+50]);
+xgrid(3);
+xarc(x-Tmax,Tmax,2*Tmax,2*Tmax,0,360*64);
+x1=[m(1,:),m(2,:)],y1=[0,0];
+x2=[Sx,Sy],y2=[Txy,-Txy];
+x3=[Sx,Sx],y3=[Txy,0];
+x4=[Sy,Sy],y4=[-Txy,0];
+x5=[x,x],y5=[Tmax,-Tmax];
+plot(x1,y1,x2,y2,'--',x3,y3,'--',x4,y4,'--',x5,y5,'--');
+datatipToggle();
+printf("\n\nclick on the point to view its coordinate on the plot");
+xtitle( 'Mohr’s circle of stress', ' σ (N/mm^2)', ' τ (N/mm^2)', boxed = 1 );
\ No newline at end of file diff --git a/564/CH1/EX1.3/1_03MohrCircle.JPG b/564/CH1/EX1.3/1_03MohrCircle.JPG Binary files differnew file mode 100755 index 000000000..101cd4dca --- /dev/null +++ b/564/CH1/EX1.3/1_03MohrCircle.JPG |