SimWe仿真论坛 » A15：iSIGHT » 请教一个关于Gradient Block 的问题

rabitzh 发表于 2008-6-16 19:27

请教一个关于Gradient Block 的问题

[b]请教一个关于Gradient Block 的问题
c6CuzB%u&ci"^*} [/b]最近做了一个课题，需要利用iSIGHT+Nastran进行优化，本来是iSIGHT带有Nastran的集成模块的，但里面存在些问题，再者，由于结构设计优化的对象不再是传统的结构重量，而是与结构重量相关的显示函数。
NSQ6c3Q)|1uha 一个解决方法就是利用iSIGHT内部的自动差分方法计算结构的敏度信息，由于数值噪声以及数值传递过程中存在误差，因此这种方法是应该极力避免的。
h9ibrV Nastran软件可以输出敏度信息，因此为我们优化提供了方便。
Mb4[m
R!`7e 现假设存在一个目标参数——[color=red]TotalMass[/color]
U5a-qy,FG n2hh'_ 一个约束参数——[color=orange]Feq1[/color]
K0R^IIjG-FO 敏度信息通过文件解析提取于数组中，数组参数为——[color=red]OBJGrad[2][24][/color] t,p_)M1G ltP
现在添加优化梯度信息模块：8r\/E8[i W#v!s@ [u
Gradient
X
LtOE _0leJ             GradProcess ObjFuncGrad
8zUR| a4Ua'q                 Outputs
yEuR&I{                     TotalMass Feq1
U#l+~i7}                 Inputs
p\&u,|"Z.KY                     PLTT PLBT PFTT PFMT PFBT SCT PLRT TTRW TTRH TMRWR+N:W,|3lFZ
                    TMRH TBRW TBRH THW THH TVW TVH TAW TAH PFBW4])Gu@ Y&RZW2_ TRW
                    PFBH PFTW PFTH AT
a|lb g fo/D9t                 Control: [
3^F+HA AJ6k[N                     TclR%ey+^7J'nzQR
                        esi_GradientObject set -array OBJGrad
i8eL$y+[@nHn@%S                         puts "Gradient block running###"'w{)~eC!q\4t
                    End Tcl
@$YZL,aZ.F?9q                 ]Q)E!l5Jbh*@E*QI
            End GradProcess ObjFuncGradv`X*|)|
            Outputs
.\?/x7A"j8Rc                 TotalMass Feq1H.Dr:H_,@
            Inputsp3Sk,_z{
                PLTT PLBT PFTT PFMT PFBT SCT PLRT TTRW TTRH TMRW
zDA3c(n&I                 TMRH TBRW TBRH THW THH TVW TVH TAW TAH PFBW
t5HM)DzE0q                 PFBH PFTW PFTH AT
^&hi9Gw             Control: [%YOC/F&G,m:w
                ObjFuncGrad;W&y(lz??8n9e i
            ]U{Byk/rB)y$l
        End GradientP'rv$R O
存在的问题，优化迭代一次就完成，显然是有错误存在。
2F w8rb&\-\]G @1aX 仔细查看其中的log文件，里面提示其中的错误：
z%x3l_E(~p!i      CONSTRAINED FUNCTION MINIMIZATION D2m1g(aGe8jRX;Z
     CONTROL PARAMETERS
['V&io7F      IPRINT  NDV    ITMAX    NCON    NSIDE  ICNDIR   NSCAL   NFDGRs%h+n?
       2      24      20       2       1      25       0       1
U!W'}iN2n+R7FO      LINOBJ  ITRM     N1      N2      N3      N4      N52Jz
e,ryB4b
       0       3      26      50       7      24      480x
Q |]
M.c0~a0F
         CT              CTMIN           CTL             CTLMINzblqa5fC+v@
     -0.10000E-03     0.40000E-03    -0.10000E-01     0.10000E-02
;G#Sm,b9] {          THETA           PHI             DELFUN          DABFUN
@[cj {h{/KK       0.10000E+01     0.50000E+01     0.10000E-02     0.10000E-021Xq-n%_#p'qm5W
         FDCH            FDCHM           ALPHAX          ABOBJ1

`!^ a*l@H3C9~       0.10000E-01     0.10000E-01     0.10000E+00     0.10000E+00kL bn3T$~&S'e2cN
     LOWER BOUNDS ON DECISION VARIABLES (VLB)\.QF eb2qOw
       1)    0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01
"C2g&g4XP;E`j        7)    0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-017D!wb:k+Lw7p6t4Ye
      13)    0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01
+XLA$uT\x       19)    0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01  0.10000E-01
4F'@`Uk|Ohx3B      UPPER BOUNDS ON DECISION VARIABLES (VUB)
;IQ,OY9C'W        1)    0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00
"Qp.Y$l8?r        7)    0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00)Q/h C1q9w:p(`5sM
      13)    0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00
![ ud
AV.p%OJ       19)    0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00?,_b*?w.g5a
     ALL CONSTRAINTS ARE NON-LINEAR+Nm.^,l(zb9b"z

w6V)_T5k      INITIAL FUNCTION INFORMATIONs!tT:~#z]0I
     OBJ =   0.110460E+04U[M7Z,A2\}@
     DECISION VARIABLES (X-VECTOR)
s4]!H4rx`&l(Vm        1)    0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00~`[%^#}
       7)    0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00
rPK!y@LH/{ T/e       13)    0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00
B4p B&q1O)j9V       19)    0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00  0.24000E+00^`#_%mY~4O*j
     CONSTRAINT VALUES (G-VECTOR)
'R(g%X!K#G        1)   -0.10526E+02 -0.10000E+01

A?3HcwK+a
WgnON;C$f)w      THE NUMBER OF ACTIVE AND VIOLATED CONSTRAINTS EXCEEDS N3-1..GD!S&h`#Y
     DIMENSIONED SIZE OF MATRICES A AND B AND VECTOR IC IS INSUFFICIENTC E5_c3[
     OPTIMIZATION TERMINATED AND CONTROL RETURNED TO MAIN PROGRAM.
y"MC&I4C PATH: /tm/Structure TYPE: Status-Descfile FipFopStatus is OKY{7DwXrmS,e
'O"tB0g,{
P3M
[color=mediumturquoise]PATH: /tm/Structure TYPE: Error-Optimization API Error : Unable to access gradient repository2V&\Z[Mp
PATH: /tm/Structure TYPE: Error-Technique Bad Optimization Run - Could not get objective gradients in CONMIN.[/color]Y5W7KM-t:G
对于提示的结果，我不是很明白，好像说是约束数量超了，这是怎么回事？A、B的矩阵维数不对？9lXXgx}*`
请问，上述问题如何解决呢？
9l#QHP|2ea 另外，esi_GradientObject set (-list{<real>...<real>}|-var<name>|-array <name> |-val <real>)+yN)zYGH
这里面是怎么应用的啊？看了帮助文档，还是不明白，有没有这方面的实例啊？
I)c2T0akh#Xr&IT 版上有没有人遇到过类似的问题呢？麻烦交流交流。

查看完整版本: 请教一个关于Gradient Block 的问题