基于损伤塑性模型的地下洞室结构地震作用分析

基于损伤力学原理,通过引入损伤变量的方法,详细推导了ABAQUS中的损伤塑性本构模型,描述混凝土后继屈服阶段及其卸荷后的损伤机制和力学行为。此本构模型可考虑循环荷载条件下材料刚度退化和应变率的影响,适用于类似混凝土材料的脆性材料,如岩石等。以某水电站地下厂房洞室结构为例,进行二维动力时程非线性分析,考虑不同的刚度恢复权重因子和应变率相关性对结构损伤的影响。多个数值算例表明：洞室混凝土结构和其周围一定深度的围岩的损伤是一个渐进积累的过程;材料的率相关性效应并不是在地震动作用开始时就立即产生,而是在一段时间后,结构经历一定程度的损伤后才慢慢体现出来;此损伤塑性模型适合地震荷载等循环荷载作用下地下洞室结构的动力时程非线性分析,尤其是强震下采用此损伤塑性模型是合理和必要的。     

基于最小耗能原理的黏弹塑性损伤模型

在基于“应变等效假定”和“余应变能等效假定”两种不同损伤力学基本假定的黏塑性损伤本构方程基础上,将最小耗能原理应用于固体介质的黏弹塑性问题,给出了损伤材料的黏塑性应变流变速率、黏塑性损伤发展速率和黏弹-塑性材料积累应变强弱化速率的非线性演化发展方程,建立了两个不同的黏弹塑性损伤模型（A和B）,并进行了试验验证。作为算例,最后运用黏塑性损伤模型A,对地震荷载作用下建于岩石基础上的某混凝土重力坝的动力损伤问题进行了分析,结果表明该模型应用于混凝土重力坝的动力损伤分析是可行的。     

鱼简河水库坝基岩体模量研究

岩体变形模量的确定是岩石力学界研究的难点问题。贵州省鱼简河水利工程是国家重点工程,该工程大坝坝型为拱坝,坝基岩体变形模量的大小直接影响其坝体结构设计。本文在坝基岩体分组的基础上,利用Monte-Carlo模拟技术,采用RMR法与岩体损伤的分形模型对岩体变形模量进行估算,最后结合原位试验和工程地质类比给出岩体参数建议值,为坝体设计提供合理的计算参数。      

金沙江阿海水电站混凝土重力坝地震抗滑稳定性的时程法分析

通过动力有限元时程法分析了金沙江阿海混凝土重力坝在设计地震动(0.344g)作用下的动态响应,计算模型考虑了动水压力作用、大坝-岩基相互作用,用黏弹性吸波边界模拟了地震动能量向无限远域逸散的地基“辐射阻尼”效应。通过将坝体惯性力时程极大值与振型分解反应谱法的惯性力极值对比,从而验证了复杂动力数值模拟结果的正确性,为后续分析打下了基础。采用沿滑面的应力积分法得到了大坝沿建基面的瞬态抗滑安全系数时程。计算结果表明：阿海重力坝在设计地震动(0.344g)作用下的安全系数时程的极小值为1.146,大坝整体稳定性良好,可不采取提高整体抗滑     

考虑初始损伤和残余强度的岩石变形过程模拟

针对现有损伤模型的不足,提出同时反映岩石初始损伤和残余强度的岩石变形过程模拟方法.将岩石材料分为未损伤部分、损伤部分和微缺陷,未损伤部分和损伤部分共同承担岩石所受的应力.通过岩石材料几何条件及微观受力分析建立未损伤部分的应变分析.再考虑岩石的形变比能只与损伤部分材料相关,从而建立损伤部分的应变分析.然后探讨岩石损伤和耗散能的关系,提出了可以反映岩石初始损伤的损伤演化方程,基于各部分应变分析建立模拟岩石变形过程的损伤本构模型,同时给出模型各参数的确定方法.结果表明:损伤演化方程不仅表现了岩石的初始损伤特征,也可以完整地体现岩石变形破坏过程的5个阶段;与前人模型相比,提出的岩石变形模拟方法与试验曲线更加吻合,体现了岩石的变形全过程,特别是更加直观地反映了岩石的初始损伤和残余强度特征.综合体现了本文模型在模拟岩石变形全过程时的优势.     

