CN109635316B - 一种螺栓仿真分析方法 - Google Patents

一种螺栓仿真分析方法 Download PDF

Info

Publication number
CN109635316B
CN109635316B CN201811272391.7A CN201811272391A CN109635316B CN 109635316 B CN109635316 B CN 109635316B CN 201811272391 A CN201811272391 A CN 201811272391A CN 109635316 B CN109635316 B CN 109635316B
Authority
CN
China
Prior art keywords
bolt
thread
finite element
model
grid
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
CN201811272391.7A
Other languages
English (en)
Other versions
CN109635316A (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.)
Ningbo University
Original Assignee
Ningbo University
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 Ningbo University filed Critical Ningbo University
Priority to CN201811272391.7A priority Critical patent/CN109635316B/zh
Publication of CN109635316A publication Critical patent/CN109635316A/zh
Application granted granted Critical
Publication of CN109635316B publication Critical patent/CN109635316B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种螺栓仿真分析方法，特点是包括以下步骤：(1)将螺栓纵向摆放，螺栓头朝下，建立螺栓有限元模型的各节点受力方向模型；(2)若对称轴方向与正方向不一致则模型对应旋转；(3)将上述模型导入Fortran编译器中，若螺栓为标准件则输入螺栓标准件代号；若螺栓为非标准件则输入螺杆信息，运行Fortran编程语言得到螺栓表面有限元网格以及各结点受力方向；(4)将螺栓表面有限元网格以及各结点受力方向坐标向量输入到HyperMesh软件中，运行TCL宏，输出螺栓输出螺栓3D二阶实体网格和2D螺纹线网格各节点受力方向，优点是网格简单，仿真精度高且耗时少。

Description

