地震条件下挡土墙主动土压力及其分布的新分析法

采用旋转挡土墙计算模型的变换法,将在地震和拟静力法条件下主动土压力的求解问题转化为在静力条件下主动土压力的求解问题。根据在静力条件下水平层分析法的主动土压力推导结果,直接获得在地震条件下主动土压力强度分布、土压力合力及其作用点位置的表达式,并运用图解法得到了临界破裂角的解析解。公式可考虑水平和垂直地震加速度、不同墙背倾角、墙背和坡面倾角与填料存在黏结力和外摩擦角、存在均布超载等诸多因素的影响,公式可以适用于在常用边界和地震条件下黏性土的主动土压力计算。旋转地震角法是将在地震和拟静力法条件下挡土墙计算模型旋转为在静力条件下挡土墙计算模型,但旋转挡土墙计算模型并不改变挡土墙和墙后填土的应力状态,按在静力条件下挡土墙主动土压力求解方法求解在地震和拟静力法条件下主动土压力,该方法大大简化了在地震和拟静力法条件下的主动土压力计算公式推导过程,统一了在拟静力法条件下的地震土压力求解,理论更加完善。     

临近地下室外墙影响下的考虑土拱效应的挡土墙主动土压力研究

当挡土墙附近存在临近建筑地下室外墙时,其挡土墙土压力与传统的Rankine理论基于无限半空间体假定不符,因而在这种新的工程背景下需要采用合适的理论来计算挡土墙土压力及其作用点高度。已有的研究表明,这种条件下土体的变形趋势可分为上、下两大部分：上部土体变形类似于Terzaghi的活动门试验,土体沿着墙体下滑,而下部土体则沿着土楔形体而变形。因而将土拱效应用于求解挡土墙土压力的计算分成了上、下两大部分考虑。假定土拱形状为圆弧,基于主应力旋转概念分别给出了上、下两部分的侧向土压力系数,运用水平微分层析法基于静力平衡思想给出了两部分的水平向主动土压力分布公式。最后通过坐标平移的方式给出了主动土压力合力及其作用点高度的表达式。算例表明,计算结果与数值计算结果较为接近,其结果对实际工程有一定的参考价值。     

考虑土拱效应的挡土墙主动土压力与被动土压力统一解

土拱效应对倾斜挡土墙下的主动土压力及被动土压力有重要的影响,但是相关计算理论研究略显不足。为了将土拱效应考虑到倾斜挡土墙下的土压力计算中,首先通过应力摩尔圆及静力平衡法分别给出了考虑土拱效应下主动土压力及被动土压力计算所需的两大因素：侧向土压力系数及竖向平均应力公式。在此基础上建立了考虑土拱效应的倾斜挡土墙主动土压力及被动土压力的统一表达式,并将其应用到求解土压力合力及其作用点高度的计算中。算例表明,土拱效应对于主动土压力与被动土压力的影响不同。随着墙体倾角的增大,主动土压力作用点高度逐渐降低,即土拱效应随着墙体倾角的增大而降低。与前述相反,随着墙体倾角的增大,被动土压力作用点高度逐渐降低,即土拱效应的影响随着墙体倾角的增大而增大。     

地震作用下挡土墙主动土压力及转动位移分析

分析地震引起的挡土墙位移及墙后土压力,对于评估挡土墙可靠性具有重要意义。基于拟动力法,考虑时效、地震波传播的相位差、超载、墙背摩擦角、填土黏聚力以及填土开裂等影响,建立地震作用下挡土墙主动土压力计算模型,获得挡土墙绕墙趾转动模式下主动土压力大小、分布形式及作用点高度。同时,考虑挡土墙本身受地震荷载作用的影响,求出挡土墙绕墙趾的转动位移。通过与Mononobe-Okabe法对比可知,文中获得的主动土压力值与Mononobe-Okabe法接近,但Mononobe-Okabe法低估了主动土压力作用点高度,表明采用Mononobe-Okabe法设计存在风险。通过算例分析了地震系数、墙背摩擦系数、超载大小、时间、填土黏聚力和内摩擦角对挡土墙转动位移的影响。     

地震条件下黏性土挡土墙土压力分析

