** Test Kamal cure model with C3D8RT element  with isootropic shrinkage
** ....................................... 
*Heading
*Physical Constants, absolute zero=-273.150
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=pCure
*Node, nset=nall
  1,  0.0, 0.0, 0.0
  2,  1.0, 0.0, 0.0
  3,  1.0, 1.0, 0.0
  4,  0.0, 1.0, 0.0
  5,  0.0, 0.0, 1.0
  6,  1.0, 0.0, 1.0
  7,  1.0, 1.0, 1.0
  8,  0.0, 1.0, 1.0
*Element, type=C3D8RT
  1, 1,2,3,4,5,6,7,8
*Else,elset=eall
  1
*Nset, nset=nfix1
 1,4,5,8
*Nset, nset=nfix2
 1,2,4,5
*Nset, nset=nfix3
 1,2,3,4
*End Part
**  
** ASSEMBLY
**
*Assembly, name=Assembly
**  
*Instance, name=iCure, part=pCure
 *solid section, material=CURE_MATERIAL, elset=eall
*End Instance
**  
*elset,elset=eall
 iCure.eall
*nset,nset=nall
 iCure.nall
*nset,nset=nfix1
 iCure.nfix1
*nset,nset=nfix2
 iCure.nfix2
*nset,nset=nfix3
 iCure.nfix3
*End Assembly
** 
** MATERIALS
** 
** ......... MATERIAL and Property and Parameter Tables
**
*Material, name=CURE_MATERIAL
*Density
 1.25000e-09,
*Conductivity
 3.20000e-04,
*Specific Heat
 1.86000e+06,
*Elastic, dependencies=1, moduli=INSTANT
 3.000e-03, 4.5000e-01, 2.00000e+01, 0.00000e+00
      19.2, 4.5000e-01, 2.00000e+02, 1.00000e+00
*Viscoelastic, Time=Prony
 0.4,,2.0E-7
 0.3,,4.0
 0.2,,5.0E8
*Trs, Definition=Tabular, Dep=1
 -3.5,  20.0, 0.0
 -3.0, 200.0, 0.0
  0.5,  20.0, 1.0
-16.0, 200.0, 1.0
*expansion,dep=1
 2.E-4,   0.0,  0.0
 2.E-4, 200.0,  0.0
 2.E-4,   0.0,  0.6
 2.E-4, 200.0,  0.6
 5.E-5,   0.0,  0.7
 1.5E-4,200.0,  0.7
 1.E-5,   0.0,  1.0
 1.E-4, 200.0,  1.0
*Cure Kinetics, Definition=Kamal
**  Zi       Ei       mi    ni   bi
 3.5E14, 12027.41,  1.05,  1.3,  0.0
*Cure Shrinkage, type=Iso
 0.014333333333333333
*Cure Heat Generation
 2.5e+07
*User Defined Field, type=specified
 1, DOC
**
** 
** BOUNDARY CONDITIONS
** 
*Boundary
 nfix1, 1,1
 nfix2, 2,2
 nfix3, 3,3
** 
** PREDEFINED FIELDS
** 
** Name: initialTemperature   Type: Temperature
*Initial Conditions, type=TEMPERATURE
 nall, 75.0
*Initial Conditions, type=CURE
 eall, 0.001
*Initial Conditions, type=FIELD, VAR=1
 nall, 0.001
** ----------------------------------------------------------------
** 
** STEP: Cure
** 
*Step, name=Cure, nlgeom=YES, inc=5000
 Curing
*Coupled Temperature-displacement, deltmx=200.
 1.0, 50., 1e-09, 1.0
** 
** OUTPUT REQUESTS
** 
** FIELD OUTPUT: fieldOutput
** 
*Output, field, freq=50
*Node Output
 NT, RF, U
*Element Output, directions=YES
 E, EVOL, HFL, S, TEMP, THE, CURESE
 EE,CE,DOC,DOCR,DDOCRDTEMP
** 
** HISTORY OUTPUT: monitorHistoryOutput
** 
*Output, history
*Element Output,Elset=eall
 DOC
*End Step