颗粒物质在等比例应变加载下的分散性失稳模式

由地下水引起的静力液化可能是边坡失稳的隐含机制之一,松砂在不排水剪切条件下可能发生静力液化,密实的颗粒集合体在特定的应变路径下也会出现相似的现象,即试样整体发生急剧的失稳,应力状态尚处于峰值强度线以内。该种失稳模式称为分散性失稳,是为了强调失稳模式中没有出现应变局部化或者剪切带。采用连续-离散耦合分析方法,研究由不规则形状颗粒组成的密实集合体在等比例应变加载路径下的力学特性。根据Hill的材料失稳理论,当试样的应力增量 和应变增量 对应的2阶功 为负时,试样即发生不可逆的整体失稳破坏。以根据不同等比例应变路径得到 曲线为界,在 平面内将试样的应力状态分为剪缩区、剪胀-稳定区和剪胀-非稳定区,连接不同围压下试样发生分散性失稳时的应力状态形成失稳线发现,峰值强度线高于临界状态线,临界状态线高于失稳线。     

基于复杂网络的岩土颗粒材料分散性失稳先兆研究

岩土颗粒材料在达到Mohr-Coulomb塑性极限之前可能已经发生了没有明显局部损伤的全域准均匀破坏,即分散性失稳。近年来,网络科学工具在颗粒材料环境中的应用为研究其失稳破坏现象提供了有趣新颖的洞见。采用离散单元法对不同初始密实度的颗粒体系进行等体积应变加载路径的数值试验,发现初始状态越松散的颗粒体系越有可能发生分散性失稳;采用网络科学的理论和方法分析颗粒接触网络的拓扑结构特性和演化规律,发现颗粒材料分散性失稳的结构根源是颗粒接触网络完全崩塌。将颗粒体系划分为强、弱接触系统以构建强接触子网络、弱接触子网络和强−弱接触子网络,研究颗粒接触网络及其3个子网络的特征量演化,结果表明：当颗粒材料开始处于不稳定状态时,易受扰动的弱接触系统的部分接触先失效,导致强接触系统失去其稳定支撑作用而降低承载能力。随着加载的继续,系统通过非局部的自组织过程导致接触结构的全面崩塌,最终触发整体分散性失稳。因此,先于整体失稳的部分弱接触的失效,可视为颗粒材料分散性失稳的关键性临界征兆。     

空心圆柱扭转仪定向剪切试验离散元模拟

为分析砂土在复杂应力条件下的剪切力学特性,采用商业离散元软件PFC3D对单粒组中密砂的空心扭剪试验进行了仿真模拟,分析了数值试样的应力-应变关系,研究了不同剪切方向下离散介质的强度、体积应变特性以及中主应力比对它们的影响,再现了力链在加载过程中的演化,并对剪切带的倾角做了深入分析。同时,从细观上看,以颗粒接触数和纯转动率变量为中心,观察了试样内部颗粒的运动状态,对比了不同剪切方向下剪切带内外颗粒接触数与纯转动位移的变化。最后,将数值试验结果与已有的室内试验结果进行了对比。此研究实现了复杂应力条件下空心扭剪试样的三维离散元模拟,加深了对空心扭剪试验过程和结果的理解和解释。     

含软弱夹层土样变形破坏过程细观数值模拟及分析

基于颗粒流理论,引入接触连接模型和滑动模型,建立了含软弱夹层试样的颗粒流模型。通过颗粒流程序(PFC)数值模型试验,对含软弱夹层试样的强度和破坏发展进行了数值模拟,分别对比了不同围压以及不同夹层参数条件下的应力-应变关系曲线,通过位移矢量场分析了破坏发展趋势。模拟结果表明,一般围压条件下试样沿软弱夹层滑动破坏,但在某特定围压下,软弱夹层的存在并不起主要控制作用;试样应力-应变关系曲线峰值随软弱夹层颗粒的摩擦系数和法向接触刚度的减小而下降,当颗粒法向接触刚度趋近0时,试样加载初始阶段呈现塑性流变特征;夹层颗粒的切向接触刚度只有低于某一特定值时才会使试样应力-应变关系曲线峰值降低。通过分析,得到了土体颗粒细观参数和宏观力学行为的内在联系,并对土体软弱夹层力学性质和渐进破坏过程有了更进一步的认识。     

