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

1、6 / 13螺栓联接的一种改进刚度模型Sayed A.Nassar,A ntoine Abboud摘要：建立了由相同和不同联接板组成的螺栓联接的改进刚度模型。一种新方法用来获得决定联接刚度有效接触面积的表达式。提供一个更加精确的螺栓刚度的计算和更可靠的联接行 为的预测，不管在初始装配中还是在后来实际应用中施加拉伸载荷。研究了联接板厚度与螺栓公称直径之比、联接板尺寸、螺栓头接触半径比、 孔间隙和板材料/厚度比对刚度的影响。 根据实验数据来决定分析模型的锥角a。有限元模型用来评价刚度模型的精确性。1.引言精确测定螺栓联接刚度对于在外载荷作用下预测螺栓总成的行为是至关重要的；联接刚度也影响到在大批量

2、生产过程中如汽车工业，达到预期的预紧力水平所需的周期。在初始装配中，螺栓预紧力Fb等于联接压力Fj都等于同一预紧力水平Fi。当外载荷Fe施加在联接上时(如图1 (a)和1 ( b),螺栓拉力的增加和相应联接压力的减小是由联接和螺栓刚度 之比kc/kb决定的。Juvinail和Marshek1提供了图1 (a)和1 (b)联接模型的螺栓拉力Fb和压紧力Fc如下式所示：Fb 二 F&kbkc)Fe(1)k“-仁斤(2)式中：Fi是初始螺栓拉力；Fe是外载荷；kb是螺栓刚度；kc为联接刚度。Motosh2提出了一种计算螺栓联接刚度的技术如图2所示；他使用锥角和球形包络线来显示高应力部分的联接

3、。一个关于r的四次多项式来表达应力在Z方向上的分布，如图 3和4所示。计算相应的联接弹性压缩量，联接刚度kc可以用公式3来表示：,Fckc-(3)0式中：是联接压缩量(联接厚度的总减少量)°Motosh2认为联接挠度 与螺栓孔边缘(r=dh/2) 的垂直线一致。这相当于在Z方向上使用一个均匀的最大压力，这会导致总的联接挠度的值偏大，螺栓联接刚度偏小。另外，这个方法对于联接刚度没有固定的解决形式；为了获得联接刚度的数值解，他十分依赖每一个几何构型的数值积分。Lehnhoff等3,4提供了一个数值模型和有限元分析( FEA )来计算联接刚度。在他们的 研究工作中，他们假设在 Z方向上有一

4、个均匀的压应力，呈锥角分布，如图3所示。他们推导的联接刚度kc如下式所示兀d tanE,.kC 二(4)2ln (Ltana +Yd -d id +d )n (Ltana "d +d X" -d )式中：d为螺栓公称直径，E为联接材料的弹性模量，L为联接长度，为决定有效联接面积所需的角度，为螺栓头接触半径比。这个假设将导致挠度值偏小和联接刚度偏大。另外，在他们的研究中，锥角不是一个固定的常数。他们不得不随着联接参数来变化角度，以此来匹配他们的 FEA的结果。Shigley等5对于Lehnoff等3,4的方程进行了修改，取和接触半径比为，0.5774 dE(5)获得了联接刚度

