- 积分
- 1
- 注册时间
- 2005-5-16
- 仿真币
-
- 最后登录
- 1970-1-1
|
;第三步:支护
;set large
;锚杆支护
def cct1;定义常数
mchdu=2.5;锚杆长度
mjju=0.8;0.5,1.0锚杆间距
yuch=0.2;0.4(多加根锚杆),0.4从半圆形拱与直墙相接处出发的圆弧长
jid=cos(45*degrad)
jd=sin(45*degrad)
end
cct1
;定义锚固起止坐标
def star
array a_x(7),a_z(7)
a_x(1)=2.5
a_z(1)=-1.4
a_x(2)=2.5
a_z(2)=-0.6
loop i (3,7)
a_x(i)=2.5*cos((yuch+(i-3)*mjju)/2.5)
a_z(i)=2.5*sin((yuch+(i-3)*mjju)/2.5)
end_loop
end
;
def rstar
array b_x(7),b_z(7)
b_x(1)=2.5+mchdu*jid
b_z(1)=-1.4-mchdu*jd
b_x(2)=2.5+mchdu
b_z(2)=-0.6
loop i (3,7)
b_x(i)=(2.5+mchdu)*cos((yuch+(i-3)*mjju)/2.5)
b_z(i)=(2.5+mchdu)*sin((yuch+(i-3)*mjju)/2.5)
end_loop
end
;
star
rstar
;
def cct2;定义常数
zhjing=0.02
jmj = pi*zhjing*zhjing/4
pper=pi*zhjing
gxj=pi*power(zhjing,4)/64
jgxj=2*gxj
end
cct2
;加锚杆
def pile_setup
loop j(1,3);第j排,每排7根
loop i(1,7)
pnum=i+(j-1)*7
yj=1.7+(j-1)*0.8;根据锚杆间距取值
command
sel pile id pnum begin a_x(i) yj a_z(i) end b_x(i) yj b_z(i) nseg 10
;这里会出错,cos(),sin()不能这样用吗
sel pile property rockbolt on
;sel pile property cs_ng on
sel pile pro emod 220e9 nu 0.3 tyield 8e8 xcarea jmj &
xciy 0 xciz 0 xcj 0 &
perimeter pper &
ydirection 0,1,0 &
CS_sK=2e10 CS_sCoh=18.8e5 CS_sFric=10.0 &
CS_nK=2e10 CS_nCoh=18.8e5 CS_nFric= 10.0
;gr_per 0.009425 gr_coh 18.8e5 yten 141.37 gr_k 1.107e9
end_command
;加载
loop k(1,10)
command
sel node apply force 50e3*(11-k)/10 0 0 cid=k
end_command
end_loop
end_loop
end_loop
end |
|