考虑膨胀力的非饱和介质热-水-应力耦合二维有限元分析

从建立应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵入手,使用Galerkin方法,将各控制方程分别在空间域和时间域进行离散,开发出了一个可考虑膨胀力的用于分析非饱和介质中热-水-应力耦合弹塑性问题的二维有限元程序.通过对一个假定的核废料地下处置库的热-水-应力耦合问题的数值计算,比较了无、有膨胀力时的情况,在定性上验证了该程序的正确性.     

脆性材料热-力耦合模型及热破裂数值分析方法

针对混凝土、岩石等脆性材料,利用热传导和热-力耦合的相关理论,并结合材料在细观尺度上的损伤演化规律,提出了一种考虑损伤的热-力耦合模型,并在原有材料破坏过程分析系统RFPA(Realistic Failure Process Analysis)模型的基础上建立了脆性材料热破裂过程分析的数值模拟方法.该方法考虑了脆性材料在细观层次上力学性质的非均匀性(包括强度、弹模、传导系数等),并通过统计分布函数建立了宏、细观力学性能之间的联系.对不同均匀程度材料的数值模拟结果表明:材料的非均匀性对热传导规律、热应力分布以及热破坏模式有较大的影响.材料热力学性质的非均匀性加剧了材料内部热应力分布的非均匀性,这是致使非均匀材料热破裂的一个重要因素.对稳态和瞬态热传导两种条件下的脆性介质破裂过程模拟分析表明,考虑瞬态热传导计算所得到的破裂区小于相同条件下稳态热传导所得到的结果,表明在热破裂过程分析中,应注重考虑瞬态热传导对破裂过程的影响.     

高速穿甲过程的热-结构耦合数值模拟

研究子弹侵彻炮弹药筒及发射药的高速穿甲问题及其热变化.利用有限元方法进行数值模拟,计算两种不同药筒情况下的结构及温度变化.观察整个穿甲侵彻过程,绘制应力、应力变化率、应变、应变率、温度等常规参数及可能引燃或引爆发射药的相应参数的变化云图及时间历程曲线,为发射药安全性分析和设计提供参考意见及一些有价值的数据.     

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

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

混凝土中化学-热-湿-力耦合过程的数值方法

提出了一个火灾下混凝土中化学-热-湿-力耦合过程分析的两级数学模型. 混凝 土模型化为充满两种非混溶孔隙流体的非饱和变形多孔多相介质. 数学模型基于控制干空 气、湿份及基质溶解物的质量守恒、混凝土介质混合体的动量守恒和焓(能量)守恒的耦合 偏微分方程组. 模型中特别考虑到了高温下的脱盐过程. 构造了一个用于数值模拟 化学-热-湿-力耦合行为的有限元求解过程的混合弱形式. 并且针对其中具有非自伴随算子特性的 双曲线控制方程的空间离散进行了特殊考虑. 数值结果例题显示所发展的数学模型和数值方 法在重现火灾条件下的混凝土中化学-热-湿-力耦合行为的有效性.     

荷载作用下多孔饱和地基的热-水-力耦合动力响应分析

基于广义热弹性理论,结合达西定律,对Biot波动方程进行修正,研究了一个受到荷载作用的多孔饱和地基的热-水-力多场耦合动态响应问题。建立了多孔饱和地基在荷载作用下的热-水-力耦合模型及控制方程,该模型可退化为热弹性耦合模型。采用正则模态法求解,得到了问题的解析解,讨论了热-水-力耦合模型和热弹性耦合模型的区别,分析了荷载频率变化对地基中各物理量的影响。最终给出了无量纲的竖向位移、超孔隙水压力、竖向应力和温度等物理量的分布规律。     

饱和颗粒土-维冰分凝模型及数值模拟

