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


博学慎思,坚持创新,求知若飢,虛心若愚                       Keyword Type
  • 博主
  • 主页
  • 日志
  • 研究
  • 论文
  • 自编程序
  • 工程
  • Dino笔记
  • ETE专题
  • DinoBox
  • OpenSEES
  • Dino公开课
  • 每周挑战
  • Dino爱编程
  • 留言板

DinoBox

DinoSap结构力学求解器V1.0的介绍:
(1)DinoSAP2D简单易用,全图形操作,不需要命令流
(2)提供 弹性模量E 截面面积A 与截面抗弯惯性矩 I 的输入
(3)显示的弯矩图等计算结果
(4)颜色,字体,弯矩图的内线都可以自已调整,可以出计算图形与内力图形
(5)计算快,出结果人性化,内力细节查询
(6)由于操作简单,可以直接看视频就学会了
(7)全过程是:结点、截面材料E A I 、单元、支座、荷载、分析得内力。
结构力学求解器
Author: dinochen1983
Time: 2020-12-04
Graphic Statics是基于Pascal语言编写的一个基于桁架图解静力法求解桁架内力的一个小程序,采用ETABS软件对桁架进行建模,导入Graphic Statics桁架图静静力法程序中,程序自动计算桁架内力并绘制桁架对应的静力求解图,图解法的静力求解图对桁架结构的优化具有很好的工程意义。为了更好的进行优化,该程序可以实现实时修改节点坐标后得到更新的静力求解图,且能得到桁架的虚功总值,通过对比虚功总值可以得到桁架的形态是否优化。
桁架图解静力法
Author: dinochen1983
Time: 2020-11-08
[程序]DXF to Xtract 任意截面生成器
DXF2Xtract
Author: dinochen1983
Time: 2020-11-01
上一篇 | 下一篇
Dino笔记 订阅所有Dino笔记的日志

[Dino笔记]遗传算法在结构优化分析中的应用

作者:dinochen1983 日期:2020-12-06

字体大小: 小 中 大

【Dino结构笔记】遗传算法在结构优化分析中的应用


【点击下载】下载相关的ETABS与OPENSEES模型文件

结构优化的方法有很多种,中间有的涉及到非常系统的工程优化算法,如之前介绍的单纯体法(算是比较简单的),在这里推荐大家一本好书《工程优化方法》陈卫东写的。
结构设计中很多时侯我们并不是要得到最最优解,相对优化对工程师就很不错了,但是太复杂的算法会让工程师很难把控,特别是计算时间长,因此遗传算法作为只是找到相对最优解的方法在工程上能很好地被利用。正所谓“遇事不决,量子力学,不知解法,遗传算法”。
 
遗传算法的背后思想就是优胜劣汰的自然法则z


在这里我们只是把遗传算法的思想应用在结构工程上,可以说是似遗传算法,并不是工程优化研究或算法研究中的复杂的遗传算法。在MARK SARKISIAN的《DESIGNING TALL BUILDINGS》(非常推荐大家去看)的书上提到这种遗传算法的应用。
 
MARK SARKISIAN的《DESIGNING TALL BUILDINGS》应用介绍

那么这一期的Dino结构笔记就通过OPENSEES, ETABS等软件结合编程实现一个遗传算法的计算,以一个桁架结构为例进行介绍吧。

实例:对一定总重量的钢桁架结构进行结构优化,使桁架在重力荷载作用下的跨中变形最小,也即是刚度最大。桁架上弦的节点的作用力为150kN,由于桁架是对称的,采用半跨建模,初始条件,全部截面的尺寸为4000mm2,统计一下,半跨的桁架的总重量约为2.285 ton。
计算简图如下图所示。通过遗传算法确定每一节截面的大小,以保证在用钢量一定的情况下总刚度最大。
 
桁架的计算简图

(1) 这个实例我打算采用OPENSEES进行分析,那么首先我们建立ETABS模型,加入支座与荷载条件,如下图所示,导出S2K文件以供导入OPENSEES所用。
 
桁架在ETABS中建模

(2) 导入ETO生成OPENSEES的分析命令流文件,导入ETO以后处到桁架的构件编号如图所示。在OPENSEES中我们只需要输入点13的竖向位移,做为遗传算法的评分标准,点13的位移越小,评分越高,代表遗传算法中的后代更加优秀。
 
