GAS-WATER TWO-PHASE DIFFUSION MODEL FOR UNSATURATED LAYERED COVER SYSTEMS

为了研究垃圾填埋场内部生物降解产生的气体(以甲烷为主) 和氧气等在非饱和覆盖层中的运移规律,建立了气体在成层非饱和覆盖层中的一维扩散模型. 该模型考虑了覆盖层含水量瞬态变化和气体扩散的耦合作用. 采用有限元多物理场耦合分析软件COMSOL Multiphysics 4.3进行求解. 其计算结果与前人实验结果能符合得很好. 研究结果表明,对典型垃圾填埋场中厚度为1m的非饱和黏土覆盖层,含水量变化对覆盖层中气体污染物的扩散运移有显著影响. 与含水量稳态变化情况相比,20 d到100 d的运移时间内同一覆盖层深度处耦合含水量作用下,气体相对浓度可以相差接近8倍. 含水量是气体在非饱和覆盖层中不稳定运移的重要影响因素.     

煤层气在非饱和水流阶段的非定常渗流摄动解

煤层甲烷由煤层的割理裂隙系统流入生产井一般经历：单相水流、非饱和流和气、水两相饱和流三个阶段，在非饱和流阶段，储层压力降至临界解吸压力之后，储存在煤基质中的吸附气体少量被解吸出来形成互不连续的气泡并阻止水的流动，含气量尚未达到饱和程度。同时煤层甲烷运移包含渗流场、变形场和应力场的动态耦合过程。本文考虑渗流过程中水-气两相不溶混流体与固体耦合作用，建立了非饱和水流阶段非定常渗流问题的流固耦合数学模型，对该强非线性一维数学模型采用摄动法和积分变换法进行解析求解，并讨论了其压力动态特性，分析了压力随饱和度S及时间t变化的规律和气相及耦合作用的影响，这些研究对煤层气、石油和天然气的开采等地下工程领域具有一定的指导意义。     

基于统计损伤理论的冻土蠕变本构模型研究

假定冻土粘塑性微元损伤符合修进的莫尔-库仑准则,损伤变量服从Weibull随机概型分布;基于热力学原理和统计损伤理论,通过推导得到了相关联流动法则下的蠕变损伤耦合本构方程.为了便于研究冻土蠕变和损伤的耦合作用,将提出的本构模型,通过用户子程序嵌入到有限元程序中.最后,对冻土中桩基承载力模型试验进行了数值模拟,结果表明:基于统计损伤理论的冻土蠕变本构模型能较好的模拟冻土蠕变变形的全过程,并且与实测蠕变变形十分吻合,该模型对冻土结构物长期稳定性分析和工程预测具有重要参考意义.     

冻土水热力耦合作用的数学模型及数值模拟

将冻土体视为空间弹性体,提出了土体在冻结过程中水分场、温度场、应力场三场耦合的一般数学模型,并给出了相应的离散方程及其解法,最后给出了数值算例,并与实测值比较,证明了该模型和算法的正确性。     

冻土的有限变形本构关系的实验研究

根据有限变形理论,给出了冻土三轴蠕变实验数据处理时所需Green应变和Kirchhoff应力的计算公式,并根据冻土三轴蠕变实验结果给出冻土的有限变形本构关系及其蠕变参数。通过对冻土的实验数据对比分析可知:在相同试验条件下,应变较小时,小应变ε1 和Green应变Ez的数值几乎相等,随着应变的增加,两种表征方式的计算结果差别越来越大,可见对于冻土如果不考虑试件变形前后长度的变化,计算得出的应变偏离实际变形情况较大,因此冻土材料的本构关系采用有限变形表征更确切。根据冻土的有限变形本构关系计算得到的冻结壁内侧最大径向位移更接近实测结果,因此冻结壁设计计算依据冻土的有限变形本构关系更为合理。对今后冻结工程的设计具有参考价值。     

非饱和粘性土中气体渗透特征

