基于结构功能衍生系数的桥式起重机金属结构优化设计
Optimization design of metal structure of bridge crane based on structural function derivative coefficients
收稿日期: 2024-02-19 修回日期: 2024-04-21
基金资助: |
|
First author contact:
Received: 2024-02-19 Revised: 2024-04-21
作者简介 About authors
董青(1989—),女,副教授,博士,从事重大装备数字孪生、疲劳寿命预测、安全评估、报废决策及其相关专业软件开发等研究,E-mail:
关键词:
Keywords:
本文引用格式
董青, 张天祥, 戚其松, 徐格宁.
DONG Qing, ZHANG Tianxiang, QI Qisong, XU Gening.
起重机械是工业、建筑等领域的关键设备,其结构设计对于整机的承载性能和安全运行至关重要。同时,随着实际工程对大吨位、高性能、小质量起重机需求的增长,起重机优化设计已成为行业关注的重点。
目前,国内外学者针对起重机结构设计及轻量化设计的研究已取得一定的成果。如:邢伟等[1]基于有限元理论和拓扑优化方法,通过理论建模、静力学计算和拓扑计算对起重机悬臂梁的质量进行优化,优化后其结构性能和经济性能都能满足工程要求;林伟等[2-4]采取改进鼠群优化等多种算法,实现了起重机主梁的轻量化设计;翁飞翔等[5-7]采取稳健优化等多种优化方法,在确保起重机结构可靠性的同时,显著减小了体积;甘盛霖等[8]的研究聚焦于改善起重机支座的应力分布;陈国雄等[9]结合有限元法与遗传算法,进行了龙门架结构多目标优化设计,取得了较高的求解效率;张磊等[10]通过Inspire和ANSYS软件对起重机伸臂上盖板和腹板的厚度进行尺寸优化;黄镇等[11-17]基于数据集训练预测模型及采用多种优化方法对起重机结构进行了优化,展现了各自方法的有效性和优越性。总体来看,这些研究为起重机结构优化设计提供了多样化的解决方案。
传统的起重机结构优化方法通常侧重于对特定设计参数的调整,而忽略了起重机各零部件之间复杂的相互作用关系。本文通过引入解释结构模型和结构功能衍生系数,提出一种能够准确捕捉起重机各零部件功能与设计参数之间关系的优化设计方法,从而求解最佳的设计参数组合,在确保起重机整机安全性的基础上,实现其金属结构的轻量化设计。
1 起重机金属结构优化设计流程
起重机金属结构设计的合理性对于设备安全运行、提质增效至关重要。为此,本文以桥式起重机(以下简称起重机)为研究对象,以轻量化其金属结构为目标,采用解释结构模型和结构耦合方法,提出基于结构功能衍生系数的起重机金属结构优化设计方法。优化设计的流程如图1所示。
图1
图1
起重机金属结构优化设计流程图
Fig.1
Optimization design flow chart of metal structure of crane
基于结构功能衍生系数的起重机金属结构优化设计流程为:首先,建立起重机的解释结构模型,分析起重机的复杂组成和运作原理,识别出影响力较大的关键结构;其次,针对关键结构进行正交试验和方差分析,建立结构影响因子矩阵,揭示不同零部件之间的耦合关系和影响程度;最后,建立传递函数,求出结构功能衍生系数,确定各个设计参数对起重机结构强度的影响趋势并进行参数优化设计,实现起重机整机的轻量化。
2 解释结构模型
图2
ISM能够全面表示复杂系统的内在结构和机制,通过图形化的方式展示系统元素之间的因果关系,使得分析结果更加直观和易于理解。将系统分解为不同层次,有助于分析系统内部的层级结构和关联关系,清晰地解释系统内各元素之间的影响关系[19]。为了分析起重机零部件之间复杂的依赖关系,确定各元素之间的层次和相互影响规律,将ISM应用于起重机的结构分析上。其具体分析步骤如下:
1)识别元素。识别和定义起重机相关元素,包括但不限于起升机构、大车运行机构、小车运行机构、桥架结构、小车架结构、钢丝绳、滑轮组、传动系统、控制系统和安全装置等,具体如表1所示。
表1 起重机元素
Table 1
符号 | 元素 | 序号 | 元素 | 序号 | 元素 | 序号 | 元素 |
---|---|---|---|---|---|---|---|
S1 | 主梁 | S8 | 小车从动车轮组 | S15 | 副横梁 | S22 | 大车主起升电机 |
S2 | 大车台车梁 | S9 | 副小车架车轮组 | S16 | 主起升定滑轮装配 | S23 | 副起升定滑轮装配 |
S3 | 大车连接梁 | S10 | 小车运行电机 (三合一) | S17 | 主起升动滑轮装配 | S24 | 副起升动滑轮装配 |
S4 | 大车台车 | S18 | 主吊钩组 | S25 | 副吊钩组 | ||
S5 | 大车车轮组 | S11 | 小车台车 | S19 | 主起升卷 筒组 | S26 | 副起升卷筒组 |
S6 | 大车运行电机 (三合一) | S12 | 小车架 | S27 | 副减速器 | ||
S13 | 副小车架 | S20 | 主减速器 | S28 | 副制动器 | ||
S7 | 小车主动车轮组 | S14 | 横梁 | S21 | 主制动器 | S29 | 大车副起升电机 |
2)构建有向图。将元素成对连接来构建有向图,以此表示元素之间的相互依赖关系和影响关系,来分析系统内部结构、组织以及元素之间的相关关系。起重机初始有向图如图3所示。图中每个节点代表一个特定的元素或组件,有向边表示元素或组件之间的直接或间接关系。如:有向边从电动机指向减速器,表示电动机对减速器有直接影响。通过观察这些节点和有向边,来深入了解起重机的组织结构以及元素之间的因果关系,找到对整个系统的性能影响最大的节点,并对其进行优化设计,从而优化起重机的整体性能。
图3
3)建立邻接矩阵。邻接矩阵用于表示有向图的矩阵形式。基于有向图,创建一个用于量化元素之间关联关系的矩阵,来表示每对元素之间的影响是否存在。如:起升电机对小车架和制动器有影响,则在矩阵中的相应位置标记为“1”,否则标记为“0”。起重机模型共有29个零部件,其邻接矩阵是一个29阶矩阵,即
4)建立可达性矩阵。根据上述矩阵建立可达性矩阵,显示关系间的传递,以此显示哪些元素直接或间接地彼此影响。通过对邻接矩阵进行重复的乘法运算,逐渐扩展其路径长度,直到获得最终的可达矩阵,用于描述从一个零部件到另一个零部件是否存在路径,即是否存在一系列有向边来连接这2个零部件。若矩阵的第i行第j列的元素值为1,表示从第i个零部件到第j个零部件存在一条路径;若不存在路径,则元素值为0。可达矩阵同样是一个29阶矩阵,即
式中:b1,1~b29,1表示b1,1,b2,1,b3,1,…,b29,1,其他同。
5)分割层次。根据可达性和影响程度,将元素划分为不同的层次,形成层次结构,其中高层元素对低层元素具有更大的影响力。由
层次分割的步骤为:①建立可达集。可达集是指从一个零部件出发,通过有向边能够到达的所有零部件的集合。若一个零部件的可达集很大,表示它能够直接或间接影响其他许多零部件,可能在系统中具有重要的影响力或控制权。起重机元素的可达集如表2所示。②建立先行集。先行集是指一个零部件所有直接或间接前驱零部件的集合,包括其本身。通过建立先行集,可了解一个零部件对哪些零部件具有依赖关系,即哪些零部件对该零部件的行为、决策具有直接或间接的影响。若一个零部件的先行集很大,表示它依赖于许多其他结构的输入或决策。起重机元素的先行集如表3所示。③根据可达集和先行集,求出可达集与先行集的共同集。④若元素的共同集与可达集相等,则该元素为第1级元素。⑤求出第1级元素后,求出去掉第1级元素后剩余元素的可达集、先行集及共同集。⑥重复上述步骤①至⑤,直至将所有元素划分完。
表2 起重机元素的可达集
Table 2
元素序号 | 可达集 | 元素序号 | 可达集 | 元素序号 | 可达集 |
---|---|---|---|---|---|
S1 | {S1, S2, S3, S11, S12, S13} | S12 | {S1, S2, S3, S11, S12, S13} | S21 | {S1, S2, S3, S11, S12, S13, S14, S16, S17, S18, S19, S20, S21} |
S2 | {S1, S2, S3,S11, S12, S13} | ||||
S3 | {S1, S2, S3, S11, S12, S13} | S13 | {S1, S2, S3, S11, S12, S13} | S22 | {S1, S2, S3, S11, S12, S13, S14, S16, S17, S18, S19, S20, S21, S22} |
S4 | {S1, S2, S3, S4, S11, S12, S13} | S14 | {S1, S2, S3, S11, S12, S13, S14, S16, S17, S18} | S23 | {S1, S2, S3, S11, S12, S13, S15, S23, S24} |
S5 | {S1, S2, S3, S4, S5, S11, S12, S13} | S15 | {S1, S2, S3, S11, S12, S13, S15, S23, S24} | S24 | {S1, S2, S3, S11, S12, S13, S15, S23, S24} |
S6 | {S1, S2, S3, S4, S5, S6, S11, S12, S13} | S16 | {S1, S2, S3, S11, S12, S13, S14, S16, S17, S18} | S25 | {S1, S2, S3, S11, S12, S13, S15, S23, S24, S25} |
S7 | {S1, S2, S3, S7, S11, S12, S13} | S17 | {S1, S2, S3, S11, S12, S13, S14, S16, S17, S18} | S26 | {S1, S2, S3, S11, S12, S13, S15, S23, S24, S26} |
S8 | {S1, S2, S3, S8, S11, S12, S13} | S18 | {S1, S2, S3, S11, S12, S13, S14, S16, S17, S18} | S27 | {S1, S2, S3, S11, S12, S13, S15, S23, S24, S26, S27} |
S9 | {S1, S2, S3, S9, S11, S12, S13, S14, S16, S17, S18} | S19 | {S1, S2, S3, S11, S12, S13, S14, S16, S17, S18, S19} | S28 | {S1, S2, S3, S11, S12, S13, S15, S23,S24, S26, S27, S28} |
S10 | {S1, S2, S3, S7, S10, S11, S12, S13} | S20 | {S1, S2, S3, S11, S12, S13, S14, S16, S17, S18, S19, S20} | S29 | {S1, S2, S3, S11, S12, S13, S15, S23, S24, S26, S27, S28, S29} |
S11 | {S1, S2, S3, S11, S12, S13} |
表3 起重机元素的先行集
Table 3
元素序号 | 先行集 | 元素序号 | 先行集 |
---|---|---|---|
S1 | {S1, S2, S3, S4, S5, S6, S7, S8, S9, S10, S11, S12, S13, S14, S15, S16, S17, S18, S19, S20, S21, S22, S23, S24, S25, S26, S27, S28, S29} | S14 | {S9, S14, S16, S17, S18, S19, S20, S21, S22} |
S15 | {S15, S23, S24, S25, S26, S27, S28, S29} | ||
S2 | {S1, S2, S3, S4, S5, S6, S7, S8, S9, S10, S11, S12, S13, S14, S15, S16, S17, S18, S19, S20, S21, S22, S23, S24, S25, S26, S27, S28, S29} | S16 | {S9, S14, S16, S17, S18, S19, S20, S21, S22} |
S17 | {S9, S14, S16, S17, S18, S19, S20, S21, S22} | ||
S3 | {S1, S2, S3, S4, S5, S6, S7, S8, S9, S10, S11, S12, S13, S14, S15, S16, S17, S18, S19, S20, S21, S22, S23, S24, S25, S26, S27, S28, S29} | S18 | {S9, S14, S16, S17, S18, S19, S20, S21, S22} |
S19 | {S19, S20, S21, S22} | ||
S4 | {S4, S5, S6} | S20 | {S20, S21, S22} |
S5 | {S5, S6} | S21 | {S21, S22} |
S6 | {S6} | S22 | {S22} |
S7 | {S7, S10} | S23 | {S15, S23, S24, S25, S26, S27, S28, S29} |
S8 | {S8} | S24 | {S15, S23, S24, S25, S26, S27, S28, S29} |
S9 | {S9} | S25 | {S25} |
S10 | {S10} | S26 | {S26, S27, S28, S29} |
S11 | {S1, S2, S3, S4, S5, S6, S7, S8, S9, S10, S11, S12, S13, S14, S15, S16, S17, S18, S19, S20, S21, S22, S23, S24, S25, S26, S27, S28, S29} | S27 | {S27, S28, S29} |
S12 | {S1, S2, S3, S4, S5, S6, S7, S8, S9, S10, S11, S12, S13, S14, S15, S16, S17, S18, S19, S20, S21, S22, S23, S24, S25, S26, S27, S28, S29} | S28 | {S28, S29} |
S13 | {S1, S2, S3, S4, S5, S6, S7, S8, S9, S10, S11, S12, S13, S14, S15, S16, S17, S18, S19, S20, S21, S22, S23, S24, S25, S26, S27, S28, S29} | S29 | {S29} |
6)绘制模型。以图表或图形显示元素的层次结构和它们之间的关联关系,使系统结构的依赖关系可视化;在此基础上,建立标注连接关系代号的ISM模型,从而确定带有连接关系的有向图。基于ISM的起重机零部件层次结构及连接关系如图4所示。
图4
图4
基于ISM的起重机零部件层次结构及连接关系
Fig.4
Hierarchy and interconnection relationships of crane parts based on ISM
7)解释和分析。对绘制的ISM模型进行解释和分析,以识别关键元素并评估元素之间的相互影响关系。起重机零部件连接关系如表4所示。
表4 起重机零部件连接关系
Table 4
连接 序号 | 装配连接关系及尺寸 | 符号 | 参数传递关系 | 符号 |
---|---|---|---|---|
S22—S21 | 主起升电机轴直径—联轴器孔径—制动器轴直径 | d1—Φ1—d2 | 主起升电机输出轴转速—联轴器转速 | n1 |
S22—S12 | 电机底座螺纹孔径—小车架固定孔径—小车架上盖板厚 | Φ2—Φ3—x14 | 主起升电机转动产生力—小车架上盖板压力 | F1 |
S29—S28 | 副起升电机轴直径—联轴器孔径—制动器轴直径 | d3—Φ4—d2 | 副起升电机输出轴转速—联轴器转速 | n2 |
S29—S13 | 电机底座螺纹孔径—小车架固定孔—小车架上盖板厚 | Φ5—Φ6—x14 | 副起升电机产生的力—小车架上盖板压力 | F2 |
S21—S20 | 联轴器孔径—减速器轴直径 | Φ1—d4 | 联轴器转速—减速器输入轴转速 | n1 |
S21—S12 | 减速器底座螺纹孔径—小车架固定孔径—小车架上盖板厚 | Φ8—Φ9—x14 | 减速器转动产生的力—小车架上盖板压力 | F3 |
S6—S5 | 大车运行电机轴直径—车轮组孔径 | d5—Φ10 | 大车运行电机输出轴转速—车轮组转速 | n3 |
S28—S27 | 副联轴器孔径—减速器孔径 | Φ11—Φ12 | 副联轴器转速—减速器输入轴转速 | n2 |
S28—S13 | 副制动器底座螺纹孔径—小车架固定孔径 | Φ13—Φ14 | 副减速器转动产生的力—小车架上盖板压力 | F4 |
S20—S19 | 主减速器输出轴直径—卷筒组孔径—卷筒轴直径 | d6—Φ15—d7 | 起升减速器输出轴扭矩—卷筒扭矩 | M1 |
S5—S4 | 大车车轮组螺纹孔直径—台车螺纹孔径—台车主副腹板厚 | Φ16—Φ17—x7 | 大车车轮组支撑力—大车台车压力 | F5 |
S27—S26 | 副减速器输出轴直径—副卷筒组孔径—副卷筒轴直径 | d8—Φ18—d9 | 副起升减速器输出轴扭矩—副卷筒扭矩 | M2 |
S18—S17 | 吊钩组孔径—主起升动滑轮组孔径—销轴直径 | Φ19—Φ20—d10 | 起重量—钢丝绳抗拉强度 | G—σ1 |
S19—S17 | 卷筒组卷筒半径—钢丝绳直经—定滑轮组滑轮半径 | r1—d11—r2 | 钢丝绳上拉力—抗拉强度—抗压强度 | F6 —σ1—σ2 |
S19—S12 | 卷筒组螺纹孔径—小车架螺纹孔径—小车架上盖板厚 | Φ21—Φ22—x14 | 卷筒组压力—小车架上盖板支撑力 | F4 |
S10—S7 | 小车运行电机轴直径—小车运行车轮组孔径 | d12—Φ23 | 小车运行电机输出轴—车轮组转速 | n6 |
S4—S2 | 大车台车连接孔径—台车梁连接孔径—销轴直径 | Φ24—Φ25—d13 | 大车台车支撑力—台车梁压力 | F5 |
S9—S13 | 副小车架车轮组螺纹孔径—副小车架螺纹孔径 | Φ26—Φ27 | 副小车架车轮组支撑力—副小车架压力 | F7 |
S26—S24 | 副卷筒组卷筒半径—钢丝绳直经—定滑轮半径 | r3—d14—r4 | 副起升钢丝绳拉力—抗拉强度—卷筒抗压强度 | F8— σ3—σ4 |
S26—S13 | 副卷筒组螺纹孔径—副小车架螺纹孔径 | Φ28—Φ29 | 副卷筒组压力—副小车架上盖板支撑力 | F9 |
S25—S24 | 副吊钩组孔径—副起升动滑轮组孔径—销轴直径 | Φ30—Φ31—d15 | 副起重量—钢丝绳抗拉强度 | G2—σ3 |
S17—S16 | 主起升动滑轮半径—钢丝绳直径—定滑轮半径 | r2—d11—r5 | 主起升和吊钩组质量—钢丝绳抗拉强度 | G3—σ1 |
S16—S14 | 主起升定滑轮孔径—滑轮轴直径—横梁孔径—横梁主腹板厚—横梁副腹板厚 | Φ32—d16—Φ33—x12—x13 | 主起升和吊钩组滑轮组质量—钢丝绳抗拉强度—横梁抗拉强度 | G4—σ1—σ5 |
S14—S12 | 横梁下盖板厚—小车基距—小车架上盖板厚 | x11—l1—x14 | 横梁压力—小车架支撑力 | F10 |
S8—S11 | 小车从动车轮螺纹孔径—台车螺纹孔径—台车主腹板厚 | Φ34—Φ35—x16 | 小车从动车轮组支撑力—小车台车压力 | F11 |
S7—S1 | 小车主动车轮组半径—大车主梁上盖板厚—轨距 | r6—x1—l2 | 小车轮压/轨道摩擦系数—轨道摩擦力 | P1/u1—f |
S2—S1 | 大车台车梁螺纹孔径—大车主梁螺纹孔径—台车梁上盖板厚—主梁下盖板厚 | Φ36—Φ37— x5—x2 | 大车台车支撑力—主梁上产生的压力 | F12 |
S3—S2 | 大车连接梁孔径—大车台车梁孔径—销轴直径—大车连接梁腹板厚—大车台车腹板厚 | Φ38—Φ39— d17—x7—x9 | 大车台车与连接梁运行过程中产生的推拉力 | F13 |
S24—S23 | 副起升动滑轮半径—钢丝绳直径—定滑轮半径 | r7—d18—r8 | 副起升和吊钩组质量—钢丝绳抗拉强度 | G5—σ3 |
S23—S15 | 副起升定滑轮孔径—滑轮轴直径—副横梁孔径 | Φ40—d19—Φ41 | 副起升量和吊钩组滑轮组质量—钢丝绳抗拉强度—横梁抗拉强度 | G6—σ3—σ6 |
S15—S13 | 副横梁下盖板厚—小车基距—副小车架上盖板厚 | x11—l1—x14 | 副横梁压力—副小车架支撑力 | F14 |
S12—S11 | 小车架孔径—小车台车孔径—销轴直径—小车架主腹板厚—小车台车主腹板厚 | Φ42—Φ43— d20—x16—x18 | 小车架压力—小车台车的支撑力 | F15 |
S11—S1 | 小车台车下盖板厚—主梁上盖板厚—轨距 | x17—x1—l2 | 从动小车轮组轮压/轨道摩擦系数—小车轨道摩擦力 | P2/u1—f |
S13—S1 | 副小车架下盖板厚—主梁上盖板厚—轨距 | x15—x1—l2 | 副小车轮组轮压/轨道摩擦系数—小车轨道摩擦力 | P3/u1—f |
S13—S11 | 副小车架孔径—小车台车孔径销轴直径 | Φ44—Φ45—d21 | 副小车架压力—小车台车的支撑力 | F16 |
起重机的机构主要负责传递和转换能量,使起重机能够进行起吊、运行等动作,而其结构承受设备自身以及货物的重量,并将其传递到地基或其他支撑结构上。在某些特殊工作场合,金属结构的稳定性和抗疲劳性尤为重要。此外,结构设计时需充分考虑起重机的使用环境、工作负荷以及安全要求等因素。若结构设计不当,要么出现起重机无法承受预期的工作负荷,导致使用过程中出现损坏,影响安全运行,要么设计裕量较大,导致资源浪费及材料利用率低。因此,本文将起重机的桥架结构和小车架结构等支撑结构作为下一步研究的主体,其余零部件以表4所示的方式作用在支撑结构上,通过零部件的结构耦合进行支撑结构的综合优化,以此实现起重机的安全性及其绿色性能的提升。
3 结构耦合方法
将结构耦合理论应用到起重机结构分析中,可以更好地分析起重机内部结构之间的相互作用,为起重机的设计、运行和维护提供一个可解释的框架。因此,从结构功能耦合系统、结构影响因子矩阵和结构功能衍生系数等三方面分析起重机结构耦合模型。其中:结构功能耦合系统表示研究策略,即将起重机整体视为一个有着许多交互耦合作用的系统;结构影响因子矩阵表示研究方法,即通过正交试验和方差分析揭示起重机零部件之间的耦合规律;结构功能衍生系数表示分析与计算方法,即根据结构功能耦合系统的传递函数建立衍生系数矩阵及其导数矩阵,反映设计参数变化对起重机结构强度的影响趋势[20]。
3.1 结构影响因子矩阵
针对起重机的桥架结构、小车架结构等支撑结构,建立结构影响因子矩阵,分析各零部件之间的相互作用和影响关系。其基本流程为:1)确定设计参数。确定影响起重机结构性能和系统关键性能的设计参数。由于起重机常常由厂家根据用户需求定制,其具体结构参数如主梁跨度、起升高度、轨距和基距等都固定,故本文重点对结构截面参数进行优化设计。起重机关键结构参数如图5所示。图中:x1, x2, x3, x4分别为主梁上翼缘板厚、下翼缘板厚、主腹板厚和副腹板厚;x5, x6, x7分别为台车梁上翼缘板厚、下翼缘板厚和腹板厚;x8, x9分别为大车连接梁翼缘板厚和腹板厚;x10, x11, x12, x13分别为小车架结构横梁上翼缘板厚、下翼缘板厚、主腹板厚和副腹板厚;x14, x15, x16分别为小车架上翼缘板厚、下翼缘板厚和腹板厚;x17, x18分别为小车走台横梁翼缘板厚和腹板厚。2)设计试验计划。采用正交试验方法,设计覆盖所有设计参数组合的一系列试验,以减少试验次数,同时保证试验的全面性。3)仿真试验。进行正交试验,记录每个试验条件下起重机的结构强度。4)分析数据。采用离差平方和分析方法,计算总离差平方和TSS、组间离差平方和BSS和组内离差平方和WSS,确定试验数据中每个设计参数对系统性能的影响程度。5)建立结构影响因子矩阵,描述不同设计参数之间的相互影响。
图5
式中:
影响因子矩阵 C 为:
式中:V1,V2,…,
3.2 结构功能衍生系数
结构功能衍生系数用于衡量起重机各零部件之间的相互作用及其对整体功能的影响程度,即通过映射函数对上述影响因子进行量化,从而了解其对起重机系统稳定运行的影响程度[21]。若某零部件的结构功能衍生系数较大,表明其对其他零部件的影响较大。因此,应根据结构功能衍生系数的计算结果,重点关注影响力较大的零部件的结构参数。
基于结构功能衍生系数的结构参数优化步骤为:
1)建立结构功能耦合系统的传递函数。针对多因素对起重机结构强度的影响,采用多元线性回归与响应面相结合的方法进行函数拟合。将起重机结构参数xi 视为函数的自变量,结构强度作为输出结果Y。为了避免在参数优化过程中拟合阶数过高而引发振荡问题,选取阶数为2的拟合函数,在保证模型简洁性的同时,可有效捕捉系统的复杂性和非线性特征。则:
式中:b0, b1, …, b3i-1为各项常系数;
2)建立衍生系数矩阵。对传递函数进行偏导数计算,得到与每个起重机结构参数相关的一组多项式
3)建立衍生系数矩阵的导数矩阵。对衍生系数矩阵求导,得到其导数矩阵,以此反映设计参数变化对起重机结构强度影响的变化率。
4)优化设计参数。通过衍生系数及其导数矩阵,分析每个结构参数对起重机性能的影响趋势,以识别关键设计参数,并确定其优化方向。
4 工程实例
以300/100 t-30 m通用桥式起重机为研究对象,利用结构功能衍生系数对起重机金属结构(即桥架结构和小车架结构)的特征参数进行优化,来验证所提出方法的有效性和可行性。该起重机整机的性能参数如下:整机工作级别为A5,跨度为30 m,大车运行速度为4~40 m/min,结构材料为Q235-B;主起升、副起升、电葫芦的起重量分别为300,100,3 t,起升速度分别为0.3~3,3,20 m/min,起升高度分别为16,16.8,18 m;小车起升/运行工作级别为M5,小车运行速度为2~20 m/min,轨距为4 200 mm,主小车基距为4 200 mm,副小车基距为2 650 mm。起重机金属结构关键参数的设计范围如表5所示。
表5 起重机金属结构关键参数的设计范围 (mm)
Table 5
结构 | 参数 | 初始值 | 设计范围 |
---|---|---|---|
桥架结构 | x1 | 24 | 20~28 |
x2 | 24 | 20~28 | |
x3 | 13 | 9~17 | |
x4 | 13 | 9~17 | |
x5 | 34 | 30~38 | |
x6 | 20 | 16~24 | |
x7 | 34 | 30~38 | |
x8 | 12 | 8~16 | |
x9 | 24 | 20~28 | |
小车架结构 | x10 | 16 | 12~20 |
x11 | 20 | 16~24 | |
x12 | 40 | 36~44 | |
x13 | 40 | 36~44 | |
x14 | 36 | 32~40 | |
x15 | 16 | 12~20 | |
x16 | 24 | 20~28 | |
x17 | 5 | 3~7 | |
x18 | 5 | 3~7 |
表6 不同结构参数下起重机金属结构的最大应力
Table 6
结构 | 序号 | x1/ mm | x2/ mm | x3/ mm | x4/ mm | 应力/MPa | 序号 | x5/ mm | x6/ mm | x7/ mm | 应力/MPa | 序号 | x8/ mm | x9/ mm | 应力/MPa | ||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
桥架结构 | 1 | 20 | 20 | 9 | 9 | 142.96 | 10 | 30 | 16 | 30 | 96.41 | 19 | 8 | 20 | 98.586 | ||
2 | 20 | 24 | 13 | 13 | 101.39 | 11 | 30 | 20 | 34 | 98.63 | 20 | 8 | 24 | 98.585 | |||
3 | 20 | 28 | 17 | 17 | 75.96 | 12 | 30 | 24 | 38 | 100.59 | 21 | 8 | 28 | 98.584 | |||
4 | 24 | 20 | 9 | 9 | 135.81 | 13 | 34 | 16 | 30 | 96.44 | 22 | 12 | 20 | 98.578 | |||
5 | 24 | 24 | 13 | 13 | 98.61 | 14 | 34 | 20 | 34 | 98.61 | 23 | 12 | 24 | 98.610 | |||
6 | 24 | 28 | 17 | 17 | 57.13 | 15 | 34 | 24 | 38 | 100.48 | 24 | 12 | 28 | 98.581 | |||
7 | 28 | 20 | 9 | 9 | 139.23 | 16 | 38 | 16 | 30 | 96.32 | 25 | 16 | 20 | 98.595 | |||
8 | 28 | 24 | 13 | 13 | 81.92 | 17 | 38 | 20 | 34 | 98.50 | 26 | 16 | 24 | 98.584 | |||
9 | 28 | 28 | 17 | 17 | 56.50 | 18 | 38 | 24 | 38 | 100.33 | 27 | 16 | 28 | 98.580 | |||
结构 | 序号 | x10/ mm | x11/ mm | x12/ mm | x13/ mm | 应力/MPa | 序号 | x14/ mm | x15/ mm | x16/ mm | 应力/MPa | 序号 | x17/ mm | x18/ mm | 应力/MPa | ||
小车架结构 | 1 | 12 | 16 | 36 | 36 | 147.88 | 10 | 32 | 12 | 20 | 158.10 | 19 | 3 | 3 | 131.08 | ||
2 | 12 | 20 | 40 | 40 | 129.30 | 11 | 32 | 16 | 24 | 141.66 | 20 | 3 | 5 | 130.99 | |||
3 | 12 | 24 | 44 | 44 | 149.95 | 12 | 32 | 20 | 28 | 147.74 | 21 | 3 | 7 | 130.79 | |||
4 | 16 | 16 | 36 | 36 | 146.50 | 13 | 36 | 12 | 20 | 148.43 | 22 | 5 | 3 | 130.91 | |||
5 | 16 | 20 | 40 | 40 | 130.56 | 14 | 36 | 16 | 24 | 130.56 | 23 | 5 | 5 | 130.56 | |||
6 | 16 | 24 | 44 | 44 | 148.96 | 15 | 36 | 20 | 28 | 147.81 | 24 | 5 | 7 | 130.53 | |||
7 | 20 | 16 | 36 | 36 | 145.33 | 16 | 40 | 12 | 20 | 134.68 | 25 | 7 | 3 | 130.85 | |||
8 | 20 | 20 | 40 | 40 | 127.26 | 17 | 40 | 16 | 24 | 120.67 | 26 | 7 | 5 | 130.75 | |||
9 | 20 | 24 | 44 | 44 | 151.95 | 18 | 40 | 20 | 28 | 147.86 | 27 | 7 | 7 | 130.44 |
表7 不同参数水平下起重机金属结构的最大应力 (MPa)
Table 7
结构 | 参数 | 参数水平 | ||
---|---|---|---|---|
V1 | V2 | V3 | ||
桥架结构 | x1 | 101.39 | 98.61 | 81.92 |
x2 | 135.81 | 98.61 | 57.13 | |
x3 | 139.23 | 81.92 | 56.50 | |
x4 | 139.23 | 81.92 | 56.50 | |
x5 | 98.63 | 98.61 | 98.50 | |
x6 | 96.44 | 98.61 | 100.48 | |
x7 | 96.32 | 98.50 | 100.33 | |
x8 | 98.58 | 98.61 | 98.58 | |
x9 | 98.58 | 98.61 | 98.58 | |
小车架结构 | x10 | 129.30 | 130.56 | 120.67 |
x11 | 146.50 | 130.56 | 147.81 | |
x12 | 145.33 | 127.26 | 151.95 | |
x13 | 145.33 | 127.26 | 151.95 | |
x14 | 141.66 | 130.56 | 120.67 | |
x15 | 148.43 | 130.56 | 147.81 | |
x16 | 158.10 | 141.66 | 147.74 | |
x17 | 130.99 | 130.56 | 130.44 | |
x18 | 130.91 | 130.56 | 130.53 |
通过计算不同参数水平组合下起重机金属结构的最大应力,可得桥架结构和小车架结构上述数据的TSS、BSS、WSS分别为:10 695.697 1,3 669.084 2,7 026.6128 和2 832.043 2,557.150 4,2 274.892 8。参数组合时给出了3组参数水平V1、V2和V3,因此所建立的桥架结构和小车架结构的影响因子矩阵 C1和 C2分别如
表8 起重机金属结构最大应力方差分析
Table 8
结构 | 参数 | 离差平方和 | 影响因子 |
---|---|---|---|
桥架结构 | x1 | 241.7 | 0.022 6 |
x2 | 3099.9 | 0.289 8 | |
x3 | 3639.5 | 0.340 1 | |
x4 | 3639.5 | 0.340 1 | |
x5 | 12.4 | 0.001 2 | |
x6 | 19.812 1 | 0.001 9 | |
x7 | 18.152 9 | 0.001 7 | |
x8 | 12.7 | 0.001 2 | |
x9 | 12.5 | 0.001 2 | |
SST | 10 695.70 | ||
小车架结构 | x10 | 372.65 | 0.131 5 |
x11 | 241.92 | 0.0854 | |
x12 | 381.26 | 0.134 6 | |
x13 | 381.26 | 0.134 6 | |
x14 | 338.97 | 0.119 6 | |
x15 | 281.36 | 0.099 3 | |
x16 | 574.44 | 0.203 3 | |
x17 | 130.18 | 0.045 9 | |
x18 | 129.97 | 0.045 8 | |
SST | 2 832.04 |
对于桥架结构,有:
由
对于小车架结构,有:
由
依据函数模型,将起重机的一个结构参数视为变量,其他结构参数视为常量,来求出所有包含该变量的单项式,从而构建关于该变量的多项式,即为结构功能衍生系数。整理结构功能衍生系数,可得到衍生系数矩阵。该矩阵的因子为:
对
在参数值按迭代步长变化时,为了使结构应力有明显的变化趋势,根据各设计参数的取值范围,设定设计参数的优化间隔为1 mm,同时优先选择结构应力影响因子大于20%的参数进行优化设计。
在初始设计参数下,桥架结构的最大应力为98.61 MPa,小车架结构的最大应力为130.56 MPa。在初始参数组合下,可求得桥架结构和小车架结构的结构功能衍生系数矩阵的导数矩阵分别为:
其单调性矩阵
由
在此基础上,采用结构耦合法对起重机金属结构进行优化设计,相关优化结果如表9所示。
表9 起重机金属结构优化结果
Table 9
结构 | 参数 | 第1次优化值/mm | 单调增减情况 | 第2次优化值/mm | 单调增减情况 | 第3次优化值/mm | 单调增减情况 | 第3次优化值/mm | 单调增减情况 |
---|---|---|---|---|---|---|---|---|---|
桥架结构 | x1 | 24 | - | 24 | - | 24 | - | 24 | - |
x2 | 23 | + | 22 | + | 21 | + | 20 | + | |
x3 | 14 | - | 14 | - | 14 | - | 14 | - | |
x4 | 14 | - | 14 | - | 14 | - | 14 | - | |
x5 | 34 | - | 34 | - | 34 | - | 34 | - | |
x6 | 20 | - | 20 | - | 20 | - | 20 | - | |
x7 | 34 | + | 34 | + | 34 | + | 34 | + | |
x8 | 12 | - | 12 | - | 12 | - | 12 | - | |
x9 | 24 | - | 24 | - | 24 | - | 24 | - | |
3.567 5 | 2.489 7 | 1.411 9 | 0.334 1 | ||||||
Y1/MPa | 102.177 5 | 104.667 2 | 106.079 1 | 106.413 2 | |||||
小车架结构 | x10 | 15 | + | 14 | 13 | 12 | |||
x11 | 20 | - | 20 | 20 | 20 | ||||
x12 | 39 | + | 38 | 37 | 36 | ||||
x13 | 39 | + | 38 | 37 | 36 | ||||
x14 | 35 | + | 34 | 33 | 32 | ||||
x15 | 16 | + | 16 | 16 | 16 | ||||
x16 | 24 | - | 24 | 24 | 24 | ||||
x17 | 5 | + | 5 | 5 | 5 | ||||
x18 | 5 | + | 5 | 5 | 5 | ||||
10.863 7 | 9.8621 | 8.860 5 | 7.858 9 | ||||||
Y2/MPa | 141.423 7 | 151.285 8 | 160.146 3 | 168.005 2 |
经过4次寻优后,x2, x10, x12, x13, x14的值分别为20,12,36,36,32 mm,即:优化后主梁下翼缘板厚减小了4 mm,使得桥架结构最大应力增大了7.91%;小车架结构横梁上翼缘板厚、主腹板厚、副腹板厚和小车架上翼缘板厚都减小了4 mm;小车架结构最大应力增大了28.68%,桥架质量和小车架质量分别减小了1 637.388,397.725 kg,起重机整机质量减小了2 035.113 kg。
5 结 论
在起重机结构优化设计中,传统的优化方法忽略了起重机各零部件之间复杂的相互作用。本文提出了基于结构功能衍生系数的起重机金属结构优化设计方法,在确保起重机整机安全性的基础上,实现了金属结构轻量化设计。本文重点如下:
1)建立了起重机零部件解释结构模型,引入了结构功能衍生系数,通过仿真分析,获得了不同设计参数组合下起重机金属结构的最大应力值,并以此确定了可近似描述设计参数与结构强度之间关系的拟合函数。通过衍生系数矩阵及其导数矩阵,量化了设计参数的变化对结构强度的影响。
2)以300/100 t-30 m通用桥式起重机为研究对象,根据设计参数对结构强度的影响趋势,通过参数优化,使得桥架结构的最大应力增大了7.91%,小车架结构的最大应力增大了28.68%,起重机整机质量减小了2 035.113 kg,在确保起重机金属结构承载能力的前提下,大幅度减小了整机质量。
该方法可以针对起重机不同的零部件进行具体分析,但存在一定的局限性,如衍生系数矩阵及其导数矩阵的建立需要大量试验数据。今后需进行更多的研究来突破其局限性,包括研究出更有效的方法来建立衍生系数矩阵及其导数矩阵,或者用更精准的数值模拟技术来减少对试验数据的依赖。
参考文献
某发射塔起重机悬臂梁结构拓扑优化设计
[J].
Topological optimization design of cantilever beam structure for a transmitter tower crane
[J].DOI:10.3969/j.issn.1001-3881.2022.21.022 [本文引用: 1]
基于改进鼠群优化算法的起重机主梁轻量化设计
[J].
Lightweight design of crane’s main girder based on improved rat swarm optimization algorithm
[J].
基于混合GSA-GA的起重机主梁优化设计
[J].
Optimization design of crane main beam based on hybrid GSA-GA
[J].DOI:10.3969/j.issn.1001-3997.2021.10.043
基于篮球联赛优化算法的通用桥式起重机主梁轻量化设计
[J].
Lightweight design of main girder of an universal bridge crane based on basketball league optimization algorithm
[J].
考虑可靠度的桥式起重机主梁金属结构稳健优化设计
[J].
Robust optimization design of metal girder structure of bridge crane considering reliability
[J].
基于稀疏网格模型与MOGA算法的起重机主梁结构优化
[J].
Optimization of the main beam of bridge crane based on sparse grid model and MOGA algorithm
[J].
基于灵敏度分析的门式起重机结构优化设计
[J].
Optimization design of gantry crane structure based on sensitivity analysis
[J].
平头塔式起重机回转总成有限元分析及结构优化
[J].
Finite element analysis and structure optimum of flat-top tower cranes
[J].DOI:10.1201/9781003078548-27 [本文引用: 1]
焊接平台龙门架结构轻量化设计
[J].
Welding platform gantry structure lightweight design
[J].
面向轻量化目标的起重机伸臂结构设计及截面尺寸优化
[J].
Crane jib structure design and section dimension optimization for lightweight target
[J].
起重机主梁截面风力系数预测及结构优化设计
[J].
Wind coefficient prediction and structural optimization design of crane girder section
[J].
Reliability-based design optimization for the lattice boom of crawler crane
[J].
Optimization of the box section of the main girder of the single-girder bridge crane by applying biologically inspired algorithms
[J].
Multi-objective design optimization using hybrid search algorithms with interval uncertainty for thin-walled structures
[J].
Structural reliability-based design optimization with non-probabilistic credibility level
[J].
Multi-objective structural optimization for the automatic member grouping of truss structures using evolutionary algorithms
[J].
Bayesian optimization-based topology optimization using moving morphable bars for flexible structure design problems
[J].
基于矩阵自乘的一种解释结构模型改进方法
[J].
An improved method of interpreting structural model based on matrix self-multiplication
[J].
基于解释结构模型和复杂网络的制造系统关键脆性因子辨识
[J].
Identification of key brittleness factor for manufacturing system based on ISM and complex network
[J].
基于零件结构耦合理论的主动再制造设计
[D].
Active remanufacturing design based on coupling theory of parts structure
[D].
基于结构功能衍生系数的主动再制造设计
[J].
Optimization design of predecisional remanufacturing based on structural function derivative coefficient
[J].DOI:10.3901/jme.2017.11.175 [本文引用: 2]
/
〈 |
|
〉 |
