Finish/Clear
/TITLE,2-D Two-Storyreinforcement problem - Maximize fundamental frequency
A=0.25 !Prepare ModelB=5C=0.375D=8E=3.75
/PREP7K,1K,2,CK,3,C+AK,4,B/2K,5,,EK,6,C,EK,7,C+A,EK,8,B/2,EKSEL,S,,,5,8KGEN,2,ALL,,,,AKSEL,S,,,9,12
KGEN,2,ALL,,,,EKSEL,S,,,13,16KGEN,2,ALL,,,,AALLSELA,1,2,6,5A,5,6,10,9A,9,10,14,13A,13,14,18,17A,2,3,7,6A,6,7,11,10A,10,11,15,14A,14,15,19,18A,3,4,8,7A,7,8,12,11A,11,12,16,15A,15,16,20,19
ET,1,82 ! Define two element type regionsET,2,82 !1 - optimized regionASEL,S,,,4,8 !2 - non-optimized regionASEL,A,,,10,12,2TYPE,2ESIZE,0.1AMESH,ALLASEL,INVETYPE,1AMESH,ALLALLSELMP,EX,1,100 ! Material of structureMP,NUXY,1,0.3MP,DENS,1,1.0E-6MP,EX,2,100 ! Material of concentrated massesMP,NUXY,2,0.3MP,DENS,2,5.0E-6
ASEL,S,,,6,8,2ESlA,S,1EMODIF,ALL,MAT,2 ! Define concentrated massesALLSELLOCAL,11,0,2.5ARSYM,X,ALL ! Full modelNUMM,KPNUMM,ELEMNUMM,NODE
LSEL,S,,,14LSEL,A,,,45NSLL,S,1D,ALL,ALLALLSEL
TOFREQ,FREQ1,SING,1 !Define single frequency function (1st)TOVAR,FREQ1,OBJ ! Define objective for topological optimizationTOVAR,VOLUME,CON,,14,ACTUAL ! Define volume constraint (upper bound = 14)TOTYPE,SCP ! Select ScP solution approach !TODEF,0.00001 !Initialize topological optimization process,accuracy = 0.00001TOLOOP,40,1 !Perform up to 40 iterationsFINISHTOPLOT,1 !Plot final pseudo-densitiesTOGRAPH,OBJ !Graph final objective (1st frequency) historyTOGRAPH,CON !Graph final constraint (volume) history
