1.一种DCR机组大变形下以稳轧为目标的轧制参数优化方法，其特征在于：其包括以下由计算机执行的步骤：(a)收集冷轧机组的主要设备与工艺参数，主要包括以下步骤：a1)收集冷轧机组的轧辊工艺参数，主要包括：工作辊半径R、表面粗糙度Rar、工作辊的弹性模量E、工作辊的泊松比ν；

a2)收集冷轧机组相关轧制工艺参数，主要包括：带材的平均变形抗力Km和屈服强度σs、带材的宽度B、来料的厚度h0、压下率ε、正常轧制速度V、轧制压力设定值P1、前张力设定值T1、后张力设定值T0；

a3)收集工艺润滑制度参数，主要包括：乳化液浓度c、初始温度t0、流量w、以及乳化液的动力粘度η0、压缩系数θ；

a4)收集冷轧机组的工艺特征参数，主要包括：许可最大轧制压力P*，临界最大前滑值S*，许用最小、最大前张力σ1min、σ1max，许用最小、最大后张力σ0min、σ0max；

(b)定义单位前张力σ1、单位后张力σ0，最佳单位前张力设定值σ1y，最佳单位后张力设定值σ0y，定义控制目标函数为G1(X)，给定单位前张力设定步长Δσ1、单位后张力设定步长Δσ0，设定目标函数初始值G0＝0；

(c)初始化单位前张力中间过程参数k1＝0；

(d)计算单位前张力σ1＝σ1min+k1Δσ1(e)初始化单位后张力中间过程参数k0＝0；

(f)计算单位后张力σ0＝σ0min+k0Δσ0(g)计算当前工况下摩擦系数 式中：a为液体摩擦影响系数；b为干摩擦影响系数；Bξ为摩擦系数衰减指数；ξ01为光辊轧制时的动态油膜厚度ψ为润滑油膜速度影响系数 ξ02为轧辊粗糙度对润滑油膜厚度影响量，主要取决于轧辊实际粗糙度；

(h)计算当前工况下轧制压力P、前滑值S，轧制压力 式中：为强 度 张 力规 格 系数 为 规 格强 度 系 数为规格压下系数 前滑值计算模型

Δh为压下量Δh＝h0ε，R'为工作辊压扁半径(i)计算控制目标函数G1(X)，计算模型为

(j)判断不等式 是否成立，如果成立，则令G0＝G1，令最佳后张力设定值σ0y＝σ0、最佳前张力设定值σ1y＝σ1，转入步骤(k)；否则，直接转入步骤(k)；

(k)判断不等式σ0＜σ0max是否成立，如果不等式成立，则令k0＝k0+1，转入步骤(f)，否则转入步骤(l)；

(l)判断不等式σ1＜σ1max是否成立，如果不等式成立，则令k1＝k1+1，转入步骤(d)，否则转入步骤(m)；

(m)输出最佳单位前张力σ1y、单位后张力σ0y，完成大变形轧制时二次冷轧机组以稳定轧制为目标的轧制工艺参数优化。