Shear strength degradation of gas hydrate-bearing sediment due to partial hydrate dissociation

Extraction of methane hydrate from subseafloor reservoir may potentially trigger seabed slides and induce subsidence.To address the problems,it is crucial to properly characterize the phase equilibrium condition of pore hydrate and the shear strength of the soil.As one of the key constitutive components,the phase equilibrium condition enforces a constraint over pore gas pressure,temperature and unhydrated water content.Such a constraint,however,has been traditionally ignored in analyzing the mechanical behavior of hydrate-bearing soil.In this paper,a series of stepwise hydrate dissociation tests was performed,and the phase equilibrium condition of pore hydrate was determined,providing an effective way to evaluate the unhydrated water content during hydrate dissociation.Meanwhile,a series of direct shear tests was also conducted to explore the shear strength characteristics of the soil.It is shown that the shear strength of the hydrate-bearing soil can be significantly influenced by pore gas pressure,unhydrated water content,hydrate saturation and several other factors.In particular,the measured shear strength depends upon the initial water content of the sample,pointing to a potential problem that the shear strength could be wrongly determined if not properly interpreted.A shear strength criterion,which enforces the equilibrium condition of pore hydrate,is developed for hydrate-bearing soil,establishing a link between the equilibrium condition and the shear strength.The proposed equation describes well the shear strength characteristics of hydrate-bearing soils,remarkably unifying the effects of pore pressure,temperature,water content and hydrate saturation.

Undrained vane shear strength of sand-foam mixtures subjected to different shear rates

The shear strength of sand-foam mixtures plays a crucial role in ensuring successful earth pressure balance(EPB)shield tunneling.Since the sand-foam mixtures are constantly sheared by the cutterhead and the screw conveyor with varied rotation speeds during tunneling,it is non-trivial to investigate the effect of shear rates on the undrained shear strength of sand-foam mixtures under chamber pressures to extend the understanding on the tunneling process.This study conducted a series of pressurized vane shear tests to investigate the role of shear rates on the peak and residual strengths of sand-foam mixtures at different pore states.Different from the shear-rate characteristics of natural sands or clay,the results showed that the peak strength of sand-foam mixtures under high vertical total stress(σ_(v)≥200 kPa)and low foam injection ratio(FIR30%)decreased with the increase in shear rate.Otherwise,the peak strength was not measurably affected by shear rates.The sand-foam mixtures in the residual state resembled low-viscous fluid with yield stress and the residual strength increased slightly with shear rates.In addition,the peak and residual strengths were approximately linear with vertical effective stress regardless of the total stress and FIR.The peak effective internal friction angle remained almost invariant in a low shear rate(γ′<0.25 s1)but decreased when the shear rate continued increasing.The residual effective internal friction angle was lower than the peak counterpart and insensitive to shear rates.This study unveiled the role of shear rates in the undrained shear strength of sand-foam mixtures with various FIRs and vertical total stresses.The findings can extend the understanding of the rate-dependent shear characteristics of conditioned soils and guide the decision-making of soil conditioning schemes in the EPB shield tunneling practice.

Anisotropic time-dependent behaviors of shale under direct shearing and associated empirical creep models

Understanding the anisotropic creep behaviors of shale under direct shearing is a challenging issue.In this context,we conducted shear-creep and steady-creep tests on shale with five bedding orientations (i.e.0°,30°,45°,60°,and 90°),under multiple levels of direct shearing for the first time.The results show that the anisotropic creep of shale exhibits a significant stress-dependent behavior.Under a low shear stress,the creep compliance of shale increases linearly with the logarithm of time at all bedding orientations,and the increase depends on the bedding orientation and creep time.Under high shear stress conditions,the creep compliance of shale is minimal when the bedding orientation is 0°,and the steady-creep rate of shale increases significantly with increasing bedding orientations of 30°,45°,60°,and 90°.The stress-strain values corresponding to the inception of the accelerated creep stage show an increasing and then decreasing trend with the bedding orientation.A semilogarithmic model that could reflect the stress dependence of the steady-creep rate while considering the hardening and damage process is proposed.The model minimizes the deviation of the calculated steady-state creep rate from the observed value and reveals the behavior of the bedding orientation's influence on the steady-creep rate.The applicability of the five classical empirical creep models is quantitatively evaluated.It shows that the logarithmic model can well explain the experimental creep strain and creep rate,and it can accurately predict long-term shear creep deformation.Based on an improved logarithmic model,the variations in creep parameters with shear stress and bedding orientations are discussed.With abovementioned findings,a mathematical method for constructing an anisotropic shear creep model of shale is proposed,which can characterize the nonlinear dependence of the anisotropic shear creep behavior of shale on the bedding orientation.

