Untitled

*Heading
** Short Fiber Fragmentation Test (SFFT)
** Glass Fiber - Epoxy Matrix with Cohesive Interface
** Axisymmetric Model
**

*Preprint, echo=NO, model=NO, history=NO, contact=NO

** ----------------------------------------------------------------
** MATERIAL DEFINITIONS
** ----------------------------------------------------------------

*Material, name=GlassFiber
*Elastic
72000., 0.22

*Material, name=Epoxy
*Elastic
3000., 0.35

*Material, name=Cohesive
*Elastic, type=TRACTION
1.E6, 1.E6, 1.E6
*Damage Initiation, criterion=MAXS
30., 15., 15.
*Damage Evolution, type=ENERGY, mixed mode behavior=BK, power=1.0
0.3, 0.3, 0.3

** ----------------------------------------------------------------
** PARTS
** ----------------------------------------------------------------

*Part, name=Fiber, type=AXISYMMETRIC
*Node
1, 0., 0.
2, 0., 1.
3, 0.005, 0.
4, 0.005, 1.

*Element, type=CAX4R
1, 1, 2, 4, 3

*Nset, nset=FiberTop, generate
2, 2, 1
*Nset, nset=FiberBottom, generate
1, 1, 1
*Elset, elset=FiberSet
1

*End Part

** Matrix Part
*Part, name=Matrix, type=AXISYMMETRIC
*Node
11, 0.005, 0.
12, 0.005, 1.
13, 0.05, 0.
14, 0.05, 1.

*Element, type=CAX4R
11, 11, 12, 14, 13

*Nset, nset=MatrixTop, generate
12, 12, 1
*Nset, nset=MatrixBottom, generate
11, 11, 1
*Elset, elset=MatrixSet
11

*End Part

** ----------------------------------------------------------------
** ASSEMBLY
** ----------------------------------------------------------------
*Assembly, name=Assembly
*Instance, name=FiberInst, part=Fiber
*End Instance
*Instance, name=MatrixInst, part=Matrix
*End Instance

** Tie the interface using cohesive behavior
*Surface, type=ELEMENT, name=FiberSurf, internal
FiberInst.FiberSet, S4
*Surface, type=ELEMENT, name=MatrixSurf, internal
MatrixInst.MatrixSet, S2

*Contact Pair, interaction=CohesiveInt
MatrixSurf, FiberSurf

*Surface Interaction, name=CohesiveInt
*Surface Behavior, pressure-overclosure=HARD
*Cohesive Behavior, coupled=YES
*Elastic, type=TRACTION
1.E6, 1.E6, 1.E6
*Damage Initiation, criterion=MAXS
30., 15., 15.
*Damage Evolution, type=ENERGY, mixed mode behavior=BK, power=1.0
0.3, 0.3, 0.3

*End Assembly

** ----------------------------------------------------------------
** STEP: APPLY TENSION
** ----------------------------------------------------------------
*Step, name=Tension, nlgeom=YES
*Static
0.01, 1., 1e-05, 1.

*CLoad
FiberInst.FiberTop, 2, 1.0

** Boundary conditions
*Boundary
** Axis of symmetry
FiberInst.1, 1, 1, 0.
** Fix bottom of matrix
MatrixInst.MatrixBottom, 2, 2, 0.
** Fix bottom of fiber
FiberInst.FiberBottom, 2, 2, 0.

** ----------------------------------------------------------------
** OUTPUT REQUESTS
** ----------------------------------------------------------------
*Output, field, variable=PRESELECT
*Output, history, variable=PRESELECT
*Node Output, nset=FiberTop
RF2, U2
*Contact Output
CSTRESS, CDISP, CSTATUS

*End Step