弹性梁动力响应分析的一种辛算法.pdfVIP

218 黄伟江等：弹性梁动力响应分析的一种辛算法 * 弹性梁动力响应分析的一种辛算法 黄伟江 罗 恩 佘 慧 中山大学应用力学与工程系, 广州 510275 摘 要 本文根据古典阴阳互补和现代对偶互补的基本思想，首次建立了线性阻尼情形下弹性 ( ) Hamilton 梁动力学的相空间 挠度、动量 非传统 型变分原理。这种变分原理不仅能反映这种动 , 力学初值-边值问题的全部特征，而且它的欧拉方程具有辛结构的特征。基于该变分原理 提出一 - 种称之为辛空间有限元 时间子域法的辛算法。这种新方法是由空间域采用有限元法与时间子域 采用 Lagrange 插值多项式插值的时间子域法相结合而成。文中用这种辛算法分析了四种支承条 件下弹性梁的动力响应问题。算例的计算结果表明，这种新方法的稳定性、收敛性、计算精度 和效率都明显高于国际上常用的 Wilson-θ 法和 Newmark-β 法。 关键词 相空间；非传统 Hamilton 型变分原理；初值-边值问题；辛算法；动力响应 中图分类号 O322 1 弹性梁动力学相空间非传统 Hamilton 型变分原理 按文[1，2] 的思路，可以建立有阻尼弹性梁动力学的相空间 （挠度、动量）非传统 Hamilton 型变分原理。对于边界条件 x = 0 : w = w , θ = θ ; x = l : M = M , Q = Q 以及初始条件 ~ ~ w(x ,0) = w0 , p (x ,0) = p 0 的情形，其泛函表达式为 2 t1 l  1 2 EI ∂ 2 w  。  Π 2 (w, p) = ∫∫ p w − p −  2  − c w w + f wdxdt 0 0  2m 2 ∂ x   2 2 t1  ∂  ∂ w  ∂w ∂ w  t1  ∂w  −∫0 (ω − ω ) −EI 2 − EI  − θ  2  dt + ∫0 Qw − M  dt  ∂x  ∂x  ∂x ∂x x =0  ∂x x =l l  。   + p (x)w(x, 0) − p (x ,t )w(x ,t ) + w(x, 0) − w (x) p (x , 0)dx (1