请教：为什么铸造温度分析中，温度加载顺序不同，结果差别这么大？

dreamlee

我做的铸造过程温度场分析，加初始温度的时候，先加高温（铸件部分）后加低温（砂型部分）和先加低温（砂型部分）后加高温（铸件部分）的两种计算结果完全不同。
原以为是自己加载错了，可是将ANSYS自带的程序换个加载顺序运行一遍，也是得到两种截然不同的结果

dreamlee

以下两段命令流除了加载初始温度的部分不同外，其余相同，请大家帮我分析一下，这样我自己的那个复杂的铸造分析估计就有思路了。
谢谢！

先加高温后加低温（照抄ANSYS里的例子）：
/TITLE,CASTING SOLIDIFICATION !Give the analysis a title      
/PREP7
K,1,0,0,0
K,2,22,0,0
K,3,10,12,0
K,4,0,12,0
/TRIAD,OFF                  
/REPLOT
!
A,1,2,3,4                    
SAVE
RECTNG,4,22,4,8              
APLOT                        
AOVLAP,1,2                  
ADELE,3,,,1                  
SAVE
!
MP,DENS,1,0.054              
MP,KXX,1,0.025
MP,C,1,0.28
!
MPTEMP,1,0,2643,2750,2875,,,
MPDATA,KXX,2,1,1.44,1.54,1.22,1.22,,,
MPDATA,ENTH,2,1,0,128.1,163.8,174.2
MPPLOT,KXX,2,,,,,            
MPPLOT,ENTH,2,,,,,           
SAVE
!
ET,1,PLANE55                 
!
SAVE
SMRT,5                       
MSHAPE,0,2D                  
MSHKEY,0                     
AMESH,5                     
!
TYPE,1                       
MAT,2                        
REAL                        
ESYS,0                       
AMESH,4                     
!
SAVE
SFL,1,CONV,0.014,,80,,      
SFL,3,CONV,0.014,,80,,
SFL,4,CONV,0.014,,80,,
SAVE
FINISH
/SOLU
!
ANTYPE,4                     
SOLCONTROL,ON,0              
!
APLOT
ASEL,S,,,4                  
NSLA,S,1                     
NPLOT                        
IC,ALL,TEMP,2875            
NSEL,INVE                    
/REPLOT                     
IC,ALL,TEMP,80               
ALLSEL,ALL                  
SAVE
!
TIME,3                       
AUTOTS,-1                    
DELTIM,0.01,0.001,0.25,1     
KBC,1                        
!
OUTRES,ALL,ALL               
SAVE
/STAT,SOLU                  
/REPLOT                     
APLOT                        
SOLVE
FINISH

先加低温后加高温（改了一下ANSYS里例子的初始温度加载顺序）：
/TITLE,CASTING SOLIDIFICATION !Give the analysis a title
/PREP7  
K,1,0,0,0   
K,2,22,0,0  
K,3,10,12,0
K,4,0,12,0  
/TRIAD,OFF                  
/REPLOT
!   
A,1,2,3,4                    
SAVE   
RECTNG,4,22,4,8              
APLOT                        
AOVLAP,1,2                  
ADELE,3,,,1                  
SAVE  
!   
MP,DENS,1,0.054              
MP,KXX,1,0.025  
MP,C,1,0.28
!   
MPTEMP,1,0,2643,2750,2875,,,   
MPDATA,KXX,2,1,1.44,1.54,1.22,1.22,,,   
MPDATA,ENTH,2,1,0,128.1,163.8,174.2
MPPLOT,KXX,2,,,,,            
MPPLOT,ENTH,2,,,,,           
SAVE   
!   
ET,1,PLANE55                 
!   
SAVE   
SMRT,5                        
MSHAPE,0,2D                  
MSHKEY,0                     
AMESH,5                     
!   
TYPE,1                          
MAT,2                           
REAL                        
ESYS,0                          
AMESH,4                       
!   
SAVE   
SFL,1,CONV,0.014,,80,,         
SFL,3,CONV,0.014,,80,,  
SFL,4,CONV,0.014,,80,,  
SAVE   
FINISH  
/SOLU   
!   
ANTYPE,4                     
SOLCONTROL,ON,0
!
APLOT   
ASEL,S,,,5  
NSLA,S,1                     
NPLOT   
IC,ALL,TEMP,80  
NSEL,INVE                    
/REPLOT
IC,ALL,TEMP,2875ALLSEL,ALL                  
SAVE
!   
TIME,3                       
AUTOTS,-1                    
DELTIM,0.01,0.001,0.25,1     
KBC,1                        
!   
OUTRES,ALL,ALL               
SAVE  
/STAT,SOLU                  
/REPLOT                     
APLOT                        
SOLVE   
FINISH  

不好意思，为了把问题说清楚，有点啰嗦，彩色部分就是改变的内容。望大家不吝赐教！

dreamlee

急呀！
大家帮帮忙吧，先谢啦！

nwpuyl

这个问题不难解决，看看如下两段命令流等价吗？（不等价）
APLOT
ASEL,S,,,4
NSLA,S,1
NPLOT
IC,ALL,TEMP,2875
NSEL,INVE
/REPLOT
IC,ALL,TEMP,80
虽然不是很明显，但是这种错误是必须注意的，小模型好查错，打模型就难了！
APLOT
ASEL,S,,,5
NSLA,S,1
NPLOT
IC,ALL,TEMP,80
NSEL,INVE
/REPLOT
IC,ALL,TEMP,2875ALLSEL,ALL
看出来了吧，主要在于NSLA,S,1这一句，它决定了两个面公用边线上的节点施加80还是2875的差别，第一段命令表示这些节点初始温度为2875，而第二段里就成了80，那么当然不等价了，结果怎么能期望相同呢！

dreamlee

谢谢斑竹
您说的这个问题我以前倒是注意到了，公共面上的节点温度应该就是后加的那个温度值
再请问一下，哪一种加载顺序比较合理呢？过程就是将熔化了的铝水浇铸到砂型里冷却。
谢谢！