Allow steady state power flow to be solved when DPLine is off.

CURENT · Mar 10, 2024 · d409483 · d409483
1 parent 5cb1a91
commit d409483
Showing 1 changed file with 6 additions and 4 deletions.
diff --git a/andes/models/line/line.py b/andes/models/line/line.py
@@ -253,14 +253,16 @@ def __init__(self, system=None, config=None):
         self.idd = State(info='real current',
                          tex_name='idd',
                          v_str='1e-8',
-                         e_str='u * ( -(x+1e-8)*iqq - (r+1e-8)*idd - v2*sin(a2) + v1*sin(a1) )',
-                         t_const=self.Leq)
+                         e_str='u * ( -(x+1e-8)*iqq - (r+1e-8)*idd - v2*sin(a2) + v1*sin(a1) ) + (1-u) * idd',
+                         t_const=self.Leq,
+                         )
 
         self.iqq = State(info='real current',
                          tex_name='iqq',
                          v_str='1e-8',
-                         e_str='u * (1*(x+1e-8)*idd - (r+1e-8)*iqq - v2*cos(a2) + v1*cos(a1) )',
-                         t_const=self.Leq)
+                         e_str='u * (1*(x+1e-8)*idd - (r+1e-8)*iqq - v2*cos(a2) + v1*cos(a1) ) + (1-u) * iqq',
+                         t_const=self.Leq,
+                         )
 
         self.a1.e_str = 'u * (idd*v1*sin(a1) + iqq*v1*cos(a1))'