螺栓固结界面接触刚度的虚拟材料建模

(1.北京科技大学机械工程学院，北京，100083；2.重庆大学机械传动国家重点实验室，重庆，400044)

螺栓联接是工程结构中最常见的连接形式之一。通过螺栓将不同的工作部件连接在一起，形成整体的机器结构，广泛应用于航空航天装备、高精度制造装备、精密仪器等机械装备[1-4]。各种类型复杂多样的螺栓联接结构的共性特征是相互接触的构件表面具有粗糙形貌，同时界面在一定预紧力下工作。螺栓固结界面的接触刚度是描述螺栓结构界面特性最重要的参数之一。研究表明，螺栓固结结构界面的接触刚度占机械结构总刚度的60%～80%，其变化直接影响界面以及机械装备系统的动力学性能[5-8]，因此，螺栓固结界面的准确动力学建模和接触刚度的准确计算，对机械结构系统的性能分析与预测至关重要。目前，人们对螺栓固结结构的动力学建模和界面刚度特性进行了大量的研究。研究人员提出了弹簧-阻尼单元模型、接触单元模型和虚拟材料模型的建模方法，来表征和描述螺栓固结界面的接触特性[9-15]。其中，弹簧-阻尼单元模型采用相互独立的离散弹簧-阻尼单元来表征螺栓固结界面的接触特性与动力学行为；接触单元建模方法采用零厚度单元、平面单元或三维单元来等效模拟螺栓固结界面的接触特性。弹簧-阻尼单元模型和接触单元模型存在无法描述实际连续的接触界面之间的相互耦合作用关系、建模精度不高、难以准确表征实际螺栓界面的接触特征，以及界面参数难以准确设置等缺点[9-10]。虚拟材料建模方法通过在固结界面之间加入一层具有一定属性参数的虚拟材料来等效接触界面的接触特性，获得了广泛的关注。例如，田红亮等[11]假设螺栓结合部是各向同性虚拟材料，并基于粗糙表面形貌的分形接触理论推导了材料参数。李小彭等[12]以组合梁为研究对象，将结合面等效为广义间隙，利用结合面接触刚度的分形模型与材料应变能等效的方法，获得了结合面广义间隙的材料常数。孙清超等[13]考虑结合面的压力分布不均匀性，提出了一种基于梯度虚拟材料的栓接结合面连接参数表征模型。但是，上述基于虚拟材料的结合面建模方法均是基于表面粗糙形貌的分形模型描述和Majumdar-Bhushan分形接触模型。对Majumdar-Bhushan接触模型，微凸体的变形过程是先塑性变形，再弹塑性变形，最后是纯弹性变形[14]。这与经典接触理论和材料本构关系中接触变形是先弹性变形后塑性变形的变形过程相矛盾[15-16]。同时，上述模型中均假设虚拟材料层的厚度为某一固定值，而不随接触界面的粗糙度而变化，不能准确表征界面粗糙度变化引起的虚拟材料厚度变化及对界面接触刚度的影响。本文作者考虑粗糙界面的弹性变形特征，基于粗糙表面形貌的Greenwood-Williamson统计学模型描述，提出一种基于统计学模型的虚拟材料建模方法，推导虚拟材料的弹性模量、泊松比、密度和厚度等材料参数的解析表达式，建立虚拟材料的材料参数与粗糙界面属性参数之间的关系式，分析螺栓固结界面的表面粗糙形貌、基体材料属性对虚拟材料的材料参数和界面接触刚度的影响特性，并采用有限元仿真和实验测试方法，验证虚拟材料建模方法的准确性。由于实际机械结构的接触界面通常承受周期性载荷的加载-卸载作用，界面粗糙微凸体的塑性变形被初始的加载-卸载作用消除，稳定工况时，界面工作在弹性接触状态，例如，GONZALEZ-VALADEZ等[17]研究发现11次加载-卸载过程可以消除界面的塑性变形。因此，本文的粗糙界面法向接触刚度计算模型仅考虑粗糙界面的弹性变形和弹性接触刚度。

1 螺栓固结界面的虚拟材料模型

1.1 模型描述

图1 粗糙界面的接触模型示意图Fig.1 Schematic of rough surface contact model

