用电网模型研究二相二维非定常渗流的一种近似方法

二相非定常渗流理论近年来得到了很大的发展。在忽略毛细管压力情况下的一维问题,原则上都已能解决。对二维问题,由于在数学上存在困难,目前还没有得到分析解。文献[3]作了在数字电子计算机上直接求二相二维问题的数值解的尝试。文献[4]在渗流过程中流管形状不变化的假设下提出了一个解决二相二维非定常渗流的近似方法,但是,这个假设只是在特定的条件下才适用,在一般情况下,流管的形状却可能随时间显著地变化。本文把已有的二相一维渗流理论和电模拟方法结合起来,提出了一个在理论上较严格、在实用上也可行的研究二相非定常渗流二维问题的近似方法。     

峰后岩石非Darcy渗流的一维非定常严格简明解析解

煤炭是中国近期的主要能源,仍需要大力研究. 而煤矿采动围岩大多处于峰后应力状态或破碎状态,其渗流一般不符合Darcy定律. 探讨非Darcy渗流系统, 对其研究既有理论创新价值,尤其在煤矿安全中更有重要的实用价值. 用近年发展的求解偏微分方程新的分离变量法------加法分离变量法,得出了Ahmed-Sunada型非Darcy渗流的3套非常简明的一维非定常严格解析解,以发展相应的渗流理论,以及作为标准解推进渗流数值计算的水平.      

????????·??????????

 阐述两种典型的计算机代数系统(CASes)------Matlab和Maple在基础 力学课程与有限元课程教学中的一些应用. CASes可以方便地将力学教学和有限元计算过程 中的模型建立、推理及其数值计算等进行演示和计算，为基础力学和有限元课程教学提供 了新的手段和途径.     

非线性渗流力学研究进展

非线性渗流力学、多重介质渗流力学、计算渗流力学(包括优化控制理论)以及渗流物理学,是现代渗流力学领域的四个重点发展方向。我国近年来多重介质渗流力学发展较快,在某些方面接近国际水平或略为领先。计算渗流力学也做了不少工作。可是,近10年来与生产实际需要关系密切的大部分内容都和非线性渗流有关。这方面国外发展很快。这应引起我国流体力学界重视,结合我国生产实际需要,将一部分科研力量转移到这个新的生长点上。      

弱不连续问题高阶有限元离散系统的GAMG法

弱不连续问题(如含夹杂问题)是固体力学计算中的一类重要问题。高阶有限元方法由于其具有更好的逼近效果,是确保数值解在界面保持较高精度的计算方法之一。但与线性元相比,高阶单元需要更多的计算机存储单元,具有更高的计算复杂性。本文利用两水平算法的思想,将高阶有限元离散系统化归于线性元离散系统的求解,为弱不连续问题高阶有限元离散系统设计了一种新的基于几何与分析信息的代数多重网格(GAMG)法,并应用于圆形求解域含单夹杂问题的高阶有限元离散系统的求解。数值试验结果表明,相比于常用GAMG法,新方法的迭代次数基本不依赖于问题规模、单元阶次以及杨氏模量的间断性,CPU计算时间得到明显改善,具有更好的计算效率和鲁棒性,可大大提高弱不连续问题有限元分析的整体效率。     

采空区流场非达西渗流一种新的迭代算法

对Bachmat非线性渗流方程求解问题，结合采空区问题的力学特点，提出基于变渗透率达西 (Darcy)渗流求解非线性渗流的迭代算法. 建立了冒落采空区非线性漏风流态的有限元数 值模型，通过对复杂边界采空区的漏风渗流的计算，得到与实际流态更接近的流动规律(风 压等值线和流函数线)和采空区漏风强度分布(速度场). 结果表明，迭代算法是振荡性收 敛的，收敛速度快；与Darcy渗流相比，采空区非线性渗流速度场趋向于平缓，在工作面附 近30\,m范围内最大风速降低到原来的0.61倍，其余大部分区域速度略有增加. 迭代方法满 足工程要求.     

裂缝性低渗透油藏流-固耦合理论与数值模拟

