很多小伙伴在后台以及文章下方留言,Matlab的动力学分析怎么编写?JYChen带大家详细的过一遍。但首先,请读者先移步阅读上面三篇文章,知道双线性计算分析原理和动力计算方法,然后我们准备发车。
本程序编写主要是单质点计算(含双线性+黏滞阻尼编写)
图 带水平粘滞阻尼器计算的双线性
图 不带水平粘滞阻尼器计算的双线性
好了,发车了~
代码讲解:
第5行:对输出界面清空,对计算保留的数据清空。
第7~10行:设置计算所需要的动力参数。
第11~12行:读取txt文本的地震波。(注:也可以其他方式读取)
第14~15行:黏滞阻尼器所需要的计算参数。
第17~19行:设置双线性滞回的关键参数。
第20行:计算系统阻尼比。
第22~23行:计算前处理所需要整理的变量。
第25行:基于地震波步进迭代开始。
第26~34行:双线性滞回的未屈服阶段判断进行确定刚度。
第35~41行:双线性滞回的正向屈服阶段判断进行确定刚度。
第42~48行:双线性滞回的反向屈服阶段判断进行确定刚度。
第49~52行:各计算迭代力提前整理与准备计算。
第53~66行:进行Newmark-β的计算分析。(详情请阅读任意一本结构动力学)
第67~69行:结果进行整理储存。
第70~72行:结果进行整理,得到了相对加速度Acc,绝对加速度Aabs,位移Dish,恢复力Fsh。
(注意:前面计算单位需保证一致,建议取:N、m、s。由于最后化单位,加速度单位为:m/s2,位移Dish单位为mm,恢复力Fsh单位为kN)
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删