具有粗糙表面形貌的2个构件通过螺栓连接实现固结，其粗糙界面接触模型示意图如图1(a)所示。图1(a)中，在螺栓预紧力P的作用下，表面均线距离为h的两粗糙表面相互接触；E1和E2分别为构件1和构件2的弹性模量；v1和v2分别为构件1和构件2的泊松比；G1和G2分别为构件1和构件2的剪切模量。图1(a)所示的双粗糙表面接触可采用图1(b)所示的粗糙表面-刚性平面接触模型进行等效。其中，等效弹性模量1/E=(1-v12)/E1+(1-v22)/E2，等效剪切模量1/G=(2-v1)/G1+(2-v2)/G2，等效表面粗糙形貌z(x)=z1(x)+z2(x)，粗糙面均线与刚性平面之间的距离为h。粗糙表面与刚性平面的接触特性可以采用虚拟材料层进行表征，如图1(c)所示，其中，h为虚拟材料层的厚度；E0为弹性模量；v0为泊松比；ρ0为密度。虚拟材料与接触基体和刚性平面刚性连接，实现位移与载荷的连续传递。

1.2 接触刚度

当粗糙接触界面采用虚拟材料层进行表征时，由法向载荷所做的虚功为[11]

式中：Pn为单位面积上的法向载荷；An为名义接触面积；kn为单位面积上的法向接触刚度。

弹性接触体的法向应变能可以表示为

法向应变能与虚功相等，故界面单位面积上的法向接触刚度可以表示为

目前我国发展形式错综复杂，国内外市场的逐步融合加大了传统产业转型升级的压力，各类不利因素和有利条件同时存在。利用“互联网+”与大数据来进行产业转型和升级，可更为便利的使用技术分析获得消费者喜好、特征，进而得到更为精准、丰富的消费者数据，如社交数据、网站浏览数据及地理追踪数据等，辅助传统产业从多个维度细分消费者，展开个性化营销等。事实上，“互联网+”大数据思维恰恰最大的优势与价值在于可以全局性的思维，针对性的细分消费者群体，并将媒体、消费者、营销等更为个性化、精准化的处理。在明确消费者和市场需求、动态走向等后，传统产业的转型与发展自然更为顺遂。

由图1可知，异步协作构成了一个设计方案的“作者链”，利益分配模式则决定了怎样在这个链上的作者之间分配利益。

2 虚拟材料参数的解析表达式

2.1 弹性模量

图1(b)所示的粗糙表面-刚性平面接触采用Greenwood-Williamson统计模型(GW模型)[18]进行描述。其中，表面粗糙体的形状为球体且各粗糙体的曲率半径均为β，粗糙微凸体的高度符合高斯分布，其高度的标准偏差为σs，微凸体的面密度(单位面积上的微凸体数量)为n，各粗糙体与刚性平面的接触符合Hertzian弹性接触理论且粗糙体之间无相互作用。

单个微凸体与刚性平面接触示意图如图2所示。

图2 单个微凸体与刚性平面接触示意图Fig.2 Geometry of single asperity contact

根据Hertzian弹性接触理论，其弹性接触力可以表示为[18]

对粗糙表面形貌的Greenwood-Williamson统计模型描述，微凸体的高度分布为标准高斯分布，微凸体高度分布在zn=[-3,3]的概率为99.73%[19]。当虚拟材料层的厚度h=6σs时，可以覆盖粗糙界面上所有的微凸体。同时，考虑到接触体之间的距离改变时，微凸体产生变形，虚拟材料的厚度也相应改变。因此，虚拟材料的厚度与接触体之间距离变化的关系可以表示为

z为单个微凸体的高度，定义为微凸体顶点与微凸体高度均线间的距离；dd为微凸体均线与粗糙面均线之间的距离，dd=1.15σs[19]。

接触面积可以表示为

考虑到单个微凸体的法向应变ε=w/β，则单个微凸体的法向接触压力可以表示为

单个微凸体的弹性模量可以表示为

采用Greenwood-Williamson统计模型(GW模型)计算粗糙界面上接触微凸体的数量，并将单个微凸体的弹性模量对微凸体数量进行积分，获得虚拟材料层的弹性模量。引入量纲一的参数：则虚拟材料层的法向接触载荷和弹性模量可以表示为

