基于绝对节点坐标法的柔性多体系统灵敏度分析

基于绝对节点坐标法的柔性多体系统灵敏度分析
王铁成;陈国平;孙东阳
【摘要】For multiple design variables,an adjoint variable method was applied in sensitivity analysis of flexible multibody systems based on the absolute node coordinate formulation.In order to verify the computational efficiency of the method,a flexible pendulum under gravity was analysed by using the direct differentiation method and the adjoint variable method respectively.The results show that the errors of the two methods are both small,and the adjoint variable method has higher computational efficiency with the increase of design variables.%针对多个设计变量情况下的柔性多体系统灵敏度分析，采用伴随变量法对基于绝对节点坐标法建立的柔性多体系统进行了研究。

为了验证的该方法的计算效率，分别采用直接微分法和伴随变量法对受重力作用的柔性单摆进行了研究，结果表明：这两种方法计算结果的误差很小，随着设计变量数量的增加，伴随变量法有更高的计算效率。

【期刊名称】《振动与冲击》
【年(卷),期】2015(000)024
【总页数】5页(P89-92,103)
【关键词】多体系统;绝对节点坐标法;灵敏度;伴随变量法
【作者】王铁成;陈国平;孙东阳
【作者单位】南京航空航天大学机械结构力学及控制国家重点实验室，南京210016;南京航空航天大学机械结构力学及控制国家重点实验室，南京 210016;重庆大学航空航天学院，重庆 400044
【正文语种】中文
【中图分类】O313.7
机械系统优化设计可以提高机械系统的精度和稳定性。

为了避免优化参数选择的盲目性，提高设计效率，进行近似分析和预测设计趋势，首先就需要对系统参数进行灵敏度分析。

因此，灵敏度分析已经成为机械系统动力学分析与优化设计的桥梁，并成为机械系统优化设计的核心问题。

灵敏度计算方法主要有直接微分法和伴随变量法。

直接微分法是通过计算系统的运动学或动力学方程对设计变量的偏导数，得到目标函数对设计变量的灵敏度。

潘振宽等[1]采用直接微分法对多刚体系统进行了灵敏度分析。

针对递推形式建立的链状多体系统，Haug等[2]采用直接微分法对其进行了灵敏度分析。

直接微分法的优点是直观，易于理解，但随着系统设计变量的增多，求解规模成倍增长，导致灵敏度求解的效率降低。

伴随变量法是通过引入伴随变量，使多设计变量系统灵敏度分析计算效率得到提高，该方法已经在多体系统的灵敏度分析中得到了重视和应用[3] 。

Etman等[4]采用伴随变量法对多体系统动态特性进行了灵敏度分析。

Li等[5]和Maly等[6] 基于隐式微分代数方程，推导了多体系统灵敏度分析的伴随变量法。

目前，多体系统灵敏度分析大多是针对刚体系统，柔性多体系统灵敏度分析的研究还比较少。

主要原因是柔性多体系统具有高维数和强非线性特性，导致分析相对困难。

虽然浮动坐标法是最常用的一种柔性多体系统建模方法，但是该方法基于小变形假设，用于分析存在大变形的柔性多体系统会带来较大误差。

Shabana等[7]提出了多柔性体绝对节点坐标建模方法，该方法的理论基础是有限元和连续介
质力学[8-9]，由该方法推导的多体系统方程的质量阵为常数矩阵，且不存在科氏力和离心力项。

为了对存在大变形的柔性多体系统进行灵敏度分析，同时考虑到多变量情况下伴随变量法有更高的计算效率，本文采用伴随变量法对基于绝对节点坐标法建立的柔性多体系统进行了灵敏度分析，并将其计算结果与直接法的计算结果进行了比较。

基于绝对节点坐标法的一维两节点梁单元如图1所示，单元上任意点的位置矢量可表示为：
根据式(1)，单元的动能可表示为
根据连续介质力学理论，单元的总应变能包含弯曲应变能Uel和轴向拉伸应变能Uet，表示为
基于虚功原理建立单元的动力学方程表示为
单元广义弹性力由单元应变能对单元坐标求偏导获得，即
则含约束的柔性体k的动力学方程为
伴随变量法计算灵敏度，需要先求解动力学方程，然后根据动力学计算结果，逆向求解伴随变量为未知量的微分方程组，其中，动力学计算结果需要拟合函数，可使在伴随变量微分方程的数值计算时，能取得在所计算时间内的任意时刻的函数值。

