用newmark法求解振动方程中参数设置有什么要求？
% NEWMARK-β法% dt为时间变化步长此处假设为0.01/120s,tend总时长假设为0.03s,n为时间片段数目。
dt = 0.01/120;
tend = 0.03;
T = 0.01/120:dt:tend;
n = tend/dt;
gama = 0.5;
beta = 0.25;
M = zeros(19,19);
= size( M );
Nalpha0 = 1./beta./dt^2 ;
Nalpha1 = gama./beta./dt ;
Nalpha2 = 1./beta./dt ;
Nalpha3 = 1./2./beta - 1 ;
Nalpha4 = gama./beta - 1 ;
Nalpha5 = dt/2*(gama./beta-2) ;
Nalpha6 = dt*(1-gama) ;
Nalpha7 = gama*dt ;
Nd = zeros( E, floor(tend/dt)) ;
Nv = zeros( E, floor(tend/dt)) ;
Na = zeros( E, floor(tend/dt)) ;
t = zeros(n,1);
for ii = 1:n
t(ii,1) = dt*ii;
end
根据查询的文章，目前已知的是gama>=0.5,beta>=0.25(0.5+gama)^2,Newmark方法是无条件稳定，但是我也可以设置gama>1e5,beta>1e10.
这里面的参数设置有什么讲究？请大牛指点。
我顶，有人有相同的疑惑吗？被困扰了好几天 楼主你联系方式多少？我目前也在做这个，可否交流一下
页:
[1]