Experimental investigation on shear strength deterioration at the interface between different rock types under cyclic loading

The shear strength deterioration of bedding planes between different rock types induced by cyclic loading is vital to reasonably evaluate the stability of soft and hard interbedded bedding rock slopes under earthquake;however,rare work has been devoted to this subject due to lack of attention.In this study,experimental investigations on shear strength weakening of discontinuities with different joint wall material(DDJM)under cyclic loading were conducted by taking the interface between siltstone and mudstone in the Shaba slope of Yunnan Province,China as research objects.A total of 99 pairs of similar material samples of DDJM(81 pairs)and discontinuities with identical joint wall material(DIJM)(18 pairs)were fabricated by inserting plates,engraved with typical surface morphology obtained by performing three-dimensional laser scanning on natural DDJMs sampled from field,into mold boxes.Cyclic shear tests were conducted on these samples to study their shear strength changes with the cyclic number considering the effects of normal stress,joint surface morphology,shear displacement amplitude and shear rate.The results indicate that the shear stress vs.shear displacement curves under each shear cycle and the peak shear strength vs.cyclic number curves of the studied DDJMs are between those of DIJMs with siltstone and mudstone,while closer to those of DIJMs with mudstone.The peak shear strengths of DDJMs exhibit an initial rapid decline followed by a gradual decrease with the cyclic number and the decrease rate varies from 6%to 55.9%for samples with varied surface morphology under different testing conditions.The normal stress,joint surface morphology,shear displacement amplitude and shear rate collectively influence the shear strength deterioration of DDJM under cyclic shear loading,with the degree of influence being greater for larger normal stress,rougher surface morphology,larger shear displacement amplitude and faster shear rate.

Anisotropic shearing mechanism of Kangding slate:Experimental investigation and numerical analysis

The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles.Direct shear tests were conducted on cubic slate samples with foliation angles of 0°,30°,45°,60°,and 90°.The effect of foliation angles on failure patterns,acoustic emission(AE)characteristics,and shear strength parameters was analyzed.Based on AE characteristics,the slate failure process could be divided into four stages:quiet period,step-like increasing period,dramatic increasing period,and remission period.A new empirical expression of cohesion for layered rock was proposed,which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments.The comparative analysis demonstrated that the new expression has better prediction ability than other expressions.The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method(FDEM),and it was found to align well with the experimental results.Considering both computational efficiency and accuracy,it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations.To balance the relationship between the number of elements and the simulation results in the direct shear simulations,the recommended element size is 1 mm.

Experimental study on the effect of unloading rate on gneiss rockburst

Rockburst are often encountered in tunnel construction due to the complex geological conditions.To study the influence of unloading rate on rockburst,gneiss rockburst experiments were conducted under three groups of unloading rates.A high-speed photography system and acoustic emission(AE)system were used to monitor the entire process of rockburst process in real-time.The results show that the intensity of gneiss rockburst decreases with decrease of unloading rate,which is manifested as the reduction of AE energy and fragments ejection velocity.The mechanisms are proposed to explain this effect:(i)The reduction of unloading rate changes the crack propagation mechanism in the process of rockburst.This makes the rockbursts change from the tensile failure mechanism at high unloading rate to the tension-shear mixed failure mechanism at low unloading rate,and more energy released in the form of shear crack propagation.Then,less strain energy is converted into kinetic energy of fragments ejection.(ii)Less plate cracking degree of gneiss has taken shape due to decrease of unloading rate,resulting in the destruction of rockburst incubation process.The enlightenments of reducing the unloading rate for the project are also described quantitatively.The rockburst magnitude is reduced from the medium magnitude at the unloading rate of 0.1 MPa/s to the slight magnitude at the unloading rate of 0.025 MPa/s,which was judged by the ejection velocity.