本文选用积分型目标函数对绝对节点坐标法建模的柔性多体系统进行灵敏度分析。

目标函数设为
状态变量和设计变量分别为
多体系统动力学方程和约束方程分别为
为了计算灵敏度，先用莱布尼兹法对目标函数求导，得到
伴随变量法的基本思想是通过选取伴随变量消除式(14)中b和Fλλb，即消去目标函数中含有对状态变量以及拉格朗日乘子的导数，而保留目标函数对设计变量的导数。

为了达到这个目的，首先对式(14)中的第三项应用分部积分法得
假设初始条件不依赖设计变量，即qb在初始时刻为零，也就是qb在t=t0时的
计算值为零，这样式(15)可以表达为
引入伴随变量μ和ν，分别乘以式(12)和式(13)并积分得
式(17)和式(18)对设计变量求导得
利用分部积分法积分对(19)式中含b和b的项进行积分得
将式(15)，式(20)和式(21)相加得
选取伴随变量使得积分式中qb和λb前的系数分别等于零，得如下关于伴随变量
的微分代数方程
选取伴随变量满足在t=t1时，式(22)中的和q前的系数为零，即下式成立
再由式(25)推出下面两式
采用绝对节点坐标法建立的柔性多体系统动力学方程，质量矩阵M为满秩常数阵，进一步可得到
根据初始条件式(28)、(29)，逆向求解关于伴随变量的微分代数方程组(23)、(24)，将计算结果代入式(16)，得目标函数灵敏度
本文以受重力作用的矩形截面柔性单摆系统为研究对象，如图2所示。

设计变量为柔性单摆长度, 密度，弹性模量，单摆横截面的长度和宽度。

在分析过程中，将柔性梁等分为4个单元，柔性单摆密度、弹性模量、长度、截面的长度
和宽度分别取为ρ=2 770 kg/m3，E=1 GN/m2，l=1 m，a=31 mm，
b=7.9 mm。

柔性单摆的截面对中性轴惯性矩I=ab3/12。

目标函数为(x-1)2dt，
式中x为单摆末端的横向坐标。

表1给出了分别采用伴随变量法和直接微分法得到的目标函数对单摆的密度、弹
性模量、长度、截面长度和宽度的灵敏度。

同时对比了两种方法计算结果的相对误差，可以发现，其相对误差都在4%以内。

而且可以发现，目标函数对单摆长度、
截面长度和宽度的变化比较敏感，而密度和弹性模量对目标函数的影响很小。

采用直接微分法计算的单摆末端横向位移对五个设计变量的灵敏度如图3~图7所示，由图可知五个设计变量中长度的变化对末端横向位移的影响最明显，这与基于目标函数的灵敏度分析结论一致。

为了分析伴随变量法和直接微分法的计算效率，对比了五个设计变量情况下两种方法的计算时间，如表2所示。

可以看出，考虑一个设计变量时，伴随变量法所用计算时间多于直接微分法，约为直接微分法计算时间的2倍；考虑两个设计变量时，两种方法的计算时间基本相当；三个设计变量时，伴随变量法所用时间约为直接微分法的75%，当五个设计变量的情况时，伴随变量法的计算时间约为直接法计算时间一半。

由此可见，随着设计变量的增加，伴随法的计算效率将高于直接微分法。

这主要是因为随着设计变量的增多，伴随变量法随之增加的计算时间主要就在灵敏度积分函数上，在此算例中大概耗时25 s左右，而直接微分法随着设计变量增加，需要再计算一组微分方程，势必增加求解时间。

本文对基于绝对节点坐标法建立的柔性多体系统进行了灵敏度分析。

以矩形截面的柔性单摆系统为例，以工程上常用的积分型函数为目标函数，使用伴随变量法和直接微分法分别计算了目标函数对柔性体长度、弹性模量、密度、截面的长度和宽度的灵敏度。

通过分析发现，两种方法计算结果的相对误差均在4%以内，因此，计算结果是可靠的。

同时可以发现随着设计变量的增加，伴随变量法比直接法有更高的计算效率。

为工程上，采用伴随变量法对存在大变形的柔性多体系统的进行灵敏度分析提供了理论依据。