Georg_Brantegger/Python-DT_Slot_3
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

probably working version (with unrealistic controller design)

Browse Source
This commit is contained in:
Brantegger Georg
2022-08-08 07:46:35 +02:00
parent ba696444bb
commit 5a790d5ca5
1 changed files with 8 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
Untertweng_mit_Pegelregler.ipynb
View File

@@ -18,7 +18,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 9, "execution_count": 2,
"metadata": {}, "metadata": {},
"outputs": [], "outputs": [],
"source": [ "source": [
@@ -40,7 +40,7 @@
" # for PI controller\n", " # for PI controller\n",
"Con_targetLevel = 8. # [m]\n", "Con_targetLevel = 8. # [m]\n",
"Con_K_p = 0.1 # [-] proportional constant of PI controller\n", "Con_K_p = 0.1 # [-] proportional constant of PI controller\n",
"Con_T_i = 10. # [s] timespan in which a steady state error is corrected by the intergal term\n", "Con_T_i = 1000. # [s] timespan in which a steady state error is corrected by the intergal term\n",
"Con_deadbandRange = 0.05 # [m] Deadband range around targetLevel for which the controller does NOT intervene\n", "Con_deadbandRange = 0.05 # [m] Deadband range around targetLevel for which the controller does NOT intervene\n",
"\n", "\n",
"\n", "\n",
@@ -99,7 +99,7 @@
"turbine.set_steady_state(flux_init,pipe.get_current_pressure_distribution()[-1])\n", "turbine.set_steady_state(flux_init,pipe.get_current_pressure_distribution()[-1])\n",
"\n", "\n",
"# influx setting turbine\n", "# influx setting turbine\n",
"turbine_in = Francis_Turbine(Tur_Q_nenn,Tur_p_nenn,Tur_closingTime/2,Pip_dt,pUnit_conv)\n", "turbine_in = Francis_Turbine(Tur_Q_nenn,Tur_p_nenn,Tur_closingTime,Pip_dt,pUnit_conv)\n",
"turbine_in.set_steady_state(flux_init,Tur_p_nenn)\n", "turbine_in.set_steady_state(flux_init,Tur_p_nenn)\n",
"\n", "\n",
"# level controll\n", "# level controll\n",
@@ -149,10 +149,11 @@
"LA_ist_vec = np.full_like(t_vec,turbine.get_current_LA())\n", "LA_ist_vec = np.full_like(t_vec,turbine.get_current_LA())\n",
"\n", "\n",
"LA_soll_vec2 = np.full_like(t_vec,turbine_in.get_current_LA())\n", "LA_soll_vec2 = np.full_like(t_vec,turbine_in.get_current_LA())\n",
"LA_soll_vec2[500:1000] = 0.\n", "LA_soll_vec2[500:] = 0\n",
"LA_soll_vec2[1000:1500] = 1. \n", "# LA_soll_vec2[500:1000] = 0.\n",
"LA_soll_vec2[1500:2000] = 0.\n", "# LA_soll_vec2[1000:1500] = 1. \n",
"LA_soll_vec2[2000:2500] = 0.5 \n" "# LA_soll_vec2[1500:2000] = 0.\n",
"# LA_soll_vec2[2000:2500] = 0.5 \n"
] ]
}, },
{ {