Hybrid finite-discrete element modelling of asperity degradation and gouge grinding during direct shearing of rough rock joints

A hybrid finite-discrete element method was implemented to study the fracture process of rough rock joints under direct shearing. The hybrid method reproduced the joint shear resistance evolution process from asperity sliding to degradation and from gouge formation to grinding. It is found that, in the direct shear test of rough rock joints under constant normal displacement loading conditions, higher shearing rate promotes the asperity degradation but constraints the volume dilation, which then results in higher peak shear resistance, more gouge formation and grinding, and smoother new joint surfaces. Moreover, it is found that the joint roughness affects the joint shear resistance evolution through influencing the joint fracture micro mechanism. The asperity degradation and gouge grinding are the main failure micro-mechanism in shearing rougher rock joints with deeper asperities while the asperity sliding is the main failure micro-mechanism in shearing smoother rock joints with shallower asperities. It is concluded that the hybrid finite-discrete element method is a valuable numerical tool better than traditional finite element method and discrete element method for modelling the joint sliding, asperity degradation, gouge formation, and gouge grinding occurred in the direct shear tests of rough rock joints.

Influence of location of large-scale asperity on shear strength of concrete-rock interface under eccentric load

The location and geometry of large-scale asperity present at the foundation of concrete gravity dams and buttress dams affect the shear resistance of the concrete-rock interface.However,the parameters describing the frictional resistance of the interface usually do not account for these asperities.This could result in an underestimate of the peak shear stre ngth,which leads to significantly conservative design for new dams or unnecessary stability enhancing measures for existing ones.The aim of this work was to investigate the effect of the location of first-order asperity on the peak shear strength of a concrete-rock interface under eccentric load and the model discrepancy associated with the commonly used rigid body methods for calculating the factor of safety(FS)against sliding.For this,a series of direct and eccentric shear tests under constant normal load(CNL)was carried out on concrete-rock samples.The peak shear strengths measured in the tests were compared in terms of asperity location and with the predicted values from analytical rigid body methods.The results showed that the large-scale asperity under eccentric load significantly affected the peak shear strength.Furthermore,unlike the conventional assumption of sliding or shear failure of an asperity in direct shear,under the effect of eccentric shear load,a tensile failure in the rock or in the concrete could occur,resulting in a lower shear strength compared with that of direct shear tests.These results could have important implications for assessment of the FS against sliding failure in the concrete-rock interface.

A Comparative Study on Asperity Peak Modeling Methods

The peak identification scheme based method(three-point definition)and the spectral moments based method(spectral moment approach)are both widely used for asperity peak modeling in tribology.To discover the differences between the two methods,a great number of rough surface profile samples with various statistical distributions are first randomly generated using FFT.Then the distribution parameters of asperity peaks are calculated for the generated samples with both methods.The obtained results are compared and verified by experiment.The variation rules of the differences between the two methods with statistical characteristics of rough surfaces are investigated.To explain for the discovered differences,the assumptions by spectral moment approach that the joint distribution of surface height,slope and curvature is normal and that the height distribution of asperities is Gaussian,are examined.The results show that it is unreasonable to assume a joint normal distribution without inspecting the correlation pattern of[z],[z′]and[z′′],and that the height distribution of asperities is not exactly Gaussian before correlation length of rough surface increases to a certain extent,20 for instance.

Influence of micro asperity contact and radial clearance on the tribological properties of crankpin bearings

A new hybrid numerical method that couples the dynamic slider-crank mechanism(SCM)and crankpin bearing(CB)lubrication models is proposed to analyze the effect of micro asperity contact on the tribological properties of a CB.In the hybrid model,the dynamic equations of the SCM are established based on the Newton method,while the lubrication equations of the CB are established on the basis of the Reynolds equation.Experimental data of the engine are also used in simulation analyses to enhance the reliability of the results.The load-bearing capacity(LBC)and friction force of the CB are selected as objective functions.Results show that the LBC has a negligible effect on the tribological properties of the CB,but the friction force greatly affects the resistance of the bearing under different radial clearances and surface roughness values.In particular,the maximum friction force in the asperity contact region accounts for 40.5%of the maximum total friction force at a radial clearance of 5μm and 77.7%of the maximum total friction of the CB with a surface roughness of 10μm.

