一种高效的涡轮叶片内部冷却结构优化方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN108595797A(43)申请公布日2018.09.28(21)申请号201810322216.8(22)申请日2018.04.11(71)申请人中国工程物理研究院化工材料研究所地址621000四川省绵阳市绵山路64号(72)发明人宋乙丹谢公南路中华陈科全陈红霞李兴隆陈翔刘柳(74)专利代理机构四川省成都市天策商标专利事务所51213代理人刘兴亮吴瑞芳(51)Int.Cl.G06F17/50(2006.01)权利要求书4页说明书10页附图2页(54)发明名称一种高效的涡轮叶片内部冷却结构优化方法(57)摘要本发明公开了一种高效的涡轮叶片内部冷却结构优化方法，包括如下步骤：步骤一、截取C3X叶片中截面，定义设计域；步骤二、获取设计变量；步骤三、坐标转换技术和建模；步骤四、有限元网格划分；步骤五、对流导热系数推导；步骤六、施加边界条件；步骤七、确定目标函数，进行有限元分析；步骤八、对变量进行扰动，记录灵敏度数据；步骤九、优化算法；步骤十、输出结果。本发明通过对C3X叶片内部冷却通道的位置、大小等进行优化，并与遗传算法比较，该优化流程明显提高了优化效率，大大缩短了优化所需的时间，并且该流程能处理应力场和温度场的耦合问题。CN108595797ACN108595797A权利要求书1/4页1.一种高效的涡轮叶片内部冷却结构优化方法，包括如下步骤：步骤一、截取C3X叶片中截面，定义设计域采用基于C3X叶片10个孔的基本拓扑结构，对设计域的基本定义条件如下表所示：表1基本设计域条件为了满足叶片边界和冷却孔边界的最小距离，将叶片边界沿边界内法向方向趋近1mm，将中心区域按照原始10孔的位置排布进行最大限度的分割，分为10块，并且尽可能的拟合叶片边界形状；第一圆孔和第二圆孔的设计域采用圆形区域，用极坐标形式来表示位置；第三圆孔～第十圆孔周围边界近似矩形，故采用矩形区域，用笛卡尔直角坐标来表示；步骤二、获取设计变量整个设计域由2个圆形和8个矩形构成，对于前两个圆域，记录其圆心位置和直径；对于矩形域，选取每个矩形左下角的点作为特征点，记录位置坐标、长、宽以及倾斜角；一个圆在空间中唯一确定需要3个变量，x，y坐标和直径，故10个冷却孔共有30个设计变量；用ρi，θi，di分别表示第1、2个孔的极坐标直径和角度、孔径；用xi，yi，di分别表示第3～10个孔的横坐标，纵坐标和直径；孔的编号顺序和原始构型相同，i∈[1,2,3…10]；步骤三、坐标转换技术和建模采用坐标转换技术，将设计域由全局坐标经过坐标系的旋转平移用与坐标轴正交的矩形域来表示，对于第一圆孔和第二圆孔，将局部坐标建立在圆域中心，由于是极坐标表示，故只需进行一次坐标平移；而对于其余矩形域，将局部坐标系的x，y轴分别与矩形的两边平行，原点建立在矩形域的特征点，从局部到全局需要进行坐标的平移和旋转；式中，x，y表示全局坐标系下的坐标；xi，yi表示在第i个冷却孔局部坐标系下的坐标；xfi，yfi表示第i个冷却孔局部坐标系的特征点坐标，θfi表示矩阵底边与全局x坐标的夹角，用来定义旋转矩阵；T为平移矩阵；R为旋转矩阵；2CN108595797A权利要求书2/4页通过该式变换后，局部坐标点将被转换成全局坐标系下点，可以在局部坐标系下规定设计变量的定常取值范围，表示出全局坐标系下几何矩阵的数学表达形式；步骤四、有限元网格划分对分析所必须的某材料属性进行赋值和规定，材料定义好后，将建好模的平面进行有限元划分，选取PLANE13耦合单元，单元尺寸设为1；步骤五、对流导热系数推导首先，对问题做出3点假设：假设气体采用常值空气的热物理性质；假设流体流量不变；假设各冷却气体的温度不变；冷却通道的对流换热系数可以看作是努赛尔数的函数:式中，λ为流体传热系数，d为冷却通道的水力直径；对于努赛尔数，在管中强制对流近似完全发展阶段的经验公式是:Nu＝0.023Re0.8Pr0.4式中，Pr和Re分别代表普朗特数和雷诺数，普朗特数取常值0.7，雷诺数用下式表示:式中，ρ代表材料密度，v代表流体流速，d代表冷却孔水力直径，μ代表粘度:将数值普朗特数和π代入式中，将表达式简化为h＝0.0242·C(λ,M,μ)d-1.8式中，M代表冷却剂流量；从该表达式中容易看出，对流换热系数h最后简化为由四个量决定的函数，分别是传热系数λ，动力粘度μ，流量M和水力直径d，前两个量是流体的性质，仅随温度变化，因为各冷却孔的温度固定，故对每一个冷却孔来说，两者固定；流量M在假设中也不发生变化，仅由具体孔来决定；实际优化中，完全未知的是水力直径d，把对流换热系数看作是自变量为直径d的表达式，而C(λ,M,μ)看作变化系数，对于一个确定的孔来说，C为一个常数；3CN