切换到宽版

 找回密码
 注册
Simdroid-非首页
Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛---(手机验证注册)»论坛 › ..:: ANSYS专区 › C03：ANSYS--结构软件（Mechanical\Stuctural\Rigid Dynamic › ANSYS optimization 例子
123下一页
返回列表 发新帖
查看: 2736|回复: 40

ANSYS optimization 例子

[复制链接]
S-Y.Chen
发表于 2004-11-16 01:00:25 | 显示全部楼层 |阅读模式 来自 台湾
之前好像有人在談到ANSYS optimization是否好用或可用的問題

我實在不想讓大家失望,可是ANSYS 的optimization,根據我自己的瞭解跟經驗,一來無法作非線性度高的題目,二來設計變數一多也很難收斂

這裏貼兩個範例,取自NASA某研究員的paper中的兩個例子. 這兩個例子ANSYS都無發收斂取得合理解. 我分兩篇回複來貼. 引用的論文在每個log檔的檔頭.

S-Y. Chen
ww.FEA-Optimization.com
ansys, Optimization, 例子

相关帖子
回复

使用道具 举报

S-Y.Chen
 楼主| 发表于 2004-11-16 01:01:22 | 显示全部楼层 来自 台湾台北

Re:ANSYS optimization 例子

Simdroid开发平台
! ========================================================
! Create by S-Y. Chen, 20001209   www.FEA-Optimization.com
! Data obtained from Dr. Patnaik at NASA @20041208
! ========================================================
! References
! ----------
! S.N.Patnaik, "Comparative Evaluation of Different Optimization Algorithms
!   for Structural Design Applications", International Journal for Numerical
!   Methods in Engineering, V39 pp1761-1774, 1996
!
! This is test problem P2 for NASA COMETBOARD
! ========================================================
! This model has shown that ANSYS has difficulty of convergence for certain problems
! ========================================================
!
*CREATE,P2_ANS
!
/FILNAME,P2_ANS
!
! ===========================================================
! INITIALIZE DV
DV001=1.0
DV002=1.0
DV003=1.0
DV004=1.0
DV005=1.0
DV006=1.0
DV007=1.0
DV008=1.0
DV009=1.0
DV010=1.0
!
/PREP7
N,    1,        0.00000,        0.00000,        0.00000
N,    2,      162.50000,        0.00000,        0.00000
N,    3,      325.00000,        0.00000,        0.00000
N,    4,      325.00000,       50.00000,        0.00000
N,    5,      162.50000,       62.50000,        0.00000
N,    6,        0.00000,       75.00000,        0.00000
!
ET,1,1
ET,2,21
MP,  EX,1,10000000
MP,NUXY,1,0.3  
MP,DENS,1,2.58799e-4
MP,ALPX,1,1.0E-7
!
!
R, 1,DV001
R, 2,DV002
R, 3,DV003
R, 4,DV004
R, 5,DV005
R, 6,DV006
R, 7,DV007
R, 8,DV008
R, 9,DV009
R,10,DV010
!
! ==================================
! TRUSS
TYPE,1
MAT,1
REAL,1
EN,    1,   1,    2,
REAL,2
EN,    2,   2,    3,
REAL,3
EN,    3,   3,    4,
REAL,4
EN,    4,   4,    5,
REAL,5
EN,    5,   5,    6,
REAL,6
EN,    6,   2,    5,
REAL,7
EN,    7,   2,    6,
REAL,8
EN,    8,   1,    5,
REAL,9
EN,    9,   3,    5,
REAL,10
EN,   10,   2,    4,
!
! ===================================
! MASS
ET,2,21
KEYOPT,2,3,4
R,102,5.0285E-4
R,103,9.0512E-4
R,104,5.0285E-4
R,105,9.0512E-4
TYPE,2
REAL,102
EN,102,2
REAL,103
EN,103,3
REAL,104
EN,104,4
REAL,105
EN,105,5
!
! ====================================
! TOTAL VOLUME
TTLVOLU=0.0
*GET,ZTMP0001,ELEM,1,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,2,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,3,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,4,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,5,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,6,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,7,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,8,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,9,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,10,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
!
!
! ========================
! MODAL ANALYSIS
/PREP7
ANTYPE,MODAL
/SOLU
EQSLV,FRONT
LSCLEAR,ALL
D,1,UX,0
D,1,UY,0
D,6,UX,0
D,6,UY,0
MODOPT,SUBSP,5
SOLVE
/POST1
SET,1,1
*GET,LS03FQ01,MODE,1,FREQ
*GET,LS03FQ02,MODE,2,FREQ
*GET,LS03FQ03,MODE,3,FREQ
!*GO,INE01
!
!
/PREP7
ANTYPE,STATIC
/SOLU
EQSLV,FRONT
!
! ========================
! Load step 1
LSCLEAR,ALL
/INPUT,P2_T00,s01,../
solve
!
! ========================
! Load step 2
LSCLEAR,ALL
/INPUT,P2_T00,s02,../
SOLVE
!
!
/POST1
SET,1
ETAB,LS1SAXL1,LS,1
*GET,LS01ES01,ELEM, 1,ETAB,LS1SAXL1
*GET,LS01ES02,ELEM, 2,ETAB,LS1SAXL1
*GET,LS01ES03,ELEM, 3,ETAB,LS1SAXL1
*GET,LS01ES04,ELEM, 4,ETAB,LS1SAXL1
*GET,LS01ES05,ELEM, 5,ETAB,LS1SAXL1
*GET,LS01ES06,ELEM, 6,ETAB,LS1SAXL1
*GET,LS01ES07,ELEM, 7,ETAB,LS1SAXL1
*GET,LS01ES08,ELEM, 8,ETAB,LS1SAXL1
*GET,LS01ES09,ELEM, 9,ETAB,LS1SAXL1
*GET,LS01ES10,ELEM,10,ETAB,LS1SAXL1
!
*GET,LS01NU01,NODE,3,U,Y
*GET,LS01NU02,NODE,4,U,Y
!
!
SET,2
ETAB,LS2SAXL1,LS,1
*GET,LS02ES01,ELEM, 1,ETAB,LS2SAXL1
*GET,LS02ES02,ELEM, 2,ETAB,LS2SAXL1
*GET,LS02ES03,ELEM, 3,ETAB,LS2SAXL1
*GET,LS02ES04,ELEM, 4,ETAB,LS2SAXL1
*GET,LS02ES05,ELEM, 5,ETAB,LS2SAXL1
*GET,LS02ES06,ELEM, 6,ETAB,LS2SAXL1
*GET,LS02ES07,ELEM, 7,ETAB,LS2SAXL1
*GET,LS02ES08,ELEM, 8,ETAB,LS2SAXL1
*GET,LS02ES09,ELEM, 9,ETAB,LS2SAXL1
*GET,LS02ES10,ELEM,10,ETAB,LS2SAXL1
!
*GET,LS02NU01,NODE,3,U,Y
*GET,LS02NU02,NODE,4,U,Y
!
SAVE
*END
!
! ======================================================
*USE,P2_ANS
!
! ======================================================
! ======================================================
! OPTIMIZATION SESSION
/OPT
!
OPVAR,TTLVOLU,OBJ
OPVAR,DV001,DV,0.01,10000,
OPVAR,DV002,DV,0.01,10000,
OPVAR,DV003,DV,0.01,10000,
OPVAR,DV004,DV,0.01,10000,
OPVAR,DV005,DV,0.01,10000,
OPVAR,DV006,DV,0.01,10000,
OPVAR,DV007,DV,0.01,10000,
OPVAR,DV008,DV,0.01,10000,
OPVAR,DV009,DV,0.01,10000,
OPVAR,DV010,DV,0.01,10000,
!
OPVAR,LS01ES01,SV,-1.0E4,+1.0E4
OPVAR,LS01ES02,SV,-1.0E4,+1.0E4
OPVAR,LS01ES03,SV,-1.0E4,+1.0E4
OPVAR,LS01ES04,SV,-1.0E4,+1.0E4
OPVAR,LS01ES05,SV,-1.0E4,+1.0E4
OPVAR,LS01ES06,SV,-1.0E4,+1.0E4
OPVAR,LS01ES07,SV,-1.0E4,+1.0E4
OPVAR,LS01ES08,SV,-1.0E4,+1.0E4
OPVAR,LS01ES09,SV,-1.0E4,+1.0E4
OPVAR,LS01ES10,SV,-1.0E4,+1.0E4
!
OPVAR,LS02ES01,SV,-1.0E4,+1.0E4
OPVAR,LS02ES02,SV,-1.0E4,+1.0E4
OPVAR,LS02ES03,SV,-1.0E4,+1.0E4
OPVAR,LS02ES04,SV,-1.0E4,+1.0E4
OPVAR,LS02ES05,SV,-1.0E4,+1.0E4
OPVAR,LS02ES06,SV,-1.0E4,+1.0E4
OPVAR,LS02ES07,SV,-1.0E4,+1.0E4
OPVAR,LS02ES08,SV,-1.0E4,+1.0E4
OPVAR,LS02ES09,SV,-1.0E4,+1.0E4
OPVAR,LS02ES10,SV,-1.0E4,+1.0E4
!
OPVAR,LS01NU01,SV,,2.2
OPVAR,LS01NU02,SV,,2.2
OPVAR,LS02NU01,SV,,2.2
OPVAR,LS02NU02,SV,,2.2
!
OPVAR,LS03FQ01,SV,26,
!
! ======================================================
! OPTIMIZATION LOOP
OPANL,P2_ANS
OPTYPE,SUBP
OPSUBP,100
OPPRNT,ON
OPEXE
!

