CURENT / andes

Python toolbox / library for power system transient dynamics simulation with symbolic modeling and numerical analysis 🔥
https://ltb.curent.org
Other
208 stars 108 forks source link

Parsing issue in dyr for HYGOV #310

Closed rodrigomha closed 2 years ago

rodrigomha commented 2 years ago

It seems that parsing dyr files with HYGOV is failing

The raw file I'm using is this one:

0,   100.00, 33, 0, 0, 60.00     / PSS(R)E 33 RAW created by rawd33  SUN, JAN 30 2022  20:47

     1,'BUS 1',  20,2,   1,   1,   1,1.03000,  27.0698,1.10000,0.90000,1.10000,0.90000
     2,'BUS 2',  20,2,   1,   1,   1,1.01000,  17.3055,1.10000,0.90000,1.10000,0.90000
     3,'BUS 3',  20,3,   2,   1,   1,1.03000,   0.0000,1.10000,0.90000,1.10000,0.90000
     4,'BUS 4',  20,2,   2,   1,   1,1.01000, -10.1917,1.10000,0.90000,1.10000,0.90000
     5,'BUS 5', 230,1,   1,   1,   1,1.00646,  20.6078,1.10000,0.90000,1.10000,0.90000
     6,'BUS 6', 230,1,   1,   1,   1,0.97813,  10.5233,1.10000,0.90000,1.10000,0.90000
     7,'BUS 7', 230,1,   1,   1,   1,0.96102,   2.1143,1.10000,0.90000,1.10000,0.90000
     8,'BUS 8', 230,1,   3,   1,   1,0.94862, -11.7551,1.10000,0.90000,1.10000,0.90000
     9,'BUS 9', 230,1,   2,   1,   1,0.97138, -25.3519,1.10000,0.90000,1.10000,0.90000
    10,'BUS 10', 230,1,   2,   1,   1,0.98347, -16.9369,1.10000,0.90000,1.10000,0.90000
    11,'BUS 11', 230,1,   2,   1,   1,1.00826,  -6.6270,1.10000,0.90000,1.10000,0.90000
0 / END OF BUS DATA, BEGIN LOAD DATA
     7,'1 ',1,   1,   1,   967.000,   100.000,     0.000,     0.000,     0.000,     0.000,   1,1,0
     9,'1 ',1,   2,   1,  1767.000,   100.000,     0.000,     0.000,     0.000,     0.000,   1,1,0
0 / END OF LOAD DATA, BEGIN FIXED SHUNT DATA
     7,'1 ',1,     0.000,   200.000
     9,'1 ',1,     0.000,   350.000
0 / END OF FIXED SHUNT DATA, BEGIN GENERATOR DATA
     1,'1 ',   700.000,   185.002,  9999.000, -9999.000,1.03000,     0,   900.000, 2.50000E-3, 2.50000E-1, 0.00000E+0, 0.00000E+0,1.00000,1,  100.0,  9999.000, -9999.000,   1,1.0000
     2,'1 ',   700.000,   234.578,  9999.000, -9999.000,1.01000,     0,   900.000, 2.50000E-3, 2.50000E-1, 0.00000E+0, 0.00000E+0,1.00000,1,  100.0,  9999.000, -9999.000,   1,1.0000
     3,'1 ',   719.083,   175.993,  9999.000, -9999.000,1.03000,     0,   900.000, 2.50000E-3, 2.50000E-1, 0.00000E+0, 0.00000E+0,1.00000,1,  100.0,  9999.000, -9999.000,   1,1.0000
     4,'1 ',   700.000,   202.038,  9999.000, -9999.000,1.01000,     0,   900.000, 2.50000E-3, 2.50000E-1, 0.00000E+0, 0.00000E+0,1.00000,1,  100.0,  9999.000, -9999.000,   1,1.0000