在学历方面，公共图书馆招聘信息中本科学历需求占比67．6%，硕士研究生及以上学历需求占比32．4%，这说明本科学历是进入公共图书馆的最低要求。

其中：ϕn(zn)是高度为标准高斯分布的概率密度函数。

太宰治：我的不幸，恰恰在于我缺乏拒绝的能力。我害怕一旦拒绝别人，便会在彼此心里留下永远无法愈合的裂痕。

2.2 泊松比

单个微凸体接触的切向应力应变关系可表示为[10]

式中：εr为切应变；qr为剪切强度，与材料的屈服强度Sy相关，对于钢材，qr=0.6Sy；f为摩擦因数；pe为法向压力。则单个微凸体的剪切模量可以表示为

虚拟材料层的泊松比可以表示为

2.3 密度

虚拟材料层的密度可以表示为

为了验证螺栓固结粗糙界面虚拟材料建模方法的准确性，进行螺栓固结结构的模态实验测试。同时，基于推导获得的虚拟材料参数，代入实验测试样品的表面粗糙形貌参数、材料参数等，建立了基于虚拟材料的螺栓固结结构有限元仿真模型，计算其固有频率与模态特性，并将有限元仿真计算结果与实验测试结果进行对比。

则虚拟材料层的密度可以表示为

鉴于金融开放与经济增长之间关系的不确定性，很多学者尝试探究其影响因素。相关影响因素主要有制度质量和宏观条件、资本流动，汇率波动等。Klein发现一国制度因素和金融机构质量影响资本账户开放与经济增长之间的互动关系[7]。Bussiere et al发现从中长期看，制度质量和资本流入构成是金融开放促进经济增长的重要因素[8]。荣晨、董瞾以资本积累率分析资本账户开放与经济增长的关系[9]。李丽玲、王曦将汇率波动的作用引入资本账户开放和经济增长的关系中，发现控制汇率的波动可以促进资本账户开放收益[10]。

2.4 厚度

式中：E为接触界面的等效弹性模量；w为弹性变形量，

式中：hn为量纲一的粗糙面均线高度。

到底给孩子看什么书？我像所有家长一样迷茫。给她买过成套的绘本，她扫了一遍，就置之不理：“这都是小毛毛看的。”假期的短暂无聊里，实在找不到童书，她居然从我书架上找了一本《吸血鬼女王》。我分析来分析去，估计她以为名叫“女王”的书和“XX公主”一样，都是童话故事。

3 模型有效性验证

式中：ρ1为接触体的密度；Ar为粗糙界面的实际接触面积，

3.1 实验测试

设计了螺栓固结界面实验装置，进行结构模态实验测试。实验装置如图3所示。实验采用自由悬挂锤击法测试，利用力锤敲击上表面产生冲击激励，在钢板下表面安装加速度传感器，测量界面的法向加速度响应。振动分析仪型号为LMS SCADAS Mobile 8CH。

图3 实验测试示意图Fig.3 Schematic diagram of experimental setup

螺栓固结界面由2块长×宽×厚为300 mm×50 mm×8 mm的45号钢板，通过14个M6螺栓联接形成。名义接触面积An=1.5×10-2m2。通过扭力扳手给每个螺栓施加5 N·m的力矩，当拧紧6个螺栓时，界面的平均压力为p=1.67 MPa。采用激光共聚焦显微镜LEXT对相互接触的2块钢板表面进行表面形貌测试，采样间隔为0.625 μm。钢板表面形貌如图4所示。

在不同的文化中，上述情况并不常见，大多数情况下，它们是不对应的，所以只有当我们弄清楚典故性成语的起源时，我们才能理解其真实含义。

表面形貌参数采用下式计算[22]：

颈源性头痛患者远比大多数学者想象的多的多。笔者门诊每天接触的头痛患者中,按笔者提的诊断标准诊断,1/3可诊断为颈源性头痛。以笔者提出的颈源性头痛的治疗标准治疗,大部分患者疗效显著,复发率低。鉴于颈源性头痛患者数目庞大,如按笔者提出的诊断、治疗标准进行,临床医生将更容易掌握该病的诊断、治疗。尤其是基层医生的诊断将变的更容易,治疗操作更简单,风险更低,患者更安全,将使更多的颈源性头痛患者免于漏诊,更多的患者将受益,更多的医疗资源将会节省,更易在基层推广。

