wangli20060913 发表于 2006-10-1 08:23

不平衡响应分析

各位高手：
      我在学习过程中遇到这样一个问题：
      JX''+KX=T
      式中，J是一个十七阶三对角质量矩阵，即
J=（J1 J2 0 0...0;J2 J3 J4 0...0;0 J3 J4 J5 0 0 ... 0;0 0 J4 J5 J6 0 0..0;...;0...0 Jn-1 Jn）
K是一个十七阶三对角刚度矩阵，即
K=(K1 K2 0 ...0;K2 K3 K4 0 ...0;0 K3 K4 K5 0...0;0 0 ...0;0 0 ... 0 Kn-1 Kn)
T是一个十七阶外界激励力列向量，即
T=（0 0 0 0 T5 0 T7 0 T9 0 T10 0 T11 0 0 0 0）
每个激励力T=P1sinwt+P2sin2wt+P3sin3wt+P4sin4wt
如何求解X的响应呢？如果用matlab,如何实现呢？
敬请各位高手指点迷津，小弟不甚感激！！

MVH 发表于 2006-10-3 22:06

可以将系统转化成34个状态方程，如果用matlab的话可以尝试一下涌ode23s看看

wangli20060913 发表于 2006-10-9 07:51

谢谢MVH!!!
能否具体解释一下呢？谢谢！！！

yang95478 发表于 2006-10-12 20:10

或者也可以先按照MVH所说将方程降阶，在Vc++中编制龙格-库塔法或基尔算法的程序进行求解，我现在就是这么做的！

wangli20060913 发表于 2006-10-14 14:50

谢谢yang95478!!!
页: [1]
查看完整版本: 不平衡响应分析