0 / END OF GENERATOR DATA, BEGIN BRANCH DATA
     5,     6,'1 ', 2.50000E-3, 2.50000E-2,   0.04375,    0.00,    0.00,    0.00,  0.00000,  0.00000,  0.00000,  0.00000,1,1,  25.00,   1,1.0000
     6,     7,'1 ', 1.00000E-3, 1.00000E-2,   0.01750,    0.00,    0.00,    0.00,  0.00000,  0.00000,  0.00000,  0.00000,1,1,  10.00,   1,1.0000
     7,     8,'1 ', 1.10000E-2, 1.10000E-1,   0.19250,    0.00,    0.00,    0.00,  0.00000,  0.00000,  0.00000,  0.00000,1,1, 110.00,   1,1.0000
     7,     8,'2 ', 1.10000E-2, 1.10000E-1,   0.19250,    0.00,    0.00,    0.00,  0.00000,  0.00000,  0.00000,  0.00000,1,1, 110.00,   1,1.0000
     8,     9,'1 ', 1.10000E-2, 1.10000E-1,   0.19250,    0.00,    0.00,    0.00,  0.00000,  0.00000,  0.00000,  0.00000,1,1, 110.00,   1,1.0000
     8,     9,'2 ', 1.10000E-2, 1.10000E-1,   0.19250,    0.00,    0.00,    0.00,  0.00000,  0.00000,  0.00000,  0.00000,1,1, 110.00,   1,1.0000
     9,    10,'1 ', 1.00000E-3, 1.00000E-2,   0.01750,    0.00,    0.00,    0.00,  0.00000,  0.00000,  0.00000,  0.00000,1,1,  10.00,   1,1.0000
    10,    11,'1 ', 2.50000E-3, 2.50000E-2,   0.04375,    0.00,    0.00,    0.00,  0.00000,  0.00000,  0.00000,  0.00000,1,1,  25.00,   1,1.0000
0 / END OF BRANCH DATA, BEGIN TRANSFORMER DATA
     1,     5,     0,'1 ',1,1,1, 0.00000E+0, 0.00000E+0,2,'TRFO1-5',1,   1,1.0000,   0,1.0000,   0,1.0000,   0,1.0000,'            '
 0.00000E+0, 1.66670E-2,   100.00
1.00000,   0.000,   0.000,   900.00,     0.00,     0.00, 0,      0, 1.10000, 0.90000, 1.10000, 0.90000,  33, 0, 0.00000, 0.00000,  0.000
1.00000,   0.000
     2,     6,     0,'1 ',1,1,1, 0.00000E+0, 0.00000E+0,2,'TRFO2-6',1,   1,1.0000,   0,1.0000,   0,1.0000,   0,1.0000,'            '
 0.00000E+0, 1.66670E-2,   100.00
1.00000,   0.000,   0.000,   900.00,     0.00,     0.00, 0,      0, 1.10000, 0.90000, 1.10000, 0.90000,  33, 0, 0.00000, 0.00000,  0.000
1.00000,   0.000
     3,    11,     0,'1 ',1,1,1, 0.00000E+0, 0.00000E+0,2,'TRFO3-11',1,   1,1.0000,   0,1.0000,   0,1.0000,   0,1.0000,'            '
 0.00000E+0, 1.66670E-2,   100.00
1.00000,   0.000,   0.000,   900.00,     0.00,     0.00, 0,      0, 1.10000, 0.90000, 1.10000, 0.90000,  33, 0, 0.00000, 0.00000,  0.000
1.00000,   0.000
     4,    10,     0,'1 ',1,1,1, 0.00000E+0, 0.00000E+0,2,'TRFO4-10',1,   1,1.0000,   0,1.0000,   0,1.0000,   0,1.0000,'            '
 0.00000E+0, 1.66670E-2,   100.00
