（一）程序简介

FEASSI程序是在QUAD-4的基础上改编而成的，它以THEPD程序进行静力计算，为动力计算提供初始条件。

THEPD是清华大学土力学教研组编制的用来分析土体和混凝土结构静力相互作用的平面应变有限元程序，可以计算高土石坝、地下连续墙（包括防渗墙）、廊道等建筑物和地基土体共同作用时的应力和变形。

THEPD设置了5种材料模型，包括线弹性模型、邓肯·张E—B模型、E—μ模型、接触面双曲线模型。FEASSI在此基础之上增加了刚塑性的接触面模型。

THEPD可以模拟分期施工和分期蓄水的实际工况，以不平衡势能原理建立平衡方程，采用变带宽一维存储直接解法，求解效率和求解精度均较高。

FEASSI的主体是动力计算部分，它在QUAD-4的基础之上，做了以下改进。

（1）增加了无厚度的三节点和四节点的接触面单元，并采用刚塑性模型来模拟混凝土与土体的相互接触作用。

（2）可以计算土与混凝土结构系统的自振特性。

（3）增加了重启动功能。在原始数据的读入、迭代计算和绘图等环节，均可以中止和重启动，以利于查错和节省机时。(www.daowen.com)

（4）增加了绘图的后处理功能。

FEASSI包含了三角形常应变单元，任意四边形等参元以及无厚度的接触面单元。整个程序由一个主程序和18个子程序组成，各子程序的主要功能简述如下：

（二）程序流程

FEASSI的程序流程示意图如图5-35所示。

图5-35　FEASSI程序流程示意图