mof healthy2.inp
para
    fct 2
    tct 2
    fbt 2
    tbt 2
    fr 200
    tw 50;
  
para
    mir 20
    mor 45;
  
para
    ex 50
    ecz 6
    ebz 4
    emz 6;
c osteophyte dimensions
para
    oh [%tbt*0.8]
    ow 1
    ox [sqrt(2)/2*(%mor+%ow)]
    oy [sqrt(2)/2*(%mor+%ow)];
  
merge
sd 1 cy 0 0 0 0 0 1 %tw ;
sd 2 plan 0 0 0 0 0 1;
sd 3 plan 0 0 [%tbt] 0 0 1;
sd 4 plan 0 0 [%tbt+%tct] 0 0 1;
sd 5 plan 0 0 0 1 0 0;
sd 6 plan 0 0 0 0 1 0;
sd 7 sp 0 0 [%tbt+%tct+%fr+1e-3] %fr ;
sd 8 sp 0 0 [%tbt+%tct+%fct+%fr+1e-3] %fr ;
sd 9 sp 0 0 [%tbt+%tct+%fct+%fbt+%fr+1e-3] %fr ;
sd 10 cy 0 0 0 0 0 1 %mir;
sd 11 cy 0 0 0 0 0 1 %mor;
sd 12 plan 0.0000000e+00 0.0000000e+00 2.9000000e+00 0 0 1
sd 13 plan 0.0000000e+00 0.0000000e+00 3.8000001e+00 0 0 1
sd 14 sp 0 0 [%tbt+%tct+0.1*%fct+%fr+1e-3] %fr ;
sd 15 sp 0 0 [%tbt+%tct+0.55*%fct+%fr+1e-3] %fr ;
sd 16 sp 0 0 [%tbt+%tct+%fr] %fr ;
sd 17 plan 0 0 [%tbt+%tct+1e-3] 0 0 1;
sd 18 cy 0 0 0 0 0 1 [0.975*%mor]
sd 19 function -1000 1000 -1000 1000 u; v; %oh*exp(-(u**2/2/%ow + v**2/2/%ow)) ; mx [%ox] my [%oy] mz [%tbt] ; ;

  
c tibia
para
    i1 1
    i2 [%i1+int(0.45*%ex+1.5)]
    i3 [%i2+int(0.45*%ex+1.5)]
    i4 [%i3+int(0.1*%ex+1.5)];
  
para
    j1 1
    j2 [%j1+int(0.45*%ex+1.5)]
    j3 [%j2+int(0.45*%ex+1.5)]
    j4 [%j3+int(0.1*%ex+1.5)];
  
para
    k1 1
    k2 [%k1+%ebz]
    k3 [%k2+int(0.45*%ecz+1.5)]
    k4 [%k3+int(0.45*%ecz+1.5)]
    k5 [%k4+int(0.1*%ecz+1.5)];
  
para
    x1 0
    x2 [%x1+0.45*%tw]
    x3 [%x2+0.45*%tw]
    x4 [%x3+0.1*%tw];
  
para
    y1 0
    y2 [%y1+0.45*%tw]
    y3 [%y2+0.45*%tw]
    y4 [%y3+0.1*%tw];
  
para
    z1 0
    z2 [%z1+%tbt]
    z3 [%z2+0.45*%tct]
    z4 [%z3+0.45*%tct]
    z5 [%z4+0.1*%tct];
  
block %i1 %i2 %i3 %i4; %j1 %j2 %j3 %j4; %k1 %k2 %k3 %k4 %k5;
%x1 %x2 %x3 %x4; %y1 %y2 %y3 %y4; %z1 %z2 %z3 %z4 %z5;
dei 1 2; 1 2;;
dei 2 4; 2 4;;
mbi -2; -4;; xy [sqrt(2)/2*0.55*%tw] [(sqrt(2)/2-1)*0.55*%tw]
mbi -2; -3;; xy [sqrt(2)/2*0.45*%tw] [(sqrt(2)/2-1)*0.45*%tw]
dei ; 3 4; 2 5;
dei 3 4; 1 2; 2 5;
bb 2 2 1 2 4 2 1;
bb 2 2 2 2 3 5 2;
bb 2 2 1 4 2 2 1;
bb 2 2 2 3 2 5 2;
sfi ; -4;;sd 1
sfi -4;;;sd 1
sfi ; -3;;sd 11
sfi -3;;;sd 11
sfi -2; 1 2;;sd 10
sfi 1 2; -2;;sd 10
eset 2 1 1 4 4 2 = etbonex
eset 1 2 1 4 4 2 = etboney
eset 2 1 2 3 4 3 = etdcx
eset 2 1 3 3 4 4 = etmcx
eset 2 1 4 3 4 5 = etscx
eset 1 2 2 4 3 3 = etdcy
eset 1 2 3 4 3 4 = etmcy
eset 1 2 4 4 3 5 = etscy
nset 1 1 1 1 4 5 = ntxm
nset 1 1 1 4 1 5 = ntym
nset 1 1 1 4 4 1 = ntzm
fset 1 1 1 4 4 1 = ftzm
nset 1 1 5 4 4 5 = ntzp
fset 1 1 5 4 4 5 = ftzp
nset 2 1 1 2 2 5 = ntfix
nset 1 2 1 2 2 5 or ntfix
nset 1 3 2 4 3 5 = ntout
fset 1 3 2 4 3 5 = ftout
nset 3 1 2 3 2 5 or ntout
fset 3 1 2 3 2 5 or ftout
interrupt

c femur
para
    i1 1
    i2 [%i1+int(0.45*%ex+1.5)]
    i3 [%i2+int(0.45*%ex+1.5)]
    i4 [%i3+int(0.1*%ex+1.5)];
  