深厚砂质覆盖层土坝的弹塑性地震反应分析

由于可液化砂质土应力-应变特性模拟的复杂性及数值计算的不稳定性,深厚砂质覆盖层土坝的弹塑性地震反应分析是土坝抗震研究中的一个尚未完全解决的课题。采用u-p完全耦合的饱和多孔介质有限元分析方法和砂土多重机构弹塑性模型,对遭受M6.7级地震的国外某深厚砂质覆盖层土坝进行弹塑性地震反应分析,研究了坝体和地基的动力反应特性及其超静孔隙水压力产生、扩散和消散的变化规律。结果表明：计算得到的坝体加速度和永久变形与实测值存在一定的差异,但基本上反映了坝体加速度与永久变形的实际分布情况,从而说明采用的本构模型和计算方法具有一定的精度;由于坝体和坝基的超静孔隙水压力较小,且坝体永久变形不大,可以不对坝体和坝基进行加固处理;坝趾附近浅层地基的超静孔隙水压力较大,有可能发生液化,因此,须采取相应的抗液化加固措施。     

Rock Mass Behavior Under Hydropower Embankment Dams: A Two-Dimensional Numerical Study

Sweden has more than 190 large hydropower dams, of which about 50 are pure embankment dams and over 100 are concrete/embankment dams. This paper presents results from conceptual analyses of the response of typical Swedish rock mass to the construction of a hydropower embankment dam and its first stages of operation. The aim is to identify locations and magnitudes of displacements that are occurring in the rock foundation and grout curtain after construction of the dam, the first filling of its water reservoir, and after one seasonal variation of the water table. Coupled hydro-mechanical analysis was conducted using the two-dimensional distinct element program UDEC. Series of the simulations have been performed and the results show that the first filling of the reservoir and variation of water table induce largest magnitudes of displacement, with the greatest values obtained from the two models with high differential horizontal stresses and smallest spacing of sub-vertical fractures. These results may help identifying the condition of the dam foundation and contribute to the development of proper maintenance measures, which guarantee the safety and functionality of the dam. Additionally, newly developed dams may use these results for the estimation of the possible response of the rock foundation to the construction.     

西南某电站河床坝基开挖卸荷条件下变形破坏的数值分析

西南某电站河床坝基开挖过程中,在河床底部出现大量表征岩体卸荷的变形破坏现象。开挖卸荷深度与卸荷强度成为大坝建设及其长期稳定性的关键问题。采用三维有限差分程序FLAC3D,分析开挖卸荷条件下坝基岩体变形破坏特征。结合现场卸荷裂隙现象调查资料,综合确定出所研究的河床坝基开挖卸荷影响深度及强度。为坝基综合治理提供必要的参考信息,也为今后此类中等地应力边坡开挖卸荷分析储备可以借鉴的资料。     

A stepwise damping‐solvent extraction method for large‐scale dynamic soil–structure interaction analysis in time domain

Time‐domain analysis of dynamic soil–structure interaction based on the substructure method plays an increasing role in practical applications as compared with the frequency‐domain analysis. Efficient and accurate modelling of the unbounded soil or rock medium has been a key issue in such an analysis. This paper presents a subregional stepwise damping‐solvent extraction formulation for solving large‐scale dynamic soil–structure problems in the time domain. Accuracy and efficiency of the formulation are evaluated in detail for a classical problem involving a rigid strip foundation embedded in a half‐space. A practical large‐scale soil–structure interaction problem, which represents a high concrete gravity dam subjected to seismic load, is then analysed using the proposed method. Various responses of the dam, including time histories of the crest displacement and acceleration and contours of the peak principal stresses within the dam body, are presented. Comparisons are also made between these results with those obtained using other models for the unbounded medium. Copyright © 2007 John Wiley & Sons, Ltd.     

The Dynamic Evaluation of Rock Slope Stability Considering the Effects of Microseismic Damage

