发明名称 |
一种直升机旋翼桨叶剖面设计方法 |
摘要 |
本发明公开了一种基于拓扑优化的直升机旋翼桨叶剖面设计方法,属于直升机技术领域。首先根据气动要求确定外形参数,指定优化目标函数、约束函数,建立优化模型;然后采用有限元软件建立桨叶几何模型并划分网格,将节点、网格和材料信息导入到MATLAB软件中,采用以每一个单元密度为设计变量的SIMP算法;再次求解当前设计状态的旋翼桨叶敏度信息,并采用敏度过滤方法以避免棋盘格问题,更新桨叶剖面单元的密度信息;最后求解更新前后的密度差值,满足收敛条件则结束优化,否则以更新后的密度为初始设计参数继续迭代计算直到收敛,最后输出桨叶结构型式。本发明提出的桨叶设计方法可以应用于概念设计阶段的桨叶结构初步设计。 |
申请公布号 |
CN102902845B |
申请公布日期 |
2015.01.28 |
申请号 |
CN201210337325.X |
申请日期 |
2012.09.12 |
申请人 |
北京航空航天大学 |
发明人 |
向锦武;任毅如;罗漳平;李道春;张连鸿;郭俊贤;张亚军;黄明其;张田田 |
分类号 |
G06F17/50(2006.01)I |
主分类号 |
G06F17/50(2006.01)I |
代理机构 |
北京永创新实专利事务所 11121 |
代理人 |
赵文利 |
主权项 |
一种直升机旋翼桨叶剖面设计方法,其特征在于,具体包括以下几个步骤:第一步,根据直升机旋翼的气动外形参数要求,确定桨叶的气动外形,根据拓扑优化设计要求,确定目标函数、约束函数,进而建立优化模型;第二步,采用有限元软件利用步骤一中确定的气动外形参数建立桨叶几何模型,并划分网格,得到n个单元,导出单元节点坐标、单元信息和材料参数信息;第三步,读入n个单元的节点坐标、单元信息和材料参数信息;第四步,根据当前桨叶剖面各个单元的密度,求解目标函数的敏度信息;采用拓扑优化中的SIMP算法,以每一个单元的密度ρ<sub>i</sub>作为设计变量,利用单元模量信息,求解当前旋翼桨叶目标函数关于各个单元设计变量的敏度信息;第五步,由步骤四得到目标函数关于各个单元设计变量的敏度信息cs<sub>j</sub>之后,采用序列线性规划算法更新每一个单元的设计变量;第六步,原结构第i个单元的密度为ρ<sub>i</sub>,通过第五步得到的更新后的单元密度为ρ<sub>i</sub><sup>new</sup>,计算更新前后每一个单元密度差值的绝对值,当最大的差值小于指定的收敛条件η后,满足收敛条件,结束优化,否则回到第四步,直到满足收敛条件;输出优化得到的各个单元的密度信息,从而得到桨叶结构的拓扑结构;所述的第一步中,以柔度或者截面刚度作为目标函数,约束条件为质量约束,设计变量为每一个单元的密度,建立优化模型如下:目标函数:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><munderover><mo>∫</mo><mi>o</mi><mi>L</mi></munderover><mi>θFθdl</mi></mrow>]]></math><img file="FDA0000608078030000011.GIF" wi="229" he="127" /></maths>或K满足约束条件:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>ρ</mi><mi>i</mi></msub><msub><mi>A</mi><mi>i</mi></msub><mo>≤</mo><msub><mi>m</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><msub><mi>ρ</mi><mi>min</mi></msub><mo>≤</mo><mi>ρ</mi><mo>≤</mo><mn>1</mn></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000608078030000012.GIF" wi="241" he="176" /></maths>其中:L表示桨叶长度,l为沿桨叶长度的变量,C为柔度,K为截面刚度,θ为截面载荷,F为截面柔度矩阵,n为单元个数,ρ<sub>i</sub>为截面第i个单元的密度,单元厚度默认为1,A<sub>i</sub>为第i个单元的面积,i为单元编号且i=1,2......n,m<sub>0</sub>为质量约束,ρ<sub>min</sub>为单元密度的最小值。 |
地址 |
100191 北京市海淀区学院路37号 |