热-结构直接耦合,接触传导热的问题
做了一个热-结构直接耦合瞬态分析,问题是这样的:两个相接触的物体,一个小,一个大,大的不动,小的在大的表面上以一定的速度移动,同时这两个物体温度不同,小的温度高,小物体可以通过接触关系向大物体传热。现在有两个问题:1.应用D命令,给物体施加移动速度vely,但物体始终不动,不知是为什么?2.帮助中说,位移命令和速度命令不能同时应用,如果同时应用,后面的将覆盖前面的命令。可是我一方面要让小物体以一定的速度向Y方向移动,另一方面要将小物体X、Z方向的自由度约束,该怎么实现呢?请各位大侠指教!
命令流如下:
fini
/clear
/filename,hotting-little
/congif,nres,1e10
blong=0.04
!接触带长
rblong=0.5
!左板长
lblong=0.5
!右板长
bwide=0.5
!板宽
bthick=0.03
!板厚
hdiameter=0.02
!火焰直径
coefficient=10000 !接触传热系数
/prep7
block,-blong/2,blong/2,-bwide/2,bwide/2,0,-bthick !建接触带
CYL4, , ,hdiameter/2, , , ,0.05
!建火焰
VGEN, ,2, , , ,bwide/2-hdiameter/2, , , ,1
!移动火焰
block,-rblong/2,rblong/2,-bwide/2,bwide/2,0,-bthick !建左板子
block,-lblong/2,lblong/2,-bwide/2,bwide/2,0,-bthick !建右板子
VGEN, ,3, , ,-rblong/2-blong/2,, , , ,1
!移动左板
VGEN, ,4, , ,lblong/2+blong/2,, , , ,1
!移动右板
VGLUE,1,3,4
ET,1,SOLID5
mptemp,1,0,200,400,800,1200
mpdata,kxx,1,1,59,54,45,28,30
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,800
MPTEMP,5,1200
MPDATA,C,1,,469
MPDATA,C,1,,519
MPDATA,C,1,,594
MPDATA,C,1,,963
MPDATA,C,1,,661
mp,dens,1,7860
MP,EX,1,1.44e10
MP,PRXY,1,0.3
TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,400e6,1e6,,,, !屈服应力,切向模量
LSEL,s, , ,23
LSEL,a, , ,37
LESIZE,all,0.01, , , , , , ,1
LSEL,s, , ,47,48
LSEL,a, , ,51,52
LESIZE,all,0.02, , , , , , ,1
LSEL,s, , ,31,34,3
LSEL,a, , ,44,45
LESIZE,all,0.01, , , , , , ,1
TYPE,
1
MAT,
1
VSWEEP,5
VSWEEP,6
LSEL,s, , ,4,5
LSEL,a, , ,2
LSEL,a, , ,7
LESIZE,all,0.01, , , , , , ,1
LSEL,s, , ,1,3,2
LSEL,a, , ,6,8,2
LESIZE,all,0.02, , , , , , ,1
TYPE,
1
MAT,
1
VSWEEP,1
LSEL,s, , ,17,20
LESIZE,all,0.004, , , , , , ,1
LSEL,s, , ,22
LESIZE,all,0.005, , , , , , ,1
TYPE,
1
MAT,
1
VSWEEP,2
!钢板的温度为20度
vsel,s,,,1
vsel,a,,,5,6
eslv,s
nsle,s
IC,all,Temp,20
!Specifies initial conditions(Temp) at V1 nodes.
!约束火焰的X、Z向;Y向施加速度
vsel,s,,,2
eslv,s
nsle,s
D,all,ux,0
D,all,uz,0
D,all,vely,0.00630
!火焰的温度为3000度
vsel,s,,,2
eslv,s
nsle,s
D,all, ,3000, , , ,TEMP, , , , ,
!约束钢板
vsel,s,,,1
vsel,a,,,5,6
eslv,s
nsle,s
D,all,ux,0
D,all,uy,0
D,all,uz,0
alls
!创建传热接触模型
MP,MU,1,
MAT,1
MP,EMIS,1,
R,3
REAL,3
ET,2,170
ET,3,174
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,10000,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,10,2
KEYOPT,3,1,3
type,2
MAT,1
ASEL,S,,,7
nsla,s,1
ESLN,S,0
esurf
type,3
MAT,1
ASEL,S,,,2
NSLA,S,1
ESLN,S,0
esurf
/solu
alls
tref,20
antype,TRANS
time,10
lnsrch,auto
autots,on
kbc,1
outres,all,all
deltime,0.1,0.01,1
solve
save