A new method based on accelerated solvent extraction (ASE) combined with response surface methodology (RSM) modeling and optimization has been developed for the extraction of four lignans in Fructus Schisandrae (the fruits of Schisandra chinensis Baill). The RSM method, based on a three level and three variable Box-Behnken design (BBD), was employed to obtain the optimal combination of extraction condition. In brief, the lignans schizandrin, schisandrol B, deoxyschizandrin and schisandrin B were optimally extracted with 87% ethanol as extraction solvent, extraction temperature of 160 ° C, static extraction time of 10 min, extraction pressure of 1,500 psi, flush volume of 60% and one extraction cycle. The 3D response surface plot and the contour plot derived from the mathematical models were applied to determine the optimal conditions. Under the above conditions, the experimental value of four lignans was 14.72 mg/g, which is in close agreement with the value predicted by the model.