1.00000,   0.000,   0.000,   900.00,     0.00,     0.00, 0,      0, 1.10000, 0.90000, 1.10000, 0.90000,  33, 0, 0.00000, 0.00000,  0.000
1.00000,   0.000
0 / END OF TRANSFORMER DATA, BEGIN AREA DATA
0 / END OF AREA DATA, BEGIN TWO-TERMINAL DC DATA
0 / END OF TWO-TERMINAL DC DATA, BEGIN VSC DC LINE DATA
0 / END OF VSC DC LINE DATA, BEGIN IMPEDANCE CORRECTION DATA
0 / END OF IMPEDANCE CORRECTION DATA, BEGIN MULTI-TERMINAL DC DATA
0 / END OF MULTI-TERMINAL DC DATA, BEGIN MULTI-SECTION LINE DATA
0 / END OF MULTI-SECTION LINE DATA, BEGIN ZONE DATA
0 / END OF ZONE DATA, BEGIN INTER-AREA TRANSFER DATA
0 / END OF INTER-AREA TRANSFER DATA, BEGIN OWNER DATA
0 / END OF OWNER DATA, BEGIN FACTS DEVICE DATA
0 / END OF FACTS DEVICE DATA, BEGIN SWITCHED SHUNT DATA
0 / END OF SWITCHED SHUNT DATA, BEGIN GNE DATA
0 / END OF GNE DATA, BEGIN INDUCTION MACHINE DATA
0 / END OF INDUCTION MACHINE DATA
Q

and dyr:

  1     'GENROU' 1    8    0.03   0.4  0.05  6.5  0  1.8  1.7   0.30  0.55  0.25    0.2  0.0      0.0 /   
  2     'GENROU' 1    8    0.03   0.4  0.05  6.5  0  1.8  1.7   0.30  0.55  0.25    0.2  0.0      0.0 /  
  3     'GENROU' 1    8    0.03   0.4  0.05  6.175  0  1.8  1.7   0.30  0.55  0.25    0.2  0.0      0.0 /   
  4     'GENROU' 1    8    0.03   0.4  0.05  6.175  0  1.8  1.7   0.30  0.55  0.25    0.2  0.0      0.0 /  
  1     'SEXS'  1    0.10000     10.000       100.00      0.10000   0.0000  5.0000  /
  2     'SEXS'  1    0.10000     10.000       100.00      0.10000   0.0000  5.0000  /
  3     'SEXS'  1    0.10000     10.000       100.00      0.10000   0.0000  5.0000  /
  4     'SEXS'  1    0.10000     10.000       100.00      0.10000   0.0000  5.0000  /
  1    'HYGOV'  1    0.04  0.5       1.0      0.05   0.20   0.10  1.0   0.00   1  1.0  0.0   0.0 / 
  2    'HYGOV'  1    0.04  0.5       1.0      0.05   0.20   0.10  1.0   0.00   1  1.0  0.0   0.0 / 
  3    'HYGOV'  1    0.04  0.5       1.0      0.05   0.20   0.10  1.0   0.00   1  1.0  0.0   0.0 / 
  4    'HYGOV'  1    0.04  0.5       1.0      0.05   0.20   0.10  1.0   0.00   1  1.0  0.0   0.0 / 

The error when doing: ss.run(raw_path, addfile = dyr_path, routine = 'eig') is the following:

line 221, in set_axis
    f"Length mismatch: Expected axis has {old_len} elements, new "
ValueError: Length mismatch: Expected axis has 14 elements, new values have 13 elements

Removing the HYGOV entries in the dyr fixes the problem.

cuihantao commented 2 years ago

Thanks for reporting the issue. Might be an issue with the yaml file defining the columns.  I will look into it.

Regards, Hantao Cui On Apr 15, 2022, 4:19 PM -0500, Rodrigo Henríquez-Auba @.***>, wrote:

Assigned #310 to @cuihantao. — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you were assigned.Message ID: @.***>

cuihantao commented 2 years ago

Fixed in https://github.com/cuihantao/andes/commit/174f8f4d11c0331d095811a31be9e45a3280ec5b. Currently available in my develop branch and will be tagged soon.