beam.h

一个计算悬臂梁的有限元vc源码

// Beam.h: interface for the CBeam class.
//
//////////////////////////////////////////////////////////////////////

#if !defined(AFX_BEAM_H__77A82520_5F23_11D5_97B8_BF31DF29CA76__INCLUDED_)
#define AFX_BEAM_H__77A82520_5F23_11D5_97B8_BF31DF29CA76__INCLUDED_

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

#include "BaseElement.h"

class CBeam : public CBaseElement  
{
public:
	virtual void MassAssemble(CSparseMatrix& smatMass);
	virtual void GetMassMatrix();
	virtual void OutputParameter(ofstream& fout);
	virtual void OutputInternalForce(ofstream& fout);
	virtual int GetElementType();
	double GetInertia();
	virtual void ReadParameter(int& iCurCharPos, CString& sData);
	void ReadNode(int& iCurCharPos, CString &sData);
	void SetInitialSectionForce(double *adQInit, double *adMInit);
	void GetSectionForce(double& dN,double *adQ, double *adM);
	void SetInitialSectionDef(double *adFyInit);
	void GetSectionDef(double *adFy,double&dU,double&dU1);
	virtual double GetLength();
	virtual void GetTransferMatrixGCToLC(CMatrix& matT);
	virtual void GetNode(int *aiNode, int& nEleNode);
	virtual void InternalForceInitial();
	virtual void CalcuInternalForce(const double* adNodeDisp);
	virtual void GetStiffness();
	virtual void CalcuGKBandWidth(unsigned long* aiGKDiagAdd);
	virtual void StiffAssemble(CSparseMatrix& smatGK);
	CBeam();
	virtual ~CBeam();
protected:
	static CMatrix m_matT;
	static CMatrix m_matTT;
	static CMatrix m_matKe;
	static CMatrix m_matNodeDisp;
	double m_dArea;
	double m_dInertia;
	int m_aiNode[2];
	double m_adM[11],m_adQ[11],m_adFy[11],m_dU,m_dU1;
	double m_dN;
	double m_dLength;
	double m_dSin,m_dCos;
};

#endif // !defined(AFX_BEAM_H__77A82520_5F23_11D5_97B8_BF31DF29CA76__INCLUDED_)