铝制翼型表面磨损破坏特性研究

铝制翼型表面磨损破坏特性研究
张涛,陈次昌,吕东莉
（西南石油大学,成都610500）
[摘要]在文丘里装置内对铝制翼型开展磨损破坏试验研究，采用SEM观测翼型表面不同位置磨损破坏微观破坏形貌，测量翼型表面从头部到尾部宏观破坏深度。

对试验工况下装置内固液两相流场进行数值模拟研究，获取翼型表面流场内液固两相速度场及固相在流场内的浓度分布等特征。

依据数值模拟结果对翼型磨损破坏微观破坏模式及宏观破坏深度进行分析。

[关键词]翼型；磨损；两相流数值模拟；破坏模式
[中图分类号]TK730,TV149[文献标识码]A[文章编号]1000-3983（2010）02-0053-04
Study of the Cha ra cter istic of Silt Abr asion on Al Hydr ofoil
ZHANG Tao,CHEN Ci-chang,LV Dong-li
（Southwest Petroleum University,Chengdu610500,China）
Abstract:the silt abrasion experiments of aluminum hydrofoil are carried out in V enturi apparatus.The microscopic appearance of hydrofoil in different position is observed by SEM.The macroscopic wear depth of hydrofoil from the leading edge to the trailing edge is measured.The liquid-particle flow at experimental conditions is studied by CFD method,and the velocity field and concentration field etc.of particle are got.The analysis of microscopic damage pattern and macro wear depth is carried out according to the result of numerical simulation.
Key wor ds:hydrofoil;silt abrasion;numerical simulation of two phase flow;wear pattern
前言
材料的磨粒磨损问题到目前已有大量的研究成果。

但在水力机械内部，其过流表面磨损破坏的机理不但与材料本身性质相关，且很大程度受到水力机械内部固液两相流流动特性，其决定了固相颗粒作用于材料表面的角度、速度等参数的影响[1]。

文献[2]对材料磨损的微观形貌进行了定义，包括：犁耕类型造成沟槽、不同角度下冲击楔入材料表面形成坑以及对表面切削去掉材料。

不同装置（体现固相的作用方式变化）下不锈钢材料磨损表现出截然不同的特性。

文献[3]对多种材料试验指出，较高速度下60°~90°的射流冲角变化对磨损破坏几乎没有影响，但角度较小时破坏程度随角度减小而减小。

J.M.Neilson等发现柔性金属材料的磨损曲线上有一个低冲角的切削高峰，而在玻璃等脆性材料只有单纯的冲击变形磨损[4]。

基金项目：国家自然科学基金项目（90410013）和教育部博士点专项基金项目（665）联合资助
目前两相流的数值方法成为预测固相运动的重要方式，文献[5]采用欧拉－拉格朗日方法加Finnie的冲刷模型计算旋转圆盘内颗粒对材料的破坏速率。

文献[6]采用
p
A
kε两相湍流模型及图尔萨大学的泥沙磨损模型分析水轮机叶片的磨损主要特征，结果也给出了泥沙相的流动特征。

本文结合固液两相流数值模拟结果对文丘里装置中铝制翼型试样磨损破坏特性进行分析。

1试验装置与方法
1.1试验装置
试验在西南石油大学水力机械磨蚀闭式试验台进行，装置中试验段为收缩-扩散型（如示意图1）；试验对象为Raf6-1翼型试件，取自轧制铝板（成分：95%Al，5%Ni），弦长为120mm，其在流动方向的右侧有较小角度的扩散角。

1.2试验方法
试验前对翼型表面用2000目的砂纸依次沿垂直于弦长方向打磨，采用XL300电子扫描显微镜（SEM）观测表面的微观形貌；试验后，对图所示测试断面
200010011
（距前侧壁面20mm 处）翼型表面位置从头部到尾部进行显微观测；同时，测量该断面试验前后翼型上表面厚度的磨损深度。

磨损破坏试验时水中添加均径为45μm 的石英砂，试验段喉部平均速度为25m/s ，其含量为15g/l ，试验时间为15h ；试验过程中同时测量试验段进出口压力。

图1试验段示意图
2流场数值模拟方法
以含沙水流为介质，开展实验装置在试验条件下
的流场数值模拟。

采用Mixture 多相流模型。

混合模型的连续方程为
()()0=+m m m v t
ρρ式中，m ρ为混合密度；m v 为质量平均速度。