Influence of strain rate on grain refinement and texture evolution under complex shear stress conditions

The effects of different complex shear stress conditions on grain refinement and texture evolution of Mg-13Gd-4Y-2Zn-0.5Zr alloy were investigated.With increasing strain rate,the average grain size of compression-shear(CS)and compression-torsion(CT)samples are decreased,and the grain size of dynamic recrystallization(DRX)grains is also decreased.This is because that the precipitation number ofβphases is increased,and the hindering effect on grain growth can be significantly enhanced.The DRX fractions of CS and CT samples are decreased with increased strain rate.The low DRX fraction at high strain rate is related to the insufficient time for grains to nucleate.The DRX process can be promoted by the PSN mechanism of second phases,and the grain growth can be restricted by the pinning effect.At the same time,the texture strength is enhanced as the strain rate increased.Besides,the kinking degree of lamellar long-period stacking ordered(LPSO)phases is increased.Under complex shear stress conditions,non-basal slip,especially pyramidal slip,is easily activated and the texture is deflected greatly.Compared with the CS samples,CT samples have smaller average grain size,higher DRX fraction,and lower texture strength for a certain strain rate.This is because that the equivalent stress of the CT sample is larger,the stress triaxiality is smaller,so more serious dislocations are piled up near grain boundaries and second phases.At the same time,since CT sample was sheared with torsion,the dislocation movement path can be called“rotational dislocation accumulation”,and the longer distribution path of the CT sample is generated,so more sub-grains and low-angle grain boundaries(LAGBs)are formed.Compared with the CS sample,more huge-angle grain boundaries(HAGBs)and DRX grains are formed from grain boundary to grain interior,so better grain refinement effect is achieved.

The Probability Density Function Related to Shallow Cumulus Entrainment Rate and Its Influencing Factors in a Large-Eddy Simulation

The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.

Failure transition of shear-to-dilation band of rock salt under triaxial stresses

Great potential of underground gas/energy storage in salt caverns seems to be a promising solution to support renewable energy.In the underground storage method,the operating cycle unfortunately may reach up to daily or even hourly,which generates complicated pressures on the salt cavern.Furthermore,the mechanical behavior of rock salt may change and present distinct failure characteristics under different stress states,which affects the performance of salt cavern during the time period of full service.To reproduce a similar loading condition on the cavern surrounding rock mass,the cyclic triaxial loading/unloading tests are performed on the rock salt to explore the mechanical transition behavior and failure characteristics under different confinement.Experimental results show that the rock salt samples pre-sent a diffused shear failure band with significant bulges at certain locations in low confining pressure conditions(e.g.5 MPa,10 MPa and 15 MPa),which is closely related to crystal misorientation and grain boundary sliding.Under the elevated confinement(e.g.20 MPa,30 MPa and 40 MPa),the dilation band dominates the failure mechanism,where the large-size halite crystals are crushed to be smaller size and new pores are developing.The failure transition mechanism revealed in the paper provides additional insight into the mechanical performance of salt caverns influenced by complicated stress states.

Prediction model for corrosion rate of low-alloy steels under atmospheric conditions using machine learning algorithms

This work constructed a machine learning(ML)model to predict the atmospheric corrosion rate of low-alloy steels(LAS).The material properties of LAS,environmental factors,and exposure time were used as the input,while the corrosion rate as the output.6 dif-ferent ML algorithms were used to construct the proposed model.Through optimization and filtering,the eXtreme gradient boosting(XG-Boost)model exhibited good corrosion rate prediction accuracy.The features of material properties were then transformed into atomic and physical features using the proposed property transformation approach,and the dominant descriptors that affected the corrosion rate were filtered using the recursive feature elimination(RFE)as well as XGBoost methods.The established ML models exhibited better predic-tion performance and generalization ability via property transformation descriptors.In addition,the SHapley additive exPlanations(SHAP)method was applied to analyze the relationship between the descriptors and corrosion rate.The results showed that the property transformation model could effectively help with analyzing the corrosion behavior,thereby significantly improving the generalization ability of corrosion rate prediction models.

