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

求助,FLAC程序,巷道锚杆支护 [复制链接]

上一主题 下一主题
离线csu瓜
关闭
个人中心可以申请新版勋章哦
立即申请知道了
 

发帖
320
土币
3633
威望
797
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2010-05-31
set log on eZ a:o1y  
set logfile mgan.log {qKxz9.y  
new eRbGZYrJ  
plot block group ^n#1<K[E  
plot show ]!:oYAm  
generate zone radcyl p0 = (0,0,0) p1 = (20,0,0) & qo+N,x9o  
p2 = (0,1,0) p3 = (0,0,20) p4 = (20,1,0) p5 = (0,1,20) & &m3.h!dq  
p6 = (20,0,20) p7 = (20,1,20) p8 = (2.4,0,0) p9 = (0,0,2.4) p10 = (2.4,1,0) & ;TQf5|R\K  
p11 = (0,1,2.4) size 10 5 10 10 dimension 2.4 2.4 2.4 2.4 ratio 1 1 1 1 fill group tunnel qZ@0]"h  
gen zone radtunnel p0 = (0,0,0) p1 = (0,0,-20) & zWw2V}U!  
p2 = (0,1,0) p3 = (20,0,0) & w)E@*h<Z  
size 10 5 10 10 dimension 2.2,2.4,2.2 fill group tunnel2 Bhp OXqg  
gen zone reflect normal (-1,0,0) origin =(0,0,0) 6Dws,_UAZ4  
group section1 range y 0 1 group tunnel any group tunnel2 any 0YH+B   
save model.sav tC8(XMVx  
pau O^LTD#}$a)  
;德鲁克-普拉格模型的材料参数 u{&B^s)k.  
model dr =9L$L|W  
prop bulk 1.41e9 shear 0.97e9  kshear 1.603e7 qvol 0.917 {-9jm%N  
;给定边界条件 iK;dU2h  
fix x range x -20.1 -19.9 +&tgJ07A  
fix x range x 19.9 20.1 4!$ M q;U  
fix z range z -20.1 -19.9 -7WW[ w  
fix y range y -0.1 0.1 HQ!Xj .y  
fix y range y 0.9 1.1 puSLqouTM  
;初始应力赋值 C2]Kc{4  
set grav 0 0 -10 B;Nl~Y|\  
ini density 2750 SEQ%'E5-'  
ini szz -17.6e6 grad 0 0 27500 T1B|w"In  
ini sxx -17.6e6 grad 0 0 27500 K.tlo^#^B[  
ini syy -17.6e6 grad 0 0 27500 "Z,q?Fc  
hist unbal kI*(V [i  
solve rh2LGuo4m  
save chushihua1_maogan.sav k'`m97B  
ini xvel 0 yvel 0 zvel 0 ,p{`pma  
ini xdis 0 ydis 0 zdis 0 .F&9.#>  
mod null range group section1 9L%I<5i  
solve MFJE6ei  
save kaiwa4.sav N\t1T(C|  
;第一排 >YW\~T  
sel cable id=1 begin 0 0 2.4 end 0 0 4.9 nseg 25     Auy".br'  
sel cable id=1 begin 0.92 0 2.2 end 1.88 0 4.53 nseg 25 '2J0>Bla  
sel cable id=1 begin 1.7 0 1.7 end 3.46 0 3.46 nseg 25 7>o .0  
sel cable id=1 begin 2.2 0 0.92 end 4.53 0 1.88 nseg 25 y#ON|c /  
sel cable id=1 begin 2.4 0 0 end 4.9 0 0 nseg 25 9D@$i<D:  
sel cable id=1 begin 2.4 0 -0.8 end 4.9 0 -0.8 nseg 25 PDx)S7+w[  
sel cable id=1 begin 2.4 0 -1.6 end 4.9 0 -1.6 nseg 25 -9P2`XQ^  
sel cable id=1 begin -0.92 0 2.2 end -1.88 0 4.53 nseg 25 ,Y_{L|:w  
sel cable id=1 begin -1.7 0 1.7 end -3.46 0 3.46 nseg 25 sfp,Lq`  
sel cable id=1 begin -2.2 0 0.92 end -4.53 0 1.87 nseg 25 9z m|Lbj  
sel cable id=1 begin -2.4 0 0 end -4.9 0 0 nseg 25 [{[N(g&d  
sel cable id=1 begin -2.4 0 -0.8 end -4.9 0 -0.8 nseg 25 k0?ZYeHC  
sel cable id=1 begin -2.4 0 -1.6 end -4.9 0 -1.6 nseg 25 Ue5O9;y]u  
sel cable pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. & QrD o|GtE  
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id 1 t$& Qv)  
solve nR \'[~+  
save maogan4.save   )!9Ifk0KH  
Itm8b4e9;  
&0N<ofYX  
请大家帮忙看看,我这个程序为什么运行了,支护没效果,谢谢!!
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      https://beian.mps.gov.cn/ 粤公网安备 44010602012919号 广州半山岩土网络科技有限公司 粤ICP备2024274469号

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