This paper introduces a model for characterizing the contact behavior of irregular asperities,transforming it into a superposition of sinusoidal asperity contact behaviors.A new sinusoidal asperity model is developed ...This paper introduces a model for characterizing the contact behavior of irregular asperities,transforming it into a superposition of sinusoidal asperity contact behaviors.A new sinusoidal asperity model is developed for bilin-ear hardening under plane strain conditions.Empirical equations are proposed,considering geometric shapes,tangent modulus,and Young’s modulus.The frequency of asperity height is extracted through Fourier transform for irregular asperities.Contact area and pressure are predicted using the sinusoidal asperity model,and the behavior of irregular asperities is obtained by superimposing those with the first three frequencies.Experimen-tal validation is conducted with milling and knurling-formed asperities,showing good alignment between the model and experimental results.In rough surface models,the proposed irregular asperity model exhibits greater accuracy in predicting contact behavior than a single sinusoidal asperity when interference exceeds 10%of the amplitude.展开更多
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 curre...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.展开更多
Fault-slip taking place in underground mines occasionally causes severe damage to mine openings as a result of strong ground motion induced by seismic waves arising from fault-slip. It is indicated from previous studi...Fault-slip taking place in underground mines occasionally causes severe damage to mine openings as a result of strong ground motion induced by seismic waves arising from fault-slip. It is indicated from previous studies that intense seismic waves could be generated with the shock unloading of fault surface asperities during fault-slip. This study investigates the shock unloading with numerical simulation. A three-dimensional (3D) numerical model with idealized asperities is constructed with the help of discrete element code 3DEC. The idealization is conducted to particularly focus on simulating the shock unloading that previous numerical models, which replicate asperity degradation and crack development during the shear behavior of a joint surface in previous studies, fail to capture and simulate. With the numerical model, static and dynamic analyses are carried out to simulate unloading of asperities in the course of fault-slip. The results obtained from the dynamic analysis show that gradual stress release takes place around the center of the asperity tip at a rate of 45 MPa/ms for the base case, while an instantaneous stress release greater than 80 MPa occurs near the periphery of the asperity tip when the contact between the upper and lower asperities is lost. The instantaneous stress release becomes more intense in the vicinity of the asperity tip, causing tensile stress more than 20 MPa. It is deduced that the tensile stress could further increase if the numerical model is discretized more densely and analysis is carried out under stress conditions at a great depth. A model parametric study shows that in-situ stress state has a significant influence on the magnitude of the generated tensile stress. The results imply that the rapid stress release generating extremely high tensile stress on the asperity tip can cause intense seismic waves when it occurs at a great depth.展开更多
By dealing with strain and acoustic emission (AE) data for two parallelling faultss instability and failure with the same slip direction including asperities, the temporal-spatial evolution of strain and AE field dist...By dealing with strain and acoustic emission (AE) data for two parallelling faultss instability and failure with the same slip direction including asperities, the temporal-spatial evolution of strain and AE field distribution on the asperity of parallelling faults is analyzed. Furthermore the failure process of asperities and interaction among the asperities, i.e., positive and negative seismicity are discussed. Results show that instability and failure for the parallelling faults is a kind of negative seismicity.展开更多
An alternative extension to the Gaussian-beam expansion technique is provided to simplify the computation of the ~esnel field integral for rectangular symmetric sources. From a known result that the circle or rectangl...An alternative extension to the Gaussian-beam expansion technique is provided to simplify the computation of the ~esnel field integral for rectangular symmetric sources. From a known result that the circle or rectangle function is approximately decomposed into a sum of Gaussian functions, the cosine function is similarly expanded by the Bessel Fourier transform. Two expansions are together inserted in this field integral, it is then expressible in terms of the simple algebraic functions. As examples, the numerical results for the sound pressure field are presented for the uniform rectangular piston transducer, in a good agreement with those directly evaluated from the Fresnel integral. A wide applicability of this approach is discussed in treatment of the ultrasonic field radiation problem for a large and important group of piston sources in acoustics.展开更多
Low dielectric constant materials/Cu interconnects integration technology provides the direction as well as the challenges in the fabrication of integrated circuits(IC) wafers during copper electrochemical-mechanical ...Low dielectric constant materials/Cu interconnects integration technology provides the direction as well as the challenges in the fabrication of integrated circuits(IC) wafers during copper electrochemical-mechanical polishing(ECMP). These challenges arise primarily from the mechanical fragility of such dielectrics, in which the undesirable scratches are prone to produce. To mitigate this problem, a new model is proposed to predict the initiation of scratching based on the mechanical properties of passive layer and copper substrate. In order to deduce the ratio of the passive layer yield strength to the substrate yield strength and the layer thickness, the limit analysis solution of surface scratch under Berkovich indenter is used to analyze the nano-scratch experimental measurements. The modulus of the passive layer can be calculated by the nano-indentation test combined with the FEM simulation. It is found that the film modulus is about 30% of the substrate modulus. Various regimes of scratching are delineated by FEM modeling and the results are verified by experimental data.展开更多
Varied slope structures have different landslide initiation mechanisms.However,the role of interfacial asperities in controlling sliding initiation is unclear.This study develops a novel analytical method for interfac...Varied slope structures have different landslide initiation mechanisms.However,the role of interfacial asperities in controlling sliding initiation is unclear.This study develops a novel analytical method for interface mechanisms in practical landslide scenarios.Based on asperity theory,theoretical derivations were conducted,yielding an equilibrium equation grounded in the maximum shear length of asperities(Llimit).A method was proposed to evaluate slope stability by integrating Llimit with monitored displacement data.On-site displacement monitoring and slope state analysis determined the distribution range of asperities,providing theoretical support for slope remediation.Between July 2023 and July 2024,a large landslide located in Shagudu Town,Zhungeer Banner,Ordos,Inner Mongolia,underwent multiple deformation stages.This landslide seriously affected the operations of the Rongwu Expressway.This study initially examined the geological structure of the landslide through geological drilling and surveys.The distribution of asperities in the landslide area was determined through physical and mechanical tests and on-site geological mapping.A layout plan for landslide disposal and monitoring points was developed based on the interface landslide asperity theory.Analysis of 20 days of displacement monitoring data was used to determine asperity distribution and analyze slope movement,informed by asperity theory.This study pioneers the systematic application of asperity theory as the primary framework for analyzing and managing a large-scale interface landslide in engineering practice.Its rationality and effectiveness are rigorously demonstrated through the successful determination of the asperity distribution range within the Rongwu Expressway landslide,leading to effective remediation and favorable monitoring outcomes.This work not only validates the asperity theory for interface landslides but,crucially,establishes a novel,practical method for their stability assessment.Furthermore,a key theoretical advancement is the development of the concept that asperity-controlled main sliding surfaces can induce secondary sliding surfaces,based on field evidence and analysis.展开更多
A new expression for contact deformation is given, and the normal contact stiff- ness between single asperities is derived according to Hooke's law. A contact model between two ellipsoidal asperities is simulated by ...A new expression for contact deformation is given, and the normal contact stiff- ness between single asperities is derived according to Hooke's law. A contact model between two ellipsoidal asperities is simulated by the FE method, the result compared with the theoretical solution. It is found that the curves of the normal contact stiffness versus the included angle in the principal curvature direction show similar trends and evolve as a cosine feature. The effects of the parameters on normal contact stiffness are found to show that normal contact stiffness increases and reaches the upper limit gradually with an increase in these parameters.展开更多
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 roug...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.展开更多
Decadal forerunning seismic activity of magnitude Mw ≥ 5.0 is mapped for all 45 mainshocks of Mw 7.7 to 9.1 at subduction zones of the world from 1993 to mid 2020. The zones of high slip in nearly all great earthquak...Decadal forerunning seismic activity of magnitude Mw ≥ 5.0 is mapped for all 45 mainshocks of Mw 7.7 to 9.1 at subduction zones of the world from 1993 to mid 2020. The zones of high slip in nearly all great earthquakes were nearly quiescent beforehand and are identified as the sites of great asperities and zones of strong seismic coupling. Much forerunning activity occurred at smaller asperities along the peripheries of the rupture zones of many great and giant mainshocks. Those sizes of great asperities as ascertained from forerunning activity generally agree with the areas of high seismic slip as determined by others from geodetic and tide-gauge data and finite-source seismic modeling. Asperities are strong, well-coupled portions of plate interfaces. Different patterns of forerunning activity on time scales of about 5 to 45 years are attributed to either the sizes and spacing of asperities (or lack of). This permits many great asperities to be mapped decades before they rupture in great and giant shocks. Several poorly coupled subduction zones such as Java, Lesser Sunda, Marianas, Tonga and Kermadec are characterized by few great thrust earthquakes and little, in any forerunning activity. Rupture zones of many great and giant earthquakes are bordered either along strike, updip, or downdip by zones of lower plate coupling. Several bordering regions were sites of forerunning activity, aftershocks, and slow-slip events. The detection of forerunning and precursory activities of various kinds should be sought on the peripheries of great asperities as well as within zones of high co-seismic slip.展开更多
Decadal forerunning seismic activity is examined for very large, shallow earthquakes along strike-slip and intraplate faults of the world. It includes forerunning shocks of magnitude Mw ≥ 5.0 for 21 mainshocks of Mw ...Decadal forerunning seismic activity is examined for very large, shallow earthquakes along strike-slip and intraplate faults of the world. It includes forerunning shocks of magnitude Mw ≥ 5.0 for 21 mainshocks of Mw 7.5 to 8.6 from 1989 to 2020. Much forerunning activity occurred at what are interpreted to be smaller asperities along the peripheries of the rupture zones of great mainshocks at transform faults and subduction zones. Several great asperities as ascertained from forerunning activity agree with the areas of high seismic slip as determined by others using geodetic, mapping of surface faulting, and finite-source seismic modeling. The zones of high slip in many great earthquakes were nearly quiescent beforehand and are identified as the sites of great asperities. Asperities are strong, well-coupled portions of plate interfaces. Different patterns of forerunning activity on time scales of up to 45 years are attributed to the sizes and spacing of asperities (or lack of). This permits at least some great asperities along transform faults to be mapped decades before they rupture in great shocks. Rupture zones of many great mainshocks along transform faults are bordered either along strike, at depth or regionally by zones of lower plate coupling including either fault creep forerunning activity, aftershocks and/or slow-slip events. Forerunning activity to transforms in continental areas is more widespread spatially than that adjacent to oceanic transforms. The parts of the San Andreas fault themselves that ruptured in great California earthquakes during 1812, 1857 and 1906 have been very quiet since 1920;moderate to large shocks have been concentrated on their peripheries. The intraplate shocks studied, however, exhibited few if any forerunning events, which is attributed to the short period of time studied compared to their repeat times. The detection of forerunning and precursory activities for various time scales should be sought on the peripheries of great asperities and not just along the major faults themselves. This paper compliments that on decadal forerunning activity to great and giant earthquakes along subduction zones.展开更多
The contact pressure of an asperity is the fundamental parameter for constructing a single asperity contact model and analyzing the contact load and contact stiffness of the mechanical joint surfaces.In response to th...The contact pressure of an asperity is the fundamental parameter for constructing a single asperity contact model and analyzing the contact load and contact stiffness of the mechanical joint surfaces.In response to the defects of existing models in contact pressure analysis that have non-monotonic changes or do not conform to physical laws,this paper proposes a novel nonlinear contact model that can achieve continuous monotonic changes in contact pressure in accordance with physical laws.This model considers the continuous deformation of the asperity after load application,and can describe three deformation states:elastic,elastic-plastic,and plastic during the loading process.For the elastic and plastic deformation stages,this paper characterizes them using the classic Hertz elastic contact theory and the complete plastic contact theory.For the elastic-plastic deformation stage,the contact pressure is characterized using an empirical pressure function,and expressions for other contact parameters are derived.Furthermore,based on the principle of probability and statistics,the solution expression for the contact parameters of the rough interface is obtained,and a novel rough interface contact model is established.By comparing with existing experimental and simulation results,it is found that:1)The model in this paper achieves monotonic and continuous changes in contact pressure during the contact process of the asperity,while complying with physical laws;2)The new model is in good agreement with experimental and simulation results,verifying the universal effectiveness and correctness of the proposed model in solving the contact parameters of the rough interface;3)The new model has simplicity in expression and high computa-tional efficiency.展开更多
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...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.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.12272227 and 11732009).
摘要This paper introduces a model for characterizing the contact behavior of irregular asperities,transforming it into a superposition of sinusoidal asperity contact behaviors.A new sinusoidal asperity model is developed for bilin-ear hardening under plane strain conditions.Empirical equations are proposed,considering geometric shapes,tangent modulus,and Young’s modulus.The frequency of asperity height is extracted through Fourier transform for irregular asperities.Contact area and pressure are predicted using the sinusoidal asperity model,and the behavior of irregular asperities is obtained by superimposing those with the first three frequencies.Experimen-tal validation is conducted with milling and knurling-formed asperities,showing good alignment between the model and experimental results.In rough surface models,the proposed irregular asperity model exhibits greater accuracy in predicting contact behavior than a single sinusoidal asperity when interference exceeds 10%of the amplitude.
基金the Research Council of Norway(Grant No.244029)the project‘Stable dams’,FORMAS(Grant No.2019e01236)+1 种基金the project‘Improved safety assessment of concrete dams’,and SVC(Grant No.VKU32019)the project‘Safe dams’,that supported the development of the research presented in this article.
摘要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.
基金financially supported by the Natural Science and Engineering Research Council of Canada(NSERC) in partnership with Vale Ltd.-Sudbury Operations,Canada,under the Collaborative Research and Development Program
摘要Fault-slip taking place in underground mines occasionally causes severe damage to mine openings as a result of strong ground motion induced by seismic waves arising from fault-slip. It is indicated from previous studies that intense seismic waves could be generated with the shock unloading of fault surface asperities during fault-slip. This study investigates the shock unloading with numerical simulation. A three-dimensional (3D) numerical model with idealized asperities is constructed with the help of discrete element code 3DEC. The idealization is conducted to particularly focus on simulating the shock unloading that previous numerical models, which replicate asperity degradation and crack development during the shear behavior of a joint surface in previous studies, fail to capture and simulate. With the numerical model, static and dynamic analyses are carried out to simulate unloading of asperities in the course of fault-slip. The results obtained from the dynamic analysis show that gradual stress release takes place around the center of the asperity tip at a rate of 45 MPa/ms for the base case, while an instantaneous stress release greater than 80 MPa occurs near the periphery of the asperity tip when the contact between the upper and lower asperities is lost. The instantaneous stress release becomes more intense in the vicinity of the asperity tip, causing tensile stress more than 20 MPa. It is deduced that the tensile stress could further increase if the numerical model is discretized more densely and analysis is carried out under stress conditions at a great depth. A model parametric study shows that in-situ stress state has a significant influence on the magnitude of the generated tensile stress. The results imply that the rapid stress release generating extremely high tensile stress on the asperity tip can cause intense seismic waves when it occurs at a great depth.
基金State Key Basic Research Development and Programming Project Mechanism and Prediction of Continental Strong Earthquake (G19980407) and Foundation of Laboratory of Tectonophysics China Seismological Bureau and State Natural Science Foundation of China (
摘要By dealing with strain and acoustic emission (AE) data for two parallelling faultss instability and failure with the same slip direction including asperities, the temporal-spatial evolution of strain and AE field distribution on the asperity of parallelling faults is analyzed. Furthermore the failure process of asperities and interaction among the asperities, i.e., positive and negative seismicity are discussed. Results show that instability and failure for the parallelling faults is a kind of negative seismicity.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11074038 and 11374051
摘要An alternative extension to the Gaussian-beam expansion technique is provided to simplify the computation of the ~esnel field integral for rectangular symmetric sources. From a known result that the circle or rectangle function is approximately decomposed into a sum of Gaussian functions, the cosine function is similarly expanded by the Bessel Fourier transform. Two expansions are together inserted in this field integral, it is then expressible in terms of the simple algebraic functions. As examples, the numerical results for the sound pressure field are presented for the uniform rectangular piston transducer, in a good agreement with those directly evaluated from the Fresnel integral. A wide applicability of this approach is discussed in treatment of the ultrasonic field radiation problem for a large and important group of piston sources in acoustics.
基金Project(50975058) supported by the National Natural Science Foundation of China
摘要Low dielectric constant materials/Cu interconnects integration technology provides the direction as well as the challenges in the fabrication of integrated circuits(IC) wafers during copper electrochemical-mechanical polishing(ECMP). These challenges arise primarily from the mechanical fragility of such dielectrics, in which the undesirable scratches are prone to produce. To mitigate this problem, a new model is proposed to predict the initiation of scratching based on the mechanical properties of passive layer and copper substrate. In order to deduce the ratio of the passive layer yield strength to the substrate yield strength and the layer thickness, the limit analysis solution of surface scratch under Berkovich indenter is used to analyze the nano-scratch experimental measurements. The modulus of the passive layer can be calculated by the nano-indentation test combined with the FEM simulation. It is found that the film modulus is about 30% of the substrate modulus. Various regimes of scratching are delineated by FEM modeling and the results are verified by experimental data.
基金supported by the National Natural Science Foundation of China(Grant No:12362034)the Science and Technology Plan Project of Inner Mongolia Autonomous Region(Grant No:2022YFSH0047)Talent Project of Inner Mongolia Autonomous Region(DC2300001439,DC2200000913).
摘要Varied slope structures have different landslide initiation mechanisms.However,the role of interfacial asperities in controlling sliding initiation is unclear.This study develops a novel analytical method for interface mechanisms in practical landslide scenarios.Based on asperity theory,theoretical derivations were conducted,yielding an equilibrium equation grounded in the maximum shear length of asperities(Llimit).A method was proposed to evaluate slope stability by integrating Llimit with monitored displacement data.On-site displacement monitoring and slope state analysis determined the distribution range of asperities,providing theoretical support for slope remediation.Between July 2023 and July 2024,a large landslide located in Shagudu Town,Zhungeer Banner,Ordos,Inner Mongolia,underwent multiple deformation stages.This landslide seriously affected the operations of the Rongwu Expressway.This study initially examined the geological structure of the landslide through geological drilling and surveys.The distribution of asperities in the landslide area was determined through physical and mechanical tests and on-site geological mapping.A layout plan for landslide disposal and monitoring points was developed based on the interface landslide asperity theory.Analysis of 20 days of displacement monitoring data was used to determine asperity distribution and analyze slope movement,informed by asperity theory.This study pioneers the systematic application of asperity theory as the primary framework for analyzing and managing a large-scale interface landslide in engineering practice.Its rationality and effectiveness are rigorously demonstrated through the successful determination of the asperity distribution range within the Rongwu Expressway landslide,leading to effective remediation and favorable monitoring outcomes.This work not only validates the asperity theory for interface landslides but,crucially,establishes a novel,practical method for their stability assessment.Furthermore,a key theoretical advancement is the development of the concept that asperity-controlled main sliding surfaces can induce secondary sliding surfaces,based on field evidence and analysis.
基金Project supported by the National Basic Research Program of China(973 Program)(No.2009CB724406)
摘要A new expression for contact deformation is given, and the normal contact stiff- ness between single asperities is derived according to Hooke's law. A contact model between two ellipsoidal asperities is simulated by the FE method, the result compared with the theoretical solution. It is found that the curves of the normal contact stiffness versus the included angle in the principal curvature direction show similar trends and evolve as a cosine feature. The effects of the parameters on normal contact stiffness are found to show that normal contact stiffness increases and reaches the upper limit gradually with an increase in these parameters.
基金support from the National Natural Science Foundation of China(Grant Nos.51991392 and 42293355).
摘要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.
摘要Decadal forerunning seismic activity of magnitude Mw ≥ 5.0 is mapped for all 45 mainshocks of Mw 7.7 to 9.1 at subduction zones of the world from 1993 to mid 2020. The zones of high slip in nearly all great earthquakes were nearly quiescent beforehand and are identified as the sites of great asperities and zones of strong seismic coupling. Much forerunning activity occurred at smaller asperities along the peripheries of the rupture zones of many great and giant mainshocks. Those sizes of great asperities as ascertained from forerunning activity generally agree with the areas of high seismic slip as determined by others from geodetic and tide-gauge data and finite-source seismic modeling. Asperities are strong, well-coupled portions of plate interfaces. Different patterns of forerunning activity on time scales of about 5 to 45 years are attributed to either the sizes and spacing of asperities (or lack of). This permits many great asperities to be mapped decades before they rupture in great and giant shocks. Several poorly coupled subduction zones such as Java, Lesser Sunda, Marianas, Tonga and Kermadec are characterized by few great thrust earthquakes and little, in any forerunning activity. Rupture zones of many great and giant earthquakes are bordered either along strike, updip, or downdip by zones of lower plate coupling. Several bordering regions were sites of forerunning activity, aftershocks, and slow-slip events. The detection of forerunning and precursory activities of various kinds should be sought on the peripheries of great asperities as well as within zones of high co-seismic slip.
摘要Decadal forerunning seismic activity is examined for very large, shallow earthquakes along strike-slip and intraplate faults of the world. It includes forerunning shocks of magnitude Mw ≥ 5.0 for 21 mainshocks of Mw 7.5 to 8.6 from 1989 to 2020. Much forerunning activity occurred at what are interpreted to be smaller asperities along the peripheries of the rupture zones of great mainshocks at transform faults and subduction zones. Several great asperities as ascertained from forerunning activity agree with the areas of high seismic slip as determined by others using geodetic, mapping of surface faulting, and finite-source seismic modeling. The zones of high slip in many great earthquakes were nearly quiescent beforehand and are identified as the sites of great asperities. Asperities are strong, well-coupled portions of plate interfaces. Different patterns of forerunning activity on time scales of up to 45 years are attributed to the sizes and spacing of asperities (or lack of). This permits at least some great asperities along transform faults to be mapped decades before they rupture in great shocks. Rupture zones of many great mainshocks along transform faults are bordered either along strike, at depth or regionally by zones of lower plate coupling including either fault creep forerunning activity, aftershocks and/or slow-slip events. Forerunning activity to transforms in continental areas is more widespread spatially than that adjacent to oceanic transforms. The parts of the San Andreas fault themselves that ruptured in great California earthquakes during 1812, 1857 and 1906 have been very quiet since 1920;moderate to large shocks have been concentrated on their peripheries. The intraplate shocks studied, however, exhibited few if any forerunning events, which is attributed to the short period of time studied compared to their repeat times. The detection of forerunning and precursory activities for various time scales should be sought on the peripheries of great asperities and not just along the major faults themselves. This paper compliments that on decadal forerunning activity to great and giant earthquakes along subduction zones.
基金Supported by Ministry of Science and Technology"the Belt and Road"Innovative Talents Exchange Foreign Experts Project(Grant No.DL2023183003L)Natural Science Basic Research Plan in Shaanxi Province of China(Grant No.2022JM-195)the National Foreign Experts Program(Grant No.H202501060).
摘要The contact pressure of an asperity is the fundamental parameter for constructing a single asperity contact model and analyzing the contact load and contact stiffness of the mechanical joint surfaces.In response to the defects of existing models in contact pressure analysis that have non-monotonic changes or do not conform to physical laws,this paper proposes a novel nonlinear contact model that can achieve continuous monotonic changes in contact pressure in accordance with physical laws.This model considers the continuous deformation of the asperity after load application,and can describe three deformation states:elastic,elastic-plastic,and plastic during the loading process.For the elastic and plastic deformation stages,this paper characterizes them using the classic Hertz elastic contact theory and the complete plastic contact theory.For the elastic-plastic deformation stage,the contact pressure is characterized using an empirical pressure function,and expressions for other contact parameters are derived.Furthermore,based on the principle of probability and statistics,the solution expression for the contact parameters of the rough interface is obtained,and a novel rough interface contact model is established.By comparing with existing experimental and simulation results,it is found that:1)The model in this paper achieves monotonic and continuous changes in contact pressure during the contact process of the asperity,while complying with physical laws;2)The new model is in good agreement with experimental and simulation results,verifying the universal effectiveness and correctness of the proposed model in solving the contact parameters of the rough interface;3)The new model has simplicity in expression and high computa-tional efficiency.
基金The first author would like to thank the supports of the NARGS, IRGS and AAS grants of Australia, and the National Science Foundation grants (No. 51574060 and No. 51079017) of China, in which the first author is the intemational collaborator. The academic visits of the third and fourth authors to the University of Tasmania are partly supported by a PhD visiting scholarship and an academic visiting scholarship, respectively, provided by the China Scholarship Council, which are greatly appreciated.
摘要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.