不等壁厚U形波纹管的轴向刚度研究

2020-06-12占丰朝张小文闫廷来

化工设备与管道 2020年2期

占丰朝，张小文，闫廷来

（1. 洛阳双瑞特种装备有限公司，河南洛阳 471001；2. 洛阳船舶材料研究所，河南洛阳 471023）

波纹管是一种具有补偿位移作用的弹性元件，广泛应用于化工、供热、航天、电力等领域。按波纹管管壁的层数，可分为单层波纹管和多层波纹管。在总厚度相同情况下，多层波纹管相对于单层波纹管具有刚度小、补偿量大、疲劳寿命长等特点，得到了广泛的应用。随着使用工况的严苛，为提高波纹管的耐腐蚀、耐压与高疲劳寿命等综合性能，设计出一种在波纹管接触环境介质侧衬一层高耐蚀合金（如Incoloy 825、Incoloy 800 等）的波纹管，并已在热力、化工产品上有所应用。EJMA—2015《膨胀节制造商协会标准》第9 章中对于多层波纹管的应用时指出为了适应管道内外不同的环境，在设计波纹管时允许使用不同的材料制作波纹管的内层和外层，用一层来防腐蚀，其他各层采用价格较为低廉的材料制作；GB/T 12777—2019《金属波纹管膨胀节通用技术条件》[1]的5.4.3 条给出了不同材料组合的多层波纹管设计温度下的许用应力计算公式。但是针对不等壁厚波纹管的应力计算、刚度计算及性能分析，较少有相关研究。

本文主要通过试验验证和有限元分析来对衬层设计方案的不等壁厚U 形波纹管的轴向刚度进行研究，以期更好地掌握不等壁厚U 形波纹管的轴向刚度，为不等壁厚波纹管的设计开发奠定基础。

1 刚度设计公式

GB/T 12777—2019 中对于无加强U 形波纹管的单波轴向刚度计算公式如式（1）所示。

式中各符号意义同GB/T 12777—2019[1]。其中n是波纹管的层数，此公式主要适用于各层壁厚相等的多层波纹管的刚度计算，当各层壁厚不相等时，基于并联弹簧刚度叠加原理，并参考等壁厚波纹管刚度计算公式得到不等壁厚波纹管轴向刚度计算公式如式（2）所示：

采用式（2）计算得到的不等壁厚波纹管的刚度数值是否能较好评估不等壁厚波纹管的轴向刚度性能，下文通过有限元分析和刚度试验做进一步验证。

2 有限元分析

2.1 有限元模型和边界条件

针对等壁厚U 形波纹管的轴向刚度性能已有多位学者进行过研究，研究结果表明采用有限元分析方法来研究波纹管的刚度性能是可行的，下文通过有限元分析软件ANSYS Workbench 对不等壁厚波纹管的轴向刚度进行分析[2]。

的两种不等壁厚方案的波纹管设计参数如表1所示。考虑到后期试验验证成本，这里所衬层波纹管的材质选择跟基层波纹管的材质一样，为耐点蚀性能较好的316L。两种方案不同之处仅在于所衬层的波纹管厚度不同，件号1 试验件所衬的一层波纹管厚度为0.5 mm，件号2 试验件所衬的一层波纹管厚度为1.0 mm，基层波纹管厚度均为2.0 mm，两个试验件的波高、波距、波根直径等其他参数保持一致。

表1 U 形波纹管设计参数Table 1 Wave form parameters of unreinforced U-shaped bellows

采用Solidworks 软件建立波纹管的三维模型，由于波纹管是轴对称结构，为了减少计算量提高计算效率，这里只建立波纹管的1/4 模型，然后将模型导入到ANSYS Workbench 中进行分析。考虑到六面体网格可以减少单元数量，加快求解收敛，这里通过扫掠方法生成六面体网格，划分网格后的模型如图1 所示。通过查看单元详细信息可以发现模型采用的单元类型对应的名字为Mesh200，这是一个特殊单元类型，并不参与实际的计算，由其对应的Abaqus 单元类型名字C3D20，可知相应的有限元单元为Solid 186 单元，是一个高阶3 维20 节点固体结构单元，该单元支持大变形、应力强化、塑性等非线性参数，能够较好模拟非线性行为，因此可以作为非线性分析的有限元单元类型[3]。

图1 波纹管有限元模型Fig.1 The finite model of bellows

300 系成型态波纹管的屈服强度与减薄率有关，不同减薄率下成型态波纹管的屈服强度按照如下公式计算[4]。

式中 Db——波纹管直边段内径；

Dm——波纹管平均直径。

根据表1 中波纹管的基本参数，可以得到成型态波纹管材料的屈服强度，成型态波纹管材料的性能参数如表2 所示。

表2 成型态波纹管材料性能参数Table 2 Material performance parameters of formed bellows