评分

1

查看全部评分

回复 赞 不支持

使用道具 举报

S-Y.Chen
 楼主| 发表于 2004-11-16 01:02:44 | 显示全部楼层 来自 台湾

Re:ANSYS optimization 例子

! ========================================================
! Create by S-Y. Chen, 20001209   www.FEA-Optimization.com
! Data obtained from Dr. Patnaik at NASA @20041208
! ========================================================
! References
! ----------
! S.N.Patnaik, "Comparative Evaluation of Different Optimization Algorithms
!   for Structural Design Applications", International Journal for Numerical
!   Methods in Engineering, V39 pp1761-1774, 1996
!
! This is test problem P9 for NASA COMETBOARD
! ========================================================
! This model has shown that ANSYS has difficulty of convergence for certain problems
! ========================================================
*CREATE,P9_ANS
!
!
/FILNAME,P9_ANS
!
DV001=1.0
DV002=1.0
DV003=1.0
DV004=1.0
DV005=1.0
DV006=1.0
DV007=1.0
DV008=1.0
DV009=1.0
DV010=1.0
!
/PREP7
N,    1,    0.00000,    0.00000,  0.00000
N,    2,  360.00000,    0.00000,  0.00000
N,    3,  720.00000,    0.00000,  0.00000
N,    4,  720.00000,  360.00000,  0.00000
N,    5,  360.50000,  360.50000,  0.00000
N,    6,    0.00000,  360.00000,  0.00000
!
ET,1,1
ET,2,21
MP,  EX,1,10000000
MP,NUXY,1,0.3  
MP,DENS,1,2.58799e-4
MP,ALPX,1,1.0E-7
!
!
R, 1,DV001
R, 2,DV002
R, 3,DV003
R, 4,DV004
R, 5,DV005
R, 6,DV006
R, 7,DV007
R, 8,DV008
R, 9,DV009
R,10,DV010
!
! ==================================
! TRUSS
TYPE,1
MAT,1
REAL,1
EN,    1,   1,    2,
REAL,2
EN,    2,   2,    3,
REAL,3
EN,    3,   3,    4,
REAL,4
EN,    4,   4,    5,
REAL,5
EN,    5,   5,    6,
REAL,6
EN,    6,   2,    5,
REAL,7
EN,    7,   1,    5,
REAL,8
EN,    8,   2,    6,
REAL,9
EN,    9,   2,    4,
REAL,10
EN,   10,   3,    5,
!
! ===================================
! MASS
ET,2,21
KEYOPT,2,3,4
R,102,5.0285E-4
R,103,9.0512E-4
R,104,5.0285E-4
R,105,9.0512E-4
TYPE,2
REAL,102
EN,102,2
REAL,103
EN,103,3
REAL,104
EN,104,4
REAL,105
EN,105,5
!
! ====================================
! TOTAL VOLUME
TTLVOLU=0.0
*GET,ZTMP0001,ELEM,1,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,2,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,3,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,4,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,5,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,6,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,7,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,8,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,9,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
*GET,ZTMP0001,ELEM,10,VOLU
TTLVOLU=TTLVOLU+ZTMP0001
!
!
! ========================
! MODAL ANALYSIS
/PREP7
ANTYPE,MODAL
/SOLU
EQSLV,FRONT
LSCLEAR,ALL
D,1,UX,0
D,1,UY,0
D,6,UX,0
D,6,UY,0
MODOPT,SUBSP,5
SOLVE
/POST1
SET,1,1
*GET,LS03FQ01,MODE,1,FREQ
*GET,LS03FQ02,MODE,2,FREQ
*GET,LS03FQ03,MODE,3,FREQ
!*GO,INE01
!
!
/PREP7
ANTYPE,STATIC
/SOLU
EQSLV,FRONT
!
! ========================
! Load step 1
LSCLEAR,ALL
/INPUT,P9_T00,s01,../,
solve
!
! ========================
! Load step 2
LSCLEAR,ALL
/INPUT,P9_T00,s02,../,
SOLVE
!
!
/POST1
SET,1
ETAB,LS1SAXL1,LS,1
*GET,LS01ES01,ELEM, 1,ETAB,LS1SAXL1
*GET,LS01ES02,ELEM, 2,ETAB,LS1SAXL1
*GET,LS01ES03,ELEM, 3,ETAB,LS1SAXL1
*GET,LS01ES04,ELEM, 4,ETAB,LS1SAXL1
*GET,LS01ES05,ELEM, 5,ETAB,LS1SAXL1
*GET,LS01ES06,ELEM, 6,ETAB,LS1SAXL1
*GET,LS01ES07,ELEM, 7,ETAB,LS1SAXL1
*GET,LS01ES08,ELEM, 8,ETAB,LS1SAXL1
*GET,LS01ES09,ELEM, 9,ETAB,LS1SAXL1
*GET,LS01ES10,ELEM,10,ETAB,LS1SAXL1
!
*GET,LS01NU01,NODE,3,U,Y
*GET,LS01NU02,NODE,4,U,Y
!
!
SET,2
ETAB,LS2SAXL1,LS,1
*GET,LS02ES01,ELEM, 1,ETAB,LS2SAXL1
*GET,LS02ES02,ELEM, 2,ETAB,LS2SAXL1
*GET,LS02ES03,ELEM, 3,ETAB,LS2SAXL1
*GET,LS02ES04,ELEM, 4,ETAB,LS2SAXL1
*GET,LS02ES05,ELEM, 5,ETAB,LS2SAXL1
*GET,LS02ES06,ELEM, 6,ETAB,LS2SAXL1
*GET,LS02ES07,ELEM, 7,ETAB,LS2SAXL1
*GET,LS02ES08,ELEM, 8,ETAB,LS2SAXL1
*GET,LS02ES09,ELEM, 9,ETAB,LS2SAXL1
*GET,LS02ES10,ELEM,10,ETAB,LS2SAXL1
!
*GET,LS02NU01,NODE,3,U,Y
*GET,LS02NU02,NODE,4,U,Y
!
SAVE
*END
!
! ======================================================
*USE,P9_ANS
!
! ======================================================
! ======================================================
! OPTIMIZATION SESSION
/OPT
!
OPVAR,TTLVOLU,OBJ
OPVAR,DV001,DV,0.01,10000,
OPVAR,DV002,DV,0.01,10000,
OPVAR,DV003,DV,0.01,10000,
OPVAR,DV004,DV,0.01,10000,
OPVAR,DV005,DV,0.01,10000,
OPVAR,DV006,DV,0.01,10000,
OPVAR,DV007,DV,0.01,10000,
OPVAR,DV008,DV,0.01,10000,
OPVAR,DV009,DV,0.01,10000,
OPVAR,DV010,DV,0.01,10000,
!
OPVAR,LS01ES01,SV,-1.0E4,+1.0E4
OPVAR,LS01ES02,SV,-1.0E4,+1.0E4
OPVAR,LS01ES03,SV,-1.0E4,+1.0E4
OPVAR,LS01ES04,SV,-1.0E4,+1.0E4
OPVAR,LS01ES05,SV,-1.0E4,+1.0E4
OPVAR,LS01ES06,SV,-1.0E4,+1.0E4
OPVAR,LS01ES07,SV,-1.0E4,+1.0E4
OPVAR,LS01ES08,SV,-1.0E4,+1.0E4
OPVAR,LS01ES09,SV,-1.0E4,+1.0E4
OPVAR,LS01ES10,SV,-1.0E4,+1.0E4
!
OPVAR,LS02ES01,SV,-1.0E4,+1.0E4
OPVAR,LS02ES02,SV,-1.0E4,+1.0E4
OPVAR,LS02ES03,SV,-1.0E4,+1.0E4
OPVAR,LS02ES04,SV,-1.0E4,+1.0E4
OPVAR,LS02ES05,SV,-1.0E4,+1.0E4
OPVAR,LS02ES06,SV,-1.0E4,+1.0E4
OPVAR,LS02ES07,SV,-1.0E4,+1.0E4
OPVAR,LS02ES08,SV,-1.0E4,+1.0E4
OPVAR,LS02ES09,SV,-1.0E4,+1.0E4
OPVAR,LS02ES10,SV,-1.0E4,+1.0E4
!
OPVAR,LS01NU01,SV,,2.2
OPVAR,LS01NU02,SV,,2.2
OPVAR,LS02NU01,SV,,2.2
OPVAR,LS02NU02,SV,,2.2
!
OPVAR,LS03FQ01,SV,26,
!
! ======================================================
! OPTIMIZATION LOOP
OPANL,P9_ANS
OPTYPE,SUBP
OPSUBP,100
OPPRNT,ON
OPEXE
!

