summaryrefslogtreecommitdiff
path: root/3863/CH3/EX3.9/Ex3_9.sce
blob: 8aacb78f16305e19061263314c4786e2c56be01b (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
clear
//

//Given
//Variable declaration
sigma1=100        //Major principal stress in N/sq.mm
sigma2=-40        //Minor principal stress in N/sq.mm
theta=90-60       //Angle of inclination in degrees

//Calculation
sigman=((sigma1+sigma2)/2)+(((sigma1-sigma2)/2)*cos((%pi/180)*(2*theta)))
sigmat=((sigma1-sigma2)/2*(sin((%pi/180)*(2*theta))))

sigmaR=(sqrt(sigman**2+sigmat**2))

sigmat_max=int((sigma1-sigma2)/2)
phi=int((180/%pi)*(atan(sigmat/sigman)))

//Result
printf("\n Resultant stress in magnitude = %0.3f  N/mm^2",sigmaR)
printf("\n Direction of resultant stress = %0.3f  degrees",phi)
printf("\n Maximum shear stress = %0.3f  N/mm^2",sigmat_max)