外文翻译译文.docVIP

丹江口重力坝后期加高的安全稳定性研究 周伟，常晓林，卢淑媛 （武汉大学 水资源与水利水电学院，中国，湖北，武汉430072） 摘要：通过采用以传统库仑摩擦模型和相应的接触元件为基础的一种大位移改善接触摩擦模型，本文对丹江口重力坝后期加高的新老混凝土结合面进行接触分析，以便研究开裂后大坝的稳定和变形状态。计算结果表明，接触面的摩擦应力和正应力的分布及大小符合一般规律，沿坝基面的抗剪断安全系数也满足重力坝设计规范要求。 关键词：安全稳定性研究;接触模型;增广拉格朗日乘子法。 1.概述 为了满足南水北调中线工程，丹江口（DJK）重力坝从162米加高至177.6米。通过温度场的模拟，新老混凝土。由于温度变化。开裂后，研究稳定性和变形状态，本文假定坝建设后缝隙，然后采用接触摩擦模型进行新老混凝土的接触分析。只有在研究整个安全能力和变形状态设计师明白如何重力坝新老混凝土，以及是否有足够的安全能力，以抵抗增加水的反作用力。 2接触模型  接触问题是高度非线性的，需要的计算机资源来解决。我们理解这个问题的物理和花时间去建立尽可能高效运行是非常重要的。接触问题两个显著的困难。首先，我们通常不知道接触的区域，直到我们已经遇到问题。根据载荷，材料，边界条件和其他因素的影响，可以和离以不可预知的和突然的方式。，接触问题需要考虑摩擦。有几种摩擦法和可供选择，并且都是非线性的。摩擦是，使得求解收敛困难。为了模拟一个接触问题，我们首先需要确定其可能的相互作用分析。如果该相互作用中的一是在一个点上，模型对应的部件是一个节点。如果该相互作用中的是在一个表面上，模型中的对应的部件是一个元素：一个梁，壳，或固态元件。本文只讨论面面接触问题和模型，其中有限元模型特定的存在识别可能的接触对。这些接触元件被叠加在正被用于相互作用模型的各个部分 在基本的库仑摩擦模型中，两个接触面可以剪应力在一定。。库仑摩擦模型（见图1 ）被定义为： (1) 其中是极限剪切应力为摩擦系数，是等效的剪应力， p为接触常压。一旦等效剪应力超过 ，和目标表面会相对于彼此滑动。这 态被称为滑动。/滑动计算决定到滑动反之亦然。  摩擦是一个复杂的物理现象，它涉及到的特性如表面粗糙度，相对速 度，应力，温度。摩擦的实际物理仍然是一个研究课题。因此，摩擦的数值模拟已经被简化为理想化模型。最流行的摩擦模型是粘着摩擦或库仑模型。本文采用了一种传统的库仑摩擦模型基础上改进的接触摩擦模型，该模型是： 其中是正应力是切向（摩擦）力，是摩擦系数， RV被理想化的滑动速度，t是在相对速度的方向上的切向矢量 是相对滑动速度。 物理上，RV值是滑动发生时的相对速度值数学模型如何密切。 一个非常大的RV有效摩擦。一个非常小的值会导致收敛不佳。建议RV值应为典型相对滑动速度1 ％或10 ％。正因为如此平滑过程，接触总是有一些打滑。摩擦模型的反正切表示这个非线性摩擦行为数学理想化，这使得解决方案收敛。  可以使用非线性弹簧进行改善接触摩擦模型的实施注意一个非线性弹簧的行为，如图2所示，由下式给出： 其中K为弹簧刚度是点1和2的位移和力，在速度方面很容易 (5) 图二 弹簧元件 由于K是相对速度的一个非线性函数， 被逐步得到解决，其中，每个增量内的迭代次数可能是必要的。对于一个典型的迭代i，则方程解决 (7) 注意对应于迭代的开始。对于变形的刚性接触，很容易看到 (8) 在静态分析，速度必须位移和时间，我们可以这样写： 其中表示位移增量迭代的修正（另见公式7） ： 式8和9代入到6中得到 (11)通过考虑先前增量方程ll第一次迭代然后9可以改写为： 可以修改 对于后续的迭代将到11 。在方程13表示在以前的增量结束时，点1和2之间的相对速度。必须注意的是，如果所述刚性体的速度相对速度差别极大13的附加术语特别重要。 图三 可能接触体之间的正常间隙 （17） 我们得到 （18） 其中K是刚度矩阵，C是约束矩阵