資源簡介
Newmark-beta法解微分方程的MATLAB程序。
function [q,v,a]=newmark(M,C,K,F,q0,v0,dt,nt)
% newmark-beta method
% [q,v,a]=newmark(M,C,K,F,q0,v0,dt,nt)
% obtain the response of the dynamic system
% M - mass matrix
% K - stiffness matrix
% C - damping matrix
%
代碼片段和文件信息
function?[qva]=newmark(MCKFq0v0dtnt)
%?newmark-beta?method
%???????????[qva]=newmark(MCKFq0v0dtnt)
%???????????obtain?the?response?of?the?dynamic?system
%???????????M?-?mass?matrix
%???????????K?-?stiffness?matrix
%???????????C?-?damping?matrix
%???????????F?-?loads?matrix(nt?columns)
%???????????q0?-?initial?displacement
%???????????v0?-?initial?velocity
%???????????dt?-?interval(time?step)
%???????????nt?-?number?of?sampling?points
%???????????[qva]?-?dispvelocityacceleration
????if?nt~=?length(F)
????????disp(‘The?size?of?force?matrix?not?match?the?number?of?sampling?points?‘)
????????return
????end
????dof=length(M);
????alpha=0
評論
共有 條評論