优化中出现的警告问题（加分题,一次性加＋4技术分,直接送你登录ftp的船票）

yajie41

请教大虾了，我在做优化时，采用的模型都不一样，但是最后都会出现这样的警告从而使优化中断，而且警告框里出现的数字都一样，而我根本没有定义过这个叫什么_z3的参数，这大概可能是什么原因？？

*** WARNING ***                        
Unknown parameter name= _Z3.  A value of 7.888609052E-31 will be used.


[ 本帖最后由 北极熊甲 于 2007-8-23 21:18 编辑 ]

北极熊甲

Z3是ANSYS内部程序的一个量，我以前好像在建立一个矩形或别的是什么几何形状时也遇到过这样的提示，就是我故意把某个尺寸给略去。所以感觉还是你的模型在优化时，可能出了某些尺寸参数冲突或矛盾，多了一个形状或少了一个形状，这样的话有可能一些尺寸就不知道了。建议把你的模型给大家说说（最好是贴几个图或动画（你把你的模型旋转一下，在屏幕录下来就行）），优化什么，几何模型可能会发生怎样的变化，你是如何控制的。

yajie41

谢谢版主，那我把我的问题具体说一下，希望版主和各位大虾帮小妹一把！
问题描述:
两块永磁铁放置在三条平行的镍带两端对其进行磁化，为了增大三条镍带中的Y向磁感应强度BY，在两块永磁铁下面放置矩形磁芯。现在就是想优化磁芯的宽度（y向）和厚度（z向）使得镍带中的BY最大。设计变量为磁芯宽度W和磁芯厚度H，目标函数为三条镍带中的BY的总和的倒数。优化方法为零阶方法，可是在优化时，总会出现下面的Warning从而使得优化中断结束，换不同的尺寸模型，连框框里的数字大小都一样，奇怪。
*** WARNING ***                        
Unknown parameter name= _Z3.  A value of 7.888609052E-31 will be used.

我也换过一阶方法，也出现同样的问题。我的镍带尺寸为8×25×0.15mm，很薄，划分网格的时候，
用的是Mesh Tool, 设置Global 中单元尺寸为0.0008m，不知道这样对于很薄的镍带怎样，可是不这样划分，又该怎样弄呢，或者是大家认为导致上面问题的原因到底是什么呢？

