CN117390788A - 一种水泵水轮机导水机构数字预装最优开档值确定方法 - Google Patents

一种水泵水轮机导水机构数字预装最优开档值确定方法 Download PDF

Info

Publication number
CN117390788A
CN117390788A CN202311372874.5A CN202311372874A CN117390788A CN 117390788 A CN117390788 A CN 117390788A CN 202311372874 A CN202311372874 A CN 202311372874A CN 117390788 A CN117390788 A CN 117390788A
Authority
CN
China
Prior art keywords
value
shift
point
optimal
gear
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.)
Pending
Application number
CN202311372874.5A
Other languages
English (en)
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.)
Huazhong University of Science and Technology
Dongfang Electric Machinery Co Ltd DEC
Original Assignee
Huazhong University of Science and Technology
Dongfang Electric Machinery Co Ltd DEC
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 Huazhong University of Science and Technology, Dongfang Electric Machinery Co Ltd DEC filed Critical Huazhong University of Science and Technology
Priority to CN202311372874.5A priority Critical patent/CN117390788A/zh
Publication of CN117390788A publication Critical patent/CN117390788A/zh
Pending legal-status Critical Current

Links

Classifications

    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/20—Hydro energy

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明公开了一种水泵水轮机导水机构数字预装最优开档值确定方法，涉及水轮机制造技术领域，本发明提供的最优开档值确定方法，用于水泵水轮机导水机构进行数字预装时分析确定导水机构开档值进行调整后的最优值，并用以来判断导水机构装配后开档值是否能满足验收要求。导水机构通过数字预装代替实物预装不仅能同等满足导水机构装配质量的判断和关键数据的验收，还能大幅度减少能耗和压缩装配周期。另外通过逆向测量得到的各零部件的实测数据和模型可用于各零部件加工的综合质量判定及数字孪生产品的实现。本发明可为水泵水轮机导水机构进行数字化预装判定装配质量提供重要的技术基础和关键支撑。

Description

