ansys载荷施加命令流求教

火电机组

各位，我现在要用瞬态分析求解，载荷的施加情况如下

      其中r为单元的响应位移，a是定值。

我现在能够做出来的命令如下

pi=2*asin(1)
ff=0.32
t=1/ff
n=20
d,3,all
*do,i,1,5*n
time,t/n*i
f,1,fy,60*sin(2*pi*ff*t/n*i)
solve
*enddo

按照这个命令流我只是施加了简谐力的部分。

我想在DO循环中施加一下选择项，从而实现K(r-a)的受力部分，可否请指点一下，主要是选择结构的判断条件如何用命令流实现以及在循环语句中的放置，

谢谢

三民

建议用matlab来做！