Thermodynamic functions in the vicinity of an ordinary critical point are expected to obey asymptotic scaling laws as $t=\frac{(T\ensuremath{-}{T}_{c})}{{T}_{c}}$ and the ordering field,
h, approach zero. However, the optimal scaling variables are the nonlinear scaling fields, ${g}_{t}=t+{b}_{t}{h}^{2}+{c}_{t}{t}^{2}+\ensuremath{\cdots}$ and ${g}_{h}=h(1+{c}_{h}t+\ensuremath{\cdots})$. The nonlinearities yield correction factors to the leading power-law (and scaling) variation of thermodynamic quantities,
L, of the form ($1+{a}_{L}t+{b}_{L}{t}^{2}+\ensuremath{\cdots}$), etc., where the correction amplitudes
aL,bL,…, are uniquely determined by the nonlinear scaling-field coefficients. It follows that "analytic" corrections to, e.g., the susceptibility, are directly related to those for the free energy and magnetization (in zero field). The term
bth2 also generates nonanalytic contributions such as an additive, energylike term, varying as ${|t|}^{1\ensuremath{-}\ensuremath{\alpha}}$ in the zero-field susceptibility, and factors like ($1+{c}_{L}{|h|}^{2\ensuremath{-}\frac{1}{\ensuremath{\Delta}}}$) on the critical isotherm,
t=0. Irrelevant scaling fields yield further, in general distinct, nonanalytic corrections, and cause shifts in
Tc and various amplitudes although "universal" ratios remain constant.
A fully anisotropic simple-cubic Ising lattice in the geometry of periodic cylinders n\ifmmode\times\else\texttimes\fi{}n\ifmmode\times\else\texttimes\fi{}\ensuremath{\infty} is investigated by the transfer-matrix finite-size-scaling method. In addition to the previously obtained critical amplitudes of the inverse correlation lengths and singular part of the free energy [M. A. Yurishchev, Phys. Rev. B 50, 13 533 (1994)], the amplitudes of the usual (``linear'') and nonlinear susceptibilities and the amplitude of the second derivative of the spin-spin inverse correlation length with respect to the external field are calculated. The behavior of critical amplitude combinations (which, in accordance with the Privman-Fisher equations, do not contain in their composition the nonuniversal metric coefficients and geometry prefactor) are studied as a function of the interaction anisotropy parameters. A universality domain for the amplitude ratios is found in the quasi-one-dimensional regime of interactions in the system. In the case of a fully isotropic three-dimensional Ising model for which the high precision values of the critical coupling and critical-point free energy are available, improved estimates are obtained for the following four universal quantities: (1) the amplitude of spin-spin inverse correlation length, (2) the amplitude of singular part of the free energy, (3) the ratio of the amplitude of a second derivative of the spin-spin inverse correlation length with respect to the external field to the usual susceptibility amplitude, and (4) the ratio of the nonlinear susceptibility amplitude to the square of the linear susceptibility amplitude (i.e., for the finite-size counterpart of the four-point renormalized coupling constant).
It is argued that the leading corrections to scaling for planar Ising models, which occur as analytic factors, arise from the quadratic terms of the nonlinear thermal and ordering fields (rather than from irrelevant variables). This yields ${{\ensuremath{\alpha}}_{\mathrm{eff}}}^{\ensuremath{'}}+2{\ensuremath{\beta}}_{\mathrm{eff}}+{{\ensuremath{\gamma}}_{\mathrm{eff}}}^{\ensuremath{'}}=2$ for the effective critical exponents, with
no leading corrections, and is confirmed by exact square-lattice results for arbitrary anisotropy
J1J2. For isotropic lattices the ratios of correction amplitudes and quadratic nonlinear-field terms appear universal.