一种螺栓仿真分析方法
技术领域
本发明涉及一种螺栓仿真分析方法。
背景技术
螺纹联接广泛存在于各类机械结构中，是最常用的紧固手段之一。例如，发动机中各个部件装配中重要的联接方式是螺栓联接。盾构隧道是地铁在设计、施工中的重要环节，其衬砌为装配式结构，接缝多以螺栓连接，从目前已进行的整环足尺试验来看，结构的破坏始于接缝处，接缝的构造形式及极限承载能力在整环结构承载中占有重要地位。对于接头的设计，一般近似地将接头面看作是连续的钢筋混凝土截面，将螺栓看作受拉钢筋，螺栓有单排直螺栓、单排弯螺栓、斜螺栓等几种连接形式。
对于螺栓部件的仿真过程中，一般会采用精细化模型和简单化模型。其中，简单化模型网格数目少，网格为直筒状，计算简单，但精度较低；而精细化模型仿真精度高，但网格数目庞大，计算困难，影响效率。在CAE仿真过程中，盾构隧道衬砌结构繁多，螺栓数目更是极多，故需要通过一种仿真分析方法，以达到仿真精度高、耗时少两方面的要求。
发明内容
本发明所要解决的技术问题是提供一种网格简单，仿真精度高且耗时少的螺栓仿真分析方法。
本发明解决上述技术问题所采用的技术方案为：一种螺栓仿真分析方法，包括以下步骤：(1)将螺栓纵向摆放，螺栓头朝下，建立螺栓有限元模型的各节点受力方向模型：
Figure BDA0001846285420000011
其中
Figure BDA0001846285420000021
为牙口角度，β螺纹线上升角度；θ为螺栓有限元网格所占角度，
Figure BDA0001846285420000022
n为一圈的网格数，n取值为12-36内的任一整数；
Figure BDA0001846285420000023
相应为坐标系x、y、z的单位方向向量；
(2)设定螺栓头部中心坐标(x0，y0，z0)，螺栓从下到上第一圈螺纹线中心坐标(xc，yc，zc)；若|xc-x0|≤0.001，|yc-y0|≤0.001，zc≥z0时，不需旋转螺栓有限元模型的各节点受力方向模型；否则需要旋转螺栓有限元模型的各节点受力方向模型，螺栓有限元模型的各节点受力方向模型旋转后如下：
Figure BDA0001846285420000024
即
Figure BDA0001846285420000025
其中：
Figure BDA0001846285420000026
Figure BDA0001846285420000027
Figure BDA0001846285420000028
Figure BDA0001846285420000029
Figure BDA00018462854200000210
Figure BDA00018462854200000211
a＝H(2,2)×H(3,3)-H(3,2)×H(2,3)
b＝H(2,3)×H(3,1)-H(2,1)×H(3,3)
c＝H(2,1)×H(3,2)-H(2,2)×H(3,1)
Figure BDA0001846285420000031
Figure BDA0001846285420000032
Figure BDA0001846285420000033
(x7，y7，z7)为(xc，yc，zc)坐标与(x，y，z)坐标连线上的任意点；
(3)将上述模型导入Fortran编译器中，并设置一圈螺纹的网格数n；若螺栓为标准件，则在Fortran编译器中输入螺栓标准件代号，螺栓头部中心坐标即(x0，y0，z0)以及螺栓从下到上第一圈螺纹线中心坐标即(xc，yc，zc)；若螺栓为非标准件，则在Fortran编译器中输入螺头高度、螺头半径、螺杆高度、螺栓半径、螺纹圈数、螺纹间距、牙口角度、螺纹线上升角度、螺栓头部中心坐标即(x0，y0，z0)以及螺栓从下到上第一圈螺纹线中心坐标即(xc，yc，zc)，运行Fortran编程语言，得到螺栓表面有限元网格以及螺栓螺纹各节点受力方向坐标向量；其中各类螺栓标准件的螺栓信息预先存储到Fortran编译器中，所述的螺栓信息包括螺头高度、螺头半径、螺杆高度、螺栓半径、螺纹圈数、螺距、牙口角度、螺纹线上升角度；
(4)将步骤(3)得到的螺栓表面有限元网格以及螺栓螺纹各节点受力方向坐标向量输入到HyperMesh软件中，在HyperMesh软件中运行TCL宏，输出螺栓3D二阶实体网格、2D螺纹线网格各节点受力方向和螺栓接触面网格。
步骤(1)中所述的螺栓有限元模型的各节点受力方向模型建立方法如下：
1)设定
Figure BDA0001846285420000034
垂直面ABCD，A的坐标为(rsinθ,rcosθ,z)，B的坐标为
Figure BDA0001846285420000035
其中P为螺栓线程，r为螺杆半径，则
Figure BDA0001846285420000036
2)D的坐标为
Figure BDA0001846285420000037
则
Figure BDA0001846285420000038
即：
Figure BDA0001846285420000039
3)
Figure BDA00018462854200000310
得
X方向
Figure BDA0001846285420000041
Y方向
Figure BDA0001846285420000042
Z方向
Figure BDA0001846285420000043
故：
Figure BDA0001846285420000044
式中：
Figure BDA0001846285420000045
为牙口角度；β螺纹线上升角度；θ为螺栓有限元网格所占角度。
与现有技术相比，本发明的优点在于：本发明一种螺栓仿真分析方法，用于有限元分析的基于螺纹几何特征的螺旋线螺纹模型，在仿真模型中兼顾计算效率与精度，根据螺纹和螺孔接触总是在螺纹的正法线方向的特点，建立合理有效的螺纹接触优化算法。基于螺纹螺旋线形状的特点，运用Hypermesh脚本语言，完成螺栓及螺孔的螺纹网格自动生成，提高效率与网格质量，减少人工前处理时间，在保证计算准确的同时，有效地提高了计算效率。主要的实施就是通过根据实际的螺栓螺纹咬合的特征，编辑Fortran语言，生成用于有限元分析的具有螺纹特征效果的圆筒形螺栓网格模型。
综上所述，本发明运用Fortran编程语言，建立一套程序，根据简单化网格节点，自动计算螺纹线上各点所受力的方向，输出螺栓有限元模型和螺纹线，从而实现简单化网格的精细化计算，达到精细化模型的仿真精度，在保证计算准确的同时，有效地提高了计算效率，具有网格简单，仿真精度高，且效率高，创建螺纹疲劳系数分布图，在螺栓紧固件结构方面运用效果显著。
附图说明
图1为公式推理螺栓简化图；
图2为导入HyperMesh软件生成的网格直观图，其中A为螺栓3D二阶实体网格，B为2D螺纹线网格各节点受力方向；
图3为导入HyperMesh软件生成的螺栓接触面网格；
图4为发动机螺栓位置图；
图5为螺栓精细化模型结果图；
图6为本发明仿真方法所得模型结果图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
具体实施例
一种螺栓仿真分析方法，包括以下步骤：
(1)将螺栓纵向摆放，螺栓头朝下，建立螺栓有限元模型的各节点受力方向模型：
Figure BDA0001846285420000051
其中
Figure BDA0001846285420000052
为牙口角度，β螺纹线上升角度；θ为螺栓有限元网格所占角度，
Figure BDA0001846285420000053
n为一圈的网格数，n取值为12-36内的任一整数；
Figure BDA0001846285420000054
相应为坐标系x、y、z的单位方向向量；
公式具体推导过程为：
如图1所示，
Figure BDA0001846285420000055
垂直面ABCD
设A的坐标为(rsinθ,rcosθ,z)
1.B的坐标为
Figure BDA0001846285420000056
其中P为螺栓线程，r为螺杆半径；
Figure BDA0001846285420000057
2.D的坐标为
Figure BDA0001846285420000058
即：
Figure BDA0001846285420000059
Figure BDA00018462854200000510
得：
X方向
Figure BDA0001846285420000061
Y方向
Figure BDA0001846285420000062
Z方向
Figure BDA0001846285420000063
故：
Figure BDA0001846285420000064
式中：
Figure BDA0001846285420000065
为牙口角度；β螺纹线上升角度；θ为螺栓有限元网格所占角度；
(2)设定螺栓头部中心坐标(x0，y0，z0)，螺栓从下到上第一圈螺纹线中心坐标(xc，yc，zc)；若|xc-x0|≤0.001，|yc-y0|≤0.001，zc≥z0时，不需旋转螺栓有限元模型的各节点受力方向模型；否则需要旋转螺栓有限元模型的各节点受力方向模型，螺栓有限元模型的各节点受力方向模型旋转后如下：
Figure BDA0001846285420000066
即
Figure BDA0001846285420000067
推导公式如下：以下设定螺栓头部中心坐标(x0，y0，z0)，螺栓从下到上第一圈螺纹线中心坐标(xc，yc，zc)；(x7，y7，z7)为(xc，yc，zc)坐标与(x，y，z)坐标连线上的任意点；5、10为选取的任意数值，可用任意自然数替换；
1.判断：|xc-x0|≤0.001,|yc-y0|≤0.001,zc≥z0时，不需旋转；否则需要旋转；
2.旋转螺栓；
当|xc-x0|≤0.001,|yc-y0|≤0.001时：
Figure BDA0001846285420000068
当|xc-x0|≤0.001,|yc-y0|＞0.001,|zc-z0|≤0.001时：
Figure BDA0001846285420000071
当|xc-x0|≤0.001，|yc-y0|＞0.001，|zc-z0|＞0.001时：
Figure BDA0001846285420000072
当|xc-x0|＞0.001，|yc-y0|≤0.001时：
Figure BDA0001846285420000073
当|xc-x0|＞0.001，|yc-y0|＞0.001，|zc-z0|≤0.001时：
Figure BDA0001846285420000074
当|xc-x0|＞0.001，|yc-y0|＞0.001，|zc-z0|＞0.001时：
Figure BDA0001846285420000075
Figure BDA0001846285420000076
Figure BDA0001846285420000077
Figure BDA0001846285420000078
Figure BDA0001846285420000079
Figure BDA00018462854200000710
Figure BDA00018462854200000711
a＝H(2，2)×H(3，3)-H(3，2)×H(2，3)
b＝H(2,3)×H(3,1)-H(2,1)×H(3,3)
c＝H(2,1)×H(3,2)-H(2,2)×H(3,1)
Figure BDA0001846285420000081
Figure BDA0001846285420000082
Figure BDA0001846285420000083
(3)将上述模型导入Fortran编译器中，并设置一圈螺纹的网格数n；若螺栓为标准件，则在Fortran编译器中输入螺栓标准件代号，螺栓头部中心坐标即(x0，y0，z0)以及螺栓从下到上第一圈螺纹线中心坐标即(xc，yc，zc)；若螺栓为非标准件，则在Fortran编译器中输入螺头高度、螺头半径、螺杆高度、螺栓半径、螺纹圈数、螺距、牙口角度、螺纹线上升角度、螺栓头部中心坐标即(x0，y0，z0)以及螺栓从下到上第一圈螺纹线中心坐标即(xc，yc，zc)，运行Fortran编程语言，得到螺栓表面有限元网格以及螺栓螺纹各节点受力方向坐标向量；其中各类螺栓标准件的螺栓信息预先存储到Fortran编译器中，螺栓信息包括螺头高度、螺头半径、螺杆高度、螺栓半径、螺纹圈数、螺距、牙口角度、螺纹线上升角度；螺头高度、螺头半径是为了生成螺头部分表面网格；螺杆高度、螺栓半径为了生成螺杆表面网格；螺纹第一圈中心坐标、螺纹圈数、螺纹间距是为了确定计算螺栓受力方向的区域；
(4)将步骤(3)得到的螺栓表面有限元网格以及螺栓螺纹各节点受力方向坐标向量输入到HyperMesh软件中，在HyperMesh软件中运行TCL宏，输出螺栓3D二阶实体网格、2D螺纹线网格各节点受力方向(如图2所示)和螺栓接触面网格(如图3所示，即螺栓受力区域与施力区域，即螺纹受力与其他部件、螺栓头部与其他部件的接触面)。
对比试验：在耐久性实验中螺纹疲劳失效的运用。主轴承盖与发动机机体连接螺栓位置如图4所示，螺栓样式为M9。运用CAE仿真方法得到的精细化模型，螺纹网格节点共1100000，仿真计算耗时7天，仿真结果如图5所示；运用该仿真方法，网格节点共680000，仿真计算耗时1天，仿真结果如图6所示。两方法仿真结果对比可知，两者应力分布图相近；精细化模型最大应力为343MPa，其次为253MPa，运用本发明螺栓仿真分析法仿真结果最大应力为312MPa，其次为226MPa，两者数值结果相近，应力分布图位置相同。
上述说明并非对本发明的限制，本发明也并不限于上述举例。本技术领域的普通技术人员在本发明的实质范围内，做出的变化、改型、添加或替换，也应属于本发明的保护范围。

