*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
 1823,1824,1825,1826,1827,1828,1829,1830
*Nset,nset=endBox
 36,49,51,54,57,59,60,61,62,63,64,65
*Elset,elset=elemBox,generate
 1, 1572, 1
*Nset,nset=nodebox,elset=elemBox
*Elset,elset=elemPin,generate
 1573, 4257, 1
*Nset,nset=nodePin,elset=elemPin
*Elset,elset=innerPin_S3
 4242,
*Elset,elset=innerPin_S4,generate
 1579, 2657, 7
*Elset,elset=innerPin_S4
 2667,2681,2698,2723,2751,2779,2805,2830,2857,2882,2913,2939,
 2969,2997,3023,3053,3081,3112,3139,3165,3195,3221,3250,3277,
 3300,3327,3354,3377,3405,3432,3457,3483,3514,3540,3564,3592,
 3619,3644,3673,3699,3723,3753,3780,3808,3834,3862,3886,3910,
 3937,3961,3983,4007,4026,4042,4057,4072,4086,4103,4117,4129,
 4140,4151,4163,4175,4189,4210,4211,4212,4214,4237,4238,4239,
 4240,4241
*Elset,elset=surfBox_S1
 1457,1371,1318,1246,1191,1135,1075,1020, 946, 882, 795, 754, 
  697, 640, 559, 492, 430, 366, 304, 244, 166,  99,  20,  17,
    2
*Elset,elset=surfBox_S2
  977, 961, 946, 962, 964, 981, 984, 985, 983, 975, 937, 915, 
  898, 882, 899, 917, 923, 925, 922, 909, 891, 853, 832, 814,
  795, 815, 834, 837, 838, 860, 862, 849, 804, 782, 768, 754,
  769, 784, 790, 792, 791, 789, 779, 745, 726, 712, 697, 713,
  728, 732, 734, 733, 731, 724, 689, 669, 654, 640, 655, 671,
  675, 677, 674, 666, 650, 614, 592, 575, 559, 576, 594, 598,
  600, 599, 597, 586, 549, 527, 508, 492, 509, 514, 515, 535,
  539, 537, 523, 483, 463, 446, 430, 447, 449, 450, 468, 472,
  470, 460, 419, 398, 382, 366, 383, 400, 406, 408, 407, 405,
  394, 357, 337, 320, 304, 321, 323, 324, 342, 346, 344, 335,
  296, 275, 259, 244, 260, 277, 281, 283, 280, 272, 256, 219,
  198, 182, 183, 200, 204, 206, 205, 203, 194, 157, 137, 115,
  166,  99,1457,1455,1454,1447,1417,1399,1385,1371,1386,1401,
 1407,1409,1405,1389,1388,1363,1344,1330,1318,1331,1346,1350,
 1352,1349,1341,1327,1296,1275,1260,1246,1261,1277,1283,1285,
 1284,1282,1272,1237,1218,1204,1191,1205,1220,1224,1226,1225,
 1223,1215,1182,1164,1150,1135,1151,1166,1170,1172,1171,1169,
 1161,1126,1106,1091,1075,1092,1108,1112,1114,1113,1111,1104,
 1067,1047,1033,1020,1034,1049,1053,1055,1052,1042,1028, 997,
  116, 119, 121, 123, 120, 118, 106,  64,  39,  28,  20,  29,
   41,  46,  49,  47,  44,  32,  12,   6,   3,   2,  31
*Elset,elset=surfBox_S3
 1468,1456,1433,1370,1408,1381,1317,1351,1312,1245,1286,1255,
 1190,1227,1199,1134,1173,1145,1074,1115,1087,1019,1013,1054
  945, 986, 957, 881, 924, 874, 794, 839, 825, 753, 763, 696,
  735, 708, 639, 676, 633, 558, 601, 569, 491, 538, 502, 429,
  471, 440, 365, 409, 376, 303, 345, 316, 243, 282, 165, 207,
  176, 793, 238,  98, 122,  86,  19,  50,  21,   1,  25
*Elset,elset=surfBox_S4
 1370,1317,1245,1190,1134,1074,1019
  945, 881, 794, 753, 696, 639, 558, 491, 429, 365, 303, 243,
  165,  98,   1,   19,  5,   8,  11
*Elset,elset=surfPin_S3
 4225,4255,4174,4088,4004,3930,3858,3800,3741,3668,3588,3512,
 3470,3395,3323,3247,3185,3108,3048,2964,2911,2778,2775
*Elset,elset=surfPin_S4
 4242,4231,4226,4221,4222,4223,4224,4225,4220,4216,4217,4218,
 4219,4232,4235,4246,4250,4254,4257,4256,4255,4251,4247,4243,
 4233,4207,4188,4162,4150,4138,4126,4115,4099,4085,4070,4056,
 4069,4083,4097,4101,4102,4100,4075,4062,4050,4036,4021,4000,
 3981,3959,3936,3958,3979,3998,4003,4005,4006,4001,3982,3976,
 3955,3932,3925,3904,3883,3861,3884,3906,3927,3933,3935,3934,
 3909,3908,3907,3896,3874,3850,3827,3802,3779,3803,3829,3852,
 3857,3859,3860,3855,3833,3832,3831,3818,3793,3770,3744,3722,
 3745,3772,3795,3805,3807,3806,3778,3776,3752,3750,3747,3738,
 3717,3695,3672,3696,3719,3740,3748,3751,3749,3720,3698,3693,
 3669,3666,3657,3636,3614,3591,3615,3638,3659,3667,3670,3671,
 3663,3642,3618,3610,3586,3580,3559,3537,3513,3538,3561,3582,
 3587,3589,3590,3584,3563,3562,3551,3528,3501,3477,3454,3431,
 3455,3479,3503,3506,3508,3510,3511,3509,3507,3505,3498,3474,
 3451,3428,3404,3429,3453,3476,3481,3482,3480,3448,3426,3403,
 3401,3398,3391,3371,3349,3326,3350,3373,3393,3399,3402,3400,
 3375,3353,3352,3345,3321,3314,3294,3272,3249,3273,3296,3316,
 3322,3324,3325,3318,3298,3276,3275,3262,3237,3214,3187,3164,
 3188,3216,3239,3244,3246,3248,3245,3241,3219,3193,3190,3181,
 3159,3136,3111,3137,3161,3183,3191,3194,3192,3163,3162,3156,
 3132,3106,3098,3072,3045,3022,3046,3074,3100,3107,3109,3110,
 3104,3080,3079,3076,3049,3042,3019,2994,2968,2995,3021,3044,
 3050,3052,3051,3038,3016,2992,2965,2962,2957,2934,2906,2881,
 2907,2936,2959,2963,2966,2967,2961,2952,2929,2901,2878,2855,
 2850,2827,2804,2828,2852,2876,2899,2910,2912,2904,2880,2874,
 2849,2826,2803,2778,2777,2776,2775,2800,2823
**
** 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,unsymm=yes
 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
*print, solve=yes
*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
*print, solve=yes
*End step
** ----------------------------------------------------------------
** 
** STEP: Step-3
** 
*Step,name=Step-3,nlgeom
 pressurize
*Static
 0.5,1.0
*Pressure penetration,master=surfPin,slave=surfBox
 36,4731,100.
*Dsload
 innerPin,P,100.
*print, solve=yes
*End step
** ----------------------------------------------------------------
** 
** STEP: Step-4
** 
*Step,name=Step-4,nlgeom
 pull
*Static
 0.1,1.0
*Boundary
 endBox,2,2,-0.1
*print, solve=yes
*End step