Dual threshold search method for asperity boundary determination based on geodetic and seismic catalog data

As an important model for explaining the seismic rupture mode,the asperity model plays an important role in studying the stress accumulation of faults and the location of earthquake initiation.Taking Qilian-Haiyuan fault as an example,this paper combines geodetic method and b-value method to propose a multi-source observation data fusion detection method that accurately determines the asperity boundary named dual threshold search method.The method is based on the criterion that the b-value asperity boundary should be most consistent with the slip deficit rate asperity boundary.Then the optimal threshold combination of slip deficit rate and b-value is obtained through threshold search,which can be used to determine the boundary of the asperity.Based on this method,the study finds that there are four potential asperities on the Qilian-Haiyuan fault:two asperities(A1 and A2)are on the Tuolaishan segment and the other two asperities(B and C)are on Lenglongling segment and Jinqianghe segment,respectively.Among them,the lengths of asperities A1 and A2 on Tuolaishan segment are 17.0 km and 64.8 km,respectively.And the lower boundaries are 5.5 km and 15.5 km,respectively;The length of asperity B on Lenglongling segment is 70.7 km,and the lower boundary is 10.2 km.The length of asperity C on Jinqianghe segment is 42.3 km,and the lower boundary is 8.3 km.

Elastohydrodynamic Lubrication Interface Stiffness and Damping Considering Asperity Lateral Contact

Elastohydrodynamic lubrication(EHL)point contact occurs between two rough surfaces at the mesoscopic level,while the interaction of rough surfaces involves contact between asperities at the microscale level.In most cases,the contact between asperities within an interface takes the form of lateral contact rather than peak contact.Regions devoid of contact asperities are filled with lubricating oil.However,conventional models often oversimplify lateral contact forms as interactions between asperities and a smooth,rigid plane.These simplifications fail to accurately represent the true contact conditions and can lead to inaccuracies in the analysis of EHL’s contact performance.To address this issue,we have developed a novel EHL interface model comprising two rough surfaces.This model allows us to explore the influence of asperity height,contact angle,and contact azimuth angle on EHL interface performance.

Coupled evolution of piston asperity and cylinder bore contour of piston/cylinder pair in axial piston pump

The wear condition of the piston/cylinder pair is crucial to the performance and reliability of the axial piston pump.The hard piston surface,the soft cylinder bore surface,and the interface oil film affects each other during the wear process.Specifically,in the mixed lubrication region,the geometry of the hard piston surface asperity directly affects the wear of soft cylinder bore surface,while the asperities may deform or even degrade when penetrating and sliding against the cylinder bore.So far,there is no suitable method to simulate their coupled evolution.This paper proposed a wear process simulation model considering the real-time interaction between the elasto-plastic deformation of the piston surface asperity,the wear contour of the cylinder bore,and the lubrication condition of the interface.An offline library of the elasto-plastic constitutive behavior of the asperity based on the finite element method(FEM)is established as a part of the simulation model to precisely analyze the deformation and degradation of the asperity and quickly invoke them in the numerical wear process simulation.The simulation and experimental results show that the piston asperity and the cylinder bore contour converge to a steady state after running-in for about 0.5 h.The distribution of the simulated asperity degradation and wear depth is also verified by the experiment.

Numerical parametric study on the influence of location and inclination of large-scale asperities on the shear strength of concreterock interfaces of small buttress dams

