基于修正刚度方法的复合材料多钉连接钉载分配预测方法与流程

本发明涉及复合材料机械连接钉载分配的技术领域，具体涉及一种基于修正刚度方法的复合材料多钉连接钉载分配预测方法，适用于航空航天飞行器中广泛使用的复合材料多钉连接结构。

背景技术：

多钉连接是复合材料重要的连接形式，也是复合材料结构的强度薄弱环节。通常，多钉连接结构失效分析的第一步是预测螺栓载荷分布，该分布通常不均匀，受参数影响大。

目前常用的复合材料机械连接钉载分配确定方法主要有解析法、有限元法和刚度法。复合材料螺栓连接结构钉载分配确定的解析方法是基于lekhnitskii等的复合材料孔板孔边应力分析发展起来的。由于钉载分配确定比孔边应力确定复杂得多，偏微分方程的求解需要采用边界配置法进行求解。但是由于当自由度增加时，边界配置法计算量增加迅速，许希武等利用保角映射及faber级数避免了上述问题。

尽管解析法具有计算简单、容易程序实现的特点，但是，忽略紧固件形状和偏心弯矩等因素对钉载分布的影响。采用有限元分析研究钉载分配不仅可以考虑紧固件形状和偏心弯矩，还可以对复杂形状的连接结构进行分析。但是采用有限元方法获得的计算精度与所需要的建模和计算工作量是矛盾的两个方面，因此，如何在一定计算规模下尽可能提高有限元分析的精度是一个重要的研究内容。

刚度法是一种基于弹性力学的钉载分配的简化计算方法。它将连接板和螺栓均简化为沿载荷方向具有一定刚度的元件，利用外载荷作用下不同元件的变形协调关系来求解每个元件的载荷，从而得到钉载分配，具有简单、形象、求解容易的特点，在规则排列多螺栓连接的钉载计算中应用广泛。刚度方法的计算式，较多的是mccarthyandgray的矩阵形式，liu等还给出了一种螺栓载荷的单参数迭代计算公式及其计算策略，计算公式如下所示：

其中，fi^a和fi^b(i＝1,···,n)分别是a、b两板在第i和第(i+1)个螺栓之间的层板内力；δi^a和δi^b分别是a、b两板在第i和第(i+1)个螺栓之间的纵向(沿加载方向)变形；fi和δi分别代表第i个螺栓的载荷和纵向变形；和分别是a、b两板在第i和第(i+1)个螺栓之间的层板等效刚度。

上述公式(1)-(3)中的螺栓载荷fi是δ1的单调函数，因此采用的是一维搜索求解的方法。方法开始迭代求解时，首先假设钉1的孔变形δ1，然后采用二分法根据公式(1)-(3)可以计算所有的钉传载荷fi。如果计算的钉传载荷的和大于外载荷f则减小δ1，否则增大δ1，当钉传载荷的和与外载荷的差别足够小时，假设的δ1为真实的孔变形，计算的钉传载荷即为真实的钉传载荷。

螺栓刚度模型是一个复杂的函数，对钉载分配的影响很大。刘等人在mccarthy和gray的三段式公式的基础上，考虑间隙对螺栓轴承弦刚度的影响，改进了螺栓刚度模型，并通过单螺栓接头的试验确定了参数，该模型为：

其中，k¹是初始拟线性区域的连接刚度；f^c是临界摩擦载荷，表示过渡区高度；表示过渡区域高度的临界摩擦载荷；c是螺栓孔间隙，表示过渡区域长度；ac是螺栓孔间隙的影响系数；k³是螺栓载荷传递区域中的净配合接头的轴承弦杆刚度。

层板等效刚度和可以由公式(5)计算得到：

其中，el是层合板的纵向均匀杨氏模量；w是层合板的宽度；h^j是层合板a或b的厚度；li是螺栓节距。

技术实现要素：

