【JY】STKO助力OpenSEES系列：结构模态分析以及动力特性（MDOF与等效SDOF验证）
相关阅读：
【JY】从一根悬臂梁说起
【JY】主成分分析与振型分解
【JY】求解器之三维结构振型分析
文/心尘轩
网站/STKO OpenSees Software (asdeasoft.net)
欢迎关注STKO官方公众号,目前许可证正在免费开放中,欢迎大家申请
STKO模态分析视频教程
一、模型案例信息
文字描述：
结构模型为6层5跨的平面框架结构（这种框架形式多见于结构抗震设计理论的学术文章），跨度为4m, 层高为3m, 约束平面外的平动自由度以防止面外变形，为简化建模过程，模型中的所有梁柱截面均采用400*400*20mm的空心箱型截面，每个节点赋予点质量10t以考虑荷载和自重，底部为固定约束可见下图。详细的建模过程可以参见相应的视频教程。
二、振型形状对比
下图给出了abaqus 和STKO 结构在第一和第二振型形状的对比，结果表明两个软件几乎一致，且呈图均具有美观欣赏性。
三、结构动力特性对比
两个软件对比的结构主要动力特性，其中包括：特征值，圆频率，模态有效质量，模态参与系数，如下系列表所示。根据对比结果可知，两个软件分析的结果几乎一致。由有效模态质量表格可知，第一阶有效模态占总质量为83%，所以对于这样的6层结构，我们可以说他的结构动力行为由其一阶模态控制，如果考虑高阶模态影响，则前两阶模态足以较精确地表达结构的动力响应。
四、欧拉梁和铁木辛柯梁对模态分析的影响
第四点是一个细节问题，总有同学问我，为什么Abaqus 和OpenSEES 做出来的结果对不上。实际上，是我们并没有考究我们的单元选择。在Abaqus 中，若不做特殊设置，默认的梁单元是考虑剪切变形的铁木辛柯梁，而在OpenSEES中默认的elastic beam column 是欧拉梁。欧拉梁因为不考虑剪切变形则会使得结构稍稍偏刚，这一点可以从下表的对比可以看出，欧拉梁的模型的特征值和圆频率均大于铁木辛柯梁。因而在这一节，一并回复询问者，做模型一定要细心，每一个参数的设置都会导致结构的悬殊。所以在Abaqus请明确使用的是B31 还是B33，在OpenSEES中可以使用section aggregator 进行添加剪切刚度以考虑剪切变形，总之模型需要保持一致性。
五、几何模型导入的快速建模
第五点内容，可参见文首视频教程，在这里不再赘述。
六、MDOF 和SDOF 等效合理性的验证
结构多自由度等效为单自由度的理论是进行结构抗震设计的一个重要内容。通过动力学的知识，我们可知，结构在弹性状态下，模态之间是正交互不影响的，于是通过对多自由度体系的推导如下：
我们可以得出以下结论： 地震输入给结构的总能量等于各个模态等效单自由度体系相应能量的之和。因而，我们可以将复杂的MDOF 问题转化为具有相应有效模态质量的单自由度问题。以双自由度系统而言，如下图：
这个双自由度系统的地震响应就等效为为两个分别具有有效模态质量的单自由度体系的叠加，也就是振型时程分析的原理。
在这里我们探讨一个有趣的问题，我们能否通过对结构的模态分析就可以根据此结构的动力特性求得结构在相应振型荷载分布下的结构初始刚度，这从原理上是可以实现的。通过模态分析，可以得出结构的有效模态质量和结构的模态圆频率，也就求出结构在相应的模态下的刚度。那我们就尝试下吧，对象为结构一阶模态下的结构刚度。操作如下：
首先提取结构一阶模态对应的振型，按照此振型分布模型对结构施加水平分布力，并进行模态推覆分析，如下图所示。
结构一阶模态推覆得到的等效单自由度的刚度和通过模态分析得到的刚度计算如下图。两者高度吻合，说明从刚度而言，MDOF 和相应模态等效的SDOF 理论是合理的。实际上，我们这节开头也推导了MDOF和相应模态等效的SDOF 从能量角度理论也是精确的。这些视角的分析和验证可以帮助我们进一步深入理解动力学的相关知识。
七、MDOF 和SDOF 等效在弹塑性状态的延申
实际上这个问题到目前为止依然没有被解决，结构一旦进入弹塑性后，结构的各阶振型是耦联的，并不再维持弹性状态下的正交关系。因此在弹性状态下的MDOF等效为相应SDOF能否继续延申至结构的弹塑性状态，是一个值得思考的问题。前辈们也付出较多的努力和探索。直到一个知名的权威专家出来说明：我们承认结构在进入弹塑性状态后结构各个振型存在耦联，也即某一个振型会对另一个振型产生影响，但我们可以假定如果我们继续按照结构弹性状态下振型理论去计算结构弹塑性的反应响应，如果得到的结构响应和结构实际的响应相比，精度满足工程需求，那么，我们就有理由相信，即使结构进入弹塑性状态，振型之间存在耦联，但这种耦联是比较弱的，忽略它并沿用结构弹性状态下的振型理论得到的结构响应是足以精确的描述结构的行为。正是在这样的指导思想下，结构的设计理论进一步得到了发展，比如MPA，以及随后的基于塑性的设计理论，也包括时下基于概率的结构性能设计理论的等等，而提出这样的假定的人就是知名结构学家乔普拉。也非常推荐大家阅读乔普拉关于结构动力学的专著，该专著中文版由谢礼立院士主持翻译。《DYNAMICS OF STRUCTURES Theory and Applications to Earthquake Engineering》。
公众号内容回顾：
关于STKO:
1.STKO for OpenSEES 安装教程
2.STKO for OpenSEES 免费许可证申请指南（修改版）
3.研究生STKO免费许可证申请邮箱范例
4.导师STKO免费许可证申请邮箱范例
5.无需TCL编程能力，STKO带你轻松玩转OpenSEES
6.STKO助力OpenSEES系列：自复位支撑框架静力循环pushover分析
7.STKO助力OpenSEES系列：平面多层多跨混凝土框架静力循环pushover分析
关于Abaqus:
1.Abaqus 基于python的plugin 开发，生成常用建模部件，可以辅助我们快速建模
2.Abaqus基于python的悬臂梁参数化分析（基础）
STKO 经典微观案例分享：
01：钢筋混凝土节点考虑钢筋滑移的微观模拟
02 钢筋混凝土节点的微观模拟
03 砌体结构的微观模拟
往期精彩
#性能分析
【JY】基于性能的抗震设计浅析（一）
【JY】基于性能的抗震设计浅析（二）
【JY】浅析消能附加阻尼比
【JY】近断层结构设计策略分析与讨论
【JY】浅析各动力求解算法及其算法数值阻尼(人工阻尼)
理念
【JY|体系】结构概念设计之(结构体系概念)
【JY|理念】结构概念设计之(设计理念进展)
【JY】有限单元分析的常见问题及单元选择
【JY】结构动力学之显隐式
【JY】浅谈结构设计
【JY】浅谈混凝土损伤模型及Abaqus中CDP的应用
【JY】结构概念设计之（隔震概念设计）
【JY】结构概念设计之（消能减震黏滞阻尼器）
#概念机理
【JY】基于Ramberg-Osgood本构模型的双线性计算分析
【JY】结构动力学初步-单质点结构的瞬态动力学分析
【JY】从一根悬臂梁说起
【JY】反应谱的详解与介绍
【JY】结构瑞利阻尼与经济订货模型
【JY】主成分分析与振型分解
【JY】浅谈结构多点激励之概念机理（上）
【JY】浅谈结构多点激励之分析方法（下）
【JY】板壳单元的分析详解
【JY】橡胶支座的简述和其力学性能计算
【JY】振型求解之子空间迭代
【JY】橡胶支座精细化模拟与有限元分析注意要点
【JY】推开土木工程振型求解之兰索斯法(Lanczos法)的大门
【JY】基于OpenSees和Sap2000静力动力计算案例分析
【JY】建筑结构施加地震波的方法与理论机理
【JY】消能减震黏滞阻尼器的力学原理与应用
#软件讨论
【JY】复合材料分析利器—内聚力单元
【JY】SDOF计算教学软件开发应用分享
【JY】Abaqus案例—天然橡胶隔震支座竖(轴)向力学性能
【JY】Abaqus6.14-4如何关联fortran？
【JY】如何利用python来编写GUI？
【JY】如何解决MATLAB GUI编程软件移植运行问题？
【JY】浅谈结构分析与设计软件
【JY|STR】求解器之三维结构振型分析
【JY】SignalData软件开发应用分享
【JY】基于Matlab的双线性滞回代码编写教程
#其他
【JY】位移角还是有害位移角？
【JY】如何利用python来编写GUI？
【JY】今日科普之BIM
~关注未来更精彩~