我是初學(xué)ansys ,麻煩大神指點(diǎn)!模型很簡(jiǎn)單,就是一個(gè)高為0.1m,直徑為1.44m的圓柱,圓柱的底面有9個(gè)半徑為0.01m小圓面,在這個(gè)9個(gè)小圓面施加全約束,圓柱只受重力。使這個(gè)圓柱的上表面的面積變形的rms。優(yōu)化設(shè)計(jì),設(shè)計(jì)變量就是這個(gè)9個(gè)小圓面的位置,按照?qǐng)A周均勻排列,里圈3個(gè),外圈6個(gè),里圈半徑為r1,外圈半徑為r2。r1,r2就是設(shè)計(jì)變量。目標(biāo)函數(shù)就是這個(gè)上表面積變形最少。單次運(yùn)行沒(méi)有警告啊能得出,但是一優(yōu)化就是循環(huán)一次就顯示不收斂結(jié)束了  ,我用的是ansys13.0 優(yōu)化設(shè)計(jì)模塊。14.0以上都沒(méi)有這個(gè)模塊了!金幣少,就先發(fā)10個(gè)了,好的話,我剩于金幣全給你了!下面是主要程序:
/PREP7
!*
ET,1,SOLID187
!*
*SET,radius1,0.3 !*設(shè)計(jì)變量
*SET,radius2,0.6 !*設(shè)計(jì)變量
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,9.1e10
MPDATA,PRXY,1,,0.24
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,2530
CYL4,0,0,0.72, , , ,0.1 !*鏡面模型
CYL4,0,radius1,0.01 !*下面為9個(gè)支撐約束面
CYL4,radius2*cos(1.0472),radius2*sin(1.0472),0.01
CYL4,-radius2*cos(1.0472),radius2*sin(1.0472),0.01
CYL4,-radius2*cos(1.0472),-radius2*sin(1.0472),0.01
CYL4,radius2*cos(1.0472),-radius2*sin(1.0472),0.01
CYL4,radius2,0,0.01
CYL4,-radius2,0,0.01
CYL4,-radius1*cos(0.5236),-radius1*sin(0.5236),0.01
CYL4,radius1*cos(0.5236),-radius1*sin(0.5236),0.01
asel,s,loc,z,0,0
asel,u,loc,x,0,0
asel,a,loc,y,radius1,radius1
asel,all
vsel,all
vsba,all,all
ESIZE,0.1,0, !*劃分網(wǎng)格
MSHAPE,1,3D
MSHKEY,0
!*
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FINISH
/SOLU
!*
ANTYPE,0
asel,none !*先選取9個(gè)面
asel,s,loc,z,0,0
asel,u,loc,x,0,0
asel,a,loc,y,radius1,radius1
/GO
DA,all,ALL, !*對(duì)9個(gè)面施加全約束
ACEL,0,0,9.8, !*只受重力
SOLVE
FINISH
/POST1
nsel,s,loc,z,0.1,0.1 !*選取上表面,并取上面變形的RMS
*get,nnub,node,,count
*get,nmin,node,,num,min
*dim,A,,nnub
*SET,n0,nmin
*SET,A(1),n0
*do,j,2,nnub
*SET,n0,ndnext(n0)
*SET,A(j),n0
*enddo
*dim,b,,nnub
*do,j,1,nnub
*SET,b(j),uz(a(j))
*enddo
*VSCFUN,zrms,RMS,b
FINISH
![ansys 優(yōu)化設(shè)計(jì)總是循環(huán)一次就結(jié)束,顯示不收斂]()
錯(cuò)誤.png |