Mathematical osteon model for examining poroelastic behaviors

An extended and reasonable stress boundary condition at an osteon exterior wall is presented to solve the model proposed by Rémond and Naili. The obtained pressure and fluid velocity solutions are used to investigate the osteonal poroelastic behaviors. The following results are obtained. (i) Both the fluid pressure and the velocity amplitudes are proportional to the strain amplitude and the loading frequency. (ii) In the physiological loading state, the key role governing the poroelastic behaviors of the osteon is the strain rate. (iii) At the osteon scale, the pressure is strongly affected by the permeability variations, whereas the fluid velocity is not.