blob: b0d6586e5b1caef56f6c358e75e8f17a0691c2df (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
syms f J K s t
T=1/(J*(s^2+(f/J)*s+(K/J))) //Q/I
omegaN=sqrt(K/J)
Mp=6 //given
zeta=0.667
omegaD=omegaN*sqrt(1-zeta^2)
Tp=%pi/omegaD
disp(Tp," Tp = ")
I=laplace('10',t,s)
Q=I*T
x=limit(s*Q,s,0);
disp(10/0.5," K = ")
disp(K/omegaN^2," J= ")
disp(zeta*(2*sqrt(K*J))," f = ")
|