本发明要解决的技术问题为：弥补现有刚度法的不足，提出了一种预测复合材料多钉连接钉载分配的方法，该方法考虑了旁路载荷引起的孔拉伸变形的影响，能够准确的预测和研究多螺栓连接结构的钉载分配问题。

本发明解决上述技术问题采用的技术方案为：一种基于修正刚度方法的复合材料多钉连接钉载分配预测方法，实现步骤如下：

步骤a，提出与孔拉伸变形相对应的附加刚度的概念；所述步骤a中提出与孔拉伸变形相对应的附加刚度的概念实现过程为：

(a1)引入旁路载荷引起的孔拉伸变形以及对应于孔拉伸变形的附加刚度i是对螺栓个数的描述，a、b分别指代螺栓连接的上下两个层合板；

(a2)得到的孔拉伸变形及其对应的附加刚度之间的关系为：

其中，fi^a和fi^b(i＝1,···,n)分别是a、b两板在第i和第(i+1)个螺栓之间的层板内力。

步骤b，提出新的刚度方法方程；所述步骤b中提出新的刚度方法方程实现过程为：

(b1)考虑孔拉伸变形存在后，第i及第i+1个螺栓之间的变形协调方程为：

其中，δi是第i个螺栓的纵向(沿加载方向)变形；δi^a和δi^b分别是a、b两板在第i和第(i+1)个螺栓之间的纵向变形；是b板第i个螺栓孔的拉伸变形；是a板第i+1个螺栓孔的拉伸变形。

(b2)代入板和孔拉伸的刚度，将板和孔拉伸变形用力表示，则第i及第i+1个螺栓之间的变形协调方程为变为：

其中，fi^a和fi^b(i＝1,···,n)分别是a、b两板在第i和第(i+1)个螺栓之间的层板内力；和分别是a、b两板在第i和第(i+1)个螺栓之间的层板等效刚度，是对应于孔拉伸变形的附加刚度。

(b3)根据(b2)中得到的变形协调方程，得到第n个螺栓的变形协调方程为：

(b4)将(b2)和(b3)中得到的两个考虑孔拉伸变形的螺栓之间的变形协调方程，替换现有刚度法计算公式中的变形协调条件，得到考虑孔拉伸变形的新的刚度方法方程为：

步骤c，提出附加刚度的计算方法。

所述步骤c提出附加刚度的计算方法实现过程为：

为了得到附加刚度的计算方法，研究一块同时存在钉传载荷和旁路载荷的孔板，对两个孔之间的部分发生的变形进行分析，得到附加刚度的计算公式为：

其中，kh为所研究孔板的刚度，没有理论公式，采用有限元求解得到；为第i和第(i+1)个螺栓之间的层板等效刚度，采用现有刚度法中对层板刚度的计算公式求解。

步骤d，采用单变量迭代的二分法预测复合材料多钉连接的钉载分配：所述步骤d中采用单变量迭代的二分法预测复合材料多钉连接的钉载分配实现过程为：

根据背景技术中现有刚度法公式(4)计算螺栓刚度，公式(5)计算层板刚度、根据本发明提出的附加刚度计算方法计算附加刚度，假设螺栓1的变形，将假设变形和三种刚度代入本发明提出的新的刚度方法方程，以计算假设变形时钉传载荷，采用现有刚度法的一维搜索求解方法，用单变量迭代的二分法获得螺栓1的真实变形及钉载分配。

本发明与现有技术相比的优点在于：本发明考虑了旁路载荷引起的孔洞拉伸变形对螺栓载荷分布的影响，得到了一种预测复合材料多钉连接钉载分配的考虑孔拉伸变形影响的刚度方法，使得复合材料多螺栓连接结构钉载分配确定更加精确。并且随着螺栓数量的增加，预测精度的提高更加显著。

附图说明

图1是本发明一种基于修正刚度方法的复合材料多钉连接钉载分配预测方法的实现流程图；

