c语言求解有限元中常应变三角形单元分析问题

1,程序设计目的

  巩固和加深对有限元基本理论的理解，掌握程序设计的基本方法,技能及程序语言的编写及编译，提升求解实际工程问题。

2,程序说明

运用有限元方法中三角形常应变单元解平面问题的计算主程序。其步骤如下：

（1）由结点位移及坐标确定出单元内部各点位移。

（2）由单元内部各点的位移确定出对应的应变。

（3）根据应力——应变的关系由单元应变确定出相应的应力。

（4）根据虚功方程由应力确定出结点力。

（5）由单元的结点位移求出单元的结点力，据此求出单元的刚度矩阵。

3,应用实例

源程序：

//求图3-6所示等腰直角三角形单元的刚度矩阵k,设厚度t=1,弹性模量E=1；

#include <stdio.h>
#include <Windows.h>
#include <math.h>
#include <stdlib.h>
#pragma warning(disable:4996)

menu(){
printf("-------------------------------------------------------------------------------\n");
printf("--------------------------Welcome to 有限元程序设计----------------------------\n");
printf("