The authors of Reference 3 formulated poroelastic governing equations by including the charge density in the electrochemistry term, and numerically solved them for cancellous bone 4:
$${\rho _{11}}\ddot {u}+{\rho _{12}}\ddot {U} - N{u_{i,jj}} - \left( {A+N} \right){u_{j,ji}} - Q{U_{j,ji}}+b\left( {{{\dot {u}}_i} - {{\dot {U}}_i}} \right)+{q_s}{\Phi _{,i}}=0$$
1
$${\rho _{22}}{\ddot {U}_i}+{\rho _{12}}{\ddot {u}_i} - {\left( {R{U_{j,j}}+Q{u_{j,j}}} \right)_{,i}} - b\left( {{{\dot {u}}_i} - {{\dot {U}}_i}} \right)+{q_s}{\Phi _{,i}}=0$$
2
The values of the parameters in the above equations are provided in References 5,6 and the calculation steps are detailed in Reference 4. N, R, Q, and A are poroelastic parameters, u and U are displacements of fluid and solid constituents, and the variable \(\text{b}\) is defined by \({\phi ^2}\mu /K\), where \(\mu\) is the fluid viscosity, K is the permeability of bone fluid, ε is the permittivity, and the \(\phi\) is the porosity. G is the constant showing the linear relationship between the electric current and electrical potential due to the applied electrical potential. The mechanical properties of Pt and water were considered for solving the equations.