summaryrefslogtreecommitdiff
path: root/1898/CH5/EX5.21/Ex5_21.sce
blob: e30573bbe1735c97f17f9673cb47efa6f6b798fa (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
clear all; clc;

disp("Scilab Code Ex 5.21 : ")

//Given:
l = 1.5; //m
G = 42*10^3; //GPa
co = 50; //mm
ci = 25;//mm
shear_y = 84; //N/mm^2
strain_y = 0.002; //rad

//Plastic Torque:
T_p = ((2*%pi)*(co^3 - ci^3)*shear_y)/3;
phi_p = (strain_y*l*10^3)/ci;

J = (%pi/2)*(co^4 - ci^4);
shear_r = (T_p*co)/J;
shear_i = (shear_r*ci)/(co);// shear = Tc/J

G = shear_y/strain_y; 

phi_dash = (T_p*l*10^3)/(J*G); //phi = TpL/JG;

phi = phi_p - phi_dash;
T_p = T_p/10^6;

//Display:


printf('\n\nThe plastic torque Tp                                =  %1.2f x 10^6 Nmm',T_p);
printf('\nThe permanent twist of the tube if Tp is removed     = %1.5f rad',phi);


//----------------------------------------------------------------END------------------------------------------------------------------------