To improve trapped field characteristics of a high temperature superconducting (HTS) bulk, a technique to implement artificial holes has been studied. The artificial holes, filled up with epoxy or metal, may provide better cooling channel and enhance mechanical strength of the HTS bulk. Although many useful researches based on experiments have been reported, a numerical approach is still limited because of several reasons that include: 1) highly non-linear electromagnetic properties of HTS; and 2) difficulty in modeling of randomly scattered "small" artificial holes. In this paper, a 2-D finite element method with iteration is adopted to analyze trapped field characteristics of HTS bulk with artificial holes. The validity of the calculation is verified by comparison between measurement and calculation of a trapped field in a $40{ imes}40;mm$ square and 3.1 mm thick HTS bulk having 16 artificial holes with diameter of 0.7 mm. The effects of sizes and array patterns of artificial holes on distribution of trapped field within HTS bulk are numerically investigated using suggested method.