应变推知位移的方法

1、结构变形分析近年来高空长航时(High-Altitude Long-Enduranee)无人预警机越来越受重视，这类飞 机普遍具有大展弦比柔性机翼，机翼同时是一部长基线天线。图1大展弦比柔性机翼大展弦比柔性机翼在飞行载荷作用下，会产生很大的弯曲和扭转变形图2机翼变形示意图结构变形估计为保障机翼的安全性及机翼长基线天线的电性能，要求对机翼的变形进行实时测量，目前主要的实时测量方法有两种：一是模态法；二是 Ko位移理论。模态法在仅考虑Y轴方向弯曲变形时，可以将机翼简化为悬臂梁，如图3所示，长为丨，半径为c的圆截面悬臂梁。图3模态法应变测量点与位移示例由结构动力学可知， 对于分布参数系统，理论上需要

2、采用偏微分方程描述其变形。但实 际上，往往采用有限元方法对其离散化。将结构离散化后，整体结构有限元节点总数为NUE，在外载荷作用P下的变形y可以由下面的运动方程来描述。My cy Ky = P式中，M , C，K，P分别表示结构的质量阵，阻尼阵，刚度阵以及外载荷向量。对于 线性系统，可以由相应振型方程得各阶模态振型。K -co 2m式中，表示结构固有频率；表示结构模态振型向量。通过特征值和特征向量的求解, 可得各阶模态振型如下所示。门22川Nue2川Xue：.：-NUEIII“NUE式中，门i表示第i阶固有频率所对应的模态振型向量；：*表示第i阶模态振型的第j个模 态位移。根据模态叠加原理，结

3、构的变形可以表示成NUEy =、 "qi(4)iV结构位移模态反映的是结构的固有振型，应变是位移的一阶导数。 因此，对应于每一阶位移位移模态，则必有其对应的固有应变分布，这种与位移模态相对应的固有应变分布状态称之为应变模态。= DC：.)式中，表示第i阶应变模态向量；表示第i阶位移模态向量； D表示位移到应变的线 性微分算子。由应变模态的定义可知， 结构在载荷作用下的变形所对应的应变分布也可以由q。各阶应变模态的线性组合来表示，且有相同的线性组合系数， 即有相同的广义模态坐标取前n阶弯曲模态，在x；处的测量应变；j可以表示为n；j 八 J-jqii=1式中，表示第i阶应变模态在Xj；

4、处的分量。m个测量应变的矩阵形式如下12HI甲1人nq1宇2HI中2A nq2I*+-m'2HI-mT n那么91甲2由模态分析得到皿m翅，由测量得到仏讪，在m xn时就可以求出qn辺，如下所T-1 t5叹=（悴4用烏）叫血(8)在求得前n阶弯曲模态所对应的广义模态坐标qn>1，并且已知前n阶弯曲位移模态的条件下，所有参与评价位移估计精度的节点的估计位移就可以表示为?i?2?nJ 爭2III応：.：NIII XN1| qi 一1慣卜炳.q ”©烏（囹异亠）秤/qn(9)简写为N 1 "DST In m mi(10)(13)(13)式中，N表示参与评价位移估计精

5、度的节点总数。Ko位移理论从整个机翼来看，变形量很大，但是逐段分析发现，段内变形量仍是一个小量，并且满足线弹性假设。因此将整个梁沿 X轴方向划分为多个小段， 在每个小段上分别应用经典梁理 论的线弹性小变形假设来推导位移。在测得第j段段首Xj；和段尾X门位置的应变；j和；j 4之后，通过一次插值函数得到第(13)(13)j段内任意位置的应变(12)? X 二；j 一 ；j _ ；j 1« 二 X1 _Xj；X1 乞 X Xj；由欧拉-伯努利梁理论可知，梁表面的弯曲应变 ；X与梁变形后的扰度 yx之间的关系 式为d2y xdX2由此，在已知第j段段首Xj；和段尾Xjl位置的转角 片；和扰

6、度yj;的条件下，通过一次积(13)分得到第j段内任意点的转角 2 X2tan»X嘗dx ta*j、.、x ? X坏一心x +tan日代x$兰x兰x曲(14)通过两次积分得到第j段内任意点的扰度y? xXdxdxtanjdx y；xz x e Xj(15)xjxx x ? X '? Xjtanxdx yj=j.Xj-显然上述积分过程是一般意义上的积分过程，但是必须是在已知 片；和y的前提下，而这在没有完成第j -1段的转角和扰度的计算之前是未知的，即第j段的位移估计必须建立在第j -1段的基础上。另外，由悬臂梁的边界条件可知，在第1段内，转角 -0，扰度y0 ,由此可知要完成整个梁的扰度估计，必须从第1段开始依次完成各个小段内的扰度估计，如此可得第 j段内任意点的转角