本文为了改进31P的全套中子核数据，据根据最新的和更加可靠的实验数据，对0.1-20MeV能区进行了理论模型计算。通过查询EXFOR库和CINDA索引库获取最新的实验数据。充分考虑了实验数据及理论参数的系统性，尤其对实验数据少或实验数据偏差大的数据进行了研究和取舍。根据最新的和更加可靠的实验数据，通过光学势调参程序APOM94，给出了一套能再现实验数据的光学势参数。核反应的直接部分用扭曲波波恩近似的DWUCK4程序进行计算，用UNF程序进行光参及相关参数调整和全套中子核数据的理论计算，计算过程中需要根据收集到的实验数据进行参数的调节，得出可靠的全套中子数据，包括各个反应道的截面，角分布，双微分截面等。计算中用到APOM94、DWUCK4、UNF三个主要程序，及其他辅助程序。模型计算结果与实验值符合的很好，各个反应道截面，角分布数据给出了与ENDF/B-6和JENDL3.3比较的图。与ENDF/B-6和JENDL3.3相比，UNF程序的n+31p模型计算有许多改进，结果更合理可信。全套n+31p的数据以ENDF/B格式给出。