! Végeslemes modellezés gyakorlat
! eromérő cella 2 - szabad háromszög háló üreg körül sűrítve
! [N,mm,MPa]

FINISH                  ! kilépés minden folyamatból
/CLEAR                  ! memória törlése

! munkakönyvtár beállítása:
/CWD, 'D:\Google_Drive\01_education\07_VEA\APDL_20211\eromero_cella_2'
! indítógomb létrehozása:
*ABBR,START,/input,eromero_cella_2,ans

! ADATBEVITEL !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Geometriai adatok
l=80        			! az erőmérő cella hossza
a=12        			! az erőmérő cella magassága
s=9         			! furatok középpontjainak távolsága
d=10        			! furatok átmérője
b=8         			! üreg kibővítése egy négyzettel(élhossz: 8mm)

! Anyagtulajdonságok
E=2e5 				    ! rugalmassági modulus 
nu=0.3 				    ! Poisson tényező

! Peremfeltételek
l0=20      			    ! megfogás hossza
l1=30       			! terhelés hossza
n=1 					! terhelési eset

!-------------------------------------------------------------------------------------------------
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! PREPROCESING !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!-------------------------------------------------------------------------------------------------

/PREP7                  ! belépés a preprocessorba (előkészítőbe)

! GEOMETRIA MEGRAJZOLÁSA !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

BLC4,0,0,l,a         	! téglalap rajzolása (bal alsó sarok, szélesség, magasság)
CM,cella,area          	! táglalap elnevezése "cella" komponensnek

CYL4,l/2-s/2,a/2,d/2   	! baloldali furat megadása (középpont, sugár)
CYL4,l/2+s/2,a/2,d/2   	! jobboldali furat megadása (középpont, sugár)
BLC5,l/2,a/2,b,b       	! téglalap üreg kibővítéséhez (középpont, szélesség, magasság)
CMSEL,u,cella          	! első téglalap (cella) elrejtése az aktuális elemek közül
CM,ureg,area           	! ami most látható (2 furat + négyzet) elnevezése "ureg"-nek

ALLSEL                	! minden kijelölése (újra minden látható lesz)
ASBA,cella,ureg			! az ureg kivonása a cellából

! Geometria előkészítése peremfeltételekhez
LSEL,s,loc,y,0          ! az alsó vonal kijelölése
LDIV,all,l0/l           ! a kijelölt vonal felosztása l0/l arányban -> 20 és 60 mm-es szakaszokra

LSEL,s,loc,y,a          ! a felső vonal kiválasztása
LDIV,all,l1/l           ! a kijelölt vonal felosztása l1/l arányban -> 30 és 50 mm-es szakaszokra

LSEL,r,loc,x,l-l1,l     ! a 30 mm-es szakasz kijelölése
LDIV,all,,,3            ! a kijelölt szakasz felosztása 3 egyenlő részre

! Geometria előkészítése hálózáshoz
LSEL,none               ! összes vonal elrejtése
*GET,p1,kp,0,num,max    ! legnagyobb sorszámú pont lekérése
K,p1+1,l/2-s/2-d/2-0.1*d/2,0  ! pontok definiálása a felület particionálásához
K,p1+2,l/2-s/2-d/2-0.1*d/2,a
K,p1+3,l/2+s/2+d/2+0.1*d/2,0
K,p1+4,l/2+s/2+d/2+0.1*d/2,a
L,p1+1,p1+2             ! vonalak definiálása a felület particionálásához
L,p1+3,p1+4

CM,cut,line             ! aktuális vonalak elnevezése "cut"-nak
ASBL,all,cut            ! cut kivonása az összes felületből (=particionálás)

ASEL,s,loc,x,l/2        ! x=l/2 középpontú felületek kijelölése (-> középső partíció)
CM,uregresz,area        ! kijelölt rész elnevezése "uregresz"-nek
ASEL,inve               ! kijelölés invertálása
CM,kulsoresz,area       ! kijelölés elnevezése "kulsoresz"-nek

ALLSEL                  ! minden kijelölése

! ANYAGTULAJDONSÁGOK !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MP,ex,1,E 				! rugalmassági modulus              
MP,nuxy,1,nu           	! Poisson tényező   

! VÉGESELEM HÁLÓ ELKÉSZÍTÉSE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ET,1,PLANE183           ! elemtípus megadása (1 -> sorszáma) (182:lineáris elem, 183: másodfokú elem)
KEYOPT,1,3,2            ! 1-es sorszámú elem 3.keyoptja: 3 (1:FSZ, 2:SA 3:ÁSF)
MSHAPE,1,2D				! elem alakja (0: négyszög, 1: háromszög)
MSHKEY,0                ! hálózási technika (0:szabad, 1:struktúrált, 2:vegyes)

! üregrész behálózása
ESIZE,0.1               ! elemméret 0,5 mm
AMESH,uregresz          ! uregresz behálózása

! külsőrész behálózása
!ESIZE,5                ! elemméret: 5 mm
! vagy:
ASEL,s,loc,x,l/2
ASEL,inve
LSLA,1
LESIZE,all,2,,,0.2


AMESH,kulsoresz		    ! külsőrész behalozasa

! PEREMFELTÉTELEK MEGADÁSA !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,s,loc,y,0          ! a befogás helyének kijelölése
LSEL,r,loc,x,0,l0
NSLL,s,1				! összes node kijelölése az aktív vonalon (1->végpontok is)
D,all,ux,0				! a kijelölt node-ok megfogása x irányban
D,all,uy,0				! a kijelölt node-ok megfogása y irányban

LSEL,s,loc,y,a			! a megadott terhelési szakasz kijelölése
LSEL,r,loc,x,l-l1+(n-1)*l1/3,l-l1+n*l1/3
SFL,all,pres,1       	! terhelés elhelyezése (p=1N/mm^2)

ALLSEL 					! minden kijelölése (újra minden látható lesz)

FINISH                  ! kilépés a preprocessorból

!-------------------------------------------------------------------------------------------------
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! SOLVING !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!-------------------------------------------------------------------------------------------------
/SOL					! belépés a solver-be
SOLVE					! számítás elvégzése
FINISH					! kilépés a solver-ből

!-------------------------------------------------------------------------------------------------
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! POSTPROCESSING !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!-------------------------------------------------------------------------------------------------

/POST1                  ! belépés a postprocessorba (kiértékelőbe)

PLNSOL,u,sum            ! elmozdulásmező ábrázolása
*ASK,next,OK,0          ! felugró ablak a következő ábrához
PLNSOL,epel,x,          ! epsilon_x fajlagos nyúlás szemléltetése
*ASK,next,OK,0          ! felugró ablak a következő ábrához
PLNSOL,s,eqv            ! redukált feszültségeloszlás ábrázolása

*GET,nn,node,0,count    ! csomópontszám lekérése