finish/clear/batch/title,HarmonicResponseoftheStructure/prep7width=1!定义分析变量length=2high=-1mass_hig=0.1et,1,63!定义单元类型et,2,4r,1,0.02!定义单元实常数r,2,2e-4,2e-8,2e-8,0.01,0.02mp,ex,1,2e11!定义材料特性mp,nuxy,1,0.3mp,dens,1,7800rect,,length,,width!定义有限元模型k,5,,,highk,6,length,,highk,7,length,width,highk,8,,width,highl,1,5*repeat,4,1,1esiz,0.1amesh,alltype,2real,2lmesh,5,8fini/soluanty,trans!选取瞬态分析nsel,s,loc,z,high!添加边界条件d,all,alloutr,all,allallsalph,5!指定质量阻尼系数!time,6!delt,0.1,0.05,0.5!auto,on!kbc,0TM_START=1e-8TM_END=6TM_INCR=0.1F1=5000freq=302*DO,TM,TM_START,TM_END,TM_INCRTIME,TMsfa,1,,pres,10000-F1*sin(2*3.1415*freq*TM)SOLVE*ENDDO!i=0!*DO,TM,TM_START,TM_END,TM_INCR!TIME,TM!sfa,1,,pres,10000-F1*sin(2*3.1415*freq*TM)!i=i+1!lswr,i!*ENDDO!lssolve,1,ifini/post26nsol,2,146,u,z!定义变量2为节点146在X方向上的位移nsol,3,146,acc,z!定义变量3为节点146在Z方向上的加速度esol,4,48,,s,eqv!定义变量4为单元48的EQV应力/grid,1plvar,2!绘制变量2的曲线plvar,3plvar,4