式中：m0，m2和m4为粗糙表面的谱矩[22]，

经过治疗,AB两组的手术时间和出血情况对比不存在统计学差异性(P>0.05)。AB两组的止痛效率不存在统计学差异性(P>0.0167),ABC三组的止痛效率、肺炎发生率、住院时间、满意度等存在统计学差异性(P<0.05)。AB两组外观满意度、肺炎发生率、住院时间不存在统计学差异性(P>0.0167)。

由图7可知：对具有不同粗糙度的表面，虚拟材料的弹性模量随着法向载荷的增大而增大，而泊松比和虚拟材料的厚度均随着法向载荷的增大而递减；在相同载荷变化范围内，虚拟材料的弹性模量和厚度随法向载荷剧烈变化，而泊松比的变化较小，即随着接触界面法向载荷变化，引起粗糙界面接触特性发生变化，在等效的虚拟材料参数中，弹性模量和厚度更显著地表征接触变化特性。

将计算获得的表面形貌参数σs=1.07 μm，n=3.09×1011m-2，β=0.315 μm以及等效弹性模量E=115.4 GPa代入式(10)，(14)和(17)，获得虚拟材料参数：E0=255 MPa，v0=0.375，ρ0=7 700 kg/m3。

3.2 基于虚拟材料的有限元仿真分析

图4 实验样品表面形貌Fig.4 Surface topograph of test specimen

基于计算获得的螺栓固结界面虚拟材料参数，采用有限元仿真计算方法对螺栓联接界面系统进行模态分析，获得系统的固有频率与模态特性。基于虚拟材料的螺栓联接结构有限元分析模型如图5所示。钢板和虚拟材料层均采用SOLID185单元离散，虚拟材料层与上下钢板共节点连接。实验测试获得的前5阶固有频率和其中3阶模态与有限元计算获得的固有频率和模态的对比，分别如图6和表1所示。

由图6和表1可知：实验测得的固有频率与基于虚拟材料模型计算值相对误差小于6%，且各阶模态一致，说明本文的虚拟材料建模方法与材料参数的计算模型是有效的。

1.2.4 妊娠结局随访 打电话对所有孕妇妊娠结局追踪随访或查阅本院电子病历，包括孕期超声筛查结果、妊娠结局、新生儿体检的外貌、结构、智力发育等。

图5 基于虚拟材料的螺栓联接结构有限元分析模型Fig.5 Finite element analysis model for bolted jointed structure based on virtual material model

4 虚拟材料参数与接触刚度特性分析

改变界面的接触特性参数，包括粗糙表面的形貌参数和接触基体材料属性，获得不同界面特性参数时螺栓固结界面虚拟材料参数和界面的接触刚度随界面压力的变化关系曲线。3种具有不同表面粗糙度的界面分别为光滑界面、中等粗糙界面和粗糙界面，其表面形貌参数如表2所示。不同接触基体的材料参数如表3所示。对不同粗糙度和不同接触基体的界面，采用式(10)，(13)，(14)，(17)，(18)和(3)计算获得螺栓固结界面的虚拟材料参数和接触刚度随界面压力的变化关系曲线。

图6 不同阶数的实验测试模态与仿真计算模态对比Fig.6 Comparison of vibration mode shapes between simulation and experiment

表1 固有频率实验值与仿真值对比Table1 Comparison between natural frequencies from experimental test and FEA

表2 不同粗糙度表面形貌参数Table2 Parameters of surfaces with different roughness

表3 不同接触基体的材料参数Table3 Material and sound parameters of different contact bodies

4.1 虚拟材料参数变化特性

当表面粗糙度不同时，虚拟材料的弹性模量、泊松比和厚度随法向载荷的变化关系曲线如图7所示。接触基体的材料属性为钢。

一杭摇头，想借着锁把板扣拧开，试了试，没成功。老太太从床底下找到一把钳子，一杭接过来，把锁拧开了。里面是一个黑漆漆的小匣子，一杭见过，装零钞用的。一杭说：“这个你也留下吧，以后可以用到。”老太太打开看了看里面的钱，犹犹豫豫地说：“钱也归我？”一杭点头。老太太赶紧把钱塞进裤兜里，又用力按了按，积极地帮一杭清理抽屉里的东西，生怕一杭反悔似的。

