計算实习报告.docxVIP

二、解决方案基于Galerkin加权余量法的有限元方法求解，同时借助AutoCAD编制网格和Matlab编程计算。三、有限元分析可知此算例所示流动满足由流函数表示的Laplace方程：式中：边界上为强制（本质）边界条件，函数值为给定值；边界上为自然边界条件，变量的法向导数值由确定。近似解：，将近似解代入上式，得到余量由迦辽金加权余量法：取权函数为插值函数，则有：其弱形式：式中：为边界，为在边界上的边界条件。 故可写成：式中：式中：为；k为单元特征矩阵或刚度矩阵，f为单元荷载列阵。三角形单元的局部插值函数为：式中：，|D|等于三角形单元的面积A的两倍。所以可得四、边界条件处理由于本算例为对称结构，故可取1/4结构分析。其边界条件表示如下：图二：算例边界条件（1）第一类边界条件（强制边界条件）由图可知，用第一类边界条件表示的边界有bc边：。ae、ed边：。对于第一类边界条件，通过改写总体系数矩阵和总体荷载矩阵得到满足，稍后会在程序中说明。（2）第二类边界条件（自然边界条件）由图可知，ab、cd边用第二类边界条件表示：。对于边界上，是给定的，取沿边界的坐标系s-n,它与x-y坐标系的关系如下，其中l为边界2－3的长度，将上式代入得到：在边界上的法向坐标为n＝0，使用Lagrange插值多项式条件，得到边界上的插值函数为：，，最终的到结果为：代入此算例的第二类边界条件，可知对应的所需改写的