2三维瞬态密封分析基本理论.docVIP

2 三维瞬态密封分析基本理论 核反应堆压力容器运行过程中存在高温高压等复杂工况，采用瞬态耦合热弹塑性接触问题有限元法对法兰进行密封分析是一种精确而有效的方法。为此本章简要介绍弹塑性小变形理论、弹塑性接触理论、稳态及瞬态温度场理论、接触传热理论以及瞬态耦合分析方法。 2.1 密封分析基本方程及迭代格式 压力容器法兰密封分析的基本方程可以表示为 (2.1) 式中 ─ 系统刚度矩阵，它是节点位移向量的函数； ─ 已知载荷向量，包括预紧力和内压载荷； ─ 待定接触力向量，它是接触点对相对位移的函数； ─ 温度载荷，它是瞬态温度场的函数； ─ 密封环反力； ─ 右端项向量。 瞬态热弹塑性接触问题中包含了表面非线性和材料非线性两种非线性，可以在同一增量加载下处理这两类非线性的耦合问题。在每一增量步中外层循环为塑性修正迭代，内层循环为接触迭代。由密封方程可知，外载荷包括外力（已知载荷、温度载荷）和支反力（接触力、密封环反力），外力在迭代过程中保持不变，支反力因受位移的影响必需通过内层循环求得。由于位移增量和密封环反力增量将影响接触状态及接触力增量，而接触状态及接触力增量的变化反过来又会影响增量位移，因此必需进行反复迭代才能求得收敛解[78]。 非线性方程常用的解法有初刚度法、变刚度法、初应力法和初应变法。对压力容器进行弹塑性分析时，由于塑性变形仅出现在密封区域，可能进入塑性的单元极少，且塑性应变较小，同时考虑到下述的接触分析中形成柔度阵需花费大量的机时，多次形成刚度矩阵是不经济的，求解时采取初刚度法较为合理。因此在增量加载及塑性迭代过程中，本文采用修正的Newton-Raphson法求解非线性方程式（2.1），迭代格式为[79] (2.2) (2.3) 式中 ─ 初始切线刚度矩阵，即时的弹性刚度矩阵； ─ 失衡力向量； ─ 接触力增量，在接触迭代中求得； ─ 密封环反力增量，根据压扁值－回弹力曲线求得。 2.2 三维弹塑性接触分析 压力容器在预紧、升压等工况下，法兰接触面附近局部应力集中部位可能出现塑性变形。由于产生塑性变形区域很小且变形量不大，平衡方程及几何方程仍为线性，只有物理方程变为非线性，因此属于小变形弹塑性范畴。 接触问题属于不定边界问题，具有表面非线性，其中既有由接触面积变化而产生的非线性及由接触压力分布变化而产生的非线性，也有由摩擦作用产生的非线性，因此接触问题的求解是一个反复迭代的过程。为了减少运算工作量，本文采用混合求解方法，即首先把整个系统的总刚度方程凝缩到接触边界，形成关于接触内力的柔度方程，整个接触迭代过程中只需要修改柔度矩阵。对于反应堆压力容器，接触区较小，柔度矩阵的阶次远比刚度矩阵的阶次低，故明显提高了迭代效率。 2.2.1 弹塑性接触问题的基本方程 对于具有足够边界位移约束的两个接触物体、，可以写出其基本方程式 （2.4） 式中 、 ─ 、的刚度矩阵； 、 ─ 、的位移向量； 、 ─ 、的外载荷向量； 、 ─ 、的接触力向量。 上式中、、、已知，但、的作用范围、分布规律及大小均为未知，无法求解出、，故需要补充接触定解条件和判定条件。 2.2.2 接触问题混合求解算法 ① 柔度方程的建立 设某一接触点对的位移分别为、，其初始间隙为，则整体坐标系下接触点对的相对间距为 （2.5） 相对间距由外载荷、及接触力、共同产生。其中由、所引起的相对间距，只要设，直接可由式（2.4）求得。由、所引起的相对间距为单位接触力形成的柔度阵与接触力的乘积。 则式（2.5）可以改写成 （2.6） 式中 ─ 柔度矩阵，通过在各接触点对分别作用一对方向相反的单位力，然后求出各接触点对相对间距的方法得到； ─ 外载荷作用下接触点对的相对间距。 接触面局部坐标下增量形式的相对间距方程为 （2.7） 当相对间距为零时，即可得到连续性方程（柔度方程） （2.8） 在求解接触问题时，我们往往首先假定所有接触点对均连续，然后根据各接触点对的实际接触状态引入定解条件作为补充方程[78]，程序中可通过修改柔度方程来实现。 ② 凝缩方程 通常物体可能没有足够的边界位移约束，这时虽然整个系统仍处于稳定平衡状态，但如单独取出物体，则可能出现刚体位移，即该物体的整体刚度矩阵是奇异的，为此必须引入刚体位移公式及整体平衡条件以消除刚度矩阵的奇异性。引入整体平衡条件后的柔度方程为 （