  In this paper, we are concerned with the global structure of radial positive
solutions of boundary value problem$$\text{div}\big(\phi_{N}(\nabla
v)\big)+\lambda f(|x|, v)=0 \text{in} B(R), v=0 \text{on} \partial B(R),
  $$where $\phi_{N}(y)=\frac{y}{\sqrt{1-|y|^{2}}}, y\in \mathbb{R}^{N}$,
$\lambda$ is a positive parameter, $B(R)=\{x\in \mathbb{R}^{N} :|x|<R\}$, and
$|\cdot|$ denote the Euclidean norm in $\mathbb{R}^{N}$. All results, depending
on the behavior of nonlinear term $f$ near 0, are obtained by using global
bifurcation techniques.
