论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 5436阅读
  • 3回复

FLAC3D5.0新手求教 [复制链接]

上一主题 下一主题
离线yang1192226
 

发帖
4
土币
34
威望
3
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2017-02-18
— 本帖被 yanzi7023豆豆 执行合并操作(2017-02-20) —

运行这段命令的时候会出现如下的提示,请问各位大神是什么情况,该怎么解决啊?
new
config fluid
gen zone radcyl p0 0 0 0 p1 5 0 0 p2 0 20 0 p3 0 0 5 &
dim 3 3 3 3 size 2 10 4 2 ratio 1.0 1.0 1.0 1.2 fill &
group tunnel
gen z b p0 0 0 5 p1 5 0 5 p2 0 20 5 p3 0 0 15 &
rat 1 1 1 size 2 10 4
gen z b p0 5 0 0 p1 25 0 0 p2 5 20 0 p3 5 0 15 &
rat 1 1 1 size 8 10 6
group soil range group tunnel not
gen  zone reflect normal 0 0 -1
fix x range x -.1 .1
fix x range x 24.9 25.1
fix y range y -.1 .1
fix y range y 19.9 20.1
fix z range z -14.9 -15.1
;cam clay
model cam-clay
prop shear 150000 bulk_bound 20e6
prop mm 0.73 lambda 0.0938 kappa 0.0234
prop mpc 0.395e6 mp1 1.0 mv_l 3.32
ini dens 1270
;fluid para
model fl_iso
prop perm 2.66e-13
ini fdens 1000
ini sat 1.0
set grav 0 0 -10
;ini stress and pp
ini szz -190500 grad 0 0 1.27e4
ini szz add -0.795e5 grad 0 0 0.53e4
ini sxx -128509.5 grad 0 0 0.857e4
ini sxx add -0.795e5 grad 0 0 0.53e4
ini syy -128509.5 grad 0 0 0.857e4
ini syy add -0.795e5 grad 0 0 0.53e4
ini pp 1.5e5 grad 0 0 -10000
fix pp 0 range z 14.9 15.1
;xiuzeng jianqiaomoxing
def camclay_ini_p
  pnt = zone_head
  loop while pnt # null
     OCR=1.2
     s1=-z_sxx(pnt)
     s2=-z_syy(pnt)
     s3=-z_szz(pnt)
     p0=(s1+s2+s3)/3.0-z_pp(pnt)
     z_prop(pnt,'cam_cp')=p0
     q0=sqrt(((s1-s2)*(s1-s2)+(s2-s3)*(s2-s3)+(s3-s1)*(s3-s1))*0.5)
     temp1=q0/(z_prop(pnt,'mm')*p0)
     pc=p0*(1.0+temp1*temp1)*OCR
     z_prop(pnt,'mpc')=pc
     pnt=z_next(pnt)
   end_loop
end
@camclay_ini_p
;chushi pingheng
set fl off
set mech on
solve
ini xdisp 0 ydisp 0 zdisp 0
save ini_equivalence
hist fltime
;kaiwamian parameter
def sup_stress
  lumda=0.9
  o_press=-208009.55
  o_grad=1.39e4
  s_press=o_press*lumda
  s_grad=o_grad*lumda
end
@sup_stress
;first huan kaiwa parameter
def excate_step1
n=1
a1=2*n-2
a2=2*n
b1=2*n-0.01
b2=2*n+0.01
n1=16*n-1
n2=16*n
t=4*3600*n
end
@excate_step1
model null range cyl end1 0 a1 0 end2 0 a2 0 rad 3.0
离线yang1192226

发帖
4
土币
34
威望
3
原创币
0
只看该作者 1楼 发表于: 2017-02-18
FLAC3D5.0求指导
运行这段命令的时候会出现如下提示,请问各位大侠是什么情况呀,该怎么解决呢?
new
config fluid
gen zone radcyl p0 0 0 0 p1 5 0 0 p2 0 20 0 p3 0 0 5 &
dim 3 3 3 3 size 2 10 4 2 ratio 1.0 1.0 1.0 1.2 fill &
group tunnel
gen z b p0 0 0 5 p1 5 0 5 p2 0 20 5 p3 0 0 15 &
rat 1 1 1 size 2 10 4
gen z b p0 5 0 0 p1 25 0 0 p2 5 20 0 p3 5 0 15 &
rat 1 1 1 size 8 10 6
group soil range group tunnel not
gen  zone reflect normal 0 0 -1
fix x range x -.1 .1
fix x range x 24.9 25.1
fix y range y -.1 .1
fix y range y 19.9 20.1
fix z range z -14.9 -15.1
;cam clay
model cam-clay
prop shear 150000 bulk_bound 20e6
prop mm 0.73 lambda 0.0938 kappa 0.0234
prop mpc 0.395e6 mp1 1.0 mv_l 3.32
ini dens 1270
;fluid para
model fl_iso
prop perm 2.66e-13
ini fdens 1000
ini sat 1.0
set grav 0 0 -10
;ini stress and pp
ini szz -190500 grad 0 0 1.27e4
ini szz add -0.795e5 grad 0 0 0.53e4
ini sxx -128509.5 grad 0 0 0.857e4
ini sxx add -0.795e5 grad 0 0 0.53e4
ini syy -128509.5 grad 0 0 0.857e4
ini syy add -0.795e5 grad 0 0 0.53e4
ini pp 1.5e5 grad 0 0 -10000
fix pp 0 range z 14.9 15.1
;xiuzeng jianqiaomoxing
def camclay_ini_p
  pnt = zone_head
  loop while pnt # null
     OCR=1.2
     s1=-z_sxx(pnt)
     s2=-z_syy(pnt)
     s3=-z_szz(pnt)
     p0=(s1+s2+s3)/3.0-z_pp(pnt)
     z_prop(pnt,'cam_cp')=p0
     q0=sqrt(((s1-s2)*(s1-s2)+(s2-s3)*(s2-s3)+(s3-s1)*(s3-s1))*0.5)
     temp1=q0/(z_prop(pnt,'mm')*p0)
     pc=p0*(1.0+temp1*temp1)*OCR
     z_prop(pnt,'mpc')=pc
     pnt=z_next(pnt)
   end_loop
end
@camclay_ini_p
;chushi pingheng
set fl off
set mech on
solve
ini xdisp 0 ydisp 0 zdisp 0
save ini_equivalence
hist fltime
;kaiwamian parameter
def sup_stress
  lumda=0.9
  o_press=-208009.55
  o_grad=1.39e4
  s_press=o_press*lumda
  s_grad=o_grad*lumda
end
@sup_stress
;first huan kaiwa parameter
def excate_step1
n=1
a1=2*n-2
a2=2*n
b1=2*n-0.01
b2=2*n+0.01
n1=16*n-1
n2=16*n
t=4*3600*n
end
@excate_step1
model null range cyl end1 0 a1 0 end2 0 a2 0 rad 3.0
*** Bad conversion of parameter number 6 (0).
Expected tokens:
     3D vector of floating point values
离线yang1192226

发帖
4
土币
34
威望
3
原创币
0
只看该作者 2楼 发表于: 2017-02-18
有没有大神呀~~
离线北伊

发帖
174
土币
16904
威望
2
原创币
0
只看该作者 3楼 发表于: 2017-02-21
建立三维网格的时候没有严格遵守右手系规则。
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      浙公网安备 33010602003799号 浙ICP备14021682号-1

      工业和信息化部备案管理系统网站