element.h

自编的有限元方法解决弹性力学平面问题的软件。有图形显示

// Element.h: interface for the CElement class.
//
//////////////////////////////////////////////////////////////////////

#if !defined(AFX_ELEMENT_H__FB8202F6_2B26_43DE_AF6D_74AFB4E4CF18__INCLUDED_)
#define AFX_ELEMENT_H__FB8202F6_2B26_43DE_AF6D_74AFB4E4CF18__INCLUDED_

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000

#include "Node.h"
class CElement : public CObject  
{
public:		//成员变量
	//单元材料性质参数
	double m_young;		//杨氏弹性模量E
	double m_possion;	//泊松比u
	double m_thickness;	//单元厚度t
	double m_weight;	//容重
	//单元结构参数
	double m_area;		//单元面积
	CNode* m_pNode[3];	//单元结点
	//应力问题,最终要求
	double m_sigmaX;		//X应力分量
	double m_sigmaY;		//Y应力分量
	double m_taoXY;			//剪切应力分量
	double m_averageStress;	//平均应力
	double m_radiusStress;	//应力圆半径
	double m_maxMainStress;	//最大主应力
	double m_minMainStress;	//最小主应力
	double m_mainPlaneAngle;//主平面角	
public:		//成员函数
	void GetStress();					//求应力
	void GetElementStiffMatrix(double** elementStiffMatrix);//求单元刚度矩阵
	void GetMatrixS(double** matrixS);	//求得应力位移转换矩阵S
	void GetMatrixB(double** matrixB);	//求得应变位移转换矩阵B
	void GetMatrixD(double** matrixD);	//求得弹性矩阵D
	void Initial();
	void SetMaterialInfo(double young,double possion,double thickness,double weight);
	void SetNode(CNode* pNode1, CNode* pNode2, CNode* pNode3);
	int m_element_id;
	double GetArea();		//求单元的面积
	CElement(int id);				//构造函数
	CElement();
	virtual ~CElement();	//析构函数

};

#endif // !defined(AFX_ELEMENT_H__FB8202F6_2B26_43DE_AF6D_74AFB4E4CF18__INCLUDED_)