Hydroxyapatite bioceramics is simulated by using finite element method （FEM）. The influences of porosity, hole shape, angle of crack and other parameters on the ceramics are analyzed. The results show that with the increase of the angle between crack and horizontal direction, the stress intensity factor KⅠ decreases gradually, but stress intensity factor KⅡ increases at first and then it decreases. The value of KⅡ reaches maximum when the angle between crack and horizontal direction is 45°. KⅠ and KⅡ rise with the increase of porosity, and they are almost the same for the circular and hexagonal holes. For elliptical holes, KⅠ and KⅡ reach maximum when the long axis of ellipse is perpendicular to the loading direction and they reach minimum when the same axis is parallel to the loading direction. Moreover, with the increase of the angle between the long axis and loading direction, KⅠ and KⅡ increase gradually.
