基于Ogden公式的充压伸缩式水封非线性计算研究.docxVIP

3基于Ogden 公式的充压伸缩式水封非线性计算研究33刘礼华王蒂欧珠光魏晓斌( 武汉大学土木建筑工程学院 ,武汉 ,430072)摘要 考虑橡胶类止水材料性能对水封设计的重要作用 ,论文采用 Moo ney2Rivli n 公式和 Ogden 公式对橡胶类止水材料的单向拉伸 、单轴压缩和纯剪实验数据进行曲线拟合与研究 . 结合工程实例将该两种公式的拟合参数应用于充压伸缩式水封自由外伸量的非线性有限元计算中 ,并将计算结果与该水封自由外伸量模型实验进行对 比研究 . 结果表明 :Ogden 公式比 Moo ney2Rivlin 公式能更好地反映橡胶类止水材料的性能 ,在闸门水封设计中应 用 Ogden 公式能更有效地提 高其非线性有限元计算 的 精 度 , 这 为 闸 门 水 封 结 构 的 设 计 计 算 提 供 了 有 益 的 参 考依据 .关键词 充压伸缩式水封 ,橡胶类止水材料 ,非线性有限元 ,Ogden 公式0 止水材料本构关系水封是闸门最重要的部件之一 ,其性能的优劣 直接影响闸门止水效果 ,特别是随着我国高水头水 电站的不断开发 ,水封的止水问题已成为高水头闸门设计的关键性技术问题 . 在水封设计与研究过程 中 ,一般要遵循计算机理论计算与模型实验两个环 节 ,即从水封止水元件的材料实验与其止水过程的 非线性有限元计算入手 ,完成较优水封材质及水封 结构型式的选择 ,然后进行水封模型实验验证 ,以最经济省时的方法得到一种止水高效 、经久耐用的水 封. 其中水封非线性有限元计算的可靠性依赖于橡 胶类止水材料性能描述方法的准确性. 而橡胶类止 水材 料 的 不 可 压 缩 性 、大 变 形 、超 弹 性 等 力 学 特 点[ 1 ,2 ] ,可以采用应变能函数 W [ 3 ] 表示 . W 既可以表示为应变不变量的多项式 : W = W ( I1 , I2 , I3 ) , 也可 以直接用主伸长比表示 : W = W (λ1 ,λ2 ,λ3 ) . 有了应变能函数 就可 以推 导 得出 材料 的 应力 2应 变关 系 ,即常用的第二 Piola2 Kirc h hoff 应力张量σij 和 Gree n应变张量εij 间的关系[ 4 ] :Rivli n 公式由于形式简单 、易于实现而被广泛应用 ,其应变能函数为 :W = C1 ( I1 - 3) + C2 ( I2 - 3)( 2)式 ( 2) 中的 C1 , C2 可认为是橡胶类止水材料的超弹性常数 , 方程 ( 2) 可得出的应力 2应变关系如下 :单向拉伸和单轴压缩 :σ= 2 (λ- λ- 2 ) C1 + 2 ( 1 - λ- 3 ) C2纯剪 :σ= 2 (λ- λ- 1 ) C1 + 2 (λ- λ- 1 ) C2( 3)( 4)Moo ney2Rivli n 公式虽然形式简单 , 应用方便 ,但是由于其应变能函数没有取足够多的展开项 , 因 此在描述材料性能时常会出现令 人 难以 接受 的 误 差. 此前 , 作者利用 Moo ney2Rivli n 公式对水封橡胶 材料的非线性进 行大 量 实验 与计 算研 究[ 5 27 ] , 发 现在计算中若采用 Moo ney2Rivli n 公式的拟 合 参数 , 计算结果误差较大 , 计算精度不能保证 . 为此 , 作者 首次尝试在水封的非线性有限元 计 算研 究中 采 用 O gde n 公式描述橡胶类止水材料的力学性能 , 其应 变能函数为 :5W 5 I1 5W 5 I2 5 W 5 I3μn6 α (λααασij = 5 I( 1)λ+- 3λ)( )51 5ε + 5 Iij 2 5ε + 5 IW =+1 n2 n3 nij 3 5εijnn式中 I1 、I2 、I3 为 Gree n 应变不变量 , 分别为 : I1 =λ2式中αn 可取任何值 ,μn 为待定参数 . 一般情况下 , 为了表示三种常见应变模式 ( 拉伸 , 压缩和纯剪) , 采用 三项表达式 , 即 O gde n 公式 , 方程 ( 5) 中 n = 1 , 2 , 3 , 对应的应力应变关系如下 :1+λ2 +λ2 , I2 =λ2λ2 +λ2λ2 +λ2λ22 2 21 2 2 3 1 3 , I3 =λ1λ2λ3 ,λ1 、λ2 、2 3λ3 为主伸长比 , 且λi = 1 +εi ( i = 1 , 2 , 3) .在闸门水封的非线性有 限 元计 算中 , Moo ney23 国家自然科学基金资助项目 ( 资助.2010201220 收到第 1 稿 ,2010205231 收到修改稿.E2mail