[发明专利]一种机床模块式刚度表征与求解方法

摘要

一种机床模块式刚度表征与求解方法，先将机床划分成多个子系统模块，然后确定载荷在各子系统模块中的传递方向，再对各子系统模块的承载位置和约束位置进行辨识，然后对每一子系统模块，分别构建并求解可以表征其刚度特性的柔度矩阵以及可以表征其载荷传递特性的载荷传递矩阵，再对任一工作位置状态下的机床整体，根据它各子系统模块的几何尺寸、柔度矩阵、载荷传递矩阵参数以及相互之间的位置关系，求解出机床整体的刚度特性，本发明可帮助设计人员发现机床系统中刚度较为薄弱的环节，便于获得机床在任意工作位置的刚度特性。

主权项

1.一种机床模块式刚度表征与求解方法，其特征在于，包括以下步骤：步骤一：将机床按照几何结构划分成两个以上的子系统模块，每个子系统模块由两个以上的零部件组合而成；步骤二：确定载荷在各子系统模块中的传递方向，选取刀具切削工件的位置为载荷的起点，机床的安装位置，即地面为载荷的终点，并经过刀具一侧和工件一侧分别由起点传递到终点；步骤三：对各子系统模块的承载位置和约束位置进行辨识，其中：承载位置是指子系统模块中靠近步骤二中定义的载荷起点的一端、与上一子系统模块相连接的位置；承载位置会承受上一子系统模块传递来的载荷，承载位置方向和位置的变化会引起与之相连的上一子系统模块方向和位置的变化；约束位置是指子系统模块中靠近步骤二中定义的载荷终点的一端、与下一子系统模块相连接的位置；子系统模块的约束位置会限制承载子系统模块的自由度，并且约束位置的约束反力的反作用力即为该子系统模块传递给下一子系统结构的载荷；机床各子系统模块间的典型结合方式分为固定结合与移动结合，其中，移动结合又分导轨-滑块结合、丝杠-螺母结合，并且这两种移动结合都是一同出现在两个可发生相对移动的子系统模块之间的；若承载位置通过固定结合来联接其他子系统模块，直接选取结合面的方向和位置来表征承载位置的方位；若承载位置通过导轨-滑块结合联接其他子系统模块，选取各导轨-滑块法相结合平面的中心点所确定的平面的方向和位置来表征承载位置的方位；若承载位置通过丝杠-螺母结合来联接其他子系统模块，因丝杠主要承受轴向载荷，只需考虑其轴向的载荷及位移量，而该位移又会使相应的导轨-滑块结合处发生同样大小的沿丝杠轴向的相对位置变化，所以同样选取各导轨-滑块法相结合平面的中心点所确定的平面的方向和位置来表征这类承载位置的方位；若约束位置通过固定结合来联接其他子系统模块，则该处的约束为全约束；若约束位置通过导轨-滑块结合联接其他子系统模块，则该处的约束只保留导轨方向的平动自由度；若约束位置通过丝杠-螺母结合来联接其他子系统模块，则该处的约束为丝杠轴向的位移约束；步骤四：对每一个子系统模块，构建可表征其刚度特性的柔度矩阵C，根据刚度矩阵的定义，刚度矩阵K与该子系统模块的柔度矩阵C互为逆矩阵，即K=[C]^-1；子系统模块的柔度矩阵C表征该子系统模块承载位置的位移量与其承受的载荷之间的关系，该矩阵的具体维数与载荷作用位置的数目n有关，即柔度矩阵C的维数为6×6n，当n=1时具体表达式为：C=c11c12c13c14c15c16c21c22c23c24c25c26c31c32c33c34c35c36c41c42c43c44c45c46c51c52c53c54c55c56c61c62c63c64c65c66]]>矩阵中的元素，即柔度系数，c_ij表示：对所表征子系统模块承载位置施加j方向的单位大小的力或力矩时，该子系统模块的承载位置会发生相对于约束位置大小为c_ij的位移量；当j=1、2、3时，相应的载荷分别为Ｘ、Ｙ、Ｚ三个方向的力，当j＝4、5、6时，相应的载荷分别为Ｘ、Ｙ、Ｚ三个方向的力距；当i=1、2、3时，相应的位移为Ｘ、Ｙ、Ｚ三个方向的平动位移量，当i＝4、5、6时，相应的位移分别为Ｘ、Ｙ、Ｚ三个方向的转动位移量，其中，X、Y、Z坐标按照机床行业内通用的坐标轴划分方法来选取；在弹性变形范围内，机床各子系统模块的变形与其所承受的载荷之间的关系是符合线性叠加原理的，对某一子系统模块，承载位置承受载荷F时的位移量E即可根据式[E]=[C][F]求得；步骤五：构建子系统模块的载荷传递矩阵B，B表征子系统模块约束位置的约束反力的反作用力，即该子系统模块对下一子系统模块的作用力F’，与该子系统模块承载位置所承受的载荷F之间的关系，表达式为[F’]=[B][F]，载荷传递矩阵B的维数视F’和F的维数而定，维数视F’和F的维数即该子系统模块所承受的约束反力与载荷的具体数目，当子系统模块所承受的约束反力与载荷的具体数目均为1时，B的表达式为：B=b11b12b13b14b15b16b21b22b23b24b25b26b31b32b33b34b35b36b41b42b43b44b45b46b51b52b53b54b55b56b61b62b63b64b65b66]]>矩阵中的元素，即载荷传递系数，b_ij表示：对所表征子系统模块承载位置施加j方向的单位大小的力或力矩时，该子系统模块相应的约束位置会产生大小为-b_ij的约束反力，即引起的对下一子系统模块的作用力为b_ij；当j=1、2、3时，相应的载荷分别为Ｘ、Ｙ、Ｚ三个方向的力，当j＝4、5、6时，相应的载荷分别为Ｘ、Ｙ、Ｚ三个方向的力距；当i=1、2、3时，相应的约束反力为Ｘ、Ｙ、Ｚ三个方向的力，当i＝4、5、6时，相应的约束反力分别为Ｘ、Ｙ、Ｚ三个方向的力矩；根据各子系统模块的载荷传递矩阵B，即可求解载荷在整个机床系统中的传递情况；步骤六：子系统模块柔度矩阵和载荷传递矩阵的获取，需要借助有限元商业软件来求解各子系统模块的柔度矩阵以及载荷传递矩阵,具体步骤如下：6-1）在三维计算机辅助设计软件中构建各子系统模块的三维模型；6-2）将步骤6-1）中构建的子系统模块的三维模型以中间文件的格式导入到有限元商业软件中，设置材料属性，按照步骤三中给定的方法，确定相应的承载位置和约束位置，并对模型添加相应的约束；6-3）在有限元商业软件中对相应子系统模块的承载位置分别施加各个方向的单位力与力距F_x、F_y、F_z、M_x、M_y、M_z；运用有限元商业软件的求解器对子系统模块在相应载荷下的有限元模型进行求解；6-4）将步骤6-3）中求解得到的结果文件通过有限元商业软件后处理模块打开之后分别提取各载荷作用下承载位置的形变位移量以及约束位置的约束反力，形变位移量包括平动和转动；6-5）根据步骤6-4）中提取出的结果计算对应的柔度系数，之后将计算出的各柔度系数组合形成表征该部件刚度特性的柔度矩阵，具体表达式如步骤四中所示；6-6）根据步骤6-4）中提取出的结果计算对应的载荷传递系数，之后将计算出的各载荷传递系数组合形成表征该部件载荷传递特性的载荷传递矩阵，具体表达式如步骤五中所示；6-7）对于承载位置不固定的子系统模块，在其所有可能的承载位置中选取离散的若干个位置作为取样，重复步骤6-2）～6-6）得到不同位置下的子系统模块的柔度矩阵和载荷传递矩阵，并根据得到的结果拟合出以承载位置的坐标为自变量的可近似表征任一承载位置下的该子系统模块的柔度矩阵和载荷传递矩阵的函数表达式；步骤七：根据上述步骤，对于处于任一工作位置的机床整机，根据其各子系统模块的间的相对位置关系即可获得各子系统模块的承载位置参数，进而可求得各子系统模块的载荷传递矩阵，通过把两个子系统模块的载荷传递矩阵按照顺序相乘，顺序是靠近步骤二中定义的载荷作用起点的子系统模块的载荷传递矩阵位于乘式的右边，即可得到表征载荷在这两个子系统模块中传递情况的新的载荷传递矩阵，进而可得到载荷在机床整体中的传递情况，即各子系统模块的承载情况；步骤八：求解机床整体的刚度特性，具体步骤如下：8-1）以机床安装面为基准构建整体坐标系，以各子系统的约束位置为基准分别建立各子系统模块的局部坐标系；8-2）根据某一工作位置下机床各子系统模块的承载位置参数确定其柔度矩阵和步骤七中求得的承载情况，求得各子系统模块的承载位置在以其约束位置为基准的局部坐标系中方向和位置的变化情况；8-3）各子系统模块承载位置方位的变化又会引起与之相连的子系统模块的局部坐标系相对于整体坐标系的方位变化，并最终引起载荷的起点，即刀具切削工件位置的方位变化，根据空间直角坐标系的齐次坐标变换的相关理论和机床各子系统模块的几何尺寸关系，即可将任一子系统模块承载位置方位的变化转化为刀具与工件在整体坐标系内方位的变化；8-4）将步骤8-3）中得到的各子系统模块在载荷作用下对刀具与工件在整体坐标系内方位变化的效果相叠加，即可得到综合各子系统模块变形情况和载荷传递情况的整机刚度特性；至此即可得到任一工作位置下的机床整体以及各子系统模块的刚度特性以及载荷传递特性。