下天雄 发表于 2010-5-7 14:15:16

请教一个各向异性材料瞬态计算问题

本帖最后由 下天雄 于 2010-5-7 17:07 编辑

一个圆柱体，导热系数是各向异性的
向其中一个单元加一个热量，看整体温度场变化情况

问题是：计算开始后，其余单元的温度会降低？
不知道问题出在哪里？
以下为命令流，哪位大侠有时间帮忙看看，命令流问题出在哪里？
谢谢大家了。

finish
/clear
/COM,Thermal
/CONFIG,FSPLIT,750
inrad=927.3e-3/2         
outrad=1022.28e-3/2      
coilh=(1.32e-3+0.01e-3)*268         
   
!/NUMBER,-1
/UIS,MSGPOP,3
/PREP7
CSWPLA,11,1,1,1,   !cyindrical local coordination is very important
ESYS,      11
ET,1,solid70
cylind,inrad,outrad,-coilh/2,coilh/2,
vsel,s,,,1
allsel,below,volu
lsel,r,length,,coilh
lesize,all,,,50
lsel,s,radius,,inrad,outrad
lesize,all,,,50
vatt,5
vsweep,all

!材料物性
MPTEMP,,,,,,,,
MPTEMP,1,4,6,8,10,15,20,
mptemp,7,25,30,35,40,50,60
mptemp,13,70,76,80,90,100,120
mptemp,19,140,160,180,200,250,300
MPDATA,KXX,5,1,0.1056,0.1502,0.1941,0.2353,0.3224,0.3902
MPDATA,KXX,5,7,0.4454,0.4929,0.5355,0.5747,0.6463,0.7119
MPDATA,KXX,5,13,0.7729,0.8076,0.8300,0.8838,0.9347,1.0287
MPDATA,KXX,5,19,1.1140,1.1922,1.2645,1.3321,1.4835,1.6141
MPDATA,KYY,5,1,240,370,470,600,850,1100
MPDATA,KYY,5,7,1200,1200,1100,1000,770,620
MPDATA,KYY,5,13,550,520,490,470,450,430
MPDATA,KYY,5,19,420,410,400,400,400,400
MPDATA,KZZ,5,1,0.1056,0.1502,0.1941,0.2353,0.3224,0.3902
MPDATA,KZZ,5,7,0.4454,0.4929,0.5355,0.5747,0.6463,0.7119
MPDATA,KZZ,5,13,0.7729,0.8076,0.8300,0.8838,0.9347,1.0287
MPDATA,KZZ,5,19,1.1140,1.1922,1.2645,1.3321,1.4835,1.6141
MPTEMP,,,,,,,,                     !Cu specific heat
MPTEMP,1,4,6,8,10,15,20,
mptemp,7,25,30,35,40,50,60
mptemp,13,70,76,80,90,100,120
mptemp,19,140,160,180,200,250,300
MPDATA,C,5,1,0.09944,0.2304,0.4642,0.8566,2.9627,7.5061
MPDATA,C,5,7,15.2593,26.4739,40.8794,57.8445,96.2689,135.8794
MPDATA,C,5,13,172.763,192.7603,205.1023,232.5038,255.3289,289.9726
MPDATA,C,5,19,314.4069,332.6329,347.0727,359.0114,380.6159,389.4015

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,5,,8900
allsel
save
/solu
ANTYPE,4
solcon,on   
NEQIT,100   
LNSRCH,1   
timestep=0.01
DELTIM,timestep,timestep,timestep   
TIME,timestep   
KBC,1         
hgent=1e9                           
esel,s,,,1923
BFe,all,HGEN,1,hgent
allsel,all
IC,all,TEMP,4.5,   
allsel,all
solve
/POST1
SET,LAST
PLNSOL, TEMP,, 0

下天雄 发表于 2010-5-7 15:21:30

急求啊！！！
程序卡在这里，进行不下去了
大家帮帮忙啊
1# 下天雄

下天雄 发表于 2010-5-8 09:44:52

没人愿意帮忙看看嘛？

下天雄 发表于 2010-5-10 09:04:15

大家帮帮忙啊
页: [1]
查看完整版本: 请教一个各向异性材料瞬态计算问题