理论力学单自由度小阻尼受迫振动阻尼系数近似公式的推到

对于单自由度小阻尼受迫振动阻尼系数近似公式的推导有许多方法:
1.可以从基本表达式出发,一点一点的推到,也就是从振幅的表达式出发,求解出来最大振幅频率,然后求出最大振幅的表达式,让振幅等于最大振幅的一半,求解两个驱动 ,频率由表达式一步步的变换就可以得到解析解,然后近似变换得到近似公式,但是过程有点繁琐,而且近似变换放在最后,没有发挥应有的功能。
2.也可以由振动方程的形式,直接从能量的角度出发,应用积分方法与公式变换加上近似变换,舍去高阶小项得到
3.从振幅表达式出发,直接推导,中间过程近似舍去高阶小项,就可以了,下面从第三个思路出发,给出示例,注释放在了图形的后面


注:B:振幅,k:固有角频率,p:振动角频率,n:阻尼系数,h:驱动力与物体质量的比值


总结:1.对于解析解的精确求解,大部分人都会得到同样的答案,就是所谓的条条道路通罗马,然而对于一些近似变换,就是仁者见仁,智者见智了,能够快速准确的找到主要矛 盾是艺术的,也就是两点论中的重点论是一个艺术的东西,找到这个艺术的当然也是艺术的。
            2.由于有些单位一有时被人们认为不是通常意义上的单位,其实这是不对的,因为不是所有的单位一都是对等的,不要被个人的习惯所困惑,比如弧度与转等




  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
自由度振动是指系统中有三个自由度(即有三个可以独立变化的参数),并且到外力的驱动而产生振动。这是一种常见的力学振动问题,可以用数学模型和matlab代码来求解。 假设有一个三自由度振动系统,其运动方程可以表示为: m1 * x1'' + c1 * x1' + k1 * x1 - k2 * (x2 - x1) + F * cos(ω * t) = 0 m2 * x2'' + c2 * x2' - k2 * (x2 - x1) + k3 * (x3 - x2) = 0 m3 * x3'' + c3 * x3' + k3 * (x3 - x2) = 0 其中,m1、m2、m3分别为质量,c1、c2、c3为阻尼系数,k1、k2、k3为刚度系数, F为外力振幅,ω为外力频率,t为时间。 下面给出一个简的matlab代码来求解上述三自由度振动问题: ```matlab % 定义系统参数 m1 = 1; m2 = 1; m3 = 1; c1 = 0.1; c2 = 0.1; c3 = 0.1; k1 = 1; k2 = 1; k3 = 1; F = 1; ω = 2; t = 0:0.01:10; % 时间范围 % 定义初始条件 x1_0 = 0; x2_0 = 0; x3_0 = 0; x1_dot_0 = 0; x2_dot_0 = 0; x3_dot_0 = 0; % 定义运动方程 F1 = @(x1, x2) -k1 * x1 + k2 * (x2 - x1) - F * cos(ω * t); F2 = @(x1, x2, x3) -k2 * (x2 - x1) + k3 * (x3 - x2); F3 = @(x2, x3) k3 * (x3 - x2); % 使用ode45函数求解运动方程 [t, solution] = ode45(@(t, x) [... x(2); ... F1(x(1), x(2)) / m1 - c1 * x(2) / m1; ... x(4); ... F2(x(1), x(2), x(3)) / m2 - c2 * x(4) / m2; ... x(6); ... F3(x(2), x(3)) / m3 - c3 * x(6) / m3], t, [x1_0, x1_dot_0, x2_0, x2_dot_0, x3_0, x3_dot_0]); % 绘制位移-时间曲线 figure; plot(t, solution(:, 1), 'r-', 'LineWidth', 1.5); hold on; plot(t, solution(:, 3), 'g-', 'LineWidth', 1.5); hold on; plot(t, solution(:, 5), 'b-', 'LineWidth', 1.5); xlabel('时间'); ylabel('位移'); title('三自由度振动'); legend('x1', 'x2', 'x3'); ``` 以上代码使用ode45函数求解了三自由度振动的运动方程,并绘制出了位移-时间曲线。可以根据具体需求修改系统参数和初始条件来求解不同的问题。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值