CN113468783B - 一种基于弹簧阻尼系统的动态铆接有限元模拟方法 - Google Patents

一种基于弹簧阻尼系统的动态铆接有限元模拟方法 Download PDF

Info

Publication number
CN113468783B
CN113468783B CN202110708465.2A CN202110708465A CN113468783B CN 113468783 B CN113468783 B CN 113468783B CN 202110708465 A CN202110708465 A CN 202110708465A CN 113468783 B CN113468783 B CN 113468783B
Authority
CN
China
Prior art keywords
riveting
rivet
model
top iron
finite element
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.)
Active
Application number
CN202110708465.2A
Other languages
English (en)
Other versions
CN113468783A (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.)
Chengdu Aircraft Industrial Group Co Ltd
Original Assignee
Chengdu Aircraft Industrial Group Co 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 Chengdu Aircraft Industrial Group Co Ltd filed Critical Chengdu Aircraft Industrial Group Co Ltd
Priority to CN202110708465.2A priority Critical patent/CN113468783B/zh
Publication of CN113468783A publication Critical patent/CN113468783A/zh
Application granted granted Critical
Publication of CN113468783B publication Critical patent/CN113468783B/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]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

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)
  • Vibration Dampers (AREA)

Abstract

本发明公开了一种基于弹簧阻尼系统的动态铆接有限元模拟方法,包括:建立风动铆接原理模型:风动铆接反铆载荷实质为铆卡在气压推动下,撞击铆钉钉头使整体向前移动,与顶铁共同挤压使铆钉变形后,铆卡回退至初始相对位置,并进行下一个循环,多次循环载荷的作用下,铆钉墩头彻底成型,并停止铆接,分析铆接过程各部件的动态响应,建立等效模型;有限元模型建立并模拟:在建模软件中建立系统的实体模型,实体模型截取铆接结构一部分,网格划分,铆钉变形较大区域的网格进行了局部细化,创建分析作业并提交分析进行计算,完成有限元计算得到铆钉变形结果。该有限元模型填补了风动铆接仿真技术的空白,简便易行、计算准确,工程实际应用效果极好。

Description

