$$ delta_s horizontal direction
$$ delta   certical direction

+prog template urs:1
kopf
sto#mesh_size 0.002
sto#dimension 0.2

sto#mat_beto 37.4
sto#f_ctm 3.42
sto#epsc -3.5
sto#dehn_ksv 'sl'

sto#steu_iter_opt 'iter_3'
sto#prob 'nonl'
sto#step 0.0001
sto#hor_load 5
sto#load_step 0.1
sto#iterationen 100
sto#elemente 32
ende
$f1ct = 0.90 * fct,sp  = 3.78*0.9
$fc = 46.19 150 mm cubes
$fc,cyl = 0.85*fc,cu,150 und f1c,cu = 0.95*fc,cyl => f1c = 0.81*fc,cu
$+16explizite alte Elemente beibehalten
$+32nur Dreieckselemente erzeugen
$+64nur Viereckselemente erzeugen
$+96gemischte Netze erzeugen

+prog aqua urs:3
kopf
norm en 199x-200x
let#fcr #mat_beto
beto nr 1 art c fcn 40 fcr #fcr typr b muec 0.2 fct #f_ctm
$beto nr 1 art c fcn 40 fcr 40 typr b muec 0.2 fct 3.3 gf 100 ec 30000 ecr 30000
$beto nr 2 art c fcn 40 typr b muec 0.2 feqr 1 feqt 0
$ fcn 40 dient hier nur als platzhalter -> verwendet wird Materialkurve die über fcr definiert wird.
arbl eps gebr
echo mat voll
ende

+prog sofimshc urs:2
kopf
syst art 3d gdiv 25000 gdir posz fix PY
steu mesh 2+#elemente
$steu fein w2 -0.5
steu hmin #mesh_size
steu deln 0
$echo quad ja
$echo syst ja
echo knot voll
sto#notch_length 0.025
sto#notch_width 0.005

let#col0 0
let#col1 #notch_length
let#col2 #dimension-#notch_length
let#col3 #dimension

let#row1 0
let#row2 -1*(0.5*#dimension-0.5*#notch_width)
let#row3 -0.5*#dimension
let#row4 -1*(0.5*#dimension+0.5*#notch_width)
let#row5 -1*(#dimension)

let#dh 0.002
$let#dh #notch_width/2
!*!Label Strukturpunkte
spt  nr  x              z
     10  #col0          #row1
     15  #col3          #row1
     20  #col0          #row2
     22  #col1-#dh        #row2
     32  #col1          #row3
     33  #col2          #row3
     23  #col2+#dh          #row2
     25  #col3          #row2
     40  #col0          #row4
     42  #col1-#dh          #row4
     43  #col2+#dh          #row4
     45  #col3          #row4
$     50  #col0          #row5
$     55  #col3          #row5

$LOOP#i 3
$sto#temp =c(#i)
$txe #temp
$ENDLOOP
$sto#node_vector(#i) #node_num

sto#num_ele_edge #dimension/#mesh_size
sto#sdiv_half 0.5*(#num_ele_edge)
sto#num_ele_top 0.5*#dimension/#mesh_size


loop#i #num_ele_top+1
sto#node_vector(#i) 500+#i
spt nr 500+#i x #col0+2*#mesh_size*#i z #row5
let#lastindex #i
txe #i
endloop


let#grp_val 10

sln nr 10 npa 10 npe 15 grp #grp_val fix PZ
sln nr 11 npa 10 npe 20 grp #grp_val
sln nr 12 npa 15 npe 25 grp #grp_val sdiv -#sdiv_half+1

let#grp_val 20
sln nr 20 npa 20 npe 22 grp #grp_val
sln nr 24 npa 23 npe 25 grp #grp_val

let#grp_val 30
let#middle 0.5*(#row2+#row4)
let#pt3_x #col1+0.005

!*!Label Kreislinien
$sln nr 31 npa 22 npe 32 grp #grp_val sdiv 0.001
$$slnb x1 #col1-#dh z1 #row2 x2 #col1 z2 #row3 r 0.005 ny 1
$slnb x1 #col1-#dh z1 #row2 x2 #col1 z2 #row3 r #notch_width/2 ny 1
$sln nr 32 npa 32 npe 42 grp #grp_val sdiv 0.001
$$slnb x1 #col1 z1 #row3 x2 #col1-#dh z2 #row4 r 0.005 ny 1
$slnb x1 #col1 z1 #row3 x2 #col1-#dh z2 #row4 r #notch_width/2 ny 1

sln nr 31 npa 22 npe 42 grp #grp_val sdiv 0.001
slnb x1 #col1-#dh z1 #row2 x2 #col1-#dh z2 #row4 r #notch_width*0.55 ny 1

sln nr 33 npa 23 npe 33 grp #grp_val sdiv 0.001
$slnb x1 #col2+#dh z1 #row2 x2 #col2 z2 #row3 r 0.005 ny -1
slnb x1 #col2+#dh z1 #row2 x2 #col2 z2 #row3 r #notch_width/2 ny -1
sln nr 34 npa 33 npe 43 grp #grp_val sdiv 0.001
$slnb x1 #col2+#dh z1 #row3 x2 #col2 z2 #row4 r 0.005 ny -1
slnb x1 #col2+#dh z1 #row3 x2 #col2 z2 #row4 r #notch_width/2 ny -1

let#grp_val 40
sln nr 40 npa 40 npe 42 grp #grp_val
sln nr 44 npa 43 npe 45 grp #grp_val

!*!Label Linien oben
let#grp_val 50

