有关newmark method的响应图

林曉維

如有XX烦请看图一

file:///C:/DOCUME~1/Jamin/LOCALS~1/Temp/msohtmlclip1/01/clip_image002.gif
请问我这样一个式子利用FEM分析,用newmark method去算时间域的响应,如果我初始值给的是位移给的是first mode,利用F=kq算出外力,再利用F=ma算出加速度,速度给零而我的file:///C:/DOCUME~1/Jamin/LOCALS~1/Temp/msohtmlclip1/01/clip_image004.gif
file:///C:/DOCUME~1/Jamin/LOCALS~1/Temp/msohtmlclip1/01/clip_image006.gif,这样跑出来的响应图会是homogeneous如图二,位移跟速度的响应会是一个漂亮的圆如图三

可是我现在要改成初始值给的是加速度,w方向的都给值,file:///C:/DOCUME~1/Jamin/LOCALS~1/Temp/msohtmlclip1/01/clip_image008.gif方向给零,位移和速度是零,file:///C:/DOCUME~1/Jamin/LOCALS~1/Temp/msohtmlclip1/01/clip_image010.gif
file:///C:/DOCUME~1/Jamin/LOCALS~1/Temp/msohtmlclip1/01/clip_image006.gif, ,这样跑出来的响应图不会是homogeneous, 位移跟速度的响应也不会是一个漂亮的圆如图四

请问我有哪个地方错了吗??请板友们帮帮忙,谢谢

[ 本帖最后由 林曉維 于 2008-12-25 16:41 编辑 ]