一种基于弹簧阻尼系统的动态铆接有限元模拟方法
技术领域
本发明涉及金属塑性成型技术领域,特别是涉及一种基于弹簧阻尼系统的动态铆接有限元模拟方法。
背景技术
铆钉连接是一种常用的机械连接方法,连接原理是将铆钉穿过连接件与被连接件后,通过恒压力或冲击力的作用墩粗铆钉钉杆,从而起到连接效果。铆钉连接是一种不可逆的连接方法。现今航空航天等制造行业中,大部分铆钉采用风动铆枪进行反铆连接。操作时,风动铆枪将气压力转换为铆枪铆卡的循环冲击力,并作用在铆钉钉盖上,与此同时,铆钉钉杆端放置有顶铁,人工握持并与钉杆贴合。在循环载荷的作用下,铆钉钉杆成型。与压铆、拉铆相比,风动铆接以其装配效率高,可施工环境多等优点,在装配连接领域具有不可替代的作用。
因风动铆接载荷变化复杂、顶持过程难以模拟、边界条件难以设定等特点,至今为止还未有出现有效的模型,用以模拟风动铆接状态下铆钉成型过程,导致很多围绕风动铆接下铆钉变形分析与应力状态分析的相关工作难以展开。实际制造过程中,采用风动铆接工艺的铆钉占铆钉总数的95%,建立风动铆接状态下,铆钉成型过程的有限元模型刻不容缓。为填补该技术的空白,本发明提出一种基于弹簧阻尼系统的动态铆接仿真技术,具有重要的工程价值。
发明内容
为了填补风动铆接有限元模拟技术的空白,本发明提出一种基于弹簧阻尼系统的动态铆接有限元模拟方法。深入研究风动铆接的动态过程,设计了铆接脉冲载荷与顶铁顶持的等效模型,并采用有限元软件实现,用以研究风动铆接状态下铆钉成型过程,为工艺参数优化奠定基础。
本发明的目的是通过以下技术方案实现的:
一种基于弹簧阻尼系统的动态铆接有限元模拟方法,其特征在于,包括以下步骤:
步骤1、建立风动铆接原理模型:风动铆接反铆载荷实质为铆卡在气压推动下,撞击铆钉钉头使整体向前移动,与顶铁共同挤压使铆钉变形后,铆卡回退至初始相对位置,并进行下一个循环,多次循环载荷的作用下,铆钉墩头彻底成型,并停止铆接,分析铆接过程各部件的动态响应,建立等效模型;
步骤2、有限元模型建立并模拟:在建模软件中建立系统的实体模型,实体模型截取铆接结构一部分,网格划分时,铆钉变形较大区域的网格进行了局部细化,创建分析作业并提交分析进行计算,完成有限元计算后得到铆钉变形结果。
优选的,所述步骤1中,等效模型包括铆卡的载荷施加方式模型和人手顶铁的状态模型。
优选的,所述铆卡的载荷施加方式模型中,铆卡的载荷施加方式等效为铆卡的往复运动。
优选的,所述往复运动的位移满足函数:,其中:ξ-分析步长,/>,/>-初始位移,/>-最大位移,/>-单次循环载荷开始时间,/>-单次循环载荷结束时间。
优选的,所述人手顶铁的状态模型中,将人对顶铁的作用等效为弹簧阻尼系统,铆钉受载,传递给顶铁时,顶铁产生冲量,并以初速度向后运动,此时人手臂肌肉绷紧,持续减速顶铁至最大位移后,推进顶铁至铆钉墩头处,并进行下一个循环。
优选的,所述人手顶铁的状态模型中,将人对顶铁的作用等效为弹簧阻尼系统,符合如下方程: ,其中:/>-顶铁质量,/>-等效阻尼系数,/>-等效弹性系数,/>-加速度,/>-速度;公式满足初始条件:/>
优选的,所述建立系统的实体模型中包括确定材料参数,材料种类包括:连接件为铝合金2024-O;被连接件为铝合金7050-T7451;铆钉为铝合金2A10;顶铁为钨钢S1;材料性能参数均通过单向拉伸试验获得。
优选的,所述建立系统的实体模型中,截取铆接结构一部分,建立等效模型,其中连接件为板件,被连接件为T型结构,铆钉模型参考HB6316-4x9结构。
优选的,所述建立系统的实体模型中还包括建立顶铁系统模型,顶铁系统模型包括顶铁系统顶铁、弹簧阻尼器与刚性板件,顶铁系统顶铁与刚性板件间通过弹簧阻尼器相连。
优选的,所述建立系统的实体模型中还包括载荷设定及边界条件:定义铆卡的循环位移,模拟铆卡施加在铆钉上的冲击载荷,顶铁系统中施加位移至刚性板件,模拟手持顶铁顶持;连接件顶端与被连接件底端设置对称约束边界条件。
本技术方案的有益效果如下:
本发明提出一种基于弹簧阻尼系统的动态铆接有限元模拟方法。通过分析风动铆接铆卡施加载荷的方式,建立了载荷的等效模型。通过分解手持顶铁顶持时,每个阶段的运动状态,建立了弹簧阻尼系统等效模型。基于以上原理,截取实体结构中的一部分,在有限元软件中进行了实现,对铆钉变形较大区域进行了网格细化,计算完成后得到铆钉变形量、铆钉残余应力与结构装配应力。该有限元模型填补了风动铆接仿真技术的空白,简便易行、计算准确,工程实际应用效果极好。
附图说明
下面将结合说明书附图和具体实施方式对本发明作进一步的详细说明,其中:
图1 为本发明基于风动铆接过程的等效原理图;
图2 为本发明的有限元模型;
图3 为本发明有限元模型网格;
图4 为本发明中铆钉局部细化网格图;
图5 为本发明中模型成型后应力云图;
图6为本发明中铆钉应力云图剖视图;
图7为本发明中铆钉位移云图;
图中标记:1、铆卡;2、连接件;3、被连接件;4、铆钉;5、顶铁;6、弹簧阻尼器;7、刚性板件。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
如图1-图4所示,本铆接模型截取了飞机进气道铆钉铆接结构的一段,用于模拟进气道内铆钉的风动铆接过程。以ABAQUS软件为例对过程进行实现。飞机进气道结构具有夹层薄、刚性弱、曲率大等特点,蒙皮的铆接质量,极大地影响进气道蒙皮连接可靠性。蒙皮为板材,连接件2为铝合金2024-O,被连接件为机加隔框,材料为铝合金7050-T7451,铆钉4材料为铝合金2A10。
(一)建立风动铆接原理模型
风动铆接反铆载荷实质为,铆卡1在气压推动下,撞击铆钉4钉头使整体向前移动,与顶铁5共同挤压使铆钉4成型后,铆卡1回退至初始相对位置,并进行下一个循环。多次循环载荷的作用下,铆钉4墩头彻底成型,并停止铆接。与此同时,单次铆接循环过程中铆钉4受载并传递给顶铁5,顶铁5产生冲量,并以初速度向后运动。此时人手持顶铁5限制其运动,持续减速顶铁5至最大位移后,推进顶铁5与铆钉4墩头接触,并进行下一个循环。
为简化模型,将铆卡1在气压作用下产生载荷的方式等效为铆卡1的往复运动,运动幅值设为。人手持顶铁5的动作可等效为顶铁的弹簧阻尼结构,即顶铁5通过弹簧阻尼器6连接至刚性板件7,铆接过程中运动刚体向铆钉4方向移动。单次载荷循环下,顶铁系统满足如下方程:/> ,其中:/>-顶铁质量,/>-等效阻尼系数,/>-等效弹性系数;公式满足初始条件:。
(二)模型参数确定
通过在试件上粘接应变片,测量风动铆枪铆接过程中实际载荷,载荷参数设置为:位移幅值;载荷脉冲宽度/>0.6/>s;载荷间隔/>0.03s。因ABAQUS动态响应计算需要,往复运动的位移函数符合ABAQUS软件中Smooth step,即:/>,其中:/>,/>-初始位移,/>-最大位移,/>-单次循环载荷开始时间,/>-单次循环载荷结束时间。
顶铁系统中,顶铁质量,正常情况下,操作者可提供的顶持力约100N/mm。风动铆接过程中顶铁5顶持的过程应属于过阻尼振动状态,满足条件:/>
(三)建立有限元模型及网格划分
基于理论模型,在ABAQUS软件中建立有限元模型,按实际尺寸绘制各个零组件数模,并在Assemble模块中组合,其中铆卡1、顶铁5初始状态与铆钉4钉头、墩头贴合。顶铁系统模型中包含顶铁5与刚性板件7,顶铁5与刚性板件7用wire命令连接,便于后续赋属性。在Property模块中设置材料参数,并赋予相应零件,材料参数见下表,所有零件模型截面属性选用实体单元Homogeneous类型。为防止计算过程中网格畸变过大,在铆钉4变形量较大区域,网格进行局部细化,单元类型为C3D8R。
(四)载荷设定及边界条件
定义铆卡1沿轴向最大位移量为0.75mm,其余五个方向的自由度均设为0。为让铆卡1实现循环脉冲位移,在Amplitudes模块下,按载荷脉冲宽度0.6/>s;载荷间隔/>0.03s,设置铆卡1位移幅值曲线。约束连接件2与被连接件3两端位移,等效为悬臂结构。在Interaction模块下定义弹簧阻尼系统连接件属性,并赋予在顶铁5与刚性板件7间wire上,同时定义刚性板件7沿轴向的最大位移量为3mm。
(五)提交分析作业及后处理
如图5-图7所示,创建分析作业并提交分析进行计算,完成有限元计算后得到循环载荷下铆钉的变形结果与应力分布,所述创建分析作业是瞬态求解器模块下创建默认分析作业,提交分析作业为ABAQUS默认操作。
本发明基于弹簧阻尼系统的动态铆接仿真技术,填补了风动铆接仿真技术的空白。模型结果不仅可以准确得出铆钉在铆接过程中各阶段的变形量,同时可以得到所有零件的应力应变分布,在飞行器强度设计环节具有重要的指导意义。
综上所述,本领域的普通技术人员阅读本发明文件后,根据本发明的技术方案和技术构思无需创造性脑力劳动而作出的其他各种相应的变换方案,均属于本发明所保护的范围。

