小弟最近做一个数值模拟岩石破坏实验，但是施加荷载的范围局限在一定的范围，即点荷载试验，但是呢，命令流我写了出来，可是模拟出来的结果和不理想，现在把结果呈上来，看哪位高手能解决这个问题：
;Z"MO@9: 命令流：
gm2|`^Xq$ ;开始模拟计算
?gU raSFU new
]7cciob rest origin_stress.sav
_ AFje ini xdis=0 ydis=0 zdis=0
=
g
& ;初始位移为0
xT_"` @ plot sketch blu
Pg 8boN]} def derive1
OblHN* s_mod1=E_mod1/(2.0*(1.0+p_ratio1))
O|+ZEBP b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))
_ hyqHvP end
Jk-WD"J6 set E_mod1=25e9 p_ratio1=0.26
=1,!EkG derive1
$*G3'G2'iS zN!yOlp5 m m
JZ`SV}\` pro bulk b_mod1 shear s_mod1 coh 6.378e6 fric 34.6 tens 1.75e6
O5vfcX4> >n3GvZ5% def get1
>!@D^3PPA zpnt1 = z_near(0,0,-10)
X Vt ;hO zid1 = z_id(zpnt1)
i a!!jK} end
vT0Op e6m get1
&"JC8 -%Vh-;Ie( def get2
8^+|I, zpnt2 = z_near(0,0,0)
)W 5g-@ zid2 = z_id(zpnt2)
(kB end
3))CD,| get2
oNe:<YT
h6?^rS8U ;apply zforce 5000 range id zid
p?>J86%[ ;fix x y z
j@(S7=^C6% app zvel -5e-6 range id zid2
ed/B.SY app zvel 5e-6 range id zid1
7Z`4K dh . ;hist n 1
2<5LQr hist gp zdisp 0,0,0
)L6
i t hist zone szz 0,1,0
esA^-$ hist zone szz 1,1,0
|(*btdqy3 plot hist -3 -2 vs 1
k^8;3#xG step 100
6]%=q)oL[ rh Q+ylt8I save tun_simulate
DvQV_D plot
:t!J
9 create midjiemian2
L`K)mCr set plane dip=90 dd=0 origin=0 0 0
^i#0aq2} ;plo add boundary behind ;在当前视图平面的后面绘出图
D((/fT)eD plot
_g/d/{-{Q add con smax plane ;在当前绘图剖平面内绘出图7
*V q'%b9 show
SF*mY=1 ----------------------上面问题的关键是：
}v2p]D5n. app zvel -5e-6 range id zid2
Uz,P^\8^$ app zvel 5e-6 range id zid1 及id1、id2的定义是否合理？
hNZ_=
<D! 得出来的中间截面应力图是：
~&p]kmwXSX 下载 (58.82 KB)
Xr Tc5V 纵截面应力分布，很不合理，下面还算差强人意，有点荷载的味道，上面完全没有施加荷载的意思，一点也体现不 ...
^_Lnqk6 S+[,\>pY 2 分钟前
JUA%l jZqa+ nG51 理想模拟结果应该是类似这种情况的应力分布图：
*\xRNgEQ 
]~dB| WB SaF0JPm4z 2 分钟前
O4f9n 即在中间产生拉应力，两侧产生应力集中，
w^&UMX} 
PSu]I?WF 希望有高手可以解决这个问题，也可以加QQ 289512820我们可以讨论这个问题。
] kmAN65c 有朋自网络来，不亦乐乎！！