Federated Learning Model for Auto Insurance Rate Setting Based on Tweedie Distribution

In the assessment of car insurance claims,the claim rate for car insurance presents a highly skewed probability distribution,which is typically modeled using Tweedie distribution.The traditional approach to obtaining the Tweedie regression model involves training on a centralized dataset,when the data is provided by multiple parties,training a privacy-preserving Tweedie regression model without exchanging raw data becomes a challenge.To address this issue,this study introduces a novel vertical federated learning-based Tweedie regression algorithm for multi-party auto insurance rate setting in data silos.The algorithm can keep sensitive data locally and uses privacy-preserving techniques to achieve intersection operations between the two parties holding the data.After determining which entities are shared,the participants train the model locally using the shared entity data to obtain the local generalized linear model intermediate parameters.The homomorphic encryption algorithms are introduced to interact with and update the model intermediate parameters to collaboratively complete the joint training of the car insurance rate-setting model.Performance tests on two publicly available datasets show that the proposed federated Tweedie regression algorithm can effectively generate Tweedie regression models that leverage the value of data fromboth partieswithout exchanging data.The assessment results of the scheme approach those of the Tweedie regressionmodel learned fromcentralized data,and outperformthe Tweedie regressionmodel learned independently by a single party.

Construction of a Computational Scheme for the Fuzzy HIV/AIDS Epidemic Model with a Nonlinear Saturated Incidence Rate

This work aimed to construct an epidemic model with fuzzy parameters.Since the classical epidemic model doesnot elaborate on the successful interaction of susceptible and infective people,the constructed fuzzy epidemicmodel discusses the more detailed versions of the interactions between infective and susceptible people.Thenext-generation matrix approach is employed to find the reproduction number of a deterministic model.Thesensitivity analysis and local stability analysis of the systemare also provided.For solving the fuzzy epidemic model,a numerical scheme is constructed which consists of three time levels.The numerical scheme has an advantage overthe existing forward Euler scheme for determining the conditions of getting the positive solution.The establishedscheme also has an advantage over existing non-standard finite difference methods in terms of order of accuracy.The stability of the scheme for the considered fuzzy model is also provided.From the plotted results,it can beobserved that susceptible people decay by rising interaction parameters.

Mechanical behavior of nanorubber reinforced epoxy over a wide strain rate loading

Nanorubber/epoxy composites containing 0,2,6 and 10 wt%nanorubber are subjected to uniaxial compression over a wide range of strain rate from 8×10^(-4) s^(-1) to~2×10^(4) s^(-1).Unexpectedly,their strain rate sensitivity and strain hardening index increase with increasing nanorubber content.Potential mechanisms are proposed based on numerical simulations using a unit cell model.An increase in the strain rate sensitivity with increasing nanorubber content results from the fact that the nanorubber becomes less incompressible at high strain,generating a higher hydro-static pressure.Adiabatic shear localization starts to occur in the epoxy under a strain rate of 22,000 s^(-1) when the strain exceeds 0.35.The presence of nanorubber in the epoxy reduces adiabatic shear localization by preventing it from propagating.

State-of-the-art on the anchorage performance of rock bolts subjected to shear load

Rock bolts are extensively utilized in underground engineering as a means of offering support and stability to rock masses in tunnels,mines,and other underground structures.In environments of high ground stress,faults or weak zones can frequently arise in rock formations,presenting a significant challenge for engineering and potentially leading to underground engineering collapse.Rock bolts serve as a crucial structural element for the transmission of tensile stress and are capable of withstanding shear loads to prevent sliding of weak zones within rock mass.Therefore,a complete understanding of the behavior of rock bolts subjected to shear loads is essential.This paper presents a state-of-the-art review of the research progress of rock bolts subjected to shear load in three categories:experiment,numerical simulation,and analytical model.The review focuses on the research studies and developments in this area since the 1970s,providing a comprehensive overview of numerous factors that influence the anchorage performance of rock bolts.These factors include the diameter and angle of the rock bolt installation,rock strength,grouting material,bolt material,borehole diameter,rock bolt preload,normal stress,joint surface roughness and joint expansion angle.The paper reviews the improvement of mechanical parameter setting in numerical simulation of rock bolt shear.Furthermore,it delves into the optimization of the analytical model concerning rock bolt shear theory,approached from the perspectives of both Elastic foundation beam theory coupled with Elastoplasticity theory and Structural mechanic methods.The significance of this review lies in its ability to provide insights into the mechanical behavior of rock bolts.The paper also highlights the limitations of current research and guidelines for further research of rock bolts.

