一种曲线纤维复合结构设计多层级优化方法

技术领域

本发明属于复合结构设计优化领域，更具体地，涉及一种曲线纤维复合结构设计多层级优化方法。

背景技术

曲线纤维复合结构是一种变刚度的先进复合材料结构形式，具有比强度高、比刚度大等优点。和直线纤维复合结构相比，具有更优的力学性能。这种结构可以通过改变纤维铺放角度来改变材料性能，因此它的可设计性极强。随着自动纤维铺放技术的发展，其广泛运用于航空、航天等领域。

通常曲线纤维复合结构设计优化主要是优化纤维铺设角度。在曲线纤维复合结构优化方法的许多组成部分中，参数化格式和优化算法对解的质量有很大影响。参数化格式应保证纤维铺设角的空间连续性以便于加工制造；优化算法应能快速收敛。因此，为发展优化纤维铺设角度的方法，需仔细考虑参数化格式和优化算法。

一类优化纤维铺设角度的设计方法基于参数化格式。这种参数化格式基于无散向量场插值，由于其采用有旋无散场，所以可保证纤维铺设角度的空间连续性。在优化算法中，采用基于敏度的移动渐近线法，效果不错，但该算法仍有进一步发展的空间。然而基于无散向量场插值的参数化格式和优化算法的结合以及优化算法计算成本的降低还需得到进一步解决。

发明内容

本发明提供一种曲线纤维复合结构设计多层级优化方法，用以解决现有曲线纤维复合结构设计优化方法存在计算成本高的技术问题。

本发明解决上述技术问题的技术方案如下：一种曲线纤维复合结构设计多层级优化方法，包括：

S1、建立参数化层级，层数j＝1表示层级结构中的最粗层，j＝m表示层级结构中的最细层；并在第j层的复合结构设计域D内均匀设置一系列离散的场中心点P

S2、根据场中心点坐标P

S3、定义复合结构优化设计模型，设计变量为场中心点P

S4、基于第j层的所述展开系数

S5、利用所述初始展开系数

上述技术方案的基础上，本发明还可以做如下改进。

进一步，在所述k次更新α

其中，ε

