We present skeleton-motivated evaluation of QCD observables. The approach can be applied in analytic versions of QCD in certain classes of renormalization schemes. We present two versions of analytic QCD which can be regarded as low-energy modifications of the "minimal" analytic QCD and which reproduce the measured value of the semihadronic τ decay ratio rτ . Further, we describe an approach of calculating the higher order analytic couplings A k (k = 2, 3, . . .) on the basis of logarithmic derivatives of the analytic coupling A1(Q 2 ). This approach can be applied in any version of analytic QCD. We adjust the free parameters of the afore-mentioned two analytic models in such a way that the skeleton-motivated evaluation reproduces the correct known values of rτ and of the Bjorken polarized sum rule (BjPSR) d b (Q 2 ) at a given point (e.g., at Q 2 = 2 GeV 2 ). We then evaluate the low-energy behavior of the Adler function dv(Q 2 ) and the BjPSR d b (Q 2 ) in the afore-mentioned evaluation approach, in the three analytic versions of QCD. We compare with the results obtained in the "minimal" analytic QCD and with the evaluation approach of Milton et al. and Shirkov. Changes in v3: the values of parameters of analytic QCD models M1 and M2 were refined and the numerical results modified accordingly; the penultimate paragraph of Sec. II and the ultimate paragraph of Sec. III are new; discussion of Figs. 4 was extended; new references were added.