CN109388868A - 一种考虑钉孔接触及钉弯曲的钉载计算方法 - Google Patents

一种考虑钉孔接触及钉弯曲的钉载计算方法 Download PDF

Info

Publication number
CN109388868A
CN109388868A CN201811114408.6A CN201811114408A CN109388868A CN 109388868 A CN109388868 A CN 109388868A CN 201811114408 A CN201811114408 A CN 201811114408A CN 109388868 A CN109388868 A CN 109388868A
Authority
CN
China
Prior art keywords
nail
calculation method
nail hole
plane
node
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
Application number
CN201811114408.6A
Other languages
English (en)
Other versions
CN109388868B (zh
Inventor
陈栎同
付益战
张新平
董凤武
朱海兵
温峰
朱远见
郭宇
朱虹
王旭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Aircraft Industry Co Ltd
Original Assignee
AVIC Shaanxi Aircraft Industry Group Corp Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by AVIC Shaanxi Aircraft Industry Group Corp Ltd filed Critical AVIC Shaanxi Aircraft Industry Group Corp Ltd
Priority to CN201811114408.6A priority Critical patent/CN109388868B/zh
Publication of CN109388868A publication Critical patent/CN109388868A/zh
Application granted granted Critical
Publication of CN109388868B publication Critical patent/CN109388868B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供了一种考虑钉孔接触及钉弯曲的钉载计算方法,可用于航天航空、机械设计制造领域。目前确定钉载的主要方法有试验测量和利用有限元方法进行钉载计算,前者耗时长,需要大量准备工作以及设备,后者对硬件要求更高,且接触分析对网格要求高,计算不易收敛,耗时很长等特点,在实际工程中往往不易实现。为解决现有技术的不足,本发明提出一种考虑钉孔接触特性以及钉弯曲的钉载计算方法。在充分利用现有二维模型解决方法的基础上,通过引入离散刚体接触面并定义节点自由度约束关系实现考虑钉孔接触特性以及连接钉弯曲的钉载计算方法。使工程计算中能更快速得到相对准确的钉载分配结果,提高了计算效率。

Description

一种考虑钉孔接触及钉弯曲的钉载计算方法
技术领域
本发明提供了一种考虑钉孔接触及钉弯曲的钉载计算方法,可用于航天航空、机械设计制造领域。
背景技术
目前确定钉载的主要方法有两种。一种为试验测量。在使用该方法时,由于受材料性能分散性、环境、试验手段等因素影响,往往需要设置多组试验件才能得到相对准确的测量结果,其耗时长,需要大量准备工作以及设备,尤其对于大尺寸连接试验,成本高,不易实现;另一种方法即利用有限元方法进行钉载计算,尤其在结构初步设计阶段,该方法被大量使用。
对于利用有限元方法进行机械连接仿真计算钉载分配,主要有三种方法:1)建立二维板壳单元模拟连接件,使用梁单元、杆单元或者fastener单元模拟连接钉。该方法建模速度快、计算周期短,但该方法无法模拟实际连接中的钉孔接触、配合间隙、摩擦力等的影响,其结果只能作为初步设计的依据;2)在1)方法的基础上,利用间隙元模拟钉孔配合间隙、接触的影响。但由于其在多钉或者大规模钉组分析中,钉载分配复杂,无法确定实际接触面,所以建模阶段需要对每个连接钉的整个圆周上设置间隙元,计算规模太大,不易建模,耗时多。一般仅用于钉孔细节分析中;3)利用三维有限元单元建立实体模型进行计算。该方法几乎可模拟接触、装配间隙、预紧力等几乎所有机械连接特性,但由于实体模型规模很大,对硬件要求更高,且接触分析对网格要求高,计算不易收敛,耗时很长等特点,在实际工程中往往不易实现。
发明内容
为解决现有技术的不足,本发明提出一种考虑钉孔接触特性以及钉弯曲的钉载计算方法。在充分利用现有二维模型解决方法的基础上,通过引入离散刚体接触面并定义节点自由度约束关系实现考虑钉孔接触特性以及连接钉弯曲的钉载计算方法。使工程计算中能更快速得到相对准确的钉载分配结果,提高了计算效率。
本发明考虑钉孔接触以及钉弯曲的钉载计算方法,用于对两个或多个连接板的钉载通过有限元计算获得,该方法将连接钉划分为被连接板中位面和连接板之间的中间平面打断的若干子线段,且每个被打断的子线段被划分为包括两个梁单元的网格;在连接板钉孔位置建立若干离散的圆柱刚性面,尺寸与钉孔一致,其中心为参考点,圆柱刚性面的参考点与线段在中位面的打断处节点重合,且两者在6个自由度进行耦合;并以线段在中位面的打断处节点作为控制节点,与对应的钉孔节点在轴向及周向的自由度进行耦合。
所述的连接板中位面为连接板的几何中位面。
对较大尺寸的连接板,选取部分区域参与计算。
将连接钉划分为若干子线段的步骤中,以最外侧的两个中位面之间的钉孔圆心连线建立线段,并将线段在中位面以及各连接板的中间平面打断,从而获得子线段。
所述对应的钉孔节点选取方法为:以钉孔中心为圆心,以钉孔为内环、以2倍的钉孔半径作外环划定中位面上的圆环区域,内环节点即为对应的钉孔节点。
将圆环区域以4节点板单元划分网格,在中位面的其他区域以板单元划分自由网格,后者的网格密度小于前者的网格密度;为圆柱刚性面以四节点划分网格。
有限元计算之前,按装配关系将连接板、连接钉、圆柱刚性面进行装配,在圆柱刚性面与钉孔内环节点之间建立接触关系,在接触属性中设置接触面摩擦系数、钉孔装配关系。
在建好的模型中,将约束端试件对应位置施加边界约束条件,在加载端施加对应载荷。
有限元计算时,若接触面出现穿透,则将初始增量步减小。
有限元计算完成后,提取刚性面与孔壁接触力或者对应位置梁单元节点力即为钉载。
有益效果:本发明提出考虑钉孔接触及钉弯曲的钉载计算方法,通过二维板单元、梁元、离散刚体的综合使用,实现了在二维模型仿真中考虑钉孔接触、摩擦力以及装配过盈的钉载分配计算方法。
充分挖掘了二维有限元模型的可能性。可操作性、应用性强,在实现同等目标前提下,建模更为简便,计算时间较短,提高了计算效率。本方法通过二维板单元、梁元、离散刚体的综合使用,并通过节点自由度约束耦合的方法模拟机械连接的方法较新颖,在不引入复杂三维模型条件下,得到了更为精确的钉载结果。
附图说明
图1是装配完成后的模型示意图;
图2是圆柱刚性面的几何外形;
图3是钉孔周围的局部细化网格。
具体实施方式
将此方法应用于复合材料与金属板双面对接机械连接分析中,分析内容为钉载计算,连接形式为单排连接,每排三个连接钉。进行对接的两个板为复合材料,上下盖板为金属材料。以下步骤基于有限元软件ABAQUS。
步骤1:按照复合材料与金属连接板几何尺寸,保留其钉孔的几何特征,取其中位面。
步骤2:以最外侧的两个中位面之间的钉孔圆心连线建立线段,并将线段在中位面以及各连接板的中间平面打断,从而获得子线段。
步骤3:建立离散的圆柱刚性面,尺寸与钉孔一致,其中心为参考点。
步骤4:在中位面钉孔处,以钉孔中心为圆心,以钉孔为内环、以2倍的钉孔半径作外环划定中位面上的圆环区域。
步骤5:定义材料属性。将圆环区域以4节点板单元划分网格,单元类型为S4R,在中位面的其他区域以板单元划分自由网格,单元类型为S4R或S3R,后者的网格密度小于前者的网格密度;将每个被打断的子线段划分为包括两个梁元的网格,并赋予梁元属性,单元类型为B31。
步骤6:按装配关系将连接板、连接钉、圆柱刚性面进行装配,将圆柱刚性面的参考点与线段在中位面的打断处节点重合。
步骤7:在interaction模块中,使用coupling命令将圆柱刚性面的参考点与线段在中位面的打断处节点的6个自由度进行耦合。
步骤8:在interaction模块中,使用coupling命令将线段在中位面的打断处节点作为控制节点,与对应的钉孔节点在轴向及周向的自由度进行耦合。
步骤9:在interaction模块中,使用create interaction命令在圆柱刚性面与钉孔内环节点之间建立接触关系,在接触属性中设置接触面摩擦系数、钉孔装配关系,此处摩擦系数定义为0.1,钉孔间隙定义为10um。
步骤10:在建好的模型中,在一侧复合材料板远端自由边施加位移约束,另一侧复合材料板远端施加载荷。
步骤11:进行有限元计算后,提取圆柱刚性面与孔壁接触力或者对应位置梁单元节点力即为钉载。

