三维梁-颗粒模型在岩石材料破坏模拟中的应用

用三维梁-颗粒模型BPM3D(beam-particlemodelinthreedimensions)对岩石类非均质脆性材料的力学性质和破坏过程进行了数值模拟。梁-颗粒模型是在离散单元法基础上,结合有限单元法中的网格模型提出的用于模拟岩石类材料损伤破坏过程的数值模型。在模型中,材料在细观层次上被离散为颗粒单元集合体,相邻颗粒单元由有限单元法中的弹脆性梁单元联结。梁单元的力学性质均按韦伯(Weibull)分布随机赋值,以模拟岩石类材料力学参数的空间变异性。材料内部裂纹通过断开梁单元来模拟。通过自动生成的非均质材料模型对岩石类材料的破坏机理进行研究。岩石类非均质脆性材料在单轴压缩状态下破坏过程细观数值模拟结果显示,岩石材料宏观破坏是由于其内部细观裂纹产生、扩展、贯通的结果。通过数值模拟结果之间的对比分析,揭示出岩石试样宏观破坏模式随细观层次上韦伯分布参数的变化而不同。与实际矿柱破坏形态的对比分析表明了模型的适用性。根据数值模拟结果对岩石类非均质材料的破坏机理进行了探讨。     

矿物粒径对花岗岩单轴压缩特性影响的试验与模拟研究

高放废物深部地质处置目前受到世界各国的高度重视。花岗岩是我国高放废物地质处置工程的候选围岩,深入了解处置库花岗岩的强度及破坏特性对于处置系统的设计及性能评价具有十分重要的意义。作为矿物颗粒的集合体,花岗岩是一种由石英、长石和黑云母等矿物组成的非均质岩石,矿物粒径对其宏观力学特性影响明显。以我国高放地质处置库预选区阿拉善花岗岩为例,选取矿物粒径差异明显的似斑状花岗岩和中粒花岗岩两类岩石,采用单轴压缩试验与数值模拟相结合的方式研究了矿物粒径对岩石力学特性的影响。单轴压缩试验在MTS815岩石力学试验系统进行,数值模拟采用基于离散元的颗粒流程序PFC2D完成。数值模拟过程中,以试件表面图像为基础,采用数字图像处理技术获取岩石内部矿物组分的实际空间分布,从而建立了精确反映花岗岩内部矿物种类及其空间位置的数值模型。利用该模型对花岗岩的单轴压缩试验进行了数值模拟,并与试验结果对比,论证了模型的可靠性。试验及模拟结果表明,阿拉善花岗岩破坏形式为脆性张拉破坏,裂纹大多平行于轴压方向,数字图像数值分析方法可真实地反映材料细观结构。矿物粒径对材料力学特性的影响主要表现为:细粒、等粒结构的岩石强度高,粗粒、不等粒结构的岩石强度低。研究成果可为掌握矿物粒径对岩石强度及变形特性的影响提供依据。     

循环荷载下砂土液化特性颗粒流数值模拟

利用PFC2D常体积循环双轴试验条件,对砂土在不排水循环荷载作用下的液化特性进行了颗粒流数值模拟,数值模拟按等应力幅加荷方式进行。颗粒流数值模拟的优点在于得到试样液化宏观力学表现的同时,通过不同循环加荷时刻试样内细观组构参量（包括配位数、接触法向分布、粒间法向接触力、粒间切向接触力）的演化规律,分析砂土液化过程中细观组构变化与宏观力学响应之间的内在联系,从而可进一步探讨砂土液化的细观力学机制。数值模拟研究结果表明,砂土液化现象在宏观力学表现上反映为超静孔隙水压力的累积上升和平均有效主应力的不断减小,在细观组构上对应于配位数的累积损失和粒间接触力的不断减小。砂土液化细观机制分析表明,试样配位数的减少与循环加荷过程中组构各向异性滞后于应力各向异性有关。     

散粒体中锚杆间距的锚固效应研究