When assessing the sliding stability of a concrete dam,the influence of large-scale asperities in the sliding plane is often ignored due to limitations of the analytical rigid body assessment methods provided by current dam assessment guidelines.However,these asperities can potentially improve the load capacity of a concrete dam in terms of sliding stability.Although their influence in a sliding plane has been thoroughly studied for direct shear,their influence under eccentric loading,as in the case of dams,is unknown.This paper presents the results of a parametric study that used finite element analysis(FEA)to investigate the influence of large-scale asperities on the load capacity of small buttress dams.By varying the inclination and location of an asperity located in the concrete-rock interface along with the strength of the rock foundation material,transitions between different failure modes and correlations between the load capacity and the varied parameters were observed.The results indicated that the inclination of the asperity had a significant impact on the failure mode.When the inclinationwas 30and greater,interlocking occurred between the dam and foundation and the governing failure modes were either rupture of the dam body or asperity.When the asperity inclination was significant enough to provide interlocking,the load capacity of the dam was impacted by the strength of the rock in the foundation through influencing the load capacity of the asperity.The location of the asperity along the concrete-rock interface did not affect the failure mode,except for when the asperity was located at the toe of the dam,but had an influence on the load capacity when the failure occurred by rupture of the buttress or by sliding.By accounting for a single large-scale asperity in the concrete-rock interface of the analysed dam,a horizontal load capacity increase of 30%e160%was obtained,depending on the inclination and location of the asperity and the strength of the foundation material.

Study of hydro-mechanical behaviours of rough rock fracture with shear dilatancy and asperities using shear-flow model

The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures,introducing complexities to fracture modelling.The present study focuses on the hydro-mechanical behaviours of rough rock fractures during shear-seepage processes to reveal how dilatancy and fracture asperities affect these phenomena.To achieve this,an improved shear-flow model(SFM)is proposed with the incorporation of dilatancy effect and asperities.In particular,shear dilatancy is accounted for in both the elastic and plastic stages,in contrast to some existing models that only consider it in the elastic stage.Depending on the computation approaches for the peak dilatancy angle,three different versions of the SFM are derived based on Mohr-Coulomb,joint roughness coefficient-joint compressive strength(JRC-JCS),and Grasselli’s theories.Notably,this is a new attempt that utilizes Grasselli’s model in shearseepage analysis.An advanced parameter optimization method is introduced to accurately determine model parameters,addressing the issue of local optima inherent in some conventional methods.Then,model performance is evaluated against existing experimental results.The findings demonstrate that the SFM effectively reproduces the shear-seepage characteristics of rock fracture across a wide range of stress levels.Further sensitivity analysis reveals how dilatancy and asperity affect hydraulic properties.The relation between hydro-mechanical properties(dilatancy displacement and hydraulic conductivity)and asperity parameters is analysed.Several profound understandings of the shear-seepage process are obtained by exploring the phenomenon under various conditions.

Forecasting the Epicenter Area of a Future Strong Earthquake

It is clear from the literature data that the geometric characteristics of the seismogenic source are directly related to the nucleation of their strong earthquakes. Our study starts from the observation of the proximity relationship between the epicenter of a strong earthquake and the surface trace of the seismogenic source. The proposed model is based on the mathematical analysis of seismicity falling within an area, through a polynomial function to determine a curve. The experimental results of our model confirm that 97% of the epicenters of strong earthquakes are located near the concavities or at the inflection points of the polynomial curve. Only 3% of strong earthquakes are located at the edges of the analyzed areas, where the number of data is not significant. The proposed model makes it possible to locate, with good approximation, the areas most likely to be the site of future strong earthquakes.

Multiscale study of the dynamic friction coefficient due to asperity plowing

A macroscopically nominal flat surface is rough at the nanoscale level and consists of nanoasperities.Therefore,the frictional properties of the macroscale-level rough surface are determined by the mechanical behaviors of nanoasperity contact pairs under shear.In this work,we first used molecular dynamics simulations to study the non-adhesive shear between single contact pairs.Subsequently,to estimate the friction coefficient of rough surfaces,we implemented the frictional behavior of a single contact pair into a Greenwood-Williamson-type statistical model.By employing the present multiscale approach,we used the size,rate,and orientation effects,which originated from nanoscale dislocation plasticity,to determine the dependence of the macroscale friction coefficient on system parameters,such as the surface roughness,separation,loading velocity,and direction.Our model predicts an unconventional dependence of the friction coefficient on the normal contact load,which has been observed in nanoscale frictional tests.Therefore,this model represents one step toward understanding some of the relevant macroscopic phenomena of surface friction at the nanoscale level.

A Model of Tangential Contact Damping Considering Asperity Interaction and Lateral Contact

