We present a numerical study of the GAMMA-X mixing in GaAs/AlAs/GaAs quantum well structures. A GAMMA-X mixing model proposed by Liu [Appl. Phys. Lett. 51, 1019 (1987)] is extended to include the effects of self-consistency and nonzero transverse momentum. In the present model, the coupled Schrodinger equations for GAMMA and X electron envelope wave functions are solved self-consistently with Poisson's equation to calculate the electron transmission probability and wave functions, which lead to the current-voltage (I-V) characteristics of single barrier and double barrier resonant tunneling diode structures. The quantum transmitting boundary method is employed in the model for numerical solution of the coupled Schrodinger equations, which proves to be very stable and efficient, even for large (> 2000 angstrom) structures. The features of GAMMA-X mixing, such as the resonance/antiresonance in the transmission probability and the virtual bound states, are clearly demonstrated. Additional physical features are observed in the transmission probability and the wave functions under applied bias conditions. Our work shows that inclusion of transverse momentum, variable effective mass, and the self-consistent potential is important in the realistic modeling of I-V characteristics for structures exhibiting GAMMA-X coupling.