!********************************* Vorbereitung *******************************************

FINISH  
/CLEAR,START

/triad,ltop  					! Anzeigen des globelen KO-Systems im Display
/dev,font,2,times*new*roman,medium,r,14,,,,,   	! Manipulation der Anzeige im Fenster (z.B. Schrift)

/PNUM,LINE,1  					! Einschalten der Elementnummerierung und Farbkodierung im Plot (hier für Line)
/PNUM,KPOI,1  					! Einschalten der Elementnummerierung und Farbkodierung im Plot (hier für Keypoint)
/PNUM,area,1  					! Einschalten der Elementnummerierung und Farbkodierung im Plot (hier für Area)

/FILNAM,4PB_kont,1
/TITLE,	4 - Punkt Biegung (Kontakt)
/TRIAD,LTOP

!********************************* Pre-Processor ******************************************

/prep7

et,1,plane82,,,2
et,2,conta172,,1,,,3,
  keyopt,2,5,3
  keyopt,2,9,0
  keyopt,2,10,2
  keyopt,2,12,0

et,3,targe169
mp,ex,1,169000					! E-Modul von Silizium 130 - 180 GPa (Hier 169000 N/mm2)
mp,dens,1,2390e-12				! Dichte in t/mm3
nuxy,1,0.226					! Querkontraktionszahl (Poisson-Konstante)
acel,,-9810					! 9810 mm/s2

mp,mu,2,0

! ---- Geometrie ----

hoehe   = 0.200
laenge  = 156
l_pruef = 55
l_auf	= 110

r_auf = 1
r_last = 1

Kraft = 10

! ---- Knotenpunkte ----

! Balken
k,1,0,0,0
k,2,laenge/2-l_pruef/2,0,0
k,3,laenge/2,0,0
k,4,laenge/2+l_pruef/2,0,0
k,5,laenge,0,0
k,6,laenge,hoehe,0
k,7,laenge/2+l_pruef/2,hoehe,0
k,8,laenge/2,hoehe,0
k,9,laenge/2-l_pruef/2,hoehe,0
k,10,0,hoehe,0

! ---- Linien ----

! Balken
l,1,2,280
l,2,3,230
l,3,4,230
l,4,5,280
l,5,6,3
l,7,6,280
l,8,7,230
l,9,8,230
l,10,9,280
l,1,10,3
l,2,9,3
l,3,8,3
l,4,7,3

! ---- Flächen ----

al,1,11,9,10		! Balken
al,2,12,8,11	
al,3,13,7,12
al,4,5,6,13

type,1
mat,1
amesh,1,6,1

! ---- Kontaktpaar Auflage links - Balken

*get,nodenum,node,,num,max
auflager_li = nodenum + 1

n,auflager_li,laenge/2-l_auf/2,-r_auf,0
tshape,circ
r,1,1,,1e-5
type,3
mat,2
real,1
e,auflager_li
esel,s,type,,3
esurf,,reve
allsel

*get,nodenum,node,,num,max
pilot_aufl_1 = nodenum + 1

n,pilot_aufl_1,laenge/2+l_pruef/2,hoehe+r_last,0
type,3
real,1
mat,2
tshape,pilo
e,pilot_aufl_1

lsel,s,,,1
nsll,s,1
!esln,s,0
type,2
mat,2
real,1
esurf,

allsel

! ---- Kontaktpaar Auflage rechts - Balken

! Auflager rechts

*get,nodenum,node,,num,max
auflager_re = nodenum + 1

tshape,circ
r,2,1,,1e-5
type,3
mat,2
real,2
n,auflager_re,laenge/2+l_auf/2,-r_auf,0
e,auflager_re

*get,nodenum,node,,num,max
pilot_aufl_2 = nodenum + 1

n,pilot_aufl_2,laenge/2+l_pruef/2,hoehe+r_last,0
type,3
real,2
mat,2
tshape,pilo
e,pilot_aufl_2

lsel,s,,,4
nsll,s,1
!esln,s,0
type,2
mat,2
real,2
esurf,

allsel

! ---- Kontaktpaar Last links - Balken

! Last links

*get,nodenum,node,,num,max
last_li = nodenum + 1

tshape,circ
r,3,1,,1e-5
type,3
mat,2
real,3
n,last_li,laenge/2-l_pruef/2,hoehe+r_last,0
e,last_li

*get,nodenum,node,,num,max
pilot_last_1 = nodenum + 1

n,pilot_last_1,laenge/2-l_pruef/2,hoehe+r_last,0
type,3
real,3
mat,2
tshape,pilo
e,pilot_last_1

lsel,s,,,8,9,1
nsll,s,1
!esln,s,0
type,2
mat,2
real,3
esurf,

allsel


! ---- Kontaktpaar Last rechts - Balken

! Last rechts

*get,nodenum,node,,num,max
last_re = nodenum + 1

tshape,circ
r,4,1,,1e-5
type,3
mat,2
real,4
n,last_re,laenge/2+l_pruef/2,hoehe+r_last,0
e,last_re

*get,nodenum,node,,num,max
pilot_last_2 = nodenum + 1

n,pilot_last_2,laenge/2+l_pruef/2,hoehe+r_last,0
type,3
real,4
mat,2
tshape,pilo
e,pilot_last_2

lsel,s,,,6,7,1
nsll,s,1
!esln,s,0
type,2
mat,2
real,4
esurf,

allsel

! ************************************************ Solution *********************************************************

/solu

nsel,s,,,pilot_aufl_1
d,all,all

allsel

nsel,s,,,pilot_aufl_2
d,all,all

allsel

nsel,s,,,pilot_last_1
d,all,ux,0
d,all,rotz,0

allsel

nsel,s,,,pilot_last_2
d,all,ux,0
d,all,rotz,0

allsel

ksel,s,,,3
nslk,s,1
d,all,ux,0

allsel

D,pilot_last_1,UY,-10
D,pilot_last_2,UY,-10
!F,pilot_last_1,FY,-Kraft
!F,pilot_last_2,FY,-Kraft

allsel

! -- Für Statikanalyse --
antype,static,new
nlgeom,on

autots,on
nsubst,50,200,50
!time,force
outres,all,all

SOLVE

fini