Topology optimization in Ansys Classic (lathe constraint) скачать в хорошем качестве

Topology optimization in Ansys Classic (lathe constraint)

Simple topology optimization written in APDL. During the optimization a lathe constraint has been applied (for other manufacturing constraints see my other videos). Achieved strain energy: 0.1685 Joule. For those who want to compare my results with other optimizer software, I provide the APDL code of the meshed model (with applied BCs and loads): finish /clear,start finish /prep7 !************************************** ! INPUT * !************************************** d1=60 d2=300 width=120 bushing_thickness=10 thickness=8 env_size=450 elem_size=5 elem_division=11 cylinder_lug_width=50 component_width=300 !************************************** ! END INPUT * !************************************** wpro,,,-90 csys,4 k,1,0,0,0 k,2,width,0,0 blc4,0-width,d1/2,width,env_size/2-d1/2 wpro,,,90 wpoff,0,0,0-thickness wpoff,0,d2/2,0 asbw,all wpro,,90 asbw,all wpoff,0,0,(d2-d1)/2-bushing_thickness asbw,all adele,6,,,1 csys,0 dsys,0 wpcsys,-1,0 wpro,,,-90 csys,4 et,1,200,6,0 et,2,185 et,3,200,0,0 et,4,185 et,5,4 et,6,188 et,7,170 et,8,174 keyopt,8,2,0 keyopt,8,4,0 keyopt,8,5,0 keyopt,8,7,0 keyopt,8,8,0 keyopt,8,9,1 keyopt,8,10,2 keyopt,8,11,0 keyopt,8,12,0 keyopt,7,5,0 r,5,99999999,99999999,99999999,1,1 r,6,,,1.0,0.1,0, rmore,,,1.0e20,0.0,1.0, rmore,0.0,0,1.0,,1.0,0.5 rmore,0,1.0,1.0,0.0,,1.0 mptemp,1,0 mpdata,dens,1,,7.8e-9 mpdata,ex,1,,210000 mpdata,prxy,1,,0.3 mptemp,1,0 mpdata,ex,5,,210000 mpdata,prxy,5,,0.3 mptemp,1,0 mpdata,dens,6,,7.8e-9 mpdata,ex,6,,210000 mpdata,prxy,6,,0.3 mp,mu,6,0.1 sectype,6,beam,csolid,,0 secoffset,cent secdata,20,0,0,0,0,0,0,0,0,0,0,0 aatt,1,,1 aesize,all,elem_size mshape,0,2d mshkey,1 amesh,all type,2 extopt,esize,elem_division,0 extopt,aclear,0 extopt,attr,1,0,0 real,1 esys,0 vrotat,all,,,,,,1,2,360,4 kdele,1,2,1 csys,0 dsys,0 wpcsys,-1,0 ! ! Find design space ! vsel,s,,,5,20,5 allsel,below,volu esel,r,type,,2 cm,topo_opt,elem allsel,below,volu esel,r,type,,1 nsle cm,tool_face,elem lsel,s,,,14 cm,tool_top,line lsel,s,,,20 cm,tool_left,line allsel,below,volu ! ! Make pin ! csys,0 dsys,0 wpcsys,-1,0 wpoff,0,0,elem_size*2 vsel,none cyl4,0,0,d1/2,,,,0-width-cylinder_lug_width/2-elem_size*2 cm,pin_solid,volu vatt,1,,4 esize,elem_size,0 vsweep,all type,5 real,5 mat,5 wpoff,0,0,elem_size*2 csys,4 *get,ar23,node,0,num,maxd n,ar23+1,0,0,0 asel,s,,,82 allsel,below,area nsla,s,1 *get,ar24,node,0,count *del,kortombtmp *dim,kortombtmp,array,ar23,2 *vfill,kortombtmp(1,1),ramp,1,1 *vget,kortombtmp(1,2),node,0,nsel *del,kortomb *dim,kortomb,array,ar24 *vmask,kortombtmp(1,2) *vfun,kortomb(1),comp,kortombtmp(1,1) nsel,a,,,ar23+1 e,ar23+1,kortomb(1:ar24) csys,0 dsys,0 wpcsys,-1,0 wpoff,0,0,component_width/2 csys,4 n,ar23+2,0,0,0 type,6 real,1 mat,1 secnum,6 e,ar23+1,ar23+2 nsel,s,,,ar23+2 cm,pin_middle,node d,all,uz d,all,rotx d,all,roty d,all,rotz asel,s,,,81 allsel,below,area nsla,s,1 cm,loaded,node csys,0 dsys,0 wpcsys,-1,0 allsel,all ! ! Define contact ! csys,0 dsys,0 wpcsys,-1,0 cswpla,21,1,1,1 vsel,all cmsel,u,pin_solid allsel,below,volu ksel,r,loc,x,d1/2-bushing_thickness/2,d1/2+bushing_thickness/2 lslk,s,1 asll,s,1 allsel,below,area nsla,s,1 cm,kon_targ,node asel,s,,,83,84,1 allsel,below,area nsla,s,1 cm,kon_cont,node allsel,all ksel,s,loc,x,env_size/2-elem_size,env_size/2+elem_size lslk,s,1 asll,s,1 allsel,below,area nsla,s,1 cm,bc,node type,7 real,6 mat,6 cmsel,s,kon_targ esln,s,0 esurf,,top cmsel,s,kon_cont type,8 esln,s,0 esurf,,top csys,0 dsys,0 wpcsys,-1,0 allsel,all