首页/文章/ 详情

【JY】基于Matlab的双线性滞回代码编写教程

3年前浏览2334

 单质点体系振动是最为简单的振动,通常在学习结构动力学中也是最开始学习这部分的知识和内容,这部分内容最为基础,也非常重要。回顾这几篇文章:

【JY】结构动力学初步-单质点结构的瞬态动力学分析

【JY】基于Ramberg-Osgood本构模型的双线性计算分析

【JY】SDOF计算教学软件开发应用分享

    很多小伙伴在后台以及文章下方留言,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)



理论科普求解技术代码&命令振动其他软件MATLAB
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2021-08-19
最近编辑:3年前
建源之光
博士 | 高级技术经理 个人主页:jycmf.cn~
获赞 137粉丝 331文章 212课程 5
点赞
收藏
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