三角形单元有限元计算程序

数据文件共含五个数据段，以－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