[Dino工程]深圳金地大百汇中心（375.6 m）

工程特点描述：带加强层的混合结构，钢管柱，立面曲线布置斜柱，
楼板大开洞，立面混合体系，伸臂桁架与腰桁架，徐变计算，皇冠钢结构

WSP服务类别：结构设计顾问(初步设计)、机电顾问
开发商：金地大百汇集团
合作单位：筑博设计股份有限公司
建筑方案公司：KPF

金地大百汇中心项目位于深圳市福田区中心地段，地上建筑面积约为16.5万m2,其中包括一个集办公、酒店公寓及商业于一体的375米主塔楼，并设有4层地下室[1]。主塔楼的屋面高度为338m，地上72层，顶部有皇冠钢结构，建筑效果图如图1所示。金地大百汇中心项目位于深圳市福田区中心地段，地上建筑面积约为16.5万m2,其中包括一个集办公、酒店公寓及商业于一体的375米主塔楼，并设有4层地下室。主塔楼的屋面高度为338m，地上72层，顶部有皇冠钢结构，建筑效果图如图所示。




CTUBH数据

塔楼结构平面呈带圆弧的矩形，标准层结构平面如图所示。

结构体系

工程位于强风荷载地区，风荷载为整体性能指标（层间位移角限制）的控制因素，结合建筑净高要求及纺锤形的立面特征，主塔楼采用钢框架+钢筋混凝土核心筒混合结构体系，其中钢框架由钢管混凝土柱与型钢梁组成。由于立面呈纺锤形，每层外框柱均为不同角度的斜柱，在重力荷载作用下，斜柱对框架梁产生的水平拉力，宜采用钢梁来抵抗这部分拉力。与钢筋混凝土框架相比，采用钢结构框架施工速度快，结构自重轻及延性好，提高斜柱建造的准确性。 

竖向抗侧力体系中低区采用框架核心筒体系，高区存在大面积楼板开洞，抗侧力基本由筒体承担，部分外框柱取消，为抵抗风荷载增设幕墙骨架。为减少地震对抗风的幕墙骨架的影响，幕墙钢骨架的底部采用铰接处理。塔楼最顶部为皇冠钢结构，立面结构布置如图所示。楼板体系采用钢梁+压型钢板组合楼板体系。部分连梁构件由于风荷载的承载力需要，设置成型钢混凝土连梁。

由于需要满足结构在风荷载作用下层间位移角满足限值要求，结构需要设置加强层，加强层设置三排伸臂桁架及X方向的腰桁架，总共设置2个加强层，分别设置于23层与51层。加强层的桁架三维布置如图所示。 

施工现场与最后建成的照片！

我与伸臂桁架与腰桁架的合影


2020年5月的照片

[Dino工程]深圳华润前海中心（291.5 m）

工程特点描述：无加强层的混凝土结构，叠合柱，结构优化初设的程序开发

WSP服务类别：结构设计顾问(初步设计)
开发商：华润集团
合作单位：广东省院深圳分院
建筑方案公司：美国GP
机电顾问单位：Meinhardt

华润前海中心项目位于深圳市南山区月亮湾大道西侧。项目包括三栋建筑面积总计340000m2的办公楼（T1、T2、T5）、一栋建筑面积56000m2的公寓（T4）、一栋建筑面积50000m2的酒店（T3）。
以下重点介绍主塔楼T1的结构设计。


T1塔楼设有4层地下室，地上62层，室外地面起算房屋高度为276.1m。考虑到项目存在下层广场，嵌固端设置在B1，从嵌固端起算，计算高度为298.5米。塔楼平面尺寸49.5mx49.5m。核心筒为矩形，塔楼层1~51平面尺寸为28.3mx26.8m,层52以上核心筒收进20.6mx25.8m。塔楼层5、20和37为机电和避难层，层高均为5.1m；办公标准层层高4.2米。

T1塔楼上部结构的抗侧力体系为框架-核心筒结构；高宽比为5.58（小于7），由于本结构高宽比合理，无需设置伸臂桁架或者腰桁架即可满足抗侧刚度要求。结构低区及高区的结构平面布置图如图所示。

低区外框柱采用钢管叠合柱构件，钢管内混凝土强度等级为C60，钢材强度等级为Q390，从嵌固端到31层，柱截面由1.4m×1.4m逐步减小至1.3m×1.3m；高区外框柱采用钢筋混凝土柱，从32层至顶层，混凝土强度等级由C60逐步减小至C40，柱截面由1.3m×1.3m逐步减小至0.8m×0.8m。

核心筒材料主要采用钢筋混凝土，从嵌固端到顶层，混凝土强度等级由C60逐步减小至C40，厚度由1.3m逐步减小至0.4m。此外，在51层以上，核心筒的尺寸由28.3mx26.8m缩小为20.6mx25.8m。底部加强部位为负一层至6层，在负一层至10层的核心筒剪力墙内设置型钢，以满足轴压比和中大震受拉受剪的要求。

广东省地区与香港地区建筑结构风荷载规范的比较研究

小编：辛展文，指导老师：陈学伟、杨易

【1】研究背景与研究思路

风荷载设计是建筑结构设计中的一项重要内容，对沿海超高层建筑尤其重要。高层建筑抗风设计方法：规范设计、风洞试验、数值风洞模拟仿真，其中规范是基准。


2019年，香港屋宇署新修订了风荷载规范。广东省风荷载设计规范按照《建筑结构荷载规范》（2012版）实行。掌握这两地的风荷载规范，对粤港两地的工程项目非常重要。在这样的背景下，确定了本研究的选题。

研究的思路如下：

【2】两地规范对比

基本计算公式对比如下：

可以看出两地顺风向风荷载的计算形式相似，但在横风方向上，广东规范通过计算各楼层风荷载来衡量风作用，而香港规范则通过基底弯矩来衡量。计算得到横风基底弯矩后，香港规范需要根据横风基底弯矩对顺风向风荷载进行修正，广东规范则没有这部分内容。对于扭转风荷载的计算，香港规范采用一个更简单的计算模型。它计算风振加速度的公式也较为简单。

值得注意的是，广东规范只需要计算两个方向（X、Y）的风作用，而香港规范需要验算四个方向（±X1、±X2），而且香港规范可以计算任意风向

部分计算参数对比如下：

可以看出两套规范下的计算参数在取值规定与测量方法上也各有差异。计算公式和计算参数的差异性最终导致两地风荷载规范计算结果有较大的差别。

【3】风荷载计算软件编写

在掌握了规范条文后，使用面向对象编程语言DELPHI针对两地规范分别编写风载计算软件。该计算软件可以在不建模的情况下，通过输入建筑尺寸参数、结构动力特性参数、风气候特征等参数，对建筑结构进行风荷载分析。

【4】计算软件准确性检验

利用编写得到的两个计算软件，对一67层的超高层建筑进行风荷载计算。该建筑总高度为295.1m，建筑宽度为46.9m，进深为46.9m，建筑平面体型为矩形。

广东风荷载计算软件计算得到的电算结果通过PKPM进行检验，对比结果显示PKPM与自编软件计算得到的顺风向风荷载相差1.86%，横风向风荷载只相差0.34%。

香港风荷载计算软件通过手算屋面层风荷载作用的方式进行校核，最后得到的手算结果也和电算结果相符合。

【5】实例计算结果分析

对两个软件计算得到的数据进行对比，可以看出香港规范的计算结果比广东规范的要大很多。它们之间的偏差会随着层高的增加而增大，在顶层的偏差达到最大值，并且可以看出，Y方向上的偏差比X方向上的小。

为了更深入地理解这些计算结果的差异性，该研究对计算过程中涉及到的一个重要参数，不同高度处的风压值，进行计算分析，取最接近香港的城市：深圳，所采用的基本风压0.90（对应回归周期100年），为分析对象，计算得到的风压值也比香港规范的要小得多，其中最接近香港规范的是A类地貌。

【6】总结

风荷载规范是建筑结构抗风设计的重要基准，涉及建筑结构的安全和经济。本文在深入研究和详细对比粤港两地规范基础上，依据规范中的风荷载计算模型，通过面向对象语言DELPHI编程，成功开发了两款风荷载计算软件。

通过对一栋超高层建筑工程实例的风荷载分析，验证了风荷载计算软件的准确性，计算精度达到商用软件PKPM标准。