A state-of-the-art microseismic monitoring system has been implemented at the left bank slope of the Jinping first stage hydropower station since June 2009. The main objectives are to ensure slope safety under continuous excavation at the left slope, and, very recently, the safety of the concrete arch dam. The safety of the excavated slope is investigated through the development of fast and accurate real-time event location techniques aimed at assessing the evolution and migration of the seismic activity, as well as through the development of prediction capabilities for rock slope instability. Myriads of seismic events at the slope have been recorded by the microseismic monitoring system. Regions of damaged rock mass have been identified and delineated on the basis of the tempo-spatial distribution analysis of microseismic activity during the periods of excavation and consolidation grouting. However, how to effectively utilize the abundant microseismic data in order to quantify the stability of the slope remains a challenge. In this paper, a rock mass damage evolutional model based on microseismic data is proposed, combined with a 3D finite element method (FEM) model for feedback analysis of the left bank slope stability. The model elements with microseismic damage are interrogated and the deteriorated mechanical parameters determined accordingly. The relationship between microseismic activities induced by rock mass damage during slope instability, strength degradation, and dynamic instability of the slope are explored, and the slope stability is quantitatively evaluated. The results indicate that a constitutive relation considering microseismic damage is concordant with the simulation results and the influence of rock mass damage can be allowed for its feedback analysis of 3D slope stability. In addition, the safety coefficient of the rock slope considering microseismic damage is reduced by a value of 0.11, in comparison to the virgin rock slope model. Our results demonstrate that microseismic activity induced by construction disturbance only slightly affects the stability of the slope. The proposed feedback analysis technique provides a novel method for dynamically assessing rock slope stability and can be used to assess the slope stability of other similar rock slopes.     

Seismic performance evaluation of dam-reservoir-foundation systems to near-fault ground motions

Ground motion records obtained in recent major strong earthquakes have provided evidence that ground motions recorded near the near-fault regions differ in many cases from those observed further away from the seismic source. As the forward directivity and fling effect characteristics of the near-fault ground motions, they have the potential to cause more considerable damage to structures during an earthquake. Therefore, understanding the influence of near-fault ground motions on the performance of structures is critical to mitigate damage and perform effective response. This paper presents results of a study aimed at evaluating the effects of near-fault and far-fault ground motions on seismic performance of concrete gravity dams including dam-reservoir-foundation interaction. Koyna gravity dam is selected as a numerical application. Four different near-fault ground motion records with an apparent velocity pulse are used in the analyses. The earthquake ground motions recorded at the same site from other events that the epicenter far away from the site are employed as the far-fault ground motions. The seismic performance evaluation method based on the demand-capacity ratio, the cumulative overstress duration and the spatial extent of overstressed regions is presented. The concrete damaged plasticity model including the strain hardening or softening behavior is employed in nonlinear analyses. Nonlinear seismic damage analyses of the selected concrete dam subjected to both near-fault and far-fault ground motions are performed. The results obtained from the analyses show the effects of near-fault ground motions on seismic performance of concrete gravity dams and demonstrate the importance of considering the near-fault ground excitations.     

考虑体积变化影响的岩石应变软硬化损伤本构模型及参数确定方法

针对现有岩石损伤模型的局限性与不足,将应力作用下的岩石抽象为空隙、损伤与未损伤材料三部分,以空隙率反映岩石体积或空隙的变化,以损伤变量或损伤因子反映岩石力学性态的改变程度,通过对各部分材料的微观受力及其与宏观应力和应变的关系进行分析,建立了可反映岩石变形过程中体积或空隙变化特性的新型岩石损伤模型;从岩石微元强度合理度量方法研究入手,引进统计损伤理论,建立了可考虑损伤阀值影响的岩石统计损伤演化模型,进而建立了能模拟应变软硬化全过程的岩石统计损伤本构模型,并提出了参数确定方法。该模型不仅能反映岩石应变软硬化特性,而且能反映岩石变形的初期压密和后期扩容的特性。通过实例分析与讨论,表明模型合理与可行。     

Continuum large cracking in a rate‐dependent plastic–damage model for cyclic‐loaded concrete structures