5、kc的简单表达式:2ln 50.5774L 0.5dJ 0.5774L + 2.5d 丿Wileman等和Norton7提出了一个基于有限元分析的联接刚度计算公式。联接挠度 为垫圈中心线最上面节点的位移。在他们的研究中主要考虑的变量有联接长度L和螺栓直径d。比值d/L从0.1到2.0，改变螺栓直径同时保持联接长度不变（25.4mm）来计算联接刚度。Rasmussen等8根据有限元分析得出的计算有效面积Ac的公式（6）。但是，当L/d>5时，公式不适用。A*= -( 1-dh2)+1(D(O21)tan(0.35“ J1+ 2L*2 -12( D。2 _)(6)D。=D°/( d

6、),dh "/( d),L 二 D / ( d)A =代/( d)2式中：Ac为压紧联接的有效面积，dh为螺孔直径，D0为联接板尺寸。冋* 1 m lypicsi nwdiel of e baited jcrfrtl ettid - ma<U Df a|0lmFlg_ 3； Hair me 血I nf boll&dl joint Gor L tar 皿*公式（6）后来被Alkatan等9修改；他们指出只有当摩擦系数为0时，它给出的结果才和FEA的结果一致。Alkatan等9在公式（7）中摩擦系数取 0.2，他们指出公式（7）对于薄板联接更加精确。Ac = 4（ 1-dh

7、2）+2（D°2-1）ta门屮驚 誉）42（D° -dh ）Juvinail和Marshek1用一个公式计算联接的有效面积Ac,如下式所示:兰fd+d4 A"2(7)(8)式中：d1 d， d2=1.5d,d 3=1.5d+Ltan30 ° 。联接的应力区域是一个等于30°的锥角，如图2和3所示。很明显，Juvinail的方法不考虑联接板尺寸 Do的影响。以前大多数的文献不能提供一个精确的公式来说明相关参数如比值L/d和接触半径比等的影响。本文提出一个测定轴对称且联接板中心和孔同轴螺栓联接模型刚度的分析程序。变量包括联接厚度（夹紧长度L），螺栓

8、直径d,联接板尺寸Do,接触半径比，孔隙，不等厚板和不同材料板。一些特殊的几何例子也被考虑，如偏心螺栓孔和方板联接。2.提出联接刚度模型提出了一个更加精确的计算螺栓联接刚度kc的分析模型。分析模型通过FEA和实验数据来证实。等效面积 A （图5）通过下式来获得：,AcE八kcLkc= cAcc（ 9）LE锥角通过关联分析结果和实验结果得到。不同变量的影响通过 FEA软件（ABAQUS来研究。Fig. 5 Equivalent elastic model or the |ointFiy. 6 model cf bolted |oinl ForIan刚度方程的建立分两步走。第一，压应力p（r,z）

9、通过压紧力Fc来表述。第二，获得加紧部分的相应总体挠度 ，代入公式（3）来得到联接刚度kc。图1展示了一个螺栓联接的轴 对称模型，上下板之间关于接触面水平对称，所以允许使用1/2模型（图3和6）。假设一个包络线角 的锥角应力分布模型；一个 r-z坐标系。从z轴到锥角边界的距离 x表示为：xzta n：（ 10）22.1压应力分布在模型厚度范围内的一个任意z值，假设压应力分布是关于r的4次多项式（r,z）,如图4所示：-(r, z)二 Ar4 B r3 C r2 Dr E( 11)式中：AI B C D 和E ，是z的函数，这些参数由边界条件和静力平衡所决定。假 设空隙很小；因此，现在它可以被忽

10、略。空隙的影响将在5.3讨论。对于每一个z值，假设应力在r=d/2和r=x的位置梯度为0,。另外，假设应力在r=x的位置 士 / j2 = 0以此来达到在锥角边界压紧区域时都是固定的零值。2 2在 r=x，设；二一:二 / ；：= -：；/ ；：r = 0产生以下三式:Ax4 Bx3 Cx2 Dx E 0（ 12）4Ax3 3Bx2 2Cx D JO（ 13）12Ax2 6Bx 2C =0（ 14）在r=d/2，设：:二/ ；：r =0产生下式：d 3d 2 d4A（ ）33B（ ）2 2C （ ）D =0（ 15）22 2静力平衡考虑在每一个 z值，联接压力Fc和结果压应力分布（）之间的关系

11、如下式所示：Fc 二."A匚（2 二 rdr） =2 二 爲（Ar5 Br4 Cr3 Dr2 Er）dr（ 16）求解公式（12）-（ 16）产生参数 A（z）、B （z）、C（z）、D（z）、E （z）如下式所示:(17a)(17b)(17c )(17d)(17e)获得，(15M)(|)6FcA (z) = 22-2-2-2_(空)6 +2(空)5 +5(空)4 _20(空)3 + 25(空)2 _14(经)+3dd dddd*44x d ,B(z) 一 1)(JA3 d 2C(zH2(2X)(2 2)(d)2Ad d 2D (z) 一2x2dA1 /2x、3/2x ，、/ d、4

12、E -七）（）（4）（ ） Ad d 2公式（17a）-（17e）代入公式（11）来获得一个指定的z值和夹紧力R（从距离x和锥角 如图3所示）的压应力分布 （r,z ）。2.2弹性联接压缩下一步是判断总的联接压缩量相当于公式（11）给的压应力分布。在弹性区域内，正dz所对应应变通过胡克定律获得。从FEA的结果来看，在联接接触面（在z=0处），r和 方向的应力 差不多是1/3和1/2的正应力，在z>0处，它可以被忽略。所以，高度的变化量的近似正应变平均值d average如下式所示:averaged zdz(d(XT？drd/2 E(18)平均正应变可能按照圆盘厚度除以初始厚度dz来表达(

13、19)2.3相同的板材料在这节，联接刚度kc和相应的有效面积 A所对应的模型是两压紧的板都是相同材料的不管每块板的厚度。利用对称z=L/2,公式（18）产生的总联接压缩量 如下式所示:12L/2 :=20X :一drdzd/2 E产生的联接变形量为:Fcd( 3)d 二 E tan :-Ltan：Y了d + Ltana 一d 丨 I 3d d( -1)(20)条件为 Do 一（Ltan：d）(21)必须声明的是方程（20）适用于图3所示的模型，联接板尺寸D0足够大，允许应力包络线在联接板厚度上完全展开（D0_（Ltan：d）。图6所示的螺栓联接模型应力包络线只是部分展开, 孔直径和螺栓头下的接

14、触面积。 它满足下面的条件：因为联接板尺寸 D0与联接板厚度L， 对于这个模型，联接板尺寸D）不符合公式（21）的条件;d : D0 : (Ltan:d)(22)根据相似的方法，在联接厚度丄ln2FcL和压紧力Fc下相应的联接压缩量 如下式所示：卩 Y+3d *4（LtanD° + Yd） Id二 D°+3d 八 丫1 丿（D°+3d）（D°-d）条件为 d :： D0 ： （Ltan："Cd）结合公式（9）、（ 20）、（ 23）给出的联接刚度 kc和有效压紧面积 Ac如下所示：(23)当 Do _ (Ltan:#d)时kc = Fc2了 +

15、 3丫？d + Ltann'd+Ltana+3d 丿(ot-厂丿(24)d 兀 tanL32詁仝| "d Ltan： 3d(25)-1当 d : D0 : (Ltan：d)时k =台2亠卩0 _d t 4(Yd .d'''D0+3dhY-1 丿L tan ： - D°)(Do 3d)(Do-d)(26)8 / 13(27)Ac =r( Y 十3、iDo d + 4(Yd +Ltana D0) 1lDo+3d 丿(丫-1 丿(D。+3d)(D°d) j(31)10 / 13Fig. 7 Two pldtee with dl'&

16、amp;glmllr materlai fortan a为了考虑空隙的影响，孔直径dh代替方程（25）和（27）中的螺栓直径d,这样有:当 D0 _ (Ltan：d)时Acdh 二 tan 一丄Yd+3dhYdLta n 3dhL ta n ： - dd -dh(28)当 d :： D0 ： (Ltan：d)时-2丄ln&d +3dh、D。-dh "/(Yd + Ltana - D0)-dhlD0 +3dh j< ；'d -dh /(Do+3dh)(D°-dh) 一二 L ta n 二(29)式中：孔直径dh等于螺栓直径d加上空隙2c，如图1（a）所示

17、。这里必须指出的是当联接板尺寸Db小于等于螺栓头外径时（D。兰Yd ），螺栓头下面的压应力分布可以认为是常数，有效面积近似等于螺母的联接面积。2.4不同联接材料当两个压紧板是由不同材料组成的，联接刚度可以表示为：kc =Fc(30)式中：、“和、：2分别为上下联接板各自的挠度。每块板的挠度都是通过公式（18）在板厚度上积分获得的，根据联接板尺寸D。，每块板的挠度可以表示为：当 D0 （Ltan：d）时，仆“m d 二 Em tan :当 d : D。:：: (Ltan： d)时,kc2Fc二 Em tan :"iD0 d I +4(2tm tana 一 D0 十了d) I丿 1 丿f

18、 (Do+3d)(D°-d) jm=1,2分别表示上下板（图7）,公式（30） - （ 32）产生联接刚度kc如下式所示:(EiE2)lnEdn： tana£ EM；鶯;d> E2ln(；f 豐;d(33)当 d : D0 : (Ltan：；©)时,(32)kE1E: tana( 34)C (Ei+E2)nf Y+3 *Dod $ + 4Ei(2t2tanaD。+Yd)+4E2(2ti tanaD°+Yd)d匕2 Do 3d -1(Do 3d)(D°d)式中：巳和E2分别表示上下板的弹性模量。假如板厚度是相等的（t, =t2 =t），联接

19、刚度表达式（33）和（34）可以表示为:当 Do _ (Ltan： d)时EE2d二 tan：k =(巳 + E2) |ln II 料-1)人Yd+2ttana+3d)jj(35)当 d :： Do : (Ltan：；''d)时,巳 E2 兀 tana甘印討缶阳+4器z(36)2.实验数据一个MTS测试仪用来测试联接压紧挠度；载荷挠度曲线的斜率就是联接刚度。图8为不同联接尺寸和联接长度的试验样本。三种不同联接尺寸D0，分别为25.4mm,40.5mm, 63.5mm联接长度分别为 30mm和60mm铝材料的联接板被采用，因为当施加相同 的压力载荷时，在相同的几何情况下，铝将比钢

20、铁产生更大的挠度。因此，在试验测试过程中联接厚度减少的百分比误差将比同样的钢铁材料小得多。施加11K N的压紧力并测得相应的变形量。联接刚度通过公式（3）算得。图9为在MTS测试仪中的测试夹具。11 / 13(37)为了确保测试数据的可重复性，每一个试验都要至少重复 10次。从公式（25）和（27）36。锥角是最相关的。分析模的分析模型数据和根据联接载荷变形关系的试验数据中得出 型联接刚度预测与实验数据的误差百分比如下：%differe nee =Ac(experimental) - Ae(analytical)100 Ac(exp erime ntal)F»Bh 8 SeinpiQ

21、joints with oriferein* atze* «ndgrip 血gthH(n) AITS电| DeLdilctllF® Q Expwlnwfitdi vetup For Jlhi eornprv«6ic?ni WeiEKprmerilal data on effective area-即ntalr=±一)>岁*.I 叱WA-C-UJH12 / 13Table 1 Experimental versus analytical versus REM results of AcJeiMricinimr】dLtr卜mm I1 mm)l-LMAi

22、kalvlicaEx peri menialI2.7TO25.42K027928012.76025.429735112.73040,阳+4412.76040.,52012 1K)fil fl畑12.76063.55S05&1585（只有第五组属于锥角完全展开）提出的分析模型和实验结果的误差从0.5%到15.9%，如图10和表1所示。除了第一个联接板尺寸D0 = 25.4mm以外，误差随着联接长度的增大而减小。对于一个固定精度的实验测 试装置，增加被测件的厚度L将自然减小实验测试的误差百分比。4. 有限元分析图11为本文有限元证实模型刚度提出的1/2轴对称联接模型。在ABAQUS商用软件

23、进行有限元建模并在x-y平面内采用轴对称单元分网。对称边界条件通过联接厚度强加。摩擦接触面，摩擦系数为0.2，用在螺栓头和联接面之间的接触。等效加紧力通过在螺栓头施加均 布载荷的方式施加。分析了不同的几何配置模型。 为了确定联接刚度， 平均的节点变形， 相对于两板之间的 接触面，使用的是螺栓头下面的接触面。联接刚度是由一个夹紧力Fc产生的连接诶变形决定的。提出的分析模型和 FEA结果的误差从0.3%到14%如表2所示。Knilt III Am-ifiiTFig. 13 Effect ot me Joint cllAmeler an meamTable 2 FEA v&rauifi an

24、alytical re«iLlti farL dJ1 jutti)Imnr-)"aiiaJyticiJ '(mm可Perce no diileeoice(ft)2.4LVXKb 5I01.H342-1i132.4M9.312 8245aB8.a149,37.64.7i.yX1173y.t4.7us血206.3114.75t223222V93.024I.V117肌6UJ2412.7W.4372924ic12.715M3门小d.O4719PJVS 72.4q a4.7117515.050.0to4.75.fj12.75UO.O5.200.32.4L9417.740731

