【讨论】关于用plane13单元进行热结构耦合的困惑？？

hsg8888

我做了一个分析，出的问题感到很困惑。不知这是不是常识问题——耦合单元不能进行耦合分析的restart?!  耦合单元的restart只能做时间为1s的（并且是不可更改的?）?!
  在进行热结构耦合中，我用的是plane13耦合单元进行热应力耦合分析，需要变边界条件的。
  
  分析的对象是一个薄壁的塑料矩形框的冷却问题。第一次的分析中的矩形框的外侧四边固定（上下两个边的uy＝0，左右两个边的ux＝0）边界的对流为sfl,all,conv,2000,,15。第二次的分析antype，tran，rest 边界的对流为sfl,all,conv,3000,,15，去掉上一个分析中的四边固定，加上节点1，2的uy＝0。
如下面的程序所示，
  第一次分析时设定分析时间为12s，第二次的分析中我设定为0.2s，结果分析出13s秒的结果。即把antype，tran，rest的分析当作1s来进行了。如果把0.2s变成其他的时间值（如2s，3s等）的话，程序也把它当作1s来分析。我用的是ansys7.0。
  
后台窗口中的这句话不知是不是问题的关键？
*** NOTE ***                            CP=   1187.438   TIME= 08:07:32
  New Solution Control Option is only valid for a single field
  structural or thermal analysis.  SOLCONTROL,ON has been de-activated.
  
如果把第一次和第二次分析的cnvtol,f,,0.1,,0.1 去掉，在错误文件里出现类似
*** WARNING ***                         CP=    113.663   TIME= 16:19:38
  The calculated reference FORCE CONVERGENCE VALUE = 5.763800812E-03 is   
  less than the default minimum of 1 (MINREF field of CNVTOL command).     
  Check results carefully.                                                
的错误信息。
如果在第一次分析中不加上cnvtol,f,,0.1,,0.1 ，程序能够分析到12s，在错误文件里出现类似如上的信息。
  
下面几种情况都是不能改变restart分析后，只能得到12s+1s＝13s的结果（当restart分析的时间设定为不为1s时,第一次分析时间为12s）
再把后一部分的命令读入进行restart分析时（有cnvtol,f,,0.1,,0.1），
1.把前后两个分析的cnvtol,f,,0.1,,0.1去掉或去掉其中之一
2.把第二个分析中的tref，190去掉
3.上述两条的组合
4.改变划分网格的大小，即改变esize,0.0002
5.进行restart分析时，如果把antype,trans,rest改成antype,static,rest即改成稳态分析，把
time,0.2  
autots,-1
deltim,0.04,0.004,0.4,1
kbc,0
删掉，程序出现错误信息，最后算的结果也为12s+1s
6.下面的命令流不作改变时。
总之似乎进行rest分析时只能作1s的分析，得到12s+1s＝13s的结果。好郁闷啊!!
  
另类郁闷的情况
把下面的命令中的两个tref，190都去掉，把restart的分析时间改为60s时，结果只得到60s的分析，而不是应该的12s+60s＝72s
这样的分析显然不对了，因为参考温度不设定就默认成0度了。
  
那位大侠想试一试下面的程序时，不妨把网格划分粗糙一点，第一次分析的时间设短一点。我想结果也是一样的。
这样的问题的确让我郁闷，百思不得其解。恳请高手指点为感!!!
  
/units,si  
/PREP7   
ET,1,PLANE13  ！定义单元类型
MP,DENS,1,1400  ！定义型材的密度
MPTEMP,,-3,27,47,67  ！定义随温度变化的型材的导热性
MPDATA,KXX,1,,0.167,0.170,0.171,0.172
MPTEMP  ！清除温度场
MPTEMP,,2,23,87,107,147  ！定义随温度变化的型材的比热
MPDATA,C,1,,865,960,1475,1569,1741
MPTEMP ！清除原来的温度场  
MPTEMP,,20,60,80,100,140,180
MPDATA,EX,1,,2.18E+9,1.30E+9,0.86E+8,7.2E6 ,4.22E6 ,8.31E5  !定义随温度变化的弹性模量
MP,PRXY,1,0.3  !定义泊松比
MPTEMP ！清除原来的温度场  
MP,ALPX,1,1e-4  !定义热膨胀系数
!这里定义了型材的各项参数，采用耦合单元PLANE13
blc4,0,0,0.06,0.06
blc4,0.002,0.002,0.056,0.056
aovlap,all
adele,2
esize,0.0002
amesh,all
finish
/solu
antype,trans,new
cnvtol,f,,0.1,,0.1 !设定以力作为收敛标准，不设定时在错误文件中会出现相关的错误信息。
tunif,190
tref,190
lsel,s,line,,1,4,1
sfl,all,conv,2000,,15
allsel,all
lsel,s,line,,1,3,2
dl,all,,uy,0
allsel,all
lsel,s,line,,2,4,2
dl,all,,ux,0
allsel,all  
timint,on
time,12
autots,-1
deltim,0.04,0.004,0.4,1
kbc,0
solve
finish
!为接着计算的程序
/solu
antype,trans,rest
cnvtol,f,,0.1,,0.1
tref,190
lsel,s,line,,1,4,1
sfldele,all,conv
allsel,all
lsel,s,line,,1,3,2
dldele,all,uy
allsel,all
lsel,s,line,,2,4,2
dldele,all,ux
allsel,all
nsel,s,node,,1
nsel,a,node,,2
d,all,uy,0
allsel,all
lsel,s,line,,1,4,1
sfl,all,conv,3000,,15
allsel,all
time,0.2         !可设定不为1s（如2s，3s或0.5s），但结果是一样的
autots,-1
deltim,0.04,0.004,0.4,1
kbc,0
solve
finish