Formulation and algorithmic treatment of a rate‐dependent plastic–damage model modified to capture large tensile cracking in cyclic‐loaded concrete structures are presented in detail for a three‐dimensional implementation. The plastic–damage model proposed by Lee and Fenves in 1998 was founded based on isotropic damaged elasticity in combination with isotropic multi‐hardening plasticity to simulate cracking and crushing of concrete under cyclic or dynamic loadings. In order that the model can capture large crack opening displacements, which are inevitable in plain concrete structures, the excessive increase in plastic strain causing unrealistic results in cyclic behaviors is prevented when the tensile plastic–damage variable controlling the evolution of tensile damage is larger than a critical value. In such a condition, the crack opening/closing mechanism becomes similar to discrete cracking. The consistent tangent operator required to accelerate convergence rate is also formulated for the large cracking state including viscoplasticity. The validation and performance of the modified algorithm implemented in a special finite element program is exemplified through several single‐element tests as well as three structural applications. The last example examines the model in the seismic fracture analysis of Koyna dam as a benchmark problem and the resulting crack profile is compared with the available experiment. The numerical experimentations well demonstrate that the developed model whose modification is necessary to properly simulate the cyclic behavior of plain concrete subjected to large tensile strains is robust and reasonably accurate. Copyright © 2012 John Wiley & Sons, Ltd.     

地下结构地震破坏静-动力耦合模拟研究

地下结构地震破坏过程模拟是研究地下结构抗震性能的有效手段,通过人工边界条件设置将围岩土体的自重应力作为动力响应模拟的初始状态,实现了地下结构地震响应静-动力耦合数值模拟。基于有限元分析软件ABAQUS,将极限剪应变作为破坏标准,采用生死单元法模拟了地下结构的地震破坏过程,通过结构埋深考虑了初始地应力对地下结构抗震性能的影响。模拟结果表明,地下结构的抗震性能存在一个最不利埋深,即最容易地震破坏的埋深,当结构的埋深超过最不利埋深时,结构的极限抗震能力增强,若结构产生地震破坏,则结构埋深越深即初始地应力越高其破坏程度也越强。     

一种改进的光滑粒子动力学法在地震条件下岩质边坡中的验证与应用

开发一种改进损伤框架的粒子流算法,被称为核断裂的光滑粒子流法（kernel-broken smoothed particle hydrodynamics,KBSPH）,用于模拟地震条件下岩质边坡的裂纹扩展和变形破坏过程.在KBSPH中,提出一种改进的损伤框架,通过引入断裂标志来改进损伤粒子的核函数,使损伤粒子的虚拟应力键直接断裂,裂纹在断裂的应力键间生成,从而模拟岩石的裂纹扩展过程.在地震边界上采用了双层边界,将动力输入边界与黏滞边界分离.首先通过薄板振动实验验证KBSPH的动力特性.其次以单裂隙岩体单轴压缩试验验证KBSPH的断裂力学特性.最后模拟地震条件下多节理岩质边坡中裂纹扩展过程和动力响应.薄板振动实验验证了KBSPH的动力特性的准确性.单裂隙岩体单轴压缩试验,证明了KBSPH可以正确模拟预制裂隙尖端的翼型裂纹.通过对比以往数值模拟方法和现场案例,表明KBSPH正确揭示了加速度放大效应以及地震条件下岩质边坡的裂纹扩展过程.KBSPH避免了传统算法的网格畸变,损伤粒子应力分量重新分配的问题,降低了编程难度,提高了运行速率,可为SPH在地震条件下岩石力学中的应用和理解岩石断裂机理提...     

锦屏一级坝基岩体质量爆破开挖损伤评价

锦屏一级水电站大坝为双曲拱坝,设计坝高305m,拱坝建基岩体以Ⅱ、Ⅲ₁级岩体为主,开挖过程中岩体质量受爆破开挖影响发生程度不等的损伤。为定量评价岩体质量损伤等级,利用开挖爆破检测成果,采用开挖后岩体松弛深度、松弛岩体开挖前后单孔声波衰减率作为岩体质量损伤分级评价的指标,探索一种快速检测、评价岩体质量损伤等级的划分标准,完成坝基开挖后岩体质量损伤的分级分区评价,提出不同损伤等级处理的建议,并评价处理效果。     

