根据变温粘弹性的一般理论及热粘弹性有限元方程,采用20结点三维等参元编写出相应的有限元程序.针对计算中要利用材料的拉伸与扭转的松弛模量函数,而现有的资料只有45＃钢的拉伸松弛与扭转蠕变实验数据,依据材料函数之间的关系,将扭转蠕变函数经拉普拉斯变换与反演变换,推算出相应的松弛函数关系式.另外,根据简单的实验拉伸松弛函数与扭转松弛函数推导出计算中要用的模量函数G<,1>(t)与G<,2>(t).编写了与之相配套的两个程序.一个用于函数曲线的拟合与数值计算,即用一组恒温的应力松弛实验数据拟合出相应的实验曲线并确定变温的应力松弛曲线,由此确定终态温度等效松弛曲线.另一个由蠕变函数计算相应的松弛函数.这样,只要知道了某种材料由实验所得的一组恒温下的应力松弛或蠕变数据,利用该文方法及相应的程序就可以进行该种材料的结构在变温粘弹性状态下的变形及应力的计算与分析.有限元程序中的整体刚度矩阵用一维变带宽存储,采用Crout分块直接求解技术,事先可估算计算时间并得到较高的计算精度,还可一次求解多种工况,以适用于微机解决工程实际中的问题.算例表明计算结果与解析解与实验数据相符.针对几种不同的变温历史分析了厚球壳的应力松弛变化情况.