控制体有限元方法中质量加权辽风格式的迁延修正

控制体有限元方法对于处理流动问题提供了一种较好的数值模拟方法，但对于大Ｒｅｙｎｏｌｄｓ数的回流流动问题，该方法在采用各种迎风格式时仍存在着严重的假扩散，为此，对二维质量加权迎风格式，给出了一种由插值单元及其一个与插值点最靠近的相邻单元所确定的多项式插值格式迁延修正的形式，通过对Ｒｅ＝４００、１０００、５０００三种情形的二维空穴流问题的计算，表明该格式可有效地减少控制体有限元方法地回流问题数值模拟的     

对流扩散方程的高精度差分有限元破开算子法

对于非定常二维对流扩散方程，将二次迎风插值方法（ＱｕａｄｒａｔｉｃＵｐｏｔｒｅａｍＩｎｔｅｒｐｏｌｓｔｉｏｎｆｆｏｒＣｏｎｕａｃｔｉｖｅＫｉｎｅｍａｔｉｃｓ，简称Ｑｕｉｃｋ格式）与集中质量九节点等参元的有限单元法结合起来，发展了一种有效的高精度差分有限元破开算子法．计算结果表明，这是一种高精度的数值模式．     

控制体有限元方法与流向型插值格式

给出在同位网格上对二维定常不可压层流数值模拟的控制体有限元方法及其流向型插值格式的计算过程,并对几个层流流动问题进行数值计算。     

一种同位网格上的迎风Galerkin 有限元方法

通过 Green定理将对流项变量从微分算子中分离出来 ,从插值函数入手引入迎风格式 ,是对强对流问题 Galerkin有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速—压同位网格公式 ,构成了一种对高 Reynolds数流体流动数值模拟比较有效的有限元方法。数值试验表明 ,采用该方法能较好地提高计算精度。     

基于间断有限元求解浅水方程

基于间断迦辽金有限元方法建立了求解二维、深度平均浅水方程的模型，该模型由浅水方程这一系列双曲方程组按照守恒律推导而得．通过在某一单元积分方程组并乘以基函数可以得到一个弱解，并假设未知量为间断分布的多项式逐个单元求解．由于间断有限元的本质具有局部特征，因此很容易通过提高插值多项式的次数来提高模拟精度．同时，间断有限元又具有局部守恒特征，可以结合数值通量来模拟高流动梯度问题．最后，给出了间断有限元的数值模拟结果。     

基于控制体的分步有限元算法及其在非正交网格上的应用

利用Newton-Raphson算法的局部收敛特性,建立了一种基于分步思想的有限元算法。该算法首先利用基于流动条件构造差值函数的迎风有限元方法获得流动问题的近似收敛解;然后,在近似解的基础上,引入混合插值格式,进一步计算降低误差。混合插值函数的形式由数值实验获得,数值结果与理论分析相一致。对方腔顶盖驱动流和倾斜空腔驱动流进行了数值模拟,数值解与基准解吻合很好。与传统的基于流动条件构造差值函数的迎风有限元方法相比较,该文算法能在较稀疏的网格条件下获得比较精确的计算结果,尤其针对存在漩涡的流动问题,能大大提高数值解的精度。     

污染物大气扩散数值模拟中伪扩散的控制

用数值模拟的方法预测污染物在大气中的扩散已经被广泛采用,但数值模拟产生的伪扩散一直是影响模拟结果精度的重要因素.为此分析了数值模拟中伪扩散产生的原因,用一阶迎风差分和对流项二次迎风插值格式对一稳态污染物扩散模型进行了计算,并比较了计算结果,得出了有实际意义的结论.     

澳门周围水域潮汐流动的数值模拟

采用有限体积法求解二维浅水流动方程,建立了一种在非规则结构化网格上数值求解平面二维浅水流动的数值方法,其中应用Roe格式计算数值对流通量,并通过MUSCL途径插值界面状态,以对水位插值代替水深插值来消除地形变化对格式的影响;对源项采用分步法及半隐式求解,防止了源项引起的计算不稳定;处理动边界采用人工小水深的方法．应用于澳门水域潮流的模拟,计算结果与实测吻合较好．     

一维和二维溃坝波的混合有限分析解

针对水深平均二维浅水方程高分辨率模型，采用非交错网格下的二维对流扩散方程的混合有限分析法，运用SIMPLE算法结合Rhie and Chow的动量插值技术计算二维溃坝水流．一维理论解和其他格式的二维数值结果符合良好，表明非交错网格下的混合有限分析法能正确有效地求解溃坝问题．进而对部分溃决情况下的溃坝波传播和绕流问题进行了数值模拟，给出了精细的数值结果，揭示了复杂的运动特性，进一步表明混合有限分析法具有较强的溃坝波的捕捉能力，是求解溃坝流动的有效方法之一．     

二维线性对流扩散方程一种新的特征差分算法及收敛分析

文章讨论了二维线性对流扩散方程,将特征线法和有限差分法相结合,借助于双线性插值,给出了求解二维线性对流扩散方程数值解的一种新的特征差分格式,并分析了该算法的收敛性。此算法表明对于一类对流扩散方程,应用此差分格式,能更有效地消除数值振荡现象,从而极大地提高数值逼近度。     

应用控制容积有限元法分析注射充模过程中塑料熔体的流动与传热

将Ｂａｌｉｇａ和Ｐａｔａｎｋａｒ提出的控制容积有限元法推广用于分析注射充模过程中熔体的流动与传热，并叙述了其特点．应用此方法计算分析了包括合不同几何形状嵌件的模腔内充模过程的流动与传热问题．得到了一些有益的结论．这将有利于指导在复杂型腔充模过程流动分析时的简化处理．     

