1:
|
initialize the algorithm:
|
2:
|
compute the worst case scenario
|
3:
|
set the material parameters of each nonunion area mesh cell to soft tissue
|
4:
|
compute the von Mises stress for the tibia
|
5:
|
search for the maximum stress value
|
6:
|
define max:=maximum stress value
|
7:
|
compute the initial scenario
|
8:
|
set the material parameters of each nonunion area mesh cell to cortical bone
|
9:
|
compute the von Mises stress for the tibia
|
10:
|
define the stress threshold
|
11:
|
define the stop_criterion
|
12:
|
repeat
|
13:
|
if von Mises stress of nonunion area mesh cell < threshold x max
|
14:
|
set material parameter of each nonunion area mesh cell from cortical bone to soft tissue
|