readme.txt

一个三角形单元有限元计算程序,是Dev C++的源代码,

数据文件共含五个数据段，以－1另起一行分隔，格式如下：

1. 数据常量
    依次为： 泊松比Mu，弹性模量E，厚度t，计算类型标志符flag，节点数nNum，单元数cNum
    （标志符为0则为平面应力问题，为1则为平面应变问题）    

2. 节点编号及坐标值
    依次为： 节点编号，横坐标x，纵坐标y
    （注意：节点编号自0开始）

3. 单元编号及其节点
    依次为： 单元编号，单元所含三节点编号（逆时针次序排列）
    （注意：单元编号自0开始）

4. 等效荷载
    依次为： 节点编号，等效荷载Fx，Fy
    （同一节点在不同单元的等效荷载可分别输入，程序自动进行累加）

5. 支撑条件
    依次为： 节点编号，x方向位移u

6. 支撑条件
    依次为： 节点编号，y方向位移v

7. 最后以多个－1结束


文件dat39.txt给出的是教材167页例3.9的算例。

文件dada.txt中给出了一个4节点2单元有限元结构的实例，数据如下：

0.0 	1.0e6	1.0	0	4	2	
-1

0	0.0	0.0
1	2.0	0.0
2	2.0	1.0
3	0.0	1.0
-1

0	0	1	3
1	1	2	3
-1

1	0.0	-0.5
2	0.0	-0.5
-1

0	0.0
3	0.0
-1

0	0.0
3	0.0
-1	-1	-1	-1	-1	

输出结果如下：


Node number: 0  (0, 0)  u = 0  v = 0
Node number: 1  (2, 0)  u = -1.38776e-006  v = -6.61224e-006
Node number: 2  (2, 1)  u = 1.38776e-006  v = -6.93878e-006
Node number: 3  (0, 1)  u = 0  v = 0

Elem number: 0    included nodes: 0  1  3	area = 1
	sgmX = -0.693878	sgmY = 0	tauXY = -1.65306
	sgmMax = 1.34214	sgmMin = -2.03601	theta = -39.0735
Elem number: 1    included nodes: 1  2  3	area = 1
	sgmX = 0.693878	sgmY = -0.326531	tauXY = -0.346939
	sgmMax = 0.800662	sgmMin = -0.433315	theta = -72.8921