最近在做一门课的课程设计想试着用FLAc3d做然后运算的时候就会如下提示
FLac3d 3.0-251里提示math error：invalid number
FLac3d 3.0-261里提示如下
math error：
invalid
file：.\Cycle.cpp
line：466
其实再solve一次的话就可以运算了但是运算出来的结果和理论分析的结果差异很大，朋友们帮帮忙~~
好吧，大体说下我做的一个东西吧，模拟的是一个盾构隧道的衬砌结构
new
title
盾构支护
;
;隧道所在地层
gen zone radc p0 0 0 0 p1 0 0 -6 p2 0 8 0 p3 6 0 0 p4 0 8 -6 p5 6 8 0 p6 6 0 -6 p7 6 8 -6 p8 0 0 -3.1 p9 3.1 0 0 p10 0 8 -3.1 p11 3.1 8 0 size 5 8 8 3 rat 1 1 1 1.3 fill group tunnel
gen zone radc p0 0 0 0 p1 6 0 0 p2 0 8 0 p3 0 0 4.28 p4 6 8 0 p5 0 8 4.28 p6 6 0 4.28 p7 6 8 4.28 p8 3.1 0 0 p9 0 0 3.1 p10 3.1 8 0 p11 0 8 3.1 size 5 8 8 3 rat 1 1 1 1.3 fill group tunnel
;隧道所在地层两边土体
gen zone bri p0 6 0 0 p1 add 0 0 -6 p2 add 0 8 0 p3 add 3 0 0 p4 add 0 8 -6 p5 add 3 8 0 p6 add 3 0 -6 p7 add 3 8 -6 size 4 8 3 rat 1 1 1
gen zone bri p0 6 0 0 p1 add 3 0 0 p2 add 0 8 0 p3 add 0 0 4.28 p4 add 3 8 0 p5 add 0 8 4.28 p6 add 3 0 4.28 p7 add 3 8 4.28 size 3 8 4 rat 1 1 1
;对隧道所在地层分组
group Stratum4 range z 1.28 4.28 group tunnel not
group Stratum5 range z -6 1.28 group tunnel not
;隧道上覆土层
gen zone bri p0 0 0 4.28 p1 add 6 0 0 p2 add 0 8 0 p3 add 0 0 3.5 p4 add 6 8 0 p5 add 0 8 3.5 p6 add 6 0 3.5 p7 add 6 8 3.5 size 4 8 5 rat 1 1 1
gen zone bri p0 0 0 7.78 p1 add 6 0 0 p2 add 0 8 0 p3 add 0 0 1 p4 add 6 8 0 p5 add 0 8 1 p6 add 6 0 1 p7 add 6 8 1 size 4 8 2 rat 1 1 1
gen zone bri p0 0 0 8.78 p1 add 6 0 0 p2 add 0 8 0 p3 add 0 0 1.5 p4 add 6 8 0 p5 add 0 8 1.5 p6 add 6 0 1.5 p7 add 6 8 1.5 size 4 8 3 rat 1 1 1
gen zone bri p0 6 0 4.28 p1 add 3 0 0 p2 add 0 8 0 p3 add 0 0 3.5 p4 add 3 8 0 p5 add 0 8 3.5 p6 add 3 0 3.5 p7 add 3 8 3.5 size 3 8 5 rat 1 1 1
gen zone bri p0 6 0 7.78 p1 add 3 0 0 p2 add 0 8 0 p3 add 0 0 1 p4 add 3 8 0 p5 add 0 8 1 p6 add 3 0 1 p7 add 3 8 1 size 3 8 2 rat 1 1 1
gen zone bri p0 6 0 8.78 p1 add 3 0 0 p2 add 0 8 0 p3 add 0 0 1.5 p4 add 3 8 0 p5 add 0 8 1.5 p6 add 3 0 1.5 p7 add 3 8 1.5 size 3 8 3 rat 1 1 1
;分组
group Stratum3 ran z 4.28 7.78
group Stratum2 ran z 7.78 8.78
group Stratum1 ran z 8.78 10.28
;镜像
gen zone ref o 0 0 0 n -1 0 0
;设置地层参数
model mohr ran y 0 8
prop b 6.25e6 s 2.88e6 coh 19.5e3 fri 26.8 ran group Stratum1
prop b 5.44e6 s 2.51e6 coh 30.25e3 fri 21.45 ran group Stratum2
prop b 8.79e6 s 4.06e6 coh 12.5e3 fri 31.18 ran group Stratum3
prop b 5.76e6 s 2.36e6 coh 12.2e3 fri 8.9 ran z 1.28 4.28
prop b 11.88e6 s 1.74e6 coh 12.1e3 fri 7.2 ran z -6 1.28
ini density 1800 ran group Stratum1
ini dens 1910 ran group Stratum2
ini dens 1870 ran group Stratum3
ini dens 1800 ran group z 1.28 4.28
ini dens 1710 ran group z -6 1.28
;边界条件
fix z ran z -6.001 -5.999
fix x ran x -9.001 -8.999
fix x ran x 8.999 9.001
fix y ran y -.001 .001
fix y ran y 7.999 8.001
set grav 0 0 -10
;因为一开始做得太厚了不考虑开挖似乎有点多余就删掉些
del range y 3 8
fix y ran y 2.999 3.001
;施加地表的超载
app szz -20e3 ran z 10.28
;地下水的设置，这里抄了书上的内容好像是考虑渗流的，一开始想做不考虑渗流但是下面的solve收敛不了
config fluid
model fl_iso range z 9.78 10.28 not
ini fdens 1000
ini pp 9.78e4 grad 0 0 -10e3 ran z 9.78 10.28 not
set mech on fluid off
ini fmod 0
hist id=1 unbal
solve
sav liner_without_tun_20151218.sav
restore liner_without_tun_20151218.sav
ini zdisp 0.0 ydisp 0.0 zdisp 0.0
;开挖
model null range group tunnel
;设置接触面，主要目的是想拿来给管片隔水
interface 1 face ran cyl end1 0 0 0 end2 0 3 0 radius 3.2 cyl end1 0 0 0 end2 0 3 0 radius 3.0 not cyl end1 0 3 0 end2 0 0 0 radius 3.2 cyl end1 0 3 0 end2 0 0 0 radius 3.0 not
app nstress 9.78e4 grad 0 0 -10e3 ran cyl end1 0 0 0 end2 0 3 0 radius 3.2 cyl end1 0 0 0 end2 0 3 0 radius 3.0 not cyl end1 0 3 0 end2 0 0 0 radius 3.2 cyl end1 0 3 0 end2 0 0 0 radius 3.0 not
int 1 prop kn 7.4e10 ks 7.4e10 f 0 c 1e20
int 1 perm off
plot create the_system
plot set rot 0 0 50
plot add bl group yellow
plot add axes black
;铺设衬砌管片
sel liner id=1 ran cyl end1 0 0 0 end2 0 3 0 rad 3.1
sel node local xdir 1 0 0 ydir 0 -1 0 ran x -0.1 0.1
sel node fix lsys ran x -0.1 0.1
sel node fix x yr zr range x -0.1 0.1
sel node local xdir 0 0 -1 ydir 0 -1 0 ran z -0.1 0.1
sel node fix lsys ran z -0.1 0.1
sel node fix x yr zr ran z -0.1 0.1
sel node fix y xr zr ran y -0.1 0.1
fix y ran y 2.999 3.001
fix y ran y -.001 .001
;衬砌参数
sel liner id=1 prop iso=(3.45e7, 0.25) thick=0.35
sel liner id=1 prop cs_nk=7.4e10 cs_sk=7.4e10 cs_scoh=1e20
sel liner id=1 prop den=2500
plot add sel geom scale 0.005 shrink=0.0 black black
;plot mov 3 2
plot add sel recover sres mx surfx 0 1 0
plot show
plot add sel geom
solve
额，希望大家可以帮我看下里面有哪些问题，第一次做很没谱，和学长交流了下改了改问题依旧存在。。。谢谢大家了！！
亲测可以收敛，３００００步左右吧。。。
谢谢！！！！！！！！
