纽马克方法单自由度体系Matlab编程.pdfVIP

%NewmarkMethod w=200*0%单位t(吨) pp=p1(j-1)+w*(6*y0(j-1)/t0+3*y00(j-1))+c*(3* k=12.2*09.8/2.54*100;%单位kN/m y0(j-1)+t0*y00(j-1)/2); o=0.05;%阻尼比，这里阻尼比可以任意调整 y(j)=y(j-1)+pp/k1; c=2*w*o*(k/w)^0.5;%阻尼