有限元收敛性与分片试验条件ξ.pdfVIP

第 14 卷第 3 期 计 算 力 学 学 报 V o l. 14 N o. 3 1997 年 8 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS A ugust 1997 有限元收敛性与分片试验条件 钟万勰 (大连理工大学工程力学研究所, 大连, 116023) 摘　要　本文就平面弹性问题给出了分片试验的条件, 并给出了其明确的物理意义。条件用公式表 达, 并证明了这就是工程式分片试验的充分必要条件; 同时还证明了这些公式也是单元收敛的充分 必要条件, 本文的方法可以推广到有限元分析的其它问题去。 关键词　有限元; 收敛性; 分片试验; 非协调元 分类号　 242 21 O 1　有限元分片试验 分片试验长期以来是有限元分析中行之有效的收敛判据[ 1, 2 ] , 然而其提法则是工程式描述 性的。数学家试图给分片试验以精确的数学定义, 并有鉴于分片试验主要用于非协调元的收敛 性问题, 将其变分试函数违反经典变分法理论的间断取名为变分犯规(V a ria tiona l crim e) , 但 [ 3 ] [ 4 ] 未曾对其间断性质作出规定 , 这是一个很大的缺点。从其数学定义出发, S tum m el 、石钟 慈[ 5, 6 ] 等进行探讨的结论认为, 对于非协调元收敛性, 分片试验既非其必要条件又非其充分条 件。这个结论将分片试验从数学角度上根本否定了。但工程师的看法则截然不同, 专著[2 ] 明 确指出, S tum m el 的质疑应予否定, 并继续认为分片试验是单元收敛的必要条件, 而且只要将 ( ) 力的边界条件也考虑在内 文[2 ] 中的T est C , 则分片试验也是单元收敛的充分条件。两种结 论尖锐对立。鉴于有限元收敛性及分片试验都是重要的概念、理论与方法, 不容回避。因此应 [ 7 ] 从其数学定义寻找原因, 例如其变分犯规所能容许的试函数间断性质等方面 。 分片试验一般是用平面弹性问题来表达的, 本文仍用此问题来表达。首先要明确工程式分 片试验的提法。设有一个多角形单元, 只具有最少的必要弹性支承约束, 而其余点则处理为力 ( ) ( ) 的边界条件 即自由边界条件 , 见图 1; 当这些节点上作用与该单元常应力解 线性位移解 所 相应的边界条件时, 不论在单元内部如何进一步细分单元, 其有限元法的解应仍旧得到该常应 力的解。如果有限元的解得到了偏离于该常应力的解则认为通不过分片试验, 否则就是通过了 分片试验。 以上分片试验的提法并不涉及无限维索波列夫空间等数学提法, 这是其优点; 但仍要采用 细分的有限元数值分析, 故被认为是通俗提法。本文则用此通俗提法, 在有限维空间给出一个 单元通过分片试验的准则, 反过来该准则也可用作分片试验的定义, 容易理解又便于实施, 同 ( ) 时又具有确切性。在本文定义 准则 下, 还可给出分片试验为有限元收敛性的充要条件的证 明。当然也要证明本文准则与工程式分片试验提法相等价, 即本文的分片试验准则为工程分片 试验的必要与充分条件。这是等价性的充要条件证明; 在此之后再证明收敛性的充要条件。这 收稿日期: 　钟万勰: 男, 1934 年生, 教授, 中科院院士 © 1995-2004 Tsinghua Tongfang Opt