5 files changed, 87 insertions, 0 deletions
diff --git a/564/CH6/EX6.1/6_1.JPG b/564/CH6/EX6.1/6_1.JPG
new file mode 100755
index 000000000..3d3c44983
--- /dev/null
+++ b/564/CH6/EX6.1/6_1.JPG
diff --git a/564/CH6/EX6.1/6_1.sce b/564/CH6/EX6.1/6_1.sce
new file mode 100755
index 000000000..29664b24c
--- /dev/null
+++ b/564/CH6/EX6.1/6_1.sce
@@ -0,0 +1,41 @@
+pathname=get_absolute_file_path('6_1.sce')
+filename=pathname+filesep()+'6_1data.sci'
+exec(filename)
+//theta=θ,lambda=λ,mu=μ
+theta12=0,lambda12=cos(theta12),mu12=sin(theta12);
+theta13=%pi/2,lambda13=cos(theta13),mu13=sin(theta13);
+theta23=%pi/2 + atan(L13/L12);lambda23=cos(theta23),mu23=sin(theta23);
+L23=(L12^2 +L13^2)^0.5;
+LL12=lambda12^2,LM12=lambda12*mu12,MM12=mu12^2;
+LL13=lambda13^2,LM13=lambda13*mu13,MM13=mu13^2;
+LL23=lambda23^2,LM23=lambda23*mu23,MM23=mu23^2;
+K220=[LL12 LM12;LM12 MM12];
+K222=(L12/L23)*[LL23 LM23;LM23 MM23];
+K210=[-LL12 -LM12;-LM12 -MM12];
+K232=(L12/L23)*[-LL23 -LM23;-LM23 -MM23];
+K=K220+K222;
+F=[FX2;FY2];
+V=(L12/AE)*inv(K)*F;
+u2=V(1,:),v2=V(2,:);
+u1=0;v1=0;u3=0;v3=0;
+F1=[K210';K232']*inv(K)*F;
+S12=(AE/L12)*[lambda12 mu12]*[u2-u1;v2-v1];
+S13=(AE/L13)*[lambda13 mu13]*[u3-u1;v3-v1];
+S23=(AE/L23)*[lambda23 mu23]*[u3-u2;v3-v2];
+x1=[-L12,-L12],y1=[0,L13];
+x2=[-L12,u2],y2=[0,v2];
+x3=[-L12,u2],y3=[L13,v2];
+plot(x1,y1,x2,y2,x3,y3);
+xgrid(3);
+datatipToggle();
+xtitle( 'framework after deformation', ' -x- ', ' -y- ',boxed = 1);
+printf("\n\nu2: %f m",u2);
+printf("\nv2: %f m",v2);
+printf("\nFx,1: %f N",F1(1,:));
+printf("\nFy,1: %f N",F1(2,:));
+printf("\nFx,3: %f N",F1(3,:));
+printf("\nFy,3: %f N",F1(4,:));
+printf("\nS12: %f N",S12);
+printf("\nS13: %f N",S13);
+printf("\nS23: %f N",S23);
+printf("\n\nclick on the point to view its coordinate on the plot");
+\ No newline at end of file
diff --git a/564/CH6/EX6.2/6_2.JPG b/564/CH6/EX6.2/6_2.JPG
new file mode 100755
index 000000000..8ef98acb3
--- /dev/null
+++ b/564/CH6/EX6.2/6_2.JPG
diff --git a/564/CH6/EX6.2/6_2.sce b/564/CH6/EX6.2/6_2.sce
new file mode 100755
index 000000000..aa9952125
--- /dev/null
+++ b/564/CH6/EX6.2/6_2.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('6_2.sce')
+filename=pathname+filesep()+'6_2data.sci'
+exec(filename)
+F=[-W;M;0;0];
+P=EI*[(27/(2*L^3)) (9/(2*L^2)) (6/L^2) -(3/(2*L^2));
+   (9/(2*L^2)) (6/L) (2/L) (1/L);
+   (6/L^2) (2/L) (4/L) 0;
+   -(3/(2*L^2)) (1/L) 0 (2/L)];
+V=inv(P)*F;
+Sy12= EI*(-(6*V(3)/(L^2))-(12*V(1)/(L^3)) -(6*V(2)/(L*L)));//Sy,12
+deff("[M12]=f(x)","M12=EI*(((-6*x/L*L)+(4/L))*V(3) + ((-12*x/L^3)+(6/L^2))*V(1) +((-6*x/L*L)+(2/L))*V(2))");
+x=[0:0.05:L];
+fplot2d(x,f);
+xgrid(3);
+xtitle( 'M12 versus x', ' -x- ', 'M12');
+datatipToggle();
+printf("\nυ2: %f m",V(1));//v2
+printf("\nθ2: %f ",V(2));//θ2
+printf("\nθ2: %f ",V(3));//θ1
+printf("\nθ2: %f ",V(4));//θ3
+printf("\nSy,12: %f",Sy12);//Sy,12
+\ No newline at end of file
diff --git a/564/CH6/EX6.4/6_4.sce b/564/CH6/EX6.4/6_4.sce
new file mode 100755
index 000000000..b40aa1f33
--- /dev/null
+++ b/564/CH6/EX6.4/6_4.sce
@@ -0,0 +1,25 @@
+pathname=get_absolute_file_path('6_4.sce')
+filename=pathname+filesep()+'6_4data.sci'
+exec(filename)
+P=[1 P1(1) P1(2) P1(1)*P1(2);
+   1 P2(1) P2(2) P2(1)*P2(2);
+   1 P3(1) P3(2) P3(1)*P3(2);
+   1 P4(1) P4(2) P4(1)*P4(2)];
+alpha1=inv(P)*u;
+alpha2=inv(P)*v;
+alpha=[alpha1;alpha2];
+deff("[Ex]=f(y)","Ex=alpha(2)+ y*alpha(4)");//εx
+deff("[Ey]=f1(x)","Ey=alpha(7)+ x*alpha(8)");//εy
+function[G]=Gxy(x,y)//γxy
+    G=x*alpha(4) +y*alpha(8) +alpha(3)+alpha(6);
+endfunction
+//at the centre
+Pc(1)=(P1(1)+P3(1))/2;
+Pc(2)=(P1(2)+P3(2))/2;
+Sy=(E/(1-V^2))*(f1(Pc(2)) +V*f(Pc(1)));
+Sx=(E/(1-V^2))*(f(Pc(1)) +V*f1(Pc(2)));
+gxy=Gxy(0,0);
+txy=(E/(1-V^2))*0.5*(1-V)*gxy;//τxy
+printf("\nσx: %f N/mm^2",Sx);
+printf("\nσy: %f N/mm^2",Sy);
+printf("\nτxy: %f N/mm^2",txy);
+\ No newline at end of file