I am modelling an unconfined compression experiment using the triphasic material definition, but seem to be running into some problems. I find that febio crashes when I run the analysis which leads me to believe I have done something very wrong! I have attached my .feb file for comment.
It is a quarter wedge of height 1.5mm and radius 2.75mm constructed following the guidlines from the example unconfined compression analysis in the users manual.
I have then input these initial conditions to satisfy a triphasic material:
1) Effective solute (Na+, Cl-) concentrations are 150mM on the free curved surface.
2) Effective fluid pressure is -0.3654MPa on the free curved surface (Pe=-R*T*(Bath Osm)).
It is a 2 step analysis:
1) Fixed charge density is decreased from 0 to -300mEq/L (t=0 to t=1).
2) Compressive displacement in the z-direction of 0.075mm (t=1 to t=3.5).
When I run this analysis but with a biphasic material, it works well and only crashes when I use the triphasic one. I would very much appreciate some feedback as I must be doing something wrong.
Thanks
Mark
It is a quarter wedge of height 1.5mm and radius 2.75mm constructed following the guidlines from the example unconfined compression analysis in the users manual.
I have then input these initial conditions to satisfy a triphasic material:
1) Effective solute (Na+, Cl-) concentrations are 150mM on the free curved surface.
2) Effective fluid pressure is -0.3654MPa on the free curved surface (Pe=-R*T*(Bath Osm)).
It is a 2 step analysis:
1) Fixed charge density is decreased from 0 to -300mEq/L (t=0 to t=1).
2) Compressive displacement in the z-direction of 0.075mm (t=1 to t=3.5).
When I run this analysis but with a biphasic material, it works well and only crashes when I use the triphasic one. I would very much appreciate some feedback as I must be doing something wrong.
Thanks
Mark
Comment