找回密码
 注册
Simdroid-非首页
查看: 59|回复: 1

*if,esel(i),eq,1,then 大家知道这句命令流什么意思吗?1代表什么?求指教

[复制链接]
发表于 2013-3-17 10:10:30 | 显示全部楼层 |阅读模式 来自 山东济南
*do,i,1,,50
*if,esel(i),eq,1,then
.
.
.
*endif
*enddo
 楼主| 发表于 2013-3-18 10:50:42 | 显示全部楼层 来自 山东济南
Simdroid开发平台
/filn,a11
/units,si
/prep7
et,1,solid5
et,2,solid5
mptemp,1,20,500,1000,1500,2000
mpdata,ex,1,1,2.06e11,1.7e11,1e11,0.5e11,0.1e11
mpdata,dens,1,1,7800,7800,7800,7800,7800
mpdata,nuxy,1,1,0.3,0.3,0.3,0.3,0.3
mpdata,kxx,1,1,16.3,16.3,16.3,16.3,16.3
mpdata,alpx,1,1,1.06e-5,1.06e-5,1.06e-5,1.06e-5,1.06e-5
mpdata,c,1,1,502,502,502,502,502
tb,bkin,1,5
tbtemp,20,1
tbdata,1,1.4e9,2.06e10
tbtemp,500,2
tbdata,1,1e9,1.7e10
tbtemp,1000,3
tbdata,1,0.5e9,0.9e10
tbtemp,1500,4
tbdata,1,0.07e9,0.2e10
tbtemp,2000,5
tbdata,1,0.007e9,0.01e10
mpdata,ex,2,1,2.07e11,1.8e11,0.9e11,0.3e11,0.1e11
mpdata,dens,2,1,7811,7811,7811,7811,7811
mpdata,nuxy,2,1,0.3,0.3,0.3,0.3,0.3
mpdata,kxx,2,1,14.9,14.9,14.9,14.9,14.9
mpdata,alpx,2,1,1.08e-5,1.08e-5,1.08e-5,1.08e-5,1.08e-5
mpdata,c,2,1,522,522,522,522,522
tb,bkin,2,5
tbtemp,20,1
tbdata,1,1.5e9,2.07e10
tbtemp,500,2
tbdata,1,1.1e9,1.8e10
tbtemp,1000,3
tbdata,1,0.6e9,1e10
tbtemp,1500,4
tbdata,1,0.08e9,0.3e10
tbtemp,2000,5
tbdata,1,0.008e9,0.02e10
k,1,0,0,0
k,2,0.15
k,3,0.137,0.025
k,4,0,0.025
k,5,0.152
k,6,0.3
k,7,0.3,0.025
k,8,0.152,0.025
k,9,0.145,0.013
a,1,2,3,4
a,5,6,7,8,9
k,14,0.145,0.0225
larc,10,11,14
l,10,13
l,12,13
l,12,9
l,11,9
k,20,0,0,0.01
l,1,20
vglue,all
nummrg,all
numcmg,all
allsel
/image,save,vmesh,bmp
/solu
antype,4
nsel,s,loc,x,0
d,all,ux,0
nsel,r,loc,y,0
d,all,uy,0
nsel,r,loc,z,0
d,all,uz,0
nsel,all
nsel,s,loc,x,0
d,all,temp,20
nsel,r,loc,x,0.3
d,all,temp,20
esel,all
*get,emax,elem,,num,max
vsel,s,volu,,3
eslv
*get,nse,elem,,count
*dim,ne,,nse
*dim,nex,,nse
*dim,ney,,nse
*dim,nez,,nse
*dim,neorder,,nse
mine=0
*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,nez(i),elem,ne(i),cent,z
*get,ney(i),elem,ne(i),cent,y
*get,nex(i),elem,ne(i),cent,x
*enddo
minz=1e20
miny=1e20
minx=1e20
*do,i,1,nse1
*if,nez(i),lt,minz,then
minz=nez(i)
miny=ney(i)
minx=nex(i)
mine=ne(i)
*else
*if,nez(i),eq,minz,then
*if,ney(i),lt,miny,then
minz=nez(i)
miny=ney(i)
minx=nex(i)
mine=ne(i)
*else
*if,ney(i),eq,miny,then
*if,nex(i),lt,minx,then
minz=nez(i)
miny=mey(i)
minx=nex(i)
mine=ne(i)
*endif
*endif
*endif
*endif
*endif
*enddo
neorder(i1)=mine
*enddo
max_tem=1600
dt1=1e-8
dt=0.1111
t=0
esel,all
eplot
*do,i,1,nse
ekill,neorder(i)
esel,s,live
*enddo
allsel,all
outres,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
esel,all
t=t+dt1
time,t
nsubst,1
*do,j,1,8
d,nelem(neorder(i),j),temp,max_tem
*enddo
solve
t=t+dt1
time,t
solve
*do,j,1,8
ddele,nelem(neorder(i),j),temp
*enddo
t=t+dt-2*dt1
time,t
nsubst,nsub1
solve
*enddo
/image,save,1.bmp
t=t+10000
time,t
nsubst,nsub2
solve
finish
/post1
set,last
plnsol,temp,,0,1.0
/image,save,temp1,bmp
plnsol,s,eqv,2,1.0
/image,save,eqv1,bmp
/post26
/axlab,x,Time(s)
/axlab,y,Seqv(Mpa)
ansol,2,1,s,eqv,n1
ansol,3,2,s,eqv,n2
ansol,4,16,s,eqv,n16
ansol,5,23,s,eqv,n23
ansol,6,36,s,eqv,n36
xvar,1
plvar,2,3,4,5,6,
/image,save,eqv2,bmp
/axlab,x,time(s)
/axlab,y,temp(c)
nsol,7,1,temp,,n1
nsol,8,2,temp,,n2
nsol,9,23,temp,,n23
nsol,10,16,temp,,n16
nsol,11,36,temp,,n36
plvar,7,8,9,10,
/image,save,temp2,bmp
finish
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-10-4 09:36 , Processed in 0.027466 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表