韩晨怡1,倪凯旋1,马辉1,2,王海3,林君哲1
(1. 东北大学 机械工程与自动化学院,沈阳 110819; 2.航空动力装备振动及控制教育部重点试验室(东北大学), 沈阳 110819;3. 三一重型装备有限公司,沈阳 110027)
摘要:
为准确分析液压支架在顶梁扭转加载试验时的强度特性,以ZY6800/08/18D型液压支架为研究对象,采用整机建模的方式,建立了液压支架的三维实体有限元模型. 根据液压支架试验的测试结果,选取仿真模型上的测点位置,分析对比顶梁和底座在8种不同的约束边界条件下,液压支架的应力分布情况,同时分析带摩擦的接触边界条件下底座和平台的接触应力分布及接触状态,指出顶梁扭转加载时的载荷不对称导致的底座不完全接触的现象. 在不改变顶梁和底座约束边界的条件下,分析4种销轴连接简化方法对液压支架整体应力分布的影响. 对比仿真结果与测试结果,分别评估各种边界条件和不同销轴连接简化方法的准确性,以及对仿真效率的影响. 结果表明,液压支架的顶梁边界选择固定支撑,底座边界选择与工作台摩擦接触,保留销轴并将销轴与轴孔绑定接触得到的仿真结果最接近试验结果,且有较高的计算效率. 研究结果可为液压支架设计过程中准确预估应力分布提供依据和参考.
关键词: 液压支架 强度特性 有限元 约束边界 销轴简化
DOI：10.11918/202004069
分类号:TD355
文献标识码:A
基金项目:辽宁省“兴辽英才计划”项目(XLYC1807008)
Effects of constraint boundaries and simplified pin joints on strength of hydraulic support
HAN Chenyi1,NI Kaixuan1,MA Hui1,2,WANG Hai3,LIN Junzhe1
(1. School of Mechanical Engineering and Automation, Northeastern University, Shenyang 110819, China; 2. Key Laboratory of Vibration and Control of Aerodynamic Equipment (Northeastern University), Ministry of Education, Shenyang 110819, China; 3. Sany Heavy Equipment Co., Ltd., Shenyang 110027, China)
Abstract:
To accurately analyze the strength characteristics of hydraulic support considering the canopy torsional loading condition, by taking ZY6800/08/18D hydraulic support as the research object, the three-dimensional solid finite element model is established. According to the results of hydraulic support test, selecting the location of measuring points on the simulation model, the stress distributions of the hydraulic support are compared under eight boundary conditions for the canopy and base. In addition, the contact stress and status are analyzed under the frictional contact boundary condition between the base and worktable. It is pointed out that the phenomenon of incomplete contact of the base is caused by the load asymmetry during the canopy torsional loading. The influence of four simplified methods for pin joint on the stress distribution is also analyzed under the constant constraint boundary condition. Finally, the stress distribution obtained from the experiment evaluates the accuracy and efficiency of various boundary and pin joint simplification methods. When the canopy boundary selects fixed support, the base boundary selects friction contact with the working table, and retaining the pin joint and bonding the pin joint to the axle hole, the simulation results are the closest to the test results, and the simulation has high-efficiency. The research results can provide the basis and reference for the accurate estimation of the hydraulic support stress distribution.
Key words: hydraulic support strength characteristics finite element constraint boundary pin joint simplification
韩晨怡, 倪凯旋, 马辉, 王海, 林君哲. 约束边界和销轴简化对液压支架强度影响[J]. 哈尔滨工业大学学报, 2021, 53(1): 78-86. DOI: 10.11918/202004069.
HAN Chenyi, NI Kaixuan, MA Hui, WANG Hai, LIN Junzhe. Effects of constraint boundaries and simplified pin joints on strength of hydraulic support[J]. Journal of Harbin Institute of Technology, 2021, 53(1): 78-86. DOI: 10.11918/202004069.
基金项目 辽宁省“兴辽英才计划”项目(XLYC1807008) 作者简介 韩晨怡(1996—)，男，硕士研究生;
马辉(1978—)，男，教授，博士生导师 通信作者 马辉，mahui_2007@163.com 文章历史 收稿日期: 2020-04-13
Abstract Full text Figures/Tables PDF
约束边界和销轴简化对液压支架强度影响
韩晨怡1, 倪凯旋1, 马辉1,2
, 王海3, 林君哲1 1. 东北大学 机械工程与自动化学院，沈阳 110819;
2. 航空动力装备振动及控制教育部重点试验室(东北大学)，沈阳 110819;
3. 三一重型装备有限公司，沈阳 110027
收稿日期: 2020-04-13
基金项目: 辽宁省“兴辽英才计划”项目(XLYC1807008)
作者简介: 韩晨怡(1996—)，男，硕士研究生; 马辉(1978—)，男，教授，博士生导师
通信作者: 马辉，mahui_2007@163.com
摘要: 为准确分析液压支架在顶梁扭转加载试验时的强度特性，以ZY6800/08/18D型液压支架为研究对象，采用整机建模的方式，建立了液压支架的三维实体有限元模型.根据液压支架试验的测试结果，选取仿真模型上的测点位置，分析对比顶梁和底座在8种不同的约束边界条件下，液压支架的应力分布情况，同时分析带摩擦的接触边界条件下底座和平台的接触应力分布及接触状态，指出顶梁扭转加载时的载荷不对称导致的底座不完全接触的现象.在不改变顶梁和底座约束边界的条件下，分析4种销轴连接简化方法对液压支架整体应力分布的影响.对比仿真结果与测试结果，分别评估各种边界条件和不同销轴连接简化方法的准确性，以及对仿真效率的影响.结果表明，液压支架的顶梁边界选择固定支撑，底座边界选择与工作台摩擦接触，保留销轴并将销轴与轴孔绑定接触得到的仿真结果最接近试验结果，且有较高的计算效率.研究结果可为液压支架设计过程中准确预估应力分布提供依据和参考.
关键词: 液压支架 强度特性 有限元 约束边界 销轴简化
Effects of constraint boundaries and simplified pin joints on strength of hydraulic support
HAN Chenyi1, NI Kaixuan1, MA Hui1,2
, WANG Hai3, LIN Junzhe1 1. School of Mechanical Engineering and Automation, Northeastern University, Shenyang 110819, China;
2. Key Laboratory of Vibration and Control of Aerodynamic Equipment (Northeastern University), Ministry of Education, Shenyang 110819, China;
3. Sany Heavy Equipment Co., Ltd., Shenyang 110027, China
Abstract: To accurately analyze the strength characteristics of hydraulic support considering the canopy torsional loading condition, by taking ZY6800/08/18D hydraulic support as the research object, the three-dimensional solid finite element model is established. According to the results of hydraulic support test, selecting the location of measuring points on the simulation model, the stress distributions of the hydraulic support are compared under eight boundary conditions for the canopy and base. In addition, the contact stress and status are analyzed under the frictional contact boundary condition between the base and worktable. It is pointed out that the phenomenon of incomplete contact of the base is caused by the load asymmetry during the canopy torsional loading. The influence of four simplified methods for pin joint on the stress distribution is also analyzed under the constant constraint boundary condition. Finally, the stress distribution obtained from the experiment evaluates the accuracy and efficiency of various boundary and pin joint simplification methods. When the canopy boundary selects fixed support, the base boundary selects friction contact with the working table, and retaining the pin joint and bonding the pin joint to the axle hole, the simulation results are the closest to the test results, and the simulation has high-efficiency. The research results can provide the basis and reference for the accurate estimation of the hydraulic support stress distribution.
Keywords: hydraulic support strength characteristics finite element constraint boundary pin joint simplification
液压支架是煤矿开采过程中重要的安全支护设备之一，其安全性和可靠性对于煤矿安全开采至关重要.为保证液压支架结构安全，在设计阶段均需进行结构强度校核.由于液压支架结构及载荷复杂，很难通过材料力学方法开展强度校核；而采用“数值近似”和“离散化”思想的有限元方法，不仅精度高，而且能适应各种各样复杂形状，因而成为液压支架强度计算行之有效的分析方法.
很多学者采用ANSYS Workbench和Abaqus等有限元软件，对各种型号的液压支架进行强度特性仿真分析[1-9]. Liu等[1]建立了ZT6500/19.5/34型液压支架三维有限元模型，分析了4种工况下支架的应力和变形分布.基于Abaqus软件，Li等[2]将各销轴连接简化为耦合铰接，垫块与顶梁和底座之间的接触简化为无摩擦接触，立柱和千斤顶之间设置为刚性连接，建立了ZF5000/16/28型液压支架有限元模型，分析了顶梁偏载、底座扭转工况下支架的应力和位移分布.基于Abaqus软件，李灿等[3]建立了某矿用液压支架有限元模型，模型中各零件之间的焊接采用节点共用法，通过带摩擦的接触模拟顶板和底板边界约束；分析了顶梁偏心加载、顶梁扭转加载和顶梁两端加载工况下液压支架的应力分布，并通过试验验证了仿真结果的有效性.基于双线性本构模型，Zhao等[4]建立了液压支架静力学有限元模型，分析了支架静强度及箱体焊接结构的疲劳特性. Wu等[5]建立了六柱液压支架的三维有限元模型，并进行了静力学仿真，分析了底座两端集中载荷和底座扭转载荷工况下各个部件的应力分布. Kong等[6]采用带摩擦的接触模拟销轴的转动，用有限元方法分析了3种危险工况下液压支架的强度特性.高耀东等[7]采用Shell单元和Solid单元建立了液压支架混合单元有限元模型，单元间连接采用绑定接触(honded contact)，采用摩擦接触处理部件间的铰接，将顶板和顶梁、底板和底座之间的接触处理为摩擦接触，分析了顶梁偏载和底座两端集中载荷、顶梁扭转和底座两端集中载荷工况下的液压支架的应力分布. Hu等[8]利用HyperWorks有限元软件对液压支架进行网格划分，并使用ANSYS计算了液压支架在顶梁偏心加载和底座扭转加载两种工况下的可靠性和灵敏度. Meng等[9]建立了基于ADAMS软件的液压支架动力学数值仿真模型，观察了不同荷载强度和不同位置下液压支架的响应，并对影响支架承载能力的因素进行了对比分析，为提高其稳定性和强度提供了理论依据.
为了验证有限元模型和仿真结果的有效性，很多学者也开展了试验验证工作[10-16].基于ANSYS Workbench软件，朱钧麟[10]对液压支架进行了静强度有限元仿真，对比了焊接采用绑定接触(bonded contact)和合成为一个零件处理(add)，顶梁加载时底座采用绑定和带摩擦接触，局部网格加密(子模型方法)对支架应力分布的影响.在文献[10]研究的基础上，洪岸柳[11]对某液压支架开展了有限元仿真和试验测试工作，对比结果表明，垫块和顶梁的接触方式为摩擦接触时的计算值与测试值最为接近.基于Abaqus软件，秦润泽[12]对ZF10500/20/36型放顶煤液压支架进行有限元分析和现场应力测试，其仿真模型考虑了销钉的转动关系，底座加载时对垫块的局部点的自由度进行约束处理.马园园等[13]针对ZF5000/16/28型放顶煤液压支架，基于有限元方法, 进行了顶梁两端加载和顶梁扭转加载工况下支架的强度仿真，其仿真模型在焊接接触面处选择绑定接触，在销轴连接处采用带摩擦的接触. Ma等[14]建立了三维有限元模型，分析了顶梁两端加载、顶梁扭转加载和底座扭转加载3种载荷工况下的支架应力分布，并通过试验验证了仿真结果的有效性.李祥松[15]对某型国内液压支架的应力进行了有限元仿真和试验测试，并与国外对应型号的支架进行了相对系数对比.
通过上述文献分析可以发现，很多研究在边界约束条件选取、不同部件之间联结选择，如销轴连接处的处理方法还存在较大的不同.而不同的边界约束和部件连接处理都会对计算精度和计算效率产生较大影响，目前的研究都是基于某种特定边界进行分析，还没有全面对比研究不同边界条件对仿真结果的影响；而且在与实验结果对比时，文献中所用到的测点也较少，无法准确评估边界条件对液压支架整体应力分布的影响.
本文以顶梁扭转加载工况为例，分析了不同边界条件和销轴不同处理方法对仿真结果精度和计算效率的影响，并通过对比试验测试结果，给出了合理的边界条件和销轴处连接的等效方法.
1 液压支架有限元模型 1.1 液压支架CAD模型简化ZY6800/08/18D型液压支架如图 1(a)所示，由图 1可知，实际液压支架结构非常复杂，存在很多辅助装置如护板、护帮等.在进行有限元仿真之前，需要对液压支架进行模型简化处理.文献[16]指出，根据液压支架的结构和原理，应忽略支撑时对力学影响较小的零部件，从而将液压支架转换为较为简单的四连杆机构.在不影响液压支架强度分析结果的前提下，对支架上的护板、吊环等不重要部分以及工艺孔、凸台等结构特征进行简化是符合工程实际的.
Fig. 1
1—顶梁；2—掩护梁；3—连杆；4—立柱；5—底座；6—垫块；7—工作台 图 1 ZY6800/08/18D型液压支架及CAD模型简化图 Fig. 1 ZY6800/08/18D hydraulic support and simplified CAD model 本文模型的具体简化原则如下：
1) 删除非承力件，如顶梁部装中的左右两侧的护板、护帮千斤顶、液压元件等，仅保留顶梁、掩护梁、前后连杆、底座这些关键承力部件[11]. 2)删除承力部件中的小零件，如顶梁部中的吊装液压软管的吊环，可以安放插销的凸台，放护帮千斤顶的套筒，小的液压油口等. 3)承力件的特征处理，将模型棱边的工艺倒角和圆角删去，将模型留有的焊缝及细微的间隙填平，并删去工艺孔[3]. 4)将立柱从模型中删去，用等效的载荷代替[7].简化后的顶梁扭转加载时的液压支架CAD模型，如图 1(b)所示.
1.2 液压支架载荷计算及有限元模型建立由文献[7]可知，液压支架的立柱在整机仿真分析中可视为二力杆，对顶梁和底座的柱窝上施加了一对大小相等、方向相反的力.立柱的液压缸额定工作压力p为42.3 MPa，液压缸活塞面积S为8.04×10-2 m2，参考GB25974.1—2010规定，单立柱分析载荷取1.2倍的工作阻力，由
$F = 1.2p \times S$ (1)
计算可得分析载荷F为4 080 kN.加载示意图如图 2(a)所示.平衡油缸对应力影响不大，本文忽略了平衡油缸的影响[11].
Fig. 2
图 2 液压支架边界、载荷及有限元模型 Fig. 2 Boundary and load of hydraulic support and finite element model 液压支架的材料牌号为HG70，弹性模量为210 GPa，泊松比为0.3.基于ANSYS Workbench建立液压支架有限元模型.为了提高分网的效率，顶梁、掩护梁、连杆和底座部装的各个零件之间的焊接采用绑定接触(bonded contact)进行处理，采用Solid186单元和Solid187单元混合的分网方法，将单元尺寸设为35 mm，其节点数为630 654，单元数为187 006(见图 2(b)).
2 仿真结果及实验验证本节主要对比在顶梁扭转加载工况下，液压支架有限元模型的约束边界和连接销轴的简化方式对静强度仿真结果的影响，通过对比试验测试结果，评估其计算效率和计算精度，给出最合适的边界设置和销轴简化方法.本节中所讨论的仿真结果均是在Intel i7-6700 CPU 3.4 GHz、16 G RAM的电脑上计算完成.
根据真实样机试验布置的测点位置，确定仿真模型的测点位置，如图 3所示.部分位置的测点由于测试得到的应力值过小，不具有参考价值，被舍弃，所以图 3中的标号并不连续.液压支架加载时的照片如图 4所示.使用LMS数据采集仪(VB8E、DB8板卡，128通道)，采用中航电测应变片(直角应变花，120 Ω)测试标定测点处的应力，通过软件合成提取Von Mises等效应力；试验共采集3组数据，数据分析时采用稳定状态下的Von Mises应力平均值，作为测试的准确结果.在试验过程中，在应力变化梯度较大的测试位置，通过增加额外的辅助测点来保证测试结果的准确性.如图 4所示，5#位置布置了5′和5″两个测点，来减少测试数据的偶然性，其余工况类似.
Fig. 3
图 3 试验和仿真提取测点位置图 Fig. 3 Positions of measuring points for simulation and experiment Fig. 4
图 4 液压支架试验现场及部分测点图 Fig. 4 Hydraulic support test site and part of measuring points 2.1 约束边界简化处理影响分析为了便于对比，销轴与轴孔之间统一采用绑定接触(bonded contact)处理，8种边界条件的工况如表 1所示.
表 1
表 1 8种边界条件的工况 Tab. 1 Eight boundary conditions 工况 顶梁约束边界 底座约束边界 1 固定支撑(fixed support) 固定支撑
2 固定支撑 无摩擦支撑(frictionless support)
3 固定支撑 不分离接触(no separation contact)
4 无摩擦支撑(frictionless support) 摩擦接触(frictional contact)，摩擦因数取0.2
5 固定支撑 摩擦接触，摩擦因数取0.2
6 固定支撑 无摩擦接触(frictionless contact)
7 固定支撑 粗糙接触(rough contact)
8 固定支撑 只压缩支撑(compression only support)加约束Z方向的位移(见图 5)
表 1 8种边界条件的工况 Tab. 1 Eight boundary conditions
Fig. 5
图 5 各工况下液压支架应力云图 Fig. 5 Stress nephograms of hydraulic support under different conditions 8种工况下液压支架应力分布云图如图 5所示；各个测点与试验结果对比如图 6和表 2所示.
Fig. 6
图 6 8种不同约束边界下的仿真及测试应力值对比 Fig. 6 Comparison of simulated stress and measured stress under eight constraint boundaries 表 2
表 2 试验测试值和8种工况下的仿真应力值 Tab. 2 Measured stress and simulated stress under eight boundary conditions ? MPa测点标号 不同工况的测点应力
试验 工况1 工况2 工况3 工况4 工况5 工况6 工况7 工况8
4′ 165.7 165.8 178.7 175.7 284.2 182.7 214.9 222.0 175.5
4″ 134.5 151.0 168.2 151.8 261.2 145.4 185.4 203.7 156.6
5′ 271.3 269.0 248.5 295.3 268.8 304.6 288.2 295.3 275.8
5″ 303.1 273.0 269.0 303.8 273.0 306.6 292.3 303.8 278.8
6′ 234.3 226.3 243.4 212.9 221.3 241.4 218.2 221.3 219.9
6″ 177.7 178.5 183.8 164.4 188.9 188.9 170.4 156.8 171.4
6″′ 159.2 163.1 175.7 149.7 168.6 149.3 153.1 144.4 155.3
8′ 329.5 328.6 263.7 341.8 317.0 342.1 321.6 356.6 314.2
8″ 303.3 305.4 274.6 302.7 305.8 302.5 313.4 317.9 308.0
8″′ 228.4 252.2 288.4 244.7 258.0 245.2 249.2 228.9 248.3
23′ 335.3 277.3 298.3 296.2 281.4 308.9 300.0 317.3 283.7
23″ 340.9 305.3 286.7 292.8 273.4 325.3 301.8 335.2 290.6
23″′ 353.3 301.3 274.1 281.2 261.5 326.4 282.2 336.5 300.5
10 102.4 14.3 92.2 90.8 142.6 100.3 91.5 11.1 90.1
14′ 114.8 5.1 11.4 14.5 194.4 117.9 179.7 15.5 171.4
14″ 114.4 4.5 12.0 16.0 189.7 113.7 177.3 15.6 167.0
15′ 165.3 3.3 12.1 15.2 226.9 166.2 258.4 19.5 233.5
15″ 165.6 3.6 13.8 16.4 258.0 153.7 277.0 19.0 236.0
18′ 105.4 3.7 9.9 9.6 160.0 104.7 161.8 9.4 147.3
18″ 105.8 3.8 9.6 8.7 154.9 99.8 166.7 9.0 151.4
19′ 106.2 1.2 6.3 6.8 169.0 116.2 175.9 9.5 158.4
19″ 107.1 1.3 6.6 5.4 177.7 117.5 178.7 11.5 172.6
20′ 142.2 3.4 14.9 12.4 267.8 157.7 227.3 22.8 210.1
20″ 146.2 3.2 15.3 14.2 246.2 165.3 238.2 23.9 214.9
表 2 试验测试值和8种工况下的仿真应力值 Tab. 2 Measured stress and simulated stress under eight boundary conditions ?
工况1的约束边界为顶梁处的垫块下表面与顶梁绑定接触，垫块上表面选择固定支撑，底座的底面选择固定支撑.仿真所得的应力云图表明顶梁上的应力较大，而其余位置应力很小，如图 5(a)所示.相比其余工况, 底座应力值最小，除柱窝部分之外，底座的应力仅有2~20 MPa，与测试得到的100~165 MPa差距甚远，这是由于底座约束过刚所致，与实际试验时的受力状况不符.
工况2的顶梁约束边界和工况1相同，底座的底面选择无摩擦支撑，该支撑只约束底面竖直方向的位移，允许其余各方向的平移，仿真结果相对于工况1有一定改善(见图 5(b)).掩护梁和连杆处的应力明显增大，除柱窝外，底座的应力为2~50 MPa，但仍远小于试验结果，这表明放宽底座约束边界，会导致底座应力在一定程度上增大.工况3建立了工作台，工作台与底座之间设置为不分离接触，不允许接触面之间有竖直方向的分离，但是允许沿着接触面的无摩擦滑动，仿真结果与工况2的仿真结果接近，底座应力仍与试验差距甚远.
工况4和工况5均在底座和工作台之间设置摩擦接触，不同点在于工况4顶梁的垫块上表面采用无摩擦支撑，两者的仿真结果显示底座边界改为带摩擦的接触均会导致底座应力大幅增加，而且工况5与试验结果吻合良好.工况4由于顶梁约束不足，导致底座测点处的应力大于试验结果，如图 5(d)和5(e)所示.工况4和工况5的底座底面位移和接触状态如图 7所示，从图 7可以明显看到顶梁扭转加载时，底座与工作台之间存在不完全接触现象，工况4的底座底面Y向最大偏移位移大于工况5，说明工况4的底座尾部翘起距离更大.
Fig. 7
图 7 8种工况的底座底面的位移和接触状态 Fig. 7 Displacement and contact status of bottom surface of base under eight boundary conditions 工况6把底座边界设为与工作台无摩擦的接触，分析仿真结果可知，该工况下底座的应力明显大于试验结果，顶梁处的应力与试验值接近，如图 6所示.底座底面的Y方向的最大偏移位移和X方向的滑移位移均大于工况5，如图 7所示，这说明底座的约束不足，导致底座的仿真结果偏大.工况7把底座边界设为与工作台粗糙接触，分析发现该工况下的仿真结果与工况1较为接近.底座处的应力均很小，与试验结果相距甚远，这表明该工况对底座的约束过刚，与实际的工作状态不符.工况8把底座边界设置为只压缩支撑并约束了底面Z方向的自由度，仿真结果与工况6类似，底座的测点应力较实验结果偏大，说明该工况对底座的约束不足，与试验时的受力状态不符.
液压支架测点处的仿真应力值与试验测得的应力值的误差不超过15%，即可视为该测点的应力与试验结果吻合.例如工况1的24个测点中只有12个测点与试验应力对比后的误差不超过15%，即吻合的测点数占比为50%.统计8种工况下的仿真计算时间和吻合的测点数占比，如表 3所示，由表 3也可以看出, 使用工况5的边界条件得到的结果具有最高的准确度，且计算效率也较高.
表 3
表 3 8种边界条件下的仿真效率和精度对比 Tab. 3 Comparison of simulation efficiency and accuracy under eight boundary conditions 工况 仿真计算时间/min 吻合的测点数占比/% 1 9 50.00
2 5 37.50
3 10 54.17
4 495 33.33
5 55 100.00
6 171 45.33
7 116 45.33
8 37 50.00
表 3 8种边界条件下的仿真效率和精度对比 Tab. 3 Comparison of simulation efficiency and accuracy under eight boundary conditions
2.2 连接销轴简化处理影响分析一些文献中，同样对约束边界进行了讨论，但是在对销轴的处理方式上有所不同.本节在工况5的约束边界下，讨论不同的销轴连接简化方法对仿真结果的影响，采用的销轴处理方式如表 4所示. 4种不同的销轴简化方法得到的液压支架应力分布云图如图 8所示；各个测点与试验结果的对比如图 9和表 5所示.计算效率和精度如表 6所示.
表 4
表 4 销轴连接简化方法 Tab. 4 Simplified methods of pin joint 工况 销轴处理方式 5 保留连接销轴，销轴和轴孔使用绑定接触(bonded contact)
9 去除连接销轴，连接位置直接采用绑定接触(bonded contact)
10 去除连接销轴，连接位置用考虑变形的旋转副(revolute joint)代替
11 保留销轴，销轴与轴孔之间使用摩擦接触(frictional contact)
表 4 销轴连接简化方法 Tab. 4 Simplified methods of pin joint
Fig. 8
图 8 4种不同销轴简化方法的应力云图 Fig. 8 Stress nephograms of four pin joint simplification methods Fig. 9
图 9 不同销轴简化方法下的仿真应力值与试验结果对比 Fig. 9 Comparison of simulated stress and measured stress under different pin joint simplification methods 表 5
表 5 试验测试值和4种工况下的仿真应力值 Tab. 5 Measured stress and simulated stress under four conditions ? MPa测点标号 不同工况的测点应力值
试验 工况5 工况9 工况10 工况11
4′ 165.7 182.7 176.2 270.7 206.1
4″ 134.5 145.4 143.1 194.3 159.4
5′ 271.3 304.6 286.2 302.9 305.4
5″ 303.1 306.6 286.9 304.1 313.7
6′ 234.3 241.4 219.3 239.1 236.3
6″ 177.7 188.9 176.9 191.9 201.4
6″′ 159.2 149.3 164.9 172.8 162.7
8′ 329.5 342.1 303.5 362.5 329.1
8″ 303.3 302.5 281.8 294.2 304.4
8″′ 228.4 245.2 243.8 241.5 252.9
23′ 335.3 308.9 289.4 288.4 320.2
23″ 340.9 325.3 304.5 295.6 327.4
23″′ 353.3 326.4 284.4 303.1 324.1
10 102.4 100.3 5.2 98.4 102.6
14′ 114.8 117.9 124.6 115.7 140.3
14″ 114.4 113.7 123.6 113.5 138.5
15′ 165.3 166.2 175.5 159.4 180.8
15″ 165.6 153.7 177.3 156.3 186.1
18′ 105.4 104.7 108.6 119.5 117.1
18″ 105.8 99.8 112.2 120.9 119.6
19′ 106.2 116.2 122.7 114.8 117.2
19″ 107.1 117.5 130.6 118.2 119.4
20′ 142.2 157.7 170.2 192.9 162.4
20″ 146.2 165.3 175.1 195.8 167.4
表 5 试验测试值和4种工况下的仿真应力值 Tab. 5 Measured stress and simulated stress under four conditions ?
表 6
表 6 不同销轴简化方法计算效率和精度对比 Tab. 6 Comparison of calculation efficiency and accuracy of different pin joint simplification methods 工况 仿真计算时间/min 吻合的测点数占比/% 5 55 100
9 118 70.83
10 382 83.33
11 818 83.33
表 6 不同销轴简化方法计算效率和精度对比 Tab. 6 Comparison of calculation efficiency and accuracy of different pin joint simplification methods
工况9的仿真应力云图表明, 在该工况下液压支架的连杆表面以及掩护梁处的应力值很小.由图 9可知，测点10位置的应力只有5.24 MPa，远小于试验结果.说明去除销轴后，将连杆的部分外表面与掩护梁和底座直接采用绑定接触，会抑制连接部分的变形，与实际工作状态不符，而且计算时间相比工况5增加了一倍左右，效率下降.
工况10同样去除连接销轴，在连接位置采用考虑变形的旋转副(revolute joint)代替，将仿真应力与试验数据对比，发现该工况下只有少部分测点的应力值与试验测得的应力相差较大，连杆与掩护梁与底座连接位置的应力与实验结果接近，说明该工况的销轴处理方式与实际的工作状态比较接近.由表 4可知工况10计算效率不及工况5.
工况11在工况5的基础上，保留销轴，把销轴与连杆轴孔之间的绑定接触改为摩擦接触.从图 9和表 5中可以看到，只有少部分测点应力值较实验值偏大，由表 6可知，工况11和工况10有相同的吻合测点数.但是，将销轴与轴孔的绑定接触改为摩擦接触之后，会导致计算效率大大降低.
3 结论本文建立了ZY6800/08/18D型支架有限元模型，针对顶梁扭转工况，对比分析不同边界约束条件和销轴连接简化方法对液压支架仿真应力分布的影响规律，同时比较了计算效率；并与试验结果对比，给出了不同约束边界的适用条件，主要结论如下：
1) 对于顶梁扭转加载工况，如果仅为了保证顶梁处的仿真应力较为准确，从而对液压支架顶梁应力做初步评估，当顶梁的垫块上表面采用固定支撑时，底座边界采用固定约束、与工作台不分离接触以及工作台粗糙接触时，均具有一定的精度，而且这些边界均是线性的，计算效率很高.
2) 在顶梁的边界采用固定支撑的情况下，底座下表面采用固定支撑，会使底座底面Y向和X向的位移为0，约束过刚，导致底座的应力偏小，与实际底座工作状态不符；底座下表面分别采用无摩擦支撑、与工作台不分离接触、与工作台粗糙接触，仿真的结果稍有改善，但得到的底座应力仍远小于试验结果；底座采用与工作台无摩擦接触或采用只压缩支撑并约束底面Z方向的位移均会导致底座的仿真应力值偏大.而底座选用与工作台摩擦接触, 仿真得到液压支架整体应力分布结果与试验吻合最好.还需指出摩擦接触存在非线性，很大程度上会降低仿真的计算效率.
3) 考虑销轴简化方法对仿真结果的影响时，去除销轴在连接位置直接使用绑定接触，会导致掩护梁处的应力偏小，与实际不符.使用考虑变形的旋转副代替销轴时，仿真结果较为接近试验结果，但是计算时间大大增加.保留销轴，将销轴和孔的绑定接触改为摩擦接触时，计算精度也较高，但添加的多个摩擦接触会大大降低计算效率.综合考虑计算效率和计算精度，建议保留销轴，并将销轴与轴孔绑定接触.
4) 不改变边界条件的情况下，局部结构的修改会导致仿真应力值的变化，若不改变顶梁上表面和底座下表面等重要位置的结构，仅对液压支架的部分零件做出修改，采用工况5的边界条件，应仍具有较好的准确性.
参考文献
[1] LIU Hongyu, LIU Junqing.Hydraulic support entire strength research under four kinds of typical load working conditions[C]//Proceedings of 2011 Second International Conference on Mechanic Automation and Control Engineering. Hohhot: IEEE, 2011: 329. DOI: 10.1109/MACE.2011.5986925
[2] LI Linghui, ZHENG Xiaowen. Three-dimensional modeling and finite element analysis of the hydraulic support[J]. Advanced Materials Research, 2013, 712/713/714/715: 1002. DOI:10.4028/www.scientific.net/AMR.712-715.1002
[3] 李灿, 苟晓明, 雷刚, 等. 矿用液压支架强度分析[J]. 重庆理工大学学报(自然科学版), 2014, 28(5): 44.
LI Can, GOU Xiaoming, LEI Gang, et al. Strength analysis of mine hydraulic support[J]. Journal of Chongqing University of Technology (Natural Science), 2014, 28(5): 44. DOI:10.3969/j.issn.1674-8425(z).2014.05.009
[4] ZHAO Xiaohui, LI Fuyong, LIU Yu, et al. Fatigue behavior of a box-type welded structure of hydraulic support used in coal mine[J]. Materials, 2015, 8(10): 6609. DOI:10.3390/ma8105325
[5] WU Jiangyu, FENG Meimei, TANG Kun. Strength analysis of six-pillar hydraulic support base[J]. Mechanics and Architectural Design, 2017, 61. DOI:10.1142/9789813149021_0008
[6] KONG Shuai, KONG Lingzan, WANG Lulu. The strength adaptabi-lity analysis of super-large mining height hydraulic support[J]. International Core Journal of Engineering, 2017, 3(11): 229.
[7] 高耀东, 周同. 用有限元法模拟液压支架强度试验[J]. 机械设计, 2018, 35(12): 66.
GAO Yaodong, ZHOU Tong. Simulation of strength test of hydraulic support by means of finite element method[J]. Journal of Machine Design, 2018, 35(12): 66. DOI:10.13841/j.cnki.jxsj.2018.12.012
[8] HU Denggao, FAN Xun. Analysis of the strength and reliability of hydraulic support with finite element method[J]. Advanced Materials Research, 2013, 619: 225. DOI:10.4028/www.scientific.net/AMR.619.225
[9] MENG Zhaosheng, ZENG Qingliang, GAO Kuidong, et al. Failure analysis of super-large mining height powered support[J]. Engineering Failure Analysis, 2018, 378. DOI:10.1016/j.engfailanal.2018.04.011
[10] 朱钧麟.液压支架应力研究及有限元分析[D].沈阳: 东北大学, 2011
ZHU Junlin. Stress research and finite element analysis of hydraulic support[D]. Shenyang: Northeastern University, 2011
[11] 洪岸柳.液压支架的结构强度与疲劳寿命分析[D].沈阳: 东北大学, 2012
HONG Anliu. Structural strength and fatigue life analysis of hydraulic support[D]. Shenyang: Northeastern University, 2012
[12] 秦润泽.放顶煤液压支架的应力测试及有限元分析[D].西安: 西安科技大学, 2012
QIN Runze. The stress testing and finite element analysis to the caving hydraulic support[D]. Xi'an: Xi'an University of Technology, 2012
[13] 马园园, 谢里阳, 洪岸柳, 等. 液压支架应力测试研究及有限元分析[J]. 机械设计与制造, 2012(11): 34.
MA Yuanyuan, XIE Liyang, HONG Anliu, et al. Finite element analysis and hydraulic support stress test[J]. Machinery Design & Manufacture, 2012(11): 34. DOI:10.3969/j.issn.1001-3997.2012.11.012
[14] MA Yuanyuan, XIE Liyang, QIN Xiaofeng. Strength and reliability analysis of hydraulic support[J]. Advanced Materials Research, 2012, 544: 19. DOI:10.4028/www.scientific.net/AMR.544.18
[15] 李祥松. 某型号液压支架的应力测试对比分析[J]. 煤炭技术, 2015, 34(5): 244.
LI Xiangsong. Comparative analysis of stress tests of certain type hydraulic support[J]. Coal Technology, 2015, 34(5): 244. DOI:10.13301/j.cnki.ct.2015.05.094
[16] WANG Peiqin, XU Zengshun, SUN Zuofeng, et al. Design and simulation of fully mechanized hydraulic support[J]. Applied Mechanics and Materials, 2014, 644: 199. DOI:10.4028/www.scientific.net/AMM.644-650.199
