initial commit / add all books

6 files changed, 82 insertions, 0 deletions

diff --git a/564/CH26/EX26.1/26_1.sce b/564/CH26/EX26.1/26_1.sce
new file mode 100755
index 000000000..5975e6471
--- /dev/null
+++ b/564/CH26/EX26.1/26_1.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('26_1.sce')

+filename=pathname+filesep()+'26_1data.sci'

+exec(filename)

+clear

+L41=L23*cos(asin((L12-L34)/L23));

+t41=t23;

+C=[0;cos(asin((L12-L34)/L23));0;1];

+S=[1;sin(asin((L12-L34)/L23));1;0];

+P=[1 0 0 0 0 -t12*G 0;

+   0 1 0 0 t23*G*cos(asin((L12-L34)/L23)) t23*G*sin(asin((L12-L34)/L23)) -t23*G*L12*cos(asin((L12-L34)/L23));

+   0 0 1 0 0 t34*G -t34*G*L23*cos(asin((L12-L34)/L23));

+   0 0 0 1 -t41*G 0 0;

+   0 1 0 -1 0 0 0;

+   L12 -L23*sin(asin((L12-L34)/L23)) -L34 0 0 0 0;

+   0 L23*L12*cos(asin((L12-L34)/L23)) L34*L23*cos(asin((L12-L34)/L23)) 0 0 0 0];

+X=[0;0;0;0;0;Load;Load*x];

+q=inv(P)*X;

+printf("\nτ12: %f N/mm^2",q(1)/t12);

+printf("\nτ23: %f N/mm^2",q(2)/t23);

+printf("\nτ34: %f N/mm^2",q(3)/t34);

+printf("\nτ41: %f N/mm^2",q(4)/t41);
\ No newline at end of file
diff --git a/564/CH26/EX26.2/26_2.JPG b/564/CH26/EX26.2/26_2.JPG
new file mode 100755
index 000000000..6d7566002
--- /dev/null
+++ b/564/CH26/EX26.2/26_2.JPG
diff --git a/564/CH26/EX26.2/26_2.sce b/564/CH26/EX26.2/26_2.sce
new file mode 100755
index 000000000..d3955550a
--- /dev/null
+++ b/564/CH26/EX26.2/26_2.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('26_2.sce')

+filename=pathname+filesep()+'26_2data.sci'

+exec(filename)

+clear

+du=0;

+P=[1 0 0 G*t*R;

+   0 1 -G*t 2*R*G*t;

+   0 0 1.79 -R

+   0 0 -0.13 R];

+X=[0;0;Load/(2*R*G*t);Load/((-0.5*%pi -6)*G*R*t)];

+q=inv(P)*X;

+printf("\nq12 = q34 = %f N/mm",q(1));

+printf("\nq41= %f N/mm",q(2));

+function[q23]=f(phi)

+    q23=-R*G*t*q(4) -G*t*sin(phi)*q(3);

+endfunction

+phi=linspace(0,%pi,%pi*10);

+q23=feval(phi,f);

+plot2d(phi,q23);

+xgrid(3),xtitle('q23 versus φ','φ','q23');

+datatipToggle();
\ No newline at end of file
diff --git a/564/CH26/EX26.3/26_3.sce b/564/CH26/EX26.3/26_3.sce
new file mode 100755
index 000000000..7d8f7d373
--- /dev/null
+++ b/564/CH26/EX26.3/26_3.sce
@@ -0,0 +1,19 @@
+pathname=get_absolute_file_path('26_3.sce')

+filename=pathname+filesep()+'26_3data.sci'

+exec(filename)

+clear

+Tr=T*L1/2;

+mu=(8*G*t/(A*E*(b+a)))^0.5;

+L=L1/2;

+k1=((T*(b-a)*10^3)/(8*a*b*G*t));

+k2=1/(mu*cosh(mu*L));

+k3=((4*(b-a))/(a*b*(b+a)));

+k4=(2*T*(10^3)/(a*b*G*t*(b+a)));

+function[th]=f(z)

+    w=(k1*((k2*sinh(mu*z))-z));

+    F=((k1*k3 +k4)*L*L*0.5 -(k1*k2*k3/mu)*cosh(mu*L));

+    th=(k1*k2*k3/mu)*cosh(mu*z) -(k1*k3 +k4)*z*z*0.5 + F;

+endfunction

+funcprot();

+printf("\nangle of twist at mid-span θ: %f rad",f(0));

+

diff --git a/564/CH26/EX26.4/26_4.JPG b/564/CH26/EX26.4/26_4.JPG
new file mode 100755
index 000000000..9babb9f38
--- /dev/null
+++ b/564/CH26/EX26.4/26_4.JPG
diff --git a/564/CH26/EX26.4/26_4.sce b/564/CH26/EX26.4/26_4.sce
new file mode 100755
index 000000000..e0a298ead
--- /dev/null
+++ b/564/CH26/EX26.4/26_4.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('26_4.sce')

+filename=pathname+filesep()+'26_4data.sci'

+exec(filename)

+B=((a*t2)/6)+((b*t1)/18)*3;

+A=((b*t1)/18)*6;

+L1=Load/4;

+x=L1/a;

+d=b/3;

+E=1,G=GE;

+mu=(((G*t1)/(d*E))*((A+B)/(A*B)))^0.5;

+D=(-x*A)/((A+B)*mu*cosh(L*0.5*mu));

+function[Sb]=f(z)

+     Sb=((D/B)*sinh(mu*z)) -(x/(A+B))*z;

+endfunction

+funcprot()

+z=linspace(0,L,10*L);

+Sb=feval(z,f);

+plot2d(z,Sb);

+xgrid(3),xtitle('σB versus z','z','σB');

+datatipToggle()

+printf("\nclick on the point to view its coordinate on the plot");
\ No newline at end of file