/CLEAR,NOSTART  

/TITLE, Optymalizacja Topologiczna 2-D 14.05.2012
/PREP7
a=5
b=0.1*a
delta=(a-b)/3

F1=-5000
F4=-1000

a_=(F2-F1)/(a-b)
b_=F1-a_*b

X2=b+delta
F2= a_*X2+b_

X3=b+2*delta
F3=a_*X3+b_

BLC4,0,0,a,1          ! Model belki o wymiarach ax1 prostokąt
  
ET,1,82               ! Elementy wykorzystane do analizy SOLI. 
                      ! Nr 1 do optymalizacji
ET,2,82               ! Nr 2 nie do optymalizacji 

MP,EX,1,118E9         ! Definicja Materiału liniowy izotropowy
MP,NUXY,1,0.3

ESIZE,a/100           ! Definicja wielkości elementu
TYPE,1                ! Wybór pierwszego elementu
AMESH,ALL             ! Tworzenie siatki FREEMESH
NSEL,S,LOC,X,0,b      ! Selekcja obszru nie do otymalizacji
ESLN
TYPE,2                ! Wybór drugiego elementu
EMODIF,ALL            ! Modyfikacja Siatki
ALLSEL     

NSEL,S,LOC,X,0        ! Selekcja węzłów do definicji utwierdzenia
D,ALL,ALL,0           ! Definicja utwierdzenia

                      !Definicja Obciążeń
NSEL,S,LOC,X,b
NSEL,R,LOC,Y,0
F,ALL,FY,F1       
ALLSEL

NSEL,S,LOC,X,X2
NSEL,R,LOC,Y,0
F,ALL,FY,F2      
ALLSEL

NSEL,S,LOC,X,X3
NSEL,R,LOC,Y,0
F,ALL,FY,F3       
ALLSEL

NSEL,S,NODE,,2
NSEL,R,LOC,Y,0
F,ALL,FY,F4      
ALLSEL


LSWRITE,1             ! Zapis konfiguracji obciążenia
FDEL,ALL

TOCOMP,MCOMP,SINGLE,1  ! Definicja Optymalizacji 
                       ! "MCOMP" Optymalizacja Topologiczna 
TOVAR,MCOMP,OBJ        ! Definicja "MCOMP" jako funkcja Celu 
TOVAR,VOLUME,CON,,50   ! Definicja ograniczeń 50% redukcja objętości
                           
TOTYPE,OC              ! Definicja solvera
TODEF                  ! Rozpoczęcie optymalizacji

/DSCALE,,OFF
/CONTOUR,,2
TOLOOP,100,1           ! Maksymalna liczba iteracji
FINISH


/POST1  
PLESOL, TOPO,, 0,1.0

TOPLOT,0