finish /clear /batch

/title, Harmonic Response of the Structure /prep7

width=1!定义分析变量 length=2 high=-1

mass_hig=0.1

et,1,63!定义单元类型 et,2,4

r,1,0.02!定义单元实常数 r,2,2e-4,2e-8,2e-8,0.01,0.02 mp,ex,1,2e11!定义材料特性 mp,nuxy,1,0.3 mp,dens,1,7800

rect,,length,,width!定义有限元模型 k,5,,,high

k,6,length,,high

k,7,length,width,high k,8,,width,high l,1,5

*repeat,4,1,1 esiz,0.1 amesh,all type,2 real,2 lmesh,5,8 fini /solu

anty,trans!选取瞬态分析

nsel,s,loc,z,high!添加边界条件 d,all,all outr,all,all alls

alph,5!指定质量阻尼系数 !time,6

!delt,0.1,0.05,0.5 !auto,on !kbc,0

TM_START=1e-8 TM_END=6 TM_INCR=0.1 F1=5000

freq=302

*DO,TM,TM_START,TM_END,TM_INCR TIME,TM

sfa,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,i fini /post26

nsol,2,146,u,z!定义变量2为节点146在X方向上的位移 nsol,3,146,acc,z!定义变量3为节点146在Z方向上的加速度 esol,4,48,,s,eqv!定义变量4为单元48的EQV应力 /grid,1

plvar,2!绘制变量2的曲线 plvar,3 plvar,4