优化设计实例

第26例优化设计实例一梁的优化设计［本例提示］介绍了优化设计的相关理论和应用，讲述了将设计问题的物理模型转化为数学模型一选取设计变量、写出目标函数、给出约束条件的方法。.优化分析文件!File:EXAMPLE26.LGWH=0.06B=0.06L=1/PREP7ET,1,BEAM3R,1,B*H,B*H*H*H/12,HMP,EX,1,2E11MP,NUXY,1,0.3K,1,0,0,0K,2,L,0,0LSTR,1,2LESIZE,1,,,50LMESH,1FINISH/SOLUDK,1,UXDK,1,UYDK,2,UXDK,2,UYSFBEAM,ALL,1,PRES,5000SOLVESAVEFINISH/POST1ETABLE,E_VOL,VOLUSSUM*GET,V_TOT,SSUM,,ITEM,E_VOLNSORT,U,Y*GET,UY_MIN,SORT,,MINUY_MAX=ABS(UY_MIN)FINISH.优化控制文件!File:EXAMPLE26_OPT.TXT/FILNAME,EXAMPLE26/OPTOPCLRFINISH/CLEAR/INPUT,EXAMPLE26,LGW/OPTOPANL,EXAMPLE26,LGWOPVAR,B,DV,0.05,0.1,0.001OPVAR,H,DV,0.05,0.1,0.001OPVAR,UY_MAX,SV,0,2E-4,1E-5OPVAR,V_TOT,OBJ,,,1E-5OPSAVE,EXAMPLE26,OPTOPTYPE,FIRSTOPFRST,30OPEXECOPLIST,ALLFINISH.求解方法先分另1J以文件名EXAMPLE26.LGW、EXAMPLE26_OPT.TXT将优化分析文件和优化控制文件存储在ANSYS的工作文件夹里，再在ANSYS的输入窗口输入/INPUT,EXAMPLE26_OPT,TXT,回车，开始求解。第27例优化设计实例一曲柄摇杆机构的优化设计本例提示］介绍了一个更复杂目标中野鼓用点的轨迹与目标位置点曲线的对比情况1.优化分析文件、更接近于实际的优化设计实例X璃方向川九T!File:EXAMPLE27.LGWPI=3.1415926A=2B=8C=7.5D=9.5AX=1AY=4BATA_OPT=PI/16K=7GAMA=PI/3!建立有限元模型/PREP7BATA=BATA_OPT-PI/2OMGA1=0.5T=60/OMGA1ET,1,COMBIN7ET,2,BEAM4MP,EX,1,2E11MP,PRXY,1,0.3MP,DENS,1,1E-14R,1,1E9,1E3,1E3,0R,2,4E-4,1.3333E-8,1.3333E-8,0.02,0.02COS_SETA1=((A+B)*(A+B)+D*D-C*C)/2/(A+B)/D*IF,COS_SETA1,GE,-1,AND,COS_SETA1,LE,1,THENSETA1=ACOS(COS_SETA1)*ELSESETA1=0*ENDIFBX=AX+A*COS(SETA1+BATA)BY=AY+A*SIN(SETA1+BATA)CX=BX+B*COS(SETA1+BATA)CY=BY+B*SIN(SETA1+BATA)DX=AX+D*COS(BATA)DY=AY+D*SIN(BATA)KX=BX+K*COS(BATA+SETA1+GAMA)KY=BY+K*SIN(BATA+SETA1+GAMA)N,1,AX,AYN,2,BX,BYN,3,BX,BYN,4,CX,CYN,5,CX,CYN,6,DX,DYN,7,BX,BY,-1N,8,CX,CY,-1N,9,KX,KYTYPE,1REAL,1E,2,3,7E,4,5,8TYPE,2REAL,2E,1,2E,3,4E,5,6E,3,9E,4,9FINISH!求解/SOLUANTYPE,TRANSNLGEOM,ONDELTIM,T/100KBC,0TIME,TOUTRES,BASIC,ALLAUTOTS,ONCNVTOL,F,1,0.1CNVTOL,M,1,0.1D,ALL,UZD,ALL,ROTXD,ALL,ROTYD,1,ROTZ,2*PI-0.001D,1,UXD,1,UYD,6,UXD,6,UYSOLVESAVEFINISH!提取数据/POST26*DIM,X_OBJ,ARRAY,10*DIM,Y_OBJ,ARRAY,10X_OBJ(1)=9.5Y_OBJ(1)=8.26X_OBJ(2)=9Y_OBJ(2)=8.87X_OBJ(3)=7.97Y_OBJ(3)=9.51X_OBJ(4)=5.65Y_OBJ(4)=9.94X_OBJ(5)=4.36Y_OBJ(5)=9.7X_OBJ(6)=3.24Y_OBJ(6)=9X_OBJ⑺=3.26Y_OBJ⑺=8.36X_OBJ(8)=4.79Y_OBJ(8)=8.11Y_OBJ(9)=8X_OBJ(10)=9.12Y_OBJ(10)=7.89X_OBJ_MIN=3.24X_OBJ_MAX=9.5NSOL,2,9,UXNSOL,3,9,UY*DIM,X9,ARRAY,100*DIM,Y9,ARRAY,100VGET,X9,2VGET,丫9,3*DO,I,1,100X9(I)=KX+X9(I)Y9(I)=KY+Y9(I)*ENDDOG2=B+C-A-DG3=C+D-A-BG4=B+D-A-CCOS_GAMA1=(B*B+C*C-(D-A)*(D-A))/2/B/CCOS_GAMA2=(B*B+C*C-(D+A)*(D+A))/2/B/C*IF,COS_GAMA1,GE,-1,AND,COS_GAMA1,LE,1,THENG5=ACOS(COS_GAMA1)-30*PI/180*ELSEG5=-1*ENDIF*IF,COS_GAMA2,GE,-1,AND,COS_GAMA2,LE,1,THENG6=150*PI/180-ACOS(COS_GAMA2)*ELSE*ENDIFF2=0*DO,I,1,10F2=F2+(Y9(I*10)-Y_OBJ(I))*(Y9(I*10)-Y_OBJ(I))+(X9(I*10)-X_OBJ(I))*(X9(I*10)-X_OBJ(I))*ENDDOFINISH.优化控制文件!File:EXAMPLE27_OPT.TXT/OPTOPCLRFINISH!/CLEAR/INPUT,EXAMPLE27,LGW/OPTOPANL,EXAMPLE27,LGWOPVAR,A,DV,0.1,10,0.001OPVAR,B,DV,0.1,10,0.001OPVAR,C,DV,0.1,10,0.001OPVAR,D,DV,0.1,10,0.001OPVAR,AX,DV,0.1,5,0.001OPVAR,AY,DV,0.1,5,0.001OPVAR,BATA_OPT,DV,0.1,PI,0.001OPVAR,K,DV,0,10,0.001OPVAR,GAMA,DV,0.1,PI,0.001OPVAR,G2,SV,0,100OPVAR,G3,SV,0,100OPVAR,G4,SV,0,100OPVAR,G5,SV,0,100OPVAR,G6,SV,0,100OPSAVE,EXAMPL