问题遇到的现象和发生背景
想用一个单子度做IDA即动力分析，但是在求频率时总是报错
问题相关代码，请勿粘贴截图
模型文件：
model
``
```FINISH !从当前处理器安全的退出
/clear !清除之前的操作
/PREP7 !进入前处理器
!TYPE OF ELEMENT*
ET,1,combin40 !弹簧单元
keyopt,1,3,3 !单元Uz方向自由度,只有z方向有自由度
MPTEMP,,,,,,,,
MPTEMP , 1, 0
MPDATA , EX , 1 , , 2.1E11 !弹性模量
MPDATA , PRXY , 1 , , 0.3 !泊松比
MPDATA , DENS , 1 , , 7800 !密度
!非线性材料属性表并定义
tb,bkin,1 !双线型随动强化材料选项,默认属性为1
tbdata,1,235e6,0.0 !给当前数据表定义数据,屈服强度235Mpa
N, 1,0,0,0 !创建节点
N, 2,0,0,1
TYPE,1 !创建弹簧单元
r,1
REAL,1 !定义实常数
mat,1
E,1,2 !弹簧单元是用节点1，2构成的
D,1,ALL !设置边界条件,约束
ET,40,mass21 !加质量源
R,40,,,1.62E7 !质量
TYPE,40
REAL,40
E,2 !加在2号节点上
allsel,all
主程序：
finish
/clear
*ulib,modelSY,txt
*use,model
/SOLU !相当于打开了Solution
ANTYPE,0 !定义，0为静力分析
NLGEOM,off !大变形开关
SOLVE !求解运算
FINISH
/solu
ANTYPE,2 !算模态
MODOPT,LANB,100 !模式选择，lanb模式，100阶
EQSLV,sparse !求解器（稀松矩阵）
MXPAND,0, , ,0 ! 模态扩展
LUMPM,0 ! 一致质量矩阵
/STATUS,SOLU
SOLVE
finish !求解运算
/post1
nfreq=100
*DIM,ffreq,,nfreq,2
*DO,I,1,nfreq
set,1,i
ffreq(I,1)=I
*get,ffreq(I,2),active,0,set,freq
*ENDDO
*mwrite,ffreq(1,1),freq,txt,,jik,2,nfreq !输出频率
(f8.0,f12.3)
finish
fre1=ffreq(1,2)
fre2=ffreq(7,2)
!/OUTPUT,TEMP,TXT
PI=3.141592
C=0.05
A1=1./(2.*2.PIFRE1)
B1=2.PIFRE1/2.
A2=1./(2.2.PIFRE2)
B2=2.PIFRE2/2.
BET=C/(FRE1+FRE2)/PI !BET阻尼系数
ALF=(C-B1BET)/A1 !ALF阻尼系数
finish
save
npoin=1 !加载点数
ntiqu=2 !一共有多少个节点
zyd=npoin*3 !自由度
ndanyuan=1 !一共有多少个单元
nstep=100 !可以改
nstep2=100
mchose=nstep
step=0.05 !对应的实际采样频率为50
*DIM,idforce,,npoin,1
*vREAD,idforce(1,1),idforce,txt,,jiK,1,npoin !(1,1)是第一行第一列,在2加载
(f8.0)
*DIM,DATA,,zyd,nstep
*vREAD,DATA(1,1),pcf,txt,,jiK,nstep,zyd !时程文件
(3943f12.2)
*DIM,idforcedanyuan,,ndanyuan,1
*vREAD,idforcedanyuan(1,1),idforcedanyuan,txt,,jiK,1,ndanyuan
(f8.0)
*DIM,idtiqu,,ntiqu,2
*vREAD,idtiqu(1,1),idtiqu,txt,,jiK,1,ntiqu
(f8.0)
finish
/CONFIG,NRES,90000 !允许存在结果文件中的最大号码90000步
/solu
anty,trans !选取瞬态分析
trnopt,full !完全瞬态分析法
NLGEOM,on !大变形
NROPT,MODI, ,OFF !Newton-Raphson，使用改进的Newton-Raphson，不要使用自适应下降（在所有其他情况下默认）
EQSLV,sparse !求解器（稀松矩阵）
ALPH,ALF !使用瑞丽阻尼
BETA,BET
autots,on !使用自动时间步
NCNV,0,0,0,0,0, !设置终止分析的键，第一个0：如果解决方案无法收敛，请勿终止分析，最大节点自由度解值（位移、温度等），迭代次数，运行时间，CPU时间
autots,on !使用自动时间步
timint,on !打开瞬态效果，包括瞬态（质量或惯性）效应。
kbc,0 !从上一加载步的值到该加载步的值，为每个子步线性插值（倾斜）载荷。
OUTPR,basic,last, !输出基本量，每个加载步骤的最后一个子步骤打印此项的解决方案。
outres,basic,last, !写入数据库的数据，每个加载步骤的最后一个子步骤写入指定的解决方案结果项。
time,0.005 !设置每一步的时间
acel,0.0,0.0,9.8 !设置加速度
SOLVE !求解
timint,on
nsubst,2,100,1
autots,on
kbc,0
fangdabeishu=0.8
DO,Istep,1,nstep
TIME,istepstep
acel,0,0.0,9.8
allsel
*DO,I,1,npoin
f,idforce(i),fx,DATA((i-1)*3+1,Istep)*fangdabeishu
f,idforce(i),fy,DATA((i-1)*3+2,Istep)*fangdabeishu
f,idforce(i),fz,DATA((i-1)*3+3,Istep)*fangdabeishu
*ENDDO
allsel
SOLVE
*ENDDO
FINISH
finish
/clear
*ulib,modelSY,txt
*use,model
/SOLU !相当于打开了Solution
ANTYPE,0 !定义，0为静力分析
NLGEOM,off !大变形开关
SOLVE !求解运算
FINISH
/solu
ANTYPE,2 !算模态
MODOPT,LANB,100 !模式选择，lanb模式，100阶
EQSLV,sparse !求解器（稀松矩阵）
MXPAND,0, , ,0 ! 模态扩展
LUMPM,0 ! 一致质量矩阵
/STATUS,SOLU
SOLVE
finish !求解运算
/post1
nfreq=100
*DIM,ffreq,,nfreq,2
*DO,I,1,nfreq
set,1,i
ffreq(I,1)=I
*get,ffreq(I,2),active,0,set,freq
*ENDDO
*mwrite,ffreq(1,1),freq,txt,,jik,2,nfreq !输出频率
(f8.0,f12.3)
finish
fre1=ffreq(1,2)
fre2=ffreq(7,2)
!/OUTPUT,TEMP,TXT
PI=3.141592
C=0.05
A1=1./(2.*2.PIFRE1)
B1=2.PIFRE1/2.
A2=1./(2.2.PIFRE2)
B2=2.PIFRE2/2.
BET=C/(FRE1+FRE2)/PI !BET阻尼系数
ALF=(C-B1BET)/A1 !ALF阻尼系数
finish
save
npoin=1 !加载点数
ntiqu=2 !一共有多少个节点
zyd=npoin*3 !自由度
ndanyuan=1 !一共有多少个单元
nstep=100 !可以改
nstep2=100
mchose=nstep
step=0.05 !对应的实际采样频率为50
*DIM,idforce,,npoin,1
*vREAD,idforce(1,1),idforce,txt,,jiK,1,npoin !(1,1)是第一行第一列,在2加载
(f8.0)
*DIM,DATA,,zyd,nstep
*vREAD,DATA(1,1),pcf,txt,,jiK,nstep,zyd !时程文件
(3943f12.2)
*DIM,idforcedanyuan,,ndanyuan,1
*vREAD,idforcedanyuan(1,1),idforcedanyuan,txt,,jiK,1,ndanyuan
(f8.0)
*DIM,idtiqu,,ntiqu,2
*vREAD,idtiqu(1,1),idtiqu,txt,,jiK,1,ntiqu
(f8.0)
finish
/CONFIG,NRES,90000 !允许存在结果文件中的最大号码90000步
/solu
anty,trans !选取瞬态分析
trnopt,full !完全瞬态分析法
NLGEOM,on !大变形
NROPT,MODI, ,OFF !Newton-Raphson，使用改进的Newton-Raphson，不要使用自适应下降（在所有其他情况下默认）
EQSLV,sparse !求解器（稀松矩阵）
ALPH,ALF !使用瑞丽阻尼
BETA,BET
autots,on !使用自动时间步
NCNV,0,0,0,0,0, !设置终止分析的键，第一个0：如果解决方案无法收敛，请勿终止分析，最大节点自由度解值（位移、温度等），迭代次数，运行时间，CPU时间
!nsubst,100,50000,100 !100为指定子步数,50000为使用自动时间跟踪时,最大的子步数(PS因为打开自动步长时不可能完全是设定的50步),100为使用自动时间跟踪时,最小的子步数
autots,on !使用自动时间步
timint,on !打开瞬态效果，包括瞬态（质量或惯性）效应。
kbc,0 !从上一加载步的值到该加载步的值，为每个子步线性插值（倾斜）载荷。
OUTPR,basic,last, !输出基本量，每个加载步骤的最后一个子步骤打印此项的解决方案。
outres,basic,last, !写入数据库的数据，每个加载步骤的最后一个子步骤写入指定的解决方案结果项。
time,0.005 !设置每一步的时间
acel,0.0,0.0,9.8 !设置加速度
SOLVE !求解
timint,on
nsubst,2,100,1
autots,on
kbc,0
fangdabeishu=0.8
DO,Istep,1,nstep
TIME,istepstep
acel,0,0.0,9.8
allsel
*DO,I,1,npoin
f,idforce(i),fx,DATA((i-1)*3+1,Istep)*fangdabeishu
f,idforce(i),fy,DATA((i-1)*3+2,Istep)*fangdabeishu
f,idforce(i),fz,DATA((i-1)*3+3,Istep)*fangdabeishu
*ENDDO
allsel
SOLVE
*ENDDO
FINISH
/post1
nhang=3*ntiqu
!nhang=npoinwm+nbeam
*DIM,disp,,nhang,mchose
!*DIM,stress,,ndanyuan,mchose
*DO,Istep,1,mchose
set,Istep,last
*DO,I,1,ntiqu
GET,disp(3i-2,Istep),node,idtiqu(i),u,x
GET,disp(3i-1,Istep),node,idtiqu(i),u,y
GET,disp(3i,Istep),node,idtiqu(i),u,z
*END DO
*END DO
*CFOPEN,dispreasult,txt
*DIM,dddisp,,mchose
*DO,I,1,nhang
*DO,J,1,mchose
dddisp(J)=disp(I,J)
*ENDDO
*VWRITE,dddisp(1)
(F30.15,$)
*VWRITE
(' ')
*ENDDO
*CFCLOSE
finish
/post1
nhang=3*ntiqu
!nhang=npoinwm+nbeam
*DIM,disp,,nhang,mchose
!*DIM,stress,,ndanyuan,mchose
*DO,Istep,1,mchose
set,Istep,last
*DO,I,1,ntiqu
GET,disp(3i-2,Istep),node,idtiqu(i),u,x
GET,disp(3i-1,Istep),node,idtiqu(i),u,y
GET,disp(3i,Istep),node,idtiqu(i),u,z
*END DO
*END DO
*CFOPEN,dispreasult,txt
*DIM,dddisp,,mchose
*DO,I,1,nhang
*DO,J,1,mchose
dddisp(J)=disp(I,J)
*ENDDO
*VWRITE,dddisp(1)
(F30.15,$)
*VWRITE
(' ')
*ENDDO
*CFCLOSE
finish
c## 运行结果及报错内容
报错：load set not found on result file 1.rst
