四接点两单元的平面静力问题dlg.cpp

三角单元平面静力问题

// 四接点两单元的平面静力问题Dlg.cpp : implementation file
//

#include "stdafx.h"
#include "四接点两单元的平面静力问题.h"
#include "四接点两单元的平面静力问题Dlg.h"
#include "elem1.h"
#include "matrix.h"
#include "node.h"
#include "math.h"
#include<cstdlib>
#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()

/////////////////////////////////////////////////////////////////////////////
// CMyDlg dialog

CMyDlg::CMyDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CMyDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CMyDlg)
	m_bosongbi = 0.0;
	m_moliang = 500000;
	m_houdu = 1;
	m_x1 = 0.0;
	m_x2 = 2.0;
	m_x4 = 0.0;
	m_x3 = 2.0;
	m_y1 = 0.0;
	m_y2 = 0.0;
	m_y3 = 1.0;
	m_y4 = 1.0;
	m_fx1 = 0.0;
	m_fx2 = 0.0;
	m_fx3 = 0.0;
	m_fx4 = 0.0;
	m_fy1 = 0.0;
	m_fy2 = -0.5;
	m_fy3 = -0.5;
	m_fy4 = 0.0;
	m_u1 = 0.0;
	m_u2 = 1;
	m_u3 = 1;
	m_u4 = 0.0;
	m_v1 = 0.0;
	m_v2 = 1;
	m_v3 = 1;
	m_v4 = 0.0;
	m_x5 = 0.0;
	m_y5 = 0.0;
	m_x6 = 0.0;
	m_y6 = 0.0;
	m_fx5 = 0.0;
	m_fy5 = 0.0;
	m_fx6 = 0.0;
	m_fy6 = 0.0;
	m_u5 = 0.0;
	m_v5 = 0.0;
	m_u6 = 0.0;
	m_v6 = 0.0;
	m_ch1 = true;
	m_ch2 = FALSE;
	m_ch3 = FALSE;
	m_ch4 = true;
	m_ch5 = FALSE;
	m_ch6 = FALSE;
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CMyDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CMyDlg)
	DDX_Control(pDX, IDC_COMBO2, m_enum);
	DDX_Control(pDX, IDC_COMBO1, m_flag);
	DDX_Text(pDX, IDC_EDIT35, m_bosongbi);
	DDX_Text(pDX, IDC_EDIT36, m_moliang);
	DDX_Text(pDX, IDC_EDIT37, m_houdu);
	DDX_Text(pDX, IDC_EDIT1, m_x1);
	DDX_Text(pDX, IDC_EDIT2, m_x2);
	DDX_Text(pDX, IDC_EDIT4, m_x4);
	DDX_Text(pDX, IDC_EDIT3, m_x3);
	DDX_Text(pDX, IDC_EDIT5, m_y1);
	DDX_Text(pDX, IDC_EDIT6, m_y2);
	DDX_Text(pDX, IDC_EDIT7, m_y3);
	DDX_Text(pDX, IDC_EDIT8, m_y4);
	DDX_Text(pDX, IDC_EDIT19, m_fx1);
	DDX_Text(pDX, IDC_EDIT20, m_fx2);
	DDX_Text(pDX, IDC_EDIT21, m_fx3);
	DDX_Text(pDX, IDC_EDIT22, m_fx4);
	DDX_Text(pDX, IDC_EDIT23, m_fy1);
	DDX_Text(pDX, IDC_EDIT24, m_fy2);
	DDX_Text(pDX, IDC_EDIT25, m_fy3);
	DDX_Text(pDX, IDC_EDIT26, m_fy4);
	DDX_Text(pDX, IDC_EDIT27, m_u1);
	DDX_Text(pDX, IDC_EDIT28, m_u2);
	DDX_Text(pDX, IDC_EDIT29, m_u3);
	DDX_Text(pDX, IDC_EDIT30, m_u4);
	DDX_Text(pDX, IDC_EDIT31, m_v1);
	DDX_Text(pDX, IDC_EDIT32, m_v2);
	DDX_Text(pDX, IDC_EDIT33, m_v3);
	DDX_Text(pDX, IDC_EDIT34, m_v4);
	DDX_Text(pDX, IDC_EDIT10, m_x5);
	DDX_Text(pDX, IDC_EDIT17, m_y5);
	DDX_Text(pDX, IDC_EDIT18, m_x6);
	DDX_Text(pDX, IDC_EDIT38, m_y6);
	DDX_Text(pDX, IDC_EDIT39, m_fx5);
	DDX_Text(pDX, IDC_EDIT40, m_fy5);
	DDX_Text(pDX, IDC_EDIT41, m_fx6);
	DDX_Text(pDX, IDC_EDIT42, m_fy6);
	DDX_Text(pDX, IDC_EDIT43, m_u5);
	DDX_Text(pDX, IDC_EDIT44, m_v5);
	DDX_Text(pDX, IDC_EDIT45, m_u6);
	DDX_Text(pDX, IDC_EDIT46, m_v6);
	DDX_Check(pDX, IDC_CHECK1, m_ch1);
	DDX_Check(pDX, IDC_CHECK2, m_ch2);
	DDX_Check(pDX, IDC_CHECK3, m_ch3);
	DDX_Check(pDX, IDC_CHECK4, m_ch4);
	DDX_Check(pDX, IDC_CHECK5, m_ch5);
	DDX_Check(pDX, IDC_CHECK6, m_ch6);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CMyDlg, CDialog)
	//{{AFX_MSG_MAP(CMyDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_CBN_SELCHANGE(IDC_COMBO1, OnSelchangeCombo1)
	ON_CBN_SELCHANGE(IDC_COMBO2, OnSelchangeCombo2)
	ON_BN_CLICKED(IDC_CHECK1, OnCheck1)
	ON_BN_CLICKED(IDC_CHECK2, OnCheck2)
	ON_BN_CLICKED(IDC_CHECK3, OnCheck3)
	ON_BN_CLICKED(IDC_CHECK4, OnCheck4)
	ON_BN_CLICKED(IDC_CHECK5, OnCheck5)
	ON_BN_CLICKED(IDC_CHECK6, OnCheck6)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CMyDlg message handlers