通过规范对比和工程算例发现，两地风荷载规范的计算模型和基本参数存在差异；由此导致粤港临近地区，即使建筑体型和结构动力特性相似的建筑，其风荷载规范取值存在显著差异。如本例中，顺风向风荷载基底总剪力按照广东省规范计算比香港地区的小64.98%。

香港风荷载程序的介绍与使用：

bilibili网站：https://www.bilibili.com/video/bv1Az4y1Q7uy

【ETE专题】ETE提取YJK与ETABS的超限性能总指标表格

[小编：李明]

简说：开发一个基于YJK与ETABS的后处理程序，可以对YJK模型与ETABS模型提取结构性能的总指标，并自动生成WORD文档的报告。，程序采用DELPHI进行编写，自动生成曲线包括超限报告的WORD文档。

[下载讲义PDF]：点击下载此文件：ETE提取超限性能总指标表格之方法

以下介绍一下ETE通过读取YJK与ETABS的数据后自动生成超限报告的功能介绍

【第一步】在YJK进行操作，导出YJK的结构模型与设计分析结果

【第二步】在ETABS操作，导出ETABS的分析结果（如有）一般高层的结构需要第二个软件进行复核

【第三步】ETE读取YJK的模型数据，包括小震，中震，大震（等效弹性）的文件夹

【第四步】导入ETABS的风模型与地震模型

【第五步】设置文件夹存档的位置

【第六步】补充部分后处理的数量，如工况的名字

【第七步】一键生成WORD文档的报告。

以上就是ETE【中国版】关于超限报告总指标的生成的方法（基于YJK的数据）

有兴趣请关注我们的微信公众号:DINOSTRU

以下是软件的操作视频：[Bilibili]B站视频

https://www.bilibili.com/video/BV1yA411e7Me/

中国知网自动生成的参考文献格式（更新于2020-07-16）

[1]陈学伟,韩小雷,王响.基于网络的结构协同弹塑性分析方法的研究[J].计算力学学报,2012,29(02):224-229.

[2]陈学伟,韩小雷,孙思为.三种非线性梁柱单元的研究及单元开发[J].工程力学,2011,28(S1):5-11.

[3]陈学伟,韩小雷.剪力墙非线性宏观单元的研究与单元开发[J].工程力学,2011,28(05):111-116+123.

[4]陈学伟. 剪力墙结构构件变形指标的研究及计算平台开发[D].华南理工大学,2011.

[5]韩小雷,陈学伟,熊爱国,何伟球.盛德大厦改建工程抗震性能研究[J].建筑结构,2011,41(04):88-92.

[6]韩小雷,陈学伟,梁崇,唐剑秋.不同模型对足尺钢框架振动台试验模拟的影响[J].地震工程与工程振动,2010,30(06):71-77.

[7]陈学伟,韩小雷,孙思为. 三种非线性梁柱单元的研究及单元开发[C]. 中国力学学会结构工程专业委员会、山东建筑大学、中国力学学会《工程力学》编委会、清华大学土木工程系、清华大学水沙科学与水利水电工程国家重点实验室、清华大学土木工程安全与耐久教育部重点实验室、山东建筑大学建筑结构鉴定加固与改造省重点实验室.第19届全国结构工程学术会议论文集（第Ⅰ册）.中国力学学会结构工程专业委员会、山东建筑大学、中国力学学会《工程力学》编委会、清华大学土木工程系、清华大学水沙科学与水利水电工程国家重点实验室、清华大学土木工程安全与耐久教育部重点实验室、山东建筑大学建筑结构鉴定加固与改造省重点实验室:中国力学学会工程力学编辑部,2010:273-280.

[8]陈学伟,韩小雷,文健,黄俊杰,何伟球,毛贵牛.中洲二期观光门架结构设计[J].建筑结构,2010,40(10):44-50.

[9]韩小雷,陈学伟,顾希敏,林生逸,吴培烽,毛贵牛,何伟球.珠海天朗海峰超限高层建筑结构设计[J].建筑结构,2010,40(10):51-54.

[10]陈学伟,季静,吴培烽,罗凡,吴爽.E-Defense振动台试验预测性分析比赛的研究综述[J].世界地震工程,2010,26(03):175-181.

[11]邓国基,陈学伟,杨穗华,李稼轩.汶川县人民医院带粘滞阻尼器结构耗能减震设计[J].广东土木与建筑,2010,17(09):3-6+9.

[12]马玉宏,金建敏,韩小雷,谭平,陈学伟,沈朝勇,黄襄云,陈建秋,任珉.高层钢-混凝土组合门式结构消能减震体系模拟地震振动台试验研究[J].振动与冲击,2010,29(06):134-139+241.

[13]陈学伟,韩小雷,林生逸.基于宏观单元的结构非线性分析方法、算例及工程应用[J].工程力学,2010,27(S1):59-67.

[14]韩小雷,陈学伟,林生逸,何伟球,郑宜,吴培烽,毛贵牛.基于纤维模型的超高层钢筋混凝土结构弹塑性时程分析[J].建筑结构,2010,40(02):13-16.

[15]陈学伟,韩小雷,林生逸,吴培烽,何伟球.中洲中心二期结构抗震性能分析[J].建筑结构学报,2010,31(01):101-109.

[16]季静,陈学伟,何慧贤,韩小雷.基于概率和位移的框支剪力墙结构抗震性能评估方法[J].地震工程与工程振动,2009,29(06):108-114.

[17]韩小雷,陈学伟,林生逸. 基于宏观单元的结构非线性分析方法、算例及工程应用[C]. 中国力学学会结构工程专业委员会、广州大学土木工程学院、中国力学学会《工程力学》编委会、清华大学土木工程系、清华大学结构工程与振动重点实验室.第18届全国结构工程学术会议论文集第Ⅱ册.中国力学学会结构工程专业委员会、广州大学土木工程学院、中国力学学会《工程力学》编委会、清华大学土木工程系、清华大学结构工程与振动重点实验室:中国力学学会工程力学编辑部,2009:40-48.

[18]陈学伟,韩小雷,毛贵牛,季静.粘滞阻尼器在连体高层结构中的抗风减振效果[J].土木建筑与环境工程,2009,31(05):74-80.

[19]韩小雷,唐剑秋,黄艺燕,季静,陈学伟.钢管混凝土巨型斜交网格筒体结构非线性分析[J].地震工程与工程振动,2009,29(04):77-84.

[20]熊爱国,季静,陈学伟. 某造粒钢塔结构温度效应专项分析[C]. 天津大学.第九届全国现代结构工程学术研讨会论文集.天津大学:全国现代结构工程学术研讨会学术委员会,2009:1130-1135.

[21]蔡静敏,陈学伟.人行天桥振动时程分析方法的研究[J].钢结构,2009,24(07):16-20.

[22]韩小雷,陈学伟,毛贵牛,郑宜,季静.基于人群行走仿真的楼板振动分析方法及反应谱公式推导[J].建筑科学,2009,25(05):4-9.

[23]何伟球,陈学伟,韩小雷.基于风洞试验的超高层建筑结构风振时程分析[J].广东土木与建筑,2009(02):14-16+28.

[24]韩小雷,陈学伟,吴培烽,郑宜,杨志强.OpenSEES的剪力墙宏观单元的研究[J].世界地震工程,2008,24(04):76-81.

[25]韩小雷,陈学伟,郑宜,彭樵斌.足尺钢框架振动台试验及动力弹塑性数值模拟[J].地震工程与工程振动,2008,28(06):134-141.

[26]韩小雷,陈学伟,戴金华,郑宜,何伟球.基于OpenSEES的剪力墙低周往复试验的数值分析[J].华南理工大学学报(自然科学版),2008,36(12):7-12.

[27]陈学伟,韩小雷,郑宜,毛贵牛.高层建筑结构风振时程分析软件的技术研究[J].华中科技大学学报(城市科学版),2008(03):136-139+145.

[28]吴培烽,季静,陈学伟. 钢筋混凝土剪力墙结构非线性分析宏观模型浅述[C]. 天津大学.庆祝刘锡良教授八十华诞暨第八届全国现代结构工程学术研讨会论文集.天津大学:全国现代结构工程学术研讨会学术委员会,2008:1760-1764.

