*Heading
 tail strike with elastic behavior of actruator (element 3) modeled via a VUEL rather than through
 connector elasticity
 The VUEL is defined in file: vuel_axial_connector3d.f
 The results match with those obtained from tail_strike_exp_fric.inp
*parameter
actruator_Emod=5000.0
actruator_RefLength=0.0 
actruator_Mass=1.0e-5
actruator_Damping=0.0
** ----------------------------------------------------------------
** PART INSTANCE: A-1
**
*Node
1, 14.95, 19.62, 0
**Element, type=MASS
**19, 1
*Element, type=ROTARYI
20, 1
*Nset, nset=ASSEMBLY_A-1_A-1-REFPT_
1
*Elset, elset=ASSEMBLY_A-1__I1_MASS_
19
*Elset, elset=ASSEMBLY_A-1__I1_ROTARYI_
20
***mass, elset=ASSEMBLY_A-1__I1_MASS_
**1E-05, 
*rotaryinertia, elset=ASSEMBLY_A-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: AIR_FRAME-1
**
*Node
2, -711, -56, 0
*Element, type=MASS
21, 2
*Element, type=ROTARYI
22, 2
*Nset, nset=ASSEMBLY_AIR_FRAME-1_AIR_FRAME
2
*Nset, nset=ASSEMBLY_AIR_FRAME-1_AIR_FRAME-1-REFPT_
2
*Elset, elset=ASSEMBLY_AIR_FRAME-1__I1_MASS_
21
*Elset, elset=ASSEMBLY_AIR_FRAME-1__I1_ROTARYI_
22
*mass, elset=ASSEMBLY_AIR_FRAME-1__I1_MASS_
699.322, 
*rotaryinertia, elset=ASSEMBLY_AIR_FRAME-1__I1_ROTARYI_
1.5E+08, 1.5E+08, 1.5E+08, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: AIR_FRAME_C-1
**
*Node
3, 15.29, 16.13, 0
*Element, type=MASS
23, 3
*Element, type=ROTARYI
24, 3
*Nset, nset=ASSEMBLY_AIR_FRAME_C-1_AIR_FRAME_C-1-REFPT_
3
*Elset, elset=ASSEMBLY_AIR_FRAME_C-1__I1_MASS_
23
*Elset, elset=ASSEMBLY_AIR_FRAME_C-1__I1_ROTARYI_
24
*mass, elset=ASSEMBLY_AIR_FRAME_C-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_AIR_FRAME_C-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: AIR_FRAME_D-1
**
*Node
4, 0, 0, 0
*Element, type=MASS
25, 4
*Element, type=ROTARYI
26, 4
*Nset, nset=ASSEMBLY_AIR_FRAME_D-1_AIR_FRAME_D-1-REFPT_
4
*Elset, elset=ASSEMBLY_AIR_FRAME_D-1__I1_MASS_
25
*Elset, elset=ASSEMBLY_AIR_FRAME_D-1__I1_ROTARYI_
26
*mass, elset=ASSEMBLY_AIR_FRAME_D-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_AIR_FRAME_D-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: B-1
**
*Node
5, 14.15, 10.24, 0
*Element, type=MASS
27, 5
*Element, type=ROTARYI
28, 5
*Nset, nset=ASSEMBLY_B-1_B
5
*Nset, nset=ASSEMBLY_B-1_B-1-REFPT_
5
*Elset, elset=ASSEMBLY_B-1__I1_MASS_
27
*Elset, elset=ASSEMBLY_B-1__I1_ROTARYI_
28
*mass, elset=ASSEMBLY_B-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_B-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: C-1
**
*Node
134, 15.29, 16.13, 0
*Element, type=MASS
141, 134
*Element, type=ROTARYI
142, 134
*Nset, nset=ASSEMBLY_C-1_C-1-REFPT_
134
*Elset, elset=ASSEMBLY_C-1__I1_MASS_
141
*Elset, elset=ASSEMBLY_C-1__I1_ROTARYI_
142
*mass, elset=ASSEMBLY_C-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_C-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: D-1
**
*Node
135, 0, 0, 0
**Element, type=MASS
**143, 135
*Element, type=ROTARYI
144, 135
*Nset, nset=ASSEMBLY_D-1_D-1-REFPT_
135
*Elset, elset=ASSEMBLY_D-1__I1_MASS_
143
*Elset, elset=ASSEMBLY_D-1__I1_ROTARYI_
144
**mass, elset=ASSEMBLY_D-1__I1_MASS_
**1E-05, 
*rotaryinertia, elset=ASSEMBLY_D-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: E-1
**
*Node
136, 13.69, -3.49, 0
*Element, type=MASS
145, 136
*Element, type=ROTARYI
146, 136
*Nset, nset=ASSEMBLY_E-1_E-1-REFPT_
136
*Elset, elset=ASSEMBLY_E-1__I1_MASS_
145
*Elset, elset=ASSEMBLY_E-1__I1_ROTARYI_
146
*mass, elset=ASSEMBLY_E-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_E-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: F-1
**
*Node
137, 24.72, -12.24, 0
*Element, type=MASS
147, 137
*Element, type=ROTARYI
148, 137
*Nset, nset=ASSEMBLY_F-1_F
137
*Nset, nset=ASSEMBLY_F-1_F-1-REFPT_
137
*Elset, elset=ASSEMBLY_F-1__I1_MASS_
147
*Elset, elset=ASSEMBLY_F-1__I1_ROTARYI_
148
*mass, elset=ASSEMBLY_F-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_F-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: LARGE_CYL-1
**
*Node
704, 18.39745, 8.583988, 0
*Element, type=MASS
1003, 704
*Element, type=ROTARYI
1004, 704
*Nset, nset=ASSEMBLY_LARGE_CYL-1_LARGE_CYL
704
*Nset, nset=ASSEMBLY_LARGE_CYL-1_LARGE_CYL-1-REFPT_
704
*Elset, elset=ASSEMBLY_LARGE_CYL-1__I1_MASS_
1003
*Elset, elset=ASSEMBLY_LARGE_CYL-1__I1_ROTARYI_
1004
*mass, elset=ASSEMBLY_LARGE_CYL-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_LARGE_CYL-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: LARGE_CYL_B-1
**
*Node
705, 14.15, 10.24, 0
*Element, type=MASS
1005, 705
*Element, type=ROTARYI
1006, 705
*Nset, nset=ASSEMBLY_LARGE_CYL_B-1_LARGE_CYL_B-1-REFPT_
705
*Elset, elset=ASSEMBLY_LARGE_CYL_B-1__I1_MASS_
1005
*Elset, elset=ASSEMBLY_LARGE_CYL_B-1__I1_ROTARYI_
1006
*mass, elset=ASSEMBLY_LARGE_CYL_B-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_LARGE_CYL_B-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: SMALLCYL_A-1
**
*Node
955, 14.95, 19.62, 0
*Element, type=MASS
1241, 955
*Element, type=ROTARYI
1242, 955
*Nset, nset=ASSEMBLY_SMALLCYL_A-1_SMALLCYL_A-1-REFPT_
955
*Elset, elset=ASSEMBLY_SMALLCYL_A-1__I1_MASS_
1241
*Elset, elset=ASSEMBLY_SMALLCYL_A-1__I1_ROTARYI_
1242
*mass, elset=ASSEMBLY_SMALLCYL_A-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_SMALLCYL_A-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: SMALLCYL-D-1
**
*Node
954, 0, 0, 0
*Element, type=MASS
1239, 954
*Element, type=ROTARYI
1240, 954
*Nset, nset=ASSEMBLY_SMALLCYL-D-1_SMALLCYL-D-1-REFPT_
954
*Elset, elset=ASSEMBLY_SMALLCYL-D-1__I1_MASS_
1239
*Elset, elset=ASSEMBLY_SMALLCYL-D-1__I1_ROTARYI_
1240
*mass, elset=ASSEMBLY_SMALLCYL-D-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_SMALLCYL-D-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
** ----------------------------------------------------------------
** PART INSTANCE: LARGECYL-E-1
**
*Node
703, 13.69, -3.49, 0
*Element, type=MASS
1001, 703
*Element, type=ROTARYI
1002, 703
*Nset, nset=ASSEMBLY_LARGECYL-E-1_LARGECYL-E-1-REFPT_
703
*Elset, elset=ASSEMBLY_LARGECYL-E-1__I1_MASS_
1001
*Elset, elset=ASSEMBLY_LARGECYL-E-1__I1_ROTARYI_
1002
*mass, elset=ASSEMBLY_LARGECYL-E-1__I1_MASS_
1E-05, 
*rotaryinertia, elset=ASSEMBLY_LARGECYL-E-1__I1_ROTARYI_
1E-05, 1E-05, 1E-05, 0., 0., 0.
*transform, nset="ASSEMBLY__T-DATUM CSYS-8"
0.173648177667034, 0.98480775301219, 0., 0.98480775301219, -0.173648177667034, 0.
*orientation, name="ASSEMBLY_DATUM CSYS-5"
0., 0., 1., 0., 1., 0.
1, 0.
*orientation, name="ASSEMBLY_DATUM CSYS-6"
-0.931691194573694, 0.363251315143998, 0., -0.363251315143998, -0.931691194573694, 0.
1, 0.
*orientation, name="ASSEMBLY_DATUM CSYS-7"
0.795403800364659, -0.606079858076027, 0., 0.606079858076027, 0.795403800364659, 0.
1, 0.
*connectorsection, elset=ASSEMBLY__A-CYL_CNSET_
JOIN, 
*connectorsection, elset=ASSEMBLY__AC_CNSET_
BEAM, 
**-----------------------
*connectorsection, elset=ASSEMBLY__AD_CNSET_, behavior=ACTRUATOR_MOD
AXIAL, 
**----------------------
*connectorsection, elset=ASSEMBLY__AIR_FRAME_C_CNSET_
BEAM, 
*connectorsection, elset=ASSEMBLY__AIR_FRAME_D_CNSET_
BEAM, 
*connectorsection, elset=ASSEMBLY__B_CYL_CNSET_
BEAM, 
*connectorsection, elset=ASSEMBLY__BB_CNSET_, behavior=HINGE_FRIC
JOIN, REVOLUTE
"ASSEMBLY_DATUM CSYS-5", 
*connectorsection, elset=ASSEMBLY__CB_CNSET_
BEAM, 
*connectorsection, elset=ASSEMBLY__CC_CNSET_, behavior=HINGE_FRIC
JOIN, REVOLUTE
"ASSEMBLY_DATUM CSYS-5", 
*connectorsection, elset=ASSEMBLY__D-CYL_CNSET_, behavior=HINGE_FRIC
JOIN, REVOLUTE
"ASSEMBLY_DATUM CSYS-5", 
*connectorsection, elset=ASSEMBLY__DD_CNSET_, behavior=HINGE_FRIC
JOIN, REVOLUTE
"ASSEMBLY_DATUM CSYS-5", 
*connectorsection, elset=ASSEMBLY__DE_CNSET_
BEAM, 
*connectorsection, elset=ASSEMBLY__E-CYL_CNSET_
JOIN, 
*connectorsection, elset=ASSEMBLY__EF_CNSET_
BEAM, 
*connectorsection, elset=ASSEMBLY__LARGECYL_CNSET_
ALIGN, 
*connectorsection, elset=ASSEMBLY__LARGE_CYL_CNSET_
SLIDE-PLANE, 
"ASSEMBLY_DATUM CSYS-6", 
*connectorsection, elset=ASSEMBLY__LARGE_CYL_E_CNSET_, behavior=ACTRUATOR
AXIAL, 
*connectorsection, elset=ASSEMBLY__SMALLCYL_CNSET_
SLIDE-PLANE, ALIGN
"ASSEMBLY_DATUM CSYS-7",
**--------------------------------
**the user element acts in parallel with the connector ACTRUATOR_MOD (element 3)
**it replaces the elasticity of the connector element
*User Element, Nodes=2, Type=VU1001, Properties=4, Coordinates=3, Variables=6
 1, 2, 3
