In this work, we propose a computational protocol enabling the simulation of mechanochromic responses in dye-dispersed polymer blends. The main objective is the modeling of the molecular-level structural changes responsible for the modulation of the photophysical properties that lead to the mechanochromic phenomenon. In this demonstrative study, we focus on predicting the changes in optical absorption displayed by a model system consisting of a dimer of a tetraphenylethylene derivative dispersed in a polyethylene matrix. The blend is subjected to an external stimulus that causes a modulation of the polymer matrix density that translates, in turn, into the emergence of specific mechanical constraints on the optically active dimers. The accurate description of this phenomenon requires the reliable sampling of the dimer configurations induced by the interaction with the matrix under stress. These molecular geometries are associated with modified electronic structures that confer novel absorption responses to the dispersed dyes. In the present contribution, the sampling of these structures is achieved through classical molecular dynamics (MD) simulations including a model element to apply an anisotropic mechanical force. This element allows the microscopic modeling of the chains' and dyes' structural rearrangements under stress. After the sampling, we compare the results of two approaches for the prediction of the optical response: (i) the calculation of a mean response from a statistical average over quantum chemical calculations on the sampled MD structures and (ii) a prediction via a more expensive hybrid scheme allowing the relaxation of the sampled molecular geometries in the presence of the matrix constraints.