裂隙岩体的弹塑性本构模型及有限元分析

首先用弹塑性理论建立起裂隙的非线性应力应变关系,然后引入表征单元体的概念,定义了平均应力及平均应变。在此基础上建立了裂隙岩体的本构模型。该模型克服了应变叠加原理的不足,同时有效地模拟了裂隙的强烈非线性变形及剪帐等特性,编制了相应的三维有限元分析程序,并给出了算例分析,计算表明,裂隙的存在对整个岩体的变形及稳定性起控制作用。模型可直接应用于裂隙岩体的渗流应力耦合分析研究。     

轴向变截面环空幂律流体内边界力的研究

在石油钻采工程中,井筒内流体多为非牛顿流体,在旋转杆管柱作用下做偏心环空螺旋流,而且偏心度沿杆管柱轴线呈任意变化.针对这一流体特性,建立了轴向变截面偏心环空幂律流体三维有限元模型,采用基于有限元的有限体积法,研究了流场分布规律及偏心度和转速对界面力的影响.计算结果表明:宽间隙侧最大轴向速度不变,旋转速度降低,窄间隙侧最大轴向速度降低,旋转速度增强.随着转速或偏心度的增加,流体对杆管柱作用的横向力和弯矩均增大,为进一步揭示和防治杆管偏磨提供了理论依据.     

油藏热流固耦合渗流问题的有限元方法

提出了一种有效实用的求解油藏热流固耦合渗流问题的数值计算方法。该方法以有限元法为主,结合有限差分法,用有限元法求解耦合温度场方程和岩石耦合变形方程,用有限差分法求解流体耦合渗流方程,发挥有限元法网格技术和单元划分灵活的特点及处理复杂的油藏边界优势,兼顾了有限差分法在流场分析方面的成熟应用,使复杂的热流固耦合数学模型得以完整求解,取得了单由有限元法或有限差分法难以取得的效果,是一种新型的油藏数值模拟方法。     

流体瞬变有限元分析

本文提出一种二维非定常微可压缩粘性液体流动计算的有限元瞬变分析方法及其计算模型,据此建立了采用全主元策略的波阵法分析软件,并以快速锻造液压机管道的流体瞬变计算实例说明其应用.     

有限元分析法确定Johson Cook本构方程材料参数

采用反演法确定了马氏体不锈钢410的Johnson-Cook本构方程材料参数.剪切实验采用的剪切速率为1~5 000 s-1,温度为20~1 000 °C.有限元模拟在Abaqus软件平台上进行.采用Johnson-Cook模型,施加与实验相同的初始条件.通过Matlab非线性最小二乘法优化程序来调整有限元模型材料参数,找到合适值,使得模拟得到的力 位移曲线与实验得到曲线相吻合.结果表明,所得参数都能很好地代表不锈钢410在高温高变形速度下的流动特性,Johnson Cook模型不能描述材料破坏阶段的行为特征.     

非对称井问题的一个解析解

运用保角变换方法寻求了一个半无限渗流域中非对称进流问题的解析解、分析了非对称井流水头势的分布规律,在此基础上,又利用井群工作时降深叠加原理,进一步研究了非对称井群协同工作情况,获得了非对称井群共同作用下的水头场解析关系,如将该解析解与数值解（有限单元法或差分法）相结合,将可实现对含有复杂密集排水系统的渗流控制问题进行快速求解,对复杂的工程渗流防渗排水设计进行优化分析研究。     

柴油机进气道全三维流动计算

用椭圆型Ｐｏｉｓｓｏｎ 方程作为控制方程的分形耦合法对形状复杂的柴油机进气道生成三维贴体网格,用时间推进有限容积法求解Ｅｕｌｅｒ方程,对柴油机进气道进行了全三维流动计算．其结果展示了进气道的结构型线对流动的影响,为进气道的设计、改进提供了理论依据．     

膜结构风致响应的分区耦合算法

采用完全隐式的分区耦合的方式,将用于计算流体流动和结构变形的软件结合起来.其中,结构和流场的几何模型分别通过有限元软件ANSYS和CFD软件建立.对流体域的模拟采用有限体积法,计算三维的、不可压缩的、湍动流动;对结构域的模拟采用有限元法,并且考虑了薄膜结构几何非线性.在每一时间步内对结构场、流体场依次求解,通过把前一个物理场的计算结果作为外荷载施加于后一个物理场来实现两个场的耦合.最后对一鞍形膜结构进行了风致响应计算,验证了这种耦合算法的可行性.     

轴向流中固支弹性薄板的大挠度流固耦合系统的数值模拟

就轴向流中两端固支大挠度弹性薄板的流固耦合振动特性,固支薄板的结构动力学方程用有限元法离散,流场采用不可压缩的二维粘性流体(N-S方程)用有限体积法离散,结合ADINA中的流体单元划分技术,建立了双向流固耦合作用下轴向流中两端固支薄板的二维仿真模型.通过模拟仿真分析研究了给定不同流速下固支板的流固耦合振动特征和大挠度系统的振动稳定性.分别得出了不同流速下固支板中点的挠度—流速曲线、挠度时程曲线及挠曲线图.结果表明:当流速小于固支板的临界流速时,板将处于稳定的直线平衡状态;当流速大于固支板的临界流速时,板将在新的位置达到弯曲平衡状态,以及在弯曲平衡位置附近发生极限环振动.