*HEADING
Cylindrical element testing under traction load(CCL18)
*PREPRINT, ECHO=YES, MODEL=YES, HISTORY=YES
**
**********************************************************
** Parameter definition
**
*PARAMETER
INRAD=1.5
EXRAD=2.0
HEIGHT=1.0
**
MIDHEIGHT=HEIGHT/2
MID=(INRAD+EXRAD)/2
**
***********************************************************
** Model definition
**
*NODE, NSET=SECTION
1, <INRAD>,0,0
2, <EXRAD>,0,0
3, <EXRAD>,0,<HEIGHT>
4, <INRAD>,0,<HEIGHT>
5,<MID>,0,0
6,<EXRAD>,0,<MIDHEIGHT>
7,<MID>,0,<HEIGHT>
8,<INRAD>,0,<MIDHEIGHT>
9,<MID>,0,<MIDHEIGHT>
*NCOPY, CHANGE NUMBER=100, OLD SET=SECTION, NEW SET=REST, SHIFT, MULTIPLE=7
0,0,0
0,0,0,0,0,1,45
*NSET, NSET=ALL
REST,SECTION
**
*NSET, NSET=INT, GENERATE
1,701,100
4,704,100
8,708,100
*NSET, NSET=EX, GENERATE
2,702,100
3,703,100
6,706,100
*NSET, NSET=TOP, GENERATE
3,703,100
4,704,100
7,707,100
*NSET, NSET=BOTTOM, GENERATE
1,701,100
2,702,100
5,705,100
*NSET, NSET=AR1, GENERATE
1,701,100
*NSET, NSET=AR2, GENERATE
2,702,100
*NSET, NSET=AR3, GENERATE
3,703,100
*NSET, NSET=AR4, GENERATE
4,704,100
**
*ELEMENT, TYPE=CCL18, ELSET=LAST
7, 1,2,3,601,602,603,701,702,703,5,6,9,605,606,609,705,706,709
8, 1,3,4,601,603,604,701,703,704,9,7,8,609,607,608,709,707,708
*ELEMENT, TYPE=CCL18, ELSET=FIRST
1, 201,202,203,1,2,3,101,102,103,205,206,209,5,6,9,105,106,109
2, 201,203,204,1,3,4,101,103,104,209,207,208,9,7,8,109,107,108
*ELGEN, ELSET=FIRST3
1,3,200,2
2,3,200,2
*ELSET, ELSET=ALL
FIRST3,LAST
*SOLID SECTION, ELSET=ALL, MATERIAL=MYMAT
*MATERIAL, NAME=MYMAT
*ELASTIC
30E6,0.3
************************************************************
** Boundary Condition
*TRANSFORM, TYPE=C, NSET=ALL
0,0,0, 0,0,1
*ORIENTATION, SYSTEM=C, NAME=CYL
0,0,0, 0,0,1
*BOUNDARY
INT,1,3,0
************************************************************
** History data
**
*STEP, NLGEOM=YES, UNSYMM=YES
*STATIC
0.1,1.0
*DLOAD,orientation=cyl
1,TRVEC4,-5e6,0,1,0
3,TRVEC4,-5e6,0,1,0
5,TRVEC4,-5e6,0,1,0
7,TRVEC4,-5e6,0,1,0
*OUTPUT, FIELD
*NODE OUTPUT, VARIABLE=PRESELECT
*END STEP