埋地大口径蝶阀基于刚度、密封性能的最大弯曲载荷研究
阀门的基础知识及选择
筑龙网 W W W .Z H U L O N G .C OM 阀门蝶阀蝶阀的蝶板安装于管道的直径方向。
在蝶阀阀体圆柱形通道内，圆盘形蝶板绕着轴线旋转，旋转角度为0°~90°之间，旋转到90°时，阀门则牌全开状态。
蝶阀结构简单、体积小、重量轻，只由少数几个零件组成。
而且只需旋转90°即可快速启闭，操作简单，同时该阀门具有良好的流体控制特性。
蝶阀处于完全开启位置时，蝶板厚度是介质流经阀体时唯一的阻力，因此通过该阀门所产生的压力降很小，故具有较好的流量控制特性。
蝶阀有弹密封和金属的密封两种密封型式。
弹性密封阀门，密封圈可以镶嵌在阀体上或附在蝶板周边。
采用金属密封的阀门一般比弹性密封的阀门寿命长，但很难做到完全密封。
金属密封能适应较高的工作温度，弹性密封则具有受温度限制的缺陷。
如果要求蝶阀作为流量控制使用，主要的是正确选择阀门的尺寸和类型。
蝶阀的结构原理尤其适合制作大口径阀门。
蝶阀不仅在石油、煤气、化工、水处理等一般工业上得到广泛应用，而且还应用于热电站的冷却水系统。
常用的蝶阀有对夹式蝶阀和法兰式蝶阀两种。
对夹式蝶阀是用双头螺栓将阀门连接在两管道法兰之间，法兰式蝶阀是阀门上带有法兰，用螺栓将阀门上两端法兰连接在管道法兰上。
阀门的强度性能是指阀门承受介质压力的能力。
阀门是承受内压的机械产品，因而必须具有足够的强度和刚度，以保证长期使用而不发生破裂或产生变形。
球阀球阀是由旋塞阀演变而来。
它具有相同的旋转90度提动作，不同的是旋塞体是球体，有圆形通孔或通道通过其轴线。
球面和通道口的比例应该是这样的，即当球旋转90度时，在进、出口处应全部呈现球面，从而截断流动。
球阀只需要用旋转90度的操作和很小的转动力矩就能关闭严密。
完全平等筑龙网 W W W .Z H U L O N G .C OM 的阀体内腔为介质提供了阻力很小、直通的流道。
通常认为球阀最适宜直接做开闭使用，但近来的发展已将球阀设计成使它具有节流和控制流量之用。
基于CFD
第８卷㊀第６期２０２３年１１月气体物理ＰＨＹＳＩＣＳＯＦＧＡＳＥＳＶｏｌ.８㊀Ｎｏ.６Ｎｏｖ.２０２３㊀㊀ＤＯＩ:１０.１９５２７/ｊ.ｃｎｋｉ.２０９６￣１６４２.１０８８基于ＣＦＤ/ＣＳＤ耦合的火箭跨声速气动阻尼特性分析李泳德ꎬ㊀郭㊀力ꎬ㊀季㊀辰(中国航天空气动力技术研究院ꎬ北京１０００７４)ＣｈａｒａｃｔｅｒｉｚａｔｉｏｎｏｆＴｒａｎｓｏｎｉｃＡｅｒｏｄｙｎａｍｉｃＤａｍｐｉｎｇｏｆＲｏｃｋｅｔｓＢａｓｅｄｏｎＣＦＤ/ＣＳＤＣｏｕｐｌｉｎｇＬＩＹｏｎｇ￣ｄｅꎬ㊀ＧＵＯＬｉꎬ㊀ＪＩＣｈｅｎ(ＣｈｉｎａＡｃａｄｅｍｙｏｆＡｅｒｏｓｐａｃｅＡｅｒｏｄｙｎａｍｉｃｓꎬＢｅｉｊｉｎｇ１０００７４ꎬＣｈｉｎａ)摘㊀要:随着新型大推力火箭的发展ꎬ弯曲模态频率的不断降低ꎬ以及流动分离和跨声速飞行时产生的激波震荡等因素ꎬ其在跨声速飞行过程中更容易出现非定常振动发散ꎮ文章以某带助推的运载火箭模型为研究对象ꎬ通过数值计算获取火箭强迫振动时的气动阻尼ꎬ并对影响火箭气动阻尼的因素进行了分析ꎮ包括结构节点位置㊁振动振幅大小㊁脉动压力等ꎮ研究表明:助推主要起到增大气动阻尼的作用ꎻ前节点主要影响收缩段的气动阻尼ꎻ振动振幅大小和脉动压力对气动阻尼的影响可忽略不计ꎮ关键词:气动阻尼ꎻ数值计算ꎻ跨声速ꎻ气动弹性ꎻ运载火箭㊀㊀㊀收稿日期:２０２３￣０９￣２５ꎻ修回日期:２０２３￣１０￣２３第一作者简介:李泳德(１９９５￣)㊀男ꎬ工学硕士ꎬ助理工程师ꎬ主要研究方向为气动弹性分析ꎮＥ￣ｍａｉｌ:５６２０６４１６９＠ｑｑ.ｃｏｍ通信作者简介:季辰(１９８２￣)㊀男ꎬ工学博士ꎬ研究员ꎬ主要研究方向为气动弹性力学ꎮＥ￣ｍａｉｌ:ｊｉｃｈｅｎ１６７＠ｈｏｔｍａｉｌ.ｃｏｍ中图分类号:Ｖ４７５.１㊀㊀文献标志码:ＡＡｂｓｔｒａｃｔ:Ｗｉｔｈｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆｎｅｗｈｉｇｈ￣ｔｈｒｕｓｔｒｏｃｋｅｔｓꎬｔｈｅｄｅｃｒｅａｓｉｎｇｆｒｅｑｕｅｎｃｙｏｆｔｈｅｂｅｎｄｉｎｇｍｏｄｅｓｏｆｔｈｅｒｏｃｋｅｔꎬａｓｗｅｌｌａｓｔｈｅｆａｃｔｏｒｓｓｕｃｈａｓｆｌｏｗｓｅｐａｒａｔｉｏｎａｎｄｓｈｏｃｋｏｓｃｉｌｌａｔｉｏｎｓｇｅｎｅｒａｔｅｄｄｕｒｉｎｇｔｒａｎｓｏｎｉｃｆｌｉｇｈｔｍａｋｅｉｔｍｏｒｅｐｒｏｎｅｔｏｎｏｎ￣ｃｏｎｓｔａｎｔｖｉｂｒａｔｉｏｎ.Ｉｎｔｈｉｓｐａｐｅｒꎬａｌａｕｎｃｈｖｅｈｉｃｌｅｍｏｄｅｌｗｉｔｈｂｏｏｓｔｗａｓｔａｋｅｎａｓｔｈｅｒｅｓｅａｒｃｈｏｂｊｅｃｔꎬａｎｄｔｈｅａｅｒｏｄｙ￣ｎａｍｉｃｄａｍｐｉｎｇｏｆｔｈｅｒｏｃｋｅｔｄｕｒｉｎｇｆｏｒｃｅｄｖｉｂｒａｔｉｏｎｗａｓｏｂｔａｉｎｅｄｔｈｒｏｕｇｈｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎ.Ｔｈｅｆａｃｔｏｒｓａｆｆｅｃｔｉｎｇｔｈｅａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇｏｆｔｈｅｒｏｃｋｅｔｗｅｒｅａｎａｌｙｚｅｄꎬｉｎｃｌｕｄｉｎｇｔｈｅｐｏｓｉｔｉｏｎｏｆｓｔｒｕｃｔｕｒａｌｎｏｄｅｓꎬｔｈｅｍａｇｎｉｔｕｄｅｏｆｖｉｂｒａｔｉｏｎａｍ￣ｐｌｉｔｕｄｅꎬｐｕｌｓａｔｉｎｇｐｒｅｓｓｕｒｅａｎｄｓｏｏｎ.Ｔｈｅｓｔｕｄｙｓｈｏｗｓｔｈａｔｔｈｅｂｏｏｓｔｍａｉｎｌｙｐｌａｙｓｔｈｅｒｏｌｅｏｆｉｎｃｒｅａｓｉｎｇａｅｒｏｄｙｎａｍｉｃｄａｍｐ￣ｉｎｇａｎｄｔｈｅｆｒｏｎｔｎｏｄｅｍａｉｎｌｙａｆｆｅｃｔｓｔｈｅａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇｏｆｔｈｅｃｏｎｔｒａｃｔｉｏｎｓｅｃｔｉｏｎ.Ｔｈｅｖｉｂｒａｔｉｏｎａｍｐｌｉｔｕｄｅｓｉｚｅａｎｄｔｈｅｐｕｌｓａｔｉｎｇｐｒｅｓｓｕｒｅｈａｖｅａｎｅｇｌｉｇｉｂｌｅｅｆｆｅｃｔｏｎｔｈｅａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇ.Ｋｅｙｗｏｒｄｓ:ａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇꎻｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎꎻｔｒａｎｓｏｎｉｃꎻａｅｒｏｅｌａｓｔｉｃｉｔｙꎻｌａｕｎｃｈｖｅｈｉｃｌｅ引㊀言通常情况下人们认为气动力对火箭的振动起到阻尼作用ꎬ即气动阻尼为正值ꎮ然而随着大推力火箭发展ꎬ火箭的长细比逐渐加大ꎬ导致弯曲刚度越来越小ꎬ同时为了满足有效载荷的外形要求ꎬ火箭头部整流罩尺寸不断加大ꎬ后续箱体的直径却保持不变ꎬ形成了典型的锤头体外形ꎮ国内外大量的火箭研制经验表明[１￣９]ꎬ对于此类锤头体外形火箭的气动设计ꎬ必须要进行动态气动载荷与动态气弹稳定性分析ꎬ否则设计的疏忽可能会导致火箭结构出现毁灭性的破坏进而导致发射失败ꎮ目前常用的衡量气弹稳定性的方法是通过风洞试验来获取气动阻尼系数ꎮ早在１９６３年ꎬ美国国家航空航天局Ａｍｅｓ研究中心(ＮＡＳＡＡｍｅｓＲｅ￣ｓｅａｒｃｈＣｅｎｔｅｒ)采用半刚性模型开展试验研究[１０]ꎬ获取火箭头部的气动阻尼来评估其稳定性ꎬ但这只能用来模拟火箭弯曲振型前节点之前部分的结构动力学特性ꎮ直到兰利研究中心(ＮＡＳＡＬａｎｇｌｅｙＲｅｓｅａｒｃｈＣｅｎｔｅｒ)开发了全弹性模型气动阻尼试验气体物理２０２３年㊀第８卷技术ꎬ其可以模拟整体的结构动力学特性以及气动外形ꎬ并应用于多款运载火箭研制[１１￣１５]ꎮ国内ꎬ中国航天空气动力技术研究院对气动阻尼问题开展过较多的研究[１６￣２０]ꎬ从模型设计方法㊁模型制作工艺㊁试验机构设计和数据处理等诸多方面ꎬ逐步改进实现了从半刚性模型到全弹性模型的过渡ꎬ并在多个型号上得到验证ꎮ然而通过风洞试验研究气动弹性问题ꎬ技术难度大ꎬ试验成本高ꎬ同时几乎不可能开展全尺寸试验ꎮ因此通过数值计算的方法开展相关研究是另一种重要的手段ꎮ刘子强等[２１]实现了通过数值计算确定气动阻尼系数的技术和方法ꎬ并与试验结果进行对比ꎬ证实了该方法的可靠性ꎮ冉景洪等[２２]通过模态数据结合准定常理论的方法分析了减阻杆加后体这一弹性结构的气动阻尼ꎬ结果表明减阻杆造成的分离流会对后体的气动阻尼系数产生影响ꎮ朱剑等[２３]针对新一代捆绑式运载火箭发展了非结构网格下的气动阻尼计算方法ꎬ并分析了攻角㊁Ｍａｃｈ数等参数对气动阻尼的影响ꎮ本文在之前的计算方法[２３]的基础上采用ＩＤＤＥＳ模型ꎬ考虑脉动压力的影响ꎬ通过强迫振动的方式ꎬ针对捆绑式运载火箭的某一特定模态进行数值计算仿真ꎬ研究前节点位置ꎬ振动振幅ꎬ脉动压力等参数对气动阻尼的影响规律ꎮ１㊀计算方法图１为本文所用的捆绑式运载火箭的计算模型ꎬ是典型的锤头体结构ꎮ在跨声速阶段ꎬ其头部会产生激波造成激波边界层干扰ꎬ而在锤头体外形的过渡段会出现气流分离ꎮ为探究各部分气动阻尼的变化ꎬ将整个箭体分为头部㊁过渡段㊁弹身３个部分ꎮ图１㊀表面网格及区域划分Ｆｉｇ.１㊀Ｓｕｒｆａｃｅｇｒｉｄａｎｄｒｅｇｉｏｎｄｉｖｉｓｉｏｎ１.１㊀流场仿真模型本文分别用Ｒｅｙｎｏｌｄｓ平均法(Ｒｅｙｎｏｌｄｓ￣ａｖｅｒａｇｅｄＮａｖｉｅｒ￣ＳｔｏｋｅｓꎬＲＡＮＳ)和改进的延迟分离涡模拟(ｉｍｐｒｏｖｅｄｄｅｌａｙｅｄｄｅｔａｃｈｅｄ￣ｅｄｄｙｓｉｍｕｌａｔｉｏｎꎬＩＤ￣ＤＥＳ)[２４￣２５]进行计算ꎬ在ＲＡＮＳ方程中ꎬ将变量分为平均值和波动值两部分ꎬ对于速度分量有ｕｉ＝ｕｉ＋ｕᶄｉ其中ꎬｉ＝１ꎬ２ꎬ３ꎬｕｉ和ｕᶄｉ分别代表平均量和波动量ꎬ对于压强和其他标量也采用类似的形式ꎬ将这种形式代入连续性方程和动量方程中ꎬ并写成张量形式∂ρ∂ｔ＋∂∂ｘｉ(ρｕｉ)＝０(１)∂∂ｔ(ρｕｉ)＋∂∂ｘｊ(ρｕｉｕｊ)＝∂ｐ∂ｘｉ＋∂∂ｘｊμ∂ｕｉ∂ｘｊ＋∂ｕｊ∂ｘｉ－２３δｉｊ∂ｕｋ∂ｘｋæèçöø÷éëêêùûúú＋∂∂ｘｊ(－ρｕᶄｉｕᶄｊ)(２)其中ꎬｉꎬｊꎬｋ可分别取１ꎬ２ꎬ３ꎻρ是密度ꎻｔ是时间ꎻ当ｉ＝ｊ时δｉｊ取０ꎬ否则取１ꎮ式(１)㊁(２)是ＲＡＮＳ方程ꎬ由方程可知ＲＡＮＳ方法将湍流脉动对平均流动的作用模化为Ｒｅｙｎｏｌｄｓ应力项即－ρｕᶄｉｕᶄｊꎬ之后采用湍流模型进行封闭ꎬ本文采用的湍流模型为ＳＳＴｋ￣ω模型ꎬ其输运方程为∂∂ｔ(ρｋ)＋∂∂ｘｉ(ρｋｕｉ)＝∂∂ｘｊΓｋ∂ｋ∂ｘｊæèçöø÷＋Ｇｋ－Ｙｋ∂∂ｔ(ρω)＋∂∂ｘｉ(ρωｕｉ)＝∂∂ｘｊΓω∂ω∂ｘｊæèçöø÷＋Ｇω－Ｙω其中ꎬｋ和ω分别代表湍流动能和湍流耗散率ꎬΓｋ和Γω分别代表ｋ和ω的有效扩散系数ꎬＧｋ和Ｇω分别代表ｋ和ω的生成率ꎬＹｋ和Ｙω分别代表ｋ和ω的耗散率ꎮ因此ＲＡＮＳ方法只能计算大尺度的平均流动ꎬ本文采用ＩＤＤＥＳ方法计算脉动压力对气动阻尼的影响ꎮＩＤＤＥＳ方法是由分离涡模拟(ｄｅｔａｃｈｅｄ￣ｅｄｄｙｓｉｍｕｌａｔｉｏｎꎬＤＥＳ)方法改进而来ꎬ其本质思想与ＤＥＳ方法相同ꎬ是想以网格尺度和模型中的特征尺度隐式划分ＲＡＮＳ和大涡模拟(ｌａｒｇｅ￣ｅｄｄｙｓｉｍｕｌａｔｉｏｎꎬＬＥＳ)区域ꎬ使其既能处理ＲＡＮＳ方法无法得到的脉动场ꎬ也能降低ＬＥＳ方法在模拟高Ｒｅｙｎｏｌｄｓ数流动时所需的计算资源ꎮ区别在于当边界层较厚或者分离区域较窄时ꎬＤＥＳ方法会出现如模型应力损耗(ｍｏｄｅｌｅｄｓｔｒｅｓｓｄｅｐｌｅｔｉｏｎꎬＭＳＤ)ꎬ网格诱导分离(ｇｒｉｄ￣ｉｎｄｕｃｅｄｓｅｐａｒａｔｉｏｎꎬＧＩＳ)以及对数层不匹配(ｌｏｇａｒｉｔｈｍｉｃ￣ｌａｙｅｒｍｉｓｍａｔｃｈꎬＬＬＭ)问题[２４]ꎬ而ＩＤＤＥＳ模型通过改良计算区域划分ꎬ结合延迟分离涡模拟(ｄｅｌａｙｅｄｄｅｔａｃｈｅｄ￣ｅｄｄｙｓｉｍｕｌａｔｉｏｎꎬＤＤＥＳ)和０３第６期李泳德ꎬ等:基于ＣＦＤ/ＣＳＤ耦合的火箭跨声速气动阻尼特性分析壁面模型大涡模拟(ｗａｌｌ￣ｍｏｄｅｌｅｄｌａｒｇｅ￣ｅｄｄｙｓｉｍｕｌａ￣ｔｉｏｎꎬＷＭＬＥＳ)ꎬ定义新的长度尺度解决了这些问题ꎬ具体公式详见文献[２５]ꎮ流场网格如图２㊁图３所示ꎬ边界层采用棱柱层结构ꎬ并调整第１层网格高度使得ｙ＋小于１ꎬ远场部分采用六面体结构网格ꎬ与边界层的过渡层采用非结构网格ꎮ整体网格单元数量为４.２ˑ１０６ꎮ图２㊀ｙ方向截面网格示意图Ｆｉｇ.２㊀Ｓｃｈｅｍａｔｉｃｄｉａｇｒａｍｏｆｃｒｏｓｓ￣ｓｅｃｔｉｏｎａｌｇｒｉｄｉｎｔｈｅｙ￣ｄｉｒｅｃｔｉｏｎ图３㊀ｘ方向截面网格示意图Ｆｉｇ.３㊀Ｓｃｈｅｍａｔｉｃｄｉａｇｒａｍｏｆｃｒｏｓｓ￣ｓｅｃｔｉｏｎａｌｇｒｉｄｉｎｔｈｅｘ￣ｄｉｒｅｃｔｉｏｎ物面边界条件为无滑移壁面条件ꎬ远场采用压力远场边界条件ꎬ湍流模型采用ＳＳＴｋ￣ω模型ꎬ采用密度基求解ꎬ气体黏性采用Ｓｕｔｈｅｒｌａｎｄ定律ꎬ空间离散采用２阶迎风格式ꎬ对流通量采用Ｒｏｅ格式ꎮ１.２ 结构分析模型结构与流场耦合分析过程中ꎬ结构部分可以采用模态方法描述ꎮ结构模态可以通过有限元方法与结构模态试验方法获得ꎮ本文采用有限元分析结果获得的模态ꎬ图４所示为结构的前３阶模态ꎬ本文只分析计算结果中气动阻尼最小的第２阶模态ꎮ(ａ)ｆ＝１.２００Ｈｚ(ｂ)ｆ＝２.４６０Ｈｚ(ｃ)ｆ＝２.９５７Ｈｚ图４㊀结构的前３阶模态Ｆｉｇ.４㊀Ｆｉｒｓｔｔｈｒｅｅｍｏｄｅｓｏｆｔｈｅｓｔｒｕｃｔｕｒｅ由于火箭结构外形简单ꎬ一般不考虑其扭转影响ꎬ因此可以将其简化为简单的梁模型ꎬ这样就可以给出其模态振动方程ｑ㊆ｉ＋２ｂｉωｉｑ˙ｉ＋ω２ｉｑｉ＝ｆｉ(３)式中ꎬｑｉ为第ｉ阶模态的广义位移ꎬｂｉ为第ｉ阶模态的结构阻尼系数ꎬωｉ为第ｉ阶模态的固有频率ꎬ１３气体物理２０２３年㊀第８卷ｆｉ为第ｉ阶模态下质量归一化的广义气动力ꎮ若将ｆｉ按照Ｔａｙｌｏｒ展开并略去高阶项ꎬ可以将其转化为气动阻尼项与气动刚度项的形式ꎬ则式(３)可写为ｑ㊆ｉ＋２(ｂｉ＋Ｂｉ)ωｉｑ˙ｉ＋(Ｋｉ＋１)ω２ｉｑｉ＝０(４)式中ꎬＢｉ为气动阻尼系数ꎬＫｉ为气动刚度系数ꎬ研究表明[２６]ꎬ气动刚度相对于结构刚度为小量可以忽略不计ꎬ而在计算中结构阻尼往往设置为０ꎬ因此气动阻尼可以直接反映其气弹稳定性ꎮ１.３㊀气动阻尼分析原理气动阻尼的分析可以采用强迫振动或者自由振动的方式进行ꎬ这两种方法获得的时域数据不同ꎬ提取气动阻尼的方式也不同ꎮ强迫振动方法初始演化过程较短ꎬ因此计算量较小ꎬ同时能够分析某一种振动形式的气动阻尼ꎬ明确该振动形式是收敛还是发散ꎮ分析过程中能够获得不同部位与部件的气动阻尼ꎮ但是对于多模态相互作用引起的发散(例如颤振)较难预测ꎮ自由振动方法需要一定的自由演化时间才能够对时域数据进行分析ꎬ不过自由振动方法能够获得最能够吸收能量的模态及其振动频率ꎮ对于本研究所关注的问题ꎬ气动载荷对结构振动的过程中气动阻尼的影响较大ꎬ而对气动刚度与气动质量影响较小ꎬ即结构的固有振动频率受到来流的影响较小ꎬ其稳定性问题主要由气动阻尼的正㊁负引起ꎬ所以采用强迫振动方法分析ꎮ强迫振动下结构做简谐模态振动ｑｉ(ｔ)＝Ａｓｉｎ(ωｉｔ)式中ꎬＡ表示振动的振幅ꎬ将其代入计算气动力的公式中[２１]并做正交积分可得Ｂｉ＝ʏｌ０Ｂｘ(ｘ)ｄｘ＝－１ＭｉＡω２ｉＴʏｌ０ʏｔ０＋Ｔｔ０Ｇ(ｘꎬｔ)ｃｏｓ(ωｉｔ)ｄｔｄｘ(５)式中ꎬＭｉ为第ｉ阶模态的模态质量ꎬＴ为整数倍周期ꎬＧ为广义气动力ꎮ根据式(５)便可以得到局部或分区域的气动阻尼ꎮ１.４㊀耦合计算流程首先进行模态分析ꎬ以确定结构的模态频率与振型ꎬ用以设计强迫振动的频率和振幅ꎮ非定常流场计算前先进行定常流场计算ꎬ来加快非定常计算的演化速度并增强收敛性ꎬ结构节点位移通过径向基函数(ＲＢＦ)插值方法[２７]映射到气动网格节点上ꎬ来进行网格的变形ꎬ这里径向基函数选用ＷｅｎｄｌａｎｄＣ２ꎬ如下所示φ(ｘ)＝(１－ｘ)４(４ｘ＋１)最后将计算出来的广义力提取出来ꎬ截取演化完毕的整数倍周期ꎬ进行气动阻尼计算ꎮ耦合计算流程图如图５所示ꎮ图５㊀耦合计算流程图Ｆｉｇ.５㊀Ｆｌｏｗｃｈａｒｔｏｆｃｏｕｐｌｅｄｃａｌｃｕｌａｔｉｏｎ２㊀结果分析与讨论２.１㊀流场分析结果计算的来流Ｍａｃｈ数范围为０.７~１.２ꎮ其中中截面的压力分布如图６所示ꎮ可以看出在头部出现了膨胀波以及跨声速激波ꎬ在过渡段存在流动分离ꎬ随着Ｍａｃｈ数的增大ꎬ头部低压区域逐渐扩张ꎬ并且能明显看到ꎬ在流动再附的位置产生了再附激波ꎮ(ａ)Ｍａ＝０.７０２３第６期李泳德ꎬ等:基于ＣＦＤ/ＣＳＤ耦合的火箭跨声速气动阻尼特性分析(ｂ)Ｍａ＝０.７５(ｃ)Ｍａ＝０.８０(ｄ)Ｍａ＝０.８５(ｅ)Ｍａ＝０.８８(ｆ)Ｍａ＝０.９０(ｇ)Ｍａ＝０.９２(ｈ)Ｍａ＝０.９６(ｉ)Ｍａ＝０.９８３３气体物理２０２３年㊀第８卷(ｊ)Ｍａ＝１.００(ｋ)Ｍａ＝１.０５(ｌ)Ｍａ＝１.１０图６㊀不同Ｍａｃｈ数下的中截面压力分布Ｆｉｇ.６㊀ＰｒｅｓｓｕｒｅｄｉｓｔｒｉｂｕｔｉｏｎｉｎｔｈｅｍｉｄｄｌｅｓｅｃｔｉｏｎａｔｄｉｆｆｅｒｅｎｔＭａｃｈｎｕｍｂｅｒｓ２.２ 气动阻尼分布通过上述流场分析ꎬ可以看出火箭不同部位流动结构并不相同ꎬ在头部与箭身上ꎬ流动主要为附着流动ꎬ而在过渡段会出现较为复杂的波系结构以及流动分离ꎮ针对不同的流动结构随流向站位ｘ的变化ꎬ设该位置上广义力与广义位移的相位差为φ(ｘ)ꎬ并且简谐振动没有引入其他模态的广义力ꎬ则广义力的表达式为Ｇ(ｘꎬｔ)＝Ｆｇｅｎ ｓｉｎ[ωｔ＋φ(ｘ)]＋Ｆ０(６)其中ꎬＦｇｅｎ为广义力的振动幅度ꎬＦ０为广义力的常数偏移量ꎮ将式(６)代入到式(５)中得到Ｂ(ｘ)＝－ＦｇｅｎＭＡω２Ｔʏｔ０＋Ｔｔ０ｓｉｎ[ωｔ＋φ(ｘ)]ｃｏｓ(ωｔ)ｄｔ其中ꎬ广义力的常数偏移量Ｆ０的积分为０ꎬ因此省略ꎮ通过将等式中的正弦函数部分进行和差化积得到Ｂ(ｘ)＝－ＦｇｅｎＭＡω２Ｔʏｔ０＋Ｔｔ０ｓｉｎ(ωｔ)ｃｏｓ[φ(ｘ)]ｃｏｓ(ωｔ)ｄｔ＋[ʏｔ０＋Ｔｔ０ｓｉｎ[φ(ｘ)]ｃｏｓ(ωｔ)ｃｏｓ(ωｔ)ｄｔ](７)式(７)中第１部分在整个周期中的积分为０ꎬ只有第２部分保留ꎬ因此得到Ｂ(ｘ)＝－Ｆｇｅｎｓｉｎ[φ(ｘ)]ＭＡω２Ｔʏｔ０＋Ｔｔ０ｃｏｓ２(ωｔ)ｄｔ(８)式(８)中积分部分恒为正值ꎬ决定整个气动阻尼的部分只有相位角φ(ｘ)的正弦值ｓｉｎ[φ(ｘ)]ꎬ为了能够更加直观地获得相位角与气动阻尼Ｂ之间的关系ꎬ须将符号转化为对应的正弦函数转角ꎬ根据正弦关系ꎬ此转角为πꎬ因此得到Ｂ(ｘ)＝－Ｆｇｅｎ(ｘ)ｓｉｎ[φ(ｘ)＋π]ＭＡω２Ｔʏｔ０＋Ｔｔ０ｃｏｓ２(ωｔ)ｄｔ(９)图７为气动阻尼变化曲线ꎬ可以看出随着Ｍａｃｈ数的增大ꎬ整体气动阻尼先增大后减少ꎬ在Ｍａｃｈ数为０.９８时达到最大值ꎬ过渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Ｍａｃｈ数的增大一直增大ꎬ只是增长速率变缓ꎮ图７㊀有助推时气动阻尼变化曲线Ｆｉｇ.７㊀Ａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇｃｈａｎｇｅｃｕｒｖｅｗｉｔｈｂｏｏｓｔ根据式(９)ꎬ得到相位角与气动阻尼Ｂ之间的关系为:当φ(ｘ)ɪ(－πꎬ０)时ꎬ相位角滞后ꎬ气动阻尼Ｂ为负值ꎻ当φ(ｘ)ɪ(０ꎬπ)ꎬ相位角提前ꎬ４３第６期李泳德ꎬ等:基于ＣＦＤ/ＣＳＤ耦合的火箭跨声速气动阻尼特性分析气动阻尼Ｂ为正值ꎻ为当φ(ｘ)＝０时ꎬ无相位角差别ꎬ气动阻尼Ｂ为０ꎮ在过渡段上ꎬ复杂的波系结构以及流动分离ꎬ使得气动力与结构位移之间会出现较为明显的迟滞现象ꎬ从而导致相位角φ(ｘ)ɪ(－πꎬ０)ꎬ由此在过渡段上产生了负的气动阻尼ꎮ计算过程中的广义力与广义位移随时间变化曲线如图８所示ꎬ可以看出所有工况计算结果都表现良好ꎬ需要注意的是在非定常计算初期ꎬ演化的不完全导致广义力存在一些突变异常的结果ꎬ计算气动阻尼时须剔除ꎬ选择后面演化完全的周期ꎮ本文计算了９个周期ꎬ剔除了第１个周期出现的错误结果ꎬ采用后８个周期进行气动阻尼分析ꎮ强迫运动振幅为芯级直径的０.５％ꎮ(ａ)Ｍａ＝０.７０㊀㊀㊀(ｂ)Ｍａ＝０.７５(ｃ)Ｍａ＝０.８０㊀㊀㊀(ｄ)Ｍａ＝０.８５(ｅ)Ｍａ＝０.８８㊀㊀㊀(ｆ)Ｍａ＝０.９０５３气体物理２０２３年㊀第８卷(ｇ)Ｍａ＝０.９２㊀㊀㊀(ｈ)Ｍａ＝０.９６(ｉ)Ｍａ＝０.９８㊀㊀㊀(ｊ)Ｍａ＝１.００(ｋ)Ｍａ＝１.０５㊀㊀㊀(ｌ)Ｍａ＝１.１０图８㊀不同工况下的广义力与广义位移随时间变化曲线Ｆｉｇ.８㊀Ｔｉｍｅｄｅｐｅｎｄｅｎｔｃｕｒｖｅｓｏｆｇｅｎｅｒａｌｉｚｅｄｆｏｒｃｅａｎｄｇｅｎｅｒａｌｉｚｅｄｄｉｓｐｌａｃｅｍｅｎｔｕｎｄｅｒｄｉｆｆｅｒｅｎｔｏｐｅｒａｔｉｎｇｃｏｎｄｉｔｉｏｎｓ２.３㊀气动阻尼影响因素２.３.１㊀有无助推对气动阻尼的影响捆绑式运载火箭相比于传统的运载火箭ꎬ最大的区别就是在尾部四周捆绑了助推器ꎬ使得其流场特性变得复杂ꎬ因此须分析其对气动阻尼的影响ꎮ图７㊁图９分别为有无助推时气动阻尼变化曲线ꎬ可以看出随着Ｍａｃｈ数的增大整体气动阻尼先增大后减少ꎬ在Ｍａｃｈ数为０.９８时达到最大值ꎬ过６３第６期李泳德ꎬ等:基于ＣＦＤ/ＣＳＤ耦合的火箭跨声速气动阻尼特性分析渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Ｍａｃｈ数的增大一直增大ꎬ只是增长速率变缓ꎮ对比两个图可知ꎬ助推主要起增大气动阻尼的作用ꎮ还可以看出有无助推情况下头部的气动阻尼变化很小ꎬ意味着在箭体尾部施加控制很难影响到头部的气动阻尼ꎬ特别是在超声速流场中ꎮ图９㊀无助推时气动阻尼变化曲线Ｆｉｇ.９㊀Ａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇｃｈａｎｇｅｃｕｒｖｅｗｉｔｈｏｕｔｂｏｏｓｔ２.３.２㊀前节点位置影响为了考察前节点位置变化对气动阻尼的影响ꎬ在保持振动频率不变㊁头部最大振型位置与振幅不变的条件下移动前节点ꎬ变化后的振型如图１０所示ꎮ(ａ)Ｆｒｏｎｔｎｏｄｅａｆｔｅｒｔｈｅｔｒａｎｓｉｔｉｏｎｒｅｇｉｏｎ(ｂ)Ｆｒｏｎｔｎｏｄｅｉｎｔｈｅｔｒａｎｓｉｔｉｏｎｒｅｇｉｏｎ(ｃ)Ｆｒｏｎｔｎｏｄｅｂｅｆｏｒｅｔｈｅｔｒａｎｓｉｔｉｏｎｒｅｇｉｏｎ图１０㊀前节点变化后的振型Ｆｉｇ.１０㊀Ｖｉｂｒａｔｉｏｎｍｏｄｅａｆｔｅｒｔｈｅｃｈａｎｇｅｏｆｆｏｒｍｅｒｎｏｄｅ根据对计算结果的分析分别获得了不同前节点位置的整体气动阻尼对比与过渡段气动阻尼对比ꎬ如图１１㊁图１２所示ꎬ可以看出前节点位置的改变并没有影响整体气动阻尼随Ｍａｃｈ数增大而增大的趋势ꎬ且前节点在过渡段上与过渡段前的整体气动阻尼相差不大ꎬ而前节点在过渡段后的整体气动阻尼要高于另两种情况ꎬ因此过渡段与头部放在同一侧有助于提高气动阻尼ꎮ过渡段的气动阻尼会随着前节点的变化发生剧烈改变ꎬ前节点在过渡段前后随Ｍａｃｈ数增大的变化规律相反ꎬ节点前后的振动相位变化导致不同节点位置过渡段的振动相位不同ꎬ进而导致气动阻尼发生变化ꎮ图１１㊀不同节点位置的整体气动阻尼Ｆｉｇ.１１㊀Ｏｖｅｒａｌｌａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇａｔｄｉｆｆｅｒｅｎｔｎｏｄｅｐｏｓｉｔｉｏｎｓ图１２㊀不同节点位置的过渡段气动阻尼Ｆｉｇ.１２㊀Ａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇｏｆｔｈｅｔｒａｎｓｉｔｉｏｎｒｅｇｉｏｎａｔｄｉｆｆｅｒｅｎｔｎｏｄｅｐｏｓｉｔｉｏｎｓ２.３.３㊀强迫振动振幅大小对气动阻尼的影响为了考察强迫振动振幅大小对气动阻尼的影响ꎬ在保证流场结构不发生改变的前提下ꎬ振动振幅分别为原来的一半和两倍ꎬ根据工程经验ꎬ如果振幅超过芯级直径的５％ꎬ则须考虑流场结构改变所造成的影响ꎮ图１３㊁图１４分别为不同振幅下的整体与头部气动阻尼ꎮ７３气体物理２０２３年㊀第８卷图１３㊀不同振幅下整体气动阻尼Ｆｉｇ.１３㊀Ｏｖｅｒａｌｌａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇａｔｄｉｆｆｅｒｅｎｔａｍｐｌｉｔｕｄｅｓ图１４㊀不同振幅下头部气动阻尼Ｆｉｇ.１４㊀Ａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇｏｆｔｈｅｈｅａｄｒｅｇｉｏｎａｔｄｉｆｆｅｒｅｎｔａｍｐｌｉｔｕｄｅｓ可以发现改变振幅无论是对整体气动阻尼还是头部气动阻尼来说变化都很小ꎬ这意味着气动阻尼的大小主要取决于气动力与结构振动的相位差ꎬ不依赖于振动幅度的大小ꎮ２.３.４㊀脉动压力对气动阻尼的影响为了模拟出脉动压力的影响ꎬ采用ＩＤＤＥＳ方法对火箭气动阻尼进行计算ꎬ计算来流Ｍａｃｈ数为０.９２ꎬ计算过程中的广义力与广义位移如图１５所示ꎬ相较于图８可以看出广义力随时间变化曲线并不光滑ꎬ脉动压力的存在导致广义力由多个频率叠加而成ꎮ由于第２阶模态的频率为２.４６Ｈｚꎬ而由分离流㊁激波振荡等引起的脉动压力频率往往远大于此频率ꎬ因此这里选择３.５Ｈｚ为分界ꎬ将高于３.５Ｈｚ的部分视为由抖振脉动压力引起的广义力ꎬ低于３.５Ｈｚ的部分视为强迫振动引起的广义力ꎬ通过低通滤波把高于３.５Ｈｚ的广义力滤掉ꎬ可以获得由强迫振动引起的广义力与广义位移变化曲线ꎬ如图１６所示ꎬ通过此广义力计算的气动阻尼为２.０８ɢꎮ同样地ꎬ进行高通滤波将低于３.５Ｈｚ的广义力滤掉ꎬ可以获得由抖振脉动压力引起的气动阻尼为(２.９４ˑ１０－３)ɢꎬ由此得到脉动压力引起的气动阻尼变化为０.１４％ꎬ可以忽略不计ꎮ同时使用ＲＡＮＳ方法计算的气动阻尼为２.０７ɢꎬ与ＩＤＤＥＳ的计算结果相比误差约为(２.９４ˑ１０－３＋２.０８－２.０７)/２.０７ʈ０.４８％ꎬ这说明针对气动阻尼的模拟ꎬ抖振引起的脉动压力对气动阻尼的计算结果影响很小ꎬ起主要作用的还是广义力的变化ꎬ该变化由强迫振动引起的结构边界变化所导致ꎮ图１５㊀基于ＩＤＤＥＳ的广义力与广义位移变化曲线Ｆｉｇ.１５㊀ＶａｒｉａｔｉｏｎｃｕｒｅｓｏｆｇｅｎｅｒａｌｉｚｅｄｆｏｒｃｅａｎｄｇｅｎｅｒａｌｉｚｅｄｄｉｓｐｌａｃｅｍｅｎｔｂａｓｅｄｏｎＩＤＤＥＳ图１６㊀滤波后的广义力与广义位移变化曲线Ｆｉｇ.１６㊀Ｖａｒｉａｔｉｏｎｃｕｒｅｓｏｆｇｅｎｅｒａｌｉｚｅｄｆｏｒｃｅａｎｄｇｅｎｅｒａｌｉｚｅｄｄｉｓｐｌａｃｅｍｅｎｔｖａｒｉａｔｉｏｎｃｕｒｖｅａｆｔｅｒｆｉｌｔｅｒｉｎｇ３㊀结论本文通过数值计算方法研究了火箭的气动阻尼特性ꎮ根据流动特征分析与理论推导ꎬ发现火箭过渡段几何外形的收缩导致该区域出现复杂的分离与激波结构ꎬ从而造成了气动力相对于结构振动８３第６期李泳德ꎬ等:基于ＣＦＤ/ＣＳＤ耦合的火箭跨声速气动阻尼特性分析相位的滞后ꎬ导致了该区域为气动负阻尼ꎬ即气动不稳定性的主要来源ꎮ在此机理的基础上ꎬ分析了前节点位置㊁振动振幅㊁脉动压力等因素对气动阻尼的影响规律ꎮ可以得出以下结论:１)助推增加了正阻尼区域的面积ꎬ从而相对于没有助推的构型起到了增加气动阻尼的作用ꎮ２)前节点位置的改变对过渡段气动阻尼影响很大ꎬ节点前后的振动方向相反ꎬ导致节点在过渡段前后的气动阻尼变化规律也截然相反ꎬ将过渡段与头部区域放在节点的同一侧有助于增加气动阻尼ꎮ３)在不改变流场结构的前提下ꎬ改变振动的振幅ꎬ气动力也会产生相应幅度的变化ꎬ因此结构振幅对气动阻尼的影响可忽略不计ꎮ４)高频部分的广义力对气动阻尼的贡献很小ꎬ即结构振动引起的广义力变化对气动阻尼起主要作用ꎬ而脉动压力对计算气动阻尼影响不大ꎬ可忽略不计ꎮ参考文献(Ｒｅｆｅｒｅｎｃｅｓ)[１]㊀ＣｏｅＣＦ.Ｓｔｅａｄｙａｎｄｆｌｕｃｔｕａｔｉｎｇｐｒｅｓｓｕｒｅｓａｔｔｒａｎｓｏｎｉｃｓｐｅｅｄｓｏｎｔｗｏｓｐａｃｅ￣ｖｅｈｉｃｌｅｐａｙｌｏａｄｓｈａｐｅ[Ｒ].ＮＡＳＡＴＭＸ￣５０３ꎬ１９６１.[２]ＣｏｌｅＳＲＪｒꎬＨｅｎｎｉｎｇＴＬꎬＲａｉｎｅｙＡＧ.ＮＡＳＡｓｐａｃｅｖｅ￣ｈｉｃｌｅｄｅｓｉｇｎｃｒｉｔｅｒｉａ[Ｒ].ＮＡＳＡＳＰ￣８００１(ＲＥＶ)ꎬ１９６４. [３]ＥｒｉｃｓｓｏｎＬＥꎬＲｅｄｉｎｇＪＰ.Ａｎａｌｙｓｉｓｏｆｆｌｏｗｓｅｐａｒａｔｉｏｎｅｆｆｅｃｔｓｏｎｔｈｅｄｙｎａｍｉｃｓｏｆａｌａｒｇｅｓｐａｃｅｂｏｏｓｔｅｒ[Ｊ].Ｊｏｕｒ￣ｎａｌｏｆＳｐａｃｅｃｒａｆｔａｎｄＲｏｃｋｅｔｓꎬ１９６５ꎬ２(４):４８１￣４９０. [４]ＲｅｄｉｎｇＪＰꎬＥｒｉｃｓｓｏｎＬＥ.Ｅｆｆｅｃｔｏｆａｅｒｏｅｌａｓｔｉｃｃｏｎｓｉｄｅｒａ￣ｔｉｏｎｓｏｎｓｅａｓａｔ￣Ａｐａｙｌｏａｄｓｈｒｏｕｄｄｅｓｉｇｎ[Ｊ].ＪｏｕｒｎａｌｏｆＳｐａｃｅｃｒａｆｔａｎｄＲｏｃｋｅｔｓꎬ１９８１ꎬ１８(３):２４１￣２４７. [５]程镇煌.宇航飞行器跨音速气动弹性问题探讨[Ｊ].上海航天ꎬ１９９７(６):１６￣２１.ＣｈｅｎｇＺＨ.Ｄｉｓｃｕｓｓｉｏｎｏｎｔｒａｎｓｏｎｉｃｐｎｅｕｍａｔｉｃｅｌａｓｔｉｃｉｔｙｏｆｓｐａｃｅｃｒａｆｔ[Ｊ].ＡｅｒｏｓｐａｃｅＳｈａｎｇｈａｉꎬ１９９７(６):１６￣２１(ｉｎＣｈｉｎｅｓｅ).[６]倪嘉敏.我国运载火箭气动设计回顾[Ｃ].近代空气动力学研讨会论文集ꎬ２００５.ＮｉＪＭ.ＲｅｖｉｅｗｏｆａｅｒｏｄｙｎａｍｉｃｄｅｓｉｇｎｏｆＣｈｉｎａᶄｓｌａｕｎｃｈｖｅｈｉｃｌｅ[Ｃ].Ｐｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅｍｏｄｅｒｎａｅｒｏｄｙｎａｍｉｃｓｓｙｍｐｏｓｉｕｍꎬ２００５(ｉｎＣｈｉｎｅｓｅ).[７]臧涛成ꎬ胡焕性.大长细比弹箭弹性效应研究综述[Ｊ].弹道学报ꎬ１９９９ꎬ１１(３):８９￣９３ꎬ９６.ＺａｎｇＴＣꎬＨｕＨＸ.Ａｒｅｖｉｅｗｏｆｇｒｅａｔｓｌｅｎｄｅｒｎｅｓｓｒａｔｉｏｐｒｏｊｅｃｔｉｌｅｅｌａｓｔｉｃｅｆｆｅｃｔｒｅｓｅａｒｃｈ[Ｊ].ＪｏｕｒｎａｌｏｆＢａｌｌｉｓｔｉｃｓꎬ１９９９ꎬ１１(３):８９￣９３ꎬ９６(ｉｎＣｈｉｎｅｓｅ). [８]吴志刚ꎬ杨超.细长体弹箭的气动弹性问题与研究方法[Ｃ].第九届(２００５年)全国空气弹性学术交流会论文ꎬ２００５.ＷｕＺＧꎬＹａｎｇＣ.Ａｅｒｏｅｌａｓｔｉｃｉｔｙｐｒｏｂｌｅｍｓａｎｄｒｅｓｅａｒｃｈｍｅｔｈｏｄｓｏｆｓｌｅｎｄｅｒｂｏｄｙｍｉｓｓｉｌｅｓ[Ｃ].９ｔｈＮａｔｉｏｎａｌＳｙｍ￣ｐｏｓｉｕｍｏｎＡｅｒｏｅｌａｓｔｉｃｉｔｙꎬ２００５(ｉｎＣｈｉｎｅｓｅ). [９]张贺ꎬ黄晓鹏.弹性细长旋转弹箭运动稳定性问题的研究进展[Ｃ].第九届(２００５年)全国空气弹性学术交流会论文ꎬ２００５.ＺｈａｎｇＨꎬＨｕａｎｇＸＰ.Ｒｅｓｅａｒｃｈｐｒｏｇｒｅｓｓｏｎｔｈｅｓｔａｂｉｌｉｔｙｏｆｅｌａｓｔｉｃｓｌｅｎｄｅｒｒｏｔａｔｉｎｇｐｒｏｊｅｃｔｉｌｅｓ[Ｃ].９ｔｈＮａｔｉｏｎａｌＳｙｍｐｏｓｉｕｍｏｎＡｅｒｏｅｌａｓｔｉｃｉｔｙꎬ２００５(ｉｎＣｈｉｎｅｓｅ). [１０]ＣｏｌｅＨＡＪｒ.Ｄｙｎａｍｉｃｒｅｓｐｏｎｓｅｏｆｈａｍｍｅｒｈｅａｄｌａｕｎｃｈｖｅｈｉｃｌｅｓｔｏｔｒａｎｓｏｎｉｃｂｕｆｆｅｔｉｎｇ[Ｒ].ＮＡＳＡＴＮＤ￣１９８２ꎬ１９６３.[１１]ＢａｒｔｅｌｓＲＥꎬＷｉｅｓｅｍａｎＣＤꎬＭｉｎｅｃｋＲＥ.ＣｏｍｐｕｔａｔｉｏｎａｌａｅｒｏｅｌａｓｔｉｃａｎａｌｙｓｉｓｏｆｔｈｅＡｒｅｓｌａｕｎｃｈｖｅｈｉｃｌｅｄｕｒｉｎｇａｓ￣ｃｅｎｔ[Ｒ].ＡＩＡＡ２０１０￣４３７４ꎬ２０１０.[１２]ＡｚｅｖｅｄｏＪＬＦ.Ａｅｒｏｅｌａｓｔｉｃａｎａｌｙｓｉｓｏｆｌａｕｎｃｈｖｅｈｉｃｌｅｓｉｎｔｒａｎｓｏｎｉｃｆｌｉｇｈｔ[Ｊ].ＪｏｕｒｎａｌｏｆＳｐａｃｅｃｒａｆｔａｎｄＲｏｃｋｅｔｓꎬ１９８９ꎬ２６(１):１４￣２３.[１３]ＳｉｎｃｌａｉｒＡꎬＦｌｏｗｅｒｓＧ.Ｌｏｗ￣ｏｒｄｅｒａｅｒｏｅｌａｓｔｉｃｍｏｄｅｌｏｆｌａｕｎｃｈ￣ｖｅｈｉｃｌｅｄｙｎａｍｉｃｓ[Ｒ].ＡＩＡＡ２０１０￣７７２５ꎬ２０１０. [１４]ＤｏｔｓｏｎＫＷ.Ｔｒａｎｓｉｅｎｔｃｏｕｐｌｉｎｇｏｆｌａｕｎｃｈｖｅｈｉｃｌｅｂｅｎｄｉｎｇｒｅｓｐｏｎｓｅｓｗｉｔｈａｅｒｏｄｙｎａｍｉｃｆｌｏｗｓｔａｔｅｖａｒｉａｔｉｏｎｓ[Ｊ].ＪｏｕｒｎａｌｏｆＳｐａｃｅｃｒａｆｔａｎｄＲｏｃｋｅｔｓꎬ２００１ꎬ３８(１):９７￣１０４.[１５]ＣｏｌｅＳＲꎬＨｅｎｎｉｎｇＴＬ.Ｂｕｆｆｅｔｒｅｓｐｏｎｓｅｏｆａｈａｍｍｅｒｈｅａｄｌａｕｎｃｈｖｅｈｉｃｌｅｗｉｎｄ￣ｔｕｎｎｅｌｍｏｄｅｌ[Ｊ].ＪｏｕｒｎａｌｏｆＳｐａｃｅｃｒａｆｔａｎｄＲｏｃｋｅｔｓꎬ１９９２ꎬ２９(３):３７９￣３８５.[１６]崔尔杰.流固耦合力学研究与应用进展[Ｃ].钱学森科学贡献暨学术思想研讨会论文集ꎬ２００１.ＣｕｉＥＪ.Ｒｅｓｅａｒｃｈａｎｄａｐｐｌｉｃａｔｉｏｎｐｒｏｇｒｅｓｓｏｆｆｌｕｉｄ￣ｓｔｒｕｃ￣ｔｕｒｅｉｎｔｅｒａｃｔｉｏｎｍｅｃｈａｎｉｃｓ[Ｃ].ＰｒｏｃｅｅｄｉｎｇｓｏｆＳｅｍｉｎａｒｏｆＱｉａｎＸｕｅｓｅｎＳｃｉｅｎｔｉｆｉｃＣｏｎｔｒｉｂｕｔｉｏｎｓａｎｄＡｃａｄｅｍｉｃＴｈｏｕｇｈｔｓꎬ２００１(ｉｎＣｈｉｎｅｓｅ).[１７]冯明溪ꎬ王志安.火箭跨音速动导数和抖振实验[Ｊ].宇航学报ꎬ１９８７(１):５５￣６３.ＦｅｎｇＭＸꎬＷａｎｇＺＡ.Ｅｘｐｅｒｉｍｅｎｔｓｏｆｔｒａｎｓｏｎｉｃｄｅｒｉｖａ￣ｔｉｖｅｓａｎｄｂｕｆｆｅｔｉｎｇｏｆｒｏｃｋｅｔ[Ｊ].ＪｏｕｒｎａｌｏｆＡｓｔｒｏｎａｕｔｉｃｓꎬ１９８７(１):５５￣６３(ｉｎＣｈｉｎｅｓｅ).[１８]白葵ꎬ冯明溪.弹性模型实验技术[Ｊ].流体力学实验与测量ꎬ１９９９ꎬ１３(１):３８￣４２.ＢａｉＫꎬＦｅｎｇＭＸ.Ａｅｒｏｅｌａｓｔｉｃｍｏｄｅｌａｎｄｔｈｅｂｕｆｆｅｔｅｘ￣ｐｅｒｉｍｅｎｔａｌｔｅｃｈｎｉｑｕｅ[Ｊ].ＥｘｐｅｒｉｍｅｎｔｓａｎｄＭｅａｓｕｒｅｍｅｎｔｓｉｎＦｌｕｉｄＭｅｃｈａｎｉｃｓꎬ１９９９ꎬ１３(１):３８￣４２(ｉｎＣｈｉｎｅｓｅ). [１９]ＪｉＣꎬＲａｎＪＨꎬＬｉＦꎬｅｔａｌ.Ｔｈｅａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇ９３气体物理２０２３年㊀第８卷ｔｅｓｔｏｆｅｌａｓｔｉｃｌａｕｎｃｈｖｅｈｉｃｌｅｍｏｄｅｌｉｎｔｒａｎｓｏｎｉｃｆｌｏｗ[Ｃ].Ｐｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅ６４ｔｈＩｎｔｅｒｎａｔｉｏｎａｌＡｓｔｒｏｎａｕｔｉｃａｌＣｏｎ￣ｇｒｅｓｓꎬ２０１３.[２０]季辰ꎬ吴彦森ꎬ何岗ꎬ等.运载火箭气动阻尼风洞试验研究[Ｃ].第十二届全国空气弹性学术交流会论文集ꎬ２０１１.ＪｉＣꎬＷｕＹＳꎬＨｅＧꎬｅｔａｌ.Ｅｘｐｅｒｉｍｅｎｔａｌｓｔｕｄｙｏｎａｅｒｏ￣ｄｙｎａｍｉｃｄａｍｐｉｎｇｗｉｎｄｔｕｎｎｅｌｏｆｌａｕｎｃｈｖｅｈｉｃｌｅ[Ｃ].Ｐｒｏ￣ｃｅｅｄｉｎｇｓｏｆｔｈｅ１２ｔｈＮａｔｉｏｎａｌＡｅｒｏｅｌａｓｔｉｃｉｔｙＣｏｎｆｅｒｅｎｃｅꎬ２０１１(ｉｎＣｈｉｎｅｓｅ).[２１]刘子强ꎬ白葵ꎬ毛国良ꎬ等.锤头体弹性振动跨音速气动阻尼系数的确定[Ｊ].宇航学报ꎬ２００２ꎬ２３(６):１￣７.ＬｉｕＺＱꎬＢａｉＫꎬＭａｏＧＬꎬｅｔａｌ.Ｔｈｅｄｅｔｅｒｍｉｎａｔｉｏｎｏｆａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇｏｎｈａｍｍｅｒｈｅａｄｌａｕｎｃｈｖｅｈｉｃｌｅｓａｔｔｒａｎｓｏｎｉｃｓｐｅｅｄｓ[Ｊ].ＪｏｕｒｎａｌｏｆＡｓｔｒｏｎａｕｔｉｃｓꎬ２００２ꎬ２３(６):１￣７(ｉｎＣｈｉｎｅｓｅ).[２２]冉景洪ꎬ刘子强ꎬ胡静ꎬ等.减阻杆气动阻尼研究[Ｊ].力学学报ꎬ２０１４ꎬ４６(４):６３６￣６４１.ＲａｎＪＨꎬＬｉｕＺＱꎬＨｕＪꎬｅｔａｌ.Ｒｅｓｅａｒｃｈｏｆａｅｒｏ￣ｄａｍｐｉｎｇｆｏｒｂｌｕｎｔｗｉｔｈｓｐｉｋｅ[Ｊ].ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＴｈｅ￣ｏｒｅｔｉｃａｌａｎｄＡｐｐｌｉｅｄＭｅｃｈａｎｉｃｓꎬ２０１４ꎬ４６(４):６３６￣６４１(ｉｎＣｈｉｎｅｓｅ).[２３]朱剑ꎬ冉景洪ꎬ吴彦森ꎬ等.捆绑式运载火箭的气动阻尼数值计算方法[Ｃ].第十三届全国空气弹性学术交流会论文集.哈尔滨:中国力学学会ꎬ中国空气动力学会ꎬ２０１３.ＺｈｕＪꎬＲａｎＪＨꎬＷｕＹＳꎬｅｔａｌ.Ｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｍｅｔｈｏｄｆｏｒａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇｏｆｂｕｎｄｌｅｌａｕｎｃｈｖｅｈｉｃｌｅｓ[Ｃ].Ｐｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅ１３ｔｈＮａｔｉｏｎａｌＡｅｒｏｅｌａｓｔｉｃｉｔｙＣｏｎ￣ｆｅｒｅｎｃｅ.Ｈａｒｂｉｎ:ＣｈｉｎｅｓｅＳｏｃｉｅｔｙｏｆＴｈｅｏｒｅｔｉｃａｌａｎｄＡｐ￣ｐｌｉｅｄＭｅｃｈａｎｉｃｓꎬＣｈｉｎａＡｅｒｏｍｅｃｈａｎｉｃｓＳｏｃｉｅｔｙꎬ２０１３(ｉｎＣｈｉｎｅｓｅ).[２４]ＳｐａｌａｒｔＰＲ.Ｄｅｔａｃｈｅｄ￣ｅｄｄｙｓｉｍｕｌａｔｉｏｎ[Ｊ].ＡｎｎｕａｌＲｅｖｉｅｗｏｆＦｌｕｉｄＭｅｃｈａｎｉｃｓꎬ２００９ꎬ４１:１８１￣２０２. [２５]ＧｒｉｔｓｋｅｖｉｃｈＭＳꎬＧａｒｂａｒｕｋＡＶꎬＳｃｈüｔｚｅＪꎬｅｔａｌ.Ｄｅｖｅｌ￣ｏｐｍｅｎｔｏｆＤＤＥＳａｎｄＩＤＤＥＳｆｏｒｍｕｌａｔｉｏｎｓｆｏｒｔｈｅｋ￣ωｓｈｅａｒｓｔｒｅｓｓｔｒａｎｓｐｏｒｔｍｏｄｅｌ[Ｊ].ＦｌｏｗꎬＴｕｒｂｕｌｅｎｃｅａｎｄＣｏｍｂｕｓｔｉｏｎꎬ２０１２ꎬ８８(３):４３１￣４４９.[２６]季辰ꎬ吴彦森ꎬ侯英昱ꎬ等.捆绑式运载火箭跨声速气动阻尼特性试验研究[Ｊ].实验流体力学ꎬ２０２０ꎬ３４(６):２４￣３１.ＪｉＣꎬＷｕＹＳꎬＨｏｕＹＹꎬｅｔａｌ.Ｅｘｐｅｒｉｍｅｎｔａｌｓｔｕｄｙｏｆａｅｒｏｄｙｎａｍｉｃｄａｍｐｉｎｇｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆａｌａｕｎｃｈｖｅｈｉｃｌｅｗｉｔｈｂｏｏｓｔｅｒｓｉｎｔｒａｎｓｏｎｉｃｆｌｏｗ[Ｊ].ＪｏｕｒｎａｌｏｆＥｘｐｅｒｉｍｅｎｔｓｉｎＦｌｕｉｄＭｅｃｈａｎｉｃｓꎬ２０２０ꎬ３４(６):２４￣３１(ｉｎＣｈｉｎｅｓｅ). [２７]ＡｌｌｅｎＣꎬＲｅｎｄａｌｌＴＣＳ.ＵｎｉｆｉｅｄａｐｐｒｏａｃｈｔｏＣＦＤ￣ＣＳＤｉｎｔｅｒｐｏｌａｔｉｏｎａｎｄｍｅｓｈｍｏｔｉｏｎｕｓｉｎｇｒａｄｉａｌｂａｓｉｓｆｕｎｃｔｉｏｎｓ[Ｒ].ＡＩＡＡ２００７￣３８０４ꎬ２００７.０４。
各种常见阀门的原理、分类及区别（蝶阀、闸阀、截止阀、球（指南）
各种常见阀门的原理、分类及区别在生产过程中，为了使介质的压力、流量等参数符合工艺流程的要求，需要安装调节机构对上述参数进行调节。
调节机构的主要工作原理，是靠改变阀门阀瓣与阀瓣与阀座间的流通面积，达到调节上述参数的目的。
属于这类阀门的统称为控制阀一、闸阀、截止阀、蝶阀、球阀各自简述1、闸阀是指关闭件（闸板）沿通道轴线垂直移动的阀门，在管路上主要作为切断介质用，即全开或全关使用，一般不用做调节流量使用，它可以适用低温低压也可以适用高温高压，并可根据阀门的不同材质适用不同的介质。
在各类阀中，闸阀是应用最广泛的一种。
闸阀通常使用于不需要经常启闭，而且保持闸板全开或全闭的工况。
不适用于作为调节或节流使用。
全开时水流直线通过，水流阻力小，宜在管径大于50mm 的管道上采用，但水中若有杂质落入筏座易产生磨损和漏水。
2、截止阀是指关闭件（阀瓣）沿阀座中心线移动的阀门，由于该阀的阀杆开启或关闭行程相对较短，而且具有非常可靠的切断功能，又由于阀座通口的变化与阀瓣行程成正比，也非常适用于对流量的调节。
截止阀是化工生产中使用最广的一种截断类阀门，是利用阀杆升降带动与之相连的圆形阀盘（阀头），改变阀盘与阀座间距离达到控制阀门的起闭。
流线式截止阀美标式截止阀特点：截止阀上部有手轮、阀杆，中部有螺纹和填料涵密封段，小型阀门阀。
截止阀主要用于开或者关，但是备有特制形式的柱塞或特殊的套环，也可用于调节流量。
其关闭严密，但水流阻力较大，因局部阻力系数与管径成正比，故只适用于管径小于等于50mm的管道上。
闸阀和截止阀是阀门的不同型号，其阀芯结构不同，但闸阀不属于截止阀。
两种阀门的区别在于启闭件相对于阀座的移动方向不同。
闸阀的启闭件为闸板，由阀杆带动阀板做升降运动的阀门称为闸阀。
截止阀启闭件为阀瓣，由阀杆带动，沿阀座轴线做升降运动3、蝶阀是用圆盘式启闭件往复回转90°左右来开启、关闭和调节流体通道的一种阀门。
是利用一可饶轴旋转的圆盘来控制管路的启闭，转角大小反映了阀门的开启程度。
蝶阀结构强度的力学特性分析
蝶阀结构强度的力学特性分析
孙浩伟
【期刊名称】《阀门》
【年(卷),期】2024()4
【摘要】蝶阀作为一种流体控制元件,其在石油化工、燃气输送等领域的应用极其广泛。
本文利用Workbench软件平台对蝶阀进行了热-结构耦合分析,这一过程结合了阀体温度场和应力场的精确计算。
研究从两个关键层面入手,首先是对阀体内部温度场进行模拟与分析,其次是在试验条件下以及额定工作压力状态下的蝶阀进行力学特性分析。
研究发现,试验压力与额定压力工况下应力值有较大变化,但都符合设计要求,阀门关闭状态下应力值稍大,这说明当蝶板处于全关位置时,阀体一侧受到的承压载荷是相对显著的,阀门总体的变形量较小。
【总页数】5页(P471-475)
【作者】孙浩伟
【作者单位】海装沈阳局
【正文语种】中文
【中图分类】TH134
【相关文献】
1.中高压法兰蝶阀阀体结构强度的有限元分析
2.轴向柱塞双联泵中间体动力学特性及结构强度校核分析
3.极限工况下三偏心蝶阀的力学特性分析
4.基于Workbench
的硬密封蝶阀结构强度与温度场分析5.基于Workbench的蝶阀结构强度与温度场分析
因版权原因，仅展示原文概要，查看原文内容请购买。
T947H调节阀技术说明
T947H调节阀技术说明概述T947H调节阀是用圆盘式启闭件往复回转90度左右来开启、关闭和调节流体通道的一种阀门。
T947H调节阀不仅结构简单、体积小、重量轻、材料耗用省，安装尺寸小，而且驱动力矩小，操作简便、迅速，并且还可同时具有良好的流量凋节功能和关闭密封特性，是近十几年来发展最快的阀门品种之一。
随着三偏心金属密封蝶阀的问世，蝶阀开始向高温、高压、大口径、高密封性、长寿命、优良的调节特性以及一阀多功能方向发展，其可靠性及其他性能指标均达较高水平，并已部分取代截止阀、闸阀和球阀。
我厂T947H调节阀采用三偏心多层次密封结构，吸收国外同类产品结构和特点，克服了原有产品在使用中存在的缺陷，使其寿命更耐久，密封更可靠。
1. T947H调节阀的设计说明蝶阀采用三维偏心设计原理，执行机构强制密封。
密封副由两个重合的圆锥体斜截面组成，蝶板旋转中心经过对阀体、阀座的偏心值后，实现开启和关闭的瞬间接触和分离，通过调节型电动装置，可实现近似等百分比调节功能,并由蝶板密封圈金属与石墨层叠结构实现迷宫式密封，即能使阀门开关扭矩小又能保证蝶阀密封性能。
1.1. 产品设计标准和所涉及的相关标准所涉及的相关标准JB/T8527 《金属硬密封蝶阀》GB/T12221 《法兰连接金属阀门结构长度》GB/T12230 《通用阀门不锈钢铸件技术条件》GB/T1220 《不锈钢棒》GB/T13927 《通用阀门压力试验》GB/T1993 《整体钢制管法兰》HG 20596 《整体钢制管法兰》GB/T12224 《钢制阀门一般要求》GB/T5300 《通用阀门材料》GB/T3766 《液压系统通用技术条件》GB/T12220 《通用阀门标志》GB/T7928 《通用阀门供货要求》JB/T7748 《阀门清洁度和测定方法》GB/T12231 《阀门铸钢件外观质量要求》GB/T12252 《通用阀门供货要求》GB/T228 《金属拉伸试验》GB/T230 《金属洛氏硬度的试验方法》GB/T232 《金属弯曲试验》JB/ZQ4286 《包装通用技术条件》1.2．T947H调节阀的特点由三偏心原理，关阀时嵌装于阀板上的多层次密封圈与锥形阀体密封面相互压紧并成正圆接触，开启时密封圈和阀座之间的整个密封面上同时脱开，即瞬间接触和瞬间脱离。
三偏心金属密封蝶阀的力学特性分析及结构优化研究的开题报告
三偏心金属密封蝶阀的力学特性分析及结构优化研究的开题报告一、选题背景随着现代工业的快速发展，阀门已成为工业生产领域中不可缺少的核心件。
在各类阀门中，蝶阀作为一种简单、经济、可靠的流体调节设备，受到广泛应用。
而在液气管路中，蝶阀的密封性能尤为重要。
传统的金属密封蝶阀在使用中容易出现泄漏、耐腐蚀性差等问题，为了解决这些问题，三偏心金属密封蝶阀在近年来得到了越来越多的应用与关注。
其采用三个偏心结构，使得阀座与阀片之间的密封作用得到了更强的保障。
为了进一步提高这种阀门的性能，有必要对其力学特性进行分析研究，并探索其结构的优化方案。
二、研究内容与目标本次研究的主要内容为：对三偏心金属密封蝶阀的力学特性进行分析，探究其内部应力分布、变形情况等，并结合实验验证，验证分析结果的准确性。
同时，通过对现有结构的优化，进一步提高阀门的性能，减少泄漏等问题的发生率。
本研究的目标为：建立三偏心金属密封蝶阀的力学测试系统，并通过实验验证和数值分析，研究不同条件下阀门的受力变形情况；运用有限元分析的方法，对阀门结构进行优化设计，提高其在不同介质、不同压力下的密封性和耐腐蚀能力。
三、研究方法和步骤研究方法：1.建立三偏心金属密封蝶阀的力学测试系统。
利用不同的荷载测试方法，对阀门进行力学测试，并测量阀门在不同介质、不同压力下的变形情况。
2.通过实验验证和数值分析，分析阀门的内部应力分布、变形情况等，并探究不同荷载条件下阀门的受力情况。
3.运用有限元分析方法，对三偏心金属密封蝶阀的结构进行优化设计，提高阀门的密封性和耐腐蚀能力。
研究步骤：1.研究现有文献，了解三偏心金属密封蝶阀的结构和特点，确定研究的重点。
2.建立三偏心金属密封蝶阀的力学测试系统，进行力学测试，并测量阀门在不同介质、不同压力下的变形情况。
3.通过数值分析方法，对阀门的内部应力分布、变形情况等进行分析，探究不同荷载条件下阀门的受力情况。
4.利用有限元分析方法，对三偏心金属密封蝶阀的结构进行优化设计，提高阀门的密封性和耐腐蚀能力。
浅谈恰甫其海水电站大口径、双平板蝶阀的特点
第２ ９卷 第 ５期
ｌ ４
水 电 站 机 电 技 术
Ｍｅｈｎ ａ＆ Ｅｅｔｃｌ ｅｈｉｕ ｏ ｙｒｐｗ ｒ ｔｔｎ ｃａｉ ｌ ． ｌ ｒａＴ ｃｎ ｅ ｆ ｄｏｏ ｅ Ｓａｉ ｃ ｃｉ ｑ Ｈ ｏ
Ｖ ｏ１ ９ ．２ ＮＯ． ５
在 ２ ０ｍ 以 下 选 用蝶 阀 。设 计 要 求 无 需 排 空 压 力 ０ 钢 管或拆 卸 阀体便 可 以检修 ， 换 工 作 密 封 和轴 颈 更 密 封 ， 去选 用 的的蝶 阀是单板 ， 过 一旦 发 生蝶 阀密封
不好 ， 更换 密封 时 ， 就要关 闭发 电洞前 的闸门 。 由于 电站一管 分两 岔 的引 水 方式 ， 闭发 电洞前 的 闸 门 关 势 必影 响另一 台机 组 发 电。为 了解 决此 问题 ， 据 根
恰 甫 其 海水 电站 为混 合 引水 式 地 面 电站 ， 机 装 规模 ４ ＯＭｗ 。主 坝 为粘土 心墙 坝 ， 长 ３ ２ｍ， ×８ 坝 ６
（ ） 输 困难 、 路长 ： 于阀 门尺寸过 大 ， ３运 线 由 进水 阀 阀板 要求 整体 到 工 地 ， 路 拒 绝运 输 。从 制 造 厂 铁
（） １ 阀体 ： 中型 阀体 多采用 铸钢 或钢板 焊接 结 大 构 ， 电站 阀体 亦 采用 钢 板 焊 接 。进水 阀公 称 直 径 本 为 ５０ ０ｍｍ， 虑 到运 输 尺 寸 的 限制 ， ０ 考 阀体 沿 阀轴 中心垂 直方 向分成 两瓣 制造 ， 用法 兰连 接 ， 有足 够 具 的强度 和刚度 ， 任 何情 况 下 都 不 会产 生漏 水 和 有 在 害变形 。在 阀体外 部 设 有 双 筋 板 复式 加 强 梁 结构 ， 将 最危 险工况 时 阀体变形 控 制在许 用 范围 内 。阀体 轴 承座采 用箱 式结 构 ， 承 受 最 危 险工 况 时 阀 轴传 能 来 的集 中荷载 。在 阀 门进水 和 出水 段 ， 采 用 机加 均 工过 渡斜 面 ， 流 阻 系 数 比 传 统 结 构 减 少 ７ 以 使 Ｏ
闸阀,截止阀,蝶阀的比较
蝶阀又叫翻板阀,是一种结构简单的调节阀，同时也可用于低压管道介质的 开关控制。蝶阀[1](英文：butterfly valve)是指关闭件(阀瓣或蝶板)为圆盘，围 绕阀轴旋转来达到开启与关闭的一种阀，在管道上主要起切断和节流用。蝶阀启闭 件是一个圆盘形的蝶板，在阀体内绕其自身的轴线旋转，从而达到启闭或调节的目 的。蝶阀全开到全关通常是小于90° ，蝶阀和蝶杆本身没有自锁能力，为了蝶板 的定位，要在阀杆上加装蜗轮减速器。采用蜗轮减速器，不仅可以使蝶板具有自锁 能力，使蝶板停止在任意位置上，还能改善阀门的操作性能。
2、缺点 1、流体阻力大，开启和关闭时所需力较大。 2、不适用于带颗粒、粘度较大、易结焦的介质。 3、调节性能较差。 4、截止阀的种类按阀杆螺纹的位置分有外螺纹式、内螺纹式。按介质
的流向分，有直通式、直流式和角式。截止阀按密封形式分，有填料密封 截止阀和波纹管密封截止阀。
5、波纹管密封截止阀采用波纹管密封的设计，完全消除了普通阀门阀 杆填料密封老化快易泄露的缺点，不但提高了使用能源效率，增加生产设
球阀 球阀是由旋塞阀演变而来。它具有相同的旋转90度提动作，不同的是旋塞体是
球体，有圆形通孔或通道通过其轴线。球面和通道口的比例应该是这样的，即当球 旋转90度时，在进、出口处应全部呈现球面，从而截断流动。
球阀只需要用旋转90度的操作和很小的转动力矩就能关闭严密。完全平等的阀 体内腔为介质提供了阻力很小、直通的流道。通常认为球阀最适宜直接做开闭使用 ，但近来的发展已将球阀设计成使它具有节流和控制流量之用。球阀的主要特点是 本身结构紧凑，易于操作和维修，适用于水、溶剂、酸和天然气等一般工作介质， 而且还适用于工作条件恶劣的介质，如氧气、过氧化氢、甲烷和乙烯等。球阀阀体 可以是整体的，也可以是组合式的。 截止阀
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1
前言
大口径蝶阀作为埋地管线的截流、 调节装置 被广泛应用于化工、 石化、 冶炼及给水排水等行业。 对埋地管线大口径蝶阀而言， 其强度、 刚度及密封 等性能需满足内部介质、 复杂环境， 如地基沉降、 ［1， 2］ 地震等附加载荷的特殊要求 。API 6D 规定， 当管道因弯曲载荷破坏时， 要求阀门具有足够的 结构强度、 密封性和动作性能［3］。
DN2800 埋地大口径蝶阀进行刚度、 密封性能分析； 建立蝶阀抗弯模型， 对 DN2800 蝶阀及系列产品最大弯曲载荷进且可承受最大弯曲载荷为 M2。研究方法和结论对埋地大口径蝶阀承受最大弯曲载荷的研究提 关键词： 埋地大口径蝶阀； 刚度； 密封； 临界弯曲载荷； 最大弯曲载荷 中图分类号： TH134 文献标志码： A
doi： 10.3969/j.issn.1005-0329.2017.10.007
The Research of Buried Large-diameter Butterfly Valve Maximum Bending Load Based on Stiffness and Sealing Performance ZHANG Ji-wei， LEI Yan， WANG Jian， HUANG Jian， JIN Shu-jun， PENG Lin， FAN Yi-lin
收稿日期： 2016-09-23 修稿日期： 2017-01-19
针 对 某 阀 门 厂 DN2800 埋 地 大 口 径 蝶 阀 在 服 役过程中发生泄漏案例， 本文利用有限元方法对 其刚度、 密封性进行分析研究， 同时对处于复杂环 境中埋地大口径蝶阀可承受弯曲载荷与刚度、 密 封性能之间的关系进行研究， 最终得到同时满足 大口径蝶阀刚度、 密封性能的可承受最大弯曲载 荷， 为大口径蝶阀实际可承受最大弯曲载荷性能 的研究提供一种参考依据。
2017 年第 45 卷第 10 期
文章编号： 1005-0329（2017）10-0037-05
流 体 机 械
37
埋地大口径蝶阀基于刚度、密封性能的 最大弯曲载荷研究
张继伟， 雷 艳， 王 剑， 黄 健， 靳淑军， 彭 林， 范宜霖
（合肥通用机械研究院， 安徽 合肥 230000） 摘 要： 基于 API 6D 对管线阀门因弯曲载荷破坏时应具有足够结构强度、 刚度及密封性的高要求， 采用有限元法对 结果表明： （1） 长期承受附加载荷的 DN2800 蝶阀密封失效先于刚度失效， 所能承受的最大弯曲载荷为密封有效临界弯曲 载荷； （2） DN ＞ 1000 mm 系列大口径蝶阀壁厚 t、 公称直径 d 与抗弯刚度 K、 刚度有效临界弯曲载荷 M1、 密封有效临界弯曲 供了参考依据。
（Hefei General Machinery Research Institute， Hefei 230031， China） Abstract： Based on API 6D pipeline valve needs enough structural strength， stiffness and effective sealing.The stiffness and effective sealing performance of DN2800 buried large-diameter butterfly body was analyzed by using the finit element method in this paper； established a flexural model of butterfly， the maximum bending load of DN2800 butterfly and series products were thickness t， nominal diameter d， flexural stiffness K， effective critical bending stiffness M1， effective critical sealing bending load load of the buried large-diameter butterfly study. Key words： buried large-diameter butterfly； stiffness； sealing； critical bending load； maximum bending load load， maximum bending load equal to the effective critical bending load； （2）it shows quadratic function relationship among studied.The results showed： （1） the sealing failure of DN2800 butterfly prior to the stiffness failure under long-term additional M2 in DN>1000 mm large-diameter butterfly.Research methods and results provide an effective evidence to maximum bending
38
FLUID MACHINERY
Vol. 45， No.10， 2017
2 2.1
分析参数及工况的确定
MPa 时， 弹性变形结束进入塑性变形阶段。
模型的确定 考虑严苛工况， 在土体作用力、 内压和外部 弯曲载荷作用下， 以有力矩理论和弯曲应力理论 5］ 为依据建立管 - 阀承受弯矩模型［4， 。以 PN10、 DN2800 大口径蝶阀为例， 将均布载荷转化为等 值集中载荷作用在以蝶阀为中心的管道上， 左右 管长均为 8000 mm， 对称集中载荷离两侧约束距 离相等， 集中载荷间管道与阀承受的弯矩值相同， 通过不断优化集中载荷的大小， 确定蝶阀刚度和 密封同时有效的最大弯矩值。简化模型如图 1 所 示。