有没有人用ANSYS做热屈曲分析的

超静定结构在温度载荷作用下可能发生失稳、热屈曲分析就是求似的结构失稳的临界载荷。目前，关于超静定杆的热屈曲分析理论和ANSYS计算都已完成，开始做板结构的热屈曲了，可是由于找到的理论方面的东西太少，不知道计算结果是否可信。想请问一下大家可有人做过这方面的东东，能否提供一些板屈曲临界温度理论计算的公式或者相关文献，或者帮我看一下我的计算过程和结果有没有问题：
以下为一块1000*1000*10的板，选择shell63单元，包括了四边简支、四边固支等六中支撑情况，命令流和结果如下：----------------------------------------------四边简支----------------------------------------------
finish
/clear,nostart
/prep7
!!!!!!!!输入材料参数
ET,1,shell63
r,1,10
mp,ex,1,2.1e5
mp,prxy,1,0.3
MP,ALPX,1,12e-6
TREF,0
RECTNG,-500,500,-500,500
LESIZE,all,50
MSHAPE,0,2D
MSHKEY,1
AMESH,all   
finish  
/SOLU   
ANTYPE,STATIC   
PSTRES,ON   

nsel,s,loc,y,-500
nsel,a,loc,y,500
d,all,ux
d,all,uy
d,all,uz
nsel,s,loc,x,-500
nsel,a,loc,x,500
d,all,ux
d,all,uy
d,all,uz

allsel
!BFUNIF,TEMP ,1
BFA,1,TEMP,1

nsel,all
outpr,,1
SOLVE
FINISH
/SOLU
ANTYPE,BUCKLE
BUCOPT,LANB,10
MXPAND,10,,,yes
SOLVE
FINISH

/post1
! 设置底色为白色
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  
/REPLOT
!set,list
!set,first
pldisp,0
*get,fcr,mode,1,freq
*status,parm


----------------------------------------------四边固支----------------------------------------------
finish
/clear,nostart
/prep7

!!!!!!!!输入材料参数
ET,1,shell63
r,1,10
mp,ex,1,2.1e5
mp,prxy,1,0.3
MP,ALPX,1,12e-6
RECTNG,-500,500,-500,500
LESIZE,all,50
MSHAPE,0,2D
MSHKEY,1
AMESH,all   
finish  
/SOLU   
ANTYPE,STATIC   
PSTRES,ON   

nsel,s,loc,y,-500
nsel,a,loc,y,500
d,all,all
nsel,s,loc,x,-500
nsel,a,loc,x,500
d,all,all


BFA,1,TEMP,1

nsel,all
outpr,,1
SOLVE
FINISH
/SOLU
ANTYPE,BUCKLE
BUCOPT,LANB,10
MXPAND,10,,,yes
SOLVE
FINISH

/post1
! 设置底色为白色
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  
/REPLOT
!set,list
!set,first
pldisp,0
*get,fcr,mode,1,freq
*status,parm


----------------------------------------------b简支\a固支----------------------------------------------
finish
/clear,nostart
/prep7

!!!!!!!!输入材料参数
ET,1,shell63
r,1,10
mp,ex,1,2.1e5
mp,prxy,1,0.3
MP,ALPX,1,12e-6
TREF,0
RECTNG,-500,500,-500,500
LESIZE,all,50
MSHAPE,0,2D
MSHKEY,1
AMESH,all   
finish  
/SOLU   
ANTYPE,STATIC   
PSTRES,ON   

nsel,s,loc,y,-500
nsel,a,loc,y,500
d,all,all
nsel,s,loc,x,-500
nsel,a,loc,x,500
d,all,ux
d,all,uy
d,all,uz

BFA,1,TEMP,1

nsel,all
outpr,,1
SOLVE
FINISH
/SOLU
ANTYPE,BUCKLE
BUCOPT,LANB,10
MXPAND,10,,,yes
SOLVE
FINISH

/post1
! 设置底色为白色
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  
/REPLOT
pldisp,0
*get,fcr,mode,1,freq
*status,parm


----------------------------------------------b简支\a一简支、一自由----------------------------------------------
finish
/clear,nostart
/prep7

!!!!!!!!输入材料参数
ET,1,shell63
r,1,10
mp,ex,1,2.1e5
mp,prxy,1,0.3
MP,ALPX,1,12e-6
TREF,0
RECTNG,-500,500,-500,500
LESIZE,all,50
MSHAPE,0,2D
MSHKEY,1
AMESH,all   
finish  
/SOLU   
ANTYPE,STATIC   
PSTRES,ON   

nsel,s,loc,y,-500
d,all,ux
d,all,uy
d,all,uz
nsel,s,loc,x,-500
nsel,a,loc,x,500
d,all,ux
d,all,uy
d,all,uz

BFA,1,TEMP,1

nsel,all
outpr,,1
SOLVE
FINISH
/SOLU
ANTYPE,BUCKLE
BUCOPT,LANB,10
MXPAND,10,,,yes
SOLVE
FINISH

/post1
! 设置底色为白色
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  
/REPLOT
set,list
set,first
pldisp,0
*get,fcr,mode,1,freq
*status,parm


----------------------------------------------b简支\a一固支、一自由----------------------------------------------
finish
/clear,nostart
/prep7

!!!!!!!!输入材料参数
ET,1,shell63
r,1,10
mp,ex,1,2.1e5
mp,prxy,1,0.3
MP,ALPX,1,12e-6
TREF,0
RECTNG,-500,500,-500,500
LESIZE,all,50
MSHAPE,0,2D
MSHKEY,1
AMESH,all   
finish  
/SOLU   
ANTYPE,STATIC   
PSTRES,ON   

nsel,s,loc,y,-500
d,all,all
nsel,s,loc,x,-500
nsel,a,loc,x,500
d,all,ux
d,all,uy
d,all,uz

BFA,1,TEMP,1

nsel,all
outpr,,1
SOLVE
FINISH
/SOLU
ANTYPE,BUCKLE
BUCOPT,LANB,10
MXPAND,10,,,yes
SOLVE
FINISH

/post1
! 设置底色为白色
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  
/REPLOT
set,list
set,first
pldisp,0
*get,fcr,mode,1,freq
*status,parm


----------------------------------------------b固支\a简支----------------------------------------------
finish
/clear,nostart
/prep7

!!!!!!!!输入材料参数
ET,1,shell63
r,1,10
mp,ex,1,2.1e5
mp,prxy,1,0.3
MP,ALPX,1,12e-6
TREF,0
RECTNG,-500,500,-500,500
LESIZE,all,50
MSHAPE,0,2D
MSHKEY,1
AMESH,all   
finish  
/SOLU   
ANTYPE,STATIC   
PSTRES,ON   

nsel,s,loc,y,-500
nsel,a,loc,y,500
d,all,ux
d,all,uy
d,all,uz
nsel,s,loc,x,-500
nsel,a,loc,x,500
d,all,all

BFA,1,TEMP,1

nsel,all
outpr,,1
SOLVE
FINISH
/SOLU
ANTYPE,BUCKLE
BUCOPT,LANB,10
MXPAND,10,,,yes
SOLVE
FINISH

/post1
! 设置底色为白色
/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  
/REPLOT
!set,list
!set,first
pldisp,0
*get,fcr,mode,1,freq
*status,parm


结果：
四边简支                           10.5289398
四边固定                           27.8712380
b简支，a固定                20.1417103
b简支，a一简支、一自由        8.35685414
b简支，a一固定、一自由        10.1283455
b固定，a简支                20.1417103