模态子空间中非线经典阻尼结构动力响应的求解.pdfVIP

·834· ‘工程力学》增刊 1999年 模态子空间中非经典阻尼结构 动力响应的求解 刘保东朱唏 (北疗变通大学上木建筑工程学院．北京1000¨) 摘要通过对非经典阻尼结构动力响应算法的分析． 提出在模态于空间中用模态叠加法求解非经典 阻尼结构动力响应时，应该首先确定非经典阻尼的程度， 根据不同的非经典阻尼程度应用j；W的求解方 法。最后通过算例说明本文所提算法的精度。 关键词非经典阻尼，模态，动力学 1引言 用有限元法求解n个自由度结构体系动力响应时，动力学方程可以写为： 阻，料+【c，酝)+k。M=识) (1) 式中阻，】、虹】、瞳。】分别为nxn阶的质量矩阵、阻尼矩阵和剐度矩阵，扛)为n×I 阶的结构位移向量，{只}为n×1阶的节点荷载向量。 用模态叠加法求解方程(1)时，可作如F的坐标变换…【j!： 扛j：眵怕} f2) 式中眵J为n×n阶的模态振型矩阵，{g}为广义坐标向量。 则方程(1)可化为： 阻，D怡)+【C，D№}+k，髟k}-把) (3) I}j眵】』左乘方程(3)可得到： ，、，吐’蜂}+1[c’k)+k’k)={F‘} ¨： 式中【M‘J、p‘J、k’J分别为n×n阶的广义质嚣矩阵、广义阻尼矩阵和广义刚度矩阵， 扩’}为n×l阶的广义节点荷载向量。 根据振型的正交性，如果振型采用规格化振型，则lM+l为对角元素为l的对角阵， 瞳’J为对角元素为n彳(f=1……刀)的对角阵，其中q为第f振型的圆频率。 对丁经典阻尼结构，方程(4)可解耦为： 可i+2善，∞，口，+山?g，=，(f=1··-···H) (5) 刘保东，男．1967 12出生．工学硕士．讲师 {工程力学》增刊 1999年 ‘B35’ 式中、，：=协}’{C}为第f振型的广义力，岳为第f振型的阻尼比。求解出各个广义坐标吼 后，代同方程(2)即可求出结构的位移向量{x}。由于在很多I：程实际问题中．结构的响 应往往只考患较低的前儿阶模态就完全能满足1．科精度要求，所以复杂结构的动力响应可 以通过随瑚模态叠加法得到简化。 但是对j二1F经典阻尼结构，Ic’I不能化为对角矩阵，m0方程(4)不能解耦，就不能 直接应H1模态叠mn法，F面讨论这种结构的动力响应求解方法。 2非经典阻尼结构的动力响应 具有两种或两种以上特性不同材料的结构一般为1|经典阻尼结构，很多1．群实际问题 如1．与结构相互作Ⅲ问题、海洋平台结构考虑动水阻尼问题、Hj阻尼器减振问题等都属丁 1r经典阻尼结构问题。研究非经典阻尼结构的动力响廊问题，寻找一种既满足l：程精度要 求义简单易行的方法是十分必要的。 求解非经典阻尼结构动力响应的方法人致有以F儿种“】：第一．直接积分法．此方法 适_【}；jr求解荐种情况的非经典阻尼结构，结果也相当精确，但对丁自由度比较多的复杂结 构，此种方法计算J：作鲑较人：第二，状态空间中的模态蘑加法，用此种方法求解时矩阵 的维数增加一倍，而且要进行人鼙的复数运算：第三，先按无阻尼情况计算结构的若干低 阶模态的振州和频率，然厉按经典阻尼模态叠加的方法进行模态分解，此时的r_义阻尼矩 阵不能化为对角阵，方程仍不能解耦，但如果忽略广义阻尼矩阵中的非对角线元素．则可 以按经典阻尼的方法求解结构动力响应，此种方法当IF经典阻尼氍度较小时，精度可满足 要求．但当。1F经典阻尼稗度较人时．误莘不可忽略。 如震把对结构响应贡献较人的棋态称为主要楼态，在求解1E经舆阻尼结构动力响膨 时，考虑与：土=要模态耦合程度较人的模态的影响，而忽略其它模态的影响，