基于PFC3D的黄土三轴试验细观参数敏感性分析

从细观角度出发,研究各细观参数对黄土三轴试验宏观力学响应的影响,基于三维离散元PFC3D软件及理论,结合室内三轴试验结果确定数值模拟基本细观参数取值范围,在50 kPa、100 kPa、150 kPa、200 kPa围压下对三轴试样进行颗粒流模拟,通过改变各细观参数值对黄土数值试样宏观力学行为影响程度及正交试验设计的方式,开展黄土三轴数值模拟过程中摩擦系数、孔隙率、颗粒刚度比k_n/k_s、颗粒粒径分布的敏感性分析,建立土体细观参数与宏观力学的关系式。结果表明:不同围压下数值试样峰值强度和剩余强度均与摩擦系数正相关,与颗粒刚度比、孔隙率负相关,而颗粒粒径分布的变化对其影响均小于2%。低围压下(50 kPa),孔隙率对数值试样峰值抗剪强度和剩余强度影响最大,颗粒刚度比次之,高围压下(200 kPa),摩擦系数对数值试样强度的影响比颗粒刚度比更为显著;对数值试样的初始线弹性模量和应变软化特性影响最大的细观参数为孔隙率,刚度比次之,颗粒粒径分布最小,200 kPa围压下正交试验分析结果与数值分析结果一致。通过改变细观参数取值范围的方式系统地分析了其对黄土三轴试验宏观力学性能的影响,为今后离散元软件用于黄土室内试验研究提供了参考。     

基于抗转模型的颗粒材料宏−细观关系研究

接触模型的宏?细观参数标定是成功使用离散元方法的关键。在离散元的接触模型中线性接触模型与抗转线性接触模型均可用于模拟砂性土的力学行为,其中抗转线性接触模型在模拟密砂的剪胀性方面具备优势。采用抗转线性接触模型对室内密实砂土三轴试验进行了离散元模拟,验证了抗转线性接触模型的可靠性;进而系统分析了颗粒间摩擦系数、刚度比和抗转动系数等细观参数与砂土峰值内摩擦角、残余内摩擦角、峰值剪胀角等宏观参数的相关关系并进行了验证;揭示了偏应力作用下,细观参数对密实砂土试样内部剪切带宽度与倾角变化的影响规律,提出了考虑剪胀角的剪切带倾角经验公式。通过研究建立了抗转线性接触模型宏?细观参数的量化关系并给出了标定参数的具体流程图,提出了快速标定宏观参数的方法并应用实例进行了验证,为采用抗转线性接触模型精准模拟密实砂土的力学特性提供依据。     

抗转动对颗粒材料组构特性的影响研究

颗粒材料大多由不规则形状的颗粒组成,如砂土、谷物等,抵抗转动是不规则形状颗粒的固有特性。已有研究表明,颗粒抗转动特性对其宏观力学特性有显著影响。因此,在颗粒材料的细观数值模拟中或采用非圆颗粒,或在圆颗粒离散元模拟中采用考虑抗转动的接触模型。采用有限元-离散元耦合方法（FDEM）和离散元方法（DEM）分别对椭球形状颗粒和具有抗转动能力的圆球颗粒进行三轴剪切数值模拟,指出了采用抗转动接触模型考虑颗粒形状影响的局限性,并基于颗粒的局部排布结构揭示了形状影响的细观来源。峰值内摩擦角和剪胀均随着转动摩擦系数和形状偏离圆球程度而单调增加,但颗粒形状对它们的影响呈现出明显的收敛趋势。细观组构分析也表明,虽然颗粒形状和转动摩擦都能显著增强组构各向异性,但是组构各向异性的演化模式有明显的区别。造成以上结果的差异在于其抵抗转动的影响机制不同。转动摩擦是通过限制颗粒转动,增强了颗粒间的稳定承载能力,而非圆颗粒是通过咬合作用形成稳定的局部排列结构。由于椭球颗粒腹部比端部能够传递更大的接触力,颗粒受剪切后发生转动,颗粒长轴倾向于正交大主应力方向,呈现交错排列,颗粒间相互锁定。     

