*HEADING
 TEMPERATURE LOADING OF COMPOSITE SHELLS, S4R 
 using User Subroutine VUTEMP
 See the Verification Manual:
 Abaqus->Verification->User Subroutines->
 User subroutine verification ->UTEMP, VUTEMP, 
 UFIELD, UMASFL, and UPRESS
*NODE  
1, 300.,0.,0.
17, 300.,0.,150.
401, 900.,0.,0.
417, 900.,0.,150.
*NGEN,NSET=LOADB
1,17
*NGEN,NSET=BND3
401,417
*NGEN,NSET=BND2
1,401,100
*NGEN,NSET=BND4
17,417,100
*NFILL,NSET=ALL
BND2,BND4,4,4
*ELEMENT,TYPE=S4R, ELSET=EALL
1,1,5,105,101
*ELGEN,ELSET=ALELL
1,4,4,1,4,100,4
*ELSET,ELSET=PRINT
1,6,11,16
*SHELL SECTION,COMPOSITE,TEMPERATURE=3, ELSET=ALELL
0.9,5,MAT1
1.9,5,MAT2
1.2,5,MAT3
*MATERIAL,NAME=MAT1
*CONDUCTIVITY
500., 
*ELASTIC
3.0E6,0.3
*EXPANSION
1E-5, 
*Density
7.800
*MATERIAL,NAME=MAT2
*CONDUCTIVITY
1500., 
*ELASTIC
3.0E6,0.3
*EXPANSION
1E-5, 
*Density
7.800
*MATERIAL,NAME=MAT3
*CONDUCTIVITY
700., 
*ELASTIC
3.0E6,0.3
*EXPANSION
1E-5, 
*Density
7.800
***
*Nset,Nset=QA_TEST_Nset
ALL
*Elset,Elset=QA_TEST_Elset
ALELL
*STEP
*Dynamic, Explicit
,1.0
*BOUNDARY
ALL,ENCASTRE
*TEMPERATURE,USER
ALL,
*Output,Field,NUM=20
*Node Output
NT
*Element Output
Temp, S, LE
**QA**
*Output,Field,NUM=1
*Node Output,Nset=QA_TEST_Nset
NT
*Element Output,elset=QA_TEST_Elset
Temp,
*Output,History,var=preselect,time Interval=1e-2
*END STEP