Claims (2)

1.一种螺栓仿真分析方法，其特征在于包括以下步骤：
(1)将螺栓纵向摆放，螺栓头朝下，建立螺栓有限元模型的各节点受力方向模型：
Figure FDA0001846285410000011
其中
Figure FDA0001846285410000012
为牙口角度，β螺纹线上升角度；θ为螺栓有限元网格所占角度，
Figure FDA0001846285410000013
n为一圈的网格数，n取值为12-36内的任一整数；
Figure FDA0001846285410000014
相应为坐标系x、y、z的单位方向向量；
(2)设定螺栓头部中心坐标(x0，y0，z0)，螺栓从下到上第一圈螺纹线中心坐标(xc，yc，zc)；若|xc-x0|≤0.001，|yc-y0|≤0.001，zc≥z0时，不需旋转螺栓有限元模型的各节点受力方向模型；否则需要旋转螺栓有限元模型的各节点受力方向模型，螺栓有限元模型的各节点受力方向模型旋转后如下：
Figure FDA0001846285410000015
即
Figure FDA0001846285410000016
其中：
Figure FDA0001846285410000017
Figure FDA0001846285410000018
Figure FDA0001846285410000019
Figure FDA00018462854100000110
Figure FDA00018462854100000111
Figure FDA00018462854100000112
a＝H(2,2)×H(3,3)-H(3,2)×H(2,3)，
b＝H(2,3)×H(3,1)-H(2,1)×H(3,3)，
c＝H(2,1)×H(3,2)-H(2,2)×H(3,1)，
Figure FDA0001846285410000021
Figure FDA0001846285410000022
Figure FDA0001846285410000023
(x7，y7，z7)为(xc，yc，zc)坐标与(x，y，z)坐标连线上的任意点；
(3)将上述模型导入Fortran编译器中，并设置一圈螺纹的网格数n；若螺栓为标准件，则在Fortran编译器中输入螺栓标准件代号，螺栓头部中心坐标即(x0，y0，z0)以及螺栓从下到上第一圈螺纹线中心坐标即(xc，yc，zc)；若螺栓为非标准件，则在Fortran编译器中输入螺头高度、螺头半径、螺杆高度、螺栓半径、螺纹圈数、螺距、牙口角度、螺纹线上升角度、螺栓头部中心坐标即(x0，y0，z0)以及螺栓从下到上第一圈螺纹线中心坐标即(xc，yc，zc),运行Fortran编程语言，得到螺栓表面有限元网格以及螺栓螺纹各节点受力方向坐标向量；其中各类螺栓标准件的螺栓信息预先存储到Fortran编译器中，所述的螺栓信息包括螺头高度、螺头半径、螺杆高度、螺栓半径、螺纹圈数、螺距、牙口角度和螺纹线上升角度；
(4)将步骤(3)得到的螺栓表面有限元网格以及螺栓螺纹各节点受力方向坐标向量输入到HyperMesh软件中，在HyperMesh软件中运行TCL宏，输出螺栓3D二阶实体网格、2D螺纹线网格各节点受力方向和螺栓接触面网格。
2.根据权利要求1所述的一种螺栓仿真分析方法，其特征在于：步骤(1)中所述的螺栓有限元模型的各节点受力方向模型建立方法如下：
1)设定
Figure FDA0001846285410000024
垂直面ABCD，A的坐标为(rsinθ,rcosθ,z)，B的坐标为
Figure FDA0001846285410000025
其中P为螺栓线程，r为螺杆半径，则
Figure FDA0001846285410000031
2)D的坐标为
Figure FDA0001846285410000032
则
Figure FDA0001846285410000033
即：
Figure FDA0001846285410000034
3)
Figure FDA0001846285410000035
得
X方向
Figure FDA0001846285410000036
Y方向
Figure FDA0001846285410000037
Z方向
Figure FDA0001846285410000038
故：
Figure FDA0001846285410000039
式中：
Figure FDA00018462854100000310
为牙口角度；β螺纹线上升角度；θ为螺栓有限元网格所占角度。
CN201811272391.7A 2018-10-30 2018-10-30 一种螺栓仿真分析方法 Active CN109635316B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811272391.7A CN109635316B (zh) 2018-10-30 2018-10-30 一种螺栓仿真分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811272391.7A CN109635316B (zh) 2018-10-30 2018-10-30 一种螺栓仿真分析方法