Claims (5)

1.一种基于弹簧阻尼系统的动态铆接有限元模拟方法,其特征在于,包括以下步骤:
步骤1、建立风动铆接原理模型:风动铆接反铆载荷实质为铆卡在气压推动下,撞击铆钉钉头使整体向前移动,与顶铁共同挤压使铆钉变形后,铆卡回退至初始相对位置,并进行下一个循环,多次循环载荷的作用下,铆钉墩头彻底成型,并停止铆接,分析铆接过程各部件的动态响应,建立等效模型;
等效模型包括铆卡的载荷施加方式模型和人手顶铁的状态模型;
所述铆卡的载荷施加方式模型中,铆卡的载荷施加方式等效为铆卡的往复运动;
所述往复运动的位移满足函数: ,其中:ξ-分析步长,/>,/>-初始位移,/>-最大位移,/>-单次循环载荷开始时间,/>-单次循环载荷结束时间;
所述人手顶铁的状态模型中,将人对顶铁的作用等效为弹簧阻尼系统,铆钉受载,传递给顶铁时,顶铁产生冲量,并以初速度向后运动,此时人手臂肌肉绷紧,持续减速顶铁至最大位移后,推进顶铁至铆钉墩头处,并进行下一个循环;
所述人手顶铁的状态模型中,将人对顶铁的作用等效为弹簧阻尼系统,符合如下方程: ,其中:/>-顶铁质量,/> -等效阻尼系数,/>-等效弹性系数,/>-加速度,/>-速度;公式满足初始条件/>
步骤2、有限元模型建立并模拟:在建模软件中建立系统的实体模型,实体模型截取铆接结构一部分,网格划分时,铆钉变形较大区域的网格进行了局部细化,创建分析作业并提交分析进行计算,完成有限元计算后得到铆钉变形结果。
2.根据权利要求1所述的一种基于弹簧阻尼系统的动态铆接有限元模拟方法,其特征在于:所述建立系统的实体模型中包括确定材料参数,材料种类包括:连接件为铝合金2024-O;被连接件为铝合金7050-T7451;铆钉为铝合金2A10;顶铁为钨钢S1;材料性能参数均通过单向拉伸试验获得。
3.根据权利要求2所述的一种基于弹簧阻尼系统的动态铆接有限元模拟方法,其特征在于:所述建立系统的实体模型中,截取铆接结构一部分,建立等效模型,其中连接件为板件,被连接件为T型结构,铆钉模型参考HB6316-4x9结构。
4.根据权利要求3所述的一种基于弹簧阻尼系统的动态铆接有限元模拟方法,其特征在于:所述建立系统的实体模型中还包括建立顶铁系统模型,顶铁系统模型包括顶铁系统顶铁、弹簧阻尼器与刚性板件,顶铁系统顶铁与刚性板件间通过弹簧阻尼器相连。
5.根据权利要求4所述的一种基于弹簧阻尼系统的动态铆接有限元模拟方法,其特征在于:所述建立系统的实体模型中还包括载荷设定及边界条件:定义铆卡的循环位移,模拟铆卡施加在铆钉上的冲击载荷,顶铁系统中施加位移至刚性板件,模拟手持顶铁顶持;连接件顶端与被连接件底端设置对称约束边界条件。
CN202110708465.2A 2021-06-25 2021-06-25 一种基于弹簧阻尼系统的动态铆接有限元模拟方法 Active CN113468783B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110708465.2A CN113468783B (zh) 2021-06-25 2021-06-25 一种基于弹簧阻尼系统的动态铆接有限元模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110708465.2A CN113468783B (zh) 2021-06-25 2021-06-25 一种基于弹簧阻尼系统的动态铆接有限元模拟方法

