GitHub - Rumi381/PyCiGen: Generate cohesive elements in ABAQUS mesh at the grain boundaries in 3D

Rumi381 / PyCiGen Public

forked from ngrilli/PyCiGen

Notifications You must be signed in to change notification settings
Fork 0
Star 1

Generate cohesive elements in ABAQUS mesh at the grain boundaries in 3D

Notifications

Name		Name	Last commit message	Last commit date
Latest commit History 3 Commits
Job-1.inp		Job-1.inp
README.txt		README.txt
TripleJunctionTest.cae		TripleJunctionTest.cae
TripleJunctionTest.jnl		TripleJunctionTest.jnl
UEL.for		UEL.for
UMAT.for		UMAT.for
__init__.py		__init__.py
abaqusparser.py		abaqusparser.py
cohelement.py		cohelement.py
element.py		element.py
interfelement.py		interfelement.py
interfnode.py		interfnode.py
main.py		main.py
mesh.py		mesh.py
newmesh.py		newmesh.py
node.py		node.py
twinnode.py		twinnode.py

Repository files navigation

Nicolo Grilli
University of Oxford
AWE project 2020
23 Giugno 2020

Code to generate cohesive elements
at the grain boundaries in 3D

Copy 'Job-1.inp' input file in this folder
Elset must be named 'GRAIN1', 'GRAIN2', et cetera

Run
python3.6 main.py

Old bulk element list from abaqus input file will be in:
'Job-1-bulk-elems.inp'

Old node list from abaqus input file will be in:
'Job-1-node.inp'

Old element sets with the generate keyword will be in:
'Job-1-elset.inp'

New node list will be in file:
'Job-1-node-new.inp'

New bulk element list will be in file:
'Job-1-elems-new.inp'

Interface (zero thickness) element list will be in file:
'Job-1-int-elems.inp'

Maximum triple junctions are handled
A node cannot belong to more than 3 grains
Only C3D8 elements are handled
It is fundamental that nodes at the grain boundary
belongs to bulk elements in two/three grains
and are not already duplicated