D0I:10.13374/j.issm1001-053x.1991.01.021 第13卷第1期 北京科技‘大学学报 Vol,13 No,1 1991年1月 Journal of University of Science and Technology Beijing J立.1991 用流函数分析平面轧制问题 汪家才· 痛要：给出了一种用浅鱼数求解理题刚塑性材科的平面应变轧制的方法。速度场 被分解为基暗速度场和附加速度场。基陆速度场满足边界上定的速度条并，附加速 度场满足齐次边界条件。与这青部分速度场相财应，有基请流函敦和附如施函数。基 础流面数可以确定地写出，附加流函数罪露助于Weierstras3定理写成完备空间的 向壹栈，即多项式。通过使全勇率帮小化，可以将多项式系数确定。用这一方法求得 了速度场、应力场、塑性区前后边界，接敏真上中性点的位置和轧制单位压力，并与 工程方法的计算结果作了比校。 尖输词：注函数，上限原理，全功率，平面轧制 Analysis of Plate Rolling Process Using Flow Function Wa排gJiacai" ◆ ABSTRACT:The velocity field is devided into the basic and additional fields. The basic velocity field should satisfy given conditions on boundary,but the additional velocity field should only satisfy the uniform boundary conditions. Corresponding to those two parts of velocity field there are basic flow function and additional fiow function.The basic flow function can be written definitly, while the additional flow function,according to the'Weierstrass theorem,has been written as a family of vectors in complete space,i.e.a polynomial.With the aid of minimization of the total power the coefficients of the polynomial have been defined.By using this method,the following results are obtained: velocity field,field of strain-rate,stress field,the back and front boundaries 1990一09一11收稿 ,机械工程系(Department of Mechanical Engincering)
户 第玲春第 1 期 i , , i年 落月 北 京 科 技 大 学 学 报 J o u r o a l o f U o i 二 e r s主t了 o f S e i o n o e 牌取 d T e c五n p! p g y B e i j i n暇 y o l . i 争 N o 一 1 J t 城 。 ` i , s一 ’ 用旅函数分析平面轧制问题 汪 家 才 ’ 勺 热熬锻然群熬 场 速 甚 的 向盆焦 即 多攻式 。 笼过使众功* 橄小化 , 率以挤薪项式拜数康定 . 用这一方法 求 得 了邀度场 、 应力场 . 塑性区贫后边界 , 接城弧 上中佳点e 位里和轧翻单位 压力 , 井 与 一 一 一 一 工穆方法的计算结来作了比技 。 哭恤润 : 卜 沈. 毅 , 上限服 , 全功串 , 平 面轧翻 人n a l y s i s p f P l a t e R o l l i n g p r o c e ` s u “ 恤 g F oI 脚 F u n ct i o n 、 甲 a 林夕 iJ a C a i . 几B s T RA c T : T h e v e l o c i t 了 f i e ld 1 5 ’ d e ` i d e d i n t o t 五e b a s i c a n d a d d i t i o 皿 a l f i e l通, · T从 b a s i e v e l o e i t 丫 f ie l d s b b u ld s a ti s f v 卫 i v e 红 e o n d i t i o n s o n b o u n d a r 了 。 b o t t h e a o d 主t l o n a l v e l o 仁 l t y i 姆 l d s . o u l往 0 巨玉y s a t 盈s l y t n e u n i 宜o r m O o u n o a r y c o n d 豆t l o n s - C o r r e s o o 妞 d i n 眨 t o t h o s e t w o o a r t s o f v e l o e i t v f i e ld t五e r e a r e b a s i e f l o w f u n e t i o n a n d a ` d 豆t i o n a l 1 10 , i u n c t i o n . l n e o a s i c t l o w t u n c t i o巨 c a皿 b e W r i t t e n 吐 e t i n i t l y - w 五i l e t h e a d d i ti o n a l f l o w f o n e t i o n , a e e o r d i n g t o t五e ’ W e i e r s t r a s ; t h e o er m , h a s 恤 e n w r i t t e n a s a f a m i l y o f 丫 , e t o r s i ; e o o P叉e t , s p a e e , i . e · a p o l; n o 二i a l · Wi t h t h e a l d o f m i n i m i z a t i on ’ o f t h e t o t a l p o w e r t h e e o e f f i e i e n t s o f t h e p o l y n o m i a l h ` v , b , e , d e f i“ e d : L By ” 5 1“ 9 t h i s 贝 e t五匕 d , t h e f o 叉1 0 , i n g r e s u l t s a r e o b t a i” e d , v e幕o e 滚t y 称。 l d , if e 互以 。 f , t r o i 。 一 r a t e , s t r e o s f i o ld , t h e b a o k a n d f r o o t b o o n d a r i e s 19 ” 一 9一1 侧有 机械工 程系 ( D c P ` r t 皿 e o t o f M e e h a o i e a l B n 名i 妞 e e r i n g ) DOI: 10. 13374 /j . issn1001 -053x. 1991. 01. 021
of plastic zone,the position of neuture point on the contact arc,and the rolling pressure acting on the contact surface of the strip. KEY WORDS:flow function,upper-bound theorem,total power,plate rolling 采用花函数按上限法求解金属成型问题，最初是应用于通过楔形模的平面应变挤压过 程，以及通过锥孔的轴对称挤压过程，后来又解决了圆柱辙粗这类非定常过程。对于平板轧 制二维花动问题的求解，4迄今未见报道。看来求解轧制问籍存在以下困准： (1)轧辊表面是柱面而不是平面，这给设定速度场带来困难。必须解决如何用一个函数 式统一表达整个轧制变形区速度场问题。 (2)待求的主要未知量是轧制接触压力P。“它沿轧辊法向作用，故在轧制过程中p不 做功，即欲求的轧制压力P不进入能量方程。 (3)在轧制接触弧上，以中性点为界，前滑区和后滑区相对于轧辊是“分泷”的。如果 不能确定接触弧上中性点的位置，就求不出速度场。 本文试国解决以上的困难，用泷函数求解二维轧制问题。 1基本假设 (1)变形视为二维流动 考虑宽展的三维轧制问题原则上也可以求解，得因三维的流函 y 数必须用两个函数式表达，问愿要复杂些，本 文不予讨论。 (2)变形金属视为理想刚塑性体不过读 者容易看出，文中提供的方法也可应用于弹塑 性强化介质，所以此假设并不降低方法的价 值。 (3)轧制接触真看作二次抛物线这是工 困1变形区儿何参数 Fig.1 The plastic sone and its 程中早已采用的一个合理假设。 parameters 变形区几何形状和参数如图1所示。 2轧制问题的上限原理 对于可变形连续体，设区域”的边界为∑，有下列普遍的能量方程： ∫xd=o5ar 此式表示沿边界Σ的表面力X,的功率，应等于区域V内的变形功率。上式左端的积分中， 主动力的功串视为正的，接触摩擦力·（星号表示“给定”）作为约束理想化的结果应看 作主动力。间断线上消耗的功率是负的。令工=S:+S2+y。+Y1,对于图1所示的无张力 轧制情况，有 2
o f P l a s ti e z o n e P r e s s肚 r e a c ti n g K EY W O RD忿 : t h e p o s i t i o n o f n e u t u r o p o i n t o牡 t五e e o n t a e t a r e , a n d t 五e r o l li n g o n t h e e o n t a e t s u r f a e e o f t h e s t r i P . f r o , 介 n o t一o n , 。 p p e r 一 b o u n d t h e o r e m , t吐。一 p o 谕e r , p l a t e r o l一i n g , 程 、 口 采用优函数按上限法录解金属成型间题 , 最初是应用于通过模形澳的平面应 变 挤 压 过 以及通过锥孔的轴对称拚压过程 , 后来又解决了 圆柱傲粗这类非定常过程 。 对于平板轧 制二维流动同题的求麟 , `辉今养助瞬 , 一 看乘求解浦荫 群在以姻难 : _ ( 1) 轧辊表面是柱面而不是乎面 , 这给瞬定速度场带来困难 , 必须解决如何用 一个函数 式统一表达整个乳制变形区速度场向题 。 ( 2) 待求的主要未知量是轧制接触压力 九 傲功 , 即欲求的轧麟压力 p 不进人能量方程 。 它沿轧辊法向作用 , 故在轧制 过程 中 户不 (3 ) 在牡制舞触弧上 , 以 中性点为界气 前借区和后指 区 一相对于轧辊聂 “ 分拢 ” 的 。 「 如果 不能确定接触弧大中性在的位置 , 就本不琳速度场 。 本文试图解决以上的困难 , 用优函数求解匕维轧制问题 。 舞 「 「 1 本 棍 译 (l) 变形视为二维沈动 考虑宽雇的三维轧翻何鹿翩上也可以录璐 , (碑脚主维的魂函 r r 下) 「 { ` _ 八 } 0zJ l l l . _ 1 数必须用两个函数式表达衬向肠要复杂些 , 本 文不予讨论 。 (2 ) 变形金属视为理想刚攀性体 不过读 图 1 变形区 几何参 数 F i g 。 I T 五e Pl a s t i e 名 o n e a n d i t s P a r a m 七t e r s 者容易看断 义中提供的方法妈可应用午弥塑 性哪化介厥 , 所以此假设井不降低 方 法 的 价 值 。 (3 ) 轧制接触弧看作二次抛杨线 这是工 程 中早已采用的一个合理假设 。 变形区几何形状和参数如 图 1 所示 。 2 轧制 一 问题 的上限原理 对于可变形连续体 , 设区域 犷 的 边界为 军 , 有下列普迫的能量方程 : I : X ` t, ; ` : 一 J v a 。 ! ` ! : ` 犷 此式表示沿边界 刃 的表面力 X ` 的功率 , 应等于区域 犷 内的 变形功率 。 上式左端的积分 中 , 主动力的功率视为正的 , 接触摩擦力 , . (星号表示 “ 给定 ” ) 作 为约束理想化的结果 应 看 作主动力 。 间断线上消耗的功 率是负的 。 令 E = S : + 5 2 + 夕。 + v : , 对于图 1 所示的无张 力 轧制情况 , 有
udro.ds-d (1) 或者写成： 它=它(W)+E(y0)+E(y1)+E(S1)-E(S,)=0 (2) 式中E称为全功率。△v是y上的间断速度，V.沿S的切向速度。由于在S和S2上T·的方向 不同，故相应的两项摩擦功率是异号的。此外还可看出，无张力的轧制问题，属于全功率的 值已知为零的一类问题，而轧制过程是靠后滑区上的摩擦力（工作剪应力）来实现的。 在式(1)中，各项积分都和速度场”，有关。上限定理证明，将任一运动可能的速度场代 入，均可使E≥0。如果某个”，使E=0,则此速度场”是真实的。如果一个运动可能的速 度场”，使E接近于零，则此”，是一个好的近似解。 3速度场的设定 3,1接触氯方程和无量纲坐标 接触抛物线方程是（见图1）： Y=h1+〔(h。-hi)112x2 (3) 采用无量纲坐标，x〔0,1)，ye〔0,1)，得 7=(1-)+ex2 (4 3,2用流西数表示速度场 根据流体力学，对任何无散又无源的二维流场，存在流函数中=中(×，y)，使 v,=3种13y,",=-3中13x (5) 3,3基础流西数°和基础速度场v: 将运动可能的速度分解为两部分〔1)：基础场和附加场。基础速度场应采取简单的形式， 并满足给定的边界速度条件，且不含未知量。基础速度场：可从基础流函数°导出。取 (见图1)〔2)： °=C,/Y 可以看出，|，=。=0，这是一条流线，即中心线。的°|，=r=C=-”h0,这是另一条流 线，即接触弧，C恰是变形区的上半流量。当采用相对坐标时，基础流函数的形式是： 中°=-vh,y/[1-e)+ex2) (6) 由基础流函数中°可按式()导出基础速度场)；。由9可按下式导出基础应变速串场 3
户 = 肛 _ , ` * 。 ` : d。 + 甲 、 !△。 ! d , , 甲 , · 。 ; d: 一 那 : , 。 : d , = J v 沙 , o , , 1 J s i J 吕 2 ( 1 ) 或者写 成 : E = E (犷 ) + E ( 夕。 ) + E ( ? i ) + E (5 1 ) 一 E ( S : ) = o ( 2 ) 式 中E 称为全功率 。 △。 是夕上的间断 速度 , t, : 沿 S 的切向速度 。 由于在 S : 和 S : 上 , * 的方向 不同 , 故相应的两项摩擦功率是异号的 。 此外还可看 出 , 无张力的轧制向题 , 属于全功率的 值已知为零的一类向题 , 而轧制过程是靠后滑区上的摩擦力 ( 工作剪应力) 来实现的 。 在式 l( ) 中 , 各项积分都和速度场 t, ` 有关 。 上限定理 证明 , 将任一运动可能的速度场代 ^, 均可使于) o0 如果某个 “ ` 使 E = ” , 则此 速度场 “ `是真东的 。 如果一个运动可 能 的 速 度场 。 . 使 E 接近于 零 , 则此 “ 。 是一个好的近似解 。 3 速度场的设定 3 . 1 接触弧方程和无皿纲坐标 接触弧 抛物线方程是 (见图 1 ) : Y = 入, + 〔 ( h 。 一 h : ) j 1 2〕 x Z ( 8 ) 采用无量纲坐标 , 二 ` 〔0 , 1〕 , 夕 ` 〔o , 1〕 , 得 Y = ( 1 一 ￡ ) + e x ( 4 ) 3 . 2 用滋函效表示璐度场 根据流体力学 , 对任何无散又无源的二 维流场 , 存在流函数 叻二 沪x( , y ) , 使 。 , = a 叻/ a 夕 , ” , 二 一 a护l a x ( 5 ) 3 。 3 甚摇 流函教沪 。 和基曲邃度场 川 将运动可能 的速度分 解为两部分 〔 又’ : 基础场和附加场 。 基础速度场应采取简单的形式 , 并满足给定的 边界速度条件 , 且不含未知量 。 基础速度场 ?z, 可从基础 流 函 数 犷 导 出 。 取 (见 图 1 ) 仁 幻: 护 o = C , I y 可以看出 , 犷 I , 二 。 = 。 , 这是一条流线 , 即中心线 。 护oI , 二 ; = C = 一 ” 。 h 。 , 这 是 另 一 条 流 线 , 即接触弧 , C 恰是 变形区的上半流量 。 当采用 相对坐 标时 , 基础流函数的形式是 : 护 O 二 一 。 。 h 。 夕 /〔l 一 。 ) + 。 二 “ 〕 ( 6 ) 由基础 流函数 功 。 可按式 (5 ) 导出基础速度场 ” 了 。 由 。 气可按下式 导出基础 应变 速率场
51 5i=(1/2)(av,/ax.+av./ax,) (7) 场5：和v?的表达式限于篇幅不予列出。 3.4附加流西救中“和附加速度场v; 附加速度场由附加流函数导出，它用来修正和补充基础速度场，且不应破坏已由基础场 满足了的边界速度条件。设定“时有以下考虑： (1)由对称条件✉(x,y)=v(,一y)和反对称条件”（￥，y)=-v,(x,-y),对照式 (5),可知·中y应为奇次。 (2)在中心层上应有5：，l-0=0,这就要求 (d中/dx2),0=(ò2中1∂y2)y=0 (3)在边界y=0和y=±了上，应有"=0，故的中必有因子y(少2-卫)。 综合以上各点，可以取： =-h7-7222。, (8) 此式的正确性可以用Weierstrass定理加以说明。对于二元函数，该定理说：在x∈〔0，l)和 yE〔0,1)上有界的任意连续函数f(x,y),可以用下列多项式来无限通近： f》=m22y m,n=0,1,2,… 考虑到前述特殊要求后，就得到式(8)，对式（8)加以整理，并取有限厦，最后得： -2会[y门 (9) 3.5完全流函敷中和剪应变速率强度H =+-[多+会多（医y全p,] (10) H=2√(5.+5.)2+（⑤9，+5，）2 (11) 4边界Y,Y1的确定 边界y。和y1可由沿边界法向的速度必须连续这一条件确定。先考察y。,由图2可见 Vo=-vi+Oj,V=v,i+v,i,n=sinai-cosaj 法向速度连续条件要求V。·n=V·n,即一v。sina=v,sina一v,cosa
考 ` . : 首 : . = (1 1 2) ( a ,t ` l a x , + a t, , / a 二 。 ) ( 7 ) 场古 ` . 和 。 弓的表达式限于篇幅不予列出 。 3 。 4 附加滋西教护 “ 和用加落度场 t, , 附加速度场由附加流 函数导出 , 它用来修正和补充基础速度场 , 满足了的边界速度条件 。 设定尹 时有以下 考虑 : (1 ) 由对称条件 。 二 ( 二 , 夕 ) 二 。 , (礼 一 y )和反对称条件 。 , ( 二 , y ) - s() , 苛知护中y 应为奇次 。 (2) 在中心层上应有 C , , , 。 。 二 0 , 这就要求 且不应破坏已 由基础场 一 ” , ( ` , 一 y) , 对 照式 ( J 乞 护 . / d 二 2 ) , . 。 二 ( J 艺 护 / J夕 “ ) , 一 。 (3 ) 在边界y = 。 和 y 二 士 Y 上 , 应有护 . = 。 , 故护 ` 中必有因子 夕 ( y Z 一 Y Z ) 。 综合以上各点 , 可以取 : · 「 ,一 、 ho掀刃 一 尹 )名 习 。 . 二 了歹 ” ’ ( 8 ) . 一 O 二 O 此式的正确性可以用 Wie 份 st ar , s 定趣加以说明 。 对于二元函数 , 该定理说 : 在二 〔 〔 0 , 1〕和 y 〔 〔。 , 1〕上有界的任意连续函数 f ( 二 , 川 , 可以用 ` 下列多项式来无限逼近 : f x( , y) = 替黑习习 “ 二“ ’ y ’ 沉 , n 二 O , 1 , 2 , ~ ’ “ 考虑到前述特殊要求后 , 就得到式 ( 8 ) , 对式 ( 8 ) 加以整理 , 并取有限项 , 最后得 : ` 六 二 广 一 , . 二 一 ” 。 ho 或 石 ` 二 L ` 一 1 2 . + 1 二1 : 一 〕 ( , , 3 。 5 完全流西救 砂和 剪应交邃率强度 H ` = *+o , 一认 浮 · 象郭 · (犷 一 ` , ” ` 乞 二 ’ “ 一 ’ ) ] (` ” , 一y 一 l 卜Y H = 2了 (犷乳 + 歹才 二 ) 2 + (多全 , + 歹少 , ) “ ( 1 1 ) 4 边界 0Y , 丫 , 的确定 边界 y 。和夕f 可 由沿 边界法向 的速度必须连续这一条件确定 。 先考察 夕 。 , 由图 2 可见 犷。 “ 一 秒。 i + Oj , F = t, x i + t, y i , ” = s i n a ` 一 e o s a j ( 法向速度连续条件要 求 v 。 · n 二 V · n , 即 一 ” 0 5 1扭 a = ” ; 5 1” a 一 ” , c o s q 弓
也就是u,/(v,+vo=tana。另一方面，沿y。有tanc=dy/dx,再考虑到式(5)，就得到： (x)dx +(/y)dy=-vody 此式左端为中的全微分，可见有中=-vy+f(x),利用条件l,0=0,故知f(x)=0,于 是：的=一uy=-v。hy(在Yo上) (12) 将式(10)代人式(12)，得人口边界的方程为 +22(,)-10 (13) 同理，可求得出口边界y1的方程为： +2总(y”-,)-。-0 m-12-2m-2 (14) 图2用流函数导出边界y0和Y1的方程，1|=1 图3勇事E(W)的计算简图 Fig.2 Derivation of the cquations of Fig.3 The calculations of the plastic boundaricsyo and yi power in vi and V2 5功率平衡计算 如果将运动可能的速度场，代入式(1)，就可以计算出全功率E(知)，但它不会小于 零，即E(u,)=E(W)+E(y)+E(y1)+E(S1)-E(S2)≥0 (15) 下面分别计算各项功率，推导过程从略。 5.1变形功率E(V)(见图3) 四=&[店了haa+JHa] (16) 5.2间断功率E(y。),E(y1)参见图2 2。)=可o+++()'() (17) 5
也就是 。 , / ( 。 二 + 。 。 ) = t a u a 。 另一方面 , 沿 y 。 有 t a n a 二 d夕 l d x , 再考虑到式 ( 5 ) , 就得到 : ( a叻/ a x ) d x + ( a吵/ a少 ) d 夕 = 一 。 。 d 夕 此式左端为 功的全微 分 , 可见有 功= 一 ” 。 y + f (二 ) , 利用 条件 州 。 二 。 = 0, 故知 f ( ` ) 二 O , 于 是 : 砂= 一 饥 y = 一 ” ho oy ( 在协上 ) 得人 口 边 界? 。 的方程为 ( 1 2 ) 将式 ( 1 0 ) 代人式 ( 1 2 ) , 梦 刀 告 + 互互… ( 又 . 一 1 2 . . 一 1 y 一 工 ) - 1 = 。 ( 1 3 ) 同理 , 可求得出 口边界 v ; 的 方程为 : 盆 对 会 + 烈买 “ , · . 一 1 2 . ( 妥 歹 一 妥 2 2 . 一 2 夕 夕 1 一 一 : : 1 一 己 ( 1 4 ) ` 已才二二 吮 咭一y 护 八 匕 厂 。 图 2 用 流 函数导出边 界 y o 和 , i 的 方程 , ! , l 二 F 19 . 2 D e r i v a t i o n o f t h e e q u a t i o n s o f b o u n d a r i e , y o a n d 下 i 图 3 功率 E ( V ) 的计 算简图 F 1 9 . 3 T h e c a l c u l a t i o n s o f t h e P l a s t i c P o w e r i n V i 注 n d V Z 5 功率平衡 计算 如 果将运动 可能的速度场 。 , 代 入式 ( 1 ) , 就可以计算出 全功 率 E ( 。 ` ) , 但它不 会 小 于 零 , 即 E ( 。 ` ) = E ( 犷 ) + E ( 夕。 ) + E ( 夕1 ) + E (夕 ; 卜 E ( S , ) ) o 下面分别计算各项 功 率 , 推导过程 从 略 。 ( 15 ) 5 . 1 变形 功率 E ( v ) (见 臼 3 ) 户` F , 一 ` · “ ` [I: 一 ` f {介 d交d夕· 上 一 。 宾 ( 了 ’ “ 翻夕〕 ( 1 6 ) 5 。 2 间断 功率 E (夕。 ) , E ( 夕; ) 参见图 2 ” ` , 。 , = ” 。 《了ha ; ` 不决丁万了万扩裔 d夕 ( 1 7 )
=,”+,加+：+（质）'()。d (18) 5.3摩擦功率E(S1),E(S2) B(s)=1了o+明s:1+4,e运2d (19) s,)=1,r人旺+，1+4a,e10 (20) 6应力场计算 6.1应力式 应力可以从平衡方程和流动方程导出。限于篇幅略去推导过程。从场内任一点(x。,y。) 出发，对上述方程积分，可以求得： =4[(登)-()，]-2，（管）y -2,号（骨）]，.dx+o(x) 0,=0-467 (21) 0,=2骨 6.2关于0：(×0，y0)的确定 为了确定式(21)中的0(x。y。),可将 初始点选在y1上，如图4所示，外力边界条 件要求： XN=-g.sina+o,,cosa 当无轧制张力时，端部无水平力的条件是： ,x=0[ 图4确定却始点(：0，J0)和压力的简图 -()].dy=0 Fig.4 Determination of the initial point on yiand the pressure on s (22) 由此式解出的0：，(x,,y,)必须同时满足后端无外力的条件，它是： [o-,()]y=0 (23) 若此式不能满足，则应在Y1上挪动点(xo,y),重复计算，直至满足，方可使c,(xoyo) 6
1 一 夕 E ( 协 ) = “ 。 、 I 了 〔和 1 + ” 二 ” + ” , , 一 了 d y ( 1 8 ) 摩擦功率 E ( S , ) , 石 ( s : ) 了口O . J J . ` E ( 5 一 1 ) 二 了 ( 。 , + 。 : , · ; 了 1 + ` (` 。 。 ` /` , · d` ( 1 9 ) E ( S : ) = ` I , 二了 ` · , + · , , 5 2 了 ` · “ ` 。 · “ ` , ’ ` 6 应 力 场 计 算 d二 ( 2 0 ) 应 力式 ù合, 吐 .5 · .6 应力可以从平衡方程和流动方程导出 。 限于篇幅略去推导过程 。 从场内任一点 ( ` 。 , y 。 ) 出 发 , 对上述方程积分 , 可以 求得 : 。 二 二 = 城 (合) 一 (合) , 。 -] 2 《 。 斋 (誓) d , _ _ _ J ` 犷 二 二 ” ” 一 U “ 一 任 伟 育洲 口 二 y 二 Z k g 二 , H 一 2 烤 `合’ 〕 ’ 。 d “ “ “ 、 “ ’ 。 ’ { 1 , { ( 2 ” 6 . 2 关于 a 二 二 ( 二 。 , y 。 ) 的确定 为 了确定式 ( 2 1) 中的 ` 二 括 。 里。 儿 可将 初始点选 在 夕 ; 上 , 如图 4 所示 , 外力边界条 件要求 : X 那 二 一 J , 二 s i n a + 『 , , e o s 。 当无轧制张力时 , 端部无 水平力的 条件是 : 夕 P S 骊 a 叭 ; 口 y 丫 劣 J 盆 翻 厂耳IJ I , I X · d , 二 { ( i 一 口 ) h o 0 一 a 二 (装) 〕 , 。 ` y 二 0 图 4 确定初始点 ( 二 。 , , 0) 和压力 户的简图 F i g . 4 D e t e r 位 i 公 a t i o 血 o f t 五e i n i t i a l P o i n t o n y i a n d t h e P r e s s 亡 r e o n , ( 2 2 ) 由此式解出的 口 二 二 (二 。 , y 。 ) 必 须同时满足后端无外力的条件 , 它是 : J h o 仁 。 二 二 一 (筹) 〕 , 。 ` y · 。 ( 2 3 ) 若此式不能满 足 , 则应在 夕; 上挪动点 x( 。 , y 。 ) , 重复计算 , 直至满足 , 方可使『 二 二 (为 , y 。 ) 芍
确定。 6.3轧制单位压力p的计算 轧制单位压力卫由下式计算，见图4。 | p=(0￥sina-0,cosa). (24) 式中“由下式确定： tana=dy/dx=(2hoe/12)x (25) 6计 算实例 下面针对在热轧薄板轧机上某一道次来进行计算。给定变形区参数是五。=6.25mm, 1=44.8mm,e=0.44。按照Sims热轧理论，设接触表面为全制动，即r·=k。变形速度对屈 服应力的影响已考虑在k中，所以可以取”。=1。设定的流函数如式(10)所示，由此流函数 按式(5)导出速度场”：，代入全功率式(1)，也就是式(15)。 如果用变分法求解泛函的极值函数一真速度，则需求解下列M×N个方程的非线性 方程组： aE(u,)/aa.m=1,2,…M,n=1,2,…N 为了避免元长的迭代计算，文中采用优化法。以E=0为目标函数，用Powell法对a.,施行 数值寻优。数值积分采用高斯法，高次方程求根采用牛顿法。求得的口m如下： -19.393 3.430 6.005 0.n=10-2 54.582 -3.186 -3.317 -18.841-11.847 4.279, 这里M×N=9,即逼近空间取为9维。在寻优过程中，B视为一坐标变量。做法是先给定 B为0.5，然后使式(15)的值尽量小，得到E=+6,接着调整B,使E更小。最后求得的结 果是B=0.387。各项功率的数值如表1所示。 表1各项功率的数值 Table 1 The values of power 功 事消耗 功率付出 ×k,mm3/s % xk,mm3/s (W) 1.3507 21.6 E(51) 4.6273 74.2 (y0) 0.0730 1.2 (y2) C.1865 3.0 (S2) 6.2375 果计 6.2375 100.0
确走 。 8 。 3 轧翻单位压力 P 的计算 轧制单位压力 P 由下式 计算 , 见图 4 。 P = ( 叮 二 , s i n a 一 叮 。 , e o s a ) 。 ( 2 4 ) 式 中 a 由下式确定 : t a n a = d y / d 二 二 ( Zh 。 忍 / 1 2 ) 劣 ( 2 5 ) 6 计 算 实 例 下面针对在热轧薄板轧机上某一道次来进行计算 。 给定变 形 区 参 数 是 五。 = 6 。 25 m m , I = 4 4 . 8 m m , 。 二 。 。 4 。 按照 iS m s 热轧理论 ,设接触表面为全制动 , 即了 . = 吞 。 变形速度对屈 服应力的影晌已考虑在 k 中 , 所以可以取 沙 。 二 1 。 设定的流 函数如式 ( 1 0 ) 所示 , 由此流函数 按式 ( 5) 导出速度场 t,l , 代入全功率式 ( 1 ) , 也就是式 ( 1 5 ) 。 如 果用 变分法 求解泛 函的极值函数— 真速度 , 则需 求解下列 M x N 个方程的非线 性 方程组 : a E ( t, ` ) / a a 二 二 = 1 , 2 , 一万 , n = 1 , 2 , … N 为 了避免冗长的迭 代计算 , 文 中采用优化法 。 以 E 二 0为 目标函数 , 用 P o w e n 法对 a 二 施行 数值寻优 。 数值积分采用 高斯法 , 高次方程求根采用牛顿法 。 求得的am 二 如下 : 一 1 9 。 3 9 3 5 4 。 5 8 2 一 1 8 。 8 4 1 3 。 4 3 0 一 3 。 1 8 6 一 1 1 。 8 4 7 6 。 0 0 5 一 3 。 3 1 7 4 。 2 7 9 厂 l e we卜l we、 . 2 n 一 曰 ù 上. 口 口 . 一一 这里 材 ` N 一 ” , 即通近空间取为 ” 维 。 舒 优过程 中 , 夕视为一坐乓变量 · 做法是先 给 定 声为。 。 5 , 然后使式 ( 15 ) 的值尽量小 , 得到 E 二 + 在 , 接着调整 凡 使 E 更小 。 最后求得的 结 果是 声= 0 . 3 8 7 。 各项功 率的数值如表 1 所示 。 表 1 各项功率的橄位 T a b l e i T 五e v a l u e s o f p o w e r 功 率 消 耗 功 率 付 出 x 圣 , m m 3 / s x 圣 , m m 3 s/ E ( V ) E ( 5 1 ) E ( y o ) E ( 护2 ) E ( 5 2 ) 黑 计 1 。 35 0 7 2 1 。 6 曰,自几ó, … 4 五, 3 4 月` 0 62 7 3 0 。 0 73 0 0 。 18 6 5 6 。 2 3 75 6 。 2 3 7 5 1 00 . 0
轧制单位压力的分布见图5。由此图算出的平均单位压力比值p/2k=1,71。用工程中 常用的方法〔3算得的平均单位压力比值是1.75。 Contant pressure Pn/2k=1.71 3 ...Basic flow lines Complete flow lines 1.0 2 0.8 B=voho -0.6 中1=-0.2B 0.4 中2=-0,4A p3=-0.66 0.2 94=-0.8B 0.20.40.60.81.0 00 0.20.4 0.60.81.0 (Exit) x/1 (Entry) 图5轧制单位压力/2业沿接触弧的分布 图6变形区的流线比较 Fig.5.Diagram of distribution of.roll Fig.6 Comparison of basic flow lines pressure and complete flow lines in the plastic zone 0.108 0.0246 0.098 .080 0.090 S 0.0221 0.024 [1/s] 0.0197 0,0221 〔1/s) 0.052 0.043 0.033 0,033 0.043 052 071 061 0.0248 0.0271 图7剪应变速强度H的等高线图 图8应变速串：：的等高线图 Fig.7 Contour line of shear strain-rate Fig.8 Contour line of strain-rate intensity H w总号 0.50 -0.044 0.65 -0.037 0.80 -0.029 -1.11 -0.95 +0.029 〔1/s〕 -0.95- 1.11 0.00 -0.008 0:022 -0.80 +0.021 0.34 -0.65 0.0003 -1.72 图9正应力0x:/2张的等高线图 图10剪应变速率的等高线图 Fig.9 Contour line of normal stress Fig.10 Contour line of shear strain- 0￥*/2k rate 图6上给出了基础场的流线，以及经附加场修正后的流线。每一流线对应着一定的流 量，即流函数值。中心线和接触当然也是流线。图中还给出了计算所得的边界y。和y1。 图7是剪应变速率强度H的等高线图。由图上看出，刷烈的塑性变形发生在人口处的上 方。在整个出口截面以内变形缓慢。在接触弧上中性点附近，几乎不发生变形。 图8和图9是应变速率5.和应力0，：的等高线图。5：在整个变形区内都是拉伸，但应 8
轧制单位压力的分布见图 5 。 由此图算出的平均单位压力比值 p . 12 掩= 1 . 71 。 用 工翟 中 常用 的方法 〔 “ 〕算得 的平 均单位压力比值是 1 . 7 50 3 0 e o nt a n矿p r se s , 么马 . 2 2、六 . 7: } . } / 户 { 入 / } 犷 : 竹 厂 l 人 }仁划 } 火 q 1 . 口 0 。 8 卜 O 。 6 O 。 4 0 . 2 ~ ` _ . B口 s i e fl o w l in es _ C o m醉e t e n o w l i n e s 日= 砂 Q人。 冲1二 一 0 . 2 月 冲2二 一 0 、 4 尽 哗夕二一 0 . 6 口 冲4二一 0 、 8日 乙n 趁望叭 1 0 。 2 0 。 4 0 . 6 (Ex i t ) 劣 / l 0 。 8 1 。 0 0 10 0 二 2 0 。 4 ( E nt r y ) 0 ·旦 万 0 。 8 1 。 0 瞥 图 6 轧制单位 压加 /朴沿接触弧的分布 F万g , 斤 D i a g , a m o r d i s t r i b u t i o 众 o f r o l l P r e s 吕 u r e , 图 6 变形 区的 流线比较 F i g . 6 c o m P a r i , o 公 o f 卜a s i e f l o 双 I i . e s a n d c o 也 P l e t 亡 f l o w l i 血 e s i n t h 亡 p l a s t i e : o n e D 。 1 08 l 。 0 98 0 9 0 / s 〕 图 7 剪应变速率强 度 H 的 等高线 图 气 F 19 . 7 c o 皿 t o u r l i n e o f o b e a r s t r a i n一 r a t e i n t e n s i t y 万 图 8 应 变速串乳 二 的等高线 图 F 19 . 8 c o n t o u r l i n e o f s t r a i n 一 r a t e 占 : 二 十 」 色+ U 一 0 . 04 4 0 3 7 图 9 正 应力 口 二 二 / 2`的等高线 图 F 1 9 . , C o 立 t o u r l i n 亡 o f 几 o r m a l s t , e , , 口二 二 / Zk 图 10 剪应 变速率的等高线图 F i g . x o C o n t o u r 1 1妞 e o f s h e a r s t r a i 红 - r a t e 雪 二 , 图 6 上给出 了基础场的流线 , 以及经附加场修正后 的流线 。 每一流线对应着 一 定 的 流 量 , 即流函数值 。 中心线和接触弧 当然也是流线 。 图中还给出 了计算所得的边界 为 和 夕: 。 图 7 是剪应变速率强度 H 的等高线 图 。 由图上看 出 , 剧烈的塑性变形发生在人 口处的上 方 。 在整个出 口截面 以内变形缓慢 。 在接触弧上 中性点附近 , 几乎不发 生变形 。 图 8 和 图 9 是应变速率氛 二 和 应力q 二 二 的等高线图 。 首 二 二 在整个变形区 内都是拉 伸 , 但应 母
力0：却到处是压缩！图10是剪应变速率5y的等高线图。5y在变形区前部和后部是异号的， 这正和平冲头压入的情况相似。 图11是剪应力0：，的等高线图。图上还给 1.0 出了计算所得的沿弧S剪应力x的分布。显然， 0.0 S上的剪应力条件只是近似地满足了事实上， 0.10.20.30.40.50.60.70.80.91.0 若Y真是在x=B1处突然改变方向，即外力 不连续，也就求不出连续的应力场和速度场。 -0.43 好在现在的分布正符合弧上应有一停滞区的 +0.33 +0.26 3 要求。为使优化出的B值对应于弧上计算剪应 0.29 +0,19 力x的零值点，看来还需提高逼近式的维数。 0.12 0.15 -0.02 8结 论 图11剪应力σ：y/2k等高线和沿弧S剪应力的分布 Fig,11 Contour linc of shear stress (1)用流函数求解轧制向题，需要设定一 ox,/2k and its distribution 个变量B(图1)，原则上没有其它困难。文 along arc s. 中算出的平均单位压力，略低于工程计算的结果，这是必然的、合理的。 (2)求解轧制全功率泛函的极值函数问题，可以作为优化问题处理。多次计算表明，此 泛函是一单峰凸函数。 (3)对于更复杂的情况，例如接触弧上给定混合摩擦条件T=下*(S),或材料为弹塑性 强化型，都可以用流函数求解。 参考文献 1 Ty I H,TeopernyeckHe OcHoB 06pa6oTKu MerannoB aaBnennex.MocKBa: nsa Merannypr,1980,99 2 Nagpal V.Jour.of Engineering for Industry,ASME,1977,99:754 3 Rowe G W,张子公等译。工业金属塑性加工原理。北京：机械工业出版社，1984， 187 符号说明 Yy1一变形区前后边界： e—一道次压缩率， S,S1,S2一接触弧，弧上的前 T—剪应力强度； 滑区和后滑区， h。,1一轧制前后轧件的半高多 H一剪应变速率强度； 1一接触弧的水平投影长度： x,y一相对坐标x=x/1,y=y/h。多 Y,亚一接触真上点的纵坐标，Y=Y1ho:°，一基础流函数和附加流函数： "。,v1一轧件的入口和出口速度多 E一全功率 v,,:一速度场，基础速度 E(S1),E(S2)一在S和S2 场和附加速度场多 上的摩擦功率： 5一应变速率张量场； E(yo,E(y1)一在y和y1上的间断 0:k一应力张量场： 功率： k一屈服剪应力。 卫，P。一轧制单位压力和平均单位压力影 9
力, 朴却到处是压缩! 图10 是剪应变速率氛 , 的等高线图 。 氛 , 在变形区前部和后部是异号的 , 这正和平冲头压入的情况相似 。 户 图 n 是剪应力 丁 二 , 的等高线 图 。 图上还给 出了 计算所得的沿 弧 S 剪应力 , 的分布 。 显然 , 泞上的剪应力 条件只是近似地满足了 `事实上 , 若 , . 真是在 二 = 声l 处突 然改 变 方向 , 即外力 不连续 , 也就求 不出 连续的应力场和 速度场 。 好在现在 , 的分布正符 合弧 上应有一 停滞区的 要 求 。 为使优化出的 声值对应于弧上 计算剪应 力 乍 的零值点 , 看来还需提高逼近式 的维数 。 一 O 。 0 2 \ 8 结 论 ( 1) 用 流函数求解轧制向题 , 需要设定一 个变量声 ( 图 1 ) , 原则 上没 有其它 困 难 。 文 图n 剪应 力叮二 , / 2益等高线和沿弧 s剪应力 的分布 F 1 9 . 1 2 C o n t o u r l i n e o f s h e a r s t r o s s 口二 , / 2 毛 a n d i t s d i s t r i b u t i o n a l o n g a r c 5 . 中算 出的平均单位压 力 , 略低于工程 计算的结果 , 这是必然的 、 合理的 。 (2 ) 求解轧制全功率泛函的极值函数问题 , 可以作为优化问题处 理 。 多次 计算表明 , 此 泛函是一单峰凸 函数 。 ( 3) 对 于更 复杂 的情况 ; 例如接触弧上给定混合靡擦条件 , 一 , * (S ) , 或材料为弹塑性 强化型 , 都可 以用流函数求解 。 参 考 文 蔽 1 ’I y a r H . T e o P e : 班 从 e e K “ e o e 万 0 . 班 0 6 P a 6 o T 盆 皿 M e T a 兀 zJ o s 及 a a 兀 e 且 皿 e 盆 。 M o e ` B a : 皿 3皿 M e T a 几 几 y P r 万 , , i , 80 , 9 9 2 N a g p a l v · J o u r · o f E n g i n e e r i n g f o r l n d “ s t r y 一 A s M E , 1 9 7 7 , 9 9 : , 7 5 4 3 R o w e G W , 张予公 等译 。 工业金属塑性加工原理 。 北京 : 机械工业出版社 , 19 84 , 1 8 7 符 号 说 夕。 vl — 变形 区前后 边界 ; S , S ; , S : — 接触弧 , 弧上的前 滑 区和后滑区 , h 。 , h , — 轧制前后轧件的半高多 I — 接触弧的水平投影 长 度 ; Y , Y — 接触弧上点的纵坐标 , Y 二 Y / h 。 ; 。 。 , Z, ; — 轧 件的入 口和出 口速度 ; ” . , z,7 , “ 言— 速度场 , 基础速度 场和附加速度场多 古 ` * — 应 变速率张量场 , a ` * — 应力张量场 ; k — 屈服 剪应力 。 明 。 — 道次压缩率 ; T — 剪应力强度 ; H — 剪应变速率强度 , 二 , y — 相对坐标 万 = 二 l! , 歹 = y l h 。 多 护 “ , 功 ` — 基础流函数和 附加流函数 , E— 全功率 , E ( S , ) , E ( 5 2 )— 在 S ,和 S : 上的摩擦功率 , 君 伽 。 ) , 君 ( 夕, ) — 在 夕。和 y : 上 的 间断 功率 , P , P 。 — 轧制单位压 力和平均单位压力;