24小时热门版块排行榜

返回列表

hmq891204

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 15
在线: 3.8小时
虫号: 3311261
注册: 2014-07-07
专业: 制造系统与自动化

[求助] 下面是一段移动热源计算力热耦合的命令流，用到生死单元已有1人参与

我想查看杀死后的结构的温度场，应力场等等，定义了一个数组，然后将结果存到了单元表中，发现移动的载荷步不增加了，一直显示第一载荷步的几个子步，时间还是在增加的？
问题1：命令流是否有错误，感觉是数组定义的问题
问题2 ：尧实现我需要的查看结果的功能，还有什么其他方法吗？
做的是切割时，力热移动载荷下的耦合分析，用到CZM模型，有相同的可以讨论下，这是我所有的金币
下面是我的命令流，求大神帮忙，毕业课题，搞不出来了

！下面是参数定义
/prep7
*set,l,0.0049 !板件宽度
*set,w,0.0098 !板件长度，切割方向
*set,h1,0.0005  !emc
*set,h2,0.0002  !cu
*set,wl_t,0.00032  !砂轮厚度
*set,vs,0.06 !砂轮切割速度
*set,f_z,0.04291    !切割道侧面切向力，z方向，141个节点
*set,f_z2,0.05216    !切割道正面切向力，z方向，116个节点
*set,pre_x,13.34    !切割法向力
*set,hflu_s,10.861e5  !热流密度
*set,hflu_s2,16.631e5  !热流密度
*set,b,0.00049 !切割砂轮宽度,即一个载荷步长
*set,bi,1  !流入工件的热量占总热量的比
*set,l1,w       !设定磨削长度
*set,loadstep,w/b !载荷步数
*set,t,w/vs    !设定时间常数
*set,tim,0             !设置初始时间
*set,es,0.0001       !设置单元划分大小
！下面是定义数组的命令
/solu
*dim,tempe,char,18  !定义数组存储单元表，计算18各载荷步，不知道是不是正确
*do,n,1,18,1
tempe(n)='temp_%n%'
*enddo

*dim,dparamr,char,18  !定义数组存储单元表
*do,m,1,18,1
dparamr(m)='dpa_%m%'
*enddo
！下面还有n个类似的定义，略去了！
allsel,all
/solu
!施加移动热源/力循环步
*do,i,1,loadstep-18,1             !循环开始
*set,es,0.0001
*set,tim,tim+t/loadstep    !设置载荷步结束时间
time,tim                !指定载荷步结束时间
!施加移动载荷
nsel,all
nsel,s,loc,y,(i-1)*5*es,i*5*es
nsel,r,loc,x,l/2,l/2
sf,all,hflux,hflu_s*bi
sf,all,pres,pre_x/(5*es*(h1+h2))
f,all,fz,f_z
nsel,all
nsel,s,loc,y,(i-1)*5*es,i*5*es
nsel,r,loc,x,l/2+wl_t,l/2+wl_t
sf,all,hflux,hflu_s*bi
sf,all,pres,pre_x/(5*es*(h1+h2))
f,all,fz,f_z
nsel,all
nsel,s,loc,x,l/2+0.01*es,l/2+wl_t-0.01*es
nsel,r,loc,y,i*5*es-0.5*es,i*5*es
sf,all,hflux,hflu_s2*bi
sf,all,pres,pre_x/(wl_t*(h1+h2))
f,all,fz,f_z2
allsel,all
solve
save
！下面将需要的结果存储到数组中
/post1 !进入post1,单元表不能超过200个，24*18，定义的太多
etable,tempe(i),temp, !将所有单元温度存入etable
etable,dparamr(i),NMISC,70,71,72,73!储存损伤参数
etable,deni(i),NMISC,140,141,142,143!存储1,2型能量释放率
etable,denii(i),NMISC,144,145,146,147
etable,dtstartr(i),NMISC,66,67,68,69!存储脱粘时程
etable,contpresr(i),cont,pres    !存储接触压力
etable,contsfricr(i),cont,sfric    !存储接触摩擦应力
etable,contgapr(i),cont,gap    !存储接触间隙
etable,sxyr(i),s,xy
etable,ssxyr(i),ss,xy
etable,seqvr(i),s,eqv
！下面删除载荷
allsel,all
nsel,s,loc,y,(i-1)*5*es,i*5*es
nsel,r,loc,x,l/2,l/2
sfedele,all,all,pres
sfedele,all,all,hflux
fdele,all,fz
allsel,all
nsel,s,loc,y,(i-1)*5*es,i*5*es
nsel,r,loc,x,l/2+wl_t,l/2+wl_t
sfedele,all,all,pres
sfedele,all,all,hflux
fdele,all,fz
allsel,all
nsel,s,loc,x,l/2+0.01*es,l/2+wl_t-0.01*es
nsel,r,loc,y,i*5*es-0.5*es,i*5*es
sfedele,all,all,pres
sfedele,all,all,hflux
fdele,all,fz
allsel,all
！下面删除切割去掉的单元
nsel,s,loc,y,i*5*es,(i+1)*5*es-0.5*es
nsel,r,loc,x,l/2+0.01*es,l/2+wl_t-0.01*es
esln,s
/solu    !重新进入求解器
ekill,all !杀死选择单元,并约束杀死单元的节点自由度
allsel,all
nsel,s,loc,y,i*5*es-0.5*es,(i+1)*5*es-0.5*es
nsel,r,loc,x,l/2+0.01*es,l/2+wl_t-0.01*es
d,all,ux,,,,,uy,uz
allsel,all
*enddo

回复此楼