复合材料无人机机翼布局的二级优化方法.pdf

  • 文档大小:354.3 KB
  • 文档格式:pdf
  • 约 5页
  • 2022-04-11 发布
  • 举报
复合材料无人机机翼布局的二级优化方法1 复合材料无人机机翼布局的二级优化方法2 复合材料无人机机翼布局的二级优化方法3 复合材料无人机机翼布局的二级优化方法4 复合材料无人机机翼布局的二级优化方法5
已阅读完毕,您还可以下载文档进行保存
文档单价:1.00 会员免费
开通会员可免费下载任意文档
  1. 1、本文档共5页,内容下载后可编辑。
  2. 2、本文档内容版权归属内容提供方,所产生的收益全部归内容提供方所有。如果您对本文有版权争议,可选择认领。
  3. 3、本文档由用户上传,本站不保证质量和数量令人满意,可能有诸多瑕疵,付费之前,请仔细先通过免费阅读内容等途径辨别内容交易风险。如存在严重挂羊头卖狗肉之情形,可联系本站下载客服投诉处理。
2016年第4期玻璃钢/复合材料31复合材料无人机机翼布局的二级优化方法刘波,陆振玉,袁东明,马经纬(中国科学院长春光学精密机械与物理研究所,长春130033)摘要:为减轻复合材料无人机机翼的结构质量,利用MSC.PATRAN和MSC.NASTRAN建立大展弦比机翼结构布局优化设计的二级优化方法:第一级以机翼的刚度为目标,采用响应面法对翼梁位置进行优化;第二级以机翼结构质量为目标函数,采用遗传算法对机翼各元件的铺层参数进行优化。通过对某型大展弦比无人机机翼进行结构布局优化设计,结果表明提出的大展弦比无人机机翼二级优化方法能够在满足强度、刚度性能设计要求的前提下,减轻约25%的结构质量,减重效果明显。关键词:复合材料;无人机;机翼;二级优化;响应面法;遗传算法中图分类号:TB332文献标识码:A文章编号:1003-0999(2016)04-0031-05无人机独有的低成本、可重复使用和机动能力强等优点,使其应用范围已扩展到军事和民用领域,具有广阔的市场应用前景¨J。一方面,无人机结构通常需要在满足强度和刚度的前提下,重量越轻越好,而复合材料具有高比强度和高比刚度等优点,是理想的飞机结构材料,在航空航天工业中发挥着越来越重要的作用,甚至已成为小型无人机的主体材料¨‘2。另一方面,对于大展弦比无人机机翼,其结构规模大,设计变量多,影响因素复杂,使得优化设计工作十分困难J。因此开展大展弦比复合材料无人机机翼结构的布局优化研究极为重要。近年来,国内外学者在机翼结构参数化建模、机翼布局优化和复合材料铺层优化方面开展了广泛的研究。罗明强等开发了一种飞机机翼结构有限元快速建模和自动化调整的方法,并通过试验验证了方法的有效性。舒恪晟等基于PATRAN的PCL语言开发了机身结构参数化建模的方法。韩伟等¨”为提高优化计算效率,利用响应面法对机翼结构布局进行优化,结果表明采用响应面法能明显提高优化效率。Tuan等将遗传算法用于复合材料层合板铺层优化设计,有效减轻了层合板结构重量。罗志军等提出了一种用于组合优化的基因编码方式和遗传算子,并将其用于层合板的铺层优化设计。晏飞等_l通过分析交叉和突变概率对遗传算法收敛性的影响,提出了一种自适应遗传算法,明显提高了算法收敛速率。丁玲等¨副提出了一种整数编码的遗传算法,并将其用于无人机机翼蒙皮的铺收稿日期:基金项目:作者简介:通讯作者:层优化设计,最后通过有限元分析和静力试验验证了优化的合理性。但这些方法均未能解决大展弦比机翼优化迭代收敛效率低、可靠性差方面的问题。基于上述考虑,本文借助PATRAN与NASTRAN的参数化建模与分析功能,建立二级机翼布局优化方法,最终实现复合材料大展弦比无人机机翼的布局优化设计。1机翼结构布局优化问题描述机翼主要由翼肋、翼梁以及蒙皮壁板等结构组成。为了尽可能提升无人机的各种性能,提高结构材料的利用率,结构优化设计应在各种载荷工况满足强度、刚度性能的要求下使结构的重量最轻,其优化数学模型可描述为:fminc(x1,X2)l—Is.t.Sl,caso…(I,)一S1.<0,case=1,,nl—jS2,sase…(1,)一S2.case<o,case1,,凡(1)…llFl一(,)一1<0,case=1,,nJ矸,≤≤…X】.XUi=1,,,≤≤…【霹JJ戳=1,,t,其中,G为机翼结构质量;case为结构承受载的工况;5为抗弯刚度指标,本文选用机翼翼梢处节点的最大位移来表征;.S为抗扭刚度指标,本文选用机翼翼展方向距翼梢75%处两侧位移差来表征,位移差越小,扭转刚度越大;X和。分别为梁的位置参数和机翼结构铺层参数;F为机翼结构失效因子,根据2015一l1-27中国科学院国防科技创新重点部署项目(Y45137FI40)刘波(1977一),男,博士,副研究员,主要研究方向为微机械系统研制、电机结构设计与分析。陆振玉(1989-),男,硕士,研究实习员,主要研究方向为复合材料结构设计,luzhenyucumt@sina.corn。骥0峨o32复合材料无人机机翼布局的二级优化方法2016年4月蔡一吴失效准则,失效因子的计算公式为:F=1√去叩z㈨+c去一)+c专一1)z+12式中,X、X为材料沿纤维方向拉伸和压缩极限强度;yf、Yc为材料沿垂直于纤维方向拉伸和压缩极限强度;s为材料在纤维铺层平面内的剪切极限强度;or、or:为材料纤维方向和垂直于纤维方向的主应力;"J'12为作用在垂直于纤维方向且沿垂直于纤维方向的剪应力。2机翼布局优化的二级优化方法大展弦比机翼布局优化问题实际上是包含形状和尺寸两类设计变量的整体布局优化设计。其中尺寸设计变量较多,而形状设计变量偏少。目前工程上处理这种问题主要有整体求解和分解处理策略,而由于复合材料设计变量多,使整体求解困难,甚至存在计算上的病态。基于上述考虑,本文利用分级优化的思想,将复合材料机翼优化模型分解成位置优化和尺寸优化:首先,利用PATRAN建立参数化模型;其次,利用响应面法建立近似模型,从而实现翼梁位置的优化;再次,在位置优化的基础上,采用自适应遗传算法对复合材料机翼结构的铺层顺序进行优化;最后,实现大展弦比机翼的布局优化设计,整个优化流程见图1。结果是I鉴壁窭錾妻lh匿蒸辔2.1参数化建模利用MSC.PATRAN的PCL语言可建立复合材料无人机机翼结构的参数化模型,模型参数包括初一姆c艟遮奄始位置参数和初始铺层参数:其中初始位置参数是指翼梁两端在翼弦上的百分比;初始铺层参数包括上、下蒙皮壁板的铺层,梁、墙缘条和腹板的铺层以及翼肋的铺层。机翼各结构全部采用四节点和三节点壳单元进行模拟。2.2基于响应面法的翼梁位置优化在第一级优化阶段,由于翼肋对结构的整体刚度贡献不大,因此应重点考虑梁的位置变量,在该级优化过程中初始铺层参数变量始终保持不变,优化模型可简化为:Imin(Sllc+wS2,。)…ccase:1j≤【.t.,X¨≤Xv¨…i=1,,,其中,W为权重系数,用于保证抗弯刚度指标和抗扭截面指标的数量级相同,取为100。由于该级优化目标无法用优化变量直接表述,因此可采用响应面法建立目标函数的显示关系。响应面法是一种通过构造一个具有明确表达形式的多项式来近似表达隐式功能函数的方法。多项式响应面是应用比较广泛的一种近似模型,工程中通常采用二次多项式响应面近似模型,其基本形式为:∑∑∑∑Y=c。+c+Ci+co.+(4)式中,Y为功能函数实际值;为设计变量;8为响应误差;n为设计变量的维数。响应面近似模型的构建精度主要取决于设计样本点的选取,因此通过使用试验设计方法构造样本点的设计矩阵可以在满足一定精度要求下大大减少计算量,本文采用拉丁超立方法来构造样本点设计矩阵。2.3基于遗传算法的铺层优化设计在第二级优化阶段,因制造工艺性要求,复合材料铺层角常由0O/45。/90。等标准角度组成,因此设计变量均为离散变量。在该级优化过程中,初始位置参数变量保持不变,优化模型可表示为:fminG(X1,X2)l—ls.t.S],c.se…(1,)一S1.c<o,case=1,,凡…{52一(,x2)一_21c<o,case=1,,(5)…llFl.(1,X2)一1<o,case=1,,n√≤≤…lX2J磺=1,,t,对于上述离散变量的优化模型,可利用整数编码的遗传算法进行计算,在计算过程中需要考虑以2016年第4期玻璃钢/复合材料33下几方面:适应度函数的选取、染色体的编码、遗传算子设计、种群规模设置等。(1)适应度函数在应用遗传算法对上述优化模型进行计算前,必须先用罚函数法将约束优化问题变为无约束优化问题,并定义问题的适应度函数:Fitness=c一;lFI一(,)<1且SI,c<S2'ca且是,。<S2,。6一一()0;其他式中,c为一个较大的正数,用来保证适应度函数永不为负值,本文取为100;G。为机翼结构的铺层,均为最大铺层数时的结构质量。(2)编码策略根据复合材料的特点,用定长基因串表示层合板的铺层,采用整数编码的方式。例如:若蒙皮铺层顺序为[空-45。/0。/45。/90。]的层合板,用编码表示为[0/1/2/3/4]。(3)遗传算子遗传算子包括选择、交叉、变异三个基本遗传算子。选择采用非线性排序选择算子;交叉采用改进的自适应两点交叉算子;变异采用改进的自适应变异算子,包括增加、删除、改变和换位操作。2.4二级协同优化在机翼布局优化设计中,第一级和第二级优化的变量之间存在耦合,为此可将第二级优化的结果作为第一级优化的初始参数,进入下一轮迭代优化,即通过二级协同优化设计实现机翼的布局优化设计。在迭代过程中,若连续两轮铺层优化后的机翼结构质量不变或在可接受的范围内,则停止迭代。3算例3.1算例模型无人机机翼采用单梁单墙多肋的结构形式,材料采用T300预浸料,预浸料单层厚度为0.15ram,其单层力学性能见表I,蒙皮壁板采用蜂窝夹芯结构,夹芯采用Nomax芳纶纸蜂窝,规格型号为AC.NH一—1.8340。机翼结构有限元模型如图2所示。机翼各结构全部采用四节点和三节点壳单元进行网格划分;机翼接头利用RBE2刚性单元模拟;机翼上下蒙皮壁板受气动力载荷,载荷工况为最大飞行过载。机翼各结构全部采用四节点和三节点壳单元进行模拟。翼梁初始位置参数为[0.25,0.25];机翼初始铺层参数如表2所示。表1T300复合材料性能参数Table1Prope ̄yparametersofT300composite密度E1E2G12S/GPa/GPa/GPa‘/g一/MPa/MPa/MPa/MPa/MPa1348.5O.345.81.6x10-s1470119038.517678注:1/X为纤维方向;2/Y为垂直纤维方向;为拉伸;c为压缩。图2机翼结构有限元模型Fig.2FiniteelementmodelofUAVwing表2机翼结构初始铺层参数Table2ParametersofUAVwingstructurebeforeoptimization结构铺层参数上壁板[-/4fe5/n9gw0/0/450//0o//00//0一/4-54/59/o/0/405//0o//90/0/0/0/405//09/o/45一/405/]0/。下壁板卜。翼肋[一45/45/0/0/90/0/0/45/-45]翼梁[一45/45/0/0/90/0/0/45/-45]墙[一45/45/0/0/90/0/0/45/一45]3.2优化过程与结果3.2.1翼梁位置优化采用拉丁超立方法安排响应面试验时,考虑到各参数变量的非线性关系和计算精度,样本点个数选取2O较适宜,试验数据及结果如表3所示。最后采用最d"--乘法对机翼刚度进行二次拟合,拟合表达式为:一r11‘(S1+wS1)=273+【一80455]ll—+209x12+[1001240]鲰蒸复合材料无人机机翼布局的二级优化方法从而该级优化的数学模型可近似为:发生在翼梁根部的9163单元第七层处。+209xtx2+[1001。嘲LX2≤≤≤≤s.t.0.110.4;0.120.4对上述模型进行计算,得到翼梁的最优位置参数为[0.39,0.1],此时,机翼翼梢处节点的最大位移为130mm,机翼翼展方向距翼梢75%处两侧位移差为1.1mm。袁3拉丁超立方试验设计结果Table3ExperimentresultsofLatinHypercubeSamplingpoints3.2.2机翼铺层优化此阶段的优化模型则变为:maxf71≤√≤…【s.t.,X2J=1,,t,设置遗传算法的运行参数:种群规模为80,最大进化代数为200,交叉概率为0.85,变异概率为0.01,最大铺层数为20。进化迭代过程如图3所示。从图中可以看出,种群进化到约50代即收敛到最优解,此时机翼铺层总的结构质量为56.7kg;最大位移130mm发生在机翼翼梢处,机翼翼展方向距翼梢75%处两侧位移差为1.1mm;最大失效因子为0_9769,啦0秘一图3机翼结构优化迭代过程Fig.3ConvergenceprocessforoptimizationofUAVwingstnlcture3.2.3二级协同优化将第一轮机翼铺层优化结果作为第二轮翼梁位置优化的初始参数,进入第二轮迭代优化,将20个样本点试验数据采用最小二乘法进行拟合,优化的数学模型可近似为:min272+…8+l19.9xlx2+[991.55矧圈≤≤≤≤s.t.0.110.4;0.120.4对上述模型进行计算,得到翼梁的最优位置参数仍为[0.39,0.1],进入第二级优化,优化迭代过程如图3所示,从图中可以看出第二轮优化后结构质量为46.6kg,较第一轮优化结果47.4kg满足收敛条件,迭代结束。优化后机翼结构质量较优化前62.4kg减重约25%,减重效果明显,优化后机翼各结构铺层如表4所示。优化前后机翼结构的位移云图如图4所示。表4机翼结构优化后各结构铺层参数Table4Parametersofwingstructureafteroptimization结构铺层参数上壁板:下壁板/fe4翼肋[-45/90/45/0/0/0/45/90/一45]翼梁[-45/90/45/0/0/45/90/一45]墙__45/90/45/0/0/45/90/一45]2016年第4期玻璃钢/复合材料35(a)Deformationofwingbe ̄reoptimization(b)Deformationofwingafteroptimization图4优化前后机翼变形云图Fig.4Deformationnephogramofwingbeforeoptimizationandafteroptimization算例结果表明提出的大展弦比无人机机翼二级优化方法能够在满足强度、刚度性能设计要求的前提下明显减轻机翼结构重量。4结论本文根据工程实际需要,提出一种大展弦比机翼结构布局优化设计的二级优化方法,该方法具有以下两个特点:(i)考虑了分级优化变量之间的耦合关系,在优化过程中,将第二级优化的结果作为下一轮第一级优化的初始参数,循环迭代计算;(2)考虑了多种工况,使得所计算的结果满足多种工况的设计要求,结果可行性高。参考文献[1]崔秀敏,王维军,方振平.小型无人机发展现状及其相关问题分—析[J].飞行力学,2005,23(1):1418.[2]杨乃斌,章怡宁.复合材料飞机结构设计[M].北京:航空工业出版社,2002:2-10.[3]薛忠民.中国玻璃钢/复合材料发展回顾与展望[J].玻璃钢/复—合材料,2015(1):512.[4]杜善义.先进复合材料与航空航天[J].复合材料学报,2007,—24(1):112.[5]RicheRL,HaftkarRT.Optimizationoflaminatestackingsequenceforbucklingloadmaximizationbygeneticalgorithm[J].AIAA,1993,3I:951-956.[6]梁猛,万志强.轻质双尾撑布局无人机气动弹性建模与分析[R].CARS,2007:25-28.[7]尹星研,冯振宇,卢翔.基于MSC.Nastran的无人机复合材料机翼有限元分析[J].玻璃钢/复合材料,2010(1):3-6.[8]王伟,杨伟,常楠.大展弦比飞翼结构形状、尺寸综合优化设计—[J].强度与环境,2007,34(5):4957.[9]罗明强,冯昊成,刘虎,等.机翼结构有限元的快速建模及自动—化调整[J].北京航空航天大学学报,2011,37(6):680684.[1O]舒恪晟,金海波,姜琬.基于Patran二次开发的机身结构参数化建模[J].飞机设计,2011,31(2):30-32.[11]韩伟,何景武.响应面法在机翼结构优化中的应用研究[J].飞机设计,2013,33(5):6-11.[12]TuanLM,JaehongL.Stackingsequenseoptimizationformaximumstrengthsoflaminatedcompositeplatesusinggeneticalgorithmandisogeometricanalysis[J].CompositeStructures,2014,116:357-363.[13]罗志军,乔新.基于遗传算法的复合材料层压板固有频率的铺层顺序优化[J].复合材料学报,1997,14(4):79-83.[14]晏飞,李为吉.基于自适应遗传算法的复合材料层合板铺层顺—序优化设计[J].西北工业大学学报,2001,19(1):156159.[15]丁玲,孙辉,贾宏光,等.应用遗传算法优化设计机翼复合材料蜂窝夹层结构蒙皮[J].光学精密工程,2014,22(12):3272.3279.LAYoUToPTIMIZATIONOFCOMPOSITEUAVWINGSTRUCTUREBASEDⅡoNTWo.LEVELoPTIZATIoNMETHoD———LIUBo,LUZhenyu,YUANDongming,MAJingwei(ChangchunInstituteofOptics,FineMechanicsandPhysics,ChineseAcademyofScience,Changchun130033,China)Abstract:Inordertoreducetheweightofwingstructureofcompositeunmannedaerialvehicle(UAV),amethodoftwo.1eve1optimizationwhichusesMSC.PATRANandMSC.NASTRANisproposedinthispaper.Inthefirstlevel,aresponsesurfacemethodisusedtooptimizethelocationofthelongitudinalwingspars.Inthesecond‘level,ageneticalgorithmisusedtooptimizethelayparameterofwingstructure.Asanexample,alayoutoptimiza—tionofonetypeofhighaspectratiowingstructureiscalculated.Theresultshowsthatabout25%weightofwingstructureisreduced.Itcanbeconcludedthatthetwo-leveloptimizationmethodcansignificantlyreducethewingstructureweight,whichcansatisfythedesignrequirements.—Keywords:composite;UAV;wing;twoleveloptimization;responsesurface;geneticalgorithm哪瓣■■渤r『闲
文库帮主
该用户很懒,什么也没介绍
文档单价:1.00 会员免费
开通会员可免费下载任意文档