本文采用Monte-Carlo数值模拟方法，研究无序钉扎和热涨落驱动的磁通格子的熔化相变。计算了磁通线格子在Z轴方向的关联函数(Cz)，用磁通线格子在Z轴方向的纠缠和退耦合来定义磁通格子的熔化相变。 首先，计算了磁通格子的3D-2D相变，定义Cz=0.45为相变点。发现在低温低磁场下1＞Cz＞0.45时，磁通线排列成稍微弯曲的没有位错的弹性格子，定义系统处于3D有长程序的布拉格玻璃相(BG)；随温度或磁场的增加Cz＜0.45时，磁通线格子纠缠，退耦合为2D相。这种3D-2D相变也可称为熔化相变到无序的磁通液体相(VL)或磁通玻璃相(VG)。调节参数使符合Bi2Sr2CaCu2O8材料数据，得到了与等离子共振实验相似的磁通线格子的3D-2D相图。得到了由磁场无序作用产生的BG-VG固-固熔化相变线与温度无关。而早先的理论工作得到的固-固熔化线随温度升高而减小，这与其他理论工作不同之处。解释其原因是由于在模型中，位于同一根磁通线而在不同平面上的磁通饼之间有很强的吸引互作用，使得在低温和中等温度区域内，由温度涨落引起的无序效应影响很小，所以计算得到的固-固熔化相变线平行于温度轴。计算结果表明此模型是合理的，而且具有一定的优越性。 然后，计算得到平面内和平面间磁通饼互作用力有异常的温度关系。定义在1.0＞Cz＞0.65时，系统处在具有准长程序的BG相。当Cz＜0.6并直线下降到0，平面间磁通无关联，定义系统处于无序的VL相。当Cz＜0.65，并保持在0.65-0.6之间，定义系统处于无序固体VG相，计算得到了无序钉扎强度—温度相图。我们发现在所得相图中，磁通格子的BG到VL一级熔化相变(FOT)线并没有在临界点TCP终止，而是继续延伸至低温，同时在临界温度点TCP附近区域，BG到VG一级相变线形成一个向上的突起。当温度大于TCP时，熔化线有负的斜率，磁通格子随温度的升高，发生热涨落驱动的BG-VL的一级熔化相变；当温度小于TCP时，熔化线斜率为正，磁通格子随温度的降低会发生由无序驱动的从有序BG相到非晶的无序VG相的转变。因此在中等温度区域，随温度的降低也会从有序相进入无序的非晶相的相变，表现出了异常的反向熔化行为，该现象与Avraham等人在最近实验上观测到的反向熔化行为一致。