From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 564/CH1/EX1.1/1_01.sce           |  18 ++++++++++++++++++
 564/CH1/EX1.2/1_02.sce           |  18 ++++++++++++++++++
 564/CH1/EX1.3/1_03.sce           |  27 +++++++++++++++++++++++++++
 564/CH1/EX1.3/1_03MohrCircle.JPG | Bin 0 -> 72564 bytes
 564/CH1/EX1.4/1_04.sce           |  13 +++++++++++++
 564/CH1/EX1.5/1_05.sce           |  34 ++++++++++++++++++++++++++++++++++
 564/CH1/EX1.5/1_05mohrcircle.JPG | Bin 0 -> 77772 bytes
 564/CH1/EX1.7/1_07.sce           |  19 +++++++++++++++++++
 8 files changed, 129 insertions(+)
 create mode 100755 564/CH1/EX1.1/1_01.sce
 create mode 100755 564/CH1/EX1.2/1_02.sce
 create mode 100755 564/CH1/EX1.3/1_03.sce
 create mode 100755 564/CH1/EX1.3/1_03MohrCircle.JPG
 create mode 100755 564/CH1/EX1.4/1_04.sce
 create mode 100755 564/CH1/EX1.5/1_05.sce
 create mode 100755 564/CH1/EX1.5/1_05mohrcircle.JPG
 create mode 100755 564/CH1/EX1.7/1_07.sce

(limited to '564/CH1')

