[地震][动力学][Dynamics][MATLAB] 将阻尼矩阵的非对角线元素取为0计算结果会怎么样?
实干、实践、积累、思考、创新。 如题,一看是一个莫名其妙的想法。实际上也是一个错误的想法,不过,当时脑子一热,就测算一下。具体看看结果。 模型: 20层的剪切层模型。 阻尼矩阵: (1)模态阻尼,20阶振型计算 (2)模态阻尼,20阶振型计算,把非对角线元素取为0。 采用MATLAB编程,采用Newmark-β积分法进行弹性时程分析,两种阻尼模型的计算结果对比如下。 其中,参考阻尼为完整的模态阻尼,对比阻尼为去掉阻尼矩阵对角线元素后的矩阵。 由图可见,采用仅保留对角线元素的阻尼矩阵,结构的位移及剪力响应远小于完整的阻尼矩阵,而楼层加速度响应似乎相差不大!!十分诧异!! 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号