plane8nodeele.h

平面8节点单元类

// Plane8NodeEle.h: interface for the CPlane8NodeEle class.
//
//////////////////////////////////////////////////////////////////////

#if !defined(AFX_PLANE8NODEELE_H__85745960_8653_11D5_97B8_8D56FCB44470__INCLUDED_)
#define AFX_PLANE8NODEELE_H__85745960_8653_11D5_97B8_8D56FCB44470__INCLUDED_

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

#include "BaseElement.h"

class CPlane8NodeEle : public CBaseElement  
{
protected:
	int m_aiNode[8];
	double m_dThickness;
	double m_adStressX[9],m_adStressY[9],m_adStressXY[9];
	double m_adMainStress[9],m_adMainStress1[9],m_adMainStressAngle[9];
	double m_adStrainX[9],m_adStrainY[9],m_adStrainXY[9];
	static CMatrix m_matKe;
	static CMatrix m_mat0301;
	static CMatrix m_mat0303;
	static CMatrix m_mat0316;
	static CMatrix m_mat1603;
	static CMatrix m_mat1602;
	static CMatrix m_mat0216;
	static CMatrix m_matNodeDisp;
	static double m_adGaussKxi[9],m_adGaussEta[9];
public:
	virtual void MassAssemble(CSparseMatrix& smatMass);
	void GetMassMatrix();
	void GetMainStress(double* adMainStress,double* adMainStress1,double* adMainStressAngle);
	void GetGaussPtGCoordinate(double* adGaussX,double* adGaussY);
	void GetNode(int *aiNode, int& nEleNode);
	virtual void CalcuGKBandWidth(unsigned long* aiGKDiagAdd);
	virtual void CalcuInternalForce(const double* adNodeDisp);
	virtual void InternalForceInitial();
	virtual void ReadParameter(int &iCurCharPos, CString &sData);
	virtual int GetElementType();
	virtual void OutputInternalForce(ofstream &fout);
	virtual void StiffAssemble(CSparseMatrix& smatGK);
	virtual void OutputParameter(ofstream &fout);
	virtual void GetD(int iGuassPoint, CMatrix &D);
	void GetB(int iGuassPoint,CMatrix& matB,double& dJacobi,double& dWeigh);
	void ReadNode(int& iCurCharPos, CString &sData);
	virtual void GetStiffness();
	CPlane8NodeEle();
	virtual ~CPlane8NodeEle();
	virtual void GetTransferMatrixGCToLC(CMatrix& matT){};
};

#endif // !defined(AFX_PLANE8NODEELE_H__85745960_8653_11D5_97B8_8D56FCB44470__INCLUDED_)