饱和土三维非轴对称Lamb问题

基于积分变换方法提出位移组合积分变换和应力组合积分变换式,成功地求解了Ｂｉｏｔ两相介质理论的饱和土三维非轴对称问题动力方程,得到了以土骨架位移和孔压为基本未知量的积分解。在此基础上系统地研究了饱和土的Ｌａｍｂ问题。考虑了表面排水和不排水两种条件,得到了三维非轴对称饱和弹性半空间在表面竖向和水平力作用下,表面径向位移、竖向位移和周向位移的积分形式解,并将解答准确地退化到弹性半空间经典Ｌａｍｂ问题解,     

横观各向同性磁电弹性体的基本解和边界积分方程

为了进行磁电弹性体的边界元法的计算, 需要相应的基本解. 为此从横观各向同性磁电弹性体的控制方程出发, 应用微分算子代数运算, 将相应于4个特征根互不相等情况下的通解用5个拟调和函数表达, 进而, 将用试凑法构造的拟调和函数代入通解, 推导得到了磁电弹性无限体的基本解. 利用这个基本解建立了横观各向同性磁电弹性体的边界积分方程, 并对两个算例进行了边界元法的计算, 与解析解作了比较, 吻合很好.     

压电材料平面问题的基本解

从横观各向同性压电介质平面问题基本方程出发，推导得玛组在体积力作用下的一般解，利用这组一般解结合体积势理论及积分的方程，得到了无限平面法集中力和点电荷作用下的位移和电势的有限形式的表达式，从而给出了边界元法中有用的基本解。     

考虑非线性多孔弹性介质热流固耦合动力响应模型

基于Biot波动方程及热弹性动力理论,Darcy定律及修正的Fourier热传导定律,考虑流固两相介质的压缩性,渗流速度对热扩散的影响及温度梯度对渗流速度的影响,建立了饱和多孔弹性介质热流固完全耦合的动力响应模型（THMD）及控制方程,该模型可退化为基于传统Fourier热传导定律的热弹性动力模型（TMD）.最后,利用Laplace变换分析了地表作用有随时间变化热/力源条件下的一维热动力响应问题,讨论了反映土体缩性的Biot模量、热水渗透系数对结果的影响,并与热弹性动力模型结果进行比较,验证了结果的正确性.     

非均匀饱和土中平面波的传播特性

基于Biot多孔介质理论,假设材料物理力学特性沿厚度方向连续变化,利用回传射线矩阵法(RRMM),建了非均匀饱和土平面波动问题的弥散方程,其中考虑了土体的非均匀性、惯性、粘滞以及固体颗和流体的可压缩性.通过数值算例,分析了非均匀饱和土中的弹性波的传播特性,结果表明:两类纵波的波数和耗散沿材料特性变化的方向变化不大,而非均匀性对横波的波数和耗散有着明显的影响.     

轴对称层合横观各向同性压电圆板的精确解：Ⅱ.弹性？…

用与压电圆板精确解相同的方法，求得横观各向同性圆板在刚性滑动和弹性简支两种边界条件下的三维精确解。对于各向同性圆板，处理了特征根为生根时的问题，最后给出了压电圆板和弹性圆板的若干数值结果。讨论了经典板理论的适用性。     

轴对称层合横观各向同性压电圆板的精确解——II.弹性圆板精确解和数值计算

用与压电圆板精确解相同的方法 ,求得横观各向同性圆板在刚性滑动和弹性简支两种边界条件下的三维精确解 .对于各向同性圆板 ,处理了特征根为重根时的问题 ,最后给出了压电圆板和弹性圆板的若干数值结果 ,讨论了经典板理论的适用性 .     

EDA介质中地震波速度、衰减与品质因子方位异性研究

通过考察EDA介质中应力应变关系 ,建立了EDA介质中考虑非弹性效应的地震波动方程组 .其中 ,裂隙引起的非弹性效应体现在波动方程组中位移分量关于时间t的三次导数项上 .通过平面波位移函数解 ,研究了EDA介质中地震波速度、衰减与品质因子Q值的方位各向异性 ,并分别对干裂隙与饱含水裂隙两种情形下EDA介质中地震波速度、衰减与Q值随波频率与传播方向的变化关系进行了探讨 .研究结果说明 ,裂隙的存在以及裂隙内所含物质内容对地震波速度、衰减与Q值及其各向异性特别是两类qS波衰减或品质因子比值具有重要的影响 ;低频时 ,EDA弹性效应占主导地位 ;高频时非弹性效应愈加不可忽略 .     

