小弟最近做一个数值模拟岩石破坏实验，但是施加荷载的范围局限在一定的范围，即点荷载试验，但是呢，命令流我写了出来，可是模拟出来的结果和不理想，现在把结果呈上来，看哪位高手能解决这个问题：
zf3v5Hk 命令流：
2VW}9O ;开始模拟计算
rt
JtK6t new
+cb6?? H rest origin_stress.sav
9rh}1eo7 ini xdis=0 ydis=0 zdis=0
v}WR+)uFQ ;初始位移为0
E .1J 2Ne plot sketch blu
a^ ,( v def derive1
F ~ Lx|)0M s_mod1=E_mod1/(2.0*(1.0+p_ratio1))
~>9_(L b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))
L(GjZAP end
a >\vUv * set E_mod1=25e9 p_ratio1=0.26
12Y derive1
P%K4[c W~ R_=6G ZH$G m m
! 8`3GX:B_ pro bulk b_mod1 shear s_mod1 coh 6.378e6 fric 34.6 tens 1.75e6
O-)-YVU IK*07h/! def get1
1CLL%\V zpnt1 = z_near(0,0,-10)
boG_f@dv( zid1 = z_id(zpnt1)
h gJ[LU| > end
6q
2_W X get1
E*i#?u iR" N13 def get2
}i!J/tJ)b zpnt2 = z_near(0,0,0)
d z Zb zid2 = z_id(zpnt2)
HSq}7S&U end
=B. F;4 0 get2
2w F8 P) !A48TgAeE ;apply zforce 5000 range id zid
X$/2[o#g ;fix x y z
Haqm^Ky$ app zvel -5e-6 range id zid2
1;VHM' app zvel 5e-6 range id zid1
hMtf.3S7c ;hist n 1
6;b ~Ht hist gp zdisp 0,0,0
+^;JS3p@\ hist zone szz 0,1,0
tBpC: SG hist zone szz 1,1,0
\ TU3rk&X plot hist -3 -2 vs 1
I}R0q step 100
J4 Dry< w3$ save tun_simulate
#:SNHM^>< plot
o.I6ulY8 create midjiemian2
Hc`A3SMR set plane dip=90 dd=0 origin=0 0 0
8V~vXnkM ;plo add boundary behind ;在当前视图平面的后面绘出图
s&j-\bOic9 plot
e1P"[|9>R add con smax plane ;在当前绘图剖平面内绘出图7
gCyW Vp show
W^"AU;^V56 ----------------------上面问题的关键是：
oY, %Iq app zvel -5e-6 range id zid2
*?M GMhE app zvel 5e-6 range id zid1 及id1、id2的定义是否合理？
<_Q1k> 得出来的中间截面应力图是：
Ha )ANAD 下载 (58.82 KB)
m^qBx A 纵截面应力分布，很不合理，下面还算差强人意，有点荷载的味道，上面完全没有施加荷载的意思，一点也体现不 ...
kwsp9 0) <I;2{*Q I2 2 分钟前
G}p\8Q}' 13Lr }M& 理想模拟结果应该是类似这种情况的应力分布图：
{5
sO 
u G<}N= Po%(~ )S> 2 分钟前
t45Z@hmcW 即在中间产生拉应力，两侧产生应力集中，
cv`~y'?D 
fum0>tff 希望有高手可以解决这个问题，也可以加QQ 289512820我们可以讨论这个问题。
PdvqD a8 有朋自网络来，不亦乐乎！！
