基于保结构方法的航天器中微电子器件布局优化*

蒋睿嵩 徐萌波 胡伟鹏†

（1.四川大学 空天科学与工程学院，成都 610065）（2.西安理工大学 土木建筑工程学院，西安 710048）

引言

航天器中的精密电子元器件容易受到航天器工作过程中的过大惯性力作用而影响其工作性能.特别是航天器的局部振动特性，如：共振等，甚至会造成特定位置的电子元器件失效，导致航天任务失败.因此，航天器中的电子元器件往往会被集中布局在航天器某一局部空间（如图1）[1].

图1 航天器中电子元器件封装示意图Fig.1 Layout of electron components in spacecraft

已有研究表明[1，2]，航天器中电子元器件的不同布局位置在一定程度上影响着航天器局部动力学特性，因此，航天器中的电子元器件布局优化问题是航天器设计的重要内容之一.

在布局优化方面，Liu等[1]将电子元器件的布局优化简化为一个NP问题，为本文发展微电子器件保结构优化算法提供了重要借鉴.陶鸿飞和崔升[3]通过对压电梁、板结构振动主动控制进行了分析，实现了压电执行器粘附于悬臂梁上的最佳布局.Jia等[4]考虑土体与结构的相互作用，利用非连续布局优化（discontinuity layout optimization，简称DLO）方法实现了块状结构的稳定性分析.Nanayakkara等[5]针对张拉结构拓扑优化问题，给出了一种简单的布局优化数值方法.Ambrozkiewicz和Kriegesmann[6]针对机械部件及关节设计问题，发展了能够同时实现结构拓扑和布局优化的数值方法.

上述布局优化方法虽然解决了部分工程问题中的布局优化难题，但是，对于航天器中的电子元器件布局优化问题，上述布局优化方法并不适用，原因是航天器中电子元器件的布局优化的目标是使得电子元器件布置位置的振动特性满足电子元器件工作条件的要求，关注的是航天器局部振动特性.无限维动力学系统的保结构分析方法[7-12]的优势正是能够很好地重现动力学系统局部动力学行为，为航天器中的电子元器件布局优化提供了新的途径.无限维系统的保结构分析方法源于冯康先生针对有限维系统建立的辛几何算法[13]和Bridges等针对无限维Hamilton系统建立的多辛算法[14，15]，已被广泛应用于航天动力学问题分析过程中.

本文将针对航天器中微电子器件的布局优化问题，通过建立微电子器件与柔性支撑薄板耦合动力学模型，发展耦合动力学模型的保结构分析方法，实现该耦合动力学问题的保结构分析，在此基础上，对电子元器件的布局进行优化，为航天器局部结构设计提供参考.

1 微电子器件与柔性薄板耦合动力学模型

本节考虑若干个微电子器件粘接在柔性薄板上（如图2），在计入微电子器件质量、忽略微电子器件尺寸假设条件下建立其动力学控制方程.如图2所示，假定微电子器件尺寸远小于圆形柔性薄板的平面尺寸，则微电子器件的尺寸可以忽略不计，其对结构的作用简化为质点，即体现在薄板面密度参数ρ（x，y）中.柔性薄板通过8个螺栓与航天器连接，忽略螺栓与航天器连接间隙，航天器结构通过8个螺栓传递载荷至薄板，引起薄板振动.

图2 微电子器件与柔性薄板简化耦合模型Fig.2 Simplified dynamic model of soft panel assembled with microelectronic devices

2 耦合动力学模型的保结构分析方法

其中，ρ0为未粘接微电子器件的柔性薄板的面密度， ml（l=1，2，…，n） 为柔性薄板区域 Ω ={（i，k） |iΔx≤xl≤（i+1）Δx，kΔy≤yl≤（k+1）Δy}内微电子器件的质量，（xl，yl）为第l个微电子器件在圆形薄板上的坐标值.值得注意的是，由于薄板面内网格尺寸较小，实际的微电子器件之间由于散热等要求，可以假定一个网格内至多能粘接一个微电子器件，这一假定将会具体体现在优化问题建模过程中.

3 微电子器件布局优化模型及数值算例

考虑n个微电子器件粘接在同时受8个冲击载荷作用的柔性薄板上的情形，以n个微电子器件粘结处薄板振动最大加速度的加权平均值最小为优化目标，建立如下微电子器件布局优化问题：

上式中，wl为第l个微电子器件的权重（描述该电子器件对精度和稳定性要求）， ∂ttu（jΔt，xl，yl） 为第l个微电子器件粘结处薄板在jΔt时刻的面外振动加速度，该加速度值由前述建立的保结构分析方法得到.

由于本文考虑的微电子器件质量较小，对薄板面外振动影响只存在于微电子器件粘结点附近的局部区域，因此，参考已有的优化理论及方法[18，19]，本文采用的布局优化算法主要包括如下两个步骤：

在数值算例中，系统的结构及材料参数假定为： c=0.05，R=0.5m，E=2.914 ×109Pa，νp=0.37，h=0.002m，ρ0=1kg/m2，ˉFζ=100N.网格剖分步长及时间步长分别取为：Δx=Δy=0.002m，Δt=0.001s.

图3 微电子器件加速度最大值加权平均值演化过程Fig.3 Evolution of the weighted average of the maximum accelerated speed of the microelectronic devices

从图3可以看出，微电子器件个数越少，布局优化需要的迭代次数也越少.以a*相对误差小于1×10-6为迭代终止条件，当n=2时，需要的迭代次数为106；当n=3时，需要的迭代次数为138；当n=4时，需要的迭代次数为165，迭代终止时，各种情形下各微电子器件位置坐标见表1（注：由于结构的对称性，微电子器件布局优化结果并非唯一，表中只是给出了其中的一种优化结果）.

表1 微电子器件布局优化结果Table 1 Layout optimization results of the microelectronic devices

4 结论

微电子器件在航天器中的布局与航天器局部动力学行为息息相关，对微电子器件的布局进行优化，不仅有利于航天器局部减振，也有利于保证微电子器件的工作性能.本文正是基于以上背景，建立了粘结有若干个微电子器件的柔性薄板在冲击荷载序列作用下的振动方程，基于广义多辛理论，发展了该柔性薄板振动问题的保结构分析方法对该薄板振动问题进行模拟.基于模拟结果，以微电子器件加速度最大值的加权平均值最小为优化目标，对微电子器件的布局位置进行优化，得到了满足约束条件的布局优化结果，为航天器中的微电子器件的布局优化设计提供了新的途径.

需要说明的是，为方便后续的动力学分析和优化设计，本文在薄板振动问题的动力学建模过程中，忽略了薄板厚度、微电子器件尺寸等因素，在布局优化设计过程中，不区分微电子器件的质量差别和权重差别，并仅让微电子器件布置位置在初始位置附近的一个很小范围内摄动，这些简化虽然能够加快优化收敛速度，但是不能从理论上证明布局优化结果的全局最优性.