桁架在OPENSEES中的建模


(3) OpenSEES的命令流如下所示,命令流分为3个部分:
第一部分:节点、材料与支座条件(分析过程中不发生改变)
第二部分:截面(分析过程中不断改变),采用引入文件  source section_1.tcl
第三部分:分析与结果提取,结果提取的命字改变,即:
recorder Node -file node1.out -time -node 13 -dof 3 disp
不同的个体分析存储的文件名不同, 分别是 node1~node8.out, 记录13号点的竖向位移

model basic -ndm 3 -ndf 6
node 1 0.000E+000 0.000E+000 0.000E+000
……….
node 17 1.000E+004 0.000E+000 0.000E+000
fix 1 1 1 1 1 1 1;
……….
fix 17 0 1 0 1 1 1;
uniaxialMaterial Elastic 1 1.999E+005
source section_1.tcl
recorder Node -file node1.out -node 13 -dof 3 disp
pattern Plain 1 Linear {
load 4 0.000E+000 0.000E+000 -1.500E+005 0.000E+000 0.000E+000 0.000E+000
……….
load 15 0.000E+000 0.000E+000 -1.500E+005 0.000E+000 0.000E+000 0.000E+000
}
constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.0e-6 200
algorithm Newton
integrator LoadControl 1
analysis Static
analyze 1

(4) 桁架截面在遗传算法的过程中是不断变化的,所以这个文件需要不断的修改,8个子代就有8个不同的文件存储,如其中一个子代的截面文件 section_1.tcl如下所示。
element truss 1 2 3 9267.198 1
………
element truss 31 17 11 1215.233 1

(5) 基本的OPENSEES的文件构造就搭好了,只要Run OpenSEES,程序就会计算出不同截面布置情况下的13号节点的竖向位移,我们的目标就是让这个位移变小。

整个简化遗传算法的流程图
(6) 接下来我们就需要通过编程实现遗传算法,并控制OPENSEES不断计算。这里我采用DELPHI进行编程,编了一个简化的遗传算法逻辑。流程图如下图所示。
遗传算法的主要内容如下:
a). 随机生成8个母代桁架结构模型,截面是随机大小分布,总重量恒定
b). 采用OPENSEES作为计算核心计算出8个母代桁架的位移作为评分标准
c). 对母代桁架进行排序,位移最小则最优,得到较大的交配概率为26.7%
d). 根据交配概率抽取2个母代桁架进行杂交,生成新的子代桁架,进行8次
e). 子代桁架需要出现变异,以保证多样性,抽取个别构件截面放大150%与缩小66.7%,以这个操作代表变异性。
f)把子代进行OPENEES计算,即执行b的操作,整个过程循环操作。

最后遗传算法的程序编好,如图所示,迭代次数与13号节点的位移曲线如图所示。

证明,随着遗传算法的不断迭代,桁架的刚度越来越大。13号节点的位移从413降至187mm。

 
 
基于OPENSEES的遗传算法控制程序
 
迭代次数与竖向位移值的关系


(7) 采用ETABS进行分析:原方案与最后优化方案的用钢量是相等
原方案:全部杆件的截面为4000mm2,13号节点竖向变形为350mm
优化方案:每个截面如下图表所示,13号节点竖向变形为185mm。刚度有了明显的提高。而这个算法的优化方向是指定的,虽然有可能不是最优,但是足够应用于工程了。

 
原方案的ETABS计算结果

 
优化方案的ETABS计算结果
 
优化方案的ETABS计算结果(虚功分布图)


总结:大自然总能给工程师与科学家启发,特别是优化算法的研究领域,如蚁群算法,狼群算法,人工神经网络算法等等。自然经过以亿年为单位的时间进行演化,从中有太多的规律可以被发现并为工程师学习与应用。作为工程师通过简单的编程实现遗传算法在工程计算中的应用,那一种喜悦是非常特别的。

             
 
我与桁架结构的合照