根据裂缝性低渗油藏的储层特征,建立适合裂缝性砂岩油藏渗流的等效连续介质模型。将渗流力学与弹塑性力学相结合,建立裂缝性低渗透油藏的流-固耦合渗流数学模型,并给出其数值解.通过数值模拟对一实际井网开发过程中孔隙度、渗透率的变化以及开发指标进行计算,并和刚性模型以及双重介质模型的计算结果进行了分析比较.     

功能梯度板的非线性动力分析

非线性材料功能梯度板件的动力分析是属于在数学方程上同时具有变系数、非线性、非定常特征的固体力学问题.文中首先将问题的变系数非线性偏微分方程组转化为各向异性常系数非线性常微分方程,然后用小参数法求得解析解,适用于各种形状、边界及功能梯度分布的板件非线性弹性振动分析.     

局部采用标志质点的流体网格法

本文提出了一个解具有大形变的二维非定常的多物质流体力学问题的数值方法,在采用流体网格法思想的同时采用了标志质点,在不同的物质的界面附近,借助于标志质点计算流量;对连续流的区域,用流体网格法计算流量。计算检验表明,与质点网格法相比,减少了数值波动。     

二维对流扩散方程的欧拉—拉格朗日分裂格式

本文在[1]基础上发展了一种有效的处理大P_e(R_e)数、非定常二维对流扩散方程的欧拉-拉格朗日(E-L)分裂格式,由于方法本质上与区域形状无关,且不需再分网格,因此是一种无网格的E-L方法,特别对于定常流动,E.-L.分裂格式可以导致比一阶迎风格式更精确的单调、无振荡格式,文中对于常系数、变系数和非线性的二维非定常和定常对流扩散方程的(初)边值问题进行了数值计算,数值结果与精确解的比较表明,本方法具有很好的精度,解是单调无振荡的,比通常一阶迎风格式具有较少的数值扩散,最大计算网格P-e(R-e)数可达100—500。     

谈工科理论力学教学中数学工具的应用

分析了当前理论力学教学由于缺乏对数学工具的充分应用而带来问题,提出在理论力学教学 中应充分利用矢量、矩阵等数学工具,重视对学生计算机数值求解方法的训练. 在具体的求 解过程中选用MATLAB或Maple软件作为数值计算的工具.     

低渗透裂缝性油气藏非稳态窜流因子研究简

深入分析了常定窜流因子和现有的非稳态因子对于低渗透基质-裂缝窜流的不适应问题,根据对23 个低渗岩心渗流曲线的拟合分析结果,引入了变渗透率系数,并在此基础上推导了非线性扩散方程,也就是低渗透基质-裂缝系统间窜流的控制微分方程;经过无量纲化处理,引入了动边界条件,分别利用积分方法和矩方法导出了非线性扩散方程在前后两个阶段的近似解析解,并在此基础上构建了非稳态窜流因子的表达形式;新的窜流因子可以有效应用于低渗透基质-裂缝系统间的窜流计算,有限元数值模拟的结果验证了新因子在计算非线性非稳态窜流方面的准确性与可靠性.     

低渗透裂缝性油气藏非稳态窜流因子研究

深入分析了常定窜流因子和现有的非稳态因子对于低渗透基质-裂缝窜流的不适应问题,根据对23 个低渗岩心渗流曲线的拟合分析结果,引入了变渗透率系数,并在此基础上推导了非线性扩散方程,也就是低渗透基质-裂缝系统间窜流的控制微分方程;经过无量纲化处理,引入了动边界条件,分别利用积分方法和矩方法导出了非线性扩散方程在前后两个阶段的近似解析解,并在此基础上构建了非稳态窜流因子的表达形式;新的窜流因子可以有效应用于低渗透基质-裂缝系统间的窜流计算,有限元数值模拟的结果验证了新因子在计算非线性非稳态窜流方面的准确性与可靠性.      

有限元双向分配法解无侧移刚架及其微机程序

本文的目的在于给出了一种新的有限元方法,用这种方法我们得到一个无侧移刚架渐近数值解. 我们编制了APPLE—Ⅱ型微机FORTRAN-Ⅳ程序.并给出了二个算例的数值结果.     

加权残数法用于固体力学近年来进展概况

