基于SEM所获取的岩石非均质数字图像，提出了一种既能考虑预制和原生缺陷，又能近似反映材料非均匀性的三维细观结构有限元模型的方法，即单元质心对应法．应用该方法进行了热力耦合下岩石的热开裂及变形破坏数值模拟．结果表明，北山花岗岩热开裂临界温度大约为82℃，与实验值68～88℃基本吻合．低温时北山花岗岩热开裂以沿颗粒热开裂为主；而高温后既有沿颗粒热开裂，又有穿颗粒热开裂．温度的升高导致花岗岩由脆性向延性转变，破坏将由低温下以载荷裂纹为主的裂纹模式转为高温下以热裂纹为主的裂纹模式．对比实验结果发现，数值模拟得到的
