sobol灵敏度分析matlab_闭口式轧机牌坊的有限元分析，轧机技术（61）

阅读本文前，您可先点击上面的蓝色字体“云轧钢”，“关注”后您就可以第一时间收到最优质的文章了。 添加小编微信(13585097569)把您想学的知识留言小编，小编尽量收集！

摘要:运用ANSYS12 软件对厚连铸板坯切割制造的闭口式轧机牌坊进行三维有限元分析，计算牌坊的应力和变形，计算结果为闭口式轧机牌坊的设计和制造提供了重要依据。

0、引言

轧钢机牌坊是轧钢机的重要零件，用于安装轧辊辊系和轧辊调整装置，承受全部轧制力。因此，轧钢机牌坊作为永久使用的不更换零件来进行设计， 并且牌坊重量大、制造复杂，为防止牌坊出现破坏， 一般给予很大的安全系数。为保证牌坊具有足够的强度和刚度，一般采用增加断面尺寸裕度、改善材料等途径来提高安全系数。但是，牌坊在使用中产生裂纹甚至发生断裂或者牌坊刚度不够，弹跳大的情况仍有发生。进一步研究牌坊的合理设计方法，保证牌坊使用的可靠性和足够的刚度仍是目前值得重视的问题。

有限单元法是一种用于连续场分析的数值模拟技术，其与弹性力学相结合进行分析研究，能了解牌坊的应力分布规律和危险断面。随着电子计算机技术的飞速发展，各种计算软件( 如: MATLAB、AN- SYS、SAP 等) 的出现，使有限元法的计算精度提高， 使用更方便，成为应力场分析的利器。

1、有限元法的基本理论

有限元法的基本思想是将连续的实际问题的求解域离散化，得到有限个单元，单元彼此之间仅靠节点相连。在单元内假设近似解的模式，通过适当的方法，建立单元内部点的待求量与单元节点量之间的关系。由于单元形状简单，易于由能量关系或平衡关系建立节点量之间的方程式，然后，将各个单元方程集合成总体线性方程组，引入边界条件后求解该线性方程组，即可得到所有节点量，进一步计算导出量后问题就得到了解决。

选择节点位移作为总体线性方程组的未知数， 位移函数是单元上点的位移对点的坐标的函数，用 单元内部点的坐标的多项式来表示，近似地表示了 单元内真实位移的分布。计算单个单元刚度矩阵， 将全部单元刚度矩阵相互结合起来构成整个结构的 总体刚度矩阵，把非节点载荷等效移置到节点上形 成结构总体载荷列阵; 引入约束条件和已知力求解线性方程组得到节点位移; 根据节点上位移计算单元的应力和应变。

2、ANSYS 软件简介

ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国 ANSYS 开发，它能与多数 CAD 软件接口，实现数据的共享和交换，如 Pro / Engineer， NASTＲAN， Alogor， I － DEAS，AutoCAD 等，是现代产品设计中的高级CAE 工具之一。ANSYS 系统含有参数化设计语言( APDL) ，它具有参数、数学函数、宏、判断分支及循环等高级语言要素，是一个理想的程序流程控制语言，很适合进行有限元计算和优化分析。ANSYS 有限元软件包是一个多用途的有限元法计算机设计程序，可以用来求解结构、流体、电力、电磁场及碰撞等问题。软件主要包括三个部分: 前处理模块，分析计算模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具，用户可以方便地构造有限元模型; 分析计算模块包括结构分析( 可进行线性分析、非线性分析和高度非线性分析) 、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力; 后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示( 可看到结构内部) 等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。

3、建立有限元模型

某轧钢厂在技术改造中选择闭口式牌坊机架作为连轧机架，结构如图 1 所示。

为简化制造工艺，选用厚连铸板坯切割下料，板坯厚度为0.245m，只对窗口部分和固定面A面进行机械精加工; 在加工过程中为简化加工工艺，只在少数部位采取圆弧过渡，不像整体铸造牌坊在多数部位均采取圆弧过渡，为避免应力集中的发生，有必要对该牌坊进行有限元分析，检查是否存在应力集中影响牌坊的刚度和安全。

在实际轧制过程中，牌坊受力情况很复杂，除轧制力外、还有轧件咬入时的冲击力、传动系统引起的轴向力、摩擦力其它力的影响，但以轧制力最大，其它力远小于轧制力，因此，忽略其它力的影响，认为牌坊只受轧制力。