25、5u522.6597JIO24J*35活阿5. 结果的讨论提出的联接刚度的分析模型已经被试验和FEA的结果所证实：所有的结果都是一致的。FEA模型被用来研究联接板尺寸Do,链接长度与螺栓直径之比L/d，接触半径比，空隙和不等厚和不同材料的联接板对刚度的影响。图12-22为FEA模型的有效面积Ac或者联接刚度kc的结果，并与分析模型的结果进行比较。提出的分析模型与 FEA的结果误差从0.3%到 14%FEA结果和试验结果的误差从0.39%到9.9%。又一次，对于第一个联接板尺寸 Do二25.4mm，联接长度越小误差越高。5.1联接板尺寸D0和联接长度L的影响：相同板材料图12-14为联接板尺寸

26、D0和联接长度L在三种螺孔直径，8mm 12.7mm 16mm寸于有效面积Ac的影响。对于一个圆柱体联接，固定的螺孔直径8mm联接长度与孔直径的比值L/d=2.4，增加联接板尺寸从 D0 = 15.2mm （应力包络线部分展开）到D 40mm （应力包络线完全展开），导致有效面积 Ac增加31%。当比值L/d=4.7时，相应的有效面积 Ac個为D0的增加）增加54%这是因为随着联接厚度的增加，比率L/d增加，从应力包络线部分展开（z0）_ （D0 - d）/ （2tan ：）到应力包络线完全展开（zo）： （D0 - d） /（2tan ：）。弭Joint INdnirtr IMmni*(hi

27、lt IH jiiUFili r 41 i-T mHHill biuii'lrr il IfmiiHi18 / 13Fig, 13 Effect of the jolni dameter on the effective areaFig. 14 Effect of ths |oint d amtor on the sftMtlvc anaa 血5.2接触半径比的影响虽然标准的六角头螺栓的接触半径比取是流行的，但螺帽系数取其他值也是常见的。图15和16分别对应的是接触半径比在应力包络线完全展开和应力包络线部分展开如公式（25）和（27）对有效接触面积Ac的影响。 分别取1.5、1.8、2

28、.1。有效面积随着的增加而增加因为增加了应力包络线（图 3和4）。包络的面积越大，联接的压应力分布越 低。这个结果导致低的联接变形，相应的更高的联接刚度kc。FEA结果与分析模型结果有很好的相关性，最大误差为7%TLJr DirMiN'En -117iianiSU9llul l l)Li nricl if di I IL? min11i A1.4M.113onhkM RiiiJii 期tknyfc- FEIFig. IS Effect of the ccrnKt radi ratio an 4 lor 広 tan a+yd)5.3孔隙的影响Rg. W EIIkL 口I llw cont

29、act radii raUo Qin Ac for vdk。.对于相同的板材料，公式（28）和（29）用来研究孔隙对联接刚度的影响，通过有效面Fig 17 Erreci ar tne hole特arana 口n AP.二=vupv l«l.«r-lk- 14411L. =J J.干 nwiTi26BW12Hk(. ICMI'MJKfi%)Fla it EiiKt or mo hale k»afonc& ot %积Ac的方式呈现。使用了 4个不同的孔隙（dh-d ）如图1 （ a）,螺栓公称直径d的0% 3% 6呀口 10%。孔隙从0%增加到10%有

30、效面积 Ac平均减少9%如图17和18所示。在应力包 络线完全展开（D0 X（Ltana +?d），对于螺栓直径分别从 8mm 12.7mm 16mm， FEA结果和分析模型结果的误差Ac在6%和10.7%之间。图18展示了小尺寸联接板在应力包络线部 分展开（d : D0 : （Lta），相应误差在0.5%和10%间。5.4不同材料联接的厚度比的影响用公式（33）-（ 36）来研究铝板（Al ）和钢板（St）在不同厚度比对刚度的影响。符号t1 /t2、Al / St、St/Al分别指上下板的厚度比和材料联接。下板厚度保持常数20mm。图19和20展示了在应力包络线完全展开 （D0 _ （Lta

31、nd）对于螺栓直径分别从 8mm 12mm 16mm，FEA结果和分析模型结果的误差在0%和 9%之间。图21和22展示了小尺寸联接板在应力包络线部分展开（d : D0 ：（Ltan.G："Fd），相应误差在0.5%和11%之间。Fig、19 E flKl of lFhs thlckirwifi Bllo on Ih* Jolnilortnn 卄ychFig. 20 EHect of me thickness raik> on tiw |dnLforO工(I ImMwtri iiil:(Si 11.)¥47 mf itr. jU|- Httiin录< -W fl

32、 -f FEA f w： il列aa 3MD 机Ii 1FMt kw Rjitifii (NiiiniiiFig- 2i Eflwl of ttnc ftieknas ratio on> me |slnT stiffness for ycfwD卡口 Kan «r+苛Tll-s.-kl inJI：irb：il：i il/St iFig. 22 E If eelth<E tiilcknes-s raHici on Hie IdIiitLItjoir1皿8仇UH5.5试样与欧洲标准VDI 2230比较表3和4是本文的分析模型与 VDI2230结果10的比较。VDI提供以下有效面

33、积 Ac的 表达式:代二(d)2® 6(d)(D0-(d)(x 仁1)( 38)式中：x=(L d)/Do)1/3公式（38）仅适用于联接板尺寸在以下范围内d m Do 乞(Ltan：d)表3和4中可以看出本文提供的刚度模型比Lehnhof等3,4和VDI223010加接近FEA的结果。(39)的结果更A. FEAfirim7)A. LAnhflff(011/tin )labile 3 Comparison wltii yu 12230 rasuHs ter slrnllar material ptalesL乩工方.丄申札心已队7.t v+'anJttllTIrtTwllMiuff*I I I I 3 2 r3 ! 5 4 5 c,JJ 444444 -7.r7Lr =丄工.工2 J iL5L5l!ry24ZJ评”常JD&.O 1KS.2 N/AP4I