求高手帮助看一下程序
本帖最后由 热爱 于 2012-3-20 13:10 编辑/PREP7
BLOCK,0,240,0,220,0,0.08, !定义集合体
BLOCK,0,240,0,220,0.08,1.42,
BLOCK,0,240,0,220,1.42,1.50,
BLOCK,50,90,130,170,1.50,1.65,
BLOCK,50,90,130,170,1.65,2.92,
NUMMRG,ALL, , , ,LOW !压缩编号
NUMCMP,ALL
TOFFST,273
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,240
MPDATA,KXX,1,,521.5
MPDATA,KXX,1,,560
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,8893
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,240
MPDATA,C,1,,356.8
MPDATA,C,1,,425.3
MPTEMP,1,0
MPDATA,KXX,2,,0.29
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,1859
MPTEMP,,,,,,,,
MPTEMP,1,30
MPTEMP,2,70
MPTEMP,3,120
MPTEMP,4,240
MPTEMP,5,250
MPDATA,C,2,,1200
MPDATA,C,2,,1380
MPDATA,C,2,,1500
MPDATA,C,2,,1650
MPDATA,C,2,,1600
MPTEMP,1,0
MPDATA,KXX,3,,50
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,1800
MPTEMP,,,,,,,,
MPTEMP,1,40
MPTEMP,2,80
MPTEMP,3,160
MPTEMP,4,170
MPTEMP,5,220
MPTEMP,6,225
MPTEMP,7,230
MPTEMP,8,250
MPDATA,C,3,,740
MPDATA,C,3,,850
MPDATA,C,3,,880
MPDATA,C,3,,960
MPDATA,C,3,,970
MPDATA,C,3,,1050
MPDATA,C,3,,920
MPDATA,C,3,,900
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,4,,50
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,4,,8350
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,4,,200 !材料参数定义
ET,1,SOLID70
ET,2,SOLID90 !单元选择
VSEL,,,,1
VATT,1,,1,0
vSEL,,,,2
VATT,2,,1,0
VSEL,,,,3
VATT,1,,1,0
VSEL,,,,4
VATT,4,,2,0
VSEL,,,,5
VATT,3,,1,0 !给不同实体添加材料属性
alls
plot
ESIZE,0,24,
vmesh,1
vmesh,2
vmesh,3
esize,0,10
vmesh,4
vmesh,5 !网格划分 出现警告
/sol
tunif,20, !初始温度20
antype,4 !设置为瞬态求解
x2=240
x3=300
tim=0
x1=0
timint,on !打开瞬态响应
*do,i,1,18 !载荷移动循环9次
autots,1
nsubst,30,,,1
kbc,1
tim=tim+3
time,tim !设置载荷步结束时间
x1=i*19.5
*if,x1,le,50,then
nsel,s,loc,z,1.5
nsel,r,loc,x,0,x1 !再选择进入炉区部分的节点
sf,all,conv,x3,x2 !在所选节点上施加对流传热系数
nsel,s,loc,z,0
nsel,r,loc,x,0,x1
sf,all,conv,x3,x2
alls !选择所有实体
!设置子步数
outres,all,all !设置输出
solve !求解运算
sfadele,all,conv !删除上述施加的载荷
*elseif,x1,le,90,then
nsel,s,loc,x,50
nsel,r,loc,y,130,170
nsel,r,loc,z,1.5,2.92
sf,all,conv,x3,x2
nsel,s,loc,y,130
nsel,r,loc,x,50,x1
nsel,r,loc,z,1.5,2.92
sf,all,conv,x3,x2
nsel,s,loc,y,170
nsel,r,loc,x,50,x1
nsel,r,loc,z,1.5,2.92
sf,all,conv,x3,x2
nsel,s,loc,z,1.5
nsel,r,loc,x,0,x1
nsel,r,loc,y,0,130
sf,all,conv,x3,x2
nsel,s,loc,z,1.5
nsel,r,loc,x,0,x1
nsel,r,loc,y,170,220
sf,all,conv,x3,x2
nsel,s,loc,z,1.5
nsel,r,loc,x,0,50
nsel,r,loc,y,130,170
sf,all,conv,x3,x2
nsel,s,loc,z,2.92
nsel,r,loc,x,50,x1
nsel,r,loc,y,130,170
sf,all,conv,x3,x2
nsel,s,loc,z,0
nsel,r,loc,x,0,x1
sf,all,conv,x3,x2
alls !选择所有实体
!设置子步数
outres,all,all !设置输出
solve !求解运算
sfadele,all,conv !删除上述施加的载荷
*else
nsel,s,loc,x,50
nsel,r,loc,y,130,170
nsel,r,loc,z,1.5,2.92
sf,all,conv,x3,x2
nsel,s,loc,y,130
nsel,r,loc,x,50,90
nsel,r,loc,z,1.5,2.92
sf,all,conv,x3,x2
nsel,s,loc,y,170
nsel,r,loc,x,50,90
nsel,r,loc,z,1.5,2.92
sf,all,conv,x3,x2
nsel,s,loc,x,90
nsel,r,loc,y,130,170
nsel,r,loc,z,1.5,2.92
sf,all,conv,x3,x2
nsel,s,loc,z,1.5
nsel,r,loc,x,0,x1
nsel,r,loc,y,0,130
sf,all,conv,x3,x2
nsel,s,loc,z,1.5
nsel,r,loc,x,0,x1
nsel,r,loc,y,170,220
sf,all,conv,x3,x2
nsel,s,loc,z,1.5
nsel,r,loc,x,0,50
nsel,r,loc,y,130,170
sf,all,conv,x3,x2
nsel,s,loc,z,2.92
nsel,r,loc,x,50,90
nsel,r,loc,y,130,170
sf,all,conv,x3,x2
nsel,s,loc,z,1.5
nsel,r,loc,x,90,x1
nsel,r,loc,y,130,170
sf,all,conv,x3,x2
nsel,s,loc,z,0
nsel,r,loc,x,0,x1
sf,all,conv,x3,x2
alls !选择所有实体
!设置子步数
outres,all,all !设置输出
solve !求解运算
sfadele,all,conv !删除上述施加的载荷
*endif
*enddo !循环结束
finish
现在遇到问题,求解结果不正确,不知道是否程序处问题,还请高手指点!扣扣:17 4 8 7 0 7 9 7.
页:
[1]