基于声-固耦合算法的储液容器内件模态分析.pdf

李小虎基于声-固耦合算法的储液容器内件模态分析李小虎（森松（江苏）重工有限公司，上海 200135）摘要：采用有限元法对比了干模态、质量附加法及湿模态方法对设备频率的影响，同时在湿模态方法中对比了液体晃动效应对设备频率计算的影响。基于某结晶器结构，通过湿模态方法对充水设备内件的频率及振型进行计算，分析设备内部液体对设备内件频率的影响，对于防止搅拌轴与设备部件发生共振能够起到很重要的作用，为大型薄壁容器及内部有很多内件的设备模态计算提供了准确的计算方法。关键词：内件；湿模态；声-固耦合；满液中图分类号：TQ 051.7文献标志码：A文章编号：2095-817X（2023）04-0059-005反应釜广泛应用于石油、化工、医药等领域，而反应釜中的搅拌装置占据着其中的核心位置 1-2。反应釜在工作中可能发生搅拌装置与容器共振，如何避免这一风险将是评估整个设备安全性的重要指标之一。在评判振动问题中，首先要进行模态分析得到容器的固有频率，然后与搅拌装置的频率进行对比分析。在工程应用上模态分析主要分为干模态和湿模态两种方法，干模态不考虑设备内流体对设备的影响，仅通过附加质量法将设备内介质的重量增加到容器本身，该方法计算简单，但结果精度较低；湿模态即考虑设备周围流体对其频率的影响，能够准确地得出设备在浸水状态下的真实频率，但该方法计算复杂，工作量大 3-5。对于大部分反应釜，容器内部附件少，因此校核设备振动时仅需考虑容器本体与搅拌装置的共振问题，对于此类问题通过简单的质量附加法即可满足计算要求，笔者之前对搅拌装置与容器共振问题、搅拌口密封面偏转角及搅拌疲劳等方面进行了全面分析 6。但对于大型薄壁设备充满水且设备内部有很多内件的情况，由于很难得到准确的附加质量分布，因此采用干模态算法已不能准确得到各部件在流体影响下的真实频率，此时就需要采用湿模态算法，可以直接考虑液体真实的质量效应。目前湿模态算法已应用于各行各业，相关的文献研究很多 7-10，这些文章多集中于计算方法的理论 收稿日期：2023-04-12作者简介：李小虎（1988），男，工程师，从事压力容器应力分析设计工作。公式解析及容器壳体模拟，模型较为简单，对于设备内件的频率分析很少。本文基于声-固耦合方法对大型满水设备的内件进行湿模态分析，并对采用干、湿模态的模拟结果进行对比分析，为设备内件在浸水状态下的频率提供了更为精确的工程计算方法。1 有限元算法简述目前在工程上计算湿模态的方法主要有声-固耦合法和虚拟质量法两种。声-固耦合法属于有限元范畴，它是在拉格朗日坐标系下研究流体、固体及其相互作用的动力学行为的一种方法 9。通过固体和液体的交界面建立耦合关系，需分别划分固体和液体网格，计算精度高但求解慢。虚拟质量法简化了流体和固体结构之间复杂的相互作用，将流体对固体的作用以固定的附加质量形式来体现，避免了液体单元网格的划分。对于单一简单结构，两种方法模拟出的结果相差无几，但对于复杂结构，由于无法准确地得到附加质量矩阵，因此采用虚拟质量法很难得到精确的结果，而声-固耦合算法将会真实考虑液体的动力特性，得到更为准确的结果。本文所分析结构复杂，将采用声-固耦合算法进行有限元模拟。假定条件如下：结构为线弹性，流体为无黏性、可压缩、无旋的理想液体。系统的结构动力学方程如式（1）：MCKFssssnnn+=po （1）对于无阻尼结构，且考虑自由振动时，式（1）可转化为式（2）：59化 工 与 医 药 工 程第44卷第4期2023年8月出版Chemical and Pharmaceutical Engineering 2023，44(4)装备应用与研究MMK0sasnn+=ph （2）式中 Ms系统的质量矩阵；Ma流体附加质量矩阵；n p 节点的加速度矢量；Cs系统的阻尼矩阵；n o 节点的速度矢量；Ks系统刚度矩阵；节点的位移矢量；Fs外界激励矢量。利用式（2）便可求得流固耦合中结构的模态及振型。根据假定条件，可利用声学单元模拟流体，因此结合声学流体及结构方程，可求得声-固耦合有限元方程 3，11，如式（3）：MMMpKKKupF000ssaasfsfsn+=pp=;GGEEE （3）式中 Msa流体的等效耦合质量矩阵；Kf流体刚度矩阵；Kfs等效刚度矩阵；p流体压力。2 计算方法及边界条件本文首先对比了干模态、质量附加法及湿模态方法对设备频率的影响，同时在湿模态方法中考虑液体晃动效应的影响，进而得到工程上既精确又快捷的计算方法。质量附加法与干模态计算方法相同，都是将质量等效到固体结构上，区别在于干模态仅考虑设备本体的质量，而质量附加法将容器内介质的质量也等效于固体结构上。其次以某结晶器简化结构为研究对象，通过湿模态方法对充水设备内件频率及振型进行计算，分析设备内部液体对设备内件频率的影响，为大型薄壁容器且内部有很多内件的设备模态计算提供了很好的计算方法。分析过程中需注意如下事项：在进行湿模态分析时，需要建立设备内部液体模型，由于液体和固体是两个单独的实体，在模型建立时需在两者交界面上保持单元节点两者共用，这样才能进行流-固耦合。流体部分采用声学单元 FLUID30，容器部分采用结构单元 SOLID186。对于流体材料，需要输入材料密度及流体声速，对于容器材料，需要输入材料弹性模量、泊松比和密度。设置边界条件时，选取流体单元中流固交界面上的节点，执行 FSI 命令，标记耦合界面。3 模态计算结果3.1 干、湿模态及液体自由晃动对设备频率的影响在对设备内件频率进行分析研究前，本节先采用一圆柱形满水容器进行模拟，分别对比附加质量法、湿模态法（分别考虑液体自由晃动和不考虑液体自由晃动两种情况）下容器固有频率的区别。设备主要几何参数如下：半径 2 500 mm，厚度 10 mm，高度 6 000 mm，容器内液体质量为 117 t，容器本体质量为 9 t。如表 1 所示，分别列出了各计算方法下容器本体前 6 阶频率。表 1 各计算方法下容器本体前 6 阶频率Tab.1 The top 6 frequencies of cylinder under each calculation method阶次 n干模态/Hz质量附加法/Hz湿模态（液体自由晃动）/Hz湿模态（忽略液体晃动）/Hz128.5427.66211.35311.325228.5467.66311.35311.326330.2808.12811.57311.554430.2828.12911.57411.555530.3378.14412.70112.651630.3408.14512.70412.654从表 1 可以看出，不考虑液体质量下的容器干模态频率要比考虑液体质量下容器频率高很多，而对于大型搅拌容器，由于搅拌频率很低，若直接按照干模态频率进行评估，不能准确地评估设备与搅拌发生共振的风险。质量附加法得到的频率要比湿模态法小，该方法虽然没有湿模态法结果准确，但该方法简单，对于大多数搅拌设备，本身的评定能够得到一个保守的结果。用声-固耦合湿模态法进行分析时，考虑液体表面自由晃动和不考虑液体晃动所得到的设备本体频率几乎是一样的，唯一的区别在于考虑液体自由晃动可以得到流体的低频率模态（如图 1所示）。若要研究容器内液体晃动问题，可以采用考虑液体自由晃动的声-固耦合法，若仅需要得到准确的容器自身模态，则无需考虑液体表面的自由晃动，这样能大大减少计算时间。3.2 湿模态法计算充水设备内件模态上节已经表明对于设备内液体质量占比很大时，60化 工 与 医 药 工 程第44卷第4期2023年8月出版Chemical and Pharmaceutical Engineering 2023，44(4)干模态算法得到的设备模态结果误差很大，而对于设备内部复杂的内件，通过质量附加法也是不现实的，因为无法准确地知道介质质量分布作用，此时就只能通过湿模态进行求解，以得到准确的结果。本文以一个大型薄壁设备作为案例，该结构主要特点为：容器尺寸大、壁厚薄、内件多且复杂、内部充满水。设备基本参数如下：直径 7 000 mm，壁厚14 mm，高度 10 000 mm，设备模型如图 2 所示。设备主体材料为 SB-688 N08367，密度为 8 060 kg/m3，内部液体密度为 1 500 kg/m3，设备自重 93 吨，满水重量约 700 吨。该模型主要分析内件频率，因此略去设备上部无内件部分，约束裙座支撑底面。其他边界条件见本文第 2 节。如图 3 所示，分析模型中液体部分也需要进行网格建立。由于模型尺寸较大，若网格数量太多，将给计算带来很大的困难。在湿模态计算中，液体部分网格的作用仅是为了将液体质量的影响传递到对应的部件上，因此可以适当减少液体部分的网格数量，保证设备本身的网格精度即可，模型总网格数量为图 1 考虑液体晃动效应的湿模态下容器及液体模态Fig.1 Cylinder and liquid modes of wet mode considering liquid sloshing effect a.容器一阶模态 f=11.353 Hzb.液体晃动模态 f=0.55 Hz图 2 设备有限元模型Fig.2 Finite element model of Equipment图 3 有限元模型网格Fig.3 Mesh of finite element model3 953 652。用湿模态算法对该设备内部充满水的工况进行计算，共计算设备的前 30 阶模态，得到各部件的频率。如图 4 所示，设备第一阶模态发生在设备内部内件位置，在满水状态下内件的频率仅为 3.09 Hz（图 4a所示），相比于空重下内件的频率 7.359 Hz（图 4b 所示），湿模态频率仅为干模态频率的42%。如图5所示，一直到 20 阶以后，模态才出现在设备本体上，在满水状态下本体的频率为 13.216 Hz，而在空重下本体的 1 阶频率为 35.09 Hz，同样的湿模态频率仅为干模态频率的 38%。从上述的计算结果可以看出，对于设备内部的复杂内件结构，在内部有介质影响时，干模态得到的结果已经严重偏离实际，同时由于无法准确地进行质量分配，质量附加法也无法使用，因此采用湿模态法才能得到最真实的结果，这对于防止搅拌轴与设备部件发生共振起到了很重要的作用。4 结论本文采用有限元法对比了干模态、质量附加法及湿模态方法对设备频率的影响，同时在湿模态方法61李小虎.基于声-固耦合算法的储液容器内件模态分析中研究了液体晃动效应对设备频率计算的影响，液体晃动效应主要得到低频的自由液面波动，对设备本体的频率几乎没有影响，因此在工程上计算设备频率时，既精确又快捷的计算方法，即为忽略液体晃动效应的声-固耦合湿模态法。文中以某结晶器简化结构为研究对象，通过湿模态方法对充水设备内件的频率及振型进行计算，分析设备内部液体对设备内件频率的影响，通过分析对比可知，对于设备内复杂内件结构，在内部有介质影响时，干模态得到的结果已经严重偏离实际，因此采用湿模态法才能得到最真实的结果，这对于防止搅拌轴与设备部件发生共振起到了很重要的作用。为大型薄壁容器及内部有很多内件的设备模态计算提供了准确的计算方法。参考文献刘庆荣，金海冰，丁全有，等.一种新型反应釜的结构设计 J.1 石油化工设备，2016，45（2）：37-40.佟立军.立式加压反应釜搅拌装置振动分析及解决方法探讨2 J.有色设备，2019，33（1）：12-16.刘瑞骏，郝志勇，郑旭，等.改进的湿模态法在流固耦合中的3 应用 J.振动与冲击，2017，36（22）：199-204.杨鸣，王辉，段玉康，等.基于声-固耦合算法的储液容器湿4 模态分析 J.四川兵工学报，2015，36（5）：152-154.程玉鑫，周力，洪明，等.舱内液体附加质量对船体振动影响5 的模型实验研究 J.中国舰船研究，2008，3（2）：64-67.李小虎.搅拌反应釜振动分析及结构优化 J.化工机械，6 2020，47（1）：37-40.Z.Wu，W.Zhou，H.Li.Modal Analysis for Filament Wound 7 Pressure Vessels Filled with Fluid J.Composite structures，2010，92：1994-1998.Kwon，Y.W.，Plessas，S.D.，Numerical Modal Analysis 8 of Composite Structures Coupled with WaterJ.Composite structures，2014，116：325-335.薛杰，何尚龙，杜大华，等.充满容器流固耦合模态仿真分析9 研究 J.火箭推进，2015，41（1）：90-97.李俊儒.基于声固耦合算法的塔釜持液塔器湿模态分析 J.化10 工设备与管道，2021，58（6):20-23.高海司，裴吉，袁寿其，等.蜗壳式混流泵叶轮部件湿模态流11 固耦合分析 J.流体机械，2018，46（1）：46-51.图 5 设备本体 1 阶模态频率Fig.5 The first mode frequency of the equipment图 4 内件 1 阶模态频率Fig.4 The first mode frequency of the internalsa.内部充满水内件的一阶模态（f=3.099 Hz）a.满水状态本体 1 阶模态（f=13.216 Hz）b.内部无水内件的一阶模态（f=7.359 Hz）b.空重下本体 1 阶模态（f=35.091 Hz）62化 工 与 医 药 工 程第44卷第4期2023年8月出版Chemical and Pharmaceutical Engineering 2023，44(4)Modal Analysis of Vessel Internals Filled With Liquid Based on Solid-Structure Coupling MethodLi Xiaohu（MORIMATSU Jjangsu Heavy Industry Co.,Ltd.,Shanghai 200135,China）Abstract:Compared the effects of dry mode,mass addition method and wet mode method on the vessel frequency by using FEA.The influence of liquid slosh effect on the vessel frequency calculation was analyzed in the wet mode method.The frequency and mode shape of the vessel internals filling with liquid were calculated by the wet mode method,and the influence of the liquid inside the vessel on the internal frequency was analyzed.It plays a very important role in preventing the resonance between the stirring shaft and the equipment,and provides an accurate calculation method for the modal calculation of the large thin-walled vessel and the equipment with many internal parts.Keywords:internal parts;wet mode;acoustic-structure coupling;filled with liquid欢迎订阅化工设备与管道（邮发代号：4669）化工设备与管道，创刊于1964年，是经国家科学技术部、国家新闻出版广电总局正式批准出版的专业技术期刊，曾用名化工设备简讯、化工设备设计，是国内同类期刊中历史最为悠久的知名期刊之一，曾多次获得上海市、原化工部、中国石油和化学工业协会等多项部省级荣誉与奖励。本刊为中文核心期刊、中国科技核心期刊，国家认定的学术期刊（A 类）。化工设备与管道由中国石油化工集团有限公司主管，中石化上海工程有限公司主办，全国化工设备设计技术中心站、中国石化集团有限公司机械技术中心站、全国锅炉压力容器标准化技术委员会压力管道分技术委员会、全国化工工业炉设计技术中心站承办。国内统一刊号：CN31-1833/TQ；国际标准刊号：ISSN1009-3281；邮发代号：4-669。化工设备与管道主要栏目有：综述、单元设备、压力容器、管道与管件、机械与密封等。读者对象为石油、化工、石化及相关行业从事设备与管道工程设计、建设、制造、采购、检测、维修工作的广大科研、技术与生产人员。本刊为双月刊，A4 开本，全彩印刷，自1964 年创刊以来已出版 300 余期，发行量大，影响面广。全国各地邮局均可办理订阅，10.00 元/期，全年定价 60 元。地 址：上海市浦东新区张杨路769号邮政编码：200120电 话：[phone]、52136992、62488580邮 箱：63李小虎.基于声-固耦合算法的储液容器内件模态分析