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

ngrilli / PyCiGen Public

Notifications You must be signed in to change notification settings
Fork 14
Star 28

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

28 stars 14 forks Branches Tags Activity

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