************************************************************************* 
** NOTE: PLOT CEF FOR ELASTO-PLASTIC FORCE RESPONSE			*
**									*
**									*
*************************************************************************
*HEADING
TEST PLASTIC CONNECTOR BEHAVIOR - 3D - RIGID PLASTICITY
*PREPRINT, MODEL=YES
*************************************************************************
**              SECTION I : TRUSS PLASTICITY (REFERENCE)
*************************************************************************
**
*************** ISOTROPIC TABULAR REFERENCE ****************
*NODE, NSET=REF1
1001,0,2,0
1002,0,3,0
*ELEMENT, TYPE=T3D2, ELSET=REF1
1001,1001,1002
*SOLID SECTION,ELSET=REF1,MATERIAL=TRUSS_ISO_TAB
1.0, 
*MATERIAL,NAME=TRUSS_ISO_TAB
*DENSITY
1e-3,
*ELASTIC
200E3, 0.3
*PLASTIC, HARDENING=ISOTROPIC
*INCLUDE, INPUT=misc_plas_data1.inp
**
************* ISOTROPIC EXPONENTIAL REFERENCE ****************
*NODE, NSET=REF2
2001,1,2,0
2002,1,3,0
*ELEMENT, TYPE=T3D2, ELSET=REF2
2001,2001,2002
*SOLID SECTION,ELSET=REF2,MATERIAL=TRUSS_ISO_EXP
1.0, 
*MATERIAL,NAME=TRUSS_ISO_EXP
*DENSITY
1e-3,
*ELASTIC
200E3, 0.3
*PLASTIC, HARDENING=COMBINED, DATATYPE=PARAMETERS
200, 0, 0
*CYCLIC HARDENING, PARAMETERS 
200.0, 104.9881192, 1.21871
**
****************** KINEMATIC HALF CYCLE REFERENCE  *************
*NODE, NSET=REF3
3001,2,2,0
3002,2,3,0
*ELEMENT, TYPE=T3D2, ELSET=REF3
3001,3001,3002
*SOLID SECTION,ELSET=REF3,MATERIAL=TRUSS_KIN_HC
1.0, 
*MATERIAL,NAME=TRUSS_KIN_HC
*DENSITY
1e-3,
*ELASTIC
200E3, 0.3
*PLASTIC, HARDENING=COMBINED, DATATYPE=HALFCYCLE
*INCLUDE, INPUT=misc_plas_data2.inp
**
****************** KINEMATIC STABILIZED REFERENCE***********
*NODE, NSET=REF4
4001,3,2,0
4002,3,3,0
*ELEMENT, TYPE=T3D2, ELSET=REF4
4001,4001,4002
*SOLID SECTION,ELSET=REF4,MATERIAL=TRUSS_KIN_STAB
1.0, 
*MATERIAL,NAME=TRUSS_KIN_STAB
*DENSITY
1e-3,
*ELASTIC
200E3, 0.3
*PLASTIC, HARDENING=COMBINED, DATATYPE=STAB
*INCLUDE, INPUT=misc_plas_data2.inp
**
****************** KINEMATIC PARAMETER REFERENCE ***********
*NODE, NSET=REF5
5001,4,2,0
5002,4,3,0
*ELEMENT, TYPE=T3D2, ELSET=REF5
5001,5001,5002
*SOLID SECTION,ELSET=REF5,MATERIAL=TRUSS_KIN_PARAM
1.0, 
*MATERIAL,NAME=TRUSS_KIN_PARAM
*DENSITY
1e-3,
*ELASTIC
200E3, 0.3
*PLASTIC, HARDENING=COMBINED, DATATYPE=PARAM
200, 1421.5, 15.008
**
*********************** COMBINED REFERENCE *****************
*NODE, NSET=REF6
6001,5,2,0
6002,5,3,0
*ELEMENT, TYPE=T3D2, ELSET=REF6
6001,6001,6002
*SOLID SECTION,ELSET=REF6,MATERIAL=TRUSS_ISO_KIN
1.0, 
*MATERIAL,NAME=TRUSS_ISO_KIN
*DENSITY
1e-3,
*ELASTIC
200E3, 0.3
*PLASTIC, HARDENING=COMBINED, DATATYPE=PARAM
200, 1421.5, 15.008
*CYCLIC HARDENING, PARAMETERS 
200.0, 104.9881192, 1.21871
**
********************************************************************
**		SECTION II : CONNECTOR MODEL DATA
********************************************************************
**
**********************  TEST 16 *****************************
*NODE, NSET=N16
161,5,0,0
162,5,0,0
*ELEMENT, TYPE=CONN3D2, ELSET=C16
161,161,162
*CONNECTOR SECTION, ELSET=C16, BEHAVIOR=B16
CARTESIAN, CARDAN
GLOB16,
*ORIENTATION, NAME=GLOB16
1.0,0.0,0.0,0.0,1.0,0.0
3,0.0
**
*CONNECTOR BEHAVIOR, NAME=B16
*INCLUDE, INPUT=misc_rig_conn3d.inp
*INCLUDE, INPUT=misc_plas_comp1.inp
*INCLUDE, INPUT=misc_hard_kin_par.inp
*INCLUDE, INPUT=misc_hard_iso_exp.inp
*INCLUDE, INPUT=misc_plas_comp4.inp
*INCLUDE, INPUT=misc_hard_kin_par.inp
*INCLUDE, INPUT=misc_hard_iso_exp.inp
**
*ELEMENT, TYPE=MASS, ELSET=M16
162,162
*MASS, ELSET=M16
1E-3
*ELEMENT, TYPE=ROTARYI, ELSET=R16
163,162
*ROTARY INERTIA, ELSET=R16
1E-4, 1E-4, 1E-4
**********************  TEST 26 *****************************
*NODE, NSET=N26
261,5,-1,0
262,5,-1,0
*ELEMENT, TYPE=CONN3D2, ELSET=C26
261,261,262
*CONNECTOR SECTION, ELSET=C26, BEHAVIOR=B26
CARTESIAN, CARDAN
GLOB26,
*ORIENTATION, NAME=GLOB26
1.0,0.0,0.0,0.0,1.0,0.0
3,0.0
**
*CONNECTOR BEHAVIOR, NAME=B26
*INCLUDE, INPUT=misc_rig_conn3d.inp
*INCLUDE, INPUT=misc_plas_coupled_quad.inp
*INCLUDE, INPUT=misc_hard_iso_exp.inp
*INCLUDE, INPUT=misc_hard_kin_par.inp
**
*ELEMENT, TYPE=MASS, ELSET=M26
262,262
*MASS, ELSET=M26
1E-3
*ELEMENT, TYPE=ROTARYI, ELSET=R26
263,262
*ROTARY INERTIA, ELSET=R26
1E-4, 1E-4, 1E-4
**********************  TEST 36 *****************************
*NODE, NSET=N36
361,5,-2,0
362,5,-2,0
*ELEMENT, TYPE=CONN3D2, ELSET=C36
361,361,362
*CONNECTOR SECTION, ELSET=C36, BEHAVIOR=B36
CARTESIAN, CARDAN
GLOB36,
*ORIENTATION, NAME=GLOB36
1.0,0.0,0.0,0.0,1.0,0.0
3,0.0
**
*CONNECTOR BEHAVIOR, NAME=B36
*INCLUDE, INPUT=misc_rig_conn3d.inp
*INCLUDE, INPUT=misc_plas_coupled_quad.inp
*INCLUDE, INPUT=misc_hard_kin_par.inp
*INCLUDE, INPUT=misc_hard_iso_exp.inp
*INCLUDE, INPUT=misc_plas_comp1.inp
*INCLUDE, INPUT=misc_hard_kin_par.inp
*INCLUDE, INPUT=misc_hard_iso_exp.inp
*INCLUDE, INPUT=misc_plas_comp4.inp
*INCLUDE, INPUT=misc_hard_kin_par.inp
*INCLUDE, INPUT=misc_hard_iso_exp.inp
**
*ELEMENT, TYPE=MASS, ELSET=M36
362,362
*MASS, ELSET=M36
1E-3
*ELEMENT, TYPE=ROTARYI, ELSET=R36
363,362
*ROTARY INERTIA, ELSET=R36
1E-4, 1E-4, 1E-4
**********************  TEST 46 *****************************
*NODE, NSET=N46
461,5,-3,0
462,5,-3,0
*ELEMENT, TYPE=CONN3D2, ELSET=C46
461,461,462
*CONNECTOR SECTION, ELSET=C46, BEHAVIOR=B46
CARTESIAN, CARDAN
GLOB46,
*ORIENTATION, NAME=GLOB46
1.0,0.0,0.0,0.0,1.0,0.0
3,0.0
**
*CONNECTOR BEHAVIOR, NAME=B46
*INCLUDE, INPUT=misc_rig_conn3d.inp
*INCLUDE, INPUT=misc_plas_coupled_lin.inp
*INCLUDE, INPUT=misc_hard_kin_par.inp
*INCLUDE, INPUT=misc_hard_iso_exp.inp
**
*ELEMENT, TYPE=MASS, ELSET=M46
462,462
*MASS, ELSET=M46
1E-3
*ELEMENT, TYPE=ROTARYI, ELSET=R46
463,462
*ROTARY INERTIA, ELSET=R46
1E-4, 1E-4, 1E-4
**********************  TEST 56 *****************************
*NODE, NSET=N56
561,5,-4,0
562,5,-4,0
*ELEMENT, TYPE=CONN3D2, ELSET=C56
561,561,562
*CONNECTOR SECTION, ELSET=C56, BEHAVIOR=B56
CARTESIAN, CARDAN
GLOB56,
*ORIENTATION, NAME=GLOB56
1.0,0.0,0.0,0.0,1.0,0.0
3,0.0
**
*CONNECTOR BEHAVIOR, NAME=B56
*INCLUDE, INPUT=misc_rig_conn3d.inp
*INCLUDE, INPUT=misc_plas_coupled_beta.inp
*INCLUDE, INPUT=misc_hard_kin_par.inp
*INCLUDE, INPUT=misc_hard_iso_exp.inp
**
*ELEMENT, TYPE=MASS, ELSET=M56
562,562
*MASS, ELSET=M56
1E-3
*ELEMENT, TYPE=ROTARYI, ELSET=R56
563,562
*ROTARY INERTIA, ELSET=R56
1E-4, 1E-4, 1E-4
********************************************************************
**		SECTION III : HISTORY DATA - DYNAMIC EXPLICIT TESTS
********************************************************************
**
*STEP, nlgeom=yes
*DYNAMIC, EXPLICIT
,1
**
**
*************** ISOTROPIC TABULAR REFERENCE ****************
*AMPLITUDE, NAME=CYCLE_REF1, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY, AMPLITUDE=CYCLE_REF1
1001,1,6,
1002,1,1,
1002,3,6,
1002,2,2,0.05
*OUTPUT, HISTORY, FREQUENCY=50
*ELEMENT OUTPUT, ELSET=REF1
PE11, S11, E11
************* ISOTROPIC EXPONENTIAL REFERENCE ****************
*AMPLITUDE, NAME=CYCLE_REF2, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY, AMPLITUDE=CYCLE_REF2
2001,1,6,
2002,1,1,
2002,3,6,
2002,2,2,0.05
*OUTPUT, HISTORY, FREQUENCY=50
*ELEMENT OUTPUT, ELSET=REF2
PE11, S11, E11
****************** KINEMATIC HALF CYCLE REFERENCE ************
*AMPLITUDE, NAME=CYCLE_REF3, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY, AMPLITUDE=CYCLE_REF3
3001,1,6,
3002,1,1,
3002,3,6,
3002,2,2,0.05
*OUTPUT, HISTORY, FREQUENCY=50
*ELEMENT OUTPUT, ELSET=REF3
PE11, S11, E11
****************** KINEMATIC STABILIZED REFERENCE ***********
*AMPLITUDE, NAME=CYCLE_REF4, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY, AMPLITUDE=CYCLE_REF4
4001,1,6,
4002,1,1,
4002,3,6,
4002,2,2,0.05
*OUTPUT, HISTORY, FREQUENCY=50
*ELEMENT OUTPUT, ELSET=REF4
PE11, S11, E11
****************** KINEMATIC PARAMETER REFERENCE ************
*AMPLITUDE, NAME=CYCLE_REF5, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY, AMPLITUDE=CYCLE_REF5
5001,1,6,
5002,1,1,
5002,3,6,
5002,2,2,0.05
*OUTPUT, HISTORY, FREQUENCY=50
*ELEMENT OUTPUT, ELSET=REF5
PE11, S11, E11
*********************** COMBINED REFERENCE ******************
*AMPLITUDE, NAME=CYCLE_REF6, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY, AMPLITUDE=CYCLE_REF6
6001,1,6,
6002,1,1,
6002,3,6,
6002,2,2,0.05
*OUTPUT, HISTORY, FREQUENCY=50
*ELEMENT OUTPUT, ELSET=REF6
PE11, S11, E11
********************** TEST 16 *****************************
*AMPLITUDE, NAME=CYCLE_16, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY
161,1,6,
*CONNECTOR MOTION, AMPLITUDE=CYCLE_16
161,1,0.05
161,2,0.05
161,3,0.05
161,4,0.05
161,5,0.05
161,6,0.05
*OUTPUT, HISTORY
*ELEMENT OUTPUT, ELSET=C16
CU, CUE, CUP, CUPEQ, CEF, CRF, CALPHAF
********************** TEST 26 *****************************
*AMPLITUDE, NAME=CYCLE_26, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY
261,1,6,
*CONNECTOR MOTION, AMPLITUDE=CYCLE_26
261,1,0.05
261,2,0.05
261,3,0.05
261,4,0.05
261,5,0.05
261,6,0.05
*OUTPUT, HISTORY
*ELEMENT OUTPUT, ELSET=C26
CU, CUE, CUP, CUPEQ, CEF, CRF, CALPHAF
********************** TEST 36 *****************************
*AMPLITUDE, NAME=CYCLE_36, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY
361,1,6,
*CONNECTOR MOTION, AMPLITUDE=CYCLE_36
361,1,0.05
361,2,0.05
361,3,0.05
361,4,0.05
361,5,0.05
361,6,0.05
*OUTPUT, HISTORY
*ELEMENT OUTPUT, ELSET=C36
CU, CUE, CUP, CUPEQ, CEF, CRF, CALPHAF
********************** TEST 46 *****************************
*AMPLITUDE, NAME=CYCLE_46, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY
461,1,6,
*CONNECTOR MOTION, AMPLITUDE=CYCLE_46
461,1,0.05
461,2,0.05
461,3,0.05
461,4,0.05
461,5,0.05
461,6,0.05
*OUTPUT, HISTORY
*ELEMENT OUTPUT, ELSET=C46
CU, CUE, CUP, CUPEQ, CEF, CRF, CALPHAF
********************** TEST 56 *****************************
*AMPLITUDE, NAME=CYCLE_56, DEFINITION=PERIODIC
1,31.4,0,0
0,1,
*BOUNDARY
561,1,6,
*CONNECTOR MOTION, AMPLITUDE=CYCLE_56
561,1,0.05
561,2,0.05
561,3,0.05
561,4,0.05
561,5,0.05
561,6,0.05
*OUTPUT, HISTORY
*ELEMENT OUTPUT, ELSET=C56
CU, CUE, CUP, CUPEQ, CEF, CRF, CALPHAF
*************************************************************
*END STEP