强震作用下损伤井架结构的力学演化规律

杨明飞，程华瑞，沙志平

(安徽理工大学土木建筑学院，安徽 淮南 232001)

井架作为凿井设备的主体结构，其安全对凿井工程至关重要。在长期的服役过程中，结构材料经过腐蚀介质的侵蚀，在表面逐渐出现局部锈蚀现象，这对结构的力学性能会产生不利影响。近年，国内外诸多学者对钢材腐蚀展开了大量研究。文献[7-8]研究了经锈蚀处理钢材后力学性能的退化，发现锈蚀率对钢材屈服强度有显著变化，锈蚀率越高，钢材强度退化越明显。文献[9]研究了在潮湿环境下，不同锈蚀位置和锈蚀时间等参数对钢管力学性能的影响，并根据试验结果，拟合出钢管在锈蚀状态下的极限承载力计算式。文献[10]研究了钢梁在无规则腐蚀作用下，其极限承载能力下降。试验研究表明，通过失重率可以反映钢材有效截面积减小导致应力集中对钢构件材料性质的弱化程度。文献[12]研究表明随着腐蚀龄期的延长，钢框架整体刚度减小，其整体力学性能减弱。文献[13]研究了锈蚀损伤对钢材的影响，发现复合缺陷损伤对结构稳定性明显不利。文献[14]通过低周疲劳试验研究了H型钢在锈蚀状态下的力学性能，试验结果表明构件的滞回性能随着锈蚀率的增加而降低。文献[15-16]研究发现可以建立钢结构损伤模型，反映锈蚀对钢框架抗震性能的影响。目前的研究针对静力情况较多，而对于动力荷载作用下损伤钢结构的力学研究偏少。

为了分析局部损伤对井架结构的力学影响，本文通过有限元软件ANSYS建立井架结构的数值模型。以改变结构局部构件的失重率来反映结构的损伤程度，接着对井架结构进行静力分析和强震作用下的动力分析，得出局部损伤下井架结构力学演化规律，其相关结论为现役井架结构的设计和维护提供参考。

1 井架结构案例

1.1 工程概况

井架结构所处区域7度设防、Ⅱ类场地，服役超过10a。由于气候复杂、环境恶劣等条件影响，井架结构受到酸性介质和积水的侵蚀，导致结构立柱的局部表面出现锈蚀点坑，进而影响到井架结构的力学性能。因此，需要对井架结构的力学性能进行评估分析。

1.2 结构平立面布置

井架结构立柱、首层和中间层平台布置如图1(a)所示，立柱为圆形钢管，一、二层截面为

φ

450mm×25mm，第三层截面为

φ

325mm×25mm。井架顶层为天轮平台，如图1(b)所示，天轮平台上、下布置主钢梁(Q345)，钢梁主要以N-0至N-5为主，其中N-0采用组合梁,N-1至N-5均采用I56a，天轮平台主梁与次梁连接方式为焊接。井架相邻层间布置斜撑构件，以增大井架结构的侧向刚度。

2 井架结构数值分析

2.1 模型荷载简化

井架结构上的主要荷载有吊绳拉力、滑架托平台自重、翻矸平台自重、井架整体自重和水平风荷载等。其中可以将吊绳拉力、滑架托平台自重、翻矸平台自重和井架整体自重等简化为集中荷载考虑，加载点如图1中符号⊕所示。其中，井架首层平台横梁上各点均竖向加载2kN，中间层平台横梁各点均竖向加载10.04kN，天轮平台横梁水平方向上共加载1 841.93kN，竖向共加载2 979.87kN，天轮平台与四根立柱交点处在竖向均载为8.76kN。水平方向风载根据《矿山井架设计规范》(GB 50385-2006)进行计算,在井架首层、中间层和顶层平台处的简化结果如图1(a)所示。

(a)井架杆件布置图

(b) 天轮平台布置图图1 井架结构布置图(单位：mm)

2.2 井架模型建立

利用ANSYS软件建立的无损井架结构模型如图2(a)所示。其中，井架结构使用非线性随动强化模型，构件均采用BEAM188单元模拟，泊松比为0.3，弹性模量为206GPa，井架结构模型中天轮平台上横梁的屈服强度为255MPa，立柱和其他构件的屈服强度为345MPa，节点采用刚接。结合图1(a)，损伤井架结构损伤部位为井架结构一、二层的四根立柱，如图2(b)所示。

(a)无损井架模型

(b) 损伤井架发生损伤部位图2 井架模型

根据立柱损伤情况，对其有效截面、弹模和屈服强度进行折减。其中，立柱的有效截面折减是通过体积相等原则进行，即材料损伤失去的体积与钢管柱外截面折减的体积相等，材料弹模和屈服强度的折减是由损伤钢管构件的轴压试验数据分析得出。具体折减公式如下

(1)

E

′=(1-0

.

75

D

)·

E

(2)

Δt

=

D

·

t

(3)

