newmark法源程序

function [accveldsp]=TransResp5（kkccmmfdbcdofntdtq0dq0）
%--------------------------------------------------------------------------
%  Purpose:
%     The function subroutine TransResp5.m calculates transient response of
%     a structural system using Newmark integration scheme.
%  Synopsis:
%     [accveldsp]=TransResp5（kkccmmfdbcdofntdtq0dq0）
%  Variable Description:
%     Input parameters
%       kk cc mm - stiffness damping and mass matrices
%       fd - Input or forcing influence matrix
%       bcdof - Boundary condition dofs vector
%       nt - Number of time steps
%       dt - Time step size
%       q0 dq0 - Initial condition vectors
%     Output parameters
%       acc - Acceleration response
%       vel - Velocity response
%       dsp - Displacement response
%--------------------------------------------------------------------------
%  （1） initial condition
%--------------------------------------------------------------------------
[sdofn2]=size（kk）;
 
dsp=zeros（sdofnt）;                                         % displacement matrix
vel=zeros（sdofnt）;                                             % velocity matrix
acc=zeros（sdofnt）;                                          % acceleration matrix
 
dsp（:1）=q0;                                               % initial displacement
vel（:1）=dq0;                                                  % initial velocity

alpha=0.5; beta=0.5;                                        % select the parameters
%--------------------------------------------------------------------------
%  （2） Newmark integration scheme for time integration
%--------------------------------------------------------------------------
acc（:1）=inv（mm）*（fd（:1）