孔型轧制中金属三维变形的模拟

D0I:10.13374/i.issn1001053x.1990.01.008 第12卷第1期 北京科技大学学报 Vol.12 No.1 1990年1月 Journal of University of Science and Technology Beijing Jan.1990 孔型轧制中金属三维变形的模拟 王 进苏逢西 (压力加工系) 摘要：将双流函数法进行推广，设定出了对任意的孔型形状严格满是速度边界条 件的速度场：并用若干相交平面代替人口刚塑性变界面使速度不连续条件得到满足：用罚 函数法使体积不变条件得到满足，找到了一个用上界法模拟复杂孔型中金屈三维变形的一 般途径，提出了一个通用的速度场以5*和2*角铜的切分孔轧制为例，进行了模拟计算和相 应的物理模拟实坠。求出了孔型的充满，变形区的形状，变形区内金属流线的形状，速度 场，应变速率场，应力场以及轧制压力，前滑，延伸等参数。模拟计算值与实验值基本吻 合。这种模拟复杂孔型中金属三维变形的一般方法，为CRD技术的进一步发展打下了基 础。 关键词：孔型轧制，金屈三维变形，金属流动，计算机辅助孔型设计(CARD) R Simulation of Three Dimensional Metal Deformation in Caliber Rolling Wang Jin Su Fengxi ABSTRCT:The dual-stream functions method is improved and applied to establish the velocity field which strictly satisfies velocity boundary conditions of arbitrary caliber configurations.The rigid-plastic interface at the entrance is re- placed by several planes so that the veloclty discontinuity condition at the entr- ance is satisfied.Penalty function algorithm is applied to make the incompressibity condition be satisfied.Then,in complex caliber rolling,a general method for the simulation of three dimensional metal deformation is found and a general velocity field which can be used in arbitrary three dimensional metal flow is formulated. KEY WORDS:caliber rolling,three d'mensicna!metal defcrmation,metal flow,CARD 1988一11-24收稿 43

、 、 第 卷第 期 年 月 北 京 科 技 大 学 学 报 孔型轧制 中金属三维变形的模拟 王 进 苏逢西 压 力加 下 系 摘 要 将双 疏 函 数法 进行推广 ， 设 定 出 了对任意 的 孔 型 形状 严格满足 速度边 界条 件的 速 度 场 多 并用若 千相 交 平 面代 替人 口 刚塑性 变界 面使速度不 连 续 条件得 到满足 多 用 罚 函 数法 使体 积不 变条 件得 到 满足 ， 找到了一个用上 界 法模拟 复杂 孔型 中金属三 维变形 的 一 般途径 ， 提 出 了一个 通 用的 速度 场 。 以 和 手 角钢的 切分 孔 轧制 为例 ， 进行 了 模拟计算和相 应 的 物 理 模拟 实验 。 求 出 了孔 型 的 充满 ， 变形 区 的形 状 ， 变形 区 内 金属流 线 的形 状 ， 速度 场 ， 应 变速率场 ， 应 力 场以 及 轧制压力 ， 前滑 ， 延伸 等 参数 。 模拟 计算 值 与实验值基 本吻 合 。 这 种模拟 复杂 孔 型 中金 属三维 变形的 一般 方 法 ， 为 技 术的 进 一步 发 展 打下 了基 础 。 关键 词 孔 型 轧 制 ， 金 属三 维 变 形 ， 金 属流 动 ， 计 算机辅 助 孔 型 设计 入 牙 儿 ” 劣 卜、 一 · 一 一 宝 ， ， 又 ， 了 。 ， ’ 切 。 毛 力 ， 环 · ， 日 一 一 收 稿 DOI ：10．13374／j ．issn1001－053x．1990．01．008