[29]毛贵牛,韩小雷,陈学伟. 高层建筑结构风振响应计算方法初步研究[C]. 天津大学.庆祝刘锡良教授八十华诞暨第八届全国现代结构工程学术研讨会论文集.天津大学:全国现代结构工程学术研讨会学术委员会,2008:1890-1894.

[30]陈学伟,韩小雷,季静,郑宜. 基于纤维模型的CFRP布加固混凝土柱有限元分析[C]. 中国土木工程学会FRP及工程应用专业委员会.第五届全国FRP学术交流会论文集.中国土木工程学会FRP及工程应用专业委员会:中国土木工程学会,2007:161-164.

【日志】2020年的端午节OPENSEES赠书活动！

2020年的端午节，DinoChen.com(我与林哲)发起了一个与《OpenSEES原理与实例》读者交流的线上活动——征集广大读者与第一版《OpenSEES原理与实例》的故事。从中抽取了一些精彩留言，如下表所示。为答谢各位小伙伴多年来的支持，特此送上第二版。《OpenSEES原理与实例》。新书已在某宝某东某多多发售，您们的支持是我们继续创作的动力！希望这本书能帮助有需要的学生、研究人员和工程师！谢谢大家！！ 这也是我们第一次在微信公众号平台举办的活动，积累了一些经验，希望下次能有更多的读者参与进来。 2020年4月20日，《结构弹塑性分析程序OPENSEES原理与实例》（第二版）终于出版了，同时京东，淘宝都有出售。不要再买一版了，都是盗版的。感谢伙伴林哲同学对本书的辛苦工作。感谢韩小雷老师，吕大刚教授，古泉教授及Tony为本书写序，感谢刘瑞霞编辑，感谢所有支持dinochen.com的人。样稿已经到了，新书厚了不少， 增加了10章的内容。

lc	学习OpenSees快两年了，入门就是看陈老师的教程，那时候一边记笔记，一边敲代码，渐渐的掌握OpenSees，很惊叹这门软件的开源性，在学习的过程中也完成了自己的课题“含屈曲约束支撑的自复位墙抗震性能研究”
吱唔猪	在学习opensees的过程中，最怕出现出现各种bug，建模不可怕，怕的是出现问题，还好有各位同学老师的指导。研究基于IDA的纤维混凝土桥墩的易损性研究
笛卡尔	使用opensees已两年多，从陈博士的第一版书开始入门再到相关opensees博客，在掌握相关操作的同时已完成两个小的课题，不得不说陈博士整理的笔记及相关讲解内容对科研工作有很大帮助。赞
魂馆	我是最近才接触opensees的，由于课题是竖向荷载对型钢混凝土剪力墙抗震性能的影响，本来是用abaqus的，但是由于不会写子程序导致滞回曲线捏缩的很不好，从而开始学习opensees，陈老师用etabs建模导出命令流对我这种初学者实在是太友好了，极大的加快了我对命令流的学习和认知，按照上面的实例一步步做，让我在学习opensees完成课题的路上少了很多的坎坷感谢陈老师
充满希望的六月	陈老师的这本书是很好的opensees 的入门教程，实例实用且全面，详细易懂。硕士阶段课题研究部分内容为老旧结构的易损性分析，曾按照此书学习搭建最简单的框架结构模型。在考虑结构的老旧时，不同龄期的混凝土强度相差不是特别大，opensees 相比其他计算软件来说对纤维截面的参数设置比较敏感，因此刚好适合做体现不同龄期混凝土的强度差异对结构抗震性能的影响，结果出来很不错。此次陈老师的书出了第二版，是一定要看一看的。
王帅	2017年硕士入学，2018年开始学习OpenSEES，乍看这程序挺吓人的，高端大气上档次但是一个按键都没有、自己建的模型也看不到。看着陈学伟博士的第一版教程学了下来，硕士期间的大小论文模型都是基于OS建的。我做的桥梁地震易损性方向，用OS计算时程是真的快，甩了MIDAS两条街，提取结果也方便，要什么有什么！而且学了OS和关注陈博后，慢慢喜欢上了编程，这东西啃着硬但用起来就是666。今年参加工作，以后打算继续深入学习OS，以后弹塑性分析少不了它！
故里	课题是钢结构半刚节点，通话教程实现了opensees入门的基础学习，尤其是局部坐标的转换关系，最重要的是为我打开了opensees系统学习的大门，没有教程只看使用手册中的代码真的很头疼，非常感谢陈博士
陈	想用openseess做钢管混凝土边框剪力墙结构和内置钢管混凝土异形柱结构的拟静力和整体结构的弹塑性分析,之前用opensees做钢管混凝土边框剪力墙模拟，主要参考是陈博的第一版书，那时候很期待第二版书，最近看了一些python的书藉，建议能不能出一版二次开发方面的书
飒鼎	在matlab对哈工大禁用后，国内研究者慢慢会更深的意识到国产软件的珍贵。Dinochen在Opensees的研究和应用方面都处于国内前列。国内很多Opensees的使用者都直接受惠于陈博士。陈博士本人也很亲民，自己工作很忙，还热心帮助回答建模问题，检查模型！第一版的书已经有很多科研工作者从中受益，第二版有很多升华，对推动国内研究者对opensees的了解功不可没！祝愿这本书越做越好。有个建议，是不是多做一些视频，比如慕课什么的，配合书，扫描二维码学习。
刘英棨Selby	初识第一版教程还是在2015年，那时怎会想到自己之后将会收到香港大学博士offer，更没有想到会在博一CIVL7008的课程上见到陈学伟老师本人。直到后来确定研究课题后才后悔，当时没能拿着书找您签字~~书中实例10、11的勾缝单元，以及实例20的隔振zero-length单元是非常好的学习资料，对我后续进行hybrid sliding-rocking预制拼装节段桥墩相关地震响应模拟提供了很好的入门指导。新版热卖在即，祝愿陈老师的工具书能造福更多的工程人~~（这学期也想去课堂上签字来着，但是被疫情耽误了，希望来年还能有机会hhh）
yz	学生阶段一直用ansys和abaqus。工作之后因为软件版权问题基本上都是用midas，现在做一条高烈度地震山区公路桥梁的弹塑性抗震验算，就找到了opensees，免费，计算速度快，加上可以像ansys一样，采用命令流的方式交互，可以兼顾工程需求和研究需要，再加上陈博的教程，完美！
@阿怪	自从去年导师推荐我们学习opensees后他就把陈博的教程推荐给了我，在书没到之前一脸懵逼不知从何处下手，有了教程后自己一步一步的跟着做，边做边记，越到后面学的越快，并且教程基本学一遍就记下了，通俗易懂，我的研究课题是框架柱约束下的再生混凝土墙体抗震性能研究，从 10月份开始学，学了差不多2个月，到12月份做完试验，然后疫情期间就把小论文写了，非常感谢陈博的教程，希望越做越好，师者传道授业解惑也。期待第二版教程。
CY	读研的时候用的opensees，课题是挡土墙在地震作用下的反应。导师是opensees的高手写了不少opensees的examples，当时只是在导师的程序里改来改去，想借陈博的书真正系统的学习opensees，BTW，京东一直缺货。

【ETE专题】基于YJK批量生成墙与柱PMM曲面

[小编：李明]

​

简说：基于ETE对YJK模型中的墙柱进行PMM曲线与曲面的自动化生成，程序采用DELPHI进行编写，自动生成曲线包括超限报告的WORD文档与EXCEL表格，还可以自动生成所有的XTRACT模型，用第三方的软件进行复核，截面验算的计算核心是DINOSEC，相关请到网站DINOCHEN.COM查看。
[下载讲义PDF]：点击下载此文件：ETE批量生成YJK墙柱PMM曲面的方法

以下介绍一下ETE通过读取YJK的截面数据与内力数据后自动生成PMM曲面的功能介绍

（1）众所周知，YJK可以生成全部数据包，包括前处理的结构信息与后处理的构件内力
通过ETE读取这个文件的位置，同时读取小震，中震与大震的数据。

（2） 基于DINOSEC的计算核心可以快速地计算矩形墙与不同形式的柱的PMM曲线与曲面
点击ETE直接可以生成曲面，通过内力组合生成YJK的构件内力点，以检查是否能包住。

（2） 剪力墙小墙肢也可以计算了，点击墙肢号，可以在平面上查号，然后生成PMM曲线或XY两个方向的曲线

最后，也可以生成XTRACT的数据，导入XTRACT进行第三方软件的验证，或直接生成EXCEL表格格式与WORD超限报告的相关章节的文档。

