一种基于最优化方法的回弹补偿系数获取方法
【专利摘要】一种基于最优化方法的回弹补偿系数获取方法，该方法将回弹补偿系数作为设计变量，以零件冲压回弹仿真后几何型面的节点分布与零件设计几何型面的节点分布之间的差值作为目标函数，将工程问题描述成数学问题，并基于最优化方法去搜寻满足设计要求的回弹补偿系数，将冲压回弹仿真分析过程嵌套到最优化算法中，实现冲压回弹仿真分析的智能化，避免了工程师范对回弹补偿系数的反复调整，可快速高效地获取冲压零件的最佳回弹补偿量，提高了工程师的分析效率，缩短了零件回弹补偿后的几何型面产生的周期。
【专利说明】一种基于最优化方法的回弹补偿系数获取方法
【技术领域】
[0001] 本发明主要涉及冲压件模具型面设计【技术领域】，尤其是指一种基于最优化方法 的回弹补偿系数获取方法。
【背景技术】
[0002] 所谓回弹，是指板料成形的最后阶段，随着变形力的消失，成形过程中存储的弹性 变形能释放出来，引发内应力的重组，进而导致零件整体形状的改变。因而在制作冲压零件 的过程中，当冲压零件从冲压模具中退出时，冲压零件会产生回弹，成型后冲压零件的尺寸 不再等于冲压模具几何型面的尺寸，如果按照零件设计的几何型面加工冲压模具的几何型 面，由于回弹的影响，那最后加工得到的零件尺寸将不满足设计的要求。近年来，随着高强 度钢板和铝合金板在汽车上的大量应用，汽车冲压零件的回弹控制便显得越来越重要。
[0003] 冲压零件的回弹控制方法大体来说有两种：一是工艺控制法，通过改变冲压工艺 参数来提高零件拉伸效应以减小回弹，但该方法无法完全消除回弹；二是模具型面补偿法， 通过对模具设计几何型面进行预修正，使得零件冲压回弹仿真后几何型面正好与零件设计 几何型面一致，这种方法可以从根本上消除回弹的影响。而在使用模具型面补偿法对冲压 零件进行回弹控制时，零件回弹补偿后几何型面的获取是补偿过程的难点，而零件回弹补 偿后几何型面与回弹补偿系数有关，零件回弹补偿后几何型面与回弹补偿系数的关系可表 示为：C = R+a (S-R)，其中，C表示零件回弹补偿后几何型面，R表示零件设计几何型面，S表 示零件冲压回弹仿真后几何型面，a表示回弹补偿系数。因而，若得知零件设计几何型面、 零件冲压回弹仿真后几何型面和回弹补偿系数，即可根据零件设计几何型面、零件冲压回 弹仿真后几何型面和回弹补偿系数获得零件回弹补偿后几何型面。由此可见，零件回弹补 偿后几何型面的获取关键在于回弹补偿系数的选取，合适地选取满足设计要求的回弹补偿 系数显得十分重要，而在现有的冲压回弹仿真过程中，工程师通常是根据实践经验以及参 考类似的工程设计来人工判断确定回弹补偿系数的，基于该回弹补偿系数，对零件设计几 何型面进行回弹补偿，生成零件回弹补偿后几何型面，然后，对零件回弹补偿后几何型面进 行冲压回弹仿真分析，得到零件冲压回弹仿真后几何型面，如果所获得的零件冲压回弹仿 真后几何型面满足设计要求，便停止对零件回弹补偿后几何型面的冲压回弹仿真分析，如 果获得的零件冲压回弹仿真后几何型面不满足设计要求，工程师需根据实践经验以及参考 类似的工程设计来人工重新调整回弹补偿系数,再次对零件回弹补偿后几何型面进行冲压 回弹仿真分析，直到获得的零件冲压回弹仿真后几何型面状满足设计要求。
[0004] 由此可见，现有获取满足设计要求的回弹补偿系数是一个工作繁琐、效率低下的 过程，由于工程师实验经验限制，很难一下子就寻找到满足设计要求的回弹补偿系数，如何 能快速高效寻找满足设计要求的回弹补偿系数是该领域技术人员急需解决的问题。
【发明内容】
[0005] 为克服上述现有技术中存有的缺陷，本发明的目的在于提供一种基于最优化方法 的回弹补偿系数获取方法，该方法可快速高效地获取到满足设计要求的回弹补偿系数，避 免了工程师人工对回弹补偿系数的反复调整，提高工作效率。
[0006] -种基于最优化方法的回弹补偿系数获取方法，该方法包括以下步骤：A.以回弹 补偿系数作为设计变量，以零件冲压回弹仿真后几何型面的节点分布与零件设计几何型面 的节点分布之间的差值作为目标函数；B.对零件设计几何型面进行网格划分，获得零件设 计几何型面的节点分布；C.运用有限元分析方法对零件设计几何型面的节点分布进行初 始的冲压回弹仿真分析，获取零件冲压回弹仿真后几何型面的节点分布；D.设定回弹补偿 系数的取值区域，并在该取值区域内随机生成回弹补偿系数取值，设定零件冲压回弹仿真 后几何型面与零件设计几何型面之间的最大允许误差；E.根据回弹补偿系数取值、零件设 计几何型面的节点分布和零件冲压回弹仿真后几何型面的节点分布获得零件回弹补偿后 几何型面的节点分布；F.对零件回弹补偿后几何型面的节点分布进行冲压回弹仿真分析， 获得新的零件冲压回弹仿真后几何型面的节点分布；G.判断零件冲压回弹仿真后几何型 面的节点分布与零件设计几何型面的节点分布之间的差值是否小于最大允许误差，若是， 此时的回弹补偿系数取值为需获取的回弹补偿系数，若否，则基于最优化方法对目标函数 的设计变量进行迭代计算，获得新的回弹补偿系数取值，再返回步骤E。
[0007] 上述方法将回弹补偿系数作为设计变量，而零件冲压回弹仿真后几何型面的节点 分布与零件设计几何型面的节点分布之间的差值作为目标函数，而由于回弹补偿系数的寻 求方向是使得零件冲压回弹仿真后几何型面等于或无限接近于零件设计几何型面，因而该 方法相当于将回弹补偿系数的寻找过程描述为一个最值问题，即回弹补偿系数的获取过程 使相当是基于最优化方法寻找零件冲压回弹仿真后几何型面与零件设计几何型面之间误 差最小的回弹补偿系数的最值问题，而为满足不同工程中不同的设计要求，因而本发明寻 求回弹补偿系数的终止标准为零件冲压回弹仿真后几何型面的节点分布与零件设计几何 型面的节点分布之间的差值满足对应的设计要求，即本发明可通过算法智能搜寻冲压回弹 仿真过程中满足设计要求的回弹补偿系数，因而利用本发明的方法可快速高效地获取满足 设计要求的回弹补偿系数，不需要工程师人工对回弹补偿系数的反复调整，提高工程师的 分析效率，缩短零件回弹补偿后几何型面产生的周期。
[0008] 所述步骤A前还包括有步骤A0 :设定零件设计几何型面的各节点的权数；为获取 零件冲压回弹仿真后几何型面的节点分布与零件设计几何型面的节点分布之间的差值，需 在整个冲压回弹仿真分析过程中只改变零件设计几何型面上网格节点的位置，而不增加或 减少网格节点的数量，因而零件设计几何型面上的节点与零件冲压回弹仿真后几何型面上 的节点、零件回弹补偿后几何型面的节点均一一对应，所对应的节点的权数也相对应；步骤 A中的目标函数结合步骤A0中的权数建立。因而可根据各节点误差要求程度不同进行设置 各节点的权数，因而可结合实际情况进行冲压回弹补偿仿真分析，更灵活。
[0009] 本发明所带来的有益效果为：
[0010] 1.通过算法智能搜寻冲压回弹仿真过程中满足设计要求的回弹补偿系数，避免了 工程师对回弹补偿系数的反复调整，快速高效地获取了冲压零件的最佳回弹补偿量，提高 了工程师的分析效率，缩短了零件回弹补偿后几何型面产生的周期；
[0011] 2.设计变量和目标函数的确定将工程问题描述为数学问题，更精确且便于实现；
[0012] 3.将冲压回弹仿真分析过程嵌套到最优化算法中，实现仿真分析的智能化，让回 弹补偿系数实现最优化方向的搜寻。
【专利附图】
【附图说明】
[0013] 图1为本发明实施例的流程图。
【具体实施方式】
[0014]本发明的一种基于最优化方法的回弹补偿系数获取方法包括以下步骤：A.以回 弹补偿系数作为设计变量，以零件冲压回弹仿真后几何型面的节点分布与零件设计几何型 面的节点分布之间的差值作为目标函数；B.对零件设计几何型面进行网格划分，获得零件 设计几何型面的节点分布；C.运用有限元分析方法对零件设计几何型面的节点分布进行 初始的冲压回弹仿真分析，获取零件冲压回弹仿真后几何型面的节点分布；D.设定回弹补 偿系数的取值区域，并在该取值区域内随机生成回弹补偿系数取值，设定零件冲压回弹仿 真后几何型面与零件设计几何型面之间的最大允许误差；E.根据回弹补偿系数取值、零件 设计几何型面的节点分布和零件冲压回弹仿真后几何型面的节点分布获得零件回弹补偿 后几何型面的节点分布；F.对零件回弹补偿后几何型面的节点分布进行冲压回弹仿真分 析，获得新的零件冲压回弹仿真后几何型面的节点分布；G.判断零件冲压回弹仿真后几何 型面的节点分布与零件设计几何型面的节点分布之间的差值是否小于最大允许误差，若 是，此时的回弹补偿系数取值为需获取的回弹补偿系数，若否，则基于最优化方法对目标函 数的设计变量进行迭代计算,获得新的回弹补偿系数取值，再返回步骤E。
[00彳5]所述步骤A前还包括有步骤A0 :设定零件设计几何型面的各节点的权数；为获取 零件冲压回弹仿真后几何型面的节点分布与零件设计几何型面的节点分布之间的差值，需 在整个冲压回弹仿真分析过程中只改变零件设计几何型面上网格节点的位置，而不增加或 减少网格节点的数量，因而零件设计几何型面上的节点与零件冲压回弹仿真后几何型面上 的节点、零件回弹补偿后几何型面的节点均一一对应，所对应的节点的权数也相对应；步骤 A中的目标函数结合步骤A0中的权数建立。因而可根据各节点误差要求程度不同进行设置 各节点的权数，因而可结合实际情况进行冲压回弹补偿仿真分析，更灵活。
[0016]步骤G中的最优化方法可以采用牛顿法、斐波纳契法、"0.618，，法、变尺度法或共 轭方向法等，本实施例以牛顿法为例对本发明的方法进行说明。
[0017]本实施例采用牛顿法获取回弹补偿系数的过程具体为：
[0018] A0.设定零件设计几何型面的各节点的权数Yi，其中，i的取值为丨，2,3，…，用 来对零件设计几何型面的各节点进行标识，仅包括工程师所关注的节点，由于零件设计几 何型面上的节点与零件冲压回弹仿真后几何型面上的节点、零件回弹补偿后几何型面的节 点均一一对应，即i即用于对零件冲压回弹仿真后几何型面上的节点及零件回弹补偿后几 何型面的节点进行标识；
[0019] A.以回弹补偿系数作为设计变量，以零件冲压回弹仿真后几何型面的节点分布与 零件设计几何型面的节点分布之间的差值作为目标函数，目标函数结合权数Yi建立；具体 为：在零件的不同位置设定相应的回弹补偿系数，将所有的回弹补偿系数作为设计变量，并 将其描述为一个矩阵 X，零件设计几何型面的各节点的坐标表示为（ΝΧκι，Ν}Γω，ΝΖ?。)，零件冲 压回弹仿真后几何型面的各节点的坐标表示为（ NXi(x)，Nyi(x)，NZi(X))，因而可得到目标 函数的具体表达式为： Nn
[0020] F ⑷= ￡^[(Νχ#)-Nxin)2+(Ny】⑷-NyJ2 +(NZi⑷，職0)2] (1) i=i
[0021] 其中，Nn为工程师关注的节点个数；
[0022] B·对零件设计几何型面进行网格划分，获得零件设计几何型面的节点分布；
[0023] C·运用有限元分析方法对零件设计几何型面的节点分布进行初始的冲压回弹仿 真分析，获取零件冲压回弹仿真后几何型面的节点分布；
[0024] D.设定回弹补偿系数的取值区域，并在该取值区域内随机生成回弹补偿系数取 值，设定零件冲压回弹仿真后几何型面与零件设计几何型面之间的最大允许误差；具体为： 在所设定的设计变量的取值区域内随机产生的回弹补偿系数取值x k,设计变量的取值区域 可由工程师根据经验进行设定，以使回弹补偿系数取值xk处于合适的范围内，减少计算量， 而设定零件冲压回弹仿真后几何型面与零件设计几何型面之间的最大允许误差则通过设 定目标函数梯度的终止标准ε实现，终止标准 ε的设定是为了对设计变量最后取值的精 度进行控制；
[0025] Ε·根据回弹补偿系数取值\、零件设计几何型面的节点分布和零件冲压回弹仿真 后几何型面的节点分布获得零件回弹补偿后几何型面的节点分布；
[0026] F.对零件回弹补偿后几何型面的节点分布进行冲压回弹仿真分析，获得新的零件 冲压回弹仿真后几何型面的节点分布；
[0027] G.判断零件冲压回弹仿真后几何型面的节点分布与零件设计几何型面的节点分 布之间的差值是否小于最大允许误差，若是，此时的回弹补偿系数取值为需获取的回弹补 偿系数，若否，则基于最优化方法对目标函数的设计变量进行迭代计算，获得新的回弹补偿 系数取值，再返回步骤Ε ;具体为：
[0028] 根据泰勒公式将目标函数在回弹补偿系数取值xk处做简化处理，其中，k的取值为 0，1，2,…，用来对迭代次数进行标识，得到
[0029] = F(xk >t ^) + - ^rk )T H (xk )(λγ - (2)
[0030] 其中，VF(jrk^ F(X)在xk点处的一阶导数矩阵，H(xk)为F(x)在xk点处的二阶 导数矩阵，由此可得目标函数的梯度为
[0031] VFpr) = ^ H(jy(x - j%) (3)
[0032] 判断是否满足1免_^)丨|<·^若是，即零件冲压回弹仿真后几何型面的节点分布 与零件设计几何型面的节点分布之间的差值小于最大允许误差，此时的xk即为需犹取的回 弹补偿系数，若否，则根据极值点处梯度为零的性质，将式（3)写成迭代公式的形式（在进 行迭代计算的过程中，根据生产需要，工程师可以根据实际需要让其中的一部分回弹补偿 系数保持不变），进而得到设计变量（即回弹补偿系数）的搜寻方向 [0033] ^ =xk ( 4 )
[0034] 根据式⑴计算和H(xk)，并进一步求得H(xk)的逆QKxjr1，在这里，由于 目标函数很难表达成具体的函数，可以采用差分进行求解，并将求得的￥1^》和[H (?) ] 1 代入式⑷，得到设计变量新的迭代点xk+1，Xk+1便是新的回弹补偿系数取值，使xk = xk+1，返 回步骤E。
[0035] 上列详细说明是针对本发明之一可行实施例的具体说明，该实施例并非用以限制 本发明的专利范围，凡未脱离本发明所为的等效实施或变更，均应包含于本案的专利范围 中。
【权利要求】
1. 一种基于最优化方法的回弹补偿系数获取方法，其特征在于：该方法包括以下步 骤： A. 以回弹补偿系数作为设计变量，以零件冲压回弹仿真后几何型面的节点分布与零件 设计几何型面的节点分布之间的差值作为目标函数； B. 对零件设计几何型面进行网格划分，获得零件设计几何型面的节点分布； C. 运用有限元分析方法对零件设计几何型面的节点分布进行初始的冲压回弹仿真分 析，获取零件冲压回弹仿真后几何型面的节点分布； D. 设定回弹补偿系数的取值区域，并在该取值区域内随机生成回弹补偿系数取值，设 定零件冲压回弹仿真后几何型面与零件设计几何型面之间的最大允许误差； E. 根据回弹补偿系数取值、零件设计几何型面的节点分布和零件冲压回弹仿真后几何 型面的节点分布获得零件回弹补偿后几何型面的节点分布； F. 对零件回弹补偿后几何型面的节点分布进行冲压回弹仿真分析，获得新的零件冲压 回弹仿真后几何型面的节点分布； G. 判断零件冲压回弹仿真后几何型面的节点分布与零件设计几何型面的节点分布之 间的差值是否小于最大允许误差，若是，此时的回弹补偿系数取值为需获取的回弹补偿系 数，若否，则基于最优化方法对目标函数的设计变量进行迭代计算，获得新的回弹补偿系数 取值，再返回步骤E。
2. 根据权利要求1所述的基于最优化方法的回弹补偿系数获取方法，其特征在于：所 述步骤A前还包括有步骤A0 :设定零件设计几何型面的各节点的权数；步骤A中的目标函 数结合步骤A0中的权数建立。
【文档编号】G06F17/50GK104298830SQ201410545323
【公开日】2015年1月21日 申请日期:2014年10月15日 优先权日:2014年10月15日
【发明者】马荣鑫, 刘红武, 丁桦 申请人:广州中国科学院工业技术研究院