Establishement of Shear Yield Stress Model of Magnetorheological Fluid
-
摘要:
基于磁偶极子理论和相邻链之间的影响,建立磁流变液的双链偶极子微观结构模型,从磁偶极子能量的角度推导了磁流变液的剪切屈服应力方程,采用双链偶极子模型模拟得到剪切角对剪切力的影响；采用双链偶极子模型计算不同质量分数（40%~70%）羰基铁粉磁流变液的剪切屈服应力,分析了磁性颗粒质量分数对磁流变液剪切屈服应力的影响,并与单链偶极子模型计算结果和试验结果进行对比。结果表明：采用双链偶极子模型模拟得到磁流变液的剪切力随着剪切角的增大基本呈先增后减的趋势,在剪切角为15°时达到最大；由单链偶极子模型、双链偶极子模型计算和试验得到的磁流变液剪切屈服应力均随羰基铁粉质量分数的增加而增大,单链偶极子模型和双链偶极子模型计算结果的平均相对误差分别为97.51%和29.58%,双链偶极子模型能较准确地模拟计算磁流变液在工作磁场下的剪切屈服应力。
Abstract:Based on the theory of magnetic dipole and the influence between adjacent chains, a microstructure model of double-chain dipole for magnetorheological fluid was established. The shear yield stress equation of magnetorheological fluid was derived from the perspective of magnetic dipole energy. The influence of shear angle on shear force was obtained by the double-chain dipole model. The shear yield stress of magnetorheological fluid with different mass fractions (40%–70%) of carbonyl iron powder was calculated by the double-chain dipole model. The influence of magnetic particle mass fraction on the shear yield stress of magnetorheological fluid was analyzed , and was compared to the calculation by the single-chain dipole model and the test results. The results show that the simulated shear force of the double-chain dipole model basically first increased and then decreased with the increase of the shear angle, and reached the maximum value at a shear angle of 15°. The shear yield stress of magnetorheological fluid obtained by single-chain dipole model, double-chain dipole model and experiments all increased with the increase of carbonyl iron powder mass fraction. The average relative errors of the calculation of the single-chain and the double-chain dipole model were 97.51% and 29.58%, respectively. The double-chain dipole model could accurately simulate the shear yield stress of magnetorheological fluid under magnetic field.
-
0. 引言
车轴是铁路机车中一个十分重要的构件,世界各国对提高车轴的可靠性均十分重视[1-5]。我国铁道车辆的提速加载对大功率机车及高速列车车轴用钢的性能提出了更高的要求。国内常用的机车车轴材料是40钢和50钢。40钢强度稍低,但韧性好,50钢强度较高,但韧性稍差,2种材料均已无法满足高速、重载铁路机车的要求[6]。EA4T钢是一种广泛使用于国外地铁动车车轴及大功率机车车轴的钢种,也是欧洲标准EN 13261规定的高速客车车轴用钢。之前我国的合金钢车轴主要靠进口,为了快速实现国家铁路机车重载技术和提速战略,通过引进、吸收、消化已实现车轴的国产化,许多新型机车都采用了国产的EA4T钢车轴,如9 600 kW牵引机车和武汉、深圳的地铁车辆等。
EA4T车轴钢的常规热处理工艺为调质热处理（淬火＋高温回火）,根据EN 13261：2009标准要求,调质热处理EA4T钢车轴轴颈1/2半径处的显微组织应全部为马氏体/贝氏体组织（以下简称为M/B组织）。当轴颈直径较小,如180 mm时,采用常规工艺调质热处理后,其1/2半径处可以全部获得M/B组织；但当轴颈直径较大,达到280 mm时,其1/2半径处在常规工艺下获得全部的M/B组织则较为困难[7-8]。在实际生产中为了保证大尺寸车轴淬火加热时心部能够淬透,常常采用提高加热温度（高于铁素体转变为奥氏体的终了温度Ac3 30~50 ℃）的方法进行淬火[9],而回火温度及方式尚需进一步研究。为了确定大尺寸国产EA4T钢车轴的热处理工艺,作者设计了EA4T钢等效车轴的热处理工艺,研究了热处理后等效车轴不同位置的显微组织与性能,据此确定车轴的调质热处理工艺并进行车轴的现场验证。
1. 试样制备与试验方法
试验材料为东北特殊钢集团股份有限公司抚顺特殊钢股份有限公司生产的尺寸280 mm×280 mm×350 mm的热轧退火态EA4T钢车轴毛坯,化学成分如表1所示,符合EN 13261：2009标准要求。在毛坯上加工出规格ϕ280 mm的圆柱体试样,将其作为实际尺寸为ϕ256 mm的EA4T钢车轴的等效车轴,并对其进行调质热处理预试验。
表 1 EA4T钢的化学成分Table 1. Chemical composition of EA4T steel项目 质量分数/% C Si Mn Cr Mo Ni P S V Al Cu Fe 测试值 0.28 0.32 0.70 1.06 0.24 0.17 0.004 0.001 0.04 0.029 0.03 余 标准值 0.22~0.29 0.15~0.40 0.50~0.80 0.90~1.20 0.15~0.30 ≤0.30 ≤0.020 ≤0.015 ≤0.06 ≤0.30 余 根据前期试验结果,EA4T钢的Ac3为840 ℃,在常规条件下其淬火温度应为870~890 ℃。但由于EN 13261：2009标准要求,EA4T钢在淬火时应具有较高的淬透性,而提高奥氏体加热温度可以获得更高的淬透性,故将等效车轴的淬火温度范围设置为890~920 ℃。等效车轴的淬火保温时间可以采用经验公式来估算,公式[10]如下：
(1) 式中：t为保温时间,min；α为加热系数,取值范围为0.9~1.1 min·mm−1；K为加热时的修正系数,取1.2；D为工件的有效厚度,取280 mm。
由式（1）计算得到,等效车轴的淬火保温时间为5~6 h。根据上述分析并结合前期研究[8],确定淬火工艺为900 ℃×5 h。为了确定大尺寸车轴坯的回火工艺,将EA4T钢等效车轴进行900 ℃×5 h水淬处理后,分别进行595,610,650 ℃保温6 h水冷回火处理[8]。
按照EN 13261：2009进行车轴的显微组织及力学性能研究。在不同调质工艺处理后的等效车轴表层、1/2半径处和心部截取金相试样,经打磨、抛光,用体积分数4%硝酸乙醇溶液腐蚀20 s后,用清水冲洗试样,并用乙醇擦拭,再用吹风机吹干,采用NEOPHOT-21型光学显微镜观察显微组织。按照GB/T 228—2002,在车轴表层、1/2半径处和心部位置截取拉伸试样,拉伸试样的尺寸为直径10 mm标准试样尺寸的10倍,在AG-250KNISMO型电子拉压试验机上进行室温拉伸试验,拉伸速度为10 mm·min−1,相同条件下测3次取平均值。按照GB/T 229—2007,在等效车轴表层、1/2半径处和心部分别沿轴向（即横向）和径向（即纵向）截取标准夏比U型冲击试样,在JXB-300型摆锤式冲击试验机上进行室温冲击试验,冲击速度为5 m·s−1,相同条件下测3次取平均值。根据等效车轴测试结果,确定符合标准要求的调质热处理工艺后进行规格ϕ256 mm EA4T钢车轴的现场验证,拉伸试样的尺寸为直径10 mm标准试样尺寸的5倍,冲击试样为5 mm缺口深度的U型冲击试样,测试设备及参数同前。采用Zeiss Supra 55型场发射扫描电镜（SEM）观察拉伸断口形貌。
2. 试验结果与讨论
2.1 调质热处理后等效车轴的组织和性能
由表2可以看出,与EN 13261：2009标准要求的力学性能（屈服强度不低于420 MPa,抗拉强度为650~800 MPa,断后伸长率不低于18%,纵向和横向冲击吸收能量分别不低于50,25 J）相比,等效车轴经900 ℃水淬后,除了595 ℃水冷回火后不同位置的抗拉强度,以及610 ℃水冷回火后表层的抗拉强度偏高之外,其余条件下的强度、断后伸长率、断面收缩率以及纵向与横向冲击韧性均符合要求,同时还存在较大的富裕量。
表 2 等效车轴经900 ℃×5 h水淬和不同温度保温6 h水冷回火后的力学性能Table 2. Mechanical properties of equivalent axle after 900 ℃×5 h water quenching and tempering at different temperatures for 6 h and water cooling回火温度/℃ 位置 屈服强度/MPa 抗拉强度/MPa 断后伸长率/% 断面收缩率/% 冲击吸收能量/J 纵向 横向 595 表层 719.0 866.0 19.0 62.5 156.0 128.0 1/2半径处 603.0 866.0 19.5 68.0 101.0 76.0 心部 694.0 861.0 18.5 62.5 107.0 93.0 610 表层 728.0 843.0 20.5 71.5 129.0 103.0 1/2半径处 628.0 763.0 27.5 65.0 101.0 155.0 心部 543.0 696.0 22.3 65.5 144.0 128.0 650 表层 629.0 781.0 19.0 67.0 141.0 164.0 1/2半径处 628.0 762.0 19.0 72.5 210.0 155.0 心部 595.0 749.0 18.3 66.5 111.0 169.0 由图1可以看出,等效车轴经900 ℃×5 h水淬和650 ℃×6 h水冷回火后,除了心部含有少量铁素体外,表层以及1/2半径处的组织基本为回火M/B组织,符合EN 13261：2009对车轴组织的要求。心部与表层、1/2半径处组织差异的原因主要在于冷却速率的不同,心部冷却速率较慢,未发生完全马氏体转变,导致生成少量铁素体,表层和1/2半径处冷却速率较快,形成了回火M/B组织。这种回火M/B组织具有较高的硬度和耐磨性,对车轴的表层性能有积极影响[11-12]。综上,确定EA4T钢车轴的调质热处理工艺为900 ℃×5 h水淬和650 ℃×6 h水冷回火。
![]()
图 1 等效车轴经900 ℃×5 h水淬和650 ℃×6 h水冷回火后不同位置的显微组织Figure 1. Microstructures of different areas in equivalent axle after 900 ℃×5 h water quenching and 650 ℃×6 h water cooling tempering: (a) surface layer；(b) 1/2 radius position and (c) core2.2 调质热处理工艺的现场验证
在热处理车间采用规格ϕ256 mm的EA4T钢车轴对前文确定的调质热处理工艺（900 ℃×5 h水淬和650 ℃×6 h水冷回火）进行现场验证。由表3可以看出,规格ϕ256 mm的EA4T钢车轴经900 ℃×5 h水淬和650 ℃×6 h水冷回火后,不同位置的力学性能均完全满足EN 13261：2009标准要求。由图2可以看出,不同部位的组织均为回火M/B组织,也满足EN 13261：2009标准要求。
位置 屈服强度/MPa 抗拉强度/MPa 断后伸长率/% 断面收缩率/% 冲击吸收能量/J 纵向 横向 表层 653.0 788.0 21.5 70.0 65.3 72.0 1/2半径处 569.0 717.0 21.5 65.0 80.0 71.3 心部 527.0 695.0 20.0 62.0 72.0 73.3 ![]()
图 2 ϕ256 mm车轴经900 ℃×5 h水淬和650 ℃×6 h水冷回火后不同位置的显微组织Figure 2. Microstructures of different areas inϕ256 mm axle after 900 ℃×5 h water quenching and 650 ℃×6 h water cooling tempering: (a) surface layer；(b) 1/2 radius position and (c) core由图3和图4可以看出,车轴不同位置所取的拉伸试样在断裂前均发生了大量的塑性变形,为韧性断裂。宏观断口表现出明显的缩颈,且仅存在表现韧性的纤维区和剪切唇,而不存在表现脆性的放射区。纤维区微观均呈韧窝特征。在拉伸应力的作用下,试样发生颈缩而在最小截面处形成三维应力,其值在轴线方向上最大,这些三维应力使晶界、缺陷等处形成显微孔洞；随着应力的提高,孔洞不断长大且相互连接,同时产生新的孔洞,从而使裂纹缓慢形成并扩展,最终在断口上留下韧窝状的区域。综上所述,国产EA4T钢车轴经900 ℃×5 h水淬和650 ℃×6 h水冷回火的调质热处理后,其力学性能和组织均符合EN 13261：2009标准要求。
![]()
图 3 经900 ℃×5 h水淬和650 ℃×6 h水冷回火后ϕ256 mm车轴不同位置拉伸试样试验后的宏观形貌Figure 3. Macromorphology of tensile samples in different areas ofϕ256 mm axle after 900 ℃×5 h water quenching and 650 ℃×6 h water cooling tempering after test![]()
图 4 经900 ℃×5 h水淬和650 ℃×6 h水冷回火后ϕ256 mm车轴不同位置拉伸试样的断口SEM形貌Figure 4. Fracture SEM morphology of tensile samples in different areas ofϕ256 mm axle after 900 ℃×5 h water quenching and 650 ℃×6 h water cooling tempering: (a) surface layer, overall morphology；(b) surface layer, micromorphology of fiber area；(c) 1/2 radius position, overall morphology；(d) 1/2 radius position, micromorphology of fiber area；(e) core, overall morphology and (f) core, micromorphology of fiber area3. 结论
（1）国产ϕ280 mm EA4T钢等效车轴在进行900 ℃×5 h水淬和595 ℃×6 h水冷回火处理后不同位置以及610 ℃×6 h水冷回火后表层的抗拉强度均偏高,650 ℃×6 h水冷回火后的强度、断后伸长率、断面收缩率以及纵向与横向冲击韧性均符合EN 13261：2009标准要求,组织也基本为回火M/B组织。确定EA4T钢车轴的调质热处理工艺为900 ℃×5 h水淬和650 ℃×6 h水冷回火。
（2）现场验证得到经900 ℃×5 h水淬和650 ℃×6 h水冷回火后,国产ϕ256 mm EA4T钢车轴表层、1/2半径处和心部处的力学性能和显微组织均符合EN 13261：2009标准要求,不同位置取样拉伸后均发生韧性断裂,断口均由纤维区和剪切唇组成,纤维区呈韧窝形貌。
-
![]()
图 2 单链和双链剪切变形模型
Figure 2. Single-chain (a) and double-chain (b) shear deformation models
![]()
图 4 垂直均匀磁场的边界条件以及双链偶极子模型网格划分
Figure 4. Boundary condition of vertically uniform magnetic field (a) and double-chain dipole model meshing (b)
![]()
图 5 双链偶极子模型磁感应矢量图
Figure 5. Magnetic induction vector diagram of double-chain dipole model
![]()
图 6 双链偶极子模型模拟得到不同剪切角下的磁感应强度云图
Figure 6. Magnetic induction cloud image under different shear angles simulated by double-chain dipole model
![]()
图 7 采用双链偶极子模型模拟得到磁流变液的剪切力随剪切角的变化曲线
Figure 7. Curve of shear force vs shear angle of magnetorheological fluid simulated by double-chain dipole model
-
[1] 杨三锋,黄向明,明阳,等. 新型磁流变-剪切增稠阻尼器的力学模型及试验研究[J]. 机械工程学报,2023,59（16）：418-426. YANG S F ,HUANG X M ,MING Y ,et al. Mechanical model and experimental study of a novel magnetorheological shear thickening damper[J]. Journal of Mechanical Engineering,2023,59（16）：418-426.
[2] 张杰,陈芳,李正贵,等. 高性能磁流变液及其在阻尼器上的应用[J]. 磁性材料及器件,2023,54（3）：113-119. ZHANG J ,CHEN F ,LI Z G ,et al. High performance magnetorheological fluid and its application in dampers[J]. Journal of Magnetic Materials and Devices,2023,54（3）：113-119.
[3] 黄禹铭,吴杰,邓兵兵. 永磁体和线圈叠加下多极式磁流变离合器的设计研究[J]. 现代制造工程,2022（10）：119-124. HUANG Y M ,WU J ,DENG B B. Design and research of multipole magnetorheological clutch with permanent magnet and coil superposition[J]. Modern Manufacturing Engineering,2022（10）：119-124.
[4] 张莉洁,朱德荣,舒云星,等. 磁流变离合器结构设计与分析[J]. 机械设计与制造,2018（11）：88-91. ZHANG L J ,ZHU D R ,SHU Y X ,et al. Structure design and analysis of magneto-rheological fluid clutch[J]. Machinery Design & Manufacture,2018（11）：88-91.
[5] 李伊伦,魏镜弢,吕彤辉,等. 基于磁性复合磨粒的磁流变抛光试验研究[J]. 兵器材料科学与工程,2024,47（2）：1-8. LI Y L ,WEI J T ,LÜ T H ,et al. Experimental study of magnetorheological polishing based on magnetic composite abrasive particles[J]. Ordnance Material Science and Engineering,2024,47（2）：1-8.
[6] 曹顺涛,陈观慈,李明春. 基于永磁体的磁流变抛光励磁装置设计与仿真[J]. 金刚石与磨料磨具工程,2023,43（4）：504-513. CAO S T ,CHEN G C ,LI M C. Design and simulation of magnetorheological polishing excitation device based on permanent magnet[J]. Diamond & Abrasives Engineering,2023,43（4）：504-513.
[7] 王其雨,申玉瑞,金康,等. 光束定向器方位轴磁流变液密封的设计与优化[J]. 润滑与密封,2023,48（4）：127-135. WANG Q Y ,SHEN Y R ,JIN K ,et al. Design and optimization of magnetorheological fluid seal for azimuth axis of beam director[J]. Lubrication Engineering,2023,48（4）：127-135.
[8] 李军强,姜与,郭士杰. 磁流变液动密封装置密封特性研究[J]. 传感器与微系统,2021,40（5）：18-21. LI J Q ,JIANG Y ,GUO S J. Study on sealing characteristics of MR fluid dynamic sealing device[J]. Transducer and Microsystem Technologies,2021,40（5）：18-21.
[9] 杨杨,徐赵东. 考虑磁性颗粒不均匀分布的磁流变液修正微观力学模型及试验验证[J]. 工程力学,2022,39（10）：173-181. YANG Y ,XU Z D. A modified micromodel for magnetorheological fluid considering the uneven distribution of magnetic particles and experimental verification[J]. Engineering Mechanics,2022,39（10）：173-181.
[10] 李小光,佟昱,赵鹏慧,等. 离子液体基低黏高屈服应力磁流变液[J]. 功能材料,2021,52（8）：8188-8191. LI X G ,TONG Y ,ZHAO P H ,et al. Ionic liquid-based magnetorheological fluid with low viscosity and high yield stress[J]. Journal of Functional Materials,2021,52（8）：8188-8191.
[11] 陈亚蒙,罗一平,任洪娟,等. 磁流变液剪切屈服应力计算模型[J]. 轻工机械,2020,38（1）：12-17. CHEN Y M ,LUO Y P ,REN H J ,et al. Shear yield stress model of magnetorheological fluid[J]. Light Industry Machinery,2020,38（1）：12-17.
[12] PENG Y B ,PEI P. Microstructure evolution based particle chain model for shear yield stress of magnetorheological fluids[J]. Journal of Intelligent Material Systems and Structures,2021,32（1）：49-64. [13] NIKA G ,VERNESCU B. Multiscale modeling of magnetorheological suspensions[J]. Zeitschrift Fur Angewandte Mathematik und Physik,2019,71（1）：14. [14] LI H P ,CHEN F ,LIU C H ,et al. Error analysis and optimization of shear yield stress model for magnetorheological fluid[J]. Arabian Journal for Science and Engineering,2019,44（9）：7779-7787. [15] 赵春伟,彭向和. 基于有限元计算的磁流变液特性分析及其与基于磁偶极子模型结果的比较[J]. 功能材料,2012,43（15）：2098-2101. ZHAO C W ,PENG X H. FE computation based modeling for MRFs and comparison with magnetic-dipoles based results[J]. Journal of Functional Materials,2012,43（15）：2098-2101.
[16] DONG X M ,DUAN C ,YU J Q. Research on chain-model transition identification of magnetic dipole theory for magneto-rheological fluid[J]. Advances in Mechanical Engineering,2016,8（12）：1-9. [17] XU Z D ,SUN C L. Single-double chains micromechanical model and experimental verification of MR fluids with MWCNTs/GO composites coated ferromagnetic particles[J]. Journal of Intelligent Material Systems and Structures,2021,32（14）：1523-1536. [18] ZHANG L Y ,LUO Y P ,REN H J ,et al. Microscopic modeling of magnetorheological fluids containing spherical and ellipsoidal ferromagnetic particles[J]. Journal of Magnetism and Magnetic Materials,2023,586：171204. [19] 刘成昊磁流变微观机理与壁面滑移特性研究徐州中国矿业大学2018刘成昊. 磁流变微观机理与壁面滑移特性研究[D]. 徐州：中国矿业大学,2018. LIU C HResearch on MR micro-mechanism and wall slip characteristicsXuzhouChina University of Mining and Technology2018LIU C H. Research on MR micro-mechanism and wall slip characteristics[D]. Xuzhou：China University of Mining and Technology,2018.
[20] ZHANG J T ,SONG W L ,PENG Z ,et al. Microstructure simulation and constitutive modelling of magnetorheological fluids based on the hexagonal close-packed structure[J]. Materials,2020,13（7）：1674. [21] 李昊鹏磁流变液本构模型及剪切传动特性研究徐州中国矿业大学2022李昊鹏. 磁流变液本构模型及剪切传动特性研究[D]. 徐州：中国矿业大学,2022. LI H PResearch on constitutive models and power transmission characteristics of magnetorheological fluid in shear modeXuzhouChina University of Mining and Technology2022LI H P. Research on constitutive models and power transmission characteristics of magnetorheological fluid in shear mode[D]. Xuzhou：China University of Mining and Technology,2022.
下载:
计量
- 文章访问数: 15
- HTML全文浏览量: 7
- PDF下载量: 5
