1
“万能”的有限元
在结构分析领域,有限元法是目前应用最为广泛的数值方法。按照未知数的特性划分,有限元法分为位移法、力法和混合法,其中位移法因其高度的灵活性和适应性、强大的模拟能力和求解能力,成为主要的结构分析方法。
桥梁结构分析属于结构分析家族中的一员,当然离不开有限元法(图1)。目前从桥梁的设计计算到施工模拟计算、从理论研究到应用研究,都高度依赖有限元法。可以说,有限元法于桥梁结构分析而言,近乎于“万能”的。有限元法对桥梁工程技术人员如此重要,如同一个忠实的仆人,每天为你服务。但你对她真的很熟悉吗?你知道她也有强烈的个性吗?你如果不尊重她的个性,她就会悄悄地给你设置陷阱,让你吃过期的食品或者给你用假冒伪劣商品。
在对桥梁结构进行有限元建模时,首先要了解桥梁结构分析的理论,了解有限元法原理和特性,只有这样才能建立一个能够反映结构实际特征、符合力学原理的有限元数值模型。忽略力学原理,只把结构计算当作一堆没有物理意义的数字对待,这是目前不少学生容易犯的毛病。本文拟对常规的桥梁结构有限元建模进行讨论,主要讨论其中一些容易被忽略的问题以及建模中的技巧,供大家参考。为了简单,以下均以平面单元和梁单元为例,对于其他单元,其原理是一样的。
(c)
图1 桥梁结构有限元分析(图片a、b引自互联网)
2
问题一:有限元法是一个精确的方法吗?
有些桥梁设计师或年轻桥梁工程师在使用有限元软件的初期,容易有一个潜意识里的错误认识,即认为有限元是一个精确的分析方法。
有限元是一个强有力的、有较高精度的数值方法,但并不是一个精确的方法。简单回顾一下有限元法理论,就能容易地知道这一点。下面以位移法有限元为例来说明。
有限元法是把连续的(具有无穷多个质点的)求解域划分为含有有限个子域(单元)的离散网格形式(见图2)。这些单元之间通过有限数量的网格节点相互连接,用这些离散的节点处的位移表征实际结构中连续的位移场。在这些节点上,位移是连续的,但位移的各阶导数并非都是连续的,取决于节点位移未知数中是否包含这些导数。这是有限元法的第一个近似之处。在单元内部,由于不知道真实的位移场函数,所以有限元法采用近似的数学表达式(如多项式)通过节点位移进行插值来描述,该表达式称为插值函数。这是第二个近似之处。在节点之间的连续边界上,位移及其各阶导数是否连续则取决于单元位移模式和插值函数的形式,一般最多使一阶导数连续。这是第三个近似之处。
图2 有限元网格
在位移法有限元中,采用能量原理建立以节点位移为基本未知数的平衡方程组,所以求解有限元问题归结为求解一个大型线性方程组的问题(即使对于非线性问题,在进行每一迭代步或者增量步求解时,也是一个线性的方程组)。该线性方程组一般为式(1)所示的形式。其中[K]为方程组的系数矩阵,称为总刚度矩阵,{δ}为方程组的基本未知数,即结构的节点位移列向量,{P}为荷载列向量。方程组可能含有成千上万乃至几百万个未知数,那么由于计算机中数字的字长(通俗地说就是有效数字位数)是有限的,因此在求解过程的大量计算中,会由于舍入误差而产生计算误差,这是第四个近似之处。这些只是有限元法本身的近似之处,还不包括力学理论、模型与实际结构差异、计算参数与实际结构差异等引起的误差。