blob: 04a7891f618ec7739cf28c10d41205072353e8b4 (
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
|
clear
//
//variable declaration
//Let the principal plane make anticlockwise angle theta with the plane of px with y-axis. Then
px=(-100) //N/mm^2
py=(-75) //N/mm^2
q=(-50) //N/mm^2
p1=(px+py)/2+sqrt((((px-py)/2)**2)+(q**2))
printf("\n p1= %0.2f N/mm^2",p1)
p2=(px+py)/2-sqrt((((px-py)/2)**2)+(q**2))
printf("\n p2= %0.2f N/mm^2",p2)
qmax=sqrt((((px-py)/2**2))+(q**2))
printf("\n qmax= %0.2f N/mm^2",qmax)
//let theta be the inclination of principal stress to the plane of px.
theta1=(atan((2*q)/(px-py))*180)/(%pi*2)
theta2=90+theta1
printf("\n theta= %0.2f ° and %0.2f °",theta1,theta2)
|
