ansys优化设计，最优解总是最初设置的值，不知问题出在哪？

wangkaka

ansys优化设计，最优解总是最初设置的值，不知问题出在哪？ 是不是设计变量，状态变量，目标函数哪块参数设计的不合适啊，有谁遇到过这样的问题，求高手了啊！！！！

wangkaka

顶起！！！！！

hanjixing

我也遇到同样的问题了，目标变量是最初值，状态变量也没变，只有设计变量变了，求高手解答

midasgalaxy

这种情况，很可能是由于没有生成或读入新模型造成的

ant820

可以把APDL发上来看看  我曾经也花了很长一段时间完成过一个ANSYS优化模型，据我的认知范围，ANSYS的常规优化问题最主要的是出在模型建立过程中，模型所有过程都需要参数化

hanjixing

ant820 发表于 2012-3-24 10:51
可以把APDL发上来看看  我曾经也花了很长一段时间完成过一个ANSYS优化模型，据我的认知范围，ANSYS的常规优 ...

我的模型是用壳单元建的（从其他有限元软件画好网格导入ansys），没有参数化建模，但是我提取了实常数即壳的厚度为设计变量，约束条件为最大应力，可以分析 但是目标变量和状态变量都没变，和楼主说的一样，请问怎么解决？

hanjixing

midasgalaxy 发表于 2012-3-24 08:51
这种情况，很可能是由于没有生成或读入新模型造成的

我的模型是用壳单元建的（从其他有限元软件画好网格导入ansys），没有参数化建模，但是我提取了实常数即壳的厚度为设计变量，约束条件为最大应力，可以分析 但是目标变量和状态变量都没变，和楼主说的一样，请问怎么解决？

midasgalaxy

hanjixing 发表于 2012-3-24 22:34
我的模型是用壳单元建的（从其他有限元软件画好网格导入ansys），没有参数化建模，但是我提取了实常数即 ...

把命令流传上来才知道问题在哪里。这个问题并不复杂哈。

hanjixing

midasgalaxy 发表于 2012-3-25 01:32
把命令流传上来才知道问题在哪里。这个问题并不复杂哈。

/BATCH  
/COM,ANSYS RELEASE 10.0    UP20050718       09:42:33    03/25/2012
RESUME,'opt','dbb','.'  
/COM,ANSYS RELEASE 10.0    UP20050718       09:42:36    03/25/2012
/OPT
FINISH  
/SOL
antype,0
solve   
FINISH  
/POST1  
plnsol,s,y  
plnsol,s,eqv
NSORT,U,Y   
*GET,D_MAX,SORT,,MAX
FINISH  
/OPT
! LGWRITE,'opt111','lgw','E:\ANSYSD~1\ANSYST~1\opt2\',COMMENT   
OPANL,'opt111','lgw',' '
!*  
!*  
OPVAR,HBZHL4,DV,2,4,0.02,   
!*  
OPVAR,QBZHL4,DV,2,4,0.02,   
!*  
OPVAR,ZL3,DV,2,3,0.01,  
!*  
OPVAR,HZL3,DV,2,3,0.01,
!*  
OPVAR,LJZJ2,DV, , , ,   
!*  
!*  
OPVAR,D_MAX,SV,1,10,0.09,   
!*  
OPVAR,SMAX,SV,0,200,2,  
!*  
OPVAR,VOLUME,SV, , , ,  
!*  
!*  
OPVAR,VOLUME,OBJ, , ,147478,
!*  
OPTYPE,SUBP
OPSUBP,30,7,
OPEQN,0,0,0,0,0,
!*  
!*  
OPEXE   
/COM,ANSYS RELEASE 10.0    UP20050718       09:46:48    03/25/2012
! OPTIMIZATION LOOPING HAS CLEARED THE INTERNAL LOG
KEYW,BETA,0
!*  
!*  
!*  
OPLIST,ALL, ,0  
!*  



opt.dbb是我昨天用的模型里面各部件的壳厚、最大应力、总体积都已经通过*GET提取，今天导入后又提取了Y向最大位移。
目标变量为总体积(原体积14747837.9)，但是计算很快一会就结束了，状态变量和目标变量都是初始设定值
结果图如下

midasgalaxy

hanjixing 发表于 2012-3-25 09:45
/BATCH  
/COM,ANSYS RELEASE 10.0    UP20050718       09:42:33    03/25/2012
RESUME,'opt','dbb','.'  ...

要先形成参数化APDL文件，这个文件中包含参数化建模过程。如果是数据库文件那么每次使用的都是同一个模型，也就是说，每次优化计算完成后的新设计变量并不能用来驱动模型。