[本日志由 dinochen1983 于 2020-12-06 03:13 AM 编辑]
文章来自: 本站原创
引用通告: 查看所有引用 | 我要引用此文章
Tags:
评论: 0 | 引用: 0 | 查看次数: 459
发表评论
昵 称:
密 码: 游客发言不需要密码.
验证码: 验证码
内 容:
选 项:
虽然发表评论不用注册,但是为了保护您的发言权,建议您注册帐号.
字数限制 200 字 | UBB代码 关闭 | [img]标签 关闭

信息(Information)

DinoChen微博与QQ空间
  
DinoChen 视频频道
  


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


dinochen.com需要你的支持

目录(Category)

[01]. 陈学伟Dino的介绍
[02]. 主页: 博客的全部内容
[03]. 日志: 生活日常记录
[04]. 研究: 结构工程研究
[05]. 论文: 结构工程论文
[06]. 程序: 自编程序分享
[07]. 工程: 结构设计经验
[08]. DinoBox:结构软件
[09]. ETE专题:结构分析
[10]. OPENSEES教程
[11]. 公开课视频教程
[12]. Dino每周挑战栏目
[13]. Dino爱编程栏目
[14]. 网站的留言板

日历(Calendar)

2021年1月
  • 日
  • 一
  • 二
  • 三
  • 四
  • 五
  • 六
  • 27
  • 28
  • 29
  • 30
  • 31
  • 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

链接(Blogs)

http://www.dinochen.com/attachments/month_19-08/920198103733.jpg
[009]. 微望建筑科技
[010]. 崔济东的博客
[011]. 陆新征教授的网站
[012]. 董翰林的个人网站
[014]. 贾苏的博客
[014]. 石岩博士研究小组
[020]. 补国斌的博客
[040]. 王斌的建筑结构博客
[050]. 种菜编程盖房子研习中医
[090]. 徐珂的博客
[100]. 鲁班人的博客
[120]. 王文达博士研究小组
[210]. 十千牛博客
[220]. 国巍结构抗震课题
[230]. 陈维晃的个人博客
[240]. Wei_Jiebin的博客
[260]. Html在线编辑器

需要交换链接的朋友
请写邮件给我
dinochen1983@qq.com

视频(MyVideo)

【BiliBili视频】  

【OPENSEES在结构中的应用】

查找(Search)

关键字
 
类 型

标签(Tags)

Loading...

评论(Comments)

圆形截面怎么不行呢 可以调成中文吗?请问软件还有渠道可以获取吗,请问前辈找到转换方法了吗?我也是学桥梁的回复:是16个柱,这个是验证模型不是工程模型,采用...文中3.1节第二段似乎应为18根柱,而不是16根柱...软件是否需要购买? 回复:软件没有更新,现在不能...http://www.dinochen.com/a...陈博士,您好。好像没有最新的ETO 的下载链接,请...陈博士,能否有施加集中时程力的方法,简单一点的

超链接(SLinks)

申请 | 查看多链接…

文章(Recent Article)

[Dino工程]深圳坪山华侨城项目(246.9m)...[分享]另类结构形式的思考—柔软的形式[Dino笔记]遗传算法在结构优化分析中的应用[DinoBox]DinoSap结构力学求解器(1...[论文]深圳大百汇主塔楼结构设计与分析[分享]Rhino+GH平台在结构设计中的应用[Dino爱编程]几何_用菱形填充正六边形问题[Dino笔记]结构工程的造山_不规则壳元的建模[程序]Graphic Statics 桁架图解静... [Dino爱编程]几何_派生的多边形

统计(Statistics)

日志: 351 篇
评论: 802 个
留言: 325 个
访问: 10676362 次
在线: 3 人
建站时间: 2008-02-08

留言(GuestBook)

陈博士好,我是研二的在读学生,在做装配式桥墩方向,...陈博士您好,我使用的是etabs9.7.4,导出的...陈博士好,我下了您的DXF to Xtract小工...陈博士您好!我想请教一下,我在做柱的低周往复模拟时...陈老师您好,麻烦老师给我发一个ETABS 9.7,...陈博士,第二版书中的例子绝大多数都是采用linea...请问一下陈博士,opensees是如何实现像aba...请问陈博士,当框架节点过多时,进行模态分析后,一次...陈老师,您好,我电脑上的Etabs 2018 没办...如果使用老版本的OpenSEES(如V2.5.0和...

用户面板(User Panel)

登录用户注册

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

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