专利名称：：冷连轧机基于机理与工况相结合的工作辊磨损预报方法
技术领域：
：本发明涉及一种冷连轧生产工艺技术，特别涉及一种适合于四辊冷连轧机基于机理与工况相结合的工作辊磨损预报技术。
背景技术：
：近年来，随着制造工业的快速发展，用户对冷轧板带及其涂镀产品尺寸精度的要求越来越苛刻。而轧辊磨损的计算精度直接关系到负载辊型、辊缝形状以及辊系变形的计算精度，与弯辊力的正确设定、分段冷却补偿系统的正确利用密切相关，直接影响到成品的板形质量。轧辊磨损是一个缓慢的累积过程，影响轧辊磨损的因素很多而且复杂，工况不同，磨损也不同。以往现场对冷轧轧制过程轧辊磨损的机理模型研究较少，并且主要偏重于热轧(主要文献1、SpuzicS，StraffordKN，SubramanianC，etal.Wearofhotrollingmillrolls:anoverview[J].Wear,1994，176(2):261-271;2、陈连生，黄传清，连家创.热带钢连轧机轧辊磨损数值测量及分析.轧钢，2002，19(1):23-25;3、陈连生，连家创.热带钢连轧机精轧轧辊磨损计算理论.河北理工学院学报，2001，23(1):24-28)，对轧制过程中轧辊磨损的测量技术研究较多，得到较多的磨损规律曲线，但是至今没有形成比较准确的轧辊磨损计算模型。特别的，现有关于冷轧过程的轧辊磨模型往往没有将机理与实际轧制工况相结合，所得出的磨损模型计算方法通用性不强(主要文献有付华高，李长生，等.冷轧轧辊磨损数学模型研究.钢铁研究，2003，(4):26-29)。为此，本发明经过大量的现场试验与理论研究，充分考虑到冷连轧过程中的设备与工况特点，将磨损机理与工业生产实际磨损数据的测量结果有机结合，建立了一套适合于四辊冷轧机的工作辊磨损预报技术，不但适用于工作辊的对称磨损而且适用于不对称磨损，大大提高了现场工作辊磨损的预报精度，为辊型设计、板形控制奠定了坚8实的理论基础，有效的保证了产品板形质量的稳定，并能给现场创造较大的经济效益。本发明方法的原理清晰明了，计算速度快，适于现场使用。
发明内容为了克服现有技术存在的上述问题，本发明提供一种适合于四辊冷连轧机基于机理与工况相结合的工作辊磨损预报技术，将基本轧制磨损原理与现场实际工况下的轧辊实际测量数据紧密结合，建立一套既考虑机理又考虑工况、不但适用于工作辊的对称磨损而且适用于不对称磨损的工作辊预报方法，提高工作辊磨损的预报精度，为板形控制与辊型设计奠定基础，最大程度的保证产品板形质量，为现场创造经济效益。为了实现上述目的，本发明采用了以下技术方案一种适合于四辊冷连轧机基于机理与工况相结合的工作辊磨损预报方法，包括以下步骤(a)构造工作辊磨损模型考虑到机理与实际工况综合影响的工作辊磨损预报模型可以采用以下结构R(x,)=G(Xi).J(x,)式中R0O—工作辊磨损量；G(x,)—工作辊磨损量中的工况影响方程，用于修正机理模型；J(x,)—工作辊磨损量中的机理影响方程；x,—工作辊上横向第^点的坐标，原点在工作辊中部。(b)构造工作辊磨损工况影响方程根据现场实际经验以及数据回归整理结果，构造出工作辊磨损量中的工况影响方程如下6k=0式中A—工况影响系数，将根据工业生产实际磨损数据的测量结果在步骤(f)中求出。(C)构造工作辊磨损量中机理影响方程根据四辊冷连轧机轧制过程的生9产工艺特点，构造出工作辊磨损量中的机理影响方程如下J")=￡[Jlm0,)+J2m0,)+J3(x,)+J4o,)]式中"一工作辊换辊周;月内总的轧制巻数；L")一在轧辊换辊周期内第m巻带材轧制过程中工作辊与带材之间的滑动磨损量；L")一在轧辊换辊周期内第m巻带材轧制过程中工作辊与带材之间的滚动磨损量；、")一在轧辊换辊周期内第m巻带材轧制过程中工作辊与支撑辊之间的滑动磨损量；、")一在轧辊换辊周期内第m巻带材轧制过程中工作辊与支撑辊之间的滚动磨损量。(d)解析计算工作辊磨损中的机理影响部分根据步骤(c)可以知道，工作辊磨损中的机理影响部分的求解实际上包括工作辊与带材之间的滑动磨损量L(x,)、工作辊与带材之间的滚动磨损量、(x,)、工作辊与支撑辊之间的滑动磨损量L(x,)以及工作辊与支撑辊之间的滚动磨损量J^(x,)的计算，基本方案如下dl)轧辊换辊周期内第m巻带材轧制过程中工作辊与轧件之间的相对滑动引起的工作辊磨损量L(x,)的计算，由以下模型得出<formula>formulaseeoriginaldocumentpage10</formula>式中^一工作辊与轧件接触的滑动磨损系数；^(x,)—在轧辊换辊周期内第m巻带材轧制过程中单位宽度轧制压力横向分布；4一在轧辊换辊周期内第m巻带材的轧制长度；A—工作辊直径；在轧辊换辊周期内第m巻带材轧制过程中轧制前滑系数；儿,一在轧辊换辊周期内第m巻带材轧制过程中工作辊辊面与轧件间每次接触的相对滑动距离；C^一在轧辊换辊周期内第m巻带材轧制过程中横向流动影响系数；Bm—在轧辊换辊周期内第m巻带材轧制过程中轧件宽度。d2)轧辊换辊周期内第m巻带材轧制过程中工作辊与轧件之间相互滚动接触引起的工作辊磨损量L(x,)的计算，由以下模型得出式中/t2—工作辊与带材接触的滚动磨损系数；/、一在轧辊换辊周期内第m巻带材轧制过程中考虑压扁后的变形区接触弧长；d3)在轧辊换辊周期内第m巻带材轧制过程中工作辊与支撑辊之间相互滑动接触引起的工作辊磨损量k(x,)的计算，由以下模型得出工—;式中^一工作辊与轧件接触的滑动磨损系数；aOO—在轧辊换辊周期内第m巻带材轧制过程中工作辊与支撑辊之间的单位宽度轧制压力横向分布；、*200—在轧辊换辊周期内第m巻带材轧制过程中工作辊、支承辊的负载半径分布值；A、w2m—在轧辊换辊周期内第m巻带材轧制过程中工作辊与支承辊转动的圆周速度；11A、i2—工作辊与支承辊的公称半径；2^—在轧辊换辊周期内第m巻带材轧制过程中工作辊与支承辊间的接触压扁宽度；d4)在轧辊换辊周期内第m巻带材轧制过程中工作辊与支撑滚之间相互滚动接触引起的工作辊磨损量L(x,)的计算，由以下模型得出w4m(x,)=va(x,)式中、一工作辊与支承辊接触的滚动磨损系数d5)工作辊磨损量中机理部分总的影响计算，综合dl)、d2)、d3)、d4)得出工作辊磨损量中机理部分总的影响为J",)=t[J,m",)+",)+、",)+、")]乂'm+C,,,,,j05B,(1+/,,,,)^0.5Bm卜4.+va")^M〔(e)测量与收集实际工作辊磨损数据选择N(为了充分反映轧辊的在轧制过程中的设备与工艺工况，N—般要大于15)个典型的轧辊周期作为轧辊磨损模型建模的背景工况，实际测量出各个轧辊周期内的轧辊实际磨损数据R;(x')"'=1，2，'"，^)。在实际磨损数据的测量中必须保证以下四个条件(1)凉辊24小时以上，使轧辊内外温度、表面温度分布趋于一致，有效消除热变形对轧后轧辊磨损数据的影响以及消除温度变化对所测轧辊轧前直径的影响；(2)在进行轧后数据测量时必须在最大、最小两个方向上进行两次测量，磨损数据取两者的平均值，以减小数据偏差。(3)在测量轧辊磨损量前，必须认真清除测量部位辊面氧化皮和其他粘附物；(4)测量必须进行轧辊两端轴瓦调平，保证测量中心线与被测轧辊轴线一致。(f)求解工况影响系数^利用实际磨损的测量数据与工作辊磨损量中机理影响计算结果，求解工况影响系数^，基本步骤如下fl)构造目标函数考虑到工作辊磨损预报模型的预报值应该最大程度的与现场的实际测量值相吻合，从而构造出工况影响目标函数F(X)(F(Z)越小，表示相关轧辊磨损模型越精确)，如下式所示式中Z=(ao,a丄,a2,"3,a4,a5,a6);&(X,)—工作辊磨损预报模型的预报值。f2)给定工况影响系数的初始值X。^a。。f3)计算出N个典型轧辊周期中每个轧辊周期工作辊磨损量中机理影响因素JjOO，基本表达式如下~「Jj&)=HJ—(x')+J—")+J河")+J41乂m)+C"my0.5B—式中一第j个典型的轧辊周期轧辊所轧制的总的钢巻数(其它各个参数的下标与此意义相同)。f4)计算出当前状态下典型轧辊周期中工作辊磨损量的理论计算值，如下式所示Rj(x,一G(Xi).Jj(x,)、u.,乂ffy.+c呵'0.5B—f5)计算出目标函数F(I。)的数值。+26—wlmj.o,).《,.+A4f6)判断powell条件是否成立，如果成立，则取最优工况影响系数值X^X。结束计算，转入步骤f7);如果不成立则调整A，转入步骤f3)，继续13寻优。f7)输出最优工况影响系数值X,(。。,,^，^,^,^A"。(g)建立基于机理与工况相结合的工作辊磨损预报模型将最优工况影响系数A—,~，~,~,,,％}代入相关方程，建立起一套机理与工况相结合的工作辊预报模型y咖T、0.5B—、4、]+26.vs"")4K少加《(h)特定轧制周期内工作辊磨损的预报利用步骤(g)所建立的相关模型，将特定工作辊换辊周期内的各个钢巻带材轧制时的轧制参数如带材轧制过程中单位宽度轧制压力横向分布值、轧制长度、轧制前滑系数等参数代入模型就可以特定轧制周期内的工作辊磨损曲线，完成工作辊磨损的预报。下面结合附图对本发明的具体实施方式做进一步详细具体的说明。图1是四辊冷连轧机基于机理与工况相结合的工作辊磨损预报总框图2是工作辊磨损机理影响的解析计算框图3是工况影响系数的求解过程图。具体实施例方式为实际验证本专利所述技术方案中轧辊磨损方法的正确性和参数优化的准确性，以某1420五机架冷连轧机为例，取一个两个轧制单位对第5机架进行模型的仿真检验。结合轧制单位的生产条件，将轧制单位的实际工艺参数和表1所示的参数代入工作辊磨损模型进行磨损量的预报计算。模型预报的磨损量和对应的实测磨损量平均相对误差低于5.0%，最大误差低于15%，满足了工程上对轧辊磨损精度的预报要求。为了进一步说明本发明专利的实施过程，现以某1420五机架冷连轧机第514机架为例说明简单说明整个模型的建模与预报过程首先，收集1420五机架冷连轧机20个典型的实际轧制周期(排除异常换辊)，测出20个典型轧制周期的轧辊磨损数据《(x,)(在测量时，每根轧辊沿着辊身长度方向每隔10mm测量一个点。上下辊都测量，并且将平均值作为实际磨损量。这部分工作是由操作工在磨辊前在数字磨床上进行)，同时记录下每个轧制周期内所轧制的钢巻巻号、钢种、规格、轧制参数。根据现场统计，平均每个轧制周期轧制长度为150公里左右，大约轧制15-25个巻左右(轧制周期内的钢巻个数取决于钢巻大小及带材厚度等因素)换辊。因此，这二十个典型轧制周期共涉及到453巻带材。其次，取工况影响系数的初始值Z。={a。。,"1Q,a2。,a3。,a4。,^。，"6。}={1,1,1,U,l};再次，构造出工作辊磨损工况影响方程6G(X,)二Z"乂=1+X,+X,2+X,3+x;+x，+<随后，利用已经存在的轧制工艺模型及实际轧制工艺参数计算每个轧制周期内工作辊磨损量中机理影响因素JjOc,.)(需要说明的是，在实际计算过程中JjOO是一个比较复杂的模型，通过编制程序计算而出，并且该部分都可以通过现有模型来计算)，根据文献，在模型计算过程中所涉及的模型系数如表1所示。表l模型参数模型参数数值模型参数数值&2.123xlO-6249.32.174xl(T*-0.24362.093xl(r61201.2&2.145xl(T40.51231.463-0.1224-0.174—一然后，利用本发明所构造的公式G(;c,)J(x,.)计算出初始工况系数{1，1，1，1，1，1}时的各个轧制周期内的轧辊磨损量计算值及/;0;15随后，计算出当前工况系数下的目标函数=么反0")-《"')]2=241.9852随后，将计算结果代入powell优化条件，判断powell条件是否成立，如果成立，则取最优工况影响系数值^=1。结束计算，如果不成立则调整X。，继续寻优。最后，根据20个典型轧制周期的实际轧辊磨损数据得到最优工况影响系数值^={,,,％,,％,%}={0.234,0.126,0.267,0.387,0.056,—0.00786,—0.000432}。在得到上述工况系数之后，就可以利用该最优工况系数值配合机理模型进行轧辊磨损预报，为了说明模型精度，选择两个典型的轧制周期对磨损进行预报，相关典型轧制周期内所轧制的钢巻情况分别如表2、表3所示表2工1作辊1的轧制周期轧制序号钢种5#入口厚度(mm)5#出口厚度(mm)5#入口张力(MP)5#出口张力(MP)5#轧制压力(KN)5#带材出口速度(m/min)5#轧制公里数(km)1MRT隱4CA0.35830.243174.5571.236696.461123.090.002MRT-4CA0.35840.243175.6267.797068.381177.137.963MRT-4CA0.31620.214166.8768.607207.881126.2215.714MRT-4CA0.31410.212172.8171.337036.841148.9124.52MRT-4CA0.31410.212164.1069.547012.561094.3330.006MRT-4CA0.31410.212166.6170.676962.461183.8035.21MRT-4CA0.31450.212164.2267.136690.881116.0640.328MRT懇4CA0.29140.202161.9768.896398.001198.7145.639MRT-4CA0.29130.202158.2969.747122.561198.9854.8010MRT-4CA0.29150.202151.5866.976954.601233.9960.2511MRT-4CA0.29130.202156.1367.016797.201221.2065.9512MRT-4CA0.32750.222147.5468.806654.661257.3975.2013MRT-4CA0.32690.222150.7067.156672.461228.6080.5914MRT-4CA0.32750.222145.2368.536662.461274.2586.1015MRT-4CA0.32770.222142.5867.956712.341255.0691.4716MRT-4CA0.32730.222143.9168.706664.761222.0096.8117MRT-4CA0.32730.222139.3569.296823.881219.52102.0418MRT-4CA0.35750.243141.2566.746335.001256.11107.4919MRT-4CA0.35780.243135.8067.706549.321236.81112.4420MRT-4CA0.31600.214134.0667.016818.161262.42119.781621MRT-4CA0.31640.214123.5370.106913.201287.60131.3722MRT-4CA0.31700.214120.2068.106789.861271.66143.20表3工'作辊2的轧制周期乳制序号钢种5#入口厚度(mm)5#出口厚度(mm)5#入口张力(MP)5#出口张力(MP)5#轧制压力(KN)5#带材出口速度(m/min)5#轧制公里数(km)1MR3T40.32650.233173.5266.287753.041028.370.002MR3T40.32580.233171.6569.487463.801005.197.893MR3T40.32580.233174.8967.996951.721011.2415.774MR3T40.32620.233169.9167.677414.361009.2423.66MR3T40-32670.233171.3168.727363.561004.2531.546MR3T40.32600.233171.6167.647316.561001.7539.477MR3T40.32690.233175.2367.417531.661003.3147.368MR3T40.32620.233173.4268.917366.20画.2855.279SPHC0.33630.228177.5758.005855.821107.1872.3810SPHC0.33550.228176.5758.075995.441228.6581.3911SPHC0.33610.228175.3556.455937.341226.0190.4912SPHC0.33670.228175.2355.066192.481217.1399.5413SPHC0.33680.228177.9557.716189.601214.24108.5914SPHC0.33510.228174.9057.166182.861262.45117.5915SPHC0.33680.228174.8657.786388.801308.78126.56根据A以及对应的机理模型J(x,)计算出相应的1号工作辊与2号工作辊的磨损量。为了比较同时列出相关实际测量值以及误差值，如表4。表4利用本发明所述方法得出的轧辊磨损预报值与测量值对比1号工作辊2号工作辊横向点磨损实际测量值磨损模型预报值误差磨损实际测量值磨损模型预报值误差10.08970.0884.956250.04640.04411.8952120.08670.0828011.67960.04420.0449424.4966630.0770.079763.730490.04750.0457283.5849540.07720.0776031.778130.04730.0464590.52175250.0810.0759292.249190.04610.0471376.2602360.08580.0745621.408710.04710.0470641.09810.07250.0734068.386770.04460.0483411.2492180.07650.0724045.095710.04650.048875.35490.07150.0715218.705150.04540.0493520.029113100.06940.0707317.075270.04650.049791.91731110.07250.0700166.098640.04730.0501853.4265120.07050.0693636.171560.04760.050538L6簡130.06880.0687638.6550.04680.0508510.05395417140.06930.0682075.848410.04830.0511251.5771150.06740.067698.817530.04720.0513620.429708160.070.0672064.803380.04920.0515633.99202170.0670.06675110.5350.04680.051730.371791180.0670.0663222.9060.05040.0518651.01162190.06710.0659177.370620.04840.0519671.76348200.06810.06553210.25420.04720.052043.77091210.06740.0651665.219840.04950.0520843.31442220.06720.0648178.541970.0480.05213.54586230.06650.0644846.090140.04910.052093.03189240.0649rv八r<i，户厂pi">"/"in厂r\A厂no厂厂1.13313250.06410.0638585.684710.04920.0519970.376862260.06630.0635648.383770.04790.0519164.12628270.06410.0632818.396360.04780.0518141.27734280.06490.0630095.276870.04910.0516912.91454290.06550.0627455.418030.0柳0.0515494.20568300.06470.0624917.962180.04760.051393.41419310.06420.0622454.731650.04890.0512143.04502320.06060.0620076.963990.04770.0510222.32175330.06440.0617762.656930.04950.0508154.07423340.06410.0615523.254820.0490.0505953.97458350.06310.0613354.268890.04830.0503622.79732360.06420.0611244.193730.04810.0501174.79189370.06090.0609186.315010.0柳0.0498620.029747380.06340.0607183.758290.04780.0495974.23012390.06340.0605232.541080.04810.0493224.5374400.06120.0603332.809220.04770.049041.41601410.06340.0601489.306140.04460.0487515.12901420.06380.0599686.964030.04530.0484556.00709430.06230.0597917.72550.04470.0481534.02729440.06250.0596196.991130.04390.0478474.61028450.06330.059456.824440.04450.0475376.0821460.06370.0692854.85290.04260.0472236.93063470.06440.0591244.89180.04230.0469078.19271480.06190.0589667.845330.04320.0465894.73992490.06270.0588116.367580.0435■0.046276.202500.06410.0586611.45230.03980.045958.48701510.06350.0585113.002480.04430.045637.85621520.06290.0583667.371530.04220.0453117.20879530.06330.0582235.616480.04260.0449938.02084540.06250.0580838.174680.04130.0446767.06781550.06440.0579457.413740.04130.04436210.0233560.06520.0578112.08730.03930.0440511.334618570.06370.0576778.00柳0.04050.0437429.45508580.06180.05754711.66350.03890.0434376.8824590.06270.01741911.27050.03710.0431368.42349600.06330.05729250.04080.042849.49065610.06330.0571683.272960.04120.0425499.6865620,06130.0570474.350780.04050.0422626.93884630.0670.0569270.6741940.04170.0411812.0351640.06320.0568085.584770.03950.04170610.1133650.06210.0566922.566520.04040.0414378.70838660.06060.0565786.118720.03880,0411746.63764670.06190.056465U，flfU10u■1上680.06270.0563549.914120.0370.04066810.1217690.0620.0562440.4296080.04060.0404269.28355700.060.0561362.787720.03910.040196.43952710.06330.056031.572370.04060.0419628.4851720.0630.0559250.6089720.03950.03974111.2302730.06190.0558222.402820.04050.0395279.81978740.05990.055722.671720.04040.0393216.97911750.0560.0556191.615260.03850.0391220.680644760.06240.055520.5959280.03870.0387311.0264770.05990.0554221.965420.0380.0387477.4767780.05820.0553251.235020.03810.03857194041790.0570.0552291.534310.0390.0384023.10677800.06020.0551350.9326420.03860.038248.41395810.05910.0550422.106110.03730.0380866.86695820.05650.054951.459160.03850.0379382.74406830.05740.0548596.472620.03550.0377984.42736840.05840.0547692.42470.03860.0376646.21768850.05660.054682.002380.03680.0375373.39202860.05560.0545920.7533580.03770.037416L81224870.05140.0545065.80520.03960.0373016.0422880.05090.054420.0213330.03720.0371926.91547890.05230.0543354.164190.03870.0370893.89145900.05660.0542515.153850.0390.036994.14942910.05290.0541691.871440.03760.0368962.39802920.05170.0540872.108770.03760.0368074.61622930.05480.0540065.599230.03890.0367221.44981940.05420.0539254.082970.03820.036640.506847950.05120.0538460.7214190.03630.0365625.16782960.05030.0537674.735920.03830.0364866.89342970.05110.053692.379040.03730.0364135.06784980.05180.0536138.230340.03960.0363413.49951990.05190.0535374.300530.03790.036273.15337<table>tableseeoriginaldocumentpage20</column></row><table>权利要求1.一种冷连轧机基于机理与工况相结合的工作辊磨损预报方法，其特征是该方法包括以下步骤(a)在考虑到机理与实际工况综合影响的条件下构造工作辊磨损模型；(b)根据现场实际经验以及数据回归整理的结果，构造出工作辊磨损量中的工况影响方程；(c)构造工作辊磨损量中机理影响方程；(d)解析计算工作辊磨损中的机理影响部分；(e)测量与收集实际工作辊磨损数据选择N个典型的轧辊周期作为轧辊磨损模型建模的背景工况，实际测量出各个轧辊周期内的轧辊实际磨损数据(xi)(j＝1，2，…，N)；(f)利用实际磨损的测量数据与工作辊磨损量中机理影响计算结果，求解工况影响系数ak；(g)建立基于机理与工况相结合的工作辊磨损预报模型；(h)特定轧制周期内工作辊磨损的预报利用步骤(g)所建立的相关模型，将特定工作辊换辊周期内的各个钢卷带材轧制时的轧制参数如带材轧制过程中单位宽度轧制压力横向分布值、轧制长度、轧制前滑系数等参数代入模型，以确定轧制周期内的工作辊磨损曲线，完成工作辊磨损的预报。2.根据权利要求1所述的冷连轧机基于机理与工况相结合的工作辊磨损预报方法，其特征是步骤(a)中所述考虑到机理与实际工况综合影响构造工作辊磨损模型可以采用以下结构R(x,"G(XiW(x,)式中，R(x,)为工作辊磨损量；G(x,)为工作辊磨损量中的工况影响方程，用于修正机理模型；J(x,)为工作辊磨损量中的机理影响方程；x,为工作辊上横向第z'点的坐标，原点在工作辊中部。3.根据权利要求1所述的冷连轧机基于机理与工况相结合的工作辊磨损预报方法，其特征是步骤(b)中所述根据现场实际经验以及数据回归整理结果，构造出工作辊磨损量中的工况影响方程可以采用以下结构<formula>formulaseeoriginaldocumentpage3</formula>式中，^为工况影响系数。4.根据权利要求1所述的冷连轧机基于机理与工况相结合的工作辊磨损预报方法，其特征是步骤(C)中所述工作辊磨损量中机理影响方程可以采用以下结构构造<formula>formulaseeoriginaldocumentpage3</formula>式中，n为工作辊换辊周期内总的轧制巻数；L(x,)为在轧辊换辊周期内第m巻带材轧制过程中工作辊与带材之间的滑动磨损量；L(x,)为在轧辊换辊周期内第m巻带材轧制过程中工作辊与带材之间的滚动磨损量；L")为在轧辊换辊周期内第m巻带材轧制过程中工作辊与支撑辊之间的滑动磨损量；L(;c,)为在轧辊换辊周期内第m巻带材轧制过程中工作辊与支撑辊之间的滚动磨损量。5.根据权利要求1所述的冷连轧机基于机理与工况相结合的工作辊磨损预报方法，其特征是步骤(d)中所述解析计算工作辊磨损中的机理影响部分实际上包括工作辊与带材之间的滑动磨损量L(x,)、工作辊与带材之间的滚动磨损量L(x,.)、工作辊与支撑辊之间的滑动磨损量L(x,)以及工作辊与支撑辊之间的滚动磨损量J^(x,)的计算，其计算步骤如下dl)轧辊换辊周期内第m巻带材轧制过程中工作辊与轧件之间的相对滑动引起的工作辊磨损量LOc,)的计算，由以下模型得出-<formula>formulaseeoriginaldocumentpage4</formula>式中A—工作辊与轧件接触的滑动磨损系数，A(^—在轧辊换辊周期内第m巻带材轧制过程中单位宽度轧制压力横向分布，;一在轧辊换辊周期内第m巻带材的轧制长度，A—工作辊直径，/;一在轧辊换辊周期内第m巻带材轧制过程中轧制前滑系数，L一在轧辊换辊周期内第m巻带材轧制过程中工作辊辊面与轧件间每次接触的相对滑动距离，C,一在轧辊换辊周期内第m巻带材轧制过程中横向流动影响系数，Bm—在轧辊换辊周期内第m巻带材轧制过程中轧件宽度；d2)轧辊换辊周期内第m巻带材轧制过程中工作辊与轧件之间相互滚动接触引起的工作辊磨损量L(x,)的计算，由以下模型得出<formula>formulaseeoriginaldocumentpage4</formula>式中^—工作辊与带材接触的滚动磨损系数，/、一在轧辊换辊周期内第m巻带材轧制过程中考虑压扁后的变形区接触弧长；d3)在轧辊换辊周期内第m巻带材轧制过程中工作辊与支撑辊之间相互滑动接触引起的工作辊磨损量k(x,)的计算，由以下模型得出<formula>formulaseeoriginaldocumentpage0</formula>式中工作辊与轧件接触的滑动磨损系数，aOO—在轧辊换辊周期内第m巻带材轧制过程中工作辊与支撑辊之间的单位宽度轧制压力横向分布，t")、足00—在轧辊换辊周期内第m巻带材轧制过程中工作辊、支承辊的负载半径分布值，w^、^^一在轧辊换辊周期内第m巻带材轧制过程中工作辊与支承辊转动的圆周速度，《、i2—工作辊与支承辊的公称半径，在轧辊换辊周期内第m巻带材轧制过程中工作辊与支承辊间的接触压扁宽度；d4)在轧辊换辊周期内第m巻带材轧制过程中工作辊与支承辊之间相互滚动接触引起的工作辊磨损量J4m(x,)的计算，由以下模型得出W4m(x,)=A:4.a")-式中工作辊与支承辊接触的滚动磨损系数；d5)工作辊磨损量中机理部分总的影响计算，综合dl)、d2)、d3)、d4)得出工作辊磨损量中机理部分总的影响为4,,""+C""''lo.5B卜/,,,,)."y....."(0.5B6.根据权利要求1所述的冷连轧机基于机理与工况相结合的工作辊磨损预报方法，其特征是步骤(f)利用实际磨损的测量数据与工作辊磨损量中机理影响计算结果，求解工况影响系数^的基本步骤如下fl)构造目标函数考虑到工作辊磨损预报模型的预报值应该最大程度的与现场的实际测量值相吻合，从而构造出工况影响目标函数F(X)(F(X)越小，表示相关轧辊磨损模型越精确)，如下式所示<formula>formulaseeoriginaldocumentpage6</formula>&(x,)—工作辊磨损预报模型的预报值；G)给定工况影响系数的初始值X。={a。。，fl1Q，a2。，a3。，"4。,,a6。};f3)计算出N个典型轧辊周期中每个轧辊周期工作辊磨损量中机理影响因素J」(x,)，基本表达式如下<formula>formulaseeoriginaldocumentpage6</formula>式中一第j个典型的轧辊周期轧辊所轧制的总的钢巻数(其它各个参数的下标与此意义相同)；f4)计算出当前状态下典型轧辊周期中工作辊磨损量的理论计算值，如下式所示<formula>formulaseeoriginaldocumentpage6</formula>f5)计算出目标函数F(X。)的数值;+26—伤)判断powell条件是否成立，如果成立，则取最优工况影响系数值A=1。结束计算，转入步骤f7);如果不成立则调整X。，转入步骤f3)，继续寻优；f7)输出最优工况影响系数值^={7.根据权利要求1所述的冷连轧机基于机理与工况相结合的工作辊磨损预报方法，其特征是步骤(7)中所述机理与工况相结合的工作辊磨损预报模型如下<formula>formulaseeoriginaldocumentpage7</formula>全文摘要本发明公开一种适合于四辊冷连轧机基于机理与工况相结合的工作辊磨损预报技术，其特征是a.构造工作辊磨损模型；b.构造工作辊磨损工况影响方程；c.构造工作辊磨损量中机理影响方程；d.解析计算工作辊磨损中的机理影响部分；e.测量与收集实际工作辊磨损数据；f.求解工况影响系数α<sub>k</sub>；g.建立基于机理与工况相结合的工作辊磨损预报模型；h.特定轧制周期内工作辊磨损的预报。采用本发明所提供的工作辊磨损预报技术，充分考虑到冷连轧过程中的设备与工况特点，将磨损机理与工业生产实际磨损数据的测量结果有机结合，适用于工作辊的对称磨损而且适用于不对称磨损，大大提高了现场工作辊磨损的预报精度，为辊型设计、板形控制奠定基础。文档编号B21B38/00GK101507978SQ20091007399公开日2009年8月19日申请日期2009年3月20日优先权日2009年3月20日发明者周莲莲,白振华申请人:燕山大学