弧形管疲劳寿命预测研究

2020-11-27刘文辉赵应龙

舰船科学技术 2020年5期

刘文辉，赵应龙，金著，李琛

(1. 海军工程大学舰船振动与噪声研究所，湖北武汉 430033；2. 船舶振动噪声重点实验室，湖北武汉 430033；3. 海军驻大连四二六厂军事代表室，辽宁大连 116000)

0 引言

随着舰船减振需求的不断提高，具有强大的位移补偿能力，小尺寸和优异的性能挠性接管被大量使用于船上[1]。其中弧形管指一种管体呈弧形的挠性接管，它在较小的空间内具有更强的位移补偿能力。疲劳破坏是弧形管损坏的一种失效方式，它指的是弧形管在外界应力或应变作用下伴随着疲劳过程的不断加深最终出现橡胶管体损伤的现象。在现有工程实践中，采用的疲劳可靠性评估方法一般为从一批成管中抽取几根进行给定次数的疲劳试验若为出现破坏则认为其具有一定可靠性，此法试验时间长并且因为其未从机理上分析疲劳失效原因，可置信性比较差。在汽车行业中，裂纹萌生法因其简单快速的特性而广泛用于橡胶疲劳寿命研究，该方法基于连续介质力学，从橡胶裂纹产生的机理出发，解释了疲劳失效的机理[2,3]，即橡胶产品本身存在很多内部缺陷，在疲劳载荷的作用下，内部缺陷逐渐扩展至一定程度出现破坏，宏观上表现为橡胶制品出现裂纹。同时该方法将有限元方法和材料级别的拉伸疲劳试验结合起来直接预测得到弧形管的疲劳寿命，大幅减少元器件级别的疲劳试验时间，具有简便快捷，易于计算的优点。

本文对以弧形管为研究对象，采用裂纹萌生法研究探索了一种弧形管疲劳寿命预测方法，研究结果可用于指导弧形管疲劳可靠性设计。

1 橡胶材料本构模型

为了建立准确可靠的有限元模型，需要对弧形管用橡胶材料的本构力学行为进行研究。针对橡胶材料复杂的力学行为，国内外学者利用应变能函数来表征橡胶材料的本构行为。在连续介质力学理论中，橡胶材料被认为是各项同性的超弹性材料，确定应变能和伸长率的关系，便可确定应力-应变关系。Ogden是一种常用的本构模型[4]：

式中：W为应变能函数；为材料常数；N为Ogden本构模型阶数；Di为材料常数，由橡胶的初始体积模量决定；J为橡胶变形前后体积比，对于主伸长率 λ，通常假设其为不可压缩的，，此时式（1）后半部分可以消去。

本文弧形管用橡胶进行了单轴条件下、等双轴条件下、平面条件下3种拉伸本构试验。

利用Abaqus有限元软件自带的拟合功能，发现Ogden三阶模型具有最佳的拟合效果（见图1），因此选择它作为本橡胶材料的本构模型。该模型的拟合参数为：μ1=1.19147， α1=0.72466， μ2=0.00015， α2=23.05513， μ3=0.00042， α3=-12.77194。

2 疲劳寿命损伤参量及其关系

在裂纹萌生法中，疲劳寿命通常是通过疲劳损伤参量表征的，疲劳损伤参量一般可以基于应变：应力和能量给出，其表征了初始裂纹的产生是否苛刻。

2.1 基于应变的疲劳损伤参量

以图2的均匀拉伸橡胶杆件为例，杆件未拉伸前的初始长度记为L，拉伸后的长度为l。由于橡胶材料的近似不可压缩性，图2所示的橡胶构件在变形之前和之后的构造可由下式描述：

图 1 三阶Odgen模型拟合结果Fig. 1 Fit result by 3 order Odgen model

图 2 单轴应力状态的橡胶杆件Fig. 2 Rubber rod in uniaxial tension deformation state

其中，Xi和xi（i=1，2，3）是橡胶杆件变形前后的位置。该方向的伸长率定义为在橡胶杆件变形前后的方向（轴向）上的长度比。当 λ＞1时，表示杆件处于拉伸状态时；当 0＜λ＜1时，表示杆件处于压缩状态；当 λ =1时，表示杆件处于不变形状态。

最大主Green-Lagrange应变 εG（Green-Lagrange应变张量E最大的主值）和最大主Almansi-Euler应变εA（Almansi-Euler应变张量e最大的主值）分别为：

工程应变 εE表示为：

2.2 基于能量的疲劳损伤参量

应变能密度的峰值也常被用来选作橡胶材料的疲劳损伤参量。应变能密度（W）的定义式为：

其中，σ和ε为一组共轭且功能完全确定一点应力和应变状态的应力张量和应变张量。

对于单轴应力状态，上式为

其中， σE为工程应力， εE为工程应变。