不同温度条件下盐岩卤水浸泡后损伤自恢复特性

为了研究地下盐穴开挖扰动区（EDZ）损伤盐岩在盐穴综合利用过程中的自恢复特性,设计了损伤盐岩在卤水温度作用下的自恢复试验。结果表明：损伤盐岩经过卤水浸泡后在不同温度下烘干时力学性质可以得到恢复,在试验温度范围35～80 ℃内,其恢复效果随着温度的增加而增强。经过恢复处理的损伤盐岩试件二次加载时的声发射特征显示,其Felicity比低于0.1,Kaiser效应基本不成立,说明损伤盐岩在自恢复处理过程中其内部损伤在细观尺度上发生愈合。根据损伤盐岩试件自恢复处理后的细观形貌特征,总结出盐岩损伤恢复3种不同的细观机制：基于扩散作用的损伤微裂纹愈合,基于NaCl晶体重结晶作用的微裂隙充填以及破碎区NaCl晶体颗粒的愈合黏结。盐岩损伤恢复过程中,裂纹中充填的饱和卤水为裂纹中的物质传递提供了通道,同时提供了大量的结晶基础物质,为盐岩损伤愈合奠定了基础。温度的升高为盐岩损伤恢复提供了更高的能量,提升了损伤愈合速度,同时使氯化钠重结晶密度增加,从而提高了损伤恢复程度。     

An anisotropic damage–plasticity model for saturated quasi‐brittle materials

This article is devoted to numerical modeling of anisotropic damage and plasticity in saturated quasi‐brittle materials such as rocks and concrete. The damaged materials are represented by an isotropic poroelastic matrix containing a number of families of microcracks. Based on previous works, a discrete thermodynamic approach is proposed. Each family of microcracks exhibits frictional sliding along crack surfaces as well as crack propagation. The frictional sliding is described by a Coulomb–Mohr‐type plastic criterion by taking into account the effect of fluid pressure through a generalized effective stress concept. The damage evolution is entirely controlled by and coupled with the frictional sliding. The effective elastic properties as well as Biot's coefficients of cracked porous materials are determined as functions of induced damage. The inelastic deformation due to frictional sliding is also taken into account. The procedure for the identification of the model's parameters is presented. The proposed model is finally applied to study both mechanical and poromechanical responses of a typical porous brittle rock in drained and undrained compression tests as well as in interstitial pressure controlled tests. The main features of material behaviors are well reproduced by the model. Copyright © 2012 John Wiley & Sons, Ltd.     

岩石材料的三维各向异性损伤破坏模型与数值模拟

基于复合材料以及连续介质损伤理论,给出了岩石材料的各向异性损伤破坏模型。通过引入与岩石材料单轴加载行为相对应的特征模态构成的四阶对称损伤张量,描述了岩石材料的损伤演化过程,其中对不同主应变方向采用不同的损伤变量,而对同一主应变方向拉压时的损伤则采用不同的损伤变量来描述。在数值模拟岩石破坏过程的程序中,采用了张量分解的方法。将该模型编写用户材料子程序,并嵌入到大型有限元分析程序ABAQUS中,通过ABAQUS/EXPLICIT SOLVER的显式有限元算法求解。利用此程序对岩石材料的单轴压缩进行了数值模拟。     

基于RNN模型的坝体和岩基区间参数反演方法研究

针对混凝土大坝坝体和岩基参数的区间不确定性,构造具有区间分析功能的RNN（粗糙神经网络）模型,并运用该模型反演坝体和岩基区间参数值。应用区间有限元对结构进行正分析,根据区间参数反演的需要选取相应的区间学习样本,利用RNN模型对样本进行模式学习直至网络收敛,最后通过网络回想和反归一法得到坝体和岩基力学参数的区间值。研究结果表明,该方法可用于反演混凝土坝坝体和岩基区间力学参数,反演得到的区间参数值是合理的。此外,基于RNN模型的区间参数反演方法经过一定的拓展和改进,理论上可应用于反演其他类型的区间参数。