首页 > 论文范文 > 工程工业论文 > 动力工程论文 > 水力机械论文 > 水轮机蝴蝶阀的蝶板截面多目标优化方案设计研究

水轮机蝴蝶阀的蝶板截面多目标优化方案设计研究

2020-07-02 348 上传者：管理员

摘要：水轮机蝴蝶阀是一种流量控制装置,蝶板是其重要零件,对流经它的液流优显著影响。为了深入研究蝶板截面最优形状,采用有限元方法进行截面拓扑优化,并通过正交阵列和逐步回归响应面法建立压力损失系数和最大应力的替代模型。采用权衡方法,求得多目标的最优解。结果表明有限元拓扑优化结果的重量减少40%,蝴蝶阀最小压力损失系数为0.4755,最大应力为134.09MPa。证明该设计是兼顾流动特性和结构安全性的蝴蝶阀蝶板多目标优化设计。

关键词：
响应面法
多目标优化
拓扑优化
水力机械
蝴蝶阀
蝶板截面
加入收藏

水轮机蝴蝶阀通常用于阀门所需压降较低的应用中,其主要功能是控制流量或压力的开启与关闭。蝶板通过一根阀轴定位在中心,阀轴使蝶板平行于或垂直于液流,并且阀轴连接到阀门外部的手柄或其他操作器和结构。与球阀不同,蝶板始终存在于流场中,因此无论蝶板位置如何,压降始终存在。由于蝴蝶阀对流经它的液流影响显著,因此各国学者对蝴蝶阀进行了大量研究。早在1995年,小川和木村[1]研究了阀门扭矩特性和蝶板引起的压降。JANUSZ和CZESLAW[2]通过使用实验和数值结合的方法研究了一种常用的蝴蝶阀,该阀具有薄而平的尖锐蝶板且堵塞比d/D为0.947。洪宇[3]通过对水轮机进水蝶阀进行仿真分析,对水轮机进水蝶阀的结构进行优化,提高其工作性能和可靠性。郭建章等[4]基于Fluent软件提供的湍流模型与离散相模型,对中线蝶阀内部流场进行了数值模拟研究,考察了阀门开度对蝶阀壁面及蝶板磨损的影响。张融等[5]利用ANSYS Work Bench对其结构强度及抗震能力进行了仿真分析并对阀轴进行结构优化改进。以上这类研究主要集中在蝴蝶阀的流体分析或者整体结构分析上,而不是聚焦在蝶板形状上。因此,本文提出了一种兼顾流动特性和结构安全性的蝴蝶阀蝶板多目标优化设计。

研究的过程包括三个部分。首先,依据阀门中的压力损失系数和最大应力进行拓扑优化以识别蝴蝶蝶板的截面形状。其次,采用正交阵列和逐步回归响应面法建立压力损失系数和最大应力的替代模型。最后采用权衡方法,实现多目标优化及设计变量优化。

1、试验方法与设计

在本研究中,设计条件要求该蝶阀既具有较小的压力损失系数(k≤0.8),又具有足够的结构系数(SF≥1.5)以用于特殊用途。这些关键因素并未根据用户的要求进行考虑。根据期望,蝶板的直径和最大厚度分别定义为100mm和20mm。

1.1正交实验和替代模型

替代模型方法使用一些基本函数来近似设计问题中涉及的高度复杂的目标和约束函数[6]。替代模型的主要步骤通常包括选择实验设计以生成采样数据;选择模型代表数据;最佳地将模型拟合到采样数据这三个步骤,其中的每一个都有多种选择,并且每个步骤都有一些优点和缺点。发现没有最佳的单一替代模型可用来适应所有工程问题。在这项工作中,选择正交阵列作为实验组,采用逐步回归的响应面法作为仿真试验设计。

根据阀门中的压力损失系数和最大应力进行拓扑优化,选择影响蝶板截面形状的七个设计变量,它们显示在图1的中间半径R,臂长L1和L2,蝶板边缘的厚度T1和T2以及中间厚度W1和W2。

图1蝶板截面主要形状变量

通常,以上变量需要使用全因子设计进行36次的分析。但是,通过使用L18(21×37)正交排列组合的阵列,实验组数可以减少到18[7]。

1.2优化方案

本文利用拓扑优化和尺寸优化来获得性能更好的蝴蝶阀蝶板。尺寸优化在拓扑优化的基础上确定蝴蝶阀的详细形状。本研究有两个需要优化的目标,一个是压力损失系数,另一个是作用在蝶板上的应力。多目标优化是同时优化两个或两个以上受某些约束冲突的目标的过程。一般的多目标优化表示为:

公式1

其中,f1(x),f2(x),...,fk(x)是k的目标函数;g(x)和h(x)分别是不等式约束,x=(x1,x2,...,xn)T是n个设计变量,x∈S定义设计空间。

为了进行多目标优化,本研究采用权衡方法。该方法允许设计人员按照重要性从1到k降序对目标进行排列。然后,每个目标函数都会受到最小化的约束,该约束不允许新函数的最小值超过先前函数的最小值的规定分数[8,9]。一旦从f1(x)(最重要)到fk(x)(最不重要)对k个目标进行了排序,则求解过程如下:

