水工结构工程计算中的组合网格算法.docxVIP

水工结构工程计算中的组合网格算法 科学、技术计算中存在许多接触问题。例如，在地质工程中会遇到在几百米的地层中分析几十厘米厚的断层，重力坝和拱坝工程中的坝踵开裂、接缝张开和闭合、边坡稳定等，接触区域常常又是工程的关键部位。对于这种问题，常规有限元计算经常采用的方法是在局部加密以代替整体加密；但是局部加密常常会带来有限元网格剖分十分困难，即便能剖分网格，粗细网格的过渡也会造成节点过多，矩阵病态的问题，使得计算精度降低或计算无法进行。特别是三维问题，在六面体网格方面想要得到一个具有局部加密的高质量的网格还是很不容易的。常规有限元解决混凝土坝的接触问题，网格尺寸通常为10～30 m，一旦一个单元发生拉开或剪切，则10～30 m内的混凝土全部拉开或发生剪切破坏。这对精确描述混凝土坝的接触行为显然是不够的。为此，我们提出了用组合网格法 (cgm) 来解决混凝土坝的接触问题。 1 采用粗网格求解 组合网格算法采用两套网格求解，在整个求解区域采用较粗网格，并且不考虑接触的影响。即认为：在整体区域是连续的，在接缝附近区域采用较细网格考虑接缝的影响，整体粗网格求解和局部细网格求解反复迭代，求得最终结果。其中，引入插值矩阵使得粗细两套有限元网格无需嵌套，并能够适用于非规则网格。即，粗细网格皆可独立生成，彼此互不制约，同时能真正地应用于工程实际问题。在整体区域和局部区域材料参数相同的条件下，许进超和周爱辉教授 1.1 组合网格法的计算格式 为了阐述方便，我们考虑如下的齐次边值问题： 式中，L 问题对应的虚功方程为： 将式（2）中两式相加得 由此可得到如下组合网格法的计算格式： 式中，u 1.2 算法步骤 初始化：n=0, u 步1：求v 步2：求u 2 整体区域收敛控制准则 基于以上讨论，本文选用北京飞箭有限元公司的有限元程序生成系统（FEPG）作为开发工具和平台，开发、设计相应的数值模拟分析程序软件。为了数值分析的高效及可靠性，本文采用模块化设计思路。程序设计分：整体区域不考虑接缝影响计算、局部区域考虑接缝影响计算、局部区域不考虑接缝影响计算、局部区域考虑接缝影响反力计算、局部区域不考虑接缝影响反力计算5个模块。收敛控制准则通过整体区域计算结果来实现。每个相对独立模块间的数据通讯由批处理指令自动执行。 程序计算流程： (1）整体区域不考虑接缝影响计算。 (2）读入（1）的计算结果并插值作为局部区域边界，局部区域考虑接缝影响计算。 (3）读入（1）的计算结果并插值作为局部区域边界，局部区域不考虑接缝影响计算。 (4）读入（2）的计算结果，局部区域考虑接缝影响反力计算。 (5）读入（3）的计算结果，局部区域不考虑接缝影响反力计算。 (6）读入（4）和（5）的计算结果并将（5）-（4）的值通过插值读入（1）中，进入下一流程，直至满足收敛要求。 软件框架（见图2）。 3 计算 3.1 静力中断 3.1.1 基于约束条件的乘子求解 包含lagrange乘子的有限元方程： 式中，A为刚度矩阵；U为位移项；F为右端载荷项；B表示位移与接触力的相关矩阵；λ表示接触力；g为点对中某端受约束的约束值。综合式（7）、（8）得： 对于等式约束，无需修正乘子λ，由式（9）直接求解乘子λ，利用求解得到的乘子λ带入U=A 对于有摩擦的不等式约束，通过整体坐标与局部坐标之间的转换矩阵T，将乘子转到局部坐标系下λ 式中，μ表示摩擦系数；c表示粘聚力；λ 迭代结束后通过λ=T 局部区域连续模块和接缝模块用一套网格，连续模块用等式约束来实现，接缝模块用不等式约束实现。 3.1.2 与组合网格法计算结果比较 计算模型示意见图3。 模型参数：弹性模量E=1.0×10 按照前面接触问题的组合网格算法，在缝的附近取1个3 m×3 m几何模型作为局部区域，整体区域不考虑缝的影响，整体区域粗网格尺寸为0.5 m，局部区域细网格尺寸为0.1 m（如图4、5所示）。 为了验证组合网格法计算结果的合理有效性，将整体加密（网格尺寸为0.1m）计算结果与组合网格法计算结果加以比较（见图6）。 从图6可以看出，两种方法得到节点对的开裂位移结果基本一致，说明组合网格法计算接触问题算法是合理有效的。 通过比较在静力开裂情况下组合网格法与整体加密方法计算的节点对开裂位移，说明组合网格法在计算接触类问题是有效可靠的，为下一步计算大型水工结构中的坝踵开裂、大坝加高后在地震作用下的反应以及横缝接触问题奠定了基础。 3.2 各网格的组合网格法 为了研究丹江口水库加高后新老混凝土接触面的地震反应，在新老混凝土接触面附近用细网格加密作为局部计算区域，而在整体区域上新老混凝土接触面认为是连续的，并不考虑接缝的影响，这样可以在整体区域采用较粗的网格。 整体区域粗网格尺寸在1 m左右，其计算模型及网格示意如图7。局部区