| 查看: 1503 | 回复: 4 | ||
[求助]
下面是一段移动热源计算力热耦合的命令流，用到生死单元 已有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 |
回复此楼» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有13人回复
-
2026国自然函评费到账
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
2楼2015-12-09 20:44:03
3楼2015-12-10 00:24:01
4楼2015-12-10 10:18:45
5楼2015-12-10 10:19:50
