《1 引言》
1 引言
概念设计过程是求解实现功能、满足各项技术和经济指标并最终确定综合最优方案的过程。现有的概念设计推理模型基本上是基于公理化设计理念的域与域间的映射机制。域的基本划分为顾客需求域、功能域、载体域和过程域
上述推理模型试图找到设计过程中隐含的逻辑关系, 属于自上而下的设计思路。这类模型和推理策略的不足之处在于:a. 模型的实现以功能、行为和结构等的完整存在集为基础, 但在概念设计阶段, 功能、行为、结构集定义不良, 因而使基于传统逻辑推理的求解方法遇到了困难;b. 推理过程中节点的匹配, 只能以当前节点的情况进行判断, 仅仅是一种局部的评价模式。节点所选元素对全局 (产品整体方案) 的影响, 尤其是当前节点与其他节点结合起来对全局的影响则无从评价。
可拓学中的菱形思维方法和物元变换技术, 既符合概念设计中人们的思维模式又具有较大的创造性
《2 产品及其组成元素物元模型》
2 产品及其组成元素物元模型
可拓学以物元作为认识和分析事物的基本逻辑单元, 物元由描述事物的概念N、特征c和关于特征的量值v组成, 记为
《2.1元素的物元模型》
2.1元素的物元模型
元素是组成产品方案的物质基础, 有关它们的操作——对它们的变换和组合, 构成了产品设计的过程。元素的物元模型除了包括描述自身的特征如材料、功能、参数大小等外, 还应该包括与其他元素物元的关系特征如它的接口特征、实际连接情况等。假定根据设计对象的领域知识确定了设计中元素应由c1, c2, …, cw等特征描述, 则它们的物元形式为
式中 vj (j = 1, 2, …, n) 是事物Rp 关于特征cj 的量值。元素属性之间的关系用物元特征之间的相关性来描述。
《2.2产品方案的物元模型》
2.2产品方案的物元模型
产品设计方案可以看成是由元素组成的有机系统, 具有外部和内部特征。内部特征描述方案的内部组成, 外部特征描述方案的对外整体特性。一旦确定了组成方案的元素及其特征, 就相应地确定了方案的外部和内部特征。
设计方案的外部特征由用户需求特征描述。若用户需求特征包含r1, r2, …, rm, 则方案的物元表示为,
式中 u1, u2, …, um 分别是r1, r2, …, rm 的量值。用户的需求可以看成是标准方案物元
其中u′1, u′2, …, u′m分别是用户对产品外部特征r1, r2, …, rm的需求值。
方案的内部特征用组成它的元素特征矩阵来描述。如果组成产品的元素为p1, p2, …, pw, 则方案的内部特征模型为
式中, vij是第i (1, 2, …, w) 个元素的第j (1, 2, …, n) 个特征值。
方案外部模型中的特征值, 由内部模型中元素的某些特征的值来确定, 即存在一个转换因子集合{Δ1, Δ2, …, Δz}, 将它们依次作用于产品的内部模型, 就可以得到方案的外部模型
方案的内部模型一旦确定, 其外部模型将相应地确定。如果将方案的内部模型物元的特征看成是系统的输入, 方案的外部模型特征是其内部模型元素相互作用而对外的输出, 物元之间的相互作用就是系统的有机构成。即系统由一系列规则知识或称方法组成, 它们将输入转换成输出。转换算法是利用领域知识直接将内部特征转换成外部特征。
应该指出, 在上面产品的内部物元模型中, 通过它们的连接特征联系起来, 每个产品方案对应一个由其组成元素物元连接起来的连接图, 用邻接矩阵表示为
若pi与ph相连接, 则lih = 1, 否则lih = 0 (i, h = 1, 2, …, w) 。
《3 物元变换与进化推理》
3 物元变换与进化推理
《3.1方案编码》
3.1方案编码
在概念设计推理过程中, 编码代表了产品的一个方案, 遗传算法的编码方案采用前面讨论过的方案模型如式 (5) 和式 (7) 所示。一个编码包含两部分:式 (5) 是物元集合的特征矩阵, 特征被赋予具体的量值;式 (7) 代表了物元集合中物元之间连接图的邻接矩阵。例如, 第k个染色体的编码为
式中 k = 1, 2, …, M 是遗传算法的群体规模, wk是第k个个体所含元素数目。
应该指出, 算法采用的编码方案的广义长度即矩阵的规模并不固定, 因为在概念设计阶段, 方案的组成结构还没有确定。如果采用固定长度的编码将会降低算法的创新性。在后面的论述中, 产品方案和个体具有相同的含义。
《3.2物元变换与进化算子》
3.2物元变换与进化算子
利用物元的分解变换生成设计对象的初始组成。如在减速器的方案生成中, 就是较大变速比物元被分成几个较小的变速比物元。而每个较小变速比物元又被聚分成能够实现该变速比的元素物元。
置换、扩缩、增加以及删减变换是解决设计对象初始组成和用户需求矛盾的有力工具。当初始组成中某一个元素物元的某特征不能满足要求时, 可以采用扩缩变换来改变该特征的值;如果该特征的值已经达到极限而不能扩大或收缩时, 可以采用置换的方法将该物元取代;如果置换变换和扩缩变换都不能解决矛盾, 就应该考虑采用增加变换来增加物元以协同解决该矛盾。
在推理过程中, 所有变换是逆向应用的。首先对初始组成进行一组物元变换, 再将变换后的结果同变换前的结果相比较, 如果有所改进就接受这一组变换;相反, 则抛弃该组变换而重新选取另一组变换, 直到结果满意为止。整个求解过程如图1所示。
《3.2.1 分解变换与交叉算子》
3.2.1 分解变换与交叉算子
分解变换包括聚分变换与组分变换。分解的目的是为了重新组合, 根据以上分解变换给出如下的交叉算子, 设w1, w2 为个体1和个体2所含物元数目, X
, (10)
, (11)
, (12)
, (13)
则交叉后产生的新个体的编码为
, (14)
, (15)
, (16)
。 (17)
《3.2.2 其他变换与变异算子》
3.2.2 其他变换与变异算子
变异操作的目的是在染色体中引入新的基因, 从而增加群体中基因物质的新颖性。依据物元的置换变换、增加变换、删减变换、数乘变换等, 设计置换变异、增加变异、删减变异以及扩缩变异。
1) 置换变换与置换变异
以置换变换为依据, 置换变异是指自元素库中随机选取一个元素的物元表示, 并让它替代所选个体中的一个物元, 新物元与其他物元的连接以及该物元的所有特征随机生成或根据领域知识确定。个体内部模型编码如式 (10) 和式 (12) 所示, 如用随机选定的元素物元p′替代个体内部模型中第i个物元pi, 则式 (10) 变为
式中l
2) 增加、删减变换
增加或删减变异指向所选个体中增加或删减一个物元, 该物元与个体中的其他物元的连接以及物元的特征值随机选取。以增加为例, 假设个体编码如式 (10) 、式 (12) 所示, 并设增加物元为p′, 则式 (10) 相应变为
式中最后一行和最后一列为新增元素, 其值随机选取0或1, 描述了新增加物元p′ 与原方案中w1个物元的连接情况。
对于式 (12) , 变换后的特征编码为
最后增加的一列即第 (w1 + 1) 列, 该列的元素依次取 p′ 的特征值。
3) 数乘变换与扩缩变异
扩缩变异指扩大或缩小个体中某一个物元的特征值。该变异操作将个体的特征矩阵中按变异概率随机选中的物元特征数乘预定范围的系数ψ, ψ ≤ 1, 为缩小变异, ψ> 1为扩大变异。
《3.3确定适应度函数》
3.3确定适应度函数
根据个体编码的结构, 适应值函数由两部分组成——方案内部关联性和方案外部关联性, 适应度即为
其中Din为方案内部关联性, Dout 为方案外部关联性, ωin, ωout为Din, Dout的权重。
适应值的计算过程如下:
1) 内部关联性Din包括方案内部模型中元素之间连接的正确性γ和元素物元的特征值与当前技术水平下该特征的极限值的关联性din。
连接正确性用邻接矩阵L中正确连接的元素数与实际连接元素数的比值定义。
在方案的内部模型中, 元素特征的取值vij 受到在当前技术水平下该特征能达到的极限值 v′ij的约束, 该约束用 vij 与 v′ij的关联值 din ij 表示。
其中ωin i, ω1, ω2 为其后一项的权重。
2) 外部关联性Dout 定义为外部模型与用户需求的关联性, 表征个体满足用户需求的程度, 利用关联函数来建立。
为进行定量计算, 首先将推理中涉及的定性指标利用模糊物元技术量值化
式中, ξλ为第 λ (λ = 1, …, m) 个用户需求指标的允许误差率。确定各指标对方案影响的权重系数 ωout λ 。最后计算总的外部关联性
这一节所涉及的权重均采用层次分析法确定, 总共分三层。
综合式 (21) 、式 (24) 和式 (26) 可得
《3.4遗传算法参数的选择》
3.4遗传算法参数的选择
群体规模30, 算法运行到150代终止。新群体的生成策略是:若新生成的个体优于当前代群体中的最差个体, 则取而代之, 否则放弃新生成的个体。
采用遗传算法模拟物元变换, 适当放大了变异概率。同时考虑到算法的收敛性, 变异概率随着算法的运行逐渐减少。交叉概率pc =0.6保持不变。变异概率 pm = 0.4 (1- gen /200) , 其中gen 为当前运行代数。
《4 实例求解》
4 实例求解
现以求解机械定比传动的减速器方案设计为例。三种用户需求见表1。
以下简要介绍需求1的模型建立和适应值函数计算, 算法的其他部分请参考论文第3节, 为了简化描述, 将方案的组成定义为传动组级。
Table 1 Requirements
《表1》
需求 | 输入转速 /r·min-1 | 输出转速 /r·min-1 | 传动功率 /kW | 中心距 /mm | 效率 |
1 | 960 | 78 | 4.3 | 500 | 0.85 |
2 | 1 500 | 78 | 4.3 | 500 | 0.85 |
3 | 960 | 78 | 4.3 | 1 000 | 0.85 |
注: 假定以上数据的允许误差为10 %。
用户要求的物元描述为
根据领域知识建立元素库, 包括减速器常用的元素。以直齿圆柱齿轮传动为例, 物元形式为
上式中功能参数可以根据概念设计的要求精度选用, 如只选齿轮齿数和模数。
设有某个体的内部模型为
应用领域知识可以确定该方案。传动比12.15, 最低传动功率 4.2 kW, 中心距 495 mm, 而传动效率 0.89。外部模型为
依据用户要求对各指标的允许范围, 利用式 (25) 计算得传动比、传动功率、中心距和效率的外部关联度分别为0.870, 0.767, 0.900, 0.530。
如果已经确定外部模型中的各元素的权重系数为 0.30, 0.25, 0.30, 0.15, 则个体的需求满足度为
Dout = 0.870×0.30 + 0.767 × 0.25 + 0.900 × 0.30 + 0.737 × 0.15 = 0.8 025。
通过计算个体内部关联性为0.8664。内部关联性和外部关联性的权重取 0.417 和 0.583, 个体的适应值为
f = 0.417 × 0.866 4 + 0.583 × 0.8 025 = 0.8 291。
三个可行方案列于表2。
Table 2 Reference schemes
《表2》
需求 | 传动形式 | 元素参数 | 元素参数 | 性能参数 |
1 | 链 | 22 | 66 | 15.87 |
圆柱齿轮 | 20 | 81 | 2.75 | |
功率/kW | 4.2 | 传动比=12.15 | ||
中心距 / mm | 495 | |||
2 | 链 | 22 | 66 | 3.175 |
圆柱齿轮 | 19 | 120 | 3.5 | |
功率 / kW | 6.15 | 传动比=18.95 | ||
中心距 / mm | 454 | |||
3 | V带 | 110 | 280 | 3根 |
链 | 18 | 94 | 15.87 | |
功率 / kW | 4.38 | 传动比=13.29 | ||
中心距 /mm | 1 001 |
《5 结论》
5 结论
将物元变换技术用于概念设计推理, 建立了求解概念设计的内部物元和外部模型。利用遗传算法模拟物元变换技术, 用关联函数建立适应值函数, 将物元理论和遗传算法有机地结合起来。算法应用于减速器的方案设计, 实验证明算法具有可行性。
无论对于概念设计还是物元理论的应用都是一种新的尝试。其中还有待进一步研究的问题, 如概念设计中的布局问题、知识复用问题以及利用遗传算法模拟物元变换时遗传算法运行的各种参数的确定策略的研究等。