** Test Kamal cure model with CPE4RT element with orthotropic 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
  2,  1.0, 0.0
  3,  1.0, 1.0
  4,  0.0, 1.0
*Element, type=CPE4RT
  1, 1,2,3,4
*Else,elset=eall
  1
*Nset, nset=nfix1
 1,4
*Nset, nset=nfix2
 1,2
*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
*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_Ortho"
 0.014333333333333333,0.01, 0.007
*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
** 
** 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