混合模型的动量方程可以通过对所有相各自的动量方程求和来获得，表示为
()()+
=
+p v v v t
m m m m m ρρ(
)[]+
++g v v m
T
m
m m
ρμ+
∑
=k
dr n
k k
dr k
k
v v F ,1
,ρα式中，n 为相数；F 为体积力；k α为第k 相的体积分数；k ρ为第k 相的密度；m μ为混合粘性；k dr v ,为第二相k 的漂移速度。

滑流速度qp v 被定义为第二相（p ）的速度相对于主相（q ）的速度：
q
p
qp v v v =则漂移速度和滑移速度可通过下式联系
q k
n
k m
k
k qp
p dr v v v ∑==1,ρρα从第二相的连续方程，可以得到第二相p 的体积
分数方程为：
()()()
p
dr p p
m p p
p
p
v v ,ραραρα=
+将单相标准ε
k 湍流模型扩展至多相湍流模型，单相SIMPLEC 算法扩展至多相流动，对装置中绕翼型的固液两相流进行数值模拟。

装置中翼型周围采用
C 型网格系统，网格数约40万，局部网格分布见图2。

边界条件：进口采用速度入口边界条件，根据实验给定来流速度、k 及ε估计值；出口采用压力边界条件，给出试验测量压力。

泥沙浓度通过给定进出口边界的固相体积浓度来实现。

壁面边界条件采用无滑
移条件，近壁处理采用标准壁面函数。

图2
装置内翼型头部位置网格
3
试验结果
3.1
材料原始表面形貌
图3为翼型表面原始微观形貌（弦长为左右方
向）。

竖直方向上有较长的打磨痕迹，其深度较浅且宽
度均在1~2μm 左右；此外，材料表面局部有一些缺陷小坑（大小约几个μm ）。

总体上讲，原始材料表面
较平整和致密。

图3翼型原始表面微观形貌
3.2
翼型表面磨损破坏形貌
以翼型头部为x 值零点，图4a~4d 分别为x 值约
2、5、10及25（mm ）处微观形貌。

图4a 中主要以大
角度冲击楔入坑（尺寸约在1~3μm 不等）为主，坑四周材料被挤压突出，表面较为粗糙；此外，材料表面还有特征不明显的小凹坑。

图4b 中主要是不同宽度和深度的中等角度冲击楔入坑（宽度可达4~5μm ），坑的后侧材料被切削后，部分淤积在坑的前侧。

图4c 与图4b 类似，但被冲击的角度较小。

图4d 中形成的冲击楔入坑很少且小，沿流动方向犁耕形成的沟槽和划痕增多，沟槽和划痕特征是窄、浅但较长且材料在
测试面
流动方向
测试面t
其两边淤积；同时，破坏的碎屑在材料表面杂乱无章地分布。

图5a~5f 分别为x 值约35、45、65、75、105及110（mm ）处微观形貌。

图5a 中破坏的形貌以犁耕形成的沟槽和小划痕为主，其宽度一般都小于远1μm ，但长度变化范围较大；同时，少量较小的小角度冲击楔入坑也存在。

图5b 与图5a 相比小角度的冲击楔入坑有所增多，尺寸特征也更为明显。

图5c 中小角度冲击切削坑开始占据了破坏类型的主导作用，部分坑的宽度可达1μm 左右，少量的划痕和沟槽仍存在，这些破坏形貌均向下侧以较小角度的偏转。

图5d 与图5c 较类似，但冲击切削坑更加明显。

图5e 及5f 中微观破坏形貌主要以大量中小等角度冲击切削坑和少量犁耕类沟槽相结合，其宽度均在1~2μm 之间，破坏后的材料使得翼型表面表现得非常粗糙。

图4翼型前侧磨损破坏微观形貌（流向由左向右）
图5翼型中后侧磨损破坏微观形貌（流向由左向右）
3.3
磨损破坏深度
15h 磨损破坏后，对测试断面位置测量翼型表面
材料失去的厚度。

由图6可以看出，翼型头部附近材
料在固相的冲击磨损下发生了整体的塑性变形，翼型在该位置的形状随之而发生了变化。

因此，本文对破坏程度的量化范围从距头部40mm 处开始。

图6翼型头部冲击变形
图7为断面2翼型中后部破坏深度值分布。

由图中看出，从x 约40~60mm 段翼型表面破坏程度较弱，且变化较小；从x 约60~110mm 位置，翼型表面破坏程度不断增大，并到x 约110mm 处达到最大破坏值；由x 约110mm 到尾部，翼型表面破坏程度有减小的趋
势，但破坏仍然较大。