牌坊受力情况是静力受力结构，该轧机允许的最大轧制力为 2 500 kN，单片牌坊的最大轧制力为1 250 kN，以均布载荷的形式作用于牌坊的支承面上，因为牌坊上下支承面面积不同，换算后的均布载荷分别为上支承面 15 MPa，下支承面 12 MPa。将牌坊置于 X，Y，Z 坐标系中，在轧机固定面处设置 X，Y，Z 方向约束，其等效力学模型图 2

在轧制力的作用下，牌坊处于三向应力状态，因 此采用三维块单元对其进行离散，厚板坯的材料是16Mn，属于结构钢，是线性各向同性材料，弹性模量为 2．1 GPa，泊松比为 0．3。SOLID45 单元一般用于构造三维各向同性材料固体结构; 单元通过8个节点来定义，每个节点有3个沿着 x，y，z 方向平移的自由度; 单元具有塑性、蠕变、膨胀、应力强化、大变形和大应变能力，但是它的缺点是计算结果与实际测量得到的数据相差较大。SOLID95 是 SOLID45( 3维 8 节点) 高阶单元形式，具有协调的位移函数并且能很好的模拟边界曲线，因此单元能够容许不规 则形状，并且不会降低精确性，特别适合边界为曲线的模型; 同时，其偏移形状的兼容性好。SOLID95 有 20 个节点定义，每个节点有 3 个自由度( X，Y，Z 方向) ，此单元在空间的方位任意，本单元具有塑性、蠕变、辐射膨胀、应力刚度、大变形以及大应变的能力，并且提供不同的输出项。SOLID95 单元的计算结果与实际测试得到的数据比较接近。但缺点是由于 SOLID95 是 20 节点，会耗费更多的时间和计算机资源; 在计算机能力能满足要求时应选择 SOL-ID95 单元。因此这里我们采用 SOLID95 单元对牌坊进行网格划分，为达到相应的计算精度，同时又要减少计算量，设定单元边长不超过 0．06 m，并开启 ANSYS 的局部细化网格控制功能对局部进行网格细化，采用自由网格方式模型进行划分网格，生成有限元模型，共生成节点128 446 个，划分单元 102 863 个，能满足分析精度要求。有限元模型图如图 3 所示。

4、有限元计算结果分析

ANSYS的计算过程通俗的讲就是解方程，通过各类求解器，求解由有限元方法建立的联立方程组。

进入求解器后，选择静力分析类型和小变形静力稳态分析进行求解。根据求解结果生成牌坊应力云图如图 4 所示，牌坊变形图如图 5 所示。最大应力出现在牌坊的上支承面处，最大应力值为 39．2 MPa，小于许用应力［σ］= 50 － 70 MPa［1］满足使用要求。牌坊的纵向最大变形出现在上支承面处，为0．166 mm，小于［f］= 0．5 mm［1］满足使用要求，牌坊的纵向刚度为:1 250 kN/0．166 mm = 7 530 kN/ mm

牌坊横向最大变形出现在立柱中间段为 0．11 mm，在最大轧制力的作用下，两立柱将向中间靠近 0．11× 2 = 0．22 mm。因此，牌坊和轴承座在加工装配时要预留一定的间隙，防止在工作时发生卡死现象。

5、结论

根据有限元计算结果分析，可以得出下列结论:

1)牌坊的纵向刚度为7530 kN/ mm，单片牌坊在1250 kN 的最大轧制力的作用下，最大变形为0．166 mm，能完全达到轧制过程中的质量要求。

2)纵观整个牌坊的应力分布均匀，没有应力集中现象，因此，采用厚连铸板坯切割后精加工制造闭口式牌坊的制造工艺是可行的。

3)在最大轧制力的作用下，牌坊窗口的横向尺寸减少022 mm，为使变形后的牌坊不会夹死轴承座，牌坊耐磨衬板与轴承座间的间隙不得小于 0．22 mm。

版权声明：本文由云轧钢整理发布，版权归原作者所有。转载请注明来源；文章内容如有偏颇，敬请各位指正；如标错来源或侵权，请跟我们联系。

现已开通云轧钢冷轧、热轧、炼钢、设备等技术交流群，进群添加小编微信(13585097569)即可入群。

欢迎关注“云轧钢”公众号，我们提供钢铁行业各类资讯、各类技术、现场案例分析、技术交流、并同步开通云轧钢技术交流群产品宣传推广服务、钢厂求职招聘咨询。