风力机圆筒的优化设计
/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 (minimumradius of curvature of 1.59674861 and a maximum thickness of-4.6329818E-03 ).This violates the assumption of a shell element. 应该怎么修改谢谢！
页:
[1]