** vim:ft=abaqus
** All values in SI units: m,s,kg,N,Pa,J
** (length,time,mass,force,stress,energy).

*HEADING
Steel pin in shell

** Geometry
*INCLUDE, INPUT=all.msh
*INCLUDE, INPUT=botdep.sur
*INCLUDE, INPUT=botindep.sur
*INCLUDE, INPUT=load.nam
*INCLUDE, INPUT=pin.nam
*INCLUDE, INPUT=shell.nam
*INCLUDE, INPUT=topdep.sur
*INCLUDE, INPUT=topindep.sur

** Rigid body constraints; should be defined before step!
*RIGID BODY, NSET=Nload, REF NODE=1

*BOUNDARY
** Restrain ref node in x and y.
1,1,2
** Restrain shell nodes
Nshell,1,3

** Material properties
*MATERIAL, NAME=steel_1.4658
*ELASTIC, TYPE=ISO
197e9,0.25,293
189e9,0.25,373
***PLASTIC
**7e+08,0.002,293
**7.36667e+08,0.00493,293
**7.73333e+08,0.01163,293
**8.1e+08,0.02636,293
**8.46667e+08,0.05762,293
**8.83333e+08,0.12187,293
**9.2e+08,0.25,293
*DENSITY
7800
*EXPANSION, TYPE=ISO
12.5e-6,293
12.5e-6,373
13.0e-6,473
13.5e-6,573

*MATERIAL, NAME=alumina92
*ELASTIC,TYPE=ISO
280e9,0.23,293
*DENSITY
3600
*EXPANSION, TYPE=ISO
7.2e-6,293

** Apply material properties to element sets.
*SOLID SECTION, ELSET=Epin, MATERIAL=steel_1.4658
** The shell needs properties, even if they're basically not used.
*SOLID SECTION, ELSET=Eshell, MATERIAL=alumina92

** Contact
*SURFACE INTERACTION,NAME=SI1
*SURFACE BEHAVIOR,PRESSURE-OVERCLOSURE=LINEAR
1000e9,,
*CONTACT PAIR,INTERACTION=SI1,TYPE=SURFACE TO SURFACE
Stopdep,Stopindep
*CONTACT PAIR,INTERACTION=SI1,TYPE=SURFACE TO SURFACE
Sbotdep,Sbotindep

** Concentrated load on rigid body **
*STEP, NLGEOM
*STATIC
0.1,1
** Force = 11772 N (2 pinnen!), 5886 N per pen
** Moment = 5279 Nm
** arm = round(5279/11772, 3) = 0.448 m
*BOUNDARY
** z = -0.01; fz = -3173 N
** z = -0.015; -5228 N
** z = -0.016; -5640 N
** z = -0.0165; -5846 N
** z = -0.01655; -5867 N
** z = -0.01656; -5871 N
** z = -0.01658; -5879 N
** z = -0.0166; -5887 N
1,3,3,-0.0166
*NODE FILE, FREQUENCY=1000
U,RF
*EL FILE, FREQUENCY=1000
ZZS,ME
*NODE PRINT, NSET=Nload, FREQUENCY=1000
RF
*END STEP