Abstract：In order to include the geometrical characteristics of jointed rock masses in its constitutive relations，a two-order tensor is adopted as damage variable. Joint connectivity at any direction of jointed rock masses is expressed as a continuum function of its normal vector and the damage tensor. Consequently，shear strength parameters of jointed rock masses，i.e.，friction factor and cohesion are weighted average of those of entire rock and entire fissures respectively according to joint connectivity，and implicit anisotropic yield criterion for jointed rock masses can be obtained through the Mohr-Coulomb yield condition for plane in every direction. In principal stresses space，the damage tensor is divided into two parts：the isotropic part and deviatoric part. According to rigor analytical solution to critical normal vector of the most disadvantage section for isotropic damage tensor，explicit yield criterion for isotropic joint rock masses is obtained，and an approximate analytical solution with one order precision to critical normal vector of the most disadvantage section for general damage tensor is obtained. The explicit anisotropic yield criterion for jointed rock masses is a form of quadratic equation of principal stresses in which its coefficients are functions of the damage tensor. The anisotropic yield criterion can be fully determined by the two-order damage tensor which has six independent variables，i.e.，the three principal values of damage tensor and the three Euler angles which represent the relations of principal vectors of damage tensor and those of principal stresses. In four examples，how the six variables affect the anisotropic properties of rock yield stresses in principal stresses space is analyzed by figuring out yield curve in deviatoric plane and meridian plane through the implicit and the explicit anisotropic strength criterion respectively. It is concluded that the explicit anisotropic yield criterion for jointed rock masses which is obtained by approximate analytical method has high precision.
[1] 田象燕，高尔根，白石羽. 岩石的应变率效应和各向异性的机理探讨[J]. 岩石力学与工程学报，2003，22(11)：1 789–1 792.(Tian Xiangyan，Gao Ergen，Bai Shiyu. Discussion about mechanism of strain-rate effect and anisotropy of saturated rocks[J]. Chinese Journal of Rock Mechanics and Engineering，2003，22(11)：1 789–1 792.(in Chinese))
[2] 余天堂. 岩土材料固有各向异性的模拟[J]. 岩石力学与工程学报，2004，23(10)：1 604–1 607.(Yu Tiantang. Modeling of inherent anisotropy for geotechnical material[J]. Chinese Journal of Rock Mechanics and Engineering，2004，23(10)：1 604–1 607.(in Chinese))
[3] Kawamoto T，Ichikawa Y，Kyoya T. Deformation and fracturing behaviour of discontinuous rock mass and damage mechanics theory[J]. International Journal for Numerical and Analytical Methods in Geomechanics，1988，12：1–30.
[4] 杨 强. 岩体损伤力学发展现状和面临的问题[A]. 见：中国岩石力学与工程学会第七次学术大会论文集[C]. 西安：中国科学技术出版社，2002. 46–50.(Yang Qiang. The state of art and problems in rock damage mechanics[A]. In：Symposium of the 7th Proceedings of the Institute of Chinese Rock Mechanics and Engineering[C]. Xi′an：China Science and Technology Press，2002. 46–50.(in Chinese))
[5] Pietruszczak S，Mroz Z. Formulation of anisotropic failure criteria incorporating a microstructure tensor[J]. Computers and Geotechnics，2000，26：105–112.
[6] Lydzba D，Pietruszczak S，Shao J F. On anisotropy of stratified rocks：homogenization and fabric tensor approach[J]. Computers and Geotechnics，2003，30：289–302
[7] Rabotnov Y N. On the equations of state for creep[J]. Progress in Applied Mechanics，1963，12：307–305.
[8] 余寿文，冯西桥. 损伤力学[M]. 北京：清华大学出版社，1997.(Yu Shouwen，Feng Xiqiao. Damage Mechanics[M]. Beijing：Tsinghua University Press，1997.(in Chinese))
[9] Swoboda G，Yang Q. An energy-based damage model of geomaterials—I. formulation and numerical results[J]. International Journal of Solids and Structures，1999，36：1 719–1 734.
[10] Lajtai E Z. Strength of discontinuous rocks in shear[J]. Geotechnique，1969，19(2)：218–233.
[11] 汪小刚，陈祖煜，孙文松. 应用蒙特卡罗法确定节理岩体的连通率和综合抗剪确定指标[J]. 岩石力学与工程学报，1992，11(4)：345–355.(Wang Xiaogang，Chen Zuyu，Sun Wensong. Determination of joint connectivity and comprehensive shear strength for jointed rock masses using Monte Carlo method[J]. Journal of Rock Mechanics and Engineering，1992，11(4)：345–355.(in Chinese))
[12] 杜景灿，汪小刚，陈祖煜. 结构面倾角对节理岩体的连通特性和综合抗剪屈服的影响[J]. 水利学报，2002，(5)：41–46.(Du Jingcan，Wang Xiaogang，Chen Zuyu. Determination of persistence and comprehensive shear strength of jointed rock mass based on waviness created by dip angle variations of discontinuities[J]. Chinese Journal of Hydraulic Engineering，2002，(5)：41–46.(in Chinese))
[13] 郑颖人，沈珠江，龚晓南. 岩土塑性力学原理[M]. 北京：中国建筑工业出版社，2002.(Zheng Yinren，Shen Zhujiang，Gong Xiaonan. The Principles of Geotechnical Plastic Mechanics[M]. Beijing：China Architecture and Building Press，2002.(in Chinese))
[14] Kanatani K. Distribution of directional data and fabric tensors[J]. Int. J. Engng. Sci.，1984，22：149–164.
[15] Yang Q，Li Z K，Tham L G . An explicit expression of second-order fabric-tensor dependent elastic compliance tensor[J]. Mechanics Research Communications，2001，28：255–260.
2005, Vol. 24 