以上就是ETE【中国版】关于PMM曲面自动生成的方法（基于YJK的数据）

有兴趣请关注我们的微信公众号:DINOSTRU

以下是软件的操作视频：[Bilibili]B站视频

https://www.bilibili.com/video/BV1T54y1q7o1/

[Dino笔记]风荷载规范中的SFL图表问题 


副标题是：横风向广义力功率谱SFL参数的计算分析

在大部分的工程师看来风荷载的计算是一件容易的事情，就是PKPM或YJK点一个按钮的功夫，但是风荷载背后
的计算公式是相当复杂的，2012年的风荷载规范又把这个难度提高了，里面都是一些复杂的公式，如果说结构工程除了有限元最难可能就是风荷载了。这里面有一个计算参数叫SFL，横风向广义力功率谱，规范附录给出的计算公式与图表是不对应的，一直想不到原因，现在帖出来，大家来研究一下，到底是什么原因导致的。 有兴趣的朋友下面可以参与评论。图片有可能不清楚，可以直接下载PDF查看。

【点击下载】关于横风向广义力功率谱SFL参数的计算分析.pdf

【小编】辛展文同学 (华南理工大学)
 
分析背景：
在2012年荷载规范中，计算矩形截面横风向风荷载时，需要使用横风向广义力功率谱SFL，该值可以通过公式法计算得到，涉及建筑高度、宽度、进深和折算频率等参数，计算过程较为复杂。
同时，规范也给出H/√BD=6下的横风向广义力功率谱图表，供设计人员根据折算频率与深宽比查表确定。



2012荷载规范图

但是，在编程计算的过程中，发现根据广东规范给出的公式计算得到的横风向广义力功率谱SFL和图表查得数值不完全相符。
通过观察发现规范图表上的曲线不是很光滑，于是猜测这和绘图的精度有关。

广义力功率谱SFL计算公式：

Matlab编程分析根据规范给出的公式，使用Matlab编程绘制等值线图来验证猜想。
首先使用Matlab对一个500*500的矩阵进行运算，得到横风向广义力功率谱如下所示，可以看出它与规范中的图形比较不一样。


Matlab计算得到的图形上的曲线更圆滑，没有出现曲线断开的情况，部分曲线与规范中的变化趋势基本一致。

Matlab 500*500 矩阵数值模拟结果 

在此基础上降低精度，当进行到20*20的矩阵计算时，终于得到一个和规范相似的图形。得到图形如下：

Matlab 20*20 矩阵数值模拟结果在20*20的矩阵运算下，得到的等值线和2012荷载规范给出的基本符合，特别是在曲率变化和曲线断开的地方，都和规范图形相一致。但值得注意的是，在某些曲线上依然有细微的差别，比如A类地貌中最右边的这条曲线，规范中的是发散曲线，而Matlab绘制得到的曲线是收敛的。

这可能是由于不同软件之间的绘图处理方法不同，或由取值不均匀造成。
请大家帮忙，找找为什么会出现比较大的差异，特别是这里。 


在0.002等值线上的差异

【总结】 

（1）2012年荷载规范给出的横风向广义力功率谱SFL图表精度较低
（2）在使用查表法确定横风向广义力功率谱SFL时，在曲线不连续、斜率突变或折算频率fL1*大于0.2的位置，宜使用公式法进行校核，以得到更加准确的结果。

最后说一句：

规范也是一本书，尽信书不如无书，对于书我们一定要保有怀疑的精神，实践出真知，自勉！

==============================

这个是我的微信公众号[Dino结构笔记]：DINOSTRU
也可以通过扫二维码进行关注一下。

[每周挑战]造物_造一个使命必达盒子(MeeseekBox)


Rick and Morty瑞克与莫蒂是我喜欢的美剧动画片，现在的第四季播完了。上次制作了Rick的黄油机器人，这次我们来造一个Rick and Morty的粉丝一定会制作的最受欢迎的外星人Meeseeks Box 使命必达先生与他的盒子。
以下是使命必达先生盒子的造型与设计图纸（左下角）。 


【Rick and Morty中使命必达先生的图纸】
 第1步，3D模型设计在3DSMAX设计出使命必达盒子【Meeseek Box】的外形，如下图所示。
 它的结构分为以下几个部分：
（1）上盖，主要用于放大按钮与纸板人，
（2）盒子内部：全部电路、电池、喇叭与舵机都在里面。
（3）盒子外部：全部是装饰性的设计 

以下就是3DSMAX出来的效果与实际的对比 


第2步,准备电子材料 与以往的风格一样，这个装置现在会用到arduino mini，
（尺寸比较小）外加一个MP3播放控制小板子（USB接口，提前存入MP3文件即可）与舵机，"
会基本电路就可以安装了。
材料如下：

(1)1个9V的电池给ARDUINO供电
(2)1个ARDUINO MINI 板，
(3)1片孔洞板，用于制作ARDUINO外接电路
(4)电焊铁，电线（杜邦线）
(5)小型开关2个,一个用来打开ARDUINO板，一个控制电池供电
(6)1个舵机，控制盒子的开合动作
(7)1个控制MP3播放的芯片
(8)1个喇叭的元件
(9)1个10K的电阻与1个1K的电阻
(10)1个小按钮，放在盒子的顶部
(11)电线若干。其它就是用3D打印机打印。材料就是相对多一点，经过训练的朋友都能做这个小型装置。

第3步,分开模块与打印
在3DSMAX将不同颜色不同部分的组件分别排版，如下图所示。


 采用了不同颜色的线材，这次采用黑色线材打印，最后采用喷漆形成效果，喷漆颜色包括深绿色，青绿色与蓝色。 


3D打印机采用 XYZprinting Mini Maker，如下图所示。 第4步,电子电路的设计与制作 以下是电路图的设计，分为三部分，基本采用ARDUINO实例教程的基本组合


（1） 普通按钮开关实现开关状态
（2） 舵机电路（ARDUINO直控）
（3） MP3播放电路（用到MP3播放的小板子）

如下图所示，通过焊接电线与孔洞板，实现了整个电路的连接，完成图如下。 

第5步,模块的组装对全部构件进行粘合，采用强力胶对部份构件进行连接，如下图所示。
 
其中，内部有一个橡皮筋保证盒子的关闭状态。

 第6步,写入ARDUINO的编程代码对ARDUINO板写入代码，这个代码就是就是ARDUINO板实现读按钮状态，实现舵机转动与播放MP3的语音功能。其中MP3文件可存在MP3的小板子中，附件有原音频文件。纸片人采用打印机打印即可。代码，如下所示。

下载方法如下：
【点击下载你需要的arduino代码文件】点击下载此文件arduino.rar



第7步,最后的测试 以上是使命必达先生与它的盒子的最后完成图，打开小开关可以启动盒子，
以下是动画的全部过程，上传到B站。

[每周挑战]造物_造一个黄油机器人(Butter Robot)

[每周挑战]造物_造一个黄油机器人(Butter Robot)
Rick and Morty瑞克与莫蒂还是我最喜欢的美剧动画片，现在的第四季下半部份又更新了。
上次制作了Rick的传送枪，这次我们来做个大一点的活，就是黄油机器人了。
我觉得一个Rick and Morty的粉丝兼一名Maker 一定会制作这一个剧中最受欢迎的机器人。
以下是黄油机器人Butter Robot的造型与设计图纸。 

【Rick and Morty中机器人的图纸】

 第1步，3D模型设计

在3DSMAX设计出黄油机器人【Butter Robot】的外形，如下图所示。 它的结构分为以下几个部分：

（1）头部，主要用于隐藏大部分的电子电路，
（2）小车部分，有两个电机，4个1.5V与一个9v的电池都在里面
（3）舵机支架，中间有舵机，连接头部与小车。 

以下就是3DSMAX出来的效果与实际的对比 

第2步,准备电子材料 与以往的风格一样，这个装置现在会用到arduino mini，
（尺寸比较小）外加一个电机控制小板子DRV8833，会基本电路就可以安装了。
材料如下：

