小弟最近做一个数值模拟岩石破坏实验，但是施加荷载的范围局限在一定的范围，即点荷载试验，但是呢，命令流我写了出来，可是模拟出来的结果和不理想，现在把结果呈上来，看哪位高手能解决这个问题：
z@iu$DZ 命令流：
=t\HtAXn[ ;开始模拟计算
@2c Gx/1# new
k;/U6,LQ* rest origin_stress.sav
p@!{Sh ini xdis=0 ydis=0 zdis=0
NitWIj[U; ;初始位移为0
,{MA90! plot sketch blu
?N oG. def derive1
ig3HPlC s_mod1=E_mod1/(2.0*(1.0+p_ratio1))
}D j W b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))
#)QR^ss)iw end
g+|1khS) set E_mod1=25e9 p_ratio1=0.26
~ \z7$9Q derive1
p "EQ 6_f TC;2K,.#k m m
yPL1(i; pro bulk b_mod1 shear s_mod1 coh 6.378e6 fric 34.6 tens 1.75e6
?a+>%uWt J=|PZ2" def get1
7KXc9:p + zpnt1 = z_near(0,0,-10)
E8j >Toz zid1 = z_id(zpnt1)
ohklLZoZ end
*U|K~dl]K get1
&[ ejxK" gZ F-zhnC def get2
fw Ooi 'jb zpnt2 = z_near(0,0,0)
MgyV {` zid2 = z_id(zpnt2)
ya8MjGo end
i;>Yx# get2
U<QO@5 B%QvFxZz ;apply zforce 5000 range id zid
H5j6$y|I|N ;fix x y z
\tH^w@j47 app zvel -5e-6 range id zid2
'F.Da#st!} app zvel 5e-6 range id zid1
musZCg$ ;hist n 1
_0p8FhNt hist gp zdisp 0,0,0
v#J 2yg hist zone szz 0,1,0
i_8v > F hist zone szz 1,1,0
.6z8fjttOC plot hist -3 -2 vs 1
HfEU[p7) step 100
tJ`tXO #lXwBfBMf save tun_simulate
S) ZcH plot
;DbEP. %u$ create midjiemian2
Ry9kGdqO set plane dip=90 dd=0 origin=0 0 0
hF|N81T ;plo add boundary behind ;在当前视图平面的后面绘出图
1S
.~Vh0Q, plot
?o307 r add con smax plane ;在当前绘图剖平面内绘出图7
2>3#/ I9Y show
*\><MXx ----------------------上面问题的关键是：
6wK>SW)#&j app zvel -5e-6 range id zid2
%>u (UmFO app zvel 5e-6 range id zid1 及id1、id2的定义是否合理？
t Gv5pe*r 得出来的中间截面应力图是：
~Ry?}5&: 下载 (58.82 KB)
I(!i"b9 纵截面应力分布，很不合理，下面还算差强人意，有点荷载的味道，上面完全没有施加荷载的意思，一点也体现不 ...
DtLga[M OOJg%y*H 2 分钟前
r&6X|2@ If[4]-dq 理想模拟结果应该是类似这种情况的应力分布图：
~~,] b 
(U bz@s^ 6j z6
2 分钟前
KG7 ~)g 即在中间产生拉应力，两侧产生应力集中，
L|O[u^ 
x{y}pH "H 希望有高手可以解决这个问题，也可以加QQ 289512820我们可以讨论这个问题。
&f48MtE 有朋自网络来，不亦乐乎！！