一种水泵水轮机导水机构数字预装最优开档值确定方法
技术领域
本发明涉及水轮机制造技术领域，更具体地说涉及一种水泵水轮机导水机构数字预装最优开档值确定方法。
背景技术
水泵水轮机导水机构各零部件在制造厂内加工完后需进行整体预装以检验各零部件之间的装配尺寸、形位公差等是否满足设计和使用要求。其中，顶盖与底环之间导叶密合平面的开档值是导水机构预装验收的关键数据，它的大小关系着导叶是否能灵活转动的同时还能尽量避免水流泄露损失。目前在现有的制造技术和模式下，导水机构统筹需要先在制造厂内进行整体的实物预装，将顶盖、底环、导叶等所有零部件进行装配。然后逐一测量顶盖底环之间的开档值，通过可调式的工装工具对顶盖或底环的支撑面的各个方向进行调节来寻找导水开档数据的最优值，并以此来判断导水机构预装后的开档值是否满足设计和使用要求。
然而随着先进制造技术的发展，尤其是数字化大尺寸三维测量系统及逆向工程技术的发展，水泵水轮机导水机构可以通过对参与预装的零部件进行精确的逆向三维测量并将逆向建模数据进行装配和计算来实现与实物预装同等的装配质量判断和预装验收，该方法可称为水泵水轮机导水机构的数字预装技术。数字预装技术是水轮机制造业的前沿发展技术，符合国家双碳战略和绿色发展理念。因此，如何利用数字预装技术，对水泵水轮机导水机构的最优开档值进行计算，是亟需解决的问题。
发明内容
为了克服上述现有技术中存在的缺陷，本发明公开了一种水泵水轮机导水机构数字预装最优开档值确定方法，本发明针对水泵水轮机导水机构数字预装技术中导水机构开档值的调整方法与最优开档数据的计算提供一种通用方法，该方法计算得到的最优开档值可到达与实物预装同等的验收结论以满足对导水机构预装质量的准确评价。
为了实现以上目的，本发明采用的技术方案：
一种水泵水轮机导水机构数字预装最优开档值确定方法，包括以下步骤：
一、理论预装三维模型建立
S1、根据顶盖、底环的设计尺寸及导水机构的装配约束和设计开档值，建立理论预装三维模型，并在理论预装三维模型中得到顶盖、底环开档值理论测点的坐标数据；
S1步骤中，根据顶盖、底环的设计尺寸及导水机构的装配约束和设计开档值，建立笛卡尔坐标系下顶盖与底环的标准化理论预装三维模型，依次分别得到顶盖、底环开档值理论测点的坐标数据。
优选的，所述S1步骤包括以下步骤：
S11、分别根据顶盖、底环的设计值建立顶盖、底环的三维模型；在底环的三维模型中，将止漏环圆心设置为坐标原点、止漏环圆心与1#轴孔圆心连线设置为+X轴、密合开档平面朝上设置为+Z轴建立笛卡尔坐标系；并根据导水机构的设计开档值及装配约束关系，将顶盖与底环进行模型装配得到导水机构开档计算的理论预装三维模型；
S12、设导水机构共有N个导叶，以底环密合开档平面为基准面，设1#导叶轴孔与2#导叶轴孔的圆心连线的垂直平分线为L1、导叶中心分布圆线为L2，L1与L2的交点为底环开档值计算的1#理论测点，其坐标为PT1＝(XPT1，YPT1，ZPT1)；然后以坐标原点为基准点，采用环形阵列的方式依次得到底环其余的开档计算理论测点PT2＝(XPT2，YPT2，ZPT2)、PT3＝(XPT3，YPT3，ZPT3)......PTN＝(XPTN，YPTN，ZPTN)；
S13、以顶盖密合开档平面为基准面，设1#导叶轴孔与2#导叶轴孔的圆心连线的垂直平分线为L3、导叶中心分布圆线为L4，L3与L4的交点为顶盖开档值计算的1#理论测点，其坐标为QT1＝(XQT1，YQT1，ZQT1)；然后以坐标原点为基准点，采用环形阵列的方式依次得到顶盖其余的开档计算理论测点QT2＝(XQT2，YQT2，ZQT2)、QT3＝(XQT3，YQT3，ZQT3)......QTN＝(XQTN，YQTN，ZQTN)。
二、过流面三维测量
S2、在生产加工完的实物模型中，对顶盖、底环的过流面进行三维测量获得原始的测量数据，并对测量数据进行预处理并匹配S1步骤的理论预装三维模型，再根据S1步骤理论测点的坐标数据，得到底环、顶盖实测点的实测坐标基础数据；
S2步骤中，对顶盖、底环的过流面进行全面的三维测量并获得原始的测量数据，对测量数据进行一定的计算处理得到可用于下一步优化计算的基础数据。
优选的，所述S2步骤包括以下步骤：
S21、在生产加工完的实物模型中，通过三维触点测量的方式获得底环过流面上关于密合开档平面、止漏环柱面、1#导叶轴孔的实测坐标数据；再采用最小二乘法将密合开档平面的实测坐标拟合成平面N1，将止漏环柱面实测坐标、1#导叶轴孔实测坐标先向N1平面进行投影后再用最小二乘法拟合为圆，得到止漏环投影圆心坐标P1＝(XP1，YP1，ZP1)、1#导叶轴孔投影圆心坐标Q1＝(XQ1，YQ1，ZQ1)；
上述步骤中，密合开档平面、止漏环柱面、1#导叶轴孔的测点数量均≥3。
S22、在生产加工完的实物模型中，通过三维触点测量的方式获得顶盖过流面上关于密合开档平面、止漏环柱面、1#导叶轴孔的实测坐标数据；再采用最小二乘法将密合开档平面的实测坐标拟合成平面N2，将止漏环柱面实测坐标、1#导叶轴孔实测坐标先向N2平面进行投影后再用最小二乘法拟合为圆，得到止漏环投影圆心坐标P2＝(XP2，YP2，ZP2)、1#导叶轴孔投影圆心坐标Q2＝(XQ2，YQ2，ZQ2)；
上述步骤中，密合开档平面、止漏环柱面、1#导叶轴孔的测点数量均≥3。
S23、在S21步骤和S22步骤的基础上，在同一个测量坐标系下采用三维连续扫描测量的方式获得顶盖、底环过流面中关于密合开档平面的点云数据；再通过S21步骤、S22步骤中的N1平面、N2平面、P1、Q1、P2、Q2坐标点，采用点、线、面的匹配方法分别将底环、顶盖三维实测数据与S1步骤中的理论预装三维模型进行匹配，获得实测扫描点云在理论预装三维模型的标准坐标系下的坐标数据；
S24、根据S12步骤和S13步骤中标准坐标系下顶盖及底环的开档值计算理论测点PTi(i＝1、2......N)、QTi(i＝1、2......N)，在三维测量软件中通过理论坐标点的引导分别得到底环、顶盖相应的最近实测点的实测坐标PS1＝(XPS1，YPS2，ZPS3)、PS2＝(XPS2，YPS2，ZPS2)......PSN＝(XPSN，YPSN，ZPSN)，以及QS1＝(XQS1，YQS2，ZQS3)、QS2＝(XQS2，YQS2，ZQS2)......QSN＝(XQSN，YQSN，ZQSN)。
本发明中，需对导水机构顶盖、底环开档调整的最优值进行说明和定义，具体如下：
优选的，顶盖与底环进行预装时导水开档高度调整的最优值为：首先在导水机构开档L有的N个测点中分别计算出对应N个初始开档值，然后在N个测点中寻找某一最优的测点并相应抬高或者降低一最优值使得原有的N个开档值中的最大值与最小值之间的差值最小，该最优值即为导水机构顶盖、底环开档高度调整的最优值。
三、确定开档调整值的比例关系
S3、根据导水机构开档值调整时的几何关系，设置其最优开档值的调整点和调整值，并确定其余各点开档调整值随最优调整值变化的比例关系；
S3步骤中，根据导水开档值调整时的几何关系，假设第i个开档测点是最优开档值的调整点，分析确定其余各点开档值的随之变化的比例关系。
优选的，所述S3步骤包括以下步骤：
S31、设在N个开档测点中最优的调整点为第i点且最优调整值为Hi，当Hi≥0时，该点的开档值抬高|Hi|；当Hi≤0时，该点的开档值降低|Hi|；
优选的，所述S31步骤中，第i点的开档调整量最大，第i点的左右两侧依次相邻的点的调整量逐渐减小，直至到第i点左右两侧完全垂直的两点的调整量为0，然后接着再往前两侧点的调整量值方向与之相反，直至第i点对称点的开档调整值与第i点大小相等、方向相反。
S32、设导叶轴孔分布圆半径为R，则第i点的抬升半径为R，其相邻的第(i-1)的抬升半径R(i-1)为：
Ri-1/R＝COS(2π/N)
分析可知，若第i点的开档调整值为Hi，则第(i-1)的开档调整值H(i-1)与Hi之间的比例关系为：
Hi-1/Hi＝COS(2π/N)
由此类推可得，第(i-m)个点的开档调整值H(i-m)与Hi之间的比例关系为：
Hi-m/Hi＝COS(2π×m/N)
其中，当N是4的整数倍时，m值的上限值为N/4；当N不是4的整数倍时，m值的上限值为(N+2)/4。
上述步骤中，m值的取值上限与导叶轴数N有关，由于导水机构中导叶数量一般为偶数个，当N是4的整数倍时，m值的上限值为N/4；当N不是4的整数倍时，m值的上限值为(N+2)/4。
本发明中，由S31步骤和S32步骤可知，当对导水机构的开档值进行调整时，若第i点为最大调整点，则其余各点的调整值的大小与方向与其在第i点半径方向上的投影半径的大小和方向有关。在第i点附近的1/4圆周测点内，其依次相邻测点开档调整值与第i的调整值之间的比例关系为：
当N是4的整数倍时：
Hi-m/Hi＝COS[(2π×m/N]
其中，i＝1、2......N，m＝1、2......N/4；
当N不是4的整数倍时：
Hi-m/Hi＝COS[(2π×m/N]
其中，i＝1、2......N，m＝1、2......(N+2)/4。
四、穷举计算开档值
S4、利用S2步骤的实测坐标基础数据计算导水机构开档的初始值，并利用初始值得到各对称测点组开档初始值差值以及初始开档调整值；再采用穷举法并结合S3步骤的比例关系，逐一计算每组对称测点作为最优开档调整点时其余开档值并建立开档值矩阵；
S4步骤中，根据理论装配模型提供的标准化坐标系计算导水机构开档的初始值，同时分析得到各对称测点开档差值的差值。根据各对称点初始开档差值确定该点及其对称点的初始开档调整值，然后采用穷举法逐一计算每组对称点作为最优开档调整点时其余开档值的变化情况。
优选的，所述S4步骤包括以下步骤：
S41、在顶盖、底环的初始装配坐标系下，设底环第个i开档测点的开档值为Ki，根据S24步骤中顶盖、底环的实测开档计算坐标数据得到初始状态下导水机构开档值：
S42、计算每个开档测点与其对称点的开档值之间的差值：
ΔKi＝KCi-KC(i+N/2)，(i＝1......2/N)
本发明中，对于任意一个开档测点的开档值，其可调整值的大小是与其对称点的开档值相关的。因此每个开档测点与其对称点的开档值之间的差值可由上式计算得到。
S43、计算由第i个点及第(i+N/2)个点组成的一组对称开档测点的初始开档的最佳调整值：
ΔHi＝ΔKi/2，(i＝1......2/N)
本发明中，由S42步骤得到的差值可知，对于第i个点及第(i+N/2)个点组成的一组对称开档测点来说，其初始开档的最佳调整值为上述ΔHi。
S44、利用S43步骤的最佳调整值，结合S3步骤，得到其余的开档测点相应的调整值；
本发明中，由前述分析可知当导水机构总计有N个开档测点时，其两两对称的开档测点组一共则N/2。当某一组的开档测点最佳调整值由S43步骤计算确定时，其余的开档测点相应的变化值可由前述的S31步骤至S32步骤分析计算得到。
S45、对每一组对称测点的最佳开档调整值对其余开档值的变化情况逐一进行穷举计算并建立开档值矩阵。
本发明中，整个导水机构的开档调整最优值是指在所有的开档测点中寻找出使得最大开档值与最小开档值的差值最小的调整点和调整值。因此若将S44步骤中每一组对称测点的最佳开档调整值对其余开档值的变化情况逐一进行穷举计算，便能寻找出整个导水机构开档调整的最优值。
本发明中，为便于清楚地展示穷举法的计算过程，以第1个开档测点及第1+N/2个开档测点组成的对称开档测点组为例，计算分析整个导水机构的开档值的调整情况，具体如下：
优选的，所述S45步骤的穷举法中，对于第1个开档测点及第1+N/2个开档测点组成的对称开档测点组：
第1个开档测点的调整值为：
第2个开档测点的调整值为：
本发明中，由S32步骤可知，第2个开档测点为第1个开档测点的相邻测点，其开档的调整是与第1个开档测点的调整值成比例的，因此第2个开档测点的调整值为上式。
同理可得当1≤m≤N/4时，第m个开档测点的开档值为：
当开档测点编号超过N/4时，由S31步骤可知其开档测点的调整值的方向将会发生改变，但其调整值的大小与前面的对应测点的调整大小相同。因此：
当m＝N/4+1时：
N/4+1＜m≤N/2时：
当m＝N/2+1时：
当N/2+1＜m≤3N/4时：
当m＝3N/4+1时：
当3N/4＜m≤N时：
由此可知，当导叶轴孔数N为4的整倍数时，开档测点的调整值：
当导叶轴孔数N不为4的整倍数时，开档测点的调整值：
由上述得出的每个开档测点的调整值得到所有开档测点调整后的开档值为：
本发明中，上述方法仅举例说明了在穷举计算中当以第1个开档测点及第1+N/2个开档测点组成的对称开档测点作为最优开档调整值时，整个导水机构开档值的变化情况。在剩余的N/2-1对称开档测点组中，可参考S31步骤的变化规律和上述的计算方法得到其余的开档调整值。
优选的，所述S45步骤中，若干对称开档测点组中，调整后的开档值为：
若干对称开档测点组调整后的开档值组成的矩阵为：
本发明中，由上述公式穷举计算得到的关于导水机构开档调整的所有值，并将其组成以上的矩阵，而整个导水机构开档最优调整值就是从其中每一行表示的调整后开档值中进行比较计算而得。
五、最优值及位置确定
S5、利用S4步骤的开档值矩阵，确定导水机构开档最优调整值、位置以及最优开档值。
S5步骤中，根据前述分析，确定导水机构开档调整的最优值及位置。
优选的，所述S5步骤包括以下步骤：
S51、在S4步骤所有穷举计算的N/2对开档调整组中，计算每一组调整后的所有开档值的最大值与最小值之差εi：
S52、对所有的εi进行比较，找到其最小值所代表的编号并将其设为Z，整个导水机构开档调整后的最优值K′m为：
上述步骤中，对所有的εi(i＝1、2......N/2)进行比较，其中εi的最小值所代表的编号即是整个导水机构开档调整的位置并设其为Z，则整个导水机构开档调整后的最优值K′m可用以上计算式表示。式中KCm、等由前述步骤计算而得。
本发明的有益效果：
本发明提供的最优开档值确定方法，用于水泵水轮机导水机构进行数字预装时分析确定导水机构开档值进行调整后的最优值，并用以来判断导水机构装配后开档值是否能满足验收要求。导水机构通过数字预装代替实物预装不仅能同等满足导水机构装配质量的判断和关键数据的验收，还能大幅度减少能耗和压缩装配周期。另外通过逆向测量得到的各零部件的实测数据和模型可用于各零部件加工的综合质量判定及数字孪生产品的实现。
目前，数字化预装技术在水轮机制造领域还属于探索研究的前沿技术，尚未得到充分的应用和推广。本发明对水泵水轮机导水机构数字预装时关键验收的导水的最优开档值的计算方法进行了分析并形成了规范化的计算方法，可为水泵水轮机导水机构进行数字化预装判定装配质量提供重要的技术基础和关键支撑。
附图说明
图1为本发明底环理论模型示意图；
图2为本发明顶盖底环理论预装模型示意图；
图3为本发明底环三维测量数据示意图；
图4为本发明顶盖三维测量数据示意图；
图5为本发明底环测量数据与理论模型匹配示意图；
图6为本发明顶盖测量数据与理论模型匹配示意图；
图7为本发明开档调整分析示意图；
图8为本发明相邻测点开档调整比例关系；
附图标记：
图1中：1、N#导叶轴孔；2、N#理论开档测点坐标(PTN)；3、1#导叶轴孔；4、1#理论开档测点坐标(PT1)；5、2#导叶轴孔；6、2#理论开档测点坐标(PT2)；7、导叶轴孔中心分布圆(L2)；8、1#、2#导叶轴孔中心垂直平分线(L1)；
图2中：9、底环；10、顶盖；
图3中：11、底环扫描测量点云；12、底环1#轴孔圆心；13、N1平面；14、底环止漏环圆心；
图4中：15、顶盖扫描测量点云；16、顶盖1#轴孔圆心；17、N2平面；18、顶盖止漏环圆心；
图7中：19、第i点开档；20、第i点对称点开档测；
图8中：21、第i个开档测点；22、第i-1个开档测点。
具体实施方式
以下将结合实施例和附图对本发明的构思、具体结构及产生的技术效果进行清楚、完整的描述，以充分地理解本发明的目的、特征和效果。
实施例1
一种水泵水轮机导水机构数字预装最优开档值确定方法，包括以下步骤：
S1、根据顶盖、底环的设计尺寸及导水机构的装配约束和设计开档值，建立理论预装三维模型，并在理论预装三维模型中得到顶盖、底环开档值理论测点的坐标数据；
S2、在生产加工完的实物模型中，对顶盖、底环的过流面进行三维测量获得原始的测量数据，并对测量数据进行预处理并匹配S1步骤的理论预装三维模型，再根据S1步骤理论测点的坐标数据，得到底环、顶盖实测点的实测坐标基础数据；
S3、根据导水机构开档值调整时的几何关系，设置其最优开档值的调整点和调整值，并确定其余各点开档调整值随最优调整值变化的比例关系；
S4、利用S2步骤的实测坐标基础数据计算导水机构开档的初始值，并利用初始值得到各对称测点组开档初始值差值以及初始开档调整值；再采用穷举法并结合S3步骤的比例关系，逐一计算每组对称测点作为最优开档调整点时其余开档值并建立开档值矩阵；
S5、利用S4步骤的开档值矩阵，确定导水机构开档最优调整值、位置以及最优开档值。
实施例2
本实施例在实施例1的基础上对S1步骤作进一步的阐述，S1步骤中，根据顶盖、底环的设计尺寸及导水机构的装配约束和设计开档值，建立笛卡尔坐标系下顶盖与底环的标准化理论预装三维模型，依次分别得到顶盖、底环开档值理论测点的坐标数据。
所述S1步骤包括以下步骤：
S11、分别根据顶盖、底环的设计值建立顶盖、底环的三维模型；在底环的三维模型中，将止漏环圆心设置为坐标原点、止漏环圆心与1#轴孔圆心连线设置为+X轴、密合开档平面朝上设置为+Z轴建立笛卡尔坐标系，如图1所示；并根据导水机构的设计开档值及装配约束关系，将顶盖与底环进行模型装配得到导水机构开档计算的理论预装三维模型，如图2所示；
S12、设导水机构共有N个导叶，以底环密合开档平面为基准面，设1#导叶轴孔与2#导叶轴孔的圆心连线的垂直平分线为L1、导叶中心分布圆线为L2，L1与L2的交点为底环开档值计算的1#理论测点，其坐标为PT1＝(XPT1，YPT1，ZPT1)；然后以坐标原点为基准点，采用环形阵列的方式依次得到底环其余的开档计算理论测点PT2＝(XPT2，YPT2，ZPT2)、PT3＝(XPT3，YPT3，ZPT3)......PTN＝(XPTN，YPTN，ZPTN)，如图1所示；
S13、以顶盖密合开档平面为基准面，设1#导叶轴孔与2#导叶轴孔的圆心连线的垂直平分线为L3、导叶中心分布圆线为L4，L3与L4的交点为顶盖开档值计算的1#理论测点，其坐标为QT1＝(XQT1，YQT1，ZQT1)；然后以坐标原点为基准点，采用环形阵列的方式依次得到顶盖其余的开档计算理论测点QT2＝(XQT2，YQT2，ZQT2)、QT3＝(XQT3，YQT3，ZQT3)......QTN＝(XQTN，YQTN，ZQTN)。
实施例3
本实施例在实施例2的基础上对S2步骤作进一步的阐述，S2步骤中，对顶盖、底环的过流面进行全面的三维测量并获得原始的测量数据，对测量数据进行一定的计算处理得到可用于下一步优化计算的基础数据。
所述S2步骤包括以下步骤：
S21、如图3所示，在生产加工完的实物模型中，通过三维触点测量的方式获得底环过流面上关于密合开档平面、止漏环柱面、1#导叶轴孔的实测坐标数据；再采用最小二乘法将密合开档平面的实测坐标拟合成平面N1，将止漏环柱面实测坐标、1#导叶轴孔实测坐标先向N1平面进行投影后再用最小二乘法拟合为圆，得到止漏环投影圆心坐标P1＝(XP1，YP1，ZP1)、1#导叶轴孔投影圆心坐标Q1＝(XQ1，YQ1，ZQ1)；
上述步骤中，密合开档平面、止漏环柱面、1#导叶轴孔的测点数量均≥3。
S22、如图4所示，在生产加工完的实物模型中，通过三维触点测量的方式获得顶盖过流面上关于密合开档平面、止漏环柱面、1#导叶轴孔的实测坐标数据；再采用最小二乘法将密合开档平面的实测坐标拟合成平面N2，将止漏环柱面实测坐标、1#导叶轴孔实测坐标先向N2平面进行投影后再用最小二乘法拟合为圆，得到止漏环投影圆心坐标P2＝(XP2，YP2，ZP2)、1#导叶轴孔投影圆心坐标Q2＝(XQ2，YQ2，ZQ2)；
上述步骤中，密合开档平面、止漏环柱面、1#导叶轴孔的测点数量均≥3。
S23、在S21步骤和S22步骤的基础上，在同一个测量坐标系下采用三维连续扫描测量的方式获得顶盖、底环过流面中关于密合开档平面的点云数据；再通过S21步骤、S22步骤中的N1平面、N2平面、P1、Q1、P2、Q2坐标点，采用点、线、面的匹配方法分别将底环、顶盖三维实测数据与S1步骤中的理论预装三维模型进行匹配，获得实测扫描点云在理论预装三维模型的标准坐标系下的坐标数据，如图5、6所示；
S24、根据S12步骤和S13步骤中标准坐标系下顶盖及底环的开档值计算理论测点PTi(i＝1、2......N)、QTi(i＝1、2......N)，在三维测量软件中通过理论坐标点的引导分别得到底环、顶盖相应的最近实测点的实测坐标PS1＝(XPS1，YPS2，ZPS3)、PS2＝(XPS2，YPS2，ZPS2)......PSN＝(XPSN，YPSN，ZPSN)，以及QS1＝(XQS1，YQS2，ZQS3)、QS2＝(XQS2，YQS2，ZQS2)......QSN＝(XQSN，YQSN，ZQSN)。
本发明中，需对导水机构顶盖、底环开档调整的最优值进行说明和定义，具体为：首先在导水机构开档已有的N个测点中分别计算出对应N个初始开档值，然后在N个测点中寻找某一最优的测点并相应抬高或者降低一最优值使得原有的N个开档值中的最大值与最小值之间的差值最小，该最优值即为导水机构顶盖、底环开档高度调整的最优值。
实施例4
本实施例在实施例3的基础上对S3步骤作进一步的阐述，S3步骤中，根据导水开档值调整时的几何关系，假设第i个开档测点是最优开档值的调整点，分析确定其余各点开档值的随之变化的比例关系。
所述S3步骤包括以下步骤：
S31、设在N个开档测点中最优的调整点为第i点且最优调整值为Hi，当Hi≥0时，该点的开档值抬高|Hi|；当Hi≤0时，该点的开档值降低|Hi|；
所述S31步骤中，第i点的开档调整量最大，第i点的左右两侧依次相邻的点的调整量逐渐减小，直至到第i点左右两侧完全垂直的两点的调整量为0，然后接着再往前两侧点的调整量值方向与之相反，直至第i点对称点的开档调整值与第i点大小相等、方向相反，如图7所示。
S32、如图8所示，设导叶轴孔分布圆半径为R，则第i点的抬升半径为R，其相邻的第(i-1)的抬升半径R(i-1)为：
Ri-1/R＝COS(2π/N)
分析可知，若第i点的开档调整值为Hi，则第(i-1)的开档调整值H(i-1)与Hi之间的比例关系为：
Hi-1/Hi＝COS(2π/N)
由此类推可得，第(i-m)个点的开档调整值H(i-m)与Hi之间的比例关系为：
Hi-m/Hi＝COS(2π×m/N)
其中，当N是4的整数倍时，m值的上限值为N/4；当N不是4的整数倍时，m值的上限值为(N+2)/4。
上述步骤中，m值的取值上限与导叶轴数N有关，由于导水机构中导叶数量一般为偶数个，当N是4的整数倍时，m值的上限值为N/4；当N不是4的整数倍时，m值的上限值为(N+2)/4。
本发明中，由S31步骤和S32步骤可知，当对导水机构的开档值进行调整时，若第i点为最大调整点，则其余各点的调整值的大小与方向与其在第i点半径方向上的投影半径的大小和方向有关。在第i点附近的1/4圆周测点内，其依次相邻测点开档调整值与第i的调整值之间的比例关系为：
当N是4的整数倍时：
Hi-m/Hi＝COS[(2π×m/N]
其中，i＝1、2......N，m＝1、2......N/4；
当N不是4的整数倍时：
Hi-m/Hi＝COS[(2π×m/N]
其中，i＝1、2......N，m＝1、2......(N+2)/4。
实施例5
本实施例在实施例4的基础上对S4步骤作进一步的阐述，S4步骤中，根据理论装配模型提供的标准化坐标系计算导水机构开档的初始值，同时分析得到各对称测点开档差值的差值。根据各对称点初始开档差值确定该点及其对称点的初始开档调整值，然后采用穷举法逐一计算每组对称点作为最优开档调整点时其余开档值的变化情况。
所述S4步骤包括以下步骤：
S41、在顶盖、底环的初始装配坐标系下，设底环第个i开档测点的开档值为Ki，根据S24步骤中顶盖、底环的实测开档计算坐标数据得到初始状态下导水机构开档值：
S42、计算每个开档测点与其对称点的开档值之间的差值：
ΔKi＝KCi-KC(i+N/2)，(i＝1......2/N)
本发明中，对于任意一个开档测点的开档值，其可调整值的大小是与其对称点的开档值相关的。因此每个开档测点与其对称点的开档值之间的差值可由上式计算得到。
S43、计算由第i个点及第(i+N/2)个点组成的一组对称开档测点的初始开档的最佳调整值：
ΔHi＝ΔKi/2，(i＝1......2/N)
本发明中，由S42步骤得到的差值可知，对于第i个点及第(i+N/2)个点组成的一组对称开档测点来说，其初始开档的最佳调整值为上述ΔHi。
S44、利用S43步骤的最佳调整值，结合S3步骤，得到其余的开档测点相应的调整值；
本发明中，由前述分析可知当导水机构总计有N个开档测点时，其两两对称的开档测点组一共则N/2。当某一组的开档测点最佳调整值由S43步骤计算确定时，其余的开档测点相应的变化值可由前述的S31步骤至S32步骤分析计算得到。
S45、对每一组对称测点的最佳开档调整值对其余开档值的变化情况逐一进行穷举计算并建立开档值矩阵。
本发明中，由S3步骤可知，整个导水机构的开档调整最优值是指在所有的开档测点中寻找出使得最大开档值与最小开档值的差值最小的调整点和调整值。因此若将S44步骤中每一组对称测点的最佳开档调整值对其余开档值的变化情况逐一进行穷举计算，便能寻找出整个导水机构开档调整的最优值。
本发明中，为便于清楚地展示穷举法的计算过程，以第1个开档测点及第1+N/2个开档测点组成的对称开档测点组为例，计算分析整个导水机构的开档值的调整情况，具体如下：
所述S45步骤的穷举法中，对于第1个开档测点及第1+N/2个开档测点组成的对称开档测点组：
第1个开档测点的调整值为：
第2个开档测点的调整值为：
本发明中，由S32步骤可知，第2个开档测点为第1个开档测点的相邻测点，其开档的调整是与第1个开档测点的调整值成比例的，因此第2个开档测点的调整值为上式。
同理可得当1≤m≤N/4时，第m个开档测点的开档值为：
当开档测点编号超过N/4时，由S31步骤可知其开档测点的调整值的方向将会发生改变，但其调整值的大小与前面的对应测点的调整大小相同。因此：
当m＝N/4+1时：
N/4+1＜m≤N/2时：
当m＝N/2+1时：
当N/2+1＜m≤3N/4时：
当m＝3N/4+1时：
当3N/4＜m≤N时：
由此可知，当导叶轴孔数N为4的整倍数时，开档测点的调整值：
当导叶轴孔数N不为4的整倍数时，开档测点的调整值：
由上述得出的每个开档测点的调整值得到所有开档测点调整后的开档值为：
本发明中，上述方法仅举例说明了在穷举计算中当以第1个开档测点及第1+N/2个开档测点组成的对称开档测点作为最优开档调整值时，整个导水机构开档值的变化情况。在剩余的N/2-1对称开档测点组中，可参考S31步骤的变化规律和上述的计算方法得到其余的开档调整值。
所述S45步骤中，若干对称开档测点组中，调整后的开档值为：
若干对称开档测点组调整后的开档值组成的矩阵为：
本发明中，由上述公式穷举计算得到的关于导水机构开档调整的所有值，并将其组成以上的矩阵，而整个导水机构开档最优调整值就是从其中每一行表示的调整后开档值中进行比较计算而得。
实施例6
本实施例在实施例5的基础上对S5步骤作进一步的阐述，S5步骤中，根据前述分析，确定导水机构开档调整的最优值及位置。
所述S5步骤包括以下步骤：
S51、在S4步骤所有穷举计算的N/2对开档调整组中，计算每一组调整后的所有开档值的最大值与最小值之差εi：
S52、对所有的εi进行比较，找到其最小值所代表的编号并将其设为Z，整个导水机构开档调整后的最优值K′m为：
上述步骤中，对所有的εi(i＝1、2......N/2)进行比较，其中εi的最小值所代表的编号即是整个导水机构开档调整的位置并设其为Z，则整个导水机构开档调整后的最优值K′m可用以上计算式表示。式中KCm、等由前述步骤计算而得。
以上对本发明的实施方式进行了具体说明，但本发明并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可作出种种等同变型或替换，这些等同或替换均包含在本发明权利要求所限定的范围内。

Claims (10)

1.一种水泵水轮机导水机构数字预装最优开档值确定方法，其特征在于，包括以下步骤：
S1、根据顶盖、底环的设计尺寸及导水机构的装配约束和设计开档值，建立理论预装三维模型，并在理论预装三维模型中得到顶盖、底环开档值理论测点的坐标数据；
S2、在生产加工完的实物模型中，对顶盖、底环的过流面进行三维测量获得原始的测量数据，并对测量数据进行预处理并匹配S1步骤的理论预装三维模型，再根据S1步骤理论测点的坐标数据，得到底环、顶盖实测点的实测坐标基础数据；
S3、根据导水机构开档值调整时的几何关系，设置其最优开档值的调整点和调整值，并确定其余各点开档调整值随最优调整值变化的比例关系；
S4、利用S2步骤的实测坐标基础数据计算导水机构开档的初始值，并利用初始值得到各对称测点组开档初始值差值以及初始开档调整值；再采用穷举法并结合S3步骤的比例关系，逐一计算每组对称测点作为最优开档调整点时其余开档值并建立开档值矩阵；
S5、利用S4步骤的开档值矩阵，确定导水机构开档最优调整值、位置以及最优开档值。
2.如权利要求1所述的最优开档值确定方法，其特征在于，所述S1步骤包括以下步骤：
S11、分别根据顶盖、底环的设计值建立顶盖、底环的三维模型；在底环的三维模型中，将止漏环圆心设置为坐标原点、止漏环圆心与1#轴孔圆心连线设置为+X轴、密合开档平面朝上设置为+Z轴建立笛卡尔坐标系；并根据导水机构的设计开档值及装配约束关系，将顶盖与底环进行模型装配得到导水机构开档计算的理论预装三维模型；
S12、设导水机构共有N个导叶，以底环密合开档平面为基准面，设1#导叶轴孔与2#导叶轴孔的圆心连线的垂直平分线为L1、导叶中心分布圆线为L2，L1与L2的交点为底环开档值计算的1#理论测点，其坐标为PT1＝(XPT1，YPT1，ZPT1)；然后以坐标原点为基准点，采用环形阵列的方式依次得到底环其余的开档计算理论测点PT2＝(XPT2，YPT2，ZPT2)、PT3＝(XPT3，YPT3，ZPT3)......FTN＝(XPTN，YPTN，ZPTN)；
S13、以顶盖密合开档平面为基准面，设1#导叶轴孔与2#导叶轴孔的圆心连线的垂直平分线为L3、导叶中心分布圆线为L4，L3与L4的交点为顶盖开档值计算的1#理论测点，其坐标为QT1＝(XQT1，YQT1，ZQT1)；然后以坐标原点为基准点，采用环形阵列的方式依次得到顶盖其余的开档计算理论测点QT2＝(XQT2，YQT2，ZQT2)、QT3＝(XQT3，YQT3，ZQT3)......QTN＝(XQTN，YQTN，ZQTN)。
3.如权利要求2所述的最优开档值确定方法，其特征在于，所述S2步骤包括以下步骤：
S21、在生产加工完的实物模型中，通过三维触点测量的方式获得底环过流面上关于密合开档平面、止漏环柱面、1#导叶轴孔的实测坐标数据；再采用最小二乘法将密合开档平面的实测坐标拟合成平面N1，将止漏环柱面实测坐标、1#导叶轴孔实测坐标先向N1平面进行投影后再用最小二乘法拟合为圆，得到止漏环投影圆心坐标P1＝(XP1，YP1，ZP1)、1#导叶轴孔投影圆心坐标Q1＝(XQ1，YQ1，ZQ1)；
S22、在生产加工完的实物模型中，通过三维触点测量的方式获得顶盖过流面上关于密合开档平面、止漏环柱面、1#导叶轴孔的实测坐标数据；再采用最小二乘法将密合开档平面的实测坐标拟合成平面N2，将止漏环柱面实测坐标、1#导叶轴孔实测坐标先向N2平面进行投影后再用最小二乘法拟合为圆，得到止漏环投影圆心坐标P2＝(XP2，YP2，ZP2)、1#导叶轴孔投影圆心坐标Q2＝(XQ2，YQ2，ZQ2)；
S23、在S21步骤和S22步骤的基础上，在同一个测量坐标系下采用三维连续扫描测量的方式获得顶盖、底环过流面中关于密合开档平面的点云数据；再通过S21步骤、S22步骤中的N1平面、N2平面、P1、Q1、P2、Q2坐标点，采用点、线、面的匹配方法分别将底环、顶盖三维实测数据与S1步骤中的理论预装三维模型进行匹配，获得实测扫描点云在理论预装三维模型的标准坐标系下的坐标数据；
S24、根据S12步骤和S13步骤中标准坐标系下顶盖及底环的开档值计算理论测点PTi(i＝1、2......N)、QTi(i＝1、2......N)，在三维测量软件中通过理论坐标点的引导分别得到底环、顶盖相应的最近实测点的实测坐标PS1＝(XPS1，YPS2，ZPS3)、PS2＝(XPS2，YPS2，ZPS2)......PSN＝(XPSN，YPSN，ZPSN)，以及QS1＝(XQS1，YQS2，ZQS3)、QS2＝(XQS2，YQS2，ZQS2)......QSN＝(XQSN，YQSN，ZQSN)。
4.如权利要求3所述的最优开档值确定方法，其特征在于，所述S3步骤包括以下步骤：
S31、设在N个开档测点中最优的调整点为第i点且最优调整值为Hi，当Hi≥0时，该点的开档值抬高|Hi|；当Hi≤0时，该点的开档值降低|Hi|；
S32、确定第(i-m)个点的开档调整值H(i-m)与Hi之间的比例关系为：
Hi-m/Hi＝COS(2π×m/N)
其中，当N是4的整数倍时，m值的上限值为N/4；当N不是4的整数倍时，m值的上限值为(N+2)/4。
5.如权利要求4所述的最优开档值确定方法，其特征在于，所述S31步骤中，第i点的开档调整量最大，第i点的左右两侧依次相邻的点的调整量逐渐减小，直至到第i点左右两侧完全垂直的两点的调整量为0，然后接着再往前两侧点的调整量值方向与之相反，直至第i点对称点的开档调整值与第i点大小相等、方向相反。
6.如权利要求4所述的最优开档值确定方法，其特征在于，所述S4步骤包括以下步骤：
S41、在顶盖、底环的初始装配坐标系下，设底环第个i开档测点的开档值为Ki，根据S24步骤中顶盖、底环的实测开档计算坐标数据得到初始状态下导水机构开档值：
S42、计算每个开档测点与其对称点的开档值之间的差值：
ΔKi＝KCi-KC(i+N/2)，(i＝1......2/N)
S43、计算由第i个点及第(i+N/2)个点组成的一组对称开档测点的初始开档的最佳调整值：
ΔHi＝ΔKi/2，(i＝1......2/N)
S44、利用S43步骤的最佳调整值，结合S3步骤，得到其余的开档测点相应的调整值；
S45、对每一组对称测点的最佳开档调整值对其余开档值的变化情况逐一进行穷举计算并建立开档值矩阵。
7.如权利要求6所述的最优开档值确定方法，其特征在于，所述S45步骤的穷举法中，对于第1个开档测点及第1+N/2个开档测点组成的对称开档测点组：
当导叶轴孔数N为4的整倍数时，开档测点的调整值：
当导叶轴孔数N不为4的整倍数时，开档测点的调整值：
由上述得出的每个开档测点的调整值得到所有开档测点调整后的开档值为：
8.如权利要求7所述的最优开档值确定方法，其特征在于，所述S45步骤中，若干对称开档测点组中，调整后的开档值为：
若干对称开档测点组调整后的开档值组成的矩阵为：
9.如权利要求1所述的最优开档值确定方法，其特征在于，所述S5步骤包括以下步骤：
S51、在S4步骤所有穷举计算的N/2对开档调整组中，计算每一组调整后的所有开档值的最大值与最小值之差εi：
S52、对所有的εi进行比较，找到其最小值所代表的编号并将其设为Z，整个导水机构开档调整后的最优值K′m为：
10.如权利要求1所述的最优开档值确定方法，其特征在于，顶盖与底环进行预装时导水开档高度调整的最优值为：首先在导水机构开档已有的N个测点中分别计算出对应N个初始开档值，然后在N个测点中寻找某一最优的测点并相应抬高或者降低一最优值使得原有的N个开档值中的最大值与最小值之间的差值最小，该最优值即为导水机构顶盖、底环开档高度调整的最优值。
CN202311372874.5A 2023-10-23 2023-10-23 一种水泵水轮机导水机构数字预装最优开档值确定方法 Pending CN117390788A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311372874.5A CN117390788A (zh) 2023-10-23 2023-10-23 一种水泵水轮机导水机构数字预装最优开档值确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311372874.5A CN117390788A (zh) 2023-10-23 2023-10-23 一种水泵水轮机导水机构数字预装最优开档值确定方法

Publications (1)

Publication Number Publication Date
CN117390788A true CN117390788A (zh) 2024-01-12

Family

ID=89440453

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311372874.5A Pending CN117390788A (zh) 2023-10-23 2023-10-23 一种水泵水轮机导水机构数字预装最优开档值确定方法

Country Status (1)

Country Link
CN (1) CN117390788A (zh)

Similar Documents

Publication Publication Date Title
CN106020147B (zh) 一种复杂精密机械产品的装配精度系统分析方法
CN105956234B (zh) 基于逆向工程的汽轮机部套虚拟检测及修配方法
CN113884024B (zh) 一种用于大型风洞收缩段的制作安装检测方法
CN112539721B (zh) 机器人用减速机的三曲柄摆线轮关键加工误差测量方法
CN112395799B (zh) 用于拱坝体形设计合理性及运行期工作状态的评价方法
CN117390788A (zh) 一种水泵水轮机导水机构数字预装最优开档值确定方法
CN112817271A (zh) 基于在机测量的铸造机匣毛坯加工余量优化方法
CN106294940B (zh) 一种空间大型弯曲圆钢管多工序累积误差控制方法
CN107617794A (zh) 一种人字齿内齿圈加工方法
CN110544207A (zh) 大面积微结构形貌的子孔径数据拼接方法、系统及介质
CN1772410A (zh) 锻造复杂零件的飞边流道的几何构造工艺
CN109746495B (zh) 一种四瓣座环单瓣加工找正的工艺方法
CN115857433A (zh) 一种极坐标平面轮廓加工轨迹补偿方法
CN114662334A (zh) 一种水泵水轮机顶盖底环数字预装的最优配对方法
CN111178425B (zh) 一种基于卷积神经网络的大型高速回转装备面型分类方法
CN110956000A (zh) 大口径检验镜重锤支撑参数设计方法
CN112308891A (zh) 一种近净成形类零件多公差约束的配准计算方法
CN112539722A (zh) 机器人用减速机的双曲柄摆线轮关键加工误差测量方法
CN116227241B (zh) 一种锻轴最优加工定位位置确定方法及锻轴加工方法
CN112577454B (zh) 机器人用摆线针轮减速机的针齿壳关键加工误差测量方法
CN221269461U (zh) 大厚板双曲面冷加工一次成型胎具
CN116736796B (zh) 一种钢材尺寸的高精度工控方法
CN110276164B (zh) 环形零部件筋板焊接错位量引起转子结构不平衡量分析法
CN207155330U (zh) 一种手机电池盖加工治具
CN116186935A (zh) 一种轴承装配高极值及相应接触角的确定方法

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