The mathematical model and boundary condition describing airflow around a mutual insert multi-greenhouse were established based on CFD theory. Wind pressure distribution was simulated and analyzed using Fluent software. The simulation results indicated that the primary effect on the greenhouse of wind was upward force. The numerical simulation results of wind pressure distribution along the greenhouse span were consistent with the wind tunnel test result. The wind pressure coefficient of windward and leeward was quadric distribution along the greenhouse's height, and correlative coefficient was 0.975 and 0.990 respectively. The simulative wind pressure coefficient was transformed to wind load shape coefficient, the coefficient provides reference for optimizing the structure of greenhouse.