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

小弟最近做一个数值模拟岩石破坏实验，但是施加荷载的范围局限在一定的范围，即点荷载试验，但是呢，命令流我写了出来，可是模拟出来的结果和不理想，现在把结果呈上来，看哪位高手能解决这个问题： y ~!Zg}o  
命令流： <s<n  
;开始模拟计算 AogVF  
new 9H`XeQ.  
rest origin_stress.sav *]/zc1Q4M  
ini xdis=0 ydis=0 zdis=0 {go;C}  
;初始位移为0 iN8zo:&Z  
plot sketch blu nwRc%C``UK  
def derive1 *kDCliL  
s_mod1=E_mod1/(2.0*(1.0+p_ratio1)) CSq4x5!_7>  
b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1)) 70tH:Z)"  
end >rKIG~P_  
set E_mod1=25e9 p_ratio1=0.26 j0evq+  
derive1 mQ26K~  
dV$gB<iS  
m m 9} M?
P  
pro bulk b_mod1 shear s_mod1 coh 6.378e6 fric 34.6 tens 1.75e6 %^GfS@t  
s/ qYa])  
def get1 9BBmw(M}  
zpnt1 = z_near(0,0,-10) ex (.=X 1  
zid1 = z_id(zpnt1) 3/e.38m|  
end cZU=o\  
get1 9oR@U
W1  
'RYIW/a  
def get2 xvl#w  
zpnt2 = z_near(0,0,0) l=)x
o@6  
zid2 = z_id(zpnt2) [<6^qla  
end /quc}"__  
get2 1.{z3_S21:  
][h%UrV  
;apply zforce 5000 range id zid ^-Kf']hU  
;fix x y z j8{i#;s!"  
app zvel -5e-6 range id zid2 ?5|>@>  
app zvel 5e-6 range id zid1 J1RJ*mo7,  
;hist n 1 I%X6T@P  
hist gp zdisp 0,0,0 Jg\zdi:t  
hist zone szz 0,1,0 ~{B7 k:  
hist zone szz 1,1,0 sRL`dEl4l  
plot hist -3 -2 vs 1 A_ N;   
step 100 2E/"hQw  
a /l)qB#  
save tun_simulate lN?qp'%H`  
plot >j(_[z|v3  
create midjiemian2 )j(7]uX`  
set plane dip=90 dd=0 origin=0 0 0 a#y;d
K  
;plo add boundary behind ;在当前视图平面的后面绘出图 0gy/:T  
plot D] jzAx  
add con smax plane ;在当前绘图剖平面内绘出图7 {V$|3m>:*  
show Tx=-Bb~;  
----------------------上面问题的关键是： E+R1 !.  
app zvel -5e-6 range id zid2 ,+vy,<e&  
app zvel 5e-6 range id zid1 及id1、id2的定义是否合理？ m=A(NKZ  
得出来的中间截面应力图是： K'Tm_"[u  
下载 (58.82 KB) %.Fi4}+O  
纵截面应力分布，很不合理，下面还算差强人意，有点荷载的味道，上面完全没有施加荷载的意思，一点也体现不 ... RJ ||}5  
}{qZ[/JwqN  
2 分钟前 [.'|_l  
Ng>5?F^v  
理想模拟结果应该是类似这种情况的应力分布图： 3kIN~/<R+7  

图片:22.jpg

(vP
N5F  
6?mibvK  
2 分钟前 zRl3KjET  
即在中间产生拉应力，两侧产生应力集中， THbh%)Zv+  

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

qL&[K>2z  
希望有高手可以解决这个问题，也可以加QQ 289512820我们可以讨论这个问题。 _8riUt  
有朋自网络来，不亦乐乎！！

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