我做的隧道全断面开挖、梅花型锚杆大家可以看看，希望大家批评建议啊！
我做的隧道全断面开挖、梅花型锚杆大家可以看看，希望大家批评建议啊！[ 本帖最后由 benjackxu 于 2007-1-5 22:05 编辑 ] 虽然没看得出程序哪里除了毛病
但是应力释放好像没有控制啊 本来就没毛病啊！！ 其他方面还的控制啊~！！ 不知道 有什么好的建议 可以发表下！我好借鉴下！！ 斑竹可以看看啊 我没做过你这个，等有时间学学，先加分鼓励一下。 定义超前锚秆
def fen_jiao_chao
jiao=float(120.0/32.0)
end
fen_jiao_chao
def cab_parm_chao
array x_b(1,33),z_b(1,33),x_2b(1,33),z_2b(1,33)
yd_mao=0
loop i(1,33)
x_b(1,i)=6.19*cos((30+jiao*(i-1))*degrad)
z_b(1,i)=6.19*sin((30+jiao*(i-1))*degrad)
x_2b(1,i)=(6.19+3.5*sin(14*degrad))*cos((30+jiao*(i-1))*degrad)
z_2b(1,i)=(6.19+3.5*sin(14*degrad))*sin((30+jiao*(i-1))*degrad)
end_loop
end
def inip_chao
x=x_b(1,i)
z=z_b(1,i)
x_2=x_2b(1,i)
z_2=z_2b(1,i)
end
cab_parm_chao
def ini_chao_mao
y2=yd_mao+3.5*cos(14)
loop i(1,33)
inip_chao
command
sel cable id i begin x,yd_mao,z endx_2,y2,z_2nseg 5
sel cable pro emod 45e9 xcarea 1.57e-3 gr_per 1. &
yten 25e4 gr_k 17.5e6 gr_c 20e4range id i
step 1
end_command
end_loop
end
;定义初期支护锚秆
def fen_jiao_chu
h=9.0
;9=180/20
end
fen_jiao_chu
def cab_parm_chu_1
array x2_b(1,54),z2_b(1,54),x2_2b(1,54),z2_2b(1,54)
loop n(34,54)
x2_b(1,n)=6.19*cos(h*(n-34)*degrad)
z2_b(1,n)=6.19*sin(h*(n-34)*degrad)
x2_2b(1,n)=9.19*cos(h*(n-34)*degrad)
z2_2b(1,n)=9.19*sin(h*(n-34)*degrad)
end_loop
end
def cab_parm_chu_2
loop n(34,54)
x2_b(1,n)=6.19*cos(h*(n-35)*degrad)
z2_b(1,n)=6.19*sin(h*(n-35)*degrad)
x2_2b(1,n)=9.19*cos(h*(n-35)*degrad)
z2_2b(1,n)=9.19*sin(h*(n-35)*degrad)
end_loop
end
def number
a=u/2.0
b=int(u/2)
if a = b then
cab_parm_chu_1
else
cab_parm_chu_2
endif
end
def inip_chu
x_d=x2_b(1,n)
z_d=z2_b(1,n)
x_d_2=x2_2b(1,n)
z_d_2=z2_2b(1,n)
end
number
def ini_chu_mao
loop n(34,54)
inip_chu
command
sel cable id n begin x_d,chumaod_y,z_d end x_d_2,chumaod_y, z_d_2 nseg 10;有问题
sel cable pro emod 45e9 xcarea 1.57e-3 gr_per 1. &
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id n
step 1
end_command
end_loop
end
;安装超前锚杆-----开挖------安装初期支护锚杆
def excavate
cut=0
chumaod_y=1.2
u=0
loop whilecut <= 48.0
cut1 =cut +2
cut2=cut-2
;安装超前锚杆
ini_chao_mao
;开挖
command
m n range group nsoil y cut cut1
step 1
end_command
yd_mao=yd_mao+2
;安装初期支护锚杆和喷混凝土
ini_chu_mao
p=cut1-1.2
if chumaod_y <= p
chumaod_y =chumaod_y+1.2
u=u+1
number
loop n(34,54)
inip_chu
command
sel cable id n begin x_d,chumaod_y,z_d end x_d_2,chumaod_y, z_d_2 nseg 10
sel cable pro emod 45e9 xcarea 1.57e-3 gr_per 1. &
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id n
step 1
end_command
end_loop
end_if
chumaod_y =chumaod_y+1.2
u=u+1
number
if cut >=2 then
command
sel shell id 55 range cyl end1 0 cut2 0 end2 0 cut 0 rad 6.19 &
cyl end1 0 cut2 0 end2 0 cut 0 rad 5.99z -2.19 6.19
sel shell propisotropic 10.5e9,0.25 thickness 0.3 density 2500 rangeid 55
step 1
超前锚秆的命令流值得借鉴；关于 shell 单元，sel shell id 55 range cyl end1 0 cut2 0 end2 0 cut 0 rad 6.19 &
cyl end1 0 cut2 0 end2 0 cut 0 rad 5.99z -2.19 6.19 是说喷射混凝土厚度6.19-5.99=0.2m
但是sel shell propisotropic 10.5e9,0.25 thickness 0.3 density 2500 rangeid 55 是说喷射混凝土厚度0.3m
我一直不是很明白，请高人指点迷津。 支持一下
我也想学习学习 这个方面还正在学习中，呵呵，正好借鉴一下 值得学习哈 仔细看了看，有很多学习的地方，非常感谢。
有几处个人观点：
1，y2=yd_mao+3.5*cos(14)应该为y2=yd_mao+3.5*cos(14*degrad)
2，在开挖过程中，对于超前锚杆，是不是y坐标应该控制一下，也就是yd_mao+3.5*cos(14*degrad)，他可能超出模型外.
3，既然定义了ini_chu_mao，在后面开挖的时候，为什么还要重复的把语句写出来，有什么特殊的地方吗？请指教
4，具体的计算发过来，让大家参考一下啊。谢谢 学习学习哈 学习中！！谢谢分享！！！ 支持支持强烈支持 先学习一下 学习一下，谢谢了！ 先学习一下，谢谢上传，:) 学习中
借鉴
谢谢 学习中:( 谢谢:) 请教，楼主模拟的应该是喷锚支护的隧道，但是，从程序里却看不到对初期喷射混凝土的模拟。
本人现在考虑用shell单元来模拟初喷混凝土，cable模拟锚杆，但是cable和shell的连接，现在还没有想好，不知道大家有没有好的建议 学习中！！！！