*HEADING
FULL AXISYMMETRIC CONICAL CRACK MODEL
 File name: conicalcrack_fullaxi_rms.inp
*SYSTEM
*NODE,SYSTEM=C
*NGEN,LINE=C,NSET=TIP
*SYSTEM
*NODE,SYSTEM=C
**
**CRACK TIP REGION
*NGEN,NSET=TIPI
*NGEN,NSET=OUTER1
*NGEN,NSET=OUTER2
*NGEN,NSET=OUTER3
*NGEN,NSET=OUTER4
*NGEN,NSET=OUTER5
*NGEN,NSET=OUTER6
*NSET,NSET=OUTER
OUTER1,OUTER2,OUTER3,OUTER4,OUTER5,OUTER6
*NFILL,SINGULAR=1,NSET=JR
TIPI,TIP,10,-50
*NFILL,NSET=JREGION,BIAS=0.75
TIP,OUTER,16,2
**
**SECTION 9
*NGEN,NSET=BOT9
*NGEN,NSET=TOP9
*NFILL,NSET=ALL9
BOT9,TOP9,8,100
**
**SECTION 25
*NGEN,NSET=BOT25
*NGEN,NSET=TOP25,LINE=C
*NFILL,NSET=ALL25,BIAS=0.8
BOT25,TOP25,16,50
**
**SECTION 41
*NGEN,NSET=BOT41
*NGEN,NSET=TOP41,LINE=C
*NFILL,NSET=ALL41,BIAS=0.8
BOT41,TOP41,16,50
**
**SECTION 925
*NFILL,NSET=ALL925
BOT9,BOT25,16,1000
**
**SECTION 2541
*NSET,NSET=BOT2533,GENERATE
*NGEN,NSET=TOP2533,LINE=C
*NSET,NSET=BOT3341,GENERATE
*NGEN,NSET=TOP3341,LINE=C
*NSET,NSET=BOT2541
BOT2533,BOT3341
*NSET,NSET=TOP2541
TOP2533,TOP3341
*NFILL,NSET=ALL2541,BIAS=0.8
BOT2541,TOP2541,16,50
**
**SECTION 4157
*NGEN,NSET=BOT57
*NFILL,NSET=ALL4157
BOT41,BOT57,16,1000
**
**SECTION 57
*NGEN,NSET=TOP57
*NFILL,NSET=ALL57,BIAS=1.0
BOT57,TOP57,16,50
**
**SECTION 5765
*NGEN,NSET=TOP5765
*NFILL,NSET=ALL5765,BIAS=1.0
OUTER6,TOP5765,16,50
**
**SECTION 19
*NGEN,NSET=TOP19
*NFILL,NSET=ALL19
OUTER1,TOP19,16,50
**
**INFINITE ELEMENT REGION
*NGEN,NSET=INF25,LINE=C
*NGEN,NSET=INF2533,LINE=C
*NGEN,NSET=INF3341,LINE=C
*NGEN,NSET=INF41,LINE=C
*NSET,NSET=INF
INF25,INF2533,INF3341,INF41
**
**CRACK TIP REGION ELEMENTS
*ELEMENT,TYPE=CAX8R
*ELGEN,ELSET=RINGS
**
**ELEMENTS SECTION 9
*ELEMENT,TYPE=CAX8R
*ELGEN,ELSET=SECT9
*ELSET,ELSET=TOP9,GENERATE
**
**ELEMENTS SECTION 25
*ELEMENT,TYPE=CAX8R
*ELGEN,ELSET=SECT25
**
**ELEMENTS SECTION 19
*ELEMENT,TYPE=CAX8R
*ELGEN,ELSET=SECT19
**
**ELEMENT SECTION 925
*ELEMENT,TYPE=CAX8R
*ELGEN,ELSET=SECT925
**
**ELEMENT SECTION 2541
*ELEMENT,TYPE=CAX8R
*ELGEN,ELSET=SECT2541
**
**ELEMENT SECTION 41
*ELEMENT, TYPE=CAX8R
*ELGEN,ELSET=SECT41
**
**ELEMENT SECTION 4157
*ELEMENT,TYPE=CAX8R
*ELGEN,ELSET=SECT4157
**
**ELEMENT SECTION 57
*ELEMENT, TYPE=CAX8R
*ELGEN,ELSET=SECT57
**
**ELEMENT SECTION 5765
*ELEMENT,TYPE=CAX8R
*ELGEN,ELSET=SECT5765
**
**INFINITE ELEMENTS
*ELEMENT,TYPE=CINAX5R
*ELGEN,ELSET=INF41
*ELEMENT,TYPE=CINAX5R
*ELGEN,ELSET=INF2541
*ELEMENT,TYPE=CINAX5R
*ELGEN,ELSET=INF25
*ELSET,ELSET=INF
INF41,INF2541,INF25
**
*ELSET,ELSET=E1
RINGS,SECT9,SECT25,SECT41,SECT57,SECT19,SECT925,
SECT2541,SECT4157,SECT5765,INF
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,0.3
*SOLID SECTION, MATERIAL=STEEL, ELSET=E1
**
**ADD BOUNDARY CONDITIONS
*NSET,NSET=R9,GENERATE
*NSET,NSET=R925,GENERATE
*NSET,NSET=R25,GENERATE
*NSET,NSET=YAXIS
R9,R925,R25
*NSET,NSET=L57,GENERATE
*NSET,NSET=L4157,GENERATE
*NSET,NSET=L41,GENERATE
*NSET,NSET=XAXIS
TOP9,L57,L4157,L41
*BOUNDARY
YAXIS,XSYMM
**
**MPC'S TO TIE REDUNDANT NODES
*NSET,NSET=TIP1,GENERATE
*NSET,NSET=TIP2,GENERATE
*MPC
TIE,TIP1,TIP2
*MPC
TIE,57865,57833
*MPC
TIE,57833,65833
*MPC
TIE,1833,9833
**
*STEP
APPLY PRESSURE LOAD
*STATIC
1.0,1.0
*DLOAD
TOP9,P2,10.
*CONTOUR INTEGRAL,CONTOURS=5,OUTPUT=BOTH
TIPI,0.707107,-0.707107
*CONTOUR INTEGRAL,CONTOURS=5,OUTPUT=BOTH,TYPE=K FACTORS
TIPI,0.707107,-0.707107
*CONTOUR INTEGRAL,CONTOURS=5,OUTPUT=BOTH,TYPE=T-STRESS
TIPI,0.707107,-0.707107
*EL PRINT,freq=0
S,
*NODE PRINT,freq=0
U,RF
*ENDSTEP