有限元中的半解析法.doc

有限元中的半解析法 在现实生活中，由于求解的问题复杂、规模较大，常规的有限单元法的费用较高，已经不适用。因此，我们希望找到其他的方法以减少计算工作量，降低费用。这时，半解析法具有其优势，它是一种离散与解析相结合的方法。目前，常用的半解析法有三种：组合条---元法，有限元线法，有限条法； 1. 组合条---元法（Combinatory Strip-Element Method） 为了克服有限条半解析法所存在的局限性，又能存有有限条的一些优点，我们又提出了所谓的组合条---元法。这是一种将有限条和有限元的特点组合起来的方法。下面我来介绍组合条---元法的构造思路。 与条元法一样，以窄长条带为单元，不同的是，节线的两端设置有结点。由于任一函数均可由完备函数集中的基函数来表达，可采用如下两步法构造单元位移场： （1）由结点位移参数，采用形函数插值构造条带单元的节线位移，这一步与有限元一样。 （2）以上述节线位移作参数，沿条带短边方向进行多项式插值，从而构造条带的位移场。 经过以上两步，即可得到 （7） 式中 -----由结点位移参数构造的位移部分； -------沿长边方向由级数构造的位移部分。 然后就可以按普通有限元进行分析。 这种方法克服了有限条法的缺陷，币有限元减少了很多未知量。使用这种方法我们解决了平面问题、薄板问题、折板与平面壳体等的线性与非线性、静力与动力分析；并联合应用了有限元、组合条元与映射无限元求解过路面力学问题。是一种可行的方法。 2有限元线法（Finite Element Method of Lines） 有限元线法（FEMOL）是袁駟提出的一种以常微分方程求解器为支撑软件的新型半解析法。有限元线法的构造思路有以下几步： （1）建立参数FEMOL的单元映射。为适应复杂形体问题的计算，可建立母单元与子单元的映射关系。 （2）构造参数FEMOL的变量场：单元上的变量场可由节线未知函数ui（η）通过ξ方向的形函数Ni(ξ)插值得到。 （8） 式中 （3）参数FEMOL的能量泛函的确定：结构中每个单元的能量为∏e，它是η的函数。则整个求解域的能量为： （9） （4）建立常微分方程体系：常微分方程建立后，经过一系列的处理后即可用求解器（Solver）来求未知节线位移函数。 有限线元法中，由于引入参数单元，是可用于不规则区域的求解；由于未知节线位移是通过解常微分方程组得到的，其自然精度要比其他方法高。也是一种很有效的半解析方法。 小结：半解析法的具体方法有多种，这里只介绍了三种方法。并对有限条法作了详细的介绍。在实际中每种方法都有其优势，也有其不足。我们应根据具体的情况和要求，采用某种合适的方法，或者联合使用多种方法进行具体分析，已达到要求的目的。 3有限条法（Finite Strip Method） 有限条法是由张佑启先生提出的一种方法，用以解决规则形体问题。本方法具有工作量小、精度高的优点。下面将以薄板为例，介绍位移场的构造方法。 如图1（a）所示，有一矩形薄板，设每条边界的支承条件相同，图中表示了三种支承情况，图1（b）用一些与边界线平行的直线将板分割成若干窄长的条带以此组成有限元分析中的单元。下面介绍这种条带单元位移场的建立思路。 图1 矩形薄板与有限条离散示意图 确定位移模式 对于薄板来说，挠度ω可用分离变量形式表示 (1) 边界条件的确定 本例中也可取Xm(x)满足条带两端的边界条件的梁振形函数，它是如下微分方程的解： （2） 式中 ---是振型参数，由边界条件确定 图1(a)所示的是一端固定一端简支的情况，则有： （3.1） 式中，振型参数由tg=tanh确定，具体取值见表1。 表1 的取值 m： 1 2 3 ≥4 3.9266 7.0685 10.2102 ≈ 其他对边约束条件情况振型函数Xm(x)为： (1) 两端简支 =mπ （3.2） (2) 两端固定 （3.3） 由cosch=1确定，取值如下表2。 表2 的取值 m： 1 2 3 ≥4 4.73004 7.8532 10.955608 ≈ （3）两端自由 （3.4） 在m≥3时，由表2中取得。 （4）一端固定一端自由 （3.5） 由cosch= -1确定，取值如下表3。 表3 的取值 m： 1 2 3 ： 1.875 4.694 ≈ （5）一端