|
导读 作者：覃雄臻1 ，莫易敏2 ，覃佳亮1 ，郭胜会2 （1.上汽通用五菱汽车股份有限公司，广西 柳州 545007；2.武汉理工大学机电工程学院，湖北 武汉 430070） 来源：《机械设计与制造》2022年5月 摘要：针对螺纹联接结构在实际工况载荷下结合面相对滑动与螺栓（螺母）旋转松动失效预测精度较低的问题，通过研
究结构参数、工况载荷与松动失效之间的内在关系和力学原理，建立松动失效判定准则；在此基础上建立实际螺纹联接
结构的有限元模型，合理设置材料与接触属性，通过施加实际工况载荷进行有限元计算并提取被连接件结合面与螺纹
面的摩擦应力与剪切应力，以此判定螺纹联接结构是否发生松动失效；最后，通过道路试验验证了所提方法的准确性
与可靠性。 1 引言 螺纹联接结构在振动、冲击及变载荷的服役环境中，易发生 松动失效［1］。松动失效分为被连接件结合面的滑动与螺栓（螺 母）的旋转松动，两种松动都不允许在工程实际中发生。（1）螺纹 联接结构在一定预紧力的条件下因受外部载荷的作用，导致结合 面的夹紧力不足、摩擦力降低，引起结合面的相对滑动，从而打破 力（力矩）传递的稳定性，螺纹联接失效；（2）当螺纹面上原有的受 力平衡被外部载荷打破，引起螺纹面的相对滑动，就会造成螺栓 （螺母）相对转动，使预紧力下降，导致螺纹联接失效。 目前针对螺纹联接的松动失效问题，常用的解决方法是理论 力学推导、试验研究与有限元分析。VDI-2230［2］在设计螺纹联接 结构时，通过计算分界面抗滑动安全余量，保证结构的防松性能， 安全余量的限制值根据经验或由用户确定。文献［3］研究指出，螺 杆头下承压面发生整体滑移是螺栓发生松动的必要条件，并由此 给出了横向载荷作用下，临界松动载荷的计算公式。文献［3-6］通 过试验和有限元仿真研究发现，在横向载荷反复作用下，螺杆头 下承压面和螺纹面的局部滑移会在接触面逐渐累积弹性应变能， 当能量累积到一定程度，接触面将发生整体滑移，从而引起螺栓 松动。而在螺栓松动中，螺纹面先发生转动，带动螺杆头下承压 面转动［3-4］。文献［7］研究了随机振动载荷下的螺栓松动行为进行 了数值仿真分析。文献［8］研究了头部结构、弹簧垫圈、预紧力与 表面处理等因素对防松性能的影响规律。 综上，现有文献主要研究单一载荷作用下螺纹联接结构的松 动行为，而在工程实际中，螺纹联接结构所承受的载荷复杂（横向载荷、轴向载荷、弯矩、扭矩等），现有的研究成果无法准确预测实 际工况载荷下的松动失效。因此，在设计过程需要综合考虑结构 参数、实际工况载荷以及松动失效之间的内在关系。采用仿真的 方式可以通过对紧固件尺寸结构、接触状态、载荷工况、边界条件 进行模拟，实现最大程度地模拟实际情况，从而准确预测螺纹联接 的松动失效。鉴于此，通过研究松动失效机理建立松动条件的判 定准则，在此基础上运用ABAQUS仿真软件建立实际螺纹联接 结构的有限元模型，并施加由Adams软件提取的实际工况载荷进 行有限元分析，从而准确判定螺纹联接的滑动与旋转松动失效。 2 松动失效判定准则 针对螺纹联接结构结合面的相对滑动与螺栓（螺母）的旋转 松动问题，采用理论力学分析其松动条件，并建立相应的失效判 定准则。 2.1 被连接件结合面相对滑动判定准则 由VDI-2230［2］可知，被连接件结合面不发生滑动的条件为 结合面的夹紧载荷FKR大于传递横向载荷FQ与扭矩MY需要的夹 紧载荷，如式（1）所示，夹紧载荷FKR的计算，如式（2）所示：
式中：FQ—横向载荷，N；MY—绕螺栓轴线的扭矩，N·mm；qF和qM— 传递力FQ和扭矩MY的分界面数量；μT—分界面摩擦系数；ra—扭矩MY所作用的夹紧部件的摩擦半径，mm。
式中：FM—装配预紧力，N；Φ* en—同心夹紧与偏心加载的载荷系 数；FA—轴向载荷，N；FZ—预紧力损失量，N。 由此，可建立基于有限元分析的结合面不发生相对滑动判定 准则，夹紧载荷FKR替换为有限元模型中的摩擦应力，横向载荷 FQ与扭矩MY替换为剪切应力，即结合面上任一节点的摩擦应力 均大于剪切应力。 2.2 螺栓（螺母）旋转松动判定准则 抑制螺纹面的局部滑移，保证螺纹面不同位置均不满足产生 滑移的力学条件，可从原理上防止螺栓（螺母）发生相对转动，由 此建立的旋转松动判定准则为外部载荷沿螺纹面任一点的切向 合应力τ小于摩擦应力f，如式（3）所示［1］。
式中：σQ、τQ—外部载荷在螺纹面产生的法向与切向应力，N/mm2 ；σA、τA—外部载荷在螺纹面产生的法向与切向应力N/mm2 ；γ—横向载荷与轴向载荷在螺纹面上投影线的夹角。 由此，可建立基于有限元分析的螺栓（螺母）不发生旋转松动 的判定准则，即螺纹面上任一节点的摩擦应力均大于剪切应力。 3 判定松动失效的FEM方法 运用ABAQUS仿真软件建立实际螺纹联接结构的有限元模 型，并施加实际工况载荷进行有限元分析，通过建立的结合面相对 滑动与螺栓（螺母）旋转松动判定准则，实现松动失效的准确预测。 3.1 螺纹联接结构有限元模型 这里以某车型上横臂与后副车架螺纹联接结构为例，如图1 所示。结构包括上横臂、后副车架、螺栓和螺母。
3.1.1 网格模型 将上横臂与后副车架几何模型导入Hypermesh软件中划分 网格模型。其中，螺栓布种密度为1，螺母布种密度为0.75，上横 臂与后副车架支架布种密度为3，接触面网格适当加密，得到的网 格模型，如图2所示。
3.1.2 材料模型 各部件材料模型设置为弹塑性材料，螺栓和螺母为高强度钢 SCM435，被连接件上横臂与后副车架为汽车结构钢 SAPH440。两种材料应力-应变曲线由材料拉伸试验获得，通过计算转换为 真实应力-真实应变曲线［9］，如图3、图4所示。材料属性参数，如 表1所示。
3.1.3 接触模型 接触属性包括接触面之间的法向作用和切向作用。法向作 用属性定义的是接触压力与接触面穿透之间的关系。法向作用 属性一般选择系统默认的接触压力与间隙间的关系—“硬接触”， 但有时为得到更精确的接触面应力分布而选择“软”接触。切向 作用属性是定义接触面之间的摩擦属性，一般选择罚函数模型。 这里研究的螺纹有限元模型定义了如下接触：螺栓头部与 上被联接件（Int-1），螺母支撑面与下被联接件（Int-2），上被联 接件与下被联接件（Int-3）、螺栓螺母螺纹接触面（Int-4）。由于 被联接件有发生分离的风险，故对其选择有限滑移设置，其余 接触面设置选择小滑移。对于螺栓和螺母的螺纹接触面设置， 在切向对其设置罚函数摩擦算法，法向设置“软”接触，其余各接 触面之间法向作用均设置“硬”接触，并且各接触面均将摩擦系 数根据实际结构提取，分别为0.13（Int-1）、0.13（Int-2）、0.22（Int3）与0.12（Int-4）。 3.1.4 预紧力加载及工况载荷加载 采用截面法［10］对螺栓施加预紧力，在螺杆上创建截面并在截 面上施加法向作用的目标预紧力。为使计算容易收敛，在第一个 分析步中施加一个10N的螺栓预紧力，第二个分析步中施加目标 预紧力。 创建被连接件中心线上的基准点，如图5所示。令该点为参 考点RP-1，并将参考点与上横臂进行耦合，在参考点RP-1上施 加极限工况载荷，如表2所示。
3.2 有限元分析与讨论 通过有限元仿真，得到被连接件结合面与螺纹面上的摩擦应 力云图，如图6所示。剪切应力云图，如图7所示。通过建立结合 面与螺纹面路径，提取周向节点的剪切应力与摩擦应力，为使得 提取结果最大程度的贴近真实仿真结果，共布点41个，周向路径 的应力曲线图，如图8、图9所示。通过比较各节点剪切应力与摩 擦应力的大小，判断是否发生相对滑动，若相对滑动，则加大预紧力重新计算；反之，则减小预紧力，直至得到临界滑动状态，即所 有节点的剪切应力均小于或等于摩擦应力，且至少有一个节点 的剪切应力等于摩擦应力。被连接件结合面、螺纹面临界滑动 状态的应力曲线图，如图 10、图11所示。对应的预紧力分别为 88kN、85kN。
4 效果验证 对螺纹联接结构进行道路试验，被连接件结合面发生相对滑 动的判断依据是汽车底盘有异响，螺栓（螺母）发生旋转松动的判 断依据是标记线发生了相对转动。分别施加不同大小的预紧力， 判断是否发生松动，得到结合面发生与不发生相对滑动时对应的 预紧力分别为83kN、86kN，螺栓（螺母）发生与不发生旋转松动时 对应的预紧力分别为 80kN、83kN，与有限元仿真结果列于表 3 中。对比发现，道路试验结果均比有限元仿真结果略小，最大误 差分别为6.02%与6.25%，分析原因为有限元仿真中采用局部节 点滑动作为松动条件偏保守，但在工程应用中，保守的设计偏安 全，可降低事故的发生概率。
5 结论 针对实际工况载荷下预测螺纹联接结构结合面相对滑动与 螺栓（螺母）旋转松动失效的难点，通过建立有限元分析方法的 失效判定准则，结合仿真软件得到了较为准确的松动失效预测 方法。 （1）通过研究结构参数、工况载荷与松动失效之间的内在关 系和力学原理，建立以被连接件结合面与螺栓（螺母）螺纹面的局 部滑动作为有限元分析松动失效的判定准则，以发挥有限元分析 方法的优势，为准确预测松动失效提供技术条件； （2）建立实际螺纹联接结构的有限元模型，合理设置材料与 接触属性，通过施加实际工况载荷进行有限元计算并在计算结 果中提取被连接件结合面与螺纹面的摩擦应力与剪切应力，以 此判定螺纹联接结构是否发生松动失效，实现了松动失效的准 确预测。 （3）通过道路试验验证所提方法的准确性与可靠性。对 比发现，在判定结合面相对滑动与螺栓（螺母）旋转松动失效 时，道路试验结果均比有限元仿真结果略小，最大误差分别为 6.02%与6.25%，说明有限元仿真结果偏保守，这有利于设计方 案的可靠性与安全性，在工程实际中具有一定的参考依据和应 用价值。 本公众号由北京艾法斯特科技发展有限公司运营，欢迎大家技术交流，有感兴趣的朋友可以添加下面的微信号，进入我们的讨论群： 本公众号由北京艾法斯特科技发展有限公司运营，欢迎大家技术交流，请添加上面微信二维码，或致电联系我们： 杜：18601022409 任：13466666822 共同探讨螺栓预紧力测量的相关技术和市场信息，共同推进此产业的发展，实现共赢。北京艾法斯特公司（www.ifast-sensor.com）专业从事超声、视觉螺栓预紧力测量技术，产品主要有螺栓预紧力测量仪系列产品、螺栓预紧力长期监测产品、智能紧固件、以螺栓轴力控制的智能紧固工具。
参考文献： ［1］张明远，鲁连涛，唐明明，等.横向载荷作用下螺栓临界松动载荷数值 计算方法研究［J］.机械工程学报，2018，54（5）：173-178. （Zhang Ming-yuan，Lu Lian-tao，Tang Ming-ming，et al.Research on numerical calculation method of critical load for bolt loosening under transverse loading［J］.Journal of Mechanical Engineering，2018，54（5）：173-178.） ［2］VDI-2230，Systematic calculation of high duty bolted joints：joints with one cylindrical bolt，Verein Deutscher Ingenieure，2014. ［3］PAI N G，HESS D P.Experimental study of loosening of threaded fasten⁃ ers due to dynamic shear loads［J］. Journal of Sound and Vibration， 2002，253（3）：585-602. ［4］王崴，徐浩，马跃，等.振动工况下螺栓连接自松弛机理研究［J］. 振动 与冲击，2014，33（22）：198-202. （Wang Wei，Xu Hao，Ma Yue，et al.Self-loosening mechanism of bolted joints under vibration［J］.Journal of Vibration and Shock，2014，33（22）：198-202.） ［5］PAI N G，HESS D P.Three dimensional finite element analysis of thread⁃ ed fastener loosening due to dynamic shear load［J］.Engineering Failure Analysis，2002（9）：383-402. ［6］Jing Xiang-jun，Zhu Yong-sheng，Hong Jun，et al.Investigation into the loosening mechanism of bolt in curvic coupling subjected to transverse loading［J］.Engineering Failure Analysis，2013，32（3）：360-373. ［7］朱光宇，郑保仲.随机振动螺栓松动数值仿真分析［J］.机械设计与制 造，2018（1）：233-235. （Zhu Guang-yu，Zheng Bao-zhong. Numerical simulation analysis of random vibration on bolt relaxing［J］.Mechanical Design and Manufac⁃ turing，2018（1）：233-235.） ［8］莫易敏，梁绍哲，晏熙，等.汽车高强度螺栓的防松性能的影响因素研 究［J］.机械设计与制造，2015（9）：89-92. （Mo Yi-min，Liang Shao-zhe，Yan Xi，et al.A study of effects on antiloosening performance of automobile high-strength threaded fasteners ［J］.Mechanical Design and Manufacturing，2015（9）：89-92.） ［9］邓小艳.轿车轮毂螺栓联接预紧力分析及其提升技术研究［D］.秦皇 岛：燕山大学，2018. （Deng Xiao-yan. Preload analysis and improving methods of bolt con⁃ nection in automotive aluminum alloys wheel［D］. Qinhuangdao：Yan⁃ shan University，2018.） ［10］刘建华.轴向激励下螺栓连接结构的松动机理研究［D］.成都：西南交 通大学，2016. （Liu Jian-hua. Research on the self-loosening mechanism of bolted joints under axial excitation［D］.Chengdu：Southwest Jiaotong Univer⁃ sity，2016.） |
|
|
