古建筑风化石质构件力学参数的确定方法

针对目前常规确定岩石力学参数的方法只能得出岩块的整体强度,不能反映古建筑石质构件风化层内力学参数逐渐变化,无法很好地满足石质古建筑稳定性、耐久性研究的不足,提出了综合运用现场穿透波CT测试、表面回弹测试和室内试验确定风化石质构件力学参数的方法。方法的实施主要包括：通过现场穿透波CT测试,获得纵波波速在石质构件深度剖面上的分布;现场选取与欲研究的石质构件类似的岩石做室内力学试验和波速测试,得出力学参数与波速的统计关系;上述两者结合得出力学参数在石质构件深度剖面上的分布,在此基础上分析风化层内力学参数随深度的变化关系;对不同风化分区代表性石构件进行风化深度研究和表面回弹测试,得出构件表面强度比与风化深度的关系。以义乌宋代古月桥为例,对其方法进行了系统的介绍,并对古月桥承载条石风化层内力学参数随深度的变化关系以及表面强度与风化深度的关系进行了研究。研究表明：义乌宋代古月桥承载条石风化层内的抗压强度比、弹性模量比随深度呈负指数关系变化;表面强度比与风化深度呈2次多项式关系。     

弹黏塑性模型及其在预压处理软基沉降计算中的应用

以准确计算预压荷载作用下软土路基的工后沉降为目的,分析了不同预压法加固软土地基时土体的应力变化。提出了用弱化应力路径变化过程的弹黏塑性模型表示软黏土的应力-应变关系。结合经典比奥固结理论,推导了预压荷载作用下土体沉降计算的有限单元法。通过具体算例,验证了所提方法。结果表明：（1）不同预压方法的加固机制和应力路径均不相同。堆载预压法加固软基是一不等向的正压固结过程,真空预压法则是在负压作用下的等向固结过程;（2）弹黏塑性模型是不受瞬时施加应力的大小及具体应力路径变化过程影响的软土流变模型,运用于预压荷载作用下的软基时具有不考虑应力路径和不划分主次固结的优越性;（3）在比奥固结理论中,用弹黏塑性模型反映软土的物理方程时,可计算软基考虑了弹黏塑变形性质的固结沉降和工后沉降;（4）弹黏塑性模型的软土体黏性参数? /V和塑性参数? /V对竖向位移都有比较明显的影响。     

考虑发育薄弱区对管涌发展影响的数学模型研究

为了揭示发育薄弱区对管涌发展的影响机制,进行了管涌砂槽模型试验和颗粒流数值模拟,建立了3组模型,以研究孔隙通道的空间发育对颗粒运移过程的影响,分析了土体内部细颗粒运移引起管涌破坏的动态过程。根据数值模拟结果,分析了管涌发展过程中试样细观及流体的变化规律,并与砂槽模型试验结果进行了比较,模拟与模型试验结果较吻合。结果表明,薄弱区越发育,压力分布调整的趋势越强,颗粒运移的趋势相应增加,达到破坏需要的时间越短。数值模拟获得的临界水力梯度与室内试验以及理论分析结果均比较接近,证明了采用颗粒流程序研究管涌问题是合理的,对研究存在薄弱区的砂土管涌现象具有明显的优势和广阔的应用前景。所揭示的结果有益于管涌机制的更深入研究。     

含双裂隙岩石裂纹演化机理的离散元数值分析

采用离散单元法探讨了预制双裂隙岩石的裂纹演化机理。用近期从试验资料提取的无胶结厚度含抗转动能力的岩石微观力学模型和相应的离散单元法商业软件,模拟了含不同预制倾角的双裂隙岩石试样在单轴压缩作用下裂纹的扩展与贯通规律,揭示了裂纹演化的宏微观机理。同时,将离散元法DEM岩石试样的裂纹的扩展和贯通规律以及强度特性与室内试验结果进行了比较分析。结果表明,预制裂隙之间以及端点处的拉应力集中是导致裂隙岩石破坏的主要原因,且DEM数值试验得到裂纹的演化规律与室内试验结果较为一致。含30°的预制裂隙的岩石试样最容易起裂,含75°的预制裂隙的岩石试样最困难起裂,造成此种现象的原因可能是裂纹在垂直于主应力方向上的长度不同导致试样受拉区域大小不同。     

