| 查看: 6011 | 回复: 2 | ||
ee522522新虫 (初入文坛)
|
[求助]
帮忙啊specified range of 0 to 0 is not permitted. the esel command is ignored！
|
|
/clear,start /filename,weld /title,weld Analysis by "Element Birth and Death" /prep7 /unit,si!采用国际单位制 ！ et,1,13,4!13号二维耦合单元，同时具有温度和位移自由度 et,2,13,4 ! mp,ex,1,2.3e11 mp,ex,2,2.05e11 ! tb,bkin,1,5 tbtemp,20,1 tbdata,1,2800e6,9.055e10 tbtemp,300,1 tbdata,1,2800e6,9.055e10 tbtemp,600,1 tbdata,1,2800e6,9.055e10 tbtemp,900,1 tbdata,1,2800e6,9.055e10 tbtemp,1200,1 tbdata,1,2800e6,9.055e10 ! tb,bkin,2,5 tbtemp,20,1 tbdata,1,335e6,8.2e10 tbtemp,300,1 tbdata,1,335e6,8.2e10 tbtemp,600,1 tbdata,1,335e6,8.2e10 tbtemp,900,1 tbdata,1,335e6,8.2e10 tbtemp,1200,1 tbdata,1,335e6,8.2e10 ! !材料密度 mptemp,1,20,300,600,900,1200 mpdata,dens,1,1,7700,7700,7700,7700,7700 mpdata,dens,2,1,7800,7800,7800,7800,7800 !热膨胀系数 mptemp,1,20,300,600,900,1200 mpdata,alpx,1,1,1.2e-5,1.2e-5,1.2e-5,1.2e-5,1.2e-5 mpdata,alpx,2,1,1.159e-5,1.31e-5,1.467e-5,1.65e-5,1.70e-5 !泊松比 mptemp,1,20,300,600,900,1200 mpdata,nuxy,1,1,0.27,0.27,0.27,0.27,0.27 mpdata,nuxy,2,1,0.25,0.25,0.25,0.25,0.25 !热传导系数 mptemp,1,20,300,600,900,1200 mpdata,kxx,1,1,7.23,7.23,10.5,13.8,17 mpdata,kxx,2,1,50.2,50.2,50.2,50.2,50.2 !比热 mptemp,1,20,300,600,755,1000 mpdata,c,1,1,3.73,3.73,3.73,3.73,3.73 mpdata,c,2,1,480,524,700,1064,602 !由于该13号单元有磁自由度，此处假设一个磁特性，但并不影响我们关心的结果 mptemp,1,20,300,600,900,1200 mpdata,murx,1,1,1,1,1,1,1 mpdata,murx,2,1,1,1,1,1,1 !假设涂层的初始温度为1000 mp,reft,1,20 mp,reft,2,1000 mp,reft,3,20 ! !下面建立几何模型 csys,0 k,1,0,0,0 k,2,0.06,0,0 k,3,0.06,0.005,0 k,4,0,0.005,0 k,5,0,0.010,0 k,6,0.06,0.010,0 a,1,2,3,4 a,3,4,5,6 !划分单元 esize,0.0005 type,2 mat,2 MSHAPE,0,2D mshkey,1 amesh,1 esize,0.00005 type,1 mat,1 amesh,2 eplot ! /solu antype,4 trnopt,full ! !在模型下面施加位移约束 nesl,all *get,miny,node,,mnloc,y nsel,s,loc,y,miny d,all,uy,0 ! !假设模型的下边界温度始终为20 nsel,all *get,miny,node,,mnloc,y nsel,s,loc,y,miny d,all,temp,20 ! ！杀死涂层区域单元 nna=2 esel,all *get,emax,elem,,num,max asel,s,area,,nna esla *get,nse,elem,,count *dim,ne,,nse *dim,nex,,nse *dim,ney,,nse *dim,neorder,,nse mine=0 !下面的do循环用于将涂层区域的单元按其形心Y坐标排序 !模拟涂层由下向上“生长”的过程 *do,i1,1,nse esel,u,elem,,mine !此句有问题？ *get,nse1,elem,,count ii=0 *do,i,1,emax *if,esel(i),eq,1,then ii=ii+1 ne(ii)=i *endif *enddo *do,i,1,nse1 *get,ney(i),elem,ne(i),cent,y *get,nex(i),elem,ne(i),cent,x *enddo miny=1e20 minx=1e20 *do,i,1,nse1 *if,ney(i),lt,miny,then miny=ney(i) minx=nex(i) mine=ne(i) *else *if,ney(i),eq,miny,then *if,nex(i),lt,minx,then miny=ney(i) minx=nex(i) mine=ne(i) *endif *endif *endif *enddo neorder(i1)=mine *enddo ! max_tem=1000 dtl=1e-3 dt=8.33*20e-8 t=0 esel,all eplot /auto,1 /replot *do,i,1,nse ekill,neorder(i) esel,s,live eplot *enddo allsel,all outress,all,all ic,all,temp,20 kbc,1 timint,0,struct timint,1,therm timint,0,mag tintp,0.005,,,1,0.5,0.2 ! ! nsub1=2 nsub2=40 *do,i,1,nse ealive,neorder(i) esel,s,live eplot esel,all !下面的求解用于建立温度的初始条件 t=t+dtl time,t nsubst,1 *do,j,1,4 d,nelem(neorder(i),j),temp,max_tem *enddo solve !下面的求解用于保证初始的升温速度为零 t=t+dtl time,t solve !下面的步骤用于求解温度分布 *do,j,1,4 ddele,nelem(neorder(i),j),temp *enddo t=t+dt-2*dtl time,t nsubst,nsub1 solve *enddo t=t+50000 time,t nsubst,nsub2 solve finish [ Last edited by ee522522 on 2013-11-8 at 10:27 ] |
回复此楼» 猜你喜欢
生物学303求调剂，一志愿华中农微生物，六级已过，有科研有文章，党员
已经有3人回复
-
求调剂
已经有8人回复
-
求调剂
已经有6人回复
-
353求调剂
已经有10人回复
-
298求调剂
已经有6人回复
-
生物学308分求调剂（一志愿华东师大）
已经有5人回复
-
一志愿电子科技大学085600材料与化工 329分求调剂
已经有9人回复
-
284求调剂
已经有6人回复
-
0817化学工程与技术求调剂，一志愿中海洋319
已经有11人回复
-
304求调剂（085602，过四级，一志愿985）
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 大神们，求帮忙啊，谢谢
已经有3人回复
- 新手求助高斯输入文件写法，求大神帮忙
已经有8人回复
- 求两个值中最大的值但编写的代码不知道哪里错了，求大神帮忙。
已经有10人回复
- 求助fortran大神，我是初学者，帮我分析一下下面的代码
已经有5人回复
- 求大神帮忙设计一个工装夹具
已经有6人回复
- 求各位大神帮忙下一篇文章
已经有3人回复
- 跪求各位大神解答~本人高三学生 高考想考建筑学 求各位大神帮忙解答一下
已经有9人回复
- 磷酸铁锂电池40度的容量居然比23度低，求大神帮忙解释！有金币感激！
已经有12人回复
- 跪求大神帮忙解决A类错误
已经有3人回复
- 求大神帮忙研究一下关于铝合金的抗弯强度的问题！！！！！！！
已经有13人回复
- 求各位大神帮忙 检测猪肝中铜铁锌铅的前处理方法有哪些！！！
已经有8人回复
- 哪位大神能帮忙总结下常见天然产物的预检识方法啊，膜拜~~~~
已经有3人回复
- 初学者请各位高手帮忙看一下GSAS精修的结果，请指点
已经有14人回复
2楼2015-05-10 16:11:24
3楼2018-04-09 09:29:43
10