(1)1个9V的电池给ARDUINO供电
(2)1个ARDUINO MINI 板，
(3)几片孔洞板，用于制作LED灯与外接接收的外接电路
(4)电焊铁，电线（杜邦线）
(5)小型开关2个,一个用来打开ARDUINO板，一个控制电池供电
(6)1个红外接收器与1个红外控制器
(7)1个LED灯，用来显示红外的接收情况
(8)1个舵机，控制头部的活动
(9)2个小黄电机，控制机器人的活动
(10)电线若干。其它就是用3D打印机打印。材料就是相对多一点，
经过训练的朋友都能做这个小型装置。

第3步,分开模块与打印 在3DSMAX将不同颜色不同部分的组件分别排版，
如下图所示。

采用了不同颜色的线材，包括银色，灰色（由于银色比较难找，所以采用喷漆）

3D打印机采用 XYZprinting Mini Maker，如下图所示。

第4步,电子电路的设计与制作 以下是电路图的设计，分为三部分，

基本采用ARDUINO实例教程的基本组合
（1） 红外接收电路（增加一个LED显示状态）
（2） 舵机电路（ARDUINO直控）
（3） 小黄电机电路（用到DRV8833板子）
如下图所示，通过焊接电线与孔洞板，实现了整个电路的连接，完成图如下。
 
第5步,模块的组装
对构件进行粘合，采用强力胶对部份构件进行连接，如下图所示。
 其它部位不需要采用胶水连接，直接采用M3X10的螺线连接就可以了。
其中，头部的两根电线是装饰没有功能。

第6步,写入ARDUINO的编程代码
对ARDUINO板写入代码，这个代码就是就是ARDUINO板实现读取红外信号，
实现电机转动与舵机动作的代码。
下载方法如下：【点击下载你需要的arduino代码文件】
点击下载此文件arduino.txt

第7步,最后的测试 以上是Butter Robot的最后完成图，

打开小开关可以启动机器人了，通过控制器可以控制它，以下是动画的全部过程。 

【Dino结构笔记】伸臂桁架的几何优化方法
结构工程分析计算中的一题多解法

最近学习【非解构】的公众号文章【超高层伸臂桁架的那些事儿】，有所心得体会，
伸臂桁架作为重要的加强层方案在结构的刚度难以满足的时会被使用。这么好的东西，
我肯定用过了，如下图所示（左侧是伸臂桁架的上部，右侧是腰桁架）。 


主要的使用场景有如下：
（1） 普通的框筒结构的刚度满足不了，相差还比较大，靠加大截面会造成不经济的情况下采用伸臂桁架，对于腰桁架效果更好
（2） 巨柱的方案，由于要把内筒的弯矩传给外筒，需要采用伸臂桁架
（3） 吊柱体系，由于要制作大空间，或外框柱间断了，局部采用伸臂桁架连接吊柱。更合理的名称也许是悬臂桁架（英文Hanger Truss）。"
（4） 无梁楼盖的框筒高层体系，因为没有内框梁与外框梁，所以需要伸臂桁架把弯矩带到外柱上，这种情况多发生于国外的高层。     


这篇笔记讨论伸臂桁架的几何优化，也就是伸臂桁架的局部布置的优化，通过以下的过程，
你会发现采用连续体的拓扑优化与基于虚功原理的桁架简化优化得出来的结果是一样了。
以下是这个方法的详细描述。

我们从伸臂桁架X支撑开始吧。

 方法一，连续体的拓扑优化 

算例：设一个连续体9mx9m（正方形比例），上端自由作用作用两个集中力P，底部两点固定支座，对该连续体进行拓扑优化，本文采用ESO方法，笔者基于OPENSEES的二次开发小程序，下载方式如下。连续体的拓扑优化的原理是简单的，经过应力分析后，把应力值最小的一部分（通过拒绝率RR进行控制）删除掉，然后进行下一步分析，到不要再删除后提高拒绝率，如此循环下去，就可以得到一个优化后的拓扑形状，该形状的特点就是在一定外力作用下，构件的每个区块的应力是饱满的，也就是把材料用尽达到优化的要求。



以下是优化的整个过程的小动画。

关于基于OPENSEES的拓扑优化的原理与程序，可以转到这个帖子去看。
[OPENSEES]实例35_基于OPENSEES拓扑优化程序开发
http://dinochen.com/article.asp?id=315  

方法二，基于虚功原理的优化 

（1） 由于伸臂桁架的轴力较大，以下实例简化了桁架构件（其实桁架的上弦下弦是压弯构件，所以应该采用梁柱单元来做的），建立ETABSS模型如下图所示。中间节点的四个杆件要划分单元，生成S2K文件后备用。计算草图如图所示。

（2） 采用ETO程序把模型导入S2K文件，如下图所示，生成OPENSEES的分析代码。

（3） 打开之前开发的OPENSEES的桁架找形程序，对中间点的位置进行找形分析，找形分析程序通过迭代后（也就是不断改变中间点的X坐标，重新提交计算）得到最优解如图表所示。

 基于OPENSEES的桁架找形程序  ETABS与ETO程序的伸臂桁架的局部建模

 基于能量的原理算出最优的位置图表所得解X=2250mm,也就是X=0.25B。
最终，两个方法是【殊途同归】的得到相同的结论，正方形比例的伸臂桁架的最优解是X=0.25B。
我再采用图解法，得到这个桁架的内力图解的图如下所示，非常用趣。 

那么换一个问题，换成V形支撑形式的伸臂桁架，虚功的最优解在哪里呢？
以下图为算例的求解过程。V形支撑的最优点在中间 
这就是局部伸臂桁架的几何优化的计算方法的简单介绍，
以上就是这篇笔记的全部内容。

==============================

这个是我的微信公众号[Dino结构笔记]：DINOSTRU
也可以通过扫二维码进行关注一下。

【Dino结构笔记】结构图解法的学习与编程
 基于桁架结构图解法的找形优化方法的开始探索

这期继续学习【结构小站】那一篇关于【题目是：建筑的“形”与“力”—形由力生，以力塑形，以形驭力】的文章里面相关的内容。这文章里面介绍了图解静力法，又可以称为桁架结构图解法等。这个方法的详细原理可以参考如下图的这本书Graphical Analysis（出版于1921年）。还有一点，这里面的图解法的插图太漂亮了，想起以前的桁架钢构桥都是用这个方法来设计的，如下图所示。

 悉尼海港大桥（Sydney Harbour Bridge） 

经典力学教科书 图解分析法

说起来也是有趣，以前的结构力学书是有这部分内容的，由于现在电算程序为主要计算手段，教材也没有介绍这部分桁架图解法，可能是因为有了电算这部分的价值不高，但是最近几年的超高层的设计发展，桁架图解法可用于大跨或高层的结构优化设计当中，这不得不重新去学习一下图解法了。


基于图解法的结构优化的工程应用   

【原理】图解静力法的基本原理
1. 作为结构，所有外力与支座反力可以构成一个环的向量和，也就是所有力向量之和是为零的。
2. 每个结点的力需要平衡，也就是一个节点多个作用力（外力支座反力与相近构件的轴力）对这个节点求和也是为零的。
3. 桁架由许多三角形组成，每个三角形三个构件的内力和不一定平衡，不平衡的力为作用于相邻近的三角形的桁架上，但是全部桁架三角形的内力传递完成后最终的合力也是为零。

综合上述的原理，就通过简单的几何原理就发展出结构图解法。


【操作】图解静力法的基本操作。对于操作而言，再多的文字与描述也是很难记录细节的，所以我通过ipad录制了以下实例的图解法操作的整个过程，有兴趣的朋友可以查看B站的视频，连接如下。 
1. 求解结构的支座反力
2. 把外力与支座反力绘制出力线部分（一般在最右端），标出由外力划分出来的线段名（大写字母）
3. 把桁架划分成许多三角形并标出号码，在桁架是面域，在图解法的图中是一个节点。
4. 根据三角形内的桁架角度关系与已标出的力线段的关系，通过直接相关找出节点位置
5. 相邻三角形可以通过相邻边的角度查找节点位置。

通过以上步骤就可以绘制出图解法的“图”，然后就可以通过测量力线的长度得到最后的构件轴力了。通过以下的草图向大家讲述求解的过程，详细可以去看看动画展示。


