[地震][动力学][Dynamics][MATLAB] 将阻尼矩阵的非对角线元素取为0计算结果会怎么样?

实干、实践、积累、思考、创新。


如题,一看是一个莫名其妙的想法。实际上也是一个错误的想法,不过,当时脑子一热,就测算一下。具体看看结果。

模型:

20层的剪切层模型。

阻尼矩阵:

(1)模态阻尼,20阶振型计算

(2)模态阻尼,20阶振型计算,把非对角线元素取为0。

采用MATLAB编程,采用Newmark-β积分法进行弹性时程分析,两种阻尼模型的计算结果对比如下。

其中,参考阻尼为完整的模态阻尼,对比阻尼为去掉阻尼矩阵对角线元素后的矩阵。

由图可见,采用仅保留对角线元素的阻尼矩阵,结构的位移及剪力响应远小于完整的阻尼矩阵,而楼层加速度响应似乎相差不大!!十分诧异!!


You already voted!

  • 微信公众号 ( Wechat Subscription)

WeChat_QRCode

欢迎关注 “结构之旅” 微信公众号

2 thoughts on “[地震][动力学][Dynamics][MATLAB] 将阻尼矩阵的非对角线元素取为0计算结果会怎么样?

  1. Reply

    博士您好,能分享一下Matlab建立非线性时程的代码吗,我最近在用matlab编写这个程序遇到了很多困难

    • CJD Post authorReply

      看这个连接:http://www.jdcui.com/?page_id=16529,新书《结构地震反应分析 编程与软件应用》对于单自由度、多自由度飞弹性及弹塑性的动力分析都给出例子和全套代码,建议直接看看就知道了。

Leave a Reply to CJD Cancel reply

Your email address will not be published. Required fields are marked *

This site uses Akismet to reduce spam. Learn how your comment data is processed.