24小时热门版块排行榜    

北京石油化工学院2025年硕士研究生招生接受调剂公告
查看: 6313  |  回复: 13
【悬赏金币】回答本帖问题,作者wzxchars将赠送您 10 个金币

wzxchars

铁杆木虫 (著名写手)

[求助] 如何用matlab求解矩阵系数的二阶微分方程已有1人参与

以一个简单的有限元动力方程为例,M*X''+K*X=F(t),其中M和K为已知质量和刚度矩阵(16*16),F(t)为已知外部载荷输入(16*1),如何求解X呢?用ode45还是用simulink?希望高手们给出详细步骤,或相关例子

[ 来自科研家族 机器人 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

静坐常思自己过,闲谈莫论他人非
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

jun_1129

铜虫 (初入文坛)

【答案】应助回帖

function z=zi(t,x)
global m
global k
z1=x(4);
z2=x(5);
z3=x(6);
z4=(k/m)*x(2)-3*(k/m)*x(1);
z5=(k/m)*x(1)-2*(k/m)*x(2)+(k/m)*x(3);
z6=(k/m)*x(2)-3*(k/m)*x(3);
z=[z1;z2;z3;z4;z5;z6];
end

……………………main.m………………………………
global m
global k
z1=1;
z2=1;
z3=-1;
z4=1;
z5=1;
z6=1;
m=2;
k=1;
t1=0;
t2=10;
[t,f]=ode45(@func,[t1 t2],[z1 z2 z3 z4 z5 z6])
plot(t,f)
8楼2013-10-30 17:59:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jun_1129

铜虫 (初入文坛)

【答案】应助回帖

上述程序为三个质量块四个弹簧的无阻尼震动波形图
前三个变量为位移后三个变量为速度
9楼2013-10-30 18:01:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wzxchars: 金币+5, 有帮助 2012-10-24 19:30:25
xiegangmai: 金币+1, 谢谢参与 2012-10-24 19:57:25
比较麻烦,你搜搜帮助里的Several Differential Equations有个用dsolve的方法,但是估计不好解出来。
showmethemoney
2楼2012-10-24 17:07:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzxchars

铁杆木虫 (著名写手)

引用回帖:
2楼: Originally posted by csgt0 at 2012-10-24 17:07:18
比较麻烦,你搜搜帮助里的Several Differential Equations有个用dsolve的方法,但是估计不好解出来。

ode45是能够解这种带有时间参数的常系数二阶微分方程的,但对于这种矩阵形式的解法我还没有用过,有一种最基本的方法是矩阵拆分后联立微分方程组求解,但是实际上我要解的是100*100的矩阵,光是函数表达就要有好几百个项,唉头大,希望能在这里遇到高手指点。谢谢版主的热心回复
静坐常思自己过,闲谈莫论他人非
3楼2012-10-24 19:30:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

eezyzyzzzee

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wzxchars: 金币+15, ★★★很有帮助, 谢谢了 2012-10-24 23:32:27
xiegangmai: 金币+1, 谢谢参与 2012-10-27 19:43:23
你搜搜——结构分析的有限元法与MATLAB程序设计_徐荣桥.pdf
看看最后一章
http://wenku.baidu.com/view/5cf31dc8050876323112124f.html

不懂你的领域,就是偶然看到过上面这个教程,有动力学响应这块的
4楼2012-10-24 21:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzxchars

铁杆木虫 (著名写手)

引用回帖:
4楼: Originally posted by eezyzyzzzee at 2012-10-24 21:25:09
你搜搜——结构分析的有限元法与MATLAB程序设计_徐荣桥.pdf
看看最后一章
http://wenku.baidu.com/view/5cf31dc8050876323112124f.html

不懂你的领域,就是偶然看到过上面这个教程,有动力学响应这块的...

看了一下比较复杂,方程组矩阵系数实在太大了,用Newmark逐步积分求解,得好好编程。谢谢应助
静坐常思自己过,闲谈莫论他人非
5楼2012-10-24 23:33:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

twoflm

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wzxchars: 金币+10, ★★★很有帮助, 谢谢啊 2012-10-26 08:22:14
xiegangmai: 金币+1, 谢谢参与 2012-10-27 19:43:36
看了上面的几个帖子,不知道是不是跟结构动力学相关的书比较相关?
《Matlab有限元结构动力学分析与工程应用》的书及程序的下载:
http://ishare.iask.sina.com.cn/f/21145871.html
http://ishare.iask.sina.com.cn/f/23700400.html?from=like
6楼2012-10-25 08:28:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

放飞心中梦想

新虫 (初入文坛)

如果质量矩阵M是奇异的,是不是就求解不出来呢?
7楼2012-12-28 10:44:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
10楼2014-10-06 16:17:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wzxchars 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 313求调剂 +20 2426399003 2025-04-13 21/1050 2025-04-15 09:53 by 2426399003
[考研] 305生物求调剂086000专硕六级 +19 路过的野人谷 2025-04-13 20/1000 2025-04-15 09:27 by taosheng130
[考研] 考研调剂 +11 Apailian 2025-04-12 11/550 2025-04-14 23:02 by 白白12138
[考研] 生物301能调去哪 +8 池恩陈_ 2025-04-10 12/600 2025-04-14 22:33 by vivegetang
[考研] 071生物学323分 +11 土豆补数学 2025-04-10 11/550 2025-04-14 10:01 by 07759hry
[考研] 大连工业大学接收调剂生 +3 造纸的神 2025-04-14 3/150 2025-04-14 08:55 by 轻工造纸
[考研] 262求调剂 +4 v1cent_tree 2025-04-12 7/350 2025-04-13 17:16 by Xiaomu121.
[考研] 求调剂(四六级已过) +7 观澜湖2002 2025-04-12 8/400 2025-04-13 14:28 by forestmumu
[考研] 材料调剂 +21 yo黎明 2025-04-10 24/1200 2025-04-13 14:16 by 齐工大轻工
[考研] 285求调剂 +5 美味轩 2025-04-12 5/250 2025-04-12 15:57 by DNDYR
[论文投稿] 本科生综述投稿 +4 kiligg 2025-04-10 5/250 2025-04-11 11:41 by nono2009
[考研] 已调剂 +12 高等数学甲 2025-04-10 13/650 2025-04-11 10:32 by Smorant
[考研] 265求调剂 +4 sqc0802 2025-04-08 4/200 2025-04-11 10:02 by shenrf
[考研] 生物与医药288求调剂 +14 yelongyu288 2025-04-10 15/750 2025-04-11 07:47 by calebby
[考研] 279求调剂 +16 Two.Dark 2025-04-08 17/850 2025-04-10 20:10 by lhy2902
[考研] 293求调剂(本科 工程管理 研究生报考交通运输,肯吃苦) +3 wcq719 2025-04-09 3/150 2025-04-10 15:42 by 褚迎松小小
[考研] 372求调剂 +3 yuiyiy 2025-04-09 3/150 2025-04-10 13:21 by orlando
[考研] 二次调剂 求稳 +14 小路 一一 2025-04-09 15/750 2025-04-09 22:43 by barlinike
[考研] 302求调剂 +3 毛毛毛, 2025-04-09 3/150 2025-04-09 17:28 by haomaier
[考研] 302求调剂 +4 毛毛毛, 2025-04-08 5/250 2025-04-09 11:01 by 毛毛毛,
信息提示
请填处理意见
濠㈣泛绉烽惁顖滄嫬閸愩劌顣�濠㈣泛绉烽惁顖滄嫬閸愩劌顣�
濠㈣泛绉烽惁顖滄嫬閸愩劌顣�濞存粌鐬煎ǎ顕€鎯嶉敓锟�闁兼儳鍢茶ぐ鍥ㄥ緞瀹ュ牏妲� & 閻犲鍟晶蹇涘箰閸パ冪闁挎冻鎷�鐎甸偊鍠曟穱濠囧箥椤愶紕顏遍柟纰夋嫹
闁哥喎妫楁慨鐔煎触閸粎绋�闁哥喎妫楁慨鐔煎触閸粎绋�
闁哥喎妫楁慨鐔煎触閸粎绋�濞存粌鐬煎ǎ顕€鎯嶉敓锟�闁哥喎妫楁慨鐔煎触閸粎绋婇柛婧垮姀椤曪拷鐎甸偊鍠曟穱濠囧箥椤愶紕顏遍柟纰夋嫹(濠㈣泛娲﹂弫鐐烘晬濮橆剚娅岄柛鏂衡偓铏€ゅù锝忔嫹)
閻庢冻闄勫﹢鎶界嵁閼艰埖褰�閻庢冻闄勫﹢鎶界嵁閼艰埖褰�
閻庢冻闄勫﹢鎶界嵁閼艰埖褰�濞存粌鐬煎ǎ顕€鎯嶉敓锟�鐎甸偊鍠曟穱濠囧箥椤愶紕顏遍柟纰夋嫹濞存粌妫滆闁哄牃鍋撻柡鍌涙緲椤掔喖寮甸姘彙閻犳劧鎷�
閻犱胶鍎ら弸鍐╂綇閸涱収鍤�閻犱胶鍎ら弸鍐╂綇閸涱収鍤�
閻犱胶鍎ら弸鍐╂綇閸涱収鍤�濞存粌鐬煎ǎ顕€鎯嶉敓锟�鐎甸偊鍠曟穱濠囧箥椤愶紕顏遍柟纰夋嫹婵烇綀顕ф慨鐐烘嚀娴e摜鐟庨柨娑樿嫰閸樸倗鎷归悷鐗堝皰閻犲浄鎷�
閻㈩垰鐗愮粊顐﹀箯濞戞ê鐝�閻㈩垰鐗愮粊顐﹀箯濞戞ê鐝