Mechanical behaviors of backfill-rock composites: Physical shear test and back-analysis

The shear behavior of backfill-rock composites is crucial for mine safety and the management of surface subsidence.For exposing the shear failure mechanism of backfill-rock composites,we conducted shear tests on backfill-rock composites under three constant normal loads,compared with the unfilled rock.To investigate the macro-and meso-failure characteristics of the samples in the shear tests,the cracking behavior of samples was recorded by a high-speed camera and acoustic emission monitoring.In parallel with the experimental test,the numerical models of backfill-rock composites and unfilled rock were established using the discrete element method to analyze the continuous-discontinuous shearing process.Based on the damage mechanics and statistics,a novel shear constitutive model was proposed to describe mechanical behavior.The results show that backfill-rock composites had a special bimodal phenomenon of shearing load-deformation curve,i.e.the first shearing peak corresponded to rock break and the second shearing peak induced by the broken of aeolian sand-cement/fly ash paste backfill.Moreover,the shearing characteristic curves of the backfill-rock composites could be roughly divided into four stages,i.e.the shear failure of the specimens experienced:stage I:stress concentration;stage II:crack propagation;stage III:crack coalescence;stage IV:shearing friction.The numerical simulation shows that the existence of aeolian sand-cement/fly ash paste backfill inevitably altered the coalescence type and failure mode of the specimens and had a strengthening effect on the shear strength of backfillrock composites.Based on damage mechanics and statistics,a shear constitutive model was proposed to describe the shear fracture characteristics of specimens,especially the bimodal phenomenon.Finally,the micro-and meso-mechanisms of shear failure were discussed by combining the micro-test and numerical results.The research can advance the better understanding of the shear behavior of backfill-rock composites and contribute to the safety of mining engineering.

Shear band evolution and acoustic emission characteristics of sandstone containing non-persistent flaws

Direct shear tests were conducted on sandstone specimens under different constant normal stresses to study the coalescence of cracks between non-persistent flaws and the shear sliding characteristics of the shear-formed fault.Digital image correlation and acoustic emission(AE)techniques were used to monitor the evolution of shear bands at the rock bridge area and microcracking behaviors.The experimental results revealed that the shear stresses corresponding to the peak and sub-peak in the stressdisplacement curve are significantly affected by the normal stress.Strain localization bands emerged at both the tip of joints and the rock bridge,and their extension and interaction near the peak stress caused a surge in the AE hit rate and a significant decrease in the AE b value.Short and curvilinear strain bands were detected at low normal stress,while high normal stress generally led to more microcracking events and longer coplanar cracks at the rock bridge area.Furthermore,an increase in normal stress resulted in a higher AE count rate and more energetic AE events during friction sliding along the shearformed fault.It was observed that the elastic energy released during the crack coalescence at the prepeak stage was much greater than that released during friction sliding at the post-peak stage.More than 75%of AE events were located in the low-frequency band(0e100 kHz),and this proportion continued to rise with increasing normal stress.Moreover,more AE events of low AF value and high RA value were observed in specimens subjected to high normal stress,indicating that greater normal stress led to more microcracks of shear nature.

Retraction to: Shear strength of extremely altered serpentinites based on degree of saturation (Ankara, Turkey)

The Editor in Chief has retracted this article.After publication,concerns were raised over potentially inconsistent data-for example,the author mentioned residual soil with silt fraction,but Fig.5 shows even the smallest grains are beyond silt grains.Further analysis showed that Fig.6 and Fig.9 show both loose and dense soils exhibiting very similar dilatance behaviour,which is not typical.The author shared raw data,but it was incomplete and could not be verified.The Editor in Chief,therefore,has lost confidence in the findings of this article.The author does not agree to this retraction.

