FINISH
/CLEAR,START

!!!!!!!!!!!!!!!!!!!!!!!Parametry (wymiary w milimetrach) !!!!!!!!!!!!!!!
dlugosc=5000
Rm=600				!!!Wytrzymałość stali
x=2				!!!Współczynnik bezp.
wielkosc_elementu=30

wysokosc_min=30
wysokosc_max=100
wysokosc_skok=30

szerokosc_min=30
szerokosc_max=100
szerokosc_skok=30

grubosc_min=10
grubosc_max=20
grubosc_skok=10


naprezenie_dop=(Rm/x)*1000000

gr_skok=grubosc_skok
szer_skok=szerokosc_skok
wys_skok=wysokosc_skok


!!!!!!*DO,grubosc,grubosc_max,grubosc_min,-gr_skok

!!!!!!!!*DO,szerokosc,szerokosc_max,szerokosc_min,-szer_skok

!!!!!!!*DO,wysokosc,wysokosc_max,wysokosc_min,-wys_skok


*DO,grubosc,grubosc_min,grubosc_max,gr_skok
*DO,szerokosc,szerokosc_min,szerokosc_max,szer_skok
*DO,wysokosc,wysokosc_min,wysokosc_max,wys_skok



!!!!!!!!!! PARAMETRYCZNY MODEL GEOMETRYCZNY !!!!!!!!

wys=wysokosc*0.001
szer=szerokosc*0.001
grub=grubosc*0.001
dlug=dlugosc*0.001

wys=wys-grub
szer=szer-(grub/2)

!!*IF,grub*2,GT,wys,CYCLE		!!!! Jeżeli podwojona grubość ścinaki jest większa od wysokości ceownika, to przestaje one być ceownikiem, a zamienia się w pełną belkę
!!*IF,grub,GT,szer*2,CYCLE	!!! Podobnie j.w.


/PREP7

K, 1, 0, 0, wys/2
K, 2, szer, 0, wys/2
K, 3, 0, 0, -wys/2
K, 4, szer, 0, -wys/2

K, 5, 0,dlug, wys/2
K, 6, szer,dlug, wys/2
K, 7, 0,dlug, -wys/2
K, 8, szer,dlug, -wys/2


L,1,2
L,1,3
L,3,4

L,5,6
L,5,7
L,7,8

L,3,7
L,4,8
L,1,5
L,2,6

AL, 1,9,4,10
AL, 3,7,6,8
AL, 7,5,9,2




!!!!!!!!!!!!! PODZIAŁ NA ELEMENTY !!!!!!!!!!!!
ET,  1,SHELL93  
R,1,grub,grub,grub,grub
MP,EX,1,2.1E11
LESIZE,ALL,wielkosc_elementu*0.001 !!!! przerobić na żadsze
!LESIZE,ALL,,,10
!ESHAPE,1
AMESH,ALL
/SOLU


!!!!!!!!!!!!!! Utwierdzenie !!!!!!!!!!!!!!!!
NSEL,S,LOC,Y,0
D,ALL,ALL
ALLSEL           


!!!!!!!!!!!!!!! Obciążenie (tylko od siły ciężkości) !!!!!!!!!!!
/PREP7
MP,DENS,1,7.81e3  ! Określenie masy(gęstości)
ACEL,0,0,9.81	! Określenie grawitacji
/SOLU


!!!!!!!!!!!!!!!! Duże przemieszczenia !!!!!!!!!!!!!!!  
!ANTYPE,0            ! analiza statyczna
!NLGEOM,ON
!KBC,0               ! liniowe narastanie warunków brzegowych
!OUTPR,ALL,ALL   
!OUTRES,ALL,ALL
!AUTOTS,ON
!NSUBST,5,6,5
!!NSUBST,5,50,5       ! podział obliczeń na podkroki obliczeniowe (początek, max, min)
!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!Tworzenie zmiennej "Masa" !!!!!!!!!!!!
/SOL
IRLF,-1                    !* Precalculation of mass
SOLVE		!!!!!! Zostawić SOLVE!
*GET,Masa,ELEM,0,MTOT,X    !* liefert die Masse aller Elemente
IRLF,0                      !* No inertia relief calculations


!!!!!!!!!!!!!!!!! Tworzenie zmiennej "vonMises" (maksymalna wartość naprężeń) !!!!!!
/POST1 
set,last          

NSORT,S,INT        
*get,vonMises,sort,,max        


!!!!!!!!!!!!!!!!!Zapis zmiennych do pliku "dane.txt"!!!!!!!!!!!!!!!
*CFOPEN, dane, txt,,APPEND
*CFWRITE,ciezar, Masa, gr:, grubosc,sze:, szerokosc,wy:, wysokosc, napr:, vonMises

!!!!!!!!!!!!!!!!Czyszczenie geometrii!!!!!!!!!!!!!!
/PREP7
ACLEAR,ALL
ADELETE,ALL,,,1



!!!!!!!!*IF,vonMises,GT,naprezenie_dop,EXIT
	

*ENDDO
!!!!!!!!*IF,vonMises,GT,naprezenie_dop*50,THEN
!!!!!!!!wys_skok=wysokosc_skok*10
!!!!!!!!*ELSEIF
!!!!!!!!wys_skok=wysokosc_skok
!!!!!!!!*ENDIF
*ENDDO
!!!!!!!!*IF,vonMises,GT,naprezenie_dop*50,THEN
!!!!!!!!gr_skok=grubosc_skok*10
!!!!!!!!*ELSEIF
!!!!!!!!gr_skok=wysokosc_skok
!!!!!!!!*ENDIF
*ENDDO


FINISH