分片实验与有限元法.docxVIP

分片实验与有限元法 摘要：本文提岀分片试验在有限元法中有着重要的作 用，它是近代有限元发展的一个主要特色。得出分片试验 对位移函数和应变函数的要求，这些要求便是一个好的有 限元法所应保证的；分析了几何方程弱形式与分片试验的 关系，借此分析了杂交元、拟协调元如何满足这些要求， 以及在满足这些要求的同时产生的对其他条件的影响；分 析了精化直接刚度法、广义协调元和双参数法如何保证分 片试验的满足；最后作为位移条件的应用例子，改进了 BCIZ 元。 关键词：分片试验，弱形式，网线函数，有限元法 1引言 连续问题极大地推动了有限元的发展，目前，成熟的 构造单元的方法有传统的位移法有限元［1 ］、应力杂交元 ⑷、杂交混合元［5］、拟协调元⑵⑶、广义协调元［6］、 双参数法［7］、精化直接刚度法［8］等多种。有些方法在数 学上已有证明，但这些方法的更为完善的证明仍是一个课 题，而且其数学证明还很难被研究力学的人们所理解。人 们仍比较普遍以事后的分片试验来验证单元的收敛性。尽 管当前仍有对分片试验的讨论，但以往的大量实践说明： 通过分片试验的单元使用起来是令人放心的。通过分片试 验是绝大多数有限元分析方法的共同点，近期有限元的发 展可以说是以分片试验为一个主要内涵的发展。 众所周知，分片试验是与单元间的位移协调性密切相 关的。人们在进行有限元分析时，不可避免的涉及了单元 间的协调关系，这种协调关系与两个单元有关，文[4] [5] 釆用了单元边界上的公共的位移插值函数，文[9]把这种位 移插值函数成为“网线函数”。正式这种所谓的“网线函数” 的采用，单元间的协调问题可以在单元内独立考虑。目前 成功解决连续问题的有限元法均有意或无意地使用了这种 网线函数。本文通过网线函数给出了分片试验对应变和位 移的要求。 目前对各种有限元法分析的方法均是在单元一级上采 用变分原理，从而得到单元的应变（或应力）的，由结点 位移为参数表达的表达式，再把它们代入最小势能原理得 到刚度阵。各种有限元法在得到应变（或应力）的做法上 不同，好的有限元法得到的应变表达式已满足了通过分片 实验所应满足的条件。 2分片检验的要求 因有限元法最终列出的是势能的方程，因此分片试验 可以看作：在常应变情况下，位移的不协调部分对势能无 贡献，在薄板弯曲问题中，可如下表达： （1） 其中,A：单元域，为位移的不协调部分，有: (2) 为位移，为位移的协调部分。 方程(1 )可以理解为：在常内力情况下，不协调位移对 应变能无贡献。把(2)式代入方程(1) (3) 对(3)式中的项应用格林公式，并应用坐标变换公式： (4) 其中、分别为位移协调部分在单元边界的法向和切向 的导数，即为文中的网线函数，、为单元边界外法线的方向 余弦。对含的项再分步积分得： (r 时)(5) r表示单元的边数，表示结点的位移参数。对(3 )中的 含项也进行分步积分并整理有： (6) 同样，对项再分步积分得： (7) ai、bi、ci为由各边的nx与ny组成的参数，表示位 移函数在结点处的值。 (4)、(5)、(6)、(7 )便是通过分片检验所需满足的方 程。 (4)、(5)是从应变的角度反映了分片试验对单元的要 求，这里称之为应变约束条件；(6)、(7)是从位移的角度 反映了分片试验对单元的要求，这里称之为位移约束条件。 成熟的有限元法都自觉或不自觉地应用了这些条件。 传统的位移法构造的协调元自动满足了上述各式，下 面对其它有限元分析方法进行分类分析。 3使用应变约束的有限元法 方程(4)、(5)是对应变的要求，没有涉及刚体位移， 同时应力和应变之间只有一个线性关系，所以，假设应变 或应力的有限元法都应满足这两个方程。 方程(4)、(5)表达的是应变与位移之间的关系，它们 必然与弹性力学的几何方程： (8) 有着密切的关系。把几何方程()写成弱形式： (9) 、、为权函数，应用两次格林公式变换上述方程： (10) 在上式中，单元边界上的、、分别以它们对应的网线 函数、、代替： (1 1) 如果方程(1 1)中、、是应力的变分，即满足了齐次的平 衡方程： (12) 则方程（12）变为： （1 3） 此即为薄板弯曲问题在单元上的最小余能原理的变分 方程。 方程（11）与（13）便是连续性方程弱形式中的两个典型 形式。在方程（11）与（13）中当、、分别取常数，另两个为零 时，便可得到方程（4）或（5）,即符合分片试验的要求。 拟协调元与杂交混合元便是采用方程（11）对应变或应 力进行离散，而应力杂交元釆用的是（13）式。不同的是应 力杂交元与杂交混合元是由假设应力出发，而拟协调元是 由假设应变入手。而应力与应变之间的关系只是一个线性 变换，如果应力与应变设在同一空间，仅是设应力与设应 变的不同是不会影响最终结果的。