Seismic performance evaluation of hybrid coupled shear wall system with shear and flexural fuse-type steel coupling beams

Replaceable flexural and shear fuse-type coupling beams are used in hybrid coupled shear wall(HCSW)systems,enabling concrete buildings to be promptly recovered after severe earthquakes.This study aimed to analytically evaluate the seismic behavior of flexural and shear fuse beams situated in short-,medium-and high-rise RC buildings that have HCSWs.Three building groups hypothetically located in a high seismic hazard zone were studied.A series of 2D nonlinear time history analyses was accomplished in OpenSees,using the ground motion records scaled at the design basis earthquake level.It was found that the effectiveness of fuses in HCSWs depends on various factors such as size and scale of the building,allowable rotation value,inter-story drift ratio,residual drift quantity,energy dissipation value of the fuses,etc.The results show that shear fuses better meet the requirements of rotations and drifts.In contrast,flexural fuses dissipate more energy,but their sectional stiffness should increase to meet other requirements.It was concluded that adoption of proper fuses depends on the overall scale of the building and on how associated factors are considered.

Photoacoustic elastography based on laser-excited shear wave

Elastography can be used as a diagnostic method for quantitative characterization of tissue hardness information and thus,differential changes in pathophysiological states of tissues.In this study,we propose a new method for shear wave elastography(SWE)based on laser-excited shear wave,called photoacoustic shear wave elastography(PASWE),which combines photoacoustic(PA)technology with ultrafast ultrasound imaging.By using a focused laser to excite shear waves and ultrafast ultrasonic imaging for detection,high-frequency excitation of shear waves and noncontact elastic imaging can be realized.The laser can stimulate the tissue with the light absorption characteristic to produce the thermal expansion,thus producing the shear wave.The frequency of shear wave induced by laser is higher and the frequency band is wider.By tracking the propagation of shear wave,Young’s modulus of tissue is reconstructed in the whole shear wave propagation region to further evaluate the elastic information of tissue.The feasibility of the method is verified by experiments.Compared with the experimental results of supersonic shear imaging(SSI),it is proved that the method can be used for quantitative elastic imaging of the phantoms.In addition,compared with the SSI method,this method can realize the noncontact excitation of the shear wave,and the frequency of the shear wave excited by the laser is higher than that of the acoustic radiation force(ARF),so the spatial resolution is higher.Compared to the traditional PA elastic imaging method,this method can obtain a larger imaging depth under the premise of ensuring the imaging resolution,and it has potential application value in the clinical diagnosis of diseases requiring noncontact quantitative elasticity.

On the calibration of a shear stress criterion for rock joints to represent the full stress-strain profile

Conventional numerical solutions developed to describe the geomechanical behavior of rock interfaces subjected to differential load emphasize peak and residual shear strengths.The detailed analysis of preand post-peak shear stress-displacement behavior is central to various time-dependent and dynamic rock mechanic problems such as rockbursts and structural instabilities in highly stressed conditions.The complete stress-displacement surface(CSDS)model was developed to describe analytically the pre-and post-peak behavior of rock interfaces under differential loads.Original formulations of the CSDS model required extensive curve-fitting iterations which limited its practical applicability and transparent integration into engineering tools.The present work proposes modifications to the CSDS model aimed at developing a comprehensive and modern calibration protocol to describe the complete shear stressdisplacement behavior of rock interfaces under differential loads.The proposed update to the CSDS model incorporates the concept of mobilized shear strength to enhance the post-peak formulations.Barton’s concepts of joint roughness coefficient(JRC)and joint compressive strength(JCS)are incorporated to facilitate empirical estimations for peak shear stress and normal closure relations.Triaxial/uniaxial compression test and direct shear test results are used to validate the updated model and exemplify the proposed calibration method.The results illustrate that the revised model successfully predicts the post-peak and complete axial stressestrain and shear stressedisplacement curves for rock joints.