sln nr 48 npa 40 npe 500 grp #grp_val sdiv -#sdiv_half+1 fix PX
sln nr 49 npa =node_vector(#lastindex) npe 45 grp #grp_val
$sln nr 50 npa 50 npe 55 grp #grp_val

loop#i #num_ele_top
sln nr 500+#i npa 500+#i npe 501+#i sdiv -1
endloop


!*!Label  sar
sar nr 1 mnr 1 t 50[mm] grp 2
sarb typ aus nl 10,11,12
sarb typ aus nl 20,24
$sarb typ aus nl 31,32,33,34
sarb typ aus nl 31,33,34
$sarb typ aus nl 40,44,48,49,50
sarb typ aus nl 40,44,48,49
loop#i #num_ele_top
sarb typ aus nl 500+#i
endloop

ende

+prog sofimsha urs:4
kopf Kopplungen
syst art rest
let#row2 -1*(0.5*#dimension-0.5*#notch_width)
knot nr x z fix
     1   #dimension*1.2  0.5*#row2   PZ
     2   #dimension*0.5  -#dimension*1.5  PX
grup 19 bez 'Kopplungen'
loop#i #sdiv_half
    let#x_pos (#i+1)*#row2/#sdiv_half
    getn x #dimension y 0 z #x_pos var node typ node tol 0.002[m]
    KNOT nr #node FIX KPPX NR1 1
endloop

sto#postion_ds_x 0.01
sto#postion_ds_z 0.05
sto#postion_d_x 0.03
sto#postion_d_z 0.0325

getn x #postion_ds_x z -#postion_ds_z var quad typ quad tol 0.002[m]
sto#ele_ds_BL #quad
getn x #dimension-#postion_ds_x z -#dimension+#postion_ds_z var quad typ quad tol 0.002[m]
sto#ele_ds_TR #quad

getn x #postion_d_x z -0.5*#dimension+#postion_d_z var quad typ quad tol 0.002[m]
sto#ele_d_BL #quad
getn x #postion_d_x z -0.5*#dimension-#postion_d_z var quad typ quad tol 0.002[m]
sto#ele_d_TL #quad

getn x #dimension-#postion_d_x z -0.5*#dimension+#postion_d_z var quad typ quad tol 0.002[m]
sto#ele_d_BR #quad
getn x #dimension-#postion_d_x z -0.5*#dimension-#postion_d_z var quad typ quad tol 0.002[m]
sto#ele_d_TR #quad

$loop#i #num_ele_edge+1
$    let#x_pos #i*#dimension/#num_ele_edge
$    getn x #x_pos y 0 z -#dimension var node typ node tol #mesh_size*0.5
$    sto#node_vector(#i) #node
$endloop

ende


+prog sofiload urs:5
kopf
lf nr 1 fakg 0
    knot nr 1 typ pxx p1 -#load_step
lf nr 10 fakg 0 bez 'H-Kraft 0'
    knot nr 1 typ pxx p1 -0.01
lf nr 11 fakg 0 bez 'H-Kraft05 voll'
    knot nr 1 typ pxx p1 -5
lf nr 12 fakg 0 bez 'H-Kraft10 voll'
    knot nr 1 typ pxx p1 -10
lf nr 15 fakg 0 bez 'H-Displacement 0 kN'
    knot nr 1 typ wxx p1 -0.001 [mm]
lf nr 16 fakg 0 bez 'H-Displacement 5 kN'
    knot nr 1 typ wxx p1 -0.009443 [mm]
lf nr 17 fakg 0 bez 'H-Displacement 10 kN'
    knot nr 1 typ wxx p1  -0.02876 [mm]
    $knot nr 1 typ wxx p1  -0.022 [mm]
lf nr 20 fakg 0
    loop#i node_vector
        let#temp_node =node_vector(#i)
        knot nr #temp_node typ wzz p1 -#step [mm]
    endloop

if #hor_load >7.5
sto#offset 2
elseif #hor_load > 1
sto#offset 1
elseif #hor_load < 1
sto#offset 0
endif
ende

-prog ase urs:13
kopf 10 - horizontal force 0.01
syst prob #prob nmat ja
steu core 12
steu solv 4
steu warn 353
STEU WARN 610
dehn ksv #dehn_ksv
steu beto w2 0 w4 #f_ctm w6 75
lf nr 10
ende

-prog ase urs:10
kopf 11 - horizontal force 5
syst prob #prob nmat ja iter 500 tol 0.0005 FMA4 0.5 FMA8 0.25
steu core 12
steu solv 4
steu warn 353
STEU WARN 610
dehn ksv #dehn_ksv
steu beto w2 0 w4 #f_ctm w6 75
lf nr 11
ende

+prog ase urs:7
kopf 20 - vertikal disp
syst prob #prob nmat ja iter 500 tol 0.0005 FMA4 0.5 FMA8 0.25
steu core 12
steu solv 4
steu warn 353
STEU WARN 610
dehn ksv #dehn_ksv
steu beto w2 0 w4 #f_ctm w6 75
lf nr 20
ende

-prog ase urs:14
kopf trag vertikale Verschiebung
syst prob #prob nmat ja iter #iterationen FMA4 0.5 FMA8 0.25
steu core 12
steu solv 4
steu warn 353
STEU WARN 610
if #steu_iter_opt=='iter_0'
STEU ITER 0
elseif #steu_iter_opt=='iter_1'
STEU ITER 1
elseif #steu_iter_opt=='iter_3'
steu iter 3
endif
dehn ksv #dehn_ksv
steu beto w2 0 w4 #f_ctm w6 75

trag ntra 150 dfak 1 pro 1 eg nein prim nein epsc #epsc fake 150 dmin 1
LF 1000
    lc nr 10+#offset trag nein
    lc nr 20 trag alle
ende
sto#loadsteps #ASE_ITER(1)
ende