其中：Nr为粗糙表面的数据点总数。

由图7还可知：在相同载荷作用下，随着界面粗糙度增大，虚拟材料的弹性模量、泊松比均减小，而虚拟材料的厚度递增。

山西黄河一线旅游资源的厚重程度与其市场竞争力不相匹配，黄河旅游产品总体知名度低，吸引东南沿海、南方和海外游客的能力低，客源主要来自省内、周边省份和环渤海地区[12]，众多游客只知壶口而不知山西沿黄其他景点。以碛口景区为例，近年来其游客量和旅游综合收入增长迅速（见表2），但通过笔者的实地调查发现，众多游客前来碛口游玩的主要动机是这里门票免费，大多数游客在游玩过后的感受平平，重游意愿不强，还未形成良好的旅游口碑效应。

当接触基体不同时，虚拟材料的弹性模量、泊松比和厚度随法向载荷的变化关系曲线如图8所示，界面为中等粗糙度。

由图8(a)可知：与钢基体相比，基体为铸铁和铝时，虚拟材料的弹性模量降低。这是因为钢材、铸铁和铝的弹性模量依次递减(见表3)。由图8(b)可知：当接触基体材料不同时，钢基体虚拟材料的泊松比最大，铝基体次之。这是因为铸铁的泊松比最小，其相应的虚拟材料的泊松比也最小；同时，随着基体材料的弹性模量减小，虚拟材料的厚度也递减：在相同载荷作用下，钢基体材料的虚拟材料厚度最大，铝基体材料的厚度最小。

C2组总体氨基酸分析，鸡胸肉、腿肉和肝脏中分别为0.75%、0.44%和9.89%，D2组鸡胸肉、腿肉和肝脏中分别为0.76%、0.35%和7.92%。表明在放养条件下至180 d，无抗养殖的腿肉和肝脏中总氨基酸含量明显更高。

4.2 接触刚度特性

不同表面粗糙度时界面接触刚度随法向载荷的变化关系曲线如图9所示，接触基体的材料属性为钢。

由图9可知：对具有不同粗糙度的表面，界面接触刚度均随着法向载荷的增加而递增，其中，光滑接触界面的刚度最大，粗糙界面的刚度最小，界面接触刚度随着表面粗糙度的增大而递减。

当接触基体不同时，虚界面接触刚度随法向载荷的变化关系曲线如图10所示，界面为中等粗糙度。

图7 不同表面粗糙度时虚拟材料的材料参数随法向载荷的变化关系Fig.7 Plots of virtual material parameters versus normal load for surfaces with different roughness

由图10可知：与钢基体相比，基体为铸铁和铝时，界面接触刚度依次递减，铸铁基体次之，铝基体最小；当法向载荷较小时(p＜50 MPa)，接触基体材料属性的差异对界面接触刚度的影响较小；随着法向载荷增大，接触基体材料属性的影响增大。

图8 不同接触基体时虚拟材料的材料参数随法向载荷的变化关系Fig.8 Plots of virtual material parameters versus normal load for different contacting bodies

图9 不同表面粗糙度时粗糙界面接触刚度随法向载荷的变化关系Fig.9 Plots of contact stiffness versus normal load for surfaces with different roughness

图10 不同接触基体时粗糙界面接触刚度随法向载荷的变化关系Fig.10 Plots of contact stiffness versus normal load for different contacting bodies

5 结论

1)提出了一种基于统计学模型的粗糙接触界面虚拟材料建模方法，推导了虚拟材料的弹性模量、泊松比、密度和厚度等材料参数的解析表达式，分析了螺栓固结界面的表面粗糙形貌、基体材料属性对虚拟材料的材料参数和界面接触刚度的影响特性。

2)通过螺栓联接结构的模态实验测试和基于虚拟材料的有限元仿真分析，验证了粗糙接触界面虚拟材料建模方法的正确性。

3)虚拟材料的材料参数与法向载荷、界面粗糙度和接触基体弹性模量相关：随法向载荷和基体弹性模量的增加，虚拟材料的弹性模量增大，泊松比和厚度减少；随界面粗糙度的增加，弹性模量和泊松比减少，而厚度增大。

4)螺栓固结界面接触刚度随法向载荷增大而递增，随界面粗糙增大而减小。