hanjixing

midasgalaxy 发表于 2012-3-25 11:35
要先形成参数化APDL文件，这个文件中包含参数化建模过程。如果是数据库文件那么每次使用的都是同一个模型 ...

可是我的有限元模型是从外部导入的呀，没有参数化过程，直接从导入的有限元模型中提取的厚度

midasgalaxy

hanjixing 发表于 2012-3-25 11:42
可是我的有限元模型是从外部导入的呀，没有参数化过程，直接从导入的有限元模型中提取的厚度 ...

没有参数化过程就是错误原因了。你可考虑如何将这一模型参数化，也就是使设计参数驱动模型。

hanjixing

midasgalaxy 发表于 2012-3-25 12:24
没有参数化过程就是错误原因了。你可考虑如何将这一模型参数化，也就是使设计参数驱动模型。 ...

模型比较复杂，我还是考虑在HYPREWORK里做了。。。参数化好像比较麻烦

midasgalaxy

hanjixing 发表于 2012-3-25 12:32
模型比较复杂，我还是考虑在HYPREWORK里做了。。。参数化好像比较麻烦

我可以确定，你的模型之需要加入不超过10行命令流，问题就能解决。

hanjixing

midasgalaxy 发表于 2012-3-25 12:39
我可以确定，你的模型之需要加入不超过10行命令流，问题就能解决。

还请版主指点迷津。。。具体怎么做？

刘英玉2088

参数都没定义你优化什么啊

deng466904825

midasgalaxy 发表于 2012-3-24 08:51
这种情况，很可能是由于没有生成或读入新模型造成的

能不能指点一下我这为题出在哪？也是同样的问题

finish
/clear
/title,the optimization of the grid structure
/replot
*set,t4,8
*set,t5,6
*set,t6,10
*set,md,7.8e-9
*set,pi,acos(-1)
*set,d1,114                 !杆1-6直径和壁厚
*set,t1,4
*set,d2,88.5
*set,d3,60
*set,t3,3.5
*set,d4,140
*set,t4,8
*set,d5,114
*set,t5,6
*set,d6,159
*set,t6,10
*set,area1,pi*(d1*d1/4-(d1/2-t1)*(d1/2-t1))  !杆1-6截面面积
*set,area2,pi*(d2*d2/4-(d2/2-t1)*(d2/2-t1))
*set,area3,pi*(d3*d3/4-(d3/2-t3)*(d3/2-t3))
*set,area4,pi*(d4*d4/4-(d4/2-t4)*(d4/2-t4))
*set,area5,pi*(d5*d5/4-(d5/2-t5)*(d5/2-t5))
*set,area6,pi*(d6*d6/4-(d6/2-t6)*(d6/2-t6))
*create,opgrid
/PREP7  !进入前处理模块: 定义模型


ET, 1, LINK8    !定义第一类单元为三维杆单元LINK8
MP, EX, 1, 207E3    !定义第一类材料弹性模量EX
MP, PRXY,1, 0.3 !材料密度
MP, DENS, 1, 7.8e-9  !材料密度
R, 1, area1 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米（114*4=1381）斜杆
R, 2, area2 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米（88.5*4=1061）上弦杆
R, 3, area3 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米（60*3.5=621）下弦杆
R, 4, area4 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米（140*8=3315）斜杆
R, 5, area5 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米（114*6=2034）上弦杆
R, 6, area6 !定义杆件第一类实常数--截面积. 该问题的长度单位为毫米（159*10=4678）斜杆