基于相平衡和力甲衡提出了一个关于饱和颗粒土的冰分凝模型，引入了冰分凝时未冻水迁移的驱动势，基于Clapeyron方程和吸附膜的观点简要叙述了相变区内未冻水含量与温度及孔隙水压的关系，数值模拟了饱和颗粒士的一维冻结过程，得到了冰分层分凝的计算结果，与已有实验所观测到的现象较为吻合，计算所得的冻胀量的变化趋势与已有实验所测结果是相符的．     

饱和多孔黏弹地基热-水-力耦合动力响应分析

天然土体由于沉积条件和应力状态不同, 往往会表现出一定的流变性. 本文研究地基上表面受外载荷作用时, 渗透系数和孔隙率变化对饱和多孔黏弹性地基热-水-力耦合动力响应问题的影响. 基于Biot波动方程、达西定律和Lord-Shulman广义热弹性理论, 并引入了考虑黏弹性松弛时间因子的Kelvin-Voigt黏弹性模型研究地基上表面受热/力源作用时, 孔隙率和渗透系数变化对均质各向同性饱和多孔黏弹性地基中所考虑的各无量纲量的影响. 根据不同的边界条件采用正则模态法推导出无量纲竖向位移、超孔隙水压力、竖向应力和温度的解析表达式, 结合算例分析了不同变量对各物理量的影响. 正则模态法是一种加权残差法, 可不经正、反积分变换将方程快速解耦并消除数值反变换的局限性. 结果表明: 无论何种载荷作用时, 载荷频率变化对所有考虑的物理量均有明显的影响; 孔隙率和渗透系数均对无量纲超孔隙水压力有明显的影响, 当仅考虑热载荷作用时, 孔隙率和渗透系数变化对无量纲温度均无影响. 正则模态法可广泛应用于岩土工程领域, 尤其适用于商业建筑、高速铁路和公路能源基础的热、力学特性研究中. 该研究结果可为工程施工奠定一定的理论基础, 具有一定的指导性意义.      

基于热-磁-弹性耦合的二维导电薄柱壳数值分析

基于几何方程、物理方程、运动方程和电动力学方程,建立了二维导电薄柱壳热磁弹性基本方程.考虑到Joule热效应,引入热平衡方程及广义Ohm定律,得出导电薄柱壳的温度场.利用变量代换方法,整理成具有10个基本未知量的标准型方程组.采用差分法和准线性化方法,给出准线性微分方程组.对于二维导电薄柱壳,导出了Lorentz力表达式、温度场积分特征值.通过实例计算,得到了二维导电薄柱壳应力、位移、温度随外加电磁参量的变化规律.研究结果可为二维导电薄壳热磁弹性问题研究提供理论参考.     

双重孔隙介质THMM耦合模型及其有限元分析

建立了一种饱和-非饱和遍有节理岩体的双重孔隙-裂隙介质热-水-应力-迁移耦合模型,其特点是应力场和温度场是单一的,但具有不同的孔隙渗流场、裂隙渗流场和孔隙浓度场、裂隙浓度场,以及可考虑裂隙的组数、间距、方向、连通率和刚度对本构关系的影响;并开发了相应的二维有限元程序.针对一个假定的高放废物地质处置库,就岩体为非饱和双重孔隙-裂隙介质和放射性核素泄漏的情况进行了数值分析,考察了岩体中的温度、负孔隙水压力、饱和度、地下水流速、核素浓度和主应力的状态.结果显示:孔隙和裂隙中的负水压力及核素浓度呈现不同的变化、分布;尽管裂隙水饱和度平均仅为孔隙水饱和度的1/10,但因裂隙的渗透系数比孔隙的渗透系数大4个数量级,故裂隙中地下水的流速约是孔隙中相应值的3倍;孔隙和裂隙中核素浓度的量值接近.     

含损伤的冻土本构模型及耦合问题数值分析

