blob: 59de3cf68970129ce323262a613041a67d5a0060 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
clear
clc
Pm=2
Pi=1
H=6
G=1
f=50
p=Pi/Pm
M=G*H/(%pi*f)
d0=asin(p)
dcc=acos(((p*(%pi - (2*d0)))- (Pi*cos(d0)))/(Pm-Pi))
mprintf("Critical Clearing angle = %.4f rad\n\n", dcc)
tcc=sqrt(2*M*(dcc-d0)/Pi)
mprintf("Critical Clearing time = %.3f sec = %.2f cycles", tcc , tcc*50)
|