Publications (2)

Publication Number Publication Date
CN113468783A CN113468783A (zh) 2021-10-01
CN113468783B true CN113468783B (zh) 2023-11-10

Family

ID=77872914

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110708465.2A Active CN113468783B (zh) 2021-06-25 2021-06-25 一种基于弹簧阻尼系统的动态铆接有限元模拟方法

Country Status (1)

Country Link
CN (1) CN113468783B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114850386B (zh) * 2022-04-22 2023-02-03 成都飞机工业(集团)有限责任公司 一种飞机拉杆装配用铆接装置及其方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103455671A (zh) * 2013-08-27 2013-12-18 西北工业大学 电磁铆接接头疲劳寿命预测方法
CN107194057A (zh) * 2017-05-19 2017-09-22 西北工业大学 壁板铆接扭翘变形预测方法
CN107590318A (zh) * 2017-08-22 2018-01-16 株洲时代新材料科技股份有限公司 一种汽车推力杆热铆接工艺仿真分析方法
CN107876685A (zh) * 2017-11-14 2018-04-06 山东大学 一种超声振动辅助铆接方法及其应用方法
CN111898291A (zh) * 2020-06-30 2020-11-06 西北工业大学 一种基于子结构的大壁板连续铆接变形预测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103455671A (zh) * 2013-08-27 2013-12-18 西北工业大学 电磁铆接接头疲劳寿命预测方法
CN107194057A (zh) * 2017-05-19 2017-09-22 西北工业大学 壁板铆接扭翘变形预测方法
CN107590318A (zh) * 2017-08-22 2018-01-16 株洲时代新材料科技股份有限公司 一种汽车推力杆热铆接工艺仿真分析方法
CN107876685A (zh) * 2017-11-14 2018-04-06 山东大学 一种超声振动辅助铆接方法及其应用方法
CN111898291A (zh) * 2020-06-30 2020-11-06 西北工业大学 一种基于子结构的大壁板连续铆接变形预测方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A Simplified Finite Element Riveted Lap Joint Model in Structural Dynamic Analysis;Marco Daniel Malheiro Dourado 等;Advanced Materials Research;第1016卷;第185-191页 *
气动锤铆自动铆接的工艺研究;刘现伟 等;有色金属材料与工程;第38卷(第03期);第154-159页 *
铆钉材料对铆接变形影响的有限元分析;李艳 等;机床与液压;第41卷(第04期);第50-52页 *
飞机薄壁件批量铆接过程的模拟方法研究;王仲奇 等;航空制造技术(第Z1期);第34-37页 *
飞机装配铆接工艺减振外骨骼关键技术研究;范曙远;中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑(第04期);第C031-48页 *

