求解中子输运问题的方法主要有确定论方法和蒙特卡罗(Monte Carlo)方法。蒙特卡罗方法具有较真实模拟粒子历史和模拟复杂几何形状等优点。过去,蒙卡方法受制于计算机配置和计算能力的限制,通常是被用作确定论方法的补充。目前随着计算机技术的飞速发展,蒙卡方法在反应堆临界安全分析、输运计算及全堆芯精细模拟等领域得到了广泛应用。在反应堆计算中,堆芯物理参数与热工参数相互影响和制约,共同调节反应堆的运行状况,这一现象被称为反应堆热工反馈效应,是分析堆芯物理现象的基础。因此,研究蒙卡方法及蒙卡程序,并对其热工反馈效应进行研究具有重要意义,研究反应堆物理热工耦合计算有助于未来更真实模拟反应堆运行状态及揭示堆芯的物理特性。本文首先调研了蒙卡方法及程序并建立适合压水堆的热工计算模型。调研蒙卡方法及程序方面,主要对蒙卡程序的临界计算流程和Keff、功率分布的计算及统计方法进行了初步研究,选取了国核软件中心与清华大学联合开发的cos RMC蒙卡程序。建立热工计算模型方面,考虑压水堆堆型及适于蒙卡程序参数反馈的目的,选取单通道模型,通过建立燃料、包壳及冷却剂的轴向及横向传热换热简化模型,选取合适的热工参数及计算公式采用直接数值模拟方法加以研究。本文还对温度反馈的截面处理进行了相关研究,选取伪材料插值方法并制作不同温度下ACE格式截面库。首先对比分析了截面随温度反馈的处理方法,综合考虑计算精度及计算效率等因素选取伪材料插值方法;其次,使用SIGACE程序制作一定温度间隔的ACE格式截面库。选取栅元基准题,计算分析了使用不同温度截面库的燃料对有效增值因子Keff的影响,并对比了不同温度下235U截面的变化趋势,初步证明了所制作截面库的合理性和可行性。然后,对蒙卡程序cosRMC的临界功能进行测试验证并基于蒙卡程序cos RMC的临界功能开发热工反馈模块及耦合程序。首先,使用cos RMC对BEAVRS基准题进行建模计算,通过对临界本征值、控制棒价值及全堆功率分布的计算与分析,初步验证cos RMC的临界计算功能。其次,在已有蒙卡程序cos RMC的基础上使用C++进行编程二次开发,采用内耦合方式对源代码进行修改,添加了热工反馈计算模块及耦合程序,可实现热工参数(燃料、包壳、冷却剂温度及冷却剂密度)与物理参数(轴向功率分布)的反馈调节,完善了蒙卡程序的临界计算功能。最后,本文选取压水堆栅元计算模型,对耦合程序计算结果进行测试与验证分析。通过计算迭代收敛规律、Keff及功率分布结果验证热工反馈效应对蒙卡程序临界计算结果的影响,同时,通过对轴向燃料温度、包壳温度、冷却剂温度和密度等热工参数的结果对比分析,验证热工模块计算的准确性。本论文的研究为未来蒙卡程序更真实、准确的模拟反应堆运行状况和更广阔的应用空间打下基础,并为蒙卡程序热工反馈计算和分析提供借鉴与参考。