杆系结构自由振动的精确有限元法及动力刚度法的等价性.pdf

杆系结构自由振动的精确有限元法与 动力刚度法的等价性1 袁驷，叶康生 清华大学土木工程系（100084 ） Email: yuans@, yeks@ 摘 要：本文论述了动力刚度法和精确有限元法的等价性，证明了等价性定理、梁氏定理及 本文导出的微分等式定理，并为新近提出的导护型 Newton 法求解振型和频率提出了一个 改进的外推公式。文中给出了算例用以表明改进公式的良好效果。 关键词：自由振动，精确有限元法，动力刚度法，杆系结构 1．引言 对于杆系结构的静力分析计算，传统的刚度法 （矩阵位移法）和现代的有限元法在某 些情况下 （如等截面杆件）殊路同归，人们常将二者混为一谈。其实，两种方法的做法是 有区别的： 刚度法：求得满足控制微分方程的解，由其解通过取导数直接计算杆端刚度系数，汇 集后得到单元刚度矩阵。 有限元法：构造单元（杆件）形函数，用能量法（虚功原理、最小势能原理）通过积 分生成单元刚度矩阵。 这里，之所以可以得到相同的单元刚度矩阵，其要点是有限元中采用了精确的形函数，即 满足控制微分方程的形函数（如等截面梁单元采用三次多项式为形函数）。换言之，两种 方法的共同点是都采用了满足控制微分方程的解，刚度法利用其解直接计算刚度系数，而 有限元法则利用其解构成形函数并用能量法计算刚度系数。 对于杆系结构的动力分析计算（自由振动问题），也相应地存在两种方法：动力刚度 法（Dynamic Stiffness Method, 简称 DSM ）和精确有限元法（FEM ）。但是，多年来，基 本上只有动力刚度法得到发展，而精确有限元法则很少为人提及。本文首先介绍动力刚度 法，然后讨论精确有限元法的列式，在此基础上进一步讨论两种方法的等价性，由此等价 性可导出梁氏定理及另一个等式定理，并为新近提出导护型 Newton 法求解振型和频率提 出了一个改进的外推公式。 1 本课题由高等学校博士学科点专项科研基金（20020003045 ）和教育部科技创新工程重大项目培育资金项 目（704003 ）联合资助。 - 1 - 2 ．模型问题 为简明起见，以杆件轴向振动作为模型问题来举例说明，其控制微分方程可表示为 [1][8] d du x ⎛ ( ) ⎞ 2 − EA m u x , ⎜ ⎟−ω ( ) 0 0 x l (1a) dx ⎝ dx ⎠ EA ω 其中， 为轴向振幅函数， 为分布质量， 为截面抗压刚度， 为自振频率， 为杆 u(x) m