Mononobe-Okabe理论是现阶段计算地震土压力的常用方法,但Mononobe-Okabe理论的诸多假设使其具有一定的局限性。针对Mononobe-Okabe理论的不足,考虑到地震作用下挡土墙偏转对土压力的影响,采用斜向条分法推导了复杂条件下黏性土地震土压力强度分布、土压力合力及其作用点位置公式,并利用图解法给出了临界破裂角的解析解。研究表明：填土黏聚力和地震系数对土压力影响显著;忽略黏性填土表面开裂与地震作用对均布超载及开裂填土等效超载的影响将使主动土压力计算结果偏小,其误差随着填土黏聚力和均布超载的增大而增大;在不同水平地震系数下土压力沿墙高呈非线性分布;所提公式适用范围更广,有效完善了Mononobe-Okabe理论。     

临近地下室外墙影响下的挡土墙空间土压力研究

目前关于临近地下室外墙影响的挡土墙空间土压力的计算理论的研究还比较少,原有的平面应变条件下的理论不能满足挡土墙的长高比B/H较小时的挡土墙土压力计算要求。通过将土拱效应原理引入顾慰慈[8]建立的空间土压力计算模型,建立了考虑土拱效应的临近地下室外墙影响的空间土压力计算模型,根据挡土墙和地下室外墙的间距与土体破裂面状态的关系将该模型分为3种情况,并将各模型划分为Ⅰ、Ⅱ、Ⅲ、Ⅳ四个区域,通过在各个区域内取水平微分单元体,建立各微分单元体的水平和竖向静力平衡方程,推导出了各区相应的挡土墙空间主动土压力计算公式,该公式可以计算出墙背任意位置的主动土压力;并提出了空间土压力合力及其合力作用点的计算方法。通过算例计算可以直观地看出挡土墙后主动土压力的空间分布,由此可以看出,当空间效应存在时,考虑土拱效应的挡土墙主动土压力沿墙长的分布与平面应变条件时有很大的不同,此时挡土墙两端附近区域的主动土压力远小于平面应变条件下计算出的主动土压力,同时可以看出,考虑空间效应的挡土墙主动土压力合力作用点要比平面应变条件下的位置要高,挡土墙长高比B/H越小,空间效应对主动土压力沿墙长的分布和主动土压力合力作用点位置的影响越大。     

地震条件下挡墙后黏性土主动土压力研究

采用水平层分析法,得到了地震条件下挡墙后黏性土主动土压力合力和作用点位置、土压力强度分布以及临界破裂角的解析解。公式考虑了水平和垂直地震加速度、挡墙墙背倾角、填料黏聚力和内摩擦角、填料与墙背的黏结力和外摩擦角、均布超载等因素,并分析了这些因素对主动土压力的影响。结果表明,朗肯和库伦理论下的主动土压力公式以及Mononobe-Okabe主动土压力公式与地震条件下的主动土压力公式完全一致。地震条件下的主动土压力强度沿墙高呈非线性分布。水平地震加速度增大了主动土压力,垂直地震加速度使得主动土压力有所减小     

挡土墙土体复合破裂面主动土压力分析

基于普朗德尔-瑞斯纳课题认为,挡土墙土体破裂面是由对数螺线与直线组成的复合滑动面的结论,采用极限平衡法和约束优化相结合的方法,重新推导了主动土压力及其合力和作用点高度的计算公式。在此基础上分析了主动土压力在纵向的分布特征,并讨论了主要参数对破裂面、主动土压力及其合力和作用点的影响。     

地震条件下悬臂式挡墙主动土压力的极限分析方法

为确定地震条件下悬臂式挡土墙主动土压力,考虑假想坦墙墙背的可能不同位置,给出了墙后填土5种可能的失稳破坏模式;在此基础上,采用拟静力法,基于极限分析上限定理,推导了作用于坦墙墙背上的地震主动土压力计算公式,包括填土性质、填方坡面倾角、踵板长度、墙体高度、水平及竖向地震影响系数等多因素,其中除填土黏聚力与竖向地震影响系数与该土压力呈线性相关性外,其余因素呈非线性影响。实例分析表明,基于本方法地震土压力而计算的墙体抗滑与抗倾稳定系数,多数情况下均比经典的Mononobe-Okabe法略偏大;在填土中存在第二破裂面情况下,以踵板下边缘作为假想墙背端点的计算模式相对略偏不安全;竖直假想墙背模式相应的土压力计算值最小,但相应的墙体稳定系数却不一定最大。     

RT模式下地震非极限主动土压力的拟动力分析