Claims (10)

1.一种考虑钉孔接触及钉弯曲的钉载计算方法,用于对两个或多个连接板的钉载通过有限元计算获得,其特征在于,该方法将连接钉划分为被连接板中位面和连接板之间的中间平面打断的若干子线段,且每个被打断的子线段被划分为包括两个梁单元的网格;在连接板钉孔位置建立若干离散的圆柱刚性面,尺寸与钉孔一致,其中心为参考点,圆柱刚性面的参考点与线段在中位面的打断处节点重合,且两者在6个自由度进行耦合;并以线段在中位面的打断处节点作为控制节点,与对应的钉孔节点在轴向及周向的自由度进行耦合。
2.根据权利要求1所述的钉载计算方法,其特征在于:所述的连接板中位面为连接板的几何中位面。
3.根据权利要求1所述的钉载计算方法,其特征在于:对较大尺寸的连接板,选取部分区域参与计算。
4.根据权利要求1所述的钉载计算方法,其特征在于:将连接钉划分为若干子线段的步骤中,以最外侧的两个中位面之间的钉孔圆心连线建立线段,并将线段在中位面以及各连接板的中间平面打断,从而获得子线段。
5.根据权利要求1所述的钉载计算方法,其特征在于:所述对应的钉孔节点选取方法为:以钉孔中心为圆心,以钉孔为内环、以2倍的钉孔半径作外环划定中位面上的圆环区域,内环节点即为对应的钉孔节点。
6.根据权利要求5所述的钉载计算方法,其特征在于:将圆环区域以4节点板单元划分网格,在中位面的其他区域以板单元划分自由网格,后者的网格密度小于前者的网格密度;为圆柱刚性面以四节点划分网格。
7.根据权利要求1所述的钉载计算方法,其特征在于:有限元计算之前,按装配关系将连接板、连接钉、圆柱刚性面进行装配,在圆柱刚性面与钉孔内环节点之间建立接触关系,在接触属性中设置接触面摩擦系数、钉孔装配关系。
8.根据权利要求7所述的钉载计算方法,其特征在于:在建好的模型中,将约束端试件对应位置施加边界约束条件,在加载端施加对应载荷。
9.根据权利要求8所述的钉载计算方法,其特征在于:有限元计算时,若接触面出现穿透,则将初始增量步减小。
10.根据权利要求8所述的钉载计算方法,其特征在于:有限元计算完成后,提取刚性面与孔壁接触力或者对应位置梁单元节点力即为钉载。
CN201811114408.6A 2018-09-25 2018-09-25 一种考虑钉孔接触及钉弯曲的钉载计算方法 Active CN109388868B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811114408.6A CN109388868B (zh) 2018-09-25 2018-09-25 一种考虑钉孔接触及钉弯曲的钉载计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811114408.6A CN109388868B (zh) 2018-09-25 2018-09-25 一种考虑钉孔接触及钉弯曲的钉载计算方法