采用细观数值模拟方法研究散粒体的锚固效应,基于随机模拟技术生成三维多面体颗粒及其在空间中的分布,在随机散粒体不连续变形模型的基础上将砾石锚固试验进行数值实现,分析加锚散粒体材料的宏观与细观力学性能,研究加锚密度及其与颗粒粒径的关系对散粒体力学性质的影响,并探讨锚杆在散粒体材料中的作用机制。分别建立不同锚杆间距和不同颗粒粒径的数值试样,数值模拟结果表明：散粒体锚固数值试验能够较好地反映不同加锚散粒体结构的变形规律与锚固效应;散粒体材料的宏观特性与其细观组构的演化密切相关;锚杆加固散粒体的作用机制为加锚散粒体内形成压缩区,挤压加固作用提高了散粒体间的接触作用力,散粒体结构的整体性得到加强并能承受一定荷载;不考虑锚杆长度的情况下,当锚杆间距小于3倍的散粒体平均粒径时,锚杆能够有效地加固散粒体形成稳定结构。     

岩石细观结构对其变形强度影响的数值分析

利用数字图像处理技术表征岩石的非均匀性,并映射到有限元网格中,建立了能比较准确地反映岩石细观结构的数值模型。通过对花岗岩数字图像的处理,利用该模型进行花岗岩单轴压缩数值试验,研究花岗岩细观结构对其变形和强度的影响。研究表明,花岗岩细观结构对应力集中影响显著,在不同方向对花岗岩进行加载时,其抗压强度差异较大,表现出较强的各向异性行为,但对花岗岩弹性模量的影响很小。基于数字图像的数值分析方法为研究岩石细观力学性质提供了一种新的方法。     

石灰岩细观力学特性的颗粒流模拟

岩体地区地质灾害的发生和发展取决于岩石细观组分的运动学行为。研究岩石运动学行为时通常将岩石作为整体研究对象较多,而直接以细观组分为对象的研究较少。以石灰岩为例,根据室内试验获得的岩石力学性质指标,使用基于非连续介质理论的颗粒流方法,将材料离散成刚性颗粒组成的模型,把颗粒细观变化与宏观力学特性联系起来,建立了石灰岩的细观结构模型,获得了颗粒接触力、颗粒接触模量、接触连接强度和连接刚度比等细观力学参数。由于文中直接以细观成分为研究对象、反映了岩石和岩体组成的本质特点,所得结论不仅对含裂隙岩石本构关系研究具有广阔的应用前景,而且对岩体工程性质和地质灾害机制研究也具有重要的理论意义。     

基于Otsu图像处理的岩石细观模型重构及数值试验

基于大津法(Otsu)多阈值分割方法,利用数字图像处理技术获取了花岗岩细观结构的表征图像,并结合颗粒流程序,重构了反映岩石非均质结构特征的细观模型。通过单轴压缩试验,分析了花岗岩的失稳破坏过程,并研究了细观结构强度对宏观强度及破坏形态的影响。研究发现:花岗岩试件失稳破坏过程中存在4种裂纹形式,加载初期首先产生界面裂纹,其次是云母裂纹,长石和石英在应力接近极限强度时才破坏;当增大胶结界面或矿物强度时,花岗岩的抗压强度随之增加,但增加梯度逐渐减小;胶结界面或3种矿物细观强度的增大会使该结构处裂纹数量减小,而增加其他结构处的裂纹分布;相对于中等强度的长石和高强度的石英,低强度的胶结界面和云母矿物对花岗岩宏观强度和破坏形态的影响较大,是关键细观结构。将Otsu图像处理技术与颗粒流程序相结合,建立反映岩石真实结构的细观模型,为研究岩石非均质性对宏观力学特性的影响提供了一种有效手段。     

应用数字图像分析法评价红砂岩渐进损伤破坏特性

