We construct nonassociative quasi-stationary solutions describing deformations of Schwarzschild black holes, BHs, to ellipsoid configurations, which can be black ellipsoids, BEs, and/or BHs with ellipsoidal accretion disks. Such solutions are defined by generic off-diagonal symmetric metrics and nonsymmetric components of metrics (which are zero on base four dimensional, 4-d, Lorentz manifold spacetimes but nontrivial in respective 8-d total (co) tangent bundles). Distorted nonassociative BH and BE solutions are found for effective real sources with terms proportional to hbar kappa (for respective Planck and string constants). These sources and related effective nontrivial cosmological constants are determined by nonlinear symmetries and deformations of the Ricci tensor by nonholonomic star products encoding R-flux contributions from string theory. To generate various classes of (non) associative /commutative distorted solutions we generalize and apply the anholonomic frame and connection deformation method for constructing exact and parametric solutions in modified gravity and/or general relativity theories. We study properties of locally anisotropic relativistic, optically thick, could and thin accretion disks around nonassociative distorted BHs, or BEs, when the effects due to the rotation are negligible. Such configurations describe angular anisotropic deformations of axially symmetric astrophysical models when the nonassociative distortions are related to the outer parts of the accretion disks.