*Heading
*preprint,model=yes,history=yes
*Node
1, 0., -0.5, 0.25
2, 0., 0.5, 0.25
3, 0., 0.5, 0.
4, 0., -0.5, 0.
5, 2., 0.5, 0.25
6, 2., 0.5, 0.125
7, 2., -0.5, 0.125
8, 2., -0.5, 0.25
9, 2., -0.5, 0.
10, 2., 0.5, 0.
11, 11., 0.5, 0.
12, 11., -0.5, 0.
13, 11., -0.5, -0.25
14, 11., 0.5, -0.25
15, 9., -0.5, 0.
16, 9., -0.5, -0.125
17, 9., 0.5, -0.125
18, 9., 0.5, 0.
19, 9., 0.5, -0.25
20, 9., -0.5, -0.25
21, 6., -0.5, 0.125
22, 6., 0.5, 0.125
23, 5., 0.5, -0.125
24, 5., -0.5, -0.125
25, 0., 0.5, 0.125
26, 0., -0.5, 0.125
27, 1., -0.5, 0.
28, 1., -0.5, 0.25
29, 1., 0.5, 0.
30, 1., 0.5, 0.25
31, 11., -0.5, -0.125
32, 11., 0.5, -0.125
33, 10., 0.5, -0.25
34, 10., 0.5, 0.
35, 10., -0.5, -0.25
36, 10., -0.5, 0.
37, 3., -0.5, 0.125
38, 4., -0.5, 0.125
39, 5., -0.5, 0.125
40, 5., 0.5, 0.125
41, 4., 0.5, 0.125
42, 3., 0.5, 0.125
43, 6., -0.5, -0.125
44, 7., -0.5, -0.125
45, 8., -0.5, -0.125
46, 8., 0.5, -0.125
47, 7., 0.5, -0.125
48, 6., 0.5, -0.125
49, 1., -0.5, 0.125
50, 1., 0.5, 0.125
51, 10., 0.5, -0.125
52, 10., -0.5, -0.125
53, 2., 0.5, 0.125
54, 2., -0.5, 0.125
55, 9., -0.5, -0.125
56, 9., 0.5, -0.125
57, 5.5, 0., 0.125
*Element, type=SC8R, elset=solid
1, 28, 30, 50, 49, 1, 2, 25, 26
2, 49, 50, 29, 27, 26, 25, 3, 4
3, 8, 5, 6, 7, 28, 30, 50, 49
4, 7, 6, 10, 9, 49, 50, 29, 27
5, 34, 36, 52, 51, 11, 12, 31, 32
6, 51, 52, 35, 33, 32, 31, 13, 14
7, 18, 15, 16, 17, 34, 36, 52, 51
8, 17, 16, 20, 19, 51, 52, 35, 33
*Element, type=S4R, elset=shell
9, 53, 54, 37, 42
10, 42, 37, 38, 41
11, 41, 38, 39, 40
12, 40, 39, 21, 22
13, 23, 24, 43, 48
14, 48, 43, 44, 47
15, 47, 44, 45, 46
16, 46, 45, 55, 56
*Elset, elset=Solid_12
1,2
*Elset, elset=SurfSolidA
3,4
*Elset, elset=SurfSolidB
7,8
*Elset, elset=EdgeShellA
9,
*Elset, elset=EdgeShellB
16,
*Elset, elset=SurfShellA
9,10,11,12
*Elset, elset=SurfShellB
13,14,15,16
*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
1, 2, 3, 4, 11, 12, 13, 14, 25, 26, 31, 32
*Nset, nset=refNodeFast
57,
*Nset, nset=loadSubstruct
13,14
*Nset, nset=fix
1,2,3,4,25,26
*Elset, elset=SolidsAndShells, generate
1,16,1
*Shell Section, elset=solid, material=Material-1, Stack direction=2
1., 5
*Shell Section, elset=shell, material=Material-1,shell thickness=Shthck
1.0, 5
*Distribution Table,name=distThickTable
length
*Distribution, Location=Element, Table=distThickTable, Name=Shthck
, 1.
shell, 0.25
** 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
0.2,
**
** 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
3,
CTSHR13,
*El print, elset=Solid_12, Totals=yes
1,2,3,4,5
CTSHR13,
STH,
SF4,
*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
3,
CTSHR13,
*El print, elset=Solid_12, Totals=yes
1,2,3,4,5
CTSHR13,
STH,
SF4,
*El Print, freq=1, elset=SolidsAndShells
S,
STH,
*Node Print, freq=1, nset=loadSubstruct
U,
*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