*Heading
*preprint,model=yes,history=yes
*Node
*Element, type=SC8R, elset=solid
*Element, type=S4R, elset=shell
*Elset, elset=Solid_12
*Elset, elset=SurfSolidA
*Elset, elset=SurfSolidB
*Elset, elset=EdgeShellA
*Elset, elset=EdgeShellB
*Elset, elset=SurfShellA
*Elset, elset=SurfShellB
*Surface, name=EdgeShellA1
 EdgeShellA, Edge
*Surface, name=EdgeShellB1
 EdgeShellB, Edge
*Surface, name=SurfShellA1
 SurfShellA, SPOS
*Surface, name=SurfShellB1
 SurfShellB, SPOS
*Surface, name=SurfSolidA1
 SurfSolidA,
*Surface, name=SurfSolidB1
 SurfSolidB,
*Nset, nset=RetainedNodes
*Nset, nset=refNodeFast
*Nset, nset=loadSubstruct
*Nset, nset=fix
*Elset, elset=SolidsAndShells, generate
*Shell Section, elset=solid, material=Material-1, Stack direction=2
 1., 5
*Shell Section, elset=shell, material=Material-1
 1.0, 5
*Distribution, Location=Element, Type=Scalar, Name=Shthck
 shell, 0.25
*Element Properties
 Thickness, Shthck
** Constraint: Constraint-1
*Shell to Solid Coupling, constraint name=Constraint-1
EdgeShellA1, SurfSolidA1
** Constraint: Constraint-2
*Shell to Solid Coupling, constraint name=Constraint-2
EdgeShellB1, SurfSolidB1
*Fastener, interaction name=inter-1, reference node set=refNodeFast, property=prop-1
 SurfShellA1, SurfShellB1
*Fastener Property, name=prop-1
**
** MATERIALS
**
*Material, name=Material-1
*Elastic
 3e+07, 0.3
*Restart, write, freq=99
**
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1, Perturbation
 Bend Beam
*Static
*Boundary
 fix,1,3
*Load Case, name = LoadCase_1
*Cload
 loadSubstruct,3,-1.
*End Load Case
*Load Case, name = LoadCase_2
*Cload
 loadSubstruct,3,-2.
*End Load Case
*Load Case, name = LoadCase_3
*Cload
 loadSubstruct,3,-3.
*End Load Case
*Output, field
*Node output
 U,
*Output, history
*Element Output, elset=solid
 CTSHR13,
*Interaction output, name=inter-1
 FTF,
*El print, elset=Solid_12, Totals=yes
 CTSHR13,
 STH,
 SF4,
*Interaction Print, name=inter-1
 FTF,
*End Step
**
** ----------------------------------------------------------------
**
** STEP: Step-2
**
*Step, name=Step-2
 Bend Beam
*Static
 0.25, 1.
*Boundary
 fix,1,3
*Cload
 loadSubstruct,3,-1.
*Output, field
*Node output
 U,
*Output, history
*Node output, nset=loadSubstruct
 U,
*Element Output, elset=solid
 CTSHR13,
*Interaction output, name=inter-1
 FTF,
*Interaction output, nset=refNodeFast
 FTF,
*El print, elset=Solid_12, Totals=yes
 CTSHR13,
 STH,
 SF4,
*El Print, freq=1, elset=SolidsAndShells
 S,
 STH,
*Node Print, freq=1, nset=loadSubstruct
 U,
*Interaction Print
 FTF,
*Interaction Print, nset=refNodeFast
 FTF,
*Interaction Print, name=inter-1
 FTF,
*End Step
**
** ----------------------------------------------------------------
**
** STEP: Step-3
**
*Step, name=Step-3
*Substructure generate, type=z1, overwrite
*Boundary, op=new
*Retained nodal dofs
 RetainedNodes,1,3
*Substructure load case, name=load
*Cload
 loadSubstruct,3,-1.
*End Step