24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1966  |  回复: 2

阿斯顿秦文

新虫 (正式写手)

[交流] 如何使用Optistruct进行应力拓扑优化或多目标、多约束优化 已有1人参与

Optistruct是一款非常优秀的商业有限元求解器、优化求解器,功能强大到炸裂,使用起来也很方便。但偶尔用起来也有一点点小麻烦,初学者经常会碰到的问题就是不知道怎么使用Optistruct进行多目标优化或应力优化这种涉及多个响应的优化。Optistruct中的响应是指要作为目标函数或约束函数的结构的性能,比如质量、体积、体积分数、应力、位移等等,其中应力和位移这种响应属于局部响应,即结构中有很多个这种响应,某点的位移或应力不能代表结构的整体性能。以应力优化为例,假如我们想要进行应力最小优化,我们实际上是要使结构中的最大应力最小,但是值得注意的是,优化过程中,具有最大应力的单元一直会变,因此不可能使某个应力值最小,而另外一方面,optistruct也只允许有一个目标函数,怎么办呢?

       通过查询各类资料,本人摸索出一个行之有效的方法,概述如下:

1. 创建一个公式,应力优化经常使用p范数凝聚所有的应力值,结构中有多少个单元,公式中即有多少个未知数

2. 创建NEL个应力响应,NEL为结构中的单元个数

3. 创建一个总的响应,类型选择为:function,勾选第一步创建的公式,然后不要着急create,先点击edit,勾选response,在数目中输入单元个数,然后挨个在弹出的NEL个框中,填入一个个响应

4. 在目标函数中,选择第三步创建的总响应作为目标函数。

        假如这么干的话,难点在于第二步和第三步,因为我们要创建NEL个应力响应,每个响应对应一个单元。第三步则要勾选一个个响应,这些响应要和公式中的未知数一一对应,NEL一般很大,手动创建太难了,为此我们很自然想到使用脚本自动创建这些数据,只要了解Optistruct的语言习惯,自动创建不是难题,本人在研究时,首先创建了基本的优化模型,并验证可运行,然后使用Python脚本自动完成了1-3步。对于下面的经典L型梁




如何使用Optistruct进行应力拓扑优化或多目标、多约束优化
其中有1600个单元,手动创建很不现实,部分响应和公式如下图所示
如何使用Optistruct进行应力拓扑优化或多目标、多约束优化-1
使用本人单独开发的小脚本,很方便完成这部分内容,然后将它们粘贴到之前创建的简单优化模型中即可。

优化过程和优化后的效果如下图所示:
如何使用Optistruct进行应力拓扑优化或多目标、多约束优化-2
拓扑优化后的等值图如下图,可见应力集中的部位基本被消除。本人使用的p范数为6,大家可以尝试其它数值得到更优异的解:
如何使用Optistruct进行应力拓扑优化或多目标、多约束优化-3
而柔度最小化优化的结果是下面这样的,显然和应力最小拓扑优化是不一样的,因为拐角还是直的,没有去除应力集中。
如何使用Optistruct进行应力拓扑优化或多目标、多约束优化-4


创建公式需要注意格式,感兴趣的同学可以在以下链接内
https://www.jishulink.com/conten ... 1-8566-244e8522b4ac
尝试下载附件,附件包括脚本和有效的应力优化模型,谢谢。

同样使用本文的方法还可以求解多目标优化和多约束优化等等,不仅仅是拓扑优化,尺寸优化,形状优化也能解决,因为这三种优化类型都需要响应。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ozfx4857

禁虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

2楼2021-01-24 08:18:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿斯顿秦文

新虫 (正式写手)

引用回帖:
2楼: Originally posted by ozfx4857 at 2021-01-24 08:18:08
难道不能直接使用最小最大化应力吗

你可以试试

发自小木虫IOS客户端
3楼2021-01-27 13:33:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿斯顿秦文 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 hvkbtfonbv 2026-05-23 4/200 2026-05-24 17:21 by 75ui6h7z2t
[教师之家] 论文撤稿了 +4 bjvtcliu 2026-05-24 6/300 2026-05-24 17:11 by jurkat.1640
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:01 by 75ui6h7z2t
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 4/200 2026-05-24 16:16 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 pmo95bazuy 2026-05-23 8/400 2026-05-24 15:56 by 1uy1ht2y9r
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15 popular289 2026-05-18 26/1300 2026-05-21 10:35 by 西葫芦炒鸡蛋
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
请填处理意见