ICRF heating in the Large Helical Device is studied applying two global simulation codes; a drift kinetic equation solver, GNET, and a wave field solver, TASK/WM. Characteristics of energetic ion distributions in the phase space are investigated changing the resonance heating position; i.e. the on-axis and off-axis heating cases. A clear peak of the energetic ion distribution can be seen in the off-axis heating case because of the stable orbit of resonant energetic ions. The simulation results are also compared with experimental results evaluating the count number of the neutral particle analyzer and a relatively good agreement is obtained.