位移约束下火箭橇结构拓扑优化方法研究摘要：兼顾平安性与经济性的高性能结构优化设计是结构设计领域的重要需求，面向火箭橇结构的拓扑优化设计对火箭橇试验具有重要意义。本文在火箭橇结构动力学仿真分析的基础上，基于变密度法建立了火箭橇结构的拓扑优化模型，利用伴随矢量求解了目标函数与约束函数关于设计变量的灵敏度，通过移动渐近线算法驱动了设计变量的更新迭代，最后对拓扑优化设计的结果进行了校验。结果说明，提出的位移约束下的火箭橇结构拓扑优化设计方法在经济性与平安性方面对火箭橇结构的设计是有效的，拓展了火箭橇结构的设计思路。关键词：火箭橇；动力学分析；变密度法；结构设计；拓扑优化中图分类号：TH122文献标识码：ADOI：10.19452/j.issnl007-5453.2022.02.015火箭橇是指采用火箭彝动机作为动力，推动载有试验件的火箭沿专门建造的高精度轨道高速滑行的地面设备，可以较好地解决各类飞行器在研制过程中有关高速度及高加速可能带来的诸多技术问题［1-3］,为军队提供强大的体系化空中力量，是国防军队现代化的重要表达［4］。火箭橇与轨道的动态性能直接影响到车载试验件的试验环境，将直接决定动态跑车试验结果的准确性和可靠性［5-6］。国内外一些学者对高速火箭橇系统动力学方面开展了研究，Gerasimov等［7］对橇轨之间的橇-轨接触变形及运动稳定性进行了分析；Laird等［8］对高速火箭橇橇-轨撞击过程进行了仿真分析；王健等［9］对火箭橇轨道不平顺功率谱密度进行了分析；顾凯旋等［10］提出了基于梁单元等效的全时段火箭橇动力学分析方法。随着航空零部件研制工作的日新月异,其设计思路也在不断拓宽和自由化［11］,拓扑优化技术可用于结构概念设计阶段，以获得具有更合理材料分布的启发式设计［12］,这已成为高性能工程结构设计领域最活跃课题之一。结构拓扑优化常用方法主要有变密度法(SIMP)［13-14］,水平集法(LSM)［15］和渐进优化法(ESO)［16］等。尽管历经30多年的研究开展，结构拓扑优化技术已经逐渐被众多工程领域所重视，但大局部研究成果只停留在学术层面，少数较为成熟的应用在构件级、简单工况的层面［17］。在整体复杂结构层面，在火箭橇车体设计中运用拓扑优化技术的研究还比拟缺乏。为满足高性能火箭橇复杂车体设计的要求,本文拟建立位移约束下火箭橇结构的拓扑优化模型，并且对拓扑优化设计结果进行校验。首先，在火箭橇动力学仿真分析的基础上，建立基于SIMP法的火箭橇结构拓扑优化模型；其次，基于移动渐近线法(MMA)求解位移约束下的火箭橇结构的拓扑优化问题；最后，对拓扑优化设计结果建立含细节的火箭橇动力学分析的有限元模型，对设计结果进行平安性校验［18］。1火箭橇动力学仿真分析高精度轨道在加工、安装、使用过程中造成的几何不平顺是引起火箭橇与轨道结构振动的重要因素。本节将介绍考虑静态几何不平顺引起的附加接触力火箭橇结构动力学仿真方法。3数值算例考虑橇车结构的服役环境，橇车结构需要与试验件、滑块等相连接，因此需要建立非设计单元。单轨三滑块火箭橇模型的设计域及非设计域如图2所示，其中设计域为蓝色单元，非设计域为紫色单元。有限元模型弹性模量为206GPa,泊松比为0.3o拓扑优化经过185次迭代终止，整个拓扑优化过程花费15min,最终获得的拓扑构型如图3所示，在保证位移滑块处不超过0.5mm的条件下实现减重超过60%o拓扑优化过程中目标函数的变化曲线如图4所示。约束函数随迭代过程的变化曲线如图5所示。针对拓扑优化的结果，建立基于梁单元等效的含滑块细节的有限元模型如图6所示。其中，同一种颜色代表同一种属性。梁单元的信息见表1,其中截面面积1代表梁两端截面面积中较大的一个，截面面积2代表两端截面面积较小的一个，由于火箭橇结构尾端需要与推进装备连接，因此建立偏置梁单元模型。基于第1节介绍的火箭橇动力学仿真分析方法对根据火箭橇拓扑优化设计结果建立的等效梁单元模型进行动力学校验，得到滑块处的动力学响应曲线，如图7〜图9所示，本文仅给出了其中一个滑块的结果，其余两个滑块的响应趋势和响应量级基本相同，三个滑块处的最大均方根加速度和最大应力见表2O由于滑车速度不断加快，单位时...