【Bilibili】结构图解法的详细绘图方法的视频
https://www.bilibili.com/video/BV1Lf4y1U7FS


 
【编程】图解静力法的程序可视化编程这一次学习图解法是为了更好地学习结构优化的老方法与新方法，所以需要大量的绘制图解法的图，那么：
手画是不可能的，这辈子都不可能手画的，手画又容易错，肯定要编程啦，编程的感觉好多了，画图的速度又快，准确率又高，超喜欢编程的。（老梗如果懂的请在下面回复） 


 (1) 编程的实现：编程的基本原理就是参考三角函数以及图解法的基本过程， 程序的界面如下图所示。

 (2) 验算算例。之前学习了SOM的马克写的《高层建筑设计-以结构为建筑》这本书，这本书提到了很多优化的方法，建议大家学习，同在SOM的Baker写了这一篇文章《Structural optimization using graphic statics》就是介绍SOM如何把图解法应用到结构找形与优化上面来。那么我就展示一个这个软件绘制图解法的算例吧。

算例1：  

算例2：  

 
算例3：  

 算例4：  

 
(3) 后续开发。

现在已经教会电脑程序自已通过计算绘制这个图解法的图了，接下来就要学习如何通过它进行结构优化，我开发了一个节点偏移的功能，随时改变节点就可以得到图解法的图的变化，这个就是个入门级别的优化手段了。这个就是这一篇DINO结构笔记的全部内容。

==============================

这个是我的微信公众号[Dino结构笔记]：DINOSTRU

【Dino结构笔记】桁架结构找形的优化方法
 基于OPENSEES的桁架找形程序的编写

【计算文件内容】相关的ETABS/SAP2000/OPENSEES模型的下载

最近结构圈多了好多优秀的公众号，推荐一下这几个，如【istructure】，【结构小站】这些优秀的公众号展示了工程师对结构工程的热爱与执着，配合最近的bilibili的宣传片【后浪】，让我这个老浪要觉得更加要好好学习，最近看到结构小站的一篇关于形与力的总结，题目是：建筑的“形”与“力”—形由力生，以力塑形，以形驭力，大家可以上网搜一下。 

                                   什么时侯有空也可以试下逆吊实验法

这篇文章真的写得不错，也启发了我，特别是介绍到前面大师们的找形作品，包括“逆吊实验法”去找形等，这个应用于西班牙圣家族大教堂的方法，我到西班牙第一次看到真实的模型，真的是震惊了。这个公众号的文章让我非常想做一下找形分析的尝试，那么我们开始吧。这次的前言写的有点多，一般我说做就做的。内容过于简单，请家长陪同观看。 

                                       去西班牙一定要拍一下这张照片

简单的原理：我通过一个两端铰接的梁来介绍找形的原理。首先以一个两端铰支的梁为例，大家都知道要提高刚度需要起拱，其实起拱的方法也可以通过找形来做，通过施加重力荷载，使梁发生变形，如图所示。记录每个节点的变形位移，最后把变形量乘于以一个负倍数得到反拱量，然后用于原结构的几何设计调整。原理是很简单的，这个一个很好的工程方法，确定这个倍数是一个知识点。如果梁的找形就会变成拱，板的找形就会变成拱壳，桁架的找形就变成拱壳（曲面）桁架等等。  

静力法找形状。本文以一个桁架的算例，如图所示。采用ETABS或SAP2000进行弹性建模可以得到以下的模型。施加重力后可以得到变形图如下图所示。采用SAP2000可以把变形值导出来，这个不是反拱值，这个是每个节点位置反拱的比例，至于总反拱要用优化算法算出来。桁架没有反拱的变形大约是16mm。 

打开ETO，ETABS TO OPENSEES，把ETABS建出来的模型导出E2K，导入到ETO程序中去，通过这个程序了解到21号节点是跨中节点，我们用这个节点作为控制节点。点击生成OPENSEES的分析代码，保存为co.tcl。

 在建模的过程中保存结构的节点坐标，节点在重力荷载下的变形，这个数据作用往后每步生成OPENSEES模型的准备。作用文本格式导入自编程序用于计算。 

调试OPENSEES使程序可以运行一个算例即可，证明OPENSEES可以进行一个荷载步的静力分析。那么打开编写的小程序，小程序太简单了，主要就是写一个小程序不断改变节点计算新的OPENSEES模型。每步生成带有新节点信息的OPENSEES模型提交运算，最后输入所有桁架构件的内力，用于计算指标。还有21号节点的位移值。 以下是程序不断调用OPENSEES程序的效果。

为了得到较为优化的结果，对于桁架结构，反拱后的结构形状必须使结构满足以下公式（我就简称力流公式吧，也可以理解成能量），代表最后满应力优化后，这个值小代表结构总的用钢量少，那么这个形状就是相对合理的形状了。如果TCL语言好的同学，这个公式的计算可以采用OPENSEES的TCL命令流计算得到。



目前这个算例的反拱值是从1m到5m，通过公式找到最小值，局部细化就可以找到最小值是2.8m的反拱量，最后我们把这个反拱量导入SAP2000就可以显示出这个相对优化的形状了。是不是很简单呢。大家可以试一试。

由于OPENSEES是一个没有界面的小程序，所以计算速度还是很快的。这个反拱后的桁架是不是好熟悉，对比一下，有点像。最后的优化结果桁架的变形从16mm优化到4.277，用钢量比例从4.17优化到3.82。 

最后有没有觉得桁架的找形后的形状有点眼熟，放在以前的古建筑中就非常之和谐了。




======================================================== 

这个是我的微信公众号[Dino结构笔记]：DINOSTRU

【Dino结构笔记】

真的好久没有在这个网站上写下Dino结构笔记了，今天总结一下以前关于弹性时程分析思考。
今天主要讲我对弹性时程分析中出现的楼层剪力R形突变的看法与研究。

【PPT下载地址】：弹性时程分析顶部剪力变大的原因研究.pdf 

【计算文件内容】： attachment.rar 包括ETABS模型，地震波与反应谱，最后表格数据 

（1）一般我们进行弹性时程分析，需要进行选波，选波后进行时程分析后一定会对比楼层剪力与反应谱的区别
往往在这个过程中，我们会发现顶部的楼层剪力比反应谱得出来的剪力还要大，似乎采用CQC或SRSS这种包络法
也不能包得住，那么一般文献会说这是由于高阶振型引起的顶部放大，其实，一直我有所怀疑，原因如下。

（2）为了更好的证明我的假定，我采用数值分析的方法比较直观，如下所示，我找了一个普通的超高层结构，
对于200米的超高层，一般会开始出现这种剪力的R形突变。

（3）以下是结构的空间三阶振型，要注意，这个是前几阶空间振型，不是X方向的前几阶振型，X方向的前几阶
振型分别是T1,T4,T7，如下图所示。

（4）以下是X方向前几阶周期值在反应谱的位置，表上是每个振型对X方向变形的参与系数，可见，第一阶振型为
70%，到第14振型，参与系数为1.122%,这个数值已经表明，高阶振型对宏观的位移及剪力的影响越来越小。
理论上不存在高阶振型影响楼层剪力产生R形突变。



（5）图为X方向每个主要振型的周期所在的位置，那么接下来为了证明R形突变是前几阶振型引起的，我们采用
扰动的方式造人工波进行时程分析。


（6）我们通过对反应谱进行扰动，如图所示，然后采用造人工波软件进行人工波的制作。GM0是不加扰动的，
GM1为对第一周期位置进行扰动，如此类推。



（7）图为扰动后的造波的反应谱与目标谱的对比。
（8）经过弹性时程分析，得到X方向的楼层剪力图如下图所示，可以发现只有第一阶扰动与第二阶扰动才对楼层
剪力有宏观的影响。

（a）从第三阶开始，扰动基本不影响剪力分布，基本可以排除高阶振型对剪力R形突变的影响
（b）第一二阶的扰动，对总剪力影响较大，第一阶扰动并不会引起R形突变
（c）第二振型是关键 ，第二振型开始出现R形突变
（d）为了进一步研究，我们增加对第二阶的扰动



（9）为了进一步证明第二周期对楼层剪力的R形突变的影响，我们增加对第二周期的扰动，操作如下：

（a）降低第一阶附近的地震影响系数
（b）增加对第二阶周期附件的扰动

在天然波的选择中，这的确是一个普遍现象，在周期大于4S的情况下，大部分地震波的影响系数小于规范值
为满足总剪力的需要，很多满足条件的地震波会在1~3s的区间有较大的影响系数值 。

