收稿日期 :2009‐07‐13 ;修改稿收到日期 :2009‐11‐12畅
基金项目 :国家自然科学重点基金 (90510018) ;
大连理工大学交叉学科建设专项 (数学 + X)项目
(M XDU T072001)资助项目 畅
作者简介 :胡志强
倡
(1972‐) ,男 ,博士 ,副教授
(E‐mail :huzhq@ dlut .edu .cn) ;
林 皋 (1929‐) ,男 ,教授 ,中国科学院院士 .
第28卷第4期
2011 年 8 月
计 算 力 学 学 报
Chinese Journal of Computational Mechanics
Vol .28 ,No .4
August 2011
文章编号 :1007‐4708(2011)04‐0510‐07
基 于 Hamilton 体 系 的 弹 性 力 学 问 题 的
比 例 边 界 有 限 元 方 法
胡志强
倡
, 林 皋 , 王 毅 , 刘 俊
(大连理工大学 建设工程学部 水利工程学院 ,大连 116024)
摘 要 :比例边界有限元方法是求解偏微分方程的一种半解析半数值解法 。 对于弹性力学问题 ,可采用基于力学
相似性 、基于比例坐标相似变换的加权余量法和虚功原理得到以位移为未知量的系统控制方程 ,属于 Lagrange
体系 。 但在求解时 ,又引入了表面力为未知量 ,控制方程属于 Hamilton 体系 。 因而 ,本文提出在比例边界有限元
离散方法的基础上 ,利用钟万勰教授提出的弹性力学对偶(辛 )体系求解方法 ,通过引入对偶变量 ,直接在 Hamil‐
ton 体系框架内建立控制方程 。 再利用区段混合能和对偶方程得到了有限域 、无限域边界静力刚度所满足的代数
Riccati 方程 ,该方程可采用特征向量展开方法和精细积分方法进行求解 。
关键词 :Hamilton 体系 ;比例边界有限元 ;弹性力学 ;边界刚度矩阵
中图分类号 :O241 文献标志码 :A
1 引 言
钟万勰教授
[1]
根据结构力学和控制理论的模
拟关系 ,在 Hamilton 变分原理基础上 ,发展了求
解弹性力学问题的对偶(辛)求解体系 。 在该体系
框架内 ,也可采用数值方法进行求解 ,如钟万勰
[2]
、
周建方等
[3]
提出的半解析半数值方法 ,选择某一坐
标方向进行解析求解 ,其余方向则采用有限元进行
离散 ,但求解域的形状受到限制 ,大多为条形区域 、
轴对称区域或球对称区域 。
由 John P .Wolf 和 Chongmin Song 两位教
授
[4‐6]
提出的比例边界有限元方法 SBFEM (Scaled
Boundary Finite Element Method) ,仅求解区域的
环向边界需要离散 ,适用于复杂形状的求解域 ,无
需基本解 ,且径向为解析解 ,特别适合于求解含有
无限域
[4]
、裂纹尖端应力奇异
[7]
的问题 。 其控制方
程可由基于力学相似性
[6]
,基于比例边界变换的加
权余量法
[6]
和虚功原理
[8]
等方法得到 ,未知量为位
移 ,属于 Lagrange 体系框架 。 但在求解时 ,引入了
边界表面力作为未知量 ,使得控制方程由以位移
为未知量的二阶线性常微分方程组变为以位移和
表面力为未知量 、系数矩阵是 Hamilton 矩阵的一
阶线性常微分方程组 ,最后在状态空间中进行求
解 。 这一求解方法与对偶体系所采用的方法相同 ,
而且对偶体系在求解含无限域
[1]
和裂纹尖端应力
奇异
[9]
问题也有很大的优势 。 因此基于这些相似
性 ,本文拟采用弹性力学的对偶体系求解框架推导
SBFEM 的控制方程 ,揭示二者之间的内在联系 ,
并为采用 SBFEM 求解偏微分方程提供一种更为
统一的求解框架 ,同时 ,由于利用 SBFEM 可适用
于复杂求解域形状的优点 ,拓宽了对偶体系的半解
析半数值方法的应用范围 。
另外 ,在采用 SBFEM 求解无限域和有限域静
力问题时 ,边界刚度矩阵满足代数 Riccati 方程 ,常
采用特征向量展开的方法来求解
[4]
。 但是 Deeks
et al .
[10]
和 Song
[11]
均指出 ,在二维无限域问题中 ,
如果在某一方向上外荷载不是自平衡的 ,则位移的
解析解中含有对数项 ,从而使得特征向量中出现约
当型 ,导致特征向量不完备 ,无法准确计算静力刚
度阵 。 为了解决这一问题 ,Deeks ,et al .
[10]
提出了
增加 与刚 体平动 位移 相关 联的对 数项 的方 法 ,
Song
[11 ]
提出了矩阵函数和 Schur 分解相结合的方
法 ,这两种方法的目的均是使特征向量完备化 。 而
钟万勰教授提出的精细积分方法
[1]
在求解代数
Riccati 方程时 ,可以避免由于特征向量中出现约
当型所带来的数值求解困难 ,从而提高求解精度和
效率 。