找回密码
 注册
Simdroid-非首页
查看: 2118|回复: 31

ansys可以做大型结构优化设计吗?

[复制链接]
发表于 2004-12-2 23:10:02 | 显示全部楼层 |阅读模式 来自 北京
看了一些资料好象ansys在结构优化方面很弱,一方面优化算法少,只提供了两中方法,但是在帮助里却说可以大多数的工程问题;另外,ansys在做优化循环的时候都要重新划分网格,但是设计参数改变后划分的网格质量还能保证吗?还有每次循环都要进行的话不是很耗费时间吗?但是,我在帮助里看到了子结构的概念,还有网格尺寸参数话的概念,如果说ansys这些功能都有的话,那么对参数不改变的部分就可以用子结构,为什么不能解决大型复杂结构的优化设计问题?请高手指教!
 楼主| 发表于 2004-12-3 10:01:29 | 显示全部楼层 来自 北京

Re:ansys可以做大型结构优化设计吗?

Simdroid开发平台
还请高手们发表一下看发啊!
 楼主| 发表于 2004-12-3 19:01:07 | 显示全部楼层 来自 北京

Re:ansys可以做大型结构优化设计吗?

怎么没有人做优化吗?大家讨论一下吧!
发表于 2004-12-4 02:20:55 | 显示全部楼层 来自 台湾

Re:ansys可以做大型结构优化设计吗?

Ansys的方法有两种,一是用SUMT (Sequential Unconstrained Minimization techniques),一是用response surface approximation. ANSYS在用这两种方法时,当遇到constraint,都是用penalty function. 也就是将有constraint的函数设法简化成没有constraint的函数,而不是将每一条的constraint都分开讨论.

所以本来的函数如果是
to minimize f(x)
subjected to g(x) <0
ANSYS会将之化成

To minimize f(x)+w*G(x)

所以,当 (1)函数非线性度高 (2)constraint多 (3)设计变量多 时,收敛就会有问题(optimization解的多半是非线性,而且是不等式,不是等式)

另一个问题就是你提的网格变化的问题了.我可以很肯定地告诉你, 如果循环之间的网格无法保持一致,收敛会容易出现问题. Paramesh有维持网格一致的东功能,但是也不知道跟ansys的整合性如何.为了解决这个问题,有所谓的”natural shape function”方法的出现.在我的网站里

www.FEA-Optimization.com > Forum, News, Resource > Resource > Technical Paper

有一篇·  S-Y. Chen, Nov 2002, "Integrating ANSYS with Modern Numerical Optimization Techniques - Part II : A Reverse Parametric Modeling Approach for Structural Shaping Optimization, 2002 Taiwan Area ANSYS Users Conference.
可参考一下

还有ansys多半是取最大应力来作constraint.这也会出问题的. 即使网格可以保持一致的情况下, 只取最大应力都会有问题,何况是网格不能一致?

以下贴两个ANSYS optimization的log file,是很简单的truss结构,但是ANSYS无法收敛.供大家参考…………

S-Y. Chen
www.FEA-Optimization.com

! ========================================================
! 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
!

! ========================================================
! 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,:LINE01
!
!
/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-4 20:51:17 | 显示全部楼层 来自 北京

Re:ansys可以做大型结构优化设计吗?

S-Y.Chen 你好,你说的循环之间的网格无法保持一致指的是什么?解释一下好吗?谢谢!
发表于 2004-12-5 11:07:12 | 显示全部楼层 来自 台湾

Re:ansys可以做大型结构优化设计吗?

所谓的”网格无法保持一致”,指的是这样:
在作优化设计时,程序会依据参数模型来变动设计. 但是在变动几何设计的同时,还要再建立网格才能进行分析. 所以如果是用free mesh来建立网格,两个设计的网格型式可能大不相同. 例如本来是10x10个网格经几合变动后可能变成10x9,或是100个节点变成101个等等. 如此在计算微分及梯度时会有很大的问题.

我的网站上 technical paper底下有几篇有提到
S-Y. Chen
www.FEA-Optimization.com
发表于 2004-12-5 11:09:18 | 显示全部楼层 来自 台湾

Re:ansys可以做大型结构优化设计吗?

所谓的”网格无法保持一致”,指的是这样:
在作优化设计时,程序会依据参数模型来变动设计. 但是在变动几何设计的同时,还要再建立网格才能进行分析. 所以如果是用free mesh来建立网格,两个设计的网格型式可能大不相同. 例如本来是10x10个网格经几合变动后可能变成10x9,或是100个节点变成101个等等. 如此在计算微分及梯度时会有很大的问题.

