热力耦合PFC模型，岩热处理后的剪切试验

使用PFC2D5.0软件，模拟含结构面的花岗岩热处理后的剪切试验。试件尺寸是50mm×50mm，矿物成分有四种，按照一定比例分组，各种矿物热膨胀系数不同，试件内含一条结构面。 首先进行热处理，将试件加热到某一指定温度，然后进行冷却，冷却条件为：结构面放入20摄氏度水中冷却，其余面在20摄氏度空气中冷却，直到试件整体降温至20摄氏度；其次进行剪切试验，在试件垂直方向上施加压力，该压力在剪切过程中保持恒定，然后以一定速率进行剪切，剪切过程中只移动上部岩块，剪切速率恒定直到达到制定剪切位移时停止剪切。 希望在加热过程中能够观察到合理的热损伤裂纹以及温度变化，在剪切过程中能够观察到合理的剪切裂纹，监测峰值抗剪强度、剪切应力-应变曲线、剪胀角、剪胀量、剪切位移等相关参数。
详细内容见文献：Shear behavior of Beishan granite fractures after thermal treatment
希望能够复现文献中的模型