附加质量法和FSI法用于弧形闸门动力特性分析的适用性.pdf

附加质量法和FSI法用于弧形闸门动力特性分析的适用性胡友一，王珏*，陆中伟，郑赟（河海大学机电工程学院，江苏常州213022）摘要：为研究附加质量法和流固耦合界面（FSI）法用于弧形闸门动力特性分析的适用性，首先通过建立弧形闸门有限元模型，采用收敛性分析的方法研究了FSI法计算时流体长度的合理取值，保证了流体-弧形闸门系统自振频率的计算精度。其次，从弧形闸门水头高度和刚度的角度出发，比较了采用附加质量法和FSI法对弧形闸门流固耦合效应的模拟结果，讨论了采用简化的附加质量法的适用性。最后，根据有限元计算的数值结果讨论了流固耦合效应对弧形闸门动力特性的影响。结果表明：1）FSI法中流体长度取2倍水头高度可代替实际流体的无界性；2）在水头较低且闸门刚度较大的情况下，采用附加质量法模拟流固耦合效应能在保证计算精度的基础上提高计算效率。随着水头增高，精度有所下降。此时，弧形闸门的自振频率明显降低，并且主要的变形由支臂转移到面板，有必要采用FSI法代替附加质量法来确保计算精度。关键词：弧形闸门；流固耦合；附加质量法；FSI法；自振频率中图分类号：U641；TV3文献标志码：A文章编号：2095-7874（2024）01-0001-06doi：10.7640/zggwjs202401001Applicability of added-mass method and FSI method for dynamiccharacteristics analysis of radial gatesHU You-yi,WANG Jue*,LU Zhong-wei，ZHENG Yun（College of Mechanical and Electrical Engineering,Hohai University,Changzhou,Jiangsu 213022,China）Abstract：To investigate the applicability of the added-mass method and the fluid-structure interaction interface(FSI)methodfor analyzing the dynamic characteristics of radial gates,firstly,by establishing a finite element model of radial gate,thereasonable value of fluid length in FSI method was studied by convergence analysis,ensuring the calculation accuracy ofnatural frequency of fluid-radial gate system.Secondly,the simulation results of the fluid-structure interaction effect of theradial gate using the added-mass method and FSI method are compared from the perspective of the head height and stiffness ofthe radial gate,and the applicability of the simplified added-mass method was discussed.Finally,the influence of fluid-structure interaction effect on the dynamic characteristics of radial gate was discussed based on the numerical results of finiteelement calculations.The results show as follows:1)the unboundedness of the fluid can be simulated by the length of twice thewater head in the FSI method;2)In the case of lower water head and greater stiffness,using the added-mass method tosimulate the fluid-structure interaction effect can improve calculation efficiency while ensuring calculation accuracy.As thewater head increases,the accuracy decreases.At this point,the natural vibration frequency of the radial gate significantlydecreases,and the main deformation is transferred from the support arm to the panel.It is necessary to use the FSI methodinstead of the added-mass method to ensure calculation accuracy.Key words：radial gate;fluid-structure interaction;added-mass method;FSI method;natural vibration frequency中国港湾建设China Harbour Engineering第44卷第1期2024年1月Vol.44No.1Jan.2024收稿日期：2023-06-15修回日期：2023-11-13基金项目：国家自然科学基金项目（52178474）；中国博士后基金（2018M632216）作者简介：胡友一（1998），男，河南鲁山人，硕士研究生，主要从事抗震减灾研究。*通讯作者：王珏，E-mail:中国港湾建设2024年第1期0引言弧形闸门因其具有启门力小、过水条件好、不设门槽等特点，是水工建筑物中运用最广泛的门型之一1。与平面闸门相比，弧形闸门具有特殊的结构形式，其支臂结构以及门叶上下部悬臂端在动水荷载作用下更易发生振动破坏。为确保闸门安全，有必要对弧形闸门进行动力分析，在设计计算中通过计算闸门的自振频率，并与作用水流的脉动频率相比较，使闸门自振频率尽量远离水流的脉动频率2。由于弧形闸门和流体之间存在相互作用，因此很有必要对考虑流固耦合效应的弧形闸门动力特性计算方法的适用性进行研究。闸门的流固耦合动力分析可采用原型观测3-4、模型试验5-6、数值模拟7-9等方法，其中数值模拟方法相比于原型观测和模型试验更加便捷、高效，近年来在闸门动力学分析中得到广泛应用。曹青等2采用Nastran软件建立了考虑弧门流固耦合效应的附加质量数值模型，分析表明流体对弧形闸门的动力特性有显著的影响；上述研究中附加质量法是一种考虑流固耦合效应的简化动力分析计算方法，计算效率相对较高，广泛应用于闸门和流体相互作用的有限元模拟中，其中模型中的附加质量通常采用Westergaard公式10来计算。而通过建立流体模型的流固耦合界面法（Fluid-SolidInterction，以下简称“FSI”）能够准确地描述流体的运动，能更真实地反映流体和结构的动力相互作用。但是由于流体的无界性，建立足够长的流体模型是保证有限元分析精度的一个关键问题。Buldgen等11-12通过比较附加质量法和FSI法计算得到的平面闸门在地震作用下的动水压力以及动力响应，发现对于不完全刚性的结构，附加质量法的计算结果往往是不精确的，有必要采用FSI法计算以保证分析精度。与平面闸门不同，弧形闸门结构更加复杂。因此，有必要在保证弧形闸门问题中FSI法流体模型长度合理取值的基础上，来验证附加质量法计算得到流体和闸门系统动力特性的计算精度，从而保证考虑流固耦合效应的闸门动力特性分析的准确性。在此基础上，进一步采用拟静力法13、振型分解反应谱法14、动力时程法15-16等对考虑流固耦合效应的闸门开展时域响应分析，研究流固耦合效应对闸门动力响应的影响规律。本文以某工程弧形闸门为例，利用大型有限元分析软件ANSYS建立数值模型。首先为避免建立流体所需要的大量单元而降低计算效率，讨论了FSI法中流体有限元模型长度的合理取值，进而从弧形闸门刚度和水头高度讨论了附加质量法的计算精度，在此基础上研究了流固耦合效应对弧形闸门动力特性的影响。1分析方法1.1附加质量法附加质量法是由Westergaard10提出来的，其假设流体是无限长且不可压缩的，基本思想是将动水压力等效成质量附加在结构上，从而简化流体对结构的作用，据此弧形闸门在动水压力作用下的控制方程可写为：（M0+Mf）俟（t）+（C0+Cf）u（t）+（K0+Kf）u（t）=F（t）（1）式中：M0、C0、K0分别为弧形闸门的质量矩阵、阻尼矩阵、刚度矩阵；Mf、Cf、Kf分别为流体的附加质量矩阵、附加阻尼矩阵、附加刚度矩阵；俟（t）、u（t）、u（t）分别为弧形闸门加速度矩阵、速度矩阵和位移矩阵；F（t）为水流脉动压力作用在结构上引起的荷载矩阵。对于附加质量矩阵Mf，可以通过在ANSYS中对弧形闸门与流体交界面上的节点施加质量单元，每个单元的数值根据Westergaard公式赋予：Mi=78籽fAiHhi姨（2）式中：Mi为对应节点的附加质量；籽f为流体密度；Ai为流体作用下对应节点关联单元的有效面积；H为总水头高度；hi为对应节点与水面之间的距离。对于结构阻尼矩阵C0，本文采用瑞利阻尼，公式如下：C0=4仔孜fifj（fi+fj）M0+孜仔（fi+fj）K0（3）式中：孜为结构阻尼比，本文取0.05；fi、fj为结构第i、j阶自振频率，本文取i=1，j=2。在进行弧形闸门考虑流固耦合的结构动力分析时，附加阻尼Cf和附加刚度Kf对闸门动力特性的影响微乎其微，因此，这里忽略水体附加阻尼Cf与附加刚度Kf的影响，将水体附加质量对闸门动力特性的影响作为考虑的重点，因此运动方程式（1）可简化为：（M0+Mf）俟（t）+C0u（t）+K0u（t）=F（t）（4）2窑窑2024年第1期1.2流固耦合界面（FSI）法FSI法是指对结构和流体分别建立运动方程，采用数值方法对其耦合求解，并将计算结果通过流固交界面进行信息交互传递的一种计算方法。FSI法需要在确定流域范围的基础上对流体进行建模，由于流体的无界性，需要建立足够长的流体模型才能保证有限元分析结果的精度，从而反映实际流体的作用。在ANSYS中采用流体单元建立流体模型，并施加边界条件。而定义流固耦合界面是弧形闸门进行流固耦合模拟仿真时的关键问题，通过对弧形闸门与流体接触的流体单元的FSI命令来实现。最后采用非对称矩阵法（UNSYMMETRIC）进行模态分析求解17。本文为避免建立流域需要划分大量单元而增加计算难度，计算了不同流体模型长度以及水头高度时弧形闸门的自振频率，当结果收敛时即可确定流域的最小取值范围，从而简化流体模型。2分析对象及其模型2.1结构布置及基本参数本文以某工程弧形闸门为例，闸门尺寸为6.9m伊9.21 m（宽伊高），半径为13.5 m，设计水头6 m，底槛和支铰高程分别为138.42 m和128.72 m。该闸门为主横梁式露顶弧形闸门，闸门采用板梁结构等高布置，水平向设置2根实腹板主横梁以及14根水平次梁（包含顶、底梁），垂直向设置5道实腹式隔板及2道边梁，据此组成梁格与面板直接焊接。闸门直支臂与主横梁采用螺栓连接构成主框架，上下支臂中间设置一系列斜撑。闸门结构和流体的材料属性如表1所示。2.2有限元模型1）弧形闸门模型将弧形闸门面板、主横梁、纵隔板、边梁、支臂划分为Shell181板单元，水平次梁、顶梁、底梁、支臂斜撑划分为Beam189梁单元，活动铰划分为Solid186实体单元，弧形钢闸门有限元计算模型如图1所示。其中节点和单元总数分别为14 575、11 922个，图中x轴（正）为顺水流方向，y轴（正）竖直向上，z轴（正）为横河向指向左岸。在弧形闸门活动铰的内表面约束x、y、z方向的移动自由度以及x、y方向的旋转自由度，仅放松可以绕支铰转动的z向旋转自由度；约束面板底边y方向移动自由度。2）流体模型流体分别采用FSI法和附加质量法建立，水头高度为6 m，如图2所示。FSI法流体模型采用Fluid30流体单元建立，模型长度根据弧形闸门水头高度按一定比例选取。该流体模型的节点和单元总数分别为23 977、51 777个，在模型底面约束y方向移动自由度，两个侧面约束z方向移动自由度。附加质量模型采用MASS21单元建立，施加在面板与流体接触面上的所有节点上，单元实常数根据公式（2）的计算结果来定义。3计算结果分析3.1FSI法流体收敛性分析采用附加质量法和FSI法计算得到弧形闸门在不同水头下的自振频率，计算结果如表2所示。其中h为水头高度，L为流体模型长度。从表2可以看出，FSI法中流体模型的长度由h增加到1.5h表1材料属性Table 1Material properties参数钢闸门铸铁闸门流体弹性模量E/（105MPa）2.101.22泊松比滋0.300.25密度籽s/（kgm-3）7 8007 300声速c/（ms-1）1 460密度籽f/（kgm-3）1 000（b）闸门背水面（a）闸门迎水面图1弧形闸门有限元模型Fig.1Finite element model of radial gate（a）FSI模型（b）附加质量模型图2流体有限元模型Fig.2Finite element model of fluidxyzxyz胡友一，等：附加质量法和FSI法用于弧形闸门动力特性分析的适用性3中国港湾建设2024年第1期3.2附加质量法适用性分析图3将采用FSI法（L=2h）计算所得弧形闸门在不同水头下的自振频率与采用附加质量法的计算结果进行了对比，结果表明：附加质量法在计算3 m和6 m水头下的弧形闸门前六阶自振频率时，能得到较高的精度。但随着水头增大到9 m，采用附加质量法计算得到的闸门高阶频率与FSI法具有一定偏差，其中最大的偏差近20%，平均偏差近10%。因此对于低水头的弧形闸门，为了提高计算效率，可以采用附加质量法模拟流体；但是对于高水头的弧形闸门，为了保证计算精度建议采用FSI法。为了研究弧形闸门刚度对附加质量法精度的影响，本算例基于FSI法（L=2h）和附加质量法计算对比了弧形铸铁闸门（弹性模量1.22伊105MPa）和钢闸门（弹性模量2.10伊105MPa）在6 m水头下的自振频率，如表3所示。结果表明，随着弧形闸门刚度减小，附加质量法相比于FSI法各阶自振频率偏差明显增大。因此，对于刚度较小、柔度较大的弧形闸门，附加质量法的计算精度较低，建议采用FSI法模拟流体。3.3弧形闸门动力特性参数分析主要从弧形闸门在不同情况下的自振频率和振型来分析流固耦合效应的作用效果。表2反映了考虑流固耦合效应下水头高度对弧形闸门自振频率的影响。通过对比自振频率的结果，可以看出：随着水头的增加，弧形闸门的各阶自振频率均有不同程度的下降。其中，闸门前两阶的频率变化相对较小，原因是由于闸门本身跨度较小、刚度较大，附加质量对自振频率影响相对较小。时，弧形闸门所得同阶自振频率相差较大，尤其是在9 m水头下第二阶的自振频率相差最大且超过了40%；而当流体长度达到2h时，前六阶频率的计算结果收敛。因此，为避免建立流体所需要划分大量单元而增加计算难度，本文建议FSI法模型流体长度取L=2h来代替实际流体的无界性。表2基于附加质量法和FSI法的不同水头下弧形闸门自振频率Table 2Natural vibration frequencies of radial gates under different water heads based on the added-mass and FSI methods阶数无水3 m水头6 m水头9 m水头FSI法附加质量法FSI法附加质量法FSI法附加质量法L=hL=1.5hL=2hL=hL=1.5hL=2hL=hL=1.5hL=2h11.4591.4011.4121.4121.4011.3991.4101.4111.3951.3971.4081.4081.38324.3153.1804.2094.2804.2763.0654.1484.2594.1372.8694.0254.1364.02636.6994.5316.3676.4296.1704.2965.6175.7585.6144.1105.5635.6205.602412.73010.22010.87011.06010.5109.99110.49010.7009.3628.7039.2789.4837.633514.05010.66010.99011.09010.95010.38010.61010.7309.8789.3509.4369.5268.029615.61012.00012.06012.07011.12011.85011.86011.86010.81011.23011.25011.2609.679图3附加质量法与FSI法计算所得频率偏差Fig.3Frequency deviation calculated by the added-massand FSI methods3 m水头6 m水头9 m水头20151050自振频率阶数254163表3不同刚度弧形闸门基于FSI和附加质量法在6 m水头下的自振频率Table 3Natural vibration frequencies of radial gate withdifferent stiffness under 6 m water head basedon the added-mass and FSI methods阶数E=2.10伊105MPaE=1.22伊105MPaFSI法/Hz附加质量法/Hz误差/%FSI法/Hz附加质量法/Hz误差/%11.4111.3951.11.2101.1812.424.2594.1372.93.2683.1025.135.7585.6142.54.7344.3488.2410.7009.36212.58.5697.23215.6510.7309.8787.99.6648.69310.0611.86010.8108.910.2109.11010.8Hz4窑窑2024年第1期该闸门从无水到9 m水头后四阶的自振频率降幅分别为14.3%、40.0%、42.9%、38.0%。由此可见，流体与结构的流固耦合作用对于弧形闸门自振的影响不可忽略。从闸门振型的角度出发，讨论了弧形闸门在无水和9 m水头工况下的前五阶振型，表4阐述了相应的振型特点，图4图8为相应的振型图。结果表明：1）在第1阶模态下，2种工况下弧形闸门均表现为整体（包括支臂和门叶）的侧向运动，流体对结构的局部变形影响较小；在高阶模态下，弧形闸门基本以支臂和面板的振动为主，这是由于面板和支臂的刚度相对较低，在动水压力作用下可能引起强烈振动；2）无水情况下闸门的低阶模态均呈现出支臂的变形，其中前三阶表现为整体运动，四到五阶表现为支臂的局部变形；随着水头增至9 m，闸门主要的变形由支臂转移到面板，究其原因在于面板上的附加质量对支臂的偏心矩减小，其中上面板（上主梁上方）相当于悬臂梁，极易引起其发生径向振动。4结语本文讨论了附加质量法和FSI法在考虑弧形闸门流固耦合动力特性分析中的计算精度问题，并从弧形闸门水头高度和刚度的角度提出了附加质量法的适用性建议。1）为避免建立流域所需要划分大量单元而增加计算难度，本文建议FSI模型流体长度取到2倍的水位高度来模拟实际流体的无界性。对于低水头且刚度较大的弧形闸门，为了提高计算效率，可以采用附加质量法考虑流固耦合效应，但对于高水头或者低刚度的弧形闸门，为了保证计算精度建议采用FSI法。表4弧形闸门振型特点Table 4Mode of vibration characteristics of radial gate阶数无水9 m水头12闸门整体扭转面板上端径向（沿逆水流方向）运动3闸门整体纵向（沿y轴正方向）运动主梁上部面板上端扭转4支臂反对称弯曲上下主梁中部面板径向运动5支臂对称弯曲上下主梁中部面板两侧径向运动闸门整体侧向（沿z轴正方向）运动闸门整体侧向（沿z轴正方向）运动图6第三阶振型图Fig.6The 3rd mode of vibration diagram图5第二阶振型图Fig.5The 2nd mode of vibration diagram图4第一阶振型图Fig.4The 1st mode of vibration diagram（b）9 m水头工况（a）无水工况图8第五阶振型图Fig.8The 5th mode of vibration diagram（b）9 m水头工况（a）无水工况（b）9 m水头工况（a）无水工况图7第四阶振型图Fig.7The 4th mode of vibration diagram（b）9 m水头工况（a）无水工况（b）9 m水头工况（a）无水工况胡友一，等：附加质量法和FSI法用于弧形闸门动力特性分析的适用性5中国港湾建设2024年第1期2）弧形闸门的自振频率随着水头的增高而降低，其振型以面板和支臂的振动为主，闸门主要的变形由支臂转移到面板，这是因为面板上的附加质量对支臂的偏心矩减小，其中上面板（上主梁上方）极易发生径向振动。参考文献：1王正中，张雪才，刘计良.大型水工钢闸门的研究进展及发展趋势J.水力发电学报，2017，36（10）：1-18.WANG Zheng-zhong,ZHANG Xue-cai,LIU Ji-liang.Advancesand developing trends in research of large hydraulic steel gatesJ.Journal of Hydroelectric Engineering,2017,36(10):1-18.2曹青，王光纶，才君眉.弧形闸门自振特性研究J.水利水电技术，2001（5）：16-19.CAO Qing,WANG Guang-lun,CAI Jun-mei.Research on the nat原ural vibrating properties of radial gatesJ.Water Resources andHydropower Engineering,2001(5):16-19.3徐建光，李长和，王延斌，等.龙羊峡水电站表孔弧形闸门原型试验研究J.水力发电，2004（5）：28-31.XU Jian-guang,LI Chang-he,WANG Yan-bin,et al.Prototypetest and research on arc steel gate in Longyangxia Hydropower Sta原tionJ.Water Power,2004(5):28-31.4严根华，陈发展，赵建平.表孔弧形闸门流激振动原型观测研究J.水力发电学报，2006（4）：45-50.YAN Gen-hua,CHEN Fa-zhan,ZHAO Jian-ping.Prototype ob原servation study on the flow-induced vibration of surface radial gateJ.Journal of Hydroelectric Engineering,2006(4):45-50.5张文远，杨帆，章晋雄，等.白鹤滩水电站泄洪洞工作弧门流激振动模型试验研究J.噪声与振动控制，2018，38（S2）：422-426.ZHANG Wen-yuan,YANG Fan,ZHANG Jin-xiong,et al.Experi原mental study on flow induced vibration of spillway tunnel in Bai原hetan Hydropower StationJ.Noise and Vibration Control,2018,38(S2):422-426.6李忠献，吴堃，石运东，等.水-振动台相互作用竖向动力试验研究J.地震工程与工程振动，2019，39（3）：1-7.LI Zhong-xian,WU Kun,SHI Yun-dong,et al.Vertical dynamictests on the interaction between water and shaking tableJ.Earth原quake Engineering and Engineering Dynamics,2019,39(3):1-7.7李桑军，秦战生.基于ANSYS的流固耦合弧形闸门振动特性研究J.水力发电，2018，44（1）：64-67.LI Sang-jun,QIN Zhan-sheng.Research on vibration characteris原tics of radial gate based on ANSYS by considering fluid-structureinteractionJ.Water Power,2018,44(1):64-67.8刘竹丽，陈赟，伊元忠.不同开度下平面钢闸门流固耦合数值模拟研究J.人民黄河，2020，42（11）：79-83.LIU Zhu-li,CHEN Yun,YI Yuan-zhong.Fluid structure interac原tion numerical simulation of plane steel gate with different openingdegreesJ.Yellow River,2020,42(11):79-83.9仵凡，王正中，张雪才，等.弧形钢闸门地震响应的数值分析C/中国力学大会论文集（CCTAM 2019）.杭州：中国力学学会，2019：13.WU Fan,WANG Zheng-zhong,ZHANG Xue-cai,et al.Numericalanalysis of seismic response of radial steel gateC/CCTAM 2019.Hangzhou:Chinese Society of Theoretical and Applied Mechanics,2019:13.10 WESTERGAARD H M.Water pressures on dams during earth原quakesJ.Transactions of the American Society of Civil Engineers,1933,98(2):0004496.11 BULDGEN L,CAPRACE J D,RIGO P,et al.Investigation of theadded mass method for seismic design of lock gatesJ.EngineeringStructures,2017,131:380-393.12 BULDGEN L,RIGO P,LE SOURNE H.A simplified analyticalmethod to evaluate the seismic pressure on plane lock gatesJ.Engineering Structures,2015,100:522-534.13 DANIELL W E,TAYLOR C A.Seismic analysis and dynamic test原ing of a spillway radial gateC/Proceeding 12th World Conferenceon Earthquake Engineering,2000.14余雄，杨怀平，王曾璇，等.挡水闸的地震响应分析J.水利水运科学研究，2000（3）：22-26.YU Xiong,YANG Huai-ping,WANG Zeng-xuan,et al.SeismicresponseanalysisofsluiceJ.Journal ofNanjingHydraulicResearchInstitute,2000(3):22-26.15 DAVEY R,HO K V,WALKER J.Seismic assessment and remedi原ation of the Aviemore Dam spillway and sluice gatesC/2007NZSEE Conference,2007.16 FARIDMEHR I,FAROKHI NEJAD A,BAGHBAN M H,et al.Numerical and physical analysis on the response of a dam忆s radialgate to extreme loading performanceJ.Water,2020,12(9):2425.17杨吉新，张可，党慧慧.基于ANSYS的流固耦合动力分析方法J.船海工程，2008，37（6）：86-89.YANG Ji-xin,ZHANG Ke,DANG Hui-hui.Dynamic analysismethod for fluid-structure interaction based on ANSYSJ.Ship&Ocean Engineering,2008,37(6):86-89.6窑窑