模拟三维裂纹问题的多尺度扩展有限元法

扩展有限元法模拟裂纹时独立于网格,因此该方法是目前求解裂纹问题最有效的数值方法。为了在计算代价不大的情况,实现大型结构分析中考虑小裂纹或提高裂纹附近精度,在裂纹附近一般采用小尺度单元,其他区域采用大尺度单元。提出了分析三维裂纹问题的多尺度扩展有限元法,在需要的地方采用小尺度单元。基于点插值构造了六面体任意节点单元。所有尺度单元都采用8节点六面体单元,这样六面体任意节点单元可方便有效地连接不同尺度单元。采用互作用积分法计算三维应力强度因子。边裂纹和中心圆裂纹算例分析结果表明,该方法是正确和有效的。     

无网格局部Petrov-Galerkin法大地电磁二维正演

将无网格局部Petrov-Galerkin算法用于大地电磁二维正演。介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法详细推导了与之对应的局部Petrov-Galerkin弱式方程,并用高斯积分法将其离散化。论述了无网格局部Petrov-Galerkin法较无单元Galerkin法及有限元法的优缺点,最后通过二维模型的计算验证了算法的有效性。      

深厚钙质黏土层冻结特征与冻结温度场数值模拟

为了解决杨村煤矿井筒在深厚黏土层层位掘砌时遇到的冻结管断裂难题,分析了施工中产生的问题,并采用ANSYS有限元软件模拟出该黏土层段温度场分布规律,细化导热系数区间,使计算结果与测温孔数据接近。结果表明,强化冻结工艺和采用短段掘砌施工是解决深厚黏土层安全施工的有效措施,对日后井筒在该层位的施工具有一定参考价值。      

安徽宁国地区岩-土系统元素迁移及其农业地质环境效应

通过多条基岩-土壤剖面和区域性土壤地球化学样品的微量元素丰度测试,分析了宁国地区不同基岩类型地球化学元素在岩石和土壤中分布的基本特点,表明同一地层的岩-土系统在不同地点有基本相同的元素迁移规律,导致土壤中的微量元素分布特征受控于其下地质单位的空间分布,植物生长必需的微量营养元素Fe、Mn、Cu、Mo、Zn、B6种元素表现出在南华纪-寒武纪地层中向土壤强烈富集的趋势．岩-土系统元素迁移特点对宁国市天目山区重要名优特产品-山核桃的生长环境产生重要影响,该区主导性名特优产品-山核桃,应该选择南沱组、兰田组、杨柳岗组、华严寺组、西阳山组等地层出露的地段种植才能产生好的经济效益．     

地壳块体动应力变化特征及其与地震关系的初步研究

在运用三维非线性动态有限元方法模拟华北块体对渤海、 黄海潮汐水位变化的动态响应的基础上,结合有关地学资料,综合研究了地壳块体动应力变化特征,分析其与地震关系,得到了对地壳动应力特征的新认识:① 动应力在传播变化过程中,能够通过应力波的干涉叠加呈现高应力值,与动压应力相比动剪应力的这种效应更为明显;② 动应力的变化影响孕震体应力环境及地震的发生;③ 从动应力变化特征看, 动剪应力叠加变化过程反映了孕震过程.       

采用显式有限元方法求解开放系统动力响应的主要环节和问题

分析了采用显式有限元方法求解开放系统动力响应的主要环节中影响精度和稳定性的因素及其相互关系,指出数值解的精度和稳定性是这些因素共同作用的结果。归纳了合理输入数字荷载应满足的条件以及系统离散化准则。基于有限元方法推导了吸收边界的一般表达式,并据此提出吸收边界和透射边界之间存在确定的函数关系,并不是两类独立人工边界的新见解。认为内域的阻尼稳定性和人工边界的近场稳定性是大型复杂开放系统动力响应数值方法研究和应用中的关键科学问题,并从应用的角度提出了解决问题的思路和看法。