【BiliBili】 【知乎】 【QQ空间】 【微博】


博学慎思,坚持创新,求知若飢,虛心若愚                       关键字     类型    
  • 博主
  • 全部
  • 日志
  • 研究
  • 论文
  • 自编程序
  • 工程
  • Dino笔记
  • DinoBox
  • OpenSEES
  • Dino公开课
  • 每周挑战
  • Dino爱编程
  • 留言板

DinoBox

ETE_Frame3D
结构分析程序
ETE_Frame3D
2022-01-19
ETE-section截面分析程序
ETE-section
2021-08-19
ETO_OpenSEES前后处理程序
程序通过邮件申请。
★注意,软件不能下载,是通过邮箱发送软件的。

申请ETO时,注意EMAIL 的写法:(格式一定要对,否则电脑收不到)
邮件题目: ETO (版本号) 下载
注册人: (真实姓名)
学校/单位: (大学或工作单位)
邮件回复,大约需要1~4天的时间。
ETO_程序
2014-03-26
上一篇 | 下一篇
Dino笔记 订阅所有Dino笔记的日志

[Dino笔记]多自由度体系模型时程分析的逐步积分法

作者:dinochen1983 日期:2022-01-08

字体大小: 小 中 大

【Dino结构笔记】多自由度体系模型时程分析的逐步积分法


【程序下载】程序: MDOF逐步积分法小程序mdof_hist2.rar
【附件下载】附件: MDOF时程分析的学习资料PDF
【附件下载】附件: SAP2000模型及EXCEL表格

上一篇有关于结构动力学的笔记的内容关于采用振型叠加法求解结构在地震波作用下的时程响应,也简称为频域方法,那么对应的就有时域方法,也就是地震动时程分析的逐步积分法,举例来说有Newmark法与Wilson法。

小记:我参考了几本不同的《结构动力学》教材,包括克拉夫《结构动力学》,还有柴田明德的《结构抗震分析》、RoyRCraig.Jr著的《结构动力学》,这三本经典教材,其中克拉夫《结构动力学》太过省略,适全高端读者,而RoyRCraig.Jr著版本,简单略过,柴田明德的公式列得不错,多自由度的例子不多。最后我在香港图书馆找到一本书《Computer-Aided Mechanical Design and Analysis》(V.Ramamurti著),这本书真的不错,公式列得好的同时,还增加了多自由度的计算实例,程序按这本书的公式编写。


无论是地震动时程分析还是风振时程分析,都有频域法与时域法,频域法速度快,但是不能解决结构变刚度(弹塑性分析)变阻尼(带阻尼器)的情况的结构动力分析(FNA方法除外)。频域法是通过把结构拆成多个模态的组合,对每个模态进行单自由度求解,叠加后得到结构的响应,而时域方法就是把时间离散为一个个时间段,如0.02s,然后每个时间更新结构的动力荷载及动力刚度,通过每一步的静力求解得到下一步的响应,所以叫做逐步积分法,上一步的计算结果影响下一步的计算结果。
注意,时域分析不需要进行模态分析也能进行,但程序中进行模态分析是因为要得到第1周期与第2周期,用来算瑞利阻尼。时域分析因为没有把模型拆成单自由度,所以阻尼不能指定给模态,必须组装整体阻尼矩阵,因此需要通地前2阶的周期计算出阻尼矩阵。

 上述程序可帮助同学学习多自由度体系的逐步积分法,整个计算的细节通过文本以计算书的形式展示出来,如展示质量矩阵,刚度矩阵及阻尼矩阵,还有逐步法的计算参数。程序提供两种算法Wilson-theta与Newmark法

本篇笔记通过一个4自由度的计算实例讲解如何通过逐步积分法来求解多自由度体系(层剪切模型)的地震动时程分析。这里面有几个概念要搞清楚。
(1) 模态分析只为了求第1第2周期,用于计算阻尼矩阵
(2) 逐步积分法要形成结构的刚度矩阵,质量矩阵与阻尼矩阵
(3) 通地Newmark, Wilson-theta法可以把结构动力分析变成静力分析,每步求解为静力求解,因此需要计算等效刚度矩阵与等效荷载向量。
(4) 静力求解得到下一步(或延长步)的位移向量,通过位移向量,根据不同的逐步积分法公式就可以求得下一步的速度向量与加速度向量。
(5) 弹性结构不存在变刚度,所以每一步只需进行一步静力求解,不需要迭代。
学习之前的准备:
(1) 层模型的简化与建模
http://www.dinochen.com/article.asp?id=260
(2) 层模型的模态分析求解方法
http://www.dinochen.com/article.asp?id=376
(3) 多自由度体系模型时程分析的振型叠加法
http://www.dinochen.com/article.asp?id=405
在学习多自由度(层剪切模型)的地震动时程分析法的原理之前,需要了解如何把多层结构简化为层剪切模型及在ETABS或其它软件如何建立剪切模型,这方面的内容可以参考以下笔记。