进一步，所述无散向量场表示为φ(x

其中，‖x

进一步，所述力平衡方程为K

进一步，所述有限元分析的实现方式为：

建立第j层各单元刚度矩阵K

则所述目标c

基于所述有限元分析得到的位移场u

进一步，所述向量v

进一步，更新设计变量α

总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：

(1)本发明提供的基于无散向量场插值的曲线纤维复合结构设计多层级优化方法，其建立参数化层次，从最粗层到最细层；对于层次结构的每一层，在结构内部均匀布局一系列离散场中心点和样本点，赋值展开系数作为设计变量，通过样本点处的纤维向量，利用无散向量场插值构建一个连续全局函数来表达整个设计域的纤维向量；利用有限元分析建立刚度矩阵与设计变量的关系；通过基于敏度的优化算法更新设计变量达到结构柔度最小的目标；得到了较粗层优化问题的解，进而计算出邻近较细层的设计初始值；重复以上步骤，得到了曲线纤维复合结构的最优纤维角度空间连续变化布局。优化算法与基于无散向量场插值的参数化格式结合，在设计变量减少优化效率较高的同时，优化过程的计算成本降低。

(2)该优化方法可以很容易地与基于无散向量场插值的参数化格式结合。

(3)通过该优化方法所得结构的性能接近于通过无散向量场插值单层直接计算得到的结果，但是计算成本较低。

附图说明

图1为本发明实施例提供的一种基于无散向量场插值的曲线纤维复合结构设计多层级优化方法流程框图；

图2为本发明实施例提供的一个曲线纤维复合材料悬臂梁结构优化设计示意图；

图3为本发明实施例提供的基于无散向量场插值的曲线纤维复合结构设计多层级优化算法流程图；

图4为图2对应的示例参数化层次中第1、2、3层关于设计点处的纤维角度值优化结果；

图5为图2对应的示例采用单级优化算法(SLO)后设计点处的纤维角度值优化结果；

图6为图2对应的示例目标函数关于迭代次数的收敛历程曲线；

图7为图2对应的示例基于SLO法优化得到的关于迭代次数的收敛历程曲线。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

实施例一

一种基于无散向量场插值的曲线纤维复合结构设计多层级优化方法，如图1所示，其包括以下步骤：

(1)建立参数化层级，层数为m，求解第j层的优化问题，其中j＝1,2,...m。j＝1表示层级结构中的最粗层，j＝m表示层级结构中的最细层。

(2)在第j层的复合结构设计域D内均匀设置一系列离散的场中心点P

(3)根据场中心点坐标P

(4)定义复合材料结构优化设计问题，设计变量为场中心点P

findα

min.c(α

s.t.K

其中，f为外力向量,K

(5)进行有限元分析：建立第j层各单元刚度矩阵K

(6)基于位移场u

(7)基于灵敏度

其中ε

(8)基于步骤(7)得到的第j层优化得到的展开系数

(9)利用步骤(8)得到的初始设计值

优选的，步骤(3)中，无散向量场φ(x

其中，‖x

优选的，步骤(3)中，各单元中心点处的向量v

优选的，步骤(7)中，更新设计变量α

本实施例提供的一种基于无散向量场插值的曲线纤维复合结构设计多层级优化算法，满足以下两种要求：第一、通过该优化算法所得结构的性能接近于通过最速下降法得到的结果，计算成本较低；第二、该优化算法可以很容易地与基于无散向量场插值的参数化格式结合。主要思路如下：第一、建立参数化层次，其每一层对纤维角度分布都有明显的参数化。从层次结构的顶层到底层，无散向量场插值的设计点的数量和密度增加，。第二、在层次结构的每一层，均制定优化问题。从最粗层到最细层，通过求解较粗层的优化问题，得到较细层优化问题的初始设计。由于求解较粗层的优化问题在计算成本方面低于求解最细层的优化问题，所以相比直接在最细层进行优化的传统算法，多级优化算法能更快地收敛于最优解。

为了更好的说明本实施例方法，现给出如下示例：

请参阅图2，本实施方式以带有面内载荷的平面悬臂梁纤维增强结构柔度最小化的优化问题为例来解释本发明。在给定的2×1矩形设计域D内给定纤维角度初始布局，本例设定初始纤维角度均为90°，区域左侧边界施加位移约束，区域右边界中心点处分布面内载荷f＝1。对所述悬臂梁纤维增强结构进行纤维角度布局优化，使其刚度最大化。

请参阅图3的流程图，本实施方式中，所述基于无散向量场插值的曲线纤维复合结构设计多层级优化方法，其包括以下步骤：

步骤一，建立参数化层次。在第一层，有限单元均匀排布为12×6，场中心点均匀排布为13×7；在第二层，有限单元均匀排布为24×12，场中心点均匀排布为25×13；在第三层，有限单元均匀排布为40×20，场中心点均匀排布为41×21。

步骤二，在第1层的复合结构设计域D内均匀设置一系列离散的场中心点P

步骤三，定义复合材料结构优化设计问题，设计变量为设计点P

findα

min.c(α

s.t.K

其中，f为外力向量,K

步骤四，在设计域D内划分12×6的有限元网格，生成72个单元，在每个单元e(e＝1,2,3,…,72)上建立各单元刚度矩阵K

步骤五，计算目标函数c

步骤六，利用基于敏度的移动渐近线优化算法更新设计变量α

步骤七，利用在第1层的无散向量场插值计算得到第2层P

步骤八，利用步骤七中得到的第2层初始设计值，重复步骤三到步骤七，得到第3层的初始设计值，进而解得第3层的纤维角度值布局。

优化结果如下：优化后参数化层次的第1、2、3层设计点处的纤维角度值布局如图4，单级优化算法(SLO)分析结果如图5，最小柔度值为10.5814，图6反映了得到最小柔度值的迭代次数为24。作为对比，采用单级优化算法(SLO)迭代25次优化得到的结构有限元分析结果如图7，柔度值为9.5537。可以看出，本发明的优化算法相较于SLO法，能在较少的迭代次数得到更小的柔度值，减少了计算时间，降低了计算成本。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。