elem1.h

三角单元平面静力问题

// elem1.h: interface for the elem class.
//
//////////////////////////////////////////////////////////////////////

#if !defined(AFX_ELEM1_H__79166E59_A0C4_4789_BBA0_25C7C4C035F5__INCLUDED_)
#define AFX_ELEM1_H__79166E59_A0C4_4789_BBA0_25C7C4C035F5__INCLUDED_

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000
#include "node.h"
#include "matrix.h"
#include<cstdlib>
class Elem
{
	

public:
	Node* ndp[3];         //节点指针 
	
	int elNum;
	
	double sgmX;   //单元应力 
	double sgmY;
	double tauXY;
	double sgmMax;   //主应力 
	double sgmMin;   //
	double theta;    //应力主角 
	
	double area;      //单元面积
 	double a[3];        //计算中间变量 
  	double b[3];
  	double c[3];
  	
  	static double Mu;   //泊松比 
  	static double E;    //弹性模量 
  	static double t;    //厚度 
  	static int ndTotNum;//总节点数 
  	static int elTotNum;//总单元数   	
  	
  	void getAbc();      //计算变量数组a[],b[],c[]及面积elArea  
  	
public:
	Elem();
	~Elem() { };
    
	void setElem(int, Node*, Node*, Node*);                 //设置成员变量值 
	void setMtxk(Mtx<double>&);     //生成整体刚度矩阵
	void setStress();                  //计算单元应力 
	void output();
	static void setStatics(double, double, double, int, int, int);
 	static int getNdTN() { return ndTotNum; }; 
  	static int getElTN() { return elTotNum; };	
	
	//friend ostream& operator<<(ostream&, const Elem&);
};




#endif // !defined(AFX_ELEM1_H__79166E59_A0C4_4789_BBA0_25C7C4C035F5__INCLUDED_)