blob: 1e76238210c910f6c1eaeea4573ab792ddcf6bd8 (
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
|
clear
//
//
//
//Variable declaration
h1=1
h2=1
k1=1
k2=1
l1=1
l2=1
l3=0
s=3.5*10**6 //stress(Pa)
//Calculation
x=sqrt(h1**2+k1**2+l1**2)
y=sqrt(h2**2+k2**2+l2**2)
z=sqrt(h2**2+k2**2+l3**2)
cos_phi=((h1*h2)-(k1*k2)+(l1*l2))/(x*y)
sin_phi=sqrt(1-(cos_phi)**2)
cos_theta=((h1*h2)+(k1*k2)+(l1*l3))/(x*z)
ss=s*cos_theta*cos_phi*sin_phi //critical resolved shear stress(Pa)
//Result
printf("\n critical resolved shear stress is %0.3f MPa",ss/10**6)
|
