We develop a new self-consistent model for simulation of the electron–cyclotron maser interaction in cylindrical structures, where expansion of the fields in transverse eigenmodes cannot be directly applied. Instead of solving the nonhomogeneous equation for the fields as a differential equation, a different approach is followed. First, the Green's function for elementary azimuthal and radial RF current sources is analytically derived by expanding the fields in longitudinal modes. Then, the total generated field is calculated by representing the perturbed electron beam as a sum of elementary RF current sources along the axis with amplitude coefficients that are found from the kinematic quantities of the electrons. The self-consistent stationary solution is found by solving the equations of motion along with the field equation in an iterative procedure. The model is useful for the full-wave simulation of lossy structures, which are frequently found in gyro-devices, such as ceramic-loaded interaction circuits of gyro-traveling-wave tubes and beam tunnels of gyrotron oscillators.