目前的结论：
（1）高阶振型并不会影响剪力的分布，特别不会引起R形突变
（2）通过修改反应谱造人工波，证明剪力的R形突变是第二振型与反应谱的形状引起的。
（3）天然波产生的基底剪力不足，经常会第一不够第二来凑，因此产生R形突变
（4）本文只用了一个算例，提供数值如附件所示，大家可以测试与测算
（5）本文基于数值分析得到小结，如果有大神可以采用解析公式法进行证明就更好了

本人才疏学浅，以上只是学习小结，以供参考，如有不对，请指正。

如果你们对这种文章感兴趣，请留言并关注我的微信公众号 DINOSTRU

重要通知
[2018-12-21]重新整理OPENSEES的汇总帖，增加视频教程、相关资料等栏目如下。
[2014-09-01]结构弹塑性分析程序OPENSEES原理与实例出版
[2011-12-30]开始写关于OPENSEES的实例教程的书

（1）程序的基础是，openSEES 程序下载， 下载地址为：
=======================================
点击下载此文件:OpenSEES v2.2.2
点击下载此文件:ActiveTcl8.5.11
点击下载此文件:OpenSEES v2.5.0
点击下载此文件:ActiveTcl8.5.18

程序的安装及入门使用，可以参考以下：
http://dinochen.com/article.asp?id=92


=======================================
（2）OpenSEES的前后处理程序(ETO)下载地址：
ETO软件下载方法：
微信搜索微信公众号: DINOSTRU  
回复 ETO 得到软件下载地址
扫以下二维码：微信公众号: DINOSTRU

【2014-03-25】DINOETO程序下载文件
【2016-08-01】DINOETO程序下载文件
【2019-07-01】DINOETO程序下载文件
【2019-07-17】DINOETO程序下载文件
【2019-08-06】DINOETO程序下载文件

微信公众号：搜索 DINOSTRU 然后回复 ETO 
=======================================
(3) OPENSEES的实例教程书上配有的光盘



附带光盘：下载地址：http://pan.baidu.com/s/1sjxGCZN

=======================================
(4) 技术讨论的相关OPENSEES的QQ群：

(5) 微信公众号：搜索 DINOSTRU
=======================================


 软件下载地址：http://dinochen.com/article.asp?id=234
 版本更新说明：http://www.dinochen.com/article.asp?id=209

程序界面如下所示。简单而言，程序是基于openSEES的分析核心，采用ETABS进行
结构建模，导入到ETO进行其它信息的补充后直接生成  openSEES 的代码文件，可
以按您的意思进行修改成为你的tcl命令。学习ETO之前，可能需要看看ETABS的书。
联系邮箱：dinochen1983@qq.com
详细的操作过程，请看各章实例：

(6) OPENSEES实例教程的每个章节内容

 

[实例00]_前后处理程序ETO的入门教程 [实例01]_桁架桥结构静力分析 [实例02]_多层框架结构静力分析 [实例03]_简支梁弹塑性分析 [实例04]_框架结构推覆分析 [实例05]_框架结构模态分析 [实例06]_框架结构弹性时程分析 [实例07]_框架结构弹塑性时程分析 [实例08]_钢结构低周往复分析 [实例09]_钢结构网壳的屈曲分析 [实例10]_单压连接单元的应用 [实例11]_缝单元的应用分析 [实例12]_杆件铰接的处理方法 [实例13]_弹性壳单元的应用分析 [实例14]_网架弹塑性分析 [实例15]_预应力梁弹塑性分析 [实例16]_桥梁多点激励动力分析 [实例17]_剪力墙构件低周往复分析 [实例18]_框架剪力墙结构推覆分析

[实例19]_带粘滞阻尼器的框架动力分析 [实例20]_带隔震的框架动力分析 [实例21]  时程曲线转化为反应谱方法 [实例22]  截面PM曲线分析方法 [实例23]  实体单元的建模及应用 [实例24]  三维钢结构节点应力分析 [实例25]  桥梁结构的影响线分析方法 [实例26]  组合梁的弹塑性分析 [实例27]  型钢混凝土柱的静力弹塑性分析 [实例28]  带防屈曲支撑的钢结构低周往复分析 [实例29]  框架结构拟倒塌试验分析 [实例30]  铰纤维单元的弹塑性分析 [实例31]  单元生死在分析当中的应用 [实例32]  基于分层壳的剪力墙弹塑性分析 [实例33]  侧向多自由度简化模型的建模 [实例34]  基于OPENSEES的桥梁游戏开发 [实例35]  基于OPENSEES拓扑优化程序开发

(6) 基于OPENSEES的小软件、小游戏开发

[程序1]基于OPENSEES的游戏BLOCK STRU GAME
[程序2]基于OPENSEES的截面分析程序
[程序3]基于OPENSEES的桥梁设计小游戏
[程序4]基于OPENSEES的结构拓扑优化小程序

(7) OPENSEES视频教程系列


[视频教程01]OPENSEES安装与基本操作
[视频教程03]单柱构件的静力弹塑性分析
[视频教程04]框架结构的静力弹塑性分析

[出书] 结构弹塑性分析程序OPENSEES
原理与实例[第二版]终于出版

书的购买方法：
【淘室】中国建筑工业出版社的淘宝店
【京东】中国建筑工业出版社的京东店

请关注我的微信公众号 dinostru, 往后的Opensees的视频教程发布

2020年4月20日，《结构弹塑性分析程序OPENSEES原理与实例》
（第二版）终于出版了，同时京东，淘宝都有出售。不要再买一版了，
都是盗版的。感谢伙伴林哲同学对本书的辛苦工作。感谢韩小雷老师，
吕大刚教授，古泉教授及Tony为本书写序，感谢刘瑞霞编辑，感谢所
有支持dinochen.com的人。样稿已经到了，新书厚了不少，
增加了10章的内容。

现在的书没有配光盘的了，大家可以扫描书上的二维码得到与以前光
盘一样的内容。


★以下是第二版的书新增加的实例

[实例26]，组合梁的弹塑性分析
[实例27]，型钢混凝土柱的静力弹塑性分析
[实例28]，带防屈曲支撑的钢结构低周往复分析
[实例29]，框架结构拟倒塌试验分析
[实例30]，塑性铰纤维单元的弹塑性分析
[实例31]，单元生死在分析当中的应用
[实例32]，基于分层壳的剪力墙弹塑性分析
[实例33]，侧向多自由度简化模型的建模
[实例34]，基于OPENSEES的桥梁游戏开发
[实例35]，基于OPENSEES拓扑优化程序开发

[每周挑战]造物_造一盏心形的七彩灯

三八女神节又来了，创客们又开始制作应节作品了。今天为了投稿instructable的小比赛。
特意制作了一个心形的彩灯，也当做对Arduino项目制作的练手。


 第1步，3D模型设计

在3DSMAX设计出心形灯【arduino heart shape Light】的外形，
如下图所示。它的结构分为以下几个部分：

（1）透明的灯壳。
（2）灯座，主要用于隐藏大部分的电子电路
（3）灯座的顶盖.

以下就是3DSMAX出来的效果图

第2步,准备电子材料

与以往的风格一样，这个装置现在会用到arduino mini，（尺寸比较小）
会基本电路就可以安装了。材料如下：

(1)1个9V的电池给ARDUINO供电
(2)1个ARDUINO MINI 板，
(3)几片孔洞板，用于制作LED灯的外接电路
(4)电焊铁，电线（杜邦线）
(5)小型开关1个,一个用来打开ARDUINO板
(6)小型按钮1个，用来触发LED灯
(8)1个三色RGB LED灯
(9)一个用于开关的电阻，100K欧电阻
(10)电线若干。其它就是用3D打印机打印。
材料就是相对多一点，经过训练的朋友都能做这个小型装置。
注意，本项目没有电机，只有LED灯，采用9V电池就可以了。
       如果这个灯是固定的，也可以用USB直接供电。



第3步,分开模块与打印

在3DSMAX将不同颜色不同部分的组件分别排版，如下图所示。
采用了不同颜色的线材，包括白色与一个白色透明颜色。
3D打印机采用 XYZprinting Mini Maker，如下图所示。
用于本项目的STL文件在这里下载：
【点击下载STL】: arduino_heart_shape_light.stl



第4步,电子电路的设计与制作

以下是电路图的设计，LED，按钮与Arduino的连接。
如下图所示，基本上用了所有的Arduino的输出与输入口了，没有进行简化。