步骤1:假定不存在所有其他目标函数,则在原始约束集的约束下找到f1(x)及其对应变量x1的最佳点。

步骤2:定义一个值为εj,然后在附加约束条件下找到第j个目标函数fj(x)的最优点:

fj−1(x)≤(1±εj−1100)fj−1(xj−1),j=2,...,k(2)

其中εj是假设的函数增量系数,%。

2、实验结果与分析

2.1拓扑优化

图2给出了蝴蝶阀蝶板截面的有限元模型,用于拓扑优化。该模型的材料是STS316,其杨氏模量E=193GPa、泊松比v=0.27。黑色区域表示无法删除的元素,浅灰色区域表示可以在拓扑优化中将其删除。在各个方向上带有箭头的短线表示已知实验数据分析获得的各种流体压力。应该注意的是,在阀门下游也有压力,但是下游作用在蝶板表面上的压力与上游表面上的压力相比太小,因此可以忽略。

图2二维蝴蝶阀蝶板的有限元模型

图3显示了拓扑优化结果的重量减少40%。圆圈D中浅灰色区域可以删除,对整体配置影响不大。也就是说,蝶板的周边可以做得很薄。

图3拓扑优化

根据拓扑优化的结果,新模型(见图4)。相应的流场分析在相同的边界条件下进行应力分析。标记为“MX”的最大应力为50.5MPa,出现在拐角区域,这表明流动的流体在前场相反的一半上产生更高的阻力。由于STS316的屈服应力为205MPa,因此该新模型的安全系数大于4.0(K1=1.5779),表明可以引用此模型以作为原型,下一步将基于该原型进行尺寸优化。

图4拓扑优化后的新有限元模型及其应力分布

2.2逐步回归的响应面法

获得实验结果后,在线性函数中构建压力损失系数的替代模型,并在二次函数中构建应力的替代模型,如下所示[10]:

表1适用于K1模型的方差分析

表2拟合σmax模型的方差分析

表1和表2分别显示了K1和σmax的方差分析。替代模型与计算机分析之间的K1和σmax的总误差R2分别为93.06%和68.53%,R2调整分别是88.09%和46.49%。因此,替代模型的准确性被认为足以进行本研究。

2.3权衡方法的多目标优化

如上所述,本研究的设计目标是最小压力损失系数以及最大应力。因此,目标和约束函数可以定义如下:

其中x=[R,L1,L2,W1,W2,T1,T2]T,数字“1”和“3”分别代表下限和上限。

基于这两个替代模型,使用权衡方法求解多目标优化。此方法将目标按重要性降序排列。然后,每个目标函数都会受到最小化的约束,该约束不允许新函数的最小值超过先前函数的最小值的规定分数。在本研究中,解决方法如下:

步骤1:

下一步是找到第二个目标函数的最优值,它受一个附加约束。

步骤2:

经过上述两步解决方法,可以获得最佳解决方案和最优结果(见表3)。x的最优水平(不是实际尺寸)在表4可以发现,T1、W1、L1的尺寸分别接近其上限。相反,T2、W2、W2的尺寸分别接近其下限。

表3最优结果

表4设计变量的最优结果

3、结论

本文提出了在多目标优化下水轮机蝴蝶阀蝶板设计过程。该过程根据初始设计给出压力损失系数,确定阀门行业中常用的系数,采用有限元方法进行拓扑优化,得到基本形状,然后采用正交阵列和逐步回归响应面法建立压力损失系数和最大应力的替代模型。最后采用权衡方法,根据公式直接求得多目标的变量最优解。权衡方法与响应面模型相结合,可以准确并有效地预测最佳条件。

需要注意的是,本文仅对蝶板的截面进行了优化,采用的是二维有限元模型。将来的文章中应考虑包括蝶板、流体在内的三维有限元分析。

参考文献：

[3]洪宇.水轮机进水蝶阀关键零件的分析与优化[D].武汉:武汉轻工大学,2018.

[4]郭建章,耿跃虎.固液两相流下蝶阀的磨损特性研究与结构优化[J].机械制造,2019,57(3):42-46,55.

[5]张融,殷德奎,徐旭升,等.基于ANSYS Work Bench的水轮机进水蝶阀的结构分析与优化[J].机械,2019,46(12):36-41.

[6]刘万刚,宋述芳,樊维超,等.模型替代方法在优化设计中的应用[J].中国工程机械学报,2017,15(2):119-124.

[7]洪新阳,王士灿,周昊奕,等.20CrMnTi倒挡中间齿轮精锻成形工艺模拟与优化[J].浙江水利水电学院学报,2019,31(4):67-73.

[8]赵冬青,张双全,高向正.基于多元统计的转轮优化分析[J].浙江水利水电学院学报,2018,30(1):24-27,35.

[9]何勇,张浩,杨晓霞,等.双排钢板桩围堰钢拉杆有限元分析[J].浙江水利水电学院学报,2016,28(2):58-63.

[10]宋志远,刘淑梅,莫壮壮,等.基于响应面法的蓄能器壳体工艺优化[J].轻工机械,2020,38(2):90-94.

牛萌萌,李宪华.水轮机蝴蝶阀蝶板截面多目标优化设计[J].浙江水利水电学院学报,2020,32(03):38-41.