Mars[6]通过研究发现应变能密度峰值与工程应变峰值之间存在着幂函数关系。本文测试了4种工程应变水平下的应力应变曲线，在经过5次循环后取最后一个循环加载段的曲线。将试验数据导入Matlab中，通过梯形数值积分函数TRAPZ对其进行积分求得对应的应变能密度峰值。基于非线性最小乘法得到的拟合模型如图3所示，其数学表达式为：

式（10）的拟合精度为0.998 5，较为接近1，说明该模型精度较高，是可取的。

3 疲劳寿命预测模型的建立

在裂纹萌生法中，为了构建疲劳损伤参量与疲劳寿命之间的关系，需要对橡胶材料哑铃型试片进行单轴拉伸疲劳试验。拟合得到疲劳损伤参量与疲劳寿命之间的数学关系，即疲劳寿命预测模型。

3.1 单轴拉伸疲劳试验

图 3 单轴拉伸试验结果及拟合结果Fig. 3 Strain energy density experimental data

本次橡胶材料的疲劳性能的试验研究参照《GT/T 1600-2008 硫化橡胶伸张疲劳的测定》，试验试样为哑铃型试片。在单轴循环应变载荷的作用下，将样品反复拉伸至断裂。样品循环的应变比为0，样品断裂时记录的循环次数为疲劳寿命。本次试验一共开展了5种不同拉伸载荷工况下的疲劳试验，应变幅值分别为80%，100%，120%，140%，160%，各个工况下测试试片数量为12片。记录每一组中哑铃型试样中间狭长部位断裂时试验的循环次数，若试件在非狭长部分断裂，则剔除此样件。

根据ASTM标准，橡胶材料的疲劳寿命服从对数正态分布，平均寿命的计算公式为：

式中：n为相同载荷工况的一组数据经去除离散数据后剩余的个数；Ni为第i片试片对应的疲劳寿命。将单轴拉伸试验得到的结果按照式（11）进行处理后，得到结果如表1所示。

表中Nf，max，Nf，min表示每种载荷条件下测量的12个橡胶试样的最大值和最小值。Nf，min/Nf，max的大小体现试验结果的分散性。从试验结果可以看出，橡胶拉伸疲劳试验的分散度在（2.38，6.71）之间，基本满足橡胶材料疲劳寿命分布要求。

3.2 疲劳寿命预测模型

根据相关文献，疲劳裂纹萌生寿命与疲劳损伤参量之间通常满足幂律[7]：

表 1 橡胶哑铃型试件片拉伸疲劳试验结果Tab. 1 Tensile fatigue results of rubber

其中：P为疲劳损伤参量：Nf是裂纹萌生寿命；K和b为材料常数。基于最小二乘法，对不同损伤参量建立橡胶试片建立寿命预测模型如表2和图4所示。

表 2 疲劳寿命模型拟合结果Tab. 2 Fitting results of fatigue life model

图 4 寿命预测模型拟合结果Fig. 4 Life prediction model fitting result

可以看出，不同的疲劳损伤参量拟合得到的寿命预测模型大相径庭，其中r2表示预测模型的拟合精度，r2的值越接近1，说明该模型的拟合效果越好。可以看出，基于最大主伸长率建立的寿命预测模型拟合效果最好。

4 弧形管疲劳特性分析

4.1 弧形管有限元模型的建立

图 5 弧形管结构和管体模型Fig. 5 Flexible pipe structure and model

弧形体挠性接管主要由內胶层，芳纶帘线骨架层，外胶层和法兰接头4个部分组成[8-9]。骨架层是由芳纶帘线组成，是主要的承载部件，包在内外胶层之间，分为多层，各层严格按照一定角度交叉包覆，呈现出各向异性的特点。各胶层使用的橡胶是含有各种填料的氯丁橡胶。

在有限元软件Abaqus中选用Rebar单元M3D4R来模拟帘线、三维应力单元C3D8RH来模拟橡胶、三维应力单元C3D8R来模拟上下盖板。帘线材料按照刚性材料定义，杨氏模量为33.882 GPa，泊松比为0.3。缠绕角度理论缠绕等效缠绕平衡角度52°进行布置，在此缠绕角度下，仅受内压时，弧形管将不会产生轴向力和位移，具有较好的平衡性。因为帘线层包络在内外橡胶层中，所以在Abaqus中采取内置区域定义橡胶层和帘线层之间的关系。橡胶材料选用超弹性本构模型，本构模型参数按选用前文建立的本构模型——三阶Odgen模型中的参数。由于本文不考虑盖板的影响，故将上下盖板都设置为刚体。

将该有限元模型的刚度仿真结果与刚度特性试验结果进行比较发现相对误差不超过5%，一致性较好，故仿真可代替试验进行下一阶段研究。

有限元仿真模型的精度决定了仿真结果是否可信，其影响因素有很多如网格划分的质量，本构模型的选用等。通过橡胶本构试验拟合得出的Odgen三阶模型大大提高了有限元模型精度，刚度试验与仿真值的对比也反映出此现象。