大型隔膜泵橡胶隔膜的有限元分析_凌学勤_图文

1、大型隔膜泵橡胶隔膜的有限元分析凌学勤1张伟21东北大学,沈阳,1100042中国有色(沈阳泵业有限公司,沈阳,110141摘要:针对隔膜腔的隔膜工作过程,建立简化的二维轴对称结构隔膜有限元模型,对其进行隐式动力学分析,同时建立三维对称结构的隔膜腔及隔膜整体有限元模型,对其进行流固耦合分析,最后针对两种不同的方法和结果展开对比,以此期望得出模拟隔膜运动中应力和位移变化的合理数值模拟方法,以指导隔膜设计。关键词:隔膜泵;隔膜;流固耦合;数值模拟;隐式动力学中图分类号:TB301文章编号:1004132X (201115181904Finite Element Analysis of Large D

2、iaphragm Pump Rubber DiaphragmLing Xueqin Zhang Wei1Northeastern University ,Shenyang ,1100042NFC Shenyang Pump Industry Co,Ltd,Shenyang ,110141Abstract :In allusion to the work process for a diaphragm of diaphragm chamber ,a simplified two di-mensional axisymmetric finite element model of membrane

3、structure was established ,and laid bare the implicit dynamics analysis ,while a three dimensional symmetrical structure whole finite element model of the dia-phragm chamber and diaphragm was established and laid bare the analysis for fluid solid couplingAt last ,for two different methods and result

4、s began comparison in order to elicit a reasonable numerical simulation method which will be used to simulate the changes on stress and displacement of diaphragm movement ,so as to guide reasonable diaphragm designKey words :diaphragm pump ;diaphragm ;fluid and solid interaction ;numerical simulatio

5、n ;implicit dy-namics收稿日期:201010110引言大型隔膜泵橡胶隔膜把液压油和输送浆体隔离开来,使输送的浆体无法外漏,同时避免了浆体中高腐蚀性、高磨砺性的颗粒性介质对泵体造成的损坏。橡胶隔膜使用寿命的长短直接关系到隔膜泵的性能和工作效率。而隔膜是在封闭的隔膜室内往复工作的,这就给实验测试隔膜工作时的受力及变形状况带来了很大的麻烦。因此有必要采用有限元分析的方法对隔膜的工作过程进行数值模拟,研究隔膜工作时的受力及变形情况以指导隔膜的设计。目前,不考虑流体对隔膜作用的单纯隔膜结构件的有限元分析已有不少研究成果,但对于考虑隔膜真实工作时候两侧流体影响的数值模拟研究却鲜有报道。

6、文献1将隔膜简化为壳单元,同时施加位移载荷作为隔膜运动的边界条件。文献2将隔膜简化为圆盘形状来进行分析研究,这和隔膜实际形状出入太大。文献3对隔膜头部密封结构的密封性能进行了静态有限元分析,得出了一些结论,然而实际隔膜工作是动态的过程。文献4对真实几何模型的隔膜做了二维轴对称结构的有限元分析。本文针对隔膜腔的隔膜工作过程,建立简化的二维轴对称结构隔膜有限元模型,对其进行隐式动力学分析,同时建立三维对称结构的隔膜腔及隔膜整体有限元模型,结合隔膜真实工况对其进行流固耦合瞬态分析。1隔膜橡胶材料的基础实验及数据拟合橡胶材料不同于金属材料仅需要几个参数描述其材料特性,橡胶的行为复杂,材料本构关系是非线

7、性的,简单依赖单向拉伸性能实验并不能完全描述材料包括压缩及剪切在内的所有力学行为,为使隔膜的有限元分析和结构模拟接近工程实际,必须对橡胶材料进行包括拉伸、压缩、剪切及体积实验等在内的全部基础实验。图1、图2所示是本文研究工作中所用到的一组橡胶材料数据,该实验在美国AXEL 实验室完成,实验橡胶试件的配方采用中国有色(沈阳泵业有限公司生产隔膜产品的配方。橡胶力学性能的描述方法主要分为两类5:一类是认为橡胶为连续介质的现象学描述;另一 图1拉、 剪实验数据图2体积压缩实验数据类是基于热力学统计的方法。对于各向同性材料,应变能可加法分解成应变偏量能和体积应变能两部分,形式如下5:U =f (I 13

8、,I 23+g (J 1(1令g =Ni =11D i(J 12i ,并且进行泰勒展开,可得U =N i+j =1C ij(I13i (I 23j+Ni =11D i(J 12i(2式中,I 1、I 2分别为第一、第二偏差应变不变量;J 为弹性体积比;C ij 为与温度有关的材料参数;i 、j 表示第i 阶和第j 阶材料参数。这种形式为多项式表示的应变能,参数N 为我们选择的多项式阶数。D i 的值决定材料是否可压:如所有D i 都为0,说明材料是完全不可压的。对于多项式模型,无论N 值如何,初始的剪切模量0和初始的体积模量k 0,都决定于多项式一阶(N =1系数:0=2(C 10+C 01k

9、 0=2/D 1(3对于完全多项式,如果N =1则只有线性部分的应变能量保留下来,即Mooney Rivlin 形式:U =C 10(I 13+C 01(I 23+1D 1(J 12(4这种形式是最简单的超弹模型,未知精确参数的材料通常使用这种形式。这种函数不能够表示应力应变曲线大应变部分的“陡升”行为,但是在小应变和中等应变时是可以较好地模拟材料特性的6-7。利用给出的这组实验数据,材料的参数可由最小二乘法确定。对于n 组应力应变的实验数据,相对误差E 取最小值:E =ni =1(1Tthi /T test i2(5式中,T testi 为实验数据中的应力值;T thi 为按照本构关系与伸长

10、率对应的应力表达式。通过线性最小二乘法对以上线性方程组进行求解,可以得到确定的弹性系数C ij 。拟合得到的材料参数见表1。表1多项式拟合得到的材料参数Mooney Rivlin 弹性系数C ij2参数C 01=04469;C 10=0 09973参数C 10=02155;C 01=04757;C 11=038125参数C 10=00350;C 01=07276;C 20=5789;C 11=1172;C 02=50429参数C 10=03389;C 01=1181;C 20=8788;C 11=1763;C 02=9088;C 30=8967;C 21=2264;C 12=1532;C 03

11、=36462橡胶隔膜隐式动力学分析21橡胶隔膜的几何结构隔膜泵橡胶隔膜的几何结构8如图3所示。从图3可以看出橡胶隔膜为轴对称结构,故在分析时可以只对其截面进行分析。图3橡胶隔膜的几何结构22橡胶隔膜隐式动力学有限元模型的建立对模型进行适当简化,隔膜头部完全约束,并在隔膜两侧施加不同压力,使其在两侧微小压差作用下产生运动。施加边界条件的橡胶隔膜平面对称有限元模型如图4所示。图4橡胶隔膜的位移及应力加载方式23计算结果及分析通过图5、图6可以看出,隔膜在加载后,所受的最大应力点在隔膜发生翻转的地方,这为隔膜设计提供了设计依据,在设计隔膜时,可以优化隔膜翻转处的尺寸,从而避免应力过大而造成的隔膜破损

12、。 图5橡胶隔膜截面的应力分布图 图6橡胶隔膜截面的位移分布图3橡胶隔膜流固耦合分析 上述基于隐式动力学分析的方法对隔膜变形及受力过程的分析与实际隔膜的工况有一定的偏差,具体偏差体现在以下几个方面:(1实际隔膜工作时两侧受到的液体压力是不断变化的,而上述分析将液体的作用力简化为恒定不变的压力直接施加在隔膜内外表面上。(2实际隔膜工作时是三维结构,而上述分析将模型简化为二维轴对称结构,与实际情况有一定的偏差。 (3实际隔膜工作时的每一个周期内均有料浆吸入和排出过程,同时进料口和出料口均存在一定的喂料压力和出料压力,而上述分析无法体现出这种实际工况对于隔膜运动和变形以及受力的影响。 鉴于上述方法存

13、在的不足和局限以及隔膜实际工作时所处的工作环境及工况,对隔膜腔整体模型进行简化处理,考虑进出料阀门的影响,以及进料压力和出料压力的作用,建立三维二分之一对称结构流固耦合模型,进行三维流固耦合分析。31隔膜腔隔膜工作原理的几何结构隔膜泵用的隔膜腔主要由隔膜、活塞、进出料阀门及金属导杆组成,如图7所示。32隔膜腔隔膜流固耦合分析有限元模型的建立建立实际三维二分之一对称结构隔膜腔整体图7隔膜泵隔膜腔示意结构流固耦合分析模型,流体部分模型和结构部分模型分别建立;不考虑金属导杆的作用;进出料处流体单元施加类似阀门功效的GAP 特殊边界条件,通过GAP 的开闭来控制流体的吸入和排除过程;在进出料处分别施加

14、进料、出料压力;将活塞位移引起的液压油推进过程用活塞缸入口的速度代替;流体及结构两侧对应耦合界面处分别定义流固耦合特殊边界条件。施加边界条件的隔膜腔隔膜流固耦合分析有限元模型如图8所示。(a 流体部分边界条件(b 结构部分边界条件图8隔膜腔整体有限元模型边界条件33计算结果及分析对整体模型进行流固耦合分析,得出一个周期不同时段的结果,如图9所示。(a 二分之一周期时刻(b 四分之三周期时刻(c 最终时刻图9流场流速矢量图最终时刻隔膜的受力云图如图10所示。 图10最终时刻隔膜受力云图对一个冲次不同时刻的分析结果进行提取、整理,分别得出隔膜位移、应力、吸入容积与时间函数的关系曲线,如图11 图1

15、3所示。 图11隔膜位移随时间变化曲线 图12隔膜应力随时间变化曲线 图13隔膜吸入体积随时间变化曲线由图11 图13可以看出,随着时间的推进,隔膜位移、应力及吸入体积均逐渐增大,同时当时间为0.65s 时隔膜应力迅速增大,到0.95s 时隔膜应力高达19.6MPa ,应力过大容易导致隔膜过早撕裂,因此隔膜运动至0.65s (对应位移为197.08mm 的时刻为一个冲次结束的时刻。按照该过程设计橡胶隔膜后进行现场工业运行,记录结果表明隔膜使用寿命从6000h 延长到了8000h ,从而说明了流固耦合分析过程的正确性。4结论(1三维流固耦合分析方法反映了隔膜实际工作时两侧液体对隔膜的压力随时间变

16、化的真实特性。而不考虑隔膜两侧流体作用力的二维隐式动力学分析方法将两侧液体的压力简化为恒定不变的压力,这导致隔膜的受力情况与实际情况有较大的偏差。(2流固耦合分析方法是三维模型,比二维轴对称结构模型更接近实际情况。(3三维流固耦合分析方法建立了整体隔膜腔结构的模型,考虑了进出料阀门开启关闭以及进料压力和出料压力对于隔膜腔内流场变化的影响,更符合实际情况。而二维轴对称结构分析方法无法体现出这种实际工况对于隔膜运动和变形以及受力的影响。(4依照三维流固耦合分析方法修改隔膜一个冲次作用时间后的现场工业测试结果表明,该方法有效地延长了隔膜的使用寿命(8000h 。因此,隔膜腔隔膜三维流固耦合分析方法能

17、充分反映隔膜运动过程中的变形、位移、受力情况,对于隔膜的形状尺寸设计以及受隔膜腔内流体的流场分布影响的隔膜腔腔型设计有一定的指导意义。参考文献:1李振河隔膜泵及其管道的技术研究D 沈阳:东北大学,20082陶自强往复式隔膜泵隔膜设计及行程控制系统的研究D 沈阳:东北大学,20083Lee B ,Rivin EFinite Element Analysis of Load de-flection and Characteristics of Compressed Rubber Components for Vibration Control Devices J Jounrnal of Mechanical Design ,1996,11(8:328-3354Charlton D ,Yang JA Review of Methods to Charac-terize Rubber Elastic Behavior for Use in Finite Ele-ment Analysis J Rubber Chemical and Technology ,1994,67(3:481-5035ABAQUSTheory Manual M Botha :ABAQUS Inc ,20036Yeoh OSome Forms of the S