发射动力系统连接结构薄层单元建模的不确定性研究

李 昊，王鸿丽，卢丙举

(1. 海军装备部装备审价中心，北京 100071；

2. 中国船舶重工集团公司第七一三研究所，河南 郑州 450015；

3. 河南省水下智能装备重点实验室，河南 郑州 450015)

0 引 言

在某发射动力系统连接结构中，螺栓连接以其连接刚性好，安装方便等诸多优点得到广泛的应用，但是螺栓连接件会影响结构局部刚度和整机的动力学性能，所以建立准确的动力学连接模型是准确计算结构动态特性与响应的前提。

螺栓连接的有限元建模方法主要分为两类：非参数化建模方法和参数化建模方法。非参数化建模将连接刚度和连接阻尼等效为对结构动力学特性作用相同的附加外力的形式[1-4]，缺点是无法揭示连接处的微观特性。与非参数化建模不同的是参数化建模是根据连接处的物理特性和连接界面的微观滑移、碰撞等特性建立连接模型，因为能够深刻揭示连接的特性而得到越来越多的研究。参数化建模分为零厚度单元（Zero thickness element）、一般单元（Generic element）和薄层单元（Thin layer element）3 种类型。零厚度单元用集中质量、弹簧和摩擦滑块来表示连接界面处的力—位移关系。Iwan[5]模型因为能够表示连接处的粘滑运动而被广泛应用；Ahmadian 和Jalali[6]采用一般单元模拟连接处的特性，并通过实验频响函数对连接参数进行了识别；薄层单元理论最早应用于岩石接触的力学分析[7]，Mayer 和Gaul[8]研究表明薄层单元能够有效模拟机械连接面的接触特性。

在连接的不确定性方面，国内外学者进行了大量的研究。Gangaharan[9]提出一种基于统计学的参数识别方法，运用系统静力学响应来进行参数识别；Qiao LS[10]采用迟滞模型研究了非自治系统的刚度不确定性；Guo Qin tao 和Zhang Ling mi[11]将连接结构的接触刚度及接触阻尼作为随机分布参数处理，然后基于响应面法的模型修正方法和分布代数方法进行参数识别；姜东和费庆国[13]基于薄层单元理论，提出了一种螺栓连接接触面不确定性参数识别方法；C.Zang[13]基于鲁棒性设计研究了两自由度结构的动力学不确定性问题。本文以某发射动力系统连接结构为研究对象，针对薄层单元连接模型，通过Taguchi 方法研究了薄层单元的厚度、弹性模量和密度不确定性。

1 薄层单元基本理论

薄层单元建模方法把2 个连接件接触面的部分等效为连续的厚度非常薄的单元，如图1 所示。

图 1 薄层单元Fig. 1 Thin layer element

对于尺寸为l1×l2×d 的薄层单元，根据虚位移原理得到虚功方程：

对式（2）运用2 节点高斯积分可得到刚度矩阵的数值表达式：

对于薄层单元，假设厚度远远小于长度和宽度，文献[15 - 16] 表明当厚度趋近于零时单元应变可以忽略。此把接触面法向和接触面切向， 分别定义为局部坐标， 和的方向；同理，应力分对应；变，分对应和，则本构方程可以转化为：

式中： 和 分别为弹性模量，剪切模量。此时 和为2 个独立的参数，如果要求两者的耦合关系，可在本构关系式（5）中加入耦合项实现。

2 Taguchi 设计方法

鲁棒设计是在不减弱或消除不确定性因素影响的情况下，寻找合理的设计值，使响应对参数变化不敏感的设计方法。结构动力学鲁棒设计问题属于不确定结构动力学分析的反问题。Taguchi 设计是一种重要的鲁棒设计方法，最早由日本学者田口谷一提出[15]。这种设计方法把统计学的实验设计方法应用到了工业设计中，通过较少的实验设计来减少设计人员的工作量，同时设计出较高质量的产品。Taguchi 方法将工程优化过程和产品设计过程分为3 步：系统设计、参数设计和容差设计，具体设计过程见文献[16]。

Taguchi 利用信噪比考虑结构性能的鲁棒性以及产品的综合经济性。根据产品的质量特性，信噪比分为3 类：望目特性的信噪比，望小特性的信噪比和望大特性的信噪比。

望目特性的信噪比：

望小特性的信噪比：

望大特性的信噪比：

通过信噪比分析，可以初步得到一组参数组合，经过进一步优化分析，最终确定最优化参数组合。

3 薄层单元的鲁棒设计

3.1 设计模型

某发射动力系统连接结构模型如图2 所示，连接件由前后2 个结构组成，2 个结构的法兰边通过薄层单元进行连接。考虑到螺栓连接对连接结构局部刚度影响比较大，而局部刚度的变化对连接结构整体影响非常大。在连接结构的建模中，薄层的厚度对建模的准确性至关重要，对连接结构的整体响应影响巨大。而目前国内外对于如何准确确定薄层厚度尚无定论。本文对连接结构的整体动力学响应为研究目标，以薄层厚度，弹性模量和密度为设计变量，对发射动力系统连接结构进行鲁棒设计。

3.2 信噪比分析

图 2 螺栓连接结构薄层单元示意图Fig. 2 Thin-layer element in thecontact surface of assembled aero-engine casing

表 1 设计参数水平设置Tab. 1 Factor levels

表 2 正交试验结果及信噪比Tab. 2 Experimental design and results

图 3 三阶振动响应的信噪比Fig. 3 The S/N graph for the third response

图 4 五阶振动响应信噪比Fig. 4 The S/N graph for the fifth response

正交试验的信噪比分析结果如图3 和图4 所示，表示各参数值在各水平位置时对输出结果的影响大小，其中表示输入参数在水平的信噪比均值，如图3中表示厚度在1 水平的信噪比均值。可知，3 个输入参数对连接结构第3 阶频率和第5 阶频率影响均比较大，并且可以初步得到一组较优的组合 A1B3C1。

3.3 优化分析

对正交试验结果进行回归分析，通过回归分析运用最小二乘法可以得到表示固有频率关于厚度、弹性模量和密度的经验公式：

以和的设计值作为优化的目标，通过权重系数法把多目标优化问题转化为单目标优化问题，目标函数如下：

式中： d 为测试的响应值， n为实验的次数。

表 3 优化前后频差对比Tab. 3 The comparison before optimization and after

4 结 语

1）薄层单元的厚度、弹性模量和密度等因素对连接结构的动力学响应影响较大。在弹性模量和密度保持不变时，随着厚度的增加，连接结构的模态频率逐渐减小；在厚度保持不变时，随着弹性模量的增加，连接结构的模态频率逐渐增大；在薄层厚度和弹性模量保持不变时，连接结构的模态频率逐渐减小；

2）通过回归分析，运用最小二乘法建立了连接结构动力学响应关于薄层厚度、弹性模量和密度的经验公式；

3）对薄层参数进行了优化分析，经过优化，建立了精确的连接模型。