为考虑挡墙位移效应对地震土压力的影响,依据前人试验研究的结论,将摩擦角表示为与挡墙位移量和位置高度相关的函数,然后基于拟动力法和水平层分析法,推导得出RT位移模式下的地震非极限主动土压力和合力作用点的计算表达式。计算模型可描述摩擦角沿着墙高逐渐发展的不同非极限位移状态工况,并建立了挡墙位移、地震动荷载和土压力之间的相互联系。参数分析讨论了振动时间、挡土墙位移状态、地震加速度参数和土体摩擦角对地震主动土压力分布、合力大小以及合力作用点高度的影响。相比于传统的极限状态地震土压力理论,所提方法更合理地描述了地震土压力随挡墙位移的发展过程,对发展非极限土压力理论和改进边坡工程中的抗震计算方法具有一定的参考意义。     

Acknowledgement for Reviewers for 2021

Mathematical Geosciences -     

Framework for soil suitability evaluation for sewage effluent renovation

Current methods of establishing suitable locations for onsite wastewater treatment systems (OWTS) are inadequate, particularly in light of the numerous cases of onsite system failure and the resulting adverse consequences. The development of a soil suitability framework for assessing soil suitability for OWTS allows a more practical means of assessment. The use of multivariate statistical analysis techniques, including Principal Component Analysis (PCA) and multi-criteria decision aids of PROMETHEE and GAIA, enabled the identification of suitable soils for effluent renovation. The outcome of the multivariate analysis, together with soil permeability and drainage characteristics permitted the establishment of a framework for assessing soil suitability based on three main soil functions: (1) the ability of the soil to provide suitable effluent renovation, (2) the permeability of the soil, and (3) the soils drainage characteristics. The developed framework was subsequently applied to the research area, Gold Coast, Queensland, Australia, and the use of standard scoring functions were utilised to provide a scoring system to signify which soils were more suitable for effluent renovation processes. From the assessment, it was found that Chromosol and Kurosol soils provided the highest level of effluent renovation, closely followed by Ferrosol and Dermosol, Kandosol and Rudosol soil types. Tennosol and Podosol soil types were found to have a significantly lower suitability, with Hydrosol soils proving the least suitable for renovating effluent from OWTS.     

湘南锡矿找矿方向探讨

湘南地区蕴藏着丰富的锡矿资源,锡及锡多金属矿矿床(点)分布广泛,类型复杂,它们的形成除与花岗质侵入岩体直接相关外,同时还受到区域地层和地质构造等因素的共同制约。通过对锡矿的地质背景、分布特征、矿床特征与矿化规律、矿床成因及找矿信息和标志的研究,把湘南地区锡矿的找矿归纳为:"划带、圈区、寻体、定位"的找矿模式,不仅从宏观上提出了或解决了锡矿找矿方向的问题,而且近几年来的地质找矿工作得到了实践应用,并取得了丰硕的找矿成果。     

A new approach for calculating strain for particulate media

Discrete element modelling is a viable alternative to conventional continuum‐based analysis for analysing problems involving localized deformations of particulate media. However, to aid in the interpretation of the results, it is useful to express the results of discrete element analyses in terms of the continuum parameters of stress and strain. A number of homogenization methods have been proposed to calculate strain in discrete systems; however, two significant limitations of these methods remain. First, none of these methods incorporate particle rotation effects satisfactorily, although significant particle rotation occurs in shear bands in both physical tests and numerical simulations of granular materials. Additionally, observations of the particle displacement fields in shear bands in granular materials indicate that the displacements within the localizations are erratic. Consequently, existing linear, local interpolation approaches produce substantial variations in the strain values calculated in adjacent elements in the region of localization, hindering clear visualization of the strain localization as it evolves. A new method of domain discretization for calculating strain is proposed. This method is capable of capturing particle rotation and employs a non‐local meshfree interpolation procedure capable of smoothing the erratic displacements in strain localizations, which better defines their evolution. The proposed method is validated for problems involving both two and three dimensions. A number of methods are compared with the proposed method and pertinent insights are made. Copyright © 2003 John Wiley & Sons, Ltd.     

ISA-Hypoplasticity accounting for cyclic mobility effects for liquefaction analysis

Acta Geotechnica - The hypoplastic model for sands proposed by Wolffersdorff (Mech Cohes Frict Mater 1: 251–271, 1996) combined with the intergranular strain anisotropy by Fuentes and...