It is widely acknowledged that the natural rock mass is anisotropic and its failing type is also non-isotropic. An orthotropic elastic damaged model has been proposed in which the elastic deformation, the damaged deformation and irreversible deformation can be identified respectively. A second rank damage tensor is employed to characterize the induced damage and damage evolution related to the propagation conditions of microcracks. A specific form of the Gibbs free energy function is used to obtain the effective elastic stiffness and the limited scopes of damage parameters are suggested. The model's parameter determination is proposed by virtue of conventional tri-axial test. Then, the proposed model is developed to simulate the coupled hydraulic mechanical responses and traction behaviors in different loading paths of porous media.anisotropy, damage modeling, saturated rock, tensile failure Citation:Lu Y F, Wu X X, Shao J F. Anisotropic damage coupled modeling of saturated porous rock.Rock geo-materials exhibit a rich series of complex behaviors. Despite many researchers have studied the fundamental mechanical comportments of brittle rock for a long time, it still remains a challenge, particularly the anisotropic porous medium. The main consequences of saturated and drained tests can be summarized as follows: 1) deterioration of elastic properties and non-linearity of stressstrain relation; 2) induced material anisotropy and significant volume dilatation; 3) irreversible deformation, etc. The fundamental characteristics of injection and nondrained experiments can also be listed as follows: 1) reduction of initial maximum linear deformation and of failure peak stress; 2) more increment of radial deformation and volume dilatancy; 3) variation of volume compressive modulus; 4) strength change under different loading type, etc.The above fundamental features have to be taken into consideration in constructing constitutive models of geomaterial like rock. The first problem to be solved is the estimation of effective elastic properties of the geo-materials with micro-cracks. The most widely-used averaging schemes are the statistical method [16], the homogenization technology [17], the self-consistent method [16] and the differential scheme [16]. The statistical method, the self-consistent method and the homogenization technology are suitable for the materials of arbitrary distribution of micro-cracks with feeble concentration. The homogenization technology works for not only the material described in the above but also the material with a periodic distribution.