initial commit / add all books

8 files changed, 129 insertions, 0 deletions

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
--- /dev/null
+++ b/564/CH1/EX1.3/1_03MohrCircle.JPG
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
--- /dev/null
+++ b/564/CH1/EX1.5/1_05mohrcircle.JPG
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