K, 1 ,-5000 ,-5000 , 0    !定义网架的所有结点,单位毫米
K, 2 ,-3000 ,-5000 , 0
K, 3 ,-1000 ,-5000 , 0
K, 4 , 1000 ,-5000 , 0
K, 5 , 3000 ,-5000 , 0
K, 6 , 5000 ,-5000 , 0
K, 7 ,-4000 ,-4000 ,-800
K, 8 ,-2000 ,-4000 ,-800
K, 9 , 0 ,-4000 ,-800
K, 10 , 2000 ,-4000 ,-800
K, 11 , 4000 ,-4000 ,-800
K, 12 ,-5000 ,-3000 , 0
K, 13 ,-3000 ,-3000 , 0
K, 14 ,-1000 ,-3000 , 0
K, 15 , 1000 ,-3000 , 0
K, 16 , 3000 ,-3000 , 0
K, 17 , 5000 ,-3000 , 0
K, 18 ,-4000 ,-2000 ,-800
K, 19 ,-2000 ,-2000 ,-800
K, 20 , 0 ,-2000 ,-800
K, 21 , 2000 ,-2000 ,-800
K, 22 , 4000 ,-2000 ,-800
K, 23 ,-5000 , -1000 , 0
K, 24 ,-3000 , -1000 , 0
K, 25 ,-1000 , -1000 , 0
K, 26 , 1000 , -1000 , 0
K, 27 , 3000 , -1000 , 0
K, 28 , 5000 , -1000, 0
K, 29 ,-4000 , 0,-800
K, 30 ,-2000 , 0 ,-800
K, 31 , 0 , 0 ,-800
K, 32 , 2000 , 0 ,-800
K, 33 , 4000 , 0 ,-800
K, 34 ,-5000 , 1000 , 0
K, 35 ,-3000 , 1000 , 0
K, 36 ,-1000 , 1000 , 0
K, 37 , 1000 , 1000 , 0
K, 38 , 3000 , 1000 , 0
K, 39 , 5000 , 1000 , 0
K, 40 ,-4000 , 2000 ,-800
K, 41 ,-2000 , 2000 ,-800
K, 42 , 0 , 2000 ,-800
K, 43 , 2000 , 2000 ,-800
K, 44 , 4000 , 2000 ,-800
K, 45 ,-5000 , 3000 , 0
K, 46 ,-3000 , 3000 , 0
K, 47 ,-1000 , 3000 , 0
K, 48 , 1000 , 3000 , 0
K, 49 , 3000 , 3000 , 0
K, 50 , 5000 , 3000 , 0
K, 51 ,-4000 , 4000 ,-800
K, 52 ,-2000 , 4000 ,-800
K, 53 , 0 , 4000 ,-800
K, 54 , 2000 , 4000 ,-800
K, 55 , 4000 , 4000 ,-800
K, 56 ,-5000 , 5000 , 0
K, 57 ,-3000 , 5000 , 0
K, 58 ,-1000 , 5000 , 0
K, 59 , 1000 , 5000 , 0
K, 60 , 3000 , 5000 , 0
K, 61, 5000 , 5000 , 0

*DO,I,1,5      !利用循环创建多段下弦杆直线
L, I , I+1   
*ENDDO
*DO,I,12,16   
L, I , I+1   
*ENDDO
*DO,I,23,27   
L, I , I+1   
*ENDDO
*DO,I,34,38   
L, I , I+1   
*ENDDO
*DO,I,45,49   
L, I , I+1   
*ENDDO
*DO,I,56,60   
L, I , I+1   
*ENDDO      
*DO,I,1,6   
L, I , I+11   
*ENDDO
*DO,I,12,17   
L, I , I+11   
*ENDDO
*DO,I,23,28   
L, I , I+11   
*ENDDO
*DO,I,34,39   
L, I , I+11   
*ENDDO
*DO,I,45,50   
L, I , I+11   
*ENDDO        !利用循环创建多段下弦杆直线


*DO,I,7,11    !利用循环创建多段斜杆直线
L, I , I-6   
*ENDDO        
*DO,I,18,22  
L, I , I-6   
*ENDDO
*DO,I,29,33   
L, I , I-6   
*ENDDO
*DO,I,40,44   
L, I , I-6   
*ENDDO  
*DO,I,51,55  
L, I , I-6   
*ENDDO
*DO,I,7,11   
L, I , I-5   
*ENDDO
*DO,I,18,22  
L, I , I-5   
*ENDDO
*DO,I,29,33   
L, I , I-5   
*ENDDO
*DO,I,40,44   
L, I , I-5   
*ENDDO  
*DO,I,51,55  
L, I , I-5   
*ENDDO
*DO,I,7,11   
L, I , I+5   
*ENDDO
*DO,I,18,22  
L, I , I+5   
*ENDDO
*DO,I,29,33   
L, I , I+5   
*ENDDO
*DO,I,40,44   
L, I , I+5   
*ENDDO  
*DO,I,51,55  
L, I , I+5   
*ENDDO        
*DO,I,7,11  
L, I , I+6   
*ENDDO
*DO,I,18,22
L, I , I+6   
*ENDDO
*DO,I,29,33  
L, I , I+6   
*ENDDO
*DO,I,40,44  
L, I , I+6   
*ENDDO
*DO,I,51,55
L, I , I+6   
*ENDDO        !利用循环创建多段斜杆直线



