We report on the results from 3-D SPH simulations of TeV binaries with Be stars. Since there is only one TeV binary (B 1259-63) where the nature of the compact companion has been established, we mainly focus on this Be-pulsar system. From simulations of B 1259-63 around periastron, we find that the pulsar wind dominates the Be-star wind and strips off an outer part of the Be-star disk, causing a strongly asymmetric, phase-dependent structure of the circumstellar material around the Be star. Such a large modulation may be detected by optical, IR, and/or UV observations at phases near periastron. We also discuss the results from simulations of another TeV binary LS I+61 303, for which the nature of the compact object is not yet known.