Abstract: An extended and reasonable stress boundary condition at an osteon exterior wall is presented to solve the model proposed by R´emond 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.

Key words: mathematical programs with equilibrium constraints (MPEC), SQP algorithm,successive approximation, global convergence, superlinear, transverse isotropy, pressure, bone fluid flow, osteon, poroelasticity