![[程序]DXF to Xtract 任意截面生成器](attachments/month_2010/kimage01.gif)
dinochen1983
48
0
Rhino+GH平台在结构设计中的应用
小编：Yifei.
结构的边界是形式。
新时代里，优秀的结构设计已经难以脱离建筑而评判。结构工程师的追求从超高层与大跨度，逐渐转向了形式的自由，并开始对除力学性能之外的其他属性关心了起来。除此之外，进入信息社会以后，人类的生产生活组织形态在看不见的地方飞速更迭。人们对空间的需求更加多样且灵活，但是更进一步，不仅住宅是居住的机器，每一个建筑甚至城市也都变成了更为复杂的机器。电讯派的Walking City从科幻逐渐降格为比喻。
Walking City
设计对象的极大复杂化使得建筑行业重新孕育了自己的系统工程。系统工程讲求加大设计阶段的投入，以尽量全面的设计考虑，避免在加工制造阶段的临时补救甚至返工。越复杂的系统，越难以通过现场的临时设计变更实现补救。因为建筑方案开始考虑更多的因素，在早期就会有大量的方案推敲和论证，而在系统工程中，需要借助基于仿真的设计方法（simulation based design），将难以实现的方案排除出去，以避免无用投入造成的浪费。因此复杂项目在方案阶段就需要引入可以及时配合的粗略分析结果来判定性能是否（难以）符合指标。
Rhino + Grasshopper因其既适合敏捷化快速配合，又有极大潜力通过开发实现深度配合，逐渐在设计阶段流行起来。设计团队可以在初期便实现对建筑的外观，热工，通风，视野，噪音，照度，能耗，力学，人流，微气候等等的分析，对方案表现进行综合性评估。同时对构造，施工过程的可行性进行推敲。其模型不仅可配合3D打印技术对原型进行论证，也可以生成可直接交付工厂加工的数据。
Rhino建模指导构件加工
无人机航拍生成场地模型辅助方案设计
场地开挖坡度分析
对于结构专业来说，置身于基于仿真的设计环境中，工程师需要能够快速对方案进行结构模型抽取分析，或对结构形态进行参数变量的批量分析，以实现敏捷配合。下面介绍一些简单工具，能够帮助工程师实现联动分析。
以下推荐几个GH插件，用于结构工程分析与设计：
(1)
Karamba3D应该是目前应用最多的结构分析插件。可以处理杆系和平面单元的有限元分析。分析方法包括线性分析，二阶分析，非线性分析，大变形分析，模态分析，屈曲分析，振动分析，截面优化，BESO拓扑优化等，分析结果的显示方式非常丰富。Karamba3D官网有中文版的使用手册可以免费下载，还有大量的分析案例和教程可以下载学习。
Karamba3D对一壳单元桥分析的教程案例
(2)
Kiwi!3D是一个应用IGA（isoGeometric Analysis）方法进行结构分析的差价。IGA方法直接以NURBS函数建立有限元模型，对于曲线构件无需进行折线重建，是有限元方法比较新的一个方向。该插件目前免费，除线性分析外，还支持非线性分析，找形，考虑施工工序的分析等。随插件会有教学案例下载，简单易学。
Kiwi!3D的一膜结构分析案例
(3)
Alpaca4D是一款调用OpenSees实现结构分析的插件，对OpenSees分析的输入和输出实现可视化的模型和分析结果。作为一款非常新的插件，目前后处理功能还比较简略。随插件下载压缩包同样内附有教学案例文件。
Alpaca4D一结构分析案例
项目进行到一定阶段，结构可通过OAPI生成Sap2000或Etabs模型进行进一步深化。Sap2000和Etabs都算是更新比较频繁的•软件，其OAPI亦随版本常常更新。所以可以想象建立模型转化的这一步工作也需要常常更新自测维护。目前虽然有一些插件可以实现Rhino模型到Sap2000或Etabs的转换，但是这里还是建议使用者学会对OPAI的使用，自己编写转换模块，以确保生成模型的正确，避免设计事故。
结构的生形逻辑多种多样，结构的布置方式要结合方案和场地，结构的设计理念需要工程师自己探索。但是俗话说“什么都不懂，就用遗传算法”，下面介绍两个遗传算法的插件。
(4)
Octopus是很多设计师在grasshopper平台应用遗传算法的启蒙，最早版本发布于2012年。随安装包附带使用说明手册以及使用培训案例。
Octopus计算模块对算法处理结果的可视化
(5)
Wallacei是一款比较新的基于遗传算法，并可进行多目标优化的插件。开发者主要为AAEmTech的毕业生。后处理及可视化更为丰富。随插件下载压缩包附带教学案例，同时官网也有视频教程。而且插件目前仍在不断更新。
Wallacei优化案例
ETH的DFAB House做了一系列研究性探索，前段时间被各媒体纷纷转载，火遍网络。他们的一个个子项目不仅展示了如何实现高度定制化的设计，同时也成为了前沿的数字建造技术的模板。在复杂形体，酷炫机械臂以及机器视觉的背后，是我们熟悉的Rhino + GH工作环境。
Smart Slab
Spatial Timber Assembly
Mesh Mould
【DINO爱编程】几何_用菱形填充正六边形问题
【数学原理】Tiling a Hexagon with Diamonds
图形程序的下载地址：
Tiling a Hexagon with Diamonds.zip
1989 年的《美国数学月刊》（American Mathematical Monthly）上有一个貌似非常困难的数学问题：下图是由一个个小三角形组成的正六边形棋盘，现在请你用右边的三种（仅朝向不同的）菱形把整个棋盘全部摆满（图中只摆了其中一部分），证明当你摆满整个棋盘后，你所使用的每种菱形数量一定相同。
文章末尾提供了一个非常帅的“证明”。把每种菱形涂上一种颜色，整个图形瞬间有了立体感，看上去就成了一个个立方体在墙角堆叠起来的样子。三种菱形分别是从左侧、右侧、上方观察整个立体图形能够看到的面，它们的数目显然应该相等。
为了证明这个原理，通过Pascal进行编程生成不同的菱形布置情况，通过不同的颜色展示，你会发现这个二维铺砖问题变成了一个三维模型生成的问题，而这些模型非常像城市规划或建筑师布置塔楼与裙房方案的试算。还有一个特点就是生成的三维模型不存在上部大而下部小的情况，都是稳定的长方体布置。一个从二维问题变成三维布置的“升维”。
生成填满菱形的正六边形的编程关键在于格子的二维化，我采用二维坐标记录这些三角形格子，有纵横坐标，三种不同填充对应三种相邻关系（奇偶行情况不同），如图所示。
这个问题及其鬼斧神工般的“证明”流传甚广，深受数学家们的喜爱。《最迷人的数学趣题——一位数学名家精彩的趣题珍集》（Mathematical Puzzles: A Connoisseur's Collection）一书的封皮上就赫然印着这个经典图形。在数学中，类似的流氓证明数不胜数，不过上面这个可能算是最经典的了。以下附上动画与相关的图案。
【Dino结构笔记】结构工程的造山_不规则壳元的建模
以下分享工程师遇到问题解决问题的过程：
现在越来越多的建筑师与结构工程师可以通过RHINO+GH对复杂曲面进行建模，通过参数化建模最后把模型导入SAP2000等软件进行结构分析。去年有幸参与一个主题乐园项目的结构分析，从中遇到了一种特别的情况，需要把不规则壳元进行有限元分析。事情是这样子的：
目前的建筑设计已经大量应用了3D打印及3D扫描技术，建筑师或设计师，只需要建立缩尺模型，采用高 精度的3D扫描技术，把模型数码化，在电脑的数码格式，可以把缩尺模型放大至实际建筑物或景观的大小，这也就是现在的人造景观的技术。目前已经可以做到“愚公造山”，结构工程师只需要设计一个钢结构 框架，如图所示，外表面可以采用GRC混凝土（带一定受拉承载力的混凝土），这种混凝土往后发展可以成为3D打印结构的主要材料。
钢结构与外蒙皮结构的关系
3D 扫描技术在建模中的应用
岩石山体扫描后数码化的效果
落在手上的问题是，需要对3D扫描的GRC混凝土壳进行有限元分析，与传统的参数化生成的壳元不同，3D扫描得到的壳元的点分布是杂乱无章的（3D扫描是得到点云，Point Cloud，这与激光的频率，点与点的距离有关，没有什么规律），需要把壳元规划划分后导入SAP2000进行计算。因为这个项目少说也有几百块壳体，批量化的划分单元是非常重要的。结构工程师遇到复杂模型的建模问题，一般通过编程解决，与传统的网格划分技术不同，这次是把网格规则化。在这个过程中，享受采用小学二年级学习的空间几何知识解决问题的快感。
以上外壳为例，在3DSMAX打开壳元的实际网格如图所示。
自编小程序用于批量划分网格 3D surface Mesh.
1）编制程序读取3D扫描得到的不规则曲面
2）生成正交规范的点坐标
3）求解点坐标到三角面内的映谢，得到正交点的Z轴坐标。
4）最后通过正交点的三维坐标形成新的规则曲面
5）另存为DXF，导入SAP2000进行静力分析，如图所示。
6）由于有很多块壳元，写个LISP就可以自动操作起来了
7）特别三轴复杂的曲面，需要人为处理一下
规划划分后的SAP2000 模型
这里用到的小学二年级的空间知识居然是已知三点坐标求得平面方程，如下图所示。有时候简单的数学拯救了工程师。这个是对去年的有限元网格划分在工程中应用的一次小总结。现在有了Rhino与GH,这种问题又有了新的解决方法，但是对于大量单元的网格划分，如有代码进行后台计算的效率，而且还可以打包成容易入手的小工具。
最后发一张照片留念一下。
以上就是这篇笔记的全部内容。
以上操作用到的AutoCAD的二次开发，SAP2000的二次开发，三维网格划分的知识，空间几何的知识等。
==============================
这个是我的微信公众号[Dino结构笔记]：DINOSTRU
也可以通过扫二维码进行关注一下。
[程序]Graphic Statics 桁架图解静力法绘图程序
程序下载：点击下载此文件：Graphic Statics 桁架图解静力法绘图程序
使用说明：点击下载此文件：
Graphic Statics 桁架图解静力法绘图程序.PDF
下载地址：http://www.dinochen.com/article.asp?id=355
Graphic Statics 桁架图解静力法使用说明
1) Graphic Statics程序的介绍
以下通过一个实例介绍这个程序的操作与建模。
2)程序操作说明
（1）打开ETABS，进行桁架结构的建模，采用的二维立面是X-Z立面，建模完成后，对桁架结构进行三角面的蒙皮操作，如图所示。模型中注意要施工节点约束。
建立ETABS模型
（2）在ETABS中，施加重力荷载，如图所示，采用点荷载输入，荷载工况名字为 DEAD 恒载工况。然后建模完成，导出S2K文件，如图所示。
施加DEAD工况的节点荷载
输出S2K文件
（3）打开Graphic Statics程序，点击按钮【Load Etabs S2k File】导入刚才ETABS生成的S2K文件，即可以导入桁架模型到程序当中。导入成功后，会出现在右图框内。
（4）点击按钮【Run Graphic Analysis】，程序自动计算桁架的构件内力并绘制图解静力法的图形，如下图所示。根据图解法（左图），每个节点代表桁架内的每个三角形，节点与节点的线，代表两个三角形共线的边的杆件的轴力（有正有负），内力之和等于外力之后，所以所有线的矢量合为0，也就是一个封闭的图形。
图解法的图与桁架布置图
（5）在右图中可以选取节点，采用键盘的WSAD四键可以上下左右偏移节点，每次偏移后，程序会计算桁架的静力图及计算虚功总值，在【Energy】的框中显示，值越小证明构件的布置越是优化。如图所示原方案的虚功总值的初始值是5.8483e6，如果移动节点得到以下图形，则虚功总值为3.7547e6，证明这个桁架更加优化。(从图解法的图来看，一开始的力线是一系列的平行线，优化后的力线是一系列的放射线，放射线越多图解法的图形就越扁，虚功总值就会越小)。
（6）关于图解静力法，可以参考我之前的帖子：
【Dino结构笔记】结构图解法的学习与编程
http://dinochen.com/article.asp?id=338
这次主要补充这一个小程序，帮助更好的理解图解静力法。想要了解如何绘制图静静力法的图，可以看我的操作视频：https://www.bilibili.com/video/BV1Lf4y1U7FS
最后分享一下SOM建筑师与工程师通地图解静力法设计的桁架结构，体验一下结构设计之美 。
【DINO爱编程】几何系列_派生的多边形
【数学原理】derived polygon
图形程序的下载地址：derived_polygon.rar
伐里农的平行四边形是通过连接任意四边形每条边的中点而派生出来的。
这个定理同样可以延伸到其它多边形上。
你可以随意地画出一个不规则的多边形，将这个多边形每条边的中点连在一起，就会派生出一个中点多边形。
让人惊讶的是，如果你继续按照这样的方式去做，那么派生出来的多边形就会变得接近正多边形，它们的每条边在长度上也越来越接近。更让人惊讶的是，同样的多边形按不同的比例去分割每条边（比例为派生derived ratio），最终派生出来的多边形也是相似的。
我举个例子，左图是一个很像三角形的四边形，经过几次迭代以后就可以得到红色的平行四边形
右图是一个各边长度相差很大的六边形，经过几次迭代后就变成了一个斜切的六边形
为了测试这个东西，我编写了一个图形小程序用来验证并绘图,界面如图所示。
大家可以绘制不同的多边形去测试，图右是一个计算各边线长度的统计，可以看出是否两两边是相等的。
如上图，经过20次迭代后，各边长度相差很少是一个正多边形了。
derived polygon 派生的多边形 其实就是一个斜切的多边形
以六边形为例，两对边是相等的六边形。四边形就是平行四边形。
其中 Node Number是多边形的节点数，Iteration是迭代数量， derived ratio的新点在线上的比例。
最后，这个现象有没有数学高手可以给出一个漂亮的证明，证明出为什么会出现这个现象呢。
[程序]DXF to Xtract 任意截面生成器
【CHANGE LOG】
*感谢网友 李昱 同学发现问题，我找到一个小bug，修复了一下，以下是更新的版本。
程序下载：
点击下载此文件：DXT to XTRACT version2020-11-06_建模小工具
程序下载：点击下载此文件：DXT to XTRACT_建模小工具
使用说明：
点击下载此文件：DXT to XTRACT操作说明PDF
下载地址：http://www.dinochen.com/article.asp?id=353
DXF to Xtract使用说明
1)Xtract CAD程序说明
Xtract 是一个非常好的截面分析程序。Xtract适合于建立非常复杂的组合截面，包括混凝土与钢材的组合材料，如果采用AUTOCAD进行前处理建模，建立复杂截面，推荐使用该项程序，由于CAD存在不同的版本，之前开发的CAD TO XTRACT不好使用，该程序是基于CAD输出的DXF文件为基础，建立复杂的构件截面，输入DXF文件，然后采用该程序就可以建立复杂的截面了，该程序需要使用AUTOCAD DXF 2000版，只要是2000年以后的AUTOCAD版本就可以使用了，免费的小程序请自行下载。
2)程序操作说明
（1）打开AUTOCAD（任意版本），建立你要分析的复杂截面，可以是多块混凝土，多块钢材截面组成。如下图所示。最后一定要另存为AUTOCAD的2000年的DXF文件格式。
注意：混凝土与大型的钢材一定要是头尾相连的POLY LINE ,多义线。而钢筋一定是圆形（Circle）。在CAD里面不需要明确图层。可以点击【帮助】采用【生成DXF模板范例】，生成一个范例的DXF文件可供练习。
为了简化操作，所有的钢筋采用圆，材料为REBAR
不需要选择与指定，钢筋的实际直径直接读取CAD图中圆的直径。
（2）打开DXF TO XTRACT 小程序，点击 【导入DXF文件】按钮，把刚才输出的文件导入。如下图所示。
（3）按【ctrl+S】或点击按钮 ，可以选取每个多边形（多义线），然后按【ctrl+B】或点击按扭 ，赋予多边形属性，如下图所示，可以赋予多边形不同的混凝土材料，保护层材料，保护层厚度及划分网格大小等属性，点【确定】退出。材料编号有3种混凝土材料，有3种钢板材料，具体参数可在XTRACT设置。
（4）最后点击【生成XTRACT】按钮 ，就可以生成（*.sec）的XTRACT文本格式文件，可以导入XTRACT程序当中。如下图所示。
（5）将SEC1.SEC导入XTRACT以后，就可以看到以下结果，表明截面导入成功。
（6）接下来可以在XTRACT程序中修改材料属性参数，最后完成截面分析，具体在XTRACT如何分析，请看XTRACT的说明书，或者看这个视频教程：
[Dino公开课]截面分析XTRACT程序介绍
http://dinochen.com/article.asp?id=244
[小记]至从在华南理工大学博士毕业以来相隔好多年没有发表工程设计相关的论文了，这篇论文总结了一个在2011年做的一个大连的连体项目的设计经验。处于7度区（安评地震加速度偏大）的连体超高层，在转换桁架的设计中要充分忽略上部混凝土框架的有利贡献，在模型中进行修正处理。这个项目同时采用的搭接柱与单向少墙，也是一个两个方向不同体系的超高层项目，虽然高度不高，但是难度还是可以的。这个项目虽说是大连的项目也经过北京全国专家的审查。另外，搭接柱虽然不是一种好的受力方式，但由于往后建筑功能的要求，以后的项目也是会大量采用的，关于搭接柱的设计方法目前来说还是以经验为主，希望往后的规范把搭接柱的设计细化。
【论文下载PDF】大连某高层建筑高位连体结构设计.pdf
| 【英文篇名】 | Design of High-rise connective structure for a high-rise buidling in Dalian |
| 【中文篇名】 | 大连某高层建筑高位连体结构设计 |
| 【作者】 | 陈学伟 黄汉华 |
| 【英文作者】 | Chen Xuewei, Henry Wong |
| 【作者单位】 | WSP Hong Kong LTD, WSP科进香港顾问有限公司 |
| 【文献出处】 | 建筑结构 , Building Structure, 编辑部邮箱 未知 网络首发 期刊荣誉：中文核心期刊要目总览 ASPT来源刊 中国期刊方阵 CJFD收录刊  |
| 【中文关键词】 | 超限高层建筑、连体结构、转换桁架、高宽比、弹塑性时程分析 |
| 【英文关键词】 | out-of-code high-rise building; connective structure; transfer truss; aspect ratio; elastoplastic time history analysis |
| 【摘要】 | 某高层酒店项目位于大连市，是包含多栋复杂高层建筑的大型综合发展项目。介绍 B 栋塔楼的连体结构设计。该塔楼结构高度为 187.25m，地上 59 层，地下 4 层。结构采用带部分钢骨的钢筋混凝土框架剪力墙体系。结构的 Y 向高宽比达到 8.9，该方向的刚度不足，通过在两端设置了“山墙”构件提高整体抗侧刚度，并通过设置钢骨提高构件的抗震抗拉性。 结构属于高位连体结构，在 29 层设置巨型转换桁架。采用性能化的分析方法对转换桁架进行重力与抗震分析，以满足关键构件性能目标的要求。最后采用大震弹塑性时程分析方法检验结构及关键构件在大震作用下的性能。 |
| 【英文摘要】 | Located in the Dalian, the project is a large-scale comprehensive development project with multiple complex high-rise building. This paper introduces the structural design of tower B which is the connected structure. The structural height of tower B is 187.25m with 59 floors above the ground and 4 floors under the ground. Tower B is reinforced concrete frame-shear wall structure system with some SRC members. The spect ratio in Y direction reaches 8.9, indicating insufficient stiffness. The gable walls are set at both ends to improve the overall lateral stiffness, and steel is set to improve the seismic tensile resistance of the members. A mega transfer trusses are set on 29h floor and it meet high location loading transfer problem. Performance-based analysis method is applied to analyze the gravity and seismic behavior of the transfer truss, make sure it meet the performance target as the key member. Finally, nonlinear time history analysis is adopted to review the performance of the structure and key members under severe earthquake level. |
| 【DOI】 | CNKI:SUN:JCJG.0.2011-04-018 |
| 【更新日期】 | 2020-10-29 |
| 【分类号】 | TU318.2 |
[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站视频
中国知网自动生成的参考文献格式（更新于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原理与实例》。新书已在某宝某东某多多发售，您们的支持是我们继续创作的动力！希望这本书能帮助有需要的学生、研究人员和工程师！谢谢大家！！ |
|
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站视频
[每周挑战]造物_造一个使命必达盒子(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站。
[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
也可以通过扫二维码进行关注一下。
[每周挑战]造物_造一个黄油机器人(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
