弹性力学中三种典型有限体积法方案的精度分析.docxVIP

弹性力学中三种典型有限体积法方案的精度分析 0 有限体积法 有限单元法是一种通用的数值解算方法。根据最小势能原理，有限单元法获得的结果是最好的。然而，有限单元法采用分散法，需要找到与微分方程对应的泛函，并分析数学模型的问题复杂。有限体积法是继有限单元法以来发展起来的新的数值计算方法，在流量计算中取得了巨大的进步。由有限体积法获得的离散方程具有统计学意义清晰、方程形式规范、易于控制性等优点。此外，在有限体积的情况下，有限体积法具有一定的灵活性。国内外科学家提出了不同的曲线曲线法。 本文研究有限体积法在弹性力学中的应用和对比三种有限体积法方案的精度,以二维平面应力问题、三角形网格剖分为例,对弹性力学中有限体积法的三种方案进行理论分析和算例应用为工程计算中合理的选择有限体积法方案提供重要理论依据. 1 角形单元 与有限单元法一样,有限体积法首先需要对整体区域进行网格剖分.但有限体积法还需圈划有限体积,如图1. 图1中实线表示三角形网格剖分,虚线表示有限体积圈划.本文称图1(a)中的四边形P0BM1A、图1(b)中的(35)P0BA等单元为体积单元.图1(a)中, 且Mi(i=1,…,6)为三角形单元重心;图1(b)中, 2 有限体积法的原理 对于弹性力学问题,有限单元法通过虚功原理或最小势能原理来建立刚度方程组;有限体积法则直接由有限体积的力平衡建立线性方程组. 2.1 面力和强化矢量 设在空间内一封闭曲面S,它所包围的区域记作V,对该区域进行力平衡分析,那么面力和体力的矢量和应为0,即 式中,FS为曲面S所受的单位面力矢量;FV为区域V所受的单位体力矢量. 对于二维平面问题,式(1)可简化为 式中,C为区域的边界曲线;S为区域面积;FC和FS分别为曲线C和面积S上所受的单位载荷矢量. 2.2 总体方程的建立 当区域剖分网格后,其边界也相应地被划分为块,本文称每个边界块为边界单元.而相对于边界单元,所有的三角形单元都被称为内部单元. 内部节点是只在内部单元上的节点,边界节点是同时在内部单元和边界单元上的节点.当节点是边界节点时,关联于节点的有限体积受到每个体积单元的应力和体力,还可能受到实际边界上的力.边界节点P0总是关联于两个边界单元,且可能受到集中力,因此根据式(2),力平衡方程为 式中,n为与该节点关联的体积单元的个数;lk为与该节点关联的第k个体积单元的外边界长度;pk为lk上所受总应力密度;lqr第r个边界单元的长度;qr为lqr上所受的载荷密度;FP为节点所受的集中力;Ak为与该节点关联的第k个体积单元的面积;fk为Ak上所受的体力密度. 当式(3)中qr和FP为零,式(3)则变成内部节点的情况,故对内部节点的讨论包含在对边界节点的讨论中.式(3)即 式(4)按x方向、y方向展开,即 式(6)中,σ为正应力,τ为切应力;α为边界外法线方向与x轴的夹角. 令Pk=pklk,Fk=fkAk,Qr=qrlqr,Jk=Pk+Fk,式(4)和式(5)即 式(7)表明以任何一个节点建立的有限体积平衡方程可分解成各个相关内部单元对平衡的贡献(即Jk)、各个相关边界单元对平衡的贡献(即Qr),以及相关点载荷对平衡的贡献(即FP).因此,在建立总体方程组时,可采用单元搜索法,如式(9),分别搜索所有内部单元计算J,搜索所有边界单元计算Q,搜索所有集中力计算FP,将所有贡献总体合成后即可得到要求解的总体方程组. 式中,k为内部单元编号,N为内部单元数;r为边界单元编号,L为边界单元数;s为集中力编号,S为集中力数;i为节点编号,M为节点数;仅当内部单元k关联于节点i时,Jik非零;仅当边界单元r关联于边界节点i且该单元受力时,Qir非零;仅当第s个集中力关联于节点i时,FPis非零. 为便于讨论,本文称图1(a)、(b)所对应的有限体积法方案为FVM1和FVM3.图1(a)的有限体积圈划方法在边界单元上受力的情况,可以使用等效载荷处理,称这种有限体积法方案为FVM2. 3 体积单元的应力 下面以FVM1为例,分析内部单元对其所关联的节点的作用.设四边形i AOB是关联于节点i的有限体积的第k个体积单元,如图2. 首先,考虑应力,因为在(35)OBA中边OA和边OB上所受的力与边AB上所受的力平衡,又根据力的相互作用,所以五边形Ajm BO对体积单元i AOB作用的应力,相当于四边形Ajm B对(35)i AB作用的应力.根据式(5)可得 由文献,式(10)可得 式中, bi=yj-ym,ci=xm-xj(i,j,m);(i,j,m)表示下标轮换,如i→j,j→m,m→i同时更换. 其次,体积力对有限体积的力平衡的贡献为 式中,Aik为关联于节点i的第k个体积单元面积,Ak为内部单元(35)ijm的面积. 最后,可以得到该内部单元对关联于节