beam-shell结构谱分析程序 哪位大侠看看出了什么问题

yuanbobeyond

/TITLE,PSD ansysis of Beam-Shell struture
/PREP7
ET,1,BEAM4
ET,2,SHELL63                     !定义单元类型
R,1,0.01,8.3E-5,8.3E-5,0.1,0.1,,
R,2,0.02,
MP,EX,1,2.0E8
MP,PRXY,1,0.3
MP,DENS,1,7.8                    !定义材料特性
K,,,,,
K,,11,,,
K,,11,,11,
K,,0,,11,
KGEN,2,1,4,,,-10,,,0
KGEN,2,5,8,,,-10,,,0
KGEN,2,9,12,,,-10,,,0
KGEN,2,13,16,,,-10,,,0  
LSTR,1,5  
LSTR,5,9  
LSTR,9,13  
LSTR,2,6  
LSTR,6,10  
LSTR,10,14  
LSTR,3,7  
LSTR,7,11  
LSTR,11,15  
LSTR,4,8  
LSTR,8,12
LSTR,12, 16
A,1,2,3,4,   
A,5,6,7,8,
A,9,10,11,12,                      !以关键点生成面
k,17,6.5,0,6.5,

AATT,1,2,2,0,                      !定义面属性
LSEL,S,LINE,,13,24,                !选择面上的线         
LESIZE,ALL,,,11,,,,,1     
MSHAPE,0,2D
MSHKEY,0
AMESH,ALL                           !划分面

LATT,1,1,1,,17,,                   !定义线的属性
LSEL,S,LINE,,1,12,                 !选择梁
LESIZE,ALL,,,10,,,,,1
LMESH,ALL                          !划分、生成梁单元
ALLSEL

/VIEW,1,1,1,1
/REPLOT
FINISH
/SOLUTION
KSEL,S,KP,,13,16
DK,ALL,ALL                !在13-16点施加约束
ANTYPE,MODAL              !选择模态分析
MODOPT,LANB,10            !制定模态提取方法、阶数
EQSLV,SPAR                !制定求解器
MXPAND,10,,,YES,          !进行10阶模态扩展
LUMPM,0                   !使用一致质量矩
PSTRES,0                  !不计算预应力
SOLVE
FINISH

/SOLUTION
ANTYPE,SPECTR               !分析为谱分析
SPOPT,PSD,10,YES        
PSDUNIT,1,ACEL,386.4,        !功率谱类型为加速度
PSDFRQ,1, ,10,20,120,121,200,240,340,500,
PSDVAL,1,0.0065,0.0065,0.0002,0.003,0.003,0.0015,0.00003,0.00015,
                             !定义频率与功率谱密度
DMPRAT,0,                    !阻尼为0
KSEL,S,KP,,13,16
DK,ALL,UY,1.0, ,1            !施加基础激励
LVSCALE,0.0003               !施加压力功率谱激励
PFACT,1,BASE,                !计算PSD激励参与系数
PSDRES,DISP,REL
PSDRES,VELO,OFF
PSDRES,ACEL,OFF              !设置输出项
SOLVE
FINISH

/SOLUTION
ANTYPE,SPECTR   
SPOPT,PSD,10,1
PSDCOM,0.005,10,                !模态合并
SOLVE  
FINISH


/POST1
/DSCALE1,10                    !定义变形比例因子为10
PLDISP,1                       !图形显示结构变形
/POST26                        !计算PSD响应
STORE,PSD,1                    !存储频率矢量到变量 1
NSOL,2,535,U,Y,UY535           !定义变量 2 (记录节点535的Y方向位移)
RPSD,3,2,,1,1,                 !计算节点 535 在Y方向上的位移响应（最后                                一个1表示取的是相对数值）
RPSD,4,2,,2,1,                 !计算节点 535 在Y方向上的速度响应
RPSD,5,2,,3,1,                 !计算节点 535 在Y方向上的加速度响应
PLTIME,0,4
/AXLAB,X,Frequency  
/AXLAB,Y,Displacement-Y
PLVAR,3, , , , , , , , , ,      !绘制节点 535 在Y 方向上的位移响应曲线
/AXLAB,X,Frequency  
/AXLAB,Y,Velocity-Y
PLVAR,4, , , , , , , , , ,
/AXLAB,X,Frequency  
/AXLAB,Y,Aacceleration-Y
PLVAR,5, , , , , , , , , ,
RESET                           !设置POST26成缺省的初始状态
NSOL,2,535,U,Y,UY535            !定义变量 2 (记录节点535的Y方向位移)
NSOL,3,507,U,Y,UY507
CVAR,4,2,3,1,2,                 !计算两节点位移协方差
*GET,U535,VARI,4,EXTREM,CVAR    !提取协方差值

hongshuaijun

模态扩展最好是单独作为一个求解步骤，然然后读入.mcom文件