一种柔性铰链的拓扑优化设计方法
【专利摘要】本发明公开了一种柔性铰链的拓扑优化设计方法，包括以下步骤：步骤1：建立柔性铰链拓扑优化设计模型，将典型切口柔性铰链轮廓设为设计域的形状，定义刚性区域(非设计域)；步骤2：建立柔性铰链拓扑优化有限元模型；步骤3：基于有限元模型建立柔性铰链拓扑优化问题数学模型；步骤4：计算柔性铰链拓扑优化问题的灵敏度；步骤5：采用优化算法求解柔性铰链拓扑优化问题，更新设计变量，得到最终拓扑结果图；步骤6：根据拓扑优化得到的最终拓扑结果图，提取其轮廓，通过适当的修改得到新型柔性铰链。本发明采用拓扑优化方法在概念层面上设计柔性铰链，可以设计出结构更复杂、性能更优越的新型柔性铰链，具有更大柔度，更高精度以及更小的最大应力。
【专利说明】
一种柔性铰链的拓扑优化设计方法
技术领域
[0001] 本发明涉及柔性铰链结构设计领域，特别涉及一种柔性铰链的拓扑优化设计方 法。
【背景技术】
[0002] 柔性铰链是利用自身的弹性变形来完成运动和力的传递与转换，通常由一块板材 料通过线切割加工而成，与传统的刚性运动副结构有很大的不同。柔性铰链是集中式柔顺 机构中最为关键的柔性单元，其性能对柔顺机构有着重要的影响。与传统的刚性运动副比 较，柔性铰链具有质量轻、无摩擦磨损、可一体化加工、精度高等众多优点，引起了国内外学 者的广泛关注。柔性铰链在精密工程、仪器测量、航空航天以及微机电系统(MEMS)等领域有 着广泛的应用。
[0003] 目前各种切口型柔性铰链已被国内外学者广泛研究，如圆型、直圆型、椭圆型、角 圆型、叶型、V型、抛物线型、双曲线型、混合型等柔性铰链。柔性铰链性能指标主要有三个： 柔度(或刚度），转动精度、应力水平。一个理想的柔性铰链应该有大的转动柔度、高的转动 精度和低的应力水平。但是，由于柔性铰链自身的弹性变形，要同时满足这三个指标是非常 困难的。例如在常见的柔性铰链中，直圆型和V型柔性铰链虽然有着很高的转动精度，但是 它们的转动柔度小且应力很大，这就大大限制了它们的运动范围。另外，传统的柔性铰链设 计方法主要依赖于设计者的经验来确定柔性铰链的构型，构型都非常简单，通常为不同形 状的切口型。此时，柔性铰链的性能很大程度上取决于设计者的经验，这样很难设计出构型 复杂、性能优越的柔性铰链。
【发明内容】
[0004] 本发明的目的在于克服现有技术存在的问题，提供一种柔性铰链的拓扑优化方法 来设计新型的柔性铰链，使其在保持较高的转动精度的同时提高转动柔度和降低最大应 力。
[0005] 为实现上述目的，本发明采用如下技术方案：
[0006] -种柔性铰链的拓扑优化设计方法，包括以下步骤：
[0007]步骤1:建立柔性铰链拓扑优化设计模型，将典型切口柔性铰链轮廓设为设计域的 形状，定义刚性区域(非设计域）；
[0008] 步骤2:建立柔性铰链拓扑优化有限元模型；
[0009] 步骤3:基于有限元模型建立柔性铰链拓扑优化问题数学模型；
[0010] 步骤4:计算柔性铰链拓扑优化问题的灵敏度；
[0011] 步骤5:采用优化算法求解柔性铰链拓扑优化问题，更新设计变量，得到最终拓扑 结果图；
[0012] 步骤6:根据拓扑优化得到的最终拓扑结果图，提取其轮廓，通过适当的修改得到 新型柔性铰链。
[0013] 进一步地，步骤1中所述的将典型切口柔性铰链轮廓设为设计域的形状的步骤具 体包括:将圆，椭圆，抛物线，双曲线，V型等常见形状设为设计域的形状，非设计域设置为矩 形，其高度和设计域高度相等，长度设为设计域高度的3倍或以上。
[0014] 进一步地，步骤2中所述的建立柔性铰链拓扑优化有限元模型的步骤具体包括:将 设计域和非设计域分别离散为N和η个有限单元，非设计域的相对密度预先设为 Xl = l，在输 入端X和Y方向和输出端X和Y方向分别添加4个虚拟弹簧用于模拟铰链与工件之间的间隙 和反作用力。
[0015] 进一步地，步骤3中所述的基于有限元模型建立柔性铰链拓扑优化问题数学模型 的步骤具体包括:最大化设计域右端中点位移且最小化设计域右端中点位移作为 目标函数，体积比为约束条件，其优化数学模型如下：
[0017] s. t:Kux=Fx
[0018] Kuy = Fy
[0019] fv(x)=vTx<V*
[0020] 0<xmin<xi<l，i = l，2.·_，N，
[0021] 其中:w是权重系数，Fx为在非设计域右端中点X方向施加的载荷，Fy为在非设计域 右端中点Y方向施加的载荷，为在载荷F x作用下产生的设计域右端中点X方向的位移， 为在载荷Fy作用设计域右端中点Y方向的位移，￡/^和％%是拓扑优化迭代中的初始 值，f为体积比，K为全局刚度矩阵，ux是由载荷Fx产生的位移场，uy是由载荷F y产生的位移 场，N为有限单元个数，Xi为各有限单元的相对密度，Xmin为最小相对密度，通常设为0.001。
[0022] 进一步地，所述的〖人斤=(Χ,)Ρ ??ΛΓ?ΜΛ..
[0023] 进一步地，步骤4中所述的计算柔性铰链拓扑优化问题的灵敏度的步骤采用伴随 矩阵法求解目标函数和约束条件的灵敏度，其目标函数敏度分析为：
[0027] 进一步地，步骤5中所述的优化算法可以采用优化准则法(0C)或者移动渐进算法 (MMA)〇
[0028] 进一步地，步骤6中所述的根据拓扑优化得到的最终拓扑结果图，提取其轮廓，通 过适当的修改得到新型的柔性铰链的步骤具体包括:柔性铰链设计域的轮廓保持不变，通 过适当修改拓扑优化得到的切口或者孔洞，设计出一种新型的柔性铰链。
[0029] 与现有的技术相比，本发明的有益效果在于：
[0030] 在典型的切口型柔性铰链的基础上去除部分材料，并且采用属于概念设计方法拓 扑优化来最佳的去除给定的材料，得到一种性能更优的新型柔性铰链。与优化之前的柔性 铰链相比，通过本发明得到的新型柔性铰链有更高的转动柔度和精度，更低的应力水平。此 外，现有的技术设计柔性铰链在很大程度上都是依靠设计者的经验来决定柔性铰链的拓扑 构型，而本发明是从拓扑层面上来改变柔性铰链的构型，给定目标函数和约束条件，通过拓 扑优化方法可以自动生成满足给定要求的铰链，它为柔性铰链的设计提供了一种全新的设 计思路和方法。
【附图说明】：
[0031] 图1是本发明方法实施例中设计域与非设计域，尺寸，边界条件以及载荷等示意 图。
[0032] 图2是本发明方法实施例不同参数下的拓扑优化结果示意图。
[0033] 图3是本发明方法实施例中得到的一种具体的新型柔性铰链平面示意图。
[0034] 图4是本发明方法实施例中得到的一种具体的新型柔性铰链三维模型示意图。
【具体实施方式】
[0035] 为更好地理解本发明，下面结合附图对本发明作进一步说明，但是本发明的实施 方式不限于此。
[0036] 参照图1-4,一种柔性铰链的拓扑优化设计方法，具体步骤如下：
[0037]步骤1:建立柔性铰链拓扑优化设计模型，将典型柔性铰链轮廓设为设计域，此实 例以直圆型柔性铰链的轮廓作为设计域1，但本发明的设计域不限于此。定义设计域1和非 设计域2(刚性构件）的尺寸，设计域1的半径R设为40mm，最小厚度t设为20mm，非设计域长度 Is为300mm，厚度均为1mm;
[0038] 步骤2:建立柔性铰链拓扑优化有限元模型，设计域1和非设计域2材料的杨氏模量 均为lGPa，泊松比为0.3，将设计域离散为80 X 100个4节点正方形有限单元，非设计域离散 为300X100个4节点正方形有限单元。每个有限单元的相对密度Xl(0< Xl<l)为一个设计变 量。由于非设计域为刚性区域，所以非设计域的相对密度预先设为1。在非设计域右端中点 施加载荷F X = 50N和Fy=lN，在设计域1右端中点添加虚拟输出弹簧，弹簧刚度设为 =? =丨，在非设计域右端中点添加虚拟输入弹簧，弹簧刚度设为％ =? = L，通过求 解得到位移场ux(由载荷Fx产生的)和uy(由载荷Fy产生的）；
[0039] 步骤3:基于有限元模型建立柔性铰链拓扑优化问题数学模型，首先定义α,,?；为 载荷Fx在设计域右端中点产生的X方向位移：
[0040] 心=(X,)W、
[0041 ]定义为载荷Fy在设计域右端中点产生的Y方向位移：
[0042] Uy Fi ^(xi)P l~l/Kuy _
[0043] 一个理想柔性铰链在转动方向应该有较大的柔度而在非期望运动的方向应该有 较大的刚度(较小的柔度），因此定义柔性铰链拓扑优化的目标函数为最大化位移UvA且 最小化位移，约束条件为材料体积比，所以柔性铰链拓扑优化问题数学模型如下：
[0045] s. t:Kux=Fx
[0046] Kuy = Fy
[0047] fv(x)=vTx<V*
[0048] 0<xmin<xi<l，i = l，2.··，N，
[0049] 其中:w是权重系数，Fx为在非设计域右端中点X方向施加的载荷，Fy为在非设计域 右端中点Y方向施加的载荷，巧,#;为在载荷F x作用下产生的设计域右端中点X方向的位移， 为在载荷Fy作用设计域右端中点Y方向的位移，巧是拓扑优化迭代中的初始 值，f为体积比，K为全局刚度矩阵，ux是由载荷Fx产生的位移场，uy是由载荷F y产生的位移 场，N为有限单元个数，Xi为各有限单元的相对密度，Xmin为最小相对密度，通常设为0.001 ;
[0050] 步骤4:优化问题敏度分析计算，采用伴随矩阵法求解目标函数和约束条件的灵敏 度，目标函数敏度分析：
[0054]步骤5:此实例采用优化准则(0C)算法求解柔性铰链拓扑优化问题，更新设计变量 ，直到迭代收敛得到最终拓扑结果。图2所示为不同体积比V*和t/R比值的拓扑结果图； [0055]步骤6:根据拓扑优化得到的最终拓扑结果图，提取其轮廓，通过适当的修改得到 新型的柔性铰链。
[0056]以下给出通过适当处理柔性铰链拓扑优化结果图得到的一种具体的新型柔性铰 链，并且和优化之前的直圆型柔性铰链对比来验证本发明的有效性。
[0057]在此实例中，由于最终的拓扑结果图中被去除材料的切口形状最接近抛物线，因 此柔性铰链左右两端的切口形状被设为抛物线，且整个铰链左右上下对称。图3所示为此实 例通过对拓扑优化结果进行后处理得到的一种新型柔性铰链的平面示意图，图4为它的三 维模型图。
[0058]将此实例得到新型柔性铰链和直圆型柔性铰链三维建模，导入到商业有限元分析 软件ANSYS 13.0进行有限元分析，计算出两个铰链的转动柔度，转动精度(数值越小精度越 高)和最大应力这三个性能，最大应力是在转角位移（1. 111 X 104rad)相同的情况下比较， 对比结果由表格1给出。
[0059]表格 1:
[0061]通过比较发现由本发明得到的新型柔性铰链比直圆型柔性铰链有着更大的转动 柔度，更高的转动精度和更低的最大应力，验证了本发明的有效性。
[0062]最后说明的是，以上优选实施例仅用以说明本发明的技术方案而非限制，尽管通 过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可以在 形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。凡 对本发明所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种柔性较链的拓扑优化设计方法，其特征在于，包括W下步骤： 步骤1:建立柔性较链拓扑优化设计模型，将典型切口柔性较链轮廓设为设计域的形 状，定义刚性区域； 步骤2:建立柔性较链拓扑优化有限元模型； 步骤3:基于有限元模型建立柔性较链拓扑优化问题数学模型； 步骤4:计算柔性较链拓扑优化问题的灵敏度； 步骤5:采用优化算法求解柔性较链拓扑优化问题，更新设计变量，得到最终拓扑结果 图； 步骤6:根据拓扑优化得到的最终拓扑结果图，提取其轮廓，通过适当的修改得到新型 柔性较链。2. 根据权利要求1所述的一种柔性较链的拓扑优化设计方法，其特征在于，步骤1中所 述的将典型切口柔性较链轮廓设为设计域的形状的步骤具体包括:将圆，楠圆，抛物线，双 曲线，V型设为设计域的形状，非设计域设置为矩形，其高度和设计域高度相等，长度设为设 计域高度的3倍或W上。3. 根据权利要求1所述的一种柔性较链的拓扑优化设计方法，其特征在于:步骤2中所 述的建立柔性较链拓扑优化有限元模型的步骤具体包括:将设计域和非设计域分别离散为 N和η个有限单元，非设计域的相对密度预先设为xi = l，在输入端X和Y方向和输出端X和Y方 向分别添加4个虚拟弹黃用于模拟较链与工件之间的间隙和反作用力。4. 根据权利要求1所述的一种柔性较链的拓扑优化设计方法，其特征在于:步骤3中所 述的基于有限元模型建立柔性较链拓扑优化问题数学模型的步骤具体包括:最大化设计域 右端中点位移Uy,A且最小化设计域右端中点位移化,?；作为目标函数，体积比为约束条件， 其优化数学模型如下：s.t:Kux = Fx KUy = Fy fv(x)=v\^V* 0<Xmin《Xi《l , i = l ,2... ,Ν， 其中:w是权重系数，Fx为在非设计域右端中点X方向施加的载荷，Fy为在非设计域右端 中点Y方向施加的载荷，为在载荷Fx作用下产生的设计域右端中点X方向的位移，吗,5. 为在载荷Fy作用设计域右端中点￥方向的位移，U、%、和^/,%、是拓扑优化迭代中的初始值，V* 为体积比，K为全局刚度矩阵，化是由载荷Fx产生的位移场，Uy是由载荷Fy产生的位移场，N为 有限单元个数，XI为各有限单元的相对密度，Xmin为最小相对密度，通常设为0.001。5. 根据权利要求4所述的一种柔性较链的拓扑优化设计方法，其特征在于：所述的6. 根据权利要求1所述的一种柔性较链的拓扑优化设计方法，其特征在于:步骤4中所 述的计算柔性较链拓扑优化问题的灵敏度的步骤采用伴随矩阵法求解目标函数和约束条 件的灵敏度，其目标函数敏度分析为：7. 根据权利要求1所述的一种柔性较链的拓扑优化设计方法，其特征在于:步骤5中所 述的优化算法可W采用优化准则法或者移动渐进算法。8. 根据权利要求1所述的一种柔性较链的拓扑优化设计方法，其特征在于:步骤6中所 述的根据拓扑优化得到的最终拓扑结果图，提取其轮廓，通过适当的修改得到新型的柔性 较链的步骤具体包括:柔性较链设计域的轮廓保持不变，通过适当修改拓扑优化得到的切 口或者孔桐，设计出一种新型的柔性较链。
【文档编号】G06F17/50GK106096158SQ201610440183
【公开日】2016年11月9日
【申请日】2016年6月16日
【发明人】张宪民, 刘敏
【申请人】华南理工大学