一种宽弦风扇叶片二维有限元建模方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN108090292A(43)申请公布日2018.05.29(21)申请号201711427365.2(22)申请日2017.12.26(71)申请人中国航发四川燃气涡轮研究院地址610500四川省成都市新军路六号(72)发明人潘容孙海鹤米春虎王春健(74)专利代理机构中国航空专利中心11008代理人杜永保(51)Int.Cl.G06F17/50(2006.01)权利要求书1页说明书3页附图2页(54)发明名称一种宽弦风扇叶片二维有限元建模方法(57)摘要本发明属于航空发动机领域，具体涉及到一种宽弦风扇叶片二维有限元建模方法。一种适用于宽弦风扇叶片的二维有限元建模方法，包括如下步骤：1二维有限元网格划分，2三维叶片厚度信息提取，3二维有限元网格厚度映射，4叶片离心力补偿修正。能真实模拟宽弦风扇叶片的三维质量和刚度分布，使二维有限元模型计算得到的风扇轮盘应力符合实际情况；其技术原理清晰，操作方法简便易行，适合于工程应用；还可推广到其余叶片类结构的二维有限元建模，提高发动机零件的强度仿真精度。CN108090292ACN108090292A权利要求书1/1页1.一种适用于宽弦风扇叶片的二维有限元建模方法，其特征在于；包括如下步骤：步骤一、二维有限元网格划分宽弦风扇转子由宽弦风扇叶片(1)和轮盘(2)组成；将宽弦风扇转子的二维结构沿轮缘流道(3)分开，其中宽弦风扇叶片(1)子午投影轮廓的四个角点从叶根前缘(4)按顺时针排列分别为叶根前缘(4)、叶尖前缘(5)、叶尖后缘(6)、叶根后缘(7)，以叶根前缘(4)为原点建立直角坐标系，X轴沿竖直方向，Y轴沿水平方向；以宽弦风扇叶片(1)子午投影轮廓建立子午平面，设定单元长度将宽弦风扇叶片(1)子午平面离散为四边形单元，其中单元长度根据模拟精度的需求设置，单元类型为带厚度的平面应力单元，每一个单元都有唯一的编号和位置信息；可得宽弦风扇叶片(1)子午面二维有限元网格中每个单元Ei的质心I，其质心坐标为(xi,yi)，单元的厚度hi为待设定的值；步骤二、三维叶片厚度信息提取将宽弦风扇转子的三维结构沿轮缘流道(3)分开，同样以叶根前缘(4)为原点建立直角坐标系，X轴沿竖直方向，Y轴沿水平方向；提取宽弦风扇叶片(1)三维结构的叶型截面(8)数据，叶型截面(8)上任一点J的坐标为(Xi，Yi，Zi)，任一叶型截面(8)在YZ平面上的投影如图4，获取任一点J沿Z方向在叶型截面范围内的长度Hi＝Zi1-Zi2，并沿Y轴方向设定步长提取叶型截面厚度；步骤三、二维有限元网格厚度映射设定二维有限元网格中所有叶片单元的质心坐标数组为[x，y]，对应的厚度信息数组为h；三维叶型截面(8)在XY平面上投影的坐标数组为[X，Y]，对应的厚度信息数组为NH，其中N为风扇转子叶片数；通过基于位置的矩阵映射将三维模型[X，Y]位置的厚度信息数组NH映射到二维模型[x，y]处的厚度信息数组h，以此方法将三维叶型的质量和刚度分布映射到二维有限元网格上；步骤四、叶片离心力补偿修正对于上述过程中产生的叶片离心力误差，需进行补偿修正；设宽弦风扇叶片(1)的材料密度为ρ1，三维叶片的体积为V1，质心半径为R1；二维有限元模型中单元Ei的质心半径为xi，通过映射已得到厚度为hi，其单元面积可测量并设为Ai，二维有限元模型的网格总数为n；通过叶片离心力等效原则，得到二维有限元模型的等效密度应为即完成了宽弦风扇叶片(1)的二维有限元建模。2.根据权利要求1所述的一种适用于宽弦风扇叶片的二维有限元建模方法，其特征在于；所述步骤一中单元长度范围为0.1mm～1mm。3.根据权利要求1所述的一种适用于宽弦风扇叶片的二维有限元建模方法，其特征在于；所述步骤二中三维叶片厚度提取步长为0.05mm～0.1mm。2CN108090292A说明书1/3页一种宽弦风扇叶片二维有限元建模方法技术领域[0001]本发明属于航空发动机领域，具体涉及到一种宽弦风扇叶片二维有限元建模方法。背景技术[0002]在航空发动机的设计过程中，转子组件有限元应力分析时经常采用二维模型，而叶片的二维有限元模型通常是采用带厚度的平面应力单元，整个叶片二维模型具有同样的厚度以模拟其质量离心力。这种转子叶片二维模型主要关注总体离心力载荷，并未考虑真实叶片的质量和刚度分布，因而转子二维模型有限元应力分析中，可能导致叶片传递的载荷分布与真实结构产生差异，甚至二维有限元模型计算得到轮盘的应力与真实三维结构出现不同，特别是对于多辐板宽弦风扇整体叶盘结构，这种差异更加明显，使得二维有限元模型计算得到各辐板的应力分布与真实结构不符。发明内容[0003]本发明是针对航空发动机设计中使用二维有限元模型分析宽弦风扇转子应力分析时，