混凝土坝坝体与坝基互馈研究.pdfVIP

混凝土坝坝体与坝基互馈研究 向衍 (水利部大坝安全管理中心，江苏 南京210029) 摘要：混凝土坝坝体与坝基系统中，存在人类活动力、岩体应力和地下水渗透力之间的耦合作用。因而，研究 人类活动力、岩体应力和地下水渗流场之问的耦合关系即为本文的互馈作用。本文基于Hertz接触力学理论和 弹塑性理论，考虑坝体与坝基互馈面在外界荷载作用下面积的变化，并借助岩体水力学知识，建立了互馈系统 中渗流场与应力场之间的耦合模型，并基于遗传算法开发了相应的计算程序。 关键词：混凝土坝坝体坝基互馈接触 1互馈的力学行为 1．1互馈面的接触理论 经典接触力学的重要假设之一就是接触表面为几何光滑面，因而表面接触应力分布也是几何光滑的。但 是近代研究表明，固体表面有粗糙度。然而，若表面粗糙度远小于接触区的尺寸，经典接触力学仍然适用， 而当表面粗糙度接近接触区尺寸时，用经典接触力学研究接触问题往往会带来相当大的误差。鉴于实际工程 的需要，近代从经典接触力学中产生了一个新的分支，即粗糙表面接触力学I¨。表面的粗糙性导致非完整接 触，从而使物体间的接触面积仅为表观面积的一部分。研究表明【¨，实际接触面积A，和作用在表面的压缩荷 载尸存在幂律关系，即 P芘群 (1) 1．2互馈面的本构关系 段云岭【21研究了接缝闭合时缝面法向和切向的应力应变关系。南摩擦粘着理论，互馈面摩擦力以及互馈 面之间相互变形的大小与两种互相接触材料的剪切变形特性有直接关系。闭合以后互馈面的切向刚度除了与 切向应变相关以外，还与互馈面的压力相关。现以有限元分析中某个互馈面上单元的一个应力Gauss点及其 附近区域进行分析。假设，在该Gauss点附近的一个较小的区域S内，互馈面法向压应力一致，则在该区域 内有一个实际接触面积S’，而S’心是法向压应力盯的函数。仃越大，S／S越大，其最大值为1，即s，-S。在该 实际接触区域y中，最大剪切应力‰ax和峰值剪切应变而的关系为 ‰=既 (2) 式中：G为较软材料峰值剪切应变对应的割线剪切刚度，可根据峰值剪应变由剪切应力应变全曲线求得，通 常小于初始剪切刚度Go。 ， r=ftr+c (3) 其中：f为摩擦系数；C为黏结强度。 区域S上抗剪力为 F=rS (4) 事实上，抗剪力F由‰。在区域S’上积分得到，对于同样的材料，钿。。相同，故由式(2)得 f=GYvs,／s (5) 由于在S上出现最大抗剪力时，S上所有点的名义剪切应变都应该是肠，即S上名义剪切应力为 基金项目：国家973项目(2002CB412707)和国家自然科学基金重点项目 ·339· f=G’儿 其中，G’为s上的名义剪切刚度．，即 ． GI-cs／s (7) 有限元分析的对象是区域s而不是区域∥，所得到的剪切应力是名义剪切应力而不是最大剪切应力‰。。， 因此有限元分析中应采用名义剪切刚度G7而不是G。本文研究互馈面实际接触面积的影响。 ‘ 1．3互馈面的应力场计算 互馈面单元用空间八节点薄层单元，坝基岩体和坝体用空间八节点块体单元模拟，薄层取为不抗拉(或低 抗拉