一种协同拓扑构型与纤维路径的柔性机构优化设计方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN112989648A(43)申请公布日2021.06.18(21)申请号202110154497.2(22)申请日2021.02.04(71)申请人西安理工大学地址710048陕西省西安市碑林区金花南路5号(72)发明人同新星元振毅高新勤魏锋涛李淑娟(74)专利代理机构西安弘理专利事务所61214代理人曾庆喜(51)Int.Cl.G06F30/23(2020.01)G06F111/04(2020.01)G06F111/06(2020.01)G06F113/26(2020.01)权利要求书3页说明书8页附图3页(54)发明名称一种协同拓扑构型与纤维路径的柔性机构优化设计方法(57)摘要本发明公开了一种协同拓扑构型与纤维路径的柔性机构优化设计方法，通过在柔性机构构型拓扑优化时，同时考虑连续纤维路径的优化，解决盲目选择纤维铺放路径的问题；通过采用多项式模拟纤维路径，间接优化多项式控制系数的设计思想，解决非凸优化及纤维路径难以加工的问题。该方法充分考虑了拓扑优化与复合材料设计优势，能从根本上改善柔性机构的变形和承载能力。CN112989648ACN112989648A权利要求书1/3页1.一种协同拓扑构型与纤维路径的柔性机构优化设计方法，其特征在于，步骤一：边界条件定义及初始化参数设定柔性机构的设计区域、边界条件、复合材料力学属性、初始纤维铺放路径，设置输入、输出的虚拟弹簧刚度、确定变量的过滤半径，将初始设计域和纤维路径进行有限元离散，设置全局收敛的移动渐近算法的初始变量；步骤二：纤维铺放参考路径的参数化描述引入多项式函数近似描述单层上纤维铺放参考路径，参数化多项式函数为：其中，Aω为第ω项多形式的系数，X为无量纲化的位置坐标，其取值范围是[‑1，1]区间，n为多项式的项数；步骤三：构造基于惩罚模型的曲线纤维铺放层合板本构关系通过离散曲线纤维路径，离散单元上可等效为直线纤维，以单元中心点处的曲线切线方向代表离散单元的纤维方向，将离散的单元等效为直线纤维铺放层合版，根据经典层合板理论，可得曲线纤维铺放层合板离散单元i的本构关系；步骤四：建立协同拓扑与纤维铺放路径的柔性机构优化模型利用单元相对密度和多项式系数分别描述结构拓扑形状和纤维铺放路径，以柔性机构系统互应变MSE为设计目标，体积为约束条件，得到协同拓扑与纤维铺放路径的柔性机构优化模型；步骤五：有限元位移场分析结合步骤三层合板本构关系和平面四边形单元，有限元分析中离散单元i的刚度矩阵ki表示为其中，B为几何矩阵，为曲线纤维铺放层合板离散单元i的本构关系，Ω为离散单元积分区域，通过对各个单元的刚度矩阵进行集成，可获得整体刚度矩阵，采用有限元编程技术可实现结构的位移场响应；步骤六：敏感性分析根据步骤四的优化模型，结合力平衡方程和伴随载荷法，推导可得到目标函数f0(x)对设计变量的灵敏度分别为体积约束函数f1(x)对设计变量的敏感度可直接求导为2CN112989648A权利要求书2/3页其中，f0(x)和f1(x)分别代表互应变能和体积约束函数，xi代表第i单元的相对密度，p为惩罚因子，qi和ui分别表示在虚拟载荷和实际载荷作用下单元节点处的位移矢量；步骤七：敏感度滤波针对优化模型中目标函数对单元密度的敏感度进行滤波，保证拓扑构型的光滑性，而对多项式系数的敏感度不需要滤波；步骤八：设计变量更新在获取柔性机构优化模型对设计变量敏感性的基础之上，利用全局收敛的移动渐近线方法同时更新单元密度设计变量和多项式系数变量，作为下一步迭代使用的初始值；步骤九：迭代终止及结果输出以相邻两次迭代设计变量的差值作为判定优化迭代的终止条件，当满足其终止条件时，优化过程收敛，迭代结束，输出最后的单元密度和多项式系数；若不满足终止条件则返回步骤二，继续优化迭代，直到满足终止条件为止。2.根据权利要求1所述的一种协同拓扑构型与纤维路径的柔性机构优化设计方法，其特征在于，步骤三曲线纤维铺放层合板离散单元i的本构关系为其中，hj为单层厚度，h是曲线纤维层合板总厚度，m为总层数，E1、E2、G12、μ12、μ21为材料j的弹性常数，Ri是第j层上i单元的转置矩阵，其表达式为其中，代表第j层上i单元中心处的曲线的斜率，A表示多项式系数矢量；基于惩罚模型，对称曲线纤维铺放层合板离散单元i的本构关系为式中，xi是第i单元上的相对密度，p为惩罚因子。3.根据权利要求2所述的一种协同拓扑构型与纤维路径的柔性机构优化设计方法，其特征在于，步骤四协同拓扑与纤维铺放路径的柔性机构优化模型为3CN112989648A权利要求书3/3页式中，x是设计变量矢量，代表单元相对密度和多项式系数，K为层合板整体刚度矩阵，F是外力矢量，L为虚拟载荷矢量，U和Q分别表示在F和L