将节理法向循环加卸载试验所得的法向应力–闭合量曲线分别采用半对数、双曲线本构方程进行拟合,得到2种类型节理法向刚度表达式,并将其引入持续屈服节理模型中,同时编写节理法向循环加卸载程序,首次实现节理非线性本构循环加卸载数值模拟分析。研究结果表明:建立的节理法向非线性本构循环加卸载数值试验方法行之有效;由双曲线方程得到的节理法向刚度值较由半对数方程得到的法向刚度值略有不同,但2种方法得到的节理法向应力–闭合量曲线数值模拟结果均与室内试验十分接近,均可有效表征节理法向循环加卸载力学特征。研究成果可为后续节理岩体动力学性质研究,边坡动力响应分析以及边坡稳定性、安全性评价等奠定基础。 The normal stress-closure curves obtained by the cyclic loading and unloading test from the joint method were respectively fitted by semi-logarithmic and hyperbolic constitutive equations to obtain two types of joint normal stiffness expressions and introduced them into the continuous yield joint In the model, we also write the method of join and unload to the loop, and for the first time realize the numerical simulation analysis of the jointed and unloaded non-linear constitutive loop. The results show that the proposed method is effective in numerical simulation of nonlinear constitutive cycle loading and unloading. The normalized stiffness of the joint obtained by hyperbolic equation is slightly different from the normalized stiffness obtained by semi-logarithmic equation. However, the numerical simulation results of the joint normalized stress-closure curves obtained by the two methods are very close to those of the laboratory tests, which can effectively characterize the mechanical properties of the cyclic loading and unloading. The research results may lay the foundation for the subsequent study of the dynamic properties of jointed rock mass, the dynamic response analysis of slope and the slope stability and safety evaluation.