在线求高手指点

liwei4201

做荷载工况的时候出现问题如下，求解。
Load step 5 Substep 0 is not on file.

局部命令流如下：
/solu
eqslv,SPARSE                          !!!!!选择 稀疏矩阵直接法求解器
nropt,full                            !!!!!     完全牛顿拉夫逊法
NLGEOM,ON                               !大变形打开

acel,,,10                       !!!!!     定义结构的线性加速度
esel,all
ekill,all                        !!!!!     生死单元
。。。。。。。
。。。。。。。
time,5

lsel,none
esel,none
cmsel,s,maincable
esll
ealive,all

lsel,none
esel,none
cmsel,s,subcable
esll
ealive,all

*dim,mcf_area,array,14
*dim,mcf_target,array,14
*dim,mcf_input,array,14
*dim,mcf_result,array,14
*dim,mcf_error,array,14
*dim,mcf_perc,array,14

*dim,scf_area,array,11
*dim,scf_target,array,11
*dim,scf_input,array,11
*dim,scf_result,array,11
*dim,scf_error,array,11
*dim,scf_perc,array,11


cmsel,s,maincable
lsel,r,loc,x,0
esll
cm,mcable_1,elem

*do,i,2,13
        cmsel,s,maincable
        lsel,r,loc,x,-4*(i-1),4*(i-1)
        lsel,u,loc,x,-4*(i-2),4*(i-2)
        esll
        cm,mcable_%i%,elem
*enddo       

cmsel,s,subcable
lsel,r,loc,x,0
esll
cm,scable_1,elem

*do,i,2,11
        cmsel,s,subcable
        lsel,r,loc,x,-4*(i-1),4*(i-1)
        lsel,u,loc,x,-4*(i-2),4*(i-2)       
        esll
        cm,scable_%i%,elem
*enddo

cmsel,s,xigan
esll
cm,mcable_14,elem

/input,cf_array2,txt


N_max=20
Toler=5

/input,smcf_input2,txt
!!!!!!!!!!!!!!!!!输入初始迭代值!!!!!!!!!!!!!!!!!!!!!
!*do,i,1,14
        !mcf_input(i)=5*mcf_target(i)
!*enddo
       
!*do,i,1,11
        !scf_input(i)=10*scf_target(i)
!*enddo

!>>>>>>>>>>>>>>>>>>>>>>开始调索<<<<<<<<<<<<<<<<<<<<<<<<!
/solu
       
        *do,isolve,1,N_max
                *do,j,1,14
                        cmsel,s,mcable_%j%
                        /solu
                        bfe,all,temp,1,-mcf_input(j)/1.2e-5/mcf_area(j)/1.95e11
                *enddo
               
        *do,j,1,11
                cmsel,s,scable_%j%
                /solu
                bfe,all,temp,1,-scf_input(j)/1.2e-5/scf_area(j)/1.95e11
        *enddo
       
allsel
solve

        *do,j,1,14
                cmsel,s,mcable_%j%
                *get,aaa,elem,,num,max
                *get,mcf_result(j),elem,aaa,smisc,1
                mcf_error(j)=mcf_target(j)-mcf_result(j)
                mcf_perc(j)=abs(mcf_error(j)/mcf_target(j))*100
        *enddo
               
        *do,j,1,11
                cmsel,s,scable_%j%
                *get,aaa,elem,,num,max
                *get,scf_result(j),elem,aaa,smisc,1
                scf_error(j)=scf_target(j)-scf_result(j)
                scf_perc(j)=abs(scf_error(j)/scf_target(j))*100
        *enddo

        /input,cf_out2,txt
       
        *VSCFUN,mcable,max,mcf_perc
        *VSCFUN,scable,max,scf_perc
        *if,mcable,lt,Toler,and,scable,lt,Toler,exit
       
       
       
        !输入下一次迭代索力
        *do,i,1,14
                mcf_input(i)=mcf_input(i)+mcf_error(i)
        *enddo
       
        *do,i,1,11
                scf_input(i)=scf_input(i)+scf_error(i)
        *enddo
       

*enddo

allsel
solve

save

/post1

allsel

lcdef,1,1
lcdef,2,2
lcdef,3,3
lcdef,4,4
lcdef,5,5