*HEADING
 Tensile load. Residual thermal stresses and no defects
**
*preprint,model=yes,history=yes
*NODE
**
*ELEMENT, TYPE=C3D8I, ELSET=FLANT
*ELEMENT, TYPE=C3D8I, ELSET=FLANB
** *********
***USER ELEMENT,TYPE=U1,NODES=8,COORDINATES=3,I PROPERTIES=2,
** PROPERTIES=6, UNSYMM, VARIABLES=8
** 1,2,3
** *********
**ELEMENT, TYPE=U1, ELSET=INTERF
*ELEMENT, TYPE=COH3D8, ELSET=INTERF
*ELEMENT, TYPE=C3D8I, ELSET=SKINT
*ELEMENT, TYPE=C3D8I, ELSET=SKINB
*ELEMENT, TYPE=C3D8I, ELSET=FT1
*ELEMENT, TYPE=C3D8I, ELSET=FT2
*ELEMENT, TYPE=C3D8I, ELSET=FT1R
**
**
** *********
*NSET, NSET=CLAMP
*NSET, NSET=DISPL, GENERATE
*NSET, NSET=EXTENS
** *********
*BOUNDARY, OP=NEW
CLAMP, 1,,          0.
CLAMP, 2,,          0.
CLAMP, 3,,          0.
DISPL, 2,,          0.
DISPL, 3,,          0.
**
** DEFINE NODE SET FOR TEMPERATURE FIELDS
**
*NSET, NSET=TFIELD, GENERATE
** *********
***UEL PROPERTY,ELSET=INTERF
**** GIc    GIIc    SI     SII    Kp     Eta (Alpha) NLGEOM, ITMAX
** 0.075,  0.547,  61.0,   68.0, 1.E+6,  1.45,         1,      8
** *********
*ORIENTATION,NAME=ZERO
*ORIENTATION,NAME=NINE
*ORIENTATION,NAME=P45
*ORIENTATION,NAME=M45
*COHESIVE SECTION,elset=INTERF,material=adhesive,THICKNESS=specified,
response=traction separation,controls=ctrs
1.0,
*section control,name=ctrs,viscosity=1.0e-3
** *********
** TOP AND BOTTOM SKIN SUBJECTED TO COMBINED
** THERMAL EXPANSION AND RESIN SHRINKAGE
** *********
*SOLID SECTION,CONTROLS=EHC,ELSET=SKINB,COMPOSITE
0.188,1,CFRPS,ZERO
0.188,1,CFRPS, P45
0.188,1,CFRPS,NINE
0.188,1,CFRPS, M45
0.188,1,CFRPS, P45
0.188,1,CFRPS, M45
0.188,1,CFRPS,ZERO
*SOLID SECTION,CONTROLS=EHC,ELSET=SKINT,COMPOSITE
0.188,1,CFRPS,ZERO
0.188,1,CFRPS, M45
0.188,1,CFRPS, P45
0.188,1,CFRPS, M45
0.188,1,CFRPS,NINE
0.188,1,CFRPS, P45
0.188,1,CFRPS,ZERO
** *********
** FLANGE SUBJECTED TO THERMAL EXPANSION ONLY
** *********
*SOLID SECTION,CONTROLS=EHC,ELSET=FLANB,COMPOSITE
0.188,1,CFRP, P45
0.188,1,CFRP,NINE
0.188,1,CFRP, M45
0.188,1,CFRP,ZERO
0.188,1,CFRP,NINE
*SOLID SECTION,CONTROLS=EHC,ELSET=FLANT,COMPOSITE
0.188,1,CFRP,NINE
0.188,1,CFRP,ZERO
0.188,1,CFRP, M45
0.188,1,CFRP,NINE
0.188,1,CFRP, P45
*SOLID SECTION,CONTROLS=EHC,ELSET=FT1,COMPOSITE,STACK DIRECTION=1
0.188,1,CFRP,NINE
0.188,1,CFRP,ZERO
0.188,1,CFRP, M45
0.188,1,CFRP,NINE
0.188,1,CFRP, P45
*SOLID SECTION,CONTROLS=EHC,ELSET=FT2,COMPOSITE,STACK DIRECTION=2
0.188,1,CFRP,NINE
0.188,1,CFRP,ZERO
0.188,1,CFRP, M45
0.188,1,CFRP,NINE
0.188,1,CFRP, P45
*SOLID SECTION,CONTROLS=EHC,ELSET=FT1R,COMPOSITE,STACK DIRECTION=1
0.188,1,CFRP,NINE
0.188,1,CFRP, P45
** 0.188,1,CFRPD,NINE
** 0.188,1,CFRPD, P45
*section controls,name=ehc,hourglass=enhanced
*MATERIAL,NAME=adhesive
*ELASTIC,type=traction
**<Emod>,<Emod>,<Emod>
1.0E6,1.0E6,1.0E6
*DAMAGE INITIATION,CRITERION=QUADS
**<ultI>, <ultII>,<ultII>
61.0,68.0,68.0
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=1.45
**<GIc>,<GIIc>,<GIIc>
0.075,  0.547,0.547
**0.075,  0.60,0.60
***UEL PROPERTY,ELSET=INTERF
**** GIc    GIIc    SI     SII    Kp     Eta (Alpha) NLGEOM, ITMAX
** 0.075,  0.547,  61.0,   68.0, 1.E+6,  1.45,         1,      8
*MATERIAL,NAME=CFRP
*ELASTIC,TYPE=ENGINEERING CONSTANTS
 144.7E+03, 9.65E+03, 9.65E+03, 0.30,0.30,0.45, 5.2E+03,5.2E+03
 3.4E+03