图2是本发明复合材料多钉连接的弹簧模型图；

图3是本发明附加刚度计算推导模型。

具体实施方式

下面结合附图以及具体实施方式方式进一步说明本发明。

如图1所示，本发明一种基于修正刚度方法的复合材料多钉连接钉载分配预测方法的具体实现为：

1、首先提出与孔拉伸变形相对应的附加刚度的概念。

引入旁路载荷引起的孔拉伸变形以及对应于孔拉伸变形的附加刚度i是对螺栓个数的描述，a、b分别指代螺栓连接的上下两个层合板。考虑旁路载荷引起的孔拉伸变形后，建立复合材料多钉连接的弹簧模型如图2所示，可以得到孔拉伸变形及其对应的附加刚度之间的关系为：

其中，fi^a和fi^b(i＝1,···,n)分别是a、b两板在第i和第(i+1)个螺栓之间的层板内力。

2、提出新的刚度方法方程。

根据图2所示的弹簧模型，可以列出，考虑孔拉伸变形后的第i及第i+1个螺栓之间的变形协调方程为：

其中，δi是第i个螺栓的纵向(沿加载方向)变形；δi^a和δi^b分别是a、b两板在第i和第(i+1)个螺栓之间的纵向变形；是b板第i个螺栓孔的拉伸变形；是a板第i+1个螺栓孔的拉伸变形。

代入板和孔拉伸的刚度，将板和孔拉伸变形用力表示，则第i及第i+1个螺栓之间的变形协调方程为变为：

其中，fi^a和fi^b(i＝1,···,n)分别是a、b两板在第i和第(i+1)个螺栓之间的层板内力；和分别是a、b两板在第i和第(i+1)个螺栓之间的层板等效刚度，是对应于孔拉伸变形的附加刚度。

根据上述含孔板刚度的新的变形协调方程，得到第n个螺栓的变形协调方程为：

将得到的两个考虑孔拉伸变形的螺栓之间的变形协调方程，替换现有刚度法计算公式中的变形协调条件，得到考虑孔拉伸变形的新的刚度方法方程为：

3、提出附加刚度的计算方法。

为了得到附加刚度的计算方法，研究一块同时存在钉传载荷和旁路载荷的孔板(如图3(a)所示)，两个螺栓孔的钉传载荷分别为fd和fe；层合板左边的旁路载荷为fc–fe；右边的旁路载荷为fc+fd；两个螺栓孔中间的旁路载荷为fc。两个螺栓孔中间的层合板变形(如图3(b)所示)由三部分组成：钉传载荷fd引起的孔挤压变形(如图3(c)所示)；旁路载荷fc引起的层合板拉伸变形(如图3(d)所示)；旁路载荷fc引起的孔拉伸变形(如图3(e)所示)。其中孔挤压变形是包含在螺栓刚度中的，不需赘述。图3(d)和图3(e)的变形的和等价于孔板(如图3(f)所示)在fc下的拉伸变形，其中一部分可分解为长度为l的无孔板的拉伸变形(如图3(g)所示)(该部分的层板刚度即为现有刚度法公式(5)中的)，则另一部分的变形为可以得到图3(f)中的变形关系为：

因此，附加刚度的计算公式为：

其中，kh为所研究孔板的刚度，没有理论公式，采用有限元求解得到；为第i和第(i+1)个螺栓之间的层板等效刚度，采用现有刚度法中对层板刚度的计算公式求解。

4、根据背景技术中现有刚度法公式(4)计算螺栓刚度，公式(5)计算层板刚度、根据本发明提出的附加刚度计算方法计算附加刚度，假设螺栓1的变形，将假设变形和三种刚度代入本发明内容中提出的新的刚度方法方程，以计算假设变形时钉传载荷，采用现有刚度法的一维搜索求解方法，用单变量迭代的二分法获得螺栓1的真实变形及钉载分配。本发明未详细阐述部分属于本领域公知技术。