summaryrefslogtreecommitdiff
path: root/3764/CH3/EX3.09/Ex3_09.sce
blob: 97ad0ba15315b4087d110baea77a2bdfd4297665 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
clc
//

//Variable declaration
T=4.60*(10**3)                                                // Torque(N.m)
L=1.2                                                         // length(m) 
G=77*(10**9)                                                  // modulus of rigidity(Pa)
J=614*(10**-9)                                                // Polar moment of inertia(m**4)
phy=8.50  
c=25*(10**-3)                                                 // radius(m)

//Calculation         
// Case(a)
phyl=((T*L)/(J*G))*(360/(2*(%pi)))                        // Lateral twist(degree)
phyp=phy-phyl                                                 // Permanent twist(degree)
// Case(b)
Tlmax=((T*c)/(J))/((10**6))                                 // Residual stresses(MPa) 

// Result
printf("\n Case(a): Permanent twist = %1f degree' ,phyp)
printf("\n Case(b): Residual stress = %1f MPa ' ,Tlmax)