Trefftz法中的自锁现象及变量减缩法①.PDF

固 体 力 学 学 报 Vol. 21 　No. 4 　　 　　　　　　 　 　　　　　　　　　 　 第21卷 第4期 2000年 12 月 ACTA MECHANICA SOLIDA SINICA December 2000 Trefftz 法中的自锁现象及变量减缩法① 金吾根 张佑启 ( 复旦大学应用力学系 ,上海 ,200433) (香港大学 ,香港) 摘 　要 　用 Trefftz 型边界解法分析了中厚板弯曲问题 ,发现了一类新的、因 Trefftz 函数溢机 引起的“自锁现象”,并提出了一种消除这类自锁问题的“变量减缩法”. 关键词 　Trefftz 法 , 自锁现象,变量减缩法 1 　引言 各类边界解法已广泛用于解中厚板问题. 考虑板内横向剪切变形的中厚板问题 ,其控制 方程为六阶的偏微分方程 ,边界条件有三个. 尽管用传统边界元法解中厚板问题已有多 年[1～4 ] ,但是以完备解为基础的 Trefftz 法求中厚板问题还是近几年的事[5 ,6 ] . Trefftz 法采用 ( ) 光滑的、非奇异的 Trefftz 完备解为权函数和 或 试函数 ,其边界积分方程是非奇异的 ,从而 避免了奇异积分. ( ) ( ) 本文以中厚板为例 ,首先介绍了 Trefftz 直接法 TDM 和间接法 TIM . 在改进的中厚板 理论基础上[7 ] ,引入胡海昌教授提出的两个广义位移函数 F 和f [ 8] , 将六阶控制方程分离为 二个低阶偏微分方程 :一个是以 F 表示的板弯曲方程, 另一个是以 f 表示的弹性地基上的 ( ) 薄膜方程. 以上述方程的完备解为权函数和 或 试函数构造相应的边界积分方程. 其次介绍 了 Trefftz 法中发现的一类新的“自锁现象”. 它与有限元曾出现过的那类“自锁现象”, 从产生 的原因, 到克服的办法都不一样. 文中提出了一种所谓的“变量减缩法”, 有效地解决了这类 新的“自锁现象”, 数值结果表明方法是有效的. 2 　基本公式 改进的中厚板理论指出, 极坐标下的平衡方程为 Qr 1 Qθ 1 r + r θ + r Qr + q = 0 M r 1 M rθ M r - Mθ ( ) r + r θ + r - Qr = 0 1 M rθ 1 Mθ 2 r + r θ + r M rθ- Qθ = 0 式中 M , Q 和 q 分别为板的弯矩、扭矩、横向剪力和横向荷载, 且忽略体积力. 根据改进的 αβ α ① ( ) 国家自然科学基金 批准号 资助项目; 同济大学固体力学国家教委开放研究实验室资助. 收到第 1 稿 收到修改稿. © 1995-2004 Tsinghua Tongfang Optical Disc Co