Publications (2)

Publication Number Publication Date
CN109388868A true CN109388868A (zh) 2019-02-26
CN109388868B CN109388868B (zh) 2023-04-07

Family

ID=65419032

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811114408.6A Active CN109388868B (zh) 2018-09-25 2018-09-25 一种考虑钉孔接触及钉弯曲的钉载计算方法

Country Status (1)

Country Link
CN (1) CN109388868B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109813587A (zh) * 2019-03-08 2019-05-28 南昌航空大学 一种用于铆钉偏心连接接头钉载分配研究的试验件
CN112800539A (zh) * 2021-01-15 2021-05-14 中国商用飞机有限责任公司北京民用飞机技术研究中心 一种钉载分布预测方法及系统

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7207448B1 (en) * 2004-07-06 2007-04-24 Marino Jr Frank Modular storage system
CN102013349A (zh) * 2010-12-17 2011-04-13 常熟开关制造有限公司(原常熟开关厂) 三位置自动转换开关的电器工作状态指示装置
CN102350686A (zh) * 2011-09-26 2012-02-15 邵文彪 一种利用u型钉绕接伞骨和伞面的装置
CN102519648A (zh) * 2011-12-22 2012-06-27 北京航空航天大学 单剪连接结构钉载矢量测量方法和测量仪
CN202292676U (zh) * 2011-09-26 2012-07-04 邵文彪 一种利用u型钉绕接伞骨和伞面的装置
CN102622472A (zh) * 2012-02-27 2012-08-01 西北工业大学 分析复合材料机械连接钉载和应力分布及应力水平的方法
CN103020350A (zh) * 2012-12-10 2013-04-03 中国飞机强度研究所 一种紧固件载荷变形曲线数值仿真方法
CN103335886A (zh) * 2013-06-25 2013-10-02 北京航空航天大学 一种基于三参数特征曲线的复合材料多钉双剪连接失效预测方法
CN103386593A (zh) * 2013-07-18 2013-11-13 哈尔滨工业大学 一种碳纤维增强树脂基复合材料与金属的连接方法
CN103593542A (zh) * 2013-12-03 2014-02-19 北京航空航天大学 一种考虑间隙和拧紧力矩的复合材料螺栓连接结构钉载分配确定方法
CN106193304A (zh) * 2016-07-11 2016-12-07 宋启宽 一种钢绞线连续悬挂体系
CN106370724A (zh) * 2016-11-22 2017-02-01 中国航空工业集团公司北京航空材料研究院 一种用于航空盘环件涡流自动检测的探头装置及扫查方法
CN106768916A (zh) * 2016-11-21 2017-05-31 北京航空航天大学 一种采用均匀钉载进行复合材料螺栓连接结构最终失效载荷预测的方法
US20170238684A1 (en) * 2016-02-20 2017-08-24 James Anthony Cureington Magnetic multipurpose wearable retainer
CN107729603A (zh) * 2017-09-04 2018-02-23 中国商用飞机有限责任公司北京民用飞机技术研究中心 复合材料壁板螺栓连接的非线性分析方法
CN108460237A (zh) * 2018-04-08 2018-08-28 大连理工大学 一种考虑螺纹实际切向力状态的螺栓连接结构松动有限元仿真方法

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7207448B1 (en) * 2004-07-06 2007-04-24 Marino Jr Frank Modular storage system
CN102013349A (zh) * 2010-12-17 2011-04-13 常熟开关制造有限公司(原常熟开关厂) 三位置自动转换开关的电器工作状态指示装置
CN102350686A (zh) * 2011-09-26 2012-02-15 邵文彪 一种利用u型钉绕接伞骨和伞面的装置
CN202292676U (zh) * 2011-09-26 2012-07-04 邵文彪 一种利用u型钉绕接伞骨和伞面的装置
CN102519648A (zh) * 2011-12-22 2012-06-27 北京航空航天大学 单剪连接结构钉载矢量测量方法和测量仪
CN102622472A (zh) * 2012-02-27 2012-08-01 西北工业大学 分析复合材料机械连接钉载和应力分布及应力水平的方法
CN103020350A (zh) * 2012-12-10 2013-04-03 中国飞机强度研究所 一种紧固件载荷变形曲线数值仿真方法
CN103335886A (zh) * 2013-06-25 2013-10-02 北京航空航天大学 一种基于三参数特征曲线的复合材料多钉双剪连接失效预测方法
CN103386593A (zh) * 2013-07-18 2013-11-13 哈尔滨工业大学 一种碳纤维增强树脂基复合材料与金属的连接方法
CN103593542A (zh) * 2013-12-03 2014-02-19 北京航空航天大学 一种考虑间隙和拧紧力矩的复合材料螺栓连接结构钉载分配确定方法
US20170238684A1 (en) * 2016-02-20 2017-08-24 James Anthony Cureington Magnetic multipurpose wearable retainer
CN106193304A (zh) * 2016-07-11 2016-12-07 宋启宽 一种钢绞线连续悬挂体系
CN106768916A (zh) * 2016-11-21 2017-05-31 北京航空航天大学 一种采用均匀钉载进行复合材料螺栓连接结构最终失效载荷预测的方法
CN106370724A (zh) * 2016-11-22 2017-02-01 中国航空工业集团公司北京航空材料研究院 一种用于航空盘环件涡流自动检测的探头装置及扫查方法
CN107729603A (zh) * 2017-09-04 2018-02-23 中国商用飞机有限责任公司北京民用飞机技术研究中心 复合材料壁板螺栓连接的非线性分析方法
CN108460237A (zh) * 2018-04-08 2018-08-28 大连理工大学 一种考虑螺纹实际切向力状态的螺栓连接结构松动有限元仿真方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴存利 等: "铆钉载荷变形曲线有限元数值分析", 《机械科学与技术》 *
张冠彪: "碳纤维层合板与铝合金板机械连接钉载分配研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109813587A (zh) * 2019-03-08 2019-05-28 南昌航空大学 一种用于铆钉偏心连接接头钉载分配研究的试验件
CN109813587B (zh) * 2019-03-08 2021-04-02 南昌航空大学 一种用于铆钉偏心连接接头钉载分配研究的试验件
CN112800539A (zh) * 2021-01-15 2021-05-14 中国商用飞机有限责任公司北京民用飞机技术研究中心 一种钉载分布预测方法及系统

