悬臂梁结构的应变模态积分法及灵敏度分析.docVIP

悬臂梁结构的应变模态积分法与灵敏度分析 邢建伟， (清华大学航天航空学院，北京100084) 摘 要：关键词：；；； 中图分类号： 文献标识码：A 文章编号： 引 言 1 应变模态积分法 任意截面非均质欧拉梁的无阻尼运动微分方程可写成 (1) 式中为轴向坐标，表示中性面的横向位移，梁单位长度的质量记为，且，其中为梁材料密度、为梁的横截面面积；梁横截面的抗弯刚度记为，且，其中为梁材料的杨氏模量、为梁横截面绕中心轴轴的转动惯量；表示梁单位长度上的横向载荷，梁长。 采用分离变量法令 (2) 分析梁的自由振动，设为自振角频率，引入无量纲参数，并令，则式(1)可写成 (3) 对各向同性材料，仅考虑梁截面轴向的正应力，则应力应变关系为 (4) 其中和分别为轴向处的正应力与正应变，应力不可直接测量，需通过公式(4)由应变换算得到，二者相互对应。 假设截面上某一点距离中性面的距离为，该点处应变与中性面曲率的关系为 (5) 可见截面上任意点的应变均正比于该截面处的中性面曲率，所以可以用欧拉梁的中性面曲率表征应变，即下文中所述的应变与曲率二者等同。 1.1 曲率微分方程 同一个模态频率（为一正整数）下，每一个位移模态振型都会对应一个应变模态振型，这两个模态振型是同一个能量平衡状态的两种不同表现形式。根据欧拉梁理论，梁截面上一点的应变与中性面曲率成正比，那么对应于应变模态，同时也存在一个等效的曲率模态，反映了结构在某一能量平衡状态下的中性面曲率状态。 由于曲率是梁横向位移的二阶导数，则有 (6) 式中和为内层积分变量。 对于悬臂梁结构，边界条件为 (7) 考虑到边界条件(7)，将式(6)代入到方程(3)中得 (8) 由高等微积分知识可知 (9) 则式(8)可化简为 (10) 上式即为以曲率为变量的自由振动微分方程。 将方程(10)等号两端在范围内依次积分两次得 (11) (12) 式中为待定系数，为积分变量，式(12)进一步化简得 (13) 将悬臂梁边界条件(7)代入到上式(11)和(13)得 (14) 求得后代回至方程(13)，最终得到 (15) 至此，对于悬臂梁结构最初的以位移为变量的运动微分方程，通过积分方法转化为以曲率为变量、适用于任意弯曲刚度和质量函数的积分方程(15)。求解该方程可直接获得梁中性面曲率，进而确定结构应变。 1.2 曲率模态方程 为求解积分方程(15)，将曲率展开成幂级数形式 (16) 式中N为一个正整数，令为幂级数系数向量，将上式代入到方程(15)可得 (17) 用乘以上式等号两边，并对等号两边在范围内积分 (18) 以上方程最终可简化为如下形式 (19) 式中分别为曲率刚度和质量矩阵，令，其中，则刚度和质量矩阵的任一元素可表示为 (20) 至此，通过级数展开求解积分方程(15)，将方程化为矩阵形式的广义特征值问题，求解特征值方程(19)得到特征频率和特征系数向量，进而根据式(16)获得曲率模态振型。 2 模态设计 结构的应变模态振型决定了激励载荷作用下的动应变分布，尤其是当激励频率处于模态频率附近的时候。通过修改结构设计参数，以优化应变模态振型，从而避免动应变集中、提高结构寿命。 在初步设计阶段，许多工程结构可以简化成梁模型，例如传动轴、输运管道、高层建筑等。这些结构的质量函数与刚度函数在轴向都是变化的，有限个设计参数决定了结构的刚度与质量函数形式，分析这些模态振型对这些参数的灵敏度是进行模态设计的关键。 2.1 参数灵敏度 由于前面已经将连续梁的曲率模态求解转化为矩阵形式的广义特征值问题，所以就可以采用Nelson法求解特征向量（曲率模态展开系数）参数灵敏度。虽然Nelson在其文献[19]中的分析对象是标准特征值问题（质量矩阵为单位矩阵），但是下面将证明其方法也同样适用于本文情况。 方程(19)中，设为第i个特征值（特征频率），为右特征向量，为左特征向量，满足 (21) (22) 由式(20)可知，刚度矩阵是对称的而质量矩阵不一定对称。将右特征向量对质量矩阵进行归一化 (23) 式中为Kronecker-符号，假设设计参数在处对应N个互异特征值，对方程(21)求参数的导数并左乘，最后可整理为 (24) 其中 (25) 方程(24)是一个非齐次代数方程组，系数矩阵的秩是N-1，不能直