由于气体的进入形成了非饱和带,饱和土层中的气体渗透实质上是非饱和渗透问题。本文借助基于多孔介质气体渗流理论并考虑气体压缩性的、描述非饱和土中气体渗流运动的本构方程,试验研究了上海地区饱和粘性土层中的气体渗透规律。结果表明,饱水粘性土层中,气体渗透存在起始压力值,数值为100～200kPa, 与上海地区粘土的非饱和进气值较为接近;气体渗透速度与外界所施加的"压力平方差梯度"之间存在明显的分段特征。后者可能与在充满结合水的微孔隙通道中,只有一部分结合水在气压梯度的作用下发生流动,从而形成气体通道,施加的外力越大,形成的气体通道也越大,气体渗流也越明显。非饱和粘性土中气体渗透规律的研究,对于饱和粘土中的气压法施工,具有重要的理论与工程实践意义。     

非饱和介电凝胶的力学性能及力电失稳行为研究

力-电耦合场作用下,介电凝胶的体积、形状会发生相应的变化,在智能传感器与促动器等领域具有广阔的应用前景．浸入溶液中,介电凝胶将吸收溶液而溶胀．当外界溶液的量足够多时,介电凝胶将吸收足量的溶液而达到饱和状态;但如果外界溶液的量较小,没有足够的溶液可被吸收时,介电凝胶将处于非饱和状态．基于介电凝胶大变形与极化理论,对非饱和介电凝胶在等双轴应力与电场耦合作用下的力学行为及力电失稳现象进行研究．结果表明,饱和度越低介电凝胶的刚性越大,材料越硬;预应力技术可以提高介电凝胶的力电稳定性,预应力越大,力电失稳出现时的临界电压值越高．     

刚塑性圆板受水下爆炸载荷时的动力响应

分析了简支刚塑性圆板受水下爆炸载荷时的塑性动力响应。考虑到流固耦合作用，应用Taylor平板理论求水下爆炸压力，用刚塑性法分析了圆板的永久变形场，引入膜力因子考虑了圆板大变形时的膜力效应。计算结果同用ABAQUS程序的计算结果进行了比较，表明不考虑空泡时所得结果较精确。最后考察了流固耦合作用及空泡对结构响应的影响。     

基于内时理论的冻土试验研究与数值分析

从内时理论出发,基于考虑了温度和应变率影响的冻土内时流变本构方程,采用Matlab拟合法确定内时材料参数,对冻土在单轴和3轴应力条件下的应变进行了求解、分析与试验验证,并将理论结果与试验所得的应力-应变曲线进行了比较,理论曲线与试验曲线基本吻合. 与传统塑性理论的研究方法相比,内时理论的研究方法可以免去确定屈服面的困难. 基于所使用的内时本构模型,对冻土路基的冻结进行了数值分析,结果表明该理论具有很好的适用性.      

非饱和结构性黄土本构模型的研究

本文在饱和土扰动概念基础上,提出了可以考虑外力和增湿两个扰动因素对非饱和土结构性影响的耦合扰动变量。根据复合材料均匀化理论的思想建立了适用于非饱和结构性黄土的本构模型,并给出了耦合扰动变量的演化方程。通过将BBM模型和本文模型的计算结果分别与已有的实验结果进行对比说明,本文模型能够更好地反映非饱和原状黄土的力学特性。通过对结构性参数β、α进行分析说明,所给出的耦合扰动变量的演化方程能够较好地描述土结构性随变形劣化的规律。     

非饱和土-隧道系统动力响应计算的波函数法

针对非饱和地基土中埋置隧道的三维动力响应计算问题, 提出了波函数法.采用无限长的Flügge薄壁圆柱壳模拟圆形隧道衬砌,采用流、固、气组成的三相介质模拟非饱和地基土体.分别采用分离变量法以及Helmholtz矢量分解定理求解薄壁圆柱壳的振动控制方程与非饱和土的波动方程.根据隧-土交界面与地表面处的应力、位移以及孔隙流体压力等边界条件,利用平面波与柱面波的转换性质,实现了隧道内作用单位简谐载荷时隧道衬砌与土体系统动力响应的耦合求解.通过与既有单相弹性介质2.5维有限元-边界元法、两相饱和多孔介质2.5维有限元-边界元法以及三相非饱和介质Pip in Pip半解析法的计算结果进行对比, 验证了本文计算方法的可靠性. 最后,基于该方法, 通过算例分析了不同饱和度下非饱和土-隧道系统的动力响应特征.结果表明, 饱和度对土体动位移与超孔隙水压力的幅值响应有较大影响.该方法的非饱和地基土参数退化后,也可用来计算和分析饱和地基土或单相弹性地基土与隧道系统的动力响应.      

