beam-shell结构谱分析程序 哪位大侠看看出了什么问题
/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 !提取协方差值 模态扩展最好是单独作为一个求解步骤，然然后读入.mcom文件
页:
[1]