基于CT图像生成晶粒模型步骤（在CT图像中，浅色区域表示高密度矿物，灰色区域表示低密度矿物；云母矿物显示为白色，长石显示为深灰色，石英矿物显示为浅灰色[24]）

下载：全尺寸图片幻灯片

基于CT图像生成晶粒模型步骤（在CT图像中，浅色区域表示高密度矿物，灰色区域表示低密度矿物；云母矿物显示为白色，长石显示为深灰色，石英矿物显示为浅灰色^[24]）

性质⁽¹⁾		ρ_m/ (kg·m^-3)	ν_m	E_m/ GPa	P_ct/ %	性质接触	φ_m/ (°)	c_m/ MPa	σ_tm/ MPa	k_n/ (GPa·m^-1)
晶体	石英	2650	0.16	37.5	44	Q-Q⁽²⁾ M-M⁽²⁾	30 15	60 40	13 7	119802 51376
	长石	2620	0.19	22.5	50	F-F⁽²⁾ Q-F⁽²⁾	25 27.5	50 55	11 12	74105 119802
	云母	3050	0.22	15	6	Q-M⁽²⁾ F-M⁽²⁾	22.5 20	50 45	10 9	119802 74105
注：（1）微观密度( ${\rho _{\text{m}}}$ )、微观杨氏模量( ${E_{\text{m}}}$ )、微观泊松比( ${v_{\text{m}}}$ )、微观内聚力( ${c_{\text{m}}}$ )和微观摩擦角( ${\varphi _{\text{m}}}$ )，数据源自文献[]；（2）石英、云母、长石晶粒内的接触分别表示为Q-Q、M-M、F-F，晶粒间的接触为Q-F、Q-M、F-M，其强度为相邻矿物晶体的平均强度；接触法向刚度 ${k_{\text{n}}} = 5(k + {\text{4}}G/3)/\Delta {z_{{\text{min}}}}$ ，其中K和G分别是体积模量和剪切模量， $\Delta {z_{{\text{min}}}}$ 是法线方向上相邻区域的最小宽度^[13]；为求简化，晶粒接触切向刚度（ ${k_{\text{s}}}$ ）均设置为 ${k_{\text{n}}}/{\text{4}}$ 。

性质⁽¹⁾		α_t⁽²⁾ / (10^-6·K^-1)	k/ (W·m^-1·K)	C_p/ (J·kg^-1·K)
晶体	石英	13.40	7.70	750.0
	长石	3.20	2.38	640.0
	云母	12.81	2.01	769.5
注：（1）热性质 ${\alpha _{\text{t}}}({\text{CLTE}})$ ，k（导热系数）， ${C}_{\text{p}}(比热)$ ，数据源自文献[]，其中k和 ${C_{\text{p}}}$ 均设置为恒定值；（2）石英在冷却期间的CLTE遵循图 1（b）中的温度依赖性曲线^{[12, 21-22]}。因为云母晶体仅占6%，并且没有经历可逆的相变，为求简化，其在冷却过程中的线性热膨胀系数被设定为加热过程中数值的1/2。