** THERMAL EXPANSION ONLY
*EXPANSION, TYPE=ORTHO, ZERO = 0.0
** 4.248E-6, -6.55E-3, -6.55E-3,
 3.77E-6, -5.81E-3, -5.81E-3,
**
*MATERIAL,NAME=CFRPS
*ELASTIC,TYPE=ENGINEERING CONSTANTS
 144.7E+03, 9.65E+03, 9.65E+03, 0.30,0.30,0.45, 5.2E+03,5.2E+03
 3.4E+03
** THERMAL EXPANSION PLUS RESIN SHRINKAGE
*EXPANSION, TYPE=ORTHO, ZERO = 0.0
 3.77E-6, -5.81E-3, -5.81E-3,
** -2.075E-5, -1.00E-2, -1.00E-2,
** -2.1E-5, -0.93E-2, -0.93E-2,
**
*MATERIAL,NAME=CFRPD
*ELASTIC,TYPE=ENGINEERING CONSTANTS, DEPENDENCIES=3
 144.7E+03, 9.65E+03, 9.65E+03, 0.30,0.30,0.45, 5.2E+03,5.2E+03
 3.4E+03,  0., 0,0,0
 144.7E+03, 9.65E+02, 9.65E+02, 0.05,0.30,0.05, 5.2E+03,5.2E+02
 3.4E+02,  0., 1,0,0
 144.7E+03, 9.65E+03, 9.65E+03, 0.05,0.05,0.05, 5.2E+01,5.2E+03
 3.4E+01,  0., 0,1,0
 144.7E+03, 9.65E+01, 9.65E+01, 0.05,0.05,0.05, 5.2E+01,5.2E+03
 3.4E+01,  0., 1,1,0
**DEPVAR
**3
**USER DEFINED FIELD
**
** THERMAL EXPANSION ONLY
*EXPANSION, TYPE=ORTHO, ZERO = 0.0
** 4.248E-6, -6.55E-3, -6.55E-3,
 3.77E-6, -5.81E-3, -5.81E-3,
** *************************************************************
**
*RESTART,WRITE,FREQ=100
**-----------------------------------------------------------------------------
**
** STEP 1- APPLY TEMPERATURE/SHRINKAGE TO MODEL
**
**-----------------------------------------------------------------------------
*STEP,NLGEOM,INC=100
*STATIC,STABILIZE
 0.5E-5,1.E-5
*TEMPERATURE, OP=NEW
 TFIELD,   1.0
*EL PRINT,FREQUENCY=0
 S
**EL PRINT,ELSET=INTERF,FREQUENCY=1
** SDV5,SDV6,SDV7,SDV8
*NODE PRINT,FREQUENCY=10,TOTALS=YES,NSET=DISPL
 RF
*OUTPUT, HISTORY, FREQ=1
*NODE OUTPUT, NSET=DISPL
U1, RF1
*NODE OUTPUT, NSET=EXTENS
 U1
*energy output,var=all
*OUTPUT, FIELD, FREQ=10
*NODE OUTPUT
 U
*ELEMENT OUTPUT, ELSET=SKINT
 E,S
*END STEP
**
**-----------------------------------------------------------------------------
**
** STEP 2-
**
**-----------------------------------------------------------------------------
*step,nlgeom,inc=2000
*STATIC,STABILIZE
 0.01,0.85,1.e-8,0.01
*BOUNDARY,TYPE=DISP
  DISPL, 1,,       0.85
*CONTROLS,PARAMETERS=FIELD,FIELD=DISPLACEMENT
** .1,1.0
**0.03,1.0
,1.0
*CONTROLS,PARAMETERS=TIME INCREMENTATION
** 1500,1500,,1500,1500
 , , , , , , , 10,
*CONTROLS,PARAMETERS=LINE SEARCH
 4,4,0.25,0.25,0.15
** 6, 5.0, 0.20
*controls, analysis=discontinuous
*EL PRINT,FREQUENCY=0
 S
**EL PRINT,ELSET=INTERF,FREQUENCY=1
** SDV5,SDV6,SDV7,SDV8
*NODE PRINT,FREQUENCY=10,TOTALS=YES,NSET=DISPL
 RF
*OUTPUT, HISTORY, FREQ=1,var=preselect
*NODE OUTPUT, NSET=DISPL
 U1, RF1
*NODE OUTPUT, NSET=EXTENS
 U1
*OUTPUT, FIELD, FREQ=1,var=preselect
*NODE OUTPUT
 U
*element output,elset=interf
s,sdeg
dmicrt
*END STEP