生活中存在很多与热传导相关的实际问题,解析法所能解决的范围是极其有限的,随着计算机的发展,用数值方法求解复杂问题是一个大的趋势。热传导问题无论在理论上还是在数值方法上都已经有很长的研究历史,常用的数值方法有有限元法、有限差分法等。然而上述方法在求解无限域问题时往往是人为确定边界并施加边界条件,导致误差较大。不过由于上述方法本身的特点,求解无限域问题时往往人为确定边界,并施加边界条件。因此,误差较大。而边界元法具有降低维数、精度高等优点,且对于无限域问题在数学上是严密成立的。因此,本文选择基于边界元法来研究三维非稳态热传导的求解问题。本文基于非稳态热传导微分方程、基本解及δ-函数推导出域内及边界积分公式；然后,对时间和空间进行离散化,推导出三维非稳态热传导的离散化方程式。利用上述边界元算法、选用Fortran语言开发出了三维非稳态热传导问题的数值模拟系统；对开发的三维非稳态热传导边界元数值系统可靠性进行验证,将该系统得出的数值解同解析解进行比较分析,同时还与有限元数值解进行比较,结果表明基于边界元方法求解三维非稳态热传导问题的精度高、速度快、系统稳定可靠。在数值系统开发中解决了多项技术难题,并进行了关于时间积分的TIME INTEGRAL法与STEP WISE法的理论推导及系统开发。