*DO,I,7,10    !利用循环创建多段上弦杆直线
L, I , I+1   
*ENDDO
*DO,I,18,21   
L, I , I+1   
*ENDDO
*DO,I,29,32   
L, I , I+1   
*ENDDO
*DO,I,40,43  
L, I , I+1   
*ENDDO
*DO,I,51,54   
L, I , I+1   
*ENDDO
*DO,I,7,11   
L, I , I+11   
*ENDDO
*DO,I,18,22   
L, I , I+11   
*ENDDO
*DO,I,29,33   
L, I , I+11   
*ENDDO
*DO,I,40,44  
L, I , I+11   
*ENDDO         !利用循环创建多段上弦杆直线


LSEL,S,,,1,60,1, !选择上弦杆指定单元属性
LATT,1,2,1
LSEL,S,,,61,160,1,  !选择斜杆指定单元属性
LATT,1,1,1
LSEL,S,,,161,200,1, !选择下弦杆指定单元属性
LATT,1,3,1
LSEL,S,,,61,136,25, !选择斜杆指定单元属性
LATT,1,6,1
LSEL,S,,,71,146,25,  !选择斜杆指定单元属性
LATT,1,6,1
LSEL,S,,,81,156,25,  !选择斜杆指定单元属性
LATT,1,6,1
LSEL,S,,,63,138,25,  !选择斜杆指定单元属性
LATT,1,4,1
LSEL,S,,,73,148,25,  !选择斜杆指定单元属性
LATT,1,4,1
LSEL,S,,,83,158,25,  !选择斜杆指定单元属性
LATT,1,4,1
LSEL,S,,,1,26,5,    !选择上弦杆指定单元属性
LATT,1,5,1
LSEL,S,,,31,55,6,    !选择上弦杆指定单元属性
LATT,1,5,1
LSEL,S,,,32,56,6,  !选择上弦杆指定单元属性
LATT,1,5,1
LSEL,S,,,2,27,5,   !选择上弦杆指定单元属性
LATT,1,5,1


ALLSEL,all       !网格划分
LESIZE,all,,,1,,,,,1
LMESH,all
!接下来可以利用命令定义边界条件

D,37,all,0   !完全固定结点37，39，41,47,49,51,57,59,61
D,39,all,0
D,41,all,0
D,47,all,0
D,49,all,0
D,51,all,0
D,57,all,0
D,59,all,0
D,61,all,0

ALLSEL$FINISH              !退出前处理模块
/SOLU$ACEL,9800,0,0$SOLVE  !自重为第一荷载步
F,1,FZ,-140000
F,2,FZ,-242500
F,3,FZ,-150000
F,4,FZ,-75000
F,5,FZ,-19000
F,6,FZ,0
F,7,FZ,-280000
F,8,FZ,-485000
F,9,FZ,-300000
F,10,FZ,-150000
F,11,FZ,-38000
F,12,FZ,0
F,13,FZ,-280000
F,14,FZ,-485000
F,15,FZ,-300000
F,16,FZ,-150000
F,17,FZ,-38000
F,18,FZ,0
F,19,FZ,-280000
F,20,FZ,-485000
F,21,FZ,-300000
F,22,FZ,-150000
F,23,FZ,-38000
F,24,FZ,0
F,25,FZ,-280000
F,26,FZ,-485000
F,27,FZ,-300000
F,28,FZ,-150000
F,29,FZ,-38000
F,30,FZ,0
F,31,FZ,-140000
F,32,FZ,-242500
F,33,FZ,-150000
F,34,FZ,-75000
F,35,FZ,-19000
F,36,FZ,0
SOLVE                 !节点荷载为第二荷载步
FINISH                !结束求解状态

/post1
etable,evolume,volu
ssum
*get,vtot,ssum,,item,evolume
etable,smax_e,ls,1
esort,etab,smax_e,0,1,,
*get,smaxe,sort,,max
*get,dzmax,node,53,u,z
dzmax=abs(dzmax)
finish
*end
*use,opgrid


/opt
opanl,opgrid,
opvar,t1,dv,1,5
opvar,t3,dv,3,5
opvar,t4,dv,5,10
opvar,t5,dv,4,8
opvar,t6,dv,7,12
opvar,dzmax,sv,5,30,0.1
opvar,smaxe,sv,180,210,1
opvar,vtot,obj,,,
optype,first
opfrst,25
opsave,,,
opexe
*status
finish

deng466904825

midasgalaxy 发表于 2012-3-25 12:24
没有参数化过程就是错误原因了。你可考虑如何将这一模型参数化，也就是使设计参数驱动模型。 ...

谢谢你的提示，我的问题解决了，我用两个文件来执行读入的问题。可以循环了