From 7aecc37219dee51cc026dc4f57110a6788a3c908 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Thu, 21 May 2020 15:42:05 +0200 Subject: [PATCH] using ordering as in DAMASK autogenerated by a small script from lattice.f90 --- processing/post/addSchmidfactors.py | 153 +++++++++++++--------------- 1 file changed, 72 insertions(+), 81 deletions(-) diff --git a/processing/post/addSchmidfactors.py b/processing/post/addSchmidfactors.py index b3c4c7500..9b694e3a1 100755 --- a/processing/post/addSchmidfactors.py +++ b/processing/post/addSchmidfactors.py @@ -15,91 +15,82 @@ scriptID = ' '.join([scriptName,damask.version]) slipSystems = { 'fcc': np.array([ - # Slip direction Plane normal - [ 0, 1,-1, 1, 1, 1, ], - [-1, 0, 1, 1, 1, 1, ], - [ 1,-1, 0, 1, 1, 1, ], - [ 0,-1,-1, -1,-1, 1, ], - [ 1, 0, 1, -1,-1, 1, ], - [-1, 1, 0, -1,-1, 1, ], - [ 0,-1, 1, 1,-1,-1, ], - [-1, 0,-1, 1,-1,-1, ], - [ 1, 1, 0, 1,-1,-1, ], - [ 0, 1, 1, -1, 1,-1, ], - [ 1, 0,-1, -1, 1,-1, ], - [-1,-1, 0, -1, 1,-1, ], - ],'f'), + [+0,+1,-1 , +1,+1,+1], + [-1,+0,+1 , +1,+1,+1], + [+1,-1,+0 , +1,+1,+1], + [+0,-1,-1 , -1,-1,+1], + [+1,+0,+1 , -1,-1,+1], + [-1,+1,+0 , -1,-1,+1], + [+0,-1,+1 , +1,-1,-1], + [-1,+0,-1 , +1,-1,-1], + [+1,+1,+0 , +1,-1,-1], + [+0,+1,+1 , -1,+1,-1], + [+1,+0,-1 , -1,+1,-1], + [-1,-1,+0 , -1,+1,-1], + ],'d'), 'bcc': np.array([ - # Slip system <111>{110} - [ 1,-1, 1, 0, 1, 1, ], - [-1,-1, 1, 0, 1, 1, ], - [ 1, 1, 1, 0,-1, 1, ], - [-1, 1, 1, 0,-1, 1, ], - [-1, 1, 1, 1, 0, 1, ], - [-1,-1, 1, 1, 0, 1, ], - [ 1, 1, 1, -1, 0, 1, ], - [ 1,-1, 1, -1, 0, 1, ], - [-1, 1, 1, 1, 1, 0, ], - [-1, 1,-1, 1, 1, 0, ], - [ 1, 1, 1, -1, 1, 0, ], - [ 1, 1,-1, -1, 1, 0, ], - # Slip system <111>{112} - [-1, 1, 1, 2, 1, 1, ], - [ 1, 1, 1, -2, 1, 1, ], - [ 1, 1,-1, 2,-1, 1, ], - [ 1,-1, 1, 2, 1,-1, ], - [ 1,-1, 1, 1, 2, 1, ], - [ 1, 1,-1, -1, 2, 1, ], - [ 1, 1, 1, 1,-2, 1, ], - [-1, 1, 1, 1, 2,-1, ], - [ 1, 1,-1, 1, 1, 2, ], - [ 1,-1, 1, -1, 1, 2, ], - [-1, 1, 1, 1,-1, 2, ], - [ 1, 1, 1, 1, 1,-2, ], - ],'f'), + [+1,-1,+1 , +0,+1,+1], + [-1,-1,+1 , +0,+1,+1], + [+1,+1,+1 , +0,-1,+1], + [-1,+1,+1 , +0,-1,+1], + [-1,+1,+1 , +1,+0,+1], + [-1,-1,+1 , +1,+0,+1], + [+1,+1,+1 , -1,+0,+1], + [+1,-1,+1 , -1,+0,+1], + [-1,+1,+1 , +1,+1,+0], + [-1,+1,-1 , +1,+1,+0], + [+1,+1,+1 , -1,+1,+0], + [+1,+1,-1 , -1,+1,+0], + [-1,+1,+1 , +2,+1,+1], + [+1,+1,+1 , -2,+1,+1], + [+1,+1,-1 , +2,-1,+1], + [+1,-1,+1 , +2,+1,-1], + [+1,-1,+1 , +1,+2,+1], + [+1,+1,-1 , -1,+2,+1], + [+1,+1,+1 , +1,-2,+1], + [-1,+1,+1 , +1,+2,-1], + [+1,+1,-1 , +1,+1,+2], + [+1,-1,+1 , -1,+1,+2], + [-1,+1,+1 , +1,-1,+2], + [+1,+1,+1 , +1,+1,-2], + ],'d'), 'hex': np.array([ - # Basal systems <11.0>{00.1} (independent of c/a-ratio, Bravais notation (4 coordinate base)) - [ 2, -1, -1, 0, 0, 0, 0, 1, ], - [-1, 2, -1, 0, 0, 0, 0, 1, ], - [-1, -1, 2, 0, 0, 0, 0, 1, ], - # 1st type prismatic systems <11.0>{10.0} (independent of c/a-ratio) - [ 2, -1, -1, 0, 0, 1, -1, 0, ], - [-1, 2, -1, 0, -1, 0, 1, 0, ], - [-1, -1, 2, 0, 1, -1, 0, 0, ], - # 2nd type prismatic systems <10.0>{11.0} -- a slip; plane normals independent of c/a-ratio - [ 0, 1, -1, 0, 2, -1, -1, 0, ], - [-1, 0, 1, 0, -1, 2, -1, 0, ], - [ 1, -1, 0, 0, -1, -1, 2, 0, ], - # 1st type 1st order pyramidal systems <11.0>{-11.1} -- plane normals depend on the c/a-ratio - [ 2, -1, -1, 0, 0, 1, -1, 1, ], - [-1, 2, -1, 0, -1, 0, 1, 1, ], - [-1, -1, 2, 0, 1, -1, 0, 1, ], - [ 1, 1, -2, 0, -1, 1, 0, 1, ], - [-2, 1, 1, 0, 0, -1, 1, 1, ], - [ 1, -2, 1, 0, 1, 0, -1, 1, ], - # pyramidal system: c+a slip <11.3>{-10.1} -- plane normals depend on the c/a-ratio - [ 2, -1, -1, 3, -1, 1, 0, 1, ], - [ 1, -2, 1, 3, -1, 1, 0, 1, ], - [-1, -1, 2, 3, 1, 0, -1, 1, ], - [-2, 1, 1, 3, 1, 0, -1, 1, ], - [-1, 2, -1, 3, 0, -1, 1, 1, ], - [ 1, 1, -2, 3, 0, -1, 1, 1, ], - [-2, 1, 1, 3, 1, -1, 0, 1, ], - [-1, 2, -1, 3, 1, -1, 0, 1, ], - [ 1, 1, -2, 3, -1, 0, 1, 1, ], - [ 2, -1, -1, 3, -1, 0, 1, 1, ], - [ 1, -2, 1, 3, 0, 1, -1, 1, ], - [-1, -1, 2, 3, 0, 1, -1, 1, ], - # pyramidal system: c+a slip <11.3>{-1-1.2} -- as for hexagonal ice (Castelnau et al. 1996, similar to twin system found below) - [ 2, -1, -1, 3, -2, 1, 1, 2, ], # sorted according to similar twin system - [-1, 2, -1, 3, 1, -2, 1, 2, ], # <11.3>{-1-1.2} shear = 2((c/a)^2-2)/(3 c/a) - [-1, -1, 2, 3, 1, 1, -2, 2, ], - [-2, 1, 1, 3, 2, -1, -1, 2, ], - [ 1, -2, 1, 3, -1, 2, -1, 2, ], - [ 1, 1, -2, 3, -1, -1, 2, 2, ], - ],'f'), + [+2,-1,-1,+0 , +0,+0,+0,+1], + [-1,+2,-1,+0 , +0,+0,+0,+1], + [-1,-1,+2,+0 , +0,+0,+0,+1], + [+2,-1,-1,+0 , +0,+1,-1,+0], + [-1,+2,-1,+0 , -1,+0,+1,+0], + [-1,-1,+2,+0 , +1,-1,+0,+0], + [-1,+1,+0,+0 , +1,+1,-2,+0], + [+0,-1,+1,+0 , -2,+1,+1,+0], + [+1,+0,-1,+0 , +1,-2,+1,+0], + [-1,+2,-1,+0 , +1,+0,-1,+1], + [-2,+1,+1,+0 , +0,+1,-1,+1], + [-1,-1,+2,+0 , -1,+1,+0,+1], + [+1,-2,+1,+0 , -1,+0,+1,+1], + [+2,-1,-1,+0 , +0,-1,+1,+1], + [+1,+1,-2,+0 , +1,-1,+0,+1], + [-2,+1,+1,+3 , +1,+0,-1,+1], + [-1,-1,+2,+3 , +1,+0,-1,+1], + [-1,-1,+2,+3 , +0,+1,-1,+1], + [+1,-2,+1,+3 , +0,+1,-1,+1], + [+1,-2,+1,+3 , -1,+1,+0,+1], + [+2,-1,-1,+3 , -1,+1,+0,+1], + [+2,-1,-1,+3 , -1,+0,+1,+1], + [+1,+1,-2,+3 , -1,+0,+1,+1], + [+1,+1,-2,+3 , +0,-1,+1,+1], + [-1,+2,-1,+3 , +0,-1,+1,+1], + [-1,+2,-1,+3 , +1,-1,+0,+1], + [-2,+1,+1,+3 , +1,-1,+0,+1], + [-1,-1,+2,+3 , +1,+1,-2,+2], + [+1,-2,+1,+3 , -1,+2,-1,+2], + [+2,-1,-1,+3 , -2,+1,+1,+2], + [+1,+1,-2,+3 , -1,-1,+2,+2], + [-1,+2,-1,+3 , +1,-2,+1,+2], + [-2,+1,+1,+3 , +2,-1,-1,+2], + ],'d'), } # --------------------------------------------------------------------