*HEADING SIMPLE CONTACT TEST, COUPLED TEMPERATURE-DIPLACEMENT *NODE 1, 1., 5. 2, 1.5, 5. 3, 2., 5. 4, 1., 7.5 5, 2., 7.5 6, 1., 10. 7, 1.5, 10. 8, 2., 10. 9, 1., 5., 0.5 10, 2., 5., 0.5 11, 1., 10., 0.5 12, 2., 10., 0.5 13, 1., 5., 1. 14, 1.5, 5., 1. 15, 2., 5., 1. 16, 1., 7.5, 1. 17, 2., 7.5, 1. 18, 1., 10., 1. 19, 1.5, 10., 1. 20, 2., 10., 1. 21, 1., 12.5 22, 1., 11.25 23, 1., 10., 24, 1.25, 12.5, 25, 1.25, 10., 26, 1.5, 12.5 27, 1.5, 11.25 28, 1.5, 10. 29, 1., 12.5, 0.5 30, 1., 10., 0.5 31, 1.5, 12.5, 0.5 32, 1.5, 10., 0.5 33, 1., 12.5, 1. 34, 1., 11.25, 1. 35, 1., 10., 1. 36, 1.25, 12.5, 1. 37, 1.25, 10., 1. 38, 1.5, 12.5, 1. 39, 1.5, 11.25, 1. 40, 1.5, 10., 1., 44, 1.75, 12.5, 45, 1.75, 10. 46, 2., 12.5 47, 2., 11.25 48, 2., 10. 51, 2., 12.5, 0.5 52, 2., 10., 0.5 56, 1.75, 12.5, 1. 57, 1.75, 10., 1. 58, 2., 12.5, 1. 59, 2., 11.25, 1. 60, 2., 10., 1. 61, 1., 12.5 62, 1.5, 12.5 63, 2., 12.5 64, 1., 13.75 65, 2., 13.75 66, 1., 15. 67, 1.5, 15. 68, 2., 15. 69, 1., 12.5, 0.5 70, 2., 12.5, 0.5 71, 1., 15., 0.5 72, 2., 15., 0.5 73, 1., 12.5, 1. 74, 1.5, 12.5, 1. 75, 2., 12.5, 1. 76, 1., 13.75, 1. 77, 2., 13.75, 1. 78, 1., 15., 1. 79, 1.5, 15., 1. 80, 2., 15., 1. *NSET,NSET=BOTTOM 1,2,3,9,10,13,14,15, *NSET,NSET=TOP 66,67,68,71,72, 78,79,80, *NSET,NSET=BACK 1,2,3 66,67,68, *ELEMENT, TYPE=C3D20HT, ELSET=SOLID 1, 1, 3, 8, 6, 13, 15, 20, 18, 2, 5, 7, 4, 14, 17, 19, 16, 9, 10, 12, 11 2, 21, 23, 28, 26, 33, 35, 40, 38, 22, 25, 27, 24, 34, 37, 39, 36, 29, 30, 32, 31 3, 26, 28, 48, 46, 38, 40, 60, 58, 27, 45, 47, 44, 39, 57, 59, 56, 31, 32, 52, 51 4, 61, 63, 68, 66, 73, 75, 80, 78, 62, 65, 67, 64, 74, 77, 79, 76, 69, 70, 72, 71 *SOLID SECTION, ELSET=SOLID, MATERIAL=SOLMAT *MATERIAL, NAME=SOLMAT *ELASTIC 30.E6,0.3 *EXPANSION 10.E-6, *CONDUCTIVITY 43.0, *DENSITY 7700., *SPECIFIC HEAT 600., *SURFACE INTERACTION, NAME=GAP 1.0, *GAP CONDUCTANCE 1.,0. 1.,1. *GAP RADIATION 3.4664e-2, 3.4664e-2 1.,0. 1.,1. ** converted from old data: 1.E-9, 1.E-9 *GAP HEAT GENERATION 0.5,0.5 ** *PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16, STEFANB=5.6697E-8 ** *** nset definitions for prescribing temperatures. ** *NSET, NSET=NELEM1 1,3,13,15, 18,20,6, 8, *NSET, NSET=NELEM23 35,40,60, 23,28,48, 73,75,63 61, 78,80,66, 68, ** ** surfaces ** *surface, NAME=ASURF 1, S5 *surface, NAME=BSURF 2, S4 3, S4 *CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP BSURF, ASURF *surface, NAME=ATIED1 2, S6 *surface, NAME=ATIED2 3, S6 *surface, NAME=BTIED 4, S3 *TIE,NAME=TIE,POSITION TOLERANCE=0.01 **CONTACT PAIR, SMALL SLIDING, TIED, ADJUST=0.01, INTERACTION=GAP ATIED1, BTIED ATIED2, BTIED ** *** Apply different temperatures to the nodes above and *** below the contact surface ** *AMPLitude,name=step 0,1,100,1 *STEP *COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE *ENERGY FILE *EL FILE S, E, EE, TEMP, *NODE FILE RF, NT,RFL *ENERGY PRINT *EL PRINT ,SUM=NO S, E, EE, TEMP, *NODE PRINT,SUM=NO RF, NT,RFL U, *BOUNDARY BACK,3,3 BOTTOM,2,2 TOP,2,2 1, 1, 1 9, 1, 1 13,1, 1 66,1, 1 71,1, 1 78,1, 1 NELEM1, 11, 11, 700.0 NELEM23, 11, 11, 500.0 *CONTACT PRINT,MAIN=ASURF,SECONDARY=BSURF CSTRESS,CDISP HFL,HFLA,HTL SFDR,SFDRA *CONTACT FILE,MAIN=ASURF,SECONDARY=BSURF CSTRESS,CDISP HFL,HFLA,HTL SFDR,SFDRA *END STEP ** **** STEP 2. Push top block to the right. ** *STEP, NLGEOM ,INC=10 *COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE 0.1,1.0,0.1,0.1 *SOLUTION TECHNIQUE,TYPE=SEPARATED *CHANGE FRICTION,INTERACTION=GAP,amplitude=step *FRICTION 0.01, *BOUNDARY, OP=NEW BACK,3,3 23, 1, 1, 0.01 30, 1, 1, 0.01 35, 1, 1, 0.01 22, 1, 1, 0.01 34, 1, 1, 0.01 61, 1, 1, 0.01 69, 1, 1, 0.01 73, 1, 1, 0.01 64, 1, 1, 0.01 76, 1, 1, 0.01 66, 1, 1, 0.01 71, 1, 1, 0.01 78, 1, 1, 0.01 NELEM1, 11, 11, 700 NELEM23, 11, 11, 500 TOP, 2, 2 BOTTOM, 2, 2 1, 1, 1 9, 1, 1 13,1, 1 4,1,1 16,1,1 6,1,1 11,1,1 18,1,1 *EL FILE S, E, *NODE FILE RF, NT,RFL, *CONTACT PRINT,MAIN=ASURF,SECONDARY=BSURF CSTRESS,CDISP HFL,HFLA,HTL, SFDR,SFDRA, *CONTACT FILE,MAIN=ASURF,SECONDARY=BSURF CSTRESS,CDISP HFL,HFLA,HTL SFDR,SFDRA *END STEP ** *** STEP 3. Push top block back into position. ** *STEP, NLGEOM, INC=10 *COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE 0.1,1.,0.1,0.1 *BOUNDARY, OP=NEW BACK,3,3 23, 1, 1, 30, 1, 1, 35, 1, 1, 22, 1, 1, 34, 1, 1, 61, 1, 1, 69, 1, 1, 73, 1, 1, 64, 1, 1, 76, 1, 1, 66, 1, 1, 71, 1, 1, 78, 1, 1, NELEM1, 11, 11, 700 NELEM23, 11, 11, 500 TOP, 2, 2 BOTTOM, 2, 2 1, 1, 1 9, 1, 1 13,1, 1 4,1,1 16,1,1 6,1,1 11,1,1 18,1,1 *EL FILE S, E, *NODE FILE RF, NT, RFL, *CONTACT PRINT,MAIN=ASURF,SECONDARY=BSURF CSTRESS,CDISP HFL,HFLA,HTL, SFDR,SFDRA, *CONTACT FILE,MAIN=ASURF,SECONDARY=BSURF CSTRESS,CDISP HFL,HFLA,HTL SFDR,SFDRA *END STEP ** **** test for linear perturbation *** *STEP,NLGEOM,PERTURBATION *STATIC *BOUNDARY,OP=NEW BACK,3,3 BOTTOM,2,2 TOP ,2,2 , 0.05 1, 1, 1 9, 1, 1 13,1, 1 66,1, 1 71,1, 1 78,1, 1 *ENERGY FILE *EL FILE S, E, *NODE FILE RF, *ENERGY PRINT *EL PRINT ,SUM=NO S, E, TEMP, *NODE PRINT,SUM=NO RF, NT, *CONTACT PRINT,MAIN=ASURF,SECONDARY=BSURF *CONTACT FILE,MAIN=ASURF,SECONDARY=BSURF *END STEP