A wide range of biological applications such as drug delivery, biosensors and hemodialysis can be provided by nanoporous biomaterials due to their uniform pore size as well as considerable pore density. In the current study, the size dependency in the nonlinear primary resonance of micro/nano-beams made of nanoporous biomaterials is anticipated. To accomplish this end, a refined truncated cube is introduced to model the lattice structure of nanoporous biomaterial. Accordingly, analytical expressions for the mechanical properties of material are derived as functions of pore size. After that, based upon a nonlocal strain gradient beam model, the size-dependent nonlinear Duffing type equation of motion is constructed. The Galerkin technique together with the multiple time-scales method is employed to obtain the nonlocal strain gradient frequency-response and amplitude-response related to the nonlinear primary resonance of a micro/nano-beam made of the nanoporous biomaterial with different pore sizes. It is indicated that the nonlocality causes to decrease the response amplitudes associated with the both bifurcation points of the jump phenomenon, while the strain gradient size dependency causes to increase them. Also, it is found that increasing the pore size leads to enhance the nonlinearity, so the maximum deflection of response occurs at higher excitation frequency. Keywords: Biomechanics, Size effect, Nonlinear oscillation, Porosity, Nonlocal strain gradient elasticity