ainite

可能你没理解重启动分析中 TIME 的含义，下面是略微修改过的命令流：
/PREP7  
ET,1,PLANE13 ！定义单元类型  
MP,DENS,1,1400 ！定义型材的密度  
MPTEMP,,-3,27,47,67 ！定义随温度变化的型材的导热性  
MPDATA,KXX,1,,0.167,0.170,0.171,0.172  
MPTEMP ！清除温度场  
MPTEMP,,2,23,87,107,147 ！定义随温度变化的型材的比热  
MPDATA,C,1,,865,960,1475,1569,1741  
MPTEMP ！清除原来的温度场  
MPTEMP,,20,60,80,100,140,180  
MPDATA,EX,1,,2.18E+9,1.30E+9,0.86E+8,7.2E6 ,4.22E6 ,8.31E5 !定义随温度变化的弹性模量  
MP,PRXY,1,0.3 !定义泊松比  
MPTEMP ！清除原来的温度场  
MP,ALPX,1,1e-4 !定义热膨胀系数  
!这里定义了型材的各项参数，采用耦合单元PLANE13  
K,1
K,2,0.002
K,3,0.058
K,4,0.06
K,5,,0.002
K,6,0.002,0.002
K,7,0.058,0.002
K,8,0.06,0.002
K,9,,0.058
K,10,0.002,0.058
K,11,0.058,0.058
K,12,0.06,0.058
K,13,,0.06
K,14,0.002,0.06
K,15,0.058,0.06
K,16,0.06,0.06
a,1,2,6,5
a,2,3,7,6
a,3,4,8,7
a,5,6,10,9
a,7,8,12,11
a,9,10,14,13
a,10,11,15,14
a,11,12,16,15
esize,0.001
MSHKEY,1
amesh,all  
finish  
  
/solu  
antype,trans,new  
trnopt,full        !full方法
cnvtol,f,,0.1,,0.1 !设定以力作为收敛标准，不设定时在错误文件中会出现相关的错误信息。  
tunif,190  
tref,190  
lsel,s,,,1  
lsel,a,,,5
lsel,a,,,8
lsel,a,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
sfl,all,conv,2000,,15  
allsel,all  
lsel,s,,,1  
lsel,a,,,5
lsel,a,,,8
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
dl,all,,uy,0  
allsel,all  
lsel,s,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
dl,all,,ux,0  
allsel,all  
timint,on  
time,12  
deltim,0.2,0.02,0.2,1
kbc,0  
solve  
finish  
  
!为接着计算的程序  
/solu  
antype,trans,rest  
trnopt,full
cnvtol,f,,0.1,,0.1  
tref,190  
lsel,s,,,1  
lsel,a,,,5
lsel,a,,,8
lsel,a,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
sfldele,all,conv
allsel,all  
lsel,s,,,1  
lsel,a,,,5
lsel,a,,,8
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
dldele,all,uy
allsel,all  
lsel,s,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
dldele,all,ux
allsel,all  
d,1,ux
d,1,uy
d,178,uy
lsel,s,,,1  
lsel,a,,,5
lsel,a,,,8
lsel,a,,,9
lsel,a,,,14
lsel,a,,,23
lsel,a,,,24
lsel,a,,,22
lsel,a,,,18
lsel,a,,,19
lsel,a,,,13
lsel,a,,,4
sfl,all,conv,3000,,15  
allsel,all  
time,12.2
deltim,0.05,0.005,0.05,1
kbc,0  
solve  
finish  
这里面对模型的修改只是为了分网的需要，时间的定义是主要修改之处。

hsg8888

非常感谢ainite的回复！！！
我明白了。
simwe的确是个好的网站，我每次上网都是必上的。论坛给了我很多的帮助。