用newmark法求解振动方程中参数设置有什么要求？

Roga

% 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);
    [E,E] = 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.
这里面的参数设置有什么讲究？请大牛指点。

Roga

我顶，有人有相同的疑惑吗？被困扰了好几天

chaostic

楼主你联系方式多少？我目前也在做这个，可否交流一下