para
    j1 1
    j2 [%j1+int(0.45*%ex+1.5)]
    j3 [%j2+int(0.45*%ex+1.5)]
    j4 [%j3+int(0.1*%ex+1.5)];
  
para
    k1 1
    k2 [%k1+int(0.1*%ecz+1.5)]
    k3 [%k2+int(0.45*%ecz+1.5)]
    k4 [%k3+int(0.45*%ecz+1.5)]
    k5 [%k4+%ebz];
  
para
    x1 0
    x2 [%x1+0.45*%tw]
    x3 [%x2+0.45*%tw]
    x4 [%x3+0.1*%tw];
  
para
    y1 0
    y2 [%y1+0.45*%tw]
    y3 [%y2+0.45*%tw]
    y4 [%y3+0.1*%tw];
  
para
    z1 [%tbt+%tct+1e-3]
    z2 [%z1+0.1*%fct+1e-3]
    z3 [%z2+0.45*%fct+1e-3]
    z4 [%z3+0.45*%fct+1e-3]
    z5 [%z4+%fbt+1e-3];
block %i1 %i2 %i3 %i4; %j1 %j2 %j3 %j4; %k1 %k2 %k3 %k4 %k5;
%x1 %x2 %x3 %x4; %y1 %y2 %y3 %y4; %z1 %z2 %z3 %z4 %z5;
dei 1 2; 1 2;;
dei 2 4; 2 4;;
mbi -2; -4;; xy [sqrt(2)/2*0.55*%tw] [(sqrt(2)/2-1)*0.55*%tw]
mbi -2; -3;; xy [sqrt(2)/2*0.45*%tw] [(sqrt(2)/2-1)*0.45*%tw]
dei 1 2; 3 4; 1 4;
dei 3 4; 1 2; 1 4;
bb 2 2 4 2 4 5 3;
bb 2 2 1 2 3 4 4;
bb 2 2 4 4 2 5 3;
bb 2 2 1 3 2 4 4;
sfi -2; 1 2;;sd 10
sfi 1 2; -2;;sd 10
sfi -3;;;sd 11
sfi ; -3;;sd 11
sfi ; -4;;sd 1
sfi -4;;;sd 1
sfi -2; 1 2;;sd 10
sfi 1 2; -2;;sd 10
sfi -3;;;sd 11
sfi ; -3;;sd 11
sfi ; -4;;sd 1
sfi -4;;;sd 1
sfi ;; -5;sd 9
sfi ;; -4;sd 8
sfi ;; -3;sd 15
sfi ;; -2;sd 14
sfi ;; -1;sd 7
eset 2 1 4 4 4 5 = efbonex
eset 1 2 4 4 4 5 = efboney
eset 2 1 3 3 4 4 = efdcx
eset 2 1 2 3 4 3 = efmcx
eset 2 1 1 3 4 2 = efscx
eset 1 2 3 4 3 4 = efdcy
eset 1 2 2 4 3 3 = efmcy
eset 1 2 1 4 3 2 = efscy
nset 1 1 1 1 4 5 = nfxm
nset 1 1 1 4 1 5 = nfym
nset 1 1 1 4 4 1 = nfzm
fset 1 1 1 4 4 1 = ffzm
nset 1 1 5 4 4 5 = nfzp
fset 1 1 5 4 4 5 = ffzp
nset 2 1 1 2 2 5 = nffix
nset 1 2 1 2 2 5 or nffix
nset 1 3 2 4 3 5 = nfout
fset 1 3 2 4 3 5 = ffout
nset 3 1 2 3 2 5 or nfout
fset 3 1 2 3 2 5 or ffout
c meniscus
para
    i1 1
    i2 [%i1+int(%ex/2+0.5)]
    i3 [%i2+int(%ex*(%mor-%mir)/%tw+0.5)];
  
para
    j1 1
    j2 [%j1+int(%ex/2+0.5)]
    j3 [%j2+int(%ex*(%mor-%mir)/%tw+0.5)];
  
para
    k1 1
    k2 [%k1+%emz];
  
para
    x1 0
    x2 [%x1+%mir]
    x3 [%x2+%mor];
  
para
    y1 0
    y2 [%y1+%mir]
    y3 [%y2+%mor];
  
para
    z1 [%tbt+%tct+1e-3]
    z2 [%z1+1];
block %i1 %i2 %i3; %j1 %j2 %j3; %k1 %k2;
%x1 %x2 %x3; %y1 %y2 %y3; %z1 %z2;
dei 2 3; 2 3;;
dei 1 2; 1 2;;
mbi -2; -3;; xy [sqrt(2)/2*%mor] [(sqrt(2)/2-1)*%mor]
bb 2 2 1 2 3 2 5;
bb 2 2 1 3 2 2 5;
sfi -1;;;sd 5
sfi ; -1;;sd 6
sfi -2; 1 2;;sd 10
sfi 1 2; -2;;sd 10
sfi ; -3;;sd 18
sfi -3;;;sd 18
sfi ;; -1;sd 17
sfi ;; -2;sd 16
nset 1 2 1 2 2 2 = nmfix
nset 2 1 1 2 2 2 or nmfix
nset 1 1 1 1 3 2 = nmxm
nset 1 1 1 3 1 2 = nmym
nset 1 1 1 3 3 1 = nmzm
fset 1 1 1 3 3 1 = fmzm
nset 1 1 2 3 3 2 = nmzp
fset 1 1 2 3 3 2 = fmzp
nset 3 1 1 3 3 2 = nmout
nset 1 3 1 3 3 2 or nmout
eset 1 2 1 2 3 2 = emeniy
eset 2 1 1 3 2 2 = emenix
merge
stp 1e-14
abaqus write