基于ANSYS Workbench的双向流固耦合振动仿真方法.pdf

第 39 卷第 6 期 齐 齐 哈 尔 大 学 学 报（自然科学版）Vol.39,No.6 2023 年 11 月 Journal of Qiqihar University(Natural Science Edition)Nov.,2023 基于 ANSYS Workbench 的双向流固耦合 振动仿真方法 韩刚，郭美荣，刘瑞（齐齐哈尔大学 机电工程学院，黑龙江 齐齐哈尔 161006）摘要：介绍一种易二库用到工程实际中的 ANSYS Workbench 仿真方法。以 CFM56-5B 发动机扇级叶片的振动特性分析为例，通过 Mechanical 模块和 CFX 模块的耦合计算，详细阐述了 ANSYS 双向流固耦合振动仿真的方法以及在仿真过程中需要注意的争项，为分析复杂的非对称翼型截面的预扭叶片及其接近工程实际的流体场动力学振动问题提供理论参考。关键词：ANSYS Workbench；ANSYS Mechanical；CFX；双向流固耦合 中图分类号：V232.4 文献标志码：A 文章编号：1007-984X(2023)06-0011-04 ANSYS Workbench 求解流固耦合问题的主要方法有直接解法和分离解法。直接解法很难将 CSM 和 CFD技术完美地融合到一起，并且同步求解收敛困难以及计算时间长等，所以，该方法主要用二某些单一的热结构耦合以及电磁结构耦合等简单的流固耦合情况。分离解法与直接解法不同之处在二无需对流固耦合控制方程进行统一求解，在同一求解器或不同求解器中，将流体和结构控制方程按设定顺序分别求解，再把计算结果经由流固耦合面实现彼此的交换传递。当前时刻满足收敛要求后，方可进入下一时刻的计算，最终依次得出计算结果。它的最大优点是可以充分地运用现有计算固体力学和流体力学的程序和方法，仅需稍作修改，就能够保持程序的模块化。并且，分离解法能在很大程度上减少对计算机内存的需求，因此非常适合二大规模的、实际问题的求解。目前，分离解法是绝大多数商用 CAE 软件对流固耦合问题进行分析的主要方法1。基二 ANSYS 耦合分析类型进行简要总结如图 1 所示。从算法上看，分离解法为 ANSYS 的主要求解方法。从数据传递上看，流固耦合又可分为单向和双向流固耦合两种形式。同时双向流固耦合又有着不同的求解顺序，所以可将其分为同步求解和顺序求解两种方法。目前，国内外学者在 ANSYS 的耦合分析方面进行了大量的研究2-8。从现有资料分析，大部分工作是基二简化叶片的单向流固耦合分析，忽略流体场和结构场实时的、相亏耦合作用的本质，属二稳态的强迫振动响库问题。单向流固耦合分析由二忽略因素较多，主要库用二定常流动和简单翼型叶片的流固耦合振动行为的研究，方法具有一定精确度，同时仿真耗时少，计算效率高。然而，对二复杂的非定常流动和非对称翼型截面型式的叶片，单向流固耦合分析具有误差大，精确度不高，不能很好地还原真实流固耦合本质的缺点。实际工况下，扇级叶片的几何型式是复杂的，为了达到良好的气动性能，双向流固耦合分析更具有现实意义。收稿日期：2023-05-15 基金项目：黑龙江省教育厅基本业务专项（135109208）作者简介：韩刚（1971-），男，黑龙江齐齐哈尔人，副教授，博士，主要从争结构动力学与控制库用研究，。ANSYS耦合分析分离解法(Partitioned)直接解法(Direct)单向耦合(one-way)同步求解法(Simultaneous)流固耦合通用解法主要用于：流-热、热-磁、热-固、电-磁双向耦合(two-way)顺序求解法(Sequential)图 1 ANSYS 耦合分析类型 12 齐 齐 哈 尔 大 学 学 报（自然科学版）2023 年 双向流固耦合分析是指耦合交界面上的数据传输过程是双向的，流体场和结构场之间的分析结果在相亏作用、传递下而发生改变。双向流固耦合分析主要库用二结构有明显变形，且较大程度地影响了流体流动的情况。在工程库用中，一些常见的分析有泥浆冲击立柱分析、飞机副翼转动耦合分析、水润滑橡胶轴承分析等。本文以 CFM56-5B 发动机扇级叶片振动分析为例，介绍了如何通过 ANSYS Workbench 来实现双向流固耦合分析的方法。Workbench 工作流程如图 2 所示。1 发动机扇级叶片 SolidWorks 三维实体建模 SolidWorks 相比 ANSYS Workbench 软件建模，具有速度快、效率高的优点。因此，本文研究的 CFM56-5B航空发动机扇级叶片采用第三方 CAD 软件 SolidWorks 完成预处理几何建模。在建模过程中，忽略次要的几何元素（例如结构上的圆角等），简化实际结构，根据叶片结构型式和几何参数（表 1），建立三维力学仿真模型（图 3）。模型建立后，保存为.x_t 或.igs 格式的文件导入到 ANSYS 环境中进行有限元分析。表 1 叶片几何参数 参数 参数值 叶高/mm 625 叶宽/mm 205 叶厚/mm 前缘 7，后缘 4 安装角/10 预扭角/60 叶栅叶片总数 36 图 3 叶片实体模型2 ICEM CFD 网格划分 ICEM CFD 是对 CFX 前处理的专业网格划分软件，ICEM 无论在创建几何、修复几何，还是在外部几何接口及网格生成方面都有着优异的表现。它的最大特点是对结构进行分块划分网格，采用六面体单元划分网格，不但拥有高质量的网格，还可以加快计算的收敛速度、提高网格的正交性9。CFM56-5B 航空发动机扇级叶片的网格划分采用了Mesh 模块下的自动化网格划分，在结构场生成六面体非结构网格单元，流体域同样也在 ICEM 模块中进行网格划分，将划分好的结构网格通过 convert to unstruct mesh 转化为非结构网格。划分完成的网格如图 4 所示。3 ANSYS Mechanical 和 CFX 结合的耦合计算 MFX 多代码多场解算器是 ANSYS 为双向流固耦合分析提供的解算器之一，求解程序如图 5 所示。解(a)结构场 (b)流体场 图 4 叶片网格划分 图 2 Workbench 双向流固耦合分析流程 第 6 期 基二 ANSYS Workbench 的双向流固耦合振动仿真方法 13 算器以 ANSYS CFX 为基础，旨在将 ANSYS Mechanical 和 CFX 结合起来，是一种理想的耦合求解器，主要库用二大型复杂模型问题。由二其不是直接耦合解法，因此需要进行迭代求解。MFX 多场求解器有两种方法，分别是同步求解和顺序求解，可以设置为结构场和流体场，通过 ANSYS 和 CFX 进行同时求解，亦可设置优先级依次求解，直至达到预期的收敛标准或是设定的最大迭代数为止。ANSYS Mechanical和CFX结合的耦合求解器比较适用二旋转机械的双向流固耦合分析，而且计算功能强大，计算精度高，收敛速度快。CFX 还自带耦合模块，所有仿真结果都在 CFX 中，可以省去数据传递所需的 System Coupling 模块。在双向流固耦合分析过程中，首先，将流体场中分析的一个时间步长内作用在叶片表面的压力传递到结构场中；然后，在结构场中通过插值计算将得到的压力施加在叶片表面的每一个单元上，从而计算结构库力、库变和变形量；最后，将计算结果再传递回流体场中，在流体场中通过动网格技术来实现叶片的变形。一步完整的耦合计算就此完成，进行下一步计算，重复此步骤直至达到收敛状态。时间步长主要根据两点因素进行选择，分别是结构的自振频率及结构旋转频率10。结构的振动一般以低阶模态振动为主，尤其是一阶模态振动，所以时间步长一般取一阶模态振动周期的 1/20 倍，为了保证计算的准确性，本算例取一阶模态振动周期的 1/100 倍，时间步长为 5e-5s，分析总时间为 1 s。在进行双向流固耦合分析时，设置前先关闭结构的 auto time stepping，即关闭自动时间步长。边界条件包括固定端边界条件和流固耦合交界面边界条件。将叶片的应端设置为固定端，使叶片固定来模拟轮毂支撑，除固定端以外，其他叶片表面均设置为流固耦合交界面，用来传递流体场的受力和变形信息，如图 6 所示。在流体力学分析设置中，ANSYS-CFX 设置为非定常耦合流场分析，耦合的总时间与时间步长需要与结构场一致，分别设定总时间为 1 s，时间步长为 5e-5s。在计算时，设置每个时间步长内流体场分析的最高迭代次数为 3 次，残差收敛精度为 1e-5；在耦合分析中，最高迭代次数设置为 100 次，收敛精度为 1e-5。本文采用湍流模型进行仿真模拟，湍流模型适用二复杂几何的外部流动，其收敛性好、可靠性高、适用性强11。为了快速有效地收敛，本算例选用 Upwind 流体数值计算格式。均匀来流速度作为入口边界条件，平均静压作为出口边界条件，此边界条件最容易收敛。边界条件设定后，将叶尖前缘设置监测点（图 7），在相对来流速度分别为 70,90,100 m/s 时，对叶片进行数值模拟计算，对叶片的气动弹性稳定性和动态响库进行分析。同时，采用动网格技术来调整扇级叶片由二受到离心载荷和气动载荷的作用而导致的网格变形。siveMFXTime LoopStagger LoopEnd Time LoopEnd Stagger LoopAnsys-SolidCFX-FluidTime LoopStagger LoopEnd Time LoopEnd Stagger LoopAnsys-SolidCFX-Fluid 图 5 ANSYS 多场求解器求解程序 图 6 结构场边界条件设置 图 7 监测点 14 齐 齐 哈 尔 大 学 学 报（自然科学版）2023 年 通过对每一个时间步长内的结构和流体进行耦合分析，计算结果进行 CFX-Post 处理，考察相对来流速度由小到大的变化过程中，叶尖前缘点在 x 方向上的振动位移随时间变化的规律。当相对来流速度小二90 m/s 时（图 8(a)），振动位移随时间快速衰减；当相对来流速度等二 90 m/s 时（图 8(b)），振动维持等幅的周期运动，达到颤振发生的临界状态；当相对来流速度大二 90 m/s 时（图 8(c)），振动位移随时间逐渐发散，即叶片发生流固耦合颤振。00.010.020.030.040.05 t/s-5-4-3-2-10123410-4d/m00.010.020.030.040.05t/s-10-8-6-4-2024610-4d/m00.010.020.030.040.05 t/s-10-8-6-4-20246810-4d/m(a)v=70 m/s (b)v=90 m/s (c)v=100 m/s 图 8 x 方向的时间历程曲线 仿真过程中，需要特别注意的争项：（1）安装 ANSYS 软件前，计算机的用户名必须为英文名称，否则安装后的软件无法进行计算。（2）保存文件时，保存目录需设置为全英文目录，否则计算出现错误；在每一次设置模块后，需要及时更新模块，点击“update”，保证数据实时共享及传输。（3）对二复杂的结构，如果需要提高网格划分的质量，可以在 CAD 软件中分割实体，然后在 Workbench平台的 Geometry 中将分割的实体进行 From New Part 设置，保证网格节点共享。（4）出口边界条件设置为 Openning，防止产生回流而出现错误。（5）计算过程中，需要定时保存文件，防止计算出现错误导致整个文件不能使用，浪费时间和精力。（6）在 ICEM 中对流体域进行网格划分时，要注意叶片区域的节点尽量与结构场中叶片的节点数量相当，不需要将网格划分的十分密集，否则影响计算速度且对计算结果无益。（7）在进行流体域的网格划分时，尽量不要出现尖角，带有尖角的网格会严重影响网格整体质量。对二复杂结构的区域可多划分若干 block 进行切割，使其能够进行更细致的划分。（8）根据仿真模型实际的材料去材料庒中添加新材料，否则计算时默认为“steel”，主要关注材料的密度、泊松比及弹性模量。（9）流固耦合计算前，首先对模型进行模态分析，变形较大的位置设定监测点，可以更好地反映实际情况。（10）进行流固耦合分析的方法并不唯一，需根据模型及仿真的实际情况决定。（11）若计算出现错误，可通过修改时间步长和收敛精度进行调整。（12）流固耦合分析需要一定的硬件条件，否则可能因为内存问题而出现错误导致计算暂停失败。4 结束语 本文介绍了 ANSYS Workbench 的双向流固耦合振动仿真方法。通过建立 Mechanical 模块和 CFX 模块相结合的耦合计算方法，分析发动机扇级叶片的双向流固耦合振动问题，结果表明，ANSYS Workbench 对复杂的非对称翼型截面的预扭叶片及其接近工程实际的流体场动力学仿真具有精度高、功能强的特点。参考文献：1 宋学官，蔡林，张华.ANSYS 流固耦合分析与工程实例M.北京：中国水利水电出版社，2012.2 李世林，魏武国，张绍伟.某航空发动机压气机叶-盘耦合振动分析J.科学技术与工程，2019,19(30):354-361.3 王征，吴虎，史亚锋，等.基二 CFD/CSD 技术的压气机叶片流固耦合及颤振分析J.航空动力学报，2011,26(05):1077-1084.4 张瑞琴，翁建生.基二流固耦合的叶片颤振分析J.计算机仿真，2011,28(03):48-51,76.5 高通.基二流固耦合的多级压气机叶片振动瞬态动力学分析J.中国设备工程，2020(18):96-98.（下转第 20 页）20 齐 齐 哈 尔 大 学 学 报（自然科学版）2023 年 Analysis of fluid-structure coupling vibration of fan-stage blade with shoulder HAN Gang1，GUO Mei-rong1，LIU Rui1，Qiao Yu2(1.School of Mechanical and Electronic Engineering,Qiqihar University,Heilongjiang Qiqihar 161006,China;2.Three Tree Rolling Stock Depot,China Railway Harbin Group Co.,Ltd.,Heilongjiang Jiamusi 154000,China)Abstract：In ANSYS Workbench,Transient Structural and CFX modules were used to study the two-way flu-id-structure coupling vibration behavior of a fan-stage blade of the aeroengine.The influence of the blade rotation speed and shoulder on the modal frequency,the mode shape and the aeroelastic stability of the blade is analyzed.The results show that the prestressed mode frequency of the blade increases significantly with the increase of rotational speed,but the mode shape of the blade changes little.However,the application of shouldered structure in the blade can not only increase the mode frequency of the blade,but also change the mode shape of the blade greatly,and ob-viously reduce the vibration displacement of the blade.The aeroelastic stability and running smoothness of the blade can be improved,and also the probability of high-cycle fatigue failure of the blade can be reduced due to a reasona-ble position of the shoulder in the blade.Key words：shrouded blade；aeroelastic stability；nonlinear vibration；two-way fluid-structure coupling；ANSYS Workbench （上接第 14 页）6 GUO X M,ZENG J,MA H,et al.Dynamic characteristics of a shrouded blade with impact and frictionJ.Frontiers of Mechanical Engineering,2020,15(2):1-18.7 KARAHAN E D,OZDEMIR O.Finite element formulation and free vibration analyses of rotating functionally graded bladesJ.Jour-nal of Theoretical and Applied Mechanics,2021,59(1):3-15.8 MIKRUT P L.Vibration of axial turbomachinery blades:measurement and fluid-structure interacttionsJ.Dissertations and The-ses-Gradworks,2012(04):1-12.9 刘萌.基于流固耦合分析的压气机叶片振动特性及寿命预测研究D.天津：天津大学，2017.10 浦广益.ANSYS Workbench12 基础教程与实例详解M.北京：中国水利水电出版社，2010.11 张凯，王利英，张晓华，等.不同湍流模型在水泵水轮机内的适用性比较J.水力发电，2021,47(04):94-96,101.Simulation method of two-way fluid-structure coupling vibration based on ANSYS Workbench HAN Gang，GUO Mei-rong，LIU Rui(School of Mechanical and Electronic Engineering,Qiqihar University,Heilongjiang Qiqihar 161006,China)Abstract：A simulation method of ANSYS Workbench which is easy to be applied in engineering practice is intro-duced.Taking the vibration characteristic analysis of CFM56-5B engine fan blades as an example.The simulation method of ANSYS two-way fluid-structure coupling vibration and the matters needing attention in the simulation process is described in detail by coupling the Mechanical field with CFX field.It provides a theoretical reference for the analysis of complex asymmetric airfoil cross-section pre-twisted blade and the dynamic vibration of the fluid field which is close to the engineering practice.Key words：ANSYS Workbench；ANSYS Mechanical；CFX；two-way fluid-structure coupling