In the present study, we propose Monte Carlo simulation that takes into consideration the interface phenomena between liquid crystal and substrate. We use rigid model molecules of liquid crystal and substrate. Interface is generated using potential field that induces decomposition of molecules. We use hard spherocylinders as model liquid crystal molecules. Substrate is modeled as region composed of shorter spherocylinders. Our results show that there is a case in which nematic order is reinforced in the vicinity of rubbed substrate.