从复合材料的细观力学机理出发,建立了含损伤的冻土弹性本构模型. 对于在不同冰体积 含量和不同温度下的冻结砂土,由该损伤本构模型计算的结果与实测的应力-应变曲线比较 吻合. 利用自主开发的程序,分别对渠道冻结和冻土路基的水分场、温度场、应力场3场耦 合进行数值模拟,得到了比较准确、详尽的符合实际的温度场与应力场、位移场、应变场耦 合的计算结果,与前人的计算及实测结果相符合,规律一致. 表明该程序能够计算冻土材料 的相关物理量,并能很好地描述这些相关物理量之间的关系.     

切削加工过程数值模拟的研究进展

近年来,有限元方法在切削加工模拟中得到了越来越广泛的应用,在研究切削工艺参数及切屑成形机理方面有着不可替代的作用。本文介绍了国外切削加工过程有限元数值模拟的研究进展,阐述了切削过程有限元模拟的关键技术,包括切屑形成、切削加工中的热力耦合、工件与刀具接触和摩擦、切屑分离和断裂准则以及工件残余应力、残余应变的模拟等技术;最后,还对切削工艺有限元数值模拟的发展方向作了探讨。      

混凝土断裂过程的力学模型与数值模拟

混凝土是一种典型的非均匀准脆性材料,其断裂过程是非常复杂的.评述了混凝土断裂以及断裂过程研究的进展状况,主要总结了一些研究混凝土断裂的宏观力学方法和模型,介绍了从混凝土的细观层次出发模拟其断裂过程的数值模型.最后,简单介绍了作者用细观数值方法模拟混凝土宏观断裂过程的研究成果.      

非局部模型与变形局部化数值模拟

基于非线性身体场论,建立了材料非局部连续模型、变分方程及相应的实时拖带系大变形有限元数值模型,设计了这一模型的数值卷积算法,由于广义函数弱收敛定理和卷积理论,证明所提出的非局部连续模型具备收敛性和稳定性。并阐明了材料特征尺度数学物理意义,统计加权函数的选择原则,数值结果表明非局部模型描述变形局部化问题是适当的。     

格构细胞机模型在岩石劈裂过程中的数值模拟研究

基于弹塑性损伤理论,结合格构细胞机以及格构模型的优点,提出了格构细胞机模型,用于模拟岩石的劈裂破坏机理并取得了较好的结果.     

NUMERICAL SIMULATION FOR A PROCESS ANALYSIS OF A COKE OVEN

A computational fluid dynamic model is established for a coking process analysis of a coke oven using PHOENICS CFD package. The model simultaneously calculates the transient composition, temperatures of the gas and the solid phases, velocity of the gas phase and porosity and density of the semi-coke phase. Numerical simulation is illustrated in predicting the evolution of volatile gases, gas flow paths, profiles of density, porosity of the coke oven charge, profiles of temperatures of the coke oven gas and the semi-coke bed. On the basis of above modeling, the flow of coke oven gas （COG） blown from the bottom of the coke oven into the porous semi-coke bed is simulated to reveal whether or not and when the blown COG can uniformly flow through the porous semi-coke bed for the purpose of desulfurizing the semi-coke by recycling the COG. The simulation results show that the blown COG can uniformly flow through the semi-coke bed only after the temDerature at the center of the semi-coke bed has risen to above 900℃.     

注塑模充模过程的数值分析

利用非等温、粘性广义Hele－Shaw流动理论,建立了三维薄壁型腔塑料熔体充填过程的数学模型．耦合利用有限元／有限差分法数值求解压力和能量方程,实现充模过程的动态模拟．     

多煤层开采覆岩破断过程的模型试验与数值模拟

首先分析了研究区煤矿山西组煤层剖面和水文地质特征,针对多煤层联合开采的特点和覆岩的工程地质特征,采用工程地质力学模型实验、数值模拟计算相结合的综合研究方法,分析了多煤层开采的采动影响及岩层动态断裂机理,得出了有关岩层移动参数和多层煤同采时应力分布状态,计算得到了多煤层开采垮落带与导水裂隙带的发育高度分别为32m和81.5m,导水裂隙带高度影响范围已经达到风化带,未形成切冒,局部穿透粘土层。