diff --git a/564/CH1/EX1.1/1_01.sce b/564/CH1/EX1.1/1_01.sce
new file mode 100755
index 000000000..7907f07d6
--- /dev/null
+++ b/564/CH1/EX1.1/1_01.sce
@@ -0,0 +1,18 @@
+pathname=get_absolute_file_path('1_01.sce')
+filename=pathname+filesep()+'1_01data.sci'
+exec(filename)
+Sx1= p*d/(4*t);
+Sy= p*d/(2*t);//σy
+printf("\nLongitudnal stress produced by internal pressure: %f N/mm^2",Sx1);
+printf("\nCircumferential stress produced by internal pressure (σy): %f N/mm^2",Sy);
+Sx2= Load/(%pi*d*t);
+printf("\ndirect stress due to the axial load: %f N/mm^2",Sx2);
+Sx= Sx1+Sx2;//σx
+printf("\nσx is: %f N/mm^2",Sx);
+printf("\nσy is: %f N/mm^2",Sy);
+Sn=(Sx)*(cos(theta))^2 + Sy*(sin(theta))^2;//σn
+T= (Sx-Sy)*sin(2*(theta)) /2;//τ
+printf("\nσn is: %f N/mm^2",Sn);
+printf("\nτ is: %f N/mm^2",T);
+Tmax= (Sx-Sy)/2;//τ max
+printf("\nτmax: %f N/mm^2",Tmax);
\ No newline at end of file
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
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
new file mode 100755
index 000000000..101cd4dca
Binary files /dev/null and b/564/CH1/EX1.3/1_03MohrCircle.JPG differ
diff --git a/564/CH1/EX1.4/1_04.sce b/564/CH1/EX1.4/1_04.sce
new file mode 100755
index 000000000..7100dce78
--- /dev/null
+++ b/564/CH1/EX1.4/1_04.sce
@@ -0,0 +1,13 @@
+pathname=get_absolute_file_path('1_04.sce')
+filename=pathname+filesep()+'1_04data.sci'
+exec(filename)
+Ex= (1/E)*(Sx-v*Sy);//εx
+Ey= (1/E)*(Sy-v*Sx);//εy
+Ez= (1/E)*(0-v*(Sx+Sy));//εz
+printf("\nεx: %f",Ex);
+printf("\nεy: %f",Ey);
+printf("\nεz: %f",Ez);
+Tmax=(Sx-Sy)/2;//τmax
+Gmax=2*(1+v)*Tmax/E;//γmax
+printf("\nγmax: %f",Gmax);
+printf("\nτmax: %f N/mm^2",Tmax);
\ No newline at end of file
diff --git a/564/CH1/EX1.5/1_05.sce b/564/CH1/EX1.5/1_05.sce
new file mode 100755
index 000000000..f44bb5b3f
--- /dev/null
+++ b/564/CH1/EX1.5/1_05.sce
@@ -0,0 +1,34 @@
+pathname=get_absolute_file_path('1_05.sce')
+filename=pathname+filesep()+'1_05data.sci'
+exec(filename)
+Ex= (1/E)*(Sx-v*Sy);//εx
+Ey= (1/E)*(Sy-v*Sx);//εy
+G=E/(2*(1+v));//shear modulus
+Gxy=Txy/G;//γxy
+printf("\n εx:  %f",Ex);
+printf("\n εy: %f",Ey);
+printf("\n G: %f",G);
+printf("\nγxy:  %f",Gxy);
+PE1=(Ex+Ey)/2 + 0.5*((Ex-Ey)^2 +Gxy^2)^0.5;//εI
+PE2=(Ex+Ey)/2 - 0.5*((Ex-Ey)^2 +Gxy^2)^0.5;//εII
+theta=0.5*atan(Gxy/(Ex-Ey)) *(180/%pi);//θ in degree
+printf("\n εI:  %f",PE1);
+printf("\nεII: %f",PE2);
+printf("\n  θ:  %f degree",theta);
+printf("\n  θ:  %f degree",theta+90);
+//plotting mohr circle
+x=abs((PE1+PE2)/2);//Centre
+r=(abs(PE1)+abs(PE2))/2;//radius
+plot2d(0,0,-1,"031"," ",[x- 1.5*r,-1.5*r,x+ 1.5*r,1.5*r]);
+xgrid(3);
+xarc(x-r,r,2*r,2*r,0,360*64);
+x1=[PE1,PE2],y1=[0,0];
+x2=[Ex,Ey],y2=0.5*[Gxy,-Gxy];
+theta1=atan(Gxy/(2*Ex)) *(180/%pi);
+plot(x1,y1);
+plot(x2,y2);
+xarc(x-r/4,r/4,r/2,r/2,0,theta1*64);
+xarc(x-r/2,r/2,r,r,180*64,(180+theta1)*64);
+datatipToggle();
+printf("\n\nclick on the point to view its coordinate on the plot");
+xtitle( 'Mohr’s circle of strain', ' ε', ' γ (gamma)', boxed = 1 );
\ No newline at end of file
diff --git a/564/CH1/EX1.5/1_05mohrcircle.JPG b/564/CH1/EX1.5/1_05mohrcircle.JPG
new file mode 100755
index 000000000..7f24c7de5
Binary files /dev/null and b/564/CH1/EX1.5/1_05mohrcircle.JPG differ
diff --git a/564/CH1/EX1.7/1_07.sce b/564/CH1/EX1.7/1_07.sce
new file mode 100755
index 000000000..52a9a57f1
--- /dev/null
+++ b/564/CH1/EX1.7/1_07.sce
@@ -0,0 +1,19 @@
+pathname=get_absolute_file_path('1_07.sce')
+filename=pathname+filesep()+'1_07data.sci'
+exec(filename)
+E1=0.5*(Ea+Ec) + (((Ea-Eb)^2 +(Ec-Eb)^2)/2)^0.5;
+E2=0.5*(Ea+Ec) - (((Ea-Eb)^2 +(Ec-Eb)^2)/2)^0.5;
+disp(E1,"εI=",E2,"εII=");
+S1=E*(E1+ v*E2)/(1-v^2);//σI
+S2=E*(E2+ v*E1)/(1-v^2);//σII
+Sx=S1+S2;//σx
+printf("\nσI: %f N/mm^2",S1);
+printf("\nσII: %f N/mm^2",S2);
+printf("\nσx: %f N/mm^2",Sx);
+P=Sx*(%pi*d^2)/4;//axial tensile load
+Txy=(((S1-S2)^2 -(S1+S2)^2)/4)^0.5;//τxy
+printf("\nτxy: %f N/mm^2",Txy);
+J=(%pi*d^4)/32;//torsion constant
+T=2*Txy*J/d;//Torque
+disp(p/10^3,"P in KN",P,"P(N)=");
+disp(T/10^6,"T in KN.m",T,"T (N.mm)=");
\ No newline at end of file
-- 
cgit