Publications (2)

Publication Number Publication Date
CN109635316A CN109635316A (zh) 2019-04-16
CN109635316B true CN109635316B (zh) 2022-12-06

Family

ID=66066829

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811272391.7A Active CN109635316B (zh) 2018-10-30 2018-10-30 一种螺栓仿真分析方法

Country Status (1)

Country Link
CN (1) CN109635316B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112115565B (zh) * 2020-09-28 2022-06-03 重庆理工大学 螺栓螺纹牙高度优化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102332047A (zh) * 2011-10-09 2012-01-25 西安交通大学 一种可实现六面体网格划分的螺栓有限元参数化建模方法
CN106202639A (zh) * 2016-06-29 2016-12-07 北京强度环境研究所 一种mj螺栓及螺母有限元网格建模方法
CN108197418A (zh) * 2018-03-14 2018-06-22 上海理工大学 一种模拟螺纹配合的六面体有限元网格划分方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102332047A (zh) * 2011-10-09 2012-01-25 西安交通大学 一种可实现六面体网格划分的螺栓有限元参数化建模方法
CN106202639A (zh) * 2016-06-29 2016-12-07 北京强度环境研究所 一种mj螺栓及螺母有限元网格建模方法
CN108197418A (zh) * 2018-03-14 2018-06-22 上海理工大学 一种模拟螺纹配合的六面体有限元网格划分方法