A precise tangential contact damping model is proposed,which includes the lateral contact of the upper-lower asperities and the interaction of adjacent asperities.The effects of the normal static preload,frequency,and amplitude of tangential displacement on the tangential contact damping were analyzed by simulation,respectively.Furthermore,the results of simulation are verified by experiment.The tangential contact damping of considering the interaction and lateral contact of asperity is very close to the experimental results.

Shear behaviours and roughness degeneration based on a quantified rock joint surface description

The asperity wear of rock joints significantly affects their shear behaviour.This study discusses the wear damage of the asperities on the joint surface,highlighting the roughness degradation characteristics during the shear process.The direct shear experiment of artificial specimens containing rock joints was conducted under different normal stresses based on three-dimensional scanning technology.These experimental results showed the contribution of joint wear to roughness degeneration,such as the height,zone,and volume of asperity degeneration.The wear coefficient of the rock joint was obtained based on the volume wear of asperities in the laboratory experiment.The functional relationship between the friction coefficient and wear coefficient is subsequently determined.To quantitatively analyse the wear damage of a joint surface,a calculation method for determining the wear depth of the rock joint after shearing was proposed based on wear theory.The relationship between the ultimate dilation and wear depth was analysed.A coefficient m,which can describe the damage degree of the joint surface,and a prediction method of joint surface roughness after shearing are established.Good agreement between analytical predictions and measured values demonstrates the capability of the developed model.Lastly,the sensitivity factors on the wear depth are explored.

经验格林函数方法模拟强地面运动的研究进展

经验格林函数方法(简称EGFM)是当前模拟强地面运动的一种主流方法。主要介绍EGFM方法的发展历程和特点,重点解释了EGFM核心理论的沿革;详细介绍了加入凹凸体模型的改进的经验格林函数方法,并指出一种划分凹凸体区域新方法的可能性。

广州市主要断裂带特征震源模型研究

在广州市目标区内的主要活动断层危险性评价的基础上,综合活断层探测的研究成果和资料,对于具有潜在发震可能的主要断层,基于凹凸体震源模型建模理论进行了特征计算模型的建模。根据断层地震危险性分析确定的断层发震震级和几何参数,进行了断层的宏观参数和微观参数的设定。为通过复核预测方法计算与合成近断层强地震动场和城市危害性评价提供了科学的依据。

沈阳市F6断层的强地面运动分布模拟

目的研究沈阳市F6断层在设定地震6.0级时的强地面运动分布,为避免城市活动断层的地震破坏和为重要建筑合理的抗震设防提供依据.方法通过断裂所处应力环境、形变特征及活动性、结构、规模、地震活动、深部构造等综合分析,确定F6断层的地震危险性.以深、浅层人工地震、钻探、多道直流电法、探地雷达和地脉动数据为基础,建立了沈阳市接近真实的三维地下速度结构模型,强地面运动计算采用了三维有限差分法和统计学格林函数法的合成方法.结果 F6断层为中更新世活动断裂,最大潜在地震震级为6.0级;F6断层的震源模型为具有两个Asperity(凹凸体)的矩形断层:一个最大18km2的Asperity和一个7km2的Asperity.在设定地震作用下,F6断层强地震动PGA的影响范围集中于F6断层两侧,最大值达到285gal;PGV、PGD的分布形态与PGA相似.结论 F6断层发生6.0级时在沈阳市能够形成带状的强地面运动分布,其幅值甚至超出大震设防水准.

大震近场地震动数值模拟的不确定性研究

