一种燃气轮机风扇叶片反扭过程的预估方法技术

本发明专利技术提供一种燃气轮机风扇叶片反扭过程的预估方法，目的是为了解决目前燃气轮机风扇叶片反扭设计方法不完善的问题，该方法对热态叶型进行ANSYS离心力静态变形分析，预测冷态叶型，再对其进行ANSYS离心力静态变形分析，过渡虚热态叶型，修正固体和流体网格，气动力变形的流固耦合分析等步骤最终得到冷态叶型。本发明专利技术为实现各个气动工况下预测的热态叶型与叶片实际工作态叶型的一致性，对精确预估燃气轮机风扇叶片反扭过程的方法进行研究是十分重要的，尤其对于发展大涵道比的跨音速风扇有着重要的工程应用价值。

【技术实现步骤摘要】

本专利技术涉及叶片叶型设计的
，具体涉及。
技术介绍
燃气轮机，特别是涡扇发动机风扇级的性能对发动机总性能有着至关重要的影响，而叶型则是决定风扇性能的最重要的参数之一，所以保证风扇叶片在工作点的形状和设计点的相同或接近是实现其设计性能指标的关键环节。燃气轮机风扇在工作过程中，叶片在气动力、离心力和热应力的共同作用下，会产生伸长、展开等复杂的变形，叶片会从加工的形状(冷态叶型)变为工作转速下的叶型(热态叶型)，而热态叶型又受转速和负载的制约，所以不同气动工况下热态叶型略有不同。在航空发动机的叶片设计中，一般根据气动性能、强度、寿命等要求设计出在设计工况下的叶片最佳气动几何形状(热叶片)，再利用结构分析对设计叶型进行反推以求得冷态加工叶型(热到冷)，即叶片反扭设计。另一方面，许多 涉及到叶片变形方面的研究，如气动弹性的强迫响应、颤振等问题，往往需要叶片在非设计气动工况下的精确叶型。以冷态加工叶型为分析起点，得到工作热态下叶片形状的预估称为“冷到热”的过程。大量研究表明，燃气轮机叶片形状的细微差别会影响叶尖间隙、气流角等，从而导致风扇叶片通道内的流动和设计工况有一定差别。特别对跨音速风扇，叶片变形可能会引起激波位置的移动，从而偏离其最佳设计点，进而影响燃气轮机的性能、效率，甚至带来气动弹性稳定性等问题。而现代发动机高推重比的需求，要求风扇叶片高负载，低重量。新型轻质复合材料、空心结构等措施的采用使得叶片相对刚性下降，在大的气动载荷下，较之常规风扇更易产生较大变形，气动力的非线性效应更明显。而在超、跨音速流动工况下叶片变形对激波位置有重要影响。叶片在气动力作用下的变形作用体现的更突出。而传统叶片反扭设计，只是根据结构强度需求，叶片安装角偏转一个角度。这种方法对新结构、新材料的设计很难得到优化的叶片制造叶型，从而导致叶片在设计转速下工作时与设计叶型不一致。无论是“冷到热”还是“热到冷”的叶型计算，都要求反推的冷态叶型具有较高精度，才能实现各个气动工况下预测的热态叶型与叶片实际工作态叶型的一致性。因此对精确预估燃气轮机风扇叶片反扭过程的方法进行研究是十分重要的，尤其对于发展大涵道比的跨音速风扇有着重要的工程应用价值。
技术实现思路
专利技术目的是为了解决目前燃气轮机风扇叶片反扭设计方法不完善的问题，进而提供了。本专利技术为解决上述技术问题采取的技术方案是，该方法由以下步骤实现的步骤一以设计时叶片几何数据为原始热态叶型，叶片材料为设计时选取的材料，获取弹性常数，泊松比性能参数，以设计点的气动工况为工作点；步骤二 根据原始设计工况下热态叶片几何型值，构建风扇叶片的几何模型；步骤三将风扇叶片进行网格划分，并结合设计时给定的材料性能参数，构建风扇叶片的有限元模型；步骤四结合风扇叶片的设计点运行工况，即转速的设计要求，叶片位移约束边界条件，进行叶片在离心力作用下的静态变形分析，获得离心负载下叶片的静变形量；步骤五使用原始热态叶型减去离心静变形得到固体域中的预测冷态叶型①首先提取叶片各节点的几何坐标即XsO，YsO, ZsO，然后根据离心静变形计算结果，在X、Y、Z轴方向上提取叶片各节点处变形值即UXl，UYl，UZ1，在步骤三中的风扇叶片几何外形的基础上实行叶片初步反扭计算，在X、Y、Z轴方向上使各节点坐标取值为XsO-UXl，YsO-UYl, ZsO-UZl,从而,确定风扇叶片的预测冷态几何形式；步骤六将燃气轮机风扇叶片的几何模型导入计算流体力学分析软件的前处理器AUT0GRID,建立流体域，完成流体网格的划分；步骤七使用原始热态叶型减去离心静变形得到流体域中的预测冷态叶型①首先提取叶片各节点的几何坐标即XfO，YfO, ZfO，然后使用所述有限元模型得到节点上的静变形量UX1，UY1，UZl插值到叶片的流体网格节点上以获得UXfl，UYfLUZfl ;然后，使用动网格技术，将基于原始叶型下的流体域网格节点以叶片固面上节点变形量为边界条件，设流体域其他边界变形为0，进行调整，并最大程度保证流体网格质量不变；步骤八使用获得的冷态叶型，加载设计转速，材料性能参数，及叶片支撑边界条件构建新的风扇叶片有限元模型；进行叶片在离心力作用下的静态变形分析，获得离心负载下叶片的静变形量；在X、Y、Z轴方向上提取叶片各节点处变形值即UX2，UY2，UZ2 ;完成自由振动模态分析，求解动力学方程，其动力学方程如下Mi/ + Cf/ + Kf/ = 0,丫、其中M为质量矩阵，C为阻尼矩阵，K为刚度矩阵T力加速度，为速度，U为位移矢量；提取风扇叶片多阶固有频率，以及各阶频率下有限元各节点的振动模态值；步骤九在流体分析模型中，使用获得的冷态叶型加上离心静变形得到一个过渡虚热态叶型首先提取叶片各节点的几何坐标即Xf，Yf, Zf ，然后使用所述叶片有限元模型节点上的静变形量UX2，UY2，UZ2插值到叶片的流体网格节点上以获得UXf2，UYf2, UZf2 ；然后，使用动网格技术，以叶片固面上节点为移动边界，流体域其他边界变形为0，并最大程度保证流体网格质量不变，调整基于冷态叶型下的流体域网格节点；步骤十结合风扇叶片的运行工况进行叶片流固耦合的瞬态响应分析，获得气动力作用下的叶片变形量；建立基于RANS方程的流体动力学非定常流动计算模型，如下式—jOdn+ f (Pt + FJdS = j Sc/il⑵ son其中Q是某一物理时刻的控制体的体积、表面面积是守恒型流体状态向量，^为对流通量，为粘性通量d是源项矢量，代表了外力和坐标系运动对守恒律的贡献；瑞流流动的模拟采用了 Sparlart-Allmas —方程瑞流模型；求解RANS方程,从而获取流动速度及作用在风扇叶片上的气动力；将有限元模型上节点的模态值插值到叶片的流体网格点上；组集动力方程所需的风扇叶片固有模态集；使用获得的固有模态集解耦并求解动力学方程，得到广义坐标下N个线性无关的常微分方程组如下，其中N为叶片模态数，MsIs-^CgIfKKjg) = {Oj、12)其中{4}为广义位移向量，是广义速度向量，是广义加速度向量，18、匕和Kg分别是归一化后的广义质量矩阵、广义阻尼矩阵和广义刚度矩阵，{Qg}为广义气动力，采用多步龙格库塔积分法分别求解各阶模态方程(12)，得到各个模态的广义位移，并线性叠加得到叶片的总物理位移；叶片在非定常气动力作用下的变形使用流固耦合方法计算在每个物理时间步，耦合求解流体运动方程和叶片变形的控制方程，并互相作为边界条件进行物理量交换；首先流体域求解了时间精确的非定常雷诺平均N-S方程，以得到每一步由于叶片变形而引起的流场变化，计算气动力作用下叶片表面受到的广义力，并求解广义坐标的结构动力学方程，得到广义位移{(}，再利用模态叠加法得到叶片表面每一点的位移{x};采用投影-形函数插值法处理流体域和固体域的信息交换，然后以叶片外其他边界静止或滑动为约束条件，采用动网格技术，移动流场中所有网格点，以满足给定的叶片变形条件；以变化后网格为初值，迭代求解N-S方程，得到下一时刻的流场；进行以上循环达到收敛性目标为止；通过模态叠加的结构非线性瞬态响应分析，得到稳态响应的值，即为非定常气动力引起的变形，此刻叶型为虚热态叶型；步骤^^一 获取的虚热态叶型，判断其是否与原始本文档来自技高网...

【技术保护点】
一种燃气轮机风扇叶片反扭过程的预估方法，其特征在于：该方法由以下步骤实现的：步骤一：以设计时叶片几何数据为原始热态叶型，叶片材料为设计时选取的材料，获取弹性常数，泊松比性能参数，以设计点的气动工况为工作点；步骤二：根据原始设计工况下热态叶片几何型值，构建风扇叶片的几何模型；步骤三：将风扇叶片进行网格划分，并结合设计时给定的材料性能参数，构建风扇叶片的有限元模型；步骤四：结合风扇叶片的设计点运行工况，即转速的设计要求，叶片位移约束边界条件，进行叶片在离心力作用下的静态变形分析，获得离心负载下叶片的静变形量；步骤五：使用原始热态叶型减去离心静变形得到固体域中的预测冷态叶型①：首先提取叶片各节点的几何坐标即Xs0，Ys0，Zs0，然后根据离心静变形计算结果，在X、Y、Z轴方向上提取叶片各节点处变形值即UX1，UY1，UZ1，在步骤三中的风扇叶片几何外形的基础上实行叶片初步反扭计算，在X、Y、Z轴方向上使各节点坐标取值为Xs0？UX1，Ys0？UY1，Zs0？UZ1，从而，确定风扇叶片的预测冷态几何形式；步骤六：将燃气轮机风扇叶片的几何模型导入计算流体力学分析软件的前处理器AUTOGRID，建立流体域，完成流体网格的划分；步骤七：使用原始热态叶型减去离心静变形得到流体域中的预测冷态叶型①：首先提取叶片各节点的几何坐标即Xf0，Yf0，Zf0，然后使用所述有限元模型得到节点上的静变形量UX1，UY1，UZ1插值到叶片的流体网格节点上以获得UXf1，UYf1，UZf1；然后，使用动网格技术，将基于原始叶型下的流体域网格节点以叶片固面上节点变形量为边界条件，设流体域其他边界变形为0，进行调整，并最大程度保证流体网格质量不变；步骤八：使用获得的冷态叶型，加载设计转速，材料性能参数，及叶片支撑边界条件构建新的风扇叶片有限元模型；进行叶片在离心力作用下的静态变形分析，获得离心负载下叶片的静变形量；在X、Y、Z轴方向上提取叶片各节点处变形值即UX2，UY2，UZ2；完成自由振动模态分析，求解动力学方程，其动力学方程如下：MU··+CU·+KU=0---(1)其中M为质量矩阵，C为阻尼矩阵，K为刚度矩阵，为加速度，为速度，U为位移矢量；提取风扇叶片多阶固有频率，以及各阶频率下有限元各节点的振动模态值；步骤九：在流体分析模型中，使用获得的冷态叶型加上离心静变形得到一个过渡虚热态叶型：首先提取叶片各节点的几何坐标即Xf，Yf，Zf，然后使用所述叶片有限元模型节点上的静变形量UX2，UY2，UZ2插值到叶片的流体网格节点上以获得UXf2，UYf2，UZf2；然后，使用动网格技术，以叶片固面上节点为移动边界，流体域其他边界变形为0，并最大程度保证流体网格质量不变，调整基于冷态叶型下的流体域网格节点；步骤十：结合风扇叶片的运行工况进行叶片流固耦合的瞬态响应分析，获得气动力作用下的叶片变形量；建立基于RANS方程的流体动力学非定常流动计算模型，如下式：其中Ω、是某一物理时刻的控制体的体积、表面面积；是守恒型流体状态向量，为对流通量，为粘性通量；是源项矢量，代表了外力和坐标系运动对守恒律的贡献；湍流流动的模拟采用了Sparlart？Allmas一方程湍流模型；求解RANS方程，从而获取流动速度及作用在风扇叶片上的气动力；将有限元模型上节点的模态值插值到叶片的流体网格点上；组集动力方程所需的风扇叶片固有模态集；使用获得的固有模态集解耦并求解动力学方程，得到广义坐标下N个线性无关的常微分方程组如下，其中N为叶片模态数，Mg{ζ··}+Cg{ζ·}+Kg{ζ}={Qg}---(12)其中{ζ}为广义位移向量，是广义速度向量，是广义加速度向量，Mg、Cg和Kg分别是归一化后的广义质量矩阵、广义阻尼矩阵和广义刚度矩阵，{Qg}为广义气动力，采用多步龙格库塔积分法分别求解各阶模态方程（12），得到各个模态的广义位移，并线性叠加得到叶片的总物理位移；叶片在非定常气动力作用下的变形使用流固耦合方法计算：在每个物理时间步，耦合求解流体运动方程和叶片变形的控制方程，并互相作为边界条件进行物理量交换；首先流体域求解了时间精确的非定常雷诺平均N？S方程，以得到每一步由于叶片变形而引起的流场变化，计算气动力作用下叶片表面受到的广义力，并求解广义坐标的结构动力学方程，得到广义位移{ζ}，再利用模态叠加法得到叶片表面每一点的位移{x}；采用投影？形函...

【技术特征摘要】

【专利技术属性】
技术研发人员：郑赟，杨慧，王彪，方韧，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：