为了分析含水状态下红砂岩渐进破坏行为,在已有的岩石细观加载仪上配置数字显微观测系统,对试件表面的裂纹扩展过程进行多角度的实时全场监测和局部显微放大监测,并将可监视到的图像进行录像和计算机自动采集处理。利用该系统进行了红山窑水利枢纽工程地基红砂岩在不同含水率条件下的单轴压缩破裂过程试验,应用数字图像相关技术处理试样表面细观裂纹萌生、扩展图像灰度分布图,去除伪裂纹,探讨了红砂岩动态损伤状态与其宏观力学响应之间的关联。为了获得试样表面裂纹萌生、扩展的定量描述,引入表观裂损度与损伤张量2个物理量,定量评价含水状态下红砂岩渐进损伤破坏特性,进一步分析单轴压缩状态下岩石裂纹损伤演化规律。     

宏细观缺陷对岩体力学特性及边坡稳定影响研究

岩体是同时含有宏、细观缺陷的天然损伤地质体,如何更好地反映该两类缺陷对岩体力学特性的共同影响是亟待解决的一个重点和难点问题。采用FLAC3D中的弹脆性模型和Null模型分别描述完整岩石和上述两类缺陷的力学行为,用超细单元剖分计算模型的方法模拟岩体破坏,同时提出一个利用岩石孔隙率反映细观缺陷的新方法,研究细观缺陷即岩石孔隙率、裂隙倾角及长度等对裂隙岩体力学特性的影响。利用该方法研究了上述两类缺陷对岩质边坡临界滑面(CFS)及安全系数(FS)的影响。结果表明,宏观缺陷对岩体单轴压缩破坏模式、峰值强度和弹性模量起控制作用,同时控制岩质边坡的破坏模式和安全系数;尽管细观缺陷并不能改变宏观缺陷对岩体力学行为的控制作用,但其对这种控制作用将产生一定影响,宏细观缺陷以不同的作用机制对岩体力学特性产生影响。     

基于BP神经网络的土体细观力学参数反演分析

利用离散元方法对颗粒材料的细观力学特性研究, 目前确定数值计算模型的细观力学参数大多数通过反复调试获取, 效率低、可重复性差。本文采用开源的颗粒离散元程序LMGC开展了土体双轴压缩数值试验, 通过25组土体细观力学参数计算得到相应的宏观力学参数, 建立了BP人工神经网络反演系统。利用土体物理试验得到的土体宏观力学参数, 输入BP神经网络, 反演得到土体的细观力学参数。将所得细观力学特性参数输入所建立的土体数值计算模型, 得到土体破坏过程中的应力-应变关系曲线, 以及土体颗粒的力链图和旋转变形云图。所建立的土体数值试验模型能够较好地模拟土体变形破坏过程, 利用BP神经网络反演细观力学参数以及数值模型计算得到的土体宏观力学参数与物理试验吻合较好, 误差在10%左右, 土颗粒间力链云图以及旋转变形云图较好地揭示了土体变形破坏的机理。     

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

在土压缩模量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.     

地震次生地质灾害风险评估集对态势分析

5．12四川汶川8级强震引发了大量次生地质灾害,加剧了灾情并严重影响抗震救灾。地震次生地质灾害受多重因素制约,其孕育与致灾过程具有不确定性。利用集对分析原理与方法,对地震次生地质灾害进行同异反态势分析,拟定了风险评估的集对分析同一度、差异度、对立度等指标体系的构建原则与赋值标准,对承灾体系统的不确定性及其作用作了刻画与分析,建立了不同风险分区代表性的集对分析联系度表达式,为地震灾害研究提供了可资借鉴的新思路和技术方法。     

下扬子区三叠纪古地理演化

下扬子区在三叠纪期间接受了早-中三叠世海相碳酸盐沉积、中-晚三叠世海陆交互相和陆相湖沼沉积.沉积相带的空间分布和古地理格局的变迁明显地受到区域构造的控制.区内盆地发育与演化是在华南板块与扬子板块、扬子板块与华北板块相互作用的背景下进行的.早-中三叠世末期的印支运动是下扬子区构造-古地理格局改变的决定性因素,使下扬子海盆闭合,沉积类型由海相沉积变为陆相沉积.下扬子区三叠纪古地理的演化特征也为板块碰撞提供了沉积证据.     

多年冻土融沉性分类研究

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