本发明公开了一种回转体截面与质量属性快速精确计算方法，该方法可根据回转体的截面图和材料密度等信息快速准确计算回转体的截面面积、质量、重心轴向位置、极转动惯量、直径转动惯量等属性参数。该方法将预处理过的回转体截面的CAD图纸导入后，以逆时针方向依次提取截面边界上各线段的起点，获得截面边界节点坐标集合，进而根据推导所得的公式计算各属性参数。本发明可以将回转体的二维面积分和三维体积积分转换至一维线积分，避免回转体形状不规则等带来的积分困难，计算更快速；而且只需要截面边界信息，避免了像素化剖分的精度损失，计算更精确。
1.一种回转体截面与质量属性快速精确计算方法，其特征在于，包括以下步骤：步骤1：对回转体截面CAD图纸进行预处理，生成DXF格式文件；
步骤2：导入并解析DXF格式的CAD图纸文件，获得截面图中线段、多线段、圆或弧、样条曲线和圆角的数据；
步骤3：将步骤2得到的多线段、圆或弧、样条曲线和圆角均转换为线段，并和步骤2得到的线段构建线段矩阵L1：其中，xi1、yi1分别代表线段i起点的x、y坐标，xi2、yi2分别代表线段i终点的x、y坐标，其中1≤i≤n；具体方法如下：步骤3-1：将步骤2得到的多线段转换为线段；
步骤3-2：利用De Boor's算法求解指定位置处的样条曲线数值S(x)＝∑iciBi,p(x)，其中，Bi,p(x)代表B样条曲线方程，ci为向量值常数，代表控制点，将步骤2得到的样条曲线转换为线段；
步骤3-4：删除步骤2得到的圆角，连接圆角首末点构成线段；
步骤4：将线段矩阵L1中的线段按照国际单位制缩放，得到线段矩阵L2：L2＝αL1
其中，α代表缩放系数，如果DXF格式的CAD图纸文件设置的长度单位为“mm”，α为0.001；
步骤5：提取截面边界线上各线段起点坐标，即为截面边界节点坐标集合；
步骤6：将截面边界节点坐标集合、回转体材料密度代入计算公式，即得到对应的截面与质量属性参数；具体方法如下：步骤6-1：格林公式及其性质推导
等式左边为面积分，P(x,y)和Q(x,y)为两个在R区域上连续可微函数；等式右边为曲线积分，b为R区域的边界线；
令b(t)＝(x(t),y(t)),t∈[t1,t2]为曲线，f(x,y)为一连续函数；则将以下两个曲线积分变换为定积分形式：曲线积分中有以下性质：
性质1：令b1(t),t∈[t1,t0]和b2(t),t∈[t0,t2]为两个曲线，b1(t0)＝b2(t0)，并且b(t)＝b1∪b2，那么b上的曲线积分由下式计算性质2：若曲线反向，如b’为b的反向曲线，那么积分结果符号改变步骤6-2：封闭曲线的应用分析；
对于封闭图形p，为了运用格林公式以及性质，首先在其边界布局节点，pi＝(xi,yi)，i∈{0,…n},且p0＝pn，从而将边界分段：边界曲率大的地方节点间距小，反之节点间距大，保证两个相邻节点之间的插值函数看作直线，那么整个图形就看作多边形；边界b分为n段，则式中，bi(t),t∈[0,1]定义如下bi(t)＝pi-1+t(pi-pi-1)那么，节点坐标之间存在如下关系
xi(t)＝xi-1+t(xi-xi-1)yi(t)＝yi-1+t(yi-yi-1)xi'(t)＝xi-xi-1
则由性质1得任意连续可微函数沿多边形p＝b(t)的曲线积分由下式计算规定节点在边界上逆时针方向排列，此时面积以及阶矩计算结果为正，若节点排列方向为顺时针，根据性质2知计算结果会与原结果正负相反；
步骤6-3：回转体截面与质量属性参数求解，计算结果均采用国际单位制；
根据回转体特性，在其截面上取微元(dx,dy)则对应整体上的一个以微元为截面、以y为半径的圆环，其体积与质量分别为dv＝2πydxdy
所述计算公式由格林公式推导得出，是材料密度与截面边界节点坐标集合的函数，具体如下：截面面积：
2.根据权利要求1中所述的回转体截面与质量属性快速精确计算方法，其特征在于，步骤1中回转体截面CAD图纸的预处理方法如下：使回转体轴线与x坐标轴重合，图纸比例缩放至1:1，去除回转体截面边界线以外的包括剖面线、尺寸、标注以及轴线，删除截面在第三、四象限内的部分。
3.根据权利要求1中所述的回转体截面与质量属性快速精确计算方法，其特征在于，步骤5中截面边界节点坐标集合的提取方法如下：调整线段矩阵L2中的线段的顺序，使线段编号增加的方向对应截面边界线逆时针方向，从L2中依次提取所有线段起始端点的x、y坐标至对应节点坐标集合(xi,yi)中，其中i＝1,
2,…,n，并增加坐标(x0,y0)，使x0＝xn,y0＝yn。
一种回转体截面与质量属性快速精确计算方法
技术领域
[0001] 本发明属于数字化机械设计技术领域，涉及一种回转体截面与质量属性快速精确计算方法。
背景技术
[0002] 机床主轴、航空发动机等机械产品的转子动力学计算是其设计分析过程中的关键环节。而在计算中，需要输入各个回转体零件的质量、重心轴向位置、极转动惯量、直径转动惯量等属性参数。目前主要CAD软件积分的方法来求解这些参数，即首先将回转体划分成微小规则形状，然后对微小形状分别进行计算后累积叠加，这种方法往往涉及到二维面积分和三维体积分的求解计算，造成形状复杂的零件运算量巨大，求解速度慢；而且由于像素剖分的原因，会使边界处微小形状的计算精度出现损失，求解结果不够精确。
发明内容
[0003] 本发明的目的在于解决回转体截面与质量属性求解过程中，求解速度慢以及求解精度损失的问题，提供一种回转体截面与质量属性快速精确计算方法，
[0004] 为达到上述目的，本发明采用以下技术方案予以实现：
[0005] 一种回转体截面与质量属性快速精确计算方法，包括以下步骤：
[0006] 步骤1：对回转体截面CAD图纸进行预处理，生成DXF格式文件；
[0007] 步骤2：导入并解析DXF格式的CAD图纸文件，获得截面图中线段、多线段、圆或弧、样条曲线和圆角的数据；
[0008] 步骤3：将步骤2得到的多线段、圆或弧、样条曲线和圆角均转换为线段，并和步骤2得到的线段构建线段矩阵L1：
[0009]
[0010] 其中，xi1、yi1分别代表线段i起点的x、y坐标，xi2、yi2分别代表线段i终点的x、y坐标，其中1≤i≤n；
[0011] 步骤4：将线段矩阵L1中的线段按照国际单位制缩放，得到线段矩阵L2：
[0012] L2＝αL1
[0013] 其中，α代表缩放系数，如果DXF格式的CAD图纸文件设置的长度单位为“mm”，α为0.001；如果长度单位为“m”，α为1；
[0014] 步骤5：提取截面边界线上各线段起点坐标，即为截面边界节点坐标集合；
[0015] 步骤6：将截面边界节点坐标集合、回转体材料密度代入计算公式，即得到对应的截面与质量属性参数。
[0016] 本发明进一步的改进在于：
[0017] 步骤1中回转体截面CAD图纸的预处理方法如下：
[0018] 使回转体轴线与x坐标轴重合，图纸比例缩放至1:1，去除回转体截面边界线以外的包括剖面线、尺寸、标注以及轴线，删除截面在第三、四象限内的部分。
[0019] 步骤3中构建线段矩阵L1具体方法如下：
[0020] 步骤3-1：将步骤2得到的多线段转换为线段；
[0021] 步骤3-2：利用De Boor's算法求解指定位置处的样条曲线数值S(x)＝∑iciBi,p(x)，其中，Bi,p(x)代表B样条曲线方程，ci为向量值常数，代表控制点，将步骤2得到的样条曲线转换为线段；
[0022] 步骤3-3：将步骤2得到的圆或弧转换为线段；
[0023] 步骤3-4：删除步骤2得到的圆角，连接圆角首末点构成线段。
[0024] 步骤5中截面边界节点坐标集合的提取方法如下：
[0025] 调整线段矩阵L2中的线段的顺序，使线段编号增加的方向对应截面边界线逆时针方向，从L2中依次提取所有线段起始端点的x、y坐标至对应节点坐标集合(xi,yi)中，其中i＝1,2,…,n，并增加坐标(x0,y0)，使x0＝xn,y0＝yn。
[0026] 步骤6中的计算公式由格林公式推导得出，是材料密度与截面边界节点坐标集合的函数，具体如下：
[0027] 截面面积：
[0028]
[0029] 质量：
[0030]
[0031] 重心轴向位置：
[0032]
[0033] 极转动惯量：
[0034]
[0035] 直径转动惯量：
[0036] Xi＝xi-xc
[0037]
[0038] 步骤6具体如下：
[0039] 步骤6-1：格林公式及其性质推导
[0040] 格林公式为
[0041]
[0042] 等式左边为面积分，P(x,y)和Q(x,y)为两个在R区域上连续可微函数；等式右边为曲线积分，b为R区域的边界线；
[0043] 进一步，等式右侧化为
[0044]
[0045] 令b(t)＝(x(t),y(t)),t∈[t1,t2]为曲线，f(x,y)为一连续函数；则将以下两个曲线积分变换为定积分形式：
[0046]
[0047]
[0048] 曲线积分中有以下性质：
[0049] 性质1：令b1(t),t∈[t1,t0]和b2(t),t∈[t0,t2]为两个曲线，b1(t0)＝b2(t0)，并且b(t)＝b1∪b2，那么b上的曲线积分由下式计算
[0050]
[0051] 性质2：若曲线反向，如b’为b的反向曲线，那么积分结果符号改变[0052]
[0053] 步骤6-2：封闭曲线的应用分析；
[0054] 对于封闭图形p，为了运用格林公式以及性质，首先在其边界布局节点，pi＝(xi,yi)，i∈{0,…n},且p0＝pn，从而将边界分段：边界曲率大的地方节点间距小，反之节点间距大，保证两个相邻节点之间的插值函数看作直线，那么整个图形就看作多边形；边界b分为n段，则
[0055]
[0056] 式中，bi(t),t∈[0,1]定义如下
[0057] bi(t)＝pi-1+t(pi-pi-1)
[0058] 那么，节点坐标之间存在如下关系
[0059] xi(t)＝xi-1+t(xi-xi-1)
[0060] yi(t)＝yi-1+t(yi-yi-1)
[0061] x'i(t)＝xi-xi-1
[0062] y'i(t)＝yi-yi-1
[0063] 则由性质1得任意连续可微函数沿多边形p＝b(t)的曲线积分由下式计算[0064]
[0065] 规定节点在边界上逆时针方向排列，此时面积以及阶矩计算结果为正，若节点排列方向为顺时针，根据性质2知计算结果会与原结果正负相反；
[0066] 步骤6-3：回转体截面与质量属性参数求解，计算结果均采用国际单位制；
[0067] 根据回转体特性，在其截面上取微元(dx,dy)则对应整体上的一个以微元为截面、以y为半径的圆环，其体积与质量分别为
[0068] dv＝2πydxdy
[0069] dm＝ρdv＝ρ2πydxdy
[0070] 具体各个参数的推导过程如下：截面面积：
[0071]
[0072] 质量：
[0073]
[0074] 重心轴向位置：
[0075]
[0076] 极转动惯量：
[0077]
[0078] 直径转动惯量：
[0079]
[0080] 与现有技术相比，本发明具有以下有益效果：
[0081] 本发明将回转体的二维面积分和三维体积积分转换至一维线积分，避免回转体形状不规则等带来的积分困难，计算更快速。本发明计算中只需要截面边界信息，避免了像素化剖分的精度损失，计算更精确。最后，本发明利用格林公式推导将回转体的二维面积分和三维体积积分转换至一维线积分形式，有助于回转体截面与质量属性的快速精准求解，也为机床主轴、航空发动机等转子动力学分析时参数的输入提供了支持。
附图说明
[0082] 图1为回转体截面与质量属性快速精确计算方法的具体流程；
[0083] 图2为回转体截面CAD图的预处理；
[0084] 图3为回转体截面转换为多边形。
具体实施方式
[0085] 下面结合附图对本发明做进一步详细描述：
[0086] 参见图1，本发明回转体截面与质量属性快速精确计算方法，包括以下步骤：
[0087] 步骤1：对回转体截面CAD图纸进行预处理；
[0088] 如图2所示，步骤1的具体方法如下：
[0089] 使回转体轴线与x坐标轴重合(也可与其他轴线重合，本发明中所有公式均参照轴线与x坐标轴重合完成)，图纸比例缩放至1：1，去除回转体截面边界线以外的包括剖面线、尺寸、标注、轴线等内容，删除截面在第三、四象限内的部分。进而生成DXF格式文件；
[0090] 步骤2：导入并解析DXF格式的CAD图纸文件，获得截面图中线段、多线段、样条曲线、圆(弧)和圆角的数据；
[0091] 步骤3：将步骤2得到的多线段、样条曲线、圆(弧)和圆角均转换为线段，并和步骤2得到的线段构建线段矩阵L1：
[0092]
[0093] 其中，xi1、yi1(1≤i≤n)分别代表线段i起点的x、y坐标，xi2、yi2(1≤i≤n)分别代表线段i终点的x、y坐标。如图3所示，此时截面图形转换为多边形。
[0094] 构建线段矩阵L1具体步骤为：
[0095] 步骤3-1：将步骤2得到的多线段转换为线段；
[0096] 步骤3-2：利用De Boor's算法求解指定位置处的样条曲线数值S(x)＝∑iciBi,p(x)，其中，Bi,p(x)代表B样条曲线方程，ci为向量值常数，代表控制点，将步骤2得到的样条曲线转换为线段；
[0097] 步骤3-3：将步骤2得到的圆(弧)转换为线段；
[0098] 步骤3-4：删除步骤2得到的圆角，连接圆角首末点构成线段；
[0099] 步骤4：将线段矩阵L1中的线段按照国际单位制缩放，得到线段矩阵L2：
[0100] L2＝αL1
[0101] 其中，α代表缩放系数，如果DXF格式的CAD图纸文件设置的长度单位为“mm”，α为0.001；如果长度单位为“m”，α为1；
[0102] 步骤5：提取截面边界线上各线段起点坐标，即为截面边界节点坐标集合。
[0103] 具体的提取方法为：调整线段矩阵L2中的线段的顺序，使线段编号增加的方向对应截面边界线逆时针方向，从L2中依次提取所有线段起始端点的x、y坐标至对应节点坐标集合(xi,yi)中，其中i＝1,2,….n，并增加坐标(x0,y0)，使x0＝xn，y0＝yn；
[0104] 步骤6：将截面边界节点坐标集合、回转体材料密度(单位为kg/m3)等代入计算公式，即可快速精确得到对应的截面与质量属性参数。
[0105] 各计算公式由格林公式推导得到，具体推导过程如下：
[0106] 步骤6-1：格林公式及其性质推导
[0107] 格林公式为
[0108]
[0109] 等式左边为面积分，P(x,y)和Q(x,y)为两个在R区域上连续可微函数。等式右边为曲线积分，b为R区域的边界线。
[0110] 进一步，等式右侧可化为
[0111]
[0112] 令b(t)＝(x(t),y(t)),t∈[t1,t2]为曲线，f(x,y)为一连续函数。则可将以下两个曲线积分变换为定积分形式
[0113]
[0114]
[0115] 曲线积分中有以下性质：
[0116] 性质1：令b1(t),t∈[t1,t0]和b2(t),t∈[t0,t2]为两个曲线，b1(t0)＝b2(t0)，并且b(t)＝b1∪b2，那么b上的曲线积分可以由下式计算
[0117]
[0118] 性质2：若曲线反向，如b’为b的反向曲线，那么积分结果符号改变[0119]
[0120] 步骤6-2：封闭曲线的应用分析。
[0121] 对于封闭图形p，为了运用格林公式以及性质，首先在其边界布局节点(pi＝(xi,yi)，i∈{0,…n},且p0＝pn)，从而将边界分段：边界曲率大的地方节点间距小，反之节点间距大，保证两个相邻节点之间的插值函数可以看作直线，那么整个图形就可以看作多边形。边界b分为n段，则
[0122]
[0123] 式中，bi(t),t∈[0,1]定义如下
[0124] bi(t)＝pi-1+t(pi-pi-1)
[0125] 那么，节点坐标之间存在如下关系
[0126] xi(t)＝xi-1+t(xi-xi-1)
[0127] yi(t)＝yi-1+t(yi-yi-1)
[0128] x'i(t)＝xi-xi-1
[0129] y'i(t)＝yi-yi-1
[0130] 则由性质1可得任意连续可微函数沿多边形p＝b(t)的曲线积分可由下式计算[0131]
[0132] 规定节点在边界上逆时针方向排列，此时面积以及阶矩计算结果为正，若节点排列方向为顺时针，根据性质2可知计算结果会与原结果正负相反。
[0133] 步骤6-3：回转体截面与质量属性参数求解，计算结果均采用国际单位制。
[0134] 根据回转体特性，在其截面上取微元(dx,dy)则对应整体上的一个以微元为截面、以y为半径的圆环，其体积与质量分别为
[0135] dv＝2πydxdy
[0136] dm＝ρdv＝ρ2πydxdy
[0137] 具体各个参数的推导过程如下：
[0138] 截面面积：
[0139]
[0140] 质量：
[0141]
[0142] 重心轴向位置：
[0143]
[0144] 极转动惯量：
[0145]
[0146] 直径转动惯量：
[0147]
[0148] 以上内容仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明权利要求书的保护范围之内。
