c output file
mof single_cell.inp
c block size
para bel 0.05;
  
c cell, pcm dimensions
para
    cr 0.005
    pr 0.0075;
  
c element edge lengths
para
    eel 0.002
    cel 0.0005
    pel 0.0005;

merge
sd 1 plan [-%bel/2] 0 0 1 0 0;
sd 2 plan [%bel/2] 0 0 1 0 0;
sd 3 plan 0 [-%bel/2] 0 0 1 0;
sd 4 plan 0 [%bel/2] 0 0 1 0;
sd 5 plan 0 0 [-%bel/2] 0 0 1;
sd 6 plan 0 0 [%bel/2] 0 0 1;
  
sd 7 sp 0 0 0 %cr;
sd 8 sp 0 0 0 %pr;
  
sd 100 sds 1 2 3 4 5 6;
c cell/pcm block definition
para i1 1
    i2 [%i1+int((%pr-%cr)/%pel+1.5)]
    i3 [%i2+int(%cr/2/%cel+1.5)]
    i4 [%i3+int(%cr/%cel+1.5)]
    i5 [%i4+int(%cr/2/%cel+1.5)]
    i6 [%i5+int((%pr-%cr)/%pel+1.5)];
para j1 1
    j2 [%j1+int((%pr-%cr)/%pel+1.5)]
    j3 [%j2+int(%cr/2/%cel+1.5)]
    j4 [%j3+int(%cr/%cel+1.5)]
    j5 [%j4+int(%cr/2/%cel+1.5)]
    j6 [%j5+int((%pr-%cr)/%pel+1.5)];
para k1 1
    k2 [%k1+int((%pr-%cr)/%pel+1.5)]
    k3 [%k2+int(%cr/2/%cel+1.5)]
    k4 [%k3+int(%cr/%cel+1.5)]
    k5 [%k4+int(%cr/2/%cel+1.5)]
    k6 [%k5+int((%pr-%cr)/%pel+1.5)];
para x1 [-%pr]
    x2 [-%cr]
    x3 [-%cr/2]
    x4 [%cr/2]
    x5 [%cr]
    x6 [%pr];
para y1 [-%pr]
    y2 [-%cr]
    y3 [-%cr/2]
    y4 [%cr/2]
    y5 [%cr]
    y6 [%pr];
para z1 [-%pr]
    z2 [-%cr]
    z3 [-%cr/2]
    z4 [%cr/2]
    z5 [%cr]
    z6 [%pr];

