[发明专利]基于流固耦合的自锁阻尼叶片气流弯应力设计方法

摘要

本发明公开了一种基于流固耦合的自锁阻尼叶片气流弯应力设计方法，其特征在于将自锁阻尼叶片与周围流体作为一个三维耦合系统，设计一套面向节点坐标一一对应的迭代方法，通过耦合求解叶片变形和流场得到叶片的气流弯应力。首先采用Turbogrid软件建立流体计算区域，将网格输入到CFX计算该叶片周围流场分布；采用Solidworks软件建立固体计算区域，将CFX计算得到的叶片表面压力通过节点坐标映射到固体计算网格上计算，在ANSYS中计算得到的气流弯应力，通过安全倍率校核准来判断是否满足要求。本发明所得到的一种基于流固耦合的自锁阻尼叶片气流弯应力设计方法，能更真实反应叶片气流弯应力分布情况，有助于优化叶片设计，减少叶片试验费用，适合于大范围推广。

主权项

一种基于流固耦合的自锁阻尼叶片气流弯应力设计方法，其特征在于：（1）将自锁阻尼叶片与周围流体作为一个三维耦合系统，设计一套面向节点坐标一一对应的迭代方法，并通过耦合求解叶片变形和流场得到叶片的气流弯应力；（2）根据自锁阻尼叶片的实际尺寸要求，采用Turbogrid软件建立流体计算区域，将网格输入到CFX计算该叶片周围的流场分布；（3）采用Solidworks软件建立固体计算区域，将CFX计算得到的叶片表面压力通过节点坐标一一映射到固体计算网格上，进行流固耦合计算，并在ANSYS中计算得到带预应力的气流弯应力；（4）通过安全倍率校核准则，判断气流弯应力分布是否满足设计要求，然后调整叶片叶型参数，循环计算得到所要设计的自锁阻尼叶片；计算详细步骤包括以下步骤：步骤一，几何物理模型的建立；根据初步设计的自锁阻尼叶片参数，先采用Turbogrid软件建立流体域，再采用Solidworks软件建立固体域，并对叶片围带与凸台进行旋转周期性分割；步骤二，网格划分；对流体域采用Turbogrid进行六面体结构化网格划分，对固体域采用Ansysmeshing进行六面体四面体混合网格划分；步骤三，流场计算；将流体网格导入CFX，并设置以下边界条件：动域采用多参考系方法处理，叶片采用周期性边界，入口给定总压、总温、干度，出口给定静压，流体采用IF97水蒸气模型，设置湍流模型RNGΚ−ε]]>，流体控制方程采用有限体积方法求解，选用二阶空间离散格式；步骤四，固体静应力计算；将固体网格导入ANSYS，采用以下循环对称周期边界条件：固定叶根接触面、给定转速载荷，然后再采用隐式动力分析方法，固体域控制方程求解采用有限元方法离散，计算离心力分布；步骤五，流固耦合数值计算；在步骤四叶片承受离心预应力的基础上，将流体计算的叶片表面压力通过节点坐标一一映射到固体域的叶片表面上，采用迭代耦合分析方法，控制方程采用完全牛顿迭代方法，计算带气流弯应力的叶片应力分布；步骤六，气流弯应力后处理；将步骤四与步骤五计算得到的叶片应力分布输出到Excel做差值运算，其结果即为叶片在承受预应力状态下的气流弯应力分布；步骤七，安全倍率校核；根据叶片材料复合疲劳强度曲线，查到叶片在平均应力下的耐振强度，且安全倍率的计算公式为Ab=k⋅σa∗σqw∗=k⋅耐振强度气流弯应力]]>，根据叶片振动强度安全准则评判叶片的性能是否满足；如果不满足，则重新修改叶片型线参数，返回第一步，循环迭代，直到获得满足设计要求的自锁阻尼叶片。