论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 1691阅读
  • 1回复

[求助]彭文斌老师课本例题问题求解答 [复制链接]

上一主题 下一主题
离线miner007
 

发帖
72
土币
3720
威望
12
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2014-06-04

new
; Uniaxial test of strain-softening material
title
应变软化材料单轴实验
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 4 0 p3 0 0 1 size 12 30 12
gen zone reflect norm 1,0,0
gen zone reflect norm 0,0,1
model ss ;应变硬化/软化模型
pro den 2500 bulk 2e8 she 1e8 co 2e6 fric 45 ten 1e6 dil 10
pro 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 ; cylinder radius = 1
end
hist n 1
hist gp ydisp 0,0,0
hist ax_str
hist gp xdisp 1,1,0
plot hist -2 vs 1 ;axial stress vs axial disp.
step 5000
save beforeplzones.sav
; Plot of plastic region as zones with strain > 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
plo crea qqq
plo add surf red range group yield
plo add axes green
plo set rot 123 313 3
plo set mag 1.5
ret
上述是彭文斌书例8.1 请教大神为什么红色区域这儿不是固定了y=0和y=4这两个面吗?为什么还可以赋值他们的速度?
离线dsaaamj

发帖
30
土币
5352
威望
1
原创币
0
只看该作者 1楼 发表于: 2014-09-09
书的 9.3有说明。不能apply,但是能ini
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      浙公网安备 33010602003799号 浙ICP备14021682号-1

      工业和信息化部备案管理系统网站