一种惯性开关双稳态结构的稳健性优化方法技术

本发明专利技术涉及微机械结构鲁棒性优化领域，具体涉及一种惯性开关双稳态结构的稳健性优化方法，用于解决由于加工工艺技术的局限，惯性开关双稳态结构的尺寸通常会存在一定的不确定性，尤其双稳态结构中的细长梁结构，对加工工艺引起的结构参数和材料性能的不确定更为敏感的不足之处。该惯性开关双稳态结构的稳健性优化方法依据惯性开关的性能要求，确定其双稳态结构稳健性优化设计目标函数和约束条件，以及双稳态结构参数设计的不确定域及工艺制备偏差不确定集，建立基于区间理论和BP神经网络代理模型的两层嵌套遗传算法，实现双稳态结构设计变量不确定性的稳健性优化。构设计变量不确定性的稳健性优化。构设计变量不确定性的稳健性优化。

【技术实现步骤摘要】
一种惯性开关双稳态结构的稳健性优化方法


[0001]本专利技术涉及微机械结构鲁棒性优化领域，具体涉及一种惯性开关双稳态结构的稳健性优化方法。

技术介绍

[0002]惯性开关是利用惯性力实现开关功能的一类传感器，双稳态惯性开关的稳健性指惯性开关对其双稳态结构的参数、材料性能和噪声干扰等不确定性的不敏感性，进行稳健性优化的主要目的是在惯性开关双稳态结构的参数的设计范围内寻求最佳的设计值，降低惯性开关的性能对不确定性因素的敏感性。
[0003]由于加工工艺技术的局限，惯性开关双稳态结构的尺寸通常会存在一定的不确定性，尤其双稳态结构中的细长梁结构，对加工工艺引起的结构参数和材料性能的不确定更为敏感。
[0004]做好前期对结构参数和材料性能不确定因素的优化，可以较好的提高惯性开关的双稳态结构一次通过的成功率，从而缩短开发周期，同时可以提高商业批量生产时的可靠性和收益性。基于区间理论和BP神经网络代理模型的两层嵌套遗传算法，利用非概率稳健性优化可以方便地解决惯性开关器件不确定性信息不足的问题，对惯性开关进行非概率稳健性优化的研究具有重要的工程意义。

技术实现思路

