Adams-Bashforth-Moulton法计算结构的动力响应

Yang J h

已于 2023-12-29 15:49:24 修改

阅读量803

点赞数 6

分类专栏： matlab数值计算方法文章标签：机器学习算法 matlab

于 2023-12-29 15:42:36 首次发布

本文链接：https://blog.csdn.net/yang_junhui/article/details/135292552

版权

matlab数值计算方法专栏收录该内容

1 篇文章

订阅专栏

function [D,V,A] = ABM4(M, C, K, F1, D0, V0, t)

% 利用4阶2.1.1.8 Adams-Bashforth-Moulton法计算结构的动力响应

%--------------------------------------------------------------------------

n_time = length(t);

A = zeros(n,n_time);

D = [D0,zeros( n, n_time-1)] ;

V = [V0,zeros( n, n_time-1)] ;

for i=4 : N1-1

D_predict = D(:,i) + (55*V(:,i) - 59*V(:,i-1) + 37*V(:,i-2) - 9*V(:,i-3))*dt/24;

V_predict = V(:,i) + (55*A(:,i) - 59*A(:,i-1) + 37*A(:,i-2) - 9*A(:,i-3))*dt/24;

A(:,i+1) = M\(F1(:,i+1)-C*V_predict-K*D_predict);

V(:,i+1) = V(:,i) + (9*A(:,i+1) + 19*A(:,i) - 5*A(:,i-1) + A(:,i-2))*dt/24;

D(:,i+1) = D(:,i) + (9*V(:,i+1) + 19*V(:,i) - 5*V(:,i-1) + V(:,i-2))*dt/24;

A(:,i+1) = M\(F1(:,i+1)-C*V(:,i+1)-K*D(:,i+1));

end



end