[发明专利]基于虚拟厚度的反击式转轮全局特性数值计算方法和系统

说明书

本发明公开了一种基于虚拟厚度的反击式转轮全局特性数值计算方法和系统。步骤包括：数据输入；坐标变换；过水断面形成线构建；轴面速度计算；轴面流线计算；S1、S2流面构建；空间相对流线构建及分割；空间相对水流角计算；空间相对流线上各点轴面速度计算；增厚系数试算；虚拟厚度叶片排挤系数计算；欧拉能量数值求解；效率数值求解；CFD数值模拟结果验证；判断精度；增厚系数插值；将计算结果输入到Excel表格进行显示和存储；计算结束。所述系统用于实现对应功能。该数值计算方法和对应的系统提高了计算过程精简性和高效性；采用模块化设计思路，各子模块既可独立运行又相互之间紧密联系，计算流程清晰；适用于工作在非压缩流体中的反击式叶轮旋转机械，具有广泛应用对象，涉及众多行业领域。

技术领域

本发明涉及反击式旋转机械工作特性数值计算领域，尤其涉及一种反击式转轮全局工况下效率和水力特性的快速预测方法。

背景技术

水力发电因其输出电能稳定、调节反应迅速已成为电网保证高质量电能不可或缺的能源结构，因此水电为了高效发挥其调节性能，要求其能量转换的核心元件——转轮工作在偏离最优工况。我国，始于上世纪70年代对转轮性能的试验研究和始于90年代的基于有限元软件的数值模拟研究经过多年发展日趋成熟，但由于与这两种研究方法相关的辅助设备未能得到充分发展，导致的高成本投入、长时间耗资，比如模型试验的长周期、数值模拟的高计算性能要求等，已成为新建水力发电项目初期阶段不可忽视的成本投入。

发明内容

本发明所要解决的技术问题是：提出基于虚拟厚度的反击式转轮全局特性数值计算方法和系统，缩短研发周期，降低成本。

本发明为解决上述技术问题采用以下技术方案：

本发明提出一种基于虚拟厚度的反击式转轮全局特性数值计算方法，进一步地：步骤包括：

S001、数据输入：读取反击式转轮的全局参数和木模图参数；

S002、坐标变换：将笛卡尔坐标系内的空间叶片投影到Z-r平面上，且遵循以下旋转投影公式：

r_i²＝x_i²+y_i²，z_i＝z_i

其中，(x_i，y_i，z_i)为点在笛卡尔坐标系内的坐标，(z_i，r_i)为对应的点在Z-r平面内的坐标；

将笛卡尔坐标系内坐标(x,y,z)变换到S1流面坐标系内坐标(m,θ)，其中m为S1流面的母线，θ为S1流面圆周方位角；流线的坐标变换并无量纲化处理遵循以下公式：

式中：x₀，y₀为流线起始端点坐标分量，(m_i，θ_i)为笛卡尔坐标系内点(x_i，y_i，z_i)在S1流面坐标系内的坐标；

S003、过水断面形成线构建：在反击式转轮轴面流道内作内切于上冠和下环流线的公切圆，圆心为O，切点为A和B，连接OA、OB和AB，作通过A、B并分别与OA、OB相切的圆弧即为过水断面形成线；

S004、轴面速度计算：求得过水断面形成线形成过水断面时的旋转半径，进一步求得对应的过水断面面积和轴面速度；针对多条过水断面形成线，求得轴面速度在轴面流道中的离散分布；所述过水断面是由过水断面形成线绕转轮转轴旋转一周形成的回转面；