声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1845|回复: 9

[综合讨论] matlab求解动力方程输入地震波的问题

[复制链接]
发表于 2011-5-3 16:09 | 显示全部楼层 |阅读模式

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

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

x
动力方程M*d2X+C*dX+K*X=M*a其中:M,C,K是n阶方阵,X是n维列向量,d2X表示X的二阶导数,dX表示X的一阶导数。式中的a是地震加速度,是一系列已知的具体的数值,我应该怎么在matlab输入呢?或者说a的向量形式应该是什么样的?

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2011-5-3 17:05 | 显示全部楼层
a虽然是一些类数值,那也应该是和时间有关的了,把a写成矩阵形式就是了!
 楼主| 发表于 2011-5-3 18:22 | 显示全部楼层
这是时间间隔0.02S的地震波加速度
-9.33E-03
-7.20E-02
-6.73E-02
-5.87E-02
-6.33E-02
-8.00E-02
-9.47E-02
-8.53E-02
-7.33E-02
-5.67E-02
-5.67E-02
-8.73E-02
-1.17E-01
-1.29E-01
-1.08E-01
-9.60E-02
-7.20E-02
-5.47E-02
-2.80E-02
-4.40E-02
-8.73E-02
-1.27E-01
-1.31E-01
-4.40E-02
2.00E-02
9.40E-02
-3.27E-02
-8.53E-02
-9.60E-02
-1.35E-01
-1.73E-01
-2.17E-01
-2.04E-01
-1.15E-01
-1.31E-01
-1.09E-01
-1.09E-01
-4.47E-02
1.67E-02
1.00E-01
1.57E-01
1.68E-01
2.24E-01
3.09E-01
3.28E-01
2.79E-01
2.39E-01
1.81E-01
1.57E-01
2.26E-01
2.75E-01
质量矩阵是9*9的矩阵。我现在迷惑的是加速度的矩阵应该怎么表示呢,才能和质量矩阵相乘?弱智问题,谢谢指点!
发表于 2011-5-3 20:18 | 显示全部楼层
那要看你的那个自由度上要添加加速度激励啊!
 楼主| 发表于 2011-5-3 21:22 | 显示全部楼层
每个自由度都加加速度激励
发表于 2011-5-3 21:46 | 显示全部楼层
哎,没有模型不好说清楚了,建议你带加速度激励一起建方程。猜测是不应该是对角阵呢,对角阵的元素就是各个自由度的质量!
 楼主| 发表于 2011-5-4 09:08 | 显示全部楼层
谢谢啊,可是我是初学者,不明白你说的意思。现在质量矩阵、刚度矩阵、阻尼矩阵都有了,都是9*9阶的,我就想知道等式右面加速度的矩阵应该怎么写呢?冒昧问一下,能不能加你的QQ?
发表于 2011-5-4 21:52 | 显示全部楼层
加速度向量,每个时刻t对应一个加速度值,九阶的,九个自由度。每个自由度对应一个加速度值,9X1的向量,同一个向量的元素都相同。t+0.02然后又是一个每个元素都相同的加速度列向量…………
 楼主| 发表于 2011-5-6 19:27 | 显示全部楼层
回复 8 # hiv5 的帖子

谢谢,现在终于弄明白了
发表于 2011-5-7 15:00 | 显示全部楼层
这个解出来的结果是什么?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-28 21:35 , Processed in 0.056720 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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