pssdlg.cpp

电力系统暂态稳定分析程序

// PSSDlg.cpp : implementation file
//

#include "stdafx.h"
#include "PSS.h"
#include "PSSDlg.h"
#include "DlgGen.h"
#include "DlgUpTran.h"
#include "Dlgxianlu.h"
#include "DlgDownTran.h"
#include "DlgLoad.h"
#include "math.h"
#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About

class CAboutDlg : public CDialog
{
public:
	CAboutDlg();

// Dialog Data
	//{{AFX_DATA(CAboutDlg)
	enum { IDD = IDD_ABOUTBOX };
	//}}AFX_DATA

	// ClassWizard generated virtual function overrides
	//{{AFX_VIRTUAL(CAboutDlg)
	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
	//}}AFX_VIRTUAL

// Implementation
protected:
	//{{AFX_MSG(CAboutDlg)
	//}}AFX_MSG
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
	//{{AFX_DATA_INIT(CAboutDlg)
	//}}AFX_DATA_INIT
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CAboutDlg)
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
	//{{AFX_MSG_MAP(CAboutDlg)
		// No message handlers
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CPSSDlg dialog

CPSSDlg::CPSSDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CPSSDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CPSSDlg)
	m_sb = 250.0f;
	m_vb = 115.0f;
	m_dt = 0.0f;
	m_t = 0.0f;
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CPSSDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CPSSDlg)
	DDX_Text(pDX, IDC_EDIT5, m_sb);
	DDX_Text(pDX, IDC_EDIT6, m_vb);
	DDX_Text(pDX, IDC_EDIT7, m_dt);
	DDX_Text(pDX, IDC_EDIT8, m_t);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CPSSDlg, CDialog)
	//{{AFX_MSG_MAP(CPSSDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_BUTTON_Gen, OnBUTTONGen)
	ON_BN_CLICKED(IDC_BUTTON_UpTran, OnBUTTONUpTran)
	ON_BN_CLICKED(IDC_BUTTON_Line, OnBUTTONLine)
	ON_BN_CLICKED(IDC_BUTTON_DownTran, OnBUTTONDownTran)
	ON_BN_CLICKED(IDC_BUTTON_Load, OnBUTTONLoad)
	ON_BN_CLICKED(IDC_Caculate, OnCaculate)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CPSSDlg message handlers

BOOL CPSSDlg::OnInitDialog()
{
	CDialog::OnInitDialog();

	// Add "About..." menu item to system menu.

	// IDM_ABOUTBOX must be in the system command range.
	ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
	ASSERT(IDM_ABOUTBOX < 0xF000);

	CMenu* pSysMenu = GetSystemMenu(FALSE);
	if (pSysMenu != NULL)
	{
		CString strAboutMenu;
		strAboutMenu.LoadString(IDS_ABOUTBOX);
		if (!strAboutMenu.IsEmpty())
		{
			pSysMenu->AppendMenu(MF_SEPARATOR);
			pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
		}
	}

	// Set the icon for this dialog.  The framework does this automatically
	//  when the application's main window is not a dialog
	SetIcon(m_hIcon, TRUE);			// Set big icon
	SetIcon(m_hIcon, FALSE);		// Set small icon
	
	// TODO: Add extra initialization here
	
	return TRUE;  // return TRUE  unless you set the focus to a control
}

void CPSSDlg::OnSysCommand(UINT nID, LPARAM lParam)
{
	if ((nID & 0xFFF0) == IDM_ABOUTBOX)
	{
		CAboutDlg dlgAbout;
		dlgAbout.DoModal();
	}
	else
	{
		CDialog::OnSysCommand(nID, lParam);
	}
}

// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.

void CPSSDlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // device context for painting

		SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);

		// Center icon in client rectangle
		int cxIcon = GetSystemMetrics(SM_CXICON);
		int cyIcon = GetSystemMetrics(SM_CYICON);
		CRect rect;
		GetClientRect(&rect);
		int x = (rect.Width() - cxIcon + 1) / 2;
		int y = (rect.Height() - cyIcon + 1) / 2;

		// Draw the icon
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

// The system calls this to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CPSSDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

void CPSSDlg::OnBUTTONGen() 
{
	CDlgGen dlg;
	dlg.DoModal();
	
}

void CPSSDlg::OnBUTTONUpTran() 
{
	CDlgUpTran dlg;
	dlg.DoModal();
	
}

void CPSSDlg::OnBUTTONLine() 
{
	CDlgXianlu dlg;
	dlg.DoModal();
	
}

void CPSSDlg::OnBUTTONDownTran() 
{
	CDlgDownTran dlg;
	dlg.DoModal();
	
}

