In adherent cells, actomyosin contractility is regulated mainly by the RhoA signaling pathway, which can be controlled by optogenetics. To model the mechanochemical coupling in such systems, we introduce a finite element framework based on the discontinuous Galerkin method, which allows us to treat cell doublets, chains of cells and monolayers within the same conceptual framework. While the adherent cell layer is modeled as an actively contracting viscoelastic material on an elastic foundation, different models are considered for the Rho-pathway, starting with a simple linear chain that can be solved analytically and later including direct feedback that can be solved only numerically. Our model predicts signal propagation as a function of coupling strength and viscoelastic time scales and identifies the conditions for optimal cell responses and wave propagation. In general, it provides a systematic understanding of how biochemistry and mechanics simultaneously contribute to the communication of adherent cells.