砂岩常规三轴的颗粒流数值模拟

以砂岩的室内试验数据为参照,使用颗粒流数值程序PFC3D,建立刚性颗粒组成的模型.通过设置颗粒接触模量、颗粒刚度比和接触粘结强度等细观力学参数,得到有合适宏观力学反应的模型.通过数值模拟获得砂岩微裂隙发展的过程,得到不同围压下砂岩的应力-应变曲线,反应出砂岩试样的弹性模量随围压增大而增大的性质.结果表明,PFC3D程序能够较好地模拟荷载作用下岩石的细观力学特性和宏观力学反应.     

基于微观力学的胶结岩土材料破损规律离散元模拟

采用离散元法（DEM）研究胶结岩土材料在不同加载条件下的结构破损规律。首先,基于微观力学理论,考虑胶结岩土材料颗粒间胶结特性,给出表征结构性损伤的破损参数式。该式具有微观物理意义,但不能直接用于建立宏观本构模型。其次,采用二维离散元源程序NS2D模拟等向压缩、等应力比压缩以及双轴压缩试验,分析破损参数在不同加载条件下随宏观力学变量（体积应变和剪应变）的演变规律。最后,根据模拟结果提出破损参数数学表达式,其为大主应变的函数。研究结果表明：胶结强度、应力比以及围压均在一定程度上影响了数值试样的结构破损规律。在等向压缩和等应力比压缩试验中,容易用函数式描述数值试样破损参数随体积应变或偏应变的演变规律;而在双轴压缩条件下,由于数值试样有剪胀特性,破损参数随体积应变的演变规律则不易描述。建议的破损参数数学表达式能够较好地描述数值试样在不同加载条件下结构破损规律。     

考虑转动阻抗的粗粒土离散元模拟

转动阻抗被定义为作用颗粒接触上的一对对称力偶,用来抵抗颗粒之间的相互转动。将转动阻抗引入到离散元模拟中是对传统离散单元法的重要改进。开发出考虑颗粒转动阻抗的接触模型,并将其嵌入到PFC2D中,利用该模型进行粗粒土的双轴剪切数值模型试验,研究剪切过程中转动阻抗对粗粒土的宏细观力学性质的影响。结果显示,在宏观方面,颗粒转动阻抗对粗粒土的宏观力学行为（应力-应变及体应变-轴应变行为）有重要的影响,随着转动阻抗的增加,粗粒土的剪切强度和最大摩擦角随之增加,这与已有的研究结果一致,证明所建模型是可靠的;在微观方面,考察转动阻抗对粗粒土内部微观结构的影响发现,随着转动阻抗的增加,粗粒土的内部的接触数目减少,而粗粒土的剪切强度增加,表明转动阻抗能够提高粗粒土力链网络的稳定性,同时发现随着转动阻抗的增加,粗粒土的各向异性增加主要是强力链各向异性的增加,说明转动阻抗增强了强力链的传递力的能力以及抵抗力链屈曲破坏能力。数值模拟结果表明,增加颗粒转动阻抗,粗粒土出现组构与轴应变非共轴的现象。     

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.     

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

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

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.     

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

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

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

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

黄土路基沉降曲线特征分析

在分析郑西客运专线试验段路基沉降实测资料的基础上,提出沉降变形常规预测方法受观测时间及单次观测结果的影响比较大,不利于用来预测最终沉降。通过分析现场沉降变形与时间关系曲线,提出路基沉降的发展主要受固结的影响,其发展是一个从沉降增加到趋缓的过程,可以明显分为固结变形及流变变形,并用非饱和土等效固结理论加以解释,提出了一种利用s- 曲线特征预测最终沉降量的方法,经检验对比发现其预测结果受观测时间及单次观测结果的影响比较小,优于其它常规预测方法     

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.