VASP计算弹性常数