TeV gamma-ray binaries are a new class of gamma-ray sources. They are rare systems: Only three binaries (B1259-63, LS I+61 303, and LS 5039) show persistent TeV emission. All of them are high-mass X-ray binaries, although the nature of the compact object has been established only in one system (B1259-63). Unexpectedly, two of three TeV binaries (B1259-63 and LS I+61 303) have Be stars, which are non-supergiant, early-type stars with a fast polar wind and a dense equatorial disk. In B1259-63, the wind from a radio pulsar collides with the wind and disk of a Be star. On the other hand, in LS I+61 303, a similar interaction is expected if the compact object is a radio pulsar, while the accretion from the Be disk becomes important otherwise. In this talk, I will report on the results from 3-D SPH simulations of the interactions in two TeV binaries, B1259-63 and LS I+61 303.