求一个计算悬臂梁简谐载荷作用下动力学响应的程序

bisifeng

一维连续系统的悬臂梁，L长为2m，其他参数E,I,截面积A,密度等取钢的默认值即可，或者直接用符号计算。在其自由端作用幅值为1N的正弦激励，频率在0-1000Hz变化。当然，可以先计算激励频率为定值（例如4.2Hz）时的响应。
计算得到结果应该是频域内的复幅值，与时间无关。最好可以算出梁的纵向位移y(x)（挠度），剪力Q(x)，弯矩M(x)。
最好给出源程序，并适当解释一下。

结构动力学书籍中关于悬臂梁受迫振动没怎么看懂。。请各位高手帮下忙啦

bisifeng

:Q