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 --- 866/CH14/EX14.1/14_1.sce | 19 +++++++++++++++++++ 866/CH14/EX14.10/14_10.sce | 13 +++++++++++++ 866/CH14/EX14.2/14_2.sce | 19 +++++++++++++++++++ 866/CH14/EX14.3/14_3.sce | 24 ++++++++++++++++++++++++ 866/CH14/EX14.5/14_5.sce | 24 ++++++++++++++++++++++++ 866/CH14/EX14.6/14_6.sce | 21 +++++++++++++++++++++ 866/CH14/EX14.7/14_7.sce | 21 +++++++++++++++++++++ 866/CH14/EX14.8/14_8.sce | 19 +++++++++++++++++++ 866/CH14/EX14.9/14_9.sce | 15 +++++++++++++++ 9 files changed, 175 insertions(+) create mode 100755 866/CH14/EX14.1/14_1.sce create mode 100755 866/CH14/EX14.10/14_10.sce create mode 100755 866/CH14/EX14.2/14_2.sce create mode 100755 866/CH14/EX14.3/14_3.sce create mode 100755 866/CH14/EX14.5/14_5.sce create mode 100755 866/CH14/EX14.6/14_6.sce create mode 100755 866/CH14/EX14.7/14_7.sce create mode 100755 866/CH14/EX14.8/14_8.sce create mode 100755 866/CH14/EX14.9/14_9.sce (limited to '866/CH14') diff --git a/866/CH14/EX14.1/14_1.sce b/866/CH14/EX14.1/14_1.sce new file mode 100755 index 000000000..7ead5d32e --- /dev/null +++ b/866/CH14/EX14.1/14_1.sce @@ -0,0 +1,19 @@ +clc +//initialisation of variables +d= 2 //m +t= 20 //mm +p= 1.5 //N/mm^2 +load= 2500 //KN +alpha= 60 //degrees +//CALCULATIONS +Cs= (p*d*10^3)/(2*t) +Ls= (p*d*10^3)/(4*t) +Ds= (load*10^3)/(%pi*d*t*10^3) +Ts= Ls+Ds +sigman= (Ts*t*(cosd(90-alpha))^2+Cs*t*(cosd(alpha))^2)/t +tab= Ts*sind(alpha)*cosd(alpha)-Cs*sind(alpha)*cosd(alpha) +tmax= (Ts-Cs)/2 +//RESULTS +printf ('direct stress= %.1f N/mm^2',sigman) +printf (' \n Shear stress=%.1f N/mm^2',tab) +printf (' \n maximum Shear stress=%.1f N/mm^2',tmax) diff --git a/866/CH14/EX14.10/14_10.sce b/866/CH14/EX14.10/14_10.sce new file mode 100755 index 000000000..9358e82f0 --- /dev/null +++ b/866/CH14/EX14.10/14_10.sce @@ -0,0 +1,13 @@ +clc +//initialisation of variables +BeamL= 250 //mm +BeamH= 500 //mm +spanlength= 4 //m +Failurestress= 1.5 //N/mm^2 +//CALCULATIONS +W1= (Failurestress*BeamL*BeamH^3)/(10^3*BeamL*12) +W2= (4*Failurestress*BeamL*BeamH)/3 +if(W1>W2) + printf ("W2= %.2e kN",W2) +else + printf ("W2= %.2e kN",W1) diff --git a/866/CH14/EX14.2/14_2.sce b/866/CH14/EX14.2/14_2.sce new file mode 100755 index 000000000..12ff51aff --- /dev/null +++ b/866/CH14/EX14.2/14_2.sce @@ -0,0 +1,19 @@ +clc +//initialisation of variables +load= 50000 //N +torque= 1200 //Nm +d= 60//mm +t= 1.5 //mm +alpha= 60 //degrees +//CALCULATIONS +BM= load*t +axialload= (load*4)/(%pi*d^2) +bendingmoment= (BM*d*64)/(%pi*d^4*2) +Ts= axialload+bendingmoment +shearstress= (torque*10^3*d*32)/(2*%pi*d^4) +sigman= -Ts*(cosd(alpha-30))^2+shearstress*cosd(alpha-30)*sind(alpha-30)+shearstress*cosd(alpha-30)*sind(alpha-30) +T= -Ts*sind(alpha)*cosd(alpha)-shearstress*(sind(alpha))^2+shearstress*(cosd(alpha)^2) +//RESULTS +printf ('direct stress= %.1f N/mm^2',sigman) +printf (' \n Shear stress=%.1f N/mm^2',T) + diff --git a/866/CH14/EX14.3/14_3.sce b/866/CH14/EX14.3/14_3.sce new file mode 100755 index 000000000..af9f5eecd --- /dev/null +++ b/866/CH14/EX14.3/14_3.sce @@ -0,0 +1,24 @@ +clc +//initialisation of variables +directstress= 80 //N.mm^2 +shearstress= 45 //N/mm^2 +//CALCULATIONS +Sx= directstress +Sy= 0 +Txy= shearstress +sigma1= ((Sx+Sy)/2)+(sqrt((Sx-Sy)^2+4*Txy^2))/2 +sigma2= ((Sx+Sy)/2)-(sqrt((Sx-Sy)^2+4*Txy^2))/2 +theta= (atand(-((2*Txy)/(Sx-Sy))))/2 +direction1= theta +direction2= theta-90 +Tmax= (sigma1-sigma2) /2 +theta1= direction1-45 +theta2= direction2-45 +//RESULTS +printf ('principal stress1= %.1f N/mm^2',sigma1) +printf (' \n principal stress2=%.1f N/mm^2',sigma2) +printf (' \n direction of principal stress1=%.1f degrees',direction1) +printf (' \n.direction of principal stress2=%.1f degrees',direction2) +printf (' \n.maximum shear stress=%.1f N/mm^2',Tmax) +printf (' \n.derection of plane1=%.1f degrees',theta1) +printf (' \n.derection of plane2=%.1f degrees',theta2) diff --git a/866/CH14/EX14.5/14_5.sce b/866/CH14/EX14.5/14_5.sce new file mode 100755 index 000000000..d2abc4e3b --- /dev/null +++ b/866/CH14/EX14.5/14_5.sce @@ -0,0 +1,24 @@ +clc +//initialisation of variables +sigmax= 60 //N/mm^2 +sigmay= -40 //N/mm^2 +Txy= 50 //N/mm^2 +E= 200000 //N/mm^2 +mu= 0.3 +//CALCULATIONS +Sx= (sigmax-mu*sigmay)/E +Sy= (sigmay-mu*sigmax)/E +G= E/(2*(1+mu)) +gamaxy= Txy/G +e1= ((Sx+Sy)/2+(sqrt((Sx-Sy)^2+(gamaxy)^2))/2) +e2= ((Sx+Sy)-(sqrt((Sx-Sy)^2+gamaxy^2)))/2 +theta= atand(-gamaxy/(Sx-Sy))/2 +theta1= -90+theta +//RESULTS +printf ('direct strain in x direction= %.2e ',Sx) +printf (' \n direct strain in y direction=%.2e ',Sy) +printf (' \n Gamaxy=%.2e ',gamaxy) +printf (' \n.principal strain1=%.9e ',e1) +printf (' \n.principal strain2=%.9e ',e2) +printf (' \n.direction of plane1=%.1f degrees',theta) +printf (' \n.direction of plane2=%.1f degrees',theta1) diff --git a/866/CH14/EX14.6/14_6.sce b/866/CH14/EX14.6/14_6.sce new file mode 100755 index 000000000..22f8b2afc --- /dev/null +++ b/866/CH14/EX14.6/14_6.sce @@ -0,0 +1,21 @@ +clc +//initialisation of variables +ea= 1000*10^-6 +eb= -200*10^-6 +ec= -300*10^-6 +sigmay= 0 //N/mm^2 +E= 70000 //N/mm^2 +mu= 0.3 +d= 50 //mm +//CALCULATIONS +e1= ((ea+ec)/2)+sqrt(((ea-eb)^2+(ec-eb)^2)/2) +e2= ((ea+ec)/2)-sqrt(((ea-eb)^2+(ec-eb)^2)/2) +sigma1= (E*(e1+mu*e2))/(1-mu^2) +sigma2= (E*(e2+mu*e1))/(1-mu^2) +sigmax= sigma1+sigma2 +Txy= sqrt((((sigma1*2)-(sigmax))^2-(sigmax^2)))/2 +P= sigmax*((%pi/4)*(d^2)) +T= Txy*((%pi/32)*(50^4)/25) +//RESULTS +printf ('T= %.2fNm',T) +printf ('\n P= %.2fN',P) diff --git a/866/CH14/EX14.7/14_7.sce b/866/CH14/EX14.7/14_7.sce new file mode 100755 index 000000000..32d194fe7 --- /dev/null +++ b/866/CH14/EX14.7/14_7.sce @@ -0,0 +1,21 @@ +clc +//initialisation of variables +sigmax= 140 //N/mm^2 +sigmay= -70 //N/mm^2 +Txy= 60 //N/mm^2 +T= 225 //N/mm^2 +//CALCULATIONS +sigma1= (sigmax+sigmay)/2+(sqrt((sigmax-sigmay)^2+4*(Txy^2)))/2 +sigma2= (sigmax+sigmay)/2-(sqrt((sigmax-sigmay)^2+4*(Txy^2)))/2 +sigmaY= sigma1-sigma2 +sigmaY1= sqrt(sigma1^2+sigma2^2-sigma1*sigma2) +if(sigmaY>T) + disp("Tresca theory failure has ocurred") +else + disp("Tresca theory failed") +end +if(sigmaY1