地震动数值模拟是研究地震发生、发展特征和规律性的基础研究,涉及震源、传播路径、场地三方面的一系列基本问题。研究结果可以直接用于强地震动基本参数的标定和地震动的工程预测以及地震危险性分析,是防震减灾中一项很重要的基础工作。当前考虑断层模型中的几个重要震源参数模拟未来地震动成为热点。例如,不管是随机有限断层方法,还是经验格林函数方法,都特别重视震源参数的选取,但当前对于强地面运动的数值模拟大多是对已发生地震的反演,反演所需的震源参数多是根据地震发生后的观测资料、经验关系推断或者人为的假定得到的。地震动特征预测的好坏,关键也在于震源参数的选取是否可靠,针对工程应用,对未来近断层地面运动的预测应当考虑参数选取的不确定性。强震动数值模拟方法已经比较成熟,比如随机有限断层法、经验格林函数法、混合方法等。但这些模拟地震动的方法多是确定性方法,对模拟过程中各个环节的不确定性因素考虑较少,尤其是断层震源参数的选取都是精确的数值,模拟结果自然不能综合表征未来地震发生时的地震动特征。另外一种地震动特性表征形式是利用地震动预测方程来表现参数的衰减关系,但是预测方程也会受观测数据、模型选取、随机因素的影响,不确定性无法避免。在地震动预测中每个步骤都带有很大的不确定性,主要包括认知的不确定性和随机的不确定性。震源模型的建立依赖各种专家的观点和意见,这种观点和意见往往相互差别很大,有很大的人为不确定性。地震资料的完整性、可靠性,以及理论模型等的建立也会造成很大的不确定性。因此,本论文的研究内容重点在预测未来大震近场地震动时震源环节的不确定性因素,系统地研究影响大震近场地震动数值模拟不确定性的影响因素,以期建立描述影响大震近场地震动数值模拟不确定性影响因素的模型,并将考虑不确定性因素的震源模型应用在破坏性大地震的地震动特征预测中去。针对地震发生时震源参数具有的随机不确定性,通过统计学方法得到震源参数间的经验关系。同时由于地球内部的不可入性以及观测技术手段的限制,地震震源模型还存在很多的认知不确定性因素,我们采用逻辑树方法处理认知不确定性因素,引入强震生成区的概念,重点分析Asperity的认知不确定性因素并在此基础上建立Asperity震源模型,取多种方案结果的最优值作为地震动的预测结果。论文具体的研究内容有以下5点:(1)在Hiroe Miyake工作的基础上,用经验格林函数法模拟了1997年日本九州鹿儿岛县的MJMA6.5地震,验证了经验格林函数法模拟强地面运动的有效性。数据使用的是日本K-NET强震台站数据。选择此次地震以及用日本K-NET数据的原因是日本的地震记录比较丰富,台站收集到的地震记录质量相对较好。同时分析了地震动模拟值中在工程领域使用较多的参数,表明模拟结果较好地反映了真实的地震动记录。证明经验格林函数预测未来强地震动是切实可行的。(2)本论文在大量地震记录及文献调研的基础上,运用统计学方法,通过统计回归分析得到地震密集区域震源参数间的经验关系,重点统计了包含龙门山断裂带在内的中国西南地区震源参数间的经验关系。为了得到更加适应于中国大陆区域的局部震源参数的经验关系,我们将地震样本局限在一定的区域内,缩小统计区域的范围,同时得到多个Asperity等局部震源参数的经验关系。(3)研究和分析各不确定性因素对地震动结果的影响程度对提高地震动预测的可靠性和准确性意义重大。本论文以2016年发生的日本熊本7.3级地震为例,通过参数敏感性分析技术筛选了对拟合结果有决定性影响的关键模型参数。未来地震动预测时会考虑在这几个不确定性因素中加入更高的权重,以提高地震动预测的可信度和准确性。(4)基于汶川地震资料,研究分析Asperity的认知不确定性对地震动预测的影响,运用逻辑树方法建立了考虑不确定性因素的Asperity震源模型。采用Hiroe Miyake提出的强震生成区的概念,将Asperity区域等同于强震生成区域,主要从Asperity的面积、数量、应力降、上升时间等方面进行了认知不确定性的分析研究,得到了针对强震的可能性较大的Asperity的面积和数量的最佳搭配形式。本文的创新点将概率地震危险性分析中不确定性因素的处理方法应用到未来大震近场地震动数值模拟中去。随机不确定性因素用统计学方法处理,认知不确定性用逻辑树方法处理。模拟结果能更好地表征未来大震近场的地震动特征,地震动参数的模拟结果在未来抗震设防及灾害预防中会有更好、更高的参考价值。考虑不确定性因素预测未来破坏性大地震,从而定量地评价发生地震时近场地面运动的综合特征,为防震减灾提供科学的依据,为工程设计决策提供参考依据。

