Numerical studies of the flux creep in superconductors show that the distri
bution of the magnetic field at any stage of the creep process can be well
described by the condition of spatial constancy of the activation energy U
independently on the particular dependence of U on the field B and current
j. This results from a self-organization of the creep process in the underc
ritical state j<j(c) related to a strong nonlinearity of the flux motion. U
sing the spatial constancy of U, one can find the field profiles B(x), form
ulate a semianalytical approach to the creep problem and generalize the log
arithmic solution for flux creep, obtained for U = U(j), to the case of ess
ential dependence of U on B. This approach is useful for the analysis of dy
namic formation of an anomalous magnetization curve ("fishtail"). We analyz
e the quality of the logarithmic and generalized logarithmic approximations
and show that the latter predicts a maximum in the creep rate at short tim
es, which has been observed experimentally. The vortex annihilation lines (
or the sample edge for the case of remanent state relaxation), where B=0, c
ause instabilities (flux-flow regions) and modify or even destroy the self-
organization of flux creep in the whole sample. [S0163-1829(98)02842-2].