机械建模与仿真.ppt

建模与仿真 北京科技大学 机械工程学院 张 杰 ZJ_me@ 6.1 板带轧制问题 四辊轧机 轧件 热轧 宽度：700～1900mm 厚度：1.0～25mm 温度：700～1200℃ 冷轧 宽度：700～1900mm 厚度：1.0～25mm 轧辊 工作辊： 直径：500～700mm 长度：1250～2250mm 支持辊 直径：1300～1500mm 长度：1250～2250mm 轧制力：4500t 四辊轧机示意图 6.2 塑性变形问题－轧制力计算 （1）杆件压缩变形载荷计算 条件 弹塑性材料 屈服极限为 σs 强度极限为σb 单向压缩 主应力σ1 方向 σ2 ＝ σ3 ＝0， 面积为S 压力P计算结果 弹性阶段：P＝S ·σ＝ S · Eε 屈服阶段：P＝S ·σs 强化阶段：P＝ ？ 问题—圆柱体的墩粗 摩擦导致的鼓形变形，变形不再是均匀的 变形较大时材料进入强化阶段 二维和三维变形问题 （2）简单轧制过程 两个轧辊都驱动 两个轧辊直径相等 两个轧辊转速相同 轧件做等速运动 轧件上除轧辊施加的力以外，无其它作用力 轧件的机械性能均匀 金属不可压缩：h0b0L0=h1b1L1 轧制力的作用方向 简单轧制：垂直 有张力时：不垂直 （2）基本假设 1）平面变形（应变）假设 描述：轧件只在厚度方向和轧制方向发生变形，宽度无变化 成立理由：轧件宽度大，横向变形阻力大，难变形 2）平面断面假设 描述：由于轧件厚度小，轧前垂直平面在轧中和轧后仍然为平面； 应力特点：忽略切应力； 平行平面上的应力视为主应力； 水平方向应力沿轧件高度均匀分布； 3）变形阻力在轧制区为常数 （3）静力平衡微分方程 取微单元， 在x面上有应力分量 sx ，t xy ， t xz 在x+dx面上，应力分量相对x 截面有一个增量 根据微分单元体x方向平衡，∑Fx=0，有 简化并且略去高阶小量，可得 （4）轧制变形区内的微分方程 前后滑区 前滑区：轧件速度大于轧辊速度 后滑区：：轧件速度小于轧辊速度 在后滑区接触弧上某点（x,y）处取一微分体积 abcd 宽度：dx， 高度：y ~ (y＋dy) 在接触ab 弧上的接触力: 单位压力 px 单位摩擦力 tx （4）轧制变形区内的微分方程 接触弧ab上的合力的水平分量为 微分体内合力的水平分量为 水平方向力的平衡条件 在接触弧上的（x，y）点，有 代入上式并略去二次微分项，得后滑区的静力平衡方程 在前滑区，摩擦力方向与后滑区相反，静力平衡方程为 两式可统一写作 该方程有4个未知数：px，σx，tx，y，如何求解？ 方法： 通过塑性变形条件，寻找 px，σx 之间的关系 通过轧辊与轧件表面的接触性质，确定px， tx 之间的关系 y(x)是接触弧曲线方程 （5）塑性方程（屈服条件） 特雷斯卡（Tresca）屈服条件（最大剪应力理论） σ1 -σ3 = σs （σ1 ≥σ2 ≥σ3 ） 米塞斯（Mises）屈服条件（形状改变比能理论） 两者可以统一写作σ1 -σ3 = βσs —— 塑性方程 当采用米塞斯（Mises）屈服条件时，β＝1～1.15 对于板带轧制，取β＝1.15，故有σ1 -σ3 = 1.15σ σ为变形抗力，与材料、 变形条件（变形温度、变形程度、变形速度）有关 （6）卡尔曼方程 设微单元abcd上的垂直和水平应力就是主应力 在水平方向：σ3 ＝ σx 在垂直方向：σ1 ＝接触弧ab上的合力的垂直分量，即 因第二项摩擦力 tx 的垂直分量比较 px 的垂直分量小得多，可忽略，故有 σ1＝ px 代入塑性方程σ1 -σ3 = βσ ，得 px - σx = βσ ＝ k 二边取微分，有 dpx = dσx 以上给出了px，σx 之间的关系 设压力px与摩擦力tx之间为干摩擦，故可给出px与tx之间的关系 tx =μ px px - σx = βσ ＝ k 将 dpx = dσx