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