pydae 0.54.0__py3-none-any.whl → 0.54.1__py3-none-any.whl

pydae/__init__.py

@@ -7,5 +7,5 @@ Differential Algebraic Equations in Python
 
 """
 
-__version__ = "0.54.0"
+__version__ = "0.54.1"
  

pydae/urisi/urisi_builder.py

@@ -70,7 +70,7 @@ class urisi:
                     'u_ini_dict':{},'u_run_dict':{},'params_dict':{},
                     'h_dict':{},'xy_0_dict':{}}
 
-        self.uz_jacs = True     
+        self.uz_jacs = False     
         self.verbose = False 
 
         self.omega_coi_numerator = 0.0

pydae/urisi/utils/model2svg.py

@@ -19,6 +19,7 @@ class model2svg(svg):
         self.set_buses_title()
         self.set_lines_title()
         self.set_transformers_title()
+        self.set_sources_title()
 
         self.set_lines_currents_colors()
         self.set_buses_voltages_colors()
@@ -53,27 +54,42 @@ class model2svg(svg):
                 U_base = bus['U_kV']*1000
                 V_base = U_base/np.sqrt(3)
                 string = 'Voltages and powers:\n'
+                p_total = 0.0
+                q_total = 0.0
                 for phn in ['an','bn','cn','ng']:
                     V = get_v(self.model,bus_name,f'V_{phn}_m')
                     V_pu = V/V_base
 
                     p = 0.0
                     q = 0.0
+
                     if bus_name in load_buses:
                         ph = phn[0]
                         if  not ph == 'n':
                             p = self.model.get_value(f'p_load_{bus_name}_{ph}')
                             q = self.model.get_value(f'q_load_{bus_name}_{ph}')
 
+                            p_total += p
+                            q_total += q
+
                     if not phn == 'ng':  
                         if q < 0.0:    
-                            string += f"V{phn} = {V:7.1f} V ({V_pu:4.2f} pu), S{phn} = {p/1e3:5.1f} - {np.abs(q)/1e3:5.1f}j kVA \n"
+                            string += f"V{phn} = {V:7.1f} V ({V_pu:4.2f} pu), S{phn[0]} = {p/1e3:5.1f} - {np.abs(q)/1e3:5.1f}j kVA \n"
                         else:
-                            string += f"V{phn} = {V:7.1f} V ({V_pu:4.2f} pu), S{phn} = {p/1e3:5.1f} + {np.abs(q)/1e3:5.1f}j kVA \n"
+                            string += f"V{phn} = {V:7.1f} V ({V_pu:4.2f} pu), S{phn[0]} = {p/1e3:5.1f} + {np.abs(q)/1e3:5.1f}j kVA \n"
 
                     else:
-                        string += f"V{phn} = {V:7.1f} V ({V_pu:4.2f} pu)\n"
-    
+                        string += f"V{phn} = {V:7.1f} V ({V_pu:4.2f} pu), St = {p_total/1e3:5.1f} + {np.abs(q_total)/1e3:5.1f}j kVA \n"
+
+                v_m_array = get_v(self.model,bus_name,'V_abcn_m')
+                v_m_min = np.min(v_m_array)
+                v_m_max = np.max(v_m_array)
+                max_dev = v_m_max - v_m_min
+                v_avg = np.sum(v_m_array)/3
+                unbalance = max_dev/v_avg
+                string += f"voltage unbalance = {unbalance*100:5.2f} %\n"
+
+
             if bus['N_nodes'] == 3:
                 bus_name = bus['name']
                 U_base = bus['U_kV']*1000
@@ -211,4 +227,28 @@ class model2svg(svg):
             string += f"Ic = {np.abs(I_1_c):6.1f} A\tIc = {np.abs(I_2_c):6.1f} A \n"
             string += f"\t\t\t  In = {np.abs(I_2_n):6.1f} A \n"
 
-            self.set_title(f'trafo_{bus_j}_{bus_k}_g',string) 
+            self.set_title(f'trafo_{bus_j}_{bus_k}_g',string) 
+
+
+    def set_sources_title(self):
+
+        for src in self.grid_data['sources']:
+            bus = src['bus']
+
+            I_a = self.model.get_value(f'i_vsc_{bus}_a_m')
+            I_b = self.model.get_value(f'i_vsc_{bus}_a_m')
+            I_c = self.model.get_value(f'i_vsc_{bus}_a_m')
+
+            p = self.model.get_value(f'p_{bus}')
+            q = self.model.get_value(f'q_{bus}')
+
+            string = 'Currents:\n'
+            string += f"Ia = {I_a:6.1f} A \n"
+            string += f"Ib = {I_b:6.1f} A \n"
+            string += f"Ic = {I_c:6.1f} A \n"
+            string += 'Powers:\n'
+            string += f"p = {p/1e3:6.1f} kW\n"
+            string += f"q = {q/1e3:6.1f} kvar \n"
+
+
+            self.set_title(f'src_{bus}',string) 

{pydae-0.54.0.dist-info → pydae-0.54.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pydae
-Version: 0.54.0
+Version: 0.54.1
 Summary: =====
 Keywords: engineering
 Author-email: Juan Manuel Mauricio <jmmauricio6@gmail.com>

{pydae-0.54.0.dist-info → pydae-0.54.1.dist-info}/RECORD

@@ -1,5 +1,5 @@
 pydae/.gitignore,sha256=c1jUQV9WXX3PWrlOzHiKvUWFomLJUrfYogbU7nYwAwc,288
-pydae/__init__.py,sha256=9UHMkOqnidL3wwAP8vNX-luL1D3WZn-9zVDzGGNOtzE,106
+pydae/__init__.py,sha256=5Xl8hDbIFcwcha4-6VIUFCHNG3cejkRWVU9qO4q8jWs,106
 pydae/b_pu.py,sha256=FdRRVieJLQDYPwYNWhbXJ0SXvf4jtPdbvIZTKWi16I0,6775
 pydae/bpu.py,sha256=FdRRVieJLQDYPwYNWhbXJ0SXvf4jtPdbvIZTKWi16I0,6775
 pydae/build.py,sha256=ECfc_doyK9PhrjEAlWJGLfirQmpT24XuYNxUgtrXI74,46151
@@ -453,7 +453,7 @@ pydae/templates/solver_template_v2.py,sha256=ciKWIwlHECY39kjBQJBnow4pA8bKTif8nfB
 pydae/templates/unit_test_template.py,sha256=tSRNAwtmcHc4QujOGzZFXDge1vTRiyIcbvemerOFROQ,1117
 pydae/templates/uz_jacs_v2.py,sha256=JzYMbn6N015656BAWrkkOnrrys9akQ1y1GRH60ZJOAY,4340
 pydae/urisi/__init__.py,sha256=r7mgVYduCfuZ_pqaNXp_Y6Y-L3wcOKT4dzBLPsr8oak,105
-pydae/urisi/urisi_builder.py,sha256=KJzDEqSUVWvWa-VdLlLRkuDXYUD2KAIaeilC66QY8_A,13870
+pydae/urisi/urisi_builder.py,sha256=VagNhnOR4XzXLPvGpwoPbGJJGRkC1Gw_1AdZYUteaLE,13871
 pydae/urisi/xy_0.json,sha256=mA-fx_pCJkcXtZzhO1zRS5GONGvgFs-Gj1gqszaGaR0,842
 pydae/urisi/ess/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 pydae/urisi/ess/bess.ipynb,sha256=Glx6-BcexpCl-j_pu_TJNYD8uN5ouX9F7kina33hOYI,54561
@@ -539,7 +539,7 @@ pydae/urisi/transformers/transformers.py,sha256=qquEa6j0HXWXqwROcZEifUhm6fZpkSx3
 pydae/urisi/transformers/xy_0.json,sha256=2Q_yStbln3nKPxHTbwAloe--eGlI0IGkqY1L7QTJWYI,881
 pydae/urisi/utils/__init__.py,sha256=byp7NSkTjT3waG_5zrUeXPPm8xlebXNGhO7a9orE66k,201
 pydae/urisi/utils/measurements.py,sha256=ubUEMkvGCe_2dEWHu7ID2z1J5ODmG2brUS8fbKpgLyM,4541
-pydae/urisi/utils/model2svg.py,sha256=G3TyJeMG29eKvkbP0r2mFSlR6GB795wXSxdkRldC6xM,9021
+pydae/urisi/utils/model2svg.py,sha256=DTTGpFMkBxT6XVKPcliArXBWTjsQl3Q5WVep4KMggM8,10404
 pydae/urisi/utils/reports.py,sha256=Gge0iIGpNTNGaK-Ex0esYlfLyH3q6VIphabyvmaEc-Y,4944
 pydae/urisi/vsc_ctrls/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 pydae/urisi/vsc_ctrls/ctrl_3ph_4w_droop.py,sha256=iW-34_Dh0SXNvq5Z-pMsNf8JdT8SQFipb5zE1sn7psQ,3638
@@ -622,8 +622,8 @@ pydae/urisi/vsgs/xy_0.json,sha256=Q6qWPIqEP-w5Qu2twd3S4PX4zMggpNpG2W62e7ry70g,85
 pydae/utils/__init__.py,sha256=kwp76qeU2mnwrCexDSF15id9ih12Q0VRL3l-_c6omcQ,184
 pydae/utils/tools.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 pydae/utils/utils.py,sha256=OrKP6NAUGdCOa6gO_MYm8FdvVuIFFSlt4jYkLsrbUtg,1069
-pydae-0.54.0.dist-info/COPYING,sha256=QWnWoslbTjQmtlFSRoGF3wxJEJuTJdlEMZ7lHtbGAeQ,1139
-pydae-0.54.0.dist-info/LICENSE,sha256=8hQe1oM4ySdliF3R-NEvR6HqrcGDKvsLFJC3gA1sNjY,1108
-pydae-0.54.0.dist-info/WHEEL,sha256=EZbGkh7Ie4PoZfRQ8I0ZuP9VklN_TvcZ6DSE5Uar4z4,81
-pydae-0.54.0.dist-info/METADATA,sha256=cJ73vHB8nZ80B2gf7iw_PibBC3Qfyt0KdSMggVk9Cwg,1120
-pydae-0.54.0.dist-info/RECORD,,
+pydae-0.54.1.dist-info/COPYING,sha256=QWnWoslbTjQmtlFSRoGF3wxJEJuTJdlEMZ7lHtbGAeQ,1139
+pydae-0.54.1.dist-info/LICENSE,sha256=8hQe1oM4ySdliF3R-NEvR6HqrcGDKvsLFJC3gA1sNjY,1108
+pydae-0.54.1.dist-info/WHEEL,sha256=EZbGkh7Ie4PoZfRQ8I0ZuP9VklN_TvcZ6DSE5Uar4z4,81
+pydae-0.54.1.dist-info/METADATA,sha256=CvIMhokj1YM76g7mbgZJQ4eDGCdks6YlQbgesZDUxj8,1120
+pydae-0.54.1.dist-info/RECORD,,