热分析*do循环中重启动出现问题

xzyzqs

研究地源热泵，间歇运行温度场的时空变化。每天先运行9h（换热管内壁热流=500），然后停运15h(换热管内壁热流=0)，运行30天（如图所示）。
这是一个载荷随时间变化的热分析问题，而且有规律。
我的做法是采用重启动分析，然后有规律的部分采用*do循环来实现。

一、部分命令流
。。。前处理省略。。。
载荷施加和求解控制如下：

/config,nres,50000    !VALUE is maximum number of results sets (substeps) allowed on the result file     
!第一天（施加初始地温场和前9h小时热流）
/SOLU                 !进入求解器
ANTYPE,TRANS          !指定分析类型p235
TIMINT,ON,THERM       !打开瞬态效应P245
TIME,9*3600            !定义计算终止时间
deltim,3600,3600,3600      !一小时计算一次
AUTOTS,ON             !打开自动时间步长
OUTRES,,ALL           !写入每个子步
bfunif,temp,19        !施加均匀温度荷载
KBC,1                 !指定载荷步阶跃方式
LSEL,S,,,9,56,1
cm,xzy,line
sf,xzy,hflux,156.1    !施加热流荷载
ALLSEL               
SOLVE
FINISH
!剩下15小时停止运行
/SOLU                 !进入求解器
ANTYPE,TRANS,rest          !指定分析类型
TIMINT,ON,THERM       !打开瞬态效应
TIME,24*3600            !定义计算终止时间
deltim,3600,3600,3600      !定义初始荷载步
AUTOTS,ON             !打开自动时间步长
OUTRES,,ALL           !写入每个子步
KBC,1                 !指定载荷步阶跃方式
LSEL,S,,,9,56,1
cm,xzy,line
sf,xzy,hflux,0    !施加热流荷载
ALLSEL               
SOLVE
FINISH

！从第2天开始采用循环计算
/SOLU                    !进入求解器
*do,i,2,30,1
ANTYPE,TRANS,rest          !指定分析类型p235
TIMINT,ON,THERM         !打开瞬态效应P245
TIME,i*86400-(24-9)*3600           !定义计算终止时间
deltim,3600,3600,3600      !一个小时计算一次
AUTOTS,ON              !打开自动时间步长
OUTRES,,ALL            !写入每个子步
KBC,1                 !指定载荷步阶跃方式
LSEL,S,,,9,56,1
cm,xzy,line
sf,xzy,hflux,0    !施加热流荷载p273
ALLSEL               
SOLVE
finish
/SOLU  
ANTYPE,TRANS,rest          !指定分析类型
TIMINT,ON,THERM       !打开瞬态效应
TIME,i*86400           !定义计算终止时间
deltim,3600,3600,3600      !定义初始荷载步
AUTOTS,ON             !打开自动时间步长
OUTRES,,ALL           !写入每个子步
KBC,1                 !指定载荷步阶跃方式
LSEL,S,,,9,56,1
cm,xzy,line
sf,xzy,hflux,0    !施加热流荷载
ALLSEL               
SOLVE
*enddo
finish

二、结果：得到如下错误信息

三、讨论
从产生的错误信息来看，是第一次进入循环后，i=2,根本没有赋给表达式：i*86400-(24-9)*360,导致等于-54000,请教各位指点下！