*HEADING ** ** Analytical contact with linear softening in 3D. ** *NODE,NSET=REFNOD 1,2.,0.,0.5 *RIGID BODY,REF NODE=1,ANALYTICAL SURFACE=SB1,ISOTHERMAL=YES *HEATCAP,ELSET=PTCAPAC 0.004, *ELEMENT,TYPE=HEATCAP,ELSET=PTCAPAC 1,1 *BOUNDARY 1,1,6 ** *NODE,NSET=NB2 101,0.,1.1,0. 102,1.,1.1,0. 103,2.,1.1,0. 104,3.,1.1,0. 105,4.,1.1,0. 111,0.,2.1,0. 112,1.,2.1,0. 113,2.,2.1,0. 114,3.,2.1,0. 115,4.,2.1,0. 121,0.,1.1,1. 122,1.,1.1,1. 123,2.,1.1,1. 124,3.,1.1,1. 125,4.,1.1,1. 131,0.,2.1,1. 132,1.,2.1,1. 133,2.,2.1,1. 134,3.,2.1,1. 135,4.,2.1,1. *ELEMENT,ELSET=EB2,TYPE=SC8RT 101,101,102,112,111,121,122,132,131 *ELGEN,ELSET=EB2 101,4,1,1 *NSET,NSET=LOAD_N 111,112,113,114,115,131,132,133,134,135 ** *ELSET,ELSET=EALL EB2, ** *SHELL SECTION,ELSET=EALL,MATERIAL=MAT1 1.0 *MATERIAL,NAME=MAT1 *ELASTIC 1.E6,0.2 *DENSITY 1., *CONDUCTIVITY 10., *SPECIFIC HEAT 0.001, *PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFAN=1.E-6 *SURFACE,TYPE=CYLINDER,NAME=SB1 -1.,1.,2.,0.,1.,2. -1.,1.,1. START,0.,0. LINE,6.,0. *SURFACE,TYPE=ELEMENT,NAME=SB2 EB2, ** *AMPLITUDE,NAME=SMOOTH,DEFINITION=SMOOTH STEP 0.,0.,1.E-3,-0.1,5.E-3,1. ** *NSET, NSET=QA_TEST_ALLNODES NB2, REFNOD ** *INITIAL CONDITIONS,TYPE=TEMP REFNOD, 100. NB2, 0. ** *STEP *DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT ,5.E-3 ** *BOUNDARY,AMPLITUDE=SMOOTH LOAD_N,2,2,0.001 ** *CONTACT PAIR,INTERACTION=LINEAR, MECHANICAL CONSTRAINT=PENALTY SB1,SB2 *SURFACE INTERACTION,NAME=LINEAR,PADTHICKNESS=0.1 *GAP CONDUCTANCE 1.,0. 1.,0.1 **Gap conductance = k_interface / t_interface *GAP HEAT GENERATION 0., *CONTACT DAMPING,DEFINITION=CRITICAL DAMPING FRACTION 0., *SURFACE BEHAVIOR,NO SEPARATION,PRESSURE-OVERCLOSURE=LINEAR 1.E6, **Pressure-overclosure slope = E_interface / t_interface *FRICTION,ROUGH,SHEAR TRACTION SLOPE=416666.7 **Shear traction slope = G_interface / t_interface ** = E_int / t_interface / (2*(1+nu_interface)) ** *OUTPUT, FIELD, NUMBER INTERVAL=1 *NODE OUTPUT, NSET=QA_TEST_ALLNODES U, NT *OUTPUT,FIELD,VARIABLE=PRESELECT *OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.005 *END STEP