We use Schwarzschild ’ s orbit-superposition technique to construct self-consistent models of the Galactic bar . Using \chi ^ { 2 } minimisation , we find that the best-fit Galactic bar model has a pattern speed \Omega _ { p } = 60 { km s ^ { -1 } kpc ^ { -1 } } , disk mass { M _ { d } = 1.0 \times 10 ^ { 11 } M _ { \odot } } and bar angle \theta _ { bar } = 20 ^ { \circ } for an adopted bar mass { M _ { bar } = 2 \times 10 ^ { 10 } M _ { \odot } } . The model can reproduce not only the three-dimensional and projected density distributions but also velocity and velocity dispersion data from the BRAVA survey . We also predict the proper motions in the range l = [ -12 ^ { \circ } , 12 ^ { \circ } ] , b = [ -10 ^ { \circ } , 10 ^ { \circ } ] , which appear to be higher than observations in the longitudinal direction . The model is stable within a timescale of 0.5 Gyr , but appears to deviate from steady-state on longer timescales . Our model can be further tested by future observations such as those from GAIA .