新手请教关于一个热分析的问题

junnerf

模型是两个粘在一起的体，上面的体有一个悬浮的横梁，下体温度固定为100，上面初始温度25.由于膨胀系数不一样，导致上面的体有变形，要求此形变对横梁的影响很小。我的命令流是下面这样子的：
/FILENAME,PiezoTransducer
/TITLE,Struct-Thermal Analysis
/UNITS,SI          !国际单位制
*set,length,0.17e-3 !定义梁的长度
*set,width,0.04e-3 !定义梁的宽度
!定义各层距底面的厚度
*set,thickness1,0.45e-3 !Layer 1
*set,thickness2,0.44e-3 !Layer 2
*set,thickness3,0.43e-3 !Layer 3
*set,thickness4,0.39e-3 !Layer 4
*set,thickness5,0.10e-3 !Layer 5

/PREP7
ET,1,SOLID90
KEYOPT,1,1,1
DOF,UX,UY,UZ
MP,DENS,1,2330  !定义材料1（Si）的材料属性：密度
MP,EX,1,130E9  !弹性模量
MP,C,1,700  !比热
MP,PRXY,1,0.278  !泊松比
MP,ALPX,1,2.6E-6 !线膨胀系数
MP,KXX,1,148  !热传导率
MP,DENS,2,7850  !定义材料2（Cu）的材料属性
MP,EX,2,105E9
MP,C,2,386
MP,PRXY,2,0.31
MP,ALPX,2,17.5E-6
MP,KXX,2,398
!创建几何模型
BLOCK,0.5e-3,-0.5e-3,0.5e-3,-0.5e-3,0,0.45e-3    !Volume1

!Keypoints in Layer 1(Volume1 includes 8 Kpts)
k,9,-0.245e-3,length/2,thickness1
k,10,-0.245e-3,0.104e-3,thickness1
K,11,,0.245e-3,thickness1
k,12,0.245e-3,0.104e-3,thickness1
k,13,0.245e-3,-0.104e-3,thickness1
k,14,,-0.245e-3,thickness1
k,15,-0.245e-3,-0.104e-3,thickness1
k,16,-0.245e-3,-length/2,thickness1
k,17,-0.245e-3-width,-length/2,thickness1
k,18,-0.245e-3-width,-0.104e-3-width*0.577,thickness1
k,19,,-0.245e-3-width*1.155,thickness1
k,20,0.245e-3+width,-0.104e-3-width*0.577,thickness1
k,21,0.245e-3+width,0.104e-3+width*0.577,thickness1
k,22,,0.245e-3+width*1.155,thickness1
k,23,-0.245e-3-width,0.104e-3+width*0.577,thickness1
k,24,-0.245e-3-width,length/2,thickness1
!Keypoints in Layer 4
k,25,-0.245e-3,length/2,thickness4
k,26,-0.245e-3,0.104e-3,thickness4
K,27,,0.245e-3,thickness4
k,28,0.245e-3,0.104e-3,thickness4
k,29,0.245e-3,-0.104e-3,thickness4
k,30,,-0.245e-3,thickness4
k,31,-0.245e-3,-0.104e-3,thickness4
k,32,-0.245e-3,-length/2,thickness4
k,33,-0.245e-3-width,-length/2,thickness4
k,34,-0.245e-3-width,length/2,thickness4
!Keypoints in Layer 5
k,35,-0.245e-3-width,-0.104e-3-width*0.577,thickness5
k,36,,-0.245e-3-width*1.155,thickness5
k,37,0.245e-3+width,-0.104e-3-width*0.577,thickness5
k,38,0.245e-3+width,0.104e-3+width*0.577,thickness5
k,39,,0.245e-3+width*1.155,thickness5
k,40,-0.245e-3-width,0.104e-3+width*0.577,thickness5

!Keypoints in Layer 2
k,41,-0.175e-3,0.028e-3,thickness2
k,42,,0.13e-3,thickness2
k,43,0.175e-3,0.028e-3,thickness2
k,44,0.175e-3,-0.028e-3,thickness2
k,45,,-0.13e-3,thickness2
k,46,-0.175e-3,-0.028e-3,thickness2

!Keypoints in Layer 3
k,47,-0.175e-3,0.028e-3,thickness3
k,48,,0.13e-3,thickness3
k,49,0.175e-3,0.028e-3,thickness3
k,50,0.175e-3,-0.028e-3,thickness3
k,51,,-0.13e-3,thickness3
k,52,-0.175e-3,-0.028e-3,thickness3

V,41,42,45,46,47,48,51,52 !Volume 2
V,42,43,44,45,48,49,50,51 !Volume 3
VADD,2,3   !Volume 4,即空腔
V,10,11,14,15,26,27,30,31 !Volume 2
V,11,12,13,14,27,28,29,30 !Volume 3
VADD,2,3   !Volume 5
VSBV,5,4   !Volume 2
V,24,9,16,17,34,25,32,33 !Volume 3,即梁
VADD,2,3   !Volume 4
V,23,22,19,18,40,39,36,35 !Volume 2
V,22,21,20,19,39,38,37,36 !Volume 3
VADD,2,3   !Volume 5
VSBV,1,5   !Volume 2
VADD,2,4   !Volume 1,即所需的上体

BLOCK,0.5e-3,-0.5e-3,0.5e-3,-0.5e-3,0,-0.4e-3    !Volume2
VGLUE,1,2   !下面的体编号变成了3，不知道为什么
VSEL,S,,,3
VATT,2,,1,0
VSEL,S,,,1
VATT,1,,1,0
VSEL,A,,,3
VPLOT,1,3
VSEL,s,,,1,3
SMRTSIZE,4
mshkey,0  !自由网格
mshape,1,3D  !三角网格
VMESH,ALL
FINISH
/SOL
ANTYPE,4
TRNOPT,FULL
LUMPM,0
!给底面增加位移限制
DA,22,UX,0
DA,22,UY,0
DA,22,UZ,0
!底座温度限制为100
VSEL,S,,,3
BFV,ALL,TEMP,100
!上面主体初始温度
VSEL,S,,,1
NSLV,S,0
IC,ALL,TEMP,25
vsel,all
vplot,all
!设置分析选项
TIME,300
AUTOTS,-1
DELTIM,10,1,20,1
KBC,1   !Stepped
OUTRES,ALL,ALL
SOLVE
FINISH

不过分析的时候出现了问题，不知道是不是分析选项出的问题，总是一求解就死机，过一会就自动退出。想了半天也找不到事哪里出错了。请高手指教。万分感谢！！

junnerf

还有一个问题，我先设置的下体固定为100度，然后选中上体，选中所有节点，设置初始条件为25度。然后list了一下所有初始温度，发现有的节点是25，有的节点是0.不知道这种设置方法是不是有问题。

schuhi

黏贴命令流的时候不要一下子全复制进去。将计算那块的单独复制进去，试试。

junnerf

schuhi 发表于 2012-5-14 13:35
黏贴命令流的时候不要一下子全复制进去。将计算那块的单独复制进去，试试。 ...

我试过了，先把命令流建模及以前的一起复制进去，然后手动mesh，手动施加载荷，手动设置分析选项，然后还是不行。刚才我试着只画出上面的一个体，进行了结构应力分析，发现那个模型应该是没问题的。