Also Published As

Publication number Publication date
CN113468783A (zh) 2021-10-01

Similar Documents

Publication Publication Date Title
CN107729603B (zh) 复合材料壁板螺栓连接的非线性分析方法
CN105426595A (zh) 一种铝合金热弹塑性变形模拟用本构模型的建立方法
KR20080019253A (ko) 부재의 설계 방법, 부재의 설계 장치 및 컴퓨터 프로그램을 기록한 컴퓨터로 판독 가능한 기록 매체
CN113468783B (zh) 一种基于弹簧阻尼系统的动态铆接有限元模拟方法
CN106407582B (zh) 一种用于柔性机构静力学和动力学分析的高效计算方法
CN107679301B (zh) 一种分段式重载横梁缩比模型设计方法
CN105912823A (zh) 一种基于n-2-1定位的飞机蒙皮和长桁装配偏差分析方法
Li et al. Dynamic modeling and simulation of percussive impact riveting for robotic automation
Repetto et al. A finite element study of electromagnetic riveting
Aljibori et al. Finite element analysis of sheet metal forming process
CN111922316B (zh) 一种压铸机机铰系统的设计方法
Reddy Evaluation of Single Point Incremental Forming Process for Parabolic AA6082 Cups
CN112464401B (zh) 一种金属材料焊点的精准建模方法
Grujicic et al. Self-piercing riveting process and joint modeling and simulations
Reddy et al. A review on finite element simulations in metal forming
Chen et al. Analysis of an equivalent drawbead model for the finite element simulation of a stamping process
JP4305645B2 (ja) 板成形のシミュレーション方法
Ahmetoglu et al. Computer simulation for tool and process design in sheet forming
Adrian et al. Curating Datasets of Flexible Assemblies to Predict Spring-Back Behavior for Machine Learning Purposes
Aljibori Non-linear elastic-plastic behaviour of aluminium sheet metal using finite element analysis
Asadi et al. A simplified model to simulate crash behavior of honeycomb
CN108664727B (zh) 一种车身钣金凹陷快速修复模型的构建方法及系统
JP2021148742A (ja) 部分構造試験方法
Anwekar et al. Finite element simulation of single stage deep drawing process for determining stress distribution in drawn conical component
Xia et al. Springback Behavior of AA6111‐T4 with Split‐Ring Test

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