A coupled model including the combination of pressure solution model and stress corrosion model is proposed to simulate the observed strength recovery of limestone fracture. Thus

the coupled model can consider both pressure solution and stress corrosion in modification of the fracture contact area

which is used to predict the frictional healing of the fracture. The simulation results show that the pressure solution model underestimates frictional healing during short hold times

and with extending the hold times

the simulated results of the pressure solution model obviously increase with the growth of temperature

whereas the results are slightly influenced by the effective confining stress. Generally

the stress corrosion model reflects the growth of fracture contact area caused by compaction. Therefore

the coupled model can simulate the frictional healing well observed from Slide-Hold-Slide (SHS) tests. In conclusion

the stress corrosion model should be considered

when the pressure solution model is employed to simulate the frictional healing

especially under short hold times and at low temperatures when the domain effects are mechanical.