关键词:分片试验,弱形式,网线函数,有限元法
1 引言
目前对各种有限元法分析的方法均是在单元一级上采用变分原理,从而得到单元的应变(或应力)的,由结点位移为参数表达的表达式,再把它们代入最小势能原理得到刚度阵。各种有限元法在得到应变(或应力)的做法上不同,好的有限元法得到的应变表达式已满足了通过分片实验所应满足的条件。
2 分片检验的要求
因有限元法最终列出的是势能的方程,因此分片试验可以看作:在常应变情况下,位移的不协调部分对势能无贡献,在薄板弯曲问题中,可如下表达:
(1) (2) (3)(5)
(6)此即为薄板弯曲问题在单元上的最小余能原理的变分方程。
精化直接刚度法虽然从设位移出发,但又对应变矩阵进行了修正。以下讨论其应变的改进作用。
上式表达了单元的平均应变所应满足的方程。可把上式写成如下矩阵形式:
(15) (16)其单元的平均应变:
(17) (18)使用位移约束方程的方式有两种:第一种是位移的广义参数的个数不增加,改变以往的采用结点参数确定各广义参数的方法,广义协调元和双参数法便是采用这种方法;第二种方法是采用增加位移中的广义参数的做法。此外两种做法也可混合使用。
4.1 广义协调元和双参数法
(19)
广义协调元与双参数法在确定位移广义参数的时候包含上述方程。这两种有限元法得到的位移插值函数在结点处的表达不一定精确,有时会有一个高阶小量的误差。而边界位移条件是直接由结点位移表示的,因此在做分片检验时会有一定的误差,即不很准确地通过分片检验。这一点可由文[8]中的算例看出。
4.2 增加位移中的广义参数
