小弟最近做一个数值模拟岩石破坏实验，但是施加荷载的范围局限在一定的范围，即点荷载试验，但是呢，命令流我写了出来，可是模拟出来的结果和不理想，现在把结果呈上来，看哪位高手能解决这个问题：
S=eY`,'#R 命令流：
F2! _Z= ;开始模拟计算
Y.[^3 new
8iqx*8} rest origin_stress.sav
]\ sBl ini xdis=0 ydis=0 zdis=0
L*D-RY W ;初始位移为0
FTtYzKX(bv plot sketch blu
"Z& { def derive1
dfO@Yo-?*' s_mod1=E_mod1/(2.0*(1.0+p_ratio1))
g5;
W6QX b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))
=5[}&W end
7'J}|m{7 set E_mod1=25e9 p_ratio1=0.26
H}dsd=yO derive1
!3kyPoq+ FxTOc@< m m
9fzbR ~s pro bulk b_mod1 shear s_mod1 coh 6.378e6 fric 34.6 tens 1.75e6
vz#rbBY*; 2%rAf8= def get1
Ae <v zpnt1 = z_near(0,0,-10)
b U-Cd zid1 = z_id(zpnt1)
N$6Rg1 end
[JO' ta get1
4%KNHeaN ne oT\HV def get2
( T VzYm
y zpnt2 = z_near(0,0,0)
HulN84 zid2 = z_id(zpnt2)
4 "@BbVYR end
PHyS^J` get2
Hd
U1gV> ujXC#r& ;apply zforce 5000 range id zid
QDzFl1\P ;fix x y z
q-KN {y/ app zvel -5e-6 range id zid2
C8m8ys app zvel 5e-6 range id zid1
j;=+5PY ;hist n 1
fM]zD/ g hist gp zdisp 0,0,0
94 H\,}i 8 hist zone szz 0,1,0
r\F`xtR( hist zone szz 1,1,0
2w $o;zz1 plot hist -3 -2 vs 1
JrQd7 step 100
)U6T]1 {-PD3 [f" save tun_simulate
|p4F^! 9 plot
R.Ao%VT create midjiemian2
2~<?E`+ set plane dip=90 dd=0 origin=0 0 0
$-|`#|CBd ;plo add boundary behind ;在当前视图平面的后面绘出图
2ZNTj u7h plot
IR;lt 3 add con smax plane ;在当前绘图剖平面内绘出图7
u `1cXL[' show
Dr^#e ----------------------上面问题的关键是：
w-Y-;*S app zvel -5e-6 range id zid2
N32!*TsWs app zvel 5e-6 range id zid1 及id1、id2的定义是否合理？
8{Bcl 5]< 得出来的中间截面应力图是：
gI~ R u8 下载 (58.82 KB)
B<S uNbR 纵截面应力分布，很不合理，下面还算差强人意，有点荷载的味道，上面完全没有施加荷载的意思，一点也体现不 ...
^G .PdX$M 3fgVvt-2 2 分钟前
; f:}gMK X-|Lg.s 理想模拟结果应该是类似这种情况的应力分布图：
bF;g.-.2 
<5?.s<
y$" %>_ZUu3M 2 分钟前
(]n^_G#-$ 即在中间产生拉应力，两侧产生应力集中，
c | 
5 6;lB$)" 希望有高手可以解决这个问题，也可以加QQ 289512820我们可以讨论这个问题。
l ~kxt2& 有朋自网络来，不亦乐乎！！
