Newmark-beta法解微分方程的MATLAB法式。
function[q,v,a]=newmark(M,C,K,F,q0,v0,dt,nt)%newmark-betamethod%[q,v,a]=newmark(M,C,K,F,q0,v0,dt,nt)%obtaintheresponseofthedynamicsystem%M-massmatrix%K-stiffnessmatrix%C-dampingmatrix%F-loadsmatrix(ntcolumns)%q0-initialdisplacement%v0-initialvelocity%dt-interval(timestep)%nt-numberofsamplingpoints%[q,v,a]-disp,velocity,acceleration
2016/11/15 17:43:41
1KB
MATLAB
1