[0005]本专利技术的目的是解决由于加工工艺技术的局限，惯性开关双稳态结构的尺寸通常会存在一定的不确定性，尤其双稳态结构中的细长梁结构，对加工工艺引起的结构参数和材料性能的不确定更为敏感的不足之处，而提供一种惯性开关双稳态结构的稳健性优化方法。
[0006]为了解决上述现有技术所存在的不足之处，本专利技术提供了如下技术解决方案：
[0007]一种惯性开关双稳态结构的稳健性优化方法，其特殊之处在于，包括以下步骤：
[0008]步骤1、依据惯性开关的闭合和反向解锁的性能要求，选取K组双稳态结构静态力学特性的有限元仿真非线性曲线数值，获取双稳态结构的设计变量组和静态力学曲线F
max
、F
min
的仿真数值结果；所述双稳态结构为细长梁结构，设计变量组为[θ，w1，l1，l2]，α为倾斜角、w1为细梁宽度、l1为细梁长度、l2为中间梁长度；F
max
为从第一稳态转变到第二稳态所需的临界力、F
min
为从第二稳态转变到第一稳态所需要的临界力；K≥50；
[0009]步骤2、建立BP神经网络，以设计变量组为输入层，F
max
、F
min
为输出层，利用BP神经网络训练得到双稳态结构的设计变量组与F
max
、F
min
之间关系的代理模型；
[0010]步骤3、基于区间理论、步骤2得到的代理模型，建立两层嵌套遗传算法，两层嵌套遗传算法分为外层遗传算法和内层遗传算法；
[0011]步骤4、在外层遗传算法中，对双稳态结构的设计变量组进行种群初始化，产生N个设计变量组的初始种群；N≥40；
[0012]步骤5、对每一个设计变量组，调用内层遗传算法进行优化；
[0013]对设计变量组的每个变量值增加不确定域，生成不确定变量及其上界约束、下界约束；通过罚函数方法，将不确定目标优化问题转换为确定性多目标优化；计算转换后的确定性多目标优化问题的目标函数F
maxI
和约束函数F
minI
，从而得到罚函数通过线性组合计算外层遗传算法的适应度函数值fitvalue，并将结果返回外层遗传算法；I为区间；
[0014]步骤6、在外层遗传算法中，适应度函数值fitvalue与外层遗传算法的设计变量组组合并作为最后一列，依据最后一列的适应度函数值fitvalue，再次进行选择、交叉和变异；若收敛或达到最大化迭代次数，则输出优化结果，否则外层遗传算法产生下一代种群，返回步骤5。
[0015]进一步地，所述步骤1具体为：
[0016]加载K组双稳态结构的设计变量组：倾斜角度θ、细梁宽度w1、细梁长度l1、中间梁长度l2，以及静态力学曲线的F
max
、F
min
的仿真数值结果，分别得到设计变量组的输入数组P和静态力学曲线的F
max
、F
min
的输出数组T；并通过mapminmax()函数对输入数组P和输出数组T进行归一化，得到归一化数组P
n
和T
n
，以及描述数据的结构体PS_P.mat和PS_T.mat。
[0017]进一步地，所述步骤2具体为：
[0018]步骤2.1、建立BP神经网络并设置网络创建参数；
[0019]步骤2.2、设置BP神经网络的训练参数，并调用train函数训练BP神经网络，得到BP神经网络的最终结构体net.mat，即代理模型；
[0020]步骤2.3、重新选取M组双稳态结构的设计变量组构成的数组P_test，利用net.mat、PS_P.mat和PS_T.mat，计算得到BP神经网络的计算结果F
max_net
、F
min_net
；M≥20；
[0021]并对所述M个设计变量组进行有限元仿真得到对应的有限元结果F
max_FEA
和F
min_FEA
；
[0022]利用均方差性能分析函数mse比较BP神经网络的计算结果F
max_net
、F
min_net
与有限元结果F
max_FEA
和F
min_FEA
，衡量BP神经网络的性能。
[0023]进一步地，步骤2.1中，所述网络创建参数包括隐含层节点数NodeNum、输出维数TypeNum、神经网络层传递函数TF，采用梯度下降自适应学习率训练函数traingdx，利用BP神经网络的创建函数newff，获得BP神经网络的初始结构体net0.mat；步骤2.2中，所述训练参数包括结果显示周期net.trainParam.show、最大训练次数net.trainParam.epochs、训练精度net.trainParam.goal、学习效率net.trainParam.lr。
[0024]进一步地，所述步骤4具体为：
[0025]步骤4.1、通过外层遗传算法加载步骤2.3生成的net.mat、PS_P.mat和PS_T.mat，设置外层遗传算法的交叉概率pc和变异概率pm；
[0026]步骤4.2、设置设计变量组的变量区间，变量区间为[θ
I
，w
1I
，l
1I
，l
2I
]；
[0027]步骤4.3、通过拉丁采样法在双稳态结构的变量区间内获得N个设计变量组的初始种群；N≥40；
[0028]步骤4.4、利用dec2bin函数将步骤4.3的初始种群转换为二进制编码，再利用bin2gray将二进制编码转换为格雷码。
[0029]进一步地，所述步骤5具体为：
[0030]步骤本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种惯性开关双稳态结构的稳健性优化方法，其特征在于，包括以下步骤：步骤1、依据惯性开关的闭合和反向解锁的性能要求，选取K组双稳态结构静态力学特性的有限元仿真非线性曲线数值，获取双稳态结构的设计变量组和静态力学曲线F
max
、F
min
的仿真数值结果；所述双稳态结构为细长梁结构，设计变量组为[θ，w1，l1，l2]，θ为倾斜角、w1为细梁宽度、l1为细梁长度、l2为中间梁长度；F
max
为从第一稳态转变到第二稳态所需的临界力、F
min
为从第二稳态转变到第一稳态所需要的临界力；K≥50；步骤2、建立BP神经网络，以设计变量组为输入层，F
max
、F
min
为输出层，利用BP神经网络训练得到双稳态结构的设计变量组与F
max
、F
min
之间关系的代理模型；步骤3、基于区间理论、步骤2得到的代理模型，建立两层嵌套遗传算法，两层嵌套遗传算法分为外层遗传算法和内层遗传算法；步骤4、在外层遗传算法中，对双稳态结构的设计变量组进行种群初始化，产生N个设计变量组的初始种群；N≥40；步骤5、对每一个设计变量组，调用内层遗传算法进行优化；对设计变量组的每个变量值增加不确定域，生成不确定变量及其上界约束、下界约束；通过罚函数方法，将不确定目标优化问题转换为确定性多目标优化；计算转换后的确定性多目标优化问题的目标函数F
maxI
和约束函数F
minI
，从而得到罚函数通过线性组合计算外层遗传算法的适应度函数值fitvalue，并将结果返回外层遗传算法；I为区间；步骤6、在外层遗传算法中，适应度函数值fitvalue与外层遗传算法的设计变量组组合并作为最后一列，依据最后一列的适应度函数值fitvalue，再次进行选择、交叉和变异；若收敛或达到最大化迭代次数，则输出优化结果，否则外层遗传算法产生下一代种群，返回步骤5。2.根据权利要求1所述的一种惯性开关双稳态结构的稳健性优化方法，其特征在于，所述步骤1具体为：加载K组双稳态结构的设计变量组：倾斜角度θ、细梁宽度w1、细梁长度l1、中间梁长度l2，以及静态力学曲线的F
max
、F
min
的仿真数值结果，分别得到设计变量组的输入数组P和静态力学曲线的F
max
、F
min
的输出数组T；并通过mapminmax()函数对输入数组P和输出数组T进行归一化，得到归一化数组P
n
和T
n
，以及描述数据的结构体PS_P.mat和PS_T.mat。3.根据权利要求2所述的一种惯性开关双稳态结构的稳健性优化方法，其特征在于，所述步骤2具体为：步骤2.1、建立BP神经网络并设置网络创建参数；步骤2.2、设置BP神经网络的训练参数，并调用train函数训练BP神经网络，得到BP神经网络的最终结构体net.mat，即代理模型；步骤2.3、重新选取M组双稳态结构的设计变量组构成的数组P_test，利用net.mat、PS_P.mat和PS_T.mat，计算得到BP神经网络的计算结果F
max_nef
、F
min_net
；M≥20；并对所述M个设计变量组进行有限元仿真得到对应的有限元结果F
max_FEA
和F
min_FEA
；利用均方差性能分析函数mse比较BP神经网络的计算结果F
max_net
、F
min_net
与有限元结果F
max_FEA
和F
min_FEA
，衡量BP神经网络的性能。4.根据权利要求3所述的一种惯性开关双稳态结构的稳健性优化方法，其特征在于：步骤2.1中，所述网络创建参数包括隐含层节点数NodeNum、输出维数TypeNum、神经网络层传
递函数TF，采用梯度下降自适应学习率训练函数traingdx，利用BP神经网络的创建函数newff，获得BP神经网络的初始结构体net0.mat；步骤2.2中，所述训练参数包括结果显示周期net.trainParam.show、最大训练次数net.trainParam.epochs、训练精度net.trainParam.goal、学习效率net.trainParam.lr。5.根据权利要求4所述的一种惯性开关双稳态结构的稳健性优化方法，其特征在于，所述步骤4具体为：步骤4.1、通过外层遗传算法加载步骤2.3生成的net.mat、PS_P.mat和PS_T.mat，设置外层遗传算法的交叉概率pc和变异概率pm；步骤4.2、设置设计变量组的变量区间，变量区间为[θ
I
，w
1I
，l
1I
，l
2I
]；步骤4.3、通过拉丁采样法在双稳态结构的变量区间内获得N个设计变量组的初始种群；步骤4.4、利用dec2bin函数将步骤4.3的初始种群转换为二进制编码，再利用bin2gray将二进制编码转换为格雷码。6.根据权利要求5所述的一种惯性开关双稳态结构的稳健性优化方法，其特征在于，所述步骤5具体为：步骤5.1、对每一个设计变量组，调用内层遗传算法，设置步骤1中的F
min
的约束条件值F
x
，保证|F
min
|≥|F
x
|，|F
x
|为双稳态结构的质量块(1)与开关底部电极接触时所要求的接触压力；步骤5.2、设置约束函数F
minI

【专利技术属性】
技术研发人员：王卫东，刘敏，朱应敏，李思新，王超，王锋，
申请(专利权)人：无锡必创数据科技有限公司，
类型：发明
国别省市：