评分

1

查看全部评分

回复 赞 不支持

使用道具 举报

光的轨迹
发表于 2004-12-11 12:27:17 | 显示全部楼层 来自 北京

Re:ANSYS optimization 例子

以前来论坛少,没看到这么有价值的帖子,谢谢!
回复 赞 不支持

使用道具 举报

cmelinda
发表于 2004-12-13 03:30:16 | 显示全部楼层 来自 美国

Re:ANSYS optimization 例子

工程优化通常是比较难的问题,尤其对复杂结构而言.但你给的优化例子,对于用ANSYS的工程师而言, 有很大帮助.
THANKS
回复 赞 不支持

使用道具 举报

zhulin 该用户已被删除
发表于 2005-1-27 16:35:10 | 显示全部楼层 来自 黑龙江哈尔滨
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

dolphin993722
发表于 2005-3-22 22:46:21 | 显示全部楼层 来自 陕西西安

Re:ANSYS optimization 例子

有价值阿·！！！
回复 赞 不支持

使用道具 举报

summer_tree 该用户已被删除
发表于 2005-4-7 11:42:49 | 显示全部楼层 来自 西安交通大学
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

leeyuhv 该用户已被删除
发表于 2005-4-9 10:22:53 | 显示全部楼层 来自 山东青岛
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

