为什么会出现bad conversion ...floating.....的错误

侯健

最近做3DEC模拟进行赋参时，总是会出现bad conversion.....floating....这是命令流
;定义材料模型
change mat 1 cons 2
change mat 2 cons 2
change mat 3 cons 2
change mat 4 cons 2
change mat 5 cons 2
change mat 6 cons 2
;定义节理模型
change jmat 1 cons 1
change jmat 2 cons 1
change jmat 3 cons 1
change jmat 4 cons 1
change jmat 5 cons 1
change jmat 6 cons 1
;赋参数
def propmod
;煤
elas_m=0.475e9
u_m=0.01
bulk_m=elas_m/(3*(1-2*u_m))
shear_m=elas_m/(2*(1+u_m))
i1=out('bulk='+string(bulk_m))
j1=out('shear='+string(shear_m))
elas_fy=0.475e9
u_fy=0.01
bulk_fy=elas_fy/(3*(1-2*u_fy))
shear_fy=elas_fy/(2*(1+u_fy))
i2=out('bulk='+string(bulk_fy))
j2=out('shear='+string(shear_fy))
;砂质泥岩第一段
elas_sny1=0.375e9
u_sny1=0.15
bulk_sny1=elas_sny1/(3*(1-2*u_sny1))
shear_sny1=elas_sny1/(2*(1+u_sny1))
i3=out('bulk='+string(bulk_sny1))
j3=out('shear='+string(shear_sny1))

;砂质泥岩第二段
elas_sny2=0.05e9
u_sny2=0.1
bulk_sny2=elas_sny2/(3*(1-2*u_sny2))
shear_sny2=elas_sny2/(2*(1+u_sny2))
i4=out('bulk='+string(bulk_sny2))
j4=out('shear='+string(shear_sny2))
;砂质泥岩第三段
elas_sny3=3.0e9
u_sny3=0.25
bulk_sny3=elas_sny3/(3*(1-2*u_sny3))
shear_sny3=elas_sny3/(2*(1+u_sny3))
i4=out('bulk='+string(bulk_sny3))
j4=out('shear='+string(shear_sny3))

;粉砂岩
elas_fs=4.0e9
u_fs=0.2
bulk_fs=elas_fs/(3*(1-2*u_fs))
shear_fs=elas_fs/(2*(1+u_fs))
i5=out('bulk='+string(bulk_fs))
j5=out('shear='+string(shear_fs))
end
@propmod
  
prop mat 1 dens 4530 bulk bulk_sny3 shear shear_sny3 bfric 30 bcoh 1.8e6 bten 1.4e6 dil 0  
prop jmat 1 jkn=1.2e10 jks=1.2e10  jfric=30  jcoh=1.8e6 jten=0
prop mat 2 dens 4530 bulk bulk_sny2 shear shear_sny2 bfric 26 bcoh 1.2e6 bten 1.2e6 dil 0
prop jmat 2 jkn=1.2e10 jks=1.2e10  jfric=26  jcoh=1.2e6 jten=0
prop mat 3 dens 4530 bulk bulk_sny1 shear shear_sny1 bfric 22 bcoh 1.0e6 bten 1.0e6 dil 0
prop jmat 3 jkn=1.2e10 jks=1.2e10  jfric=22 jcoh=1.0e6 jten=0

prop mat 4 dens 4530 bulk bulk_m shear shear_m bfric 30 bcoh 1.8e6 bten 0.9e6 dil 0
prop jmat 4 jkn=1.2e10 jks=1.2e10  jfric=30  jcoh=1.8e6 jten=0
prop mat 5 dens 4530 bulk bulk_m shear shear_m bfric 45 bcoh 2.2e6 bten 2.0e6 dil 0
prop jmat 5 jkn=1.2e10 jks=1.2e10  jfric=45  jcoh=2.2e6 jten=0
prop mat 6 dens 4530 bulk bulk_fs shear shear_fs bfric 35 bcoh 2.6e6 bten 1.35e6 dil 0
prop jmat 6 jkn=1.2e10 jks=1.2e10  jfric=35  jcoh=2.6e6 jten=0