blob: bc96d913a60747cbe9e5c70407cd29fd36c17c20 (
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
|
// Given:-
omega = 100.0 //motor rotation speed in rad/s
tau = 18.0 //torque applied by shaft in N.m
Welecdot = -2.0 //electric power input in KW
Wshaftdot = (tau*omega)/1000 //shaft work rate in KW
Wdot = Welecdot + Wshaftdot //net work rate in KW
//function [Qdot]=f(t)
//Qdot = (-0.2)* [1-2**(-0.05*t)]
//function [Edot]=f1(t) //function for rate of change of energy
//Edot = (-0.2)*[1-2**(-0.05*t)] - Wdot
//function [deltaE] =f2(t) //function for change in energy
t = linspace(0,120,100);
for i = 1:100
Qd(i) = i
Wd(i) = i
dltaE(i) = i
Qd(i) = (-0.2*(1-%e^(-0.05*t(i))))
Wd(i) = Wdot
dltaE(i) = 4*(1 - %e^(-0.05*t(i)))
end
subplot(2,2,1)
plot(t,Qd)
xlabel("Time (s)")
ylabel("Qdot (KW)")
subplot(2,2,2)
plot(t,Wd)
xlabel("Time (s)")
ylabel("Wdot (KW)")
subplot(2,2,3)
plot(t,dltaE)
xlabel("Time (s)")
ylabel("deltaE (KJ)")
|
