基于ansysworkbench的密炼机转子端面静力分析.docxVIP

基于ansysworkbench的密炼机转子端面静力分析 密度机是橡胶加工过程中必要的设备。密度机的混合和分散效果，确保了橡胶和各种配合材料的良好混合和分散效果。 本工作通过ANSYS Workbench软件中的结构静力（Static Structure）分析模块，对密炼机转子端面无介质密封（混炼物料颗粒不进入端面）和介质密封（混炼物料颗粒进入端面）2种情况进行有限元分析，研究介质（混合物料）泄漏对动静密封环使用寿命的影响。 1 振动环模型的构建 1.1 安装压紧导柱和压紧封环 密炼机转子端面密封结构如图1所示，处理前动静密封环结构如图2所示。 图1所示的密封结构采用单端密封，静密封环固定端均布3个直径4 mm的小孔，用来安装压紧导柱。在有限元模拟中，孔的存在会产生应力集中，且使得网格划分复杂，由于3个孔仅用来安装压紧导柱，对接触瞬态分析影响不大，因此为了提高网格划分质量，减少计算量，动静密封环简化模型去掉了孔，如图3所示。 1.2 材料参数的配置 动静密封环皆为45 1.3 网络划分中的输入和中断 1.3.1 接触的区域范围 弹性力学中，接触问题是指2个物体在受压状态下所发生的力学行为，包括局部应力和应变等；在实际工程中2个物体相互接触并产生相对运动时接触表面即产生摩擦；在有限元中，以节点来定义接触，理论上2个节点在某一个区域或线上恰好重合即被认为发生了接触。但是在计算中不能精确描述节点的接触，需要设置接触区域。实际计算过程中，当某一个时间增量或载荷步中接触物体移动超出了接触距离，就会被认为发生了穿透 接触模型可描述为求区域内位移场，使系统势能在接触边界条件约束下达到最小，即 式中，K,U和F分别为有限元位移法中的刚度矩阵、未知结点位移向量和结点荷载向量。 实际计算中，由于相互接触的物体材料属性、力学性能等不同，有限元将接触划分为两类：刚性体-柔性体接触和柔性体-柔性体接触。在密炼机密封装置的动静密封环中，动密封环在密炼室内部，静密封环在密炼室端面，更换静密封环较为容易，动静密封环材料均为45 1.3.2 接触模型的建立 穿透在实际接触物体之间不会发生，因此在有限元分析软件中采用程序强制接触协调来控制穿透，以期模拟真实接触状态。罚函数法接触模型如图4所示，当接触协调性不被强制时就会发生穿透。Workbench Mechanical模块中提供了罚函数法 (1）罚函数法。罚函数法要求有限的穿透量在交界面处产生接触压力，保持系统平衡需要此接触压力，即接触压力必须大于零。然而，实际的接触物体相互不穿透，因此为了得到最高精度，应使发生在接触界面处的穿透量最小，这意味着理想的接触刚度应非常大，但是刚度值太大会引起收敛困难。如果接触刚度太大，微小的穿透即会产生很大的接触压力，在下一次迭代中可能会将接触面推开，采用过大的接触刚度通常会导致收敛振荡且发散。 (2）拉格朗日乘子法。拉格朗日乘子法接触模型如图5所示，其实质是增加一个附加自由度（接触压力）来满足不可穿透条件。 (3）增广拉格朗日乘子法。将以上两种方法结合起来，就是增广拉格朗日乘子法，其接触模型如图6所示。接触协调性基于惩罚刚度确定，一旦达到平衡，检查穿透容差，如果有必要，接触压力增大，迭代继续。 1.3.3 接触因数和介质因数 在密封端面摩擦磨损过程中，理论上金属材料的摩擦因数与表面接触面积及相对运动速度无关，但是在密炼机密封装置中摩擦因数随着介质环境发生变化，且介质环境影响较大，甚至改变端面接触应力和应变的大小。端面密封副的摩擦因数包括接触摩擦因数和介质摩擦因数。摩擦因数（f）的计算公式为 式中，F 在平衡情况下，W为 式中，W 混合摩擦状态下密封副的f为 式中，x 微凸体F 式中，R 综上所述，在此次模拟中，静密封环表面设置为接触面，动密封环表面设置为目标表面，端面接触摩擦因数为0.2，端面介质摩擦因数为0.3，默认增广拉格朗日乘子法，使结果更容易收敛。模型采用正六面体划分网格，网格总数为44 800，网格划分模型如图7所示。 1.4 边境条件的定义 动密封环随着转子旋转，静密封环固定且受横向压力，两表面接触起到密封作用，因此在动密封环上施加转动惯性负荷，转动速度为120 r·s 2 结果与讨论 2.1 浮动密封环接触分析 在动静密封环端面直接接触和介质接触2种状态下，对动静密封环接触端面进行应力分析，结果如图8所示。 从图8可以看出，在端面直接接触和介质接触时，在均匀密封压力作用下，动静密封环的受力均沿圆周方向均匀分布，在接触端面压力最大，摩擦力也最大，磨损较为剧烈。 对比图8(a）与（b）可以看出，动静密封环端面直接接触时，其最大应力为1.373 8 MPa，端面介质接触时，其最大应力为2.108 2 MPa，相较于直接接触增大了53%。因此，密炼过程中，混合物料泄漏会导