Also Published As

Publication number Publication date
CN109635316A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN110580376B (zh) 一种基于Revit和Dynamo的桩基础模型创建方法
CN103914590A (zh) 一种电力铁塔三维实体模型生成方法
CN107133383B (zh) 一种抛物线双曲拱坝拱圈编绘方法
US11549490B2 (en) Reinforcement assembly for wind turbine tower
CN106971046B (zh) 一种高位收水冷却塔中央竖井设计方法
CN112081285B (zh) 一种索结构预应力态拉索长度确定方法
CN109635316B (zh) 一种螺栓仿真分析方法
CN112986731B (zh) 计及地震不确定性的电气互联系统韧性评估与提升方法
CN105354382B (zh) 一种建立输电线路塔线体系有限元模型的方法
Bhooshan et al. Use of sub-division surfaces in architectural form-finding and procedural modelling
CN112836318A (zh) 基于代理模型的海上风机支撑结构优化设计方法及系统
CN107679332B (zh) 基于Tekla平台的输电铁塔BIM模型构建方法
Liu et al. A digital twin-based framework for simulation and monitoring analysis of floating wind turbine structures
CN111832185A (zh) 基于蒙皮模型的精确装配仿真方法及系统
CN114626273B (zh) 一种六面体网格划分的螺栓及螺母有限元参数化建模方法
CN110414142B (zh) 浓密机的参数化建模方法
CN112502527B (zh) 一种杆塔原位升高方法
Gu et al. Equivalent static wind loads for stability design of large span roof structures
CN107944100A (zh) 综合生产脐带缆的线型设计方法
CN109313819A (zh) 线路模型实现方法、装置和计算机可读存储介质
CN112507445A (zh) 一种拓扑优化结果的几何结构重构方法
CN116305415A (zh) 一种基于3DExperience的铰接柱脚参数化建模方法
CN116305396A (zh) 基于三维计算机视觉的复杂异形构件节点有限元建模方法
CN113761617A (zh) 矩形风管建筑信息模型的建模方法
CN113722950B (zh) 一种风力机t型法兰锚栓疲劳强度校核方法

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