基于动力学和材料软化特性的边坡渐进破坏特征研究

边坡的失稳是一个从量变到质变的动态渐进破坏过程,此问题也是边坡领域研究的重点与难点之一。在考虑岩土材料软化特性和动力学求解的基础上,建立了边坡渐进破坏仿真的理论框架;利用ABAQUS软件的动力显式求解模块实现了边坡的渐进破坏仿真;根据塑性应变揭露了剪切带的扩展过程,由软化本构确定了滑面材料的分区演化规律,根据等效塑性应变确定了边坡的滑面,通过滑面位置将边坡分为滑体、滑带、滑床,并分别研究了边坡各分区内部特征点运动学变量的发展过程,从而揭示了边坡的渐进破坏过程;基于材料参数沿滑面的时空分布,利用矢量和法得到了边坡不同演化阶段的安全系数。对比该方法与Bishop法确定的滑面位置与安全系数,发现两种方法峰值和残余强度对应的安全系数比较接近,该方法搜索所得滑面位于Bishop法自动搜索的滑面之间,验证了此方法的合理性及可靠性。最后分析了材料软化特征对边坡稳定性的影响,在保持其他参数不变的条件下,增大残余黏聚力,边坡的滑面位置加深,安全系数的初始值减小,安全系数的快速减小阶段有所推迟,并且快速减小阶段经历的时间有所延长,稳定后的安全系数有所增大。保持其他参数不变,增大残余黏聚力对应的等效塑性应变阈值,边坡的滑面位置加深,安全系数的初始值减小,安全系数的快速减小阶段有所推迟,但快速减小阶段经历的时间基本不变,达到稳定的时间有所推迟,同时稳定后的安全系数略微有所增大。     

基于应变软化的多阶边坡稳定分析

多阶边坡除了最危险滑面需要关注,每级台阶的潜在滑面及其他不满足规范要求的滑面亦需要得到重视。以双台阶边坡和复杂多台阶边坡为例,基于应变软化理论,利用FLAC3D软件分析边坡的剪应变发展规律,根据剪应变的集中带获取边坡的多个潜在滑面。分析滑面的剪入、剪出和滑面中部点的强度参数变化规律,并结合矢量和法分析滑面在应变软化过程中安全系数的演化规律。研究结果表明：采用该方法可获取边坡的最危险滑面及次级滑动面;随着软化的发展,滑面上点的强度参数从峰值强度逐渐降低至残余强度,然而演化的速度不同;随着应变软化过程的发展,所有滑面的安全系数均逐渐降低,最后趋于一个定值。     

基于矢量和分析方法的边坡滑面搜索

基于力是矢量这个基本概念,给出了边坡矢量和法安全系数的矢量表达式。在边坡荷载已知的情况下,运用有限元法计算得到边坡真实的应力分布状态,使用矢量和法安全系数并采用优化算法中的直接法搜索边坡的临界滑面。计算结果表明,使用矢量和分析法搜索得到的临界滑面位置与极限平衡法得到的滑面位置基本吻合,且安全系数相对误差最大为2.09 %,验证了该方法在边坡稳定分析中的合理性与可靠性。矢量和分析方法与基于强度储备安全系数的分析方法不同,具有明确的物理和力学意义,且不需引入过多假设,公式简洁,易于进行三维边坡稳定分析,不需进行迭代计算,可在工程实际中得到广泛的运用。     

边坡局部和整体稳定性评价方法的辩证统一

