应变软化材料的单轴试验小例子中最后group yield的yield无法识别怎么回事？求助
;应变软化材料的单轴试验new
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 4 0 p3 0 0 1 size 12 30 12
plot add surface yellow
plot add axes
plot show
gen zone reflect norm 1 0 0
gen zone reflect norm 0 0 1
mo ss
prop den 2500 bulk 2e8 she 1e8 co 2e6 fric 45 ten 1e6 dil 10
prop ftab 1 ctab 2 dtab 3
table 1 0 45 .05 42 .1 40 1 40
table 2 0 2e6 .05 1e6 .1 5e5 1 5e5
table 3 0 10 .05 3 .1 0
fix x y z range y -.1 .1
fix x y z range y 3.9 4.1
ini yvel 2.5e-5 range y -.1 .1
ini yvel -2.5e-5 range y 3.9 4.1
def ax_str
str=0
pnt=gp_head
loop while pnt # null
if gp_ypos(pnt)<0.1 then
str=str + gp_yfunbal(pnt)
endif
pnt=gp_next(pnt)
endloop
ax_str=str/pi
end
hist n 1
hist gp ydisp 0 0 0
hist ax_str
hist gp xdis 1 1 0
plot hist -2 vs 1
step 100
;绘出切应变率等值线图
save beforeplzones.sav
plot create shcrate
plot add contour ssr outline on
plot set rot 123 313 3
plot set mag 0.64
plot show
;绘出塑性区应变>0.2的单元体
def ShowPlasticZones
zp=zone_head
loop while zp # null
if z_prop(zp,'es_plastic')>0.2
z_group(zp)= 'yield'
else
z_group(zp)= 'other'
endif
zp=z_next(zp)
endloop
end
ShowPlasticZones
plot crea qqq
plot add surf red range group yield
运算到此提示unrecgnized parameter 6 (yidld)请问是那里有问题吗？
plot add axes green
plot set rot 123 313 3
plot set mag 1.5
plot show
ret 1# GODSAVEME
我也做这个例子了 也是很纳闷啊希望高手指点 原因是加载过程中没有出现塑性剪切应变。可以通过查单元属性得到。也就是ShowPlasticZones没找到相应的单元 3# hh199422
那怎么改呢 ？？
页:
[1]