基于热力学的临界状态模型的返回映射算法及实现

相对于其他临界状态模型,基于热力学的临界状态（TCS）模型不需要引入塑性势假设,能够自动满足热力学定律。通过对TCS的修正,使其能够模拟初始K0固结的影响,利用了返回映射算法进行了TCS模型的ABAQUS本构二次开发,通过与ABAQUS内嵌的修正剑桥模型（MCC）计算结果的对比,证实了该程序的可靠性。对模型中的2个控制屈服面形状的参数进行了讨论,分析了它们对应力-应变关系与剪胀关系的影响,修改TCS模型参数可以实现非椭圆的屈服面,进而拓展了模型的适用范围,不同的参数对于屈服面形状和大小的影响也有所不同。同时还比较了采用考虑K0固结及旋转硬化的TCS模型与不考虑K0固结及旋转硬化的MCC模型所得到的应力-应变关系与剪胀关系的差异。对于真实土体,K0固结及旋转硬化是土体的基本力学特性。最后证明了TCS模型较MCC模型（模型中没有考虑旋转硬化与K0固结对硬化规律的影响）可以很好地模拟该特性,所以更加有效。     

抗拔桩极限平衡方程及其应用

将抗拔桩侧阻力分解为与桩侧正压力不相关的桩-土黏结强度 、与桩侧有效正压力成正比的摩擦力 两部分,采用摩擦定律计算摩擦力 。基于轴对称条件,假定土体为半无限弹性体,以Mindlin公式积分计算分析极限平衡状态下桩-土共同作用,依据平面应变条件下柱状孔扩张的弹性力学解建立桩-土界面位移协调方程,推导出抗拔桩极限平衡方程,给出了求解方法及计算参数确定方法。该方程能反映桩与土的材料特性、桩体尺寸、桩顶埋深、群桩效应、卸荷效应等多因素对抗拔桩极限承载力的影响。结合海上风电大直径超长抗拔钢管桩足尺试验进行验证。对比分析结果表明,该方法计算的抗拔极限承载力与实测值接近,计算精度远高于现行规范推荐方法,其结果可为工程应用及抗拔桩承载力机制研究提供参考。     

单一粒径组颗粒材料破碎规律及破碎极限研究

岩土材料颗粒破碎演化规律是颗粒破碎过程复杂性的集中体现,对单一粒径组颗粒材料破碎规律的研究是进行多粒径组颗粒材料破碎规律的基础。从概率统计角度入手,先对前期研究中所提出的描述单一粒径组颗粒材料破碎规律的3个参数-破碎概率p、破碎状态参数a、b的演化规律进行了研究。通过对3组不同类型（不同试验材料、试验类型）的试验数据的研究,发现随着破碎的发展3个参数趋于定值,表明破碎极限的存在。基于上述研究,最后提出了计算单一粒径组颗粒材料破碎极限的方法,分别得到了3组试样的最终分形维数分别为2.3、2.34、2.28。     

基于扰动状态概念的结构性土压缩特性分析

由于土结构性的影响天然沉积土的结构是亚稳定的,屈服后的压缩变形阶段必然伴随着结构的破损。将附加孔隙比?e作为表征原状土的一个结构状态参数,基于扰动状态概念,引入定量分析附加孔隙比?e随固结压力变化关系的方法,得到一维扰动状态概念(DSC)压缩模型,以描述结构性土压缩损伤现象,并通过结构破损指数b来刻画压缩过程中的结构破损速率。根据该模型,对太湖湖沼相典型的天然沉积软黏土、粉质黏土和硬黏土不扰动试样的一维压缩试验结果进行模拟和分析,并结合试验数据提出了模型参数的测定方法。分析结果验证了该模型能够描述具有不同结构形式土样的压缩特性及其实际变形过程中表现出的非线性,这给结构性土非线性固结与沉降计算提供了一定的理论基础。     

考虑位移效应的土压力计算理论

基于朗肯土压力理论,假定填土的内摩擦角与该点土体位移呈非线性关系,进而提出挡土墙的主动和被动土压力的计算模式。该模式随位移变化是连续的,且能考虑墙土间的摩擦的影响。又将其与砂土模型试验结果进行了对比分析,吻合较好,从而证明用该计算模式计算其主动或被动土压力是合理的。此外,根据该计算模式,导出了一种有效估算其静止土压力系数的计算方法。     

