风力机圆筒的优化设计

lvgang637

/units,si
/FILNAM,SCRATCH
/TITLE,OPTIMIZATION OF TOWER
/PREP7
*SET,D1,6
*SET,D2,5
*SET,D3,5
*SET,D4,4.5
*SET,D5,4.5
*SET,t1,0.025
*SET,t2,0.025
*SET,t3,0.020
*SET,t4,0.020
*SET,t5,0.015
*SET,P,1.29                 !定义空气的密度
*SET,R,20                   !定义风轮半径
*SET,V,10                   !定义额定风速
*SET,V1,11                  !定义风轮顶端风速
*SET,V2,9.5                 !定义风轮底端风速
*SET,S_XUYONG,2.35e8        !定义许用应力
*SET,value,0.75             !定义挠度           
MP,EX,1,2.07e11
MP,PRXY,1,0.3   
MP,DENS,1,7850  
ET,1,SHELL93
K,1,D1/2,,,  
K,2,D2/2,,20.5,  
K,3,D3/2,,41,
K,4,D4/2,,61.5,  
K,5,D5/2,,82,
l,1,2   
l,2,3   
l,3,4   
l,4,5   
LPLOT
K,100,0,,0,     
K,101,0,,82,   
Arotat,all,,,,,,100,101
ASEL,S,LOC,Z,0,20.5        !第一节
R,1,t1
AMESH,ALL
EPLOT
ASEL,S,LOC,Z,20.5,41       !第二节
AMESH,ALL
MZNODE=NDINQR(0,14)
*DIM,THICK,,MZNODE
*DO,NODE,1,MZNODE
  *IF,NDINQR(NODE,1),EQ,1,THEN
     THICK(NODE)=t1-(t1-t2)/20.5*NZ(NODE)
   *ELSE
      THICK(NODE)=0
   *ENDIF
*ENDDO
NODE=
MZNODE=
RTHICK,THICK(1),1,2,3,4
ASEL,S,LOC,Z,41,61.5      !第三节
R,2,t3
AMESH,ALL
EPLOT
ASEL,S,LOC,Z,61.5,82       !第四节
AMESH,ALL
MZNODE2=NDINQR(0,14)
*DIM,THICK2,,MZNODE2
*DO,NODE,1,MZNODE2
  *IF,NDINQR(NODE,1),EQ,1,THEN
     THICK2(NODE)=t4-(t4-t5)/20.5*NZ(NODE)
   *ELSE
      THICK2(NODE)=0
   *ENDIF
*ENDDO
NODE=
MZNODE2=
RTHICK,THICK2(1),1,2,3,4
EPLOT
ALLSEL,ALL
CSYS,0               
K,9999,0,0,82
ET,2,MASS21
KSEL,S,LOC,Z,78,82
KMESH,ALL   
NSEL,S,LOC,Z,80,82  
CERIG,NODE(0,0,82),ALL,ALL,,,,!刚性区域的划分
/SOL
ANTYPE,0
LSEL,S,LOC,Z,0,0
DL,ALL,,ALL
ACEL,0,0,9.8,
nsel,s,loc,x,0
nsel,r,loc,y,0
nsel,r,loc,z,82
F,ALL,FY,1/8*4/9*P*3.14*(R**2)*(V**2)       !轴向推力
F,ALL,FY,1/2*1/9*P/3*3.14*(R**2)*(V1-V2)**2 !脉动力
F,ALL,MX,1/2*8/27*P/3*3.14*(R**3)*(V1**2-V2**2)!俯仰力矩
nsel,all
SOLVE                  
FINISH
/POST1
PLNSOL, U,SUM,0,1   
PLNSOL, S,EQV,2,1
/POST1
ETABLE,VOLU,VOLU
ETABLE,SMAX_I,NMISC,1
ETABLE,SMAX_J,NMISC,3
ESORT,ETAB,SMAX_I,,1
*GET,SMAXI,SORT,,MAX
ESORT,ETAB,SMAX_J,,1
*GET,SMAXJ,SORT,,MAX
SMAX=ABS(SMAXI)>ABS(SMAXJ)
SSUM
*GET,VOLUME,SSUM,,ITEM,VOLU
*STATUS,Parm             ! 列表显示内存变量的值
FINISH                   ! 退出后处理模块
LGWRITE,SCRATCH,lgw      ! 宏程序定义结束
/OPT
OPANL,SCRATCH,lgw        ! 定义优化设计文件SCRATCH
OPVAR,t1,DV,0.02,0.04
OPVAR,t2,DV,0.02,0.04
OPVAR,t3,DV,0.01,0.03
OPVAR,t4,DV,0.01,0.03
OPVAR,t5,DV,0.01,0.03   
OPVAR,SMAX,SV,0,S_XUYONG !定义应力状态变量（约束）
opvar,DMAX,sv,,value     !定义位移状态变量
OPVAR,VOLUME,OBJ         !定义目标函数
OPTYPE,SUBP              !使用子模型近似优化设计方法
OPSUBP,15                ! 最大优化设计15个循环      
OPPRNT,ON                ! 显示优化过程中的详细情况
OPEXE                    ! 用子模型近似优化设计方法开始优化设计
OPLIST,ALL,,1            ! 显示设计集中的参数
/AXLAB,Y,VOLUME          ! 将目标函数设置为Y轴
PLVAROPT,VOLUME          ! 显示体积收敛过程
*STATUS,Parm             ! 列表显示内存变量的值
请问：我在进行优化循环的过程中，为什么出现了这个错误Element 6003 has a radius/thickness ratio of -344.648151 (minimum  radius of curvature of 1.59674861 and a maximum thickness of  -4.6329818E-03 ).  This violates the assumption of a shell element.   应该怎么修改谢谢！