*HEADING THERMAL STRESS IN A CYLINDRICAL SHELL [S4RT] The end temperatures should be: NT11=100; NT12=125;NT13=150;NT14=175;NT15=200 The Analytical stress value is 171.43 MPa (the /explicit solution is about 170.3 MPa at SNEG and 172.5 at SPOS); see the verification Manual *NODE,NSET=ALLN,SYSTEM=C 1, 0.1,0.0,0.0 2, 0.1,0.0,0.02 3, 0.1,11.25,0.02 4, 0.1,11.25,0.0 *ELEMENT,TYPE=S4RT,ELSET=ALLE 1, 1,2,3,4 *SHELL SECTION,ELSET=ALLE,MATERIAL=STEEL 0.001,5 *Elset,elset=qa_test alle *nset,nset=qa_test alln *MATERIAL,NAME=STEEL *DENSITY 7800., *CONDUCTIVITY 52., *SPECIFIC HEAT 586., *ELASTIC 200.E9, 0.3 *EXPANSION 1.2E-5, *NSET,NSET=ROTATE 3,4 *NSET,NSET=SYMM 1,2,3,4 *TRANSFORM,NSET=ROTATE,TYPE=C 0.,0.,0.,0.,0.,1. *Amplitude,name=ramp 0,0,0.01,1 *STEP *DYNAMIC TEMPERATURE-DISPLACEMENT, Explicit ,0.06 *BOUNDARY,amp=ramp ALLN, 11,11, 100. ALLN, 15,15, 200. *BOUNDARY SYMM, 2 SYMM, 4 SYMM, 6 1, 3 4, 3 1, 5 2, 5 3, 5 4, 5 ***** *OUTPUT,FIELD,NUMBER INTERVAL=4 *ELEMENT OUTPUT 1,2,3,4,5 Temp,S *NODE OUTPUT NT, U **** *OUTPUT,FIELD,NUMBER INTERVAL=1 *ELEMENT OUTPUT, ELSET=QA_TEST MisesMax *NODE OUTPUT,NSET=QA_TEST NT *** *OUTPUT,HISTORY,TIME INTERVAL=0.003 *node output,nset=alln NT, *ENERGY OUTPUT,VAR=ALL *** *ENDSTEP