多自由度体系MDOF的地震动时程分析的逐步积分法的主要原理就是:通过不同的积分方法,把动力分析转换为静力分析。
(4) 关于如何得到简化模型及相关矩阵,可以看《层模型的简化与建模》
http://www.dinochen.com/article.asp?id=260
 
求解过程如下:
(1) 形成结构的刚度矩阵,质量矩阵,如上图所示。
 
(2) 通过模态分析求得结构的第1与第2周期,再通过周期求解阻尼矩阵,采用瑞利阻尼,计算过程如下。计算公式参考克拉夫《结构动力学》书上公式12-41(第185页)。
 

(3) 逐步积分法,以Wilson-theta法为例整个计算过程如下所示:


详细资料可下载PDF。
 (a) 形成刚度矩阵,质量矩阵,阻尼矩阵
(b) 计算Wilson-theta法的参数c0~c3
(c) 计算等效刚度矩阵
(d) 根据上一步的位移,速度与加速度计算等效荷载向量
(e) 静力求解当前步(或延长步)的位移向量
(f) 根据不同的积分方法求解下一步的位移,速度与加速度向量,用于生成下一步的等效荷载向量,迭代进行至全部荷载步算完。

(4) 最终求得的结构位移时程分下图所示,与SAP2000进行对比非常吻合。SAP2000模型可以附件下载。时程分析对比,可见EXCEL表格。
如果采用SAP2000进行逐步积分法,瑞利阻尼的设置,如下图所示。
 
在SAP2000中设置瑞利阻尼的设置方法

 
Wilson-theta方法程序与SAP2000的结果对比
 
Newmark方法程序与SAP2000的结果对比

其中,Wilson-theta, theta = 1.4, Newmark方法,Gamma = 0.5, beta = 0.25
以上就是本节笔记的全部内容。

====本节完====

(微信公众号用户,点击【阅读原文】后下载)


更多结构动力学的东西,可以关注我的微信公众号: DINOSTRU




[本日志由 dinochen1983 于 2022-01-08 10:58 PM 编辑]
文章来自: 本站原创
引用通告: 查看所有引用 | 我要引用此文章
Tags: 动力学动力学
评论: 0 | 引用: 0 | 查看次数: 3348
发表评论
昵 称:
密 码: 游客发言不需要密码.
验证码: 验证码
Dino:如果看不到验证码,再刷几下就可以了
内 容:
选 项:
虽然发表评论不用注册,但是为了保护您的发言权,建议您注册帐号.
字数限制 200 字 | UBB代码 关闭 | [img]标签 关闭

信息(Information)

DinoChen 视频频道

My GitHub Page



扫描以下二维码,加入公众号!


dinochen.com需要你的支持


日历(Calendar)

2022年7月
  • 日
  • 一
  • 二
  • 三
  • 四
  • 五
  • 六
  • 26
  • 27
  • 28
  • 29
  • 30
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6

查找(Search)

关键字
 
类 型

文章(Recent Article)

[OpenSEES]OpenSEES实例教程[汇总...[Dino笔记]粒子群算法在桁架结构中的应用[获奖]WSP亚洲2021年卓越奖获得一等奖与创新...[Dino笔记]基于遗传算法对剪力墙结构优化设计方...[Dino笔记]基于OpenSEES虚功原理的钢结...[OpenSEES] 新版ETO更新及相关下载方法...[新闻]祝韩老师荣获广东省科技进步一等奖[程序]ETE-section功能更新:可修改划分...[App开发]采用Swfit语言编制经典的推箱子游...[程序]ETE-section更新_考虑混凝土保护...

评论(Comments)

陈博士您好,这个小工具好像有点问题,有效翼板宽度b...打开的时候显示“你没有许可证”,怎么弄啊陈博士您好,这个是需要购买的软件还是可以试用的怎么查看文档计算过程呢,陈博士一直显示出现一个错误字符无法计算是什么原因呢?为什么一直出错计算不了?陈博士,能共享下文件嘛?第二版没有光盘陈博士您好,我换了好几台电脑ETE还是需要许可,陈博士您好,现在程序还是打不开,需要许可请问如何获得会员许可证?

留言(GuestBook)

Thanks regarding furnishi...陈博士您好,想请问一下零长度单元用在桥梁的支座上,...陈博,新版ETO下载以后界面打不开是什么原因呢,旧...陈博士您好,请问可以出一期画易损性曲线的教程吗陈博士您好,我在学习您4.33侧向多自由度简化模型...陈博士您好,計算wind2019的program ...陈博士,您好,我下载使用了您编写的压型钢板组合楼板...陈博士您好,我是opensee初学者,想向您请教两...陈博士你好,请问开洞混凝土截面,ETE-secti...陈博士您好,感谢您的回复,我的刚性隔板设置了转动惯...

TagBar

Loading...

用户面板(User Panel)

登录用户注册

统计(Statistics)

日志: 408 篇
评论: 893 个
留言: 413 个
访问: 14023740 次
在线: 9 人
建站时间: 2008-02-08

Powered By PJBlog3 V2.7 Build 05 CopyRight 2005 - 2011, 陈学伟的博客

Processed in 0.281250 second(s) , 6 queries 粤ICP备15071789号-1