谱分析问题请教

mlovehan

我是ANSYS新手，想做PSD分析，如图的简易结构
振动谱施加的面是A20,同时A20是约束固定面，
想得到的结果是H型叉边界点的UZ（位移）谱响应
振动谱值如下
频率(Hz)        PSD值（g^2/Hz）
20        0.01
80        0.04
350        0.04
2000        0.0056
按参考书以及例子，知道PSD分析步骤该是
建模分网
模态求解
模态扩展
获得谱解
模态叠加

按部就班照例子用GUI方式输入命令，POST26结果读不出来任意节点UZ的响应（PSD响应曲线），不知道是谱响应值没有叠加上，还是模态扩展方式有问题？
我用的是SUBPLACE模态求解器。因为对Reduced，不知道定义MDOF，而选择所有节点来定义，后续求解会出错。
因为原来结构相对比较复杂，所以没有采用批命令方式
所以有老师，高手如果有时间，可否帮忙过一遍下面的批命令，帮忙看下相应问题出在哪里.
非常感谢，敬待佳音。

批命令如下：
/TITLE,psd analysis of simple H-folk
/PREP7
BLOCK,0,1,0,0.05,0,0.05,
BLOCK,0,1,0.25,0.3,0,0.05,  
BLOCK,0.45,0.55,0.05,0.25,0,0.05,  
BLOCK,0.45,0.55,0.1,0.2,0.05,0.1,  !简单建模
VGLUE,all
ET,1,SOLID98
MP,DENS,1,2649                     ! density
MP,PERX,1,4.51                     ! permittivity (X direction)
MP,PERY,1,4.51                     ! permittivity (Y direction)
MP,PERZ,1,4.62                     ! permittivity (Z direction)
TB,ANEL,1                          ! define elastic table   
TBDATA,1,8.674E10,0.699E10,1.191E10,0,-1.791E10     ! [C] matrix
TBDATA,7,8.674E10,1.191E10,0,1.791E10   
TBDATA,12,10.72E10  
TBDATA,16,3.9875E10,0,-1.791E10
TBDATA,19,5.794E10  
TBDATA,21,5.794E10                 !定义单元类型和材料属性
mshape,1,3d     
mshkey,0         
allsel           
vmesh,all
eplot                              !分网，SMART1
FINISH  
/SOL
ANTYPE,2                           !模态求解，要求扩展，利用subspace求解器
MODOPT,SUBSP,10
EQSLV,FRONT
MXPAND,10, , ,1
LUMPM,0
PSTRES,0
MODOPT,SUBSP,10,0,5000, ,OFF
RIGID,  
SUBOPT,8,4,14,0,0,ALL               
FLST,2,1,5,ORDE,1                   !固支面约束   
FITEM,2,20  
/GO
DA,P51X,ALL,0
/STATUS,SOLU                        !模态求解
SOLVE                           
FINISH                             

/SOL                                 !模态扩张
EXPASS,1
MXPAND,10,0,0,1,0.005,  
/STATUS,SOLU
SOLVE   
FINISH

/SOLUTION                            !获得频解
ANTYPE,8
SPOPT,PSD,10,1  
PSDUNIT,1,ACCG,386.4,                !施加振动谱
PSDFRQ,1, ,20,80,350,2000,           !定义频率——谱值表中的频率
PSDVAL,1,0.01,0.04,0.04,0.0056,      !定义对应的谱值

ASEL,S,,,20
APLOT   
NSLA,R,1
NPLOT   
FLST,2,97,1,ORDE,2  
FITEM,2,98  
FITEM,2,-194
D,P51X,UZ,1.0                                 !施加节点振动谱
ALLSEL,ALL

PFACT,1,BASE,                          !计算参与系数
PSDRES,DISP,ABS                        !定义输出控制  
PSDRES,VELO,ABS
PSDRES,ACEL,ABS
SOLVE  
FINISH

/SOLU
ANTYPE,SPECTR
PSDCOM,0.005,10                   !定义PSD模态叠加
SOLVE
FINISH

/POST26
FILE,'psd-simple','rst','.'
STORE,PSD,1,                          !存储矢量到变量1
NSOL,2,15,U,Z,                        !定义变量2（记录节点15的Z向位移）
RPSD,3,2,,1,1，                       
/AXLAB,X,freq
/AXLAB,Y,UZ
PLVAR,3                                 !绘制RPSD曲线