请楼主说明和介绍资料啊，规范发帖啊~详细请参考： http://bbs.co188.com/list/viewthread.php?tid=1609707&extra=page%3D1%26amp%3Bfilter%3D0%26amp%3Borderby%3Ddateline%26amp%3Bascdesc%3DDESC&ofid=0

function [dis_new,vel_new,acc_new] = Newmark(M,C,K_inv,R_new,dis_old,vel_old,acc_old,a0,a1,a2,a3,a4,a5,a6,a7)
%--(1)--OUTPUT VARIBALES
%-------dis_new---new displcement vector at time (t + delt t);
%-------vel_new---new velocity vector at time (t + delt t);
%-------acc_new---new acceloration vector at time (t + delt t);
%--(2)--INPUT VARIBALES
%-------M---------mass matrix;
%-------C---------damping matrix;
%-------K---------stiffness matrix;
%-------dis_old---old displcement vector at time t;
%-------vel_old---old velocity vector at time t;
%-------acc_old---old acceloration vector at time t;
%-------K_inv-----inverse of effective stiffness matrix;
%-------R_new-----new load vector at time (t + delt t);

%---(3)---Calculate effective loads at time (t + delt t)
R_eff = R_new + M*(a0*dis_old + a2*vel_old + a3*acc_old)+ C*(a1*dis_old + a2*vel_old + a5*acc_old );
%---(4)----Calculate new displacement at time (t + delt t)
dis_new = K_inv*R_eff;
%---(5)----Calculate new velocity at time (t + delt t)
vel_new = vel_old + a6*acc_old + a7*acc_new;
%---(6)----Calculate new acceleration at time (t + delt t)
acc_new = a0*(dis_new - dis_old )- a2*vel_old - a3 * acc_old ;
就这些东西

[ 本帖最后由 lc850407 于 2010-4-10 17:48 编辑 ]

下载下来研究研究，正在研究matlab，研究其在土木工程中的应用

谢谢，楼主，好人一生平安！！！！！！