24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ee522522

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 21
帖子: 14
在线: 33.6小时
虫号: 2700944
注册: 2013-10-06
性别: GG
专业: 数值预报与数值模拟

[求助] 帮忙啊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 ]

回复此楼