CONSOLIDATION HEORY OF UNSATURATED SOIL BASED ON THE THEORY OF MIXTURE(Ⅰ)

Unsaturated soil is a three-phase media and is composed of soil grain,water andgas.In this paper,the consolidation problem of unsaturated soil is investigated basedon the theory of mixture.A theoretical formula of effective stress on anisotropicporous media and unsaturated soil is derived.The principle of effective stress and theprinciple of Curie symmetry are taken as two fundamental constitutive principles ofunsaturated soil.A mathematical model of consolidation of unsaturated soil isproposed,which consists of25 partial differenfial equations with25 unknowns.Withthe help of increament linearizing method,the model is reduced to5 governingequations with5 unknowns,i.e.,the three displacement components of solid phase,thepore water pressure and the pore gas pressure.7 material parameters are involved inthe model and all of them can be measured using soil tests.It is convenient to use themodel to engineering practice.The well known Biot’s theory is a special case of themodel.     

土力学有效应力及其作用的讨论

讨论了饱和土有效应力的作用及其使用中应注意的问题,以及土力学预测中存在的不确定性的几种情况.从土的各相平衡方程出发推导得到了非饱和土的总应力和各相应力之间的关系, 以及非饱和土的平均骨架应力的表达式,并定义该式为非饱和土的有效应力;该表达式与变形功推导得到的有效应力公式是一致的. 文中还就非饱和土有效应力的一些问题进行了探讨.      

A Meshless Solution Technique for the Solution of 3D Unsaturated Zone Problems,Based on Local Hermitian Interpolation with Radial Basis Functions

Recent developments in meshless numerical methods have led to algorithms that can be used to solve arbitrarily large problems without the support of a connected mesh, and without the computational cost and numerical ill-conditioning issues usually associated with such solution techniques. This work applies the Local Hermitian Interpolation (LHI) method, based on local interpolation with Radial Basis Functions (RBFs), to the solution of 3D unsaturated porous media problems. The proposed implementation is capable of handling real soil properties, provided either as an analytical function or as a series of pointwise measurements. The technique is implemented with implicit and explicit timestepping, and is validated against two transient Richards’ equation models, of which one has a known analytical solution. In addition, a real-world infiltration problem based on a saturated–unsaturated formulation is modelled, using a realistic variation of soil properties with water-pressure.     

Numerical and experimental investigation of convective drying in unsaturated porous media with bound water

The heat and mass transfer in an unsaturated wet cylindrical porous bed packed with quartz particles was investigated theoretically for relatively low convective drying rates. Local thermodynamic equilibrium was assumed in the mathematical model describing the multi-phase flow in the unsaturated porous media using the energy and mass conservation equations to describe the heat and mass transfer during the drying. The drying model included convection and capillary transport of the free water, diffusion of bound water, and convection and diffusion of the gas. The numerical results indicated that the drying process could be divided into three periods, the temperature rise period, the constant drying rate period and the decreasing drying rate period. The numerical results agreed well with the experimental data verifying that the mathematical model can evaluate the drying performance of porous media for low drying rates. The effects of drying conditions such as the ambient temperature, the relative humidity, and the velocity of the drying air, on the drying process were evaluated by numerical solution.     

非饱和土-结构动力响应的多耦合周期性有限元法

真实的地基土体-隧道系统中土体及结构性质往往沿线路纵向变化.为考虑土体与结构沿纵向的变化特性,提出了一种非饱和土-结构系统动力响应分析的多耦合周期性有限元法.首先基于非饱和土的实用波动方程,采用Galerkin法推导了单节点5个自由度的非饱和土ub-pl-pg格式有限元表达式,相比于单节点9个自由度的ub-v-w格式有...     

