pfc2d 5.0热力破岩

首先我有单种成分的热力破岩的代码，我的需求是帮我把成样和胶结改成近似花岗岩的组成（三种成分按照一定比例随机分布，或者网上又通过图片扫描到cad然后再导入到这个pfc中的也行）就是希望可以满足三种主要成分建的初始模型，然后我有给胶结模型设置一定初始温度的代码（我只有热传递的代码），但是我的实物实验是等离子破岩它有一定的力但是主要还是热破岩，如果可以给我的热代码加一个单独的力等效成等离子破岩（力加热）是最好的希望大神们踊跃接单咨询。
在我的需求图片中前两个是我做的岩石模型都是只有单一成分，需要升级，第三个是热传递希望可以搭配一个力来等效成高温粒子冲击加热共同破岩的效果。最后就是代码给我写出来需要一定的注释，并且最后要指导我一下怎么将仿真模拟的结果导出（我的变量是初始温度和加热时间，结果是裂纹数量）