图7翼型前侧到尾部破坏深度
4
两相流数值模拟结果
数值模拟完成后，取测试断面位置流场进行分
析，图8为断面2上沙粒速度分布矢量图。

在翼型头部处沙粒速度约20m/s 左右，沙粒以与翼型表面成接近90°的角度对其冲击磨损破坏；绕过翼型头部达到
约翼型最大厚度位置附近，沙粒运动速度达最大约25m/s 左右，其运动方向与翼型表面近平行；在沙粒体向翼型尾部运动过程中，随着流道的扩散，其运动速度由中后部约20m/s 左右减小到尾部的16m/s 左右。

此时，在流道扩散及翼型尾部上侧顺时针涡流的
影响下，翼型尾部上侧沙粒运动方向以约30°大小的夹角作用在翼型表面。

图9为测试断面流场固相含率的分布。

翼型头部
附近固相浓度最高，其原因在于固相以较高速度从喉
b
c
d
e
f
a
a b
c d
4050607080901001101200.00
0.010.020.030.04距翼型头部距离x/mm
单
位时间破坏深度h /
(m m /
h )
40
5060708090100110120
0.040.030.020.010.00
部流出后沿直线直接冲击翼型，形成沙粒的聚集；过头部后，两相流体作绕流运动，沙粒在离心力作用下沙粒偏离翼型表面，靠近翼型上部壁面浓度降低；接近翼型后部位置，较高浓度沙粒群在后侧回流区的影
响下向翼型表面偏转。

图8断面2
流场沙粒速度矢量
图9
断面2流场沙粒浓度分布
5结果分析
结合铝制翼型磨损破坏试验结果及相应两相流数
值模拟结果看出，翼型头部附近（x 从0~20mm ）沙粒含率最高，并以约20m/s 的速度大角度冲击翼型表面，头部微观破坏形貌为各种大角度的冲击楔入坑及冲击变形破坏，宏观上也表现出了整体的塑性变形；在翼型厚度最大位置附近（x 约30~60mm ），由于绕流效应的影响，虽然颗粒运动速度较大（约25m/s ），但其方向基本平行于翼型表面（图8）。

翼型表面微观破坏特征大多为刮擦或犁耕类的小沟槽，其破坏能力较弱，这在宏观的破坏深度上体现得较为明显；在翼型中后部，在涡流的作用下，沙粒以小角度对翼型壁面产生冲切作用，形成冲击楔入坑，利用颗粒进一步切削破坏，且在宏观上表现除了磨损破坏的峰值。

6结论
结合两相流场的数值模拟结果，本文对铝制翼型
的磨损破坏研究有下述结论：
（1）固相颗粒冲击角度反应出完全不同的磨损破坏模式，较大角度以冲击变形为主，较小角度以冲切破坏为主，而接近0°以刮擦划痕为主；
（2）铝翼型磨损破坏在大角度时会出现宏观上的变形，而在较小角度出现材料损失的峰值。

材料的性能在磨损破坏时成为重要的影响因素；
（3）水力机械内部的磨损破坏问题需要以两相流流动的准确预测为基础。

[参考文献]
[1]
张涛,陈次昌,等.空蚀磨损及联合作用研究进展[J].排灌机械,2006,（3）:47-53.[2]
J.Sato ，ami ，T .Okamura ，K.Oshima,et parison of Silt Erosion Characteristics of Materials Through Different Types of T est Methods.The 3rd Japan －China Joint Conference on Fluid Machinary .Osaka.1990.[3]
Kenichi Sugiyama,et al.Slurry Wear and cavitation erosion of thermal －sprayed cermets[J].wear.2005,
258:768-775.
[4]邓军,杨永全,等.水流含沙量对磨蚀的影响[J].
泥沙研究,2000,(4):65-68.[5]
徐姚,张政,等.旋转圆盘上液固两相流冲刷磨损数值模拟研究[J].北京化工大学学报,2002,29(3):12-16.
[6]唐学林,唐宏芬,吴玉林.水轮机转轮内固液两
相紊流场数值模拟及磨蚀预估[J].工程热物理学报,2001,22(1):51-54.
[收稿日期]2008-12-09
[作者简介]
张涛（1978-），2009年毕业于西南石油大学油气储运工程专业，博士，主要从事流体机械及工程、油气储运工程等领域内的科研工
作，讲师。