边坡稳定性评价方法有局部和整体两类,两者之间存在差异有其历史的局限性和存在的合理性。以极限平衡法和有限元强度折减法为例,分析了典型整体与局部稳定性评价方法在滑面搜索、平衡条件、安全系数的定义、材料参数和失稳等5个方面存在的异同。利用PFC算例揭示了真实边坡带有强烈空间和时间特征的渐进破坏特性,发现边坡失稳是自然地形的重塑、再造和稳定的过程,也是一个动势能转换和耗散的过程。针对边坡稳定问题及其特点,提出了基于广义软化的边坡稳定矢量和法的理论构架和技术路线,并进行了程序实现,克服了传统强度折减脱离实际材料劣化特征的不足,安全系数求解仅与应力状态相关,假设条件少,可以实现局部与整体稳定评价的辩证统一,可以反映边坡失稳的渐进特性和参数演化的时空效应。与传统极限平衡法(LEM)和强度折减法(SRM)在滑面和安全系数方面的综合对比发现,局部稳定揭示的滑面更稳定,获取的安全系数是全周期的,文中方法所得安全系数比整体方法的结果稍高,差异是整体方法的假设条件造成的。     

考虑张拉-剪切渐进破坏的边坡矢量和法研究

针对边坡渐进破坏过程中岩土体强度参数不断劣化现象,采用应变软化模型代替传统的理想塑性模型,提出一种考虑张拉?剪切渐进破坏的边坡矢量和分析方法。从安全系数表达式与最危险滑动面搜索两方面进行研究,首先基于滑动和力的矢量特性,改进矢量和法滑动趋势与安全系数表达式,进而基于改进自适应遗传算法搜索边坡最危险滑动面。算例表明,采用应变软化模型时的矢量和法安全系数与最危险滑动面形态具有动态变化与调整的特性。该基于张拉-剪切应变软化模型的计算方法进一步扩展了边坡矢量和分析法。     

基于瑞典条分法的应变软化边坡稳定性评价方法

引入简化的应变协调方程,基于瑞典条分法的基本原理,提出一个新的能够考虑土的剪应力-应变关系的应变软化边坡稳定性简化分析方法。推导了该方法的计算公式,提出了具体算法并编制了计算程序,分析了不同应变软化特性的边坡稳定性。实际应用表明,该方法简单可靠、效率较高,计算所采用的分条数量和初始滑裂面位置对计算结果精度的影响均很小。应变软化边坡的安全系数介于基于峰值强度和残余强度的极限平衡条分法计算结果之间,具体数值不仅取决于土的抗剪强度指标,还与其剪应力-应变关系密切相关。     

基于聚类分析的边坡稳定性研究

在边坡稳定性分析中,滑面搜索方法对于得到边坡的稳定性系数至关重要。现今两种常用的边坡稳定性分析方法,极限平衡法以及强度折减法,可以得到极限应力状态下边坡的潜在滑面。而为了评价边坡在当前应力场下的稳定性,提出一种基于聚类分析概念的滑面搜索方法。这种方法综合考虑了边坡各点的位移以及点安全系数,采用K均值聚类法将具有相近性质的点进行分类,把边坡体分为潜在危险区以及稳定区。潜在危险区的边界即为边坡在当前应力状态下的潜在滑面,并运用矢量和方法计算得到其稳定性系数。采用上述方法对3个算例进行验算,并与极限平衡法以及强度折减法进行对比。结果表明：（1）基于聚类分析的边坡稳定性评价方法具有较好的合理性和可靠性;（2）当边坡处于临界状态附近,通过聚类分析得到的潜在滑面与其他两种方法得到的滑面基本一致;当边坡处于较稳定状态时,当前应力状态下的潜在滑面要比强度折减法得到的滑动带范围要更深更广,但求出的稳定性系数也较高;当存在软弱夹层时,矢量和计算出的滑面更为陡峻,稳定性系数比其他两种方法的结果要小。     

基于岩土材料软化特性的滑坡多级滑动面分析方法研究

