关于施加荷载的命令流---你如果能耐心解决的话，你的水平已经是大牛了
本帖最后由 sun99shine2008 于 2009-9-28 10:59 编辑小弟最近做一个数值模拟岩石破坏实验，但是施加荷载的范围局限在一定的范围，即点荷载试验，但是呢，命令流我写了出来，可是模拟出来的结果和不理想，现在把结果呈上来，看哪位高手能解决这个问题： 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 60]
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的定义是否合理？
得出来的中间截面应力图是：
http://www.yantubbs.com/attachment-2/Fid_75/75_20963_9a1bcc31636fbbd.jpg
纵截面应力分布，很不合理，下面还算差强人意，有点荷载的味道，上面完全没有施加荷载的意思，一点也体现不出来。
理想模拟结果应该是类似这种情况的应力分布图：
http://www.yantubbs.com/attachment-2/Fid_75/75_20963_aabfef04dd141db.jpg
即在中间产生拉应力，两侧产生应力集中，
希望有高手可以解决这个问题，也可以加QQ 289512820我们可以讨论这个问题。 ,"
有朋自网络来，不亦乐乎！！ 你这个问题很简单，fish也很简单，既然都知道要施加荷载的节点在哪，干嘛不去，打开看看，那个节点具体是多少呢，建议你把获得节点号的fish去掉，直接去看那个节点号是多少。 谢谢楼上的关注
关键我施加点荷载的时候，是施加一个范围，不可能就在一个节点上施加荷载；
我的意思是，用fish获得一定的节点号范围，从而在这个范围内施加荷载 有时即使是同一个面上，节点号可能会出现不连续，你这样应当用坐标来控制。 zpnt1 = z_near(0,0,-10)
zid1 = z_id(zpnt1)
app zvel 5e-6 range id zid1
我是用坐标来获取周围的节点啊，莫非哥们认为这个语句有什么问题？
还是你有什么好的建议，要不直接把命令写出来看看？
真心谢谢你一直关注这个贴:handshake ??请这位大哥帮忙哈
4# cloud328 版大
参考了你的建议直接施加在节点上,通过查看节点号，上下两个面中心的节点号为1&621
app zvel 5e-6 range gridpointid 621
app zvel 5e-6 range gridpointid 1
不知道为什么，还是效果还是不好。
页:
[1]