CN101042712A - 实现面向大规模定制的产品族设计的方法 - Google Patents
实现面向大规模定制的产品族设计的方法 Download PDFInfo
- Publication number
- CN101042712A CN101042712A CN 200710039622 CN200710039622A CN101042712A CN 101042712 A CN101042712 A CN 101042712A CN 200710039622 CN200710039622 CN 200710039622 CN 200710039622 A CN200710039622 A CN 200710039622A CN 101042712 A CN101042712 A CN 101042712A
- Authority
- CN
- China
- Prior art keywords
- product
- design
- variable
- product family
- platform
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种实现面向大规模定制的产品族设计的方法,步骤为:确定决策变量及其各种约束条件,即确定出单个产品的表现型X和问题的解空间;根据确定的决策变量和约束条件建立优化模型;采用实数编码表示可行解的染色体编码;在可行解的编码方式确定后选择相应的解码方法,作为由单个产品基因型X到单个产品表现型X的对应关系或转换方法;确定单个产品适应度的量化评价方法,确定出由目标函数值f(X)到单个产品适应度F(X)的转换规则;确定出选择、交叉、变异遗传算子的具体操作方法;进行具体操作,确定产品的全部设计变量x的值,最终获得产品族。本发明改变了以往需要预先指定平台变量的不足,有效的支持大规模定制。
Description
技术领域
本发明涉及的是一种计算机应用技术领域的方法,具体涉及一种实现面向大规模定制的产品族设计的方法。
背景技术
大规模定制是一个集计算机技术、机械原理、计算方法等诸多技术于一体的复杂的系统工程。其中产品族设计是一项最为关键的技术,在大规模定制不断普及的今天,产品族的设计具有举足轻重的作用。人们进行了许多关于产品族设计方法的研究,得到了很多产品族设计方法。
经对现有技术的文献检索发现,王云霞等人在《机械设计》里撰文《面向大规模定制的产品族设计方法研究综述》(第21卷第6期,2004年6月)中指出:根据设计过程的不同产品族设计可分为top-down与bottom-down两种设计方法;根据产品族中产品特征的实现方式的不同可分为基于模块化的产品族设计和基于可调节变量的设计方法。基于可调节变量的设计方法主要通过改变一个或多个产品平台参数来形成产品族。产品平台参数主要指组成产品族的产品设计变量,包含公共变量和可调节变量。产品族中的产品特性通过调节变量的不同取值而变化。在可调节变量方法中,公共设计变量的确定是问题解决的关键和难点,也就是确定采用何种方式获取产品间的最大设计共性。传统的产品族设计方法需要预先指定平台变量,而这主要取决于产品族设计者的从业经验,当遇到设计者不熟悉的情况时,他所给出的产品族信息是不可靠的,这是传统方法缺陷之所在。
发明内容
本发明的目的在于克服现有技术中的不足,提供一种实现面向大规模定制的产品族设计的方法,使其可以自动改变平台共性并搜索共性与产品性能之间的最佳平衡点,在满足性能需求与约束的前提下获得合适的产品族,从而改变了以往需要预先指定平台变量的不足,有效的支持大规模定制。
本发明是通过以下技术方案实现的,步骤如下:
第一步:确定决策变量及其各种约束条件,即确定出单个产品的表现型X和问题的解空间。
第二步:根据第一步确定的决策变量和约束条件建立优化模型,即确定出目标函数的类型(是求目标函数的最大值还是求目标函数的最小值)及其数学描述形式或量化方法,也就是最终获取的产品族。
第三步:建立优化模型之后就要确定表示可行解的染色体编码方法,也即确定出单个产品的基因型X及遗传算法的搜索空间。
第四步:在可行解的编码方式确定之后选择相应的解码方法,即确定出由单个产品基因型X到单个产品表现型X的对应关系或转换方法。
第五步:确定单个产品适应度的量化评价方法,即确定出由目标函数值f(X)(也就是最终的产品族)到单个产品适应度F(X)的转换规则。
第六步:设计遗传算子,即确定出选择运算、交叉运算、变异运算等遗传算子的具体操作方法。
第七步:根据第六步所确定的选择、交叉、变异等具体操作方法,进行具体操作,进而确定产品的全部设计变量x的值,最终获得产品族。
以下对本发明方法作进一步的描述,具体实现过程如下:
1.1问题描述
产品族是共享公共特征、零部件与模块等通用技术并定位于一系列细分市场应用的一组产品,其中产品之间通用的这些参数、特征与零部件等称为产品平台,由产品平台衍生的一族个性化产品个体称为产品变量。
一个产品族可以通过一组设计变量表示如下:
x={x1,x2,…,xm}
式中xj——设计变量列向量,包括平台变量与个性化变量
j——设计变量
m——描述产品族的设计变量的总数量
对于产品变量i=1,2,…,n
xj={xj,1,xj,2,…,xj,n}′
式中xj,i——产品变量i的一个设计变量j的值
i——产品变量
n——产品族内产品总数量
1.2建立优化模型
产品族优化问题可用公式描述如下:
Fori=1,2,…,n
Find: xc,xv,i
Minimize: f(xc,xv,i)
Subject to:g(xc,xv,i)≤0
h(xc,xv,i)=0
xl≤x≤xu
式中c——平台设计变量的数量
v——个性化设计变量的数量
xc——平台变量的值
xv,i——个性化设计变量的值
1.3确定编码方式
实数编码适合在遗传算法中表示范围较大的数,满足精度要求较高的要求,便于搜索较大的产品族解空间,提高了产品族设计的效率,综合以上因素,在本发明中采用的是实数编码。
1.4适应度函数的确定
为了体现染色体的适应能力,引入了对问题中的每一个染色体都能进行度量的函数,叫适应度函数。通过适应度函数来决定染色体的优、劣程度,它体现了自然进化中的优利劣汰原则。对优化问题,适应度函数就是目标函数。
f(x)——适应度值;
g(x)——最大值问题下的适应度值。
在进行产品族设计的时候,往往需要实现多个目标,即为多目标优化问题。多目标最优化的数学描述如下:
1.5设定遗传算子
对遗传算子进行选择、交叉、变异等具体操作结合控制基因,并确定操作终止原则,以产生产品族的可行解空间。
1.6获得产品族
在可行解的解空间中,按照具体要求,获得产品族最优解。
本发明采用简单的交叉方法,即首先对群体中的个体进行随机配对,然后在配对个体中随机设定交叉点,配对个体彼此交换部分基因。交叉操作的示例如下:
为了尽可能挖掘可行解,以免搜索过早的陷于局部最优解,采用变异操作可取得较好的效果,既要尽量不破坏已搜索得到的较优解,又要充分考虑到多样性,应该根据实际情况确定变异率。在产品族优化设计过程中,对于产品的某个变量,如果允许变异,则在所规定的范围内,再随机产生一个赋值。
相对于现有其它方法,本发明优点如下:(1)引入了产品族惩罚函数(PFPF)解决多平台问题:以往的产品族设计方法在处理多平台问题时未能提出一种高效的方法,本发明引入了PFPF,PFPF在优化过程中确定哪个变量适合作平台变量及其最佳类数,不需要借用其他方法,在优化的过程中自动平衡平台化和定制化。产品族惩罚函数(Product Family Penalty Function,PFPF)针对未指定平台变量的情况,寻找对共性性能影响最大的变量,从而确定平台变量。
其中:
xij是第i个产品的第j个变量值,PFPF即是把每个变量的标准差之和。
PFPF越小则产品族共性越大,PFPF越大则共性越小。把PFPF加入到适应度函数中,来决定产品每个变量的最佳平台数,从而代替了由设计者指定平台变量的方法。
(2)转化处理强约束问题:在普通的产品族设计方法中,对于约束处理采取的方法一般是罚函数法。结合后面所举的实例,由于电机产品对扭矩T相比其他约束有着更高的要求,而且该产品族中每个产品所要求的扭矩T值都不一样。在研究过程中发现,如果仅仅采用罚函数来处理这些约束条件,发现只有很少一部分产品的扭矩T能满足条件。然而,在把该强约束(对T的要求)转化到目标函数中去,然后调整系数,发现每个产品的扭矩T都能很好的满足,而且其他约束条件及变量值都能令人满意。
附图说明
图1基于0/1型两层结构的产品族设计流程图
图2基于0/m型两层结构的产品族设计流程图(m=3为例)
图3适应度函数流程图
图4轮盘赌选择法操作流程图
图5交叉操作函数流程图
图6变异操作函数流程图
图7控制基因控制函数流程图(m=3为例)
图8每代最大适应度值示意图
图9每代最大适应度值示意图
图10每代最大适应度值示意图
具体实施方式
下面对本发明的实施例作详细说明:本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
2.1问题描述:
通用电机因其可在交流和直流电下运行而得名,具有很高的性能与适应性,目前得到了广泛应用,尤其在家用领域,例如,电钻、电锯与缝纫机。本实例引用于Simpson等人的论文,设计一个具有10个电机的产品族,满足不同的扭矩与功率性能。有关电机的构造请参考相关专业文献。在本例中,设计变量包括:转子与定子的绕线的横截面积及其匝数,电机半径、厚度与堆栈长度,发电电流,其符号与值域如下:
转子绕线匝数Nc,(100≤Nc≤1500匝)
场磁极绕线匝数,Ns(1≤Ns≤500匝)
转子绕线横截面积,Awa(0.01≤Awa≤1.0mm2)
定子绕线的横截面积,Awf(0.01≤Awf≤1.0mm2)
电机半径,ro(0.01≤ro≤0.10m)
定子厚度,t(0.0005≤t≤0.10m)
电流,I(0.1≤I≤6.0Amp)
堆栈长度,L(0.001≤L≤0.10m)
两个设计目标为:
效率(η,%)最大化,每个电机比较理想的目标效率是70%;
重量(M,kg)最小化,目标重量是0.5kg。
同时,遵循6个约束:
每个电机的扭矩,Ti(i=1,…,10):Ti={0.05,0.10,0.125,0.15,0.20,0.25,0.30,0.35,0.40,0.50}Nm;
每个电机的功率P=300W;
ro/t>1以确保几何结构可行;
每个电机的磁通强度,H<5000A·turns/m(安培·匝/米);
每个电机的重量小于2kg;
每个电机的效率大于15%。
2.2数学建模
①功率P
电机的功率为输入的功率减去铜线和电刷上损失的功率。由于机械和内核上的功率损失较少,因此忽略不计。
P=Pin-Ploss
其中:Pin=VtI
Ploss=Pcopper+Pbrush
有:Pcopper=I2(Ra+Rs)
Pbrush=2I
其中:Ra=ρ[2L+4(ro-t-lgap)]Nc/Awa
Rs=ρ2[2L+4(ro-t)]Ns/Awf
这里常数ρ是铜线的电阻系数,ρ=1.69×10-8Ω·m.
②效率η
效率η是输出功率和输入功率的比
③质量M
电机的质量由三部分组成:定子,电枢,以及电场内和电枢上的线圈
M=M1+M2+M3
其中:M1=πL[ro 2-(ro-t)2]ρsteel
M2=πL(ro-t-lgap)2ρsteel
M3=ρcopper{[2L+4(ro-t-lgap)]NcAwa+2[2L+4(ro-t)]NsAwf}
常数ρcopper=8960,ρcopper=7850,且π=3.14159.
④扭矩T
最后,电机扭矩是由常数K,磁通量φ和电流I组成。
T=KφI
其中:
其中:lc=π(2ro+t)/2
lr=2(ro-t-lgap)
As=tL
Ar=lrL
Aa=Ar
μ由磁化强度曲线得到,该曲线分为三个区域,为了更好的适应这个模型,每个部分都有一个适当的公式,如下:
μsteel=-0.22791·H2+52.411·H+3115.8,H≤220
μsteel=11633.5-1486.33·ln(H),220<H≤1000
μsteel=1000,H>1000
自由空间渗透性μ0是一个常数,μ0=4π×10-7(unit:Henrys/m)。
到此,电机优化问题的数学建模完成。
3.1参数设置
经过多次试验,对比所得到的数据结果,不断尝试,以及考虑遗传算法的一般要求,设定运行参数如下:
群体大小POPSIZE:200
遗传代数MAXGEN:100或200(根据不同情况具体确定)
交叉概率Pc:0.8
变异概率Pm:0.05
3.2适应度函数的设定
(1)目标转化
由于在轮盘赌选择法(图4)中,适应度值大的个体被选择进入下一代的几率也大,因此把目标都转化为使目标最大化的情况。同时考虑以下几点要求:
1)目标函数和罚函数应该在数量级上相差不大,要具有可比性,从而使罚函数发挥作用
2)∑Mi与∑ηi不能简单的以权重分配,应该具有一定的独立性
3)为了目标函数和罚函数的可比性,因此把该值归一化,因此,目标函数设为:(∑Mi/∑ηi)/(1+∑Mi/∑ηi)
(2)约束设定
总共有6个约束条件,用数组Cstr[6]来表示。凡是有不满足约束的情况,都要处以一定的罚函数。总的规则是,以超过界限的百分比来表示不满足程度,具体的惩罚系数要根据试运行的结果进行调整,以使每个约束都尽量满足。
①
Ti={0.05,0.10,0.125,0.15,0.20,0.25,0.30,0.35,0.40,0.50};
②
③如果Hi-5000>0,则
④因为根据单独设计和两阶段法已经得到质量M应该都小于1kg,因此可以缩小范围空间,以提高搜索效率和最终的优化结果。
如果Mi-5000>0,则
⑤根据单独设计和两阶段法,发现较优解的效率都大于40%,因此可以加强约束,以进一步提高最终效果。
如果0.4-ηi>0,则
⑥如果1-ri/ti>0,则
其中:i=1,2,…,10
(3)适应度函数确定(图3)
对于电机产品,一个重要的指标就是扭矩T必须满足条件,而功率P只要在规定的误差内即可。在本案例中,产品族中的10个产品扭矩都不同,因此如果只用罚函数处理扭矩这个约束,在运行程序的过程中发现所得结果中扭矩T大多都不能满足要求,如果仅仅调整该项的惩罚系数,仍然不能确保每个产品的T约束都能很好。为解决这个问题,创造性的把该项约束当作目标函数处理,即以1/Cstr[0]表示,发现最终效果非常好。把目标函数和罚函数结合起来,可以得到适应度函数。最后通过不断各项的系数,能得到满意的结果。
3.3输入设置
输入分为文件输入和屏幕提示输入。文件输入为产品族中8个变量各自的取值范围。由于已经由单独设计和两阶段设计法得到较优的结果,因此可以根据上述结果适当减小变量取值范围,这样可以提高本方法的搜索效率,无需耗费大量机器运行时间。
屏幕提示输入分为两种:
①基于0/1型两层结构遗传算法的产品族设计方法,运行时产生如下提示是否指定基因,若输入’y’,则提示输入’0’或’1’,回车后即可执行;若输入’n’,则随机产生控制基因:
Do you want to designate controlling genes?<y or n>y
Please input the 8 controlling genes<0 or 1>:
②基于0/m型两层结构遗传算法的产品族设计方法,本例中m=3,因此运行时产生如下提示是否指定基因,若输入’y’,则提示输入[0,3]内的整数,回车即可执行;若输入’n’,则随机产生控制基因:
Do you want to designate controlling genes?<y or n>y
Please input the 8 controlling genes<integer between 0 and 3>:
3.4输出设置
输出以下几项:
①每代中最大适应度值,以显示不断优化的过程;
②最终最优个体中10个产品的8个变量值Nc,Ns,Awa,Awf,t,l,I,r,以及6个约束结果T,P,H,M,r/t;
③输出群体中每个个体的约束结果,以作整体上的观察比较。
4设计问题模型构造及优化结果
4.10/1型2层结构GA用于单平台产品族优化(预先指定平台变量)
这种情况预先指定了平台变量,一次性求解整个产品族。预先指定的平台变量是根据观察单独设计和两阶段设计的结果由设计者依据设计经验指定的。此案例中,预先指定的平台变量为:Nc,Ns,Awa,t,l,即控制基因为(11100101)。
数学模型为:
求解:xc={Nc,Ns,Awa,t,l},xv,i={Awf,r,I}-xc为平台变量
目标:Min f1(x)=∑Mi (尽量以0.5kg为目标)
Max f2(x)=∑ηi
约束:Ti(x)={0.05,0.10,0.125,0.15,0.20,0.25,0.30,0.35,0.40,0.50}Nm
Pi(x)=300W
Hi(x)<5000A·turns/m
Mi(x)≤2.0kg
ηi(x)≥15%
ri/ti>1
其中:i=1,2,…,10.
4.1.1适应度函数
经过不断的调整系数进行试运行,以得到满意的结果,确定适应度函数为:
Fitness=9/Cstr[0]+(∑Mi/∑ηi)/(1+∑Mi/∑ηi)
-(120*Cstr[1]+Cstr[2]+Cstr[3]+Cstr[4]+Cstr[5])/5
4.1.2输出结果(运行200代)
4.1.2.1每代最大适应度值,如下:
代数 | 最大适应度值 | 代数 | 最大适应度值 | 代数 | 最大适应度值 | 代数 | 最大适应度值 |
1 | 22.12 | 51 | 134.99 | 101 | 246.88 | 151 | 363.61 |
2 | 24.20 | 52 | 135.34 | 102 | 249.17 | 152 | 363.61 |
3 | 33.01 | 53 | 144.01 | 103 | 249.17 | 153 | 363.61 |
4 | 33.01 | 54 | 144.01 | 104 | 285.00 | 154 | 363.61 |
5 | 33.01 | 55 | 151.63 | 105 | 285.00 | 155 | 363.62 |
6 | 33.02 | 56 | 166.58 | 106 | 285.00 | 156 | 363.62 |
7 | 39.33 | 57 | 166.58 | 107 | 285.00 | 157 | 363.63 |
8 | 40.84 | 58 | 166.58 | 108 | 285.00 | 158 | 363.63 |
9 | 40.84 | 59 | 166.59 | 109 | 285.04 | 159 | 363.63 |
10 | 40.84 | 60 | 166.62 | 110 | 285.05 | 160 | 363.65 |
11 | 44.36 | 61 | 166.62 | 111 | 285.07 | 161 | 363.65 |
12 | 44.36 | 62 | 166.93 | 112 | 285.07 | 162 | 366.64 |
13 | 58.19 | 63 | 177.54 | 113 | 291.55 | 163 | 366.64 |
14 | 58.20 | 64 | 177.54 | 114 | 297.87 | 164 | 366.64 |
15 | 58.20 | 65 | 187.35 | 115 | 307.74 | 165 | 366.64 |
16 | 58.20 | 66 | 187.35 | 116 | 307.74 | 166 | 366.66 |
17 | 58.20 | 67 | 199.90 | 117 | 307.74 | 167 | 366.66 |
18 | 58.21 | 68 | 199.90 | 118 | 307.74 | 168 | 366.66 |
19 | 59.17 | 69 | 199.96 | 119 | 309.18 | 169 | 378.63 |
20 | 65.19 | 70 | 205.25 | 120 | 309.19 | 170 | 378.63 |
21 | 65.19 | 71 | 205.25 | 121 | 309.19 | 171 | 378.63 |
22 | 65.19 | 72 | 205.25 | 122 | 309.19 | 172 | 405.97 |
23 | 65.19 | 73 | 205.25 | 123 | 319.75 | 173 | 405.97 |
24 | 65.19 | 74 | 207.46 | 124 | 319.75 | 174 | 405.97 |
25 | 66.80 | 75 | 207.48 | 125 | 319.76 | 175 | 405.99 |
26 | 67.27 | 76 | 207.48 | 126 | 319.77 | 176 | 405.99 |
27 | 68.94 | 77 | 207.48 | 127 | 319.77 | 177 | 405.99 |
28 | 74.66 | 78 | 210.33 | 128 | 319.77 | 178 | 405.99 |
29 | 74.66 | 79 | 210.33 | 129 | 319.80 | 179 | 405.99 |
30 | 82.00 | 80 | 210.72 | 130 | 328.41 | 180 | 406.00 |
31 | 86.05 | 81 | 214.20 | 131 | 328.41 | 181 | 406.02 |
32 | 86.05 | 82 | 214.21 | 132 | 328.41 | 182 | 406.02 |
33 | 87.23 | 83 | 214.21 | 133 | 328.41 | 183 | 406.02 |
34 | 87.23 | 84 | 217.00 | 134 | 328.41 | 184 | 406.02 |
35 | 87.23 | 85 | 217.00 | 135 | 328.41 | 185 | 406.02 |
36 | 108.82 | 86 | 217.00 | 136 | 328.41 | 186 | 406.02 |
37 | 108.82 | 87 | 221.61 | 137 | 343.29 | 187 | 406.90 |
38 | 108.82 | 88 | 237.00 | 138 | 343.29 | 188 | 406.90 |
39 | 108.82 | 89 | 246.87 | 139 | 343.29 | 189 | 406.90 |
40 | 108.82 | 90 | 246.87 | 140 | 343.29 | 190 | 406.90 |
41 | 108.82 | 91 | 246.87 | 141 | 343.29 | 191 | 406.90 |
42 | 108.95 | 92 | 246.88 | 142 | 343.34 | 192 | 406.90 |
43 | 112.77 | 93 | 246.88 | 143 | 343.34 | 193 | 406.90 |
44 | 112.96 | 94 | 246.88 | 144 | 343.34 | 194 | 406.90 |
45 | 127.76 | 95 | 246.88 | 145 | 343.34 | 195 | 406.90 |
46 | 127.76 | 96 | 246.88 | 146 | 343.36 | 196 | 406.91 |
47 | 127.96 | 97 | 246.88 | 147 | 343.36 | 197 | 414.51 |
48 | 127.96 | 98 | 246.88 | 148 | 343.38 | 198 | 414.51 |
49 | 134.99 | 99 | 246.88 | 149 | 346.93 | 199 | 414.51 |
50 | 134.99 | 100 | 246.88 | 150 | 363.61 | 200 | 414.52 |
通过观察上表,可以发现在程序执行的每代过程中,每代最大适应度的值是不断增加了,表现在不断寻找着最优点。在200时终止了程序,此时观察所输出的结果,发现各项都能很好的满足,因此选择此项输出作为最终结果。
4.1.2.2每个产品的8个变量值,如下:
控制基因值 | 1 | 1 | 1 | 0 | 0 | 1 | 0 | 1 |
电机号 | 变量 | |||||||
Nc | Ns | Awa(m2) | Awf(m2) | r(m) | t(m) | I(A) | l(m) | |
1 | 985 | 93 | 0.00000026 | 0.00000023 | 0.02008 | 0.011993 | 3.01 | 0.01702 |
2 | 985 | 93 | 0.00000026 | 0.00000022 | 0.02482 | 0.011993 | 3.162 | 0.01702 |
3 | 985 | 93 | 0.00000026 | 0.00000023 | 0.03054 | 0.011993 | 3 | 0.01702 |
4 | 985 | 93 | 0.00000026 | 0.00000024 | 0.03186 | 0.011993 | 3.2 | 0.01702 |
5 | 985 | 93 | 0.00000026 | 0.00000025 | 0.0356 | 0.011993 | 3.448 | 0.01702 |
6 | 985 | 93 | 0.00000026 | 0.00000021 | 0.03312 | 0.011993 | 4.024 | 0.01702 |
7 | 985 | 93 | 0.00000026 | 0.00000023 | 0.03558 | 0.011993 | 4.22 | 0.01702 |
8 | 985 | 93 | 0.00000026 | 0.00000024 | 0.03542 | 0.011993 | 4.602 | 0.01702 |
9 | 985 | 93 | 0.00000026 | 0.00000023 | 0.03646 | 0.011993 | 4.814 | 0.01702 |
10 | 985 | 93 | 0.00000026 | 0.00000025 | 0.03998 | 0.011993 | 4.996 | 0.01702 |
4.1.2.3每个产品的6个约束值,即0/1型指定平台变量性能参数,如下表:
电机 | 性能参数 |
号 | T(Nm) | P(W) | H(A turns/m) | M(kg) | η | r/t |
1 | 0.05697221 | 294.3432617 | 2853.640625 | 0.33339331 | 0.85033447 | 1.67430999 |
2 | 0.09935479 | 291.2606201 | 2401.187744 | 0.46839705 | 0.80098069 | 2.06954051 |
3 | 0.12470806 | 263.4855957 | 1837.015381 | 0.65877485 | 0.76372635 | 2.54648529 |
4 | 0.15033506 | 272.6925049 | 1875.665527 | 0.7101301 | 0.74101222 | 2.6565494 |
5 | 0.20035274 | 272.8192139 | 1802.564209 | 0.85444969 | 0.68803394 | 2.96839809 |
6 | 0.25001302 | 303.1043091 | 2266.121338 | 0.74950254 | 0.65499246 | 2.76161078 |
7 | 0.29991841 | 299.8554077 | 2207.43042 | 0.84997034 | 0.61787641 | 2.96673036 |
8 | 0.35481173 | 312.1981506 | 2418.440674 | 0.84585464 | 0.58991009 | 2.95338945 |
9 | 0.40131646 | 307.5742798 | 2455.650635 | 0.88579589 | 0.55557936 | 3.04010674 |
10 | 0.47639486 | 284.9743347 | 2318.347656 | 1.03716493 | 0.49600437 | 3.33361103 |
从上表可以看出,每个产品的强约束扭矩T能够很好的满足条件,弱约束P值在误差允许的范围内,其他约束全部满足。对于0/1型预先指定平台变量的这种情况,此组数据比较理想。
4.2 0/1型2层结构GA用于单平台产品族优化(不指定平台变量)
对于未指定平台变量的情况,采用一种共性指数,即产品族惩罚函数(Product Family Penalty Function,PFPF),来寻找对共性性能影响最大的变量,从而确定平台变量。
其中:
xij是第i个产品的第j个变量值,PFPF即是把每个变量的标准差之和。
PFPF越小则产品族共性越大,PFPF越大则共性越小。把PFPF加入到适应度函数中,来决定产品每个变量的最佳平台数,从而代替了由设计者指定平台变量的方法。
该情况下电机产品族设计的数学模型为:
求解:x={Nc,Ns,Awa,Awf,r,t,I,l},xcc={xcc,j}-xcc为平台变量的控制基因值
目标:Min f1(x)=∑Mi (尽量以0.5kg为目标)
Max f2(x)=∑ηi
Min f3(x)=PFPF
约束:Ti(x)={0.05,0.10,0.125,0.15,0.20,0.25,0.30,0.35,0.40,0.50}Nm
Pi(x)=300W
Hi(x)<5000A·turns/m
Mi(x)≤2.0kg
ηi(x)≥15%
ri/ti>1
xcc,j={0,1}
其中:i=1,2,…,10.
4.2.1适应度函数
Fitness=2.5/Cstr[0]+(∑Mi/∑ηi)/(1+∑Mi/∑ηi)-PFPF
-(120*Cstr[1]+Cstr[2]+50*Cstr[3]+Cstr[4]+Cstr[5])/5
4.2.2输出结果(运行200代)
4.2.2.1每代最大适应度值,即0/1型不指定平台变量每代最大适应度值,如下表:
代数 | 最大适应度值 | 代数 | 最大适应度值 | 代数 | 最大适应度值 | 代数 | 最大适应度值 |
1 | 5.54 | 51 | 41.79 | 101 | 562.11 | 151 | 823.75 |
2 | 5.60 | 52 | 43.96 | 102 | 562.11 | 152 | 823.75 |
3 | 5.70 | 53 | 58.10 | 103 | 562.11 | 153 | 823.75 |
4 | 5.83 | 54 | 104.70 | 104 | 562.11 | 154 | 823.75 |
5 | 5.94 | 55 | 104.70 | 105 | 562.11 | 155 | 823.75 |
6 | 6.04 | 56 | 104.70 | 106 | 562.11 | 156 | 823.75 |
7 | 6.24 | 57 | 104.70 | 107 | 562.11 | 157 | 833.03 |
8 | 6.40 | 58 | 110.71 | 108 | 562.11 | 158 | 952.98 |
9 | 6.95 | 59 | 115.34 | 109 | 562.16 | 159 | 952.98 |
10 | 6.95 | 60 | 115.34 | 110 | 562.21 | 160 | 952.98 |
11 | 7.00 | 61 | 115.36 | 111 | 562.21 | 161 | 952.98 |
12 | 7.42 | 62 | 151.17 | 112 | 562.21 | 162 | 952.98 |
13 | 7.42 | 63 | 151.17 | 113 | 562.21 | 163 | 952.98 |
14 | 7.42 | 64 | 151.17 | 114 | 562.21 | 164 | 952.98 |
15 | 7.42 | 65 | 153.28 | 115 | 562.21 | 165 | 952.98 |
16 | 7.48 | 66 | 171.77 | 116 | 562.21 | 166 | 952.98 |
17 | 7.48 | 67 | 171.77 | 117 | 562.21 | 167 | 952.98 |
18 | 7.48 | 68 | 171.77 | 118 | 562.21 | 168 | 952.98 |
19 | 7.48 | 69 | 171.77 | 119 | 562.21 | 169 | 952.98 |
20 | 7.85 | 70 | 171.80 | 120 | 562.21 | 170 | 952.98 |
21 | 7.85 | 71 | 171.80 | 121 | 562.21 | 171 | 952.98 |
22 | 8.86 | 72 | 180.84 | 122 | 562.21 | 172 | 952.98 |
23 | 11.11 | 73 | 180.84 | 123 | 562.34 | 173 | 952.98 |
24 | 11.11 | 74 | 180.84 | 124 | 562.34 | 174 | 952.98 |
25 | 11.11 | 75 | 180.84 | 125 | 606.81 | 175 | 952.98 |
26 | 12.59 | 76 | 180.84 | 126 | 606.81 | 176 | 952.98 |
27 | 12.59 | 77 | 180.84 | 127 | 606.81 | 177 | 952.98 |
28 | 12.59 | 78 | 181.39 | 128 | 750.75 | 178 | 952.98 |
29 | 12.78 | 79 | 355.55 | 129 | 750.75 | 179 | 952.98 |
30 | 12.78 | 80 | 355.55 | 130 | 750.75 | 180 | 952.98 |
31 | 12.78 | 81 | 355.55 | 131 | 750.75 | 181 | 952.98 |
32 | 14.38 | 82 | 355.55 | 132 | 750.75 | 182 | 952.98 |
33 | 14.38 | 83 | 355.57 | 133 | 750.75 | 183 | 952.98 |
34 | 14.38 | 84 | 355.57 | 134 | 750.75 | 184 | 952.98 |
35 | 14.38 | 85 | 438.47 | 135 | 750.92 | 185 | 952.98 |
36 | 14.38 | 86 | 518.23 | 136 | 750.92 | 186 | 952.98 |
37 | 20.09 | 87 | 518.23 | 137 | 823.73 | 187 | 952.98 |
38 | 20.47 | 88 | 518.23 | 138 | 823.73 | 188 | 952.98 |
39 | 23.86 | 89 | 526.74 | 139 | 823.73 | 189 | 952.98 |
40 | 23.86 | 90 | 526.84 | 140 | 823.73 | 190 | 952.98 |
41 | 23.86 | 91 | 526.85 | 141 | 823.73 | 191 | 952.98 |
42 | 23.86 | 92 | 526.85 | 142 | 823.73 | 192 | 952.98 |
43 | 24.93 | 93 | 526.85 | 143 | 823.74 | 193 | 952.98 |
44 | 36.96 | 94 | 526.96 | 144 | 823.75 | 194 | 952.98 |
45 | 41.45 | 95 | 561.99 | 145 | 823.75 | 195 | 952.98 |
46 | 41.45 | 96 | 561.99 | 146 | 823.75 | 196 | 952.98 |
47 | 41.79 | 97 | 562.01 | 147 | 823.75 | 197 | 952.98 |
48 | 41.79 | 98 | 562.01 | 148 | 823.75 | 198 | 952.98 |
49 | 41.79 | 99 | 562.11 | 149 | 823.75 | 199 | 952.98 |
50 | 41.79 | 100 | 562.11 | 150 | 823.75 | 200 | 952.98 |
由图8看出,适应度值在不断增大,虽然在200代前似乎已经找到最优点,但通过观察前面的曲线,发现有多处平坦的区域,表明在这些代内并没有找到更优的点,但这并不表明已经得到最优点了,因此在200代处最后获得的解不一定是最优点。但此时观察输出结果,发现各个变量和约束值都已很好的满足条件,这已足够说明方法的优点。
4.2.2.2每个产品的8个变量值,即0/1型不指定平台变量设计结果,如下表:
控制基因值 | 0 | 0 | 1 | 0 | 1 | 0 | 0 | 0 |
电机号 | 变量 | |||||||
Nc | Ns | Awa(m2) | Awf(m2) | r(m) | t(m) | I(A) | l(m) | |
1 | 803 | 41 | 0.00000021 | 0.00000022 | 0.02506 | 0.005133 | 3.498 | 0.0142 |
2 | 803 | 44 | 0.00000021 | 0.00000022 | 0.02506 | 0.005133 | 3.77 | 0.02281 |
3 | 788 | 61 | 0.00000021 | 0.00000022 | 0.02506 | 0.005133 | 3.746 | 0.02113 |
4 | 803 | 65 | 0.00000021 | 0.00000022 | 0.02506 | 0.005133 | 3.412 | 0.0283 |
5 | 803 | 100 | 0.00000021 | 0.00000022 | 0.02506 | 0.005133 | 3.65 | 0.02113 |
6 | 803 | 65 | 0.00000021 | 0.00000022 | 0.02506 | 0.008983 | 3.992 | 0.03883 |
7 | 803 | 73 | 0.00000021 | 0.00000022 | 0.02506 | 0.00549 | 4.862 | 0.02491 |
8 | 803 | 99 | 0.00000021 | 0.00000022 | 0.02506 | 0.006372 | 4.684 | 0.0235 |
9 | 803 | 88 | 0.00000021 | 0.00000022 | 0.02506 | 0.005133 | 4.732 | 0.02902 |
10 | 807 | 87 | 0.00000021 | 0.00000024 | 0.02506 | 0.005133 | 4.562 | 0.03925 |
4.2.2.3每个产品的6个约束值,即0/1型不指定平台变量性能参数,如下表:
电机号 | 性能参数 | |||||
T(Nm) | P(W) | H(A turns/m) | M(kg) | η | r/t | |
1 | 0.04989308 | 304.5837402 | 1132.439209 | 0.38869148 | 0.75716245 | 4.88213543 |
2 | 0.09990544 | 302.6721802 | 1309.800781 | 0.54681307 | 0.69812524 | 4.88213543 |
3 | 0.12430941 | 302.4458008 | 1804.300415 | 0.52110088 | 0.70207244 | 4.88213543 |
4 | 0.14998446 | 270.4446716 | 1751.191406 | 0.65824634 | 0.68924171 | 4.88213543 |
5 | 0.19715768 | 286.2072144 | 2882.067139 | 0.54350138 | 0.68185163 | 4.88213543 |
6 | 0.24993716 | 286.8616943 | 2075.960205 | 0.82947874 | 0.62486207 | 2.78971384 |
7 | 0.29981852 | 326.3250122 | 2805.916992 | 0.59778273 | 0.58362997 | 4.56466309 |
8 | 0.34973904 | 322.1720581 | 3676.97583 | 0.57870787 | 0.59809911 | 3.93283106 |
9 | 0.4004958 | 300.1734619 | 3288.051514 | 0.68413937 | 0.55160695 | 4.88213543 |
10 | 0.50021398 | 266.5535278 | 3133.904297 | 0.87798488 | 0.50807905 | 4.88213543 |
强约束扭矩T非常满足,弱约束P在允许的误差范围内,其他约束条件都很好的满足,因此这组数据比较令人满意。
4.30/m型2层结构GA用于多平台产品族优化(预先指定平台变量)
预先指定的平台变量为:Nc,Ns,Awa,t,l,根据单独设计和两阶段设计法结果设定最大类数为3,即m=3,并指定控制基因为(23100203)。
对于m取3的补充说明:
根据单独设计和两阶段设计法结果设定最大类数为3,即m=3。原因出于如下两个考虑:(1)在具体实例中产品数量为10,若分为太多类则聚类的效果就不明显了;(2)m=3足以说明本方法的各种操作机制,因为m不为3时机制跟m=3时的机制是相同的。
数学模型为:
求解:xc={Nc,1,Nc,2;Ns,1,Ns,2,Ns,3;Awa;t1,t2;l1,l2,l3},xv,i={Awf,r,I}-xc为平台变量
目标:Min f1(x)=∑Mi (尽量以0.5kg为目标)
Max f2(x)=∑ηi
约束:Ti(x)={0.05,0.10,0.125,0.15,0.20,0.25,0.30,0.35,0.40,0.50}Nm
Pi(x)=300W
Hi(x)<5000A·turns/m
Mi(x)≤2.0kg
ηi(x)≥15%
ri/ti>1
其中:i=1,2,…,10.
4.3.1适应度函数
Fitness=3/Cstr[0]+2*(∑Mi/∑ηi)/(1+∑Mi/∑ηi)
-(250*Cstr[1]+Cstr[2]+200*Cstr[3]+Cstr[4]+Cstr[5])/5
4.3.2输出结果(运行100代)
4.3.2.1每代最大适应度值,即0/m型指定平台每代最大适应度值(m=3),
代数 | 最大适应度值 | 代数 | 最大适应度值 | 代数 | 最大适应度值 | 代数 | 最大适应度值 |
1 | 5.66 | 26 | 49.89 | 51 | 75.26 | 76 | 247.69 |
2 | 6.41 | 27 | 51.54 | 52 | 93.29 | 77 | 247.69 |
3 | 6.58 | 28 | 51.54 | 53 | 93.29 | 78 | 261.10 |
4 | 9.31 | 29 | 54.38 | 54 | 93.29 | 79 | 267.85 |
5 | 9.80 | 30 | 54.38 | 55 | 93.29 | 80 | 267.85 |
6 | 11.13 | 31 | 54.38 | 56 | 93.29 | 81 | 267.85 |
7 | 11.13 | 32 | 54.38 | 57 | 96.78 | 82 | 269.48 |
8 | 11.13 | 33 | 54.38 | 58 | 96.78 | 83 | 269.48 |
9 | 11.24 | 34 | 54.38 | 59 | 96.78 | 84 | 269.48 |
10 | 17.15 | 35 | 54.38 | 60 | 96.78 | 85 | 269.48 |
11 | 17.15 | 36 | 61.68 | 61 | 96.78 | 86 | 269.48 |
12 | 18.22 | 37 | 61.68 | 62 | 99.37 | 87 | 269.48 |
13 | 28.08 | 38 | 61.68 | 63 | 99.37 | 88 | 269.48 |
14 | 28.99 | 39 | 61.68 | 64 | 106.93 | 89 | 269.48 |
15 | 28.99 | 40 | 67.07 | 65 | 107.00 | 90 | 269.48 |
16 | 29.45 | 41 | 67.07 | 66 | 196.43 | 91 | 269.48 |
17 | 34.31 | 42 | 72.17 | 67 | 196.43 | 92 | 269.48 |
18 | 34.31 | 43 | 72.17 | 68 | 196.43 | 93 | 269.48 |
19 | 39.62 | 44 | 72.24 | 69 | 196.43 | 94 | 269.48 |
20 | 39.62 | 45 | 72.24 | 70 | 218.70 | 95 | 269.48 |
21 | 39.62 | 46 | 72.24 | 71 | 245.72 | 96 | 269.48 |
22 | 39.62 | 47 | 72.24 | 72 | 245.72 | 97 | 269.48 |
23 | 39.62 | 48 | 72.24 | 73 | 245.72 | 98 | 269.48 |
24 | 49.89 | 49 | 72.24 | 74 | 245.72 | 99 | 269.48 |
25 | 49.89 | 50 | 72.24 | 75 | 247.69 | 100 | 269.56 |
由图9可见,随着遗传算法遗传代数的增加,每代的最大适应度值不断增大,表明在连续搜索最优点的过程。在60~80代之间,图形波动幅度较大,说明通过交叉或变异,找到了比本代最优点优异的多的另外一点。当代数大于80时,适应度值几乎不变。通过观察相关输出数据,发现很符合所求目标和约束条件。
4.3.2.2每个产品的8个变量值,即0/m型指定平台设计结果(m=3),如下表:
控制基因值 | 2 | 3 | 1 | 0 | 0 | 2 | 0 | 3 |
电机号 | 变量 | |||||||
Nc | Ns | Awa(m2) | Awf(m2) | r(m) | t(m) | I(A) | l(m) | |
1 | 773 | 88 | 0.00000023 | 0.00000024 | 0.02032 | 0.010537 | 3.162 | 0.015121 |
2 | 773 | 69 | 0.00000023 | 0.00000023 | 0.03046 | 0.006372 | 3.472 | 0.015121 |
3 | 1001 | 47 | 0.00000023 | 0.00000022 | 0.0306 | 0.010537 | 4.294 | 0.015121 |
4 | 773 | 88 | 0.00000023 | 0.0000002 | 0.03134 | 0.006372 | 3.712 | 0.015121 |
5 | 1001 | 88 | 0.00000023 | 0.00000021 | 0.0214 | 0.010537 | 3.866 | 0.028477 |
6 | 773 | 88 | 0.00000023 | 0.0000002 | 0.025 | 0.010537 | 4.26 | 0.028477 |
7 | 1001 | 69 | 0.00000023 | 0.00000022 | 0.02724 | 0.010537 | 4.338 | 0.028477 |
8 | 1001 | 88 | 0.00000023 | 0.0000002 | 0.02812 | 0.006372 | 3.76 | 0.028477 |
9 | 773 | 88 | 0.00000023 | 0.00000022 | 0.02748 | 0.010537 | 4.996 | 0.028477 |
10 | 1001 | 88 | 0.00000023 | 0.00000021 | 0.02554 | 0.006372 | 4.786 | 0.028477 |
4.3.2.3每个产品的6个约束值,即0/m型指定平台性能参数(m=3),如下表:
电机号 | 性能参数 | |||||
T(Nm) | P(W) | H(A turns/m) | M(kg) | η | r/t | |
1 | 0.0502375 | 311.20279 | 2783.822754 | 0.2821775 | 0.8558226 | 1.9284431 |
2 | 0.0991432 | 293.06857 | 1556.871216 | 0.569061 | 0.7339926 | 4.7799918 |
3 | 0.1234558 | 326.22147 | 1320.710205 | 0.5842041 | 0.6606214 | 2.9040531 |
4 | 0.1478007 | 294.80725 | 2062.190186 | 0.600402 | 0.6906092 | 4.9180873 |
5 | 0.2030594 | 309.60162 | 3224.489258 | 0.5483173 | 0.6963756 | 2.030939 |
6 | 0.2494355 | 336.47418 | 3022.804443 | 0.641996 | 0.6868221 | 2.3725922 |
7 | 0.2985629 | 300.0416 | 2208.456787 | 0.7911885 | 0.6014425 | 2.5851765 |
8 | 0.33795 | 249.77809 | 2332.661865 | 0.8745927 | 0.5776552 | 4.4127829 |
9 | 0.3948761 | 352.43036 | 3214.542725 | 0.7546485 | 0.6134131 | 2.6079534 |
10 | 0.4997971 | 279.23337 | 3275.501709 | 0.7564876 | 0.5073373 | 4.0079117 |
由上表可见,强约束扭矩T非常的符合,弱约束P值分布在300附近,而且偏差幅度在允许的范围内,其他约束都很好的满足了条件,因此该结果令人满意,可用于进一步的结果分析。
4.4 0/m型2层结构GA用于多平台产品族优化(不指定平台变量)
PFPF在优化过程中确定哪个变量适合作平台变量及其最佳类数。
数学模型为;
求解:x={Nc,Ns,Awa,Awf,r,t,,I,l},xcc={xcc,j}-
xcc为平台变量的控制基因值
目标:Min f1(x)=∑Mi (尽量以0.5kg为目标)
Max f2(x)=∑ηi
Min f3(x)=PFPF
约束:Ti(x)={0.05,0.10,0.125,0.15,0.20,0.25,0.30,0.35,0.40,0.50}Nm
Pi(x)=300W
Hi(x)<5000A·turns/m
Mi(x)≤2.0kg
ηi(x)≥15%
ri/ti>1
xcc,j={0,1,2,3}
其中:i=1,2,…,10.
4.4.1适应度函数
Fitness=3.5/Cstr[0]+3.5*(∑Mi/∑ηi)/(1+∑Mi/∑ηi)-PFPF
-(310*Cstr[1]+Cstr[2]+200*Cstr[3]+Cstr[4]+Cstr[5])/5
4.4.2输出结果(运行100代)
4.4.2.1每代最大适应度值,即0/m型不指定平台每代最大适应度值(m=3),如下表:
代数 | 最大适应度值 | 代数 | 最大适应度值 | 代数 | 最大适应度值 | 代数 | 最大适应度值 |
1 | 9.11 | 26 | 12.89 | 51 | 16.53 | 76 | 31.48 |
2 | 9.11 | 27 | 13.32 | 52 | 16.53 | 77 | 32.05 |
3 | 9.95 | 28 | 13.34 | 53 | 18.45 | 78 | 40.91 |
4 | 10.36 | 29 | 13.34 | 54 | 18.45 | 79 | 40.91 |
5 | 10.47 | 30 | 13.34 | 55 | 18.45 | 80 | 40.91 |
6 | 10.47 | 31 | 14.18 | 56 | 18.45 | 81 | 41.17 |
7 | 10.55 | 32 | 14.18 | 57 | 18.45 | 82 | 46.17 |
8 | 10.55 | 33 | 14.18 | 58 | 18.45 | 83 | 46.17 |
9 | 10.55 | 34 | 14.18 | 59 | 18.45 | 84 | 46.17 |
10 | 10.55 | 35 | 14.61 | 60 | 18.45 | 85 | 49.77 |
11 | 11.39 | 36 | 14.82 | 61 | 18.45 | 86 | 62.69 |
12 | 11.39 | 37 | 14.82 | 62 | 18.45 | 87 | 62.69 |
13 | 11.39 | 38 | 14.82 | 63 | 19.07 | 88 | 62.69 |
14 | 11.39 | 39 | 16.21 | 64 | 19.07 | 89 | 78.04 |
15 | 11.39 | 40 | 16.21 | 65 | 19.07 | 90 | 87.88 |
16 | 11.39 | 41 | 16.21 | 66 | 19.35 | 91 | 98.76 |
17 | 11.39 | 42 | 16.21 | 67 | 19.35 | 92 | 116.45 |
18 | 11.39 | 43 | 16.21 | 68 | 19.35 | 93 | 116.45 |
19 | 11.91 | 44 | 16.21 | 69 | 20.86 | 94 | 116.45 |
20 | 12.24 | 45 | 16.44 | 70 | 20.86 | 95 | 117.95 |
21 | 12.24 | 46 | 16.44 | 71 | 20.86 | 96 | 117.95 |
22 | 12.24 | 47 | 16.52 | 72 | 26.79 | 97 | 118.66 |
23 | 12.25 | 48 | 16.52 | 73 | 26.79 | 98 | 126.00 |
24 | 12.26 | 49 | 16.52 | 74 | 26.79 | 99 | 126.00 |
25 | 12.89 | 50 | 16.53 | 75 | 31.48 | 100 | 157.28 |
由图10看出,每代最大适应度值在不断增大,反映了本发明不断搜索最优点的过程。虽然到100代时,最大适应度值仍有大幅度增加的趋势,但此时通过观察输出结果,发现各项指标都能很好的满足,因此这种情况下,可以采取这个较优结果,来与其他方法所得结果进行比较分析。
4.4.2.2每个产品的8个变量值,即0/m型不指定平台设计结果(m=3),如下表:
控制基因值 | 3 | 3 | 3 | 3 | 1 | 3 | 3 | 3 |
电机号 | 变量 | |||||||
Nc | Ns | Awa(m3) | Awf(m2) | r(m) | t(m) | I(A) | l(m) | |
1 | 903 | 52 | 0.0000002 | 0.00000023 | 0.0251 | 0.010489 | 3.4914 | 0.01207 |
2 | 903 | 52 | 0.0000002 | 0.00000024 | 0.0251 | 0.008675 | 3.4914 | 0.021362 |
3 | 1008 | 52 | 0.00000027 | 0.00000021 | 0.0251 | 0.010489 | 3.955 | 0.021362 |
4 | 762 | 89 | 0.00000027 | 0.00000024 | 0.0251 | 0.006278 | 3.4914 | 0.021362 |
5 | 903 | 69 | 0.00000027 | 0.00000024 | 0.0251 | 0.006278 | 3.4914 | 0.032673 |
6 | 1008 | 89 | 0.00000027 | 0.00000024 | 0.0251 | 0.010489 | 3.4914 | 0.032673 |
7 | 1008 | 89 | 0.00000027 | 0.00000024 | 0.0251 | 0.006278 | 3.4914 | 0.032673 |
8 | 903 | 89 | 0.00000027 | 0.00000021 | 0.0251 | 0.006278 | 3.955 | 0.032673 |
9 | 1008 | 89 | 0.00000027 | 0.00000021 | 0.0251 | 0.010489 | 4.4642 | 0.032673 |
10 | 1008 | 89 | 0.00000027 | 0.00000024 | 0.0251 | 0.006278 | 4.4642 | 0.032673 |
4.4.2.3每个产品的6个约束值,即0/m型不指定平台性能参数(m=3),如下表:
电机号 | 性能参数 | |||||
T(Nm) | P(W) | H(A turns/m) | M(kg) | η | r/t | |
1 | 0.04946928 | 312.4021606 | 457.6501465 | 0.3279295 | 0.77806324 | 2.39305698 |
2 | 0.09619252 | 286.4941406 | 1448.59668 | 0.51500428 | 0.71353722 | 2.89337124 |
3 | 0.12540966 | 337.3135376 | 1651.195557 | 0.58277684 | 0.74163371 | 2.39305698 |
4 | 0.1512866 | 309.930481 | 2459.140869 | 0.57718325 | 0.77190733 | 3.99822482 |
5 | 0.21258827 | 283.5084229 | 1906.524902 | 0.83158255 | 0.70610106 | 3.99822482 |
6 | 0.2558434 | 283.1493225 | 2494.824219 | 0.83587229 | 0.70520669 | 2.39305698 |
7 | 0.30609271 | 267.7033081 | 2459.140869 | 0.87890756 | 0.66673708 | 3.99822482 |
8 | 0.35186058 | 294.5454712 | 2785.663086 | 0.83730495 | 0.64760178 | 3.99822482 |
9 | 0.41828144 | 318.2096252 | 3189.974121 | 0.83027625 | 0.61982197 | 2.39305698 |
10 | 0.5004347 | 297.1102295 | 3144.3479 | 0.87890756 | 0.57872373 | 3.99822482 |
强约束扭矩T非常符合要求,弱约束P也在误差允许的范围内,其他约束条件都能很好的满足,因此这组数据比较理想。
5结果比较与分析
多平台与单平台比较
5.1 0/m型(m=3)指定平台与0/1型指定平台:
电机号 | 0/1型(指定平台变量) | 0/m型m=3(指定平台变量) | 目标比较(%) | |||
η(%) | M(kg) | η(%) | M(kg) | η | M | |
1 | 85.033447 | 0.33 | 77.81 | 0.28 | -8.49 | -0.15 |
2 | 80.098069 | 0.47 | 71.35 | 0.57 | -10.92 | 0.21 |
3 | 76.372635 | 0.66 | 74.16 | 0.58 | -2.90 | -0.11 |
4 | 74.101222 | 0.71 | 77.19 | 0.60 | 4.17 | -0.15 |
5 | 68.803394 | 0.85 | 70.61 | 0.55 | 2.63 | -0.36 |
6 | 65.499246 | 0.75 | 70.52 | 0.64 | 7.67 | -0.14 |
7 | 61.787641 | 0.85 | 66.67 | 0.79 | 7.90 | -0.07 |
8 | 58.991009 | 0.85 | 64.76 | 0.87 | 9.78 | 0.03 |
9 | 55.557936 | 0.89 | 61.98 | 0.75 | 11.56 | -0.15 |
10 | 49.600437 | 1.04 | 57.87 | 0.76 | 16.67 | -0.27 |
平台变量 | 11100101 | 23100203 | 3.81 | -0.12 | ||
平台变量个数 | 5 | 5 | ||||
平均变化百分比 |
对比0/m型(m=3)指定平台与0/1型指定平台两种情况发现:
●前者产品族产品平均效率比后者提高3.81%,平均质量减小了0.12kg;
●前者的产品平台变量与后者相同,都是5个。
这两种情况都是预先指定了平台变量,然后用面向大规模定制的产品族设计方法进行优化。两者平台变量都是5个,虽然后者是单一平台,通用性更好,但前者平均效率比其提高了3.81%,而且平均质量也减少了0.12kg,因而前者的性能大大好于后者。相比之下,后者牺牲了很大的性能只换来通用性上的一点优越性,因此从总体上看,0/m型(m=3)指定平台这种情况在设计产品族问题上比0/1型指定平台要好。
5.2 0/m型(m=3)不指定平台与0/1型不指定平台:
电机号 | 0/1型(不指定平台变量) | 0/m型m=3(不指定平台变量) | 目标比较(%) | |||
η(%) | M(kg) | η(%) | M(kg) | η | M | |
1 | 75.72 | 0.39 | 77.81 | 0.33 | 2.76 | -0.16 |
2 | 69.81 | 0.55 | 71.35 | 0.52 | 2.21 | -0.06 |
3 | 70.21 | 0.52 | 74.16 | 0.58 | 5.63 | 0.12 |
4 | 68.92 | 0.66 | 77.19 | 0.58 | 11.99 | -0.12 |
5 | 68.19 | 0.54 | 70.61 | 0.83 | 3.56 | 0.53 |
6 | 62.49 | 0.83 | 70.52 | 0.84 | 12.86 | 0.01 |
7 | 58.36 | 0.60 | 66.67 | 0.88 | 14.24 | 0.47 |
8 | 59.81 | 0.58 | 64.76 | 0.84 | 8.28 | 0.45 |
9 | 55.16 | 0.68 | 61.98 | 0.83 | 12.37 | 0.21 |
10 | 50.81 | 0.88 | 57.87 | 0.88 | 13.90 | 0.00 |
平台变量 | 00101000 | 33331333 | 8.78 | 0.15 | ||
平台变量个数 | 2 | 8 | ||||
平均变化百分比 |
对比0/m型(m=3)不指定平台与0/1型不指定平台两种情况发现:
●前者产品族产品平均效率比后者提高8.78%,平均质量增加了0.15kg;
●前者的产品平台变量有8个,而后者只有2个。
从目标指标上看,前者虽然平均质量增加了0.15kg,但平均效率却大幅提高了8.78%,因此从总体上看,前者的性能比后者有一定的提高。而且,前者平台变量数大大多于后者,产品族共性非常高,这在大规模定制中会有非常广泛的应用。综合考虑,0/m型(m=3)不指定平台的情况效果大大好于0/1型不指定平台的情况。
5.3 0/m型(m=3)不指定平台变量与0/m型(m=3)指定平台变量:
电机号 | 0/m型m=3(指定平台变量) | 0/m型m=3(不指定亚台变量) | 目标比较(%) | |||
η(%) | M(kg) | η(%) | M(kg) | η | M | |
1 | 85.58226 | 0.28 | 77.81 | 0.33 | -9.09 | 0.16 |
2 | 73.39926 | 0.57 | 71.35 | 0.52 | -2.79 | -0.09 |
3 | 66.06214 | 0.58 | 74.16 | 0.58 | 12.26 | 0.00 |
4 | 69.06092 | 0.60 | 77.19 | 0.58 | 11.77 | -0.04 |
5 | 69.63756 | 0.55 | 70.61 | 0.83 | 1.40 | 0.52 |
6 | 68.68221 | 0.64 | 70.52 | 0.84 | 2.68 | 0.30 |
7 | 60.14425 | 0.79 | 66.67 | 0.88 | 10.86 | 0.11 |
8 | 57.76552 | 0.87 | 64.76 | 0.84 | 12.11 | -0.04 |
9 | 61.34131 | 0.75 | 61.98 | 0.83 | 1.04 | 0.10 |
10 | 50.73373 | 0.76 | 57.87 | 0.88 | 14.07 | 0.16 |
平台变量 | 23100203 | 33331333 | 5.43 | 0.12 | ||
平台变量个数 | 5 | 8 | ||||
平均变化百分比 |
对比0/m型(m=3)不指定平台变量与0/m型(m=3)指定平台变量两种情况发现:
●前者产品族产品平均效率比后者提高5.43%,平均质量增加了0.12kg;
●前者的产品平台变量为8个,而后者有5个。
两种方法区别在于前者的平台变量及其所分类数都是在遗传算法的优化过程中不断寻找最优解所得,而后者则是先指定平台变量,然后运用改进的遗传算法确定所有变量的值。两者对比,前者虽然平均质量增加了0.12kg,但平均效率却提高了5.43%。综合看来,两者总体性能相差不大,但前者的共性比后者好很多,所有变量都可以作为平台变量,更适合于大规模定制。因此,0/m型(m=3)不指定平台变量的情况效果更好。
通过以上的结果对比和分析,可以定量的总结出在产品族设计问题上,运用本发明方法相比于其他设计方法效果更好,特别是0/m型两层结构方法,在不指定平台变量的情况下,能自动寻找出最优的产品平台变量个数及其分类方法,并且确定所有变量的值。不仅在整体性能上有了一定的提高,而且能找到最佳平台。
Claims (4)
1、一种实现面向大规模定制的产品族设计的方法,其特征在于,步骤如下:
第一步:确定决策变量及其各种约束条件,即确定出单个产品的表现型X和问题的解空间;一个产品族通过一组设计变量表示如下:
x={x1,x2,…,xm}
式中xj为设计变量列向量,包括平台变量与个性化变量;j为设计变量;m为描述产品族的设计变量的总数量;
对于产品变量i=1,2,…,n
xj,i={xj,1,xj,2,…,xj,n}′
式中xj,i为产品变量i的一个设计变量j的值;i为产品变量;n为产品族内产品总数量;
第二步:根据第一步确定的决策变量和约束条件建立优化模型,确定出目标函数的类型及其数学描述形式或量化方法,也就是最终获取的产品族;建立优化模型用公式描述如下:
For i=1,2,…,n
Find:xc,xv,i
Minimize:f(xc,xv,i)
Subject to:g(xc,xv,i)≤0
h(xc,xv,i)=0
xl≤x≤xu
式中c为平台设计变量的数量,v为个性化设计变量的数量,xc为平台变量的值,xv,i为个性化设计变量的值;
第三步:建立优化模型之后采用实数编码表示可行解的染色体编码,确定出单个产品的基因型X及遗传算法的搜索空间;
第四步:在可行解的编码方式确定后选择相应的解码方法,作为由单个产品基因型X到单个产品表现型X的对应关系或转换方法;
第五步:确定单个产品适应度的量化评价方法,由目标函数值f(X)到单个产品适应度F(X)的转换规则;适应度函数就是目标函数:
f(x)为适应度值;g(x)为最大值问题下的适应度值;
第六步:设计遗传算子,确定出选择、交叉、变异遗传算子的具体操作方法;
第七步:根据第六步所确定的选择、交叉、变异具体操作方法,进行具体操作,进而确定产品的全部设计变量x的值,最终获得产品族。
2、根据权利要求1所述的实现面向大规模定制的产品族设计的方法,其特征是,所述第四步中,进行产品族设计的时候,当实现多个目标时,即为多目标优化问题,多目标最优化的数学描述如下:
其中f(x)为适应度值,x为要设计的产品族,x={x1,x2,…,xm},xj为设计变量列向量,包括平台变量与个性化变量,X指产品族设计的可行方案,Rm为所有产品族设计方案。
3、根据权利要求1所述的实现面向大规模定制的产品族设计的方法,其特征是,所述第六步中,交叉方法是指:首先对群体中的个体进行随机配对,然后在配对个体中随机设定交叉点,配对个体彼此交换部分基因。
4、根据权利要求1所述的实现面向大规模定制的产品族设计的方法,其特征是,所述第六步,在产品族优化设计过程中,对于产品的某个变量,如果允许变异,则在所规定的范围内,再随机产生一个赋值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100396225A CN100492376C (zh) | 2007-04-19 | 2007-04-19 | 实现面向大规模定制的产品族设计的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100396225A CN100492376C (zh) | 2007-04-19 | 2007-04-19 | 实现面向大规模定制的产品族设计的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101042712A true CN101042712A (zh) | 2007-09-26 |
CN100492376C CN100492376C (zh) | 2009-05-27 |
Family
ID=38808217
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007100396225A Expired - Fee Related CN100492376C (zh) | 2007-04-19 | 2007-04-19 | 实现面向大规模定制的产品族设计的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100492376C (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101276436B (zh) * | 2008-04-25 | 2011-08-17 | 浙江大学 | 一种利用基于关联规则的遗传算法进行仪表设计的方法 |
CN103295063A (zh) * | 2013-05-17 | 2013-09-11 | 浙江大学 | 基于遗传算法对suv产品族基因进化的方法 |
CN105825291A (zh) * | 2016-03-09 | 2016-08-03 | 同济大学 | 一种基于食物链生态进化算法的产品优化更新方法 |
CN105844398A (zh) * | 2016-03-22 | 2016-08-10 | 武汉大学 | 一种基于plm数据库面向dpipp产品族的挖掘算法 |
CN106295004A (zh) * | 2016-08-12 | 2017-01-04 | 江苏大学 | 一种考虑扰动型设计变量区间的永磁电机优化设计方法 |
CN106528964A (zh) * | 2016-10-21 | 2017-03-22 | 浙江大学城市学院 | 一种面向用户偏好的饮料瓶产品智能设计方法 |
-
2007
- 2007-04-19 CN CNB2007100396225A patent/CN100492376C/zh not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101276436B (zh) * | 2008-04-25 | 2011-08-17 | 浙江大学 | 一种利用基于关联规则的遗传算法进行仪表设计的方法 |
CN103295063A (zh) * | 2013-05-17 | 2013-09-11 | 浙江大学 | 基于遗传算法对suv产品族基因进化的方法 |
CN103295063B (zh) * | 2013-05-17 | 2015-09-30 | 浙江大学 | 基于遗传算法对suv产品族基因进化的方法 |
CN105825291A (zh) * | 2016-03-09 | 2016-08-03 | 同济大学 | 一种基于食物链生态进化算法的产品优化更新方法 |
CN105844398A (zh) * | 2016-03-22 | 2016-08-10 | 武汉大学 | 一种基于plm数据库面向dpipp产品族的挖掘算法 |
CN106295004A (zh) * | 2016-08-12 | 2017-01-04 | 江苏大学 | 一种考虑扰动型设计变量区间的永磁电机优化设计方法 |
CN106295004B (zh) * | 2016-08-12 | 2020-03-06 | 江苏大学 | 一种考虑扰动型设计变量区间的永磁电机优化设计方法 |
CN106528964A (zh) * | 2016-10-21 | 2017-03-22 | 浙江大学城市学院 | 一种面向用户偏好的饮料瓶产品智能设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100492376C (zh) | 2009-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1253790C (zh) | 指令调度方法和指令调度设备 | |
CN101042712A (zh) | 实现面向大规模定制的产品族设计的方法 | |
CN1237502C (zh) | 生成声音模型的方法、装置和生成声音模型的计算机程序 | |
CN1845021A (zh) | 指令生成装置 | |
CN1934570A (zh) | 文本挖掘装置和其方法以及程序 | |
CN1292081C (zh) | 热处理装置、热处理方法及记录热处理程序的记录介质 | |
CN1695138A (zh) | 从mpc(模型预测控制器)的模型中移除pid动态值的方法 | |
CN1331449A (zh) | 用于将粘着法构成的文本或文档分段成词的字符串划分或区分的方法及相关系统 | |
CN1842702A (zh) | 声音合成装置和声音合成方法 | |
CN1265355C (zh) | 音源矢量生成装置及语音编码/解码装置 | |
CN1063381A (zh) | 电压或无功控制方法及其控制设备 | |
CN101046682A (zh) | 一种预测热轧含Nb带钢组织及力学性能的方法 | |
CN1447261A (zh) | 特定要素、字符串向量生成及相似性计算的装置、方法 | |
CN1873625A (zh) | 一种基于功能覆盖率的随机激励的自动产生方法 | |
CN1776621A (zh) | 程序变换方法 | |
CN1755663A (zh) | 信息处理设备、信息处理方法和程序 | |
CN1718802A (zh) | 镍基高温合金成分的优化设计方法 | |
CN1319275A (zh) | 用于控制电动机的设备 | |
CN1121447A (zh) | 轧机中控制形状的方法 | |
CN1889039A (zh) | 基于xml实现通用浏览器工具条的方法 | |
CN1624698A (zh) | 高阶合成方法和高阶合成装置 | |
CN1595752A (zh) | 大型发电机定子绕组内部故障主保护配置优选方法 | |
Ueda et al. | Calculating four-loop massless propagators with Forcer | |
CN1920868A (zh) | 信息处理装置及信息处理方法 | |
CN101046867A (zh) | 工作流决定方法以及工作流决定系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090527 Termination date: 20120419 |