fini
/clear,start

filename='Doppel-Elem-Tetra-4'
/filnam,filename
/title,%filename%

*abbr,Doppel-Elem-Tetra-4,Doppel-Elem-Tetra-4.txt


/prep7
/vup,,z
/triad,rbot 								! turn off the global triad display
/pbc,f,1


ET,1,285					! Tretraeder ohne Mittelknoten

et,2,205

young1	= 1000
MP,ex,1,young1
MP,PRXY,1,0.4



*do,i,1,2
sigma_max%i%  = 1
maxtau%i%     = 1
choesion%i%   = 0.01
ucn%i%	  = sigma_max%i%/young1*2
ucnt%i%       = sigma_max%i%/young1*2
u_n%i%        = sigma_max%i%/young1
stabfact%i%   = 1e-3
beta%i%=1
*enddo


! 1. Korngrenzen
! 2. Matrix
! 3. Korn

*do,i,1,2

ET,100+i*2-1,170   				! Target
ET,100+i*2,173   					! Kontakt
keyopt,100+i*2,2,1
keyopt,100+i*2,5,1
keyopt,100+i*2,7,1  
keyopt,100+i*2,9,1
keyopt,100+i*2,10,2
keyopt,100+i*2,12,5

tb,czm,100+i*2-1,,,cbdd
tbdata,1,sigma_max%i%,ucn%i%,maxtau%i%,ucnt%i%,stabfact%i%,beta%i%
r,100+i*2-1,,,-sigma_max%i%/(u_n%i%),,,0		! 0= auto Pinball
RMODIF,100+i*2-1,13,choesion%i%   		! Cohesion
RMODIF,100+i*2-1,21,0.01 			! dynamische Kontaktdämpfung
mp,mu,100+i*2-1,0.3				! Reibung Mörtel Mauerwerk			

*enddo


!maßgebende Verschiebung
weg = max(ucn1,ucn2)




block,,2,,1,,1
vgen,10,all,,,2
vgen,5,all,,,,1
vgen,5,all,,,,,1

esize,4
vmesh,all


!modmsh,detach
!aslv,s
!lsla,s
!ksll,s
!vdel,all
!adel,all
!ldel,all
!kdel,all



alls
numm,kp
numm,node
numm,kp
numm,node
cm,alt,elem
cm,alt1,elem
numc,elem
numc,node
*get,anze,elem,,count
*get,anzn,node,,count

mat,2
type,2

/VIEW,1,1,1,1
eplo

imme,0
*do,i,1,anze
*get,mine,elem,mine,NXTH
esel,s,,,mine
nsle,s,1
NGEN,2,10000+anzn*i,all
EGEN,2,10000+anzn*i,all
cmsel,s,alt,elem
*enddo
imme,1

/eof

cmsel,s,alt
cmdele,alt
edel,all
ekill,all

alls
type,101
mat,101
real,101
alls
esurf

type,102
esurf
alls

type,103
mat,103
real,103
ESEL,S,TYPE,,102
NSLE,S  
ESLN,S,0
ESURF   
TYPE,104
ESEL,S,TYPE,,101
NSLE,S  
ESLN,S,0
ESURF   

alls


*get,minx,node,,mnloc,x
*get,mixx,node,,mxloc,x
nsel,s,loc,x,minx
d,all,ux
*get,minn,node,,num,min
nsel,r,,,minn
d,all,all
nsel,s,loc,x,mixx
d,all,ux
*get,minn,node,,num,min
nsel,r,,,minn
d,all,all
alls




!---------------------
/com, Lösungsteil
!---------------------

/solu

esel,s,live
nsle,s,1

time,1e-6
outres,all,all
sstif,on
neqit,50
RESCONTRL,DEFINE,ALL,1,1
NSUBST,10,10000,5
AUTOTS,on
solv


nsel,s,loc,x,mixx
d,all,ux,weg*2*mixx
alls
esel,s,live
nsle,s,1

NLHIST,PAIR,101,CONT,PRES,101 

nlgeom,on
solcontrol,on,on
time,1
AUTOTS,on
NSUBST,160,2000,30

solv

save



!---------------------
/com, Auswertung und Bestimmung des Kriteriums
!---------------------


/POST1  
esel,s,ename,,174
plns,cont,press