24小时热门版块排行榜

返回列表

秋水0718

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 4135.8
帖子: 643
在线: 102.1小时
虫号: 1572340
注册: 2012-01-10
性别: GG
专业: 电能储存与节电技术

[求助] 切割仿真加移动热源力热耦合，杀死单元后处理时改怎么处理不同载荷步下的结果的问题

图片是查到的一个解决方法，没有具体说怎么做？
对不同载荷步做后处理，一定要先确信数据库中有和该载荷步生死状态相匹配的所有单元的生死状态（对于改变生死状态的每一载荷步应作一数据库副本）
但是没有具体的解决方法，求教？？

有类似的高手可以私聊请教下，270954803

我想查看杀死后的结构的温度场，应力场等等，定义了一个数组，然后将结果存到了单元表中，发现移动的载荷步不增加了，一直显示第一载荷步的几个子步，时间还是在增加的？
问题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

bw178h3311261_1449665012_204.jpg

回复此楼