| > | restart; |
| > | M0:=13:k:=13: nu:=40: mu:=12: |
| > | I1:=19: m2:=17: m3:=35: m4:=27: |
| > | R1:=0.36: r1:=0.25: R3:=0.26: |
| > | R4:=0.12: |
| > | I4:=m4*R4^2/2: w1:=0.3: f:=1.4: |
| > | I3:=m3*R3^2/2: w3:=w1*R1/R3: T3:=I3*w3^2/2: |
| > | T1:=I1*w1^2/2: |
| > | A:=(T1+T3)/(w1^2/2); |
| > | Va:=r1*w1: |
| > | T2:=m2*Va^2/2: B2:=T2/(w1^2/2); |
| > | T4:=3/4*m4*Va^2: B4:=T4/(w1^2/2); |
| > | B:=B2+B4; |
| > | w3:=w1*R1/R3:w4:=w1*sin(f)*r1/R4: |
| > | Q3:=(M0-k*w1);Qm:=-mu*w4^2/w1; |
| > | v5:=w1*2*r1*sin(f): |
| > | Qf:=-nu*v5^2/w1; |
| > | Q:=Qf+Qm+Q3; |
| > | eps:=(Q-0.5*B*sin(2*f)*w1^2)/(A+B*sin(f)^2); |