有限元上机报告.docVIP

弹性力学及有限元基础上机实践报告

PAGE

PAGE1

在以常应变三角形单元计算弹性力学平面问题原理的基础上，编写可计算平面应力问题，亦可计算平面应变问题的通用程序。给出程序设计的流程图，程序源代码和程序应用的工程算例及其结果分析。

一、程序流程图

此程序采用本人比较熟悉C语言编写。C语言相比于FORTRAN语言拥有更强大的实用编程环境，同时由于提供了指针变量、动态内存分配函数和结构体变量等，使其编写的软件在维护性、可读性和内存利用效率方面具有明显的优势。

程序总框图

图1程序总框图

整个程序由一个主调主程序(主函数main（）)和七个子程序组成，其中数据子程序DATA()用来接受输入的参数和变量。单刚子程序ELEST()为含有参数的函数，功能控制参数iask可取1、2、3，分别计算相应单元的面积、应力矩阵和单元刚度矩阵，主程序不直接调用它，而通过总刚子程序、载荷子程序和求应力子程序间接调用。总刚子程序TOTSTI()用来合成总刚矩阵，载荷子程序LOAD()用来计算合成载荷，支承子程序SUPPOR()用来引入约束，解方程子程序SOLVEQ()用来求解并输出各个结点的位移，求应力子程序STRESS()用来计算和输出应力、主应力及主平面角。

2.主函数及各子程序流程图

图2voidmain()和voidDATA()

图3voidELEST