最近做仿真的时候需要施加一个简谐力，就这么一个简单的需求，网上找了半天居然没找到合适的。。。

瞬态分析

! 瞬态分析
/solution
antype, trans
ff=16.350
T=1/ff
n=15                        ! 每周期取点
nn=40*n                        ! 总采样点
*do, i, 1, nn
    time, i*T/n
    allsel, all
    fk, 3, fy, 10*cos(2*pi*ff*i*T/n)
    outres, all, all
    solve
*enddo
finish

/post26
nsol, 2, 3, u, y, uy1
nsol, 3, 1, u, y, uy2

prcplx, 1                   ! 相位幅值输出
plcplx, 0
plvar, 2, 3
finish

很暴力的加一个循环。。。

谐分析

! 谐分析
/solution

allsel, all
fk, 3, fy, 10*cos(0), 10*sin(0)
allsel, all

antype, harmic
dmprat, 0
harfrq, 10, 50
nsubst, 200
solve
finish

/post26
nsol, 2, 3, u, y, uy1
nsol, 3, 1, u, y, uy2

prcplx, 1                   ! 相位幅值输出
plcplx, 0
plvar, 2, 3
finish

这个只能算稳态响应，但是优势在于计算快。