在实际工程中,滑坡经常呈现出多级破坏的特征,而一般计算和设计中仅关注最危险的滑动面及对应的最小安全系数,这往往遗留安全隐患。在考虑岩土材料软化特性的基础上,借助FLAC3D和Matlab软件平台构建了无需人工干预便可有效模拟与评价滑坡多级破坏的理论框架。以旬阳县党家坝廉租房小区滑坡为例,通过对塑性剪应变、塑性拉应变、剪切应变增量等特征变量的渐进发展规律分析发现：滑坡多级滑动面的形成,其时间和空间顺序并不一定相对应。各级滑动面在时间上按产生的先后顺序依次为第1级主滑面、第2级主滑面、次级滑面;而在空间上按照从前到后的顺序依次为第1级主滑面、次级滑面、第2级主滑面。从计算所得的滑坡最终破坏形态中发现,坡体内部出现的滑动面条数与现场采集的拉裂缝数基本保持一致;第1级主滑面的入口位置与滑坡前缘拉裂缝位置基本吻合;而第2级主滑面及次级滑面的位置与现场勘查到的拉裂缝位置出现偏差。从材料参数的渐进发展规律中发现,滑面上强度参数的分布及大小均随着滑面的产生、发展而逐渐变化,其变化区间为峰值强度到残余强度,这是此方法可以有效模拟多级滑动面形成过程与各级滑动面之间相互影响的核心所在。借助矢量和法成功地实现了基于滑面上强度参数渐进发展规律的安全系数演化过程的确定。通过对安全系数演化过程的分析发现,在滑坡多级滑动面的形成过程中安全系数大小顺序呈现出3种不同的状态,这很好地揭示了各级滑动面在形成过程中的主、被动关系。     

LU Kun-lin,ZHU Da-yong,YANG Yang 均质边坡准三维安全系数实用计算曲线

针对三维极限平衡法计算边坡安全系数存在着分析过程繁琐、较难普及推广等诸多不足,建立了一套实用的三维边坡安全系数计算曲线。假设三维滑面为中间圆柱和两端为部分球面的组合滑面,采用基于滑面正应力修正的极限平衡分析方法,计算36 864个均质边坡,根据计算结果绘制该套曲线。依据该套曲线可以得到,（1）根据已知边坡体的相关指标快速获得三维安全系数,能够得到较二维分析方法更准确、更客观的稳定性评价;（2）根据实际坡体的稳定状况反分析获得滑面土体强度参数,有效地克服二维反分析得到滑面土体强度参数偏高的弊端。研究成果可为均质边坡三维稳定性初步评价提供简单、实用的参考依据。     

求解边坡矢量和安全系数的条分法

矢量和安全系数物理、力学意义明确,通常采用有限元计算得到的应力场来求解。提出一种简便的求解方法,在通用条分法（GLE）的基础上,满足力和力矩平衡的严格条件,推导了基于条分法的矢量和安全系数。采用广义简约梯度法（GRG）搜索临界滑面,在极限抗滑合力矢方向上求解矢量和安全系数。整个计算过程都在电子表格软件Excel内实现,易于在工程实践中推广。对于圆弧滑面边坡,其计算结果表明：①矢量和安全系数与代数和安全系数非常接近,对应的滑面位置也相近。②不同的条间力函数假设（常数和半正弦函数）对边坡的矢量和安全系数影响很小。③极限抗滑合力矢方向? 小于边坡坡角,其值与假定条间力平行时的条间力方向非常一致。在极限平衡分析框架内,计算方法可进一步简化为：假设条间力平行,将条间力方向作为矢量和安全系数的计算方向。对于非圆弧滑面边坡,其计算结果表明：①均质边坡的矢量和安全系数与代数和安全系数相对误差小于3 %,滑面位置接近,而软弱层边坡两种方法的计算结果有一定差别。②不同的条间力函数假设对边坡的矢量和安全系数影响很小。     

土压缩模量预测回归分析方法

在土压缩模量Es的分析计算中,常用的分析方法是根据土的e-p曲线试验资料建立公式求解,本文引入了数学上的回归分析方法来探讨Es与其它各项物理力学指标参数之间的相关关系,旨在找出对Es影响最显著的指标参数,并建立最优回归方程式,以此来预测Es。研究结果表明:该方法通用性强,应用范围广,亦可推广应用于其它相关领域。     

Geochemistry of Thermal Waters of Continental Margin of Far East of Russia