BOOL CMyDlg::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
	m_flag.SetCurSel(0);
    m_enum.SetCurSel(0);
	GetDlgItem(IDC_EDIT10)->EnableWindow(false);
	GetDlgItem(IDC_EDIT17)->EnableWindow(false);
	GetDlgItem(IDC_EDIT39)->EnableWindow(false);
	GetDlgItem(IDC_EDIT40)->EnableWindow(false);
	GetDlgItem(IDC_EDIT43)->EnableWindow(false);
	GetDlgItem(IDC_EDIT44)->EnableWindow(false);

	GetDlgItem(IDC_EDIT18)->EnableWindow(false);
	GetDlgItem(IDC_EDIT38)->EnableWindow(false);
	GetDlgItem(IDC_EDIT41)->EnableWindow(false);
	GetDlgItem(IDC_EDIT42)->EnableWindow(false);
	GetDlgItem(IDC_EDIT45)->EnableWindow(false);
	GetDlgItem(IDC_EDIT46)->EnableWindow(false);
	
	GetDlgItem(IDC_EDIT28)->EnableWindow(false);
	GetDlgItem(IDC_EDIT32)->EnableWindow(false);
	GetDlgItem(IDC_EDIT29)->EnableWindow(false);
	GetDlgItem(IDC_EDIT33)->EnableWindow(false);
	
	kk=0;
	ek=0;
	return TRUE;  // return TRUE  unless you set the focus to a control
}

void CMyDlg::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 CMyDlg::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 CMyDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

void CMyDlg::OnOK() 
{
	// TODO: Add extra validation here
	
	nptr = 0; //节点指针 
	eptr = 0; //单元指针 
	kptr = 0; //整体刚度矩阵指针，初始值为0 	
	
	readData(nptr, eptr, kptr);//读取节点数、单元数及常量Mu, E, t 
	solveEq(nptr, kptr);
	stress(eptr);
 	MessageBox(output(nptr, eptr));
	
	//eptr[1].output();
 	UpdateData(false);
  	delete[] nptr;
 	delete[] eptr;
 	delete kptr;
	
}
void CMyDlg::readData(Node*& ndp, Elem*& elp, Mtx<double>*& kp)
{
    UpdateData(true);
  	
	Elem::setStatics(m_bosongbi, m_moliang, m_houdu, kk, ek+4, ek+2); 

	const int nNum=Elem::getNdTN();
 	const int eNum=Elem::getElTN();
	
	ndp = new Node[nNum];                  //定义节点数组 
	elp = new Elem[eNum];                  //定义单元数组 
	kp = new Mtx<double>(2*nNum, 2*nNum);  //定义整体刚度矩阵 
  	
   // CString ss;
	//ss.Format("%d",nNum);
	//MessageBox(ss);
	// 读取节点数据 
     double x[6],y[6];
	 {x[0]=m_x1;
	  x[1]=m_x2;
	  x[2]=m_x3;
	  x[3]=m_x4;
	  x[4]=m_x5;
	  x[5]=m_x6;
	 
	  y[0]=m_y1;
	  y[1]=m_y2;
	  y[2]=m_y3;
	  y[3]=m_y4;
	  y[4]=m_y5;
	  y[5]=m_y6;
	 }
    for(int i=0;i<nNum;i++)
	 {  ndp[i].setPoint(i, x[i], y[i]);
		
	 }	
		  //读取单元数据 
	
	     for(int q=0;q<eNum;q++)