1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
#ifndef ELECTROMECHANICAL_H
#define ELECTROMECHANICAL_H
#include "ElectricCalculation.h"
#include <wx/progdlg.h>
class ControlElementSolver;
/**
* @brief Synchronous machine data for different models
*
* This struct have the parameters and calculated data that depends on each machine model.
*/
struct SyncMachineModelData {
/** Direct-axis reactance (transient for models 1, 2 and 3; subtransient for models 4 and 5) */
double xd;
/** Quadrature-axis reactance (transient for models 1, 2 and 3; subtransient for models 4 and 5) */
double xq;
/** Direct-axis internal voltage (transient for models 1, 2 and 3; subtransient for models 4 and 5) */
double ed;
/** Quadrature-axis internal voltage (transient for models 1, 2 and 3; subtransient for models 4 and 5) */
double eq;
};
/**
* @class Electromechanical
* @author Thales Lima Oliveira <thales@ufu.br>
* @date 23/09/2017
* @file Electromechanical.h
* @brief Calculates the electromechanical transient based on disturbances (e.g. system fault).
*/
class Electromechanical : public ElectricCalculation
{
public:
Electromechanical(wxWindow* parent, std::vector<Element*> elementList, SimulationData data);
~Electromechanical();
bool RunStabilityCalculation();
wxString GetErrorMessage() const { return m_errorMsg; }
std::vector<double> GetTimeVector() const { return m_timeVector; }
std::vector<double> m_wErrorVector;
std::vector<double> m_deltaErrorVector;
std::vector<double> m_transEdErrorVector;
std::vector<double> m_transEqErrorVector;
std::vector<double> m_numItVector;
protected:
void SetEventTimeList();
bool HasEvent(double currentTime);
void SetEvent(double currentTime);
inline bool EventTrigger(double eventTime, double currentTime);
// double GetPowerValue(double value, ElectricalUnit unit);
void InsertSyncMachinesOnYBus();
std::complex<double> GetSyncMachineAdmittance(SyncGenerator* generator);
bool InitializeDynamicElements();
bool CalculateMachinesCurrents();
void CalculateIntegrationConstants(SyncGenerator* syncGenerator, double id, double iq, double k = 1.0);
bool SolveSynchronousMachines();
void SetSyncMachinesModel();
SyncMachineModelData GetSyncMachineModelData(SyncGenerator* syncMachine);
double CalculateSyncMachineIntVariables(SyncGenerator* syncGenerator,
double id,
double iq,
double sd,
double sq,
double pe,
double k = 1.0);
void CalculateSyncMachineNonIntVariables(SyncGenerator* syncGenerator,
double& id,
double& iq,
double& sd,
double& sq,
double& pe,
double k = 1.0);
void CalculateReferenceSpeed();
bool CalculateSyncMachineSaturation(SyncGenerator* syncMachine, double k = 1.0);
void SaveData();
wxWindow* m_parent = NULL;
wxString m_errorMsg = _("Unknown error");
double m_systemFreq = 60.0;
double m_refSpeed = 2.0 * M_PI * 60.0;
bool m_useCOI = false;
std::vector<std::vector<std::complex<double> > > m_yBus;
std::vector<std::vector<std::complex<double> > > m_yBusU;
std::vector<std::vector<std::complex<double> > > m_yBusL;
std::vector<std::complex<double> > m_vBus;
std::vector<std::complex<double> > m_iBus;
double m_powerSystemBase = 100e6;
double m_simTime = 10.0;
double m_plotTime = 1e-2;
double m_timeStep = 1e-2;
double m_ctrlTimeStepMultiplier = 0.1;
double m_tolerance = 1e-8;
int m_maxIterations = 100;
double m_saturationTolerance = 1e-4;
std::vector<double> m_eventTimeList;
std::vector<bool> m_eventOccurrenceList;
std::vector<double> m_timeVector;
// tests
double m_wError = 0.0;
double m_numIt = 0;
};
#endif // ELECTROMECHANICAL_H
|
