** Test Kamal cure model with C3D8RT element  with isootropic shrinkage
** ....................................... 
*Heading
*Physical Constants, absolute zero=-273.150, universal gas constant= 8.31434e+00
*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="ABQ_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
** 
** ......... Property and Parameter Types
**
*include,input=ABQ_Cure_types.inp
**
** ......... MATERIAL and Property and Parameter Tables
**
*Material, name="ABQ_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
*Eigenstrain, user
*Parameter Table, type="ABQ_Cure_ReactionKinetics_Kamal"
**  Zi       Ei     mi    ni   bi
 3.5E14, 1.0E5,  1.05,  1.3,  0.0
*Property Table, type="ABQ_Cure_ShrinkageCoeff_Iso"
 0.014333333333333333
*Heat Generation
*Parameter Table, type="ABQ_Cure_HeatGeneration"
 2.5e+07, 1.25000e-09
*User Defined Field, type=specified
 1, SDV1
*Depvar
 3
 1,ALPHA,"Cure state"
 2,ALPHAR,"Rate of curing"
 3,DALPHARDT,"Derivative of rate of curing w.r.t. temperature"
**
** 
** 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=SOLUTION
 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, SDV, TEMP, THE, EEIG
 EE,CE,
** 
** HISTORY OUTPUT: monitorHistoryOutput
** 
*Output, history
*Element Output,Elset=eall
 SDV1
*End Step