如下图所示，通过焊接电线与孔洞板，实现了整个电路的连接，完成图如下。




第5步,模块的组装

对构件进行粘合，采用强力胶对连接环与顶盖进行连接，如下图所示。其它部位
不需要采用胶水连接，直接拼掉就可以了。

如下图所示，把主要电路在灯座中。

如下图所示是未拼接之前的完成图。拼接后的完成图如下图所示。



第6步,写入ARDUINO的编程代码

对ARDUINO板写入代码，这个代码就是就是ARDUINO板实现读取按钮的值
实现对LED进行控制。

代码，如下图所示。下载方法如下：

点击下载此文件arduino.txt

第7步,最后的测试





以上是arduino heart shape light的最后完成图，
打开左侧的小开关可以启动ARDUINO板。
按下小按钮，灯可以发出不同颜色的光。
注意，采用同样的结构，可以制作出不同效果的灯，如闪灯，火灯及渐变灯等。



 

[每周挑战]造物_造一把瑞克的时空传送枪

Rick and Morty瑞克与莫蒂是最近我最喜欢的美剧动画片，现在
已经是第四季了，已经成功路转粉了。我觉得一个Rick and Morty
的粉丝兼一名Maker 一定会制作一把属于自已的Arduino枪。
再加上，这是一个很好的用于Arduino的入门教程，包括四位数码管
滑动变阻器，开关与LED灯都有了。

以下是传送枪在动画中的图形。

【Rick and Morty与时空传送枪的图片】

 第1步，3D模型设计

在3DSMAX设计出时空传送枪【Portal Gun】的外形，如下图所示。
它的结构分为以下几个部分：

（1）枪身，主要用于隐藏大部分的电子电路，
（2）枪头，有三个LED灯
（3）顶部，上部的大灯管，四位数码管也在上面
（4）手柄，带有一个小按钮
（5）旋转按钮，颜色是黑色的。

以下就是3DSMAX出来的效果图

第2步,准备电子材料

与以往的风格一样，这个装置现在会用到arduino mini，（尺寸比较小）
会基本电路就可以安装了。材料如下：

(1)1个9V的电池给ARDUINO供电
(2)1个ARDUINO MINI 板，
(3)几片孔洞板，用于制作LED灯与数码管的外接电路
(4)电焊铁，电线（杜邦线）
(5)小型开关1个,一个用来打开ARDUINO板
(6)小型按钮1个，用来触发LED灯
(7)一个扭动变阻器，用于控制数码管的数字
(8)7个三角LED灯，（去掉左右线变成纯绿色的LED灯）
(9)一个用于开关的电阻，100K欧电阻
(10)电线若干。其它就是用3D打印机打印。
材料就是相对多一点，经过训练的朋友都能做这个小型装置。
注意，本项目没有电机，只有LED灯，采用9V电池就可以了。

第3步,分开模块与打印

在3DSMAX将不同颜色不同部分的组件分别排版，如下图所示。
采用了不同颜色的线材，包括白色，黑色与一个白色透明颜色。
3D打印机采用 XYZprinting Mini Maker，如下图所示。

第4步,电子电路的设计与制作

以下是电路图的设计，分为两部分，一部分是四位数码管与Arduino的连接。
另一部分介绍变阻器，LED，按钮与Arduino的连接。
如下图所示，基本上用了所有的Arduino的输出与输入口了，没有进行简化。

如下图所示，通过焊接电线与孔洞板，实现了整个电路的连接，完成图如下。



第5步,模块的组装

对构件进行粘合，采用强力胶对手柄与枪身部件进行连接，如下图所示。其它部位
不需要采用胶水连接，直接拼掉就可以了。

如下图所示，把主要电路身入枪当中。

如下图所示是未拼接之前的完成图。拼接后的完成图如下图所示。



第6步,写入ARDUINO的编程代码

对ARDUINO板写入代码，这个代码就是就是ARDUINO板实现读取电阻器的电阻值
实现对数码管显示的控制，通过读取按钮状态，确定是否打开或闪一下LED灯组。

代码，如下图所示。下载方法如下：

点击下载此文件arduino.txt

第7步,最后的测试



以上是Portal Gun的最后完成图，打开左侧的小开关可以启动ARDUINO板。

按下小按钮，枪发出绿光，数码管启动，可以调节变阻器，实现数码管的数字
的显示，从C-000宇宙可以穿越到 C-255宇宙了。当然包手了C-137了，如果
你能看懂上面的梗的话。

测试装置，打开开关，它就实现了上述的功能了，如下图所示。

 

【DINO爱编程】编写一个FlappyBird的游戏 

【编程工具】采用delphi进行编程

【资源网站】https://opengameart.org/

【游戏介绍】https://en.wikipedia.org/wiki/Flappy_Bird

游戏的设计理念
1. FlappyBird 是一只小鸟，要通过飞行穿过一系列的水管
2.点击屏幕或点空格键，小鸟向小飞一下
3.穿过一个水管得一分。如果撞到水管就会GameOver



【技术要点】

(1) 声音的播放方式。

procedure PlaySoundFile(FileName: string);
begin
 if FileExists(FileName)
 then PlaySound(pchar(FileName), 0, SND_ASYNC or SND_FILENAME);
 { Flags are:
 SND_SYNC =0 = Start playing, and wait for the sound to finish
 SND_ASYNC =1 = Start playing, and don't wait to return
 SND_LOOP =8 = Keep looping the sound until another sound is played }
end;

2. 对象之间的碰撞的分析与判断,如下图所示

if (((bird.x+20>pile[i].x)and(bird.x+20<pile[i].x+pile[i].w))
or ((bird.x-20>pile[i].x)and(bird.x-20<pile[i].x+pile[i].w)))
and (((bird.y-15>0) and (bird.y-15<pile[i].up))
or ((bird.y+15>pile[i].down-15) and (bird.y+15<image1.height)) ) then
begin
   collision:=true;
end;

3.按空格键对游戏的响应。

procedure TForm1.FormKeyPress(Sender: TObject; var Key: Char);
begin
if Key=#32 THEN
 begin
 end;
end;

 【FlappyBird的源代码】点击下载此文件：FlappyBirdexe.zip
【FlappyBird游戏的EXE文件】点击下载此文件：FlappyBird.zip


delphi 常用按键对照表

【DINO爱编程】编写一个扫雷 Mine Sweeper 小游戏 

【编程工具】采用delphi进行编程

游戏原理：(应该很多人会玩的，不用太多的介绍)
（1）点击鼠标左键扫描，如果碰到雷，全部雷爆炸，GAME OVER
（2）如果没有碰到雷，显示数字，数字代表9个格子内的炸弹数量
（3）可以标旗子，代表雷的位置，把全部雷标出就可以完成任务

 

【技术要点】

（1）BRCC32的方法，把文件放入程序的资源区
具体的方法如下：

创建一个MYRC.RC的文本文件，写下一段代码后，运行BRCC32.EXE就可以了

BMP1 Bitmap bmp1.bmp
BMP2 Bitmap bmp2.bmp
BMP3 Bitmap bmp3.bmp
WAV1 Wave k1.wav
WAV2 Wave k2.wav
WAV3 Wave k3.wav

（2）通过读取资源文件得到图片BITMAP与声音音效的文件

具体的方法如下（播放音效与LOAD入图片）：

BMP1.LoadFromResourceName(HInstance,'BMP1');
PlaySound('WAV1', 0, SND_RESOURCE);

（3）在开发游戏的过程中，应用到FloodFill的功能（图形学）
代码如下，也就是区域填充的功能，用到一点点的递归

PROCEDURE TFORM1.FILL(I,J:INTEGER);
BEGIN
IF (I<>0) AND (I<>SIZE+1) AND (J<>0) AND (J<>SIZE+1) THEN
BEGIN
IF (GRID[I,J].BOMB=0) AND (GRID[I,J].COVER=0)
AND (GRID[I,J].NUMBER<2) THEN
BEGIN
GRID[I,J].COVER:=1;
FILL(I,J+1);
FILL(I,J-1);
FILL(I+1,J);
FILL(I-1,J);
END;
IF (GRID[I,J].NUMBER>=2) THEN
BEGIN
GRID[I,J].COVER:=1;
END;
END;
END;

【小游戏的源代码】点击下载此文件：source code.rar
【小游戏的EXE文件】点击下载此文件：mine sweeper.rar