混凝土高坝仿真应力分析新方法.pdfVIP

混凝1-高坝仿真应力分析新方法 朱伯芳 许 平 (中国水利水电科学研究院) 摘 要 施工过程对坝体应力有重要影响，但混凝土高坝的仿真应力分析需要大量机时和存储容量。 本文作者提出了混凝土高坝仿真应力分析的一套新的计算方法，包括井层算法、并层接缝单元、分 区异步长算法扣考虑水管冷却的等效热传导方程，它们使仿真应力分析的效率得到了极大的提高。 关键词 混凝土坝仿真计算 并层算法并层接缝单元分区异步长算法水管冷却 等效热传导方程 一、前 言 混凝土坝是逐层浇筑的，施工过程对坝体温度场和应力场有重要影响。模拟施工过程进 行仿真计算是困难的，因为①高坝往往有lOO～200层，而且施工期长达数年；②在新浇混 凝土中沿层厚方向的温度梯度和应力梯度都比较大，每一浇筑层内至少要划3～5层单元． 并且需要采用很短的时间步长。为了真实地模拟施工过程并保证较高的计算精度，需要庞大 的存储量和很长的机时，限于目前的计算机硬件，除了巨型机外，一般的计算机难以完成计 算。采用作者提出的一套新的计算方法，在普通微机上即可完成混凝土高坝的仿真计算。 二、并层算法 如图1所示，当坝体逐步升高时，将坝体从上到下划分为四个区域如下： R1：密集单元 在上面新的2～3个混凝土浇筑层中，由于在厚度方向的 温度梯度和应力梯度都比较大，每个浇筑层都分为3～5层单 兀o R2：单层单元 由于在每个浇筑层内，温度增量和应力增量的梯度已趋于 线性化，因此每个浇筑层内只需采用一层单元。 R3：复合并层单元 在区域3中，温度增量和应力增量在铅直方向已进一步趋 于线性变化，但各浇筑层的弹性模量和徐变度仍有一定差异， 把几个浇筑层加以合并，并层后，按复合单元计算，单元内各 图l坝段 层仍保持原来的弹性模量、徐变度和绝热温升，计算单元刚度矩阵[K。]、徐变引起的结点 荷载{pc．}和应力增量{△“}如下(图2a)：。7】 [K。]=…[B]’[西。][B]dxdydz (1) 胛 }△P；}= (2) …【B]T[西。]{％：dxdydz (3) {△“}=[D。](}Ae。}一{‰}一}△Ej}一1△￡!}一{△E；}) 式中}r／．}是与混凝土徐变有关的向量，见文献[7]。 图2并层单Ⅱ (a)复台并层单元；(b)均质并层单元 R4：均质并层单元 在最下面的区域4中，属于老混凝土，各层的材料参数已充分接近，可把若干个浇筑层 合并为一个均质层，并层后的单元采用平均的材料参数，即平均弹性模量E。(r)和平均徐变 度C。(t，r)。 经过这样处理，坝体有限单元的层数，可从100--500层减少到20--30层，计算得到极 大的简化【1，“。 三、并层坝块接缝单元 在常规的混凝土坝施工中，常用纵横接缝把坝体分成许多柱状块体，各坝块的上部为新 混凝土．要采用密集网格，下部可以并层。但采用有限元方法计算时，相邻坝块的网格分层 高度在同一高程最好保持一致，这就限制了并层 的效果。例如高坝块的下部本来可以并层，但因 与相邻坝块的新混凝土接触，必须相应地划分为 许多很薄的单元。另外，在岸坡坝段，基岩本来