Substituent effects in Nd2Fe14B-based alloys interpreted on a crystal field model