cele.h

有限元计算,可以计算一型梁单元的软件. 我们可以共同提高

#include "CMatrix.h"

class CEle
{
private:

public:
	CEle();
	CEle(double dL,double dB,double dH,double dQ,int nN1,int nN2,double dE,double dMu);
	~CEle();
	void GetN(double mx,double my);
	void GetB(int IM,double mx,double my);
	void GetD();
	void GetKe(int IM);
	void GetPe(int IM);
	void GetK();
	void GetDise(int IM);
	void GetSe(int IM);
	void GetS();
	void GetJ(int IM,double mx,double my);
	void Solve();
	void Solve2();
	void GetN2(double mx,double my);
	void SetUp(double dL,double dB,double dH,double dQ,int nN1,int nN2,double dE,double dMu);
	void restrict(char flag1,char flag2);	

	CMatrix m_mCOORD;//结点坐标
	CMatrix m_mMP;//单元的结点编号
	CMatrix m_mKe;//单元刚度矩阵
	CMatrix m_mPe;//单元的结点力
	CMatrix m_mP;//整体结点力矩阵
	CMatrix m_mK;//整体刚度矩阵
	CMatrix m_mDise;//单元位移矩阵
	CMatrix m_mDis;//整体位移矩阵
	CMatrix m_mSe;//单元应力
	CMatrix m_mS;//整体应力矩阵
	CMatrix m_mB;//应变矩阵
	CMatrix m_mD;//弹性矩阵
	CMatrix m_mN;//型函数矩阵
	CMatrix m_mML;//结点荷载矩阵
	CMatrix m_mJ;//jacob矩阵
	CMatrix m_mN2;//四结点型函数矩阵

	int m_nNJ;//总结点数
	int m_nNM;//总单元数
	double m_dDetJ;//矩阵的绝对值
	double m_dE;//弹性模量
	double m_dMu;//泊淞比
	int m_nN1;//划分的格数
	int m_nN2;//竖向格数
	double m_dL;//长度
	double m_dH;//高度
	double m_db;//宽度
	double m_dQ;//均布力
};