基于SIMP方法的爬杆机器人结构优化与分析.pdf

第 30 卷第 3 期2023 年 6 月 工程设计学报 Chinese Journal of Engineering DesignVol.30 No.3Jun.2023基于SIMP方法的爬杆机器人结构优化与分析郑小飞，黄镇海，马小龙，王建新，王斌锐（中国计量大学 机电工程学院，浙江 杭州 310018）摘 要：为了实现爬杆机器人的结构紧凑和轻量化设计，以及克服垂直爬杆运动困难的问题，可通过减小零件厚度的方法来进行减重，但这会造成局部应力集中，导致机器人结构刚度不足。为此，从尺蠖的形态出发，对爬杆机器人的主要结构进行SIMP（solid isotropic material with penalization，具有惩罚作用的固体各向同性材料）变密度拓扑优化设计，在保证机器人整体性能的基础上实现轻量化设计。首先，选取爬杆机器人的夹持器与T关节等主要部件为对象，通过ADAMS仿真软件对3种运动步态进行动态仿真，模拟了在极限工况下3种运动步态引起的载荷变化。然后，根据极限工况载荷等信息设置边界条件，通过ANSYS Workbench软件中的Topology Optimization模块对爬杆机器人进行SIMP变密度拓扑优化，去除机器人结构中受载荷影响微小的冗余材料，对优化后的机器人模型进行重构并与优化前对比。结果表明，在爬杆机器人整机质量减小了7.6%（由13.60 kg减小为12.57 kg）的同时，受力较大的夹持器和T关节的综合性能有所提高。能耗测试实验结果表明，优化后爬杆机器人的运行能耗比优化前降低了7.0%。所提出的SIMP变密度拓扑优化方法在仿生爬杆机器人的结构设计中具有较高的参考价值。关键词：机器人；运动步态；结构分析；拓扑优化；变密度中图分类号：TH 164；TP 241 文献标志码：A 文章编号：1006-754X（2023）03-0342-11Structure optimization and analysis of pole-climbing robot based on SIMP methodZHENG Xiaofei,HUANG Zhenhai,MA Xiaolong,WANG Jianxin,WANG Binrui(College of Mechanical and Electrical Engineering,China Jiliang University,Hangzhou 310018,China)Abstract：In order to achieve the structural compactness and lightweight design of pole-climbing robots,and overcome the difficulty of vertical pole movement,the weight reduction can be realized by reducing the thickness of parts.However,it can cause local stress concentration,resulting in insufficient stiffness of robot structure.Therefore,starting from the morphology of inchworms,the SIMP(solid isotropic material with penalization)variable density topology optimization design was carried out for the main structure of the pole-climbing robot,realizing the lightweight design of the robot while ensuring its overall performance.Firstly,the main components such as the gripper and T-joint of the pole-climbing robot were selected as the objects,and three motion gaits were dynamically simulated by ADAMS simulation software to simulate the load changes caused by the three motion gaits under extreme working conditions.Then,the boundary conditions were set based on information such as extreme working condition loads.The Topology Optimization module in ANSYS Workbench software was used to perform SIMP variable density topology optimization on the pole-climbing robot to remove the redundant materials that were slightly affected by loads from the robot structure.The optimized robot doi:10.3785/j.issn.1006-754X.2023.00.042收稿日期：20230110 修订日期：20230215本刊网址在线期刊：http:/ 3 期郑小飞，等：基于SIMP方法的爬杆机器人结构优化与分析model was reconstructed and compared with before optimization.The results showed that while the overall weight of the pole-climbing robot was reduced by 7.6%(from 13.60 kg to 12.57 kg),the comprehensive performance of the gripper and T-joint with larger force was improved.The experimental results of energy consumption test showed that the operating energy consumption of the optimized pole-climbing robot decreased by 7.0%compared with before optimization.The proposed SIMP variable density topology optimization method has a high reference value in the structural design of biomimetic pole-climbing robots.Key words：robot;motion gait;structural analysis;topology optimization;variable density近年来，提高能源效率、降低材料成本和提升机器人结构紧凑性1推动了对轻量化机器人的需求，机器人的轻量化设计得到了广泛研究。爬杆机器人在高空杆上作业是一项重要的功能，这对机器人的负载提出了更高的要求。研究人员采用了大量的方法来使机器人在性能降低较小的同时实现轻量化。目前，轻量化技术已成功应用于各个领域，如汽车、航空航天等领域。应用于机器人领域的2种典型的轻量化技术为：材料轻量化和结构轻量化。材料轻量化能够获得较为良好的减重效果，因此机器人轻量化首选材料轻量化，即在保证机器人零部件力学性能的基础上优先选取轻质材料，达到轻量化的目的。具有材料轻量化代表的是工业机器人四大家族之一的德国KUKA机器人公司，其研究设计的轻型机器人LWR-32的躯干均采用轻型材料（碳纤维），负重比接近11，是当今负重比较高的机器人。材料轻量化的优势是可直接替换原有质量较大的金属材料以获得良好的减重效果，但弊端是轻质材料的成本较高且耐磨性比金属材料差，零部件反复拆装时容易磨损甚至发生断裂，这给机器人的后期维护带来了极大困难。而结构轻量化弱化了材料轻量化的弊端，在减重的同时保证机器人零部件的综合力学性能不变，降低了机器人的运行能耗，从而使机器人的综合性能得到提升，该方法是当前应用最广泛的机器人轻量化方法之一。与材料轻量化相比，结构轻量化只改变材料分布方式，通过去除冗余材料的方式直接达到轻量化目的，具有计算速度快、准确性高和易实现等优势。目前，常用的结构优化方法可归类为以下3种：尺寸优化、形状优化、拓扑优化3。这3种方法适用的设计阶段和应用场合有所区别。其中：尺寸优化是在不改变零部件形状和材料的情况下，对内部结构尺寸进行优化设计，其设计变量通常选择厚度、孔径和深度等，但是该方法忽略了结构的形状和连接性；形状优化是指确定零部件的拓扑结构（即孔的数量和位置），并优化每个孔的形状；拓扑优化是指在给定设计域内，基于施加的约束条件（通常为质量保留百分比、刚度和变形量等）对材料分布与结点连接方式进行优化，该方法只研究材料的最佳分布情况。与尺寸优化和形状优化相比，拓扑优化不仅允许更大的设计自由度，而且能改变材料的分布状态。因此，采用拓扑优化方法进行机器人结构优化设计具有更高的层次。随着计算机软件技术的发展及结构优化理论的逐渐成熟，拓扑优化方法已被广泛应用于各类机器人的结构轻量化设计4。赖欣等5利用变密度拓扑优化法，先将爬壁机器人分成若干个零部件，再根据各零部件的受力情况逐个进行优化，解决了机器人自身质量过大、结构复杂等问题，达到了轻量化目的。结果表明，优化后爬壁机器人的整机质量从1.14 kg减小为0.58 kg，轻量化效果显著，且优化后结构更为精简，为后续的装配环节提供了操作空间。Tran等6采用拓扑优化方法对Delta机器人的上臂进行了设计，在保证机器人刚度和耐用性的同时节省了材料和成本，重新设计后上臂的总质量为0.23 kg，相较于重新设计前减小了48.8%。Chu等7提出了一种以机器人结构的变形量、振动频率、质量为约束条件的多目标拓扑优化方法，优化后的机器人在实现减重的同时减小了整体变形量，提高了固有频率。黄永平等8利用变密度拓扑优化法对建筑重载机器人的主体结构进行了优化，使得其大臂的质量减小了 29.6%，小臂的质量减小了 33.4%，并进一步通过动力学仿真实验分析得到，优化后机器人的能耗明显改善。综上所述，拓扑优化方法非常适用于机器人的结构优化设计。为实现爬杆机器人的结构紧凑和轻量化设计，克服其垂直爬杆困难等问题，笔者以自制仿尺蠖爬杆机器人为研究对象，以质量最小与柔度最小为目标函数，以3种运动步态的实际极限工况为边界约束条件，采用ANSYS Workbench软件对机器人关键零部件的结构进行变密度拓扑优化设计。同时，对优化后的机器人进行力学性能校验，以验证其减 343工程设计学报第 30 卷 重降耗效果以及综合性能提升效果。1 拓扑优化数学方法与分析过程 1.1变密度拓扑优化问题的目标函数及约束一般来说，基于连续体的变密度拓扑优化问题的基本数学模型包含目标函数、边界条件约束和设计变量。基于线性静力有限元分析的变密度拓扑优化问题的一般公式可以描述为：min=f(),Ust:K()U=F()g(),U 0 0 1(1)式中：f为目标函数；为材料相对密度，=1表示满材料状态，=0表示空白状态，即无材料；U为全局位移矢量；F为载荷矢量；g为约束函数；K为全局结构刚度矩阵。在上述广义描述中，可以考虑各种目标和约束，包括柔度、应力、频率和位移等。1.2变密度拓扑优化数学模型拓扑优化9的目的是在满足一定约束条件的前提下，根据载荷的分布方式将材料最佳地分布在约束支撑位置，使得结构性能达到最优。拓扑优化的一般流程为：根据定义的设计变量以及优化目标构建优化目标函数，在满足约束的条件下进行迭代求解，获得最优解后停止迭代。在连续体拓扑优化方法中，典型的方法包括均匀化方法10、变密度法11、变厚度法12、峰值函数插值法13和渐进优化法14等。通常变密度材料插值模型可分为 2类：SIMP（solid isotropic material with penalization，具有惩罚作用的固体各向同性材料）模型15和 RAMP（rational approximation of material properties，材料属性有理近似）模型 16。其中，采用SIMP变密度法可使材料相对密度取0,1中的任意值，当弹性模量和材料相对密度之间的关系恒定时，SIMP变密度法的插值模型可表示为：E(xi)=Emin+xi(E0-Emin),i 1,N(2)K=i=1N(Emin+xpiE)ki(3)式中：xi为设计变量，即设计单元 i的相对密度，xi0,1；p为中间密度材料的惩罚因子，可使设计变量xi趋近于0或1；E(xi)为插值后得到的设计单元i的弹性模量；Emin为空白单元的弹性模量，通常情况下Emin取极小值，以确保优化过程中不会出现刚度矩阵奇异；E0为实体单元的初始弹性模量，E=E0-Emin；ki为设计单元i的刚度矩阵。本文以质量最小和柔度最小作为目标函数17，以应力条件和载荷为边界约束条件，建立基于SIMP变密度法的结构拓扑优化数学模型：Find:x=x1 x2 xNT NMinimize:C()x=UTKU=i=1NE()xiuTik0uis.t:V()x=i=1Nvixi V0 F=KU 0 xmin xi 1 max(4)式中：x为设计变量矩阵；为连续体设计域；C(x)为结构柔度；N为设计单元个数；k0为单位弹性模量单元的刚度矩阵；ui为设计单元i的位移矢量；V(x)为优化后的结构体积；vi为设计单元i的相对体积；V0为设计域的总体积；为预设的体积约束分数；xmin为单元相对密度的最小下限；max为结构的最大应力；为结构的许用应力。2 爬杆机器人动力学分析与结构优化 2.1爬杆机器人拓扑优化流程本文以 5-DOF（degree of freedom，自由度）串联爬杆机器人为研究对象。该机器人由2个夹持器、3个T关节（由摆动件和连接座组成）和2个I关节串联组成，其结构如图1所示。经分析，爬杆机器人主要靠夹持器与T关节支撑。选取负载最大的零部件（夹持器和T关节）作为优化对象，采用SIMP变密度拓扑优化方法对爬杆机器人进行结构优化，具体流程如图2所示。根据图2，基于SIMP变密度法的爬杆机器人拓扑优化的主要步骤如下：1）根据爬杆机器人各零部件的连接关系，在三维软件中建立机器人整机三维模型。图1自制5-DOF串联爬杆机器人Fig.1Self-made 5-DOF series pole-climbing robot 344第 3 期郑小飞，等：基于SIMP方法的爬杆机器人结构优化与分析2）在ADAMS软件中对爬杆机器人的3种运动步态进行动力学仿真，分析其主要零部件在极限运动工况下的受力情况，获得极限工况信息。3）根据极限工况信息，在ANSYS Workbench软件中对夹持器及T关节设置边界约束条件并进行力学分析，以得到所受的载荷、扭矩、等效应力和变形量的分布情况，从而确定这2个零部件是否可以进行拓扑优化。4）利用SIMP变密度法对夹持器及T关节进行拓扑优化，并根据拓扑优化结果，对原夹持器及T关节的模型进行改进，并重新建立新的模型。5）在ANSYS Workbench软件中对优化后的夹持器及T关节进行结构强度校核，判断其是否满足要求：若不满足，则根据拓扑优化结果重新对结构进行改进；若满足，则视新结构为最终设计方案。2.2ADAMS动力学仿真根据爬杆机器人的机构特点，设计3种运动步态，分别为蠕动步态、翻转步态和旋转步态，如图3所示。由于爬杆机器人（整机质量为13.60 kg）主要在高空环境中进行作业，其外部负载主要来自自身重量。利用ADAMS软件对爬杆机器人的3种运动步态进行动力学仿真，得到这3种运动步态下夹持器和T关节分别在t=3.5，1.5，0.5 s时的受力较大，即机器人运行的极限工况。在爬杆机器人攀爬过程中，各夹持器和T关节均受到1个六维力（Fx、Fy、Fz，Tx、Ty、Tz）的作用18，其中Fx、Fy、Fz分别表示沿x、y、z方向的载荷，Tx、Ty、Tz分别表示绕x、y、z轴的扭矩。通过ADAMS仿真分析得到夹持器的载荷和扭矩变化曲线，如图4所示。蠕动步态、翻转步态和旋转步态下3个T关节的载荷和扭矩变化曲线如图5所示。从图4中可以看出，夹持器所受的最大载荷为750 N左右，符合设计要求。由图5可得，T关节在翻转步态极限工况下所受的载荷最大，其峰值分别为Fx=-974 N、Fy=-43 N、Fz=-713 N；对 应 的 扭 矩分别为Tx=-14 637 Nmm、Ty=-296 Nmm、Tz=-58 677 Nmm。将获得的极限工况载荷、扭矩数据作为约束条件，对夹持器和T关节进行拓扑优化分析。图4夹持器的载荷和扭矩变化曲线Fig.4Load and torque change curves of gripper图2基于SIMP变密度法的爬杆机器人拓扑优化流程Fig.2Topology optimization process of pole-climbing robot based on SIMP variable density method图3爬杆机器人的3种运动步态Fig.3Three motion gaits of pole-climbing robot 345工程设计学报第 30 卷 2.3ANSYS有限元分析利用ANSYS Workbench软件中的Static Structual模块对爬杆机器人的夹持器和T关节进行静力学分析19。将通过ADAMS动力学仿真得到的翻转步态极限工况下的载荷、扭矩施加在夹持器的实际接触位置处。在实际工作过程中，夹持器的接触为点或线接触。为模拟真实工况，在 ANSYS Workbench软件中将夹持器前端的2个销孔作固定支撑；根据夹持器在实际夹持工况下所受的载荷，可将其接触类型视为线接触；在前处理中，在夹持接触面处添加2条辅助直线，用于载荷施加。通过运算求解，得到该极限工况下夹持器的等效应力及等效变形云图，结果如图6所示。由图6(a)可知，夹持器的最大等效应力主要分布于夹持器两端连接处，为49.992 MPa；对应的最大等效变形主要集中在夹持器末端，为0.25 mm。在分析T关节摆动件时，根据实际工况将摆动件前端的固定孔作固定支撑，在轴承孔处施加由ADAMS动力学仿真得到的翻转步态极限工况下的载荷及扭矩。通过运算求解，得到该极限工况下T关节摆动件的等效应力及等效变形云图，结果如图7所示。由图7可知，T关节摆动件的最大等效应力为37.923 MPa；等效变形集中在T关节摆动件末端处，最大值为0.23 mm。同理，按照ADAMS动力学仿真结果，将T关节连接座前端的螺丝定位孔固定，并在轴承孔处施加翻转步态极限工况下的载荷及扭矩。通过求解运算，得到该极限工况下T关节连接座的等效应力及等效变形云图，结果如图8所示。由图8可知，T关节连接座的等效应力集中在定位孔周边，最大值为77.897 MPa，因此建议增加此部位的定位孔，以减少应力集中；等效变形同样集 中 在 T 关 节 连 接 座 末 端 远 处，最 大 值 为0.03 mm。图5不同运动步态下T关节的载荷和扭矩变化曲线Fig.5Load and torque change curves of T-joint under different motion gaits 346第 3 期郑小飞，等：基于SIMP方法的爬杆机器人结构优化与分析本文爬杆机器人的夹持器及T关节均选用航空金属材料7075-T6铝合金，其密度为2.85 g/cm3，泊松比为 0.33，屈服强度为 504 MPa。综合上文的ANSYS有限元仿真结果可知，夹持器及T关节的最大等效应力远小于材料的屈服强度，因此可对其进行SIMP变密度拓扑优化。2.4结构优化分析通过前文对爬杆机器人运行工况的综合分析，可知夹持器和T关节这2个关键零部件承受了大部分的载荷（包括自身重力以及负载）。基于爬杆机器人的ADAMS模型，利用ANSYS Workbench软件对夹持器和关节进行SIMP变密度拓扑优化，优化目标是在体积约束下使柔度最小。按照约束条件，对设计变量进行寻优求解，求解器选用典型的优化准则法；设求解器最优迭代数为200次，最小相对密度为0.001，计算精度为0.15%。根据实际夹持工况，在Topology Optimization模块中设置优化区域为夹持器的整个材料区域，选择夹持面和支撑夹持器的2个销孔为优化排除域，响应约束设置为以质量最小为目标，保留质量百分比设置为45%。材料相对密度在01内变化，当相对密度接近1时，表明材料应保留；相反，当相对密度接近0时，表明材料可被去除。通过21次迭代计算，得到夹持器的变密度拓扑优化结果，如图9所示。由图9可知，夹持器在受载荷微小的材料区域出现明显桁架空穴，说明原夹持器具有拓扑性质。图6夹持器等效应力及等效变形云图Fig.6Cloud maps of equivalent stress and equivalent deformation of gripper图7T关节摆动件的等效应力及等效变形云图Fig.7Cloud maps of equivalent stress and equivalent deformation of T-joint swinging component图8T关节连接座的等效应力及等效变形云图Fig.8Cloud maps of equivalent stress and equivalent deformation of T-joint connection seat 347工程设计学报第 30 卷 根据拓扑优化云图，重新构建夹持器模型，去除原夹持器中受载荷微小的材料区域，使夹持器在保持力学性能不变的同时，达到减重和提高综合性能的设计目的。优化后夹持器的模型如图10所示。采用同样的方法对T关节进行SIMP变密度拓扑优化。根据实际运行工况，在Topology Optimization模块中设置优化区域为T关节的整个材料区域，选择底部固定孔和轴承孔作为优化排除域，边界条件约束设置为以质量最小为目标，保留质量百分比设为45%。T关节摆动件和T关节连接座的变密度拓扑优化结果分别如图11和图12所示。通过图11和图12可以看出，T关节出现了明显孔洞，说明其具有拓扑性质。根据图11和图12所示的拓扑优化云图，重新构建T关节模型，在保证其力学性能不变的基础上，去除受载荷不大的材料区域，达到轻量化目的。优化后T关节摆动件和T关节连接座的模型分别如图13和图14所示。3 爬杆机器人结构优化结果分析 利用ANSYS Workbench软件对优化后的夹持器和T关节进行有限元分析，以校验其强度是否满足要求，并分析存在问题的区域并进行修改。经过重新校验得到优化后夹持器的等效应力及等效变形云图，如图15所示。优化前后夹持器的质量和力学性能对比如表1所示。结果表明，优化后夹持器的质量由 0.35 kg 减小为 0.25 kg，减小了28.6%；最大等效应力出现在夹持器连接端之间的尖角处，为 59.758 MPa，远小于材料许用应力；等效变形量增大了0.07 mm，对整个爬杆机器人结构的影响微小；整个夹持器的力学性能良好。综上，建议将夹持器连接端之间的尖角修改为圆角，以减少应力集中。图9夹持器拓扑优化云图Fig.9Topology optimization cloud map of gripper图10优化后夹持器模型Fig.10Optimized gripper model图11T关节摆动件拓扑优化云图Fig.11Topology optimization cloud map of T-joint swinging component图12T关节连接座拓扑优化云图Fig.12Topology optimization cloud map of T-joint connection seat图13优化后T关节摆动件模型Fig.13Optimized T-joint swinging component model图14优化后T关节连接座模型Fig.14Optimized T-joint connection seat model 348第 3 期郑小飞，等：基于SIMP方法的爬杆机器人结构优化与分析优化后T关节摆动件的等效应力及等效变形云图如图16所示，优化前后T关节摆动件的质量和力学性能如表2所示。由图16可以看出，T关节摆动件的最大等效应力出现在轴承的两端连接处，为88.911 MPa，相比于优化前的 37.923 MPa 增大了50.988 MPa，等效应力分布规律与变密度拓扑优化云图中的材料密度分布规律基本一致；最大等效变形量基本不变。由表2可知，优化后T关节摆动件的质量由0.27 kg减小为0.16 kg，减小了40.7%，轻量化效果较好。优化后T关节连接座的等效应力及等效变形云图如图17所示，优化前后T关节连接座的质量和力学性能对比如表3所示。对比图8和图17可知，经过变密度拓扑优化，T关节连接座的最大等效应力从77.897 MPa增大到142.880 MPa，后端固定孔区域出现了较大的等效应力。建议增加后端固定孔，以减少应力集中。T关节连接座的最大等效变形极其微小，可以忽略不计。由表3可知，优化后T关节连接座的质量由0.34 kg减小为0.23 kg，减小了32.4%。上述拓扑优化结果表明，在保证结构强度的前提下，采用SIMP变密度拓扑优化方法可显著减小爬杆机器人的质量，并提高其综合性能。4 实验验证 根据拓扑优化结果，重新设计爬杆机器人的夹持器和T关节，并完成加工和装配。最后，在自制的攀爬实验平台上测试爬杆机器人的性能，对其质量和能耗等进行对比。优化后爬杆机器人的质量为12.57 kg，相较于优化前减小了7.6%。图15优化后夹持器的等效应力及等效变形云图Fig.15Cloud maps of equivalent stress and equivalent deformation of optimized gripper表 1优化前后夹持器的质量和力学性能对比Table 1Comparison of quality and mechanical properties of gripper before and after optimization优化前后优化前优化后变化率/%质量/kg0.350.25-28.6最大等效应力/MPa49.99259.75819.5最大等效变形量/mm0.250.32图16优化后T关节摆动件的等效应力及等效变形云图Fig.16Cloud maps of equivalent stress and equivalent deformation of optimized T-joint swinging component表2优化前后T关节摆动件的质量和力学性能对比Table 2Comparison of quality and mechanical properties of T-joint swinging component before and after optimization优化前后优化前优化后变化率/%质量/kg0.270.16-40.7最大等效应力/MPa37.92388.911134.4最大等效变形量/mm0.230.49 349工程设计学报第 30 卷 4.1爬杆机器人能耗仿真实验将优化后的爬杆机器人导入ADAMS软件，并添加相关约束和设置材料参数，其能耗仿真模型如图18所示。对优化后的爬杆机器人进行动力学仿真分析，并以蠕动步态工况为例，对比分析优化前后机器人T关节电机的功率，以分析其能耗情况，结果如图19所示。由图可知，相较于优化前，优化后爬杆机器人T关节电机的功率有所下降（约降低了7%），达到了能耗降低的目的。4.2爬杆机器人样机能耗实验在爬杆机器人样机中，采用型号为QDD Pro-PR60-100-90 的谐波减速电机，其额定扭矩为47 Nm，输出功率为147 W。使用STM32VET6主控板接收 QT5.9.0 版本的上位机指令并通过 CAN（controller area network，控制器局域网）接口接收电机的驱动信号。在攀爬实验平台上开展爬杆机器人样机能耗实验（以蠕动步态工况为例），如图20所示。结果表明，蠕动步态符合仿生学设计，步态较柔顺，该步态下夹持器夹持稳定，未出现上下滑落现象。由于缺乏有效的能耗测量方案，本文通过读取爬杆机器人各关节电机内置函数值，并利用单片机将其转换为电压数据发送到上位机，经MATLAB软件处理得到功率曲线，结果如图21所示。由于能耗仿真实验中采用step函数来添加驱动，爬杆机器人的仿真轨迹与实际轨迹存在差异，故其关节电机功率的仿真结果与实验结果不完全吻合，但从总体上看，两者的相对误差在允许范围内，验证了仿真结果的正确性，表明优化后机器人可实现能耗降低的目标。图17优化后T关节连接座的等效应力及等效变形云图Fig.17Cloud maps of equivalent stress and equivalent deformation of optimized T-joint connection seat表3优化前后T关节连接座的质量和力学性能对比Table 3Comparison of quality and mechanical properties of T-joint connection seat before and after optimization优化前后优化前优化后变化率/%质量/kg0.340.23-32.4最大等效应力/MPa77.897142.88083.4最大等效变形量/mm0.030.05图18优化后爬杆机器人的能耗仿真模型Fig.18Energy consumption simulation model of optimized pole-climbing robot图19优化前后爬杆机器人能耗仿真结果Fig.19Simulation results of energy consumption of pole-climbing robot before and after optimization图20爬杆机器人样机蠕动步态能耗实验Fig.20Energy consumption experiment on peristaltic gait of pole-climbing robot prototype 350第 3 期郑小飞，等：基于SIMP方法的爬杆机器人结构优化与分析5 结 论 本文以自行开发的仿尺蠖爬杆机器人为研究对象，对其结构的拓扑优化进行分析。实验结果表明，所采用的SIMP变密度拓扑优化方法对机器人结构优化有显著效果。主要研究过程及结论如下：1）构建了爬杆机器人的动力学模型，通过ADAMS软件仿真得到了该机器人在3种运动步态极限工况下的载荷和扭矩，为机器人结构拓扑优化的边界条件设置提供了基本依据。2）利用ANSYS Workbench软件对爬杆机器人的力学性能进行有限元分析，得到了夹持器和T关节的等效应力和等效变形分布情况。3）根据有限元分析结果，采用 SIMP 变密度拓扑优化方法对夹持器和 T 关节进行拓扑优化，去除了冗余材料，从而实现爬杆机器人的轻量化和间接提高负载自重比的目标，并通过实验验证了爬杆机器人可在减小质量的同时降低能耗。本文研究方法可为高空爬杆机器人的结构优化设计提供参考。参考文献：1 VERTUT J，LIEGEOIS A.General design criteria for manipulators J.Mechanism and Machine Theory，1981，16（1）：65-70.2 POVSE B，HADDADIN S，BELDER R，et al.A tool for the evaluation of human lower arm injury：approach，experimental validation and application to safe roboticsJ.Robotica，2016，34（11）：2499-2515.3 杨姝.复杂机械结构拓扑优化若干问题研究 D.大连：大连理工大学，2007：1-3.YANG S.Study on topology optimization of complex mechanics structure D.Dalian：Dalian University of Technology，2007：1-3.4 姚屏，林源，汤勇，等.结构拓扑优化方法及其机器人轻量化应用现状及发展 J.机电工程技术，2018，47（11）：117-122.YAO P，LIN Y，TANG Y，et al.Current situation and improving trends of structural topology optimization method and the corresponding application in robot light-weight designJ.Mechanical&Electrical Engineering Technology，2018，47（11）：117-122.5 赖欣，师靖远，彭天宇，等.爬壁机器人变密度拓扑优化吸附结构研究J.机械科学与技术，2021，40（6）：821-827.LAI X，SHI J Y，PENG T Y，et al.Research on adsorption structure of wall climbing robot with variable density topology optimizationJ.Mechanical Science and Technology for Aerospace Engineering，2021，40（6）：821-827.6 TRAN T H T，NGUYEN D S，VO N T，et al.Design of Delta robot arm based on topology optimization and generative design methodC/2020 5th International Conference on Green Technology and Sustainable Development（GTSD），Ho Chi Minh City，Nov.27-28，2020.7 CHU X Y，XU H H，SHAO G F，et al.Multi-objective topology optimization for industrial robot C/2016 IEEE International Conference on Information and Automation（ICIA），Ningbo，Aug.1-3，2016.8 黄永平，郭秋明.基于拓扑优化的起重机转台的轻量化设计 J.机电工程技术，2022，51（10）：211-214.HUANG Y P，GUO Q M.Lightweight design of the crane turntable based on topology optimazation J.Mechanical&Electrical Enginee