blob: 73f126afbab482a3c1343c5877888c547ec5c4fa (
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
36
37
38
39
40
41
42
43
44
45
46
47
48
|
clc
ti=0.075 //m
to=0.015//m
a=0.05//m
b=0.25//m
delta=0.05 //mm
w=6900*(2*%pi/60) //rpm
s=1
row=7.8*10^3//Ns^2/m^4
E=200 //GPa
//solution a:
t1=ti*a^s
disp(t1,"t1 is=")
t1=to*b^2
disp(t1,"t1 is=")
//(ti/to)=(t1*a^-s)/(t1*b^-s)=(b/a)^s
c=(b/a)^s
(ti/to)==c
disp(ti/to,"ti/t0 is=")
m1=-0.5+sqrt((0.5)^2+(1+0.3*1))
disp(m1,"m1 is=")
m2=-0.5-sqrt((0.5)^2+(1+0.3*1))
disp(m2,"m2 is=")
//sigmar=0=(c1/t1)*(0.05)^m1+(c2/t1)*(0.05)^(m2)-0.00176*row*w^2 // r=0.05
//sigmar=0=(c1/t1)*(0.25)^m1+(c2/t1)*(0.25)^(m2)-0.0439*row*w^2 // r=0.25
c1=t1*0.12529*row*w^2
disp(c1,"c1 is=")
c2=t1*-6.272*10^-5*row*w^2
disp(c2,"c2 is=")
r=0.05
sigmar=(0.12529*r^0.745-6.272*10^-5*r^(-1.745)-0.70*r^2)//*row*w^2
disp(sigmar,"sigmar is= ")
sigmatheta=(0.09334*r^0.745+1.095*10^-4*r^(-1.745)-0.40*r^2)//*row*w^2
disp(sigmatheta,"sigmatheta is= ")
//solution b:
r=0.05
//ur=(r*sigmatheta)/E
ur=(r*sigmatheta)
disp(ur,"ur is= ")
|
