In Fracture phenomenon of orthotropic materials, generally in crack tip vicinity, an area called damaged zone appears that in quasi-brittle materials, is known as fracture process zone. This area contains a multitude of micro cracks which due to various reasons, failure analysis and fracture process of these materials has been difficult. Determination of Mechanical properties in this region can help to predict the value or even the direction of crack growth in orthotropic materials. So far, several models have been proposed to determine the mechanical properties of this region, but due to the immense complexity of this region, the results have not been expressed the behavior of this region properly. Moreover, the existence methods have not been verified with new experimental and numerical data, yet. the present paper, it was attempted to present a new numerical model based on viscoelastic theory (the linear springs and damping) and, according to Perony series, the mechanical properties of the damaged zone simulated. In new approach, that is based on experimenal and FEM results, the mechanical behavior of damaged zone can accurately simulated.