CN107526854A - 一种三辊行星轧机轧制力的计算方法 - Google Patents
一种三辊行星轧机轧制力的计算方法 Download PDFInfo
- Publication number
- CN107526854A CN107526854A CN201610459111.8A CN201610459111A CN107526854A CN 107526854 A CN107526854 A CN 107526854A CN 201610459111 A CN201610459111 A CN 201610459111A CN 107526854 A CN107526854 A CN 107526854A
- Authority
- CN
- China
- Prior art keywords
- rolling
- roll
- force
- roller
- rolled piece
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000005096 rolling process Methods 0.000 title claims abstract description 115
- 238000000205 computational method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000009826 distribution Methods 0.000 claims abstract description 19
- 239000000463 material Substances 0.000 claims description 10
- 238000005520 cutting process Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 abstract description 8
- 238000004458 analytical method Methods 0.000 abstract description 7
- 230000008859 change Effects 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 6
- 229910000831 Steel Inorganic materials 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 229910052802 copper Inorganic materials 0.000 description 2
- 239000010949 copper Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 229910000881 Cu alloy Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005242 forging Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000003913 materials processing Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Control Of Metal Rolling (AREA)
Abstract
本发明涉及三辊行星轧机领域，具体为一种三辊行星轧机轧制力的计算方法。在确定一定轧制工艺参数下一个轧辊与轧件接触区周向分布范围夹角的基础上，计算轧制变形接触区的面积，建立了轧制力的解析模型和计算方程组，通过数值解析获得轧制变形区轧制力的分布，进而获得轧制力。本发明方法可以基于轧辊辊型曲线和轧制基本工艺条件，包括轧辊的偏转角、倾斜角、轧制材料强度和摩擦系数等，定量计算轧件在轧制区内的轧制力变化，较好地反映了工艺条件对轧制力的影响规律。
Description
技术领域
本发明涉及三辊行星轧机领域，具体为一种三辊行星轧机轧制力的计算方法。
背景技术
三辊行星轧机发明之初是用于钢棒材的轧制，后期也向钢管材扩展，可以达到高效低耗短流程的目的。自上世纪九十年代，芬兰奥托昆普公司将三辊行星轧机引入铜管材的连续生产中，显著提高了加工效率，奠定了当前世界范围内主流的铸轧法工艺(“水平连铸-行星轧制”)基础，有关三辊行星轧机结构、运动特点、力能参数、轧件组织等研究不断深入。
轧制力作为轧机的关键参数，对于三辊行星轧机工作过程中承力部件(如：齿轮、轴承等)的稳定性、使用寿命有重要影响，通过建立轧机轧制力的解析模型，分析结构与工艺参数对于轧制力的影响规律，对于轧机的设计与优化具有重大的应用价值。
由于轧件受力状况的复杂性，而且轧辊与轧件接触区域不是规则形状，加上整个轧制过程的压力与摩擦力的变化，使得对于轧制力的计算需要做各种简化与近似。这包括接触面形状的规则化、接触压力的均匀分布等，在简化计算的同时带来较大的误差。文献1[康祖立，秦念祖，项德。行星斜轧机特征参数的理论分析及轧制力的确定。北京钢铁学院学报，1984(s1)：1-10]通过图解法给出了轧制时轧辊与轧件的接触面积的计算方法，乘以单位平均压力后可以得到总的轧制力。但由于三辊行星轧机轧制空间关系的复杂性，三维转二维的图解的过程中有简化，且力的作用区域与方向难以确定。文献2[黄勇，三辊行星轧机轧制过程力分析。昆明理工大学硕士学位论文，2008]和文献3[孙明伟，三辊行星轧机的辊形研究。昆明理工大学硕士学位论文，2012]，利用单位压力微分求解原理推导了总体轧制力的计算公式，并分析计算了一组轧辊设计方案条件下的轧制力和轧制功率的数值。更多的轧制力计算采用了数值模拟的方法，如文献4[Shih C K，Hung C H.Experimental andNumerical Analyses on Three-roll Planetary Rolling Process.Journal ofMaterials Processing Technology,2003,142(3)：702-709]、文献5[李冰、李章刚、张士宏等，铜合金管坯旋轧成形的三维热力耦合有限元模拟。科学技术与工程，2005，17：1294-1297]和文献6[陈生平，铜管三辊行星轧制过程轧制力及轧辊轴力学行为分析。山东大学硕士学位论文，2014]。
发明内容
本发明的目的在于提出一种三辊行星轧机轧制力的计算方法，通过数值解析可以获得轧制变形区轧制力的分布，进而获得轧制力。
为实现以上目的，本发明的技术方案是：
一种三辊行星轧机轧制力的计算方法，包括如下步骤：
步骤一、计算给定轧制工艺参数下一个轧辊与轧件接触区周向分布范围的夹角；
步骤二、根据轧辊辊型曲线确定轧辊辊面展开扇形面积；
步骤三、利用切块法建立轧制力的解析方程组：
式中，r和为柱坐标系下的坐标量，θout为轧出线的周向角，ρ表示轧辊外轮廓上各点到轧件轴线(x轴)的距离；α为偏转角；β为倾斜角；A为行星轮公转半径，k为轧制材料强度；μ为轧辊与轧件间的摩擦系数，pθ为轧制力；
步骤四、给定轧制材料强度k和摩擦系数μ后数值求解以上方程组，获得轧制力pθ的分布情况；
步骤五、取轧制变形区的轧制力平均值作为所求结果。
本发明的设计过程如下：
基于文献7[刘海、张士宏、程大勇等，三辊行星轧制中轧辊与轧件接触区分布。锻压技术，2011，36(3)：46-49]对三辊行星轧制中轧辊与轧件接触区分布的研究，本发明理论推导了轧制力的解析模型。采用该模型，利用本发明方法可以获得轧制力分布的解析公式，进而采用数值方法求解。
一、参照文献7的方法计算接触区的咬入线与轧出线分布，确定给定轧制工艺参数下一个轧辊与轧件接触区周向分布范围的夹角ψ。
二、轧辊辊面展开图形为扇形。根据轧辊辊型曲线确定扇形的两组半径值r1(mm)和r2(mm)。计算轧制变形接触区的面积S为
三、利用切块法进行轧制力的解析推导，如附图1所示。x方向平衡方程为：
式中，σθ+dσθ.为作用在切块左侧x方向上的应力分量(MPa)；(y+dy)·dz为切块左侧的投影面积(mm2)；σθ为作用在切块右侧x方向上的应力分量(MPa)；y·dz为切块右侧的投影面积(mm2)；为法向轧制力在y方向上的分量(MPa)；为切向摩擦力在y方向上的分量(MPa)；为轧制力在切块在x方向上的投影面积(mm2)。
其中，
整理得
σθ·dy+y·dσθ+dσθ·dy-pθ·dy+tθ·dx＝0 (4)
略去二次微分项，则有
对于前滑区，摩擦力方向相反，则有
且摩擦力tθ(MPa)大小为：
tθ＝μ·pθ(7)
平面问题条件下的屈服条件为：
σ1-σ3＝σs (8)
式中，σ1为第一主应力；σ3为第三主应力；σs为许用应力；
轧制过程中，单元体上的垂直应力为第一主应力σ1，水平应力为第二主应力σ2，轴向(即z方向)应力为第三主应力σ3。而平面问题中，三个主应力满足如下关系：
联立式(8)和式(9)得
根据附图1所示可得
σ2＝σθ (12)
式(11)中第一分项为轧辊的压力项，第二分项为摩擦力项。摩擦力在轧制中主要起到咬入的作用，相比于压力要小得多，可忽略。则
σ1＝pθ (13)
根据式(10)、(11)和(12)得
式(14)和(7)代入式(5)和(6)得
在式(15)中表征了轧件的减径变形量。根据附图1与文献7中图1、图2所示坐标系存在如下关系：
其中r和为柱坐标系下的坐标量(mm)，θout为轧出线的周向角(°)，可由文献7中式(5)求解，即
因此，文献7中的式(4)就可以写为
式中，ρ表示轧辊外轮廓上各点到轧件轴线(x轴)的距离(mm)；α为偏转角(°)；β为倾斜角(°)；A为行星轮公转半径(mm)。
式(15)可以写为
综上，联立方程组如下：
根据边界条件
代入式(14)得
这样，轧制区轧制力pθ满足以下方程组：
式中，k为轧制材料强度(MPa)；μ为轧辊与轧件间的摩擦系数，pθ为轧制力(MPa)。
四、给定轧制材料强度k和摩擦系数μ值后通过数值方法求解方程组(22)，可以获得轧制力的分布情况。
五、取轧制变形区的轧制力平均值作为所求结果。
本发明的优点及有益效果是：
本发明基于轧辊辊型曲线和轧制基本工艺条件包括轧辊的偏转角、倾斜角、轧制材料强度和摩擦系数等，就可以定量计算轧件在轧制区内的轧制力变化。
附图说明
图1为本发明的三辊行星轧制过程微元受力分析。
图2实施例1中的轧辊辊型曲线。其中，横坐标代表轧辊轴向坐标，纵坐标代表轧辊周向轮廓线坐标。
图3实施例2中的轧辊辊型曲线。其中，横坐标代表轧辊轴向坐标，纵坐标代表轧辊周向轮廓线坐标。
具体实施方式
下面结合实施例对本发明做进一步说明。
实施例1
轧辊辊型轮廓如图2所示，曲线方程如下：
一、参照文献7的方法计算接触区的咬入线与轧出线分布，确定给定轧制工艺参数下一个轧辊与轧件接触区周向分布范围的夹角ψ＝22°。
二、轧辊辊面展开图形为扇形。根据轧辊辊型曲线确定扇形的两组半径值r1＝127mm和r2＝95mm，计算轧制变形接触区的面积为1363.87mm2。
三、利用切块法建立轧制力的解析方程如下：
四、给定轧制材料强度k＝58MPa和摩擦系数μ＝0.5值后通过数值方法求解方程组(22)，获得轧制力pθ的分布情况。
五、取轧制变形区的轧制力平均值66MPa作为所求结果。
实施例2
轧辊辊型轮廓如图3所示，曲线方程如下：
一、参照文献7的方法计算接触区的咬入线与轧出线分布，确定给定轧制工艺参数下一个轧辊与轧件接触区周向分布范围的夹角ψ＝25°。
二、轧辊辊面展开图形为扇形。根据轧辊辊型曲线确定扇形的两组半径值r1＝135mm和r2＝28mm，计算轧制变形接触区的面积为3805mm2。
三、利用切块法建立轧制力的解析方程如下：
四、给定轧制材料强度k＝170MPa和摩擦系数μ＝0.3值后通过数值方法求解方程组(22)，获得轧制力pθ的分布情况。
五、取轧制变形区的轧制力平均值210MPa作为所求结果。
实施例结果表明，本发明在确定一定轧制工艺参数下一个轧辊与轧件接触区周向分布范围夹角的基础上，计算轧制变形接触区的面积，建立了轧制力的解析模型和计算方程组，通过数值解析获得轧制变形区轧制力的分布，进而获得轧制力。本发明方法可以基于轧辊辊型曲线和轧制基本工艺条件，包括轧辊的偏转角、倾斜角、轧制材料强度和摩擦系数等，定量计算轧件在轧制区内的轧制力变化，较好地反映了工艺条件对轧制力的影响规律。
Claims (1)
1.一种三辊行星轧机轧制力的计算方法，其特征在于，包括如下步骤：
步骤一、计算给定轧制工艺参数下一个轧辊与轧件接触区周向分布范围的夹角；
步骤二、根据轧辊辊型曲线确定轧辊辊面展开扇形面积；
步骤三、利用切块法建立轧制力的解析方程组：
式中，r和为柱坐标系下的坐标量，θout为轧出线的周向角，ρ表示轧辊外轮廓上各点到轧件轴线x轴的距离；α为偏转角；β为倾斜角；A为行星轮公转半径，k为轧制材料强度；μ为轧辊与轧件间的摩擦系数，pθ为轧制力；
步骤四、给定轧制材料强度k和摩擦系数μ后数值求解以上方程组，获得轧制力pθ的分布情况；
步骤五、取轧制变形区的轧制力平均值作为所求结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610459111.8A CN107526854B (zh) | 2016-06-22 | 2016-06-22 | 一种三辊行星轧机轧制力的计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610459111.8A CN107526854B (zh) | 2016-06-22 | 2016-06-22 | 一种三辊行星轧机轧制力的计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107526854A true CN107526854A (zh) | 2017-12-29 |
| CN107526854B CN107526854B (zh) | 2020-08-21 |
Family
ID=60733735
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610459111.8A Active CN107526854B (zh) | 2016-06-22 | 2016-06-22 | 一种三辊行星轧机轧制力的计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107526854B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109376444A (zh) * | 2018-11-07 | 2019-02-22 | 太原科技大学 | 三辊孔型轧制棒材横断面面积计算方法 |
| CN113579129A (zh) * | 2021-07-29 | 2021-11-02 | 武汉理工大学 | 基于分割累加的复杂异型截面环件轧制力计算方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101062509A (zh) * | 2006-04-29 | 2007-10-31 | 中国科学院金属研究所 | 一种三辊行星轧机轧辊参数自动计算方法 |
| CN101201871A (zh) * | 2007-12-18 | 2008-06-18 | 东北大学 | 一种预测热轧过程轧制力的刚塑性有限元方法 |
| CN101927267A (zh) * | 2009-06-22 | 2010-12-29 | 宝山钢铁股份有限公司 | 一种精轧带钢辊缝的控制方法及装置 |
| CN102303051A (zh) * | 2011-06-15 | 2012-01-04 | 天津商业大学 | 使用七机架连轧机轧制无缝钢管的管形控制方法 |
| CN103191919A (zh) * | 2012-01-05 | 2013-07-10 | 鞍钢股份有限公司 | 一种带钢轧制在线控制摩擦系数模型优化方法 |
-
2016
- 2016-06-22 CN CN201610459111.8A patent/CN107526854B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101062509A (zh) * | 2006-04-29 | 2007-10-31 | 中国科学院金属研究所 | 一种三辊行星轧机轧辊参数自动计算方法 |
| CN101201871A (zh) * | 2007-12-18 | 2008-06-18 | 东北大学 | 一种预测热轧过程轧制力的刚塑性有限元方法 |
| CN101927267A (zh) * | 2009-06-22 | 2010-12-29 | 宝山钢铁股份有限公司 | 一种精轧带钢辊缝的控制方法及装置 |
| CN102303051A (zh) * | 2011-06-15 | 2012-01-04 | 天津商业大学 | 使用七机架连轧机轧制无缝钢管的管形控制方法 |
| CN103191919A (zh) * | 2012-01-05 | 2013-07-10 | 鞍钢股份有限公司 | 一种带钢轧制在线控制摩擦系数模型优化方法 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109376444A (zh) * | 2018-11-07 | 2019-02-22 | 太原科技大学 | 三辊孔型轧制棒材横断面面积计算方法 |
| CN109376444B (zh) * | 2018-11-07 | 2022-04-19 | 太原科技大学 | 三辊孔型轧制棒材横断面面积计算方法 |
| CN113579129A (zh) * | 2021-07-29 | 2021-11-02 | 武汉理工大学 | 基于分割累加的复杂异型截面环件轧制力计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107526854B (zh) | 2020-08-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103302094B (zh) | 一种横向楔形轧制变厚度钢板的生产方法 | |
| Pater et al. | Complex numerical analysis of the tube forming process using Diescher mill | |
| CN106903166A (zh) | 一种铝合金板材异步轧制翘曲预报和优化的方法 | |
| CN107526854B (zh) | 一种三辊行星轧机轧制力的计算方法 | |
| Joun et al. | An axisymmetric forging approach to preform design in ring rolling using a rigid–viscoplastic finite element method | |
| CN105458012A (zh) | 一种与高次曲线工作辊配合的支撑辊辊形通用设计方法 | |
| CN104438493B (zh) | 等厚度曲面的连续辊压成形方法 | |
| Jiang et al. | Rotation mechanics and numerical simulation of hot rolling process under asymmetric rolls | |
| CN105013835B (zh) | 冷连轧机组极薄带轧制中基于热凸度的原始辊缝设定方法 | |
| Wang et al. | Investigation and improvement of 3D rolling process for 3D surface parts | |
| CN102553945A (zh) | 一种适合于四辊轧机的非常态板形预报方法 | |
| Lee et al. | Control method for centering rolls in radial-axial ring rolling process | |
| Shatalov et al. | SHEET MILL CONTROL IN STEEL STRIP HOT ROLLING. | |
| Han et al. | An innovative eccentric ring rolling method for fabricating eccentric rings | |
| Sun et al. | Study on asymmetrical deformation and curvature of heavy cylinder rolling | |
| Sun et al. | Numerical and experimental investigations on the effect of mandrel feeding speed for high-speed rail bearing inner ring | |
| Wang et al. | Analysis of metal forming in two-roll cross wedge rolling process using finite element method | |
| WO2018163930A1 (ja) | クロス角同定方法、クロス角同定装置、及び圧延機 | |
| CN202356455U (zh) | 多辊矫直机辊系辊缝调整装置 | |
| Cui et al. | Modelling analysis of coupling deformation between strip steel and roller system | |
| Zhang et al. | Research on wall thickness uniformity of hollow axles by three-roll skew rolling | |
| Soszyński et al. | A review of contemporary solutions for cold rolling that allow quality improvement | |
| JP4603707B2 (ja) | 継目無管の製造方法 | |
| CN105880299B (zh) | 一种确定冷轧机工作辊水平移动距离的方法 | |
| Guo et al. | Numerical modelling and simulation of radial-axial ring rolling process |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |