structure.h

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

// Structure.h: interface for the CStructure class.
//
//////////////////////////////////////////////////////////////////////

#if !defined(AFX_STRUCTURE_H__C4F258F7_BC51_4758_BE85_7FFA911E9409__INCLUDED_)
#define AFX_STRUCTURE_H__C4F258F7_BC51_4758_BE85_7FFA911E9409__INCLUDED_

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

class CStructure : public CObject  
{
public:
	int m_halfband;				//半带宽
	double** m_aStiffMatrix;	//刚度矩阵
	double*	 m_aLoadVector;		//荷载向量
	CString m_title;
	//成员变量========================
	//结构控制参数
	int m_element_num;		//结构单元总数
	int m_node_num;			//结构节点总数
	int m_supportnode_num;	//支承节点数
	int m_loadnode_num;		//有荷载节点数
	int m_question_id;		//问题类型码
	//材料性质参数
	double m_young;			//杨氏弹性模量	
	double m_possion;		//泊松比
	double m_thickness;		//厚度
	double m_weight;		//容重
	//元素对象
	CObArray m_nodeArray;	//结点数组
	CObArray m_elementArray;//单元数组
	//成员函数=========================
	CStructure();
	virtual ~CStructure();
	void Initial();				//初始化所有数据成员
	//-----------------------------------------
	int DataInput(CString &sData);	//取得各项数据
	double ReadDouble(int& iCurCharPos, const CString &sData);	//从指定位置读取一个双精度数据
	int ReadInt(int& iCurCharPos, const CString &sData);		//从指定位置读取一个整数
	int FindKeyWordsPos(CString& sKeyWords, CString& sData);	//在缓冲区中定位字符，返回位置
	//------------------------------------------
	void AnalysisStructure();				//对结构进行分析，取得必要的数据
	int GetHalfBand();				//取得半带宽
	void GetLoadVector(double* loadvector);			//生成荷载向量
	void GetStiffMatrix(double** stiffmatrix);//生成整体刚度矩阵
	void SetSupport(double** stiffmatrix);		//引入支承条件
	void SolveEquationGroup(double** matrixKE, double* vectorP);//解方程组等半带宽存储
	void GetStress();			//求应力，并且存到单元对象中去。
};

#endif // !defined(AFX_STRUCTURE_H__C4F258F7_BC51_4758_BE85_7FFA911E9409__INCLUDED_)