波纹管的刚度分析采用的边界条件为一端固支，另一端施加位移载荷。多层波纹管的非线性分析除了考虑材料的非线性行为，还涉及几何非线性和接触非线性。几何非线性通过打开大变形开关来实现，而接触非线性通过设置层与层之间的接触方式来实现，根据现有文献研究波纹管层与层之间比较合适的接触参数可设置为：层间间隙取0.1，法向接触刚度系数取0.1，层间摩擦系数取0.12，通过上述参数的设定来进行非线性分析[5]。

3 刚度试验

为了更好地了解无加强U 形波纹管在不等壁厚设计方案下的刚度性能，同时对有限元分析和理论公式计算的结果进行验证，进行了两种设计方案下波纹管的轴向刚度试验。

3.1 试验装置

对于无加强U 形波纹管的刚度试验，按GB/T 12777—2019 规定，实测刚度一般在无压力状态下进行测量，试验前要求波纹管以自由长度置于试验装置内，试验时通过一端固支，另一端施加轴向位移，到力与位移的曲线，由最终的力除以最大位移得到波纹管的实测刚度。

试验装置如图2 所示，试验时膨胀节一端固定在试验底座上，底座固定在100T 波纹管液压试验机底盘上，另一端通过拉压传感器与试验机上滑块连接，采用液压试验机对试验件施加位移，在试验件上方对称的两个位置各放置一个百分表来记录试验过程中的位移量，通过拉压传感器测量出波纹管的弹性反力反馈给称量显示仪。

图2 刚度试验Fig.2 Stiffness testing

3.2 试验方法

试验时注意保证传感器位于试验件上方中心，为了便于与有限元分析结果作对比，加载过程采用与有限元分析相同的过程，即加载位移为每次增加5mm 直至加载至波纹管的设计补偿位移40 mm，稳定后记录百分表数值和记录仪读数，由此得到一组力与位移的曲线，为了使试验数据更有说服力，每个设计方案的试验件进行3 次测量[6]，同时为了更好了解不等壁厚波纹管的轴向刚度性能，分别进行了拉伸位移和压缩位移时的轴向刚度测试。

4 结果分析

两种设计方案下的试验件测量所得力与位移关系分别如下图3 所示，其中图3a 所示为内部衬层为0.5 mm 时，在逐步拉伸到设计补偿位移40 mm 得到的力与位移关系曲线图，图3b 为内部衬层为1.0 mm时在逐步拉伸到设计补偿位移40 mm 得到的力与位移关系曲线图，从图中可以看出，随着位移增大，力逐渐增大，力与位移基本成线性关系。

图3 力与位移曲线（拉伸试验）Fig.3 Curves of force versus displacement （stretch test）

图4 为两个试验件在压缩到设计位移40mm 过程中的力与位移关系曲线。从图4 中可以看到，压缩时波纹管的力和位移的关系为非线性，从曲线斜率可以看出压缩过程中刚度逐渐减小，此现象与压缩时波纹管的层与层之间作用力增大有关，所衬层的厚度不影响压缩时力与位移变化趋势。

对拉伸得到的实测刚度值和理论刚度值进行对比，结果如表3 所示，从表中可以看出理论计算刚度和实测刚度较接近，偏差最大不到20%，由于试验测试时无法保证试验件整体的垂直度，测量刚度值应比实际刚度值偏大，按此原则实际刚度值和理论计算值的偏差会进一步缩小，因此可认为不等壁厚设计时可以按叠加原则进行总的刚度计算，按文中给出的刚度计算公式可以指导工程应用。

图4 力与位移曲线（压缩试验）Fig.4 Curves of force versus displacement (compression test)

表3 理论刚度和实测刚度对比Table 3 Comparison of the oretical stiffness and measured stiffness

图5 为两个试验件通过试验和有限元分析所得到的拉伸到设计位移40 mm 时力和位移关系对比曲线，其中试验数值为三次试验平均值，从图中可以看到拉伸时有限元分析得到的力和位移的关系与试验测试得到的变化趋势接近。

图6 为两个试验件通过试验和有限元分析所得到的压缩到设计位移40 mm 时力和位移关系对比曲线，试验数值为三次试验平均值，从图中可以看到压缩时有限元得到的力和位移的关系跟试验测试得到的变化趋势接近，均为非线性，曲线的斜率逐渐减小。

图5 拉伸时试验测试与有限元计算的力与位移曲线Fig.5 Force and displacement curve of experimental test and finite element calculation when stretch

将拉伸时分别经过有限元软件分析计算、理论公式计算及试验测试得到的实测刚度平均值进行对比，对比结果如表4 所示，序号1 为衬层0.5 mm 的设计方案，序号2 为衬层1.0 mm 的设计方案，两种设计方案下有限元分析得到的刚度值与实测计算刚度值最大偏差为10.4%，相比较于推导出的理论公式计算得到的刚度值更接近实测刚度值，这应与有限元分析时考虑了层与层之间的作用有关，衬层壁厚的增加会使层间作用力增大，因此可以认为采用有限元分析方法来计算不等壁厚波纹管的刚度是可行的，可以一定程度减少试验的数量。