为了得出损伤对结构力学的影响规律，失重率设置为0、15%、30%、45%、60%。钢材属性折减后数值如表1所示。

表1 腐蚀折减后的钢材属性

3 井架结构力学分析

3.1 静力分析

对模型进行静力分析，得到不同失重率下的井架结构在损伤处最大位移和顶层最大应力曲线，如图3所示。

图3 井架损伤区域最大位移和顶层最大应力

从图3得知， 井架结构在静力荷载作用下， 随失重率的增加， 损伤处最大应力和井架顶层最大位移均增加， 呈现出二次函数关系。 无损井架结构在一、 二层立柱部位最大应力为19.9MPa，发生在立柱底部，同时无损井架在井架顶层最大位移为2.287mm。当失重率达到60%时，立柱损伤区域最大应力为32.2MPa，和无损伤结构相比增加61.3%，最大应力出现在井架一、二层立柱与中间层平台横梁交接处，顶层最大位移为2.841mm，和无损伤结构相比增加24.2%。

3.2 模态分析

对无损伤井架结构和失重率为60%的损伤井架结构进行模态分析，分析结果如表2所示。

表2 无损伤井架/损伤井架结构前10阶模态

由表2可以看出，损伤井架结构的模态频率发生了改变，说明结构出现损伤。对前10阶井架结构模态进行对比，发现由于井架结构的刚度较大，因此对前四阶模态频率的影响较小，同时振型没有明显变化。

3.3 强震作用分析

在强震分析中，选用无损伤井架结构和失重率分别为60%的损伤井架模型进行分析。地震峰值加速度按照现有规范进行选取(见表3)。本文选取三条地震波，分别为Elcentro波、Taft波和一条Ⅱ类场地的人工波对井架结构进行三向激励，峰值加速度的

x

∶

y

∶

z

=1∶0.85∶0.65。现以Taft波为例进行说明。

表3 地震波加速度峰值 cm·s-2

在不同峰值加速度的作用下，得到无损伤和有损伤井架结构在损伤处最大应力和井架结构顶层最大位移曲线，如图4所示。

(a)井架一、二层立柱处最大应力 (b) 井架顶层最大位移图4 各地震波峰值加速度下井架结构损伤处最大应力和顶层最大位移

通过图4可知，随地震波峰值加速度的增加，井架结构损伤处最大应力和顶层最大位移均呈现出线性增加的趋势。

当地震波峰值加速度为35cm/s时，无损伤井架一、二层立柱最大应力为1.23MPa，井架结构最大应力为13.0MPa；损伤井架损伤立柱的最大应力为1.36MPa， 井架结构最大应力为14.8MPa。地震波峰值加速度为620cm/s时，无损伤井架一、二层立柱处最大应力是24.9MPa， 井架结构最大应力为111.0MPa；损伤井架损伤立柱的最大应力为 25.1MPa，而井架结构最大应力达到141.0MPa。表明随地震波峰值加速度的增加，无损伤井架和损伤井架的损伤处应力以及井架结构最大应力也随之增大。但是文中所设置的立柱失重率对井架损伤处最大应力影响不明显，表明井架结构刚度较大，仍处于安全范围。由图4(b)知，不同失重率下的井架顶层最大位移均在地震波峰值加速度为620cm/s的时候出现，最大位移达到33.936mm，和无损井架结构顶层最大位移为19.185mm相比，增幅达到76.89%，说明失重率对井架结构顶层最大位移影响显著。

通过计算分析，得到井架结构层间位移角(按弹性方法计算的风荷载或多遇地震标准值作用下的楼层层间最大水平位移与层高之比)随地震波峰值加速度增加的变化曲线，如图5所示。

图5 不同强震作用下的井架最大层间位移角

由图5知，发现强震荷载作用下，井架层间位移角随地震波峰值加速度的增加而线性增长。无损井架最大层间位移角为1/926，发生在井架一、二层立柱底部；而失重率为60%的损伤井架最大层间位移角为1/752，发生在井架一、二层立柱中间部位，层间位移角数值仍满足《构筑物抗震设计规范》要求。

4 结论

(1)在静力荷载作用下，井架结构随失重率的增加，其顶层最大位移和损伤处的最大应力均明显增大，并且呈现出二次函数关系。当失重率达到60%的情况下，井架结构在损伤处最大应力发生的位置出现上移。

(2)模态分析显示，损伤井架结构的模态频率出现变化，可以证实结构受损，但文中设置的立柱局部损伤对井架结构的模态影响不明显。

(3)在地震作用下，损伤井架的顶层最大位移和损伤处最大应力均呈现出线性增加。随着地震峰值加速度的增加，失重率对井架结构的损伤处最大应力影响较小，但是对井架顶层最大位移影响显著。

(4)井架层间位移角随地震波峰值加速度增加而出现线性增长。在失重率为60%、地震波峰值加速度为620cm/s时，井架层间位移角达到1/752，仍满足《构筑物抗震设计规范》要求。