*HEADING
Axisymmetrie, Mohr-Coulomb, CAX4
Modell mit Pfahlkopfelement über GOK
Statische Berechnung, kraftgesteuert
Primärspannungszustand -> GEOSTATIC
mit KONTAKT
************
***** Knoten
************
*NODE, NSET=top
90001,0,0.03
90002,0.0125,0.03
90003,0.026,0.03
*NODE, NSET=fuermaster1
50001,0,0
50002,0.0125,0
50003,0.026,0
*NODE, NSET=knoten1
1,0,0
3,0.025,0
*NODE, NSET=knoten1a
4,0.03,0
*NODE, NSET=knoten1b
14,0.5,0
*NGEN, NSET=knoten1
1,3,1
*NFILL,NSET=knoten1, BIAS=0.78
knoten1a,knoten1b,10,1
*NCOPY,CHANGE NUMBER=10000, OLD SET=knoten1, SHIFT, NEW SET=knoten2
0,-1
0,0
*NCOPY,CHANGE NUMBER=10100, OLD SET=knoten1, SHIFT, NEW SET=knoten3
0,-1.005
0,0
*NCOPY,CHANGE NUMBER=11100, OLD SET=knoten1, SHIFT, NEW SET=knoten4
0,-2
0,0
*NCOPY,CHANGE NUMBER=10000, OLD SET=fuermaster1, SHIFT, NEW SET=fuermaster2
0,-1.001
0,0
*NFILL, NSET=knoten1
knoten1,knoten2,10,1000
*NFILL, NSET=knoten1
knoten2,knoten3,1,100
*NFILL, NSET=knoten1, BIAS=0.68
knoten3,knoten4,10,100
*NFILL, NSET=fuermaster1
fuermaster1,fuermaster2,10,1000
*******************
***** Pfahlelemente
*******************
*ELEMENT,TYPE=CAX4,ELSET=z
50001,50001,51001,51002,50002
*ELGEN, ELSET=z
50001,2,1,1,10,1000,1000
*ELEMENT,TYPE=CAX4,ELSET=z
90001,90001,50001,50002,90002
90002,90002,50002,50003,90003
***************
***** Scherfuge
***************
*ELEMENT,TYPE=CAX4,ELSET=c
3,3,1003,1004,4
*ELGEN, ELSET=c
3,1,1,1,10,1000,1000
**********************************
***** Boden rechts neben Scherfuge
**********************************
*ELEMENT,TYPE=CAX4,ELSET=d
4,4,1004,1005,5
*ELGEN, ELSET=d
4,10,1,1,10,1000,1000
******************
***** Spitzendruck
******************
*ELEMENT,TYPE=CAX4,ELSET=e
10001,10001,10101,10102,10002
*ELGEN, ELSET=e
10001,2,1,1,1,1,1
***********************************************
***** unter Scherfuge, neben Spitzendruck, Rest
***********************************************
*ELEMENT,TYPE=CAX4,ELSET=f
10003,10003,10103,10104,10004
*ELGEN, ELSET=f
10003,11,1,1,11,100,100
************************
***** unter Spitzendruck
************************
*ELEMENT,TYPE=CAX4,ELSET=g
10101,10101,10201,10202,10102
*ELGEN, ELSET=g
10101,2,1,1,10,100,100
**************
***** Sections
**************
*SOLID SECTION, ELSET=z, MATERIAL=beton
1.,
*SOLID SECTION, ELSET=c, MATERIAL=sand
1.,
*SOLID SECTION, ELSET=d, MATERIAL=sand
1.,
*SOLID SECTION, ELSET=e, MATERIAL=sand
1.,
*SOLID SECTION, ELSET=f, MATERIAL=sand
1.,
*SOLID SECTION, ELSET=g, MATERIAL=sand
1.,
**************************************************
***** Material, Sand
***** gamma = 17 kN/m3, E = 12000 kN/m2, nü = 0,33
***** phi = 39°, psi = 0°, c = 0.1 kN/m2
**************************************************
*MATERIAL, NAME=sand
*DENSITY
1.7,
*ELASTIC, TYPE=ISO
12000,0.33
*MOHR COULOMB
39,0
*MOHR COULOMB HARDENING
0.1
*************************************************
***** Material, Beton
***** gamma = 25 kN/m3, E = 30000 MN/m2, nü = 0,2
*************************************************
*MATERIAL, NAME=beton
*DENSITY
2.5,
*ELASTIC, TYPE=ISO
3E7,0.2
*********************
***** Randbedingungen
*********************
*NSET, NSET=rand, GENERATE
1,11101,100
14,11114,100
50001,60001,1000
*NSET, NSET=rand
30001,90001
*NSET, NSET=unten, GENERATE
11101,11114,1
***********************
***** Pfahlkopfelemente
***********************
*ELSET, ELSET=kopf
90001,90002
**************************
***** Interaktionselemente
**************************
*ELSET, ELSET=elemaster1
90002
*ELSET, ELSET=elemaster2
59001
*ELSET, ELSET=elemaster3, GENERATE
50002,59002,1000
*ELSET, ELSET=elemaster4
59002
*ELSET, ELSET=eleslave1
4
*ELSET, ELSET=eleslave2, GENERATE
3,9003,1000
*ELSET, ELSET=eleslave3
3
*ELSET, ELSET=eleslave4
10001,10002
**************************************
***** Boden für Primärspannungszustand
**************************************
*ELSET, ELSET=init
c, d, e, f, g
*******************************************
***** Knoten k1 für Verschiebung des Pfahls
*******************************************
*NSET, NSET=k1
90001
*******************************************
***** Interaktion
*******************************************
*SURFACE, NAME=Master, TYPE=ELEMENT, TRIM=YES
elemaster1,s3
elemaster2,s2
elemaster3,s3
elemaster4,s2
*SURFACE, NAME=Slave, TYPE=ELEMENT
eleslave1,s4
eleslave2,s1
eleslave3,s4
eleslave4,s4
*CONTACT PAIR, ADJUST=0.004, INTERACTION=Kontakt
Slave, Master
*SURFACE INTERACTION, NAME=Kontakt
** mü = tan deltas = 2/3 phi = 0.49
** taumax = 1,3 KN?
***FRICTION, TAUMAX=
*FRICTION
0.49
***********************
******  Randbedingungen
***********************
*BOUNDARY
rand,1,1,0
unten,2,2,0
************************
***** Initial conditions
************************
*INITIAL CONDITIONS, TYPE=STRESS, GEOSTATIC
init, 0, 0, -36, -2, 0.5
z, 0, 0.03, -25.75, -1.03, 1
**************************
******  Berechnung, STEP 1
**************************
*STEP, UNSYMM=YES
*GEOSTATIC
***STATIC
*DLOAD
init, GRAV, 10, 0, -1 
z, GRAV, 10, 0, -1
*MONITOR, NODE=k1, DOF=2
*END STEP
************************
***** Berechnung, STEP 2
************************
*STEP, INC=10000, AMPLITUDE=RAMP, UNSYMM=YES
*STATIC
0.0001,1,1E-10,1E-1
*DLOAD
kopf, p4, 1000
*PRINT, PLASTICITY=YES
***CONTROLS, PARAMETERS=FIELD, FIELD=DISPLACEMENT
**0.02, 0.05
*OUTPUT, FIELD, VARIABLE=PRESELECT, FREQUENCY=1
*OUTPUT, HISTORY, VARIABLE=PRESELECT, FREQUENCY=1
*END STEP