Power System Platform  2026w23a-beta
Loading...
Searching...
No Matches
EMTElementForm.cpp
1#include "EMTElementForm.h"
2#include "../elements/powerElement/EMTElement.h"
3#include "../editors/ChartView.h"
4#include "../editors/Workspace.h"
5#include "../extLibs/fftw/fftw3.h"
6#include "ATPFileEditorForm.h"
7
8#include <wx/wfstream.h>
9#include <wx/datstrm.h>
10#include <wx/dir.h>
11#include <wx/process.h>
12
13EMTElementForm::EMTElementForm(wxWindow* parent, EMTElement* emtElement, Workspace* workspace)
14 : EMTElementFormBase(parent), m_parent(parent), m_emtElement(emtElement), m_workspace(workspace)
15{
16 m_staticTextNode->SetLabel(m_staticTextNode->GetLabel() + wxT(" (3\u0278)"));
17
18 auto data = m_emtElement->GetEMTElementData();
19 m_textCtrlName->SetValue(data.name);
20 m_filePickerATPFile->SetPath(data.atpFile.GetFullPath());
21 LoadNodes(data.atpFile.GetFullPath());
22 m_choiceNodes->SetSelection(m_choiceNodes->FindString(data.atpNodeName));
23 m_textCtrlStepSize->SetValue(wxString::Format(wxT("%e"), data.stepSize));
24 m_spinCtrlCyclesToSS->SetValue(data.cyclesToSS);
25 m_textCtrlRecFreq->SetValue(wxString::Format(wxT("%d"), data.recordFrequency));
26 m_checkBoxUseFilter->SetValue(data.useMedianFilter);
27 m_spinCtrlMaxNumHarmonics->SetValue(data.numMaxHarmonics);
28 m_textCtrlHarmThreshold->SetValue(wxString::Format(wxT("%.2f"), data.harmonicsThreshold));
29
30 m_properties = workspace->GetProperties();
31
32 SetSize(DoGetBestSize());
33}
34
35EMTElementForm::~EMTElementForm()
36{
37}
38
39bool EMTElementForm::ValidateData()
40{
41 wxFileName fileName(m_filePickerATPFile->GetPath());
42 if (!fileName.IsOk()) {
43 wxMessageDialog msgDialog(m_parent, _("The ATP file doesn't exist. Insert the correct path."), _("Error"), wxOK | wxCENTRE | wxICON_ERROR);
44 msgDialog.ShowModal();
45 return false;
46 }
47 double stepSize = 0.0;
48 if (!m_textCtrlStepSize->GetValue().ToDouble(&stepSize)) {
49 wxMessageDialog msgDialog(m_parent, _("The step size must be a number."), _("Error"), wxOK | wxCENTRE | wxICON_ERROR);
50 msgDialog.ShowModal();
51 return false;
52 }
53 long recFreq = 0;
54 if (!m_textCtrlRecFreq->GetValue().ToLong(&recFreq)) {
55 wxMessageDialog msgDialog(m_parent, _("The steps to record data must be an integer number."), _("Error"), wxOK | wxCENTRE | wxICON_ERROR);
56 msgDialog.ShowModal();
57 return false;
58 }
59 if (m_choiceNodes->GetSelection() == wxNOT_FOUND) {
60 wxMessageDialog msgDialog(m_parent, _("No connection node selected."), _("Error"), wxOK | wxCENTRE | wxICON_ERROR);
61 msgDialog.ShowModal();
62 return false;
63 }
64 double harmThreshold = 0.0;
65 if (!m_textCtrlHarmThreshold->GetValue().ToDouble(&harmThreshold)) {
66 wxMessageDialog msgDialog(m_parent, _("The threshold to detect harmonics must be a number."), _("Error"), wxOK | wxCENTRE | wxICON_ERROR);
67 msgDialog.ShowModal();
68 return false;
69 }
70
71 auto data = m_emtElement->GetEMTElementData();
72 data.name = m_textCtrlName->GetValue();
73 data.atpFile = m_filePickerATPFile->GetPath();
74 data.atpNodeName = m_choiceNodes->GetString(m_choiceNodes->GetSelection());
75 data.stepSize = stepSize;
76 data.cyclesToSS = m_spinCtrlCyclesToSS->GetValue();
77 data.recordFrequency = recFreq;
78 data.useMedianFilter = m_checkBoxUseFilter->GetValue();
79 data.numMaxHarmonics = m_spinCtrlMaxNumHarmonics->GetValue();
80 data.harmonicsThreshold = harmThreshold;
81
82 m_emtElement->SetEMTElementData(data);
83
84 return true;
85}
86
87void EMTElementForm::OnCancelClick(wxCommandEvent& event)
88{
89 EndModal(wxID_CANCEL);
90}
91void EMTElementForm::OnOKClick(wxCommandEvent& event)
92{
93 if (ValidateData()) EndModal(wxID_OK);
94}
95
96void EMTElementForm::LoadNodes(wxString atpFilePath)
97{
98 wxArrayString atpFile;
99 wxTextFile file(atpFilePath);
100 if (!file.Exists()) return;
101
102 if (file.Open(atpFilePath)) {
103 atpFile.Add(file.GetFirstLine());
104 while (!file.Eof())
105 {
106 atpFile.Add(file.GetNextLine());
107 }
108 file.Close();
109 }
110 wxArrayString nodeList = m_emtElement->GetATPNodes(atpFile);
111 m_choiceNodes->Clear();
112 m_choiceNodes->Append(nodeList);
113}
114
115void EMTElementForm::OnEditATPFileClick(wxCommandEvent& event)
116{
117 wxFileName fileName(m_filePickerATPFile->GetPath());
118 if (!fileName.IsOk()) {
119 wxMessageDialog msgDialog(m_parent, _("The ATP file doesn't exist. Insert the correct path."), _("Error"), wxOK | wxCENTRE | wxICON_ERROR);
120 msgDialog.ShowModal();
121 return;
122 }
123 ATPFileEditorForm atpFileEditorForm(m_parent, m_filePickerATPFile->GetPath());
124 atpFileEditorForm.SetSize(800, 600);
125 atpFileEditorForm.ShowModal();
126}
127
128void EMTElementForm::OnTestClick(wxCommandEvent& event)
129{
130 if (ValidateData()) {
131 wxString errorMsg;
132 m_emtElement->UpdateData(m_properties);
133 if (!m_emtElement->CalculateCurrent(errorMsg, Mode::Both, true, true)) {
134 wxMessageDialog msgDialog(m_parent, errorMsg, _("Error"), wxOK | wxCENTRE | wxICON_ERROR);
135 msgDialog.ShowModal();
136 return;
137 }
138 std::vector<double> time, in, timeSamp, inSamp, outR, outI, outM, outP, freq;
139 std::vector <double> timeThreePhase, currentA, currentB, currentC;
140 std::vector <double> magIa, magIb, magIc, angIa, angIb, angIc;
141 //auto atpData = m_emtElement->GetEMTElementData().atpData;
142 auto inFFTData = m_emtElement->GetEMTElementData().inFFTData;
143 auto outFFTData = m_emtElement->GetEMTElementData().outFFTData;
144 //for (auto data : atpData) {
145 // time.emplace_back(data.first);
146 // in.emplace_back(data.second);
147 //}
148 for (auto& data : m_emtElement->GetEMTElementData().atpSampleData)
149 {
150 time.emplace_back(data.t);
151 currentA.emplace_back(data.current[0]);
152 currentB.emplace_back(data.current[1]);
153 currentC.emplace_back(data.current[2]);
154 magIa.emplace_back(data.mag[0]);
155 magIb.emplace_back(data.mag[1]);
156 magIc.emplace_back(data.mag[2]);
157 angIa.emplace_back(data.angle[0] * 180.0 / M_PI);
158 angIb.emplace_back(data.angle[1] * 180.0 / M_PI);
159 angIc.emplace_back(data.angle[2] * 180.0 / M_PI);
160
161 }
162 for (auto data : inFFTData) {
163 timeSamp.emplace_back(data.first);
164 inSamp.emplace_back(data.second);
165 }
166 for (auto data : outFFTData) {
167 outR.emplace_back(data.second.real());
168 outI.emplace_back(data.second.imag());
169 outM.emplace_back(abs(data.second));
170 outP.emplace_back(arg(data.second) * 180.0 / M_PI);
171 freq.emplace_back(data.first);
172 }
173
174
175 std::vector<double> currMag, currPha, harmOrder;
176 auto harmData = m_emtElement->GetEMTElementData().currHarmonics;
177 //auto harmOrderData = m_emtElement->GetEMTElementData().currHarmonicsOrder;
178 //for (size_t i = 0; i < harmData.size(); i++) {
179 // currMag.emplace_back(abs(harmData[i]));
180 // currPha.emplace_back(arg(harmData[i]) * 180.0 / M_PI);
181 // harmOrder.emplace_back(static_cast<double>(harmOrderData[i]));
182 //}
183 for (auto const& data : harmData) {
184 currMag.emplace_back(abs(data.second));
185 currPha.emplace_back(arg(data.second) * 180.0 / M_PI);
186 harmOrder.emplace_back(static_cast<double>(data.first));
187 }
188
189 std::vector<ElementPlotData> plotDataList;
190 //ElementPlotData plotATPData(_("ATP Output"), ElementPlotData::CurveType::CT_TEST);
191 ElementPlotData plot3PATPData(_("ATP Three-phase output"), ElementPlotData::CurveType::CT_TEST);
192 ElementPlotData plotFFTData(_("Fourier Analysis"), ElementPlotData::CurveType::CT_TEST);
193 ElementPlotData plotCurrData(_("Current Phasor"), ElementPlotData::CurveType::CT_TEST);
194
195 //plotATPData.AddData(in, _("Data"));
196
197 plot3PATPData.AddData(currentA, _("Phase A Current"));
198 plot3PATPData.AddData(currentB, _("Phase B Current"));
199 plot3PATPData.AddData(currentC, _("Phase C Current"));
200 plot3PATPData.AddData(magIa, _("Magnitude A Current"));
201 plot3PATPData.AddData(magIb, _("Magnitude B Current"));
202 plot3PATPData.AddData(magIc, _("Magnitude C Current"));
203 plot3PATPData.AddData(angIa, _("Angle A Current"));
204 plot3PATPData.AddData(angIb, _("Angle B Current"));
205 plot3PATPData.AddData(angIc, _("Angle C Current"));
206 plot3PATPData.AddData(timeThreePhase, _("Time 3f"));
207
208 plotFFTData.AddData(inSamp, _("Input"));
209 plotFFTData.AddData(timeSamp, _("Time Sampled"));
210 plotFFTData.AddData(outR, _("Real"));
211 plotFFTData.AddData(outI, _("Imaginary"));
212 plotFFTData.AddData(outM, _("Magnitude"));
213 plotFFTData.AddData(outP, _("Phase"));
214 plotFFTData.AddData(freq, _("Frequency"));
215
216 plotCurrData.AddData(currMag, _("Magnitude"));
217 plotCurrData.AddData(currPha, _("Phase"));
218 plotCurrData.AddData(harmOrder, _("Harmonic Order"));
219
220 //plotDataList.push_back(plotATPData);
221 plotDataList.push_back(plot3PATPData);
222 plotDataList.push_back(plotFFTData);
223 plotDataList.push_back(plotCurrData);
224 //ChartView* cView = new ChartView(m_parent, plotDataList, time, m_properties->GetGeneralPropertiesData().plotLib);
225 //cView->Show();
226 auto& cwMngr = m_workspace->GetChartWindowManager();
227 cwMngr.Open(plotDataList, time);
228 }
229}
230void EMTElementForm::OnATPFileChange(wxFileDirPickerEvent& event)
231{
232 LoadNodes(m_filePickerATPFile->GetPath());
233}
EMTElement
Element to connect ATP-EMTP.
Definition EMTElement.h:83
Workspace
This class manages the graphical and power elements. It is responsible for handling the user's intera...
Definition Workspace.h:105