一种生成30的制作方法

一种生成30
°
楔形防松螺纹精密有限元模型的方法
技术领域
1.本发明属于螺纹有限元模型技术领域，具体涉及一种生成30
°
楔形防松螺纹精密有限元模型的方法。


背景技术：

2.目前，在螺栓精密有限元研究领域，国内外的研究基本集中于对标准米制螺纹的研究。由于有限元模型技术限制且国内外文献中暂无30
°
楔形防松螺纹轮廓公式的相关研究。轮廓公式的缺失导致30
°
楔形防松螺纹的精密有限元计算暂无人研究。


技术实现要素：

3.本发明的目的是为了解决30
°
楔形防松螺纹无相应轮廓公式导致精密有限元模型无法生成的问题，提出了一种生成30
°
楔形防松螺纹精密有限元模型的方法。
4.本发明的技术方案是：一种生成30
°
楔形防松螺纹精密有限元模型的方法包括以下步骤：
5.s1：获取螺母物理参量、轮廓分段点数据和螺纹有限元网格节点三维坐标信息，并根据待生成30
°
楔形防松螺纹轮廓确定螺母牙型的基本结构参数；
6.s2：根据螺母物理参量和螺纹有限元网格节点三维坐标信息计算螺纹轴向节点最小间距与螺纹周向节点最小间距，并根据螺纹轴向节点最小间距与螺纹周向节点最小间距确定容差值；
7.s3：获取螺母有限元网格节点，并对螺母有限元网格节点进行坐标转换；
8.s4：根据坐标转换后的螺母有限元网格节点进行轮廓偏移操作，根据螺母牙型的基本结构参数生成30
°
楔形防松螺纹有限元模型，并确定容差值。
9.进一步地，步骤s1中，螺母物理参量包括内螺纹公称直径d、截距p、螺距等分数n
p
、周向等分数nr、密集网格层数nm和螺纹旋向；
10.轮廓分段点数据包括第一轮廓分段点θ1、第二轮廓分段点θ2、第三轮廓分段点θ3和第四轮廓分段点θ4；其中，
[0011][0012]
螺母牙型的基本结构参数包括螺纹原始三角形高度h和内螺纹基本小径d1。
[0013]
进一步地，步骤s2中，螺纹轴向节点最小间距am的计算公式为：
[0014][0015]
其中，hm表示密集网格层高，p表示螺纹截距，n
p
表示螺距等分数；
[0016]
步骤s2中，将笛卡尔坐标系下编号为1的螺纹有限元网格节点的极坐标(r
p1
,θ
p1
)转换为直角坐标(x
p1
,y
p1
)，将笛卡尔坐标系下编号为2的螺纹有限元网格节点的极坐标(r
p2
,θ
p2
)转换为直角坐标(x
p2
,y
p2
)，并根据直角坐标计算螺纹周向节点最小间距cm，其计算
公式为：
[0017][0018]
其中，θ
p1
＝0x
p1
表示笛卡尔坐标系下编号为1的螺纹有限元网格节点的x轴坐标，y
p1
表示笛卡尔坐标系下编号为1的螺纹有限元网格节点的y轴坐标，x
p2
表示笛卡尔坐标系下编号为2的螺纹有限元网格节点的x轴坐标，y
p2
表示笛卡尔坐标系下编号为2的螺纹有限元网格节点的y轴坐标，r
p1
表示笛卡尔坐标系下编号为1的螺纹有限元网格节点的极径，θ
p1
表示笛卡尔坐标系下编号为1的螺纹有限元网格节点的极角，r
p2
表示笛卡尔坐标系下编号为2的螺纹有限元网格节点的极径，θ
p2
表示笛卡尔坐标系下编号为2的螺纹有限元网格节点的极角，d表示内螺纹公称直径，nr表示周向等分数；
[0019]
步骤s2中，若螺纹轴向节点最小间距am小于螺纹周向节点最小间距cm，则容差值err的计算公式为：
[0020]
err＝errp
×am
[0021]
其中，errp表示容差系数；
[0022]
若螺纹轴向节点最小间距am大于等于螺纹周向节点最小间距cm，则容差值err的计算公式为：
[0023]
err＝errp
×cm
。
[0024]
进一步地，步骤s3中，进行坐标转换的具体方法为：在柱坐标系的θ轴坐标为0-2πrad的范围内，对所有螺母有限元网格节点进行笛卡尔坐标到柱坐标转换。
[0025]
进一步地，步骤s4包括以下子步骤：
[0026]
s41：根据坐标转换后的螺母有限元网格节点，确定多层密集网格节点；
[0027]
s42：对多层密集网格节点进行轮廓偏移操作；
[0028]
s43：设置多层密集网格节点的径向方向插值系数，并更新进行轮廓偏移操作后多层密集网格节点的柱坐标；
[0029]
s44：将更新后多层密集网格节点的柱坐标转换为笛卡尔坐标，并利用转换为笛卡尔坐标的多层密集网格节点替换螺母有限元网格节点，生成30
°
楔形防松螺纹有限元模型并确定容差值。
[0030]
进一步地，步骤s41中，将柱坐标的r轴坐标值在到的螺母有限元网格节点作为多层密集网格节点，其中，d表示内螺纹公称直径，nm表示密集网格层数，hm表示密集网格层高。
[0031]
进一步地，步骤s42中，进行轮廓偏移操作的计算公式为：
[0032][0033]
其中，r'(θ,z)表示在柱坐标系下以θ和z两个坐标轴数据为自变量且r坐标为因变量的函数，p表示螺纹截距，d表示内螺纹公称直径，d1表示内螺纹原始基本小径，h表示螺纹原始三角形高度，θ1表示第一轮廓分段点，θ2表示第二轮廓分段点，θ3表示第三轮廓分段点，θ4表示第四轮廓分段点。
[0034]
进一步地，步骤s43中，多层密集网格节点的径向方向插值系数kr的计算公式为：
[0035][0036]
其中，ri表示编号为i的多层密集网格节点柱坐标系下的r轴坐标值，r
min
表示进行轮廓偏移操作后多层密集网格节点的最内层节点的r轴坐标值，r
max
表示进行轮廓偏移操作后多层密集网格节点的最外层节点的r轴坐标值；
[0037]
步骤s43中，更新后多层密集网格节点的坐标表达式为i”(r
i”,θi,zi)，柱坐标的r轴坐标r
i”的表达式为r
i”＝krri(θi,zi)+(1-kr)r
min
，其中，r
i”表示经过两次变换操作后的r轴坐标，θi表示序号为i的节点的柱坐标系θ轴坐标，zi表示序号为i的节点的柱坐标系z轴坐标。
[0038]
本发明的有益效果是：
[0039]
(1)本发明提出的30
°
楔形防松螺纹轮廓公式可以快速生成30
°
楔形防松螺纹轮廓。该轮廓公式是针对30
°
楔形防松螺纹防松性能有限元研究的必要条件与前提。
[0040]
(2)本发明可以在全自动计算的情况下高效率、高精准度且大批量的获取各种尺寸参数下的30
°
楔形防松螺纹精密有限元模型。将复杂的30
°
楔形防松轮廓网格绘制工作从无法绘制推进到了输入计算轮廓公式后纯算法绘制，为后续30
°
楔形防松轮廓的应力分布、防松性能研究与改进奠定了技术基础。
[0041]
(3)本发明采用的螺纹轮廓网格生成方式为将未变形螺栓模型中的密集网格通过节点偏移的方式变形为螺纹轮廓。节点偏移操作仅对密集网格的形状产生影响，不会对螺栓的任何机械参数产生影响。同时，控制最外层网格不变，最内层网格贴合轮廓后，再对中间层节点进行统一的偏移操作，此算法的优点在于，密集网格层数可以任意设定，无需对算法本身进行修改，适用性强。
[0042]
(4)本发明使用径向插值系数对内部多层网格进行节点坐标更新。此算法的关键点在于，插值系数的选取对控制内部网格的畸变有很大影响，而网格的畸变程度决定了此模型的收敛性。通过改变插值系数的线性与非线性，可以得到更加理想的网格。
附图说明
[0043]
图1为生成30
°
楔形防松螺纹精密有限元模型的流程图；
[0044]
图2为30
°
楔形防松螺母剖面结构示意图；
[0045]
图3为螺母密集网格变换为内螺纹网格示意图。
具体实施方式
[0046]
下面结合附图对本发明的实施例作进一步的说明。
[0047]
如图1所示，本发明提供了一种生成30
°
楔形防松螺纹精密有限元模型的方法，包括以下步骤：
[0048]
s1：获取螺母物理参量、轮廓分段点数据和螺纹有限元网格节点三维坐标信息，并根据待生成30
°
楔形防松螺纹轮廓确定螺母牙型的基本结构参数；
[0049]
s2：根据螺母物理参量和螺纹有限元网格节点三维坐标信息计算螺纹轴向节点最小间距与螺纹周向节点最小间距，并根据螺纹轴向节点最小间距与螺纹周向节点最小间距确定容差值；
[0050]
s3：获取螺母有限元网格节点，并对螺母有限元网格节点进行坐标转换；
[0051]
s4：根据坐标转换后的螺母有限元网格节点进行轮廓偏移操作，根据螺母牙型的基本结构参数生成30
°
楔形防松螺纹有限元模型，并确定容差值。
[0052]
在描述计算机程序流程前，首先需要提出30
°
楔形防松螺纹的轮廓公式。相比普通内螺纹牙底为圆弧形结构，30
°
楔形防松螺纹的牙底结构为30
°
楔形，如图2所示。30
°
楔形防松螺纹轮廓的分段点与普通螺纹基本相同，仅第三轮廓分段点θ3略有不同。
[0053]
在本发明实施例中，步骤s1中，螺母物理参量包括内螺纹公称直径d、截距p、螺距等分数n
p
、周向等分数nr、密集网格层数nm和螺纹旋向；
[0054]
轮廓分段点数据包括第一轮廓分段点θ1、第二轮廓分段点θ2、第三轮廓分段点θ3和第四轮廓分段点θ4；其中，
[0055][0056]
螺母牙型的基本结构参数包括螺纹原始三角形高度h和内螺纹基本小径d1。
[0057]
在本发明实施例中，步骤s2中，螺纹轴向节点最小间距am的计算公式为：
[0058][0059]
其中，hm表示密集网格层高，p表示螺纹截距，n
p
表示螺距等分数；
[0060]
步骤s2中，将笛卡尔坐标系下编号为1的螺纹有限元网格节点的极坐标(r
p1
,θ
p1
)转换为直角坐标(x
p1
,y
p1
)，将笛卡尔坐标系下编号为2的螺纹有限元网格节点的极坐标(r
p2
,θ
p2
)转换为直角坐标(x
p2
,y
p2
)，并根据直角坐标计算螺纹周向节点最小间距cm，其计算公式为：
[0061][0062]
其中，θ
p1
＝0，x
p1
表示笛卡尔坐标系下编号为1的螺纹
有限元网格节点的x轴坐标，y
p1
表示笛卡尔坐标系下编号为1的螺纹有限元网格节点的y轴坐标，x
p2
表示笛卡尔坐标系下编号为2的螺纹有限元网格节点的x轴坐标，y
p2
表示笛卡尔坐标系下编号为2的螺纹有限元网格节点的y轴坐标，r
p1
表示笛卡尔坐标系下编号为1的螺纹有限元网格节点的极径，θ
p1
表示笛卡尔坐标系下编号为1的螺纹有限元网格节点的极角，r
p2
表示笛卡尔坐标系下编号为2的螺纹有限元网格节点的极径，θ
p2
表示笛卡尔坐标系下编号为2的螺纹有限元网格节点的极角，d表示内螺纹公称直径，nr表示周向等分数；
[0063]
步骤s2中，若螺纹轴向节点最小间距am小于螺纹周向节点最小间距cm，则容差值err的计算公式为：
[0064]
err＝errp
×am
[0065]
其中，errp表示容差系数；
[0066]
若螺纹轴向节点最小间距am大于等于螺纹周向节点最小间距cm，则容差值err的计算公式为：
[0067]
err＝errp
×cm
。
[0068]
在本发明实施例中，步骤s3中，进行坐标转换的具体方法为：在柱坐标系的θ轴坐标为0-2πrad的范围内，对所有螺母有限元网格节点进行笛卡尔坐标到柱坐标转换。
[0069]
在本发明实施例中，步骤s4包括以下子步骤：
[0070]
s41：根据坐标转换后的螺母有限元网格节点，确定多层密集网格节点；
[0071]
s42：对多层密集网格节点进行轮廓偏移操作；
[0072]
s43：设置多层密集网格节点的径向方向插值系数，并更新进行轮廓偏移操作后多层密集网格节点的柱坐标；
[0073]
s44：将更新后多层密集网格节点的柱坐标转换为笛卡尔坐标，并利用转换为笛卡尔坐标的多层密集网格节点替换螺母有限元网格节点，生成30
°
楔形防松螺纹有限元模型并确定容差值。
[0074]
在本发明实施例中，如图3所示，步骤s41中，将柱坐标的r轴坐标值在到的螺母有限元网格节点作为多层密集网格节点，其中，d表示内螺纹公称直径，nm表示密集网格层数，hm表示密集网格层高。
[0075]
在本发明实施例中，步骤s42中，进行轮廓偏移操作的计算公式为：
[0076][0077]
其中，r'(θ,z)表示在柱坐标系下以θ和z两个坐标轴数据为自变量且r坐标为因变
量的函数，r'的角标含义为经过单次变换操作后的r坐标，θ表示柱坐标系下的θ坐标轴自变量的数据，p表示螺纹截距，z表示柱坐标系下的z坐标轴自变量的数据，d表示内螺纹公称直径，d1表示内螺纹原始基本小径，h表示螺纹原始三角形高度，θ1表示第一轮廓分段点，θ2表示第二轮廓分段点，θ3表示第三轮廓分段点，θ4表示第四轮廓分段点。
[0078]
在螺纹旋向为右旋时，角度自变量为θ，当螺纹旋向为左旋时，需要将θ替换为2π-θ。r'(θ,z)在周向和轴向上均具有周期性：
[0079]
r'(θ,z)＝r'(θ+2mπ,z+np)m,n＝1,2,3,
…
[0080]
在本发明实施例中，步骤s43中，多层密集网格节点的径向方向插值系数kr的计算公式为：
[0081][0082]
其中，ri表示编号为i的多层密集网格节点柱坐标系下的r轴坐标值，r
min
表示进行轮廓偏移操作后多层密集网格节点的最内层节点的r轴坐标值，r
max
表示进行轮廓偏移操作后多层密集网格节点的最外层节点的r轴坐标值；
[0083]
步骤s43中，更新后多层密集网格节点的坐标表达式为i”(r
i”,θi,zi)，柱坐标的r轴坐标r
i”的表达式为r
i”＝krri(θi,zi)+(1-kr)r
min
，其中，r
i”表示经过两次变换操作后的r轴坐标，θi表示序号为i的节点的柱坐标系θ轴坐标，zi表示序号为i的节点的柱坐标系z轴坐标。
[0084]
本发明的有益效果为：
[0085]
(1)本发明提出的30
°
楔形防松螺纹轮廓公式可以快速生成30
°
楔形防松螺纹轮廓。该轮廓公式是针对30
°
楔形防松螺纹防松性能有限元研究的必要条件与前提。
[0086]
(2)本发明可以在全自动计算的情况下高效率、高精准度且大批量的获取各种尺寸参数下的30
°
楔形防松螺纹精密有限元模型。将复杂的30
°
楔形防松轮廓网格绘制工作从无法绘制推进到了输入计算轮廓公式后纯算法绘制，为后续30
°
楔形防松轮廓的应力分布、防松性能研究与改进奠定了技术基础。
[0087]
(3)本发明采用的螺纹轮廓网格生成方式为将未变形螺栓模型中的密集网格通过节点偏移的方式变形为螺纹轮廓。节点偏移操作仅对密集网格的形状产生影响，不会对螺栓的任何机械参数产生影响。同时，控制最外层网格不变，最内层网格贴合轮廓后，再对中间层节点进行统一的偏移操作，此算法的优点在于，密集网格层数可以任意设定，无需对算法本身进行修改，适用性强。
[0088]
(4)本发明使用径向插值系数对内部多层网格进行节点坐标更新。此算法的关键点在于，插值系数的选取对控制内部网格的畸变有很大影响，而网格的畸变程度决定了此模型的收敛性。通过改变插值系数的线性与非线性，可以得到更加理想的网格。
[0089]
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。