轧制过程的变形模拟是计算机辅助孔型设计(CARD)技术是否成功的关键之一，但是复 杂孔型中金属的三维变形模拟到现在还没有解决。以往在研究复杂孔型中金属的流动时，常 常将孔型几何形状进行转化，这种转化大大影响了模拟结果的精度。 本文在对任何复杂的孔型形状不进行丝转化的条件下，提出了一个用改进的上界法来 模拟孔型轧制中金属三维变形的-·般途径。这种方法既可以用于CARD技术，又可以用于指 导现场生产。 1改进的上界法及一般的三维变形速度场 上界法的基本原理是刚塑性体的第一变 y Interscetion line of work 分原理，在此基础上Kobayashi【1J进行了改 piece and upper roll 进，指出：使能串泛函取得极值的运动许可 , Upper roll 的速度场对某一给定形状是稳定的，那么这 ( F(x2) 个形状和速度场就是这一包含自由边界问题 (x2) 的解。解的精度和计算时间与所设定的速度 IVk(=) Lower roll 场有关。 G(x.2) 为了使速度场适合于任意孔型形状，设 Intersectlon line I of workpiece and 变形区中坐标系如图1示：其中，x轴与下 lower roll 辊轴线重合，上、下辊面方程分别为y= F(x,2)及y=G(x,2)。 Roll direetion Rolling direction 1.1速度边界条件 图1变形区中坐标系 将V,Nagpal提出的双流函数法【2】进 Fig.1 Coordinate in deformation zone 行推广，应用到型钢轧制中来。先设定一个 流函数X(x,y,2),Y(x,y,2)=(y-G)1(F-G),使上、下辊面同时为流表面，从而使 速度边界条件得到满足；再设定另外一个流函数中(x,y,2)使速度场对设定的变形区左、 右边界W4(x,2)=W:(2)-x=0及Wn(x,z)=WR(2)-x=0是稳定的；即： v.grad(W)=0 在x=W(2)上 (1) v·grad(Wa)=0在x=WR(2)上 其中：v为速度矢量。 这样求出的速度场为： .=[(WR-x)dwv/dz+(x-WL)dw a!dz]v.(x,y,2)/(W w-WL) Uy=-Ux (y-F)OGlx(G-y)oF/0xv(y-F)0Glz+(G-y)oFioz F-G F-G v,=0,(x,y,2) (2) 44

轧 制过 程的变形 模拟 是 计算机辅 助孔 型设计 技 术是否 成功 的关 键之 一 ， 但是 复 杂 孔 型中金属 的三维变形 模拟到 现 在还没 有解决 。 以 往在研 究复杂 孔 ‘ 型 中金属 的流动时 ， 常 常将孔 型 几何 形状 进行转化 ， 这种转化 大大 影响 了 模拟 结果 的 精 度 。 本文 在对 任何 复杂 的孔 型 形状 不进行丝毫转化 的条件下 ， 提 出了 一个 用 改 进 的上界法来 模拟孔 型 轧 制中金 属三维 变形的一搬途径 。 这种 方法既 可以 用于 技术 ， 又可 以用 于指 导 现 场生 产 。 改进的 上 界法及 一 般的 三维变形速度场 上界 法 的基 本原 理是 刚塑性体 的第一 变 分原理 ， 在此 基础 上 了 ’ 进行 了改 进 ， 指 出 使能率 泛 函取 得 极值 的运 动许可 的速度场 对某一 给定 形状 是稳定 的 ， 那 么这 个 形状和 速度场 就 是 这 一 包含 自由边界问题 的解 。 解 的精度和 计 算时 间与所 设定的速度 场有关 。 为 了使速度场 适 合于 任意孔 型形状 ， 设 变形 区中坐标 系如 图 示 其 中 ， 二 轴 与下 辊轴线重 合 ， 上 、 下辊面 方 程 分 别 为 二 ， 及 二 二 ， 的 。 份、 ‘ ， 卜 · 户 是乏 ，， 馆 速 度边 界条 件 将 提 出的双流函数 法 ’ 进 行推广 ， 应 用到型 钢轧 制 中来 。 先 设定一个 图 变形 区 中坐 标系 「 一 若 屯 流 函数 二 ， 夕， “ ， ， ， “ 一 ‘ 一 ， 使上 、 下 辊面 同时为流 表面 ， 从 而 使 速 度边界 条件得到 满 足 再 设定 另 外 一个流 函数 功 ， 夕， 使 速度场 对 设定 的 变形 区左 、 右边 界才 ， 习 附 一 二 。 及 附 。 ， 班 一 二 二 。 是稳定 的 即 石 ‘ 砰 ， ‘ 牙 ， 在 牙 上 在 不于 尹 上 其中 为速度矢 量 。 这样求 出的速度 场为 二 〔 环 二 一 、 否犷 一 不 才 〕 二 二 ， 夕 ， 才 ， 一 不不 一 一 十 一 一 一 二 一 一 一 二 戈 ， ， 之 阳 、尸

1.2速度不连续条件 就孔型轧制的-·般情形而言，变形区的入、出口各是一个速度不连续面。 轧件出口刚塑性交界面上法向速度连续的条件要求：v,l,=o=常数。将Kato提出的v,的 设定方法「81加以改进、推广，设： =,4-(4,+4'()-(÷)(✉)门 A,=a1+a前+a() (3) 4=:+a成+，() 式中，v。入1测刚区速度，。，TL来料最大宽度和高度；y,为轧件与下辊最先接触 点y座标，1最大接触弧长；(x),,(x)分别为入口刚塑性交界面与上、下银交线在水平 面内的没影；A为由体积不变条件决定的常数，a1~a为待定参数。 入口刚塑性交界面的速度不连续条件难以精确满足。为此在之(x),2。(x)上各选3点， 以过某点的微小平面的法向做为该点的法向，在这6点上列出法向速度连续条件所要求满足 的方程，解以？1~a:为未知量的六元一次方程组，求出a1~a。根据Kto【81的观点，该面 上法向速度连续的条件可以认为近似满足了。 1.3体积不变条件 将体积不变条件￡r=0等价地转化为约束： J∬(r1八ar=0 (4) 其中：V为变形区体积，ε为体积应变速率。 将约束(4)下求能率泛函极小值的问题用外罚函数转化为无约束问题，并用Hestenes乘 子法「4防止病态收敛。目标函数为： F=jlo:&dr+j儿：mkAo,ds+4 kdvds j(n1v)dr+(e12)(∬，(Vdv) (5) 其中：c,u为根据Hestenes方法决定的常数。 对上式变分，求出球应力分量： o.=（3是)(e1s(cj∬，(rlV)ndv-) (6) 这样，优化后的速度场就成为运动许可的了。上述方法，提供了一个模拟孔型轧制中金 属三维变形的一般途径，并给出了一个通用速度场。 45

速 度不连 续条 件 就孔 型轧 制 的 一 役情形而 言 ， 变形区 的入 、 出 口 各是一个 速度不连 续面 。 轧 件 出 口 刚塑性 交界面上法 向速度连续 的条件要 求 。 。 ， 。 二 常数 。 将 提出 的公 的 设定方法 〔 “ 加以 改 迸 、 推 广 ， 设 。 。 。 〔 ‘ 一 。 夕 一 夕 。 了 。 万 一 之 一 千 “ 。 工 。 戈 工 、 禹 一 “ ‘ 十 “ 瓦 十 “ “ 又瓦 〕 、 、 过 、 万 十 。 一 昨 。 了 些 、 奋 式 中 ， 。 。 入 汉刚 区速 度 ， 砰 。 ， 。 来 料最 大宽 度和 高度 夕 。 ，为轧件与 下辊最 先接触 点 夕 座标 ， 最大接触弧长 二 ， 、 劝 分 别为入 口刚塑性 交界面与 上 、 下辊 交线在水 平 面 内的 设影 左为 由体积不 变 条件 决定的常数 ， 。 ， 一 。 。 为待定 参数 。 人 口 刚塑性 交界面 的 速度 不连 续 条件难 以精 矶满 足 。 为 此 在之 。 ， 。 幼 上各选 点 ， 以 过某 点的 微小 平面 的 法 向 做为 该点 的 法 向 ， 在 这 点 上列 出 法向速度连续 条件所要 求 满足 的 方程 ， 解 以 。 一 。 。 为未 知 量的六 元 一 次 方程组 ， 求 出 ， 一 。 。 。 根据 。 〔 “ ’ 的观点 ， 该面 上 法向速度连 续 的条件可 以认为近 似满 足 了 。 体 积 不变条 件 将 体积 不变 条 件“ 。 等价地 转化 为约束 二 ‘ 二“ 厂， ‘ ’ 犷 “ 其 中 为 变形区体 积 ， 。 为体 积应 变速 率 。 将 约束 下 求 能率泛 函 极小 值的问题 用 外 罚函数转 化为无约束问题 ， 并 用 乘 子 法 〔 ‘ ’ 防 止病 态 收 敛 。 目标 函数为 饭 、 二 飞 一 犷 阴 “ 一‘ ￡ “ ‘ · ‘ · 一 。 。 二 · 厂 ， ， 厂 ‘ · ， 丁工 ‘ 二 · 犷， ‘ ’ 厂 其 中 。 ， 。 为根据 方法决定 的常数 。 对 上式 变分 ， 求 出球应 力分 量 碑 、 、 久粼 ‘六‘犷， ， 。 ‘六‘犷，“ “ ‘ 一 这 样 ， 优化后的速度场 就成为运 动许可 的 了 。 上 述方法 ， 提供 了一个 模拟 孔 型轧 制中金 属 三维 变形 的 一般 途径 ， 并 给出 了一个 通 用速度场

2角钢切分孔轧制变形区的简化及速度场 作为CARD技术的模拟，主要是模拟孔型的充满，将不充满部分以直线代替，如图2所 示，则不充满以×，与W.,表示。 不考虑轧件的宽展，根据纯几何关系，先求出变形区的形状，做为模拟计算中变形区设 定和简化的依据。 3(g,o (Weg.0) Workpicee Caliber B (g产g 。r) (o,6f》 (x) (x) P (x3g,22g) (0,2) . (,2uD 图2孔型允满 图3变形区形状 Fig.?Filling of caliber Fig.3 Configuration of deformation zone 2.1变形区左、右端面的简化 如图3，求于对称性，左端面： x=W,(z)=0 (7) 参照不考虑宽展时的计算结果，右端面取三阶多项式： x=W.(z)=W(z)=W。十S(22,-2) +(-) +(叠-2”，)- (8) 式中：S,=-dW1d2l,2,其余见图3。 2,2入口刚塑性交界面的简化 由不考虑宽展时的计算结果，设入口刚塑性交界面与上、下辊交线在x02内的投影2.(x) 与26(x)为： 21+22(x-x)0xx 2.(x）= (x (9) 21 0xx 26(x）= 239一261 {21+(i-（x-1)21xW 46

角钢切 分孔 轧制变形 区 的简化及速度场 作为 尺 技 术的 模拟， 主 要 是模拟 孔 型的 充满 ， 将不充满部分以 直 线代 替 ， 如 图 所 示 ， 则不充满 以 ， 与 附 ， 表示 。 不考虑轧 件的 宽展 ， 根 据纯几何关系 ， 定和 简化 的依 据 。 先求 出 变形区的形状 ， 做为 模拟 计算中变形区 设 一 一 以 场，〕 弘 叨 ‘ ，分 甲 甲劣 味溉勺枷 ︸︸妇‘口。 目︸曾 ， ， 井石 吸 沪 图 于 型 充 满 弓 图 变 形 区形 状 厂 变形 区 左 、 右端面 的简化 如 图 ， 求 于对 称性 ， 左端 面 戈 二 牙 二 参照 不 考虑宽展 时 的计 算结 果 ， 右端面取 三阶 多项式 才 之 之 二 。 ， 之 ， 一 之 不 。 ， 一 平 。 ， 之 ‘ 。 一 “ ， 。 详 十 ‘ ，乒 一 － ， 。 一 “ ’ 犷 “ 了一 式中 ， 一 甲 。 ， ， 其余 见 图 。 。 入 口 附塑 性 交 界面 的简 化 尸口 由不 考虑宽 展 时 的计算结果 ， 设入 口 刚塑 性 交界面与 上 、 下 辊 交线在 内的 投影 劝 与 。 二 为 口 一 之 。 ， ， 一 一 玉万厂 气工 一 峋 ’ 一 ” 义 一 火 气 ‘ 之 日 一 ， 一 一矛 义 一 。 一 、 火 ， ， 之 矛‘ 十 ﹄ ﹃ 一 、尹 之 浑 一 工 一 、 。 了 之 ， 一 之 ‘ 。 ， 一 卜 艺 一 。 。 厂 一 羊 一 才 厂白之， 、 几 一︸

式中各符号意义见图3，P、P,为轧件与上、下辊最先接触点。 入口的刚塑性交界面由过P。P:P.P2P,P的4个相交平面代替。 2.3拉缩变形区形状的假定 拉缩变形区不与上辊面接触，设其上端面为一个与x轴平行，过AB与CD的平面，参照 图2,3。其方程为： y=Fp(x,2)=y,-(y.,-y。/I)2 (10) 其中：yg=y.(x.),y。为CD的y座标。 拉缩与接触变形区的界面设为一个与y轴平行过y=F(x,)与上辊面交线的柱面。该 面是一个速度不连续面。 在(8)、(9)、(10)式中包含了6个独立的参数x:,W,,S,21,22,,23；它们 决定了整个变形区的形状，称为儿何参数。 2.4切分孔中的速度场 分别以x,W,S,21,22,2代替(8)、(9)、(10)中的x,W.,Sg,21g,22 ?:然后将(7)、(8)、(9)、(10)代人(2)、(3)可以求出切分孔中的速度场。其中包 括了7个参数，x,W,S,21,22,2a,o,称为速度场参数。 3优化计算 设定了变形区的形状与速度场，可以求出目标函数(5)式，它包括了6个几何参数与T 个速度场参数。对(5)式求极小值，直到达到稳定状态，即x。=x,W。=W,S=S,,21 =21,22=22,23=2g,这时的变形区形状与速度场为真实解【1门。在求目标函数时采用 了高斯数值积分法【5J;在优化计算中采用了Powell-一Sargent.与抛物线插值法。 4物理模拟实验 为了验证模拟研究的合理性，在实验室用铅试样、在红治钢厂用普碳钢分别模拟了2与 5角钢切分孔型的变形，得到了孔型充满x,与W。延伸4，前滑S和轧制压力P。此外，在 实验室实验中，通过在铅试样表面与剖分面上刻网格，考察了变形区中金属的流动与流线形 状。实验中所用材料的变形抗力为： 对工业纯铅：0s=49.3e25e004 (MN/m2) (11) 对锅(BY,F):os=9,8exp(片+a:）(/10)r4〔u6(e/0.4)s-(“。-1)0.4) (MN/m2) (12) 其中：T=(t+273)/1000,t变形温度°C, 47

式 中各 符号意义 见 图 ， 尸 。 、 。 为轧件与 上 、 下辊最 先接触点 。 入 口 的 刚塑 性 交界面由过尸。 尸 尸 。 尸 。 。 的 个相 交 平面代 替 。 拉 缩 变形 区形 状 的假 定 拉缩 变形 区 不与上辊面接触 ， 设其上端面为一个与 轴 平行 ， 过 与 刀 的平面 ， 参照 图 ， 。 其方程为 少 尸 工 ， 夕 ， ， 一 少， ， 一 夕 。 、 、 、 其 中 。 ， 少 “ ， ， 。 为 刀 的 座标 。 拉缩与接触 变形区 的界面设为 一个与 轴 平行过 夕 二 ， ， 习 与上辊面 交 线的柱面 。 该 面 是一个 速度不连 续面 。 在 、 、 式 中包含 了 个独立 的参数 二 。 ， 冲 ， ， ， ， ， ， ， ， ， ， 它 扩 决定 了整个变形 区 的 形状 ， 称 为 几何 参数 。 切 分 孔 中的 速度场 分别以二 。 ， 。 然 后将 括 了 个 参数 ， 班 。 ， ， ， 。 ， 。 ， ， ， 代 替 代人 一 ， ， ， ， 、 中的二 ， ， 才 ， ， ， ， ， ， ， 、 可 以求 出切 分孔 中的速度场 。 其 中包 称 为 速度场 参数 。 优 化 计 算 设定 了变 形区 的形状 与 速度场 ， 可 以求 出 目标 函数 式 ， 它 包括 了 个几何 参数与 了 个 速度场 参数 。 对 式 求 极小值 ， 直到 达到 稳定状 态 ， 即二 一 二 ， ， ， 牙 。 研 。 。 ， ， ， 二 ，， ， ， 习 ， 。 二 。 。 ， 这 时 的变 形区形状与速度场为真实 解 〔 ” 。 在求 目标 函数 时采 用 了高斯数值积分法 〔 ’ 在 优化计算中采 用 了 一 与抛物 线插 值法 。 物理模拟实验 为 了验证 模拟 研 究的 合理性 ， 在实验室 用铅试样 、 在红冶 钢厂 用普 碳 钢分别模拟 了 井与 角钢切 分 孔 型 的变 形 ， 得到 了孔 型充满 ， ，与 附 ， ，， 延伸拼， 前滑 、 和 轧制压 力 此 外 ， 在 实验室实验 中 ， 通过 在铅试样 表面与剖分面上刻 网 格 ， 考察 了变 形 区中金属 的流 动与流 线 形 状 。 实验 中所用材 料 的变形抗 力为 对 工业纯铅 。 。 “ “ 。 。 对 钢 一 一 今 一 二 。 ，一 〔 一 “ ‘ ‘ ， “ 一 一 ‘ ，六〕 其中 了 ， 变形温度 “

41=4.0994,42=-0.5411, 4g=0.1896,44=-0.1107, u5=0.4528,4。=1.5897 5模拟计算结果及其与实验值比较 本研究用FORTRAN-77编制了切分孔模拟的全部程序，在装有8087协处理器的IBM- PC/XT机上对2铅试样与5角钢的切分孔轧制进行了模拟计算。 5,1孔型充满 将优化出的x,,W,值与实验值比较，可知x,最大误差小于15%，W,最大误差小于 10%。 固定坯料，变化辊缝，求出x,,W,随延伸的变化规律如图4。可知对一定的坯料存在 一个“，使孔型顶部充满最好。 5,2变形区内的速度场、应变速事场 将优化参数代入(2)、(3)式，可以求出速度场，进而可以求出应变速串场。模拟出的 变形区内速度分布可以合理地解释实验现象。 图5是与上银接触表面上的轧向速度分布，其中点划线为相应2截面上轧辊速度沿轧向2 的分布。两者的交点O.B~G为相应截面上的 -180 中性点，其连线为中性线，中性线左上方为前 y=d-fld-y)2-z2 0.0 滑区，右下方为后滑区。该图反应了孔型中不 千18.0 均匀变形的情况，还可以定性分析孔型磨损。 -170 -20.0—9 11r 24.5 s9,0 0 -1602 24.0 12.0 23.5 ￥15,0 G 23.0 18.8 -150 22.5 Workpieie 2.0 1.01.11.21.31.41.51.01.7 -Upper 'roll -140 6 81012 Elongation, mm 图4孔型充满陆延伸的变化 图5轧向速度Vz沿宽向x分布(2·角钢，铅试样) Fig.4 Variation of filling of caliber Fig.5 Distribution of velocity of rolling with clongation direction vx in lateral direction (No.2 steel angle,a lead sample) 5.3变形区内应力场 变形区内应力为： e:1十0m01 (13) 48

， 。 ， “ 一 。 ， 。 。 ， ‘ 一 。 。 ， 。 。 模拟计算结果及其与 实验值比较 本研 究用 一 编 制 了切分孔 模拟 的全部程序 ， 在装有 协处 理 器 的 机 上对 雄铅 试 样 与 角 钢 的 切分孔 轧 制 进 行 了 模拟计 算 。 。 孔型 充满 将 优化 出的 ， ， 才 ， ，值与 实验值 比较 ， 可 知 ，最 大 误差小 于 写 ， 砰 ， 最 大 误差 小 于 。 固定坯料 ， 变 化 辊 缝 ， 求 出劣 ， ，， 附 ， ， 随 延伸的变 化 规律如 图 。 可 知对一定的 坯料 存在 一个 产 ， 使孔 型 顶部充满最 好 。 。 变形 区 内的 速 度场 、 应变 速 率场 将 优化 参数 代人 、 式 ， 可 以求 出速度场 ， 进 而 可以求 出应 变速率场 。 模拟 出的 变 形 区内速 度分 布可 以 合理 地 解 释 实验现 象 。 图 是与 上棍接触 表面 上 的轧 向速度分布 ， 其 中点 划线 为相应 截面上轧辊速度沿轧 向 的 分 布 。 两 者 的 交点 为 相应截面上 的 中性 点 ， 其连 线 为 中性线 ， 中性 线 左 上方为 前 滑区 ， 右 下 方为后滑 区 。 该 图反应 了孔 型 中不 均 匀变 形的情况 ， 还可以定性 分 析孔型磨 损 。 戮 一 ，，匕 ‘ 口口尸 、 ，、 、 … 一 砂 声 三的‘ ︸︸ 牙 蕊 一 丫竿 厂 户口， 铲 婚 丫 尸丫 、 曰尸 二二奋 一 一 分 云 从 ， 降 图 孔型 充满随延 伸的 变化 戈 ， 石 图 轧 向速 度犷 沿 宽向 二 分布 角钢 ，铅 走样 ， 。 变形 区 内应 力场 变形 区 内应 力为 一 吐 ￡ ‘ 一 占

其中￡为等效应变速率，0s为变形抗力，0，见(6)式。 2铅试样轧向正应力G,沿宽向的分布，在对称轴处为较强的拉应力，然后减小，在孔型 肩部为较强的压应力，然后又逐新增大，在边部又成为拉应力。0，的这种分布解释了刻槽试 样在对称轴处和边部被期开的现象。 5.4轧制压方 参照图3，轧制压力为： P=2s,w,|,ods (14) 计算与实验结果相比，误差小于1$%。 5.5延伸4与前滑S。 延神值由入轧件面积和模拟的出口轧件面积乙比求出。 型钢轧制中沿宽向冬点前滑值不同，取平均系数： S4=(w,-ur)/vHe)×100i (15) 式中.为轧件出口速度，)r为轧：平均线逃度。 计算与实验结果的误差4小于7%，S小于18%。 5.6变形区中金属的流线 变形区中金属的流线满是下式： (dx(2)id2=v/v. dy(2)/dz=v,/v. (16) （x(2。)=x0 y(2a）=y0 其中：(x,y,2。)为流线上任一点。 变形区内金属流线的求出，采用了定步长的Runge-Kutta方法。计算值与实验值基本 吻合【1。 6结 论 (1)采用双流函数法设定速度场：用若干相交平面代替入口刚塑性交界面使速度不连续 条件近似满足，用罚函数法在优化中使体积不变条件得到满足的一般模拟方法是可行的。 (2)模拟计算结果与实验结果相比误差不大，并可以合理地解释实验现象。 (3)用罚函数法求变形区内的球应力分量0是可行的。 参考文献 1 Oh S I and Kobayashi S.Int.J.Mech.Sci,1975;17:293 2 Nagpal V.J.of Engineering for Industry,Transactions of the ASME, 49

共 中 ￡ 为等效 应 变速率 ， 为 变形抗 力 ， 。 见 式 。 铅试 样轧 向正应 力。 二 沿 宽向的分 布 ， 在对 称轴 处为较强 的 拉应 力 ， 然 后减小 ， 在孔型 肩部为较 强 的 压应 力 ， 然 后又 逐渐增 大 ， 在边部又 成 为拉应 力 。 的这 种分 布解 释 了刻槽试 样在对称轴处 和边 部 被 撕 开 的现 象 。 。 轧 制 压 力 参照 图 ， 轧 制 压 力为 二 一 、 ， ， 。 ， 。 “ 尸 了 · 。 二 ， 二 计算与 实验 结 果相 比 ， 误 差小 于 。 延 伸 与前滑 、 延 伸 谊由人 一 】车 件面 积和 模拟 的 出 口 轧件 面 积之 比求 出 。 型 钢 轧 制 中沿 宽 向 各点 前 滑值不 同 ， 取平均 系数 。 一 口 。 。 少言 式 中。 。 为 礼 件 出 口速 度 ， 。 。 为 轧棍 平均 线速度 。 计 算与 实 验结 果 的 误差 那 小 于 了 ， 小 于 。 叹 变形 区 中金 属的流 线 变形 区 中金属 的流 线满足 下式 ︺一﹄曰二 亡口匕八 二 。 二 义 。 工 。 。 ， ，‘ 。 二 又 。 。 其 中 二 。 ， 。 ， 。 为流 线上任一点 。 变形 区 内金属流 线 的求 出 ， 采 用 了定步长 的 一 方法 。 计 算值 一 与实验值墓 本 吻 合 工“ 。 结 论 采 用双流 函数 法设定 速度场 用若 干相 交平面代 替人 口 刚塑性交界面 使速度不连续 条 件近 似满足 ， 用 罚 函数 法在优化 中使体积不 变条件得到 满 足 的 一般 模拟 方法是可 行 的 。 模拟 计算结 果与 实验结 果 相 比 误差 不大 ， 并可 以 合理地 解 释 实验 现 象 。 用 罚 函数 法求 变形区 内的 球应 力分 量 ， 是可 行 的 。 参 考 文 献 」 关 · ， ，‘ 夕 少

1977,99:624 3 Kato K,Murota T and Kawabata K.J.Jap.Soc.Tech.Plasticity, 1985;26:174 4席少霖，赵凤治.最优化计算方法，上海科学技术出版社，1984 5李庆杨等，数值分析，武汉，华中工学院出版社，1984 6王进.北京钢铁学院硕土学位论文，1987 洲水xw以*x※Nw以wwww州州州w州%xXXX以X州x 冶金人才数据管理系统 为了更好地适应冶金生产发展和技术进步的需要，进行治金专门人才的需求跟踪顶测， 开发了冶金人才数据管理系统。系统的主要功能： (1)对原始数据的收集、整理、校对、检索、修改、统计、汇总、计算和存贮等。 (2)提供治金系统13种企业的68个单位，连同治金部本身在内的69个单位，整体的和抽 样单位的基本信息。 (3)提供各单位各企业类型按年度统计的专业结构学历分布数据。 (4)提供各单位各企业类型按年度统计的专门人才的学历结构数据。 (5)提供各单位或各企业类型专门人才各学历的相对比例数。 (6)提供各种输出报表。 主要特点： (1)汉字菜单选择和操作提示，提供了良好的用户友好性。 (2)作为汉字数据管理系统，无论什么操作都可不必输入-·个议字。方便用户，加快了 操作速度，避免了操作错误。 (3)提供了3种给入方式，为成批数据的输入提供了方使。 (4)多数据库管理，加快库中的操作。 (⑤)采用了编译系统，增加了系统的安全性和保密性。 (6)提供了两级保密，加强了系统的安全性和保密性。 (7)考虑了系统数据的完整性，避免数据兀余，增强了系统的完整性。 综上所述，该系统选题有-定的社会实用价值，总体设计是成功的，实现法是可行 的，逻辑和物理设计是合理的，设计思想和采用的技术比较先进，为进一步研制决策支持系 统打下了基础。 *068601660660666603W州wX州00:6664986266806840668010X-010060800108 50

， ’ ， · · · · 少 ， 席 少 霖 ， 赵凤 治 最 优化计 算 方 法 ， 上海科 学技 术 出版社 ， 李庆杨等 数值分析 ， 武汉 ， 华 中工学 院 出版 社 ， 王 进 北 京钢铁学院硕 士 学 位论文 ， 冶 金 人 才 数 据 管 理 系 统 为 了更 好地 适应 冶 金 生 产发 展 和 技 术进步 的需要 ， 进 行冶 金 专门人 才 的 需 求跟 踪 预 测 ， 开发 了冶 金人 才数 据管理 系统 。 系 统 的 主要 功能 对 原始数据的收集 、 整理 、 校 对 、 检索 、 修 改 、 统计 、 汇总 、 计算和 存贮 等 。 提供 冶 金 系统 种企业 的 “ 个 单位 ， 连 同冶 金部本身 在内的 个单位 ， 整体 的和 抽 样单位 的基 本 信息 。 提供 各单位各 企 业 类型按 年度统计的 专业结 构 学历分 布数据 。 提 供 各 单位各 企 业 类型按 年度 统 计的 专门人 才的学历 结 构数据 。 提 供 各单位或 各企业类 型 专门人 才各 学历 的相 对 比例数 。 〔 提 供各 种输 出报表 。 主要 特 点 汉字菜单选择和 操作提 示 ， 提 供 了良好的用 户 友好性 。 作为 汉 字数据 管理 系统 ， 无论 什么操 作 都 可不必输 入 一 个 汉 字 。 方便 用户 ， 加快 了 操作速 度 ， 避 免 了 操作错 误 。 提 洪了 种输 入方 式 ， 为 成批数 据 的输 入提 供 了方便 。 多数据 库管理 ， 加快库中的操作 。 的 采 用 了编译系 统 ， 增加 了 系统 的 安全性和 保密性 。 提 供 了两 级保 密 ， 加强 了 系统 的 安全性和 保密性 。 考虑了 系统数 据 的完整性 ， 避免数据冗 余 ， 增强 了系 统 的完整性 。 综 上 所 述 ， 该系统 选题 有一定 的 社 会实用价 值 ， 总 沐 设 计是 成 功的 ， 实现 方 法 是 可 行 的 ， 逻 辑和 彻理 设计是 合理 的 ， 设 计思想 和 采用 的技 术 比较 先迸 ， 为 进 一步 研 制决策 支 持系 统 打 下 了基 础