基于旅游市场竞争态及其转移模型的入境旅游市场研究——以广西为例

以广西入境旅游市场为对象,引用旅游市场竞争态模型,采用定性定量相结合的方法,分析广西入境旅游主要客源市场与区内各城市目的地市场的竞争态。结果表明,港澳台、马来西亚、越南、日本占据着广西入境旅游客源市场的主导地位,市场发展水平不平衡;区内各城市目的地市场发展水平不合理,瘦狗市场与幼童市场数量众多,市场仍然处于开发潜力阶段;广西入境旅游市场基本上遵循了市场竞争态转移的规律和方向,但脆弱性依然明显,必须从根本上推动旅游业改革,实现旅游业的可持续发展。     

A comparative simulation study of coupled THM processes and their effect on fractured rock permeability around nuclear waste repositories

This paper presents an international, multiple-code, simulation study of coupled thermal, hydrological, and mechanical (THM) processes and their effect on permeability and fluid flow in fractured rock around heated underground nuclear waste emplacement drifts. Simulations were conducted considering two types of repository settings (1) open emplacement drifts in relatively shallow unsaturated volcanic rock, and (2) backfilled emplacement drifts in deeper saturated crystalline rock. The results showed that for the two assumed repository settings, the dominant mechanism of changes in rock permeability was thermal–mechanically induced closure (reduced aperture) of vertical fractures, caused by thermal stress resulting from repository-wide heating of the rock mass. The magnitude of thermal–mechanically induced changes in permeability was more substantial in the case of an emplacement drift located in a relatively shallow, low-stress environment where the rock is more compliant, allowing more substantial fracture closure during thermal stressing. However, in both of the assumed repository settings in this study, the thermal–mechanically induced changes in permeability caused relatively small changes in the flow field, with most changes occurring in the vicinity of the emplacement drifts.     

Influence of interaction between small asperities on various types of slow earthquakes in a 3-D simulation for a subduction plate boundary

Recently, the occurrence of slow earthquakes such as low-frequency earthquakes and very low-frequency earthquakes have been recognized at depths of about 30 km in southwest Japan and Cascadia. These slow earthquakes occur sometimes in isolation and sometimes break into chain-reaction, producing tremor that migrates at a speed of about 5–15 km/day and suggesting a strong interaction among nearby small asperities. In this study, we formulate a 3-D subduction plate boundary model with two types of small asperities chained along the trench at the depth of 30 km. Our simulation succeeds in representing various types of slow earthquakes including low-frequency earthquakes and rapid slip velocity in the same asperity, and indicates that interaction between asperities may cause the very low-frequency earthquakes. Our simulation also shows chain reaction along trench with propagation speed that can be made consistent with observations by adjusting model parameters, which suggests that the interactions also explain the observed migration of slow earthquakes.     

莫尔包络线对数模型方程的建立与应用

从研究岩石破坏机理着手,结合室内相关岩石力学试验,建立莫尔包络线对数模型方程,并对其进行参数率定,得出应用范围、变化趋势和岩石的破坏机制,简化岩石破坏强度的研究,可推广到其它的岩石力学实验中。     

金沙江某水电站上下坝址岩体卸荷差异研究

在河流下切、斜坡成坡过程中,斜坡应力场不断调整和适应坡体的变化,表浅部岩体在应力变化过程中往坡外发生卸荷,在同一地区岩体卸荷深度应当相差不大。而在西南某水电站相隔仅1.4km上下坝址勘察过程中发现,两坝址岩体卸荷深度却存在较大差异,上坝址岩体卸荷深度在20m以上,最大深度达到90m;下坝址岩体卸荷深度仅为5~20m。,通过现场调查及数值模拟成果表明,上下坝址岩层方位是两坝址岩体卸荷差异的决定性因素,上坝址为顺向坡,岩层走向与边坡走向平行,倾角较陡的岩层形成拉裂面,易于风化营力的侵入,岩体卸荷深度大;下坝址为横向坡,岩层走向与边坡走向近于垂直,岩体强度远大于岩层面的强度,且平行坡向的结构面数量少,岩体卸荷深度在20m以内。