考虑相变的近场动力学热?力耦合模型及多孔介质冻结破坏模拟

多孔介质的传热传质现象广泛存在于自然界和工业领域中. 低温条件可能导致多孔介质中的组分发生相变, 并由此诱发材料损伤, 甚至导致结构失效破坏. 对这类破坏现象的预测需要精细化建模, 以能够反映物质的相变过程和材料的破坏特征. 本文采用热焓法改写经典的热传导方程, 在近场动力学框架下, 建立了一种考虑物质相变的热?力耦合模型, 发展了交错显式求解的数值计算方法, 进行了方板角冻结、热致变形和多孔介质冻结破坏等问题的模拟, 得到了方板的冻结特征、温度场和变形场的分布规律以及多孔介质的冻结破坏过程, 与试验和其他数值方法的结果具有较好的一致性. 研究表明, 本文所建立的考虑物质相变的近场动力学热?力耦合模型能够反映材料的非局部效应和物质相变潜热的影响, 准确捕捉相变过程中液固界面的演化特征, 再现多孔介质中材料相变、基质热致变形和冻结破坏过程, 突破了传统连续性模型求解这类破坏问题时面临的瓶颈, 为深入研究多孔介质冻融破坏过程和破坏机理提供了有效途径.      

CONSOLIDATION HEORY OF UNSATURATED SOIL BASED ON THE THEORY OF MIXTURE(Ⅰ)

Unsaturated soil is a three-phase media and is composed of soil grain, water and gas. In this paper, the consolidation problem of unsaturated soil is investigated based on the theory of mixture. A theoretical formula of effective stress on anisotropic porous media and unsaturated soil is derived. The principle of effective stress and the principle of Curie symmetry are taken as two fundamental constitutive principles of unsaturated soil. A mathematical model of consolidation of unsaturated soil is proposed, which consists of 25 partial differenfial equations with 25 unknowns. With the help of increament linearizing method, the model is reduced to 5 governing equations with 5 unknowns, i.e., the three displacement components of solid phase, the pore water pressure and the pore gas pressure. 7 material parameters are involved in the model and all of them can he measured using soil tests. It is convenient to use the model to engineering practice. The well known Biot’s theory is a special case of the model.     

A hyperbolic mathematical modeling for describing the transition saturated/unsaturated in a rigid porous medium

This work proposes a mathematical model to study the filling up of an unsaturated porous medium by a liquid identifying the transition from unsaturated to saturated flow and allowing a small super saturation. As a consequence the problem remains hyperbolic even when saturation is reached. This important feature enables obtaining numerical solution for any initial value problem and allows employing Glimm’s scheme associated with an operator splitting technique for treating drag and viscous effects. A mixture theory approach is used to build the mechanical model, considering a mixture of three overlapping continuous constituents: a solid (porous medium), a liquid (Newtonian fluid) and a very low-density gas (to account for the mixture compressibility). The constitutive assumption proposed for the pressure gives rise to a continuous function of the fluid fraction. The complete solution of the Riemann problem associated with the system of conservation laws, as well as four examples, considering all the four possible connections, namely, 1-shock/2-shock, 1-rarefaction/2-rarefaction, 1-rarefaction/2-shock and 1-shock/2-rarefaction are presented.     

非饱和土力学理论的研究进展

回顾了非饱和土有效应力的发展,目前普遍认同采用两个应力变量来建立本构模型,且对基质吸力中毛细和粘吸两部分作用进行了阐述。分析了非饱和土强度问题,包括抗剪强度和抗拉强度。讨论了非饱和土的本构模型问题,包括基于净应力和基质吸力的弹塑性模型,基于Bishop有效应力和基质吸力的水力力学耦合弹塑性模型,以及双孔隙结构的模型。最后探讨了热力学方法和多孔介质理论在非饱和土中的应用,基于多孔介质理论在多场耦合条件下土体复杂的行为是当前值得研究的问题。