简支梁的随机振动分析

简支梁的随机振动分析

/PREP7 !进入前处理模块

/TITLE, EX 8.4(1) by Zeng P, Lei L P, Fang G

ET,1,BEAM4 !设置1号单元

MP,EX,1,2E11 $MP,NUXY,1,0.3 $MP,DENS,1,8000 !设定弹性模量泊松比密度L=10 !设置几何参数

R,1,4,1.333,1.333,2,2,0 !定义1号实常数AREA, IZZ, IYY , TKZ, TKY,

RMORE,0,2.2496,1.177,1.177 !1号实常数的7~10项为ISTRN, IXX, SHEARZ, SHEARY

K,1,0,0 $K,2,L,0 !生成两个关键点

L,1,2 !由1号点和2号点生成1号线

LESIZE,1, , ,10 !对1号线设置10段的网格划分

LMESH, 1 !对1号线划分网格

D,1,UX,,,,,UY,UZ,ROTX ! 在节点1处施加约束

D,2,UY,,,,,UZ !在节点2处施加约束

M,3,UY,11,1 !对3~11号节点的UY，定义主自由度UY

FINISH !结束前处理模块

/SOLU !进入求解模块

ANTYPE,MODAL ! 设定模态分析方式

MODOPT,REDUC !设置缩减算法

MXPAND,9,,,YES ! 设定模态扩展的阶数为9，并计算单元及支反力结果SOLVE !进行求解

*GET,FREQ,MODE,1,FREQ !提取第一阶模态的频率值，赋予参数FREQ

FINISH !结束

/SOLU !进入求解模块

ANTYPE,SPECTR !设置谱分析

SPOPT,PSD,9,ON !设定功率谱密度方式(针对9阶模态)

PSDUNIT,1,FORCE !设定功率谱密度分析的类型为作用力谱,为1号表DMPRAT,0.02 !定义阻尼比为0.02

F,1,FY,-0.5E6 !在端节点1处施加激励力1/2的FY力

F,2,FY,-0.5E6 !在端节点2处施加激励力1/2的FY力

F,3,FY,-1E6,,11,1 !在节点3~11处施加激励力FY

PSDFRQ,1,1,0.1,70. !设定频率范围0.1~70，前两位数值为数据表的编号PSDV AL,1,1,1 !设定功率谱值，前一位数为数据表的编号

PFACT,1,NODE !计算针对节点激振的参与系数

PSDRES,DISP,REL !设定针对位移的输出

PSDCOM !设定功率谱密度的模态合并方法

SOLVE $FINISH !进行求解，结束求解模块

/POST26 !进入时间历程后处理模块

STORE,PSD,10 !存储功率谱密度分析的数据

P=NODE(L/2,0,0) !获取位置(L/2,0,0)处的节点编号，赋给P

ELE=ENEARN(P) !获取节点P的单元编号

NSOL,2,P,U,Y !将P节点处的位移UY设为2号变量

RPSD,3,2 !计算针对2号变量的响应功率谱密度，存为3号变量PRTIME,42.640,42.641 !定义所列数据的时间间隔

PLV AR,3 !图形给出3号变量的曲线

*GET,P_UY,V ARI,3,RTIME,42.64 !提取3号变量在时间42.64处的值，赋予P_UY PM= P_UY*1000000 !进行单位换算

FINISH !结束

/POST26 !进入时间历程后处理模块

STORE,PSD,10 !存储功率谱密度的数据

ESOL,4,ELE,P,LS,7 !将单元ELM处节点P的应力设定为4号变量RPSD,5,4 !计算针对4号变量的响应功率谱密度，存为5号变量PRTIME,42.640,42.641 !设定一个时间(频率)范围

PLV AR,5 !图形显示5号变量的曲线

*GET,P2,V ARI,5,RTIME,42.64

PN=P2/(1E12) !进行单位换算

*STATUS !列显所有参数的内容

FINISH