[热点探讨]关于施加荷载的命令流---你如果能耐心解决的话，你的水平已经是大牛了 [复制链接]

小弟最近做一个数值模拟岩石破坏实验，但是施加荷载的范围局限在一定的范围，即点荷载试验，但是呢，命令流我写了出来，可是模拟出来的结果和不理想，现在把结果呈上来，看哪位高手能解决这个问题： u8JH~b  
命令流： Z}0{FwW"4  
;开始模拟计算 hC"'cUrcN  
new %`s#p` Ol1  
rest origin_stress.sav Xa>'DO2  
ini xdis=0 ydis=0 zdis=0 L /
PAC  
;初始位移为0 P-T@'}lW  
plot sketch blu T$0)un  
def derive1 d`V.i6u  
s_mod1=E_mod1/(2.0*(1.0+p_ratio1)) 2dHO!A$RF  
b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1)) qm/>\4eLt  
end 0sw;h.VY  
set E_mod1=25e9 p_ratio1=0.26 *A48shfO  
derive1 $2*_7_Qb  
M {'(+a[  
m m )P9]/y  
pro bulk b_mod1 shear s_mod1 coh 6.378e6 fric 34.6 tens 1.75e6 jcJ@A0]  
TtrO_D  
def get1 {ly<%Q7j  
zpnt1 = z_near(0,0,-10) /8]
K}yvR  
zid1 = z_id(zpnt1) S)ipkuj X  
end <L8FI78[*  
get1 9xJtDdy-O  
1!)'dL0mI  
def get2 m.ka%h$  
zpnt2 = z_near(0,0,0) AfO.D?4x  
zid2 = z_id(zpnt2) M]Vi]s  
end Tz6I
7S-w  
get2 dKEy6
C"@  
j|[$P4w}U  
;apply zforce 5000 range id zid F|+B8&-v  
;fix x y z 2K/t[.8  
app zvel -5e-6 range id zid2 Yg5o!A  
app zvel 5e-6 range id zid1 go=xx.WJ  
;hist n 1 2Jiy`(P  
hist gp zdisp 0,0,0 V!. Y M)B  
hist zone szz 0,1,0 q
0ab]g+  
hist zone szz 1,1,0 ?m5"|f\  
plot hist -3 -2 vs 1 'z}9BGR!  
step 100 /0k'w%V{n  
m#[c]v{  
save tun_simulate sd+_NtH  
plot %e25Z.Se$  
create midjiemian2 !!ZGNZ_  
set plane dip=90 dd=0 origin=0 0 0 H~#$AD+H  
;plo add boundary behind ;在当前视图平面的后面绘出图 JT<JS6vw#  
plot ex~"M&^  
add con smax plane ;在当前绘图剖平面内绘出图7 C[Q4OAFG  
show hFycSu  
----------------------上面问题的关键是： Xsanc@w)^C  
app zvel -5e-6 range id zid2 &?p(UY7'"  
app zvel 5e-6 range id zid1 及id1、id2的定义是否合理？ k9H}nP$F  
得出来的中间截面应力图是： O251. hXK  
下载 (58.82 KB) &m36h`tM  
纵截面应力分布，很不合理，下面还算差强人意，有点荷载的味道，上面完全没有施加荷载的意思，一点也体现不 ... P
Ol-S<QV  
]b;a~Y0  
2 分钟前 1jKj'7/K  
$; _{|{Yj  
理想模拟结果应该是类似这种情况的应力分布图： {I0!q"sF  

图片:22.jpg

&x*l{s[  
l{3zlXk3z  
2 分钟前 \rh+\9(  
即在中间产生拉应力，两侧产生应力集中， =ePwGm1:c  

图片:硕士毕业全文终稿_叶剑红.jpg

C[TjcHoA  
希望有高手可以解决这个问题，也可以加QQ 289512820我们可以讨论这个问题。 yMkR)HY  
有朋自网络来，不亦乐乎！！

高手啊，俺才学flac，向你学习。