Abstract: The Monte Carlo （MC）-discrete ordinates （S_N） bidirectional coupling method is an efficient approach to solve shielding calculation of the large complex nuclear facility. The test calculation was taken by the application of the MC-S_N bidirectional coupling method on the shielding calculation of the large PWR nuclear facility. Based on the characteristics of NUREG/CR-6115 PWR benchmark model issued by the NRC, 3D Monte Carlo code was employed to accurately simulate the structure from the core to the thermal shield and the dedicated model of the calculation parts locating in the pressure vessel， while the TORT was used for the calculation from the thermal shield to the second down-comer region. The transform between particle probability distribution of MC and angular flux density of S_N was realized by the interface program to achieve the coupling calculation. The calculation results were compared with MCNP and DORT solutions of benchmark report and satisfactory agreements were obtained. The preliminary validity of feasibility by using the method to solve shielding problem of a large complex nuclear device was proved.