孔隙和裂缝分形油藏的瞬态压力分析

研究了孔隙与裂缝分形油藏中流体的流动．通过引进流动路径长度的微分关系式和单位横截面面积得出分形渗透率、孔隙度和渗流速度三个基本公式．阐明了分形孔隙介质与分形裂缝介质中渗透率和孔隙度表达式的区别．建立起分形油藏中单相流体渗流的压力扩散方程,给出相关方程的解,绘制了用于瞬态压力分析的典型曲线．讨论了孔隙分形油藏和裂缝分形油藏试井解释的方法和步骤．对于多维各向异性分形介质,在直角坐标和圆柱坐标系中也给出了相应的渗流压力扩散方程．对分形油气藏勘探和开发具有重要意义．     

页岩气直井非稳态非线性渗流的数值计算及产能预测

结合我国南方海相页岩气藏的实际,针对页岩气藏的特殊孔隙结构及多尺度流动的特征,首先由Knudsen数判断出储层主体流态为滑移流,并指出页岩气在页岩多孔介质中的传输过程是一个等温过程.然后修正了Beskok-Karniadakis二阶近似模型,建立了致密页岩气直井平面径向非稳态非线性渗流问题的数学模型.进一步应用玻尔兹曼变换将问题简化,并给出了求解压力场的数值离散显式迭代格式.最后数值计算得到了内边界定压条件下直井的压力场随时间、空间分布曲线,预测了产量随时间变化的规律,并分析了解吸、滑移和扩散效应对产量的影响.     

非饱和土中功的表述以及有效应力与相分离原理的讨论

基于多相孔隙介质理论的质量、线动量、能量平衡方程和非饱和土力学的假定，推导出总变形功W的一种表达式。在总变形功的具体方程基础上给出了建议与土骨架位移在功上对偶的有效应力表达式．它和Wheeler（2003年）以及其他一些研究者根据直觉和经验给出的有效应力具有相同的形式，但本文给出的有效应力是从多相孔隙介质理论的基本原理出发，经推演而得到的，具有坚实的理论基础和科学依据。文中根据总变形功的具体表达式，指出非饱和土的物理力学性质复杂，其影响因素较多，其固体骨架的变形不可能像饱和土那样由有效应力唯一确定。另外，利用总变形功的表达式以及各相自由能A^a的关系式，论证并得出了Passman（1984年）给出的相分离原理不严格适用于非饱和土的结论。     

炭／炭复合材料微观孔隙结构演化的分形特征

提出了一种用分形理论研究炭／炭（C／C）复合材料制备过程中孔隙演化特征的新方法．基于压汞测试数据,根据海绵分形模型及热力学关系模型导出的多孔介质分形维数计算公式,计算了致密化各阶段C／C复合材料孔隙的分形维数,并研究了分形维数随孔隙演化过程的变化规律．结果表明C／C复合材料属于多孔分形介质,孔隙分形维数随孔隙率的减小而增大,但同时受到热解炭织构形态的影响,从各向同性到高织构热解炭,分形维数减小．分形维数综合反映了C／C复合材料内部孔隙的复杂程度和热解炭的形貌特征,是监控C／C复合材料致密化过程中孔隙演化的一个有效参数．     

EDA介质中地震波速度,衰减与品质因子方位异位研究

通过考察ＥＤＡ介质中应力应变关系,建立了ＥＤＡ介质中考虑非弹性效应的地震波动方程组。其中,裂隙引起的非弹性效应体现在波动方程组中位移分量关于时间ｔ的三次导数项上。     

基于Reissner理论的硬夹心夹层板的弯曲研究

传统的夹层板理论都是将夹心视为软夹心,不考虑其面内应力分量和抗弯刚度.本文在修正了Reissner理论的基础上,研究了硬夹心夹层板的弯曲问题,根据最小势能原理得到了具有3个广义位移的弯曲基本方程和边界条件,并对其进行了简化.同时还研究了硬夹心夹层板弯曲问题的有限元格式,推导了4节点四边形等参单元单元刚度矩阵和等效节点载荷列阵,并给出了矩形填充多孔夹层板的ANSYS实体模型解,计算结果表明位移的有限元解收敛于解析解且具有足够的精度,与ANSYS实体模型解也有较好的一致性.     

