关于施加荷载的命令流---你如果能耐心解决的话，你的水平已经是大牛了

sun99shine2008

小弟最近做一个数值模拟岩石破坏实验，但是施加荷载的范围局限在一定的范围，即点荷载试验，但是呢，命令流我写了出来，可是模拟出来的结果和不理想，现在把结果呈上来，看哪位高手能解决这个问题： e[Ul"pMvS`  
命令流：
;开始模拟计算 ,  
new os0"haOI9h  
rest origin_stress.sav 1
ini xdis=0 ydis=0 zdis=0

;初始位移为0
plot sketch blu
def derive1

s_mod1=E_mod1/(2.0*(1.0+p_ratio1))
b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))
end

set E_mod1=25e9 p_ratio1=0.26

derive1


m m

pro bulk b_mod1 shear s_mod1 coh 6.378e6 fric 34.6 tens 1.75e6

   
def get1  
zpnt1 = z_near(0,0,-10)
zid1 = z_id(zpnt1)

end
get1  

def get2
zpnt2 = z_near(0,0,0)
zid2 = z_id(zpnt2
end

get2  
;apply zforce 5000 range id zid  
;fix x y z
app zvel -5e-6 range id zid2  
app zvel 5e-6 range id zid1

;hist n 1
hist gp zdisp 0,0,0
hist zone szz 0,1,0
hist zone szz 1,1,0

plot hist -3 -2 vs 1
step 1000

save tun_simulate
plot

create midjiemian2
set plane dip=90 dd=0 origin=0 0 0

;plo add boundary behind ;在当前视图平面的后面绘出图

plot 2 _Z 6  0]  
add con smax plane ;在当前绘图剖平面内绘出图7
show ,U}8( D~:  
----------------------上面问题的关键是：
app zvel -5e-6 range id zid2 9
app zvel 5e-6 range id zid1 及id1、id2的定义是否合理？
得出来的中间截面应力图是：

纵截面应力分布，很不合理，下面还算差强人意，有点荷载的味道，上面完全没有施加荷载的意思，一点也体现不出来。  
理想模拟结果应该是类似这种情况的应力分布图：

即在中间产生拉应力，两侧产生应力集中，
希望有高手可以解决这个问题，也可以加QQ 289512820我们可以讨论这个问题。 ,"
有朋自网络来，不亦乐乎！！

cloud328

你这个问题很简单，fish也很简单，既然都知道要施加荷载的节点在哪，干嘛不去，打开看看，那个节点具体是多少呢，建议你把获得节点号的fish去掉，直接去看那个节点号是多少。

sun99shine2008

谢谢楼上的关注
   关键我施加点荷载的时候，是施加一个范围，不可能就在一个节点上施加荷载；
   我的意思是，用fish获得一定的节点号范围，从而在这个范围内施加荷载

cloud328

有时即使是同一个面上，节点号可能会出现不连续，你这样应当用坐标来控制。

sun99shine2008

zpnt1 = z_near(0,0,-10)
zid1 = z_id(zpnt1)
app zvel 5e-6 range id zid1
   我是用坐标来获取周围的节点啊，莫非哥们认为这个语句有什么问题？
  还是你有什么好的建议，要不直接把命令写出来看看？
  真心谢谢你一直关注这个贴

sun99shine2008

??请这位大哥帮忙哈
4# cloud328

sun99shine2008

版大  
  参考了你的建议  直接施加在节点上,通过查看节点号，上下两个面中心的节点号为1&621
app zvel 5e-6 range gridpoint  id 621
app zvel 5e-6 range gridpoint  id 1
  不知道为什么，还是效果还是不好。