Also Published As

Publication number Publication date
CN109388868B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
CN107330146B (zh) 一种考虑转动效应的岩质边坡极限承载力分析上限法
CN108153981A (zh) 一种基于有限元分析的复合材料机身加筋壁板结构后屈曲分析方法
JP6255491B2 (ja) コンピュータ支援設計(cad)モデルをシミュレートするためのファスナの自動作成
JP2022521907A (ja) ハイブリッド繊維複合材料の板巻きシェル構造に対する高速協調最適化方法
CN105653869B (zh) 一种超临界水堆堆芯三维瞬态性能分析方法
CN111460568B (zh) 一种混凝土重力坝运行期裂缝扩展判别方法
Meakin et al. Unsteady aerodynamic simulation of static and moving bodies using scalable computers
CN103267507B (zh) 基于有限元分析提取机械结构平面的平面度误差的方法
CN107958124B (zh) 一种基于应力分布模式优化的混合尺度耦合方法
CN102819632A (zh) 一种钢框架结构节点多尺度有限元模型建模方法
CN109388868A (zh) 一种考虑钉孔接触及钉弯曲的钉载计算方法
CN104298814B (zh) 一种基于参数误差累积的齿轮系统性能可靠度计算方法
CN110222420A (zh) 一种粘铆连接结构的失效判断方法
CN109726435B (zh) 一种接头钉群载荷计算的方法
CN104298823A (zh) 一种高、低温天平分析方法及分析系统
CN108491612A (zh) 为复合管液压胀形工艺提供选材方案的有限元模拟方法
CN104143215A (zh) 索膜结构的信息化模型设计方法
Suri et al. Quick and easy manufacturing systems analysis using MANUPLAN.
CN116776678A (zh) 一种基于混合单元的结构应力计算方法
CN113673023B (zh) 一种飞机复杂块体结构疲劳载荷谱确定方法
Swoboda et al. Simulation of arch dam–foundation interaction with a new friction interface element
Slater CFD methods for computing the performance of supersonic inlets
CN106156537A (zh) 基于麦夸特算法的径向基无网格软组织数据的力反馈模型建模方法
Lawrence An overview of three approaches to multidisciplinary aeropropulsion simulations
Gatzke et al. MACGS-A zonal grid generation system for complex aero-propulsion configurations

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 723213 Liulin Town, Chenggu County, Hanzhong City, Shaanxi Province

Patentee after: Shaanxi Aircraft Industry Co.,Ltd.

Address before: 723213 box 34, Hanzhong City, Shaanxi Province

Patentee before: Shaanxi Aircraft INDUSTRY(GROUP) Co.,Ltd.