The seasonal variation of Titan's atmospheric structure with emphasis on the stratosphere is simulated by a three-dimensional general circulation model. The model includes the transport of haze particles by the circulation. The likely pattern of meridional circulation is reconstructed by a comparison of simulated and observed haze and temperature distribution. The GCM produces a weak zonal circulation with a small latitudinal temperature gradient, in conflict with observation. The direct reason is found to be the excessive meridional circulation. Under uniformly distributed opacity sources, the model predicts a pair of symmetric Hadley cells near the equinox and a single global cell with the rising branch in the summer hemisphere below about z = 230 km and a thermally indirect cell above the direct cell near the solstice. The interhemispheric circulation transports haze particles from the summer to the winter hemisphere, causing a maximum haze opacity contrast near the solstice and a smaller contrast near the equinox, contrary to observation. On the other, if the GCM is run under modified cooling rate in order to account for the enhancement in nitrites and some hydrocarbons in the northern hemisphere near the vernal equinox, the meridional cell at the equinox becomes a single cell with rising motions in the autumn hemisphere. A more realistic haze opacity distribution can be reproduced at the equinox. However, a pure transport effect (without particle growth by microphysics, etc.) would not be able to cause the observed discontinuity of the global haze opacity distribution at any location. The stratospheric temperature asymmetry can be explained by a combination of asymmetric radiative heating rates and adiabatic heating due to vertical motion within the thermally indirect cell. A seasonal variation of haze particle number density is unlikely to be responsible for this asymmetry. It is likely that a thermally indirect cell covers the upper portion of the main haze layer. An artificial damping of the meridional circulation enables the formation of high-latitude jets in the upper stratosphere and weaker equatorial superrotation. The latitudinal temperature distribution in the stratosphere is better reproduced.