The crack propagation under external load is generally non-uniform in space

which causes derivative anisotropies of the material. In this context

rock materials can be discretized into a large number of randomly distributed material points connected by force bonds. Based on the directionality of the force bonds

the local elasticity tensor is discretized into a certain number of directional tensors

and the relationship between the force bond modulus and the macroscopic parameters is deduced theoretically. An anisotropic elastic damage constitutive model is developed by considering the fracture effect of the force bonds

In order to simulate the mechanical behavior of brittle to ductile transition of medium porosity rocks in conventional triaxial compression tests

a damage suppression function is introduced into force bond fracture effect. The rationality and effectiveness of the model are verified by simulating the conventional triaxial experiments of Tennessee marble and Indiana limestone and comparing with the experimental data.