mof healthy.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;
  
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;

  
c tibia
para
    i1 1
    i2 [%i1+int(%ex/2+0.5)]
    i3 [%i2+int(%ex/2+0.5)];
  
para
    j1 1
    j2 [%j1+int(%ex/2+0.5)]
    j3 [%j2+int(%ex/2+0.5)];
  
para
    k1 1
    k2 [%k1+%ebz]
    k3 [%k2+int(0.45*%ecz+0.5)]
    k4 [%k3+int(0.45*%ecz+0.5)]
    k5 [%k4+int(0.1*%ecz+0.5)];
  
para
    x1 0
    x2 [%x1+%tw/2]
    x3 [%x2+%tw/2];
  
para
    y1 0
    y2 [%y1+%tw/2]
    y3 [%y2+%tw/2];
  
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; %j1 %j2 %j3; %k1 %k2 %k3 %k4 %k5;
%x1 %x2 %x3; %y1 %y2 %y3; %z1 %z2 %z3 %z4 %z5;
dei 2 3; 2 3;;
dei 1 2; 1 2;;
mbi -2; -3;; xy [sqrt(2)/2*(%tw/2)] [(sqrt(2)/2-1)*(%tw/2)]
bb 2 2 1 2 3 5 1;
bb 2 2 1 3 2 5 1;
sfi -1;;;sd 5
sfi ; -1;;sd 6
sfi ; -3;;sd 1
sfi -3;;;sd 1
sfi ;; -1;sd 2
sfi ;; -2;sd 3
sfi ;; -5;sd 4
sfi ;; -3;sd 12
sfi ;; -4;sd 13
sfi -2; 1 2;;sd 10
sfi 1 2; -2;;sd 10
unifm 1 2 1 2 3 1 & 1 2 5 2 3 5 & 2 1 5 3 2 5 & 2 1 1 3 2 1 20 0 1
unifm 1 2 1 2 3 5 & 2 1 1 3 2 5 20 0 1
nset 1 1 1 1 3 5 = ntxm
nset 1 1 1 3 1 5 = ntym
nset 1 1 1 3 3 1 = ntzm
fset 1 1 1 3 3 1 = ftzm
nset 1 1 5 3 3 5 = ntzp
nset 2 1 1 2 2 5 = ntfix
nset 1 2 1 2 2 5 or ntfix
fset 1 1 5 3 3 5 = ftzp
eset 2 1 1 3 2 2 = etbonex
eset 2 1 2 3 2 3 = etdcx
eset 2 1 3 3 2 4 = etmcx
eset 2 1 4 3 2 5 = etscx
eset 1 2 1 2 3 2 = etboney
eset 1 2 2 2 3 3 = etdcy
eset 1 2 3 2 3 4 = etmcy
eset 1 2 4 2 3 5 = etscy
c femur
para
    i1 1
    i2 [%i1+int(%ex/2+0.5)]
    i3 [%i2+int(%ex/2+0.5)];
  
para
    j1 1
    j2 [%j1+int(%ex/2+0.5)]
    j3 [%j2+int(%ex/2+0.5)];
  
para
    k1 1
    k2 [%k1+int(0.1*%ecz+0.5)]
    k3 [%k2+int(0.45*%ecz+0.5)]
    k4 [%k3+int(0.45*%ecz+0.5)]
    k5 [%k4+%ebz];
  
para
    x1 0
    x2 [%x1+%tw/2]
    x3 [%x2+%tw/2];
  
para
    y1 0
    y2 [%y1+%tw/2]
    y3 [%y2+%tw/2];
  
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; %j1 %j2 %j3; %k1 %k2 %k3 %k4 %k5;
%x1 %x2 %x3; %y1 %y2 %y3; %z1 %z2 %z3 %z4 %z5;
dei 2 3; 2 3;;
dei 1 2; 1 2;;
mbi -2; -3;; xy [sqrt(2)/2*(%tw/2)] [(sqrt(2)/2-1)*(%tw/2)]
bb 2 2 1 2 3 5 2;
bb 2 2 1 3 2 5 2;
sfi -1;;;sd 5
sfi ; -1;;sd 6
sfi ; -3;;sd 1
sfi -3;;;sd 1
sfi ;; -5;sd 9
sfi ;; -4;sd 8
sfi ;; -3;sd 15
sfi ;; -2;sd 14
sfi ;; -1;sd 7
sfi -2; 1 2;;sd 10
sfi 1 2; -2;;sd 10
unifm 1 2 1 2 3 1 & 1 2 5 2 3 5 & 2 1 5 3 2 5 & 2 1 1 3 2 1 20 0 1
unifm 1 2 1 2 3 5 & 2 1 1 3 2 5 20 0 1
eset 2 1 4 3 2 5 = efbonex
eset 2 1 3 3 2 4 = efdcx
eset 2 1 2 3 2 3 = efmcx
eset 2 1 1 3 2 2 = efscx
eset 1 2 4 2 3 5 = efboney
eset 1 2 3 2 3 4 = efdcy
eset 1 2 2 2 3 3 = efmcy
eset 1 2 1 2 3 2 = efscy
nset 1 1 1 1 3 5 = nfxm
nset 1 1 1 3 1 5 = nfym
nset 1 1 1 3 3 1 = nfzm
nset 1 1 5 3 3 5 = nfzp
nset 2 1 1 2 2 5 = nffix
nset 1 2 1 2 2 5 or nffix
fset 1 1 1 3 3 1 = ffzm
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 3;
bb 2 2 1 3 2 2 3;
sfi -1;;;sd 5
sfi ; -1;;sd 6
sfi -2; 1 2;;sd 10
sfi 1 2; -2;;sd 10
sfi ; -3;;sd 11
sfi -3;;;sd 11
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