Ansys-载荷步练习

!问题描述
!一块长10m宽2m厚0.1m的长方形板，左侧一端完全固定，右侧一端的外侧角点上施加一个横向集中力载荷
!改载荷按正弦规律变化，幅值为3000N。计算的是板在横向集中力的作用下的横向位移响应


finish
/clear,all
/prep7
et,1,shell63
r,1,0.1 !定义壳单元厚度实常数
mp,ex,1,2e11
mp,prxy,1,.3
mp,dens,1,7800
rectng,0,10,0,2
aatt,1,1,1,0 !指定矩形面的网格属性，材料，单元，实常数，单元坐标系
aesize,all,0.8
mshape,0,2d
mshkey,0
amesh,all
finish
save



/solu
antype,4 !瞬态分析
trnopt,full !选择完全法



dl,4,,all,
fz_kp=kp(10,2,0) !提取载荷施加关键点号
*SET,dt,0.1 !设置参数dt=0.1
*do,j,1,40 !执行40次重复计算

time,dt*j !载荷步时间
nsubst,1,,,1 !设置子部数目，格式 NSUBST，NSBSTP，NSBMX，NSBMN，Carry（直接指定子步数）？？？？
kbc,0 !载荷增量方式为渐变方式--kbc,0为斜坡载荷 kbc,1为阶跃载荷 对于瞬态分析，用户必须指定载荷类型
outres,all,all !载荷步输出所有子步的所有结果项
fk,fz_kp,fz,3000*sin(2*3.14*dt*j) !载荷按三角函数变化
!格式：FK，KPOL，Lab，VALUE1，VALUE2（在关键点上施加集中载荷）

solve
*enddo


time,40*dt+2 !激励撤消后的持续震动求解
nsubst,20,,,1
kbc,0
solve
finish


/post26
n_fz=node(10,2,0) !提出节点号，并存储到变量n_fz中
nsol,2,n_fz,u,z,uz !定义节点n_fz上的uz位移记录变量2
!格式：nsol，nvar，node，item，comp，name
plvar,2 !绘制变量2位移随时间的响应曲线