发明内容
针对现有技术的不足,本发明提供了一种用于防风抗振的高耸塔器疲劳寿命设计方法,该方法可为沿海地区高耸塔器的合理设计、安全运行提供参考。
为实现上述目的,本发明采用了以下技术方案:一种用于防风抗振的高耸塔器疲劳寿命设计方法,包括以下步骤:
1、一种用于防风抗振的高耸塔器疲劳寿命设计方法,其包括如下步骤:
S1、以离地面10m高度处的平均风速作为参考风速,设定参考风速为Vr,r=1,2,…,13;最小参考风速V1取4m/s,最大参考风速V13取28m/s,中间相邻的参考风速取值间隔为2m/s。
根据高耸塔器所在的地区,选择风速风向联合分布函数,所述风速风向联合分布函数的表达式如下:
式中,u为风速风向联合分布函数的风速自变量,单位为m/s;
Ur为风速风向联合分布函数的风速值,单位为m/s;
θj为0°到360°之间平均分布的16个风向角,且按顺时针依次增大,θj=(j-1)×22.5°,各个风向角θj代表不同的风向:如θ1=0°代表正北(N)风向,θ5=90°代表正东(E)方向,θ9=180°代表正南(S)风向,θ11=225°代表西南(SW)风向;
fj代表风向角θj的频度参数,aj、rj分别是风向角θj上风速分布的尺度参数、形状参数,fj、aj、rj采用下表中的数据:
风速风向联合分布函数中fj、aj、rj的默认参数值
j |
风向 |
fj |
aj |
rj |
|
j |
风向 |
fj |
aj |
rj |
1 |
N |
0.123 |
4.724 |
1.592 |
|
9 |
S |
0.050 |
4.374 |
1.267 |
2 |
NNE |
0.084 |
3.222 |
1.160 |
|
10 |
SSW |
0.050 |
5.140 |
1.927 |
3 |
NE |
0.063 |
2.630 |
1.073 |
|
11 |
SW |
0.035 |
4.215 |
1.448 |
4 |
ENE |
0.031 |
2.057 |
0.822 |
|
12 |
WSW |
0.016 |
3.868 |
1.377 |
5 |
E |
0.076 |
3.975 |
1.033 |
|
13 |
W |
0.006 |
3.810 |
1.315 |
6 |
ESE |
0.168 |
3.154 |
1.009 |
|
14 |
WNW |
0.003 |
3.622 |
1.596 |
7 |
SE |
0.147 |
2.351 |
0.908 |
|
15 |
NW |
0.023 |
3.597 |
1.133 |
8 |
SSE |
0.055 |
3.923 |
1.380 |
|
16 |
NNW |
0.069 |
4.347 |
1.593 |
将风速风向联合分布函数离散成与各参考风速Vr相关的多个连续区间,使得风速风向联合分布函数的每个区间[Ur,Ur+1)的中值为Vr,r=1,2,…,13;参考风速Vr在第j个风向角θj上的概率Prj为
Prj=P(u<Ur+1,θj)-P(u<Ur,θj)
S2、利用ANSYS软件,建立高耸塔器壳单元有限元模型,在塔顶施加弯矩进行静力学求解,高耸塔器轴向应力最大值所在的水平截面为高耸塔器的危险截面;
S3、按照无损检测及容许缺陷要求,设计若干种达到无损检测及容许缺陷要求的危险截面部位结构形式,并分别建立危险截面部位的局部实体单元ANSYS有限元模型,求解危险截面部位在上述若干种结构形式下的应力集中因子;
S4、模拟单位时间内高耸塔器在各参考风速Vr下的危险截面部位的不同方向上的名义应力时程响应,包括如下步骤:
1)、模拟风载荷
(1)将高耸塔器沿竖直方向划分为n段塔节,划分标准参考高耸塔器的设计文件,确定每段塔节的长度lh、外直径Dh、等效直径Deh以及每段塔节的中心点距地面高度zh,其中h=1,2,3,…,n;各段塔节的中心点即风载荷的模拟点,lh、Dh、Deh和zh的单位为m;
(2)计算各参考风速Vr下各模拟点的平均风速Vrh,其表达式如下:
(3)计算各参考风速Vr下的顺风向风载荷时程样本
首先计算参考风速Vr下的顺风向脉动风速谱矩阵中的上标a代指顺风向,其表达式如下
式中,矩阵元素
n为高耸塔器的塔节数;
Xr为脉动风速谱的无量纲莫宁坐标,
ω为圆频率,单位为rad/s;
为第h1、h2段塔节的中心点距地面高度,即第h1、h2个模拟点的高度,单位为m;
为参考风速Vr下第h1、h2个模拟点的平均风速,单位为m/s;
则各参考风速Vr下的顺风向脉动风速时程样本由谐波叠加法获得:
式中,h=1,2,…,n;
Δω=(ωup-ωdown)/N;
N为频率等分数,N≥1024;
ωup为截止频率,ωdown为起始频率,单位为rad/s;
为均匀分布于0~2π之间的随机相位角,单位为rad/s;
ωl=ωdown+(l-1/2)Δω,l=1,2,…,N;
为的下三角矩阵第h行、m列元素,单位为m·s1/2/rad1/2;
时间t的增量Δt必须满足:Δt≤π/ωup,t=Δt,2Δt,3Δt,…,T,单位为s;
T为单位时间,T为增量Δt的整数倍,即每次模拟样本的采样时间,单位为s;
于是,各参考风速Vr下的顺风向风载荷时程样本可写成如下形式:
式中,为高耸塔器第h段塔节的顺风向风载荷时程样本,单位为牛顿(N);
μs为体型系数,μs=0.7;
ρ为空气密度,取ρ=1.25Kg/m3;
lh为高耸塔器第h段塔节的长度,单位为m;
Deh为高耸塔器第h段塔节的等效直径,单位为m;
(4)计算各参考风速Vr下的横风向风载荷时程样本
首先计算参考风速Vr下的横风向风载荷谱矩阵中的c代指横风向,其表达式如下:
式中,矩阵元素
n为高耸塔器的塔节数;
ω为圆频率,单位为rad/s;
为第h1、h2个模拟点的高度,单位为m;
为参考风速Vr下的旋涡脱落力互谱密度,单位为N2·s/rad;
为参考风速Vr下的来流紊流力互谱密度,单位为N2·s/rad;
由以下二式获得:
式中,为参考风速Vr下的旋涡脱落力自谱密度,单位为N2·s/rad;;
R为距离参数,
ρ为空气密度,取ρ=1.25Kg/m3;
σCL为升力系数的方差,取0.14;
B是谱曲线带宽经验参数,取0.18;
ωsh为高耸塔器h段塔节的旋涡脱落圆频率,ωsh=Vrh/(10πDh),单位为rad/s;;
分别为第h1、h2段塔节的外直径,单位为m;
由以下两式获得:
式中,为来流紊流谱的无量纲莫宁坐标,
为来流紊流力自谱密度,单位为N2·s/rad;
Crh为阻力系数, 中间值按线性插值计算;
Rerh为雷诺数,Rerh=69000VrhDh;
u*r为摩擦速度,u*r=0.064Vr,单位为m/s;
则各参考风速Vr下的横风向单位长度的风载荷时程样本由谐波叠加法公式获得:
式中,h=1,2,…,n;
Δω=(ωup-ωdown)/N;
N为频率等分数,N≥1024;
ωup为截止频率,ωdown为起始频率,单位为rad/s;
为均匀分布于0~2π之间的随机相位角,单位为rad/s;
ωl=ωdown+(l-1/2)Δω,l=1,2,…,N;
为的下三角矩阵第h行、m列元素,单位为m·s1/2/rad1/2;
T为单位时间,即每次模拟样本的采样时间,单位为s;
时间t的增量Δt必须满足:Δt≤π/ωup,t=Δt,2Δt,3Δt,…,T,T为增量Δt的整数倍,单位为s;
于是,各参考风速Vr下各点横风向风载荷时程样本为:
式中,为高耸塔器第h段塔节的横风向风载荷时程样本,单位为牛顿(N);
lh为高耸塔器第h段塔节长度,单位为m;
2)、模拟名义应力时程响应
利用ANSYS有限元软件建立高耸塔器梁单元有限元模型,分别加载各参考风速Vr下的顺风向风载荷时程样和横风向风载荷时程样本进行FULL瞬态动力学模块求解;利用ANSYS软件APDL编程,提取各参考风速下危险截面的顺风向名义应力时程响应和横风向名义应力时程响应
求解高耸塔器顺风向风振响应时,高耸塔器的阻尼比为ζa=0.01;
求解高耸塔器横风向风振响应时,参考风速Vr下高耸塔器的阻尼比为
式中,为求解横风向风振响应时高耸塔器的阻尼比;
M*为高耸塔器的广义质量,
Mh为高耸塔器第h段塔节的质量,由设计文件给出,单位为Kg;
lh为高耸塔器第h段塔节的长度,单位为m;
为高耸塔器第h段塔节振型函数值,计算公式如下:
H为高耸塔器总高度,单位为m;
Arh为高耸塔器第h段塔节的气动力系数;
为第h段塔节的临界风速,单位为m/s;
ns为高耸塔器的一阶自振频率,可由ANSYS软件模态求解得出;
ah为第h段塔节的气动力参数,ah=a1ha2ha3h,h=1,2,3,…,n;
a1h、a2h、a3h分别由以下三式确定:
Rerh为雷诺数,Rerh=69000VrhDh;
3)、叠加名义应力时程响应
当参考风速为Vr、风向角为θj时,高耸塔器的危险截面在方向θk(风向角的角度同时也是危险截面上的方向位置)上单位时间T内的名义应力时程响应σrjk(t)为
式中,为顺风向名义应力时程响应,单位为MPa;
为横风向名义应力时程响应,单位为MPa;
S5、计算高耸塔器在步骤S3中不同结构形式和缺陷评定等级下的风致疲劳寿命,包括如下步骤:
1)、统计应力循环
(1)将步骤S4中的叠加后的名义应力时程响应σrjk(t)修正为
σirjk(t)=Ki·σrjk(t),i=1,2,…,p
式中,Ki为步骤S3中所求的第i种危险截面部位结构形式的应力集中因子;
p为S3中设计的危险截面部位结构形式个数;
σirjk(t)为修正后的应力时程响应,单位为MPa;
(2)利用雨流法统计σirjk(t)对应的应力循环;
2)、计算疲劳损伤
根据不同的焊接形式和焊接缺陷所对应的疲劳曲线计算单位时间T内上述应力循环对应的疲劳损伤Dirjk;
3)、计算疲劳寿命
结合离散后的风速风向联合分布,计算单位时间T内危险截面θk方向上引起的疲劳损伤,其表达式为
取上述疲劳损伤的最大值,记为则高耸塔器在第i种危险截面部位结构形式下的风致疲劳寿命为
在实际工作中,还可以增加以下三步:
S6、根据S5的结果判定是否存在满足设计寿命的评定区域,且该区域应包含制造允许的容许缺陷要求;
S7、若步骤S6判定结果为“否”,则需重新设计高耸塔器整体结构,并返回步骤S2;
S8、若步骤S6判定结果为“是”,则确定结构设计方案。
实际工作中,在步骤S5中计算疲劳损伤Dirjk时,针对不同的焊接形式和焊接缺陷,英国标准BSI PD5500:2009给出了D、E、F、F2、G、W六种评定等级下的疲劳曲线;使用时根据高耸塔器具体的焊接形式与焊接缺陷确定评定等级,再由评定等级下的疲劳曲线计算单位时间T内上述应力循环对应的疲劳损伤Dirjk即可。
所述不同的焊接形式和焊接缺陷所对应的疲劳曲线是反应自然规律的人工绘制的曲线,前述按照英国标准BSI PD5500:2009给出的D、E、F、F2、G、W六种评定等级下的疲劳曲线进行计算是为了提高计算效率,而不是必须按照英国标准BSI PD5500:2009进行计算。
本发明的主要优点在于:
1)、本发明模拟了高耸塔器在危险截面各方向上的风振动态响应过程,这与标准JB/T4710-2005《钢制塔式容器》与GB50135-2006《高耸结构设计规范》只能计算单一的等效风振响应值相比,能够再现高耸塔器的风振响应随时间变化的过程;
2)、本发明计算疲劳损伤时可以单独考虑应力集中因子的影响,使得该方法能够仅在计算一次风振响应的条件下,就能分析危险截面不同结构对风致疲劳寿命带来的影响,这对基于寿命预测的危险截面局部优化设计带来了很大的便利;
3)、危险截面的应力循环是根据雨流法对应力幅和平均应力进行直接统计,这比其他计算疲劳损伤的间接等效计算方法更加精确,此外还可以提供应力幅和平均应力的循环周次分布情况;
4)、该发明考虑了风速风向对高耸塔器风致疲劳寿命的影响,可以根据高耸塔器所在地的气象条件选择适合的风速风向联合分布函数,这使得高耸塔器风致疲劳寿命分析更加科学准确;
5)、本发明可以根据不同的无损检测要求和容许缺陷等级对其结构疲劳寿命进行校核,设计结果的预知性保证了塔器的安全运行。
具体实施方式
本实施例中的具体计算部分没有给出,在实际计算时请参考现有技术中的计算规则。
本发明步骤S5中的疲劳曲线、焊接形式与焊接缺陷评定等级的确定,可参考英国标准BSI5500-2009《Specification for Unfired fusion welded pressure vessels》与BS7910-2005《Guidance on methods for assessing the acceptability of flaws inmetallic structures》;技术方案中的具体计算过程需通过MATLAB软件编程实现,涉及有限元法求解部分需利用ANSYS有限元软件实现。
以下结合附图和实施例对本方法做进一步说明。
对沿海的某台高耸塔器进行用于抗风防振的疲劳寿命设计,该高耸塔器高约52.5m,筒体内直径2.6m,封头壁厚14mm,上筒体壁厚12mm,下筒体壁厚14mm,材料为Q345R。
步骤S1,首先设定参考风速为Vr=4,6,8,…,26,28m/s,相邻的参考风速间隔为2m/s;其次,步骤S1中的频度参数fj、尺度参数aj、形状参数rj除选用表格中的默认参数值以外,还可以选用该地区的风速风向联合分布函数的统计数据,本实施例中,由于高耸塔器所在地区尚无风速风向联合分布函数的统计数据,故频度参数fj、尺度参数aj、形状参数rj均选用默认参数值,并将其离散为与参考风速对应的风速分布区间。风向与风向角的对应关系如图2所示,其中NNE和E两个风向上的风速分布离散化结果如图3a、3b所示。
步骤S2,使用ANSYS有限元软件SHELL93壳单元确定危险截面位置,水平截面A-A为该高耸塔器的危险截面,参见图4a、4b所示。
步骤S3,首先设计危险截面部位不同的结构形式。标准JB/T4710-2005规定该处的焊脚高度ht应不小于1.7倍的裙座名义厚度δns(参见图5a),即焊脚高度ht不小于23.8mm,为此设定25mm、30mm、35mm、40mm、45mm五种可行的焊脚高度值,图5b所示的为焊脚高度ht的上限值,图5c所示的为焊脚高度ht的下限值。
其次,利用ANSYS有限元软件SOLID95实体单元,如图6a、6b、6c、6d所示,建立高耸塔器危险截面部位的局部有限元模型,分别计算上述五种焊脚高度值所对应的结构的应力集中因子,其结果如表1所示。
表1有限元法求解应力集中因子
焊脚高度ht/mm |
45 |
40 |
35 |
30 |
25 |
应力集中因子Ki |
3.393 |
3.858 |
4.391 |
5.009 |
5.751 |
步骤S4,分以下几步获得该高耸塔器各风向上的名义应力时程响应。
1)模拟风载荷
(1)将该塔器分为10段塔节,塔节的划分和模拟点的分布见图7,每段塔节的lh、Dh、Deh和zh参数如表2所示。
表2塔节的相关参数值(单位:m)
h |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
lh |
5.65 |
7.35 |
5.00 |
5.00 |
5.00 |
5.00 |
5.00 |
5.00 |
5.00 |
4.55 |
Dh |
3.50 |
2.70 |
2.70 |
2.70 |
2.70 |
2.70 |
2.70 |
2.70 |
2.70 |
2.70 |
Deh |
5.13 |
4.33 |
4.33 |
4.33 |
4.33 |
4.33 |
4.33 |
4.33 |
4.33 |
4.33 |
zh |
2.83 |
9.33 |
15.5 |
20.5 |
25.5 |
30.5 |
35.5 |
40.5 |
45.5 |
50.3 |
(2)计算各参考风速下各模拟点的平均风速,以参考风速Vr=16m/s为例,各模拟点的平均风速Vrh如表3所示。
表3各模拟点的平均风速(单位:m)
h |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
Vrh |
13.1 |
15.8 |
17.2 |
17.9 |
18.6 |
19.1 |
19.6 |
20.0 |
20.4 |
20.7 |
(3)利用MATLAB编程,模拟了该高耸塔器在各参考风速下顺风向风载荷时程样本。谐波叠加法模拟参数如下:起始频率为ωdown=0.001rad/s,截止频率为ωup=8rad/s,频率等分数为N=2048,单位时间为T=600s,采样时距为t=0.1s。以上参数值适用于各参考风速,图8示出了参考风速Vr=20m/s时,第4模拟点前200s的顺风向风载荷时程样本。
(4)利用MATLAB编程,模拟了该高耸塔器在各参考风速下横风向风载荷时程样本。谐波叠加法的模拟参数如下:起始频率为ωdown=0.001rad/s,截止频率为ωup=12rad/s,频率等分数为N=2048,单位时间为T=600s,采样时距为t=0.1s。以上参数值适用于各参考风速,但不同参考风速还需给出特定的参数值,如摩擦速度u*r、旋涡脱落圆频率ωsrh和阻力系数Crh。以参考风速Vr=24m/s为例,计算所得的摩擦速度为u*r=1.536m/s,参数ωsrh和Crh的值如表4所示,模拟获得的第4点的横风向风载荷时程样本如图9所示。
表4各模拟点的ωsrh、Crh参数值
2)模拟名义应力时程响应
首先,求解各参考风速下的顺风向风振响应。利用ANSYS软件建立该高耸塔器的梁单元有限元模型,阻尼比设为0.01,分别加载各参考风速下顺风向风载荷时程样本,对模型进行FULL瞬态动力学求解,并提取顺风向名义应力时程响应图10示出了参考风速为28m/s时高耸塔器危险截面的顺风向名义应力时程响应。
其次,求解各参考风速下的横风向风振响应。各参考风速下的高耸塔器横风向风振响应阻尼比的计算结果如表5所示,据此设定各参考风速下的阻尼比,分别加载各参考风速下横风向风载荷时程样本,对模型进行FULL瞬态动力学求解,并提取横风向名义应力时程响应图11示出了参考风速为28 m/s时高耸塔器危险截面的横风向名义应力时程响应。
表5 各参考风速下的横风向阻尼比
3)叠加名义应力时程响应
根据公式 将各参考风速下的顺风向名义应力时程响应与横风向名义应力时程响应,叠加为不同风向角和危险截面方向上的名义应力时程。图12a、12b示出了参考风速为28 m/s时两种夹角下危险截面的名义应力时程响应。
步骤S5,针对不同的结构设计形式,分别利用雨流法统计了各参考风速下危险截面各方向上的应力循环。参考风速为28 m/s、焊脚高度为45 mm时,危险截面某一夹角下T=600s内的应力循环统计情况如图13所示。最后,计算该高耸塔器在不同结构形式和缺陷评定等级下的风致疲劳寿命,计算结果如表6所示。
表6不同焊脚高度及评定等级下的疲劳寿命(单位:年)
步骤S6,判定是否满足设计条件。参考标准BSI 5500-2009中缺陷等级的评定部分,其规定危险截面处的单面焊焊接接头在无明显缺陷时应定为E级,而F、F2、G、W均代表不同等级的缺陷,W代表最严重的缺陷等级。已知该塔器的设计寿命为20年,则表3中灰色背景的区域均可以满足设计寿命要求,且部分结构容许F、F2级缺陷。故判定结果为“是”,即经初次结构设计后不必执行步骤S7与返回步骤S2。
步骤S8,确定结构设计方案。当焊脚高度ht为30、35mm时,其疲劳寿命在无明显缺陷的E等级下可以满足要求,但是这两种结构形式没有容许的缺陷等级,若该塔器在制造时产生了缺陷或者无损检测时未能发现存在的缺陷,其使用寿命可能达不到设计要求;当焊脚高度分别ht为40、45mm时,其疲劳寿命可以满足要求,并且该塔器在缺陷等级分别为F、F2时,其疲劳寿命仍可以满足设计要求。因此,应将危险截面处的焊接结构设计为焊脚高度ht在40mm以上,使之其风致疲劳寿命不仅可以满足设计要求,而且还有一定的缺陷容许等级,从而保证该塔器的安全运行。