我在做一个轴的优化的时候,其中设计变量定义的是轴的第一段和第二段的半径和长度,状态变量定义的是模态分析的固有频率和绕度,目标变量定义的是轴的体积,是不是定义的有问题呀,还是命令流有问题?为什么求解不出来?还出现以下提示?
命令流如下:
*SET,theta,120
*SET,H1, 90
*SET,R1, 21
*SET,R2, 14
*SET,R3, 24
*SET,R4, 14
*SET,R5, 9.75
*SET,R6, 12
*SET,R7, 9.7
*SET,H2,128
*SET,H3,15
*SET,H4,427
*SET,H5,6
*SET,H6,38
*SET,H7,14
/PREP7
! Element types
et,1,SOLID92
! Material properties
mp,ex,1,2.059E11
mp,dens,1,7.8E3
mp,prxy,1,0.3
CYL4,,,R1,360,,,H1
wpoff,0,0,90
CYL4,,,R2,360,,,H2
wpoff,0,0,128
CYL4,,,R3,360,,,H3
wpoff,0,0,15
CYL4,,,R4,360,,,H4
wpoff,0,0,427
CYL4,,,R5,360,,,H5
wpoff,0,0,6
CYL4,,,R6,360,,,H6
wpoff,0,0,38
RPR4,6,0,0,r7,theta,h7
VADD,ALL
MSHAPE,1,3D
MSHKEY,0
SMRTSIZE,5,,,,,,,,,,,
VMESH,ALL
/SOL
DA,16,ALL
DA,15,ALL
SFA,4,1,PRES,1000
! /STATUS,SOLU
SOLVE
FINISH
/post1
ETABLE,EV,VOLU,
SSUM
*GET,VTOT,SSUM,,ITEM,EV
*GET,RHO,DENS,1,TEMP
WT=RHO*VTOT
U=ABS(U)
!
/SOLU
ANTYPE,2
MODOPT,SUBSP,5
EQSLV,FRONT
MXPAND,5,,,0
/SOLU
SOLVE
/POST1
SET,FIRST
*GET,F1,MODE,1,FREQ
nsort,u,sum
*get,umax,sort,,max
FINISH
! opt.inp
anfile='opt'
anext='inp'
/opt
opanl,anfile,anext
opvar,h1,dv,50,120
opvar,r1,dv,20,50
opvar,h2,dv,100,150
opvar,r2,dv,20,50
opvar,f1,sv,,100
opvar,umax,sv,,1.0
opvar,vtot,obj
optype,subp
opkeep,on
opsave,anfile,opt0
opexe |