我的网站上 technical paper底下有几篇有提到
S-Y. Chen
www.FEA-Optimization.com
fox-hu 该用户已被删除
发表于 2004-12-6 09:24:26 | 显示全部楼层 来自 江苏南京
提示: 作者被禁止或删除 内容自动屏蔽
发表于 2004-12-9 00:42:59 | 显示全部楼层 来自 台湾

Re:ansys可以做大型结构优化设计吗?

你说的是subproblem approximation跟first order method吧?first order method就是所谓的linear programming,subproblem approximation是所谓的response surface approximation (ANSYS老是用不一样的名辞……….).我建议是到google去找( 用英文找英文网站),或是找本教科书读吧.

要不找找美国大学的网站,有的会把上课奖义贴上.

评分

1

查看全部评分

发表于 2005-3-15 16:34:22 | 显示全部楼层 来自 江苏南京

Re:ansys可以做大型结构优化设计吗?

我用ANSYS里的优化功能做过一个优化设计,但是我没有用free mesh,而是人工事先划分好网格,然后不管参数如何变化,我的风格划分数量是不变的,只是大小改变,这样也可以得到结果。提出来供大家讨论,欢迎指正。
wuyibin 该用户已被删除
发表于 2006-4-5 18:41:24 | 显示全部楼层 来自 福建泉州
提示: 作者被禁止或删除 内容自动屏蔽
发表于 2006-4-7 18:38:56 | 显示全部楼层 来自 广西南宁

Re:ansys可以做大型结构优化设计吗?

我导师一般就是用ANSYS提取敏度,而用他自己的优化程序进行优化的,但ANSYS有个缺点,就是只能设置20个设计变量。
发表于 2007-4-28 15:32:31 | 显示全部楼层 来自 湖北武汉

回复 #13 moonyu 的帖子

请问一下用ANSYS怎么提取灵敏度?用哪个命令?
发表于 2007-4-29 18:28:26 | 显示全部楼层 来自 湖北武汉
关注中....
发表于 2007-4-29 19:14:54 | 显示全部楼层 来自 北京
Ansys的尺寸优化还是不错的,关键是需要搞清楚优化目标等,我建过一个非常复杂的参数化模型,完全利用APDL,优化结果不错的,以前还做过非线性优化,优化结果和实验误差在1%以内,绝对强到底,可惜我现在不做有限元了,要是ansys公司招我去做优化项目,我绝对要重新走有限元的不归路。

[ 本帖最后由 jyq110 于 2007-4-29 19:18 编辑 ]
发表于 2007-4-29 19:17:54 | 显示全部楼层 来自 北京
不过ansys的拓扑优化实在不敢恭维,还考虑用Altair吧,拓扑优化是十分有前途的优化工具,当时我在给我们师兄弟、老板介绍优化算例的时候,大家都感觉十分不错,可惜啊,我还是没干上我最喜欢的CAE模拟优化工作

[ 本帖最后由 jyq110 于 2007-4-29 19:19 编辑 ]
发表于 2007-4-29 19:20:48 | 显示全部楼层 来自 北京
原帖由 moonyu 于 2006-4-7 18:38 发表
我导师一般就是用ANSYS提取敏度,而用他自己的优化程序进行优化的,但ANSYS有个缺点,就是只能设置20个设计变量。

可以做到60的,我用8.0的可以
发表于 2007-4-29 21:58:08 | 显示全部楼层 来自 四川成都
真的能学会不少,海洋平台应该算大型结构吧
发表于 2007-5-31 21:44:07 | 显示全部楼层 来自 浙江杭州
原帖由 S-Y.Chen 于 2004-12-4 02:20 发表
Ansys的方法有两种,一是用SUMT (Sequential Unconstrained Minimization techniques),一是用response surface approximation. ANSYS在用这两种方法时,当遇到constraint,都是用penalty function. 也就是将有cons ...



弓虽,
形式简单,但程序好长,仰慕一下这位朋友.
回复 不支持

使用道具 举报

发表于 2007-5-31 21:48:41 | 显示全部楼层 来自 浙江杭州
原帖由 S-Y.Chen 于 2004-12-5 11:07 发表
所谓的”网格无法保持一致”,指的是这样:
在作优化设计时,程序会依据参数模型来变动设计. 但是在变动几何设计的同时,还要再建立网格才能进行分析. 所以如果是用free mesh来建立网格,两个设计的网格型式可能大不 ...


能否通俗一点,网格数量差几个会很影响计算精度?
这个观点我不是很赞同.

只要在目标函数、状态变量、设计变量、的提取时,应回避对网格等的依赖性。
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

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

GMT+8, 2024-3-29 18:02 , Processed in 0.071325 second(s), 19 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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