基于偶应力理论涂层半平面的二维滑动摩擦接触分析

本文基于偶应力理论,研究了平面应变状态下刚性压头与均匀涂层半平面之间的二维滑动摩擦接触问题.该理论通过引入材料特征长度来描述微结构材料的尺寸效应.采用傅里叶积分变换方法将此尺寸依赖的接触问题转化为第二类Cauchy奇异积分方程.然后,通过数值求解该类方程来确定涂层表面法向应力和面内应力的分布.详细讨论了尺度参数、剪切模量比和摩擦系数对表面法向接触应力和面内应力的影响.结果表明基于偶应力理论得到的接触应力明显偏离基于经典弹性理论预测的结果,而且表现出对材料特征长度的强烈依赖性.     

多孔介质内火焰传播与驻定机理的研究

采用多孔介质内的一维反应流模型和详细的化学反应机理GRI3.0,数值模拟多孔介质内的预混燃烧过程。通过对能量方程中各项大小的考察,以及对火焰传播速度随火焰位置、相间对流换热系数、固体导热系数、消光系数等参数的敏感性分析,深入浅出地阐述了多孔介质中浸没火焰和表面火焰的传播和驻定机理。结果表明,多孔介质中浸没火焰的传播机理与有预热的层流预混自由火焰传播机理相似,表面火焰传播机理与反应区有热损失的层流预混自由火焰传播机理相似。     

不同渗透性饱和土中三类体波的传播性质研究与参数分析

饱和土动力问题中,渗透系数决定了流-固耦合作用的强弱,是影响饱和土中体波传播特性的重要参数.针对渗透系数为有限值、零及无穷大3种情况,基于Zienkiewicz提出的u-w-p模型,推导获得了P1,P2和S波的波动方程及相应的名义波速.在此基础上,得到势函数表示的三类体波的波动方程的解,并结合解中各项的物理意义,给出了饱和土中波的实际传播波速及衰减系数.最后,逐一讨论了渗透系数、荷载频率、孔隙率和固相可压缩性系数对波速及其衰减特性的影响.     

波浪作用下海床响应的解析解

利用分离变量法,将Biot固结模型的求解问题转化为复系数线性齐次常微分方程的求解,进而引入复波数的概念,构造水平位移的通解,最终得到了Biot固结模型的显式通解.将本文给出的解析解与Yamamoto等人的试验结果、解析解及Hsu和Jeng的解析解进行了比较,验证了本文解析解的有效性和优越性.基于解析解,通过对海床深度影响问题的分析,得出了波浪动水压力在海床内的影响深度为一倍波长的结论.     

多孔介质发动机压燃着火的数值研究

多孔介质发动机能够实现发动机内的均质和稳定燃烧．为加深对多孔介质发动机特性的了解,用改进的KIVA-3V对一种形式的多孔介质发动机燃用气体燃料的工作过程进行了模拟,并讨论了多孔介质初始温度和多孔介质结构及燃料喷射时刻对多孔介质发动机实现压燃着火的影响．针对Zhdanok等人的甲烷预．空气预混合气过滤燃烧的实验,用改进的KIVA-3V进行了数值计算以检验模型的合理性,燃烧波的计算结果与实验结果吻合得很好．计算结果表明在压缩比一定时,多孔介质初始温度是决定多孔介质发动机能否压燃着火的重要因素;改变多孔介质的结构会影响多孔介质内气固两相的换热和多孔介质的弥散作用,多孔介质孔隙大小是影响多孔介质发动机压燃着火的主要因素;在上止点附近起喷燃料不能实现多孔介质发动机的压燃着火．     

基于格子Botlzmann方法的三维多孔介质中多相导热过程的研究

建立了模拟三维多孔介质中多相导热过程的格子Boltzmann模型,结合三维多孔介质重构算法,从孔隙尺度对CPL／LHP毛细芯中的瞬态导热过程进行了模拟,计算了不同条件下毛细芯内温度分布,分析了热负荷、孔隙率等因素对导热过程的影响,并计算了三维多孔介质的有效导热系数．