我的命令流：
/BATCH  
/COM,ANSYS RELEASE 10.0    UP20050718      
/TITLE,optimization of cores
*SET,W,0.004           磁芯初始宽度
*SET,H,0.004           磁芯初始厚度
*SET,HC,895000         永磁铁矫顽力
/PREP7  
!*  
ET,1,SOLID96           定义单元
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0             定义空气特性
MPDATA,MURX,1,,1
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,2,,250     定义镍带材料特性
TB,BH,3,1,24,          定义磁芯B-H曲线
TBPT,,405,0.80  
TBPT,,468,0.90  
TBPT,,550,1.0   
TBPT,,670,1.1   
TBPT,,840,1.2   
TBPT,,1060,1.3  
TBPT,,1223,1.35
TBPT,,1415,1.40
TBPT,,1722,1.45
TBPT,,2130,1.50
TBPT,,2671,1.55
TBPT,,3480,1.60
TBPT,,4496,1.65
TBPT,,5952,1.70
TBPT,,7650,1.75
TBPT,,10105,1.80
TBPT,,13000,1.85
TBPT,,15905,1.90
TBPT,,26296,2.0
TBPT,,32901,2.05
TBPT,,42700,2.10
TBPT,,61700,2.15
TBPT,,84295,2.20
TBPT,,110000,2.25   
TB,BH,4,,30               定义永磁铁B-H曲线
tbpt,,355,.7
tbpt,,405,.8
tbpt,,470,.9
tbpt,,555,1.0   
tbpt,,673,1.1   
tbpt,,836,1.2   
tbpt,,1065,1.3  
tbpt,,1220,1.35
tbpt,,1420,1.4  
tbpt,,1720,1.45
tbpt,,2130,1.5  
tbpt,,2670,1.55
tbpt,,3480,1.6  
tbpt,,4500,1.65
tbpt,,5950,1.7  
tbpt,,7650,1.75
tbpt,,10100,1.8
tbpt,,13000,1.85
tbpt,,15900,1.9
tbpt,,21100,1.95
tbpt,,26300,2.0
tbpt,,32900,2.05
tbpt,,42700,2.1
tbpt,,61700,2.15
tbpt,,84300,2.2
tbpt,,110000,2.25   
tbpt,,135000,2.3
tbpt,,200000,2.41   
tbpt,,400000,2.69   
tbpt,,800000,3.22   
!*  
*CSET,1,2,       4,5,   
MPCOPY, ,4,5
TBCOPY,ALL,4,5  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MGXX,4,,0          定义永磁铁矫顽力
MPDATA,MGYY,4,,0
MPDATA,MGZZ,4,,HC   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MGXX,5,,0
MPDATA,MGYY,5,,0
MPDATA,MGZZ,5,,-HC  
BLOCK,-0.008,0.008,0.0125,0.0125+W,0,H,                       建模磁芯1
BLOCK,-0.008,0.008,-0.0125,-0.0125-W,0,H,                     建模磁芯2
BLOCK,-0.008,0.008,0.0125,0.0175,H,H+0.004,                   建模永磁铁1
BLOCK,-0.008,0.008,-0.0125,-0.0175,H,H+0.004,                 建模永磁铁2
BLOCK,-0.008,-0.004,-0.0125,0.0125,0,0.00015,                 建模镍带1
BLOCK,-0.002,0.002,-0.0125,0.0125,0,0.00015,                   建模镍带2
BLOCK,0.004,0.008,-0.0125,0.0125,0,0.00015,                    建模镍带3
BLOCK,-0.018,0.018,-0.028,0.028,-0.008,0.02,                   建空气
FLST,2,8,6,ORDE,2   
FITEM,2,1   
FITEM,2,-8  
VOVLAP,P51X
FLST,5,3,6,ORDE,3   
FITEM,5,6   
FITEM,5,9   
FITEM,5,-10
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       2, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ESIZE,0.0008,0,                  定义镍带单元大小
MSHAPE,1,3D
MSHKEY,0
!*  
FLST,5,3,6,ORDE,3   
FITEM,5,6   
FITEM,5,9   
FITEM,5,-10
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
/VIEW,1,1   
/ANG,1  
/REP,FAST   
FLST,5,2,6,ORDE,2   
FITEM,5,13  
FITEM,5,-14
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       3, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ESIZE,0,0   
ESIZE,0.001,0,                          定义磁芯单元大小
FLST,5,2,6,ORDE,2   
FITEM,5,13  
FITEM,5,-14
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
CM,_Y,VOLU  
VSEL, , , ,      11
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       4, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ESIZE,0.001,0,  
CM,_Y,VOLU  
VSEL, , , ,      11
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
CM,_Y,VOLU  
VSEL, , , ,      12
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       5, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,VOLU  
VSEL, , , ,      12
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
/VIEW,1,,,1
/ANG,1  
/REP,FAST   
CM,_Y,VOLU  
VSEL, , , ,      15
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ESIZE,0,0   
ESIZE,0.005,0,                           定义空气单元大小
CM,_Y,VOLU  
VSEL, , , ,      15
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FINISH  
/SOL
!*  
ANTYPE,0
MAGSOLV,2, , ,0.001,25,0
FINISH  
/POST1  
EPLOT   
ALLSEL,ALL  
FLST,5,3,6,ORDE,3   
FITEM,5,6   
FITEM,5,9   
FITEM,5,-10
VSEL,S, , ,P51X
ESLV,S  
ETABLE,EBY,B,Y                             提取镍带中磁感应强度BY
!*  
SSUM
!*  
*GET,BYTOT,SSUM, ,ITEM,EBY  
ALLSEL,ALL  
*SET,OBJECT,1/BYTOT
! LGWRITE,'core','lgw','C:\DOCUME~1\yajie\',COMMENT
FINISH  
/OPT
OPANL,'core','lgw',' '  
!*  
!*  
OPVAR,H,DV,0.001,0.01, ,                 定义设计变量H
!*  
OPVAR,W,DV,0.001,0.01, ,                  定义设计变量W
!*  
/OPT
OPSAVE,'corevar','opt',' '  
!*  
!*  
OPVAR,OBJECT,OBJ, , , ,                 定义目标函数
!*  
OPTYPE,SUBP                               零阶优化
OPSUBP,30,7,                                迭代次数
OPEQN,0,0,0,0,0,   
!*  
!*  
OPEXE   

yajie41

版主加各位高手啊，发表一下意见吧

yajie41

我刚才又把镍带厚度改成0.001mm，单元大小为0.001，算了一遍，结果还是会出现这样的Warning，这就证明应该不是镍带太薄单元划分不好导致的吧，那又是什么原因呢？？？急死了！

menyueqi

*** WARNING ***                        
Unknown parameter name= _Z3.  A value of 7.888609052E-31 will be used.
我认为是你这个参数没有定义好，你查查是否定义了此参数，或者此参数的获取值的过程中出现了定义错误。

lilulilu

我用你的命令算了几个循环，没有遇到你的问题。我的版本是9.0。
出了这个问题，可以用*stat检查参数设置。
根据我计算的结果，发现B出现负值，obj设置为B倒数最小，于是负值绝对值越大，反而结果越优，但这是否是你要的结果？
矫顽力的方向是否验证过？是不是你要的方向？
同意北极熊的看法，检查参数是否冲突。

xkg520

我运行了一遍，出错，我也是ansys9.0版本的
感觉你 在之前就没有定义优化文件，无法执行循环

rodman260

兄弟们这个问题怎么解决啊！我初学，今天也碰到这样的问题了，我做结构 分析的。

dahai1268

这个问题真是个难题啊！
我发现一个现象，我最近也在做优化，就是上面那个错误的值是一样啊！不管是_Z1，_Z2，_Z3都是7.888609052E-31，版主到底怎么回事啊？

shi120058

自己写个优化代码吧。省去了用ansys的优化程序的麻烦。