复杂水闸闸室结构分析方法研究.docVIP

复杂水闸闸室结构分析方法研究 （扬州大学 水利科学与工程学院，江苏 扬州 225009） 摘 要：本文主要介绍采用Marc软件求解复杂水闸闸室结构在固定水位下的空间位移和主应力，利用程序提供的接触分析功能，构造闸室结构、灌注桩和地基结构的三维有限元模型来进行计算分析，得出整体结构的位移场和应力场，判断其整体强度和稳定性，并对在采用不同水闸底板结构情况下的水闸工作性态进行分析评价。实例分析表明，对于复杂的闸室结构，利用三维有限元分析可更直观、精确反映建筑物整体强度与位移状况，采用整体式结构情况下水闸适应地基不均匀沉降能力较强，采用分离式结构情况下底板更能体现其结构的优越性。 关键词 有限元 水闸作为水工建筑物的重要组成部分之一，近年来随着工农业生产的发展也发展迅速，水闸建设积累了丰富的经验，其结构越来越多样化。在复杂水闸闸室结构中，上下游段底板及闸墩可均为不同结构，且结构不同的闸墩均与上部撑梁相连，下游段柱状闸墩间也由撑梁连接，底板与闸墩间可以采用整体式结构浇筑，也可以采用分离式结构，设置施工缝。若闸墩下部有灌注桩，则反拱底板将不承受上部荷载作用，所有上部荷载重量由灌注桩承受。对于这样的水闸结构受力复杂，目前还没有精确的理论分析方法以进行计算。下的地基反力分布与底板的尺寸、刚度、弹性模量、荷载大小分布等因素地基的变形特性密切相关，在结点上，作用有两种力，其一为变形单元对结点的作用力；其二为作用于结点的外荷载。根据平衡条件，作用于结点上的合力应为零，即在结点，所有相邻单元的点结点力之和与结点荷载平衡，列出全部结点的平衡条件，得出总体平衡方程组，即： （1） （2） 式中：为总体刚度矩阵，由单元刚度矩阵组合而成，是一个稀疏的正定对称矩阵；为全部结点的位移列向量；为全部结点的荷载列向量，包括直接作用于结点上的集中荷载和所有分布荷载、初应力或初应变转换而成的等效结点荷载；称为单元刚度矩阵，它的元素表示当该单元发生一定的结点位移时所对应的结点力，单元刚度矩阵决定于该单元的形状、大小、方向和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平移而改变。 在Marc程序中，对于接触面的处理，一般是通过定义接触体和接触表来描述物体间的接触关系，接触关系类型包括接触、不接触和粘合三种。在接触表中定义接触提之间的摩擦系数、接触后分离所需的分离力、接触容差及可能的过盈配合值[1]。对于整体式底板结构，定义底板与闸墩间接触为粘合，对于分离式底板结构则定义底板与闸墩间的摩擦系数为0。在复杂水闸闸室结构中，通过无相对滑动的Glue功能和施加很大的分力把地基结构粘连起来如一个整体，之间无相对滑动速度，接触表面变形满足连续条件。离散成八节点六面体单元等参单元，单元之间通过有限个点连接起来[5]。 3 实例分析 本文结合江苏省黄楼闸实例进行有限元计算分析研究，并分析比较该水闸底板采取整体式和分离式两种结构下，水闸整体位移及应力情况。 3.1 概况 黄楼闸位于复新河上游宋楼镇境内，是丰县复新河流域梯级水面三级控制工程，承担116.8km2防洪和排涝任务。该闸建成于1982年6月，规模为中型，闸室为开敝式，闸门共5孔，每孔净宽4.2m，闸门为钢丝网双曲扁壳结构，启闭机为2×8T手电两用螺杆式，设计行洪流量122.6 m3/s，设计上游蓄水位40.0m。交通桥设在闸室靠下游一侧，桥面净宽7m。 （一）黄楼闸结构计算模型 取黄楼闸整个底板、闸墩和一定范围内的地基(沿水平方向距底板两端取24.0m，沿深度方向取至25.5m)作为计算结构模型，剖分网格[6]。闸顶以上排架、启闭机等设备作为荷载作用在闸墩顶上，拱内混凝土填料作为荷载作用在底板顶面上，墙后填土作为土压力作用在墙背和底板上，周边填土作为边荷载作用在地基上。由于考虑到了地基模型的尺寸范围的选择，故采用全约束[7]。坐标系取为：零点设置在底板面和闸墩的接触处，x轴垂直于水流方向指向右边，y轴垂直指向上方，z轴顺水流方向指向下游。图1为黄楼闸结构空间有限元网格图。该黄楼闸生成131786个八节点六面体等参单元，24根灌注桩共生成44880个八节点六面体等参单元，地基划分成38454个八节点六面体等参单元。  图0.05g，相当于地震基本烈度6度。可不进行地震计算。 对于黄楼闸计算模型取以下合计固定荷载+回填土荷载+设计水位荷载+浪压力 对黄楼闸闸室结构进行空间有限元计算，求出了各种工况下结构在荷载作用下的各点位移和应力。由此分析闸室位移与应力变化规律。 3.3计算结果分析及比较 分别对整体式和分离式底板模型进行材料参数的设置和边界条件的加载，进而对其位移场和应力场比较，得出两种结构的不同特征。 在整体式底板情况下，黄楼闸首在设计工况下的位移分布图示意如图2、3、4。 图图 图