block %i1 %i2 %i3 %i4 %i5 %i6; %j1 %j2 %j3 %j4 %j5 %j6; %k1 %k2 %k3 %k4 %k5 %k6;
%x1 %x2 %x3 %x4 %x5 %x6; %y1 %y2 %y3 %y4 %y5 %y6; %z1 %z2 %z3 %z4 %z5 %z6;
dei 1 3 0 4 6; 1 3 0 4 6; 1 3 0 4 6;
dei 1 3 0 4 6; 1 3 0 4 6; 1 6;
dei 1 3 0 4 6; 3 4; 1 3 0 4 6;
dei 3 4; 1 3 0 4 6; 1 3 0 4 6;
bb 1 4 3 3 4 4 1;
bb 4 4 3 6 4 4 2;
bb 4 3 3 6 3 4 3;
bb 1 3 3 3 3 4 4;
bb 1 3 3 3 4 3 5;
bb 4 3 3 6 4 3 6;
bb 4 3 4 6 4 4 7;
bb 1 3 4 3 4 4 8;
sd 9 plan 0 0 0 1 1 0
sd 10 plan 0 0 0 1 -1 0
sd 11 plan 0 0 0 -1 0 1
sd 12 plan 0 0 0 1 0 1
sfi -1 0 -6;;; sd 8
sfi -2 0 -5;;;sd 7
sfi ; -1 0 -6;;sd 8
sfi ; -2 0 -5;;sd 7
sfi ;; -1 0 -6;sd 8
sfi ;; -2 0 -5;sd 7
sfi 1 3; -4;;sd 9
sfi 4 6; -3;;sd 9
sfi 4 6; -4;;sd 10
sfi 1 3; -3;;sd 10
sfi 1 3; 3 4; -3;sd 11
sfi 4 6; 3 4; -4;sd 11
sfi 4 6; 3 4; -3;sd 12
sfi 1 3; 3 4; -4;sd 12
bb 3 4 3 3 6 4 1;
bb 4 4 3 4 6 4 2;
bb 4 1 3 4 3 4 3;
bb 3 1 3 3 3 4 4;
bb 3 3 1 3 4 3 5;
bb 4 3 1 4 4 3 6;
bb 4 3 4 4 4 6 7;
bb 3 3 4 3 4 6 8;
bb 3 4 3 4 6 3 9;
bb 3 4 4 4 6 4 10;
bb 3 1 4 4 3 4 11;
bb 3 1 3 4 3 3 12;
bb 3 4 1 4 4 3 9;
bb 3 4 4 4 4 6 10;
bb 3 3 4 4 3 6 11;
bb 3 3 1 4 3 3 12;
bb 1 3 3 1 4 4 13;
bb 6 3 3 6 4 4 14;
bb 3 1 3 4 1 4 15;
bb 3 6 3 4 6 4 16;
bb 3 3 1 4 4 1 17;
bb 3 3 6 4 4 6 18;
unifm 3 1 3 4 6 4 & 1 3 3 3 4 4 & 4 3 3 6 4 4 & 3 3 1 4 4 3 & 3 3 4 4 4 6 20 0 1
eset 2 2 2 5 5 5 = ecell
eset 1 1 1 6 6 6 = epcm
eset 2 2 2 5 5 5 - epcm
nseti -1 0 -6;;; = npcm
fseti -1 0 -6;;; = fpcm
nseti ; -1 0 -6;; or npcm
fseti ; -1 0 -6;; or fpcm
nseti ;; -1 0 -6; or npcm
fseti ;; -1 0 -6; or fpcm
nseti -2 0 -5;;; = ncell
fseti -2 0 -5;;; = fcell
nseti ; -2 0 -5;; or ncell
fseti ; -2 0 -5;; or fcell
nseti ;; -2 0 -5; or ncell
fseti ;; -2 0 -5; or fcell

c ecm block definition
para i1 1
    i2 [%i1+int((%bel-2*%cr-2*%pr)/2/%eel+1.5)]
    i3 [%i2+int(%cr/%cel+1.5)]
    i4 [%i3+int((%bel-2*%cr-2*%pr)/2/%eel+1.5)];
para j1 1
    j2 [%j1+int((%bel-2*%cr-2*%pr)/2/%eel+1.5)]
    j3 [%j2+int(%cr/%cel+1.5)]
    j4 [%j3+int((%bel-2*%cr-2*%pr)/2/%eel+1.5)];
para k1 1
    k2 [%k1+int((%bel-2*%cr-2*%pr)/2/%eel+1.5)]
    k3 [%k2+int(%cr/%cel+1.5)]
    k4 [%k3+int((%bel-2*%cr-2*%pr)/2/%eel+1.5)];
  
para x1 [-%bel/2]
    x2 [-(%cr+%pr)/2]
    x3 [(%cr+%pr)/2]
    x4 [%bel/2];
  
para y1 [-%bel/2]
    y2 [-(%cr+%pr)/2]
    y3 [(%cr+%pr)/2]
    y4 [%bel/2];
  
para z1 [-%bel/2]
    z2 [-(%cr+%pr)/2]
    z3 [(%cr+%pr)/2]
    z4 [%bel/2];
  
block %i1 %i2 %i3 %i4; %j1 %j2 %j3 %j4; %k1 %k2 %k3 %k4;
%x1 %x2 %x3 %x4; %y1 %y2 %y3 %y4; %z1 %z2 %z3 %z4;
  
dei 2 3; 2 3; 2 3;
sfi -1 0 -4;;;sd 100
sfi ; -1 0 -4;;sd 100
sfi ;; -1 0 -4;sd 100
bb 2 2 2 2 3 3 13;
bb 3 2 2 3 3 3 14;
bb 2 2 2 3 2 3 15;
bb 2 3 2 3 3 3 16;
bb 2 2 2 3 3 2 17;
bb 2 2 3 3 3 3 18;

eset 1 1 1 4 4 4 = eecm
nseti -1 0 -4;;; = nout
fseti -1 0 -4;;; = fout
nseti ; -1 0 -4;; or nout
fseti ; -1 0 -4;; or fout
nseti ;; -1 0 -4; or nout
fseti ;; -1 0 -4; or fout
merge
stp 1e-14
abaqus write