B.Lee
Student
- Jun 9, 2022
- 1
Hello everyone, I am a student learning about the cohesive element in Abaqus. I download an example and try to build a model by using the same material property, model scale and boundary condition. However, my result of the model is different and the cohesive element does not fail. Both input files are shown below and attached in the zip file as well:
This is my .inp file.
This is the example input file.
This is my .inp file.
Code:
*Heading
** Job name: Job0607 Model name: Model-1
** Generated by: Abaqus/CAE 2021
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Part-1-1, part=Part-1
*Node
1, 1., 0., 0.50999999
2, 1., 1., 0.50999999
3, 1., 1., 1.00999999
4, 1., 0., 1.00999999
5, 0., 1., 0.50999999
6, 0., 1., 1.00999999
7, 0., 0., 1.00999999
8, 0., 0., 0.50999999
9, 1., 0., 0.5
10, 1., 1., 0.5
11, 0., 0., 0.5
12, 0., 1., 0.5
13, 1., 0., 0.
14, 1., 1., 0.
15, 0., 0., 0.
16, 0., 1., 0.
*Element, type=C3D8R
1, 8, 5, 6, 7, 1, 2, 3, 4
3, 12, 11, 15, 16, 10, 9, 13, 14
*Element, type=COH3D8
2, 5, 8, 11, 12, 2, 1, 9, 10
*Nset, nset=Set-1, generate
1, 8, 1
*Elset, elset=Set-1
1,
*Nset, nset=Set-2, generate
9, 16, 1
*Elset, elset=Set-2
3,
*Nset, nset=Set-3
1, 2, 5, 8, 9, 10, 11, 12
*Elset, elset=Set-3
2,
** Section: Steel_Sec
*Solid Section, elset=Set-1, material=Steel
,
** Section: Coh_Sec
*Cohesive Section, elset=Set-3, controls=EC-1, material=Coh, response=TRACTION SEPARATION
,
** Section: Steel_Sec
*Solid Section, elset=Set-2, material=Steel
,
*End Instance
**
*Nset, nset=Set-2, instance=Part-1-1, generate
13, 16, 1
*Elset, elset=Set-2, instance=Part-1-1
3,
*Nset, nset=Set-3, instance=Part-1-1
3, 4, 6, 7
*Elset, elset=Set-3, instance=Part-1-1
1,
*Elset, elset=_Surf-1_S5, internal, instance=Part-1-1
1,
*Surface, type=ELEMENT, name=Surf-1
_Surf-1_S5, S5
*End Assembly
**
** ELEMENT CONTROLS
**
*Section Controls, name=EC-1, ELEMENT DELETION=YES
1., 1., 1.
**
** MATERIALS
**
*Material, name=Coh
*Damage Initiation, criterion=MAXS
90.,66.,66.
*Damage Evolution, type=DISPLACEMENT
0.04,
*Density
1e-09,
*Elastic, type=TRACTION
4000.,700.,700.
*Material, name=Steel
*Density
7.75e-09,
*Elastic
2.05e+08, 0.3
**
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Displacement/Rotation
*Boundary
Set-2, 1, 1
Set-2, 2, 2
Set-2, 3, 3
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES, inc=10000
*Static
0.001, 0.1, 1e-06, 0.001
**
** BOUNDARY CONDITIONS
**
** Name: BC-2 Type: Displacement/Rotation
*Boundary
Set-3, 1, 1
Set-3, 2, 2
Set-3, 3, 3, 0.1
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-2
**
*Output, field
*Node Output
U,
*Element Output, directions=YES
E, S, SDEG, STATUS
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End StepThis is the example input file.
Code:
*Heading
** Job name: Job-2 Model name: Example_01
** Generated by: Abaqus/CAE 2021
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=BLOCK
*Node
1, 0., 0., 1.
2, 1., 0., 1.
3, 1., 0., 0.
4, 0., 0., 0.
11, 0., 0.49000001, 1.
12, 1., 0.49000001, 1.
13, 1., 0.49000001, 0.
14, 0., 0.49000001, 0.
21, 0., 0.5, 1.
22, 1., 0.5, 1.
23, 1., 0.5, 0.
24, 0., 0.5, 0.
31, 0., 1., 1.
32, 1., 1., 1.
33, 1., 1., 0.
34, 0., 1., 0.
*Element, type=C3D8R
1, 1, 2, 3, 4, 11, 12, 13, 14
2, 21, 22, 23, 24, 31, 32, 33, 34
*Element, type=COH3D8
3, 11, 12, 13, 14, 21, 22, 23, 24
*Nset, nset=TOP, generate
31, 34, 1
*Nset, nset=BOTTOM, generate
1, 4, 1
*Elset, elset=TOPBLOCK
2,
*Elset, elset=BOTBLOCK
1,
*Elset, elset=ADHESIVE
3,
** Section: Section-1-BOTBLOCK
*Solid Section, elset=BOTBLOCK, material=STEEL
,
** Section: Section-2-TOPBLOCK
*Solid Section, elset=TOPBLOCK, material=STEEL
,
** Section: Section-3-ADHESIVE
*Cohesive Section, elset=ADHESIVE, material=MAT, response=TRACTION SEPARATION
,
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=BLOCK, part=BLOCK
*End Instance
**
*End Assembly
**
** MATERIALS
**
** ********************************
*Material, name=MAT
*Damage Initiation, criterion=MAXS
90.,66.,66.
*Damage Evolution, type=DISPLACEMENT
0.04,
*Density
1e-09,
*Elastic, type=TRACTION
4000.,700.,700.
** ********************************
*Material, name=STEEL
*Density
7.75e-09,
*Elastic
2.05e+08, 0.3
**
** BOUNDARY CONDITIONS
**
** Name: Disp-BC-1 Type: Displacement/Rotation
*Boundary
BLOCK.BOTTOM, 1, 1
** Name: Disp-BC-2 Type: Displacement/Rotation
*Boundary
BLOCK.BOTTOM, 2, 2
** Name: Disp-BC-3 Type: Displacement/Rotation
*Boundary
BLOCK.BOTTOM, 3, 3
** ----------------------------------------------------------------
**
** STEP: loadit
**
*Step, name=loadit, nlgeom=YES, inc=1000
*Static
0.001, 0.1, 1e-05, 0.001
**
** BOUNDARY CONDITIONS
**
** Name: Disp-BC-4 Type: Displacement/Rotation
*Boundary
BLOCK.TOP, 1, 1
** Name: Disp-BC-5 Type: Displacement/Rotation
*Boundary
BLOCK.TOP, 2, 2, 0.1
** Name: Disp-BC-6 Type: Displacement/Rotation
*Boundary
BLOCK.TOP, 3, 3
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-2
**
*Output, field
*Element Output, directions=YES
E, S, SDEG, STATUS
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history
*Node Output, nset=BLOCK.TOP
RF2, U2
**
** HISTORY OUTPUT: H-Output-2
**
*Output, history, variable=PRESELECT
*End Step