声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1960|回复: 2

[动力学和稳定性] 用newmark法求解振动方程中参数设置有什么要求?

[复制链接]
发表于 2018-1-17 20:27 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
% 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.
这里面的参数设置有什么讲究?请大牛指点。

回复
分享到:

使用道具 举报

 楼主| 发表于 2018-1-19 16:58 | 显示全部楼层
我顶,有人有相同的疑惑吗?被困扰了好几天
发表于 2018-11-7 20:30 | 显示全部楼层
楼主你联系方式多少?我目前也在做这个,可否交流一下
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-3-29 10:06 , Processed in 0.072215 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表