Studied waters belong to warm(T=30-50℃),alkaline(pH=8.9-9.3),low mineralized(TDS235 mg/1)Na-HCO_3 or Na-SO_4-HCO_3 thermal waters with high content of SiO_2(up to 81 mg/l)and F(up to 3.9 mg/1),occur on modern volcano-tectonic rejuvenated areas of Eastern Sikhote-Alin orogenic belt.Low~3He concentration as well as N_2/O_2 and N_2/Ar ratios exclude influence of deep mantle fluid.New rare earth element data constrain our understandmg of water-rock interaction occurring in the water source region.Meteoric origin of waters is proved by stable isotope values varying from-71‰to-136.1‰and from-10.8‰to-18.8‰forδ~2U andδ~(18)O respectively.REE patterns reflect high pH,resultfing from water-rock interaction and oxidative conditions.Calculations of deep aquifer temperature using Na-K and quartz geothermometers show 116.8-131.1°C and 82.2-125.8℃respectively.Presence of deep faults both with abnormal thermal gradient(~45-50 K/km)define unique geochemical shape of thermal waters of Sikhote-Alin,area,where no present volcanic activity is registered.     

岩石坚硬程度对核电厂地基-基础地震响应特征的影响分析

考虑核电厂地基-基础的动力相互作用,应用显式动力有限差分法分析了地震作用下极软岩、较软岩、坚硬岩上核电厂建筑结构基础的地震响应特征,比较了岩石坚硬程度对基础加速度反应谱的影响。研究表明：随着岩石坚硬程度的提高,核电厂建筑物结构基础的地震响应有增加的趋势;在周期轴上,基础处的加速度反应谱曲线会随着岩石坚硬程度的提高逐渐向短周期（高频段）方向移动。在高频段,建造于较坚硬岩石上基础结构的加速度反应谱值偏大;在中等频段,建造于较软岩石上基础结构的加速度反应谱值偏大;在低频段,岩石坚硬程度对加速度反应谱的影响不显著。     

Causes of salinization of the Gulf of Taganrog

Using the database of automatic hydrometeorological stations, installed in the Don RIver delta and Taganrog Bay seashore, the sources of the anomalois scale water negative surge and salinization of the Azov Sea under conditions of low river flow in 2015–2016 are studied. The new schemes of stratification and advection of salty sea waters in the Don River mouth under different weather conditions, water discharge and levels are given.     

电信网络两类逻辑图自动生成算法

为解决电信手工绘制的逻辑图不易维护的问题,研究了电信带方向逻辑图和不带方向逻辑图的自动生成算法.算法是从数据库提取缆线的路由信息,用深度优先的方式遍历缆线经过的所有坐标,采取回溯的方式,计算缆线在每个点的偏移值,调整缆线顺序,自动绘制逻辑图.产生的逻辑图没有不必要的交叉,避免了传统的手工绘制.结果表明,对于常用的电信逻辑图,自动生成算法是可以满足要求的.      

Outline of metallogeny of copper deposits of Cuba

Copper deposits of Cuba belong to the Jurassic lead-zinc-copper, Cretaceous copper and Paleogene copper complexes, which differ substantially in their formation conditions. The island is differentiated into a number of metallogenic zones formed at the initial and early stages of the Cuba folded province. Pinar del Rio, Pinos, Trinidad and Oriente zones, of an intrageanticlinal nature, are prospective with respect to copper vein deposits of the lead-Zinc-copper complex. The Sasa and Cauto zones - intrageosynclinal - are prospective with regard to copper-pyrite and copper vein deposits of the Cretaceous complex. A metallogenic study of the island is indicated in planning the exploration of its copper prospects. -- AGI Staff.     

多年冻土融沉性分类研究

在冻土工程中,冻土的融沉性评价是工程地质勘察的主要内容之一,融沉性分类是冻土地基基础设计施工的重要依据．根据345个冻土原状样品融沉压缩试验数据,提出了细砾、砂土、粉土、黏性土、泥炭化黏性土和泥炭质土等6类土的融沉系数一含水量或融沉系数一超塑含水量线性回归方程式,得到与各融沉性分级相应的界限含水量或界限超塑含水量．最后...