*Heading
*Node, nset=nall
*Include,Input=threadedconnector_cax4_n.inp
*Element,type=cax4,elset=eall
*Include,Input=threadedconnector_cax4_e.inp
*Solid Section,Elset=eall,Material=ss95
**
** node and element set definitions
**
*Nset,nset=endPin
*Nset,nset=endBox
*Elset,elset=elemBox,generate
*Nset,nset=nodebox,elset=elemBox
*Elset,elset=elemPin,generate
*Nset,nset=nodePin,elset=elemPin
*Elset,elset=innerPin_S3
*Elset,elset=innerPin_S4,generate
*Elset,elset=innerPin_S4
*Elset,elset=surfBox_S1
*Elset,elset=surfBox_S2
*Elset,elset=surfBox_S3
*Elset,elset=surfBox_S4
*Elset,elset=surfPin_S3
*Elset,elset=surfPin_S4
**
** surface definitions
**
*Surface,name=innerPin
 innerPin_S3,S3
 innerPin_S4,S4
*Surface,name=surfBox
 surfBox_S1,S1
 surfBox_S2,S2
 surfBox_S3,S3
 surfBox_S4,S4
*Surface,name=surfPin
 surfPin_S3,S3
 surfPin_S4,S4
**
** Interaction And Interaction Properties
**
*Contact pair,interaction=slide
 surfBox,surfPin
*Surface interaction,name=slide
*Friction
 0.0,
**
** Material
**
*Material,name=ss95
*Elastic
 30.E6,0.3
*Plastic
 92510,0.00
 94794,0.00176067
 97303,0.00467793
 109980,0.03573816
 122925,0.08264343
*Restart,Write,Freq=100
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step,nlgeom
 resolve contact interference
*Static
 0.01,1.0,0.001,0.5
*Contact interference,type=contact pair,shrink
 surfBox,surfPin
**
** BOUNDARY CONDITIONS
**
*Boundary
 endPin,2
 endBox,2
**
** OUTPUT REQUESTS
**
*El print,freq=0
*Node print,freq=0
*Contact print,freq=1
 CMS,
 CAREA,
 CTRQ,
*Output,history,freq=1
*Contact output
 CMS,CAREA,CTRQ
*Output,field,freq=99
*Element output
 S,PEEQ
*Node output
 U,
*Contact output
 CSTRESS,CDISP,PPRESS
*End step
** ----------------------------------------------------------------
**
** STEP: Step-2
**
*Step,name=Step-2,nlgeom
 change friction coefficient
*Static
 1.0,1.0
*Change friction,interaction=slide
*Friction
 0.1
*End step
** ----------------------------------------------------------------
**
** STEP: Step-3
**
*Step,name=Step-3,nlgeom
 pressurize
*Static
 0.5,1.0
*Pressure penetration,master=surfPin,slave=surfBox
*Dsload
 innerPin,P,100.
*End step
** ----------------------------------------------------------------
**
** STEP: Step-4
**
*Step,name=Step-4,nlgeom
 pull
*Static
 0.1,1.0
*Boundary
 endBox,2,2,-0.1
*End step