model configure thermal
zone cmodel assign mohr-coulomb
zone property density 2650 bulk 1.89e10 shear 1.04e10 friction 53 cohesion 7e6 tension 4.3e6
zone face apply velocity-normal 0 range group 'North'
zone face apply velocity-normal 0 range group 'Bottom'
zone face apply velocity-normal 0 range group 'West'
zone face apply stress-normal -2.69e7 range group 'East'
zone face apply stress-normal -5.3e7 range group 'South'
zone face apply stress-normal -2.8e7 range group 'Top'
zone thermal cmodel assign isotropic
zone thermal property conductivity=3.49
zone thermal property expansion=1.1e-6
zone thermal property specific-heat=920
zone initialize density 2650
zone thermal active on
zone gridpoint initialize temperature 293.15
;zone thermal cmodel assign null range group 'A'
zone gridpoint fix temperature 293.15 range group ...
'B' or 'V' or 'N1' or 'I1' or 'W' or 'S' or 'R' or 'Q' or 'L' or 'K'
zone gridpoint fix temperature 303.15 range group ...
'M' or 'N' or 'G' or 'F' or 'O' or 'H' or 'D' or 'E' or 'G1' or 'X'
zone gridpoint fix temperature 313.15 range group ...
'Y' or 'J1' or 'J' or 'I' or 'P' or 'A1' or 'Z' or 'L1' or 'K1' or 'B1'
zone gridpoint fix temperature 323.15 range group ...
'U' or 'T' or 'E1' or 'C1' or 'D1' or 'F1' or 'M1' or 'O1' or 'H1' or 'C'
上面这个代码，运行开挖，能够收敛， ratio 1e-3，能够达到。但是修改成下面这样，开挖就不收敛了，各位大佬知道怎么回事么。
model configure thermal
zone cmodel assign mohr-coulomb
zone property density 2650 bulk 1.89e10 shear 1.04e10 friction 53 cohesion 7e6 tension 4.3e6
zone face apply velocity-normal 0 range group 'North'
zone face apply velocity-normal 0 range group 'Bottom'
zone face apply velocity-normal 0 range group 'West'
zone face apply stress-normal -2.69e7 range group 'East'
zone face apply stress-normal -5.3e7 range group 'South'
zone face apply stress-normal -2.8e7 range group 'Top'
zone thermal cmodel assign isotropic
zone thermal property conductivity=3.49
zone thermal property expansion=1.1e-6
zone thermal property specific-heat=920
zone initialize density 2650
;zone thermal active on
zone gridpoint initialize temperature 293.15
;zone thermal cmodel assign null range group 'A'
zone gridpoint fix temperature 298.15 range group ...
'B' or 'V' or 'N1' or 'I1' or 'W' or 'S' or 'R' or 'Q' or 'L' or 'K' or 'M' or 'N' or 'G' or 'F' or 'O' or 'H' or 'D' or 'E' or 'G1' or 'X' or 'Y' or 'J1' or 'J' or 'I' or 'P' or 'A1' or 'Z' or 'L1' or 'K1' or 'B1' or 'U' or 'T' or 'E1' or 'C1' or 'D1' or 'F1' or 'M1' or 'O1' or 'H1' or 'C'
model thermal active on