有限元法自本世纪50年代发轫以来,迄今已有1/4世纪。大量的计算实践证明,它是一种应用广泛、能处理复杂问题的卓有成效的数值计算方法。它的重大贡献在于,以前在数学、力学、结构或非结构的各个技术领域中看来是无法获解的问题,现在几乎都能解决,例如塑性力学问题,广泛的非线性问题,复杂结构的动力问题等。目前有限元方法已发展得相当完备,清醒地重新估计,无可否认地有限元法也有比较重大的缺点。首先这是一个很费钱的方法。其次这个方法在一定情况下是不准确的。E.L.Wilson也指出:“现在      

低渗透多孔介质渗流动边界模型的解析与数值解

考虑启动压力梯度的低渗透多孔介质非达西渗流模型属于强非线性动边界问题, 分别利用相似变量变换方法和基于空间坐标变换的有限差分方法, 对内边界变压力情况下、考虑启动压力梯度的一维低渗透多孔介质非达西渗流动边界模型进行了精确解析与数值求解研究. 研究结果表明:该动边界模型存在唯一的精确解析解, 且所求得的精确解析解可严格验证数值解的正确性;且当启动压力梯度值趋于零时, 非达西渗流动边界模型的精确解析解将退化为达西渗流情况下的精确解析解. 由求解结果作出的非零无因次启动压力梯度下的地层压力分布曲线表现出紧支性特点, 其与达西渗流模型的有显著不同. 因此, 研究低渗透多孔介质中非稳态渗流问题时, 应该考虑动边界的影响. 研究内容完善了低渗透多孔介质的非达西渗流力学理论, 为低渗透油气藏开发的试井解释与油藏数值模拟技术提供了理论基础.      

机械荷载与热荷载共同作用下变厚度梁的热弹性力学解

本文研究了两端简支变厚度梁受机械荷载与热荷载共同作用下的热弹性力学解.温度场用调和级数展开,通过求解热传导方程可首先确定温度沿梁厚度方向的非线性分布情况.从二维热弹性力学理论的基本方程出发,导出满足控制微分方程和两端简支边界条件的位移函数的一般解,对上下表面的边界方程作傅立叶正弦级数展开确定待定系数,数值结果与商业有限元软件ANSYS进行了比较,显示出很高的精度.本文方法可直接应用于对应力和位移分析要求较高的工程问题,如航空航天和微型机械的设计.     

非线性计算力学的某些新发展——参加第二次国际非线性力学计算方法会议后的述评

七十年代以来,由于计算技术的发展,在理性连续介质力学研究获得进展之外,形成了非线性连续介质的数值分析,或非线性计算力学的重要领域,这是现代连续介质理论和现代数值分析方法相互作用的产物。J.T.Oden教授的著作“非线性连续介质的有限元”在一定程度上给出了这一研究的理论基础,也反映了这一研究的发展趋向。     

求解三维弹性力学问题的快速耦合方法

将多种数值方法耦合,充分利用各种方法的优点建立新的数值方法,是求解三维复杂问题的有效途径之一．本文将无单元Galerkin (Element-Free Galerkin, EFG)方法、有限元法和维数分裂法耦合,提出了求解三维弹性力学问题的快速耦合方法(Fast Hybrid Method, FHM)．将三维弹性力学问题分裂为若干个二维平面问题,对于每个二维问题采用罚函数法施加边界条件,并推导其相应的积分弱形式,引入Shepard基函数的移动最小二乘法建立形函数,进而推导二维平面问题的离散方程．第三个方向上采用有限元法将这些二维离散方程进行耦合,可以得到原三维弹性力学问题的快速耦合方法数值解的求解公式．通过数值算例验证了本文快速耦合方法求解三维弹性力学问题的收敛性,将数值解与解析解对比,说明了本文方法求解三维弹性力学问题的有效性．     

非饱和多孔介质中热-渗流-力学耦合的混合元法

提出了一个非饱和多孔介质中热-渗流-力学耦合分析的混合有限元 方法. 固相位移、应变和净应力；孔隙水和气的压力、压力空间梯度和Darcy速度；多相混 合介质的温度、温度空间梯度和热流量在单元内均为独立变量分别插值. 基于胡海 昌-Washizu 三变量广义变分原理给出的多孔介质中热-渗流-力学耦合问题控制方程的单元弱形式，导 出了单元公式. 采用共旋公式进行几何非线性分析. 数值结果证明了所提出的单元模拟以 应变局部化为特征的渐进破坏的能力