简支梁的地震响应分析

简支梁的地震响应分析

/PREP7 !进入前处理模块 /TITLE, EX 8.4(3) by Zeng P, Lei L P, Fang G ET,1,BEAM3 !设定1号单元

L=240 $A=273.9726 $H=14 $I=1000/3 !设定几何参数 R,1,273.9726,(1000/3),14 !设定1号实常数(梁单元)

MP,EX,1,3E7 $MP,PRXY,1,0.3 $MP,DENS,1,73E-5 !设定弹性模量, 泊松比, 密度

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

ESIZE,,8 !设定单元网格划分的分段数 LMESH,1 !对1号线划分单元网格 NSEL,S,LOC,X,0 !选择位置x=0的节点

D,ALL,UY !对所选择的节点施加位移约束UY=0 NSEL,S,LOC,X,L !选择位置x=L的节点

D,ALL,UX,,,,,UY !对所选择的节点施加位移约束UX=UY=0 NSEL,ALL !选择所有节点 FINISH !结束前处理模块 /SOLU !进入求解模块

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

MODOPT,REDUC,,,,3 !设置缩减算法，提取3阶模态

MXPAND,1,,,YES ! 设定模态扩展的阶数为1，并计算单元及支反力结果

M,ALL,UY !对所有节点定义主自由度UY OUTPR,BASIC,1 !设置输出结果的方式 SOLVE !进行求解

*GET,F1,MODE,1,FREQ !提取第一阶模态频率，赋给F1 FINISH !结束

/SOLU !进入求解模块 ANTYPE,SPECTR !设定谱分析方式 SPOPT,SPRS !设定单点激励谱分析

SED,,1, !设定单点激励的方向为Y轴 SVTYP,3 !指定单点响应谱类型为地震位移谱 FREQ,.1,10 !设定频率数据表格的频率点

SV,,.44,.44 !设定频率数据表格的对应于频率点的激励值 SOLVE !进行求解

*GET,F1_COEF,MODE,1,MCOEF !提取模态1的谱分析结果的模态系数 FINISH !结束求解

/POST1 !进入一般性后处理模块

SET,1,1,F1_COEF !调出第1阶模态的结果，并乘以模态系数 PRNSOL,DOF !打印节点结果 PRESOL,ELEM !打印单元结果 PRRSOL,F !打印支反力结果

ETABLE,STS,LS,3 !建立单元数据表格(梁单元的轴向应力)，名称为STS

*GET,P_SX,ELEM,5,ETAB,STS !在单元表STS中提取5号单元(梁的中间位置)的应力值，赋给P_SX

*GET,P_UY,NODE,6,U,Y !提取6号节点(梁的中间位置)的位移UY，赋给P_UY

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