新手请教关于一个热分析的问题
模型是两个粘在一起的体,上面的体有一个悬浮的横梁,下体温度固定为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
不过分析的时候出现了问题,不知道是不是分析选项出的问题,总是一求解就死机,过一会就自动退出。想了半天也找不到事哪里出错了。请高手指教。万分感谢!!
还有一个问题,我先设置的下体固定为100度,然后选中上体,选中所有节点,设置初始条件为25度。然后list了一下所有初始温度,发现有的节点是25,有的节点是0.不知道这种设置方法是不是有问题。 黏贴命令流的时候不要一下子全复制进去。将计算那块的单独复制进去,试试。 schuhi 发表于 2012-5-14 13:35 static/image/common/back.gif
黏贴命令流的时候不要一下子全复制进去。将计算那块的单独复制进去,试试。 ...
我试过了,先把命令流建模及以前的一起复制进去,然后手动mesh,手动施加载荷,手动设置分析选项,然后还是不行。刚才我试着只画出上面的一个体,进行了结构应力分析,发现那个模型应该是没问题的。
页:
[1]