Why are operators in quantum mechanics always linear?
The straightforward (but perhaps disappointing) answer is that linearity is specified in the postulates/axioms of quantum mechanics.
So if you are doing something that you call "quantum mechanics" then you have already admitted that you are working with linear Hermitian operators that represent observables.
Historically, Bohr worked out the "old quantum theory," which described atoms pretty well and then Schrodinger worked out the standard formalism that we teach nowadays. Effectively, Schrodinger introduced an equation like:
$$
i\frac{\partial \Psi(\vec x, t)}{\partial t} = \frac{-\hbar^2}{2m}\frac{\partial^2}{\partial \vec x^2}\Psi(\vec x, t) + V(x)\Psi(\vec x,t)\;,
\tag{A}$$
which worked really well to describe atomic-scale phenomena.
Why did he want to write such an equation in the first place? Probably because he knew that solving such an equation could result in discrete eigenvalues, and could thus explain the discrete spectra of atoms.
Note that Eq. (A) is linear. So, a linear equation for the wavefunction evolution worked really well. Linear equations for the wave evolution continued to work really well for describing other phenomena. So, why fix what ain't broke?
If you want to describe more than one particle, you again have to solve a linear equation:
$$
i\frac{\partial \Psi(\vec x_1,\ldots,\vec x_N t)}{\partial t} = \left(\sum_{i=1}^N\frac{-\hbar^2}{2m}\frac{\partial^2}{\partial \vec x_i^2}\right)\Psi(\vec x_1,\ldots,\vec x_N t) + V(\vec x_1,\ldots,\vec x_N)\Psi(\vec x_1,\ldots,\vec x_N, t)\;,
\tag{B}$$
Also note that although Eq. (A) and Eq. (B) are both linear that doesn't mean that we only care about linear equations. For example, if you try to solve an N-particle Schrodinger equation (where N>1) like Eq. (B), you usually make a bunch of approximations and reduce it down to solving a non-linear single particle equation (e.g., via density functional theory or whatever).