lbt18872
发表于 2005-4-16 19:43:19 | 显示全部楼层 来自 江苏徐州

Re:ANSYS optimization 例子

非常感谢了
回复 赞 不支持

使用道具 举报

ZHOUYUANKAI 该用户已被删除
发表于 2005-5-8 15:27:56 | 显示全部楼层 来自 黑龙江哈尔滨
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

好市民 该用户已被删除
发表于 2005-9-5 21:51:42 | 显示全部楼层 来自 辽宁阜新
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

bubu 该用户已被删除
发表于 2005-9-5 22:03:07 | 显示全部楼层 来自 陕西西安
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

lggzgg 该用户已被删除
发表于 2005-9-7 21:18:17 | 显示全部楼层 来自 河北秦皇岛
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

wubangdi
发表于 2005-9-8 09:10:20 | 显示全部楼层 来自 湖北武汉

Re:ANSYS optimization 例子

顶！！！
回复 赞 不支持

使用道具 举报

xiadan0924 该用户已被删除
发表于 2005-9-9 21:17:31 | 显示全部楼层 来自 山东聊城
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

caozuozhong 该用户已被删除
发表于 2005-11-1 03:11:46 | 显示全部楼层 来自 四川凉山州西昌
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

zch226 该用户已被删除
发表于 2006-3-13 10:34:10 | 显示全部楼层 来自 东南大学
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

zch226 该用户已被删除
发表于 2006-3-13 10:36:35 | 显示全部楼层 来自 东南大学
提示: 作者被禁止或删除 内容自动屏蔽
回复 赞 不支持

使用道具 举报

qianzd
发表于 2006-3-14 15:48:23 | 显示全部楼层 来自 广东深圳

Re:ANSYS optimization 例子

楼主给我们推荐一个做结构优化的好软件吧，正打算做结构优化问题。
回复 赞 不支持

使用道具 举报

123下一页
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-16 13:35 , Processed in 0.064586 second(s), 18 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表