void CPSSDlg::OnBUTTONLoad() 
{
	CDlgLoad dlg;
	dlg.DoModal();
	
}

extern float t1_k,t1_s,t1_v,t1_vn_u,t1_vn_d;
extern float t2_k,t2_s,t2_v,t2_vn_u,t2_vn_d;
extern float s_s,s_v,s_xd,s_xq,s_xdd,s_x2,s_tj;
extern float l_l,l_x,l_x0,l_v;
extern float load_v,load_p,load_cos;
void CPSSDlg::OnCaculate() 
{
 UpdateData(TRUE);
float vb_2=m_vb*t2_k;
float vb_1=vb_2*t1_k;
float pu_xd=s_xd*s_v*s_v*m_sb/(s_s*vb_1*vb_1);
float pu_xq=s_xq*s_v*s_v*m_sb/(s_s*vb_1*vb_1);
float pu_xdd=s_xdd*s_v*s_v*m_sb/(s_s*vb_1*vb_1);
float pu_x_t1=t1_v*t1_vn_u*t1_vn_u*m_sb/(100*t1_s*vb_2*vb_2);
float pu_x_t2=t2_v*t2_vn_u*t2_vn_u*m_sb/(100*t2_s*vb_2*vb_2);
float pu_x_l=l_l*l_x*m_sb/(vb_2*vb_2);
float pu_x_tl=pu_x_t1+pu_x_t2+0.5*pu_x_l;
float xd_sum=pu_xd+pu_x_tl;
float xq_sum=pu_xq+pu_x_tl;
float xdd_sum=pu_xdd+pu_x_tl;
float pu_v0=load_v/m_vb;
float a0=acos(load_cos)*180/3.1415;
float pu_p=load_p/m_sb;
float pu_q=pu_p*tan(a0*3.1415/180);
float e_Q0=sqrt((pu_v0+pu_q*xq_sum/pu_v0)*(pu_v0+pu_q*xq_sum/pu_v0)+(pu_p*xq_sum/pu_v0)*(pu_p*xq_sum/pu_v0));
float as0=atan((pu_p*xq_sum/pu_v0)/(pu_v0+pu_q*xq_sum/pu_v0))*180/3.1415;
float e_q0=	e_Q0*xd_sum/xq_sum+(1-xd_sum/xq_sum)*pu_v0*cos(as0*3.1415/180);
float e_qq0=e_q0*xdd_sum/xd_sum+(1-xdd_sum/xd_sum)*pu_v0*cos(as0*3.1415/180);
float ee0=sqrt((pu_v0+pu_q*xdd_sum/pu_v0)*(pu_v0+pu_q*xdd_sum/pu_v0)+(pu_p*xdd_sum/pu_v0)*(pu_p*xdd_sum/pu_v0));
float ass0=atan((pu_p*xdd_sum/pu_v0)/(pu_v0+pu_q*xdd_sum/pu_v0))*180/3.1415;
float pu_x2=s_x2*m_sb*s_v*s_v/(s_s*vb_1*vb_1);
float pu_x_l0=l_x0*pu_x_l;
float pt=pu_p;
float x_normal=pu_xdd+pu_x_t1+0.5*pu_x_l+pu_x_t2;
float x2_sum=(pu_x2+pu_x_t1)*(0.5*pu_x_l+pu_x_t2)/((pu_x2+pu_x_t1)+(0.5*pu_x_l+pu_x_t2));
float x0_sum=pu_x_t1*(0.5*pu_x_l0+pu_x_t2)/(pu_x_t1+(0.5*pu_x_l0+pu_x_t2));
float x_add=x2_sum*x0_sum/(x2_sum+x0_sum);
float x_short=x_normal+(pu_xdd+pu_x_t1)*(0.5*pu_x_l+pu_x_t2)/x_add;
float pm_short=ee0*pu_v0/x_short;
float x_after=pu_xdd+pu_x_t1+pu_x_l+pu_x_t2;
float pm_after=ee0*pu_v0/x_after;
float as_cr=180-asin(pt/pm_after)*180/3.1415;
float as_lim=acos((load_p/m_sb*3.1415/180*(as_cr-as0)+pm_after*cos(as_cr*3.1415/180)-pm_short*cos(as0*3.1415/180))/(pm_after-pm_short))*180/3.1415;
float k=18000*m_dt*m_dt/(s_tj*s_s/m_sb); 
int n=(int)(m_t/m_dt);
for(int i=0;i<n;i++)
{

}
m_sb=as_lim;
      UpdateData(FALSE);
}