lxpingguoxl 发表于 2011-5-3 16:09

matlab求解动力方程输入地震波的问题

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

meiyongyuandeze 发表于 2011-5-3 17:05

a虽然是一些类数值，那也应该是和时间有关的了，把a写成矩阵形式就是了！

lxpingguoxl 发表于 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的矩阵。我现在迷惑的是加速度的矩阵应该怎么表示呢，才能和质量矩阵相乘？弱智问题，谢谢指点！

meiyongyuandeze 发表于 2011-5-3 20:18

那要看你的那个自由度上要添加加速度激励啊！

lxpingguoxl 发表于 2011-5-3 21:22

每个自由度都加加速度激励

meiyongyuandeze 发表于 2011-5-3 21:46

哎，没有模型不好说清楚了，建议你带加速度激励一起建方程。猜测是不应该是对角阵呢，对角阵的元素就是各个自由度的质量！

lxpingguoxl 发表于 2011-5-4 09:08

谢谢啊，可是我是初学者，不明白你说的意思。现在质量矩阵、刚度矩阵、阻尼矩阵都有了，都是9*9阶的，我就想知道等式右面加速度的矩阵应该怎么写呢？冒昧问一下，能不能加你的QQ？

hiv5 发表于 2011-5-4 21:52

加速度向量，每个时刻t对应一个加速度值，九阶的，九个自由度。每个自由度对应一个加速度值，9X1的向量，同一个向量的元素都相同。t+0.02然后又是一个每个元素都相同的加速度列向量…………

lxpingguoxl 发表于 2011-5-6 19:27

回复 8 # hiv5 的帖子

谢谢，现在终于弄明白了

yzsldj 发表于 2011-5-7 15:00

这个解出来的结果是什么？
页: [1]
查看完整版本: matlab求解动力方程输入地震波的问题