A method for estimating throughput of a production line that comprises multiple work stations including industrial robots and/or dedicated machines, characterized in that the method comprises: 101) obtaining a process time parameter of each of the multiple work stations, a travel time parameter between two neighboring work stations, as well as the number of buffers between the two neighboring work stations by simulations on a single work station; 102) selecting a bottleneck work station in the production line according to a preset criterion that is defined by the process time parameter and travel time parameter of the work station; 103) calculating bottleneck sub-throughput ρ BN of a bottleneck subline, upstream sub-throughput ρ up of an upstream subline and downstream sub-throughput ρ down of a downstream subline according to the process time parameter, the travel time parameter and the number of buffers, wherein the bottleneck subline includes the bottleneck work station, at least one upstream work station of the bottleneck work station, and at least one downstream work station of the bottleneck work station, the upstream subline comprises a first work station of the bottleneck subline and at least one upstream work station of the first work station, the downstream subline comprises a last work station of the bottleneck subline and at least one downstream work station of the last work station; 104) calculating a corrected process time T' i -1,pro of the first work station of the bottleneck subline according to the process time parameter T i -1, pro of the first work station, the bottleneck sub-throughput ρ BN and the upstream sub-throughput ρ up using the equation below: T ′ i − 1, pro = T i − 1, pro + 3600 ρ up − 3600 ρ BN; 105) calculating a corrected process time T' i +1, pro of the last work station of the bottleneck subline according to the process time parameter T i +1,pro of the last work station, the bottleneck sub-throughput ρ BN and the downstream sub-throughput ρ down using the equation below: T ′ i + 1, pro = T i + 1, pro + 3600 ρ down − 3600 ρ BN; and 106) re-calculating the bottleneck sub-throughput as the throughput of the production line according to the corrected process time of the first work station and the corrected process time of the last work station to accurately estimate the throughput of the production line where only a small amount of simple data input is provided, and use the re-calculated sub-throughput as the estimated throughput of the entire production line. 101) obtaining a process time parameter of each of the multiple work stations, a travel time parameter between two neighboring work stations, as well as the number of buffers between the two neighboring work stations by simulations on a single work station; 102) selecting a bottleneck work station in the production line according to a preset criterion that is defined by the process time parameter and travel time parameter of the work station; 103) calculating bottleneck sub-throughput ρ BN of a bottleneck subline, upstream sub-throughput ρ up of an upstream subline and downstream sub-throughput ρ down of a downstream subline according to the process time parameter, the travel time parameter and the number of buffers, wherein the bottleneck subline includes the bottleneck work station, at least one upstream work station of the bottleneck work station, and at least one downstream work station of the bottleneck work station, the upstream subline comprises a first work station of the bottleneck subline and at least one upstream work station of the first work station, the downstream subline comprises a last work station of the bottleneck subline and at least one downstream work station of the last work station; 104) calculating a corrected process time T' i -1,pro of the first work station of the bottleneck subline according to the process time parameter T i -1, pro of the first work station, the bottleneck sub-throughput ρ BN and the upstream sub-throughput ρ up using the equation below: T ′ i − 1, pro = T i − 1, pro + 3600 ρ up − 3600 ρ BN; 105) calculating a corrected process time T' i +1, pro of the last work station of the bottleneck subline according to the process time parameter T i +1,pro of the last work station, the bottleneck sub-throughput ρ BN and the downstream sub-throughput ρ down using the equation below: T ′ i + 1, pro = T i + 1, pro + 3600 ρ down − 3600 ρ BN; and 106) re-calculating the bottleneck sub-throughput as the throughput of the production line according to the corrected process time of the first work station and the corrected process time of the last work station to accurately estimate the throughput of the production line where only a small amount of simple data input is provided, and use the re-calculated sub-throughput as the estimated throughput of the entire production line.