螺栓联接的一个改进刚度模型

螺栓联接的一种改进刚度模型SayedA.Nassar,AntoineAbboud摘要：建立了由相同和不同联接板组成的螺栓联接的改进刚度模型。一种新方法用来获得决定联接刚度有效接触面积的表达式。提供一个更加精确的螺栓刚度的计算和更可靠的联接行为的预测，不管在初始装配中还是在后来实际应用中施加拉伸载荷。研究了联接板厚度与螺栓公称直径之比、联接板尺寸、螺栓头接触半径比、孔间隙和板材料/厚度比对刚度的影响。根据实验数据来决定分析模型的锥角α。有限元模型用来评价刚度模型的精确性。1.引言精确测定螺栓联接刚度对于在外载荷作用下预测螺栓总成的行为是至关重要的；联接刚度也影响到在大批量生产过程中如汽车工业，达到预期的预紧力水平所需的周期。在初始装配中，螺栓预紧力Fb等于联接压力Fc；都等于同一预紧力水平Fi。当外载荷Fe施加在联接上时（如图1（a）和1（b）），螺栓拉力的增加和相应联接压力的减小是由联接和螺栓刚度之比kc/kb决定的。Juvinall和Marshek[1]提供了图1（a）和1（b）联接模型的螺栓拉力Fb和压紧力Fc如下式所示：()bbiebckFFFkk（1）()cciebckFFFkk（2）式中：Fi是初始螺栓拉力；Fe是外载荷；kb是螺栓刚度；kc为联接刚度。Motosh[2]提出了一种计算螺栓联接刚度的技术如图2所示；他使用锥角和球形包络线来显示高应力部分的联接。一个关于r的四次多项式来表达应力在Z方向上的分布，如图3和4所示。计算相应的联接弹性压缩量，联接刚度kc可以用公式3来表示：ccFk（3）式中：δ是联接压缩量（联接厚度的总减少量）。Motosh[2]认为联接挠度δ与螺栓孔边缘（r=dh/2）的垂直线一致。这相当于在Z方向上使用一个均匀的最大压力，这会导致总的联接挠度δ的值偏大，螺栓联接刚度偏小。另外，这个方法对于联接刚度没有固定的解决形式；为了获得联接刚度的数值解，他十分依赖每一个几何构型的数值积分。Lehnhoff等[3,4]提供了一个数值模型和有限元分析（FEA）来计算联接刚度。在他们的研究工作中，他们假设在Z方向上有一个均匀的压应力，呈锥角分布，如图3所示。他们推导的联接刚度kc如下式所示tantan2lntancdEkLddddLdddd（4）式中：d为螺栓公称直径，E为联接材料的弹性模量，L为联接长度，∝为决定有效联接面积所需的角度，γ为螺栓头接触半径比。这个假设将导致挠度值偏小和联接刚度偏大。另外，在他们的研究中，锥角∝不是一个固定的常数。他们不得不随着联接参数来变化角度∝，以此来匹配他们的FEA的结果。Shigley等[5]对于Lehnoff等[3,4]的方程进行了修改，取∝=30°和接触半径比为γ=1.5，获得了联接刚度kc的简单表达式：0.57740.57740.52ln50.57742.5cdEkLdLd（5）Wileman等[6]和Norton[7]提出了一个基于有限元分析的联接刚度计算公式。联接挠度为垫圈中心线最上面节点的位移。在他们的研究中主要考虑的变量有联接长度L和螺栓直径d。比值d/L从0.1到2.0，改变螺栓直径同时保持联接长度不变（25.4mm）来计算联接刚度。Rasmussen等[8]根据有限元分析得出的计算有效面积Ac的公式（6）。但是，当L/d5时，公式不适用。22212202000.351211=-d(1)tan()422()/(),/(),/(),/()chohhhccLLADDdDDddddLDdAAd（1）+（6）式中：Ac为压紧联接的有效面积，dh为螺孔直径，D0为联接板尺寸。公式（6）后来被Alkatan等[9]修改；他们指出只有当摩擦系数为0时，它给出的结果才和FEA的结果一致。Alkatan等[9]在公式（7）中摩擦系数取0.2，他们指出公式（7）对于薄板联接更加精确。22122010.75(0.2)=-d(1)tan()42()chohLADDd（1）+（7）Juvinall和Marshek[1]用一个公式计算联接的有效面积Ac，如下式所示：22321=42cddAd（8）式中：d1≈d，d2=1.5d,d3=1.5d+Ltan30°。联接的应力区域是一个等于30°的锥角，如图2和3所示。很明显，Juvinall的方法不考虑联接板尺寸D0的影响。以前大多数的文献不能提供一个精确的公式来说明相关参数如比值L/d和接触半径比γ等的影响。本文提出一个测定轴对称且联接板中心和孔同轴螺栓联接模型刚度的分析程序。变量包括联接厚度（夹紧长度L），螺栓直径d，联接板尺寸D0，接触半径比γ，孔隙，不等厚板和不同材料板。一些特殊的几何例子也被考虑，如偏心螺栓孔和方板联接。2.提出联接刚度模型提出了一个更加精确的计算螺栓联接刚度kc的分析模型。分析模型通过FEA和实验数据来证实。等效面积Ac（图5）通过下式来获得：=ccccAEkLkALE（9）锥角∝通过关联分析结果和实验结果得到。不同变量的影响通过FEA软件(ABAQUS)来研究。刚度方程的建立分两步走。第一，压应力p(r,z)通过压紧力Fc来表述。第二，获得加紧部分的相应总体挠度δ，代入公式（3）来得到联接刚度kc。图1展示了一个螺栓联接的轴对称模型，上下板之间关于接触面水平对称，所以允许使用1/2模型（图3和6）。假设一个包络线角∝的锥角应力分布模型；一个r-z坐标系。从z轴到锥角边界的距离x表示为：tan2dxz（10）2.1压应力分布在模型厚度范围内的一个任意z值，假设压应力分布是关于r的4次多项式σ（r,z）,如图4所示：432(,)rzArBrCrDrE（11）式中：A、B、C、D和E，是z的函数，这些参数由边界条件和静力平衡所决定。假设空隙很小；因此，现在它可以被忽略。空隙的影响将在5.3讨论。对于每一个z值，假设应力在r=d/2和r=x的位置梯度为0，∂σ∂𝑟=0。另外，假设应力在r=x的位置22/0r以此来达到在锥角边界压紧区域时都是固定的零值。在r=x，设22//0rr产生以下三式：4320AxBxCxDxE（12）324320AxBxCxD（13）212620AxBxC（14）在r=d/2，设/0r产生下式：324()3()2()0222dddABCD（15）静力平衡考虑在每一个z值，联接压力Fc和结果压应力分布σ（𝑟）之间的关系如下式所示：5432/2(2)2()xcdFdArdrArBrCrDrErdr（16）求解公式（12）-（16）产生参数()Az、()Bz、()Cz、()Dz、()Ez如下式所示：6654322(15/)()()222222()2()5()20()25()14()3cFdAzxxxxxxdddddd（17a）44()()(1)()32xdBzAd（17b）222()2()(2)()2xxdCzAdd（17c）2()2DzxdA（17d）34122()()()(4)()32xxdEzAdd（17e）公式（17a）-（17e）代入公式（11）来获得一个指定的z值和夹紧力Fc（从距离x和锥角∝获得，如图3所示）的压应力分布σ（r,z）。2.2弹性联接压缩下一步是判断总的联接压缩量δ相当于公式（11）给的压应力分布。在弹性区域内，正应变通过胡克定律获得。从FEA的结果来看，在联接接触面（在z=0处），r和θ方向的应力差不多是1/3和1/2的正应力，在z0处，它可以被忽略。所以，高度的变化量dz所对应的近似正应变平均值dεaverage如下式所示：/21()2xaverageddzddrddzEx（18）平均正应变ε可能按照圆盘厚度除以初始厚度dz来表达（/dzdz）。2.3相同的板材料在这节，联接刚度kc和相应的有效面积Ac所对应的模型是两压紧的板都是相同材料的不管每块板的厚度。利用对称z=L/2,公式（18）产生的总联接压缩量δ如下式所示：/20/212()2LxddrdzdEx（19）产生的联接变形量为：(3)tan2lntantan3(1)CFddLddEdLdd（20）0(tan)DLd条件为（21）必须声明的是方程（20）适用于图3所示的模型，联接板尺寸D0足够大，允许应力包络线在联接板厚度上完全展开（0(tan)DLd）。图6所示的螺栓联接模型应力包络线只是部分展开，因为联接板尺寸D0与联接板厚度L，孔直径和螺栓头下的接触面积。对于这个模型，联接板尺寸D0不符合公式（21）的条件；它满足下面的条件：0(tan)dDLd（22）根据相似的方法，在联接厚度L和压紧力Fc下相应的联接压缩量δ如下式所示：0000002134(tan)lntan31(3)()(tan)CFDdLDdEdDdDdDddDLd条件为（23）结合公式（9）、（20）、（23）给出的联接刚度kc和有效压紧面积Ac如下所示：当0(tan)DLd时tan=3tan2lntan31ccFdEkdLddLd（24）tan3tan2lntan31cdLAdLddLd（25）当0(tan)dDLd时00000tan=134(tan)2ln31(3)()ccFEkDddLDdDdDdDd（26）00000tan134(tan)2ln31(3)()cLADddLDdDdDdDd（27）为了考虑空隙的影响，孔直径dh代替方程（25）和（27）中的螺栓直径d，这样有：当0(tan)DLd时tan3tan2lntan3hchhhhdLAdddLddLddd（28）当0(tan)dDLd时00000tan134(tan)2ln3(3)()chhhhhhhLAddDddLDdDdddDdDd（29）式中：孔直径dh等于螺栓直径d加上空隙2c，如图1（a）所示。这里必须指出的是当联接板尺寸D0小于等于螺栓头外径时（0Dd），螺栓头下面的压应力分布可以认为是常数，有效面积近似等于螺母的联接面积。2.4不同联接材料当两个压紧板是由不同材料组成的，联接刚度可以表示为：12=ccFk（30）式中：1和2分别为上下联接板各自的挠度。每块板的挠度都是通过公式（18）在板厚度上积分获得的，根据联接板尺寸D0，每块板的挠度可以表示为：当0(tan)DLd时，32tanlntan2tan3(1)CmmmmFdtddEdtd（31）当0(tan)dDLd时，000002134(2tan)lntan31(3)()CmmmFDdtDdEdDdDdDd