ABAQUS user subroutines for the simulation of viscoplastic - loicz

7.2 Appendix 2: Plate containing a hole — BODNER-PARTOM model coupled with damage
*HEADING
plate containing a crack
a= 5.0, r0= .50, N-alfa=10, N-r=15, K=12 D= .20
*PREPRINT, ECHO=NO, MODEL=NO, HISTORY=NO
*NODE
1, .20, .49846, .03923
5001, .00, .49846, .03923
101, .20, .53761, .04231
5101, .00, .53761, .04231
...
1380, .20, 3.47804, .00000
6380, .00, 3.47804, .00000
1480, .20, 4.16651, .00000
6480, .00, 4.16651, .00000
1580, .20, 5.00000, .00000
6580, .00, 5.00000, .00000
*ELEMENT, TYPE=C3D8I, ELSET=solid
*** 0 - PI/4
80, 5080, 80, 180, 5180, 5001, 1, 101, 5101
180, 5180, 180, 280, 5280, 5101, 101, 201, 5201
280, 5280, 280, 380, 5380, 5201, 201, 301, 5301
380, 5380, 380, 480, 5480, 5301, 301, 401, 5401
480, 5480, 480, 580, 5580, 5401, 401, 501, 5501
580, 5580, 580, 680, 5680, 5501, 501, 601, 5601
...
1477, 6477, 1477, 1577, 6577, 6478, 1478, 1578, 6578
1478, 6478, 1478, 1578, 6578, 6479, 1479, 1579, 6579
1479, 6479, 1479, 1579, 6579, 6480, 1480, 1580, 6580
**
** define node-set
**
*NSET, NSET=zlager0, GENERATE
6550, 6570, 1
*NSET, NSET=zlager1, GENERATE
1550, 1570, 1
*NSET, NSET=ylager0, GENERATE
6530, 6550, 1
*NSET, NSET=ylager1, GENERATE
1530, 1550, 1
**
** define element-set
**
*ELSET, ELSET=zlast, GENERATE
1410, 1429, 1
*ELSET, ELSET=ylast, GENERATE
1401, 1409, 1
1470, 1480, 1
**
** define materials and UMA
**
*SOLID SECTION, ELSET=SOLID, MATERIAL=VISCOPLA
1.,
*MATERIAL, NAME=VISCOPLA
*USER MATERIAL, CONSTANTS=19
**
** IN738LC, 850 deg C
******* BP-Model, units: h, MPa
********========================================
38