*Element, type=vu1001, elset=ASSEMBLY__AD_CNSET_A
 333, 1, 135
*Uel Property, Elset=ASSEMBLY__AD_CNSET_A
 <actruator_Emod>, <actruator_Damping>, <actruator_RefLength>, <actruator_Mass> 
**------------------------------- 
*Element, type=CONN3D2
1, 1, 955
2, 1, 134
3, 1, 135
4, 2, 3
5, 2, 4
6, 705, 704
7, 5, 705
8, 134, 5
9, 3, 134
10, 135, 954
11, 4, 135
12, 135, 136
13, 136, 703
14, 136, 137
15, 704, 703
16, 704, 136
17, 704, 136
18, 954, 955
**
** Assembly Nsets
**
*Nset, nset="ASSEMBLY__T-DATUM CSYS-8"
137
*Nset, nset=QA_TEST_NSET
134, 137, 704, 703, 705, 2, 955, 954
*Nset, nset=ASSEMBLY__PICKEDSET48
2
*Nset, nset=ASSEMBLY__PICKEDSET50
137
*Nset, nset=ASSEMBLY__PICKEDSET51
1, 5, 704, 135, 137, 705, 2, 4, 136, 134, 3
**
** Assembly Elsets
**
*Elset, elset=ASSEMBLY_CONNECTORS
7, 9, 10, 11
*Elset, elset=ASSEMBLY__A-CYL_CNSET_
1
*Elset, elset=ASSEMBLY__AC_CNSET_
2
**---------------------------------
*Elset, elset=ASSEMBLY__AD_CNSET_
3
**---------------------------------
*Elset, elset=ASSEMBLY__AIR_FRAME_C_CNSET_
4
*Elset, elset=ASSEMBLY__AIR_FRAME_D_CNSET_
5
*Elset, elset=ASSEMBLY__BB_CNSET_
7
*Elset, elset=ASSEMBLY__BC-1_CNSET_
3
*Elset, elset=ASSEMBLY__B_CYL_CNSET_
6
*Elset, elset=ASSEMBLY__CB_CNSET_
8
*Elset, elset=ASSEMBLY__CC_CNSET_
9
*Elset, elset=ASSEMBLY__D-CYL_CNSET_
10
*Elset, elset=ASSEMBLY__DD_CNSET_
11
*Elset, elset=ASSEMBLY__DE_CNSET_
12
*Elset, elset=ASSEMBLY__E-CYL_CNSET_
13
*Elset, elset=ASSEMBLY__EF_CNSET_
14
*Elset, elset=ASSEMBLY__LARGECYL_CNSET_
15
*Elset, elset=ASSEMBLY__LARGE_CYL_CNSET_
16
*Elset, elset=ASSEMBLY__LARGE_CYL_E_CNSET_
17
*Elset, elset=ASSEMBLY__SMALLCYL_CNSET_
18
*connectorbehavior, name=ACTRUATOR
*connectorelasticity, component=1
5000., 
*connectordamping, component=1
1000., 
*connectorstop, component=1
8.95922, 
**----
**connector ACTRUATOR_MOD has only damping and lock behavior defined;
**elasticity is defined via separate user defined element (VUEL) in parallel
*connectorbehavior, name=ACTRUATOR_MOD
*connectordamping, component=1
1000., 
*connectorstop, component=1
8.95922, 
**----
*connectorbehavior, name=HINGE_FRIC
*connectorfriction, predefined
1.2, 1.25, 2.0
*friction
0.12, 
*boundary
ASSEMBLY__PICKEDSET48, 1, 1
ASSEMBLY__PICKEDSET48, 2, 2
ASSEMBLY__PICKEDSET48, 3, 3
ASSEMBLY__PICKEDSET48, 4, 4
ASSEMBLY__PICKEDSET48, 5, 5
ASSEMBLY__PICKEDSET48, 6, 6
*initialconditions, type=ROTATINGVELOCITY
ASSEMBLY__PICKEDSET51, -0.0976907, 0., 0., 0., 
-711., -56., 0., -711., -56., 1., 
*Step, name=Step-1
*dynamic, explicit, directusercontrol
1.0E-05, 0.2
*connectormotion
ASSEMBLY__BC-1_CNSET_, 1
*boundary, op=NEW
ASSEMBLY__PICKEDSET48, 3, 3
ASSEMBLY__PICKEDSET48, 4, 4
ASSEMBLY__PICKEDSET48, 5, 5
*boundary, op=NEW
ASSEMBLY__PICKEDSET50, 1, 1
*output, field, variable=PRESELECT
*output, field
*nodeoutput, nset=QA_TEST_NSET
U, 
*output, history
*elementoutput, elset=ASSEMBLY__LARGE_CYL_E_CNSET_
CTF1, CEF1, CVF1, CP1
*elementoutput, elset=ASSEMBLY_CONNECTORS
CUR1, CASUR1, CSM1, CNM1
*endstep