CN106295024A - 一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法 - Google Patents
一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法 Download PDFInfo
- Publication number
- CN106295024A CN106295024A CN201610670979.2A CN201610670979A CN106295024A CN 106295024 A CN106295024 A CN 106295024A CN 201610670979 A CN201610670979 A CN 201610670979A CN 106295024 A CN106295024 A CN 106295024A
- Authority
- CN
- China
- Prior art keywords
- bolt
- gap
- stiffness matrix
- vector
- represent
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
本发明是一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法,涉及复合材料螺栓连接技术。本方法建立结构节点单元模型,输入几何参数、材料参数和边界条件,完成螺栓连接模型的建立;再基于建立的螺栓连接模型,计算节点位移,在此基础上计算螺栓载荷分配。本方法基于EXCEL边框生成螺栓连接模型结构;对螺栓的单元刚度矩阵和相应的等效荷载向量考虑三种情况:不考虑摩擦和间隙的影响,考虑间隙的影响,考虑摩擦和间隙的影响。本发明结合弹簧模型计算效率高和有限元方法建模效率高的优点,避免了传统的有限元方法计算效率较低,考虑间隙和摩擦影响建模复杂的问题,相对传统弹簧模型具有较高的建模效率,对模型修改花费时间较少。
Description
技术领域
本发明涉及复合材料螺栓连接的技术,具体是一种考虑间隙及摩擦影响的复合材料载荷分布确定计算方法,实现复合材料螺栓连接模型的快速建立和载荷分布计算。
背景技术
复合材料具有比强度、比刚度高、热膨胀系数小、耐疲劳、抗腐蚀、制造周期短和维修方便等优点,被广泛应用于航空航天结构,特别是在主要承载部件中得到越来越多的应用。螺栓连接通常用于复合材料层板之间的连接,或复合材料和金属之间的连接,其可靠性高,承载能力强,拆卸方便,因此得到广泛的应用。然而由于螺栓连接会导致应力集中,连接可能因此破坏。因此,在初始设计阶段,进行螺栓连接复合结构的载荷分布的可靠预测是十分有必要的。
复合材料螺栓连接的分析方法有数值模拟和理论分析两种,虽然有限元方法可以模拟螺栓连接所有的因素的影响,但所需时间过多,特别是当摩擦和损伤的影响考虑在内时。弹簧方法是一种高效便于实施的螺栓载荷分配理论分析方法,在复合材料螺栓连接分析领域得到广泛的应用,同时可以考虑间隙、摩擦等多种因素。但由于初始设计阶段,复合材料螺栓连接模型通常根据实际需要进行调整,该方法虽对特定的模型计算效率较高,但不易重新建立计算模型。因此,需发展一种能够快速建立螺栓连接计算模型的方法。
发明内容
本发明针对现有的螺栓连接理论分析方法不易建立计算模型的缺陷,提出一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布确定方法,使既包含现有的理论分析方法计算效率高的优点,同时能够快速建立螺栓连接计算模型。
本发明提供的考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法,包括以下步骤:
步骤1,建立结构节点单元模型;
具体地,对待分析的螺栓连接结构,建立等效弹簧模型,然后基于EXCEL边框生成螺栓连接模型结构;在基于EXCEL边框生成螺栓连接模型结构时,水平边界表示板弹簧单元,垂直边框代表螺栓弹簧单元;
步骤2,输入几何参数、材料参数和边界条件,完成螺栓连接模型的建立;
所述的几何参数包括各螺栓间板的长度、宽度和厚度,所述的材料参数包括各板和螺栓的力学参数;所述的边界条件包括被约束的节点编号、施加力的节点编号以及力的大小,当考虑间隙和摩擦的影响时,还需输入各螺栓的间隙以及摩擦力的大小。
步骤3,基于建立的螺栓连接模型,计算节点位移,在此基础上计算螺栓载荷分配。
(3.1)确定所有的弹簧单元对应的节点和单元位置向量;(3.2)计算单元刚度矩阵和相应的等效荷载向量;(3.3)将所有螺栓和板的单元刚度矩阵以及荷载向量组合,形成系统的刚度矩阵和等效载荷向量;(3.4)将约束应用到系统的刚度矩阵和等效载荷向量;(3.5)在(3.4)得到的系统的刚度矩阵和等效载荷向量基础上建立刚度方程,得到节点位移。
本发明提供的方法结合弹簧模型计算效率高和有限元方法建模效率高的优点,同时可以考虑间隙及摩擦影响,其优点在于:
(1)与传统的有限元方法相比,计算效率高;
(2)与传统的有限元方法相比,无需建立复杂的有限元模型,较为容易的考虑间隙及摩擦的影响;
(3)与已有螺栓载荷分布理论分析方法相比,能够快速建立螺栓连接计算模型。
附图说明
图1是本发明的复合材料螺栓连接载荷分布计算方法的整体步骤示意图;
图2是一个三螺栓单剪搭接连接结构示意图;
图3是对图1所示的三螺栓单剪搭接连接结构建立的等效弹簧模型图;
图4是本发明方法对图1所示结构基于EXCEL边框生成的模型结构图;
图5是本发明的载荷分布计算方法的流程图;
图6是本发明算例模型结构图;
图7是本发明算例的有限元模型示意图。
具体实施方式
下面结合附图来对本发明的实现方案进行具体说明。
本发明提供的考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法,基于弹簧模型和有限元方法,实现了螺栓连接计算模型的快速建立,并在此基础上计算了螺栓载荷分配。本发明方法整体步骤如图1所示,下面对各步骤详细说明。
步骤1,建立结构节点单元模型。结构模型通过定义EXCEL表格中的单元格边框进行。
如针对图2所示的三螺栓单剪搭接连接结构,其等效弹簧模型如图3所示。图2中,Skin_plate表示上板,Splice_plate表示下板,Bolt表示螺栓,skin load表示外力P。图3中,建立7个节点,标号1~6表示三个螺栓与板的6个连接点,标号7表示外力作用点。KB1,KB2,KB3分别代表螺栓1、2、3的刚度。KSPL1代表螺栓1和螺栓2之间下板的刚度,KSPL2表示螺栓2和螺栓3之间下板的刚度,KSKIN1表示螺栓1和螺栓2之间上板的刚度,KSKIN2表示螺栓2和螺栓3之间上板的刚度。KSPL_END表示连接螺栓3和固支点之间的下板的刚度。KSKIN_END表示连接螺栓1和外力作用点之间的上板的刚度,P表示外力。
基于EXCEL边框生成的三螺栓单剪搭接连接模型结构,如图4所示,水平边界表示板弹簧单元,垂直边框代表螺栓弹簧单元。板弹簧单元可通过使用VBA函数“Cells.Borders(xlEdgeBottom).LineStyle”识别。螺栓弹簧单元可通过使用VBA函数“Cells.Borders(xlEdgeRight).LineStyle”识别。因此,连接板和螺栓的节点可以通过结合上述两个函数识别。其中,Cells表示单元格,Borders表示边框,xlEdgeBottom表示单元格下边框,xlEdgeRight表示单元格右边框,LineStyle表示返回边框的线型。
如图4所示的模型结构图中,板弹簧单元为节点1-2,2-3,3-4,5-6,6-7,7-8之间的单元,螺栓弹簧单元为节点2-5,3-6,4-7之间的单元,约束和负载被施加到节点8和节点1上。
步骤2,输入几何参数、材料参数和边界条件。需要输入的几何参数有各螺栓间板的长度、宽度和厚度,需要输入的材料参数有各板和螺栓的力学参数。需要输入的边界条件有被约束的节点编号,施加力的节点编号以及力的大小。如果考虑间隙和摩擦的影响,同时还需输入各螺栓的间隙以及摩擦力的大小。至此,建立了螺栓连接模型。
步骤3,基于上述螺栓连接模型,计算节点位移,在此基础上计算螺栓载荷分配。
步骤3包括如下子步骤3.1~3.5,如图5所示,下面具体说明各步骤。
步骤3.1,获取有限元节点和单元位置向量。
所述的有限元节点为螺栓连接模型结构中的节点,是根据等效弹簧模型中各弹簧单元确定的节点。首先确定所有的弹簧单元对应的节点数。对图3所示三螺栓单剪搭接连接结构,确定共8个有限元节点,如图4所示。然后,基于各弹簧单元对应的节点编号确定弹簧单元刚度矩阵在初始系统刚度矩阵中的位置,构成单元刚度矩阵的位置向量。
步骤3.2,计算单元刚度矩阵和相应的等效荷载向量。
如果不考虑摩擦和间隙的影响,螺栓的单元刚度矩阵和相应的等效荷载向量可以计算并存储为以下形式。
其中,kb表示螺栓刚度,前两列表示螺栓刚度向量Kb,后一列表示对应的等效荷载向量Pb。
如果考虑间隙的影响,则螺栓的单元刚度矩阵和相应的等效载荷向量存储为如下形式。
其中,cb表示螺栓间隙。
如果考虑间隙和摩擦的影响,则螺栓的单元刚度矩阵和相应的等效载荷向量存储为如下形式。
其中,ufb表示由于摩擦导致的位移,PFRICb表示摩擦力。
螺栓间板的单元刚度矩阵和相应的元素等效荷载向量可以计算并存储为以下形式。
其中,kp表示该段板的刚度,前两列表示该段板的刚度矩阵Kp,后一列表示对应的等效力向量Pp,Pp为零向量。
步骤3.3,将螺栓和板的单元刚度矩阵和荷载向量组合成系统的刚度矩阵和等效载荷向量。
步骤3.4,将约束应用到系统的刚度矩阵和等效载荷向量。删除原系统的刚度矩阵和荷载向量的对应的行和列,得到约束应用后的系统的刚度矩阵和系统等效荷载向量。
步骤3.5,通过求解刚度矩阵和系统等效荷载向量的基础上建立的刚度方程,得到节点位移。螺栓或板的内力是相应的螺栓或板的刚度和对应的节点位移差的乘积。
本发明提供的算例,如图6所示,为三螺栓双剪搭接连接结构,根据对称性,可以取对称的一半进行分析,则将其转化为等效单剪搭接模型如图3所示,对应的基于EXCEL边框生成的模型结构图如图4所示,图7为对应的用于验证的有限元模型。初始参数值如表1、表2所示,基于相关公式可以计算得到螺栓等效刚度为36.76KN/mm。
表1几何参数
表2材料参数
最终螺栓载荷分配结果及有限元对比如表3所示,最大误差为0.42%,证明本发明方法是有效的。
表3该方法与有限元结果的对比
Claims (4)
1.一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法,其特征在于,实现步骤如下:
步骤1,建立结构节点单元模型,具体是:对待分析的螺栓连接结构,建立等效弹簧模型,然后基于EXCEL边框生成螺栓连接模型结构;在基于EXCEL边框生成螺栓连接模型结构时,水平边界表示板弹簧单元,垂直边框代表螺栓弹簧单元;
步骤2,输入几何参数、材料参数和边界条件,完成螺栓连接模型的建立;
所述的几何参数包括各螺栓间板的长度、宽度和厚度,所述的材料参数包括各板和螺栓的力学参数;所述的边界条件包括被约束的节点编号、施加力的节点编号以及力的大小,当考虑间隙和摩擦的影响时,还需输入各螺栓的间隙以及摩擦力的大小;
步骤3,基于建立的螺栓连接模型,计算节点位移,在此基础上计算螺栓载荷分配;
具体包括:(3.1)确定所有的弹簧单元对应的节点和单元位置向量;(3.2)计算单元刚度矩阵和相应的等效荷载向量;(3.3)将所有螺栓和板的单元刚度矩阵以及荷载向量组合,形成系统的刚度矩阵和等效载荷向量;(3.4)将约束应用到系统的刚度矩阵和等效载荷向量;(3.5)在(3.4)得到的系统的刚度矩阵和等效载荷向量基础上建立刚度方程,得到节点位移。
2.根据权利要求1所述的一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法,其特征在于,所述的在基于EXCEL边框生成螺栓连接模型结构时,板弹簧单元通过使用VBA函数“Cells.Borders(xlEdgeBottom).LineStyle”识别;螺栓弹簧单元通过使用VBA函数“Cells.Borders(xlEdgeRight).LineStyle”识别;连接板和螺栓的节点通过结合所述的两个函数识别;其中,Cells表示单元格,Borders表示边框,xlEdgeBottom表示单元格下边框,xlEdgeRight表示单元格右边框,LineStyle表示返回边框的线型。
3.根据权利要求1所述的一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法,其特征在于,所述的(3.2)中,在计算螺栓的单元刚度矩阵时,分为三种情况:
(1)不考虑摩擦和间隙的影响,螺栓的单元刚度矩阵和相应的等效荷载向量计算并存储为以下形式:
其中,kb表示螺栓刚度,Kb表示螺栓刚度向量,Pb表示相应的等效荷载向量;
(2)考虑间隙的影响,螺栓的单元刚度矩阵和相应的等效载荷向量存储为如下形式:
其中,cb表示螺栓间隙;
如果考虑间隙和摩擦的影响,螺栓的单元刚度矩阵和相应的等效载荷向量存储为如下形式:
其中,ufb表示由于摩擦导致的位移,PFRICb表示摩擦力。
4.根据权利要求1所述的一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法,其特征在于,所述的(3.2)中,螺栓间板的单元刚度矩阵和相应的等效荷载向量计算并存储为以下形式:
其中,kp表示板的刚度,前两列表示该段板的刚度矩阵Kp,后一列表示对应的等效载荷向量Pp,Pp为零向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610670979.2A CN106295024B (zh) | 2016-08-15 | 2016-08-15 | 一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610670979.2A CN106295024B (zh) | 2016-08-15 | 2016-08-15 | 一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106295024A true CN106295024A (zh) | 2017-01-04 |
CN106295024B CN106295024B (zh) | 2019-10-01 |
Family
ID=57671457
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610670979.2A Active CN106295024B (zh) | 2016-08-15 | 2016-08-15 | 一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106295024B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111209707A (zh) * | 2020-02-21 | 2020-05-29 | 山东交通学院 | 承受压剪组合作用的摩擦型螺栓连接节点及方法及系统 |
CN112084616A (zh) * | 2020-05-19 | 2020-12-15 | 北京航空航天大学 | 一种考虑几何非线性的预测复合材料螺旋结构压缩刚度与压缩强度的方法 |
CN112528393A (zh) * | 2020-11-05 | 2021-03-19 | 中国船舶工业集团公司第七0八研究所 | 一种船舶轴系连接法兰动力学分析建模方法 |
CN113282003A (zh) * | 2021-05-19 | 2021-08-20 | 南京航空航天大学 | 一种考虑界面接触的磁悬浮轴承-转子系统建模方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102063552A (zh) * | 2011-01-07 | 2011-05-18 | 西安交通大学 | 一种螺栓连接件刚度的确定方法 |
CN102735411A (zh) * | 2012-07-10 | 2012-10-17 | 株洲时代新材料科技股份有限公司 | 汽车横置复合材料板弹簧静态刚度测试装置及测试方法 |
CN103218483A (zh) * | 2013-04-02 | 2013-07-24 | 北京航空航天大学 | 一种基于梁—弹簧模型的螺纹连接强度计算方法 |
CN103593542A (zh) * | 2013-12-03 | 2014-02-19 | 北京航空航天大学 | 一种考虑间隙和拧紧力矩的复合材料螺栓连接结构钉载分配确定方法 |
WO2016065391A1 (en) * | 2014-10-30 | 2016-05-06 | Garock Pty Ltd | Ground support apparatus |
-
2016
- 2016-08-15 CN CN201610670979.2A patent/CN106295024B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102063552A (zh) * | 2011-01-07 | 2011-05-18 | 西安交通大学 | 一种螺栓连接件刚度的确定方法 |
CN102735411A (zh) * | 2012-07-10 | 2012-10-17 | 株洲时代新材料科技股份有限公司 | 汽车横置复合材料板弹簧静态刚度测试装置及测试方法 |
CN103218483A (zh) * | 2013-04-02 | 2013-07-24 | 北京航空航天大学 | 一种基于梁—弹簧模型的螺纹连接强度计算方法 |
CN103593542A (zh) * | 2013-12-03 | 2014-02-19 | 北京航空航天大学 | 一种考虑间隙和拧紧力矩的复合材料螺栓连接结构钉载分配确定方法 |
WO2016065391A1 (en) * | 2014-10-30 | 2016-05-06 | Garock Pty Ltd | Ground support apparatus |
Non-Patent Citations (1)
Title |
---|
姚熊亮 等: "船用复合材料层合板螺栓孔处疲劳寿命实验研究", 《中国造船》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111209707A (zh) * | 2020-02-21 | 2020-05-29 | 山东交通学院 | 承受压剪组合作用的摩擦型螺栓连接节点及方法及系统 |
CN112084616A (zh) * | 2020-05-19 | 2020-12-15 | 北京航空航天大学 | 一种考虑几何非线性的预测复合材料螺旋结构压缩刚度与压缩强度的方法 |
CN112528393A (zh) * | 2020-11-05 | 2021-03-19 | 中国船舶工业集团公司第七0八研究所 | 一种船舶轴系连接法兰动力学分析建模方法 |
CN112528393B (zh) * | 2020-11-05 | 2022-09-16 | 中国船舶工业集团公司第七0八研究所 | 一种船舶轴系连接法兰动力学分析建模方法 |
CN113282003A (zh) * | 2021-05-19 | 2021-08-20 | 南京航空航天大学 | 一种考虑界面接触的磁悬浮轴承-转子系统建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106295024B (zh) | 2019-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106295024A (zh) | 一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法 | |
CN106096257B (zh) | 一种非线性索单元分析方法及系统 | |
CN108153981A (zh) | 一种基于有限元分析的复合材料机身加筋壁板结构后屈曲分析方法 | |
CN107451307B (zh) | 一种多尺度计算复杂复合材料结构等效刚度矩阵的方法 | |
US10824778B2 (en) | Computer-implemented method for simulating behavior of a screw connection between metal components of engineering structures | |
Szekrényes | Application of Reddy's third-order theory to delaminated orthotropic composite plates | |
Overgaard et al. | A methodology for the structural analysis of composite wind turbine blades under geometric and material induced instabilities | |
Bai et al. | Improved design of extended end-plate connection allowing for prying effects | |
Liang et al. | A smeared stiffener based reduced-order modelling method for buckling analysis of isogrid-stiffened cylinder | |
Bhaskar et al. | An elasticity approach for simply-supported isotropic and orthotropic stiffened plates | |
Peng et al. | Numerical investigation of the effects of structural geometric and material nonlinearities on limit-cycle oscillation of a cropped delta wing | |
Kovvali et al. | Verification of the variational-asymptotic sectional analysis for initially curved and twisted beams | |
Shen et al. | Modeling of peridynamic beams and shells with transverse shear effect via interpolation method | |
Pan et al. | Modal interactions of a dual-joint cylindrical shell system via Nonlinear Normal Modes | |
Ramaprasad et al. | Exact solutions for thin-walled composite open section beams using a unified state space coupled field formulation | |
Stapleton et al. | Adaptive shape functions and internal mesh adaptation for modeling progressive failure in adhesively bonded joints | |
CN116663362B (zh) | 一种航空航天复合材料薄壁结构装配力学建模方法 | |
Ordonneau et al. | A simplified modal analysis of a single lap bonded joint using the macro-element technique | |
Abdelrahman et al. | Generalized line-element formulations for geometrically nonlinear analysis of nonsymmetric tapered steel members with warping and Wagner effects | |
Yang et al. | A semi-analytical method for determining the strain energy release rate of cracks in adhesively-bonded single-lap composite joints | |
He et al. | A high order finite strip transfer matrix method for buckling analysis of single-branched cross-section thin-walled members | |
Xu-xi et al. | A trigonometric analytical solution of simply supported horizontally curved composite I-beam considering tangential slips | |
Blachowski et al. | Graph based discrete optimization in structural dynamics | |
CN116702381B (zh) | 一种螺栓连接非线性的等效和振动响应计算方法 | |
CN110852015B (zh) | 一种陶瓷基复合材料模态的非线性计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170104 Assignee: Beijing northern sky long hawk UAV Technology Co.,Ltd. Assignor: BEIHANG University Contract record no.: X2021990000039 Denomination of invention: A calculation method for load distribution of Composite Bolted Connections considering the effect of clearance and friction Granted publication date: 20191001 License type: Exclusive License Record date: 20210119 |
|
EE01 | Entry into force of recordation of patent licensing contract |