编辑：杭州那泰科技    www.nataid.com

有限元分析时利用式可以得到如下单盘板的载荷分布：进一步可以得到浮顶内外液面高度hw和hs分别为从式可以看出，浮顶内外液面高度hw与hs以及单盘板载荷分布q(r)都与单盘板挠度f(r)有关，且一般情况下q(r)随f(r)沿单盘板径向变化而呈非均匀分布。在沉没试验状态下有 (即灌内及单盘上均装水)，并取R=R1(相当于整个浮顶上承受h0的名义降雨量)，则q(r)简化为此时q(r)为均布载荷(亦称为当量载荷)，但与单盘板的变形有关(因hc与f(r)有关)；浮顶内外液面高度hw和hs分别简化为由式(9)和式(11)可以看出，浮顶所受载荷和单盘板的挠度是相关的，且二者是一个非线性的关系，本文通过载荷迭代修正的方法获得这一问题的解。

求解的基本思路是：首先假设单盘板有一个初始挠度(一般初始状态取为单盘板没有发生变形的状态，即f(r)=0)，然后根据这一挠度求出载荷分布并作为初始载荷利用有限元方法对浮顶结构进行几何非线性分析；得到新的挠度解后对载荷进行修正并开始下一次迭代分析，直到根据挠度修正后的载荷与本次迭代的初始载荷之差满足精度要求时为止。这一分析过程可以概括如下：

(1)根据浮顶结构的力学模型建立有限元分析模型，包括选取单元类型(如用壳单元和梁单元分别模拟板和梁等结构)、划分网格、施加载荷等。

(2)根据式(7)分别计算Hg和ΔH。

(3)取单盘板未发生变形时为初始状态，即f(r)=0。

(4)根据式(4)计算hc(可采用普通的数值积分方法求得，初始状态时hc=0)。

(5)根据式(9a)或式(9b)计算单盘板上的等效载荷分布q(r)。

(6)根据式(11a)或式(11b)计算浮舱底面载荷分布。

(7)在有限元模型中进行加载求解。

(8)根据有限元求解得到的挠度f(r)按照上述步骤(4)～(6)修正q(r)和pb。

(9)检查本次计算修正结果和上一次计算修正结果是否一致，即载荷q(r)是否满关系式 (其中i为计算修正次数，ε为给定的精度要求)。