声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2267|回复: 11

[结构振动] matlab编程多自由度系统对激励的响应

[复制链接]
发表于 2012-12-8 12:05 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 ME! 于 2012-12-8 15:48 编辑

4.jpg
求助:这个多自由度系统怎么求响应啊求解

m1=[m11,0,0,m14,0,0,0,0,0,0],
m2=[0,m22,m23,0,m25,m26,0,0,0,0],
m3=[0,m32,m33,0,m35,m36,0,0,0,0],
M=[m1;m2;m3],   
k1=[k11,0,0,k14,0,0,0,0,0,0],
k2=[0,k22,k23,0,k25,k26,0,0,0,0],
k3=[0,k32,k33,0,k35,k36,0,0,0,0],
K=[k1;k2;k3],   
F(t)=25000 * (stepfun(t, 0) - stepfun(t, 0.1));
就是转化为一阶微分方程时,有两个二阶导数未知量,我不知道怎么求解,望高手帮忙?


回复
分享到:

使用道具 举报

发表于 2012-12-11 15:26 | 显示全部楼层
方程数量不足
 楼主| 发表于 2012-12-11 16:19 | 显示全部楼层
我是举个例子,实际上我的M K都为51阶的方阵,且为非对角阵,我不知道怎么用微分方程求解
例如第一个方程m(1,1)x(1)''+m(1,4)x(4)''+k(1,1)x(1)+k(1,4)x(4)=f(t),像这样的方程有很51个非线性方程组
发表于 2013-1-22 15:16 来自手机 | 显示全部楼层
楼主,我觉得你的系统方程可能错了,因为系统的刚度矩阵一般都是对称的!建议你检查一下
发表于 2013-1-23 09:20 | 显示全部楼层
龙格库塔方法
 楼主| 发表于 2013-1-23 14:13 | 显示全部楼层
我已经解出来了!
发表于 2013-1-24 11:45 来自手机 | 显示全部楼层
ME! 发表于 2013-1-23 14:13
我已经解出来了!

解出来了,不代表你就是对的!关键在于结果是否合理!
 楼主| 发表于 2013-1-28 16:37 | 显示全部楼层
和例题结果一样啊
发表于 2013-2-21 22:11 | 显示全部楼层
matlab编程多自由度系统对激励的响应
这个例题你有么
给我传一个啊
我的qq416015617
发表于 2013-2-22 14:52 | 显示全部楼层
先辈们已经把振动问题分成频域和时域两种解答方法。。。为什么到现在还要去纠结于解方程……
 楼主| 发表于 2013-2-27 22:19 | 显示全部楼层

去看看徐荣桥和徐斌的书,上面有例子的,资料很全,只是有些小错误
发表于 2013-2-28 08:58 | 显示全部楼层
书上还有小错误啊 呵呵
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-5-17 19:17 , Processed in 0.073664 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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