密炼机二维流场流固耦合分析.docVIP

密炼机二维流场的流固耦合分析 边慧光 张俊岭 汪传生 (青岛科技大学机电工程学院，山东 青岛 266061) 摘 要：本文利用ADINA软件对密炼机同步二棱和同步四棱转子的二维流场进行了模拟研究，实现真正意义上的流固耦合的模拟计算，并在后处理中对压力场、速度场、剪切应力场进行了理论分析，可以直观地观察到在不同的瞬时压力场、速度场、剪切应力场的变化情况，为转子优化设计提供数据参考。 关键词：同步转子 有限元 网格划分 流场 流固耦合 流体、固体动力耦合广泛存在于水利、船舶、海洋以及航空航天等多个工程领域中，是一个十分重要、多学科交叉的研究课题。近一个世纪以来，国内外学者对流体与固体的动力耦合理论和方法开展了广泛的研究，取得了许多的成果，近阶段随着数学方法的逐渐成熟，越来越多地应用于实际生产和大型工程中。 ADINA[1]可以到流固耦合(Real-Time Two-Way Fluid-Structure Interaction)功能,可以在流固各使用不同的,且在界用不吻合,程依然可快速流固耦合除此之外,在流固, ADINA可以重建,以避免因而造成扭曲,ADINA中的FSI是流固耦合,而非一般CFD的其重建也非一般的ALE方法,而是真正的重建 图1 两棱同步转子示意图 图2 四棱同步转子示意图 转子和密炼室的几何参数列表如表1所示： 表1转子和密炼室的几何参数列表 单位 mm 参数 名称 基圆直径 剪切棱半径 非剪切棱半径 中心距 顶隙 棱顶宽度 容积 （L） 同步四棱 70 70 100 120 2 6 1.42 同步二棱 70 70 100 120 2 6 1.65 1.2数学模型 本文主要研究在密闭环境下完全充满胶料的流场模型，根据假设，流场为等温度场，即各点温度相同，故不考虑能量转换。首先给出广义的等温和不可压缩流体的偏微分控制方程。 （1） （2） （3） （4） 式中：D/Dt——物质导数 V——速度 P——压力 ——密度 G——重力 ——额外应力张量 为了既描述在高剪切速率下的假塑性流体的流变性质，又描述在低剪切速率下牛顿流体的流变性质，本文采用Carreau[2]模型： （5） 式中：—初始粘度，也就是零剪切速率时的粘度 —粘弹性的特征时间 —无穷剪切粘度 n—幂律指数 —混炼胶密度 其中：，=10，，n =0.23，＝1066Kg/m3 这些参数是110℃的条件下，测量混炼胶的各种参数而得到的。由n =0.23得出，可以得知本文假设胶料的流体为塑性流体。 其中转子的材料采用45号钢，其物性参数为：泊松比v=0.29,弹性模量E=2.07e11N/m2,密度=7800Kg/m3。2.有限元模型和边界条件 2.1 有限元模型 2.1.1 坐标系的建立 如图3所示，选择笛卡儿坐标系，在YZ平面内进行建模、模拟。 图3坐标系示意图 2.1.2 网格划分 对安装有同步四棱转子和同步二棱转子的密炼机流场进行模拟分析，取整个流场的中间截面。考虑到模拟过程为混炼后期，可以假设胶料趋于一种均匀的流态，转子受力情况将非常小，转子中间的水道结构将不会影响到转子的整体强度，故采取省略转子中间水道的做法，两类转子流场的有限元模型如图4所示。  （a）同步四棱流场 （b） 同步二棱流场 图4两转子流场网格模型 2.2边界条件 2.2.1压力边界条件 密炼机混炼过程中密炼室是完全密闭的，而且有限元模型假设是完全充满的，这样胶料流体就没有进出口。在这样一个边界条件下，由于密炼室内的转子高速转动，从物理模型上来理解，就会产生一个巨大的压力场，导致模型不稳定，所以必须对压力低的区域进行压力零点设置，首先设定约束边界条件，设置压力零点也就是约束压力自由度，这样整个压力场有了基准进行对比，相应的求解模型就会变得稳定。 2.2.2速度边界条件 转子流场模型不需要施加载荷，只要对转子的固体模型加载一个绕中心旋转的角速度，两转子进行相对回转即可。 （6） 转子旋转的速度为60r/min, 故=6.28 r/ s。 2.2.3边界条件加载方式 流固耦合很注重模型的稳定性，所以载荷不能在一个瞬时完全加载到模型上，需要从零开始逐步加载，而且计算的每个时间步不能太大，这就需要对步序时间进行严格控制。瞬时计算的时候，步序时间为真实的物理时间，每个步序时间定为0.5S。载荷加载曲线如图5所示： 图5速