镉胁迫对续断菊Sonchus asper L.Hill.根系分泌物的影响

通过盆栽试验,研究了不同Cd质量分数(0、50、100、200 mg.kg-1)对续断菊Sonchus asper L.Hill.根系分泌总有机酸、游离氨基酸、可溶性糖的影响,旨在探明根系分泌物对续断菊超积累Cd的影响。结果表明:Cd胁迫下续断菊根系分泌总有机酸、游离氨基酸和可溶性糖的质量浓度显著增加,同时,总有机酸、可溶性糖和游离氨基酸又促进了植株对Cd的吸收。随着Cd处理质量分数的增加,续断菊地上部和根部镉质量分数显著增加,90 d时续断菊地上部镉质量分数与可溶性糖、游离氨基酸的质量浓度呈极显著正相关,相关系数分别为0.999(P<0.01)和0.995(P<0.01),根部镉质量分数与可溶性糖、游离氨基酸的质量浓度也呈显著正相关,相关系数分别为0.998(P<0.01)和0.987(P<0.05);Cd对续断菊根系可溶性糖的分泌、游离氨基酸的合成有刺激作用,根系分泌的可溶性糖和游离氨基酸可能在续断菊累积镉的过程中有重要作用。

Effect of rock joint roughness on its cyclic shear behavior

Rock joints are often subjected to dynamic loads induced by earthquake and blasting during mining and rock cutting. Hence, cyclic shear load can be induced along the joints and it is important to evaluate the shear behavior of rock joint under this condition. In the present study, synthetic rock joints were prepared with plaster of Paris(Po P). Regular joints were simulated by keeping regular asperity with asperity angles of 15°-15° and 30°-30°, and irregular rock joints which are closer to natural joints were replicated by keeping the asperity angles of 15°-30° and 15°-45°. The sample size and amplitude of roughness were kept the same for both regular and irregular joints which were 298 mm×298 mm×125 mm and 5 mm, respectively. Shear test was performed on these joints using a large-scale direct shear testing machine by keeping the frequency and amplitude of shear load under constant cyclic condition with different normal stress values. As expected, the shear strength of rock joints increased with the increases in the asperity angle and normal load during the first cycle of shearing or static load. With the increase of the number of shear cycles, the shear strength decreased for all the asperity angles but the rate of reduction was more in case of high asperity angles. Test results indicated that shear strength of irregular joints was higher than that of regular joints at different cycles of shearing at low normal stress. Shearing and degradation of joint asperities on regular joints were the same between loading and unloading, but different for irregular joints. Shear strength and joint degradation were more significant on the slope of asperity with higher angles on the irregular joint until two angles of asperities became equal during the cycle of shearing and it started behaving like regular joints for subsequent cycles.

Assessing current faulting behaviors and seismic risk of the Anninghe-Zemuhe fault zone from seismicity parameters

Using the data of regional seismic network, this paper analyzes the current faulting behaviors of different segments of the Anninghe-Zemuhe fault zone, western Sichuan, and identifies the likely risky segments for potential large earthquakes. The authors map the probable asperities from the abnormally low b-value distribution, develop and employ a method for identifying current faulting behaviors of individual fault segment from the combinations of multiple seismicity parameter values, and make an effort to estimate the average recurrence intervals of character-istic earthquakes by using the parameters of magnitude-frequency relationship of the asperity segment. The result suggests that the studied fault zone contains 5 segments of different current faulting behaviors. Among them, the Mianning-Xichang segment of the Anninghe fault has been locked under high stress, its central part is probably an asperity with a relatively large scale. The Xichang-Puge segment of the Zemuhe fault displays very low seismicity under low stress. Both the locked segment and the low-seismicity segment can be outlined on the across-profile of relocated hypocenter depths. The Mianning-Xichang segment is identified to be the one with potential large earth-quake risk, for which the average recurrence interval between the latest M = 6.7 earthquake in 1952 and the next characteristic event is estimated to be 55 to 67 years, and the magnitude of the potential earthquake between 7.0 and 7.5. Also, it has been preliminarily suggested that for a certain fault segment, its faulting behaviors may change and evolve with time gradually.